Calculation method, system, device and storage medium for erosion rate during debris flow movement

Abstract

The present invention discloses a calculation method, system, device and storage medium for the erosion rate during the movement of debris flow, belonging to the field of debris flow prevention and control engineering. The method includes obtaining the characteristic parameters of the debris flow fluid, performing dimensionless analysis based on the characteristic parameters of the debris flow fluid, and calculating the dimensionless parameters Froude number, Savage number, and Bagnold number; obtaining the characteristic parameters of the channel through which the debris flow passes, and constructing a calculation model for 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 channel through which the debris flow passes and the erosion rate; solving the calculation model of the dimensionless erosion rate. The present invention can meet the needs of actual engineering design of debris flow, has high feasibility and practicability, can more accurately predict the erosion situation of debris flow, thereby more effectively guiding the design of debris flow prevention and control engineering, and providing technical support for realizing 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;

[0006] 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;

[0007] Solve the calculation model of dimensionless erosion rate and obtain the erosion rate.

[0008] 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:

[0009] ;

[0010] In the formula, is the Froude number, a dimensionless parameter, is the debris flow velocity; is the debris flow depth; , is the acceleration due to gravity;

[0011] ;

[0012] In the formula, is the Savage number, a dimensionless parameter, , is the density of the particles, is the density of water, is the diameter of the particles, is the average shear rate of the debris flow, is the debris flow velocity, , is the acceleration due to gravity, is the debris flow depth, is the channel slope, = 30 is the friction angle between particles;

[0013] ;

[0014] In the formula, is the Bagnold number, a dimensionless parameter; is the density of the particles; is the average shear rate of the debris flow; is the debris flow velocity; is the debris flow depth; is the diameter of the particles; is the volume fraction of solid phase material in the debris flow fluid, a dimensionless parameter; is the debris flow density; is the density of the particles; is the density of water; is the dynamic viscosity coefficient of the debris flow liquid phase slurry, taking the value of 0.03 Pa·s for dilute type, 0.05 Pa·s for transitional type, and 0.1 Pa·s for viscous type.

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

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

[0017] Optionally, in some embodiments of the present application, the calculation model of the dimensionless erosion rate of debris flow in the channel is:

[0018] ;

[0019] In the formula, is the dimensionless erosion rate, and is a dimensionless parameter; is the volumetric water content of the bottom bed, and is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, and is a dimensionless parameter; is the Froude number of debris flow, and is a dimensionless parameter; is the Savage number of debris flow, and is a dimensionless parameter; is the Bagnold number of debris flow, and is a dimensionless parameter; is the volume fraction of solid phase substances in the debris flow, a dimensionless parameter; is the density of the debris flow; is the density of the particles; is the density of water; ; ; ; ; ; ; .

[0020] Optionally, in some embodiments of the present application, the erosion rate is calculated through the dimensionless erosion rate , and the calculation formula is as follows:

[0021] ;

[0022] In the formula, E is the erosion rate; is the dimensionless erosion rate, and is a dimensionless parameter; g = 9.81 is the acceleration due to gravity; is the channel slope; is the corresponding debris flow depth.

[0023] Correspondingly, the embodiment of the present application also provides a calculation system for the erosion rate during the movement of debris flow, including: a dimensionless analysis module, configured to obtain the characteristic parameters of the debris flow fluid, perform dimensionless analysis according to the characteristic parameters of the debris flow fluid, and calculate the dimensionless parameters Froude number, Savage number, and Bagnold number;

[0024] A calculation model module for constructing a dimensionless erosion rate, configured to obtain the characteristic parameters of the debris flow passing through the channel, and construct 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;

[0025] A 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.

[0026] Correspondingly, an embodiment of the present application further provides a computer device, including a storage and a processor. The storage stores a computer program. When the computer program is executed by the processor, the processor executes the steps of the above method.

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

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

[0029] In addition to considering the debris flow fluid characteristic parameters, this application also considers the influence of the debris flow passing through the channel characteristic parameters on the erosion rate of the debris flow. A calculation method for the erosion rate of the debris flow in the channel is proposed under the influence of the debris flow fluid and the debris flow channel on the erosion ability of the debris flow. The method of this application is derived based on empirical formulas to solve the influence of the channel bed and the debris flow movement characteristics on the debris flow erosion rate, and the parameters required for calculation can be obtained through field investigations, sampling, and experimental analyses. It 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 actual debris flow engineering design, has high feasibility and practicability. The debris flow erosion rate calculation method provided by this application can more accurately predict the debris flow erosion situation, thus 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

[0030] Figure 1 It is a flowchart of the calculation method for the erosion rate during the movement of the debris flow of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

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

