A calculation method and device for hydraulic parameters of closure by end dumping

Abstract

The present invention provides a method and device for calculating the hydraulic parameters of closing a gap by vertical damming, which relates to the technical field of water conservancy projects. The method includes: determining the bottom width of the gap, and determining the assumed upstream water depth and the assumed downstream water depth of the gap according to the bottom width of the gap; obtaining the diversion flow rate of the diversion structure according to the assumed upstream water depth of the gap; constructing a first constraint condition according to the designed closure flow rate, the known upstream river channel storage flow rate and the diversion flow rate of the diversion structure; determining the discharge of the gap according to the assumed upstream water depth of the gap, the assumed downstream water depth of the gap and the obtained dike parameters; constructing a second constraint condition according to the dike parameters and the discharge of the gap; performing nonlinear programming solution according to the first constraint condition and the second constraint condition; when the solution is successful, obtaining the hydraulic parameters according to the bottom width of the gap, the upstream water depth of the gap and the downstream water depth of the gap, which are used for accurately and automatically calculating the hydraulic parameters of closing the gap by the vertical damming method.

Description

Technical Field

[0001] The present invention relates to the technical field of water conservancy engineering, and in particular to a method and device for calculating hydraulic parameters of a vertical cut-off flow. Background Art

[0002] Hydraulic calculation of diversion is an important part of the construction organization of vertical blocking and diversion. Through hydraulic calculation, the hydraulic parameters upstream of Longkou are obtained, providing a basis for diversion design.

[0003] In related technologies, interception hydraulic calculations are mainly based on the principle of flow balance and are calculated through graphical methods, manual trial calculations, and simplified methods. Among them, the graphical method draws a water level flow relationship curve on a graph and solves the water depth upstream of Longkou by finding the intersection of the curves. This method has low accuracy and is difficult to meet the requirements of modern engineering. Manual trial calculations require repeated manual trial calculations by constantly assuming independent variables, which takes a lot of time. Due to the high complexity of the calculations, the probability of obtaining an accurate solution is low. The simplified method takes simplified measures, ignores changing parameters such as dike seepage and submergence coefficient, or approximates the average water-passing width of the Longkou axis section to a constant value, resulting in certain errors in specific situations.

[0004] In this regard, there is no good method in the relevant technology to build an automated calculation system for the hydraulic parameters of the vertical plugging method. Summary of the Invention

[0005] The problem solved by the present invention is how to accurately and automatically calculate the hydraulic parameters of the vertical plugging method.

[0006] In order to solve the above problems, the present invention provides a method and device for calculating the hydraulic parameters of a vertical cut-off flow.

[0007] In a first aspect, the present invention provides a method for calculating hydraulic parameters of a vertical cut-off flow, comprising:

[0008] Determine the bottom width of the dragon mouth, and determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth;

[0009] According to the assumed water depth upstream of the dragon mouth, piecewise linear interpolation is performed in the first water level-discharge relationship curve of the diversion building to obtain the diversion flow of the diversion building;

[0010] According to the known interception design flow and the known upstream river channel storage flow, a first constraint condition is established based on the diversion flow of the diversion building and the second water level flow relationship curve downstream of the dragon mouth;

[0011] Determining the discharge of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth and the acquired dike parameters;

[0012] Constructing a second constraint condition based on a numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth, and the assumed dragon mouth downstream water depth;

[0013] Performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition;

[0014] When the solution is successful, the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth are used as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth;

[0015] Hydraulic parameters are obtained according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth, and the water depth downstream of the dragon mouth.

[0016] Optionally, the method for calculating hydraulic parameters of vertical cut-off flow further includes:

[0017] Determine different dragon mouth bottom widths, determine the corresponding assumed upstream dragon mouth water depth, assumed downstream dragon mouth water depth, first constraint condition and second constraint condition according to each dragon mouth bottom width, and perform nonlinear programming solution;

[0018] When the solution is successful, the corresponding dragon mouth bottom width, water depth upstream of the dragon mouth and water depth downstream of the dragon mouth are calculated to obtain the corresponding hydraulic parameters under each dragon mouth bottom width.

[0019] Optionally, obtaining the diversion flow of the diversion building by piecewise linear interpolation in the first water level-flow relationship curve of the diversion building according to the assumed water depth upstream of the dragon mouth includes:

[0020] Obtaining a water level interval including the assumed water depth upstream of the dragon mouth in the first water level-flow relationship curve;

[0021] Obtaining a flow rate interval corresponding to the water level interval in the first water level-flow rate relationship curve;

[0022] Linear interpolation is performed on an array consisting of the water level interval and the flow interval to obtain the diversion flow of the diversion building.

[0023] Optionally, the constructing of the first constraint condition based on the known interception design flow and the known upstream river channel storage flow, and based on the diversion flow of the diversion building and the second water level-flow relationship curve downstream of the dragon mouth includes:

[0024] Obtaining the downstream flow of Longkou according to the interception design flow, the upstream river channel storage flow and the diversion flow of the diversion building;

[0025] According to the flow rate downstream of the dragon mouth, in the second water level-flow relationship curve, the constrained water depth downstream of the dragon mouth is obtained by piecewise linear interpolation;

[0026] The first constraint condition is constructed according to the first constrained water depth downstream of the dragon mouth, wherein the first constraint condition includes: the first constrained water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth.

[0027] Optionally, the dike parameters include dike upstream slope ratio, dike downstream slope ratio, dike height, dike top width, dike entry length and permeability coefficient;

[0028] Determining the discharge volume of the Longkou according to the assumed water depth upstream of the Longkou, the assumed water depth downstream of the Longkou, and the acquired dike parameters includes:

[0029] Determining the length of the infiltration path according to the upstream slope ratio of the dike, the downstream slope ratio of the dike, the dike height, and the dike top width;

[0030] Determining the dike seepage rate based on the infiltration path length, the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the dike occupation length, and the permeability coefficient;

[0031] The Longkou discharge volume is determined based on the dike seepage volume, the diversion design flow volume, the upstream river channel storage volume and the diversion volume of the diversion building.

[0032] Optionally, the dike parameters include dike head slope ratio;

[0033] The second constraint condition is constructed according to the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth, and the assumed dragon mouth downstream water depth, including:

[0034] When a first ratio between the assumed water depth downstream of the dragon mouth and the assumed water depth upstream of the dragon mouth is greater than or equal to a preset ratio, obtaining a submergence coefficient by piecewise linear interpolation in the Pavlovsky submergence coefficient table according to the first ratio;

[0035] Determine the average water-passing width of the Longkou axis section according to the submergence coefficient, the known discharge coefficient, the assumed water depth upstream of the Longkou, and the Longkou discharge;

[0036] Determine the second constrained water depth downstream of the dragon mouth according to the average water-passing width of the dragon mouth axis section, the dragon mouth bottom width and the dike head slope ratio;

[0037] The second constraint condition is constructed according to the second constrained water depth downstream of the dragon mouth, wherein the second constraint condition includes: the second constrained water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth.

[0038] Optionally, constructing the second constraint condition based on the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed water depth upstream of the dragon mouth, and the assumed water depth downstream of the dragon mouth further includes:

[0039] When a first ratio between the assumed water depth downstream of the dragon mouth and the assumed water depth upstream of the dragon mouth is less than a preset ratio, the average water-passing width of the dragon mouth axis section is determined according to the first ratio using a non-submerged discharge formula;

[0040] Determine the third constrained water depth downstream of the dragon mouth according to the average water width of the dragon mouth axis section;

[0041] When the third constrained water depth downstream of the dragon mouth is greater than the assumed water depth downstream of the dragon mouth, the second constraint condition includes that the third constrained water depth downstream of the dragon mouth is equal to the critical flow water depth;

[0042] When the third constrained water depth downstream of the dragon mouth is less than or equal to the assumed water depth downstream of the dragon mouth, the second constraint condition includes that the third constrained water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth.

[0043] Optionally, after performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition, the method further includes:

[0044] When the solution fails, assuming a new water depth upstream of the dragon mouth and a new water depth downstream of the dragon mouth under the determined dragon mouth bottom width;

[0045] Perform nonlinear programming solution on the new assumed water depth upstream of Longkou and the new assumed water depth downstream of Longkou until a successful solution is obtained.

[0046] Optionally, the hydraulic parameters include head height, water depth, single-width flow rate and flow velocity.

[0047] In a second aspect, the present invention provides a device for calculating hydraulic parameters of a vertical cut-off flow, comprising:

[0048] An assumption module is used to determine the bottom width of the dragon mouth, and to determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth;

[0049] A piecewise linear interpolation module is used to obtain the diversion flow of the diversion building by piecewise linear interpolation in the first water level-flow relationship curve of the diversion building according to the assumed water depth upstream of the dragon mouth;

[0050] A first constraint module is used to construct a first constraint condition based on the known interception design flow and the known upstream river channel storage flow, and based on the diversion flow of the diversion building and the second water level flow relationship curve downstream of the dragon mouth;

[0051] a discharge flow determination module, configured to determine the discharge flow of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, and the acquired dike parameters;

[0052] A second constraint module is used to construct a second constraint condition according to the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth and the assumed dragon mouth downstream water depth;

[0053] a nonlinear programming solution module, configured to perform nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition;

[0054] An association module is configured to, when the solution is successful, use the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth;

[0055] A solution module is used to obtain hydraulic parameters according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth.

[0056] The beneficial effects of the method for calculating hydraulic parameters of vertical cut-off flow of the present invention are:

[0057] Given a dragon mouth bottom width, and based on this width, set the initial dragon mouth upstream and downstream water depths. By fixing the dragon mouth bottom width and setting the initial assumed dragon mouth upstream and downstream water depths, basic parameters are provided for subsequent calculations. The piecewise linear interpolation method is used to obtain the diversion flow of the diversion structure that matches the assumed dragon mouth upstream water depth from the existing first water level flow relationship curve, which improves the calculation efficiency and ensures the rationality of the interpolation result. The first and second constraints constructed according to the interception design flow, upstream river channel storage flow, diversion flow of the diversion structure, dike parameters and the dragon mouth discharge flow ensure that the dragon mouth bottom width, dragon mouth upstream water depth, dragon mouth downstream flow and dragon mouth downstream water depth conform to natural laws. Through nonlinear programming solution, the optimal solution that meets all constraints is automatically found, which realizes the efficient solution of complex hydraulic problems, avoids the uncertainty caused by manual adjustment of parameters, simplifies the calculation process, and can quickly generate a complete hydraulic parameter set under a given dragon mouth bottom width, thereby accurately and automatically calculating the hydraulic parameters of the vertical plugging method interception. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 Schematic diagram of the flow chart of the method for calculating hydraulic parameters of a vertical cut-off flow according to an embodiment of the present invention;

[0059] Figure 2 This is a flow chart of a refinement of step S200 of the method for calculating hydraulic parameters of a vertical cut-off flow according to an embodiment of the present invention;

[0060] Figure 3 1 is a flowchart of a method for calculating hydraulic parameters of a vertical cut-off flow according to an embodiment of the present invention;

[0061] Figure 4is an exemplary diagram of an electronic device according to an embodiment of the present invention;

[0062] Figure 5 This is an elevation view of the upstream of Longkou according to an embodiment of the present invention;

[0063] Figure 6 It is a cross-sectional view of the dragon mouth of an embodiment of the present invention. DETAILED DESCRIPTION

[0064] To make the above-mentioned objects, features, and advantages of the present invention more clearly understood, specific embodiments of the present invention are described in detail below with reference to the accompanying drawings. Although certain embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be construed as being limited to the embodiments described herein. Instead, these embodiments are provided to provide a more thorough and complete understanding of the present invention. It should be understood that the drawings and embodiments of the present invention are for illustrative purposes only and are not intended to limit the scope of protection of the present invention.

[0065] It should be understood that the various steps described in the method embodiments of the present invention may be performed in different orders and / or in parallel. In addition, the method embodiments may include additional steps and / or omit the steps shown. The scope of the present invention is not limited in this respect.

[0066] The term "including" and its variations used in this document are open inclusions, that is, "including but not limited to"; the term "based on" means "based at least in part on"; the term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one other embodiment"; the term "some embodiments" means "at least some embodiments"; the term "optionally" means "optional embodiments". The relevant definitions of other terms will be given in the following description. It should be noted that the concepts of "first", "second", etc. mentioned in the present invention are only used to distinguish different devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these devices, modules or units.

[0067] It should be noted that the modifications of "one" and "multiple" mentioned in the present invention are illustrative rather than restrictive. Those skilled in the art should understand that unless otherwise clearly indicated in the context, it should be understood as "one or more".

[0068] The names of the messages or information exchanged between multiple devices in the embodiments of the present invention are only used for illustrative purposes and are not used to limit the scope of these messages or information.

[0069] In response to the problems existing in the above-mentioned related technologies, this embodiment provides a method and device for calculating the hydraulic parameters of a vertical shut-off flow.

[0070] likeFigure 1 As shown, an embodiment of the present invention provides a method for calculating hydraulic parameters of a vertical cut-off flow, comprising:

[0071] Step S100: determining the bottom width of the dragon mouth, and determining an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth.

[0072] Before calculations begin, a preliminary estimate or assumption is made for the upstream and downstream water depths based on the determined dam bottom width. These assumed values form the basis for subsequent calculations, used to determine parameters such as the diversion flow rate of the diversion structure and the dam discharge. These values are gradually adjusted throughout the calculation process using nonlinear programming until all constraints are met.

[0073] Specifically, it is assumed that the water depth H upstream of Longkou refers to the height of the water surface upstream of Longkou relative to a certain reference surface; and it is assumed that the water depth h downstream of Longkou refers to the height of the water surface upstream of Longkou relative to a certain reference surface. s Refers to the height of the water surface downstream of the dragon mouth relative to the same reference plane. The two constitute important initial conditions in the calculation process. In actual operation, it may be necessary to try different H and h multiple times. s to find a solution that suits the actual situation.

[0074] Step S200: According to the assumed water depth upstream of the dragon mouth, piecewise linear interpolation is performed in the first water level-flow relationship curve of the diversion building to obtain the diversion flow of the diversion building.

[0075] In one embodiment, diversion structures include structures such as diversion tunnels and open diversion channels used to direct or change the direction of water flow and distribute flow in water conservancy and hydropower projects. The first water level-flow relationship curve is used to characterize the correspondence between different water levels and flow rates under given conditions. The flow rate diverted by the diversion structure represents the flow rate of the portion of water flow directed through the diversion structure. The first water level-flow relationship curve serves as one of the known boundary conditions.

[0076] By combining several functions in Excel, piecewise linear interpolation is performed on the first water level-discharge relationship curve to obtain the diversion flow of the diversion building under the assumed upstream water depth, which provides the necessary input for the subsequent nonlinear programming solution.

[0077] Step S300: constructing a first constraint condition based on the known interception design flow and the known upstream river channel storage flow, and based on the diversion flow of the diversion building and the second water level flow relationship curve downstream of the dragon mouth.

[0078] The design flow rate of the closure Q represents the maximum flow rate allowed through the dragon mouth during the closure construction period. r It represents the amount of water added in the upstream river channel due to the rising water level. The second water level-discharge relationship curve is one of the known boundary conditions.

[0079] Based on the relationship curve between the interception design flow, the upstream channel storage flow, the diversion flow from the diversion structure, and the second water level flow curve downstream of Longkou, the constrained water depth downstream of Longkou was obtained. This was then combined with the initially assumed water depth downstream of Longkou to form the first constraint. This first constraint ensured that the relationship between the flow downstream of Longkou and its corresponding water depth conformed to natural laws, thereby verifying the correctness of the calculation model.

[0080] Step S400: determining the discharge rate of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth and the acquired dike parameters.

[0081] Since the parameters of the dike affect the way and amount of water flowing out of the dragon mouth, the obtained dike parameters can be used together with the assumed water depth upstream and downstream of the dragon mouth to determine the dragon mouth discharge under assumed conditions, thereby affecting the design size, material selection and construction method of the intercepting structure (such as the dike).

[0082] Step S500: constructing a second constraint condition based on the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth, and the assumed dragon mouth downstream water depth.

[0083] The discharge type of Longkou discharge can be determined based on the numerical relationship between the assumed water depth upstream of Longkou and the assumed water depth downstream of Longkou. According to the discharge type, the constrained water depth downstream of Longkou under given conditions can be solved according to the dike parameters and the Longkou discharge volume. Then, a second constraint condition is constructed based on the constrained water depth downstream of Longkou. According to the second constraint condition, it is ensured that the calculated constrained water depth downstream of Longkou is consistent with the assumed water depths upstream and downstream of Longkou and conforms to specific physical laws, thereby ensuring the validity of the assumption and improving the reliability of the calculation results.

[0084] Step S600: performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition.

[0085] In one embodiment, a nonlinear programming problem is solved by calling Excel's built-in solver. The nonlinear programming problem is solved using the assumed water depths upstream and downstream of the Longkou as independent variables, and the first and second constraints as two constraint equations. The goal is to find the optimal solution that satisfies all given constraints, that is, the water depths upstream and downstream of the Longkou that accurately reflect the actual project conditions.

[0086] Step S700: When the solution is successful, the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth are used as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth.

[0087] In one embodiment, ifFigure 5 As shown, when the solution is successful, it means that the assumed water depth upstream of Longkou and the assumed water depth downstream of Longkou satisfy both the first constraint and the second constraint. Then, the assumed water depth upstream of Longkou and the assumed water depth downstream of Longkou are used as the water depth upstream of Longkou and the water depth downstream of Longkou corresponding to the Longkou bottom width determined in step S100. The Longkou bottom width is associated with its corresponding water depth upstream of Longkou and water depth downstream of Longkou and stored. In this process, the dike seepage amount Q corresponding to the water depth upstream of Longkou and water depth downstream of Longkou is also stored. s , Longkou discharge flow Q g and other parameters.

[0088] Step S800, obtaining hydraulic parameters according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth, and the water depth downstream of the dragon mouth.

[0089] The hydraulic parameters under the dragon mouth bottom width are calculated based on the dragon mouth bottom width, the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth.

[0090] Optionally, the hydraulic parameters include head height, water depth, single-width flow rate, flow rate, etc.

[0091] In this embodiment, a dragon mouth bottom width is given, and the initial dragon mouth upstream and downstream water depths are set based on this width. By fixing the dragon mouth bottom width and setting the initial assumed dragon mouth upstream and downstream water depths, basic parameters are provided for subsequent calculations. The piecewise linear interpolation method is used to obtain the diversion flow of the diversion building that matches the assumed dragon mouth upstream water depth from the existing first water level flow relationship curve, which improves the calculation efficiency and ensures the rationality of the interpolation result. The first constraint condition and the second constraint condition constructed according to the interception design flow, the upstream river channel storage flow, the diversion flow of the diversion building, the dike parameters and the dragon mouth discharge flow ensure that the dragon mouth bottom width, the dragon mouth upstream water depth, the dragon mouth downstream flow and the dragon mouth downstream water depth conform to natural laws. Through nonlinear programming solution, the optimal solution that meets all constraints is automatically found, which realizes the efficient solution of complex hydraulic problems, avoids the uncertainty caused by manual adjustment of parameters, simplifies the calculation process, and then can quickly generate a complete hydraulic parameter set under a given dragon mouth bottom width, thereby accurately and automatically calculating the hydraulic parameters of the vertical plugging method interception.

[0092] In one embodiment, the scheme formed by steps S100-S800 of the present invention can automatically assign initial values multiple times to perform nonlinear programming solutions, avoiding the situation where a single evaluation failure requires manual intervention. Fully automated solutions can be performed, and various hydraulic parameters can be iteratively obtained without manual intervention. Graphs can be drawn based on the hydraulic parameters to show a schematic diagram of the changes in various hydraulic parameters during the vertical blockage process.

[0093] Optionally, the method for calculating hydraulic parameters of vertical cut-off flow further includes:

[0094] Determine different dragon mouth bottom widths, determine the corresponding assumed dragon mouth upstream water depth, assumed dragon mouth downstream water depth, first constraint condition and second constraint condition according to each dragon mouth bottom width, and perform nonlinear programming solution.

[0095] When the solution is successful, the corresponding dragon mouth bottom width, water depth upstream of the dragon mouth and water depth downstream of the dragon mouth are calculated to obtain the corresponding hydraulic parameters under each dragon mouth bottom width.

[0096] In one embodiment, steps S100-S800 are used to obtain the upstream and downstream water depths and various hydraulic parameters corresponding to different given dam bottom widths, so as to realize automatic calculation of hydraulic parameters for vertical plugging method.

[0097] Alternatively, as Figure 2 、 Figure 5 and Figure 6 As shown, according to the assumed water depth upstream of Longkou, obtaining the diversion flow of the diversion building by piecewise linear interpolation in the first water level-flow relationship curve of the diversion building includes:

[0098] Step S210: obtaining a water level interval including the assumed water depth upstream of the dragon mouth in the first water level-flow relationship curve.

[0099] Step S220: Obtain the flow rate interval corresponding to the water level interval in the first water level-flow rate relationship curve.

[0100] Step S230: performing linear interpolation on the array formed by the water level interval and the flow interval to obtain the diversion flow of the diversion building.

[0101] In one embodiment, in Excel, the first water level-flow relationship curve of the diversion building is a data table with two adjacent columns, the left column is the water level, and the right column is the corresponding flow. The water level increases from top to bottom in the data table. The piecewise linear interpolation formula is expressed as:

[0102] =TREND(OFFSET(H1:H2,MATCH(H+E,H1:Hn,1)-1,1,2,1),OFFSET(H1:H2,MATCH(H+E,H1:Hn,1)-1,0,2,1),H+E),

[0103] Among them, H represents the assumed water depth upstream of Longkou, E represents the elevation of Longkou bottom, H1 represents the first water level (cell form) in the first water level-flow relationship table, H2 represents the second water level (cell form) in the first water level-flow relationship table, and Hn represents the last water level (cell form).

[0104] The MATCH function in the formula can find the position of the interval to be piecewise linearly interpolated in the data column based on the assumed water depth H+E upstream of Longkou; the OFFSET function can obtain the two adjacent water level cells and flow cells (in array form) of the interpolation interval based on the position obtained by the MATSET function; the TREND function interpolates and solves the diversion flow of the diversion building based on the array obtained by OFFSET. .

[0105] Optionally, the constructing of the first constraint condition based on the known interception design flow and the known upstream river channel storage flow, and based on the diversion flow of the diversion building and the second water level-flow relationship curve downstream of the dragon mouth includes:

[0106] The flow rate downstream of Longkou is obtained according to the diversion design flow rate, the storage flow rate of the upstream river channel and the diversion flow rate of the diversion building.

[0107] According to the flow rate downstream of Longkou, the constrained water depth downstream of Longkou is obtained by piecewise linear interpolation in the second water level-flow relationship curve.

[0108] The first constraint condition is constructed according to the first constrained water depth downstream of the dragon mouth, wherein the first constraint condition includes: the first constrained water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth.

[0109] In one embodiment, the upstream channel storage flow is the fixed value obtained, the flow rate downstream of Longkou can be expressed as According to the known natural water level-discharge relationship curve, the first constraint water depth downstream of Longkou is obtained by piecewise linear interpolation in the second water level-discharge relationship curve. , wherein the piecewise linear interpolation algorithm is the same as the piecewise linear interpolation algorithm in step S200, and the first constraint water depth downstream of Longkou With the assumed water depth h downstream of Longkou s Equality is the first constraint.

[0110] Alternatively, as Figure 5 and Figure 6 As shown, the dike parameters include the upstream slope ratio of the dike, the downstream slope ratio of the dike, the dike height, the dike top width, the dike occupation length and the permeability coefficient.

[0111] Determining the discharge volume of the Longkou according to the assumed water depth upstream of the Longkou, the assumed water depth downstream of the Longkou, and the acquired dike parameters includes:

[0112] The length of the infiltration path is determined according to the upstream slope ratio of the dike, the downstream slope ratio of the dike, the dike height and the dike top width.

[0113] The dike seepage volume is determined according to the length of the seepage path, the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the dike occupation length and the permeability coefficient.

[0114] The Longkou discharge volume is determined based on the dike seepage volume, the diversion design flow volume, the upstream river channel storage volume and the diversion volume of the diversion building.

[0115] The formula for determining the penetration path length is expressed as:

[0116] ,

[0117] The formula for determining the seepage rate of the dike is expressed as:

[0118] ,

[0119] in, Indicates the upstream slope ratio of the dike, Downstream slope ratio of the dike, P Height of dike, a The top of the dike is wide.

[0120] The formula for determining the dragon mouth discharge is expressed as:

[0121] ,

[0122] in, Indicates the discharge flow from the dragon mouth, Q indicates the interception design flow, Indicates the upstream channel storage flow.

[0123] Optionally, the dike parameters include the dike head slope ratio.

[0124] Alternatively, as Figure 3 and Figure 5 As shown, the second constraint condition is constructed according to the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth and the assumed dragon mouth downstream water depth, including:

[0125] When a first ratio between the assumed water depth downstream of the dragon mouth and the assumed water depth upstream of the dragon mouth is greater than or equal to a preset ratio, the submergence coefficient is obtained by piecewise linear interpolation in the Pavlovsky submergence coefficient table according to the first ratio.

[0126] The average water-passing width of the Longkou axial section is determined based on the flooding coefficient, the known flow coefficient, the assumed water depth upstream of the Longkou and the Longkou discharge.

[0127] The second constrained water depth downstream of the dragon mouth is determined according to the average water-passing width of the dragon mouth axis section, the dragon mouth bottom width and the dike head slope ratio.

[0128] The second constraint condition is constructed according to the second constrained water depth downstream of the dragon mouth, wherein the second constraint condition includes: the second constrained water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth.

[0129] When the first ratio of the assumed water depth downstream of Longkou to the assumed water depth upstream of Longkou is greater than or equal to the preset ratio, it indicates that the Longkou discharge is a submerged flow. The submerged coefficient σ is obtained by piecewise linear interpolation in the Pavlovsky submerged coefficient table using the first ratio. The piecewise linear interpolation method is consistent with the piecewise linear interpolation algorithm in step S200. Then, the average water width of the Longkou axis section is solved according to the submerged flow discharge formula. ,Depend on Dragon mouth bottom width B and dike head slope ratio Solve the second constraint water depth h downstream of Longkou, the second constraint water depth h downstream of Longkou and the assumed water depth h downstream of Longkou s Equality is the second constraint.

[0130] In one embodiment, the value of the first ratio is determined according to actual needs.

[0131] Preferably, the first ratio is 0.7.

[0132] In this embodiment, when When ≥0.7, the dragon mouth discharge is submerged flow, The inundation coefficient σ is obtained by piecewise linear interpolation from the Pavlovsky inundation coefficient table. The formula for determining the average water-passing width of the Longkou axis section is expressed as follows:

[0133] ,

[0134] in, Indicates the discharge flow from the dragon mouth, represents the flooding coefficient, m represents the discharge coefficient, which is a fixed value, H represents the assumed water depth upstream of Longkou, g represents the acceleration of gravity, It indicates the average water width of the Longkou axis section.

[0135] The formula for determining the second constraint water depth downstream of Longkou is expressed as:

[0136] ,

[0137] Among them, B represents the width of the dragon mouth bottom, represents the slope ratio of the dike head, and h represents the second constrained water depth downstream of Longkou.

[0138] The second constraint is expressed as h = h s , that is, the second constrained water depth downstream of Longkou is equal to the assumed water depth downstream of Longkou.

[0139] Alternatively, as Figure 3As shown, the second constraint condition constructed according to the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth and the assumed dragon mouth downstream water depth also includes:

[0140] When the first ratio between the assumed water depth downstream of the dragon mouth and the assumed water depth upstream of the dragon mouth is less than the preset ratio, the average water-passing width of the dragon mouth axis section is determined according to the first ratio by the non-submerged flow rate formula.

[0141] The third constrained water depth downstream of the Longkou is determined based on the average water-passing width of the Longkou axis section.

[0142] When the third constrained water depth downstream of the dragon mouth is greater than the assumed water depth downstream of the dragon mouth, the second constraint condition includes that the third constrained water depth downstream of the dragon mouth is equal to the critical flow water depth.

[0143] When the third constrained water depth downstream of the dragon mouth is less than or equal to the assumed water depth downstream of the dragon mouth, the second constraint condition includes that the third constrained water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth.

[0144] When the first ratio of the assumed water depth downstream of Longkou to the assumed water depth upstream of Longkou is less than the preset ratio, the Longkou discharge is a non-submerged flow, and the average water width of the Longkou axis section is solved according to the non-submerged flow discharge formula. ,Depend on Dragon mouth bottom width B and dike head slope ratio Solve for the second constraint water depth h downstream of Longkou. Construct the second constraint condition based on the second constraint water depth downstream of Longkou.

[0145] In one embodiment, the preset ratio is set to 0.7, then When <0.7, the Longkou discharge is non-submerged flow, and the formula for non-submerged flow discharge is:

[0146] ,

[0147] in, represents the discharge flow from Longkou, m represents the discharge coefficient, H represents the assumed water depth upstream of Longkou, g represents the acceleration of gravity, It indicates the average water width of the Longkou axis section.

[0148] The formula for determining the second constraint water depth downstream of Longkou is expressed as:

[0149] ,

[0150] Among them, B represents the width of the dragon mouth bottom, represents the slope ratio of the dike head, and h represents the second constrained water depth downstream of Longkou.

[0151] When h>h sWhen the second constraint condition is that the second constraint water depth h downstream of Longkou is equal to the critical flow depth h k .

[0152] When h≤h s When the second constraint condition includes the second constraint water depth h downstream of Longkou equal to the assumed water depth h downstream of Longkou s .

[0153] The critical water depth judgment formula is expressed as:

[0154] ,

[0155] in, Indicates the cross-sectional area corresponding to the critical water depth, Indicates the width of the water surface corresponding to the critical water depth.

[0156] Alternatively, as Figure 3 As shown, after performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition and the second constraint condition, the method further includes:

[0157] When the solution fails, a new water depth upstream of the dragon mouth and a new water depth downstream of the dragon mouth are assumed based on the determined dragon mouth bottom width.

[0158] Perform nonlinear programming solution on the new assumed water depth upstream of Longkou and the new assumed water depth downstream of Longkou until a successful solution is obtained.

[0159] In one embodiment, the nonlinear programming problem is solved by calling the slover solver in Excel through VBA. The independent variables include the assumed water depth H upstream of Longkou and the assumed water depth h downstream of Longkou. s The constraint equation includes the equation formed by the first constraint condition, and the second constraint condition is determined based on the relationship between the preset ratio and the first ratio, the second constraint water depth downstream of Longkou, and the assumed water depth downstream of Longkou, forming the constraint equation. Set the independent variables and constraint equations in the Excel planning and solving dialog box; then add a reference to solver.xlam in the development tool window; then create a new macro in the VBA code window, enter the code solversolve(True) in the macro, and you can call the solver to solve H and .

[0160] If the solution is successful, go to step S700. If the solution fails, return to step S100, assume a new upstream water depth and a new downstream water depth, and repeat steps S200-S600 until the successfully solved H and h are obtained. s .

[0161] In one embodiment, first use vba code to give H and Call the solver to solve the problem. According to the return value of the solversolve function, if the return value is 0, it means the solution is successful and the next step is entered. Otherwise, the program jumps to the next step. Assign a new initial value and continue solving. The initial value range can be determined based on H and The physical meaning is to select in increments of 2m within the range of the dike height.

[0162] Through the VBA program loop, the bottom width B of each dragon mouth is calculated in the same way as above, and the H and After each successful solution, H and , and the corresponding Furthermore, formulas are used to calculate the hydraulic parameters (drop, water depth, flow rate per unit width, velocity, etc.) for the dam bottom width B. Excel's built-in charting function then displays a schematic diagram of the changes in these parameters during the vertical flow blockage. This can be done in the form of a two-dimensional scatter plot, with the dam widths plotted on the horizontal axis and the hydraulic parameter data series plotted on the vertical axis. This scatter plot can be pre-set in the Excel menu window or created using Visual Basic (VBA) drawing functions.

[0163] In one embodiment, the solution comprised of steps S100 through S800 of the present invention can automatically assign initial values multiple times for nonlinear programming solutions, avoiding the need for manual intervention due to single evaluation failures. This solution supports fully automated solutions, iteratively calculating various hydraulic parameters without human intervention. Based on these parameters, graphs are plotted to display the changes in hydraulic parameters during the vertical flow blockage process. This approach not only improves the efficiency and stability of the solution process but also facilitates intuitive understanding of the dynamic changes in hydraulic parameters throughout the entire process.

[0164] An embodiment of the present invention provides a device for calculating hydraulic parameters of a vertical cut-off flow, comprising:

[0165] An assumption module is used to determine the bottom width of the dragon mouth, and to determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth;

[0166] A piecewise linear interpolation module is used to obtain the diversion flow of the diversion building by piecewise linear interpolation in the first water level-flow relationship curve of the diversion building according to the assumed water depth upstream of the dragon mouth;

[0167] A first constraint module is used to construct a first constraint condition based on the known interception design flow and the known upstream river channel storage flow, and based on the diversion flow of the diversion building and the second water level flow relationship curve downstream of the dragon mouth;

[0168] a discharge flow determination module, configured to determine the discharge flow of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, and the acquired dike parameters;

[0169] A second constraint module is used to construct a second constraint condition according to the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth and the assumed dragon mouth downstream water depth;

[0170] a nonlinear programming solution module, configured to perform nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition;

[0171] An association module is configured to, when the solution is successful, use the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth;

[0172] A solution module is used to obtain hydraulic parameters according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth.

[0173] like Figure 4 As shown, an electronic device 400 provided by an embodiment of the present invention includes a memory 410 and a processor 420; the memory 410 is used to store computer programs; the processor 420 is used to implement the above-mentioned method for calculating hydraulic parameters of vertical cutoff flow when executing the computer program.

[0174] In other words, an electronic device 400 includes a memory 410 and a processor 420 coupled to the memory 410; the memory 410 is configured to store a computer program; and the processor 420 is configured to perform the following operations when executing the computer program:

[0175] Determine the bottom width of the dragon mouth, and determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth;

[0176] According to the assumed water depth upstream of the dragon mouth, piecewise linear interpolation is performed in the first water level-discharge relationship curve of the diversion building to obtain the diversion flow of the diversion building;

[0177] According to the known interception design flow and the known upstream river channel storage flow, a first constraint condition is established based on the diversion flow of the diversion building and the second water level flow relationship curve downstream of the dragon mouth;

[0178] Determining the discharge of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth and the acquired dike parameters;

[0179] Constructing a second constraint condition based on a numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth, and the assumed dragon mouth downstream water depth;

[0180] Performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition;

[0181] When the solution is successful, the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth are used as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth;

[0182] Hydraulic parameters are obtained according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth, and the water depth downstream of the dragon mouth.

[0183] An embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the above-mentioned method for calculating hydraulic parameters of a vertical cut-off flow is implemented.

[0184] In other words, a non-volatile computer-readable storage medium stores a computer program, which, when executed by a processor, causes the processor to perform the following operations:

[0185] Determine the bottom width of the dragon mouth, and determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth;

[0186] According to the assumed water depth upstream of the dragon mouth, piecewise linear interpolation is performed in the first water level-discharge relationship curve of the diversion building to obtain the diversion flow of the diversion building;

[0187] According to the known interception design flow and the known upstream river channel storage flow, a first constraint condition is established based on the diversion flow of the diversion building and the second water level flow relationship curve downstream of the dragon mouth;

[0188] Determining the discharge of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth and the acquired dike parameters;

[0189] Constructing a second constraint condition based on a numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth, and the assumed dragon mouth downstream water depth;

[0190] Performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition;

[0191] When the solution is successful, the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth are used as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth;

[0192] Hydraulic parameters are obtained according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth, and the water depth downstream of the dragon mouth.

[0193] An electronic device 400 that can serve as a server or client of the present invention will now be described, which is an example of a hardware device that can be applied to various aspects of the present invention. The electronic device 400 is intended to represent various forms of digital electronic computer devices, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device 400 can also represent various forms of mobile devices, such as personal digital assistants, cellular phones, smart phones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0194] Electronic device 400 includes a computing unit that can perform various appropriate actions and processes based on a computer program stored in a read-only memory (ROM) or a computer program loaded from a storage unit into a random access memory (RAM). The RAM can also store various programs and data required for device operation. The computing unit, ROM, and RAM are interconnected via a bus. An input / output (I / O) interface is also connected to the bus.

[0195] Those skilled in the art will appreciate that all or part of the processes in the above-described method embodiments can be implemented by instructing the relevant hardware through a computer program. The program can be stored in a computer-readable storage medium. When executed, the program can include the processes in the above-described method embodiments. The storage medium can be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM). In this application, the units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they may be located in one location or distributed across multiple network elements. Some or all of these units can be selected based on actual needs to achieve the objectives of the embodiments of the present invention. Furthermore, the functional units in the various embodiments of the present invention can be integrated into a single processing unit, each unit can exist physically separately, or two or more units can be integrated into a single unit. These integrated units can be implemented in either hardware or software functional units.

[0196] Although the present invention is disclosed as above, the scope of protection of the present invention is not limited thereto. Those skilled in the art can make various changes and modifications without departing from the spirit and scope of the present invention, and these changes and modifications will all fall within the scope of protection of the present invention.

Claims

1. A calculation method for hydraulic parameters of vertical cutoff closure, characterized in that, include: Determine the bottom width of the dragon mouth, and determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth; According to the assumed water depth upstream of the dragon mouth, piecewise linear interpolation is performed in the first water level-flow relationship curve of the diversion building to obtain the diversion flow of the diversion building, including: obtaining a water level interval including the assumed water depth upstream of the dragon mouth in the first water level-flow relationship curve; obtaining a flow interval corresponding to the water level interval in the first water level-flow relationship curve; performing linear interpolation on an array consisting of the water level interval and the flow interval to obtain the diversion flow of the diversion building; According to the known interception design flow and the known upstream river channel storage flow, a first constraint condition is constructed based on the diversion flow of the diversion building and the second water level-flow relationship curve downstream of Longkou, including: obtaining the flow downstream of Longkou according to the interception design flow, the upstream river channel storage flow, and the diversion flow of the diversion building; according to the flow downstream of Longkou, obtaining a first constrained water depth downstream of Longkou through piecewise linear interpolation in the second water level-flow relationship curve; constructing the first constraint condition based on the first constrained water depth downstream of Longkou, wherein the first constraint condition includes: the first constrained water depth downstream of Longkou is equal to the assumed water depth downstream of Longkou; Determining the discharge of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth and the acquired dike parameters; A second constraint condition is constructed based on the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed water depth upstream of the dragon mouth, and the assumed water depth downstream of the dragon mouth, wherein the dike parameters include the dike head slope ratio; the second constraint condition is constructed based on the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed water depth upstream of the dragon mouth, and the assumed water depth downstream of the dragon mouth, including: When a first ratio between the assumed water depth downstream of the Longkou and the assumed water depth upstream of the Longkou is greater than or equal to a preset ratio, the inundation coefficient is obtained by piecewise linear interpolation in the Pavlovsky inundation coefficient table according to the first ratio; the average water-passing width of the Longkou axis section is determined according to the inundation coefficient, the known flow coefficient, the assumed water depth upstream of the Longkou, and the Longkou discharge; the second constrained water depth downstream of the Longkou is determined according to the average water-passing width of the Longkou axis section, the Longkou bottom width, and the dike head slope ratio; the second constraint condition is constructed according to the second constrained water depth downstream of the Longkou, wherein the second constraint condition includes: the second constrained water depth downstream of the Longkou is equal to the assumed water depth downstream of the Longkou; When the first ratio between the assumed water depth downstream of the Longkou and the assumed water depth upstream of the Longkou is less than a preset ratio, the average water-passing width of the Longkou axis section is determined according to the first ratio by using the non-submerged discharge formula; the third constrained water depth downstream of the Longkou is determined according to the average water-passing width of the Longkou axis section; when the third constrained water depth downstream of the Longkou is greater than the assumed water depth downstream of the Longkou, the second constraint condition includes that the third constrained water depth downstream of the Longkou is equal to the critical flow water depth; when the third constrained water depth downstream of the Longkou is less than or equal to the assumed water depth downstream of the Longkou, the second constraint condition includes that the third constrained water depth downstream of the Longkou is equal to the assumed water depth downstream of the Longkou; Performing nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition; When the solution is successful, the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth are used as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth; Hydraulic parameters are obtained according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth, and the water depth downstream of the dragon mouth.

2. The calculation method of hydraulic parameters for vertical cutoff closure according to claim 1, characterized in that Also includes: Determine different dragon mouth bottom widths, determine the corresponding assumed upstream dragon mouth water depth, assumed downstream dragon mouth water depth, first constraint condition and second constraint condition according to each dragon mouth bottom width, and perform nonlinear programming solution; When the solution is successful, the corresponding dragon mouth bottom width, water depth upstream of the dragon mouth and water depth downstream of the dragon mouth are calculated to obtain the corresponding hydraulic parameters under each dragon mouth bottom width.

3. The calculation method of the hydraulic parameters of the vertical cutoff closure according to claim 1, characterized in that, After performing nonlinear programming on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition, the method further includes: When the solution fails, assuming a new water depth upstream of the dragon mouth and a new water depth downstream of the dragon mouth under the determined dragon mouth bottom width; Perform nonlinear programming solution on the new assumed water depth upstream of Longkou and the new assumed water depth downstream of Longkou until a successful solution is obtained.

4. The hydraulic parameter calculation method for closing the gap by vertical cutoff according to claim 1, characterized in that, The hydraulic parameters include head height, water depth, single-width flow rate and flow velocity.

5. A calculation device for hydraulic parameters of vertical cutoff, characterized in that, include: An assumption module is used to determine the bottom width of the dragon mouth, and to determine an assumed water depth upstream of the dragon mouth and an assumed water depth downstream of the dragon mouth according to the bottom width of the dragon mouth; A piecewise linear interpolation module is configured to obtain the diversion flow of the diversion building by piecewise linear interpolation in the first water level-flow relationship curve of the diversion building according to the assumed water depth upstream of the dragon mouth, comprising: obtaining a water level interval including the assumed water depth upstream of the dragon mouth in the first water level-flow relationship curve; obtaining a flow interval corresponding to the water level interval in the first water level-flow relationship curve; and performing linear interpolation on an array consisting of the water level interval and the flow interval to obtain the diversion flow of the diversion building; A first constraint module is configured to construct a first constraint condition based on a known interception design flow and a known upstream channel storage flow, the diversion flow of the diversion structure, and a second water level-flow relationship curve downstream of the Longkou, including: obtaining a flow downstream of the Longkou according to the interception design flow, the upstream channel storage flow, and the diversion flow of the diversion structure; obtaining a first constrained water depth downstream of the Longkou by piecewise linear interpolation in the second water level-flow relationship curve based on the flow downstream of the Longkou; and constructing the first constraint condition based on the first constrained water depth downstream of the Longkou, wherein the first constraint condition includes: the first constrained water depth downstream of the Longkou is equal to the assumed water depth downstream of the Longkou; a discharge flow determination module, configured to determine the discharge flow of the dragon mouth according to the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, and the acquired dike parameters; The second constraint module is used to construct a second constraint condition based on the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth and the assumed dragon mouth downstream water depth, wherein the dike parameters include the dike head slope ratio; the second constraint condition constructed based on the numerical relationship between the dike parameters, the dragon mouth discharge, the assumed dragon mouth upstream water depth and the assumed dragon mouth downstream water depth includes: when the first ratio between the assumed dragon mouth downstream water depth and the assumed dragon mouth upstream water depth is greater than or equal to a preset ratio, according to the first ratio, in the Pavlovsky flooding coefficient table, obtain the flooding coefficient by piecewise linear interpolation; determine the average water-passing width of the dragon mouth axis section according to the flooding coefficient, the known flow coefficient, the assumed dragon mouth upstream water depth and the dragon mouth discharge; determine the dragon mouth axis section according to the average water-passing width, the dragon mouth bottom width and the dike head slope ratio a second constraint water depth downstream of the mouth; constructing the second constraint condition based on the second constraint water depth downstream of the dragon mouth, wherein the second constraint condition includes: the second constraint water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth; when the first ratio between the assumed water depth downstream of the dragon mouth and the assumed water depth upstream of the dragon mouth is less than a preset ratio, determining the average water-passing width of the dragon mouth axis section according to the non-submerged discharge volume formula based on the first ratio; determining the third constraint water depth downstream of the dragon mouth according to the average water-passing width of the dragon mouth axis section; when the third constraint water depth downstream of the dragon mouth is greater than the assumed water depth downstream of the dragon mouth, the second constraint condition includes that the third constraint water depth downstream of the dragon mouth is equal to the critical flow water depth; when the third constraint water depth downstream of the dragon mouth is less than or equal to the assumed water depth downstream of the dragon mouth, the second constraint condition includes that the third constraint water depth downstream of the dragon mouth is equal to the assumed water depth downstream of the dragon mouth; a nonlinear programming solution module, configured to perform nonlinear programming solution on the assumed water depth upstream of the dragon mouth, the assumed water depth downstream of the dragon mouth, the first constraint condition, and the second constraint condition; An association module is configured to, when the solution is successful, use the assumed water depth upstream of the dragon mouth and the assumed water depth downstream of the dragon mouth as the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth corresponding to the bottom width of the dragon mouth; A solution module is used to obtain hydraulic parameters according to the bottom width of the dragon mouth, the water depth upstream of the dragon mouth and the water depth downstream of the dragon mouth.

Images