CN116663223A - Dam break flood evolution prediction method based on wave breaking principle - Google Patents
Dam break flood evolution prediction method based on wave breaking principle Download PDFInfo
- Publication number
- CN116663223A CN116663223A CN202310087495.5A CN202310087495A CN116663223A CN 116663223 A CN116663223 A CN 116663223A CN 202310087495 A CN202310087495 A CN 202310087495A CN 116663223 A CN116663223 A CN 116663223A
- Authority
- CN
- China
- Prior art keywords
- wave
- flow
- section
- dam break
- upstream
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 230000008569 process Effects 0.000 claims abstract description 33
- 238000011144 upstream manufacturing Methods 0.000 claims description 101
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 50
- 239000011159 matrix material Substances 0.000 claims description 40
- 238000004364 calculation method Methods 0.000 claims description 32
- 238000012545 processing Methods 0.000 claims description 13
- 238000011156 evaluation Methods 0.000 claims description 7
- 238000013441 quality evaluation Methods 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 description 9
- 239000002245 particle Substances 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 239000002689 soil Substances 0.000 description 3
- 239000000758 substrate Substances 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 230000010429 evolutionary process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
The invention discloses a dam break flood evolution prediction method based on a wave breaking principle, which comprises the following steps: constructing a dam break flood evolution model; acquiring parameter information of flood evolution and flow process information generated at a breach position when dam break occurs; inputting the parameter information and the flow process information generated by the position of the breach when the dam break occurs into a dam break flood evolution model, and outputting a flood evolution result. The method can simulate the burst flood evolution in a simple and efficient manner, and provides technical support for downstream emergency decision.
Description
Technical Field
The invention relates to the technical field of flood evolution prediction, in particular to a dam break flood evolution prediction method based on a wave breaking principle.
Background
In the downstream flowing process of the burst flood, the flood peak is continuously reduced, and the flow process line is gradually flattened. The study of the dam break flood evolution model essentially reveals the change rule of hydraulic parameters such as flow, water depth, flow velocity and the like along with time and evolution distance.
At present, the simulation of Hong Shuixia game evolution process mainly adopts a numerical calculation method, and according to different theoretical basis, the numerical calculation model is divided into two types, namely a continuum model and a particle flow model. In the continuum model, the flow of flood is described by an N-S equation, and the basal resistance is expressed by a Manning equation when water flow is simulated. The particle flow model is mainly used for simulating the surface movement of the high-concentration mud sand flow, discrete bodies are scattered into particle sets, and different forms of particles are assumed to contact the model, so that the movement forms of different particles forming substances are simulated.
The two models are different in basic principle and have essential differences, and the particle flow model can better simulate the process and characteristics of rock and soil damage, but is difficult to simulate the erosion effect of flood and mud and sand flow on a substrate in the motion process, and is suitable for simulating the whole process of landslide dam-break rock and soil body breaking and evolution. The continuum model can simulate the erosion effect of a moving substance on a substrate, the physical meaning of each parameter of an equation is clear, the solving efficiency is far higher than that of a particle flow model, and the continuum model is suitable for simulating the movement characteristics of various materials by adopting different substrate friction resistance models and soil pressure coefficient expressions. However, the N-S equation is derived based on a constant gradual flow, and when a dam break occurs, the flow rate changes greatly in a short time, and the induced flow is necessarily a non-constant rapid flow. Dam break waves are typically discontinuous waves, i.e. broken waves, but no research has been done to apply the broken wave principle to breaking flood evolution.
Disclosure of Invention
The present invention aims to solve, at least to some extent, one of the technical problems in the above-described technology. Therefore, the invention aims to provide the dam break flood evolution prediction method based on the wave breaking principle, which can simulate the dam break flood evolution in a simple and efficient way and provide technical support for downstream emergency decision.
In order to achieve the above objective, an embodiment of the present invention provides a dam break flood evolution prediction method based on a wave breaking principle, including:
constructing a dam break flood evolution model;
acquiring parameter information of flood evolution and flow process information generated at a breach position when dam break occurs;
inputting the parameter information and the flow process information generated by the position of the breach when the dam break occurs into a dam break flood evolution model, and outputting a flood evolution result.
According to some embodiments of the invention, constructing a dam-break flood model includes:
constructing an initial model;
and adding an operation rule into the initial model to construct a dam-break flood evolution model.
According to some embodiments of the invention, the operation rule includes: continuity equation and wave flow, momentum equation, wave velocity of the wave break and basic equation of the wave break.
According to some embodiments of the invention, a method of determining a continuity equation and a wave flow rate includes:
setting the flow velocity of the unaffected area of the wave breaking crest as v by adopting fixed coordinates 0 Corresponding cross-sectional area A 0 The method comprises the steps of carrying out a first treatment on the surface of the The flow velocity of the area affected by the wave crest of the broken wave is v, the corresponding cross-sectional area is A, and the wave velocity of the broken wave is w; the velocity is the dynamic coordinate of the wave breaking velocity w parallel to the wave propagation direction, and the flow velocity before and after wave crest influence is v 0 -w, v-w, according to a constant flow continuity equationObtaining:
(v-w)A＝(v 0 -w)A 0
vA-v 0 A 0 ＝w(A-A 0 )
introducing a flow Q at section h0 0 And the flow rate Q at section h, can be obtained:
Q-Q 0 ＝w(A-A 0 )
let Δq=q-Q 0 Is wave flow and the area difference deltaa=A-A before and after the peak 0 Determining a relation expression between the flow rate and the wave speed and the wave height of the broken wave:
ΔQ＝wΔA＝ζB′w
in the method, in the process of the invention,B 0 the initial water surface width of the section before the occurrence of the wave breaking; b is the water surface width of the section after the occurrence of the wave breaking; ζ=h-h 0, wave height.
According to some embodiments of the invention, a method of determining a momentum equation includes:
analyzing the dynamic coordinate with the speed of the broken wave w in the broken wave unsteady flow, and taking the water body between the first section and the second section as a control body; the total dynamic water pressure distribution of the first section and the second section is respectively P and P 0 ，
Because the control body is smaller, the gravity component and the boundary friction force acting on the control body are ignored, and a constant flow momentum equation can be applied to the control body along the water flow direction, so that the following results:
wherein gamma is the weight of water, 9.8kN/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, 9.8m/s 2 。
According to some embodiments of the invention, a method of determining a breakup wave speed includes:
the wave velocity of the broken wave can be obtained according to the continuity equation and the momentum equation;
from the equation of continuity,
(v-w)A＝(v 0 -w)A 0
consider wave flow Δq=q-Q 0 And Δa=A-A 0 ，Δv＝v 0 V, the above formula can be written as:
(A 0 +ΔA)[w-(v 0 +Δv)]＝A 0 (w-v 0 )
after finishing, the method comprises the following steps:
according to the momentum equation,
wherein: Δp=p-P 0
Substituting the continuity equation into the above equation:
finishing to obtain wave velocity formula
In the above, the positive sign before the root corresponds to forward wave field, and the negative sign corresponds to backward wave field;
approximately considering that the dynamic water pressure intensity on the first section and the second section meets the distribution rule of static water pressure intensity, obtaining:
ΔP＝P-P 0 ＝γ(Ay c -A 0 y c0 )
wherein y is c Is the depth of the centroid of the area A under water, y c0 Is area A 0 The centroid of (c) is at the depth of the water; let y be c ' is the depth of the centroid of the area deltaa under water, and ζ is the height difference of the water surface before and after the peak, i.e. the wave height, according to the area moment property, for the area moment Ay c The method can obtain:
Ay c ＝(A 0 +ΔA)y c ＝A 0 y c0 +ΔAy c ′
at this time, the dynamic water pressure difference is as follows:
ΔP＝γ(A 0 y c0 +ΔAy c ′+A 0 ζ-A 0 y c0 )＝γ(ΔAy c ′+A 0 ζ)
substituting the above formula into a wave velocity formula to obtain:
assume that: Δa=ζb',
the above can be simplified as:
average water depth is introduced:obtaining:
in the case of a constant cross-sectional shape, v is 0 、h 0 、B 0 Equal known, i.e. the wave velocity is a single value of wave heightA function.
According to some embodiments of the invention, a method of determining a broken wave basis equation comprises:
dam break wave belongs to forward wave, calculates based on forward wave:
and obtaining a wave breaking basic equation according to the calculation.
According to some embodiments of the invention, the parameter information includes upstream section geometry parameters, downstream section geometry parameters, river parameters, and control calculation parameters;
the geometric parameters of the upstream section comprise the bottom elevation, the bottom width and the section slope ratio of the upstream section;
the geometric parameters of the downstream section comprise the bottom elevation, the bottom width and the section slope ratio of the downstream section;
the river parameters include: calculating time, natural runoff and Manning coefficients;
the control calculation parameters comprise an initial time interval, an evolution distance and a resistance coefficient.
According to some embodiments of the present invention, obtaining flow process information generated by a breach position when a dam break occurs includes:
dividing the flow process into n disconnected wave element waves according to a time average mode, wherein the numbers are respectively 1-n;
the number i of the broken element wave is that the average flow rate in the period is:can obtain the corresponding wave height ζ i ；
The total flow of the broken wave element wave with the number of i is the product of the average flow in the period and the time length, namely:
wherein: q i For the ith wave-breaking elementTotal flow, m 3 ；Q i For the average flow of the ith broken element wave, m 3 S; dt is the time interval, s.
According to some embodiments of the invention, obtaining the upstream cross-section geometry includes:
shooting a plurality of upstream section images based on the unmanned aerial vehicle;
performing image quality evaluation on a plurality of upstream section images, and determining a target upstream section image according to an evaluation result;
performing image recognition on the target upstream section image, and determining an upstream section geometric parameter according to a recognition result;
image quality evaluation is carried out on a plurality of upstream section images, and a target upstream section image is determined according to an evaluation result, wherein the method comprises the following steps:
performing standardization processing on the size of the upstream section image to obtain a processed upstream section image;
acquiring pixel points and pixel values on each processed upstream section image, and constructing a first pixel matrix;
calculating a difference value of the first pixel matrix and the second pixel matrix; the second pixel matrix is a matrix constructed based on pixel points and pixel values on the standard upstream section image;
wherein F is i A difference value between the second pixel matrix and the first pixel matrix constructed based on the ith processed upstream cross-sectional image;first pixel matrix A constructed for processing upstream cross-sectional image based on ith sheet i The pixel value of row t column s; c (C) s,t I=1, 2, 3 … … G for the pixel value of the s-th row t column of the second pixel matrix; g is the number of processed upstream sectional images; m is the number of pixels of the standard upstream cross-section image and the processed upstream cross-section image in the transverse direction; n is standard upstream cross-section image and processing upstream cross-section imageLike the number of pixels in the vertical direction;
and determining the processing upstream section image corresponding to the first pixel matrix with the minimum difference value as a target upstream section image.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention may be realized and attained by the structure particularly pointed out in the written description and drawings.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention. In the drawings:
fig. 1 is a flowchart of a dam break flood prediction method based on a wave breaking principle according to an embodiment of the present invention;
FIG. 2 is a schematic representation of a dam break wave according to one embodiment of the present invention;
FIG. 3 is a graph of conservation of wave mass analysis according to one embodiment of the invention;
FIG. 4 is a channel cross-sectional view according to one embodiment of the invention;
FIG. 5 is a graph of conservation of momentum analysis according to one embodiment of the invention;
FIG. 6 is a profile hydrodynamic pressure profile according to one embodiment of the invention;
FIG. 7 is a flow process division diagram according to one embodiment of the invention;
fig. 8 is a schematic diagram of dam break flood input-output according to one embodiment of the invention;
FIG. 9 is a schematic diagram of dam break flood calculation according to one embodiment of the invention;
FIG. 10 is a schematic illustration of north Sichuan station calculation parameters according to one embodiment of the invention;
FIG. 11 is a schematic illustration of North Sichuan station calculation according to one embodiment of the invention;
FIG. 12 is a schematic illustration of calculating parameters for a through-port station according to one embodiment of the invention;
FIG. 13 is a schematic diagram of the calculation result of a through-port station according to one embodiment of the present invention;
FIG. 14 is a schematic representation of calculated parameters of a Fujiang bridge in accordance with one embodiment of the present invention;
figure 15 is a schematic diagram of the calculation of the Fu Jiang Qiaoce station in accordance with one embodiment of the invention.
Detailed Description
The preferred embodiments of the present invention will be described below with reference to the accompanying drawings, it being understood that the preferred embodiments described herein are for illustration and explanation of the present invention only, and are not intended to limit the present invention.
According to the embodiment of the invention, a dam break flood evolution prediction method based on a wave breaking principle is provided according to the attached figures 1-15.
As shown in fig. 1, the embodiment of the invention provides a dam break flood evolution prediction method based on a wave breaking principle, which comprises the following steps of S1-S3:
s1, constructing a dam break flood evolution model;
s2, acquiring parameter information of flood evolution and flow process information generated by a breach position when dam break occurs;
and S3, inputting the parameter information and the flow process information generated by the position of the breach when the dam break occurs into a dam break flood evolution model, and outputting a flood evolution result.
The working principle of the technical scheme is as follows: constructing a dam break flood evolution model; acquiring parameter information of flood evolution and flow process information generated at a breach position when dam break occurs; inputting the parameter information and the flow process information generated by the position of the breach when the dam break occurs into a dam break flood evolution model, and outputting a flood evolution result.
The beneficial effects of the technical scheme are that: the method is characterized in that a wave-breaking principle is applied to break flood evolution, a flood evolution model is established based on the wave-breaking principle, an EXCEL is used as a platform, and a corresponding numerical calculation program is developed by adopting VBA language, so that break flood evolution can be simulated succinctly and efficiently, and technical support is provided for downstream emergency decision.
The basic principle of wave breaking is as follows: when the flow rate of the open cell changes significantly in a short period of time, a non-constant rapid flow is induced. The non-constant rapid flow is characterized by a stepped surge or fall front over a short distance, the instantaneous water surface gradient being very steep and its hydraulic components no longer being a continuous function of time and flow path. Known as discontinuous waves, also known as broken waves. Such as tidal bore, tidal bore generated in the channel by the change of the input flow of the hydraulic turbine of the hydropower station, and tidal bore propagated downstream after the dam break, etc. These peaks are steep and the resulting stepped leading edge shape can be maintained for a period of time as shown in FIG. 2.
According to some embodiments of the invention, constructing a dam-break flood model includes:
constructing an initial model;
and adding an operation rule into the initial model to construct a dam-break flood evolution model.
According to some embodiments of the invention, the operation rule includes: continuity equation and wave flow, momentum equation, wave velocity of the wave break and basic equation of the wave break.
According to some embodiments of the invention, a method of determining a continuity equation and a wave flow rate includes:
as shown in FIG. 3 (a), the flow velocity of the unaffected region of the broken wave crest is set to be v 0 Corresponding cross-sectional area A 0 The method comprises the steps of carrying out a first treatment on the surface of the The flow velocity of the area affected by the wave crest of the broken wave is v, the corresponding cross-sectional area is A, and the wave velocity of the broken wave is w; the velocity is the dynamic coordinate of the wave breaking velocity w parallel to the wave propagation direction, as shown in FIG. 3 (b), where the flow velocity before and after the influence of the wave peak is v 0 -w, v-w, according to a constant flow continuity equation, obtainable:
(v-w)A＝(v 0 -w)A 0
vA-v 0 A 0 ＝w(A-A 0 )
introducing a flow Q at section h0 0 And the flow rate Q at section h, can be obtained:
Q-Q 0 ＝w(A-A 0 )
let Δq=q-Q 0 Is wave flow and the area difference deltaa=A-A before and after the peak 0 Determining a relation expression between the flow rate and the wave speed and the wave height of the broken wave:
ΔQ＝wΔA＝ζB′w
in the method, in the process of the invention,B 0 the initial water surface width of the section before the occurrence of the wave breaking; b is the water surface width of the section after the occurrence of the wave breaking; ζ=h-h 0, wave height. As shown in fig. 4.
The technical scheme has the working principle and beneficial effects that: when the wave break occurs, the water flow is changed only at the wave peak of the wave break or after the wave peak passes. I.e. the water flow reflects a non-constant characteristic only in the vicinity of the peak. If flow coordinate analysis is used, the non-constant problem can be converted into a constant problem, so the non-constant flow problem can be analyzed using a constant flow equation. The calculation process is simplified, and the data calculation efficiency is improved.
According to some embodiments of the invention, a method of determining a momentum equation includes:
for the method shown in figure 3, the dynamic coordinate with the speed of the wave breaking speed w is taken from the water flow with the unsteady wave breaking flow to analyze, and the water body between the first section and the second section is taken as a control body; the total dynamic water pressure distribution of the first section and the second section is respectively P and P 0 As shown in fig. 5.
Because the control body is smaller, the gravity component and the boundary friction force acting on the control body are ignored, and a constant flow momentum equation can be applied to the control body along the water flow direction, so that the following results:
wherein, gammaFor the water weight, 9.8kN/m was taken 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, 9.8m/s 2 。
The technical scheme has the working principle and beneficial effects that: the first section is 1-1, and the second section is 2-2. The momentum equation is convenient to accurately determine.
According to some embodiments of the invention, a method of determining a breakup wave speed includes:
the wave velocity of the broken wave can be obtained according to the continuity equation and the momentum equation;
from the equation of continuity,
(v-w)A＝(v 0 -w)A 0
consider wave flow Δq=q-Q 0 And Δa=A-A 0 ，Δv＝v 0 V, the above formula can be written as:
(A 0 +ΔA)[w-(v 0 +Δv)]＝A 0 (w-v 0 )
after finishing, the method comprises the following steps:
according to the momentum equation,
wherein: Δp=p-P 0
Substituting the continuity equation into the above equation:
finishing to obtain wave velocity formula
In the above, the positive sign before the root corresponds to forward wave field, and the negative sign corresponds to backward wave field;
as shown in fig. 6, the water flow conditions approximately consider that the dynamic water pressure on the first section and the second section satisfy the static water pressure distribution rule, and the following results:
ΔP＝P-P 0 ＝γ(Ay c -A 0 y c0 )
wherein y is c Is the depth of the centroid of the area A under water, y c0 Is area A 0 The centroid of (c) is at the depth of the water; let y be c ' is the depth of the centroid of the area deltaa under water, and ζ is the height difference of the water surface before and after the peak, i.e. the wave height, according to the area moment property, for the area moment Ay c The method can obtain:
Ay c ＝(A 0 +ΔA)y c ＝A 0 y c0 +ΔAy c ′
at this time, the dynamic water pressure difference is as follows:
ΔP＝γ(A 0 y c0 +ΔAy c ′+A 0 ζ-A 0 y c0 )＝γ(ΔAy c ′+A 0 ζ)
substituting the above formula into a wave velocity formula to obtain:
assume that: Δa=ζb',
the above can be simplified as:
average water depth is introduced:obtaining:
in the case of a constant cross-sectional shape, v is 0 、h 0 、B 0 Equal are known, i.e. the wave speed is a single-valued function of the wave height.
According to some embodiments of the invention, a method of determining a broken wave basis equation comprises:
dam break wave belongs to forward wave, calculates based on forward wave:
that is to say,the right end is positively marked, and a wave breaking basic equation is obtained according to calculation:
according to some embodiments of the invention, the parameter information includes upstream section geometry parameters, downstream section geometry parameters, river parameters, and control calculation parameters;
the geometric parameters of the upstream section comprise the bottom elevation, the bottom width and the section slope ratio of the upstream section;
the geometric parameters of the downstream section comprise the bottom elevation, the bottom width and the section slope ratio of the downstream section;
the river parameters include: calculating time, natural runoff and Manning coefficients;
the control calculation parameters comprise an initial time interval, an evolution distance and a resistance coefficient.
The technical scheme has the working principle and beneficial effects that: as shown in fig. 8-9, the calculation of the parameters to be input includes four modules, respectively: (1) upstream section geometry; (2) downstream section geometry; (3) river parameters; (4) and controlling the calculated parameters.
The calculation is divided into two cases: firstly, flood/debris flow evolution analysis from a crumple to a downstream section is performed, wherein the bottom elevation, the bottom width and the crumple slope ratio of the crumple are input into the module (1); secondly, flood/debris flow evolution analysis is carried out between two sections, and at the moment, the bottom elevation, the bottom width and the river slope ratio of the upstream section are input into the module (1).
And (3) inputting geometric parameters of an upstream section of the flow process to be solved in the module (2), wherein the geometric parameters comprise bottom elevation, bottom width and river slope ratio.
In the module (3), natural runoff is filled in according to the actual perennial flow of the river to be analyzed, a Manning coefficient takes an empirical value according to the roughness of an actual river bed, the calculation time refers to the duration of an input flow process, the natural runoff is used when equal-interval wave breaking element division is carried out in an initial calculation step, and the program is self-determined.
In the module (4), the initial time interval is the time interval used when the time interval wave breaking element division is performed in the initial step, and the recommended value is 1/10-1/30 of the calculated time in the module (3); the evolution distance is the length of the axis of the river between the calculated section and the section of the burst or the section of the known flow process (upstream section); the physical meaning of the resistance coefficient is shown in an along-path head loss formula of the broken wave element wave, and plays a key role in the calculation result.
As shown in fig. 7, according to some embodiments of the present invention, acquiring flow process information generated by a breach position when a breach occurs includes:
dividing the flow process into n disconnected wave element waves according to a time average mode, wherein the numbers are respectively 1-n;
the number i of the broken element wave is that the average flow rate in the period is:can obtain the corresponding wave height ζ i ；
The total flow of the broken wave element wave with the number of i is the product of the average flow in the period and the time length, namely:
in the middle of：q i For the total flow of the ith broken element wave, m 3 ；Q i For the average flow of the ith broken element wave, m 3 S; dt is the time interval, s.
In a specific embodiment, the developed program is adopted to calculate the evolution of Tang Gushan barrier lake burst flood in an inversion way, and when Tang Gushan barrier lake burst, a North Sichuan hydrological station at 7km downstream of a dam address, a 33.5km through-port hydrological station and a Fujiang bridge hydrological station at 77km measure flow processes, and the flow processes are used as the basis of model calculation.
(1) Crumbling to north Sichuan station
The calculation is carried out by taking a trapezoid section according to the actual conditions of the crumple and the North Sichuan section, the roughness is 0.035 according to the actual measurement topography comparison and experience, the evolution distance is 7km, the initial time interval is 300s, and the specific values of all parameters are shown in figure 10. Through inversion analysis, when the resistance coefficient is 0.00007, the calculated result and the actual measurement result are better matched, as shown in fig. 11. The actual measured flood peak flow of the North Sichuan station is 6540m 3 And/s, obtaining the peak flood 6434.24m by calculation 3 And/s, the flood peak arrival time is also consistent.
(2) North Sichuan station to through port station
The calculation is carried out by taking a trapezoid section according to the actual conditions of the North Sichuan section and the through port section, the roughness is 0.035 according to the actual measurement topography comparison and experience, the evolution distance is 33.5km, the initial time interval is 300s, and the specific values of all parameters are shown in figure 12. Through inversion analysis, when the resistance coefficient is 0.00007, the calculated result and the actual measurement result are better matched, as shown in fig. 13. The actual measured flood peak flow of the through-port station is 6210m 3 And/s, obtaining the peak flood 6560.15m by calculation 3 And/s, the flood peak arrival time is also consistent.
(3) Through-port measuring station to Fu Jiang Qiaoce station
The actual conditions of the through-hole section and the Fujiang bridge section are taken as trapezoid sections during calculation, the roughness is taken as 0.035 according to actual measurement terrain comparison and experience, the evolution distance is 77km, the initial time interval is taken as 300s, and the specific values of all parameters are shown in figure 14. Through inversion analysis, when the resistance coefficient is 0.00007, the calculated result and the actual measured resultThe agreement is better, as shown in figure 15. The actual measured flood peak flow of the river bridge is 6100m 3 And/s, obtaining the peak flood 6106.27m by calculation 3 And/s, the flood peak arrival time is also consistent.
The developed program is adopted to simulate and calculate the flood evolutionary process of the Tangjia mountain engineering example, and the result shows that when the resistance coefficient is 0.00007, the calculation result is more consistent with the measured data, and compared with the model shown in the table 2-1, the model can simulate the whole evolutionary process better.
Table 1 flood evolution calculation result analysis
The flood evolution model is compiled based on a Mi crosoft Exce l platform, the interface is simple and easy to operate, the model parameters are few, the physical meaning is clear, the calculation speed is high, the result stability is good, the precision is high, and the flood evolution model has good application value to practical engineering.
According to some embodiments of the invention, obtaining the upstream cross-section geometry includes:
shooting a plurality of upstream section images based on the unmanned aerial vehicle;
performing image quality evaluation on a plurality of upstream section images, and determining a target upstream section image according to an evaluation result;
performing image recognition on the target upstream section image, and determining an upstream section geometric parameter according to a recognition result;
image quality evaluation is carried out on a plurality of upstream section images, and a target upstream section image is determined according to an evaluation result, wherein the method comprises the following steps:
performing standardization processing on the size of the upstream section image to obtain a processed upstream section image;
acquiring pixel points and pixel values on each processed upstream section image, and constructing a first pixel matrix;
calculating a difference value of the first pixel matrix and the second pixel matrix; the second pixel matrix is a matrix constructed based on pixel points and pixel values on the standard upstream section image;
wherein F is i A difference value between the second pixel matrix and the first pixel matrix constructed based on the ith processed upstream cross-sectional image;first pixel matrix A constructed for processing upstream cross-sectional image based on ith sheet i The pixel value of row t column s; c (C) s,t I=1, 2, 3 … … G for the pixel value of the s-th row t column of the second pixel matrix; g is the number of processed upstream sectional images; m is the number of pixels of the standard upstream cross-section image and the processed upstream cross-section image in the transverse direction; n is the number of pixels of the standard upstream section image and the processed upstream section image in the longitudinal direction;
and determining the processing upstream section image corresponding to the first pixel matrix with the minimum difference value as a target upstream section image.
The technical scheme has the working principle and beneficial effects that: acquiring the geometric parameters of the upstream section, including: shooting a plurality of upstream section images based on the unmanned aerial vehicle; performing image quality evaluation on a plurality of upstream section images, and determining a target upstream section image according to an evaluation result; and carrying out image recognition on the target upstream section image, and determining the upstream section geometric parameters according to recognition results. The unmanned aerial vehicle can obtain a plurality of upstream section images in the shooting process, selects the upstream section image with the best image quality as a target upstream section image, accurately determines the geometric parameters of the upstream section based on the image recognition technology, and saves time and labor. When the target upstream section image is determined, the size of the upstream section image is subjected to standardized processing to obtain a processed upstream section image, so that the size of the processed upstream section image is consistent with that of the standard upstream section image, namely the number of pixels on the two images is equal. Acquiring pixel points and pixel values on each processed upstream section image, and constructing a first pixel matrix; calculating a difference value of the first pixel matrix and the second pixel matrix; the second pixel matrix is a matrix constructed based on pixel points and pixel values on the standard upstream section image; the difference value is a quality difference representing the processed upstream cross-sectional image and the standard upstream cross-sectional image, and is represented by color, brightness, and the like. And determining the processing upstream section image corresponding to the first pixel matrix with the minimum difference value as a target upstream section image. And accurately screening out the processed upstream section image corresponding to the first pixel matrix with the smallest difference with the standard upstream section image, so that the accuracy of the subsequent image recognition is improved. Based on the formula, the difference value of the first pixel matrix and the second pixel matrix is accurately calculated, so that the processing upstream section image corresponding to the first pixel matrix with the minimum difference value is conveniently determined and used as the target upstream section image.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (10)
1. A dam break flood evolution prediction method based on a wave breaking principle is characterized by comprising the following steps:
constructing a dam break flood evolution model;
acquiring parameter information of flood evolution and flow process information generated at a breach position when dam break occurs;
inputting the parameter information and the flow process information generated by the position of the breach when the dam break occurs into a dam break flood evolution model, and outputting a flood evolution result.
2. The dam break flood evolutionary prediction method based on the wave breaking principle as claimed in claim 1, wherein constructing the dam break flood evolutionary model comprises:
constructing an initial model;
and adding an operation rule into the initial model to construct a dam-break flood evolution model.
3. The dam break flood prediction method based on the wave breaking principle as claimed in claim 2, wherein the operation rule comprises: continuity equation and wave flow, momentum equation, wave velocity of the wave break and basic equation of the wave break.
4. The dam break flood prediction method based on the wave breaking principle as claimed in claim 3, wherein the method for determining the continuity equation and the wave flow comprises the following steps:
setting the flow velocity of the unaffected area of the wave breaking crest as v by adopting fixed coordinates 0 Corresponding cross-sectional area A 0 The method comprises the steps of carrying out a first treatment on the surface of the The flow velocity of the area affected by the wave crest of the broken wave is v, the corresponding cross-sectional area is A, and the wave velocity of the broken wave is w; the velocity is the dynamic coordinate of the wave breaking velocity w parallel to the wave propagation direction, and the flow velocity before and after wave crest influence is v 0 -w, v-w, according to a constant flow continuity equation, obtainable:
(v-w)A＝(v 0 -w)A 0
vA-v 0 A 0 ＝w(A-A 0 )
introducing a flow Q at section h0 0 And the flow rate Q at section h, can be obtained:
Q-Q 0 ＝w(A-A 0 )
let Δq=q-Q 0 Is wave flow and the area difference deltaa=A-A before and after the peak 0 Determining a relation expression between the flow rate and the wave speed and the wave height of the broken wave:
ΔQ＝wΔA＝ζB′w
in the method, in the process of the invention,B 0 the initial water surface width of the section before the occurrence of the wave breaking; b is the water surface width of the section after the occurrence of the wave breaking; ζ=h-h 0, wave height.
5. The method for dam break flood evolutions prediction based on the wave breaking principle according to claim 4, wherein the method for determining the momentum equation comprises the following steps:
unsteady flow at breakAnalyzing a dynamic coordinate with the speed of the wave breaking speed w in water flow, and taking a water body between a first section and a second section as a control body; the total dynamic water pressure distribution of the first section and the second section is respectively P and P 0 ，
Because the control body is smaller, the gravity component and the boundary friction force acting on the control body are ignored, and a constant flow momentum equation can be applied to the control body along the water flow direction, so that the following results:
wherein gamma is the weight of water, 9.8kN/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, 9.8m/s 2 。
6. The dam break flood prediction method based on the wave breaking principle according to claim 5, wherein the method for determining the wave breaking speed comprises the following steps:
the wave velocity of the broken wave can be obtained according to the continuity equation and the momentum equation;
from the continuity equation:
(v-w)A＝(v 0 -w)A 0
consider wave flow Δq=q-Q 0 And Δa=A-A 0 ，Δv＝v 0 V, the above formula can be written as:
(A 0 +ΔA)[w-(v 0 +Δv)]＝A 0 (w-v 0 )
after finishing, the method comprises the following steps:
according to the momentum equation,
wherein: Δp=p-P 0
Substituting the continuity equation into the above equation:
finishing to obtain wave velocity formula
In the above, the positive sign before the root corresponds to forward wave field, and the negative sign corresponds to backward wave field;
approximately considering that the dynamic water pressure intensity on the first section and the second section meets the distribution rule of static water pressure intensity, obtaining:
ΔP＝P-P 0 ＝γ(Ay c -A 0 y c0 )
wherein y is c Is the depth of the centroid of the area A under water, y c0 Is area A 0 The centroid of (c) is at the depth of the water; let y be c ' is the depth of the centroid of the area deltaa under water, and ζ is the height difference of the water surface before and after the peak, i.e. the wave height, according to the area moment property, for the area moment Ay c The method can obtain:
Ay c ＝(A 0 +ΔA)y c ＝A 0 y c0 +ΔAy c ′
at this time, the dynamic water pressure difference is as follows:
ΔP＝γ(A 0 y c0 +ΔAy c ′+A 0 ζ-A 0 y c0 )＝γ(ΔAy c ′+A 0 ζ)
substituting the above formula into a wave velocity formula to obtain:
assume that: Δa=ζb',
the above can be simplified as:
average water depth is introduced:obtaining:
in the case of a constant cross-sectional shape, v is 0 、h 0 、B 0 Equal are known, i.e. the wave speed is a single-valued function of the wave height.
7. The dam break flood prediction method based on the principle of wave breaking as claimed in claim 6, wherein the method for determining the basic equation of wave breaking comprises the following steps:
dam break wave belongs to forward wave, calculates based on forward wave:
and obtaining a wave breaking basic equation according to the calculation.
8. The dam break flood evolutionary prediction method based on the wave breaking principle as claimed in claim 1, wherein the parameter information comprises an upstream section geometric parameter, a downstream section geometric parameter, a river channel parameter and a control calculation parameter;
the geometric parameters of the upstream section comprise the bottom elevation, the bottom width and the section slope ratio of the upstream section;
the geometric parameters of the downstream section comprise the bottom elevation, the bottom width and the section slope ratio of the downstream section;
the river parameters include: calculating time, natural runoff and Manning coefficients;
the control calculation parameters comprise an initial time interval, an evolution distance and a resistance coefficient.
9. The method for predicting dam break flood evolutions based on the wave breaking principle as claimed in claim 1, wherein the step of obtaining the flow process information generated by the position of the break when the dam break occurs comprises the steps of:
dividing the flow process into n disconnected wave element waves according to a time average mode, wherein the numbers are respectively 1-n;
the number i of the broken element wave is that the average flow rate in the period is:can obtain the corresponding wave height ζ i ；
The total flow of the broken wave element wave with the number of i is the product of the average flow in the period and the time length, namely:
wherein: q i For the total flow of the ith broken element wave, m 3 ；Q i For the average flow of the ith broken element wave, m 3 S; dt is the time interval, s.
10. The method for predicting dam break flood evolutions based on the wave breaking principle as claimed in claim 8, wherein the step of obtaining the geometric parameters of the upstream section comprises the steps of:
shooting a plurality of upstream section images based on the unmanned aerial vehicle;
performing image quality evaluation on a plurality of upstream section images, and determining a target upstream section image according to an evaluation result;
performing image recognition on the target upstream section image, and determining an upstream section geometric parameter according to a recognition result;
image quality evaluation is carried out on a plurality of upstream section images, and a target upstream section image is determined according to an evaluation result, wherein the method comprises the following steps:
performing standardization processing on the size of the upstream section image to obtain a processed upstream section image;
acquiring pixel points and pixel values on each processed upstream section image, and constructing a first pixel matrix;
calculating a difference value of the first pixel matrix and the second pixel matrix; the second pixel matrix is a matrix constructed based on pixel points and pixel values on the standard upstream section image;
wherein F is i A difference value between the second pixel matrix and the first pixel matrix constructed based on the ith processed upstream cross-sectional image;first pixel matrix A constructed for processing upstream cross-sectional image based on ith sheet i The pixel value of row t column s; c (C) s,t I=1, 2, 3 … … G for the pixel value of the s-th row t column of the second pixel matrix; g is the number of processed upstream sectional images; m is the number of pixels of the standard upstream cross-section image and the processed upstream cross-section image in the transverse direction; n is standard upstream section image and processingThe number of pixels of the upstream cross-sectional image in the longitudinal direction;
and determining the processing upstream section image corresponding to the first pixel matrix with the minimum difference value as a target upstream section image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310087495.5A CN116663223B (en) | 2023-02-09 | 2023-02-09 | Dam break flood evolution prediction method based on wave breaking principle |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310087495.5A CN116663223B (en) | 2023-02-09 | 2023-02-09 | Dam break flood evolution prediction method based on wave breaking principle |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116663223A true CN116663223A (en) | 2023-08-29 |
| CN116663223B CN116663223B (en) | 2024-05-03 |
Family
ID=87717730
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310087495.5A Active CN116663223B (en) | 2023-02-09 | 2023-02-09 | Dam break flood evolution prediction method based on wave breaking principle |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116663223B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102034001A (en) * | 2010-12-16 | 2011-04-27 | 南京大学 | Design method for distributed hydrological model by using grid as analog unit |
| CN102445239A (en) * | 2010-09-29 | 2012-05-09 | 联宇工程技术(武汉)有限公司 | Novel water metering method for open channel based on multi-point water level |
| CN104091065A (en) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | Intermittent flow numerical simulation method for solving shallow water problem |
| CN110219277A (en) * | 2018-03-26 | 2019-09-10 | 中国水利水电科学研究院 | The dam bursting flood analysis system of dam system |
| CN111089625A (en) * | 2019-12-13 | 2020-05-01 | 国网浙江省电力有限公司紧水滩水力发电厂 | Binocular vision-simulated river flow real-time monitoring system and method |
| CN114858226A (en) * | 2022-07-05 | 2022-08-05 | 武汉大水云科技有限公司 | Unmanned aerial vehicle torrential flood flow measuring method, device and equipment |
| CN114897242A (en) * | 2022-05-10 | 2022-08-12 | 中国长江电力股份有限公司 | Method for adaptively segmenting rainfall runoff in time period field |
| CN115409741A (en) * | 2022-11-01 | 2022-11-29 | 长江勘测规划设计研究有限责任公司 | Machine vision recognition algorithm for measuring sediment content by using river surface color difference |
-
2023
- 2023-02-09 CN CN202310087495.5A patent/CN116663223B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102445239A (en) * | 2010-09-29 | 2012-05-09 | 联宇工程技术(武汉)有限公司 | Novel water metering method for open channel based on multi-point water level |
| CN102034001A (en) * | 2010-12-16 | 2011-04-27 | 南京大学 | Design method for distributed hydrological model by using grid as analog unit |
| CN104091065A (en) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | Intermittent flow numerical simulation method for solving shallow water problem |
| CN110219277A (en) * | 2018-03-26 | 2019-09-10 | 中国水利水电科学研究院 | The dam bursting flood analysis system of dam system |
| CN111089625A (en) * | 2019-12-13 | 2020-05-01 | 国网浙江省电力有限公司紧水滩水力发电厂 | Binocular vision-simulated river flow real-time monitoring system and method |
| CN114897242A (en) * | 2022-05-10 | 2022-08-12 | 中国长江电力股份有限公司 | Method for adaptively segmenting rainfall runoff in time period field |
| CN114858226A (en) * | 2022-07-05 | 2022-08-05 | 武汉大水云科技有限公司 | Unmanned aerial vehicle torrential flood flow measuring method, device and equipment |
| CN115409741A (en) * | 2022-11-01 | 2022-11-29 | 长江勘测规划设计研究有限责任公司 | Machine vision recognition algorithm for measuring sediment content by using river surface color difference |
Non-Patent Citations (2)
| Title |
|---|
| 姚志坚 等: "溃坝洪水演进数学模型研究及其在大镜山水库的应用", 《人民珠江》, pages 1 - 3 * |
| 曹伟: "平面二维数学模型在土石坝溃坝洪水演进中的应用", 《山西水利》, pages 1 - 4 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116663223B (en) | 2024-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110287571B (en) | River dangerous work scouring safety analysis and bank slope stability determination method | |
| CN112464584B (en) | Method for estimating water level and flow of free surface flow | |
| Salmasi et al. | Neuro-fuzzy approach for estimating energy dissipation in skimming flow over stepped spillways | |
| CN111241757A (en) | Dam break three-dimensional numerical simulation method for uranium tailing pond based on computational fluid mechanics | |
| Torabi et al. | The effect of geometry parameters and flow characteristics on erosion and sedimentation in channels junction using finite volume method | |
| CN114781136A (en) | Method for simulating landslide and surge based on DEM-CFD-SPH modeling | |
| CN110263449B (en) | Numerical simulation method for landslide surge of landslide deformation | |
| KR101106548B1 (en) | Method for analyzing shallow water flow using the two-dimensional river flow model with tensor-type eddy viscosity | |
| CN116663223B (en) | Dam break flood evolution prediction method based on wave breaking principle | |
| CN114357903A (en) | Alluvial river resistance calculation and simulation method and device based on riverbed morphological parameters | |
| Martino et al. | Double-average mean flow and local turbulence intensity profiles from PIV measurements for an open channel flow with rigid vegetation | |
| Patel et al. | Flow and bedform dynamics in an alluvial channel with downward seepage | |
| Yeow et al. | Wind tunnel analysis of the detachment bubble on Bolund Island | |
| CN110847112B (en) | River flood discharge early warning method based on hydraulics simulation | |
| Zhou et al. | Experimental investigations of propagation characteristics and wave energy of dam-break waves on wet bed | |
| CN116562043A (en) | Data assimilation research method for three-dimensional hydrodynamic simulation of natural river reach | |
| CN108827871B (en) | Method for determining shear stress of sediment surface in tubular sediment erosion test device | |
| Mulahasan et al. | Effect of blockage ratio on flow characteristics in obstructed open channels | |
| CN116305731A (en) | Land and stone dam break analysis method based on cellular automaton principle | |
| CN115358166A (en) | HLL format-based one-dimensional open channel hydrodynamic process simulation method | |
| CN115391717A (en) | Prediction method for overtopping and bursting flood peak flow of earth and rockfill dam | |
| CN116384266B (en) | Mud-rock flow evolution prediction method based on wave-breaking principle | |
| Li et al. | Lagrangian modeling of tidal bores passing through bridge piers | |
| Behera et al. | Scour at the downstream of Ghatakeswar spillway using non-cohesive hydraulic model | |
| Dal et al. | Validation of dam-break problem over dry bed using SPH |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |