Geological disaster risk assessment method and device and storage medium

Abstract

The invention provides a geological disaster risk assessment method and device and a storage medium, and belongs to the technical field of geological assessment, the method comprises the following steps: importing geological data at the current moment, and carrying out infiltration depth analysis on the geological data to obtain a target infiltration depth; performing stability analysis on the geological data and the target infiltration depth to obtain a stability coefficient; performing kinetic analysis on the geological data according to the stability coefficient to obtain a target kinetic parameter set; and evaluating the stability coefficient and the target kinetic parameter set to obtain a geological disaster risk evaluation result. According to the method, the risk level of the geological disaster can be evaluated more accurately, the risk evaluation precision is improved, the possible influence range and degree of the disaster can be predicted more accurately, and a scientific basis is provided for disaster prevention and reduction work.

Description

Technical Field

[0001] This invention relates to the field of geological assessment technology, specifically to a geological hazard risk assessment method, device, and storage medium. Background Technology

[0002] Geological disasters refer to hazardous phenomena caused by changes in the geological environment due to natural factors or human activities, such as landslides, debris flows, and ground subsidence. These disasters are characterized by their suddenness, destructive power, and wide-ranging impact, posing a significant threat and challenge to human society.

[0003] To effectively address geological disasters, accurate risk assessments are necessary. However, traditional geological disaster risk assessment methods mainly include historical disaster case analysis, expert scoring, and fuzzy comprehensive evaluation. While these methods are simple, easy to implement, and highly operable, they suffer from drawbacks such as high subjectivity and low accuracy, making them insufficient for addressing geological disasters. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a geological disaster risk assessment method, device and storage medium to address the shortcomings of the prior art.

[0005] The technical solution of this invention to solve the above-mentioned technical problems is as follows: A geological disaster risk assessment method, comprising the following steps: Import the geological data at the current moment, perform infiltration depth analysis on the geological data, and obtain the target infiltration depth; Stability analysis was performed on the geological data and the target infiltration depth to obtain the stability coefficient; Based on the stability coefficient, a dynamic analysis is performed on the geological data to obtain a target dynamic parameter set; The stability coefficient and the target dynamic parameter set are evaluated to obtain the geological hazard risk assessment results.

[0006] Another technical solution of the present invention to solve the above-mentioned technical problems is as follows: A geological disaster risk assessment device, comprising: The import module is used to import geological data at the current time. The infiltration depth analysis module is used to perform infiltration depth analysis on the geological data to obtain the target infiltration depth; The stability analysis module is used to perform stability analysis on the geological data and the target infiltration depth to obtain the stability coefficient; The dynamic analysis module is used to perform dynamic analysis on the geological data based on the stability coefficient to obtain a target dynamic parameter set; The evaluation result acquisition module is used to evaluate the stability coefficient and the target dynamic parameter set to obtain the geological hazard risk assessment result.

[0007] Based on the above-mentioned geological disaster risk assessment method, the present invention also provides a geological disaster risk assessment system.

[0008] Another technical solution of the present invention to solve the above-mentioned technical problems is as follows: a geological disaster risk assessment system, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the geological disaster risk assessment method as described above.

[0009] Based on the above-mentioned geological disaster risk assessment method, the present invention also provides a computer-readable storage medium.

[0010] Another technical solution of the present invention to solve the above-mentioned technical problems is as follows: a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the geological disaster risk assessment method as described above.

[0011] The beneficial effects of this invention are as follows: by analyzing the infiltration depth of geological data to obtain the target infiltration depth, by analyzing the stability of geological data and the target infiltration depth to obtain the stability coefficient, by analyzing the dynamics of geological data based on the stability coefficient to obtain the target dynamic parameter set, and by evaluating the stability coefficient and the target dynamic parameter set to obtain the geological hazard risk assessment result, the risk level of geological hazards can be assessed more accurately, the accuracy of risk assessment is improved, and it helps to more accurately predict the possible scope and degree of impact of disasters, providing a scientific basis for disaster prevention and mitigation work. Attached Figure Description

[0012] Figure 1 A schematic flowchart of the geological hazard risk assessment method provided in an embodiment of the present invention; Figure 2 A block diagram of a geological disaster risk assessment device provided in an embodiment of the present invention. Detailed Implementation

[0013] The principles and features of the present invention are described below with reference to the accompanying drawings. The examples given are only for explaining the present invention and are not intended to limit the scope of the present invention.

[0014] Figure 1 This is a flowchart illustrating a geological hazard risk assessment method provided in an embodiment of the present invention.

[0015] like Figure 1 As shown, a geological hazard risk assessment method includes the following steps: S1: Import the geological data at the current moment, perform infiltration depth analysis on the geological data, and obtain the target infiltration depth; S2: Perform stability analysis on the geological data and the target infiltration depth to obtain the stability coefficient; S3: Perform dynamic analysis on the geological data based on the stability coefficient to obtain the target dynamic parameter set; S4: Evaluate the stability coefficient and the target dynamic parameter set to obtain the geological hazard risk assessment results.

[0016] It should be understood that the geological data may be geological data of a certain area or a designated area.

[0017] In the above embodiments, the target infiltration depth is obtained by analyzing the infiltration depth of geological data, the stability coefficient is obtained by analyzing the stability of geological data and the target infiltration depth, the target dynamic parameter set is obtained by dynamic analysis of geological data based on the stability coefficient, and the geological hazard risk assessment result is obtained by evaluating the stability coefficient and the target dynamic parameter set. This can more accurately assess the risk level of geological hazards, improve the accuracy of risk assessment, help to more accurately predict the possible impact range and degree of disasters, and provide a scientific basis for disaster prevention and mitigation work.

[0018] Optionally, as an embodiment of the present invention, the geological data includes saturated water content of the wet zone, initial soil water content, slope gradient, saturated permeability coefficient, depth of the wetting front, average matrix suction head of the wetting front, and rainfall intensity. The process of performing infiltration depth analysis on the geological data to obtain the target infiltration depth includes: The target infiltration depth is obtained by calculating the saturated water content of the humid area, the initial water content of the soil, the slope gradient, the saturated permeability coefficient, the depth of the wetting front, the average matrix suction head of the wetting front, and the rainfall intensity using the first formula: , in, , in, For the target infiltration depth, and All are coefficients. , , This refers to the slope infiltration depth. For the current moment, Rainfall intensity, For slope gradient, This represents the saturated water content of the humid region. This represents the initial moisture content of the soil. The time when water accumulation begins. The saturated permeability coefficient, This refers to the depth at the moistened front. The average substrate suction head is the wetting front.

[0019] It should be understood that parameters such as the infiltration coefficient (i.e., geological data) are input into the GA model to obtain the saturated and unsaturated infiltration results (i.e., the target infiltration depth) at each time step.

[0020] It should be understood that the depth of the wet front and the rainfall intensity are both obtained based on the rainfall process.

[0021] Specifically, the model uses an improved GA method to simulate the spatiotemporal dynamics of surface runoff and subsurface infiltration under rainfall conditions.

[0022] The GA slope infiltration model considering the effect of slope can be expressed as: , z w The depth at the moistening front; α is the slope gradient; h The water head at the surface; k s The saturated permeability coefficient; φ w The average substrate suction head is the wetting front.

[0023] The formula for calculating the relationship between slope infiltration depth and time is: , t For a certain moment; θ s This represents the saturated water content of the humid region. θ i This refers to the initial moisture content of the soil. t p This is the time when water began to accumulate.

[0024] Considering that in reality, the moisture content curve after water accumulation is not a single regular rectangle, but rather a segment of length... z s The rectangular saturation region and a length of z us For the elliptical unsaturated region, the GA model is modified to obtain the modified infiltration depth, as shown in the following formula: , a , b Let be the coefficient, and a <0, 0< b<1 can be obtained by fitting experimental data.

[0025] In the above embodiments, infiltration depth analysis of geological data is performed to obtain the target infiltration depth, and the spatiotemporal dynamic changes of surface runoff and underground seepage under rainfall conditions are simulated. This can more accurately assess the risk level of geological disasters, improve the accuracy of risk assessment, and help to more accurately predict the possible scope and extent of disaster impact.

[0026] Optionally, as an embodiment of the present invention, the geological data further includes effective internal friction angle, effective cohesion, saturated soil unit weight, and water unit weight; The process of performing stability analysis on the geological data and the target infiltration depth to obtain the stability coefficient includes: The stability coefficient is obtained by calculating the target infiltration depth, the effective internal friction angle, the slope gradient, the effective cohesion, the saturated soil unit weight, and the water unit weight using the second formula. The second formula is: , in, The stability coefficient, For the target infiltration depth, For the current moment, For the effective internal friction angle, For slope gradient, For effective cohesion, For water density, This is the unit weight of saturated soil.

[0027] It should be understood that Stable adjusts the choice of stability calculation method, controls parameters such as cohesion c, internal friction angle φ and unit weight, and adopts the stability calculation method of slope element to obtain the stability results (i.e. stability coefficients) of grids and slope elements at different time steps.

[0028] Specifically, the unstable elements are determined using infinite slope simulation or three-dimensional stability calculation methods.

[0029] The stability calculation uses an infinite slope model, combined with traditional slope stability calculation formulas, to obtain stability coefficients at different depths, as shown in the following formula: , γ d It is the bulk density of dry soil. γ s It is the unit weight of saturated soil. γ w α is the specific weight of water, and α is the slope gradient.

[0030] In the above embodiments, stability coefficients are obtained by performing stability analysis on geological data and target infiltration depth, which improves the accuracy of risk assessment and helps to more accurately predict the possible scope and extent of disaster impact.

[0031] Optionally, as an embodiment of the present invention, the process of performing dynamic analysis on the geological data based on the stability coefficient to obtain the target dynamic parameter set includes: Import the historical dynamic parameter set, determine whether the stability coefficient is less than a first preset stability threshold. If so, perform dynamic calculations on the geological data to obtain the original dynamic parameter set, and combine the original dynamic parameter set and the historical dynamic parameter set to obtain the target dynamic parameter set. If not, use the historical dynamic parameter set as the target dynamic parameter set.

[0032] Preferably, the first preset stability threshold can be 1.

[0033] It should be understood that dynamic calculations are performed on slope elements with stability Fs (i.e., stability coefficient) ≤ 1 to obtain the complete dynamic motion process of each slope element (i.e., the original set of dynamic parameters).

[0034] In the above embodiments, the target dynamic parameter set is obtained by performing dynamic analysis on geological data based on the stability coefficient. This can more accurately assess the risk level of geological disasters, improve the accuracy of risk assessment, and help to more accurately predict the possible scope and extent of disaster impact, thus providing a scientific basis for disaster prevention and mitigation work.

[0035] Optionally, as an embodiment of the present invention, the geological data further includes water density and block parameters of multiple blocks, the block parameters including block side length, block thickness, location coordinates and block velocity; The process of performing dynamic calculations on the geological data to obtain the original set of dynamic parameters includes: The third formula is used to calculate the side length of each block, the thickness of each block, and the position coordinates of each block, respectively, to obtain multiple target volumes corresponding to each block. , in, For the first The block moves towards the first The target volume of each block transfer. For the first The thickness of each block. For the first The side length of each block. For the first The position coordinates of each block For the first The position coordinates of each block; The fourth equation is used to calculate the water density, the side length of each block, the thickness of each block, the velocity of each block, and the position coordinates of each block, respectively, to obtain multiple target momentum corresponding to each block. The original dynamic parameter set is then obtained by combining all target volumes and all target momentum. The fourth equation is: , in, For the first The block moves towards the first The target momentum of the block transfer For the density of water, For the first The thickness of each block. For the first The side length of each block. For the first The position coordinates of each block For the first The position coordinates corresponding to each block.

