Method, system and equipment for calculating erosion rate in debris flow movement process and storage medium

Abstract

The invention discloses a method, system and device for calculating the erosion rate in the debris flow movement process and a storage medium, and belongs to the field of debris flow prevention and control engineering, and the method comprises the steps: obtaining debris flow fluid characteristic parameters, carrying out the dimensionless analysis according to the debris flow fluid characteristic parameters, and calculating the Froude number, the Savage number and the Bagnold number of the dimensionless parameters; feature parameters of the debris flow passing through the channel are obtained, and a calculation model of the dimensionless erosion rate of the debris flow in the channel is constructed according to the feature parameters of the debris flow fluid and the relation between the feature parameters of the debris flow passing through the channel and the erosion rate; and solving the calculation model of the dimensionless erosion rate. The method can meet the requirements of actual debris flow engineering design, has high feasibility and practicability, and can more accurately predict the debris flow erosion condition, thereby more effectively guiding the debris flow prevention and treatment engineering design and providing technical support for disaster reduction and prevention.

Description

Technical Field

[0001] The present invention belongs to the field of debris flow prevention and control engineering, and in particular relates to a method, system, equipment and storage medium for calculating the erosion rate during the movement of debris flow. Background Art

[0002] Debris flow is a high-speed, highly destructive natural disaster, often accompanied by mountain torrents, landslides, rock collapses and other disasters, posing a huge threat to people's lives and property. In recent years, the research on debris flow erosion at home and abroad has focused on comprehensively considering factors such as the properties of debris flow fluids and the composition of gully bed materials. The trend is to explore the relationship and mechanism of these influencing factors in erosion characteristics. For example, Smerdon and Beasley et al. believed that the shear strength of soil is related to the median particle size. d 50 There is an exponential relationship. The debris flow erosion capacity and shear resistance are used to judge the occurrence of debris flow erosion. The fine particle content also has a significant effect on soil strength. Iverson USGS Large-scale debris flow experiments have shown that the water content of the bed material has a positive feedback effect on the erosion process and momentum growth. A high-water content bed will increase pore water pressure and promote further erosion, while a low-water content bed will produce a negative feedback. The loss of soil strength in the ditch bed caused by rapid undrained loading is also a factor in the instability of the ditch bed soil.

[0003] At present, the feasibility, practicality and accuracy of predicting debris flow erosion need to be further improved. Summary of the invention

[0004] The purpose of the present invention is to provide a method, system, device and storage medium for calculating the erosion rate during debris flow movement, so as to solve the problem that the feasibility, practicality and accuracy of predicting debris flow erosion need to be further improved.

[0005] The embodiment of the present application is implemented as follows: a method for calculating the erosion rate during debris flow movement includes: obtaining characteristic parameters of debris flow fluid, performing dimensionless analysis based on the characteristic parameters of debris flow fluid, and calculating dimensionless parameters Froude number, Savage number, and Bagnold number; Obtain characteristic parameters of debris flow passing through the channel, and build a dimensionless erosion rate calculation model of debris flow in the channel according to the relationship between characteristic parameters of debris flow fluid and characteristic parameters of debris flow passing through the channel and erosion rate; Solve the calculation model of dimensionless erosion rate and obtain the erosion rate.

[0006] Optionally, in some embodiments of the present application, the calculation formulas of the dimensionless parameters of debris flow, Froude number, Savage number, and Bagnold number, are as follows: ; In the formula, is the Froude number, is a dimensionless parameter, is the debris flow velocity; For debris flow depth; , is the acceleration due to gravity; ; In the formula, is the Savage number, a dimensionless parameter, , is the density of the particle, is the density of water, is the diameter of the particle, is the average shear rate of debris flow, is the debris flow velocity, , is the acceleration due to gravity, The mudslide is deep. is the channel slope, =30 is the friction angle between particles; ; In the formula, is the Bagnold number, which is a dimensionless parameter; is the density of the particle; is the average shear rate of debris flow; is the debris flow velocity; For debris flow depth; is the diameter of the particle; is the volume fraction of solid matter in debris fluid, which is a dimensionless parameter; is the density of debris flow; is the density of the particle; is the density of water; It is the dynamic viscosity system of the liquid phase slurry of debris flow, with the dilution value being 0.03 Pa·s, the transition value being 0.05 Pa·s, and the viscosity value being 0.1 Pa·s.

[0007] Optionally, in some embodiments of the present application, the debris flow fluid characteristic parameters include debris flow velocity, debris flow depth and debris flow density.

[0008] Optionally, in some embodiments of the present application, characteristic parameters of debris flow through the channel include the volume moisture content of the bottom bed and the volume percentage of fine particles in the bottom bed.

[0009] Optionally, in some embodiments of the present application, the calculation model of the dimensionless erosion rate of debris flow in the channel is: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow; is the density of the particle; is the density of water; ; ; ; ; ; ; .

[0010] Optionally, in some embodiments of the present application, the dimensionless erosion rate Calculate the erosion rate , the calculation formula is as follows: ; Where, E is the erosion rate; is the dimensionless erosion rate, is a dimensionless parameter; g = 9.81 is the gravitational acceleration; is the channel slope; To correspond to the depth of debris flow.

[0011] Accordingly, the embodiment of the present application also provides a system for calculating the erosion rate during debris flow movement, including: a dimensionless analysis module for obtaining characteristic parameters of debris flow fluid, performing dimensionless analysis based on the characteristic parameters of debris flow fluid, and calculating dimensionless parameters Froude number, Savage number, and Bagnold number; Construct a dimensionless erosion rate calculation model module to obtain characteristic parameters of debris flow passing through the channel, and construct a dimensionless erosion rate calculation model of debris flow in the channel according to the relationship between characteristic parameters of debris flow fluid and characteristic parameters of debris flow passing through the channel and erosion rate; The calculation model module for solving the dimensionless erosion rate is used to solve the calculation model of the dimensionless erosion rate and obtain the erosion rate.

[0012] Accordingly, an embodiment of the present application also provides a computer device, including a storage device and a processor, wherein the storage device stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the above method.

[0013] Accordingly, an 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 executes the steps of the above method.

[0014] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are: In addition to considering the characteristic parameters of debris flow fluid, this application also considers the influence of the debris flow through the channel characteristic parameters on the erosion rate of the debris flow, and proposes a method for calculating the erosion rate of the debris flow in the channel under the influence of the debris flow fluid and the debris flow channel on the debris flow erosion capacity. The method of this application is based on the derivation of empirical formulas to solve the influence of the channel bed and the movement characteristics of the debris flow on the debris flow erosion rate, and can obtain the parameters required for calculation through field investigation, sampling and experimental analysis, and can reasonably determine the influence of the channel characteristics on the debris flow erosion and calculate the erosion rate during its movement, providing a basis for the design of debris flow prevention and control projects. In addition, the method of this application can meet the needs of the actual engineering design of debris flow, and has high feasibility and practicality. The debris flow erosion rate calculation method provided by this application can more accurately predict the debris flow erosion situation, thereby more effectively guiding the design of debris flow prevention and control projects, and providing technical support for disaster reduction and prevention. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 The present invention is a flow chart of a method for calculating the erosion rate during debris flow movement. DETAILED DESCRIPTION

[0016] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with 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: See also Figure 1 The present application embodiment provides a method for calculating the erosion rate during debris flow movement, including: S01. Obtain characteristic parameters of debris flow fluid, perform dimensionless analysis based on the characteristic parameters of debris flow fluid, and calculate dimensionless parameters such as Froude number, Savage number, and Bagnold number; S02, obtaining characteristic parameters of debris flow passing through the channel, and constructing a calculation model of the dimensionless erosion rate of the debris flow in the channel according to the relationship between the characteristic parameters of the debris flow fluid and the characteristic parameters of the debris flow passing through the channel and the erosion rate; S03. Solve the calculation model of the dimensionless erosion rate to obtain the erosion rate.

