Calculation method for sorting erosion of barrier dam burst and related equipment

Abstract

The invention provides a calculation method for sorting erosion of barrier dam outburst and related equipment. The calculation method comprises the following steps: determining breach bottom erosion critical shear stress and side slope erosion critical shear stress of an erosion surface; according to the shear stress, the breach bottom erosion critical shear stress and the side slope erosion critical shear stress, the thickness of an erosion layer is determined; determining the particle size distribution of the erosion layer; determining the proportion of particles in each particle size group in the coarsening layer; determining the thickness of the coarsening layer according to the thickness of the erosion layer and the ratio of the thickness of the coarsening layer to the thickness of the active layer; determining the thickness of the updated active layer according to the thickness of the erosion layer and the thickness of the coarsening layer; and determining the particle size distribution of the updated active layer according to the proportion of each particle size group of particles in the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer and the thickness of the updated active layer. The evolutionary process of the coarsening layer in the dam burst process can be effectively simulated, and the accuracy of simulating the overtopping burst key parameters of the dam is improved.

Description

Technical Field

[0001] This application relates to the technical field of geological disaster prevention and control, and particularly relates to a calculation method and related equipment for selective erosion of dam-break of barrier dams. Background Art

[0002] Barrier lakes are caused by the accumulation of materials to form barrier dams blocking river channels, and usually occur in high mountain and canyon areas. As the water level of the barrier lake rises, the barrier dam may be damaged within a short period of time, releasing potentially destructive floods, thus causing catastrophic losses. Predicting and simulating the dam-break of barrier dams is crucial. Summary of the Invention

[0003] In view of this, the purpose of this application is to propose a calculation method and related equipment for selective erosion of dam-break of barrier dams.

[0004] Based on the above purpose, this application provides a calculation method for selective erosion of dam-break of barrier dams, including:

[0005] Determine the critical shear stress of bottom erosion and the critical shear stress of side slope erosion of the breach of the erosion surface;

[0006] Determine the thickness of the erosion layer according to the shear stress, the critical shear stress of bottom erosion of the breach, and the critical shear stress of side slope erosion;

[0007] Determine the particle size distribution of the erosion layer according to the spatial average rolling probability of soil particles in the erosion layer;

[0008] Determine the proportion of particles in each particle size group in the coarsening layer;

[0009] Determine the thickness of the coarsening layer according to the thickness of the erosion layer and the ratio of the thickness of the coarsening layer to the thickness of the active layer;

[0010] Determine the thickness of the updated active layer according to the thickness of the erosion layer and the thickness of the coarsening layer;

[0011] Determine the particle size distribution of the updated active layer according to the proportion of particles in each particle size group in the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer, and the thickness of the updated active layer.

[0012] In some of these embodiments, the determination of the proportion of particles in each particle size group in the coarsening layer includes:

[0013] By formula Determine the proportion of particles in each particle size group in the coarsening layer; where P di c is the proportion of particles with a particle size range of d i-1 to d i in the coarsening layer; d max is the maximum particle size of the active layer; dmin is the minimum particle size of the active layer; is the spatial average rolling probability of soil particles with a diameter of d i ; e is the base of the natural logarithm; x is the integration variable; is the product of the reciprocal of the Shields parameter and the shielding-exposure coefficient; p w is the density of water; G s is the particle specific gravity; τ is the shear stress of the water flow; d 50 is the median grain size; is the proportion of particles in the active layer.

[0014] In some embodiments, determining the thickness of the roughened layer according to the thickness of the erosion layer and the ratio of the thickness of the roughened layer to the thickness of the active layer includes:

[0015] Determine the thickness of the roughened layer by the formula ; wherein, T c is the thickness of the roughened layer; T cmin is the minimum thickness of the roughened layer; P c is the ratio of the thickness of the roughened layer to the thickness of the active layer; d c is the size of the smallest immovable particle in the dam body material of the barrier dam; d max is the maximum particle size of the active layer; d min is the minimum particle size of the active layer; is the 0 rolling probability; is the proportion of particles in the active layer; dH is the thickness of the erosion layer; Re* is the sand grain Reynolds number; Re* = d i ·u * / υ, u * =(t / r w ) 0.5 , υ is the kinematic viscosity; r w is the density of water; d 50 is the median grain size; G s is the particle specific gravity; τ is the shear stress of the water flow; ρ is the density of the overall soil mass.

[0016] In some embodiments, determining the thickness of the erosion layer according to the shear stress, the critical shear stress for erosion at the bottom of the breach, and the critical shear stress for erosion of the side slope includes:

[0017] Determine the thickness of the erosion layer by the formula dH = E·dt; wherein, dH is the thickness of the erosion layer; dt is the time step; E = k d (τ - τ cb ) or E = k d (τ - τcs ); E is the erosion rate per unit area; k d is the erodibility coefficient; and or and d 50 is the median grain size; d 10 is the grain size at which 10% of the particles are finer; d 60 is the grain size at which 60% of the particles are finer; e s is the void ratio; is the angle of internal friction; n is the Manning coefficient; A n is a constant; r w is the density of water; Q is the flow rate in the breach; A b is the cross-sectional area of flow; R is the hydraulic radius; θ is the slope of the bottom of the downstream section of the breach; α is the side slope angle of the breach; ps is the particle density of the soil; PI is the plasticity index; P f is the fine-grained content.

[0018] In some embodiments, determining the particle size distribution of the erosion layer according to the spatial average rolling probability of soil particles in the erosion layer includes:

[0019] Determining the particle size distribution of the erosion layer through the formula ; where P di e is the proportion of particles in the erosion layer with a particle size range from d i-1 to d i .

[0020] In some embodiments, determining the thickness of the updated active layer according to the thickness of the erosion layer and the thickness of the coarsening layer includes:

[0021] Determining the sum of the thickness of the erosion layer and the thickness of the coarsening layer to obtain the thickness of the updated active layer.

[0022] In some embodiments, determining the particle size distribution of the updated active layer according to the proportion of particles in each particle size group in the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer, and the thickness of the updated active layer includes:

[0023] Determining the particle size distribution of the updated active layer through the formula ; where FP di a is the proportion of particles with a particle size less than d i ; P di o is the proportion of particles in the original material layer with a particle size range from di-1 to d i the particle ratio; T c is the thickness of the roughened layer; P di c is the particle ratio in the roughened layer with a particle size range from d i-1 to d i ; T is the thickness of the updated active layer.

[0024] An embodiment of the present application further provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the method described in any one of the previous items is implemented.

[0025] An embodiment of the present application further provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the method described in any one of the previous items.

[0026] An embodiment of the present application further provides a computer program product including computer program instructions that, when run on a computer, cause the computer to execute the method described in any one of the previous items.

[0027] As can be seen from the above, the calculation method of sorting erosion during the breach of a barrier dam provided by the present application determines the critical shear stress for bottom erosion of the breach and the critical shear stress for side slope erosion of the erosion surface; determines the thickness of the erosion layer based on the shear stress, the critical shear stress for bottom erosion of the breach, and the critical shear stress for side slope erosion; determines the particle size distribution of the erosion layer based on the spatial average rolling probability of soil particles in the erosion layer; determines the proportion of particles in each particle size group within the roughened layer; determines the thickness of the roughened layer based on the thickness of the erosion layer and the ratio of the thicknesses of the roughened layer and the active layer; determines the thickness of the updated active layer based on the thicknesses of the erosion layer and the roughened layer; and determines the particle size distribution of the updated active layer based on the proportion of particles in each particle size group within the roughened layer, the particle size distribution of the original active layer, the thickness of the roughened layer, and the thickness of the updated active layer. It can effectively simulate the evolution process of the roughened layer during the breach of a barrier dam and improve the accuracy of key parameters for simulating the overtopping breach of a barrier dam. Description of the Drawings

[0028] To more clearly illustrate the technical solutions in the present application or related technologies, the following will briefly introduce the drawings required for use in the description of the embodiments or related technologies. Obviously, the drawings in the following description are only embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0029] Figure 1Flow chart of the calculation method for sorting erosion of barrier dam breach in the embodiments of the present application;

[0030] Figure 2 Schematic diagram of particle size - selective erosion in the embodiments of the present application;

[0031] Figure 3 Schematic diagram of the electronic device in the embodiments of the present application. Detailed implementation manners

[0032] To make the objectives, technical solutions and advantages of the present application clearer and more understandable, the present application will be further described in detail below with reference to specific embodiments and the accompanying drawings.

[0033] It should be noted that unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present application should have the ordinary meanings understood by those with ordinary skills in the field to which the present application belongs. The words such as "including" or "comprising" used in the embodiments of the present application mean that the elements or objects appearing before this word cover the elements or objects listed after this word and their equivalents, without excluding other elements or objects. The words such as "connected" or "coupled" are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect.

[0034] Barrier dams are usually composed of materials of various sizes, with particle sizes ranging from clay to gravel, and even including giant gravels with diameters sometimes reaching several meters. In the case of wide - graded dam materials, particle size - selective erosion may occur and a coarsening layer will gradually form, which in turn restricts erosion. Especially when the water flow intensity is higher than the critical intensity of movable particles, smaller particles will be continuously eroded and transported, while larger and heavier particles cannot be carried by the water flow and remain in place, eventually covering the erosion surface to form a coarsening layer and preventing further scouring. This kind of coarsening phenomenon has also been observed in the barrier dam model tests.

[0035] As Figure 1 shown, a calculation method for sorting erosion of barrier dam breach provided by an embodiment of the present invention may include:

[0036] S100, determining the critical shear stress of bottom erosion of the breach of the erosion surface and the critical shear stress of side - slope erosion;

[0037] S200, determining the thickness of the erosion layer according to the shear stress, the critical shear stress of bottom erosion of the breach and the critical shear stress of side - slope erosion;

[0038] S300, determining the particle size distribution of the erosion layer according to the spatial average rolling probability of soil particles in the erosion layer;

[0039] S400. Determine the proportion of particles in each particle size group within the coarsening layer;

[0040] S500. Determine the thickness of the coarsening layer based on the thickness of the erosion layer and the ratio of the thickness of the coarsening layer to the thickness of the active layer;

[0041] S600. Determine the thickness of the updated active layer based on the thickness of the erosion layer and the thickness of the coarsening layer;

[0042] S700. Determine the particle size distribution of the updated active layer based on the proportion of particles in each particle size group within the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer, and the thickness of the updated active layer.

[0043] The calculation method for sorting erosion during the breach of a barrier dam provided by the embodiments of the present application may include: determining the critical shear stress for erosion at the bottom of the breach and the critical shear stress for side slope erosion of the erosion surface; determining the thickness of the erosion layer based on the shear stress, the critical shear stress for erosion at the bottom of the breach, and the critical shear stress for side slope erosion; determining the particle size distribution of the erosion layer based on the spatial average rolling probability of soil particles in the erosion layer; determining the proportion of particles in each particle size group within the coarsening layer; determining the thickness of the coarsening layer based on the thickness of the erosion layer and the ratio of the thickness of the coarsening layer to the thickness of the active layer; determining the thickness of the updated active layer based on the thickness of the erosion layer and the thickness of the coarsening layer; and determining the particle size distribution of the updated active layer based on the proportion of particles in each particle size group within the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer, and the thickness of the updated active layer; which can effectively simulate the evolution process of the coarsening layer during the breach of the barrier dam and improve the accuracy of key parameters for simulating the overtopping breach of the barrier dam.

[0044] See Figure 2 As shown, the soil mass / dam in contact with the water flow is divided into three layers: the erosion layer, the coarsening layer, and the original layer. Among them, the erosion layer is usually composed of particles washed away by the water flow, while the coarsening layer is composed of particles remaining in the active layer. The erosion layer and the coarsening layer are collectively referred to as the active layer, that is, the area affected by the water flow. And the original layer is the remaining part under the active layer. After erosion is completed, the thickness and particle size distribution (GSD) of the new active layer (i.e., the updated active layer) may change and be uniformly mixed internally.

[0045] In some embodiments, in step S100, determining the critical shear stress for erosion at the bottom of the breach and the critical shear stress for side slope erosion of the erosion surface can be obtained according to the modified shear stress equation. Specifically, it can be obtained by adjusting the shear stress equation and the critical shear stress equation under the condition of having a slope.

[0046] In some of these embodiments, the shear stress equation for the shear stress τ of the water flow can be in the form of the water depth and flow rate, as shown in Equation (1). The shear stress in Equation (1) can actually be understood as being based on the average flow velocity in open-channel flow, that is, the average shear stress. Generally, the shear stress acting on the side walls of the breach is lower than the shear stress acting on the bottom of the breach.

[0047] In some of these embodiments, the critical shear stress equation can be as shown in Equation (2) or Equation (4). Among them, Equation (2) is the critical shear stress equation in the actual scenario, and the corresponding soil erodibility coefficient can be as shown in Equation (3). And Equation (4) is the critical shear stress equation in the laboratory scenario, and the corresponding soil erodibility coefficient can be as shown in Equation (5). During specific detection, it can be selected according to requirements. Among them, where τ c is the critical shear stress. k d is the soil erodibility coefficient. e s is the void ratio. d 10 and d 60 respectively represent the particle sizes (mm) with 10% and 60% finer particles. ρ s is the particle density of the soil. ρ w is the particle density of the water. is the angle of internal friction (°). PI is the plastic index. P f is the fine-grained content (<0.063 mm). d 50 is the median particle size (mm).

[0048] In some of these embodiments, for the critical shear stress at the bottom of the breach on the erosion surface, it can be as shown in Equation (6). The critical shear stress for side slope erosion can be as shown in Equation (7). In the erosion of the barrier dam, the erosion surface is usually a slope. In this way, the change in the force acting on the soil particles by the slope can be better considered, and its influence on the critical shear stress can be considered.

[0049] In some of these embodiments, in step S200, determining the thickness of the erosion layer according to the shear stress, the critical shear stress at the bottom of the breach, and the critical shear stress for side slope erosion may include:

[0050] Determining the thickness of the erosion layer through Equation (8). dH = E·dt, Equation (8). Among them, dH is the thickness of the erosion layer. dt is the time step; E = k d (τ - τ cb ) is the erosion rate per unit area at the bottom of the breach or E = k d (τ - τ cs) is the erosion rate per unit area of the side slope. It should be understood that, based on the erosion rates in different directions, the thicknesses of the erosion layers in different erosion directions can be calculated, that is, the thickness of the erosion layer at the bottom of the breach and the erosion thickness of the side slope can be obtained. k d is the erodibility coefficient.

[0051] and or and d 50 is the median grain size. d 10 is the grain size at which 10% of the particles are finer. d 60 is the grain size at which 60% of the particles are finer. e s is the void ratio. is the angle of internal friction. τ is the average shear stress, as in the aforementioned formula (1). n is the Manning coefficient. A n is a constant. ρ w is the density of water. Q is the flow rate in the breach; A b is the cross-sectional area of flow. R is the hydraulic radius. θ is the slope at the bottom of the downstream section of the breach; α is the angle of the side slope of the breach; ps is the particle density of the soil; PI is the plasticity index; P f is the fine-grained content. It should be understood that the particle density of the soil does not include voids.

[0052] In some embodiments, in step S300, the grain size distribution (GSD) of the erosion layer can be determined by the erosion ratio of each grain size group, usually through moving probability calculation. Among them, the critical condition for particle rolling requires the lowest critical shear force compared with sliding or jumping. Thus, rolling is the most likely movement to occur. When the external torque is greater than the torque to keep the particle stationary, it can be considered that a single particle rolls. However, for wide-graded dam materials, there are at least two uncertainties in particle rolling: (1) For any particle on the dam surface, its size and position relative to other particles are random. Specifically, coarse particles are exposed to the water flow, while fine particles are usually affected by the shielding effect of coarse particles, and both of these factors will affect the particle initiation speed / stress. (2) The flow at the bottom of the channel is pulsating, resulting in random forces on the particles at specific positions. Therefore, the rolling probability of particles in each grain size group is regarded as the proportion of eroded particles in each grain size group.