[0036] It should be understood that the dynamic process employs the Tsunami Squares (TS) numerical simulation method, applicable to various fluids and fluid-like motions. Breaking away from the traditional model of establishing and solving numerical equations, it adopts a method of moving and dividing blocks, with each time step following volume and momentum conservation for iterative iteration. Therefore, it eliminates the need to solve equations and handle complex boundary conditions and mesh breakage. Utilizing the geometric properties of squares, the conventional calculation of traversing N squared blocks is reduced to 4N blocks, decreasing the computational load by N orders of magnitude and significantly improving computational efficiency.

[0037] Specifically, it incorporates multiple forms of motion mechanisms, as shown in the following formula: , a It is acceleration.

[0038] The entire process satisfies the laws of conservation of mass and momentum, as shown in the following formulas: , H The thickness of the block; M Momentum; D The length of the block's side; v For speed; V For volume; ρ w This is the density of water.

[0039] In conventional calculations, N 2 Adding 1 square becomes adding 4N, greatly reducing the amount of calculation. The formula is as follows: , .

[0040] In the above embodiments, dynamic calculations are performed on geological data to obtain the original set of dynamic parameters, which improves calculation efficiency, greatly reduces the amount of calculation, improves the accuracy of risk assessment, helps to more accurately predict the possible scope and degree of disaster impact, and provides a scientific basis for disaster prevention and mitigation work.

[0041] Optionally, as an embodiment of the present invention, the process of evaluating the stability coefficient and the target dynamic parameter set to obtain the geological hazard risk assessment result includes: The minimum value of the target dynamic parameter set is selected to obtain the minimum dynamic parameter; The maximum value of the target dynamic parameter set is selected to obtain the maximum dynamic parameter; The range of dynamic parameters is obtained by taking the minimum dynamic parameter as the left boundary and the maximum dynamic parameter as the right boundary; When the stability coefficient is less than or equal to the first preset stability threshold, or when the range of the dynamic parameters meets the first preset dynamic conditions, the preset first assessment result is taken as the geological disaster risk assessment result. When the stability coefficient is greater than the first preset stability threshold and less than or equal to the second preset stability threshold, or when the range of the dynamic parameters meets the second preset dynamic conditions, the preset second assessment result is taken as the geological disaster risk assessment result. When the stability coefficient is greater than the second preset stability threshold and less than or equal to the third preset stability threshold, or when the range of the dynamic parameters meets the third preset dynamic conditions, the preset third assessment result is taken as the geological hazard risk assessment result. When the stability coefficient is greater than the third preset stability threshold, or when the range of the dynamic parameters meets the fourth preset dynamic condition, the preset fourth assessment result is taken as the geological hazard risk assessment result.

[0042] Preferably, the second preset stability threshold can be 1.05, and the third preset stability threshold can be 1.25.

[0043] Specifically, according to the stability of the slope unit (i.e., the stability coefficient) and the dynamic calculation results (i.e., the range of dynamic parameters), a risk grading standard is established. The risk levels of geological disaster-affected bodies are divided into four levels, namely: extremely high risk, high risk, medium risk, and low risk. Specifically: ① For houses within slope units where Fs (i.e., the stability coefficient) ≤ 1.0 (i.e., the first preset stability threshold), or houses within the range of dynamic calculation results (i.e., the first preset dynamic condition), the risk level is extremely high (i.e., the preset first evaluation result); ② For houses within slope units where 1.0 (i.e., the first preset stability threshold) < Fs (i.e., the stability coefficient) ≤ 1.05 (i.e., the second preset stability threshold), or houses within 5 meters outside the range of dynamic calculation results (i.e., the second preset dynamic condition), the risk level is high (i.e., the preset second evaluation result); ③ For houses within slope units where 1.05 (i.e., the second preset stability threshold) < Fs (i.e., the stability coefficient) ≤ 1.25 (i.e., the third preset stability threshold), or houses within 5 - 10 meters outside the range of dynamic calculation results (i.e., the third preset dynamic condition), the risk level is medium (i.e., the preset third evaluation result); ④ For houses within slope units where Fs (i.e., the stability coefficient) ≥ 1.25 (i.e., the third preset stability threshold), or houses within 10 - 20 meters outside the range of dynamic calculation results (i.e., the fourth preset dynamic condition), the risk level is low (i.e., the preset fourth evaluation result).

[0044] Substitute the calculated landslide stability (i.e., the stability coefficient) and dynamic results (i.e., the range of dynamic parameters) into the above judgment criteria to determine the risk level of regional landslide disasters.

[0045] In the above embodiments, by evaluating the stability coefficient and the target dynamic parameter set, the geological disaster risk assessment result is obtained, which improves the accuracy of risk assessment, helps to more accurately predict the possible influence range and degree of disasters, and provides a scientific basis for disaster prevention and mitigation work.

[0046] Optionally, as another embodiment of the present invention, with the progress of science and technology, the regional geological disaster risk assessment technical method considering the dynamic process can more accurately evaluate the risk level of geological disasters by simulating the occurrence, development, and evolution process of geological disasters.

[0047] Optionally, as another embodiment of the present invention, the present invention mainly includes: 1. Currently, there is little hourly geological disaster risk prediction for real-time rainfall. To address this problem, the present invention can simulate the occurrence location and influence range according to the predicted short-term rainfall in the future at a custom time interval; 2. Traditional regional stability calculation methods that consider rainfall characteristics only calculate the stability of each calculation unit, without simultaneously performing subsequent dynamic calculations of unstable regions. This invention can perform dynamic simulations of unstable units after stability calculations and provide calculation results that simultaneously include both stability and dynamic motion paths.

[0048] 3. In performing dynamic calculations, conventional methods mostly perform dynamic simulations on unstable elements within a fixed time interval after the calculation is completed. This invention can perform dynamic simulation calculations on unstable elements within a defined time step throughout the entire calculation process.

[0049] 4. By introducing numerical simulation of dynamic processes, the occurrence and development of geological disasters can be simulated more accurately, thereby improving the accuracy of risk assessment. This helps to more accurately predict the potential scope and extent of disaster impact, providing a scientific basis for disaster prevention and mitigation efforts.

[0050] Optionally, as another embodiment of the present invention, the present invention specifically includes: (1) Layer input; (2) Parameter input; (3) Calculation of grid cells (4) Calculation method of the Main control module, terrain file, rainfall input and total calculation time; (5) Adjust the landslide modeling method in the LS file, and use it as the main control file for stability and dynamics calculations; (6) The Outline and Thickness files correspond to the slope element outline and thickness, respectively; The above documents together complete the landslide modeling and begin the next step of calculation; (7) Input parameters such as the infiltration coefficient into the GA model to obtain the saturated and unsaturated infiltration results at each time step; (8) Stable adjusts the stability calculation method selection, controls parameters such as cohesion c, internal friction angle φ and unit weight, and adopts the stability calculation method of slope element to obtain the stability results of grids and slope elements at different time steps; (9) Dynamic reads in the friction and acceleration parameters, performs dynamic calculations for the slope element with stability Fs≤1, and obtains the complete dynamic motion process of each slope element.

[0051] (10) Based on the stability and dynamic calculation results of the slope unit, a risk classification standard is established. The risk level of the geological hazard-bearing body is divided into four levels: extremely high risk, high risk, medium risk, and low risk, specifically: ① Houses within the slope unit where Fs ≤ 1.0 or houses within the range of the dynamic calculation results have an extremely high risk level; ② Houses within the slope unit where 1.0 < H ≤ 1.05 or houses within 5 meters outside the range of the dynamic calculation results have a high risk level; ③ Houses within the slope unit where 1.05 < H ≤ 1.25 or houses within 5 - 10 meters outside the range of the dynamic calculation results have a medium risk level; ④ Houses within the slope unit where H ≥ 1.25 or houses within 10 - 20 meters outside the range of the dynamic calculation results have a low risk level.

[0052] Substitute the landslide stability and dynamic results obtained from steps 8 - 9 into the above judgment criteria to determine the risk level of regional landslide disasters.

[0053] (11)Repeat steps 7 - 9 until the entire rainfall duration is completed, and save the calculation results for each time step.

[0054] Figure 2 This is a module block diagram of a geological disaster risk assessment device provided by an embodiment of the present invention.

[0055] Optionally, as another embodiment of the present invention, as Figure 2 shown, a geological disaster risk assessment device includes: An import module for importing geological data at the current moment; An infiltration depth analysis module for performing infiltration depth analysis on the geological data to obtain the target infiltration depth; A stability analysis module for performing stability analysis on the geological data and the target infiltration depth to obtain the stability coefficient; A dynamic analysis module for performing dynamic analysis on the geological data based on the stability coefficient to obtain the target dynamic parameter set; An evaluation result acquisition module for evaluating the stability coefficient and the target dynamic parameter set to obtain the geological disaster risk assessment result.

[0056] Optionally, another embodiment of the present invention provides a geological disaster risk assessment system, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the above - mentioned geological disaster risk assessment method is implemented. This system can be a computer and other systems.

[0057] Optionally, another embodiment of the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the geological hazard risk assessment method as described above.

[0058] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0059] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of the above-described apparatus and unit can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0060] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed.

[0061] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of the embodiments of the present invention, depending on actual needs.

[0062] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0063] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0064] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for assessing geological hazard risk, characterized in that, Includes the following steps: Import the geological data at the current moment, perform infiltration depth analysis on the geological data, and obtain the target infiltration depth; Stability analysis was performed on the geological data and the target infiltration depth to obtain the stability coefficient; Based on the stability coefficient, a dynamic analysis is performed on the geological data to obtain a target dynamic parameter set; The stability coefficient and the target dynamic parameter set are evaluated to obtain the geological hazard risk assessment results.

2. The geological hazard risk assessment method according to claim 1, characterized in that, The geological data include saturated water content of the humid area, initial soil water content, slope gradient, saturated permeability coefficient, depth of the wetting front, average matrix suction head of the wetting front, and rainfall intensity. The process of performing infiltration depth analysis on the geological data to obtain the target infiltration depth includes: The target infiltration depth is obtained by calculating the saturated water content of the humid area, the initial water content of the soil, the slope gradient, the saturated permeability coefficient, the depth of the wetting front, the average matrix suction head of the wetting front, and the rainfall intensity using the first formula: ， in, , in, For the target infiltration depth, and All are coefficients. , , This refers to the slope infiltration depth. For the current moment, Rainfall intensity, For slope gradient, This represents the saturated water content of the humid region. This represents the initial moisture content of the soil. The time when water accumulation begins. The saturated permeability coefficient, This refers to the depth at the moistened front. The average substrate suction head is the wetting front.

3. The geological hazard risk assessment method according to claim 2, characterized in that, The geological data also includes effective internal friction angle, effective cohesion, saturated soil unit weight, and water unit weight; The process of performing stability analysis on the geological data and the target infiltration depth to obtain the stability coefficient includes: The stability coefficient is obtained by calculating the target infiltration depth, the effective internal friction angle, the slope gradient, the effective cohesion, the saturated soil unit weight, and the water unit weight using the second formula. The second formula is: ， in, The stability coefficient, For the target infiltration depth, For the current moment, For the effective internal friction angle, For slope gradient, For effective cohesion, For water density, This is the unit weight of saturated soil.

4. The geological hazard risk assessment method according to claim 2, characterized in that, The process of performing dynamic analysis on the geological data based on the stability coefficient to obtain the target dynamic parameter set includes: Import the historical dynamic parameter set, determine whether the stability coefficient is less than a first preset stability threshold. If so, perform dynamic calculations on the geological data to obtain the original dynamic parameter set, and combine the original dynamic parameter set and the historical dynamic parameter set to obtain the target dynamic parameter set. If not, use the historical dynamic parameter set as the target dynamic parameter set.

5. The geological hazard risk assessment method according to claim 4, characterized in that, The geological data also includes water density and block parameters of multiple blocks, including block side length, block thickness, location coordinates, and block velocity; The process of performing dynamic calculations on the geological data to obtain the original set of dynamic parameters includes: The third formula is used to calculate the side length of each block, the thickness of each block, and the position coordinates of each block, respectively, to obtain multiple target volumes corresponding to each block. ， in, For the first The block moves towards the first The target volume of each block transfer. For the first The thickness of each block. For the first The side length of each block. For the first The position coordinates of each block For the first The position coordinates of each block; The fourth equation is used to calculate the water density, the side length of each block, the thickness of each block, the velocity of each block, and the position coordinates of each block, respectively, to obtain multiple target momentum corresponding to each block. The original dynamic parameter set is then obtained by combining all target volumes and all target momentum. The fourth equation is: ， in, For the first The block moves towards the first The target momentum of the block transfer For the density of water, For the first The thickness of each block. For the first The side length of each block. For the first The position coordinates of each block For the first The position coordinates corresponding to each block.

6. The geological hazard risk assessment method according to claim 1, characterized in that, The process of evaluating the stability coefficient and the target dynamic parameter set to obtain the geological hazard risk assessment result includes: The minimum value of the target dynamic parameter set is selected to obtain the minimum dynamic parameter; The maximum value of the target dynamic parameter set is selected to obtain the maximum dynamic parameter; The range of dynamic parameters is obtained by taking the minimum dynamic parameter as the left boundary and the maximum dynamic parameter as the right boundary; When the stability coefficient is less than or equal to the first preset stability threshold, or when the range of the dynamic parameters meets the first preset dynamic conditions, the preset first assessment result is taken as the geological disaster risk assessment result. When the stability coefficient is greater than the first preset stability threshold and less than or equal to the second preset stability threshold, or when the range of the dynamic parameters meets the second preset dynamic conditions, the preset second assessment result is taken as the geological disaster risk assessment result. When the stability coefficient is greater than the second preset stability threshold and less than or equal to the third preset stability threshold, or when the range of the dynamic parameters meets the third preset dynamic conditions, the preset third assessment result is taken as the geological hazard risk assessment result. When the stability coefficient is greater than the third preset stability threshold, or when the range of the dynamic parameters meets the fourth preset dynamic condition, the preset fourth assessment result is taken as the geological hazard risk assessment result.

7. A geological hazard risk assessment device, characterized in that, include: The import module is used to import geological data at the current time. The infiltration depth analysis module is used to perform infiltration depth analysis on the geological data to obtain the target infiltration depth; The stability analysis module is used to perform stability analysis on the geological data and the target infiltration depth to obtain the stability coefficient; The dynamic analysis module is used to perform dynamic analysis on the geological data based on the stability coefficient to obtain a target dynamic parameter set; The evaluation result acquisition module is used to evaluate the stability coefficient and the target dynamic parameter set to obtain the geological hazard risk assessment result.

8. The geological hazard risk assessment device according to claim 7, characterized in that, The geological data include saturated water content of the humid area, initial soil water content, slope gradient, saturated permeability coefficient, depth of the wetting front, average matrix suction head of the wetting front, and rainfall intensity. The infiltration depth analysis module is specifically used for: The target infiltration depth is obtained by calculating the saturated water content of the humid area, the initial water content of the soil, the slope gradient, the saturated permeability coefficient, the depth of the wetting front, the average matrix suction head of the wetting front, and the rainfall intensity using the first formula: ， in, , in, For the target infiltration depth, and All are coefficients. , , This refers to the slope infiltration depth. For the current moment, Rainfall intensity, For slope gradient, This represents the saturated water content of the humid region. This represents the initial moisture content of the soil. The time when water accumulation begins. The saturated permeability coefficient, This refers to the depth at the moistened front. The average substrate suction head is the wetting front.

9. A geological hazard risk assessment device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the geological hazard risk assessment method as described in any one of claims 1 to 6.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the geological hazard risk assessment method as described in any one of claims 1 to 6.

Images