[0032] The technical solution of this application is as follows:

[0033] Please refer to Figure 1 , an embodiment of the present application provides a calculation method for the erosion rate during the movement of the debris flow, including:

[0034] S01. Obtain the debris flow fluid characteristic parameters, conduct dimensionless analysis based on the debris flow fluid characteristic parameters, and calculate the dimensionless parameters Froude number, Savage number, and Bagnold number;

[0035] S02. Obtain the debris flow channel characteristic parameters, and construct a calculation model for the dimensionless erosion rate of debris flow in the channel based on the relationship between the debris flow fluid characteristic parameters and the debris flow channel characteristic parameters and the erosion rate;

[0036] S03. Solve the calculation model of the dimensionless erosion rate to obtain the erosion rate.

[0037] In the debris flow erosion formula proposed by the prior art, it is necessary to consider calculating 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 mountains and deep valleys, it is difficult to measure the shear ability of the debris flow and the shear strength between the channel. The process of debris flow movement - erosion is a dynamic result jointly affected by both the debris flow fluid characteristics and the physical conditions of the bottom channel. In addition to considering the debris flow fluid characteristic parameters, this application also considers the influence of the debris flow channel characteristic parameters on the erosion rate of the debris flow under their combined action, and proposes a calculation method for the erosion rate of debris flow in the channel under the influence of the debris flow fluid and the debris flow channel on the debris flow erosion ability. The method of this application is derived based on empirical formulas to solve the influence of the channel bed and the debris flow movement characteristics on the debris flow erosion rate, and can obtain the parameters required for calculation through field investigations, 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 actual debris flow engineering design, 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.

[0038] The method of this application can be used to predict the damage and erosion of the channel bed by the debris flow during its movement in the channel. The measured erosion rate of the debris flow can be obtained by measuring the dynamic parameters during the debris flow movement process; at the same time, the debris flow erosion rate calculation method provided by this application can also provide a scientific basis and technical support for debris flow prevention and control, and has important theoretical significance and practical application value.

[0039] In the above S01:

[0040] In some embodiments, the debris flow fluid characteristic parameters include the debris flow velocity , the debris flow depth and the debris flow density .

[0041] In some embodiments, the debris flow fluid characteristic parameters are obtained through field investigations.

[0042] In some embodiments, the calculation formulas for the dimensionless parameters of debris flow, namely the Froude number, Savage number, and Bagnold number, are as follows:

[0043] ;

[0044] In the formula, is the Froude number, a dimensionless parameter; is the debris flow velocity, with the dimension: m / s; is the debris flow depth, with the dimension: m; , is the acceleration due to gravity, with the dimension: m / s 2 .

[0045] ;

[0046] In the formula, is the Savage number, a dimensionless parameter; , is the density of the particles, with the dimension: kg / m 3 ; , is the density of water, with the dimension: kg / m 3 ; is the diameter of the particles, with the dimension: m; is the average shear rate of the debris flow, with the dimension: s -1 ; is the debris flow velocity, with the dimension: m / s; , is the acceleration due to gravity, with the dimension: m / s 2 ; is the debris flow depth, with the dimension: m; is the channel slope, with the dimension: º; = 30 is the inter-particle friction angle, with the dimension: º;

[0047] ;

[0048] In the formula, is the Bagnold number, a dimensionless parameter; is the density of the particles, with the dimension: kg / m 3 ;

[0049] is the average shear rate of the debris flow, with the dimension: s -1 ; is the debris flow velocity, with the dimension: m / s; is the debris flow depth, with the dimension: m; is the diameter of the particles, with the dimension: m; is the volume fraction of solid phase substances in the debris flow, a dimensionless parameter; is the density of the debris flow, dimension: kg / m 3 ; is the density of the particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; is the dynamic viscosity coefficient of the debris flow liquid-phase slurry, dimension: Pa·s, taking 0.03 Pa·s for dilute type, 0.05 Pa·s for transitional type, and 0.1 Pa·s for viscous type.

[0050] In the said S02:

[0051] In some embodiments, the channel characteristic parameters of the debris flow include the bottom bed volume water content and the bottom bed fine particle volume percentage .

[0052] In some embodiments, the calculation model of the dimensionless erosion rate of the debris flow in the channel is:

[0053] ;

[0054] In the formula, is the dimensionless erosion rate, a dimensionless parameter; is the bottom bed volume water content, a dimensionless parameter; is the bottom bed fine particle volume percentage, a dimensionless parameter; is the debris flow Froude number, a dimensionless parameter; is the debris flow Savage number, a dimensionless parameter; is the debris flow Bagnold number, a dimensionless parameter; is the volume fraction of solid phase substances in the debris flow, a dimensionless parameter; is the density of the debris flow, dimension: kg / m 3 ; is the density of the particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0055] In the said S03:

