CN110197030B - Regulating and controlling method for controlling brake flooding outflow - Google Patents

Regulating and controlling method for controlling brake flooding outflow Download PDF

Info

Publication number
CN110197030B
CN110197030B CN201910457604.1A CN201910457604A CN110197030B CN 110197030 B CN110197030 B CN 110197030B CN 201910457604 A CN201910457604 A CN 201910457604A CN 110197030 B CN110197030 B CN 110197030B
Authority
CN
China
Prior art keywords
gate
variable
formula
flow
solution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910457604.1A
Other languages
Chinese (zh)
Other versions
CN110197030A (en
Inventor
崔巍
穆祥鹏
陈文学
刘丰
乔雨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201910457604.1A priority Critical patent/CN110197030B/en
Publication of CN110197030A publication Critical patent/CN110197030A/en
Application granted granted Critical
Publication of CN110197030B publication Critical patent/CN110197030B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D9/00—Level control, e.g. controlling quantity of material stored in vessel
    • G05D9/12—Level control, e.g. controlling quantity of material stored in vessel characterised by the use of electric means
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Barrages (AREA)
  • Control Of Non-Electrical Variables (AREA)

Abstract

The invention relates to a check gate flooding deviceThe method for regulating and controlling the outflow comprises the following steps: preparing data; calculating a first round of correlation coefficient; calculating a second round of correlation coefficient; determining a gate overflowing formula form; calibrating a gate overflowing formula; determining the calibration precision of a gate overflowing formula; regulating and controlling overcurrent: and (4) calculating the opening degree of the gate by using a calibrated gate overflowing formula, and regulating and controlling the gate. The invention is based on a large amount of engineering case experience and providese、H u －H d 、e/H u Ande/H d to construct the factors, a check gate flooding outflow formula version of a regression parabolic parameter calibration equation is employed. The technical scheme of rapidly calibrating the gate overflowing formula through two rounds of correlation calculation and comparison is provided. The gate outflow formula is provided based on the hydraulic characteristics of the submergence outflow, the structure is flexible, the form is various, the applicability is wide, the fitting precision is high, and the regulation and control of the submergence outflow of the check gate based on the formula are more accurate, simple and convenient and easy to realize.

Description

Regulating and controlling method for controlling brake flooding outflow
Technical Field
The invention relates to a regulation and control method for controlling the flooding outflow of a sluice, which is a regulation and control method for hydraulic facilities and is a regulation and control method based on gate parameter regression calculation.
Background
The check gate is used as a flow and water level adjusting structure and is widely applied to water delivery channels and open channel water-adjusting projects. The gate overflowing formula describes the functional relation between the passing gate flow and the water level before the gate, the water level after the gate and the gate opening, and is the basis for gate regulation. According to whether the gate overflowing capacity is influenced by the water level behind the gate, the gate outflow flow state can be divided into free outflow and submerged outflow.
The regulation of the flooding outflow is relatively complex, since it is influenced by various factors,key to check brake flooding outflow regulation: the accuracy of the inundation outflow formula is difficult to determine. Firstly, the gate overflowing is a complex three-dimensional flow state, particularly, the gate passing contraction section has concentrated head loss, and accurate estimation is difficult. Secondly, under the submerged outflow flow state, along with the rising of the water level behind the floodgate, the curvature of the flow curve increases rapidly, and is very steep, easily causes big calibration deviation. Finally, the flow capacity of the check gate is affected by many environmental factors, such as the smoothness of the entrance section of the chamber, the shape of the pier nose, the sealing form and material of the gate edge, etc. Taking the north-south water regulation central line engineering regulating gate as an example, model tests show that the maximum operation regulation range of the submerged outflow reaches e/HuNear 1, far more conventional 0.65.
The traditional regulating mode of the check gate is based on the form of a free outflow formula, and the free outflow formula is multiplied by a submerging outflow coefficient representing submerging characteristics to obtain a submerging outflow formula for regulation. The inundation outflow coefficient adopts a fixed structural form, and parameters of the inundation outflow coefficient are obtained by prototype observation or test data calibration. In recent years, large and medium-sized gate projects are increased, environmental conditions and structural types are increasingly complex, the traditional formula with a fixed structure and a single form cannot fully reflect the actual project, and the calibration effect of the overcurrent relationship is poor. Therefore, it is urgently needed to develop a check gate inundation outflow calculation formula with flexible structure, various forms and wider applicability, and improve the regulation and control effect. In addition, the selection of the throttle valve overcurrent formula belongs to a process of time and labor waste due to repeated trial and comparison, and a technical scheme which is convenient and practical and can help quickly select the optimal throttle valve regulation and control form and effect is also needed.
Disclosure of Invention
In order to overcome the problems of the prior art, the invention provides a regulating method for controlling the inundation outflow of a damper. The method adopts a check gate submerging outflow formula type of a regression parabolic parameter calibration equation, and establishes a gate overflowing formula through two rounds of correlation calculation and comparison, so as to obtain more accurate gate regulation.
The purpose of the invention is realized as follows: a method of regulating a throttle gate flood effluent, the method comprising the steps of:
step 1, data preparation: at least n groups of measured data under different flow rates are arranged, including the flow rate Q of the passing gate and the water level H before the gateuWater level H behind gatedOpening e of the gate; setting single width B and fan number N of the gate, and enabling B to be Nb;
step 2, calculating a first round correlation coefficient: respectively calculate Hu－HdAnd
Figure BDA0002077074240000021
and H, andu－Hdand
Figure BDA0002077074240000022
a correlation coefficient between;
setting:
Figure BDA0002077074240000023
and 3, calculating a second round correlation coefficient:
if | R1|>|R2Respectively calculate
Figure BDA0002077074240000024
Between e and,
Figure BDA0002077074240000025
And e/HuA middle part of the upper part,
Figure BDA0002077074240000026
And e/HdA correlation coefficient between;
setting:
Figure BDA0002077074240000027
if | R1|≤|R2Respectively calculate
Figure BDA0002077074240000028
Between e and,
Figure BDA0002077074240000029
And e/HuA middle part of the upper part,
Figure BDA00020770742400000210
And e/HdA correlation coefficient between;
setting:
Figure BDA00020770742400000211
step 4, determining the gate overflowing formula form:
if | R1|>|R2I, solving: max (| R)1|，|R3|，|R4|，|R5|)；
If the solution is | R1L, with its corresponding variable Hu－HdAs the independent variable, there is a variable,
Figure BDA00020770742400000212
for dependent variables, the gate flow formula is determined in the form of
Figure BDA00020770742400000213
Wherein f (H)u－Hd) Is a regression parabolic function;
if the solution is | R3Taking the corresponding variable e as an independent variable,
Figure BDA00020770742400000214
determining gate flow formula form for dependent variable
Figure BDA00020770742400000215
Wherein f (e) is a regressive parabolic function;
if the solution is | R4L, with its corresponding variable e/HuAs the independent variable, there is a variable,
Figure BDA00020770742400000216
determining gate flow formula form for dependent variable
Figure BDA0002077074240000031
Wherein f (e/H)u) Is a regression parabolic function;
if the solution is | R5L, with its corresponding variable e/HdAs the independent variable, there is a variable,
Figure BDA0002077074240000032
determining gate flow formula form for dependent variable
Figure BDA0002077074240000033
Wherein f (e/H)d) Is a regression parabolic function;
if | R1|≤|R2I, solving: max (| R)2|，|R6|，|R7|，|R8|)；
If the solution is | R2L, with its corresponding variable Hu－HdAs the independent variable, there is a variable,
Figure BDA0002077074240000034
determining gate flow formula form for dependent variable
Figure BDA0002077074240000035
Wherein f (H)u-Hd) Is a regression parabolic function;
if the solution is | R6Taking the corresponding variable e as an independent variable,
Figure BDA0002077074240000036
determining gate flow formula form for dependent variable
Figure BDA0002077074240000037
Wherein f (e) is a regressive parabolic function;
if the solution is | R7L, with its corresponding variable e/HuAs the independent variable, there is a variable,
Figure BDA0002077074240000038
determining gate flow formula form for dependent variable
Figure BDA0002077074240000039
Wherein f (e/H)u) Is a regression parabolic function;
if the solution is | R8L, with its corresponding variable e/HdAs the independent variable, there is a variable,
Figure BDA00020770742400000310
determining gate flow formula form for dependent variable
Figure BDA00020770742400000311
Wherein f (e/H)d) Is a regression parabolic function;
step 5, calibrating a gate overflow formula:
using cubic polynomial f (x) ═ (ax)3+bx2+ cx + d) rating coefficients a, b, c, d:
1) for solution to | R1Case of |:
let Hu-HdAs an independent variable x, with
Figure BDA00020770742400000312
Drawing a y-x scatter diagram for the dependent variable y, fitting a regression type parabola trend line, and calibrating coefficients a, b, c and d of the polynomial;
the gate flow formula is:
Figure BDA00020770742400000313
2) for solution to | R3Case of |:
let e be an argument x, to
Figure BDA00020770742400000314
Drawing a y-x scatter diagram for a dependent variable y, fitting a regression parabola trend line, and respectively determining coefficients of a polynomial as a, b, c and d;
the gate flow formula is:
Figure BDA00020770742400000315
3) for solution to | R4Case of |：
Let e/HuAs an independent variable x, with
Figure BDA00020770742400000316
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure BDA0002077074240000041
4) for solution to | R5Case of |:
let e/HdAs an independent variable x, with
Figure BDA0002077074240000042
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure BDA0002077074240000043
5) for solution to | R2Case of |:
let Hu-HdAs an independent variable x, with
Figure BDA0002077074240000044
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure BDA0002077074240000045
6) for solution to | R6Case of |:
with e as the argument x, with
Figure BDA0002077074240000046
Drawing a scatter diagram of y-x for the dependent variable yFitting a regression parabola trend line, and respectively determining coefficients of the polynomial as a, b, c and d;
the gate flow formula is:
Figure BDA0002077074240000047
7) for solution to | R7Case of |:
at e/HuAs an independent variable x, with
Figure BDA0002077074240000048
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure BDA0002077074240000049
8) for solution to | R8Case of |:
at e/HdAs an independent variable x, with
Figure BDA00020770742400000410
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure BDA00020770742400000411
step 6, determining the calibration precision of the gate overflowing formula: calculating the brake-passing flow by using a calibrated gate overflowing formula, subtracting the measured flow from the brake-passing flow, and dividing the difference by the measured flow to obtain a calibrated relative percentage error;
step 7, regulating and controlling overcurrent: the opening degree of the gate is calculated by using a calibrated gate overflowing formula, and the gate is regulated and controlled: firstly, according to the gate opening e (t), the passing gate flow Q (t) and the flow control target value Q (t +1) at the current moment, determining a gate opening control increment:
Figure BDA0002077074240000051
wherein
Figure BDA0002077074240000052
The derivative is obtained by adopting the gate overflowing formula type determined in the step 5, and the current H is inputu(t)、HdAnd (t), determining a gate control target opening e (t +1) ═ e (t) + △ e based on the current gate opening e (t), and finally inputting a gate opening adjusting command to a mechanical system and adjusting the gate opening to e (t + 1).
The invention has the following beneficial effects: based on a great deal of engineering case experience, the invention provides the method by e and Hu-Hd、e/HuAnd e/HdFor constructing factors, a check gate submerging outflow formula type of a regression parabolic parameter calibration equation is adopted, and a gate overflow formula is calibrated rapidly through two rounds of correlation calculation and comparison. The gate outflow formula is provided based on the hydraulic characteristics of the submergence outflow, the structure is flexible, the form is various, the applicability is wide, the fitting precision is high, and the regulation and control of the submergence outflow of the check gate based on the formula are more accurate, simple and convenient and easy to realize. The formula provided by the invention can be applied to a regulation and control brake and can also be applied to the form of selecting the overflowing formula of the gate in engineering design.
Drawings
The invention is further illustrated by the following figures and examples.
FIG. 1 is a schematic illustration of gate parameters of a process according to an embodiment of the present invention;
FIG. 2 is a flow chart of a method according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of regression parabolic parameter fitting of the application example 1 according to the embodiment of the present invention;
fig. 4 is a schematic diagram of regression parabolic parameter fitting of the application example 2 according to the embodiment of the present invention.
Detailed Description
Example (b):
this embodiment is a check gateAnd (4) a regulation and control method of flooding outflow. The control of the throttling gate mainly uses an overflow formula to calculate the opening of the gate, so that the accurate calibration of the overflow formula is the key of the control according to actual operation. The over-current formula describes the over-current flow Q and the water level H before the gateuWater level H after gatedAnd the opening e of the gate, as shown in fig. 1, these parameters are the basis for the gate control.
Based on a great deal of engineering case experience, the invention provides the method by e and Hu－Hd、e/HuAnd e/HdFor constructing the factors, a throttle valve flooding outflow formula type of a regression parabolic parameter calibration equation is adopted, and the process is shown in fig. 2 and specifically comprises the following steps:
step 1, data preparation:
n groups (n) of measured data under different flow rates are sorted>10, and uniformly covering the flow interval of the submerged flow operation), including the flow Q of the passing gate and the water level H before the gateuWater level H behind gatedThe gate opening e is set to B equal to Nb and set to g equal to 9.81, with the given gate width B and gate fan number N.
Step 2, calculating a first round correlation coefficient (also called Pearson correlation coefficient):
respectively calculate (H)u－Hd) And
Figure BDA0002077074240000061
and H, andu－Hdand
Figure BDA0002077074240000062
the correlation coefficient between them.
Setting:
Figure BDA0002077074240000063
and 3, calculating a second round correlation coefficient:
if | R1|>|R2Respectively calculate
Figure BDA0002077074240000064
Between e and,
Figure BDA0002077074240000065
And e/HuA middle part of the upper part,
Figure BDA0002077074240000066
And e/HdThe correlation coefficient between them.
Setting:
Figure BDA0002077074240000067
if | R1|≤|R2Respectively calculate
Figure BDA0002077074240000068
Between e and,
Figure BDA0002077074240000069
And e/HuA middle part of the upper part,
Figure BDA00020770742400000610
And e/HdThe correlation coefficient between them.
Setting:
Figure BDA00020770742400000611
step 4, determining the gate overflowing formula form:
if | R1|>|R2I, solving: max (| R)1|，R3|，|R4|，|R5|)
If the solution is | R1I, |, with its corresponding variable, Hu－HdAs the independent variable, there is a variable,
Figure BDA00020770742400000612
for dependent variables, the gate flow formula is determined in the form of
Figure BDA00020770742400000613
Wherein f (H)u-Hd) Is a regression parabolic function.
If the solution is | R3Taking the corresponding variable e as an independent variable,
Figure BDA00020770742400000614
determining gate flow formula form for dependent variable
Figure BDA00020770742400000615
Wherein f (e) is a regressive parabolic function.
If the solution is | R4L, with its corresponding variable e/HuAs the independent variable, there is a variable,
Figure BDA00020770742400000616
determining gate flow formula form for dependent variable
Figure BDA00020770742400000617
Wherein f (e/H)u) Is a regression parabolic function.
If the solution is | R5L, with its corresponding variable e/HdAs the independent variable, there is a variable,
Figure BDA00020770742400000618
determining gate flow formula form for dependent variable
Figure BDA0002077074240000071
Wherein f (e/H)d) Is a regression parabolic function.
If | R1|≤|R2I, solving: max (| R)2|，|R6|，|R7|，|R8|)。
If the solution is | R2L, with its corresponding variable Hu－HdAs the independent variable, there is a variable,
Figure BDA0002077074240000072
determining gate flow formula form for dependent variable
Figure BDA0002077074240000073
Wherein f (H)u-Hd) Is a regression parabolic function.
If the solution is | R6Taking the corresponding variable e as an independent variable,
Figure BDA0002077074240000074
determining gate flow formula form for dependent variable
Figure BDA0002077074240000075
Wherein f (e) is a regressive parabolic function.
If the solution is | R7L, with its corresponding variable e/HuAs the independent variable, there is a variable,
Figure BDA0002077074240000076
determining gate flow formula form for dependent variable
Figure BDA0002077074240000077
Wherein f (e/H)u) Is a regression parabolic function.
If the solution is | R8L, with its corresponding variable e/HdAs the independent variable, there is a variable,
Figure BDA0002077074240000078
determining gate flow formula form for dependent variable
Figure BDA0002077074240000079
Wherein f (e/H)d) Is a regression parabolic function.
Step 5, calibrating a gate overflow formula:
for the previous step, solve to | R1In the case of |, let x be Hu-HdThen f (x) becomes (ax)3+bx2+ cx + d), a, b, c, d are coefficients to be calibrated. To be provided with
Figure BDA00020770742400000710
And drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively determining coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA00020770742400000711
for the previous step, solve to | R3If x is equal to e, f (x) is equal to ax3+bx2+ cx + d), a, b, c, d are coefficients to be calibrated. To be provided with
Figure BDA00020770742400000712
And drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively determining coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA00020770742400000713
for the previous step, solve to | R4In case of l, in e/HuAs an independent variable, with
Figure BDA00020770742400000714
And drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively setting coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA00020770742400000715
for the previous step, solve to | R5In case of l, in e/HdAs an independent variable x, with
Figure BDA00020770742400000716
And drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively setting coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA00020770742400000717
for the previous step, solve to | R2In case of |, by Hu-HdAs an independent variable x, with
Figure BDA0002077074240000081
Is a dependent variableAnd y, drawing a y-x scatter diagram, fitting a regression parabola trend line, and respectively determining coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA0002077074240000082
for the previous step, solve to | R6In case of l, e is an independent variable, and
Figure BDA0002077074240000083
and drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively setting coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA0002077074240000084
for the previous step, solve to | R7In case of l, in e/HuAs an independent variable x, with
Figure BDA0002077074240000085
And drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively setting coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA0002077074240000086
for the previous step, solve to | R8In case of l, in e/HdAs an independent variable x, with
Figure BDA0002077074240000087
And drawing a y-x scatter diagram for the dependent variable y, fitting a regression parabola trend line, and respectively setting coefficients of the polynomial as a, b, c and d. The gate flow formula is:
Figure BDA0002077074240000088
step 6, determining the calibration precision of the gate overflowing formula:
and calculating the brake-passing flow by using a calibrated gate overflowing formula, subtracting the actual measurement flow from the brake-passing flow, and dividing the difference by the actual measurement flow to obtain a calibrated relative percentage error.
Step 7, regulating and controlling overcurrent: and (4) calculating the opening degree of the gate by using a calibrated gate overflowing formula, and regulating and controlling the gate. Firstly, according to the gate opening e (t), the passing gate flow Q (t) and the flow control target value Q (t +1) at the current moment, determining a gate opening control increment:
Figure BDA0002077074240000089
wherein
Figure BDA00020770742400000810
The derivative is obtained by adopting the gate overflowing formula type determined in the step 5, and the current H is inputu(t)、HdAnd (t), determining a gate control target opening e (t +1) ═ e (t) + △ e based on the current gate opening e (t), and finally inputting a gate opening adjusting command to a mechanical system and adjusting the gate opening to e (t + 1).
Application example 1:
the method comprises the following steps of regulating a central line engineering ice cream inverted siphon outlet regulating gate in south-to-north water, setting a gate overflowing formula and regulating and controlling overflowing:
a) preparing data: as shown in Table 1, the (1) th to (6) th rows are the water level before the gate HuWater level H behind gatedGate opening e, passing gate flow Q, single gate width b, and number of parallel gates N.
TABLE 1 Ice cream inverted siphon outlet throttle gate flow formula calibration table
Figure BDA0002077074240000091
b) The first round of correlation coefficient calculation: (H) was calculated separately using the CORREL function of Execlu－Hd) And
Figure BDA0002077074240000092
and H, andu－Hdand
Figure BDA0002077074240000093
the correlation coefficient between them was found to be 0.75 and-0.63, respectively.
c) And calculating a second round of correlation coefficient: since |0.75>Respectively calculating with CORREL function of Excel to 0.63 |)
Figure BDA0002077074240000094
Between e and,
Figure BDA0002077074240000095
And e/HuA middle part of the upper part,
Figure BDA0002077074240000096
And e/HdThe results are 0.09, -0.98 and 0.07 respectively.
d) Determining the form of a gate overflowing formula: since the result of Max (|0.75|, |0.09|, | -0.98|, |0.07|) is 0.98, with its corresponding variable e/HuAs the independent variable, there is a variable,
Figure BDA0002077074240000101
determining the form of gate flow equation for dependent variable
Figure BDA0002077074240000102
e) Calibrating a gate overflowing formula; to be provided with
Figure BDA0002077074240000103
(i.e., column (8)) as a dependent variable y in e/Hu(i.e., column (9)) is an independent variable x, a scatter plot is plotted and a regression parabolic trend line is fitted (see fig. 3), and the coefficients of the polynomial are defined as a, b, c, and d, i.e., a-7.2152, b-13.35, c-9.4253, and d-2.8672, respectively. The gate flow formula determined by the ratio is as follows:
Figure BDA0002077074240000104
f) determining the rating precision of a gate overflowing formula: calculating the gate flow Q by using a calibrated gate overflowing formulavSee column (12), calculate its phaseFor percent error, see (13), the average is 2.59%. Better than the conventional formula of 10% -20%.
g) And (4) calculating the opening degree of the gate by using a calibrated gate overflowing formula, and regulating and controlling the gate. Firstly, according to the opening e (t) of the gate at the current moment, the passing gate flow Q (t), the water level H before the gateu(t), post-gate water level Hd(t) and a flow control target value Q (t +1), determining a gate opening control increment:
Figure BDA0002077074240000105
wherein:
Figure BDA0002077074240000106
Figure BDA0002077074240000107
and finally, inputting a gate opening adjusting command to a mechanical system, and adjusting the gate opening to e (t + 1).
Application example 2:
a south-to-north water regulation central line project Niuhe Henan Branch regulating gate comprises the following steps of:
a) data preparation, as shown in Table 2, the (1) th to (6) th columns are the pre-gate water levels HuWater level H behind gatedGate opening e, passing gate flow Q, single gate width b, and number of parallel gates N.
TABLE 2 rating table for brake control overcurrent formula of Henan Branch of Niuhuan province
Figure BDA0002077074240000108
Figure BDA0002077074240000111
b) The first round of correlation coefficient calculation. (H) was calculated separately using the CORREL function of Execlu－Hd) And
Figure BDA0002077074240000112
and H, andu－Hdand
Figure BDA0002077074240000113
the correlation coefficient between them was found to be 0.75 and-0.96, respectively.
c) And calculating the correlation coefficient in the second round. Because | -0.96|>Respectively calculating with CORREL function of Excel to |0.75 |)
Figure BDA0002077074240000114
Between e and,
Figure BDA0002077074240000115
And e/HuA middle part of the upper part,
Figure BDA0002077074240000116
And e/HdThe results were 0.68, 0.69 and 0.43, respectively.
d) And determining the form of a gate overflowing formula. Since max (| -0.96|, |0.75|, |0.69|, |0.68|) results in 0.96, with its corresponding variable Hu－HdAs the independent variable, there is a variable,
Figure BDA0002077074240000117
determining gate flow formula form for dependent variable
Figure BDA0002077074240000121
Wherein f (H)u-Hd) Is a regression parabolic function.
e) And calibrating a gate overflowing formula. To be provided with
Figure BDA0002077074240000122
(i.e., column (8)) as a dependent variable y, and Hu－Hd(i.e., column (9)) is an independent variable x, a scatter plot is plotted, and a regression formula parabola trend line (see fig. 4) is fitted, and the coefficients of the polynomial are respectively designated as a, b, c, and d, i.e., a is 0.4545, b is 2.0648, c is 2.8611, and d is 0.4082. The gate flow formula determined by the ratio is as follows:
Figure BDA0002077074240000123
f) and determining the calibration precision of the gate overflowing formula. Calculating the gate flow Q by using a calibrated gate overflowing formulavSee column (12) and calculate the relative percent error, see (13), with an average of 1.85%. Better than the conventional formula of 10% -20%.
g) And (4) calculating the opening degree of the gate by using a calibrated gate overflowing formula, and regulating and controlling the gate. Firstly, according to the opening e (t) of the gate at the current moment, the passing gate flow Q (t), the water level H before the gateu(t), post-gate water level Hd(t) and a flow control target value Q (t +1), determining a gate opening control increment:
Figure BDA0002077074240000124
wherein:
Figure BDA0002077074240000125
then, based on the current gate opening e (t), determining a gate control target opening: e (t +1) ═ e (t) +. ae; and finally, inputting a 15-gate opening adjusting command to a mechanical system, and adjusting the gate opening to e (t + 1).
Finally, it should be noted that the above is only for illustrating the technical solution of the present invention and not for limiting, and although the present invention has been described in detail with reference to the preferred arrangement, it should be understood by those skilled in the art that modifications or equivalent substitutions can be made on the technical solution of the present invention (such as the form of the check gate, the application of various formulas, the sequence of steps, etc.) without departing from the spirit and scope of the technical solution of the present invention.

Claims (1)

1. A method of regulating a throttle gate flood effluent, the method comprising the steps of:
step 1, data preparation: arranging at least the flow rates of different flowsn groups of measured data including the flow Q of passing gate and the water level H before gateuWater level H behind gatedOpening e of the gate; setting single width B and fan number N of the gate, and enabling B to be Nb;
step 2, calculating a first round correlation coefficient: respectively calculate Hu－HdAnd
Figure FDA0002562714530000011
and H, andu－Hdand
Figure FDA0002562714530000012
a correlation coefficient between;
setting:
Figure FDA0002562714530000013
wherein: r () is a functional relational expression of the correlation coefficient; g is the acceleration of gravity;
and 3, calculating a second round correlation coefficient:
if | R1|＞|R2Respectively calculate
Figure FDA0002562714530000014
Between e and,
Figure FDA0002562714530000015
And
Figure FDA0002562714530000016
a middle part of the upper part,
Figure FDA0002562714530000017
And
Figure FDA0002562714530000018
a correlation coefficient between;
setting:
Figure FDA0002562714530000019
if | R1|≤|R2Respectively calculate
Figure FDA00025627145300000110
Between e and,
Figure FDA00025627145300000111
And
Figure FDA00025627145300000112
a middle part of the upper part,
Figure FDA00025627145300000113
And
Figure FDA00025627145300000114
a correlation coefficient between;
setting:
Figure FDA00025627145300000115
step 4, determining the gate overflowing formula form:
if | R1|＞|R2I, solving: max (| R)1|，|R3|，|R4|，|R5|)；
If the solution is | R1L, with its corresponding variable Hu－HdAs the independent variable, there is a variable,
Figure FDA00025627145300000116
for dependent variables, the gate flow formula is determined in the form of
Figure FDA00025627145300000117
Wherein f (H)u－Hd) Is a regression parabolic function;
if the solution is | R3Taking the corresponding variable e as an independent variable,
Figure FDA00025627145300000118
determining gate flow formula form for dependent variable
Figure FDA0002562714530000021
Wherein f (e) is a regressive parabolic function;
if the solution is | R4L, in its corresponding variable
Figure FDA0002562714530000022
As the independent variable, there is a variable,
Figure FDA0002562714530000023
determining gate flow formula form for dependent variable
Figure FDA0002562714530000024
Wherein
Figure FDA0002562714530000025
Is a regression parabolic function;
if the solution is | R5L, in its corresponding variable
Figure FDA0002562714530000026
As the independent variable, there is a variable,
Figure FDA0002562714530000027
determining gate flow formula form for dependent variable
Figure FDA0002562714530000028
Wherein
Figure FDA0002562714530000029
Is a regression parabolic function;
if | R1|≤|R2I, solving: max (| R)2|，|R6|，|R7|，|R8|)；
If the solution is | R2In its correspondingVariable Hu－HdAs the independent variable, there is a variable,
Figure FDA00025627145300000210
determining gate flow formula form for dependent variable
Figure FDA00025627145300000211
Wherein f (H)u-Hd) Is a regression parabolic function;
if the solution is | R6Taking the corresponding variable e as an independent variable,
Figure FDA00025627145300000212
determining gate flow formula form for dependent variable
Figure FDA00025627145300000213
Wherein f (e) is a regressive parabolic function;
if the solution is | R7L, with its corresponding variable e/HuAs the independent variable, there is a variable,
Figure FDA00025627145300000214
determining gate flow formula form for dependent variable
Figure FDA00025627145300000215
Wherein
Figure FDA00025627145300000216
Is a regression parabolic function;
if the solution is | R8L, in its corresponding variable
Figure FDA00025627145300000217
As the independent variable, there is a variable,
Figure FDA00025627145300000218
determining gate flow formula form for dependent variable
Figure FDA00025627145300000219
Wherein
Figure FDA00025627145300000220
Is a regression parabolic function;
step 5, calibrating a gate overflow formula:
using cubic polynomial f (x) ═ (ax)3+bx2+ cx + d) rating coefficients a, b, c, d:
1) for solution to | R1Case of |:
let Hu-HdAs an independent variable x, with
Figure FDA0002562714530000031
Drawing a y-x scatter diagram for the dependent variable y, fitting a regression type parabola trend line, and calibrating coefficients a, b, c and d of the polynomial;
the gate flow formula is:
Figure FDA0002562714530000032
2) for solution to | R3Case of |:
let e be an argument x, to
Figure FDA0002562714530000033
Drawing a y-x scatter diagram for a dependent variable y, fitting a regression parabola trend line, and respectively determining coefficients of a polynomial as a, b, c and d;
the gate flow formula is:
Figure FDA0002562714530000034
3) for solution to | R4Case of |:
order to
Figure FDA0002562714530000035
As an independent variable x, with
Figure FDA0002562714530000036
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure FDA0002562714530000037
4) for solution to | R5Case of |:
order to
Figure FDA0002562714530000038
As an independent variable x, with
Figure FDA0002562714530000039
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure FDA00025627145300000310
5) for solution to | R2Case of |:
let Hu-HdAs an independent variable x, with
Figure FDA00025627145300000311
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure FDA00025627145300000312
6) for solution to | R6Case of |:
with e as the argument x, with
Figure FDA00025627145300000313
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure FDA00025627145300000314
7) for solution to | R7Case of |:
to be provided with
Figure FDA00025627145300000315
As an independent variable x, with
Figure FDA00025627145300000316
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure FDA0002562714530000041
8) for solution to | R8Case of |:
to be provided with
Figure FDA0002562714530000042
As an independent variable x, with
Figure FDA0002562714530000043
Drawing a y-x scatter diagram for the dependent variable y, and fitting a regression type parabola trend line, wherein coefficients of a polynomial are respectively defined as a, b, c and d;
the gate flow formula is:
Figure FDA0002562714530000044
step 6, determining the calibration precision of the gate overflowing formula: calculating the brake-passing flow by using a calibrated gate overflowing formula, subtracting the measured flow from the brake-passing flow, and dividing the difference by the measured flow to obtain a calibrated relative percentage error;
step 7, regulating and controlling overcurrent: calculating the opening of the gate by using a calibrated gate overflowing formula, and regulating and controlling the gate; firstly, according to the gate opening e (t), the passing gate flow Q (t) and the flow control target value Q (t +1) at the current moment, determining a gate opening control increment:
Figure FDA0002562714530000045
wherein:
Figure FDA0002562714530000046
the derivative is obtained by adopting the gate overflowing formula type determined in the step 5, and the current H is inputu(t)、Hd(t) calculating; then, based on the current gate opening e (t), determining a gate control target opening: e (t +1) ═ e (t) + Δ e; and finally, inputting a gate opening adjusting instruction to the mechanical system, and adjusting the gate opening to e (t + 1).
CN201910457604.1A 2019-05-29 2019-05-29 Regulating and controlling method for controlling brake flooding outflow Active CN110197030B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910457604.1A CN110197030B (en) 2019-05-29 2019-05-29 Regulating and controlling method for controlling brake flooding outflow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910457604.1A CN110197030B (en) 2019-05-29 2019-05-29 Regulating and controlling method for controlling brake flooding outflow

Publications (2)

Publication Number Publication Date
CN110197030A CN110197030A (en) 2019-09-03
CN110197030B true CN110197030B (en) 2020-08-21

Family

ID=67753378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910457604.1A Active CN110197030B (en) 2019-05-29 2019-05-29 Regulating and controlling method for controlling brake flooding outflow

Country Status (1)

Country Link
CN (1) CN110197030B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874622A (en) * 2017-03-13 2017-06-20 中国水利水电科学研究院 One kind series connection Duo Qu ponds gate lock flow coefficient rating method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2600180A1 (en) * 1986-06-13 1987-12-18 Alsthom Sluice gate for automatic regulation of an upstream or downstream level
CN108181940B (en) * 2017-12-27 2020-08-04 中国水利水电科学研究院 Gate regulating and controlling method suitable for downstream emergency water cut-off condition of series channel
CN108534843B (en) * 2018-03-02 2019-12-24 武汉大学 Single gate flow calibration method and device for open channel water delivery channel
CN108867574B (en) * 2018-05-29 2020-10-16 广东省水利水电科学研究院 Calibration model sluice for numerical simulation and generalization method thereof
CN109783934A (en) * 2019-01-15 2019-05-21 广东省水文局惠州水文分局 A kind of mean velocity in section fitting rating method based on H-ADCP

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874622A (en) * 2017-03-13 2017-06-20 中国水利水电科学研究院 One kind series connection Duo Qu ponds gate lock flow coefficient rating method

Also Published As

Publication number Publication date
CN110197030A (en) 2019-09-03

Similar Documents

Publication Publication Date Title
CN106527519B (en) Control method for variable speed pressure flow field of large direct current temporary impulse type supersonic wind tunnel
CN110285403B (en) Main steam temperature control method based on controlled parameter estimation
WO2021196997A1 (en) Method for predicting and controlling water level of series of water conveyance canals on basis of fuzzy neural network
CN107966902A (en) A kind of constraint 2D tracking and controlling methods of uncertainty batch process
CN105425612B (en) A kind of method for optimizing of Adaptive System of Water-Turbine Engine control parameter
CN108181940B (en) Gate regulating and controlling method suitable for downstream emergency water cut-off condition of series channel
CN102520617B (en) Prediction control method for unminimized partial decoupling model in oil refining industrial process
CN110197030B (en) Regulating and controlling method for controlling brake flooding outflow
CN107918283A (en) A kind of batch injection moulding process is jamproof to mix stability controller design method
CN102520616B (en) Partial decoupling unminimized model prediction function control method in oil refining industrial process
CN103257658A (en) Feeding tank temperature control method and system
CN109711013B (en) Electrohydraulic actuator model for water turbine adjusting system
CN113267314A (en) Supersonic flow field total pressure control system of temporary-impulse wind tunnel
CN103616815A (en) Control method for waste plastic oil refining cracking furnace chamber temperature based on dynamic matrix control optimization
CN105344720A (en) Online control method for finish rolling temperature of precision rolling strip steel
CN103605284B (en) The cracking waste plastics stove hearth pressure control method that dynamic matrix control is optimized
CN110673482A (en) Power station coal-fired boiler intelligent control method and system based on neural network prediction
Clemmens et al. Computational procedures used for radial gate calibration in WinGate
CN107942669A (en) The limited rolling time horizon of batch injection moulding process mixes tracking and controlling method
CN108803342B (en) Unit unit load quick response prediction control method
CN105318312A (en) Method of controlling steam temperature of power station boiler through multi-step advance prediction of full working conditions
CN110705109A (en) Method for determining equal-percentage throttling window molded line of sleeve regulating valve
CN115167308A (en) Water transfer engineering multi-target prediction control algorithm for guiding gate regulation and control
CN105138041A (en) Method for controlling main steam temperature of thermal power generating unit for implementing DCS
CN115045854A (en) Parallel fan balance control method based on machine learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant