Patent Publication Number: US-2022222393-A1

Title: Dynamic prediction method and system for initiation volume of debris flow slope source

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
     This application claims the benefit and priority of Chinese Patent Application No. 202110030645.X, titled Dynamic Prediction Method and System for Initiation Volume of Debris Flow Slope Source, filed on Jan. 11, 2021, the entire contents of which is hereby incorporated by reference herein in its entirety. 
     TECHNICAL FIELD 
     The present disclosure relates to the technical field of prevention and control of debris flow disasters, in particular to a dynamic prediction method and system for an initiation volume of a debris flow slope source. 
     BACKGROUND 
     In recent years, debris flows occurred frequently due to extreme weather (such as heavy rainfall and storms) and active earthquakes induced by global warming. This seriously affects the safety of life and property of people in mountainous areas, and restricts the sustainable development of the economy and society in China&#39;s mountainous areas. The accurate prediction of the initiation volume of the debris flow slope source directly affects the accuracy of early warning and forecasting of the debris flow and the reliability of prevention and control engineering design. The difficulty lies in the real-time prediction of changes in the initiation volume of the debris flow slope source over time based on rainfall data and the physical and mechanical mechanisms of debris flow slope source initiation. 
     Under the action of rainfall, the initiation process of the debris flow slope source is usually affected by seepage, erosion or the combined effect of seepage and erosion. Due to the complex and diverse influencing factors, the debris flow has the characteristics of suddenness, randomness and unpredictability, making it hard to accurately predict the initiation volume of the debris flow slope source. At present, most of the calculation methods for the initiation volume of the debris flow slope source at home and abroad are limited to empirical models and based on a large number of simplified methods. They are quite different from the actual initiation volume and initiation process characteristics of the debris flow slope source, resulting in large errors in the prediction of the debris flow initiation volume. 
     The existing prediction methods of the initiation volume of the debris flow slope source are essentially based on four single models, which have the following problems: 
     1. The calculation based on the stochastic process model takes into account the random distribution of rainfall, runoff erosion, rainfall infiltration, gravity and soil parameters, but it is mainly based on statistical formulas, lacking physical meaning and certainty, and the calculation results fluctuate greatly. 
     2. The calculation based on the infinite slope model divides the slope into several segments, and calculates the sliding force and anti-sliding force of the segments according to the limit equilibrium. However, it believes that the debris flow slope source is initiated at one time, and the initiation depth is consistent with the source thickness, which is inconsistent with the actual gradual initiation process of the debris flow slope source, and it does not consider runoff erosion, random distribution of soil parameters and rainfall infiltration. 
     3. The calculation based on the slope erosion model essentially considers the erosion effect of the slope runoff caused by rainfall on the debris flow slope source, without considering the rainfall infiltration, the gravity and the random distribution of soil parameters. 
     4. The calculation based on the soil column model discretizes the slope into soil column elements, taking into account the gravity-induced interaction between the soil columns. However, it believes that the debris flow slope source represented by each soil column is initiated at one time and the slope initiation depth is consistent with the source thickness, without considering the gradual initiation process of the slope source caused by factors such as runoff erosion and random distribution of soil parameters. 
     In addition to these four models, there are some empirical formulas and simplified methods for predicting the initiation volume of the debris flow slope source. However, they are subject to various restrictions, resulting in a large difference between the predicted and real initiation volumes of the debris flow slope source, and making it hard to meet the requirements of debris flow warning and prevention for the calculation accuracy. 
     SUMMARY 
     The present disclosure aims to provide a dynamic prediction method and system for an initiation volume of a debris flow slope source. The present disclosure can accurately predict a debris flow initiation volume. 
     To achieve the above objective, the present disclosure provides a solution as follows: 
     A prediction method for an initiation volume of a debris flow slope source, where the prediction method includes: 
     selecting a debris flow source slope to be predicted, and dividing the debris flow source slope to be predicted into regular hexagonal soil columns arranged in i rows and j columns; 
     calculating a most unfavorable sliding surface of each soil column according to an upper bound theorem of a limit analysis, and calculating an unbalanced force on the most unfavorable sliding surface; 
     determining whether each soil column is unstable according to the unbalanced force on the most unfavorable sliding surface; stopping the prediction of the soil column if a determination result indicates no, and if the determination result indicates yes, then 
     acquiring centers of gravity and elevations of a central soil column and six adjacent soil columns around, determining an instability direction of the central soil column, and determining a force mode of the six surrounding adjacent soil columns on the central soil column; 
     simulating the force mode of the six surrounding adjacent soil columns on the central soil column through a fiber bundle model, and determining a break status of connection bonds; 
     stopping the prediction of the soil column if a determination result indicates that connection bonds of a lateral tensile stress are not all broken, and if the determination result indicates that the connection bonds of the lateral tensile stress are all broken, then 
     determining whether the soil column is fluidized according to a fluidization criterion; stopping the prediction of the soil column if a determination result indicates no, and if the determination result indicates yes, predicting that the soil column is about to initiate a debris flow, and predicting an initiation volume of the debris flow slope source according to a fluidization status. 
     A prediction system for an initiation volume of a debris flow slope source, where the prediction system includes: 
     a slope division module, for selecting a debris flow source slope to be predicted, and dividing the debris flow source slope to be predicted into regular hexagonal soil columns arranged in i rows and j columns; 
     a sliding surface and unbalanced force calculation module, for calculating a most unfavorable sliding surface of each soil column according to an upper bound theorem of a limit analysis, and calculating an unbalanced force on the most unfavorable sliding surface; 
     an instability determination module, for determining whether each soil column is unstable according to the unbalanced force on the most unfavorable sliding surface; 
     stopping the prediction of the soil column if a determination result indicates no, and if the determination result indicates yes, connecting a force mode determination module; 
     the force mode determination module, for acquiring centers of gravity and elevations of a central soil column and six adjacent soil columns around, determining an instability direction of the central soil column, and determining a force mode of the six surrounding adjacent soil columns on the central soil column; 
     a connection bond break determination module, for simulating the force mode of the six surrounding adjacent soil columns on the central soil column through a fiber bundle model, and determining a break status of connection bonds; stopping the prediction of the soil column if a determination result indicates that connection bonds of a lateral tensile stress are not all broken, and if the determination result indicates that the connection bonds of the lateral tensile stress are all broken, connecting a fluidization determination module; and 
     the fluidization determination module, for determining whether the soil column is fluidized according to a fluidization criterion; stopping the prediction of the soil column if a determination result indicates no, and if the determination result indicates yes, predicting that the soil column is about to initiate a debris flow, and predicting an initiation volume of the debris flow slope source according to a fluidization status. 
     Based on the specific embodiments provided by the present disclosure, the present disclosure has the following technical effects: 
     By digitizing a theoretical calculation model for the debris flow initiation volume, the present disclosure can accurately predict the initiation flow of a potential debris slope flow, and solve the long-standing problems of difficulty and low accuracy in calculating the debris flow initiation volume. By simulating a failure and interaction process between the soil columns through a fiber bundle model, the present disclosure can determine whether the soil column is about to fail according to a break status of the fiber bundle before a debris flow is initiated. In this way, the present disclosure can perform early prediction and analysis on the initiation of the debris slope flow, facilitating the planning of disaster prevention measures in advance. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIG. 1  is a flowchart of a prediction method for an initiation volume of a debris flow slope source according to the present disclosure. 
         FIG. 2  shows numbering of soil columns. 
         FIG. 3  shows a positional relationship between the soil columns. 
         FIG. 4  shows stress transfer between the soil columns. 
         FIG. 5  shows a schematic diagram of a prediction system for an initiation volume of a debris flow slope source according to the present disclosure. 
         FIG. 6  shows a time-dependent curve of a cumulative initiation volume of a debris flow at a mouth of a valley. 
         FIG. 7  shows a time-dependent curve of a recorded initiation volume of the debris flow at the mouth of the valley. 
     
    
    
     DETAILED DESCRIPTION 
     The technical solutions of the embodiments of the present disclosure are clearly and completely described below with reference to the accompanying drawings. Apparently, the described embodiments are merely a part rather than all of the embodiments of the present disclosure. All other embodiments derived from the embodiments of the present disclosure by a person of ordinary skill in the art without creative efforts should fall within the protection scope of the present disclosure. 
     An objective of the present disclosure is to provide a prediction method and system for an initiation volume of a debris flow slope source. The present disclosure can accurately predict the initiation flow of a potential debris slope flow, and solve the long-standing problems of difficulty and low accuracy in calculating the debris slope flow initiation volume. 
     To make the above objectives, features and advantages of the present disclosure clearer and more comprehensible, the present disclosure is described in further detail below with reference to the accompanying drawings and specific implementations. 
       FIG. 1  is a flowchart of the prediction method for an initiation volume of a debris flow slope source according to the present disclosure. The prediction method includes: 
     Step  101 : Select a debris flow source slope to be predicted, and divide the debris flow source slope to be predicted into regular hexagonal soil columns arranged in i rows and j columns. 
     Step  102 : Calculate a most unfavorable sliding surface of each soil column according to an upper bound theorem of a limit analysis, and calculate an unbalanced force on the most unfavorable sliding surface. 
     Step  103 : Determine whether each soil column is unstable according to the unbalanced force on the most unfavorable sliding surface; stop the prediction of the soil column if a determination result indicates no, and if the determination result indicates yes, perform Step  104 . 
     Step  104 : Acquire centers of gravity and elevations of a central soil column and six adjacent soil columns around, determine an instability direction of the central soil column, and determine a force mode of the six surrounding adjacent soil columns on the central soil column. 
     Step  105 : Simulate the force mode of the six surrounding adjacent soil columns on the central soil column through a fiber bundle model, and determine a break status of connection bonds; stop the prediction of the soil column if a determination result indicates that connection bonds of a lateral tensile stress are not all broken, and if the determination result indicates that the connection bonds of the lateral tensile stress are all broken, perform Step  106 . 
     Step  106 : Determine whether the soil column is fluidized according to a fluidization criterion; stop the prediction of the soil column if a determination result indicates no, and if the determination result indicates yes, predict that the soil column is about to initiate a debris flow, and predict an initiation volume of the debris flow slope source according to a fluidization status. 
     When the determination results in Steps  103 ,  105  and  106  indicate no, the initiation volume of the debris flow slope source may be recorded as zero. 
     Step  102  specifically includes: 
     Calculate a depth H sd  and an angle α+β of the most unfavorable sliding surface according to the upper bound theorem of the limit analysis: 
     total external power: 
     
       
         
           
             
               
                 W 
                 . 
               
               ext 
             
             = 
             
               
                 
                   ( 
                   
                     
                       H 
                       sd 
                     
                     - 
                     
                       
                         1 
                         2 
                       
                       ⁢ 
                       
                         L 
                         g 
                       
                       ⁢ 
                       
                         tan 
                         ⁡ 
                         
                           ( 
                           
                             α 
                             + 
                             β 
                           
                           ) 
                         
                       
                     
                   
                   ) 
                 
                 ⁢ 
                 
                   L 
                   g 
                 
                 ⁢ 
                 
                   ρ 
                   r 
                 
                 ⁢ 
                 g 
                 ⁢ 
                 v 
                 ⁢ 
                 
                   sin 
                   ⁡ 
                   
                     ( 
                     
                       α 
                       + 
                       β 
                     
                     ) 
                   
                 
               
               - 
               
                 
                   τ 
                   T 
                 
                 ⁢ 
                 
                   
                     L 
                     g 
                   
                   
                     cos 
                     ⁡ 
                     
                       ( 
                       
                         α 
                         + 
                         β 
                       
                       ) 
                     
                   
                 
                 ⁢ 
                 
                   v 
                   . 
                 
               
             
           
         
       
     
     internal energy dissipation on the sliding surface: 
     
       
         
           
             
               
                 D 
                 . 
               
               int 
             
             = 
             
               c 
               ⁢ 
               
                 
                   L 
                   g 
                 
                 
                   cos 
                   ⁡ 
                   
                     ( 
                     
                       α 
                       + 
                       β 
                     
                     ) 
                   
                 
               
               ⁢ 
               v 
               ⁢ 
               
                   
               
               ⁢ 
               cos 
               ⁢ 
               
                   
               
               ⁢ 
               
                 φ 
                 . 
               
             
           
         
       
     
     safety factor of the soil column K: 
     
       
         
           
             
               K 
               = 
               
                 
                   
                     
                       D 
                       . 
                     
                     int 
                   
                   
                     
                       W 
                       . 
                     
                     ext 
                   
                 
                 = 
                 
                   c 
                   
                     
                       
                         ( 
                         
                           
                             H 
                             sd 
                           
                           - 
                           
                             
                               1 
                               2 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               tan 
                               ⁡ 
                               
                                 ( 
                                 
                                   α 
                                   + 
                                   β 
                                 
                                 ) 
                               
                             
                           
                         
                         ) 
                       
                       ⁢ 
                       
                         ρ 
                         r 
                       
                       ⁢ 
                       g 
                       ⁢ 
                       
                         sin 
                         ⁡ 
                         
                           ( 
                           
                             α 
                             + 
                             β 
                           
                           ) 
                         
                       
                       ⁢ 
                       
                         cos 
                         ⁡ 
                         
                           ( 
                           
                             α 
                             + 
                             β 
                           
                           ) 
                         
                       
                     
                     - 
                     
                       τ 
                       T 
                     
                   
                 
               
             
             . 
           
         
       
     
     corresponding minimum safety factor of the sliding surface: 
     
       
         
           
             
               { 
               
                 
                   
                     
                       
                         
                           
                             ∂ 
                             K 
                           
                           
                             ∂ 
                             η 
                           
                         
                         = 
                         0 
                       
                       , 
                       
                         η 
                         = 
                         
                           α 
                           + 
                           β 
                         
                       
                     
                   
                 
                 
                   
                     
                       
                         
                           ∂ 
                           K 
                         
                         
                           ∂ 
                           
                             H 
                             
                               s 
                               ⁢ 
                               d 
                             
                           
                         
                       
                       = 
                       0 
                     
                   
                 
               
             
             . 
           
         
       
     
     The calculating an unbalanced force on the most unfavorable sliding surface specifically includes: 
     Calculate a sliding force of the most unfavorable sliding surface according to 
         W   d   =H   sd [θρ w +(1−φ)ρ r ] g  sin(α+β)cos(α+β)−τ T +τ rf +γ w   h   s   J.  
 
     Calculate an anti-sliding force of the most unfavorable sliding surface according to 
     
       
         
           
             
               τ 
               s 
             
             = 
             
               { 
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 c 
                                 + 
                                 
                                   { 
                                   
                                     
                                       
                                         
                                           H 
                                           sd 
                                         
                                         ⁡ 
                                         
                                           [ 
                                           
                                             
                                               θρ 
                                               w 
                                             
                                             + 
                                             
                                               
                                                 ( 
                                                 
                                                   1 
                                                   - 
                                                   φ 
                                                 
                                                 ) 
                                               
                                               ⁢ 
                                               
                                                 ρ 
                                                 r 
                                               
                                             
                                           
                                           ] 
                                         
                                       
                                       ⁢ 
                                       g 
                                       ⁢ 
                                       
                                           
                                       
                                       ⁢ 
                                       
                                         
                                           cos 
                                           2 
                                         
                                         ⁡ 
                                         
                                           ( 
                                           
                                             α 
                                             + 
                                             β 
                                           
                                           ) 
                                         
                                       
                                     
                                     - 
                                   
                                 
                               
                             
                           
                           
                             
                               
                                 
                                   
                                     
                                       χ 
                                       ⁢ 
                                       
                                           
                                       
                                       ⁢ 
                                       h 
                                       ⁢ 
                                       
                                           
                                       
                                       ⁢ 
                                       
                                         ρ 
                                         w 
                                       
                                       ⁢ 
                                       g 
                                     
                                     ) 
                                   
                                   } 
                                 
                                 ⁢ 
                                 tan 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 γ 
                               
                             
                           
                         
                         , 
                       
                     
                     
                       
                         
                           
                             H 
                             1 
                           
                           - 
                           
                             H 
                             2 
                           
                           - 
                           
                             H 
                             sd 
                           
                         
                         &gt; 
                         
                           H 
                           d 
                         
                       
                     
                   
                   
                     
                       
                         
                           
                             
                               
                                 c 
                                 + 
                                 
                                   { 
                                   
                                     
                                       
                                         
                                           H 
                                           sd 
                                         
                                         ⁡ 
                                         
                                           [ 
                                           
                                             
                                               θ 
                                               ⁢ 
                                               
                                                   
                                               
                                               ⁢ 
                                               
                                                 ρ 
                                                 w 
                                               
                                             
                                             + 
                                             
                                               
                                                 ( 
                                                 
                                                   1 
                                                   - 
                                                   φ 
                                                 
                                                 ) 
                                               
                                               ⁢ 
                                               
                                                 ρ 
                                                 r 
                                               
                                             
                                           
                                           ] 
                                         
                                       
                                       ⁢ 
                                       g 
                                       ⁢ 
                                       
                                           
                                       
                                       ⁢ 
                                       
                                         
                                           cos 
                                           2 
                                         
                                         ⁡ 
                                         
                                           ( 
                                           
                                             α 
                                             + 
                                             β 
                                           
                                           ) 
                                         
                                       
                                     
                                     - 
                                   
                                 
                               
                             
                           
                           
                             
                               
                                 
                                   
                                     
                                       h 
                                       s 
                                     
                                     ⁢ 
                                     
                                       ρ 
                                       w 
                                     
                                     ⁢ 
                                     g 
                                     ⁢ 
                                     
                                         
                                     
                                     ⁢ 
                                     
                                       
                                         cos 
                                         2 
                                       
                                       ⁡ 
                                       
                                         ( 
                                         
                                           α 
                                           + 
                                           β 
                                         
                                         ) 
                                       
                                     
                                   
                                   } 
                                 
                                 ⁢ 
                                 tan 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 γ 
                               
                             
                           
                         
                         , 
                       
                     
                     
                       
                         
                           
                             H 
                             1 
                           
                           - 
                           
                             H 
                             2 
                           
                           - 
                           
                             H 
                             sd 
                           
                         
                         &lt; 
                         
                           H 
                           d 
                         
                       
                     
                   
                 
                 . 
               
             
           
         
       
     
     Calculate an unbalanced force on the most unfavorable sliding surface based on the sliding force and the anti-sliding force of the most unfavorable sliding surface. 
     
       
         
           
             
               τ 
               T 
             
             = 
             
               { 
               
                 
                   
                     
                       
                         
                           
                             ∑ 
                             
                               n 
                               = 
                               0 
                             
                             
                               n 
                               = 
                               5 
                             
                           
                           ⁢ 
                           
                             τ 
                             Tn 
                             ′ 
                           
                         
                         , 
                       
                     
                     
                       
                         t 
                         &gt; 
                         0 
                       
                     
                   
                   
                     
                       
                         0 
                         , 
                       
                     
                     
                       
                         t 
                         = 
                         0 
                       
                     
                   
                 
                 , 
               
             
           
         
       
     
     where τ T  represents a total stress of the six surrounding adjacent soil columns to the central soil column, and τ′ Tn  represents a stress of the surrounding adjacent soil column to the central soil column. 
     L g  represents a spacing between the soil columns; c represents a cohesion of a root-soil complex of a debris flow source to be predicted; H sd  represents a depth of the most unfavorable sliding surface; θ represents a soil moisture content of the debris flow source slope to be predicted; ρ w  represents a density of water; φ represents a porosity of the root-soil complex on the debris flow source slope to be predicted; ρ r  represents a density of the root-soil complex on the debris flow source slope to be predicted; g represents an acceleration of gravity; α represents an angle of a bedrock; β represents an angle between the most unfavorable sliding surface and the bedrock; χ represents a parameter related to a water-soil characteristic curve; h represents a capillary pressure head; γ represents an internal friction angle of the root-soil complex on the debris flow source slope to be predicted; τ rf  represents a runoff erosivity of the debris flow source slope to be predicted; γ w  represents a unit weight of the water; s represents a phreatic level head in an area where the debris flow source slope to be predicted is located; h s  represents a phreatic level hydraulic gradient in the area where the debris flow source slope to be predicted is located; H 1  represents a slope elevation of the central soil column; H 2  represents a bedrock elevation of the central soil column; H d  represents a groundwater level in the area where the debris flow source slope to be predicted is located. 
     Among the above parameters, the cohesion c of the root-soil complex of the debris flow source to be predicted and the internal friction angle γ of the root-soil complex on the debris flow source slope to be predicted obey a Weibull distribution: 
     
       
         
           
             
               
                 ( 
                 
                   c 
                   , 
                   γ 
                 
                 ) 
               
               ∼ 
               
                 f 
                 ⁡ 
                 
                   ( 
                   
                     
                       x 
                       
                         c 
                         , 
                         γ 
                       
                     
                     ; 
                     
                       λ 
                       
                         c 
                         , 
                         γ 
                       
                     
                     ; 
                     
                       k 
                       
                         c 
                         , 
                         γ 
                       
                     
                   
                   ) 
                 
               
             
             = 
             
               { 
               
                 
                   
                     
                       
                         
                           
                             k 
                             
                               c 
                               , 
                               γ 
                             
                           
                           
                             λ 
                             
                               c 
                               , 
                               γ 
                             
                           
                         
                         ⁢ 
                         
                           ( 
                           
                             
                               x 
                               
                                 c 
                                 , 
                                 γ 
                               
                             
                             
                               λ 
                               
                                 c 
                                 , 
                                 γ 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         
                           e 
                             
                           
                             - 
                             
                               
                                 ( 
                                 
                                   
                                     x 
                                     
                                       c 
                                       , 
                                       γ 
                                     
                                   
                                   / 
                                   
                                     λ 
                                     
                                       c 
                                       , 
                                       γ 
                                     
                                   
                                 
                                 ) 
                               
                               
                                 k 
                                 
                                   c 
                                   , 
                                   γ 
                                 
                               
                             
                           
                             
                             
                           
                             
                               k 
                               
                                 c 
                                 , 
                                 γ 
                               
                             
                             - 
                             1 
                           
                         
                       
                     
                     
                       
                         
                           x 
                           
                             c 
                             , 
                             γ 
                           
                         
                         ≧ 
                         0 
                       
                     
                   
                   
                     
                       0 
                     
                     
                       
                         
                           x 
                           
                             c 
                             , 
                             γ 
                           
                         
                         &lt; 
                         0 
                       
                     
                   
                 
                 . 
               
             
           
         
       
     
     The soil moisture content θ of the debris flow source slope to be predicted and the porosity φ and density ρ of the root-soil complex on the debris flow source slope to be predicted obey a normal distribution: (θ, φ, ρ)˜N(μ θ , μ φ , μ ρ , σ θ   2 , σ φ   2 , σ ρ   2 ). 
     μ θ , μ φ , μ ρ  respectively represent an expectation of the normal distribution obeyed by the moisture content, the porosity and the density; σ θ   2 , σ φ   2 , σ ρ   2  respectively represent a variance of the normal distribution of the moisture content, the porosity and the density; the expectation and the variance are determined according to a measured value of the slope. 
     The parameters of the root-soil complex on the debris flow source slope to be predicted, such as density ρ, porosity φ, cohesion c and internal friction angle γ, are measured by conducting an indoor geotechnical test of the root-soil complex on the selected debris flow source slope to be predicted. 
     The soil moisture content θ of the debris flow source slope to be predicted, the parameter λ related to a water-soil characteristic curve, the phreatic level head h s  and phreatic level hydraulic gradient J in the area where the debris flow source slope to be predicted is located are calculated based on related theories or empirical equations such as Darcy&#39;s law and Richard equation. 
     Step  104  specifically includes: 
     The tendency of the instability and failure movement of the soil column determines the direction of force distribution between the soil columns. According to the principle of least action, the movement tendency of the soil column may be determined according to the center of gravity of a failure soil column.  FIGS. 2 and 3  show numbering and positional relationship of the soil columns, where (i,j) indicate that the position of the soil column is in row i and column j. 
     As shown in  FIG. 2 , the surrounding soil columns are numbered. D 0,1 , D 1,2 , D 2,3 , D 3,4 , D 4,5  and D 5,6  represent six tendencies of instability and failure of the six soil columns. Taking D 5,6  as an example, if central soil column (i,j) is unstable along direction D 5,6 , surrounding soil columns (i,j+2) and (i+1,j+1) in front of soil column (i,j) along direction D 5,6  are subjected to a compressive stress of soil column (i,j), and other soil columns are subjected to a tensile stress of soil column (i,j). 
     Step  104  specifically includes: 
     Acquire a position (i,j) of the central soil column, and determine the center of gravity of the central soil column as (x i,j , y i,j , z i,j ). 
     Determine, according to the position coordinates of the soil columns and a geometric principle, the center of gravity of the surrounding adjacent soil column as (x k , y k , z k ) (k=0, 1, 2, 3, 4, 5,6): 
     
       
         
           
             { 
             
               
                 
                   
                     
                       
                         
                           x 
                           k 
                         
                         = 
                         
                           
                             
                               
                                 3 
                               
                               3 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     k 
                                     3 
                                   
                                   ⁢ 
                                   π 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             x 
                             ij 
                           
                         
                       
                       , 
                     
                   
                   
                     
                         
                     
                   
                 
                 
                   
                     
                       
                         
                           y 
                           k 
                         
                         = 
                         
                           
                             
                               
                                 3 
                               
                               3 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     k 
                                     3 
                                   
                                   ⁢ 
                                   π 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             y 
                             ij 
                           
                         
                       
                       , 
                     
                   
                   
                     
                         
                     
                   
                 
                 
                   
                     
                       
                         
                           z 
                           0 
                         
                         = 
                         
                           
                             z 
                             6 
                           
                           = 
                           
                             z 
                             
                               i 
                               , 
                               
                                 j 
                                 + 
                                 1 
                               
                             
                           
                         
                       
                       , 
                       
                         
                           z 
                           1 
                         
                         = 
                         
                           z 
                           
                             
                               i 
                               - 
                               1 
                             
                             , 
                             j 
                           
                         
                       
                       , 
                       
                         
                           z 
                           2 
                         
                         = 
                         
                           z 
                           
                             
                               i 
                               - 
                               1 
                             
                             , 
                             
                               j 
                               - 
                               1 
                             
                           
                         
                       
                       , 
                     
                   
                   
                     
                         
                     
                   
                 
                 
                   
                     
                       
                         
                           z 
                           3 
                         
                         = 
                         
                           z 
                           
                             i 
                             , 
                             
                               j 
                               - 
                               1 
                             
                           
                         
                       
                       , 
                       
                         
                           z 
                           4 
                         
                         = 
                         
                           z 
                           
                             
                               i 
                               + 
                               1 
                             
                             , 
                             
                               j 
                               - 
                               1 
                             
                           
                         
                       
                       , 
                       
                         
                           z 
                           5 
                         
                         = 
                         
                           z 
                           
                             
                               i 
                               + 
                               1 
                             
                             , 
                             j 
                           
                         
                       
                     
                   
                   
                     
                         
                     
                   
                 
               
               ⁢ 
               
                 
                   ( 
                   
                     
                       k 
                       = 
                       0 
                     
                     , 
                     1 
                     , 
                     2 
                     , 
                     3 
                     , 
                     4 
                     , 
                     5 
                     , 
                     6 
                   
                   ) 
                 
                 . 
               
             
           
         
       
     
     Acquire, according to the center-of-gravity coordinates of the central soil column and the surrounding adjacent soil column, a vector of the surrounding adjacent soil column relative to the central soil column, (x′ k , y′ k , z′ k ), 
     
       
         
           
             { 
             
               
                 
                   
                     
                       
                         
                           x 
                           k 
                           ′ 
                         
                         = 
                         
                           
                             
                               3 
                             
                             3 
                           
                           ⁢ 
                           
                             L 
                             g 
                           
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               
                                 
                                   k 
                                   3 
                                 
                                 ⁢ 
                                 π 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           y 
                           k 
                           ′ 
                         
                         = 
                         
                           
                             
                               3 
                             
                             3 
                           
                           ⁢ 
                           
                             L 
                             g 
                           
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               
                                 
                                   k 
                                   3 
                                 
                                 ⁢ 
                                 π 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           z 
                           k 
                           ′ 
                         
                         = 
                         
                           
                             z 
                             k 
                           
                           - 
                           
                             z 
                             
                               i 
                               , 
                               j 
                             
                           
                         
                       
                     
                   
                 
                 ⁢ 
                 
                   ( 
                   
                     
                       k 
                       = 
                       0 
                     
                     , 
                     1 
                     , 
                     2 
                     , 
                     3 
                     , 
                     4 
                     , 
                     5 
                     , 
                     6 
                   
                   ) 
                 
               
               , 
             
           
         
       
     
     where L g  represents a spacing between adjacent soil columns. 
     Determine a sum vector of two adjacent soil columns among the surrounding soil columns as (x′ k,k+1 , y′ k,k+1 , z′ k,k+1 ): 
     
       
         
           
             { 
             
               
                 
                   
                     
                       
                         
                           x 
                           
                             k 
                             , 
                             
                               k 
                               + 
                               1 
                             
                           
                           ′ 
                         
                         = 
                         
                           
                             
                               
                                 3 
                               
                               3 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     k 
                                     3 
                                   
                                   ⁢ 
                                   π 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             
                               
                                 3 
                               
                               3 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     k 
                                     3 
                                   
                                   ⁢ 
                                   π 
                                 
                                 ) 
                               
                             
                           
                         
                       
                       , 
                     
                   
                 
                 
                   
                     
                       
                         
                           y 
                           
                             k 
                             , 
                             
                               k 
                               + 
                               1 
                             
                           
                           ′ 
                         
                         = 
                         
                           
                             
                               
                                 3 
                               
                               3 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     k 
                                     3 
                                   
                                   ⁢ 
                                   π 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             
                               
                                 3 
                               
                               3 
                             
                             ⁢ 
                             
                               L 
                               g 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     k 
                                     3 
                                   
                                   ⁢ 
                                   π 
                                 
                                 ) 
                               
                             
                           
                         
                       
                       , 
                     
                   
                 
                 
                   
                     
                       
                         
                           z 
                           
                             k 
                             , 
                             
                               k 
                               + 
                               1 
                             
                           
                           ′ 
                         
                         = 
                         
                           
                             z 
                             k 
                           
                           + 
                           
                             z 
                             
                               k 
                               + 
                               1 
                             
                           
                         
                       
                       , 
                     
                   
                 
               
               ⁢ 
               
                 
                   ( 
                   
                     
                       k 
                       = 
                       0 
                     
                     , 
                     1 
                     , 
                     2 
                     , 
                     3 
                     , 
                     4 
                     , 
                     5 
                   
                   ) 
                 
                 . 
               
             
           
         
       
     
     Determine a projected vector of the sum vector in a two-dimensional plane as D k,k+1  (k=0, 1, 2, 3, 4, 5) 
     Determine a sum vector with a lowest center of gravity min(z′ k,k+1 ) (k=0, 1, 2, 3, 4, 5) as a movement direction vector of the soil column according to a principle of least action. 
     Calculate, by taking a value of k in the case of min(z′ k,k+1 ) (k=0, 1, 2, 3, 4, 5), the projected vector of the movement direction of the central soil column in the two-dimensional plane as (x″ i,j , y″ i,j ): x″ i,j =x′ k,k+1 , y″ i,j =y′ k,k+1 . 
     Determine a projected vector of a movement direction of the surrounding soil column in the two-dimensional plane as (x″ k , y″ k ). 
     Determine an angle θ k  between adjacent soil columns, 
     
       
         
           
             
               θ 
               k 
             
             = 
             
               arc 
               ⁢ 
               
                   
               
               ⁢ 
               cos 
               ⁢ 
               
                 
                   
                     
                       
                         x 
                         
                           i 
                           , 
                           j 
                         
                         ″ 
                       
                       * 
                       
                         x 
                         k 
                         ″ 
                       
                     
                     + 
                     
                       
                         y 
                         
                           i 
                           , 
                           j 
                         
                         ″ 
                       
                       * 
                       
                         y 
                         k 
                         ″ 
                       
                     
                   
                   
                     
                       
                         
                           
                             x 
                             
                               i 
                               , 
                               j 
                             
                             ″ 
                           
                           2 
                         
                         + 
                         
                           
                             y 
                             
                               i 
                               , 
                               j 
                             
                             ″ 
                           
                           2 
                         
                       
                     
                     ⁢ 
                     
                       
                         
                           
                             x 
                             k 
                             ″ 
                           
                           2 
                         
                         + 
                         
                           
                             y 
                             k 
                             ″ 
                           
                           2 
                         
                       
                     
                   
                 
                 . 
               
             
           
         
       
     
     Determine, when an angle between movement direction vectors of adjacent soil columns is an acute angle, that a force between the soil columns is a compressive stress; determine, when the angle between the movement direction vectors of the adjacent soil columns is an obtuse angle, that the force between the soil columns is a tensile stress. 
     If 
     
       
         
           
             
               
                 θ 
                 k 
               
               &lt; 
               
                 π 
                 2 
               
             
             , 
           
         
       
     
     then τ′ Tk =−τ Tk . 
     If 
     
       
         
           
             
               
                 θ 
                 k 
               
               = 
               
                 π 
                 2 
               
             
             , 
           
         
       
     
     then τ′ Tk =0. 
     If 
     
       
         
           
             
               
                 θ 
                 k 
               
               &gt; 
               
                 π 
                 2 
               
             
             , 
           
         
       
     
     then τ′ Tk =τ Tk . 
     τ Tk  represents a force exerted on the central soil column by a surrounding soil column that is unstable; 
     τ′ Tk  represents a force exerted on the central soil column by the surrounding soil column, which is used to calculate the sliding force. 
     In the fiber bundle model, a load between the most unfavorable sliding surfaces or contact surfaces of the soil column is: 
     
       
         
           
             
               F 
               = 
               
                 
                   
                     N 
                     I 
                   
                   ⁢ 
                   
                     σ 
                     F 
                   
                 
                 + 
                 
                   
                     
                       
                         f 
                         R 
                       
                       ⁢ 
                       
                         N 
                         F 
                       
                     
                     2 
                   
                   ⁢ 
                   
                     
                       
                         σ 
                         2 
                       
                       ⁢ 
                       F 
                     
                     
                       
                         σ 
                         max 
                       
                       ⁡ 
                       
                         ( 
                         θ 
                         ) 
                       
                     
                   
                 
               
             
             ; 
           
         
       
     
     a strength of the fiber bundle is: 
     
       
         
           
             
               
                 
                   σ 
                   max 
                 
                 ⁡ 
                 
                   ( 
                   θ 
                   ) 
                 
               
               = 
               
                 
                   
                     2 
                     ⁢ 
                     
                       ( 
                       
                         2 
                         - 
                         
                           f 
                           
                             R 
                             ) 
                           
                         
                       
                     
                   
                   
                     N 
                     F 
                   
                 
                 ⁢ 
                 
                   A 
                   c 
                 
                 ⁢ 
                 
                   
                     τ 
                     h 
                   
                   ⁡ 
                   
                     ( 
                     θ 
                     ) 
                   
                 
               
             
             ; 
           
         
       
     
     a count of intact fiber bundles is: 
     
       
         
           
             
               N 
               I 
             
             = 
             
               
                 
                   
                     
                       ( 
                       
                         1 
                         - 
                         
                           f 
                           R 
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       N 
                       F 
                     
                     ⁢ 
                     
                       
                         σ 
                         max 
                       
                       ⁡ 
                       
                         ( 
                         θ 
                         ) 
                       
                     
                   
                   + 
                   
                     
                       
                         N 
                         F 
                       
                       ⁢ 
                       
                         
                           
                             σ 
                             max 
                           
                           ⁡ 
                           
                             ( 
                             θ 
                             ) 
                           
                         
                         ⁡ 
                         
                           [ 
                           
                             
                               2 
                               ⁢ 
                               
                                 Ff 
                                 R 
                               
                             
                             - 
                             
                               4 
                               ⁢ 
                               F 
                             
                             + 
                             
                               
                                 N 
                                 F 
                               
                               ⁢ 
                               
                                 
                                   σ 
                                   max 
                                 
                                 ⁡ 
                                 
                                   ( 
                                   θ 
                                   ) 
                                 
                               
                             
                           
                           ] 
                         
                       
                     
                   
                 
                 
                   
                     ( 
                     
                       2 
                       - 
                       
                         f 
                         R 
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       σ 
                       max 
                     
                     ⁡ 
                     
                       ( 
                       θ 
                       ) 
                     
                   
                 
               
               . 
             
           
         
       
     
     N I  represents a count of intact fiber bundles; N F  represents a total count of fiber bundles on the contact surface; f R  represents a residual strength coefficient after the fiber bundle is broken (f R =residual tensile strength after the fiber bundle is broken/tensile strength before the fiber bundle is broken); σ F  represents a force exerted on each intact fiber bundle; σ max (θ) represents a maximum tensile strength that a fiber bundle can withstand, which is related to the moisture content (measured by experiment); A c  represents an area of the contact surface; τ h (θ) represents a tensile strength of the root-soil complex, which is related to the moisture content (measured by experiment). 
     N I  represents a count of intact fiber bundles on the sliding surface or contact surface of the soil column; N F  represents a total count of fiber bundles on the sliding surface or contact surface of the soil column. It can be found from the above equation that the size of N I  is affected by a contact force, and the size of N I  directly reflects the failure status of a foundation of the soil column. When N I =N F , it means that the fiber bundle on a failure surface of the soil column is not broken, that is, the soil column is intact without any damage. When N I &lt;N F , it means that the soil column fails along the failure surface. When a ratio of N I  to N F  reaches a failure threshold (which is set according to a specific slope), it can be considered that the soil column fails, that is, the source will initiate. By recording the value of N I , the failure level of the soil column can be predicted, that is, a potential failure-prone zone of the slope can be predicted. The purpose of disaster prevention and mitigation can be achieved by taking corresponding engineering measures before the slope fails. 
     The fiber bundle model is used to simulate the force mode of the six surrounding adjacent soil columns on the central soil column. According to the fiber bundle model, it is determined whether a ratio of the total count of fibers on the contact surface of adjacent soil columns to the count of intact fiber bundles is greater than the set threshold. If a determination result indicates no, it means that the connection bonds of a lateral tensile stress are not all broken. If the determination result indicates yes, it means that the connection bonds of the lateral tensile stress are all broken. The stress distribution of the soil column in a failure process is calculated according to a break status of the fiber bundle and the contact area between the soil columns. 
     When a soil column element of the slope is unstable along the most unfavorable sliding surface, under the action of root system, cementation, pore water, etc., the soil column distributes the compressive stress to the front soil column along the failure movement direction of the soil column through the connection bond, and distributes the tensile stress to the rear and side soil columns along the failure movement direction of the soil column. 
     As shown in a stress transfer process in  FIG. 4 , the calculation is as follows: 
     
       
         
           
             
               
                   
               
               ⁢ 
               
                 
                   A 
                   cn 
                 
                 = 
                 
                   [ 
                   
                     
                       min 
                       ⁡ 
                       
                         ( 
                         
                           
                             
                               H 
                               1 
                             
                             - 
                             
                               h 
                               rf 
                             
                           
                           , 
                           
                             
                               H 
                               
                                 1 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 n 
                               
                             
                             - 
                             
                               h 
                               rfn 
                             
                           
                         
                         ) 
                       
                     
                     - 
                     
                       ( 
                       
                         
                           H 
                           1 
                         
                         - 
                         
                           h 
                           rf 
                         
                         - 
                         
                           H 
                           sd 
                         
                       
                       ) 
                     
                   
                   ) 
                 
               
               ] 
             
             ⁢ 
             
               L 
               g 
             
           
         
       
       
         
           
             
               τ 
               
                 T 
                 ⁢ 
                 
                     
                 
                 ⁢ 
                 0 
               
             
             = 
             
               
                 
                   N 
                   0 
                 
                 
                   
                     ( 
                     
                       
                         N 
                         0 
                       
                       + 
                       
                         N 
                         1 
                       
                       + 
                       
                         N 
                         2 
                       
                       + 
                       
                         N 
                         3 
                       
                     
                     ) 
                   
                   + 
                   
                     
                       ( 
                       
                         
                           N 
                           4 
                         
                         + 
                         
                           N 
                           5 
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       ( 
                       
                         
                           K 
                           2 
                         
                         / 
                         
                           K 
                           1 
                         
                       
                       ) 
                     
                   
                 
               
               ⁢ 
               
                 
                   
                     A 
                     H 
                   
                   ⁡ 
                   
                     ( 
                     
                       
                         W 
                         d 
                       
                       - 
                       
                         F 
                         f 
                       
                       - 
                       
                         τ 
                         R 
                       
                     
                     ) 
                   
                 
                 
                   
                     cos 
                     ⁡ 
                     
                       ( 
                       
                         α 
                         + 
                         β 
                       
                       ) 
                     
                   
                   ⁢ 
                   
                     A 
                     cn 
                   
                 
               
             
           
         
       
       
         
           
             
               τ 
               
                 T 
                 ⁢ 
                 
                     
                 
                 ⁢ 
                 4 
               
             
             = 
             
               
                 
                   
                     N 
                     4 
                   
                   ⁡ 
                   
                     ( 
                     
                       
                         K 
                         2 
                       
                       / 
                       
                         K 
                         1 
                       
                     
                     ) 
                   
                 
                 
                   
                     ( 
                     
                       
                         N 
                         0 
                       
                       + 
                       
                         N 
                         1 
                       
                       + 
                       
                         N 
                         2 
                       
                       + 
                       
                         N 
                         3 
                       
                     
                     ) 
                   
                   + 
                   
                     
                       ( 
                       
                         
                           N 
                           4 
                         
                         + 
                         
                           N 
                           5 
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       ( 
                       
                         
                           K 
                           2 
                         
                         / 
                         
                           K 
                           1 
                         
                       
                       ) 
                     
                   
                 
               
               ⁢ 
               
                 
                   
                     A 
                     H 
                   
                   ⁡ 
                   
                     ( 
                     
                       
                         W 
                         d 
                       
                       - 
                       
                         F 
                         f 
                       
                       - 
                       
                         τ 
                         R 
                       
                     
                     ) 
                   
                 
                 
                   
                     cos 
                     ⁡ 
                     
                       ( 
                       
                         α 
                         + 
                         β 
                       
                       ) 
                     
                   
                   ⁢ 
                   
                     A 
                     cn 
                   
                 
               
             
           
         
       
     
     In the same way, calculation equations of τ T1 , τ T2 , τ T3  and τ T5  are derived, where n=0, 1, 2, 3, 4, 5. 
     A cn  represents a contact area between the central soil column and the six surrounding soil columns; τ T0  τ T1 , τ T2  and τ T3  respectively represent a tensile stress of soil column (i,j) to surrounding soil columns; τ T4  and τ T5  respectively represent a compressive stress of soil column (i,j) to surrounding soil columns; N n  represents a count of intact fibers between soil column (i,j) and a soil column numbered n; K 2 /K 1  represents a distribution ratio between the compressive stress and the tensile stress; h rf  represents an erosion depth of the central soil column (calculated based on an erosion model proposed by Iverson et al.); h rfn , represents an erosion depth of the surrounding soil column numbered n; H 1  represents a slope elevation of the central soil column; A H  represents a cross-sectional area of the soil column; L g  represents a spacing between adjacent soil columns; H 1n , represents a slope elevation of the surrounding soil column numbered n; τ R  represents a residual strength of the soil column after a failure along the most unfavorable sliding surface. 
     
       
         
           
             
               
                 τ 
                 c 
               
               &gt; 
               
                 
                   
                     2 
                     ⁢ 
                     c 
                     ⁢ 
                     
                       cos 
                       ⁡ 
                       
                         ( 
                         γ 
                         ) 
                       
                     
                   
                   
                     1 
                     - 
                     
                       sin 
                       ⁡ 
                       
                         ( 
                         γ 
                         ) 
                       
                     
                   
                 
                 + 
                 
                   
                     
                       2 
                       ⁢ 
                       
                         sin 
                         ⁡ 
                         
                           ( 
                           γ 
                           ) 
                         
                       
                     
                     
                       1 
                       - 
                       
                         sin 
                         ⁡ 
                         
                           ( 
                           γ 
                           ) 
                         
                       
                     
                   
                   ⁢ 
                   
                     
                       τ 
                       h 
                     
                     ⁡ 
                     
                       ( 
                       θ 
                       ) 
                     
                   
                 
               
             
             , 
           
         
       
     
     A fluidization criterion of the soil column is: where τ c  represents a compressive strength of the central soil column; τ h (θ) represents a shear strength of the central soil column related to the moisture content; γ represents an internal friction angle of the root-soil complex on the debris flow source to be predicted. 
       FIG. 5  shows a schematic diagram of the prediction system for an initiation volume of a debris flow slope source according to the present disclosure. The prediction system includes: a slope division module  201 , a sliding surface and unbalanced force calculation module  202 , an instability determination module  203 , a force mode determination module  204 , a connection bond break determination module  205 , a fluidization determination module  206 , a recording module  207  and a stopping module  208 . 
     The slope division module  201  is used for selecting a debris flow source slope to be predicted, and dividing the debris flow source slope to be predicted into regular hexagonal soil columns arranged in i rows and j columns. 
     The sliding surface and unbalanced force calculation module  202  is used for calculating a most unfavorable sliding surface of each soil column according to an upper bound theorem of a limit analysis, and calculating an unbalanced force on the most unfavorable sliding surface. 
     The instability determination module  203  is used for determining whether each soil column is unstable according to the unbalanced force on the most unfavorable sliding surface; stopping the prediction of the soil column and connecting the stopping module  208  if a determination result indicates no, and if the determination result indicates yes, connecting the force mode determination module  204 . 
     The force mode determination module  204  is used for acquiring centers of gravity and elevations of a central soil column and six adjacent soil columns around, determining an instability direction of the central soil column, and determining a force mode of the six surrounding adjacent soil columns on the central soil column. 
     The connection bond break determination module  205  is used for simulating the force mode of the six surrounding adjacent soil columns on the central soil column through a fiber bundle model, and determining a break status of connection bonds; stopping the prediction of the soil column and connecting the stopping module  208  if a determination result indicates that connection bonds of a lateral tensile stress are not all broken, and if the determination result indicates that the connection bonds of the lateral tensile stress are all broken, connecting the fluidization determination module  206 . 
     The fluidization determination module  206  is used for determining whether the soil column is fluidized according to a fluidization criterion; stopping the prediction of the soil column and connecting the stopping module  208  if a determination result indicates no, and if the determination result indicates yes, predicting that the soil column is about to initiate a debris flow, and predicting an initiation volume of the debris flow slope source according to a fluidization status. 
     The prediction system of the present disclosure further includes a recording module  207 , which is connected to the fluidization determination module  206 , and is used to record an initiation position and a corresponding initiation volume of the debris flow, so as to facilitate subsequent reference and analysis. The recording module  207  may further record that the debris flow initiation volume of a non-fluidized soil column is zero. The instability determination module  203  and the connection bond break determination module  205  may further be connected to the recording module  207  to record that the debris flow initiation volume of a soil column whose prediction is stopped is zero. 
     In a specific implementation process, the prediction method and system provided by the present disclosure were used to simulate and invert the possibility and failure level of a debris flow in a valley in southwest China under the action of rainfall. In this example, the terrain parameters, hydrological parameters and soil parameters in the valley required for modelling were acquired through ArcGIS. The method and system were programmed by C′ language to calculate the time-dependent initiation volume and cumulative total initiation volume of the debris flow in the valley, and the calculation results were imported into ArcGIS to simulate the gradual initiation process of the debris flow in the valley. 
     In the present disclosure, the most unfavorable sliding surface of the soil column is calculated by an upper bound theorem of a limit analysis. Compared with the landslide Hydro-mechanical Triggering (LHT) model assuming that the failure surface of the soil is on the bedrock surface, the calculation accuracy of the present disclosure is higher, and it is more in line with the actual initiation process of the debris flow on the slope.  FIG. 6  shows cumulative time-dependent curves of the measured and inverted initiation volumes of a debris flow at a mouth of the valley, and  FIG. 7  shows time-dependent curves of the measured and inverted initiation volumes of the debris flow at the mouth of the valley. By comparing the initiation volume and the cumulative initiation volume of the debris flow, it is found that a greater error is led in the prediction result of the LHT model. The inversion result of the present disclosure and the actual measured initiation volume of the debris flow have strong volatility, and have basically consistent peak values. The result of the total initiation volume of the debris flow derived by the simulation and inversion of the present disclosure is basically consistent with the actual observation value, and the error is small. This indicates that the prediction accuracy and precision of the present disclosure are higher. 
     Application effect 1: A valley where debris flows are not frequent but a large debris flow has occurred in history was selected, and the topographical parameters, soil parameters and historical hydrological parameters of the valley were acquired through on-site investigation, literature analysis and real-time monitoring. According to historical rainfall and hydrological data, the present disclosure predicted the initiation volume and source position of a debris flow that may be initiated in the valley under historical rainfall conditions. The peak initiation volume and total initiation volume of the debris flow can be used to determine the engineering prevention level, rationally design the size and quantity of the blocking dam, drainage channel and other projects, so as to reduce the waste of resources caused by unreasonable engineering design. The present disclosure can achieve the effect of preventing and controlling debris flows by minimizing the use of resources. 
     Application effect 2: According to the present disclosure, real-time monitoring is required for a valley with frequent debris flows. A debris flow early warning and forecasting platform can be established, which integrates functions such as real-time rainfall acquisition, initiation volume calculation and debris flow impact range division. According to the real-time acquisition of rainfall data, the present disclosure can be used to predict the initiation volume of the debris flow in real time, set a debris flow initiation volume-induced disaster threshold or level, realize the early warning and forecast of the debris flow before the initiation reaches the disaster threshold or level, and evacuate the affected people in time. The present disclosure can also quickly determine the location and scale of the debris flow after the debris flow occurs, which is helpful for planning rescue implementation plans and carrying out rescue and dredging operations in the disaster-affected area. 
     The present disclosure can meet the needs of various debris flow prevention and mitigation requirements, and is applied to different disaster prevention and mitigation scenarios. 
     For a system disclosed in the embodiments, since the system corresponds to the method disclosed in the embodiments, the description is simple, and reference may be made to the method description. 
     Several embodiments are used for illustration of the principles and implementations of the present disclosure. The description of the embodiments is used to help illustrate the method and its core principles of the present disclosure. In addition, those skilled in the art can make various modifications in terms of specific embodiments and scope of application in accordance with the teachings of the present disclosure. In conclusion, the content of the present specification should not be construed as a limitation to the present disclosure.