[0056] In some embodiments, taking the logarithm of both sides of the equation of the calculation model of the dimensionless erosion rate, we get:

[0057] ;

[0058] In the formula, is the dimensionless erosion rate, and is a dimensionless parameter; is the volumetric water content of the bottom bed, and is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, and is a dimensionless parameter; is the debris flow Froude number, and is a dimensionless parameter; is the debris flow Savage number, and is a dimensionless parameter; is the debris flow Bagnold number, and is a dimensionless parameter; is the volume fraction of solid phase substances in the debris flow fluid, a dimensionless parameter; is the density of the debris flow, dimension: kg / m 3 ; is the density of the particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0059] Based on the properties of logarithms, taking the logarithm of both sides of the equation , we get the calculation model of the debris flow erosion rate:

[0060] ;

[0061] In the formula, is the dimensionless erosion rate, and is a dimensionless parameter; is the volumetric water content of the bottom bed, and is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, and is a dimensionless parameter; is the debris flow Froude number, and is a dimensionless parameter; is the debris flow Savage number, and is a dimensionless parameter; is the debris flow Bagnold number, and is a dimensionless parameter; is the volume fraction of solid phase substances in the debris flow fluid, a dimensionless parameter; is the density of the debris flow, dimension: kg / m 3 ; is the density of the particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; e = 2.718 is the base of the natural logarithm.

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

[0063] Furthermore, through the dimensionless erosion rate the erosion rate is calculated and obtained , and the calculation formula is as follows:

[0064] ;

[0065] In the formula, E is the erosion rate, dimension: m / s; is the dimensionless erosion rate, and is a dimensionless parameter; g = 9.81 is the acceleration of gravity, dimension: m / s 2 ; is the channel slope, dimension: º; is the corresponding debris flow depth, dimension: m.

[0066] In the second aspect, an embodiment of the present application provides a calculation system for the erosion rate during the movement of debris flow, including:

[0067] A dimensionless analysis module, configured to obtain debris flow fluid characteristic parameters, perform dimensionless analysis according to the debris flow fluid characteristic parameters, and calculate dimensionless parameters such as the Froude number, Savage number, and Bagnold number;

[0068] A calculation model module for constructing a dimensionless erosion rate, configured to obtain debris flow through-channel characteristic parameters, and construct a calculation model of the dimensionless erosion rate of debris flow in the channel according to the relationship between the debris flow fluid characteristic parameters and the debris flow through-channel characteristic parameters and the erosion rate;

[0069] A calculation model module for solving the dimensionless erosion rate, configured to solve the calculation model of the dimensionless erosion rate to obtain the erosion rate.

[0070] In the dimensionless analysis module:

[0071] In some embodiments, the debris flow fluid characteristic parameters include the debris flow velocity , the debris flow depth and the debris flow density .

[0072] In some embodiments, the debris flow fluid characteristic parameters are obtained through field investigations.

[0073] In some embodiments, the calculation formulas for the dimensionless parameters of debris flow, namely the Froude number, the Savage number, and the Bagnold number, are as follows:

[0074] ;

[0075] In the formula, is the Froude number, a dimensionless parameter; is the debris flow velocity (dimension: m / s); is the debris flow depth (dimension: m); , is the acceleration due to gravity (dimension: m / s 2 );

[0076] ;

[0077] In the formula, is the Savage number, a dimensionless parameter; , is the density of the particles, dimension: kg / m 3 ; , is the density of water, dimension: kg / m 3 ; is the diameter of the particles, dimension: m; is the average shear rate of the 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 debris flow depth, dimension: m; is the channel slope, dimension: º; = 30 is the friction angle between particles, dimension: º;

[0078] ;

[0079] In the formula, is the Bagnold number, a dimensionless parameter; is the density of the particles, dimension: kg / m 3 ;

[0080] is the average shear rate of the debris flow, dimension: s -1 ; is the debris flow velocity, dimension: m / s; is the debris flow depth, dimension: m; is the diameter of the particles, dimension: m; is the volume fraction of solid phase substances in the debris flow, a dimensionless parameter; is the debris flow density, dimension: kg / m 3 ; is the density of particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; is the dynamic viscosity coefficient of the debris flow liquid-phase slurry, dimension: Pa·s, the value for dilute debris flow is 0.03 Pa·s, the value for transitional debris flow is 0.05 Pa·s, and the value for viscous debris flow is 0.1 Pa·s.

[0081] In the module for constructing the calculation model of the dimensionless erosion rate:

[0082] In some embodiments, the debris flow channel characteristic parameters include the bottom bed moisture content and the content of fine particles in the bottom bed .

[0083] In some embodiments, the calculation model of the dimensionless erosion rate of the debris flow in the channel is:

[0084] ;

[0085] In the formula, is the dimensionless erosion rate, which is a dimensionless parameter; is the bottom bed volume moisture content, 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 debris flow Bagnold number, which is a dimensionless parameter; is the volume fraction of solid-phase substances in the debris flow fluid, a dimensionless parameter; is the debris flow density, dimension: kg / m 3 ; is the density of particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0086] In the module for solving the calculation model of the dimensionless erosion rate:

[0087] In some embodiments, take the logarithm of both sides of the equation of the calculation model of the dimensionless erosion rate, and we get:

[0088] ;

[0089] In the formula, is the dimensionless erosion rate, and is a dimensionless parameter; is the volumetric water content of the bottom bed, and is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, and is a dimensionless parameter; is the debris flow Froude number, and is a dimensionless parameter; is the debris flow Savage number, and is a dimensionless parameter; is the debris flow Bagnold number, and is a dimensionless parameter; is the volume fraction of solid phase substances in the debris flow fluid, a dimensionless parameter; is the density of the debris flow, dimension: kg / m 3 ; is the density of the particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; ; ; ; ; ; ; .

[0090] Based on the properties of logarithms, take the logarithm of both sides of the equation , to obtain the calculation model of the debris flow erosion rate:

[0091] ;

[0092] In the formula, is the dimensionless erosion rate, and is a dimensionless parameter; is the volumetric water content of the bottom bed, and is a dimensionless parameter; is the volume percentage of fine particles in the bottom bed, and is a dimensionless parameter; is the debris flow Froude number, and is a dimensionless parameter; is the debris flow Savage number, and is a dimensionless parameter; is the debris flow Bagnold number, and is a dimensionless parameter; is the volume fraction of solid phase substances in the debris flow fluid, a dimensionless parameter; is the density of the debris flow, dimension: kg / m 3 ; is the density of the particles, dimension: kg / m 3 ; is the density of water, dimension: kg / m 3 ; e = 2.718 is the base of the natural logarithm.

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

[0094] Furthermore, through the dimensionless erosion rate the erosion rate is calculated to obtain , and the calculation formula is as follows:

[0095] ;

[0096] In the formula, E is the erosion rate, dimension: m / s; is the dimensionless erosion rate, and is a dimensionless parameter; g = 9.81 is the acceleration due to gravity, dimension: m / s 2 ; is the channel slope, dimension: º; is the corresponding debris flow depth, dimension: m.

[0097] 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 method for calculating the erosion rate during the movement of debris flow as described above.

[0098] 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 interact with the user through a keyboard, a mouse, a remote control, a touchpad or a voice control device and other means.

[0099] 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 calculation method of the erosion rate during the debris flow movement process. In addition, the memory can also be used to temporarily store various data that have been output or will be output.

[0100] In some embodiments, the processor may be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chips. 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 calculation method of the erosion rate during the debris flow movement process.

[0101] 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 calculation method of the erosion rate during the debris flow movement process as described above.

[0102] Among them, 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 calculation method of the erosion rate during the debris flow movement process as described above.

[0103] 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 to enable a terminal device (which can be a mobile phone, computer, server, or network device, etc.) to execute the calculation method of the erosion rate during the debris flow movement process described in the embodiments of the present application.

[0104] The invention will be further described below in combination with application cases.

[0105] Application Example 1

[0106] After the "5·12" Wenchuan earthquake, the landslide geological disasters in the Zoumalinggou area increased, exacerbating the disaster risk in this area. The accumulation amount of loose materials increased, and when encountering inducing 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.

[0107] This application example is based on a detailed investigation and research of this extremely large debris flow disaster in Zoumaling, and verifies the debris flow erosion rate calculation formula of the present application. The specific steps are as follows:

[0108] The first step is to actually sample and measure that the density of the debris flow in the gully is 2240 m 3 / kg, and it is known from on-site data that the debris flow velocity is 3.45 m / s; through on-site observation data, it is determined that the water content of the bottom bed of the Wenjiagou gully is 15%, the content of fine particles in the bottom bed is 60%, the depth of the debris flow is 4 m, The median particle size of the debris flow particles is 75 mm, the particle density is 2680 kg / m, and the gully slope is 15º, The dynamic viscosity coefficient of the liquid-phase slurry is taken as 0.1 Pa·s, =30 is the inter-particle friction angle. The calculation methods of the dimensionless parameters Froude number, Savage number, Bagnold number, and the volume fraction of solid-phase substances in the debris flow are as follows:

[0109] ;

[0110] ;

[0111] ;

[0112] ;

[0113] The calculated shear rate of the debris flow is 1.29 s -1 , the solid-phase volume fraction is 0.74, the dimensionless parameter Froude number is 0.55, the Savage number is 0.00017, and the Bagnold number is 549.63;

[0114] Second, calculate the dimensionless erosion rate of the debris flow through the following formula:

[0115] ;

[0116] Obtain the dimensionless erosion rate of the debris flow = 0.0106;

[0117] Third, calculate the erosion rate of the debris flow from the dimensionless erosion rate of the debris flow. The calculation formula is as follows:

[0118] ;

[0119] Calculate and obtain the erosion rate = 0.0538 m / s.

[0120] Application Example 2

[0121] A debris flow event occurred in a mountainous area. The debris flow rushed from the gully into the river channel. Now, the debris flow river-blocking determination method of the present invention is used to judge this event. The specific steps are as follows:

[0122] First, through actual sampling, the measured density of the debris flow in the gully is 1850 m 3 / kg. Through on-site data, the flow velocity of the debris flow is known to be 5 m / s. Through on-site observation data, it is determined that the water content of the bottom bed of Wenjiagou gully is 10%, the fine particle content of the bottom bed is 50%, the flow depth of the debris flow 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 gully slope is 15º. is the dynamic viscosity coefficient of the liquid-phase slurry, and the value is 0.05 Pa·s. = 30 is the inter-particle friction angle. The calculation methods of the dimensionless parameters Froude number, Savage number, Bagnold number, and the volume fraction of solid-phase substances in the debris flow are as follows:

[0123] ;

[0124] ;

[0125] ;

[0126] ;

[0127] The calculated shear rate of the debris flow is 2.34 s -1 , the solid-phase volume fraction of the debris flow is 0.51, the dimensionless parameter Froude number is 0.894, the Savage number is 0.000559, and the Bagnold number is 577.52.

[0128] Second, calculate the dimensionless erosion rate of the debris flow through the following formula:

[0129] ;

[0130] Obtain the dimensionless erosion rate of the debris flow = 0.00181;

[0131] Third, calculate the erosion rate of the debris flow from the dimensionless erosion rate of the debris flow. The calculation formula is as follows:

[0132] ;

[0133] Calculate and obtain the erosion rate = 0.102 m / s.

[0134] Application Example 3

[0135] A debris flow event occurred in a mountainous area. The debris flow rushed from the gully into the river channel. Now, the calculation method of the erosion rate during the movement of the debris flow of the present invention is used to calculate this event. The specific steps are as follows:

[0136] First, through actual sampling, the measured density of the debris flow in the gully is 1550 m 3 / kg. From on-site data, the flow velocity of the debris flow is known to be 6.2 m / s. Through on-site observation data, it is determined that the water content of the bottom bed of Wenjiagou gully is 5%, the fine particle content of the bottom bed is 40%, and the flow depth of the debris flow is 5.6 m, the median particle size of the debris flow particles is 62 mm, the particle density is 2680 kg / m, the gully slope is 15º, is the dynamic viscosity coefficient of the liquid-phase slurry, with a value of 0.03 Pa·s, =30 is the inter-particle friction angle. The calculation methods of the dimensionless parameters Froude number, Savage number, Bagnold number, and the volume fraction of solid-phase substances in the debris flow are as follows:

[0137] ;

[0138] ;

[0139] ;

[0140] ;

[0141] The calculated debris flow shear rate is 1.66 s -1 , the solid volume fraction of the debris flow is 0.327, the dimensionless parameter Froude number is 0.838, the Savage number is 0.00014, and the Bagnold number is 277.5.

[0142] Second, calculate the dimensionless erosion rate of the debris flow through the following formula:

[0143] ;

[0144] Obtain the dimensionless erosion rate of the debris flow = 0.014;

[0145] Third, calculate the erosion rate of the debris flow from the dimensionless erosion rate of the debris flow. The calculation formula is as follows:

[0146] ;

[0147] The calculated erosion rate = 0.059 m / s.

[0148] 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. 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 to obtain the erosion rate; 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; ; ; ; ; ; ; ; 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.

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. 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; A calculation model module for solving dimensionless erosion rate is used to solve the calculation model of dimensionless erosion rate and obtain the erosion rate; 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; ; ; ; ; ; ; ; 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.

6. Computer equipment, 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 4.

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

Images