[0018] The debris flow erosion formula proposed in the prior art needs to consider the calculation of the shear force between the debris flow and the channel and the shear strength of the main channel. Since debris flows often occur in high mountain canyons, the shear capacity of debris flows and the shear strength between channels are difficult to measure. The debris flow movement-erosion process is a dynamic result of the combined action of debris flow fluid characteristics and the physical conditions of the bottom bed channel. In addition to considering the characteristic parameters of the debris flow fluid, the present application also considers the influence of the debris flow on the erosion rate of the debris flow under the combined action of the channel characteristic parameters, and proposes a method for calculating the erosion rate of the debris flow in the channel under the influence of the debris flow fluid and the debris flow channel on the debris flow erosion capacity. The method of the present application is based on the derivation of empirical formulas to solve the influence of the channel bottom bed and the debris flow movement characteristics on the debris flow erosion rate, and can obtain the parameters required for calculation through field investigation, sampling and experimental analysis, can reasonably determine the influence of the channel characteristics on the debris flow erosion and calculate its erosion rate during movement, and provide a basis for the design of debris flow prevention and control projects. In addition, the method of the present application can meet the needs of actual engineering design of debris flow, and has high feasibility and practicality. The debris flow erosion rate calculation method provided in the present application can more accurately predict the debris flow erosion situation, thereby more effectively guiding the design of debris flow prevention and control projects, and providing technical support for achieving disaster reduction and prevention.

[0019] The method of the present application can be used to predict the damage and erosion of the channel bottom bed when debris flow moves in the channel. The debris flow erosion rate can be obtained by measuring the dynamic parameters of the debris flow during the movement of the debris flow. At the same time, the debris flow erosion rate calculation method provided by the present application can also provide a scientific basis and technical support for debris flow prevention and control, which has important theoretical significance and practical application value.

[0020] In the S01: In some embodiments, the debris flow fluid characteristic parameters include the debris flow velocity , the depth of the debris flow and debris flow density .

[0021] In some embodiments, debris flow fluid characteristic parameters are obtained through field surveys.

[0022] In some embodiments, the calculation formulas of the dimensionless parameters of debris flow, Froude number, Savage number, and Bagnold number are as follows: ; In the formula, is the Froude number, which is a dimensionless parameter; is the debris flow velocity, dimension: m / s; is the depth of debris flow, dimension: m; , is the acceleration due to gravity, dimension: m / s2 .

[0023] ; In the formula, is the Savage number, which is a dimensionless parameter; , is the density of the particle, dimension: kg / m 3 ; , is the density of water, dimension: kg / m 3 ; is the diameter of the particle, dimension: m; is the average shear rate of debris flow, dimension: s -1 ; is the debris flow velocity, dimension: m / s; , is the acceleration due to gravity, dimension: m / s 2 ; is the depth of debris flow, dimension: m; is the channel slope, dimension: º; =30 is the friction angle between particles, dimension: º; ; In the formula, is the Bagnold number, which is a dimensionless parameter; is the density of the particle, dimension: kg / m 3 ; is the average shear rate of debris flow, dimension: s -1 ; is the debris flow velocity, dimension: m / s; is the depth of debris flow, dimension: m; is the diameter of the particle, dimension: m; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; is the dynamic viscosity system of the liquid phase slurry of debris flow, dimension: Pa·s, the dilution value is 0.03Pa·s, the transition value is 0.05 Pa·s, and the viscosity value is 0.1 Pa·s.

[0024] In the S02: In some embodiments, the characteristic parameters of the debris flow through the channel include the volumetric water content of the bottom bed. and volume percentage of fine particles in the bottom bed .

[0025] In some embodiments, the dimensionless erosion rate of debris flow in the channel is calculated as: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0026] In said S03: In some embodiments, taking logarithms on both sides of the equation of the dimensionless erosion rate calculation model, we obtain: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0027] Based on the properties of logarithms, take the logarithm of both sides of the equation , and the calculation model of debris flow erosion rate is obtained: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; e = 2.718 is the base of natural logarithms.

[0028] It can be understood that based on the properties of logarithms, taking logarithms on both sides of the equation , a debris flow erosion rate calculation model that comprehensively considers the debris flow fluid properties and channel physical characteristics is obtained.

[0029] Furthermore, the dimensionless erosion rate Calculate the erosion rate , the calculation formula is as follows: ; Where, E is the erosion rate, dimension: m / s; is the dimensionless erosion rate, which is a dimensionless parameter; g = 9.81 is the gravitational acceleration, dimension: m / s 2 ; is the channel slope, dimension: º; Corresponding to the depth of debris flow, dimension: m.

[0030] In a second aspect, an embodiment of the present application provides a system for calculating the erosion rate during debris flow movement, including: The dimensionless analysis module is used to obtain the characteristic parameters of debris flow fluid, perform dimensionless analysis based on the characteristic parameters of debris flow fluid, and calculate the dimensionless parameters Froude number, Savage number, and Bagnold number; Construct a dimensionless erosion rate calculation model module to obtain characteristic parameters of debris flow passing through the channel, and construct a dimensionless erosion rate calculation model of debris flow in the channel according to the relationship between characteristic parameters of debris flow fluid and characteristic parameters of debris flow passing through the channel and erosion rate; The calculation model module for solving the dimensionless erosion rate is used to solve the calculation model of the dimensionless erosion rate and obtain the erosion rate.

[0031] In the dimensionless analysis module: In some embodiments, the debris flow fluid characteristic parameters include the debris flow velocity , the depth of the debris flow and debris flow density .

[0032] In some embodiments, debris flow fluid characteristic parameters are obtained through field surveys.

[0033] In some embodiments, the calculation formulas of the dimensionless parameters of debris flow, Froude number, Savage number, and Bagnold number are as follows: ; In the formula, is the Froude number, which is a dimensionless parameter; is the debris flow velocity (dimension: m / s); is the depth of debris flow (dimension: m); , is the acceleration due to gravity (dimension: m / s 2 ); ; In the formula, is the Savage number, which is a dimensionless parameter; , is the density of the particle, dimension: kg / m 3 ; , is the density of water, dimension: kg / m 3 ; is the diameter of the particle, dimension: m; is the average shear rate of debris flow, dimension: s -1 ; is the debris flow velocity (dimension: m / s); , is the acceleration due to gravity (dimension: m / s 2 ); is the depth of debris flow, dimension: m; is the channel slope, dimension: º; =30 is the friction angle between particles, dimension: º; ; In the formula, is the Bagnold number, which is a dimensionless parameter; is the density of the particle, dimension: kg / m 3 ; is the average shear rate of debris flow, dimension: s -1 ; is the debris flow velocity, dimension: m / s; is the depth of debris flow, dimension: m; is the diameter of the particle, dimension: m; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; is the dynamic viscosity system of the liquid phase slurry of debris flow, dimension: Pa·s, the dilution value is 0.03Pa·s, the transition value is 0.05 Pa·s, and the viscosity value is 0.1 Pa·s.

[0034] In the module for constructing the calculation model of dimensionless erosion rate: In some embodiments, the characteristic parameters of the debris flow through the channel include the bottom water content and the fine particle content of the bottom bed .

[0035] In some embodiments, the dimensionless erosion rate of debris flow in the channel is calculated as: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0036] In the calculation model module for solving the dimensionless erosion rate: In some embodiments, taking logarithms on both sides of the equation of the dimensionless erosion rate calculation model, we obtain: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0037] Based on the properties of logarithms, take the logarithm of both sides of the equation , and the calculation model of debris flow erosion rate is obtained: ; In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the density of debris flow, dimension: kg / m 3 ; is the density of the particle, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; e = 2.718 is the base of natural logarithms.

[0038] It can be understood that based on the properties of logarithms, taking logarithms on both sides of the equation , a debris flow erosion rate calculation model that comprehensively considers the debris flow fluid properties and channel physical characteristics is obtained.

[0039] Furthermore, the dimensionless erosion rate Calculate the erosion rate , the calculation formula is as follows: ; Where, E is the erosion rate, dimension: m / s; is the dimensionless erosion rate, which is a dimensionless parameter; g = 9.81 is the gravitational acceleration, dimension: m / s 2 ; is the channel slope, dimension: º; Corresponding to the depth of debris flow, dimension: m.

[0040] In a third aspect, the present application provides a computer device including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method for calculating the erosion rate during debris flow movement as described above.

[0041] The computer device may be a desktop computer, a notebook, a PDA, a cloud server, etc. The computer device may interact with the user through a keyboard, a mouse, a remote control, a touch pad, or a voice control device.

[0042] The memory includes at least one type of readable storage medium, and the readable storage medium includes flash memory, hard disk, multimedia card, card-type memory (for example, 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, disk, optical disk, etc. In some embodiments, the memory can be an internal storage unit of the computer device, such as a hard disk or memory of the computer device. In other embodiments, the memory can also be an external storage device of the computer device, such as a plug-in hard disk equipped on the computer device, a smart memory card (SmartMediaCard, SMC), a secure digital (SecureDigital, SD) card, a flash card (FlashCard), etc. Of course, the memory can also include both the internal storage unit of the computer device and its external storage device. In this embodiment, the memory is often used to store the operating system and various application software installed on the computer device, such as the program code of the calculation method of the erosion rate during the movement of the debris flow, etc. In addition, the memory can also be used to temporarily store various types of data that have been output or are to be output.

[0043] The processor may be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chip in some embodiments. 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 method for calculating the erosion rate during the movement of the debris flow.

[0044] In a fourth aspect, the present application provides 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 method for calculating the erosion rate during debris flow movement as described above.

[0045] The computer-readable storage medium stores an interface display program, and the interface display program can be executed by at least one processor so that the at least one processor executes the steps of the method for calculating the erosion rate during debris flow movement as described above.

[0046] Through the description of the above implementation methods, those skilled in the art can clearly understand that the above-mentioned embodiment method can be implemented by means of software plus a necessary general hardware platform, and of course it can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, disk, CD), and includes a number of instructions for a terminal device (which can be a mobile phone, computer, server or network device, etc.) to execute the method for calculating the erosion rate during the movement of debris flow described in the embodiment of the present application.

[0047] The invention is further described below in conjunction with application cases.

[0048] Application Example 1 After the "5.12" Wenchuan earthquake, the number of landslide geological disasters in the Zoumalinggou area increased, exacerbating the disaster risk in the area. The accumulation of loose materials increased, and when encountering factors such as heavy rain, the scale of debris flow disasters further expanded. On August 13, 2010, a large-scale debris flow disaster occurred in Zoumalinggou, Qingping Township, Mianzhu City, Sichuan Province.

[0049] This application example is based on a detailed investigation and study of the Zoumaling debris flow disaster, and verifies the debris flow erosion rate calculation formula of this application. The specific steps are as follows: The first step is to measure the density of debris flow in the channel by actual sampling, which is 2240 m 3 / kg, and the debris flow velocity is 3.45 m / s according to the on-site data; through the field observation data, it is determined that the water content of the bottom bed of Wenjiagou is 15%, the content of fine particles in the bottom bed is 60%, and the debris flow depth is 4 m. The median particle size of debris flow particles is 75 mm, the particle density is 2680 kg / m, and the channel slope is 15°. is the dynamic viscosity system of the liquid slurry, with a value of 0.1 Pa·s. =30 is the friction angle between particles. The calculation method of dimensionless parameters Froude number, Savage number, Bagnold number and volume fraction of solid phase in debris fluid is as follows: ; ; ; ; The calculated debris flow shear rate is 1.29 s -1, the solid volume fraction is 0.74, the dimensionless parameters Froude number is 0.55, Savage number is 0.00017, and Bagnold number is 549.63; The second step is to calculate the dimensionless erosion rate of debris flow by the following formula: ; Obtain the dimensionless erosion rate of debris flow = 0.0106; The third step is to calculate the debris flow erosion rate through the dimensionless erosion rate of debris flow. The calculation formula is as follows: ; Calculate the erosion rate = 0.0538 m / s.

