Landslide Simulation Method, System, Device and Storage Medium Based on Kinematic Wave Equation

Abstract

The present invention discloses a landslide simulation method, system, device and storage medium based on the kinematic wave equation, belonging to the technical field of landslide simulation. The method includes: 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; and performing special processing on the landslide movement speed through the terrain slope. Based on the terrain slope data obtained by remote sensing, this application simulates the landslide movement using the kinematic wave equation, obtains the parameter required for the simulation - the terrain slope, calculates the landslide movement speed using a friction model, and performs special processing on the speed through the terrain slope, which can improve the stability of the simulation to simulate the landslide accumulation depth and flow state at different times. By combining the kinematic wave theory and the friction model, this application can effectively simulate the movement characteristics of landslides under different conditions through accurate calculation, providing a scientific basis for disaster prevention and mitigation of landslides.

Description

Technical Field

[0001] The present invention belongs to the technical field of landslide simulation, and in particular 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 relatively fast 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, which limits 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:

[0005] Simulate the landslide movement using the kinematic wave equation;

[0006] Obtain the parameters required for the simulation, and the parameters include the terrain slope;

[0007] Calculate the landslide movement speed using the friction model;

[0008] Perform special processing on the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0009] Optionally, in some embodiments of the present application, the continuity equation of the kinematic wave equation for simulating the landslide movement is:

[0010] ;

[0011] In the formula, is the landslide depth, unit: ;

[0012] is the landslide movement speed along the direction, unit: ;

[0013] is the movement speed of the landslide along the direction, unit: ;

[0014] is the time of the landslide movement, unit: .

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

[0016] Optionally, in some embodiments of the present application, the friction model is used to calculate the landslide movement speed:

[0017] ;

[0018] ;

[0019] In the formula, is the landslide depth, unit: ;

[0020] is the depth of the landslide accumulation under constant conditions, unit: ;

[0021] is the movement speed of the landslide along the direction, unit: ;

[0022] is the movement speed of the landslide along the direction, unit: ;

[0023] is the terrain slope along the x direction, unit: dimensionless;

[0024] is the terrain slope along the y direction, unit: dimensionless;

[0025] is the flow efficiency, unit: .

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

[0027] Optionally, in some embodiments of the present application, the depth of the landslide accumulation under constant conditions is:

[0028] ;

[0029] In the formula, is the characteristic length, unit: ;

[0030] is the terrain slope, , unit: dimensionless;

[0031] is the critical value of the friction coefficient at zero shear rate, unit: dimensionless;

[0032] is the limit value of the friction coefficient at high inertia number, unit: dimensionless.

[0033] Optionally, in some embodiments of the present application, the terrain slope is specially processed for the landslide movement speed as follows:

[0034] ;

[0035] In the formula, is the landslide movement speed, , unit: ;

[0036] is the terrain slope, , unit: dimensionless;

[0037] is the critical value of the friction coefficient at zero shear rate, unit: dimensionless;

[0038] is the limit value of the friction coefficient at high inertia number, unit: dimensionless;

[0039] is a constant, its value is 0.0001, unit: dimensionless.

[0040] Correspondingly, the embodiment of the present application also provides a landslide simulation system based on the kinematic wave equation, including:

[0041] a module for simulating the landslide movement, which simulates the landslide movement using the kinematic wave equation;

[0042] a parameter acquisition module for acquiring the parameters required for the simulation, and the parameters include the terrain slope;

[0043] a landslide movement speed calculation model for calculating the landslide movement speed using the friction model;

[0044] A special processing model that specially processes the landslide movement speed based on the terrain slope to simulate the landslide accumulation depth and flow state.

[0045] Correspondingly, the embodiment of the present application also provides a computer device, including a storage and a processor. The storage stores a computer program, and when the computer program is executed by the processor, the processor is caused to execute the steps of the above method.

[0046] Correspondingly, the embodiment of the present application also 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 above method.

[0047] In summary, due to the adoption of the above technical solutions, the beneficial effects of the present invention are as follows:

