Landslide simulation method, system and equipment based on kinematic wave equation and storage medium

Abstract

The invention discloses a landslide simulation method, system and device based on a kinematic wave equation and a storage medium, and belongs to the technical field of landslide simulation. Parameters required by simulation are obtained, wherein the parameters comprise the terrain gradient; calculating a landslide movement speed by using a friction model; and the landslide movement speed is specially processed through the terrain gradient. The landslide movement is simulated by adopting the kinematic wave equation based on the terrain gradient data obtained by remote sensing, the parameter-terrain gradient required by simulation is obtained, the landslide movement speed is calculated by using the friction model, and the speed is specially processed through the terrain gradient, so that the simulation stability can be improved; the landslide accumulation depth and the flowing state at different moments are simulated; according to the method, the motion wave theory and the friction model are combined, the motion characteristics of the landslide under different conditions can be effectively simulated through accurate calculation, and a scientific basis is provided for disaster prevention and reduction of the landslide.

Description

Technical Field

[0001] The present invention belongs to the technical field of landslide simulation, and particularly relates to a landslide simulation method, system, device and storage medium based on the kinematic wave equation. Background Art

[0002] The simulation and prediction of landslide disasters are important research directions in the field of geological disaster prevention and control. Traditional landslide simulation methods mainly include numerical models based on physical processes and simplified models based on statistical experience. However, numerical models based on physical processes often have problems such as low computational efficiency and slow simulation speed in practical applications, making it difficult to meet the needs of rapid risk assessment and emergency response. Simplified models based on statistical experience usually have a faster calculation speed and are suitable for large-scale landslide risk assessment. However, such methods rely on historical data, and the prediction accuracy decreases when there is a lack of sufficient training samples or significant changes in geological conditions. In addition, statistical models usually lack a clear physical mechanism and are difficult to effectively simulate the dynamic evolution process of landslides, limiting their application in refined disaster simulation. Summary of the Invention

[0003] The purpose of the present invention is to provide a landslide simulation method, system, device and storage medium based on the kinematic wave equation to solve the problem that the existing technology is difficult to effectively simulate the dynamic evolution process of landslides.

[0004] The embodiments of the present application are implemented as follows. The landslide simulation method based on the kinematic wave equation includes the following steps: Simulate the landslide movement using the kinematic wave equation; Obtain the parameters required for the simulation, and the parameters include the terrain slope; Calculate the landslide movement speed using the friction model; Perform special processing on the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0005] Optionally, in some embodiments of the present application, the continuity equation of the kinematic wave equation for simulating the landslide movement is: ; In the formula, is the landslide depth, unit: ; is the landslide movement speed along the direction, unit: ; is the landslide movement speed along the direction, unit: ; is the time of the landslide movement, unit: .

[0006] Optionally, in some embodiments of the present application, the terrain slope includes the terrain slope in the x direction and the terrain slope in the y direction .

[0007] Optionally, in some embodiments of the present application, the friction model is used to calculate the landslide movement speed: ; ; wherein, is the landslide depth, unit: ; is the depth of landslide accumulation under constant conditions, unit: ; is the movement speed of the landslide along the direction, unit: ; is the movement speed of the landslide along the direction, unit: ; is the terrain slope in the x direction, unit: dimensionless; is the terrain slope in the y direction, unit: dimensionless; is the flow efficiency, unit: .

[0008] Optionally, in some embodiments of the present application, the calculation formula of is: , where g is the acceleration due to gravity, unit:

[0009] Optionally, in some embodiments of the present application, the depth of landslide accumulation under constant conditions is: ; wherein, is the characteristic length, unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless.

[0010] Optionally, in some embodiments of the present application, the terrain slope has an impact on the landslide movement speed and is specially processed as follows: ; In the formula, is the landslide movement speed, , unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless; is a constant, whose value is 0.0001, unit: dimensionless.

[0011] Correspondingly, the embodiments of the present application also provide a landslide simulation system based on the kinematic wave equation, including: A module for simulating landslide movement, which simulates landslide movement using the kinematic wave equation; A parameter acquisition module for obtaining the parameters required for simulation, and the parameters include the terrain slope; A landslide movement speed calculation model for calculating the landslide movement speed using a friction model; A special processing model for specially processing the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0012] Correspondingly, the embodiments of the present application also provide a computer device, including a storage and a processor. When the computer program stored in the storage is executed by the processor, the processor executes the steps of the above method.

[0013] Correspondingly, the embodiments of the present application also provide a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the processor executes the steps of the above method.

[0014] In summary, due to the adoption of the above technical solutions, the beneficial effects of the present invention are: Based on the terrain slope data obtained by remote sensing, this application uses the kinematic wave equation to simulate landslide movement, obtains the parameters required for simulation - terrain slope, calculates the landslide movement speed using a friction model, and performs special processing on the speed through the terrain slope to improve the stability of the simulation, so as to simulate the landslide accumulation depth and flow state at different times; this application combines the kinematic wave theory and the friction model, and through precise calculation, can effectively simulate the movement characteristics of landslides under different conditions, providing a scientific basis for landslide disaster prevention and mitigation. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 is a flowchart of the landslide simulation method based on the kinematic wave equation of the present invention; Figure 2 is the accumulation situation diagram of the landslide provided by the embodiment of the present invention at different times (unit: ) DETAILED DESCRIPTION OF THE EMBODIMENTS

[0016] In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0017] The technical solution of this application is as follows: Please refer to Figure 1 , on the first aspect, the embodiment of this application provides a landslide simulation method based on the kinematic wave equation, including the following steps: S01. Simulate the landslide movement using the kinematic wave equation; S02. Obtain the parameters required for simulation, and the parameters include terrain slope; S03. Calculate the landslide movement speed using a friction model; S04. Perform special processing on the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0018] Based on the terrain slope data obtained by remote sensing, this application uses the kinematic wave equation to simulate landslide movement, obtains the parameters required for simulation - terrain slope, calculates the landslide movement speed using a friction model, and performs special processing on the speed through the terrain slope to improve the stability of the simulation, so as to simulate the landslide accumulation depth and flow state at different times; this application combines the kinematic wave theory and the friction model, and through precise calculation, can effectively simulate the movement characteristics of landslides under different conditions, providing a scientific basis for landslide disaster prevention and mitigation.

[0019] According to the landslide simulation method proposed in this application, key parameters such as the movement speed and accumulation depth of the landslide can be obtained, and the corresponding departments can take disaster response measures in advance according to the simulation results; the landslide simulation method of this application fully considers the depth, speed, and terrain slope of the landslide, which helps to reveal the spatio-temporal distribution law of landslide activities and provides technical support for regional landslide risk management and early warning system construction; the landslide simulation method of this application simplifies the landslide dynamic model based on physical processes, and greatly accelerates the calculation efficiency while maintaining high accuracy, providing technical support for regional landslide risk management and early warning system construction.

[0020] In the real world, the influence of terrain slope on speed is usually continuous and smooth. By specially processing the landslide movement speed through the terrain slope in this application, the simulation can be made closer to the real situation, avoiding unrealistic physical behaviors caused by excessive speed, thereby improving the stability and credibility of the simulation.

[0021] In the said S01: In some embodiments, the continuity equation of the kinematic wave equation for simulating the landslide movement is: ; In the formula, is the landslide depth, unit: ; is the movement speed of the landslide along the direction, unit: ; is the movement speed of the landslide along the direction, unit: ; is the time of the landslide movement, unit: .

[0022] In the said S02: In some embodiments, the terrain slope includes the terrain slope along the x direction and the terrain slope along the y direction.

[0023] In the said S03: Furthermore, the friction model is used to calculate the landslide movement speed: ; ; In the formula, is the landslide depth, unit: ; is the depth of landslide accumulation under constant conditions, unit: ; is the movement speed of the landslide along direction, unit: ; is the movement speed of the landslide along direction, unit: ; is the terrain slope along the x direction, unit: dimensionless; is the terrain slope along the y direction, unit: dimensionless; is the flow efficiency, unit: .

[0024] Furthermore, The calculation formula of is: , g is the acceleration of gravity, unit: .

[0025] Furthermore, the depth of landslide accumulation under constant conditions is: ; In the formula, is the characteristic length, unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless.

[0026] It can be understood that is the lower limit value of the friction coefficient, representing the friction characteristics of granular materials in a static or extremely slow flow state, and plays an important role in determining the yield criterion and flow initiation conditions of granular flows.

[0027] It can be understood that reflects the friction characteristics of granular materials at high shear rates or high inertia states, and together with defines the range of variation of the friction coefficient with the inertia number.

[0028] In the above S04: In some embodiments, the terrain slope is specially processed for the landslide movement speed as follows: ; In the formula, is the landslide movement speed, , unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia number, unit: dimensionless; is a constant, and its value is 0.0001, unit: dimensionless.

[0029] In a second aspect, an embodiment of the present application provides a landslide simulation system based on the kinematic wave equation, including: A module for simulating landslide movement, which simulates landslide movement using the kinematic wave equation; A parameter acquisition module for obtaining the parameters required for simulation, where the parameters include the terrain slope; A landslide movement speed calculation model for calculating the landslide movement speed using a friction model; A special processing model for specially processing the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0030] In the module for simulating landslide movement: In some embodiments, the continuity equation of the kinematic wave equation for simulating landslide movement is: ; In the formula, is the landslide depth, unit: ; is the movement speed of the landslide along the direction, unit: ; is the movement speed of the landslide along the direction, unit: ; is the time of landslide movement, unit: .

[0031] In the parameter acquisition module: In some embodiments, the terrain slope includes the terrain slope along the x direction and the terrain slope along the y direction.

[0032] In the landslide movement speed calculation model: Furthermore, the friction model is adopted to calculate the landslide movement speed: ; ; In the formula, is the landslide depth, unit: ; is the depth of landslide accumulation under constant conditions, unit: ; is the movement speed of the landslide along the direction, unit: ; is the movement speed of the landslide along the direction, unit: ; is the terrain slope along the x direction, unit: dimensionless; is the terrain slope along the y direction, unit: dimensionless; is the flow efficiency, unit: .

[0033] Furthermore, the calculation formula of is: , g is the acceleration of gravity, unit: .

[0034] Furthermore, the depth of landslide accumulation under constant conditions is: ; In the formula, is the characteristic length, unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia number, unit: dimensionless.

[0035] It can be understood that is the lower limit value of the friction coefficient, representing the friction characteristics of granular materials in a static or extremely slow flow state, and plays an important role in determining the yield criterion and flow initiation conditions of granular flows.

[0036] It can be understood that reflects the friction characteristics of granular materials under high shear rates or high inertia conditions, and together define the range of variation of the friction coefficient with the inertia number.

[0037] In the special treatment model: In some embodiments, the terrain slope has a special treatment on the landslide movement speed as follows: ; where, is the landslide movement speed, , unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless; is a constant, whose value is 0.0001, unit: dimensionless.

[0038] In a third aspect, the present application provides a computer device, including a memory and a processor. When the computer program stored in the memory is executed by the processor, the processor is caused to execute the steps of the landslide simulation method based on the kinematic wave equation as described above.

[0039] Among them, the computer device can be a computing device such as a desktop computer, a notebook, a palm computer, and a cloud server. The computer device can perform human-computer interaction with the user through a keyboard, a mouse, a remote control, a touchpad, or a voice control device, etc.

[0040] The memory at least includes one type of readable storage medium, and the readable storage medium includes flash memory, hard disk, multimedia card, card-type memory (such as SD or D-interface display memory, etc.), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, magnetic disk, optical disc, etc. In some embodiments, the memory may be an internal storage unit of the computer device, such as the hard disk or memory of the computer device. In other embodiments, the memory may also be an external storage device of the computer device, such as a plug-in hard disk, SmartMedia Card (SMC), Secure Digital (SD) card, Flash Card, etc. equipped on the computer device. Of course, the memory may also include both the internal storage unit and the external storage device of the computer device. In this embodiment, the memory is commonly used to store the operating system installed on the computer device and various application software, such as the program code of the landslide simulation method based on the kinematic wave equation. In addition, the memory may also be used to temporarily store various data that have been output or will be output.

[0041] In some embodiments, the processor may be a central processing unit (CPU), controller, microcontroller, microprocessor, or other data processing chip. The processor is generally used to control the overall operation of the computer device. In this embodiment, the processor is used to run the program code stored in the memory or process data, such as running the program code of the landslide simulation method based on the kinematic wave equation.

[0042] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium storing a computer program, and when the computer program is executed by a processor, the processor is caused to execute the steps of the landslide simulation method based on the kinematic wave equation as described above.

[0043] Wherein, the computer-readable storage medium stores an interface display program, and the interface display program can be executed by at least one processor to cause the at least one processor to execute the steps of the landslide simulation method based on the kinematic wave equation as described above.

[0044] Through the description of the above embodiments, those skilled in the art can clearly understand that the above-described embodiment methods can be implemented by means of software plus a necessary general hardware platform. Of course, they can also be implemented by hardware, but in many cases, the former is a better implementation method. Based on such an understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) and includes several instructions for causing a terminal device (which can be a mobile phone, computer, server, or network device, etc.) to execute the landslide simulation method based on the kinematic wave equation described in the embodiments of the present application.

[0045] The invention will be further described below in conjunction with application cases.

[0046] Application Example Please refer to Figure 2 , taking a landslide in a certain place as an example, and analyzing it using the landslide simulation method based on the kinematic wave equation of the present application.

[0047] In this application example, the gravitational acceleration g is taken as 9.81 , the characteristic length is 0.2 , the critical value of the friction coefficient at zero shear rate = , the limiting value of the friction coefficient at high inertia number = , is taken as 0.0001, the grid size is , , the calculation time step is . Under the above parameter conditions, the accumulation depth distribution of the landslide movement is successfully simulated. The simulation results are as shown in Figure 2 , Figure 2 showing the accumulation situation of the landslide at different times (unit: ). In Figure 2 , the red-yellow band is used to visually display the accumulation depth (unit: m). Among them, the yellow area represents a smaller accumulation depth, and the red area represents a larger accumulation depth. The color gradient intuitively reflects the spatial variation law of the accumulation depth.

[0048] The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent replacements, and improvements made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims

1. Landslide simulation method based on kinematic wave equation, characterized in that, Including the following steps: Simulating the landslide movement using the kinematic wave equation; Obtaining the parameters required for the simulation, where the parameters include the terrain slope; Calculating the landslide movement speed using a friction model; Performing special processing on the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

2. The landslide simulation method based on the kinematic wave equation according to claim 1, characterized in that The continuity equation of the kinematic wave equation for simulating the landslide movement is: ； In the formula, is the landslide depth, unit: ; is the moving speed of the landslide along direction, unit: ; is the movement speed of the landslide along direction, unit: ; is the time of landslide movement, unit: .

3. The landslide simulation method based on the kinematic wave equation according to claim 1, characterized in that The terrain slope includes the terrain slope in the x direction and the terrain slope in the y direction .

4. The landslide simulation method based on the kinematic wave equation according to claim 1, characterized in that, Adopt The friction model is used to calculate the landslide movement speed: ； ； In the formula, is the landslide depth, unit: ; is the depth of landslide accumulation under constant conditions, unit: ; is the movement speed of the landslide along direction, unit: ; is the movement speed of the landslide along direction, unit: ; is the terrain slope in the x direction, unit: dimensionless; is the terrain slope in the y direction, unit: dimensionless; is the flow efficiency, unit: .

5. The landslide simulation method based on the kinematic wave equation according to claim 4, wherein The calculation formula of is: , where g is the acceleration due to gravity, unit:

6. The landslide simulation method based on the kinematic wave equation according to claim 4, characterized in that Depth of landslide accumulation under constant conditions is as follows: ； In the formula, is the characteristic length, unit: ; is the terrain slope, , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless.

7. The landslide simulation method based on the kinematic wave equation according to claim 1, characterized in that Terrain slope on landslide movement speed Special treatment is carried out as follows: ； In the formula, is the landslide movement speed, , unit: ; is the terrain slope , unit: dimensionless; is the critical value of the friction coefficient at zero shear rate, unit: dimensionless; is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless; is a constant with a value of 0.0001 and the unit: dimensionless.

8. Landslide simulation system based on kinematic wave equation, characterized in that, Including: A module for simulating the landslide movement, which simulates the landslide movement using the kinematic wave equation; A parameter obtaining module, which obtains the parameters required for the simulation, where the parameters include the terrain slope; A landslide movement speed calculation model, which calculates the landslide movement speed using a friction model; A special processing model, which performs special processing on the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

9. Computer device, characterized in that, Including a memory and a processor, where the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, Storing a computer program, and when the computer program is executed by a processor, the processor executes the steps of the method according to any one of claims 1-7.

Images