[0050] Application Example 2 A debris flow event occurs in a certain mountainous area, and the debris flow rushes into the river from the ditch. The debris flow blocking river determination method of the present invention is used to determine the time of this event. The specific steps are as follows: The first step is to measure the density of debris flow in the channel by actual sampling, which is 1850 m 3 / kg, and the debris flow velocity is 5 m / s according to the field data. According to the field observation data, the water content of the bottom bed of Wenjiagou is 10%, the fine particle content of the bottom bed is 50%, the debris flow depth is 3.2 m, the median particle size of the debris flow particles is 67 mm, the particle density is 2680 kg / m, and the channel slope is 15º. is the dynamic viscosity system of the liquid slurry, with a value of 0.05 Pa·s. =30 is the friction angle between particles. The calculation method of dimensionless parameters Froude number, Savage number, Bagnold number and volume fraction of solid phase in debris fluid is as follows: ; ; ; ; The calculated debris flow shear rate is 2.34 s -1 The solid volume fraction of the debris flow is 0.51, the dimensionless parameters Froude number is 0.894, Savage number is 0.000559, and Bagnold number is 577.52.

[0051] The second step is to calculate the dimensionless erosion rate of debris flow by the following formula: ; Obtain the dimensionless erosion rate of debris flow = 0.00181; The third step is to calculate the debris flow erosion rate through the dimensionless erosion rate of debris flow. The calculation formula is as follows: ; Calculate the erosion rate = 0.102 m / s.

[0052] Application Example 3 A debris flow event occurs in a mountainous area, and the debris flow rushes into the river from the ditch. The erosion rate calculation method of the debris flow movement process of the present invention is used to calculate this event. The specific steps are as follows: The first step is to measure the density of debris flow in the channel by actual sampling, which is 1550 m 3 / kg, and the debris flow velocity is 6.2 m / s according to the field data. According to the field observation data, the water content of the bottom bed of Wenjiagou is 5%, the fine particle content of the bottom bed is 40%, and the depth of the debris flow is 1.3 m / s. is 5.6 m, The median particle size of debris flow particles is 62 mm, the particle density is 2680 kg / m, and the channel slope is 15°. is the dynamic viscosity system of the liquid slurry, with a value of 0.03 Pa·s. =30 is the friction angle between particles. The calculation method of dimensionless parameters Froude number, Savage number, Bagnold number and volume fraction of solid phase in debris fluid is as follows: ; ; ; ; The calculated debris flow shear rate is 1.66 s -1 The solid volume fraction of the debris flow is 0.327, the dimensionless parameters Froude number is 0.838, Savage number is 0.00014, and Bagnold number is 277.5.

[0053] The second step is to calculate the dimensionless erosion rate of debris flow by the following formula: ; Obtain the dimensionless erosion rate of debris flow = 0.014; The third step is to calculate the debris flow erosion rate through the dimensionless erosion rate of debris flow. The calculation formula is as follows: ; Calculate the erosion rate = 0.059m / s.

[0054] 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 and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. The method for calculating the erosion rate during debris flow movement is characterized by: include: Obtain the characteristic parameters of debris flow fluid, perform dimensionless analysis based on the characteristic parameters of debris flow fluid, and calculate the dimensionless parameters Froude number, Savage number, and Bagnold number; Obtain characteristic parameters of debris flow passing through the channel, and build a dimensionless erosion rate calculation model of debris flow in the channel according to the relationship between characteristic parameters of debris flow fluid and characteristic parameters of debris flow passing through the channel and erosion rate; Solve the calculation model of dimensionless erosion rate and obtain the erosion rate.

2. The method for calculating the erosion rate during debris flow movement according to claim 1, characterized in that: The calculation formulas for the dimensionless parameters of debris flow, Froude number, Savage number, and Bagnold number, are as follows: ； In the formula, is the Froude number, is a dimensionless parameter, is the debris flow velocity; For debris flow depth; , is the acceleration due to gravity; ； In the formula, is the Savage number, a dimensionless parameter, , is the density of the particle, is the density of water, is the diameter of the particle, is the average shear rate of debris flow, is the debris flow velocity, , is the acceleration due to gravity, The mudslide is deep. is the channel slope, =30 is the friction angle between particles; ； In the formula, is the Bagnold number, which is a dimensionless parameter; is the density of the particle; is the average shear rate of debris flow; is the debris flow velocity; For debris flow depth; is the diameter of the particle; is the volume fraction of solid matter in debris fluid, which is a dimensionless parameter; is the debris flow density; is the density of the particle; is the density of water; It is the dynamic viscosity system of the liquid phase slurry of debris flow, with the dilution value being 0.03 Pa·s, the transition value being 0.05 Pa·s, and the viscosity value being 0.1 Pa·s.

3. The method for calculating the erosion rate during debris flow movement according to claim 1, characterized in that: The characteristic parameters of debris flow flow include debris flow velocity, debris flow depth and debris flow density.

4. The method for calculating the erosion rate during debris flow movement according to claim 1, characterized in that: The characteristic parameters of debris flow through the channel include the volume moisture content of the bottom bed and the volume percentage of fine particles in the bottom bed.

5. The method for calculating the erosion rate during debris flow movement according to claim 1, characterized in that: The calculation model of the dimensionless erosion rate of debris flow in the channel is: ； In the formula, is the dimensionless erosion rate, is a dimensionless parameter; is the volumetric water content of the bottom bed, which is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, which is a dimensionless parameter; is the debris flow Froude number, which is a dimensionless parameter; is the debris flow Savage number, which is a dimensionless parameter; is the Bagnold number of debris flow, which is a dimensionless parameter; is the volume fraction of solid matter in debris fluid, a dimensionless parameter; is the debris flow density; is the density of the particle; is the density of water; ; ; ; ; ; ; .

6. The method for calculating the erosion rate during debris flow movement according to claim 5, characterized in that: By dimensionless erosion rate Calculate the erosion rate , the calculation formula is as follows: ； Where, E is the erosion rate; is the dimensionless erosion rate, is a dimensionless parameter; g = 9.81 is the gravitational acceleration; is the channel slope; To correspond to the depth of debris flow.

7. A system for calculating the erosion rate during debris flow movement, characterized in that: include: The dimensionless analysis module is used to obtain the characteristic parameters of debris flow fluid, perform dimensionless analysis based on the characteristic parameters of debris flow fluid, and calculate the dimensionless parameters Froude number, Savage number, and Bagnold number; Construct a dimensionless erosion rate calculation model module to obtain characteristic parameters of debris flow passing through the channel, and construct a dimensionless erosion rate calculation model of debris flow in the channel according to the relationship between characteristic parameters of debris flow fluid and characteristic parameters of debris flow passing through the channel and erosion rate; The calculation model module for solving the dimensionless erosion rate is used to solve the calculation model of the dimensionless erosion rate and obtain the erosion rate.

8. Computer device, characterized in that The invention comprises a storage and a processor, wherein the storage 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 to 6.

9. A computer-readable storage medium, characterized in that A computer program is stored, and when the computer program is executed by a processor, the processor is caused to execute the steps of the method according to any one of claims 1 to 6.

Images