CN116911215A

CN116911215A - Dam break flood simulation method and system based on reservoir flood regulation coupling hydrodynamic force

Info

Publication number: CN116911215A
Application number: CN202311154591.3A
Authority: CN
Inventors: 廖添; 未荣莲; 闵皆昇; 周璐; 吴健明
Original assignee: Zhejiang Yuansuan Technology Co ltd
Current assignee: Zhejiang Yuansuan Technology Co ltd
Priority date: 2023-09-08
Filing date: 2023-09-08
Publication date: 2023-10-20
Anticipated expiration: 2043-09-08
Also published as: CN116911215B

Abstract

The invention discloses a dam break flood simulation method and system based on reservoir flood control coupling hydrodynamic force, and belongs to the technical field of reservoir flood control. The existing burst flood analysis scheme does not relate to how to obtain model parameters of a flood evolution analysis module, and influences accurate simulation of burst flood. According to the dam-break flood simulation method based on reservoir flood control coupling hydrodynamic force, disclosed by the invention, the dam-break flood simulation based on reservoir flood control coupling hydrodynamic force is realized by constructing the terrain generalization model, the reservoir flood control model, the reservoir parameter calibration model and the dam-break coupling simulation model, so that the accuracy of a dam-break flood simulation scheme is effectively improved, a user can accurately simulate the dam-break flood, and popularization and use of the dam-break flood simulation scheme are facilitated. The invention can realize the model parameter calibration of the reservoir flood control model, and can obtain different model parameters according to different reservoirs and different reservoir generalized topographic information, thereby effectively reducing the calculation error of the dam break simulation scheme.

Description

Dam break flood simulation method and system based on reservoir flood regulation coupling hydrodynamic force

Technical Field

The invention relates to a dam break flood simulation method and system based on reservoir flood control coupling hydrodynamic force, and belongs to the technical field of reservoir flood control.

Background

Chinese patent (CN 110219277B) discloses a dam-breaking flood analysis system comprising: the system comprises a topology analysis module, a data collection module, an analysis module of a single dam, a flood evolution analysis module and a prediction module of the dam system; the topology analysis module of the dam system is used for determining the grading of the dams in the dam system; the data collection module is used for collecting information of dams in the dam system; the analysis module of the single dam is used for establishing a water level model of the dam and determining a bursting state model of the dam according to the water level model of the dam and building parameter data of the dam body; the flood evolution analysis module is used for analyzing the process of flood flowing from the low-grade dam to the high-grade dam, and calculating the inflow rate of the downstream high-grade dam according to the outflow rate of the low-grade dam; the prediction module is used for predicting whether the downstream dam is broken and predicting the breaking time of the downstream broken dam and the flood peak flow of broken flood if the dam is broken in the evolution process.

However, according to the above-mentioned burst flood analysis scheme, the flood evolution analysis module is set to calculate the inflow of the downstream high-grade dam, but the model parameters of the flood evolution analysis module are not obtained, and how to rate the model parameters is not disclosed, so that a large error may exist in the burst flood analysis scheme, and further, a user cannot accurately simulate burst flood, which is not beneficial to popularization and use of the burst flood simulation scheme.

Disclosure of Invention

Aiming at the problems or one of the problems, the invention aims to provide the dam break flood simulation method based on the reservoir flood control coupling hydrodynamic force, which is used for realizing the dam break flood simulation based on the reservoir flood control coupling hydrodynamic force by constructing a terrain generalization model, a reservoir flood control model, a reservoir parameter rating model and a dam break coupling simulation model, thereby effectively improving the accuracy of a dam break flood simulation scheme, facilitating the accurate simulation of the dam break flood by users and being beneficial to popularization and use of the dam break flood simulation scheme.

Aiming at the problems or one of the problems, the invention aims to provide a dam break flood simulation system based on the reservoir flood control coupling hydrodynamic force, which is used for realizing the dam break flood simulation based on the reservoir flood control coupling hydrodynamic force by arranging a terrain generalization module, a reservoir simulation module, a reservoir flood control module, a reservoir parameter rating module and a dam break coupling simulation module, thereby effectively reducing the calculation error of the dam break simulation scheme, improving the accuracy of the dam break flood simulation scheme and facilitating the accurate simulation of the dam break flood by a user.

Aiming at the problems or one of the problems, the invention aims to provide a dam break flood simulation method and system based on reservoir flood control coupling hydrodynamic force, which can realize model parameter calibration of a reservoir flood control model and obtain different model parameters according to different reservoirs and different reservoir generalized topographic information, so that calculation errors of a dam break simulation scheme can be effectively reduced, and the scheme is simple and practical.

In order to achieve one of the above objects, a first technical solution of the present invention is:

a dam break flood simulation method based on reservoir flood regulation coupling hydrodynamic force comprises the following steps:

acquiring a certain reservoir parameter to be dam-break simulated;

processing reservoir parameters by utilizing a pre-constructed topographic generalized model to obtain reservoir generalized topographic information;

based on a reservoir flood regulation model constructed in advance, carrying out flood regulation calculation on reservoir generalized topography information according to a reservoir scheduling principle and an initial roughness rate to obtain a flood regulation simulation result;

the model parameters of the reservoir flood control model are calibrated according to the pre-constructed reservoir parameter calibration model, the flood control simulation result and reservoir flood control historical data, and the optimal model parameters of the reservoir flood control model are obtained;

According to the optimal model parameters, the reservoir flood regulating model takes a storage flow process line brought by dam break flood as a boundary condition, and flood regulating calculation is carried out to obtain a total reservoir outlet flow process line;

and carrying out two-dimensional continuous dam-break hydrodynamic coupling simulation according to the total reservoir outlet flow process line by utilizing a pre-constructed dam-break coupling simulation model to obtain a dam-break simulation result, thereby realizing dam-break flood simulation based on reservoir flood-regulating coupling hydrodynamic force.

Through continuous exploration and test, reservoir parameters are processed by constructing a topographic generalized model, a reservoir flood regulating model, a reservoir parameter rating model and a dam break coupling simulation model, so that reservoir generalized topographic information is obtained; according to reservoir dispatching principle and initial roughness, flood dispatching calculation is carried out on reservoir generalized topographic information, and flood dispatching simulation results are obtained; thirdly, according to the flood control simulation result and reservoir flood control historical data, calibrating model parameters of a reservoir flood control model to obtain optimal model parameters of the reservoir flood control model; according to the optimal model parameters, taking a warehouse entry flow process line brought by dam break flood as a boundary condition, and carrying out flood regulating calculation to obtain a total reservoir outlet flow process line; and then, according to the total outlet flow process line of the reservoir, carrying out two-dimensional continuous dam-break hydrodynamic coupling simulation to obtain a dam-break simulation result, and realizing dam-break flood simulation based on reservoir flood-regulation coupling hydrodynamic, thereby effectively improving the accuracy of a dam-break flood simulation scheme, facilitating accurate simulation of dam-break flood by users, facilitating popularization and use of the dam-break flood simulation scheme, and being scientific, reasonable and practical.

Furthermore, the invention realizes the model parameter calibration of the reservoir flood regulating model by setting and constructing the terrain generalization model and the reservoir parameter calibration model, and can obtain different model parameters according to different reservoirs and different reservoir generalized terrain information, thereby effectively reducing the calculation error of a dam break simulation scheme, and having simple scheme and practicality.

As a preferred technical measure:

the construction method of the terrain generalization model comprises the following steps:

step 1: acquiring a certain reservoir parameter, wherein the certain reservoir parameter at least comprises peripheral contour data of a certain reservoir area;

step 2: establishing a reservoir polygon according to the peripheral outline data, and calculating the area of the reservoir polygon to obtain the reservoir area under the normal water level;

step 3: generating a water level array according to the water level resolution and the elevation range of the reservoir region terrain;

step 4: based on the reservoir water level and reservoir capacity curves, solving reservoir capacities corresponding to each water level in the water level array; calculating the area of the reservoir corresponding to each water level by adopting a prismatic table volume formula in combination with the area of the reservoir under the normal water level in the step 2 to form a water level and area mapping table;

step 5: traversing the water level array, and performing offset processing on the polygon in the reservoir area until the area of the processed polygon is consistent with the corresponding area in the water level and area mapping table in the step 4; at this time, the contour line of the polygon is taken as the contour line of the current water level, and the generalized topography of the reservoir is obtained.

According to the invention, the water level and area mapping table is calculated by using the prismatic table volume formula through the reservoir region polygon and the reservoir water level and reservoir capacity curve, the nonlinear equation is solved in an iterative manner, and the reservoir two-dimensional generalized topography is established through the offset curve algorithm, so that the scheme is simple and practical.

As a preferred technical measure:

the method for obtaining reservoir generalized topographic information by using the topographic model comprises the following steps:

step 21, acquiring a certain reservoir parameter, wherein the certain reservoir parameter at least comprises peripheral contour data of a certain reservoir area;

establishing a reservoir polygon according to the peripheral outline data, wherein the reservoir polygon is a two-dimensional array of point coordinates sequentially arranged in a clockwise or anticlockwise direction; calculating the area of the polygon in the reservoir area by adopting a Gaussian area formula; taking the area of the polygon of the reservoir area as the area of the reservoir area under the normal water storage level;

step 22, obtaining water level resolution as the water level difference of two adjacent contour lines in generalized topography;

acquiring an elevation range of the topography of the reservoir area as a minimum elevation value and a maximum elevation value of the generalized topography;

then, the minimum elevation value is used as the first element of the water level array, the water level difference is used as the increment, and the maximum elevation value is used as the last element, so that a water level array is generated;

Step 23, calculating reservoir capacity corresponding to each water level in the water level array based on the reservoir water level and the reservoir capacity curve; calculating the area of the reservoir corresponding to each water level by adopting a prismatic table volume formula in combination with the area of the reservoir under the normal water level in the step 21 to form a water level and area mapping table;

step 24, traversing the water level array, and performing offset processing on the polygon in the reservoir area until the area of the processed polygon is consistent with the corresponding area in the water level and area mapping table in step 23; at the moment, the polygonal contour line can be used as the contour line of the current water level to obtain reservoir generalized topographic information; the method comprises the following specific steps:

s41: traversing the water level array, and solving the offset of the contour polygon corresponding to each water level relative to the reservoir polygon; the area of the polygon after offset is a function of the original polygon area and the enlarged distance;

the function value corresponding to each amplifying distance is obtained by an offset curve algorithm and a Gaussian area formula; traversing the water level array, and iteratively solving a nonlinear equation about the amplified distance to obtain the offset of each water level;

s42: according to the offset of each water level, using an offset curve algorithm to offset the polygon in the reservoir area to obtain a contour line of the water level;

S43: extracting vertex coordinates of each contour polygon; and according to the vertex coordinates, taking the water level value corresponding to the polygon as the elevation value of the vertex, and constructing three-dimensional terrain elevation data of the reservoir, thereby obtaining generalized terrain information of the reservoir.

As a preferred technical measure:

the system also comprises a reservoir simulation model; the reservoir simulation model is used for updating reservoir generalized topography information according to hydrological data to obtain reservoir generalized new topography;

the method for obtaining the reservoir generalized new topography by the reservoir simulation model comprises the following steps:

step 31: according to a certain reservoir parameter, an irregular triangular grid with a plurality of grid nodes is established;

step 32: mapping the three-dimensional terrain elevation data onto grid nodes using bilinear interpolation;

step 33: calculating the model water quantity when the model water level is the characteristic water level according to the grid nodes;

step 34: comparing the model water quantity with the water level and the reservoir capacity in the reservoir capacity curve to obtain an error result;

step 35: when the error result is smaller than a given threshold value, the inventory rechecking is completed;

when the error result is greater than or equal to a given threshold value, calculating elevation adjustment quantity according to a prism volume formula and the model water quantity and the prism volume formula by using a prism volume formula, and generating new three-dimensional terrain elevation data;

Step 36: and updating the reservoir generalized topography information according to the new three-dimensional topography elevation data to obtain the reservoir generalized new topography.

As a preferred technical measure:

the method for obtaining the flood control simulation result by the reservoir flood control model is as follows:

step 41: acquiring the surface type and the condition of the reservoir based on the reservoir generalized new topography;

step 42: setting initial roughness according to the surface type and condition of the reservoir, and taking reservoir design information and a checking and warehousing flow time sequence process line as upstream boundary conditions;

step 43: establishing a model downstream boundary condition according to a reservoir dispatching principle;

the reservoir dispatching principle comprises the following contents:

determining a delivery condition according to the reservoir delivery flow condition and the dam front water level condition;

if the ex-warehouse condition is to realize constant ex-warehouse flow, setting a fixed flow boundary;

if the ex-warehouse condition is that the reservoir storage flow is equal to the ex-warehouse flow, setting a fixed water level boundary;

if the delivery conditions are that the delivery flow is calculated according to the discharge curve of the reservoir flood discharge building, setting a water level-flow boundary;

step 44: according to the downstream boundary condition of the model, flood diversion simulation calculation is carried out, and a flood diversion simulation result is obtained;

the flood regulating simulation result at least comprises a dam front water level time sequence process line and a reservoir lower discharge flow process.

As a preferred technical measure:

the method for obtaining the optimal model parameters of the reservoir flood control model by the reservoir parameter calibration model comprises the following steps:

step 51: obtaining a flood control simulation result and reservoir flood control historical data;

the flood regulating simulation result at least comprises a dam front water level time sequence process line and a reservoir lower discharge flow process;

step 52: calculating maximum value error, peak time error and maximum lower discharge flow error of the dam front water level according to the flood regulating simulation result and reservoir flood regulating historical data;

step 53: normalizing the maximum value error, the peak time error and the maximum drainage flow error of the dam front water level, and taking the square sum of the relative errors as an objective function;

step 54: and solving the minimum value of the objective function by using an iteration method and taking the minimum value as the optimal model parameter of the reservoir flood regulating model.

According to the method, the objective function is established according to the peak value of the water level in front of the dam, the peak current time and the maximum drainage flow of the reservoir, the minimum value of the flood regulating model of the water reservoir is calculated in an iterative mode, the calibration and optimization of the model parameters are realized, and the scheme is practical.

As a preferred technical measure:

the method for carrying out flood control calculation by the reservoir flood control model comprises the following steps:

step 61: according to the reservoir generalized new topography, the reservoir scheduling principle and the optimal model parameters, and taking a warehousing flow process line caused by upstream dam break flood as an upstream boundary condition, carrying out flood control calculation to obtain dam front water level information;

Step 62: judging whether a top-overflow condition or a continuous burst condition occurs according to the front water level information of the dam;

if the overtaking situation occurs, calculating to obtain the overtaking flow by using a wide top weir formula, and executing a step 63;

if continuous burst occurs, ending the flood regulating operation;

step 63: calculating the total reservoir outlet flow according to the flood peak flow and the reservoir outlet boundary flow;

step 64: and obtaining a reservoir total outlet flow process line according to the reservoir total outlet flow.

As a preferred technical measure:

the method for obtaining the dam break simulation result by the dam break coupling simulation model comprises the following steps:

when the reservoirs are not burst, calculating dam site flow according to the total outlet flow process line of one or more reservoirs;

when the reservoir reaches the breaking condition, dam site flow is calculated in real time according to the local dam front water level, the dam lower water level and the topography, and dam breaking flood simulation based on reservoir flood regulating coupling hydrodynamic force is realized.

In order to achieve one of the above objects, a second technical solution of the present invention is:

firstly, acquiring a reservoir parameter to be dam-break simulated;

secondly, processing reservoir parameters by utilizing a pre-constructed topographic generalized model to obtain reservoir generalized topographic information;

Thirdly, updating reservoir generalized topography information by a reservoir simulation model built in advance and combining hydrologic data to obtain new reservoir generalized topography;

fourthly, based on a reservoir flood control model constructed in advance, carrying out flood control calculation on the reservoir generalized new topography according to a reservoir scheduling principle and an initial roughness rate to obtain a flood control simulation result;

fifthly, utilizing a pre-constructed reservoir parameter calibration model, and calibrating model parameters of the reservoir flood control model according to a flood control simulation result and reservoir flood control historical data to obtain optimal model parameters of the reservoir flood control model;

sixthly, according to the optimal model parameters, the reservoir flood control model takes a storage flow process line brought by dam break flood as a boundary condition, and performs flood control calculation to obtain a total reservoir outlet flow process line;

and seventhly, carrying out two-dimensional continuous dam-break hydrodynamic coupling simulation according to the total reservoir outlet flow process line of the reservoir by utilizing a pre-constructed dam-break coupling simulation model to obtain a dam-break simulation result, and realizing dam-break flood simulation based on reservoir flood-regulating coupling hydrodynamic.

Through continuous exploration and experiments, reservoir parameters are processed by constructing a terrain generalization model, a reservoir simulation model, a reservoir flood regulation model, a reservoir parameter rating model and a dam break coupling simulation model, so that reservoir generalized terrain information is obtained; and combining with hydrologic data to update the reservoir generalized topography information to obtain a new reservoir generalized topography; further, according to a reservoir dispatching principle and an initial roughness, flood dispatching calculation is carried out on the reservoir generalized new topography to obtain a flood dispatching simulation result; thirdly, according to the flood control simulation result and reservoir flood control historical data, calibrating model parameters of a reservoir flood control model to obtain optimal model parameters of the reservoir flood control model; according to the optimal model parameters, taking a warehouse entry flow process line brought by dam break flood as a boundary condition, and carrying out flood regulating calculation to obtain a total reservoir outlet flow process line; and then, according to the total outlet flow process line of the reservoir, carrying out two-dimensional continuous dam-break hydrodynamic coupling simulation to obtain a dam-break simulation result, and realizing dam-break flood simulation based on reservoir flood-regulation coupling hydrodynamic, thereby effectively improving the accuracy of a dam-break flood simulation scheme, facilitating accurate simulation of dam-break flood by users, facilitating popularization and use of the dam-break flood simulation scheme, and being scientific, reasonable and practical.