[0048] Based on the terrain slope data obtained by remote sensing, the present application uses the kinematic wave equation to simulate the landslide movement, obtains the parameter required for simulation - the terrain slope, calculates the landslide movement speed using the friction model, and specially processes the speed through the terrain slope, which can improve the stability of the simulation to simulate the landslide accumulation depth and flow state at different times; the present 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 disaster prevention and mitigation of landslides. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 is a flowchart of the landslide simulation method based on the kinematic wave equation of the present invention;

[0050] 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

[0051] 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.

[0052] The technical solution of the present application is as follows:

[0053] Please refer to Figure 1 , in the first aspect, the embodiment of the present application provides a landslide simulation method based on the kinematic wave equation, including the following steps:

[0054] S01. Simulate the landslide movement using the kinematic wave equation;

[0055] S02. Obtain the parameters required for simulation, where the parameters include terrain slope;

[0056] S03. Calculate the landslide movement speed using a friction model;

[0057] S04. Perform special processing on the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0058] This application is based on the terrain slope data obtained by remote sensing, uses the kinematic wave equation to simulate the landslide movement, obtains the parameter - terrain slope required for simulation, calculates the landslide movement speed using a friction model, and performs special processing on the speed through the terrain slope, which can improve the stability of the simulation 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.

[0059] According to the proposed landslide simulation method of this application, key parameters such as the movement speed and accumulation depth of the landslide can be obtained, and the corresponding departments can make advance disaster response measures 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 retaining high accuracy, providing technical support for regional landslide risk management and early warning system construction.

[0060] In the real world, the influence of terrain slope on speed is usually continuous and smooth. By performing special processing on the landslide movement speed through the terrain slope, this application can make the simulation closer to the real situation, avoid unrealistic physical behaviors caused by excessive speed, and thus improve the stability and credibility of the simulation.

[0061] In the said S01:

[0062] In some embodiments, the continuity equation of the kinematic wave equation for simulating the landslide movement is:

[0063] ;

[0064] In the formula, is the landslide depth, unit: ;

[0065] is the movement speed of the landslide along the direction, unit: ;

[0066] is the movement speed of the landslide along direction, unit: ;

[0067] is the time of the landslide movement, unit: .

[0068] In the above S02:

[0069] In some embodiments, the terrain slope includes the terrain slope along the x direction and the terrain slope along the y direction .

[0070] In the above S03:

[0071] Furthermore, the friction model is used to calculate the landslide movement speed:

[0072] ;

[0073] ;

[0074] In the formula, is the landslide depth, unit: ;

[0075] is the depth of the landslide accumulation under constant conditions, unit: ;

[0076] is the movement speed of the landslide along direction, unit: ;

[0077] is the movement speed of the landslide along direction, unit: ;

[0078] is the terrain slope along the x direction, unit: dimensionless;

[0079] is the terrain slope along the y direction, unit: dimensionless;

[0080] is the flow efficiency, unit: .

[0081] Even further, the calculation formula of is: , where g is the acceleration due to gravity, unit:

[0082] Furthermore, the depth of landslide accumulation under constant conditions is as follows:

[0083] ;

[0084] In the formula, is the characteristic length, unit: ;

[0085] is the terrain slope, , unit: dimensionless;

[0086] [[ID=DNA]]is the critical value of the friction coefficient at zero shear rate, unit: dimensionless;

[0087] is the limit value of the friction coefficient at high inertia number, unit: dimensionless.

[0088] 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 condition of granular flow.

[0089] 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. [[ID=4l]]

[0090] In the said S04:

[0091] In some embodiments, the terrain slope is specially processed for the landslide movement speed as follows:

[0092] ;

[0093] In the formula, is the landslide movement speed, ;unit: ;

[0094] is the terrain slope, , unit: dimensionless;

[0095] is the critical value of the friction coefficient at zero shear rate, unit: dimensionless;

[0096] is the limit value of the friction coefficient at high inertia number, unit: dimensionless;

[0097] is a constant, and its value is 0.0001, unit: dimensionless.

[0098] In a second aspect, an embodiment of the present application provides a landslide simulation system based on the kinematic wave equation, including:

[0099] A module for simulating landslide movement, which simulates landslide movement using the kinematic wave equation;

[0100] A parameter acquisition module for acquiring the parameters required for the simulation, where the parameters include the terrain slope;

[0101] A landslide movement speed calculation model for calculating the landslide movement speed using a friction model;

[0102] A special processing model for specially processing the landslide movement speed through the terrain slope to simulate the landslide accumulation depth and flow state.

[0103] In the module for simulating landslide movement:

[0104] In some embodiments, the continuity equation of the kinematic wave equation for simulating landslide movement is:

[0105] ;

[0106] In the formula, is the landslide depth, unit: ;

[0107] is the movement speed of the landslide along the direction, unit: ;

[0108] is the movement speed of the landslide along the direction, unit: ;

[0109] is the time of landslide movement, unit: .

[0110] In the parameter acquisition module:

[0111] In some embodiments, the terrain slope includes the terrain slope along the x direction and the terrain slope along the y direction .

[0112] In the landslide movement speed calculation model:

[0113] Furthermore, the friction model is used to calculate the landslide movement speed:

[0114] ;

[0115] ;

[0116] In the formula, is the landslide depth, unit: ;

[0117] is the depth of landslide accumulation under constant conditions, unit: ;

[0118] is the movement speed of the landslide along the direction, unit: ;

[0119] is the movement speed of the landslide along the direction, unit: ;

[0120] is the terrain slope along the x direction, unit: dimensionless;

[0121] is the terrain slope along the y direction, unit: dimensionless;

[0122] is the flow efficiency, unit: .

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

[0124] Furthermore, the depth of landslide accumulation under constant conditions is:

[0125] ;

[0126] In the formula, is the characteristic length, unit: ;

[0127] is the terrain slope, , unit: dimensionless;

[0128] is the critical value of the friction coefficient at zero shear rate, unit: dimensionless;

[0129] is the limit value of the friction coefficient at high inertia numbers, unit: dimensionless.

[0130] 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 condition of granular flow.

[0131] It can be understood that reflects the friction characteristics of granular materials under high shear rate or high inertia state, and together with defines the range of variation of the friction coefficient with the inertia number.

[0132] In the special processing model:

[0133] In some embodiments, the terrain slope has a special treatment on the landslide movement speed as follows:

[0134] ;

[0135] In the formula, is the landslide movement speed, , unit: ;

[0136] is the terrain slope, , unit: dimensionless;

[0137] is the critical value of the friction coefficient at zero shear rate, unit: dimensionless;

[0138] is the limit value of the friction coefficient at high inertia number, unit: dimensionless;

[0139] is a constant, and its value is 0.0001, unit: dimensionless.

[0140] In a third aspect, the present application provides 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 landslide simulation method based on the kinematic wave equation as described above.

[0141] Among them, the computer device can be a desktop computer, a notebook, a palm computer, a cloud server and other computing devices. 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 and other means.

[0142] 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 and various application software installed on the computer device, 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.

[0143] 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.

[0144] In a fourth aspect, 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.

[0145] 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.

[0146] Through the description of the above embodiments, those skilled in the art can clearly understand that the above-described example 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. 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.

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

[0148] Application Example

[0149] Please refer to Figure 2 , taking a landslide in a certain place as an example, the landslide simulation method based on the kinematic wave equation of the present application is used for analysis.

[0150] 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 limit 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.

[0151] 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 It 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; 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 direction, unit: ; is the movement speed of the landslide along direction, unit: ; is the time of landslide movement, unit: ; Adopt The friction model is used to calculate the landslide movement speed: ； ； In the formula, is the landslide depth, unit: ; is the depth of the 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: ; 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.

2. 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 .

3. The landslide simulation method based on the kinematic wave equation according to claim 1, characterized in that The calculation formula is: , where g is the acceleration due to gravity, unit: .

4. The landslide simulation method based on the kinematic wave equation according to claim 1, 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.

5. Landslide simulation system based on kinematic wave equation, characterized in that, It includes: 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, and the parameters include the terrain slope; A landslide movement speed calculation model, which calculates the landslide movement speed using the 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; 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 direction, unit: ; is the movement speed of the landslide along direction, unit: ; The time of landslide movement, unit: ; 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: ; 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.

6. A computer device, characterized in that, It includes a memory and a processor. The memory stores a computer program. When the computer program is executed by the processor, the processor executes the steps of the method according to any one of claims 1-4.

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

Images