[0053] In some embodiments, the determining the grain size distribution of the erosion layer according to the spatial average rolling probability of the soil particles in the erosion layer includes: determining the grain size distribution of the erosion layer through formula (9). where P di e is the grain size range in the erosion layer as di-1 to d i The particle ratio. d max is the maximum particle size of the active layer. d min is the minimum particle size of the active layer. When d i < d cmax , it is assumed that the instantaneous underflow velocity follows a normal distribution and the shelter-exposure effect is considered. The spatial average rolling probability e of soil particles with a diameter of d i That is, di . That is, is the spatial average rolling probability of soil particles with a diameter of 0 and d i as shown in Equation (10). Usually, when the particle size in the soil layer exceeds the critical diameter d of immovable particles cmax (d i > d cmax ), it will not move in the water flow, indicating that its rolling probability is zero When the Shields parameter is set to 0.023, d cmax。 d cmax can be as shown in Equation (11). is the particle ratio in the active layer.

[0054] Among them, Among them, e is the base of the natural logarithm. x is the integration variable. G s is the particle specific gravity. is the product of the reciprocal of the Shields parameter and the shelter-exposure coefficient, as shown in Equation (12). Among them, p w is the density of water. τ is the shear stress of the water flow, as shown in Equation (1). d 50 is the median grain size.

[0055] In some of these embodiments, in step S400, in each erosion event, the un-eroded soil mass usually integrates into the coarsening layer. Usually, the ratio of the thickness of the coarsening layer to the thickness of the active layer can be set to be proportional to the ratio of the remaining particle mass to the total particle mass.

[0056] In some of these embodiments, the determination of the proportion of particles in each particle size group in the coarsening layer may include: determining the proportion of particles in each particle size group in the coarsening layer through Equation (13). Among them, P di c is the particle ratio in the coarsening layer with a particle size range of d i-1 to d i . d max is the maximum particle size of the active layer. d min is the minimum particle size of the active layer. is the diameter of di Spatial average rolling probability of soil particles. e is the base of the natural logarithm. x is the integration variable. is the product of the reciprocal of the Shields parameter and the shielding-exposure coefficient, as shown in Equation (12).

[0057] In some of these embodiments, in step S500, the weighted average method can be used to determine the ratio of the thickness of the erosion layer to the thickness of the roughening layer. Among them, the particle size can be used as the weighting factor. Determining the thickness of the roughening layer according to the thickness of the erosion layer and the ratio of the thickness of the roughening layer to the thickness of the active layer can include: determining the ratio of the thickness of the roughening layer to the thickness of the active layer through Equation (14). Determining the thickness of the roughening layer through Equation (15). That is, through Determining the thickness of the roughening layer. Where T c is the thickness of the roughening layer; T cmin is the minimum thickness of the roughening layer, usually related to the size d c of the smallest immovable particles in the dam body material, as shown in Equation (16). P c is the ratio of the thickness of the roughening layer to the thickness of the active layer, as shown in Equation (14). dH is the thickness of the erosion layer, as described in Equation (8).

[0058]

[0059] d c is the size of the smallest immovable particles in the dam body material of the barrier dam, usually determined by the sand grain Reynolds number. d max is the maximum particle size of the active layer. d min is the minimum particle size of the active layer. is the rolling probability. is the particle proportion in the active layer; dH is the thickness of the erosion layer. Re* is the sand grain Reynolds number. Re* = d i ·u * / υ, u * =(t / r w ) 0.5 , υ is the kinematic viscosity; r w is the density of water; d 50 is the median particle size; G s is the particle specific gravity; τ is the shear stress of the water flow; ρ is the density of the overall soil mass. It should be understood that the density of water is generally a fixed value of 1000 kg / m 3 . Since the thickness of the erosion layer includes the thickness of the erosion layer at the bottom of the breach and the erosion thickness of the side slope in these two directions, the thickness of the roughening layer can also correspondingly include the thickness of the roughening layer at the bottom of the breach and the thickness of the roughening layer of the side slope in these two directions.

[0060] In some of these embodiments, the maximum thickness of the roughened layer can be as shown in Equation (17). Specifically, where T cmax is the maximum thickness of the roughened layer. e s is the soil porosity, and P dc is the proportion of armor particles.

[0061] In some of these embodiments, in step S600, the roughened layer and part of the original layer form a new active layer. Determining the thickness of the updated active layer according to the thickness of the eroded layer and the thickness of the roughened layer may include: determining the sum of the thickness of the eroded layer and the thickness of the roughened layer to obtain the thickness of the updated active layer. Specifically, the thickness of the updated active layer can be calculated by the formula T = dH b +T c It should be understood that since the thickness of the eroded layer includes the thickness of the eroded layer at the bottom of the breach and the erosion thickness of the side slope, and the thickness of the roughened layer also includes the thickness of the roughened layer at the bottom of the breach and the thickness of the roughened layer of the side slope in these two directions, the thickness of the updated active layer can include the thickness of the updated active layer at the bottom of the breach and the updated thickness of the side slope in these two directions.

[0062] In some of these embodiments, in step S700, the roughened layer and the original layer in the updated active layer are fully mixed to form a uniform updated active layer. Determining the particle size distribution of the updated active layer according to the proportion of particles in each particle size group in the roughened layer, the particle size distribution of the original active layer, the thickness of the roughened layer, and the thickness of the updated active layer may include:

[0063] By the formula to determine the particle size distribution of the updated active layer. Where FP di a is the proportion of particles with a particle size less than d i ; P di o is the proportion of particles in the original material layer with a particle size range from d i-1 to d i ; T c is the thickness of the roughened layer; P di c is the proportion of particles in the roughened layer with a particle size range from d i-1 to d i ; and T is the thickness of the updated active layer.

[0064] It should be noted that the method of the embodiment of the present application can be executed by a single device, such as a computer or a server. The method of this embodiment can also be applied to a distributed scenario and completed by multiple devices cooperating with each other. In such a distributed scenario, one of the multiple devices can only execute one or more steps of the method of the embodiment of the present application, and these multiple devices will interact with each other to complete the described method.

[0065] It should be noted that some embodiments of the present application have been described above. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be performed in a different order than in the above embodiments and still achieve the desired results. Additionally, the processes depicted in the drawings do not necessarily require the particular order or sequential order shown to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0066] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the calculation method of the selective erosion of the landslide dam breach described in any one of the above embodiments.

[0067] Figure 3 A more specific schematic diagram of the hardware structure of the electronic device provided in this embodiment is shown. The device may include: a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. Among them, the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040 are communicatively connected to each other inside the device through the bus 1050.

[0068] The processor 1010 can be implemented in ways such as a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0069] The memory 1020 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage devices, dynamic storage devices, etc. The memory 1020 can store the operating system and other application programs. When implementing the technical solutions provided in the embodiments of this specification through software or firmware, the relevant program codes are stored in the memory 1020 and are called and executed by the processor 1010.

[0070] The input / output interface 1030 is used to connect to the input / output module to achieve information input and output. The input / output module can be configured as a component in the device (not shown in the figure) or can be externally connected to the device to provide corresponding functions. Among them, the input devices can include keyboards, mice, touchscreens, microphones, various sensors, etc., and the output devices can include displays, speakers, vibrators, indicator lights, etc.

[0071] The communication interface 1040 is used to connect to the communication module (not shown in the figure) to achieve communication interaction between this device and other devices. Among them, the communication module can achieve communication through wired means (such as USB, network cable, etc.) or can also achieve communication through wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0072] The bus 1050 includes a path for transmitting information between the various components of the device (such as the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040).

[0073] It should be noted that although the above device only shows the processor 1010, the memory 1020, the input / output interface 1030, the communication interface 1040, and the bus 1050, in the specific implementation process, the device may also include other components necessary for normal operation. In addition, those skilled in the art can understand that the above device may also only include the components necessary to implement the solutions of the embodiments of this specification and does not necessarily include all the components shown in the figure.

[0074] The electronic device in the above embodiment is used to implement the calculation method for the selective erosion of the dam-break of the barrier dam corresponding to any one of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0075] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to execute the calculation method for the selective erosion of the dam-break of the barrier dam as described in any of the above embodiments.

[0076] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette tapes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0077] The computer instructions stored in the storage medium of the above embodiment are used to cause the computer to execute the calculation method of the sorting erosion of the dam-break of the barrier dam described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0078] Based on the same inventive concept, corresponding to the calculation method of the sorting erosion of the dam-break of the barrier dam described in any of the above embodiments, the present disclosure also provides a computer program product, which includes computer program instructions. In some embodiments, the computer program instructions can be executed by one or more processors of the computer to cause the computer and / or the processor to execute the calculation method of the sorting erosion of the dam-break of the barrier dam. Corresponding to the execution subject of each step in each embodiment of the calculation method of the sorting erosion of the dam-break of the barrier dam, the processor executing the corresponding step can belong to the corresponding execution subject.

[0079] The computer program product of the above embodiment is used to cause the computer and / or the processor to execute the calculation method of the sorting erosion of the dam-break of the barrier dam described in any of the above embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0080] Those of ordinary skill in the art should understand that: the discussion of any of the above embodiments is only exemplary, and is not intended to imply that the scope of the present application (including the claims) is limited to these examples; under the idea of the present application, the technical features in the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations in different aspects of the embodiments of the present application as described above, and they are not provided in detail for the sake of brevity.

[0081] In addition, for simplicity of explanation and discussion, and in order not to make the embodiments of the present application difficult to understand, known power / ground connections to integrated circuit (IC) chips and other components may or may not be shown in the provided drawings. Further, the devices may be shown in block diagram form in order to avoid making the embodiments of the present application difficult to understand, and this also takes into account the fact that details of the implementation of these block diagram devices are highly dependent on the platform on which the embodiments of the present application are to be implemented (i.e., these details should be entirely within the understanding of those skilled in the art). In cases where specific details (e.g., circuits) are set forth to describe exemplary embodiments of the present application, it will be apparent to those skilled in the art that the embodiments of the present application may be practiced without these specific details or with variations of these specific details. Accordingly, these descriptions should be considered illustrative rather than restrictive.

[0082] Although the present application has been described in connection with specific embodiments thereof, many alternatives, modifications, and variations of these embodiments will be apparent to those of ordinary skill in the art in light of the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0083] Embodiments of the present application are intended to cover all such alternatives, modifications, and variations that fall within the broad scope of the appended claims. Accordingly, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principle of the embodiments of the present application shall be included within the protection scope of the present application.

Claims

1. A method for calculating the sorting erosion of dam breach, characterized in that: include: Determine the critical shear stress of the breach bottom and the critical shear stress of the side slope of the erosion surface; Determine the thickness of the erosion layer based on the shear stress, the critical shear stress of erosion at the bottom of the breach, and the critical shear stress of erosion on the side slope; Determining the particle size distribution of the erosion layer according to the spatial average rolling probability of soil particles in the erosion layer; Determine the proportion of particles of each size group in the coarsening layer; According to the thickness of the erosion layer and the ratio of the thickness of the coarsening layer to the thickness of the active layer, the thickness of the coarsening layer is determined; According to the thickness of the erosion layer and the thickness of the coarsening layer, the thickness of the updated active layer is determined; The particle size distribution of the updated active layer is determined based on the proportion of particles of each particle size group in the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer and the thickness of the updated active layer.

2. The calculation method of the sorting erosion of the dam breach according to claim 1 is characterized in that: Determining the proportion of particles of each particle size group in the coarsening layer includes: Pass-through Determine the proportion of particles in each size group in the coarsening layer; where P di c The particle size range in the coarsening layer is d i-1 to d i The particle ratio; d max is the maximum particle size of the active layer; d min is the minimum particle size of the active layer; = 0, diameter is d i The spatial average rolling probability of soil particles; e is the base of natural logarithm; x is the integral variable; It is the product of the inverse of the Shields parameter and the shielding-exposure coefficient; ρ w is the density of water; G s is the particle specific gravity; τ is the shear stress of the water flow; d 50 For the middle d = value d Particle (1-diameter ε; d is the proportion of particles in the active layer.

3. The calculation method of the sorting erosion of the dam breach according to claim 1 is characterized in that: Determining the thickness of the coarsening layer according to the thickness of the erosion layer and the ratio of the thickness of the coarsening layer to the thickness of the active layer comprises: Pass-through Determine the thickness of the roughened layer; Among them, T c is the thickness of the roughened layer; T cmin is the minimum thickness of the coarse layer; P c is the ratio of the thickness of the coarsening layer to the thickness of the active layer; d c The smallest immovable particle size in the dam material; d max is the maximum particle size of the active layer; d min is the minimum particle size of the active layer; =Rolling probability is 0; is the particle ratio in the active layer; dH is the thickness of the erosion layer; Re* is the sand Reynolds number; Re* = d i ·u * / υ,u * =(t / r w ) 0.5 ,υ is the kinematic viscosity; r w is the density of water; d 50 is the median particle size; G s is the particle specific gravity; τ is the shear stress of water flow; ρ is the density of the entire soil.

4. The method for calculating the sorting erosion of dam breach according to claim 1, characterized in that: Determining the thickness of the erosion layer according to the shear stress, the critical shear stress of erosion at the bottom of the breach, and the critical shear stress of erosion at the side slope comprises: The thickness of the erosion layer is determined by the formula dH=E·dt; wherein dH is the thickness of the erosion layer; dt is the time step; E=k d (τ-τ cb ) or E = k d (τ-τ cs ); E is the erosion rate per unit area; k d is the erodibility coefficient; and or and d 50 is the median particle size; d 10 There are 10% finer particles; d 60 60% of the particles are finer in size; s is the void ratio; is the internal friction angle; n is the Manning coefficient; A n is a constant; r w is the density of water; Q is the flow rate in the breach; A b is the water flow area; R is the hydraulic radius; θ is the slope of the bottom of the downstream section of the breach; α is the side slope angle of the breach; ρ s is the particle density of soil; PI is the plasticity index; P f For fine particle content.

5. The calculation method of the sorting erosion of the dam breach according to claim 2 is characterized in that: Determining the particle size distribution of the erosion layer according to the spatial average rolling probability of soil particles in the erosion layer comprises: Pass-through Determine the particle size distribution of the eroded layer; wherein P di e The particle size range in the eroded layer is d i-1 to d i The particle ratio.

6. The method for calculating the sorting erosion of dam breach according to claim 5, characterized in that: According to the thickness of the erosion layer and the thickness of the coarsening layer, the thickness of the updated active layer is determined by: The sum of the thickness of the eroded layer and the thickness of the roughened layer is determined to obtain the thickness of the updated active layer.

7. The method for calculating the sorting erosion of dam breach according to claim 5, characterized in that: Determining the particle size distribution of the updated active layer according to the proportion of particles of each particle size group in the coarsening layer, the particle size distribution of the original active layer, the thickness of the coarsening layer and the thickness of the updated active layer includes: Pass-through Determine the updated particle size distribution of the active layer; where FP di a For particles smaller than d i The particle ratio of P di o The particle size range in the original material layer is d i-1 to d i The particle ratio of T c is the thickness of the coarse layer; P di c The particle size range in the coarsening layer is d i-1 to d i The particle ratio; T is the thickness of the updated active layer.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method according to any one of claims 1 to 7 when executing the program.

9. A non-transitory computer-readable storage medium storing computer instructions, wherein the computer instructions are used to cause a computer to execute the method according to any one of claims 1 to 7.

10. A computer program product, comprising computer program instructions, which, when executed on a computer, cause the computer to execute the method according to any one of claims 1 to 7.

Images