Still further, the invention establishes the reservoir simulation model and updates the reservoir generalized topography information by combining with hydrological data to obtain the reservoir generalized new topography, further reduces the calculation error of the dam break simulation scheme and improves the accuracy of the dam break flood simulation scheme.

In order to achieve one of the above objects, a third technical solution of the present invention is:

the dam break flood simulation system based on reservoir flood control coupling hydrodynamic force adopts the dam break flood simulation method based on reservoir flood control coupling hydrodynamic force, and comprises a terrain generalization module, a reservoir simulation module, a reservoir flood control module, a reservoir parameter calibration module and a dam break coupling simulation module:

the terrain generalization module is used for processing the reservoir parameters to obtain reservoir generalized terrain information;

the reservoir simulation module is used for updating reservoir generalized topography information by combining hydrologic data to obtain reservoir generalized new topography;

The reservoir flood regulating module is used for carrying out flood regulating calculation according to the reservoir generalized new topography, the reservoir scheduling principle and the model parameters to obtain a reservoir total delivery flow process line;

the reservoir parameter calibration module is used for calibrating the model parameters of the reservoir flood regulating module to obtain the optimal model parameters of the reservoir flood regulating module;

and the dam break coupling simulation module is used for carrying out two-dimensional continuous dam break hydrodynamic coupling simulation according to the total reservoir outlet flow process line of the reservoir to obtain a dam break simulation result and realize dam break flood simulation based on reservoir flood regulating coupling hydrodynamic force.

Through continuous exploration and experiments, the dam-break flood simulation based on reservoir flood-regulation coupling hydrodynamic force is realized by arranging the terrain generalization module, the reservoir simulation module, the reservoir flood-regulation module, the reservoir parameter rating module and the dam-break coupling simulation module, so that the calculation error of a dam-break simulation scheme is effectively reduced, the accuracy of the dam-break flood simulation scheme is improved, a user can accurately simulate dam-break flood, popularization and use of the dam-break flood simulation scheme are facilitated, and the scheme is scientific, reasonable and practical.

Compared with the prior art, the invention has the following beneficial effects:

Through continuous exploration and test, reservoir parameters are processed by constructing a topographic generalized model, a reservoir flood regulating model, a reservoir parameter rating model and a dam break coupling simulation model, so that reservoir generalized topographic information is obtained; according to reservoir dispatching principle and initial roughness, flood dispatching calculation is carried out on reservoir generalized topographic information, and flood dispatching simulation results are obtained; thirdly, according to the flood control simulation result and reservoir flood control historical data, calibrating model parameters of a reservoir flood control model to obtain optimal model parameters of the reservoir flood control model; according to the optimal model parameters, taking a warehouse entry flow process line brought by dam break flood as a boundary condition, and carrying out flood regulating calculation to obtain a total reservoir outlet flow process line; and then, according to the total outlet flow process line of the reservoir, carrying out two-dimensional continuous dam-break hydrodynamic simulation calculation to obtain a dam-break simulation result, and realizing dam-break flood simulation based on reservoir flood-regulating coupling hydrodynamic force, thereby effectively improving the accuracy of a dam-break flood simulation scheme, facilitating accurate simulation of dam-break flood by users, facilitating popularization and use of the dam-break flood simulation scheme, and being scientific, reasonable and practical.

Furthermore, through continuous exploration and experiments, the dam-break flood simulation based on the reservoir flood-regulating coupling hydrodynamic force is realized by arranging the terrain generalization module, the reservoir simulation module, the reservoir flood-regulating module, the reservoir parameter rating module and the dam-break coupling simulation module, so that the calculation error of the dam-break simulation scheme is effectively reduced, the accuracy of the dam-break flood simulation scheme is improved, the accurate simulation of dam-break flood is facilitated for users, the popularization and the use of the dam-break flood simulation scheme are facilitated, and the scheme is scientific, reasonable and practical.

Drawings

FIG. 1 is a schematic flow chart of a dam break flood simulation method of the invention;

FIG. 2 is a schematic flow chart of another dam-break flood simulation method of the invention;

FIG. 3 is a schematic flow chart of a dam break flood simulation method of the present invention;

FIG. 4 is a schematic view of an irregular triangular mesh of a downstream reservoir according to an embodiment of the present invention;

fig. 5 is a schematic diagram of a flood diversion calculation result for a typical flood process according to an embodiment of the present invention;

fig. 6 is a schematic diagram of a flood control calculation result for a dam break flood according to an embodiment of the present invention.

Detailed Description

The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.

On the contrary, the invention is intended to cover any alternatives, modifications, equivalents, and variations as may be included within the spirit and scope of the invention as defined by the appended claims. Further, in the following detailed description of the present invention, certain specific details are set forth in order to provide a better understanding of the present invention. The present invention will be fully understood by those skilled in the art without the details described herein.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The term "or/and" as used herein includes any and all combinations of one or more of the associated listed items.

As shown in fig. 1, the dam break flood simulation method based on reservoir flood control coupling hydrodynamic force of the invention is a first specific embodiment:

acquiring a certain reservoir parameter to be dam-break simulated;

according to the optimal model parameters, the reservoir flood regulating model takes a storage flow process line brought by upstream dam break flood as a boundary condition, and flood regulating calculation is carried out to obtain a reservoir total delivery flow process line;

and carrying out two-dimensional continuous dam-break hydrodynamic simulation calculation according to the total reservoir outlet flow process line by utilizing a pre-constructed dam-break coupling simulation model to obtain a dam-break simulation result, thereby realizing dam-break flood simulation based on reservoir flood-regulating coupling hydrodynamic force.

The invention relates to a dam break flood simulation method based on reservoir flood control coupling hydrodynamic force, which comprises the following steps:

firstly, acquiring a reservoir parameter to be dam-break simulated;

sixthly, according to the optimal model parameters, the reservoir flood control model uses a storage flow process line brought by upstream dam break flood as a boundary condition to carry out flood control calculation, and a reservoir total outlet flow process line is obtained;

and seventhly, performing two-dimensional continuous dam break hydrodynamic simulation calculation according to a pre-constructed dam break coupling simulation model and the total reservoir outlet flow process line to obtain a dam break simulation result, and realizing dam break flood simulation based on reservoir flood control coupling hydrodynamic force.

As shown in fig. 2, the dam break flood simulation method based on reservoir flood control coupling hydrodynamic force of the invention is a third specific embodiment:

a dam-break flood simulation method based on reservoir flood-regulating coupling hydrodynamic force comprises the steps of carrying out topography generalization on reservoirs in a convection area, then establishing a flood-regulating model of each reservoir, carrying out reservoir capacity rechecking and parameter calibration of the model, carrying out flood-regulating calculation by using a flood-regulating model by taking a storage flow process line brought by upstream dam-break flood as an upstream boundary condition, and carrying out continuous dam-break hydrodynamic force simulation calculation based on the reservoir total discharge flow process line. The method specifically comprises the following steps:

s1: and establishing a reservoir polygon according to the peripheral outline of the reservoir, and calculating the area of the reservoir polygon to obtain the reservoir area under the normal water storage level.

S2: and generating a water level array according to the water level resolution and the elevation range of generalized terrain.

S3: and solving the reservoir capacity corresponding to each water level in the water level array based on the reservoir water level and the reservoir capacity curve. And (3) calculating the area of the reservoir corresponding to each water level by adopting a prismatic table volume formula in combination with the area of the reservoir under the normal water level in the step (S1) to form a water level and area mapping table.

S4: and traversing the water level array, and performing offset processing on the polygon in the reservoir area until the area of the processed polygon is consistent with the water level in the S3 and the corresponding area in the area mapping table. At this time, the polygon can be used as the contour line of the current water level to obtain the generalized topography of the reservoir.

S5: and (3) establishing an irregular triangular grid in the reservoir range, endowing the grid node topography values, constructing a reservoir simulation model, and rechecking the model reservoir capacity under the characteristic water level of the generalized topography to update the reservoir generalized topography.

S6: and (3) setting an initial roughness rate based on the surface type and the condition of the reservoir, and carrying out flood control calculation according to a reservoir scheduling principle to obtain a flood control simulation result.

S7: and comparing the reservoir flood regulating result data with the flood regulating simulation result, and performing model parameter calibration of a reservoir flood regulating model.

S8: and taking a warehouse-in flow process line brought by upstream dam break flood as an upstream boundary condition, and performing flood control calculation by using the reservoir flood control model to obtain a reservoir total warehouse-out flow process line.

S9: and (3) performing two-dimensional continuous dam break hydrodynamic simulation calculation based on the total reservoir outlet flow process line of the reservoir in the step S8.

The invention obtains a specific embodiment of the area of the reservoir:

And S1, establishing a reservoir polygon according to the peripheral outline of the reservoir, and calculating the area of the reservoir polygon. The reservoir region polygon is specifically a two-dimensional array of point coordinates sequentially arranged in a clockwise or anticlockwise direction, and the expression is as follows:

，

wherein%) Is->Coordinates of individual points, +.>，/>Is the number of points. The calculation of the area of the polygon in the reservoir area adopts a Gaussian area formula, namely:

，

wherein, the liquid crystal display device comprises a liquid crystal display device,，/>。

the area of the polygon of the reservoir area is taken as the reservoir area under the normal water storage level。

The present invention generates a specific embodiment of a water level array:

s2, according to the water level resolution and the elevation range of the generalized topography, namely the water level differences of two adjacent contour lines in the generalized topographyMinimum elevation value of generalized topography +.>And maximum elevation value->. To->As the first element of the water level array, +.>As an increment, +.>As the last element, a water level array is generated, and the expression is as follows:

，

wherein the method comprises the steps ofIs the number of elements of the water level array.

The present invention forms one embodiment of a water level and area mapping table:

s3, a method for forming a water level and area mapping table comprises the following steps:

s31, judging the normal water storage levelAnd water level array->Is used for the relative position of the two parts.

If it isIs positioned at->Inside, then->For dividing the dots, will ∈ ->Dividing into two new arrays->And->The expression is as follows:

，

wherein, the liquid crystal display device comprises a liquid crystal display device,is composed of all->Middle greater than or equal to->An array of elements, < >>Is the number of these elements; />Is composed of all->Less than or equal to%>An array of elements, < >>Is the number of these elements.

If it isIs greater than->Then a new array is generated:

；

if it isLess than->Then a new array is generated:

。

the method comprises the following stepsLess than->The case of (2) is described as an example, and the processing method of the other cases is the same.

S32: and calculating the storage capacity corresponding to each water level in the water level array, wherein the function expression of the water level and storage capacity curve is as follows:

，

wherein, the liquid crystal display device comprises a liquid crystal display device,is a reservoir capacity>Is a function of the water level and reservoir capacity curve;

and further calculating the reservoir capacity corresponding to each water level in the water level array, wherein the calculation formula is as follows:

，

wherein the method comprises the steps ofIs corresponding to the water level array->The element of the reservoir capacity array is the reservoir capacity corresponding to each water level in the water level array.

If the curve is discrete point data, linear interpolation is carried out on the original data, so that reservoir capacity corresponding to each water level is obtained.

S33: combining the area of the reservoir area under the normal water storage level obtained in the step S1 To->And traversing the new water level group as a starting point, and sequentially calculating the area of the reservoir area corresponding to each water level by adopting a prismatic table volume formula. The volume formula of the prismatic table is as follows:

，

wherein, the liquid crystal display device comprises a liquid crystal display device,is the reservoir capacity between two adjacent water levels, < >>And->Respectively the higher water level and the lower water level of two adjacent water levels, +.>And->Is respectively corresponding to->And->Is a pool area of the container. To->The expression of the array is as follows, which is the starting point traversal array:

。

in the first step of the loop, the expression to the right of the equal sign is as follows:

，/>，/>，

the expression to the left of the equal sign is as follows:

，

further, the above volume formula becomes:

，

wherein the unknown valueIs->Corresponding reservoir area. The above equation is about->Is a unitary quadratic equation of (a). Two roots can be obtained according to the solution of the unitary quadratic equation, and the two roots are squared and then added with +.>A comparison is made. Because the higher the water level in the reservoir, the larger the water area in the corresponding reservoir, so it is less than +.>Is->Discarding the error solution, and obtaining the right solution>. And then sequentially solving the areas of the reservoir areas of the residual water levels to form a water level and area mapping table.

The invention obtains a concrete embodiment of reservoir generalized topography:

s4, a method for obtaining reservoir generalized topography comprises the following steps:

S41: traversing the water level array, and solving the offset of the contour polygon corresponding to each water level relative to the reservoir polygon. Area of polygon after offsetIs the original polygon area->And enlarge distance->The expression of which is as follows:

。

function ofIs of unknown form, but each +.>The corresponding function value can be obtained by an offset curve algorithm and a Gaussian area formula. Traversing water level array->Iterative solution about +.>Nonlinear equation>The offset of each water level is obtained.

S42: according to the previous step, each water levelOffset below->Shifting the polygon of the reservoir area by using a shifting curve algorithm to obtain the water level +.>Is a contour of (c).

S43: and extracting the vertex coordinates of each contour polygon, and taking the water level value corresponding to the polygon as the elevation value of the vertices, thereby constructing the three-dimensional terrain elevation data of the reservoir.

The invention relates to a concrete embodiment for updating reservoir generalized topography, which comprises the following steps:

s5, the method for updating the reservoir generalized topography comprises the following steps:

s51: and establishing a reservoir simulation model of the irregular triangular grid. The three-dimensional terrain elevation data obtained in S4 is mapped onto grid nodes using bilinear interpolation.

S52: the model reservoir capacity under the characteristic water level is rechecked for the generalized topography of the reservoir so as to update the generalized topography of the reservoir, and the model reservoir capacity rechecked specifically comprises the following contents:

calculating the water level of the model as the characteristic water levelWater quantity at the time->。

If the water quantity is the model water quantityStorage capacity corresponding to the water level in S32 +.>If the error of the database is smaller than the given threshold value, the inventory rechecking is completed;

if the water quantity is the model water quantityStorage capacity corresponding to the water level in S32 +.>Is greater than a given thresholdAnd according to the model water quantity and the terrain elevation adjustment quantity below the water level, the prismatic volume formula is used, and the calculation formula is as follows:

，/>

wherein, the liquid crystal display device comprises a liquid crystal display device,is the characteristic water level->Lower reservoir area. Thereby generating new three-dimensional terrain elevation data.

And S6, inputting reservoir design data, checking flood process information and initial roughness, and carrying out flood regulation calculation according to a reservoir dispatching principle. The method comprises the following specific steps:

s61: and (3) setting initial roughness based on the surface type and the condition of the reservoir, and taking reservoir design information and a time sequence process line for checking the flow of the reservoir as upstream boundary conditions.

S62: and setting up a downstream boundary condition of the model according to a reservoir dispatching principle, namely determining a delivery condition according to the reservoir delivery flow condition and the dam front water level condition.

if the delivery condition is that the delivery flow is calculated according to the discharge curve of the reservoir flood discharge building, a water level-flow boundary is established.

S63: and performing simulation calculation on the model to obtain a dam front water level time sequence process line and a reservoir drainage flow process.

The invention relates to a specific embodiment for model parameter calibration of a reservoir flood control model, which comprises the following steps:

s7, a model parameter calibration method of a reservoir flood control model is carried out, and the method comprises the following steps:

based on reservoir flood regulating result data, calculating model calculation result and flood regulating formationMaximum error of dam front water level, peak time error and maximum drain flow error in the fruit data. Because the three are different physical units, normalization processing is carried out, and the square sum of relative errors is taken as an objective functionPerforming simulation model calibration, namely:

，

wherein the method comprises the steps ofIs the true highest dam front water level +.>For the real peak present time, < >>Is the true maximum leakage flow; />Maximum dam front water level calculated for simulation, < ->Peak present time calculated for simulation, +. >The maximum leakage flow calculated for the simulation. Manning roughness coefficient of model as rating parameter +.>. I.e. < ->Iterative solution equationObtain->As the optimal roughness coefficient of the simulation model.

The invention obtains a specific embodiment of a total reservoir outlet flow process line:

s8, a method for obtaining a total reservoir outlet flow process line of the reservoir comprises the following steps:

and taking a warehousing flow process line brought by upstream dam break flood as an upstream boundary condition, and carrying out flood regulating calculation by using the simulation model. Monitoring the water level in front of the dam, and if the overflow condition occurs, calculating the overflow flow by using a wide top weir formula:

，

wherein the method comprises the steps ofTo take into account the coefficient of the influence of the downstream water level on the leakage flow,/->Is of side shrinkage factor, < >>For the flow coefficient of the weir->For the water clear width of the weir crest, < > a>Acceleration of gravity, ++>Is the total head in front of the weir including the flow rate head. At this time, the total outlet flow of the reservoir is equal to the sum of the flood peak flow and the outlet boundary flow. If continuous burst occurs, the flood regulating operation is ended.

The invention relates to a specific embodiment for carrying out two-dimensional continuous dam break hydrodynamic simulation calculation:

s9, a two-dimensional continuous dam break hydrodynamic force simulation calculation method comprises the following steps:

And (3) performing two-dimensional continuous dam break hydrodynamic simulation calculation based on the total outlet flow process of each reservoir in the step (S8). When the reservoir is not burst, calculating dam site flow by using the corresponding reservoir outlet flow process calculated by the flood regulating model; when the reservoir reaches the breaking condition, the dam site flow is calculated in real time according to the local dam front water level, the dam bottom water level and the topography.

The invention relates to a dam break flood simulation system based on reservoir flood control coupling hydrodynamic force, which comprises the following steps:

a dam break flood simulation system based on reservoir flood regulation coupling hydrodynamic force comprises a terrain generalization module, a reservoir simulation module, a reservoir flood regulation module, a reservoir parameter calibration module and a dam break coupling simulation module.

And the dam break coupling simulation module is used for carrying out two-dimensional continuous dam break hydrodynamic simulation calculation according to the total reservoir outlet flow process line of the reservoir to obtain a dam break simulation result and realize dam break flood simulation based on reservoir flood control coupling hydrodynamic force.

The invention relates to a second concrete embodiment of a dam break flood simulation system based on reservoir flood control coupling hydrodynamic force:

a dam-break flood simulation system based on reservoir flood regulation coupling hydrodynamic force comprises a reservoir topography generalization module, a reservoir flood regulation module and a two-dimensional continuous dam-break module.

The reservoir topography generalization module draws a reservoir polygon based on a satellite map of a geographic information system, calculates the reservoir area under a normal water level, generates a water level array according to reservoir engineering characteristic data, calculates the reservoir area corresponding to each water level, performs offset processing on the reservoir polygon, generates generalized topography contour lines, and further obtains three-dimensional reservoir generalized topography data.

The reservoir flood regulating module is based on a reservoir flood regulating model and reservoir flood regulating result data generated in the reservoir topography generalization module, a dispatching function capable of reflecting a reservoir flood regulating principle is added in the model, the flood regulating model is calibrated, and the optimal roughness coefficient of the reservoir flood regulating model is determined.

The two-dimensional continuous dam break module calculates the dam front water level process of the downstream reservoir based on the total reservoir outlet flow calculated by the reservoir flood regulating module, and synchronously judges the dam break occurrence condition. If continuous burst occurs, dam site flow is calculated in real time according to the local dam front water level, the dam lower water level and the topography. And obtaining the influence results of dam-break flood on the downstream flooding area, wherein the influence results comprise a maximum submerged deep field, a maximum flow velocity field and a submerged duration field.

As shown in fig. 3, a specific embodiment of the dam break flood simulation method of the present invention is applied:

and drawing a reservoir area polygon of a certain reservoir based on a satellite map of a geographic information system, calculating the reservoir area under a normal water level, generating a water level array according to reservoir engineering characteristic data, calculating the reservoir area corresponding to each water level, performing offset processing on the reservoir area polygon, generating generalized topography contour lines, and further obtaining three-dimensional reservoir generalized topography data. The specific flow is as follows:

1. and collecting reservoir engineering characteristic data, including a water level reservoir capacity curve, a reservoir scheduling principle and a reservoir engineering characteristic table.

2. And drawing a peripheral outline of a reservoir area on a satellite map of the geographic information system, and taking the outline as a reservoir area polygon under a normal water storage level. The area of the polygon in the reservoir area is calculated according to a Gaussian area formula, and the specific calculation formula is as follows:

，

3. Generating a water level array, and generalizing the elevation range of the terrain to be from the elevation of a river channel under a reservoir damTo dam crest elevation +.>The water level resolution is 2m. Normal water level->Within this range, the expression of the water level array is thus as follows: />

，

。

And two reservoir characteristic water levels are added into the water level array to respectively design flood levels for reservoirsAnd check flood level->。

4. Based on the reservoir level-reservoir capacity table and the discrete point data, the method is based on linear interpolationAnd->Corresponding to the water levels. And solving a unitary quadratic equation by combining the prismatic table volume formula to obtain a water level and area mapping table.

5. Traversing the water level array, iteratively calculating the offset corresponding to each water level in the array, and offsetting the reservoir area polygon under the normal water storage level in the step 1 by an offset curve algorithm to obtain the contour line of each water level, namely the reservoir generalized topography. And extracting vertex coordinates of each contour line, taking the corresponding water level as the terrain elevation value of the points to obtain three-dimensional terrain elevation data, wherein the expression is as follows:

，

wherein the method comprises the steps of、/>Are respectively->The (th) on the personal contour line>The abscissa of the individual points, +. >Is->Elevation corresponding to the contour lines, +.>，/>，/>Is->The number of vertices on the contour lines.

6. An irregular triangular mesh with a mesh scale of 30m is established, which comprises the following contents:

in this embodiment, an irregular triangular mesh of a downstream reservoir is established, see fig. 4, and three-dimensional terrain elevation data is mapped onto mesh nodes by using bilinear interpolation to construct a reservoir simulation model. Respectively calculating the model water level as normal water storage levelDesign flood level->And check flood level->And the water quantity is compared with the actual reservoir capacity in the water level and reservoir capacity curve, and the error is determined. Judging whether the error is higher than a threshold (5%), and when the judgment result is yes, calculating elevation adjustment according to a prism volume formula to generate new three-dimensional terrain elevation data:

，

and when the judgment result is no, namely the error is not higher than the threshold value (5%), adding a scheduling function capable of reflecting the reservoir flood scheduling principle into the model based on the reservoir simulation model and reservoir flood scheduling result data, calibrating the simulation model, and determining the optimal roughness coefficient. The specific flow is as follows:

1. the flood flow rate is designed and checked in 72 hours (h) of the reservoir as an upstream boundary condition, and the initial roughness is selected to be 0.030 of the normal roughness value of a mountain large river (no grass tree in the river channel, steeper river bank and submerged when a bank slope is used for tree cluster flood).

The expression of the reservoir 72h design and check flood flow process is as follows:

，

wherein the method comprises the steps ofTime (unit hour) is +.>Is->Time to->Average flow in each moment, +.>。

Determining a downstream boundary condition based on reservoir flood regulation principles:

when the storage flow is smaller than the maximum lower flow of the reservoir flood discharge building under the normal water storage levelWhen the warehouse-out flow is equal to the warehouse-in flow, the downstream boundary condition is set as the fixed water level boundary.

When the storage flow is greater than or equal to the maximum lower discharge flow of the reservoir flood discharge building under the normal water storage levelIn the time of this, according to the flood discharge curve of the reservoir flood discharge building +.>Calculating the reservoir outlet flow, and setting the downstream boundary condition as the water level flow boundary.

The simulation model is calculated, and the maximum value of the dam front water level, the peak current time and the maximum reservoir drainage flow are output, in this embodiment, the flood control calculation result of the downstream reservoir model for the typical flood process can be seen in fig. 5.

2. Comparing the reservoir flood regulating result data, normalizing the three calculated values in the step 7, taking the square sum of relative errors, and establishingRated objective functionAnd (3) iteratively solving the minimum value to obtain the optimal roughness coefficient of the simulation model, wherein the calculation formula is as follows:

，

3. And taking a warehousing flow process line brought by upstream dam break flood as an upstream boundary condition, and carrying out flood regulating calculation by using the simulation model. Monitoring the water level in front of the dam, judging whether a top overflow condition occurs, and when the judgment result is yes, calculating the top overflow flow by using a wide top weir formula:

，

and when the judgment result is no, the overtaking flow is zero, and the total ex-warehouse flow process is output.

In this embodiment, the calculation result of the downstream water reservoir model for flood control of a dam break flood can be seen in fig. 6, which specifically includes the following contents:

before dam break flood arrives, the water level in front of the dam is calculated to be kept at the normal water storage level of 52.50 meters (m), and the discharge flow under the reservoir is kept at 2000 cubic meters per second (m) ³ S), no overtopping condition occurs, so that the overtopping flow is zero; when dam break flood reaches a downstream water reservoir, calculating the rapid rising of the water level in front of the dam, calculating the discharge flow under the reservoir according to a flood discharge curve, and after the water level in front of the dam reaches the 54.90 meters (m) of the dam top, starting to rapidly rise the flow on top of the dam; after a period of time (h), along with the passing of a dam break flood peak, calculating the front water level of the dam to start falling, returning the top flood flow to zero after calculating the front water level to return to the dam top elevation, and calculating the reservoir bottom discharge flow according to the reservoir storage flow after calculating the front water level to return to the normal water storage level, so as to realize the reservoir storage flow equal to the reservoir discharge flow and ensure that the water level is stabilized on the normal water storage level.

And calculating the dam front water level process of the downstream reservoir based on the total reservoir outlet flow, performing two-dimensional continuous dam break calculation, and synchronously judging the dam break occurrence condition. If continuous burst occurs, dam site flow is calculated in real time according to the local dam front water level, the dam lower water level and the topography. And analyzing the influence on the downstream flooding area to obtain an influence result of dam-break flood on the downstream flooding area, wherein the influence result comprises a maximum submerged deep field, a maximum flow velocity field and a submerged duration field.

Therefore, the invention can realize generalization of reservoir topography and establishment and calibration of reservoir flood control models, and is suitable for reservoir flood control two-dimensional simulation calculation of continuous dam break analysis.

Furthermore, the invention is based on the coupling of reservoir flood control model and continuous dam break hydrodynamic simulation, can realize the flood control simulation of the reservoir before dam break under the continuous dam break working condition, and has detailed and practical scheme.

An embodiment of a device for applying the method of the invention:

an electronic device, comprising:

one or more processors;

a storage means for storing one or more programs;

and when the one or more programs are executed by the one or more processors, the one or more processors are caused to implement a dam break flood simulation method based on reservoir flood control coupling hydrodynamic force.

A computer medium embodiment to which the method of the application is applied:

a computer readable storage medium having stored thereon a computer program which when executed by a processor implements a dam break flood simulation method based on reservoir flood control coupled hydrodynamic forces as described above.

It will be appreciated by those skilled in the art that embodiments of the present application may be provided as methods, systems, computer program products. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present application is described in terms of methods, apparatus (systems), computer program products, flowcharts, and/or block diagrams in accordance with embodiments of the present application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

Finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.

Claims

1. A dam break flood simulation method based on reservoir flood regulation coupling hydrodynamic force is characterized by comprising the following steps of:

the method comprises the following steps:

acquiring a certain reservoir parameter to be dam-break simulated;

2. The dam break flood simulation method based on reservoir flood control coupling hydrodynamic force as claimed in claim 1, wherein the dam break flood simulation method is characterized by comprising the following steps:

3. The dam break flood simulation method based on reservoir flood control coupling hydrodynamic force as claimed in claim 1, wherein the dam break flood simulation method is characterized by comprising the following steps:

4. A dam break flood simulation method based on reservoir flood control coupling hydrodynamic force as claimed in claim 3, wherein the dam break flood simulation method is characterized in that:

5. The dam break flood simulation method based on reservoir flood control coupling hydrodynamic force as claimed in claim 4, wherein the dam break flood simulation method is characterized in that:

the reservoir dispatching principle comprises the following contents:

6. The dam break flood simulation method based on reservoir flood control coupling hydrodynamic force according to claim 5, wherein the dam break flood simulation method is characterized in that:

7. The dam break flood simulation method based on reservoir flood control coupling hydrodynamic force as claimed in claim 6, wherein the dam break flood simulation method is characterized in that:

if continuous burst occurs, ending the flood regulating operation;

8. The dam break flood simulation method based on reservoir flood control coupling hydrodynamic force according to any one of claims 1-7, wherein the dam break flood simulation method is characterized in that:

9. A dam break flood simulation method based on reservoir flood regulation coupling hydrodynamic force is characterized by comprising the following steps of:

the method comprises the following steps:

firstly, acquiring a reservoir parameter to be dam-break simulated;

10. Dam break flood simulation system based on reservoir flood control coupling hydrodynamic force, characterized in that:

a dam break flood simulation method based on reservoir flood control coupling hydrodynamic force according to any one of claims 1-9, comprising a terrain generalization module, a reservoir simulation module, a reservoir flood control module, a reservoir parameter calibration module and a dam break coupling simulation module: