CN112561213A - Reservoir flood control scheduling method, system, equipment and medium based on improved POA - Google Patents

Reservoir flood control scheduling method, system, equipment and medium based on improved POA Download PDF

Info

Publication number
CN112561213A
CN112561213A CN202110199694.6A CN202110199694A CN112561213A CN 112561213 A CN112561213 A CN 112561213A CN 202110199694 A CN202110199694 A CN 202110199694A CN 112561213 A CN112561213 A CN 112561213A
Authority
CN
China
Prior art keywords
reservoir
period
time interval
discharge
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110199694.6A
Other languages
Chinese (zh)
Other versions
CN112561213B (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.)
Zhejiang Guiren Information Technology Co ltd
Original Assignee
Zhejiang Guiren Information Technology Co ltd
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 Zhejiang Guiren Information Technology Co ltd filed Critical Zhejiang Guiren Information Technology Co ltd
Priority to CN202110199694.6A priority Critical patent/CN112561213B/en
Publication of CN112561213A publication Critical patent/CN112561213A/en
Application granted granted Critical
Publication of CN112561213B publication Critical patent/CN112561213B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Computational Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Operations Research (AREA)
  • General Business, Economics & Management (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Primary Health Care (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Sewage (AREA)
  • Flow Control (AREA)

Abstract

The invention relates to a reservoir flood control dispatching method, a reservoir flood control dispatching system, reservoir flood control dispatching equipment and a reservoir flood control dispatching medium based on improved POA, wherein the method comprises the following steps: firstly, determining the starting water level of each reservoir in a flow area; secondly, solving a time interval discharge interval of the reservoir according to a water quantity balance equation and a constraint condition; then, calculate the firstnAverage warehousing flow of time intervals in each reservoir scheduling period; then, judgenWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge flow interval or not is determined according to the judgment resultnInitial solution of each time interval discharge of each reservoir; finally, according tonObtaining an initial dispatching line by initial solution of the initial water level of each reservoir, the discharge of each period of the reservoir and the period of the warehousing flow, and determining a total target value of a dispatching period by solving a total target value of the dispatching period based on the initial dispatching line and a water balance equationTo make surenA dispatch line for each reservoir. The technical scheme of the invention simplifies the solving process of the initial solution, can quickly generate the initial solution meeting the conditions, and improves the solving efficiency of the initial solution.

Description

Reservoir flood control scheduling method, system, equipment and medium based on improved POA
Technical Field
The invention relates to the technical field of reservoir dispatching, in particular to a reservoir flood control dispatching method, system, equipment and medium based on improved POA.
Background
China has few water resources per capita, the rainfall space-time distribution is seriously uneven, the water resource utilization difficulty is high, the water shortage is serious, and the flood and drought disasters are prominent problems in the flood management and flood prevention scheduling work in China. Reservoir engineering is an important engineering means for regulating runoff space-time distribution by human beings, and the amount of discharged water is scientifically controlled by jointly scheduling reservoir groups, so that the flood control pressure of middle and downstream areas is remarkably reduced. The river has been built into the river basin flood control project system mainly including dikes, flood storage areas and main and branch reservoirs, the whole flood control capacity of the river basin is greatly improved, and the river basin flood control project system plays an important role in defending the flood and the extra flood all the time. Meanwhile, with the rapid development of economic society in China, the continuous increase of population and the influence of global climate change, the frequent occurrence of extreme rainfall events, the increase of drainage basin extreme flood events, and the new problems and challenges of drainage basin flood control, the development of library (group) joint scheduling or drainage basin flood control joint scheduling is a necessary choice.
The river basin flood control combined dispatching refers to unified flood control coordinated dispatching for a group of reservoirs and related engineering facilities which are in hydrological, hydraulic and water conservancy connection with each other in a river basin, and the flood loss of the whole river basin is reduced to the minimum while the benefit requirements are considered, and the flood control benefit is maximized. Although the basin flood control combined dispatching is based on the theory and the method of single reservoir flood control dispatching, the flood control requirements of all protection objects of the system are required to be met through the combined regulation and control of all flood control projects according to the flood control standards of reservoirs, upstream and downstream protection objects and the safe discharge of riverways of flood control points. The flood control in different areas is different in flood combination, the flood control engineering in different areas is required to bear different flood control functions, the problem of river basin flood control combined dispatching is far more complicated than that of single reservoir flood control dispatching, but if unified combined flood control dispatching is not implemented, contradictions can be generated inevitably when flood control tasks are shared between upstream and downstream branches, the phenomenon of flood aggravation occurs, and flood control and water ecological environment safety are affected. Therefore, a suitable flood control combined scheduling scheme is constructed according to the flood control operation characteristics of the drainage basin and the actual demands of interest, a scientific and reasonable scheduling operation mode is formed, and the method has important significance for guiding safe and stable reservoir and optimizing operation. It is more obvious necessary and urgent to implement the river basin flood control joint scheduling under the current new economic situation.
Disclosure of Invention
Technical problem to be solved
In view of the above disadvantages and shortcomings of the prior art, the invention provides a reservoir flood control scheduling method, system, device and medium based on improved POA, which solves the technical problems of excessive dependence on initial solutions, complex steps in solution and low solution efficiency in the existing reservoir scheduling optimization method.
(II) technical scheme
In order to achieve the purpose, the invention adopts the main technical scheme that:
in a first aspect, an embodiment of the present invention provides a reservoir flood control scheduling method based on improved POA, which includes:
s1, determining the adjusted water level of each reservoir in the river according to the flood season water level limit of each reservoir;
s2, solving according to the water balance equation and the constraint conditioniIn the first periodnTime interval discharge interval of each reservoir;
s3 based oniIn the first periodnThe time interval warehousing flow in each reservoir scheduling period is calculated, and the time interval average warehousing flow is calculated;
s4, judgmentnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge flow interval or not is determined according to the judgment resultnInitial solution of each time interval discharge of each reservoir;
s5, according to thenStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
Optionally, the water balance equation is:
Figure 326046DEST_PATH_IMAGE001
,
wherein the content of the first and second substances,
Figure 774345DEST_PATH_IMAGE002
is composed ofiIn the first periodnThe discharge capacity of each reservoir is controlled by the flow rate of the water,
Figure 288503DEST_PATH_IMAGE003
for giving upiIn the first periodnThe storage flow of each reservoir is measured,
Figure 214871DEST_PATH_IMAGE004
is composed ofiIn the first periodnThe water storage capacity of each water reservoir is,
Figure 825981DEST_PATH_IMAGE005
is composed ofiLast of a periodnThe water storage capacity of each water reservoir is,
Figure 179602DEST_PATH_IMAGE006
in order to convert the coefficients of the image,
Figure 977793DEST_PATH_IMAGE007
is composed ofiThe time length of the time period is in seconds.
Optionally, the constraint conditions comprise a flood discharge capacity constraint and a reservoir water storage capacity constraint;
the flood discharge capacity constraints are:
Figure 911114DEST_PATH_IMAGE008
the reservoir water storage capacity constraint is as follows:
Figure 111151DEST_PATH_IMAGE009
wherein the content of the first and second substances,
Figure 901253DEST_PATH_IMAGE010
is composed ofiIn the first periodnThe maximum flood discharge flow of each reservoir,
Figure 186741DEST_PATH_IMAGE011
is composed ofiIn the first periodnThe water storage capacity of each water reservoir is,
Figure 658173DEST_PATH_IMAGE012
is composed ofiIn the first periodnThe dead water level of each reservoir is determined,
Figure 712717DEST_PATH_IMAGE013
is composed ofiIn the first periodnCheck flood levels of individual reservoirs.
Optionally, step S2 includes:
s21, using the water level storage capacity curvenThe characteristic water level of each reservoir is converted into corresponding water storage capacity; the water level reservoir capacity curve is a curve representing the relation between the water level of the reservoir and the corresponding reservoir capacity, and the characteristic water level comprises a dead water level and a check flood level;
s22 for meeting the reservoir water storage capacity constraintiIn the first periodnThe water storage capacity of each reservoir is subjected to the following discrete treatment:
Figure 670790DEST_PATH_IMAGE014
wherein the content of the first and second substances,Sin the form of discrete parts, the amount of the particles,
Figure 709153DEST_PATH_IMAGE015
Nthe number of the reservoirs is the number of the reservoirs,
Figure 515435DEST_PATH_IMAGE016
Figure 893327DEST_PATH_IMAGE017
is as followsnThe water storage amount corresponding to the dead water level of each reservoir,
Figure 759652DEST_PATH_IMAGE018
is as followsnThe individual reservoir checks the water storage capacity corresponding to the flood level,
Figure 285311DEST_PATH_IMAGE019
is composed ofiIn the first periodnDiscrete storage capacity of individual reservoirs;
s23, mixing a plurality ofiIn the first periodnDiscrete water storage capacity of individual reservoir
Figure 895284DEST_PATH_IMAGE020
Substituting into the water quantity balance equation to obtainNAniIn the first periodnDischarge of individual reservoir
Figure 127682DEST_PATH_IMAGE021
FromNDischarge volume
Figure 430488DEST_PATH_IMAGE021
Selecting a maximum value that meets said flood discharge capacity constraint
Figure 177864DEST_PATH_IMAGE010
And minimum value
Figure 794790DEST_PATH_IMAGE022
And then determineiIn the first periodnDischarge interval of individual reservoir
Figure 678432DEST_PATH_IMAGE023
Optionally, in step S3, the time-interval average warehousing traffic
Figure 417718DEST_PATH_IMAGE024
The following formula is used to obtain:
Figure 386811DEST_PATH_IMAGE025
Tis the total number of periods within the scheduling period.
Optionally, step S4 includes:
s41, judgmentnTime interval average warehousing flow of individual reservoir
Figure 338587DEST_PATH_IMAGE026
Whether or not to discharge the flow rate interval in time period
Figure 545577DEST_PATH_IMAGE023
Within;
s42a, ifnTime interval average warehousing flow of individual reservoir
Figure 190185DEST_PATH_IMAGE026
Interval of time interval discharge
Figure 912153DEST_PATH_IMAGE023
In, then it isnInitial solution of each time interval discharge of individual reservoir
Figure 667620DEST_PATH_IMAGE027
Time-interval-based average warehousing flow
Figure 463537DEST_PATH_IMAGE026
A value of (d);
s42b, ifnTime interval average warehousing flow of individual reservoir
Figure 279047DEST_PATH_IMAGE026
Interval of discharge quantity not in time period
Figure 491241DEST_PATH_IMAGE023
In, then it isnInitial solution of each time interval discharge of individual reservoir
Figure 50398DEST_PATH_IMAGE028
Interval of time interval and discharge
Figure 966402DEST_PATH_IMAGE023
Inner most proximal tonAverage warehousing flow of individual reservoir time interval
Figure 952812DEST_PATH_IMAGE026
The value of (c).
Optionally, step S5 includes:
s51, using the water level storage capacity curvenStarting and regulating water level of individual reservoir to be converted into reservoirThe water storage amount of the 1 st period of time, thereafter according to the firstnThe water storage capacity of the 1 st period of the individual reservoirnDetermining an initial scheduling line by the initial solution of the reservoir discharge flow of each period of each reservoir and the period warehousing flow;
s52, after determining the initial dispatching line of the reservoirs, according to the discrete storage amount obtained in the step S22, keeping the discrete storage amount of each reservoir
Figure 383793DEST_PATH_IMAGE029
Under the condition that the values of two adjacent time intervals are not changed, the water storage capacity in each discrete state is calculated
Figure 684325DEST_PATH_IMAGE029
Comparing the target values in each discrete state, and selecting the optimal water storage amount from the target values; the optimal water storage capacity of each time interval forms a water storage capacity value sequence;
s53, according to the water storage quantity value sequence and the water balance equation, combiningiIn the first periodnAnd (3) solving a total target value of a dispatching period by the inflow of intervals between each reservoir and the flood control point of the reservoir:
Figure 517152DEST_PATH_IMAGE030
wherein the content of the first and second substances,
Figure 674463DEST_PATH_IMAGE031
is composed ofiIn the first periodnThe interval between each reservoir and its flood control point is entered,Nthe number of the reservoirs is the number of the reservoirs,
Figure 592741DEST_PATH_IMAGE016
,
Figure 696963DEST_PATH_IMAGE032
s54, recording the water storage amount corresponding to the scheduling period total target value of each period after the scheduling period total target value is converged, and then connecting the water storage amounts corresponding to the scheduling period total target value of each period in turn to obtain the secondnOf individual reservoirsAnd scheduling the lines.
In a second aspect, an embodiment of the present invention provides an improved POA-based cascade reservoir flood control optimal scheduling system, which includes:
the starting water level adjusting determining module is used for determining the starting water level of each reservoir in the flow field according to the flood season limiting water level of each reservoir;
a time interval discharge rate interval calculation module for calculating according to the water balance equation and the constraint conditioniIn the first periodnTime interval discharge interval of each reservoir;
an average warehousing flow calculation module based oniIn the first periodnThe time interval warehousing flow in each reservoir scheduling period is calculated, and the time interval average warehousing flow is calculated;
a judging module for judgingnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge interval or not;
an initial solution determining module for determining the first solution according to the judgment resultnInitial solution of each time interval discharge of each reservoir;
a dispatch line calculation module for calculating a dispatch line based onnStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
In a third aspect, an embodiment of the present invention provides an electronic device, including:
a processor;
and the storage is used for controlling the steps of the reservoir flood control scheduling method based on the improved POA by the processor.
In a fourth aspect, embodiments of the present invention provide a computer-readable medium, on which computer-executable instructions are stored, and when executed by a processor, the steps of the method for reservoir flood control based on improved POA are implemented.
(III) advantageous effects
The invention has the beneficial effects that: in the optimized dispatching of the cascade reservoir, the average warehousing flow in each reservoir dispatching period is firstly solved, then an initial dispatching line is obtained, the requirement of a maximum peak clipping criterion objective function is better met, meanwhile, a step-by-step optimization method (POA) applied later can quickly converge, and finally, a reasonable flood control optimized dispatching scheme is obtained. Compared with the prior art, the technical scheme of the invention has the advantages that the solving steps are simpler and clearer, the understanding is easy, the solving process of the initial solution is simplified, the initial solution meeting the conditions can be quickly generated, and the solving efficiency of the initial solution is improved.
Drawings
Fig. 1 is a schematic flow chart of a reservoir flood control scheduling method based on improved POA according to the present invention;
fig. 2 is a schematic flow chart of step S2 of the reservoir flood control scheduling method based on improved POA according to the present invention;
fig. 3 is a schematic flow chart of step S4 of the reservoir flood control scheduling method based on improved POA according to the present invention;
fig. 4 is a schematic flow chart of step S5 of the reservoir flood control scheduling method based on improved POA according to the present invention;
fig. 5 is a schematic diagram illustrating the reservoir 1 and reservoir 2 flood control scheduling effect comparison of the reservoir 1 and reservoir 2 flood control scheduling method based on the improved POA;
fig. 6 is a schematic composition diagram of a cascade reservoir flood control optimized dispatching system based on improved POA provided by the invention;
FIG. 7 is a schematic structural diagram of a computer system of an electronic device according to the present invention;
fig. 8 is an algorithm flowchart of a reservoir flood control scheduling method based on improved POA according to the present invention.
[ description of reference ]
300: reservoir flood control optimizing and dispatching system; 301: a starting water level adjusting determining module; 302: a time interval discharge rate interval calculation module; 303: an average warehousing flow calculation module; 304: a judgment module; 305: an initial solution determination module; 306: a dispatch line calculation module;
400: a computer system; 401: a CPU; 402: a ROM; 403: a RAM; 404: a bus; 405: an I/O interface; 406: an input section; 407: an output section; 408: a storage section; 409: a communication section; 410: a driver; 411: a removable media.
Detailed Description
For the purpose of better explaining the present invention and to facilitate understanding, the present invention will be described in detail by way of specific embodiments with reference to the accompanying drawings.
Fig. 1 is a schematic flow chart of a reservoir flood control scheduling method based on an improved POA according to the present invention, and as shown in fig. 1, a reservoir flood control scheduling method based on an improved POA according to an embodiment of the present invention includes: firstly, determining the starting water level of each reservoir in a flow area; secondly, according to the water quantity balance equation and the constraint condition, the method calculatesiIn the first periodnTime interval discharge interval of each reservoir; then, calculate the firstnAverage warehousing flow of time intervals in each reservoir scheduling period; then, judgenWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge flow interval or not is determined according to the judgment resultnInitial solution of each time interval discharge of each reservoir; finally, according tonStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
In the optimized dispatching of the cascade reservoir, the average warehousing flow in each reservoir dispatching period is firstly solved, then an initial dispatching line is obtained, the requirement of a maximum peak clipping criterion objective function is better met, meanwhile, a step-by-step optimization method (POA) applied later can quickly converge, and finally, a reasonable flood control optimized dispatching scheme is obtained. Compared with the prior art, the technical scheme of the invention has the advantages that the solving steps are simpler and clearer, the understanding is easy, the solving process of the initial solution is simplified, the initial solution meeting the conditions can be quickly generated, and the solving efficiency of the initial solution is improved.
For a better understanding of the above-described technical solutions, exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention can be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Specifically, the invention provides a reservoir flood control scheduling method based on improved POA, which comprises the following steps:
and S1, determining the adjusted water level of each reservoir in the river according to the flood season water level limit of each reservoir. Known asnThe flood season of each reservoir is limited to water level
Figure 384296DEST_PATH_IMAGE033
Of 1 atnStarting and regulating water level of individual reservoir
Figure 712510DEST_PATH_IMAGE034
Is lower than or equal to the flood season limiting water level (determined according to the running condition of the reservoir). The initial water level of the first stage of each reservoir is the water level
Figure 55766DEST_PATH_IMAGE035
The last stage is the control water level
Figure 760417DEST_PATH_IMAGE036
Figure 567836DEST_PATH_IMAGE032
Figure 66950DEST_PATH_IMAGE016
NThe number of the reservoirs is the number of the reservoirs,Tis the total number of periods within the scheduling period.
S2, solving according to the water balance equation and the constraint conditioniIn the first periodnThe time interval discharge rate interval of each reservoir.
Fig. 2 is a schematic detailed flowchart of step S2 of the reservoir flood control scheduling method based on improved POA according to the present invention, as shown in fig. 2, step S2 includes:
s21, using the water level storage capacity curvenThe characteristic water level of each reservoir is converted into corresponding water storage capacity; the water level and storage capacity curve is a curve representing the relation between the water level of the reservoir and the corresponding storage capacity, the storage capacity curve is usually drawn by taking the water level as a vertical coordinate and the storage capacity as a horizontal coordinate, and the mutual conversion between the water level of the reservoir and the storage capacity can be realized according to the storage capacity curve of the reservoir. The characteristic water level comprises a dead water level and a check flood level.
S22 method for meeting reservoir water storage capacity constraintiIn the first periodnThe water storage capacity of each reservoir is subjected to the following discrete treatment:
Figure 897503DEST_PATH_IMAGE037
wherein the content of the first and second substances,Sin the form of discrete parts, the amount of the particles,
Figure 405845DEST_PATH_IMAGE015
the first of each periodnWater storage capacity of individual reservoir
Figure 802191DEST_PATH_IMAGE038
All get through under all conditions
Figure 675469DEST_PATH_IMAGE039
The value of the one or more of,Nthe number of the reservoirs is the number of the reservoirs,
Figure 790056DEST_PATH_IMAGE016
Figure 836509DEST_PATH_IMAGE040
is as followsnThe water storage amount corresponding to the dead water level of each reservoir,
Figure 556204DEST_PATH_IMAGE041
is as followsnThe individual reservoir checks the water storage capacity corresponding to the flood level,
Figure 394191DEST_PATH_IMAGE020
is composed ofiIn the first periodnDiscrete storage capacity of individual reservoirs.
S23, mixing a plurality ofiIn the first periodnDiscrete water storage capacity of individual reservoir
Figure 996074DEST_PATH_IMAGE042
Substituting into the water quantity balance equation to obtainNAniIn the first periodnDischarge of individual reservoir
Figure 580639DEST_PATH_IMAGE043
FromNDischarge volume
Figure 217156DEST_PATH_IMAGE044
Selecting the maximum value meeting the flood discharge capacity constraint
Figure 432237DEST_PATH_IMAGE045
And minimum value
Figure 255837DEST_PATH_IMAGE046
And then determineiIn the first periodnDischarge interval of individual reservoir
Figure 909672DEST_PATH_IMAGE023
Further, the water balance equation is:
Figure 135117DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 521099DEST_PATH_IMAGE047
is composed ofiIn the first periodnThe discharge capacity of each reservoir and the length of each time interval are determined according to actual needs, and can be 1h, 3h, 1 day, 1 ten days or 1 month and the like,
Figure 363153DEST_PATH_IMAGE003
for giving upiIn the first periodnThe storage flow of each reservoir is measured,
Figure 289521DEST_PATH_IMAGE011
is composed ofiIn the first periodnThe water storage capacity of each water reservoir is,
Figure 572734DEST_PATH_IMAGE048
is composed ofiLast of a periodnThe water storage capacity of each water reservoir is,
Figure 191934DEST_PATH_IMAGE006
for conversion factors (e.g. storage capacity units of large reservoirs are usually in billions of meters)3Flow rate in m3S, at which time the conversion factor may be determined to be 100000000),
Figure 990126DEST_PATH_IMAGE049
is composed ofiThe time length of the time period is in seconds.
Furthermore, the constraint conditions comprise flood discharge capacity constraint, reservoir water storage capacity constraint, water level constraint in a dispatching period and ex-warehouse flow amplitude constraint.
The flood discharge capacity constraints are:
Figure 720185DEST_PATH_IMAGE008
the reservoir water storage capacity constraint is as follows:
Figure 920222DEST_PATH_IMAGE009
the water level constraint in the scheduling period is as follows:
Figure 710323DEST_PATH_IMAGE050
the variation of amplitude of flow out of the warehouse is restricted as follows:
Figure 199074DEST_PATH_IMAGE051
wherein the content of the first and second substances,
Figure 467244DEST_PATH_IMAGE052
is composed ofiIn the first periodnOne reservoir to the greatestThe flow rate of the large flood discharge is high,
Figure 787367DEST_PATH_IMAGE038
is composed ofiIn the first periodnThe water storage capacity of each water reservoir is,
Figure 751299DEST_PATH_IMAGE053
is composed ofiIn the first periodnThe dead water level of each reservoir is determined,
Figure 461766DEST_PATH_IMAGE054
is composed ofiIn the first periodnThe flood level of each reservoir is checked,
Figure 533628DEST_PATH_IMAGE055
is as followsnThe water level at the end of the dispatching period of each reservoir,
Figure 973836DEST_PATH_IMAGE056
is as followsnThe water level of each reservoir is finally controlled,
Figure 840161DEST_PATH_IMAGE057
is composed ofiLast of a periodnThe discharge capacity of each reservoir is controlled by the flow rate of the water,
Figure 100241DEST_PATH_IMAGE058
is as followsnThe maximum amplitude of variation allowed by the delivery flow of each reservoir in adjacent time intervals. The scheduling end-of-period water level constraint and ex-warehouse flow amplitude variation constraint are used in the subsequent application of a POA algorithm, and each kind of discrete constraint is usediIn the first periodnWater storage capacity of individual reservoir
Figure 913476DEST_PATH_IMAGE020
Calculating the target value, and if any constraint condition is not satisfied during the calculation process, discarding the state value
Figure 208191DEST_PATH_IMAGE020
S3 based oniIn the first periodnAnd (4) calculating the average warehousing flow of each time period in the scheduling period of each reservoir.
Further, in step S3, the time-interval average warehousing traffic is
Figure 245417DEST_PATH_IMAGE059
The following formula is used to obtain:
Figure 196056DEST_PATH_IMAGE060
Tis the total number of periods within the scheduling period.
S4, judgmentnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge flow interval or not is determined according to the judgment resultnInitial solution of discharge of each reservoir in each period.
Fig. 3 is a schematic specific flowchart of step S4 of the reservoir flood control scheduling method based on improved POA according to the present invention, as shown in fig. 3, step S4 includes:
s41, judgmentnTime interval average warehousing flow of individual reservoir
Figure 609720DEST_PATH_IMAGE026
Whether or not to discharge the flow rate interval in time period
Figure 758941DEST_PATH_IMAGE023
Within;
s42a, ifnTime interval average warehousing flow of individual reservoir
Figure 232648DEST_PATH_IMAGE026
Interval of time interval discharge
Figure 405003DEST_PATH_IMAGE023
In, then it isnInitial solution of each time interval discharge of individual reservoir
Figure 622358DEST_PATH_IMAGE061
Time-interval-based average warehousing flow
Figure 626086DEST_PATH_IMAGE026
A value of (d);
s42b, ifnTime interval average warehousing flow of individual reservoir
Figure 5115DEST_PATH_IMAGE026
Interval of discharge quantity not in time period
Figure 664766DEST_PATH_IMAGE023
In, then it isnInitial solution of each time interval discharge of individual reservoir
Figure 685812DEST_PATH_IMAGE061
Interval of time interval and discharge
Figure 544047DEST_PATH_IMAGE023
Inner most proximal tonAverage warehousing flow of individual reservoir time interval
Figure 562818DEST_PATH_IMAGE026
The value of (c).
S5, according to thenStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
Fig. 4 is a schematic detailed flowchart of step S5 of the reservoir flood control scheduling method based on improved POA according to the present invention, as shown in fig. 4, step S5 includes:
s51, using the water level storage capacity curvenThe starting water level of the individual reservoir is converted into the water storage capacity of the reservoir in the 1 st period of time and then according to the 1 st period of timenThe water storage capacity of the 1 st period of the individual reservoirnAnd determining an initial scheduling line by the initial solution of the reservoir discharge flow of each period of each reservoir and the period warehousing flow.
S52, after determining the initial dispatching line of the reservoirs, according to the discrete storage amount obtained in the step S22, keeping the discrete storage amount of each reservoir
Figure 506503DEST_PATH_IMAGE062
Two adjacent time periodsUnder the condition that the value of (A) is not changed, the water storage amount in each discrete state is calculated
Figure 62731DEST_PATH_IMAGE062
Comparing the target values in each discrete state, and selecting the optimal water storage amount from the target values; the optimal water storage amount of each time interval forms a water storage amount value sequence. Preferably, the following is also applicable: in maintaining the quantity of water dispersed from each reservoir
Figure 775472DEST_PATH_IMAGE039
Under the condition that the values of other time intervals are not changed, only the first time in the initial scheduling line is changed every timenA water reservoiriThe water storage amount of each discrete state is calculated
Figure 965145DEST_PATH_IMAGE062
And comparing the target values in each discrete state to select the optimal water storage amount.
In the above step, hold
Figure 396126DEST_PATH_IMAGE063
And
Figure 493395DEST_PATH_IMAGE064
unchanged, calculating the second after dispersionnValue of individual reservoir at time 1
Figure 326222DEST_PATH_IMAGE065
Comparing the target values in each discrete state to select the optimal storage capacity of the reservoir
Figure 686796DEST_PATH_IMAGE066
To make
Figure 605074DEST_PATH_IMAGE067
(ii) a Then hold
Figure 506034DEST_PATH_IMAGE068
And
Figure 193367DEST_PATH_IMAGE069
unchanged, by discretenValue of individual reservoir at time 2
Figure 724842DEST_PATH_IMAGE070
Calculating at each discrete state
Figure 130416DEST_PATH_IMAGE071
Comparing the target values of each state to select the optimal storage capacity of the reservoir
Figure 835067DEST_PATH_IMAGE072
To make
Figure 376907DEST_PATH_IMAGE073
. And by analogy, the optimal water storage capacity of each time interval is solved and a water storage capacity value sequence is formed.
S53, according to the water storage quantity value sequence and the water balance equation, combiningiIn the first periodnAnd (3) solving a total target value of a dispatching period by the inflow of intervals between each reservoir and the flood control point of the reservoir:
Figure 876021DEST_PATH_IMAGE030
wherein the content of the first and second substances,
Figure 706574DEST_PATH_IMAGE021
is composed ofiIn the first periodnThe discharge capacity of each reservoir is controlled by the flow control valve,
Figure 949336DEST_PATH_IMAGE031
is composed ofiIn the first periodnThe interval between each reservoir and its flood control point is entered,Nthe number of the reservoirs is the number of the reservoirs,
Figure 611262DEST_PATH_IMAGE016
,
Figure 484540DEST_PATH_IMAGE032
s54, recording the water storage amount corresponding to the scheduling period total target value of each period after the scheduling period total target value is converged, and then connecting the water storage amounts corresponding to the scheduling period total target value of each period in turn to obtain the secondnA dispatch line for each reservoir. The criterion for convergence is that the target value is stable after several program calculations.
From the equation of water balance, it is known
Figure 864706DEST_PATH_IMAGE074
Figure 645580DEST_PATH_IMAGE075
Figure 365274DEST_PATH_IMAGE076
Figure 209121DEST_PATH_IMAGE077
In combination with the firstnThe water storage capacity of the reservoir corresponding to the starting water level of the individual reservoir can be obtainednVariation process of water storage capacity of individual reservoir
Figure 811004DEST_PATH_IMAGE078
By using the firstnThe curve drawn by the water storage capacity point corresponding to the total target value of the dispatching period of each period of the reservoir is the firstnA dispatch line for each reservoir.
In a specific embodiment, a certain Dongting lake basin is characterized in that the coastal area is provided with more Chongshan and steep mountains and plateaus, the slope is large, the average slope of a river channel is 0.254 per thousand, the valley is more, the beach risk is more, the water flow is turbulent, the downstream area is supported by the top of the Dongting lake, the flood control situation is complex, the encountering probability of the flood of the basin and the flood of the adjacent basin is higher, and the combined flood control scheduling of the downstream reservoir 1 and the cascade reservoir 2 in the basin is necessary for realizing scientific and effective peak staggering and optimal water quantity scheduling of the cascade reservoir. The characteristic parameters of the reservoir 1 and the reservoir 2 are shown in table 1. In the research, the design flood process of 50-year warehousing of two reservoirs is selected as the input condition of flood control combined optimization scheduling calculation.
TABLE 1 characteristic parameters table of reservoir
Figure DEST_PATH_IMAGE080
The invention uses an improved stepwise optimization method (POA) to carry out the combined flood control optimization scheduling calculation on two step reservoirs in a certain sub-watershed of the Dongting lake, and the calculation results are shown in the table 2 and the figure 5.
TABLE 2 flood control scheduling computation results
Water reservoir Initial water level/m of dispatching period Scheduling end of term water level/m Maximum water level limit/m Maximum amplitude/(m) allowed for delivery3.s-1) Peak reduction rate/%)
Reservoir 1 198.5 198.5 209.56 2000 42
Reservoir 2 108 108 114.53 8000 20
In the above method, the present invention provides a method for generating an initial solution of POA based on the maximum peak clipping criterion in the flood control scheduling problem: the average warehousing flow in each reservoir scheduling period is firstly solved, and then an initial scheduling line is obtained, so that the requirement of a maximum peak clipping criterion target function is better met. Meanwhile, the size of the required flood regulation storage capacity does not need to be considered in the solving process of the initial solution, and only the maximum flood regulation storage capacity corresponding to the check flood level needs to be considered; the reservoir dispatching line does not need to be adjusted repeatedly, the solving process of the initial solution is simplified, the solving efficiency of the initial solution is improved, and the scheme is simple and clear.
In addition, fig. 6 is a schematic composition diagram of a cascade reservoir flood control optimized dispatching system based on the improved POA, as shown in fig. 6, the invention further provides a cascade reservoir flood control optimized dispatching system based on the improved POA, which includes:
and the starting water level adjusting determining module 301 is used for determining the starting water level of each reservoir in the flow field according to the flood season water level limit of each reservoir.
A time interval discharge rate interval calculation module 302 for calculating the water balance equation and the constraint conditioniIn the first periodnThe time interval discharge rate interval of each reservoir.
An average warehousing traffic calculation module 303 for calculating a traffic based oniIn the first periodnAnd (4) calculating the average warehousing flow of each time period in the scheduling period of each reservoir.
A judging module 304 for judgingnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge interval or not.
An initial solution determining module 305 for determining the second solution according to the determination resultnInitial solution of discharge of each reservoir in each period.
A dispatch line calculation module 306 fornStarting to adjust water level of individual reservoirnReservoir discharge of individual reservoir at each time intervalThe initial solution of the quantity and the time interval warehousing flow are obtainednThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
The functions of each module in the system are described in the above method embodiments, and are not described herein again.
In another aspect, the present invention further provides an electronic device, which includes a processor and a memory, wherein the memory stores the steps of the method for analyzing the cause of the engine abnormality, which is used for controlling the processor.
Referring now to FIG. 7, shown is a block diagram of a computer system 400 suitable for use in implementing the electronic device of an embodiment of the present application. The electronic device shown in fig. 7 is only an example, and should not bring any limitation to the functions and the scope of use of the embodiments of the present application.
As shown in fig. 7, the computer system 400 includes a Central Processing Unit (CPU) 401 that can perform various appropriate actions and processes in accordance with a program stored in a Read Only Memory (ROM) 402 or a program loaded from a storage section 408 into a Random Access Memory (RAM) 403. In the RAM 403, various programs and data necessary for the operation of the system 400 are also stored. The CPU 401, ROM 402, and RAM 403 are connected to each other via a bus 404. An input/output (I/O) interface 405 is also connected to bus 404.
The following components are connected to the I/O interface 405: an input section 406 including a keyboard, a mouse, and the like; an output section 407 including a display device such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, and a speaker; a storage section 408 including a hard disk and the like; and a communication section 409 including a network interface card such as a LAN card, a modem, or the like. The communication section 409 performs communication processing via a network such as the internet. A driver 410 is also connected to the I/O interface 405 as needed. A removable medium 411 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 410 as necessary, so that a computer program read out therefrom is mounted into the storage section 408 as necessary.
In particular, according to an embodiment of the present invention, the processes described above with reference to the flowcharts may be implemented as computer software programs. For example, embodiments of the invention include a computer program product comprising a computer program embodied on a computer-readable medium, the computer program comprising program code for performing the method illustrated in the flow chart. In such an embodiment, the computer program may be downloaded and installed from a network through the communication section 409, and/or installed from the removable medium 411. The above-described functions defined in the system of the present application are executed when the computer program is executed by a Central Processing Unit (CPU) 401.
It should be noted that the computer readable medium shown in the present application may be a computer readable signal medium or a computer readable medium or any combination of the two. A computer readable medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples of the computer readable medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the present application, a computer readable medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In this application, however, a computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: wireless, wire, fiber optic cable, RF, etc., or any suitable combination of the foregoing.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams or flowchart illustration, and combinations of blocks in the block diagrams or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units described in the embodiments of the present application may be implemented by software or hardware. The units described may also be provided in a processor, where the names of the units do not in some cases constitute a limitation of the units themselves.
In another aspect, the present invention also provides a computer-readable medium, which may be contained in the apparatus described in the above embodiments; or may be separate and not incorporated into the device. The computer readable medium carries one or more programs which, when executed by a device, cause the device to include the method steps of:
and S1, determining the adjusted water level of each reservoir in the river according to the flood season water level limit of each reservoir.
S2, solving according to the water balance equation and the constraint conditioniIn the first periodnThe time interval discharge rate interval of each reservoir.
S3 based oniIn the first periodnAnd (4) calculating the average warehousing flow of each time period in the scheduling period of each reservoir.
S4, judgmentnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge flow interval or not is determined according to the judgment resultnInitial solution of discharge of each reservoir in each period.
S5, according to thenStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
In summary, the invention provides a reservoir flood control scheduling method, system, device and medium based on improved POA. Fig. 8 is an algorithm flow chart of the reservoir flood control scheduling method based on the improved POA provided by the present invention, and as shown in fig. 8, the algorithm flow of the technical scheme of the present invention is as follows: firstly, determining the starting water level of each reservoir, and secondly, determining the time interval discharge interval of each reservoir; then, judging whether the time-interval average leakage flow is within a time-interval leakage flow interval or not, if so, obtaining a time-interval average leakage flow value by the initial solution of the time-interval leakage flow, and if not, giving the time-interval average leakage flow value which is the closest to the time-interval average leakage flow value within the time-interval leakage flow interval to the initial solution of the time-interval leakage flow; then, after the value of the initial solution of the discharge in the time period is determined, the water storage capacity of the reservoir in the current time period is obtained according to the water balance equation in sequence, the water storage capacity of two adjacent time periods is fixed, and the water storage capacity in the current time period is optimized based on the water balance equation; and obtaining a total target value of the scheduling period, and converging the total target value.
Compared with the background technology, the technical scheme of the invention has simpler and clearer solving steps, is easy to understand, can generate an initial solution meeting the conditions, and can quickly converge by a gradual optimization method (POA) in the optimal dispatching of the cascade reservoir to obtain a satisfactory and reasonable optimal dispatching scheme for flood control.
Since the system/apparatus described in the above embodiments of the present invention is a system/apparatus used for implementing the method of the above embodiments of the present invention, a person skilled in the art can understand the specific structure and modification of the system/apparatus based on the method described in the above embodiments of the present invention, and thus the detailed description is omitted here. All systems/devices adopted by the methods of the above embodiments of the present invention are within the intended scope of the present invention.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions.
It should be noted that in the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the terms first, second, third and the like are for convenience only and do not denote any order. These words are to be understood as part of the name of the component.
Furthermore, it should be noted that in the description of the present specification, the description of the term "one embodiment", "some embodiments", "examples", "specific examples" or "some examples", etc., means that a specific feature, structure, material or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, the claims should be construed to include preferred embodiments and all changes and modifications that fall within the scope of the invention.
It will be apparent to those skilled in the art that various modifications and variations can be made in the present invention without departing from the spirit or scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention should also include such modifications and variations.

Claims (10)

1. A reservoir flood control dispatching method based on improved POA is characterized by comprising the following steps:
s1, determining the adjusted water level of each reservoir in the river according to the flood season water level limit of each reservoir;
s2, solving according to the water balance equation and the constraint conditioniIn the first periodnTime interval discharge interval of each reservoir;
s3 based oniIn the first periodnThe time interval warehousing flow in the individual reservoir scheduling period is calculated, and the time interval average warehousing flow is calculatedAn amount;
s4, judgmentnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge flow interval or not is determined according to the judgment resultnInitial solution of each time interval discharge of each reservoir;
s5, according to thenStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
2. The improved POA-based reservoir flood control scheduling method of claim 1, wherein the water balance equation is as follows:
Figure 59204DEST_PATH_IMAGE001
,
wherein the content of the first and second substances,
Figure 430142DEST_PATH_IMAGE002
is composed ofiIn the first periodnThe discharge capacity of each reservoir is controlled by the flow rate of the water,
Figure 707540DEST_PATH_IMAGE003
for giving upiIn the first periodnThe storage flow of each reservoir is measured,
Figure 531139DEST_PATH_IMAGE004
is composed ofiIn the first periodnThe water storage capacity of each water reservoir is,
Figure 122658DEST_PATH_IMAGE005
is composed ofiLast of a periodnThe water storage capacity of each water reservoir is,
Figure 351032DEST_PATH_IMAGE006
in order to convert the coefficients of the image,
Figure 799331DEST_PATH_IMAGE007
is composed ofiThe time length of the time period is in seconds.
3. The improved POA-based reservoir flood control scheduling method of claim 2, wherein the constraint conditions comprise a flood discharge capacity constraint and a reservoir water storage capacity constraint;
the flood discharge capacity constraints are:
Figure 375806DEST_PATH_IMAGE008
the reservoir water storage capacity constraint is as follows:
Figure 505436DEST_PATH_IMAGE009
wherein the content of the first and second substances,
Figure 585388DEST_PATH_IMAGE010
is composed ofiIn the first periodnThe maximum flood discharge flow of each reservoir,
Figure 204588DEST_PATH_IMAGE011
is composed ofiIn the first periodnThe water storage capacity of each water reservoir is,
Figure 268359DEST_PATH_IMAGE012
is composed ofiIn the first periodnThe dead water level of each reservoir is determined,
Figure 732838DEST_PATH_IMAGE013
is composed ofiIn the first periodnCheck flood levels of individual reservoirs.
4. The improved POA-based reservoir flood control scheduling method of claim 3, wherein the step S2 comprises:
s21, using the water level storage capacity curvenThe characteristic water level of each reservoir is converted into corresponding water storage capacity; the water level reservoir capacity curve is a curve representing the relation between the water level of the reservoir and the corresponding reservoir capacity, and the characteristic water level comprises a dead water level and a check flood level;
s22 for meeting the reservoir water storage capacity constraintiIn the first periodnThe water storage capacity of each reservoir is subjected to the following discrete treatment:
Figure 136138DEST_PATH_IMAGE014
wherein the content of the first and second substances,Sin the form of discrete parts, the amount of the particles,
Figure 926239DEST_PATH_IMAGE015
Nthe number of the reservoirs is the number of the reservoirs,
Figure 211727DEST_PATH_IMAGE016
Figure 745476DEST_PATH_IMAGE017
is as followsnThe water storage amount corresponding to the dead water level of each reservoir,
Figure 3282DEST_PATH_IMAGE018
is as followsnThe individual reservoir checks the water storage capacity corresponding to the flood level,
Figure 698706DEST_PATH_IMAGE019
is composed ofiIn the first periodnDiscrete storage capacity of individual reservoirs;
s23, mixing a plurality ofiIn the first periodnDiscrete water storage capacity of individual reservoir
Figure 737069DEST_PATH_IMAGE020
Substituting into the water quantity balance equation to obtainNAniIn the first periodnDischarge of individual reservoir
Figure 808930DEST_PATH_IMAGE021
FromNA discharge flowMeasurement of
Figure 717981DEST_PATH_IMAGE021
Selecting a maximum value that meets said flood discharge capacity constraint
Figure 53147DEST_PATH_IMAGE022
And minimum value
Figure 47648DEST_PATH_IMAGE023
And then determineiIn the first periodnDischarge interval of individual reservoir
Figure 188779DEST_PATH_IMAGE024
5. The improved POA-based reservoir flood control scheduling method of claim 4, wherein in step S3, the time interval average warehousing flow rate
Figure 214985DEST_PATH_IMAGE025
The following formula is used to obtain:
Figure 455474DEST_PATH_IMAGE026
Tis the total number of periods within the scheduling period.
6. The improved POA-based reservoir flood control scheduling method of claim 5, wherein the step S4 comprises:
s41, judgmentnTime interval average warehousing flow of individual reservoir
Figure 202850DEST_PATH_IMAGE027
Whether or not to discharge the flow rate interval in time period
Figure 882093DEST_PATH_IMAGE024
Within;
s42a, ifnTime interval average warehousing flow of individual reservoir
Figure 765735DEST_PATH_IMAGE027
Interval of time interval discharge
Figure 177125DEST_PATH_IMAGE024
In, then it isnInitial solution of each time interval discharge of individual reservoir
Figure 146218DEST_PATH_IMAGE028
Time-interval-based average warehousing flow
Figure 629152DEST_PATH_IMAGE027
A value of (d);
s42b, ifnTime interval average warehousing flow of individual reservoir
Figure 367301DEST_PATH_IMAGE027
Interval of discharge quantity not in time period
Figure 215171DEST_PATH_IMAGE024
In, then it isnInitial solution of each time interval discharge of individual reservoir
Figure 671560DEST_PATH_IMAGE029
Interval of time interval and discharge
Figure 692606DEST_PATH_IMAGE024
Inner most proximal tonAverage warehousing flow of individual reservoir time interval
Figure 550841DEST_PATH_IMAGE027
The value of (c).
7. The improved POA-based reservoir flood control scheduling method of claim 6, wherein the step S5 comprises:
s51 using the curve of reservoir capacityWill be firstnThe starting water level of the individual reservoir is converted into the water storage capacity of the reservoir in the 1 st period of time and then according to the 1 st period of timenThe water storage capacity of the 1 st period of the individual reservoirnDetermining an initial scheduling line by the initial solution of the reservoir discharge flow of each period of each reservoir and the period warehousing flow;
s52, after determining the initial dispatching line of the reservoirs, according to the discrete storage amount obtained in the step S22, keeping the discrete storage amount of each reservoir
Figure 569612DEST_PATH_IMAGE030
Under the condition that the values of two adjacent time intervals are not changed, the water storage capacity in each discrete state is calculated
Figure 247718DEST_PATH_IMAGE030
Comparing the target values in each discrete state, and selecting the optimal water storage amount from the target values; the optimal water storage capacity of each time interval forms a water storage capacity value sequence;
s53, according to the water storage quantity value sequence and the water balance equation, combiningiIn the first periodnAnd (3) solving a total target value of a dispatching period by the inflow of intervals between each reservoir and the flood control point of the reservoir:
Figure 72455DEST_PATH_IMAGE031
wherein the content of the first and second substances,
Figure 785196DEST_PATH_IMAGE032
is composed ofiIn the first periodnThe interval between each reservoir and its flood control point is entered,Nthe number of the reservoirs is the number of the reservoirs,
Figure 974869DEST_PATH_IMAGE016
,
Figure 405850DEST_PATH_IMAGE033
s54, recording the total dispatching period of each period after the convergence of the total dispatching period target valueThe water storage amount corresponding to the target value is connected in turn with the total target value of the dispatching period of each time interval to obtain the second water storage amountnA dispatch line for each reservoir.
8. A step reservoir flood control optimization and dispatching system based on improved POA is characterized by comprising:
the starting water level adjusting determining module is used for determining the starting water level of each reservoir in the flow field according to the flood season limiting water level of each reservoir;
a time interval discharge rate interval calculation module for calculating according to the water balance equation and the constraint conditioniIn the first periodnTime interval discharge interval of each reservoir;
an average warehousing flow calculation module based oniIn the first periodnThe time interval warehousing flow in each reservoir scheduling period is calculated, and the time interval average warehousing flow is calculated;
a judging module for judgingnWhether the time interval average warehousing flow of each reservoir is within the corresponding time interval discharge interval or not;
an initial solution determining module for determining the first solution according to the judgment resultnInitial solution of each time interval discharge of each reservoir;
a dispatch line calculation module for calculating a dispatch line based onnStarting to adjust water level of individual reservoirnObtaining the initial solution of reservoir discharge and time-interval warehousing flow of each time interval of each reservoirnThe initial dispatching line of each reservoir is determined by solving the total target value of the dispatching period based on the initial dispatching line and the water quantity balance equationnA dispatch line for each reservoir.
9. An electronic device, comprising:
a processor;
a memory for storing instructions for the processor to control the steps of a method for improved POA-based reservoir flood control scheduling according to any of claims 1-7.
10. A computer readable medium having stored thereon computer executable instructions, which when executed by a processor, perform the steps of the method for improved POA-based flood control dispatch for reservoirs of any of claims 1 to 7.
CN202110199694.6A 2021-02-23 2021-02-23 Reservoir flood control scheduling method, system, equipment and medium based on improved POA Active CN112561213B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110199694.6A CN112561213B (en) 2021-02-23 2021-02-23 Reservoir flood control scheduling method, system, equipment and medium based on improved POA

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110199694.6A CN112561213B (en) 2021-02-23 2021-02-23 Reservoir flood control scheduling method, system, equipment and medium based on improved POA

Publications (2)

Publication Number Publication Date
CN112561213A true CN112561213A (en) 2021-03-26
CN112561213B CN112561213B (en) 2021-05-25

Family

ID=75034548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110199694.6A Active CN112561213B (en) 2021-02-23 2021-02-23 Reservoir flood control scheduling method, system, equipment and medium based on improved POA

Country Status (1)

Country Link
CN (1) CN112561213B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113742637A (en) * 2021-08-16 2021-12-03 中国水利水电科学研究院 Method and device for calculating annual average silt loss rate of reservoir, electronic equipment and storage medium
CN113934777A (en) * 2021-12-16 2022-01-14 长江水利委员会水文局 Method and system for quantifying influence of backwater jacking on water level change
CN114118554A (en) * 2021-11-18 2022-03-01 武汉大学 Typhoon path-based dynamic control method and system for operating water level of reservoir in flood season

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109886543A (en) * 2019-01-16 2019-06-14 河南大学 Reservoir optimizing and dispatching method based on uniform design
CN112232659A (en) * 2020-10-15 2021-01-15 华中科技大学 Cascade reservoir power generation scheduling method and system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109886543A (en) * 2019-01-16 2019-06-14 河南大学 Reservoir optimizing and dispatching method based on uniform design
CN112232659A (en) * 2020-10-15 2021-01-15 华中科技大学 Cascade reservoir power generation scheduling method and system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李雨: "水库防洪和蓄水优化调度方法及应用", 《中国博士学位论文全文数据库工程科技II辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113742637A (en) * 2021-08-16 2021-12-03 中国水利水电科学研究院 Method and device for calculating annual average silt loss rate of reservoir, electronic equipment and storage medium
CN113742637B (en) * 2021-08-16 2024-05-07 中国水利水电科学研究院 Calculation method and device for annual average silt loss rate of reservoir, electronic equipment and storage medium
CN114118554A (en) * 2021-11-18 2022-03-01 武汉大学 Typhoon path-based dynamic control method and system for operating water level of reservoir in flood season
CN114118554B (en) * 2021-11-18 2024-05-14 武汉大学 Dynamic control method and system for reservoir flood period running water level based on typhoon path
CN113934777A (en) * 2021-12-16 2022-01-14 长江水利委员会水文局 Method and system for quantifying influence of backwater jacking on water level change

Also Published As

Publication number Publication date
CN112561213B (en) 2021-05-25

Similar Documents

Publication Publication Date Title
CN112561213B (en) Reservoir flood control scheduling method, system, equipment and medium based on improved POA
WO2021109848A1 (en) Ecological flow determination method considering lift amount
Azari et al. Multi-objective optimization of the reservoir system operation by using the hedging policy
CN112633578A (en) Optimized dispatching method for lower-grade reservoir group under influence of diversion project
CN110555239B (en) River surge sluice optimization scheduling method based on drainage pipe network model
CN112149984B (en) Reservoir flood regulation multidimensional uncertainty risk analysis method based on Bayesian network
CN111126847A (en) Cascade reservoir short-term optimization scheduling method and system coupled with riverway water power process
WO2021238170A1 (en) Method and system for calculation and transaction of binary water rights
CN112989538A (en) Control method and control device for urban drainage system
CN104636834B (en) A kind of improved joint probability plan model system optimization method
CN107038151A (en) The step Optimized Operation drawing drawing method and dispatching method controlled based on accumulation of energy
Chen et al. Applying genetic algorithm and neural network to the conjunctive use of surface and subsurface water
Kangrang et al. Optimal Reservoir Rule Curves Considering Conditional Ant Colony Optimization with
CN106779468A (en) A kind of user power utilization demand response dynamic modeling and the uncertain appraisal procedure of response
CN113050430A (en) Drainage system control method based on robust reinforcement learning
Gupta et al. Optimization techniques for leakage management in urban water distribution networks
CN110334851B (en) Hybrid cascade reservoir combined flood control optimal scheduling method considering branch flood storage civil application
CN109002946B (en) Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes
CN114357868A (en) Multi-target cooperative scheduling method and device for complex flood control system
CN109033517A (en) A kind of Engineering Optimizational Method improving Storm Sewer Network design scheme applicability
Zhu et al. Optimal operation of a parallel multireservoir system for flood control using a stagewise compensation method
CN108330896A (en) A kind of reservoir dispatching method and system
CN115983733B (en) Electric power market clear data processing method and device based on water level control
Fang et al. Parallel improved DPSA algorithm for medium-term optimal scheduling of large-scale cascade hydropower plants
JP2001160043A (en) Simulation apparatus for water pipe network

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