CN112561213B - 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
CN112561213B
CN112561213B CN202110199694.6A CN202110199694A CN112561213B CN 112561213 B CN112561213 B CN 112561213B CN 202110199694 A CN202110199694 A CN 202110199694A CN 112561213 B CN112561213 B CN 112561213B
Authority
CN
China
Prior art keywords
reservoir
period
nth
discharge
interval
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
CN202110199694.6A
Other languages
Chinese (zh)
Other versions
CN112561213A (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

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • 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 initial dispatching lines by initial solutions 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 the first target value of the dispatching period by solving the total target value of the dispatching period based on the initial dispatching lines and the water balance equationnA 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 a time interval discharge interval of the nth reservoir in the i time interval according to a water quantity balance equation and a constraint condition;
s3, calculating the average warehousing flow of the time periods based on the warehousing flow of the time periods in the nth reservoir scheduling period of the i time period;
s4, judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge interval, and determining an initial solution of the discharge of the nth reservoir in each time interval according to the judgment result;
s5, obtaining an initial dispatching line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time interval of the nth reservoir and the time interval warehousing flow, and determining the dispatching line of the nth reservoir by solving a total target value of a dispatching period based on the initial dispatching line and the water quantity balance equation.
Optionally, the water balance equation is:
Figure GDA0003012484460000021
wherein,
Figure GDA0003012484460000022
the discharge quantity of the nth reservoir in the period i,
Figure GDA0003012484460000023
for forecasting the storage flow of the nth reservoir in the period ii nThe water storage capacity of the nth reservoir in the period i,
Figure GDA0003012484460000024
the water storage capacity of the nth reservoir at the end of the previous period of the period i, zh is a conversion coefficient, and tiIs the time length of the i period 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 GDA0003012484460000031
the reservoir water storage capacity constraint is as follows:
Figure GDA0003012484460000032
wherein,
Figure GDA0003012484460000033
maximum flood discharge flow of nth reservoir in period i, Vi nThe water storage capacity of the nth reservoir in the period of i, ZnDead level of the nth reservoir at time i, FZnThe flood level of the nth reservoir is checked in the period i.
Optionally, step S2 includes:
s21, converting the characteristic water level of the nth reservoir into corresponding water storage capacity by using the water level reservoir capacity curve; 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, carrying out the following discrete processing on the storage capacity of the nth reservoir in the period i meeting the reservoir storage capacity constraint:
Figure GDA0003012484460000034
wherein, S is discrete number of parts, S1, 2nThe amount of water stored, VF, corresponding to the dead water level of the nth reservoirnChecking the water storage amount corresponding to the flood level for the nth reservoir,
Figure GDA0003012484460000035
discrete water storage capacity of the nth reservoir in the period i;
s23, discrete water storage amount of the nth reservoir in a plurality of i periods
Figure GDA0003012484460000036
Substituting the water quantity into the water quantity balance equation to obtain the discharge quantity of the nth reservoir in N i time periods
Figure GDA0003012484460000037
From N discharge flows
Figure GDA0003012484460000038
Selecting a maximum value that meets said flood discharge capacity constraint
Figure GDA0003012484460000039
And minimum value
Figure GDA00030124844600000310
Further determining the discharge interval of the nth reservoir in the period i
Figure GDA00030124844600000311
Optionally, in step S3, the time-interval average warehousing flow rate qe is obtained by the following formula:
Figure GDA00030124844600000312
t is the total number of periods within the scheduling period.
Optionally, step S4 includes:
s41, judging the time interval average warehousing flow qe of the nth reservoirnWhether or not to discharge the flow rate interval in time period
Figure GDA0003012484460000041
Within;
s42a, if the time interval average warehousing flow qe of the nth reservoirnInterval of time interval discharge
Figure GDA0003012484460000042
Within the time interval, the initial solution of the discharge of the nth reservoir in each time interval
Figure GDA0003012484460000043
Time-interval-based average warehousing flow qenA value of (d);
s42b, if the time interval average warehousing flow qe of the nth reservoirnInterval of discharge quantity not in time period
Figure GDA0003012484460000044
Within the time interval, the initial solution of the discharge of the nth reservoir in each time interval
Figure GDA0003012484460000045
Interval of time interval and discharge
Figure GDA0003012484460000046
Average warehousing flow qe of the period closest to the nth reservoirnThe value of (c).
Optionally, step S5 includes:
s51, converting the starting water level of the nth reservoir into the water storage capacity of the 1 st time period of the reservoir by using a water level reservoir capacity curve, and then determining an initial scheduling line according to the water storage capacity of the 1 st time period of the nth reservoir, the initial solution of the reservoir discharge capacity of the nth reservoir in each time period and the time 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 GDA0003012484460000047
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 GDA0003012484460000048
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 quantity balance equation, and by combining the interval inflow between the nth reservoir in the period i and the flood control point, the total target value of the dispatching period is obtained:
Figure GDA0003012484460000049
wherein,
Figure GDA00030124844600000410
the method comprises the steps that inflow flows from the nth reservoir to a flood control point of the nth reservoir in an i period, N is the number of reservoirs, i is 1, 2, T, N is 1, 2, and N;
and 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 sequentially connecting the water storage amounts corresponding to the scheduling period total target value of each period to obtain the scheduling line of the nth reservoir.
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;
the time interval discharge interval calculation module is used for solving a time interval discharge interval of the nth reservoir in the i time interval according to a water quantity balance equation and a constraint condition;
the average warehousing flow calculation module is used for calculating the average warehousing flow of the time periods based on the warehousing flow of the time periods in the nth reservoir scheduling period of the time period i;
the judgment module is used for judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge flow interval or not;
the initial solution determining module is used for determining an initial solution of the discharge capacity of the nth reservoir in each time interval according to the judgment result;
and the scheduling line calculation module is used for obtaining an initial scheduling line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time interval of the nth reservoir and the time interval warehousing flow, and determining the scheduling line of the nth reservoir by solving a total target value of the scheduling period based on the initial scheduling line and the water quantity balance equation.
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, solving a time interval discharge interval of the nth reservoir at the time interval i according to a water quantity balance equation and a constraint condition; then, calculating the average warehousing flow of the time period in the nth reservoir scheduling period; judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge interval or not, and determining an initial solution of the discharge of the nth reservoir at each time interval according to the judgment result; and finally, obtaining an initial dispatching line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time period of the nth reservoir and the time period warehousing flow, and determining the dispatching line of the nth reservoir by solving a total target value of the dispatching period based on the initial dispatching line and a water quantity balance equation.
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. The limited flood season water level of the known nth reservoir is ZXnStarting water level of nth reservoir
Figure GDA0003012484460000071
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 GDA0003012484460000081
Final stage with final water level as control water level
Figure GDA0003012484460000082
Figure GDA0003012484460000083
N is the number of reservoirs, and T is the total number of periods in the scheduling period.
And S2, solving the time interval discharge rate interval of the nth reservoir in the i time interval according to the water quantity balance equation and the constraint condition.
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, converting the characteristic water level of the nth reservoir into corresponding water storage capacity by using the water level reservoir capacity curve; 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, carrying out the following discrete processing on the water storage quantity of the nth reservoir in the period i meeting the reservoir water storage quantity constraint:
Figure GDA0003012484460000084
wherein, S is discrete number of parts, S is 1, 2i nAll get through under all conditions
Figure GDA0003012484460000085
The value N is the number of reservoirs, i 1, 2, T, VSnThe amount of water stored, VF, corresponding to the dead water level of the nth reservoirnChecking the water storage amount corresponding to the flood level for the nth reservoir,
Figure GDA0003012484460000086
discrete water storage amount of the nth reservoir in the period i.
S23, discrete water storage amount of the nth reservoir in a plurality of i periods
Figure GDA0003012484460000087
Substituting into a water quantity balance equation to obtain the discharge quantity of the nth reservoir in N i time periods
Figure GDA0003012484460000088
From N discharge flows
Figure GDA0003012484460000089
Selecting the maximum value meeting the flood discharge capacity constraint
Figure GDA00030124844600000810
And minimum value
Figure GDA00030124844600000811
Further determining the discharge interval of the nth reservoir in the period i
Figure GDA00030124844600000812
Further, the water balance equation is:
Figure GDA00030124844600000813
wherein,
Figure GDA00030124844600000814
the discharge capacity of the nth reservoir in the period i is determined according to actual needs, the length of each period can be 1h, 3h, 1 day, 1 ten days or 1 month and the like,
Figure GDA00030124844600000815
for forecasting the storage flow of the nth reservoir in the period ii nThe water storage capacity of the nth reservoir in the period i,
Figure GDA00030124844600000816
the storage capacity of the nth reservoir at the end of the period before the period i, and zh is a conversion coefficient (for example, the storage capacity of a large reservoir is generally billion m3Flow rate in m3/s, where the conversion factor may be determined to be 100000000), tiIs the time length of the i period 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 GDA0003012484460000091
the reservoir water storage capacity constraint is as follows:
Zn≤Vi n≤FZn，
the water level constraint in the scheduling period is as follows:
Zendn＝Zen，
the variation of amplitude of flow out of the warehouse is restricted as follows:
Figure GDA0003012484460000092
wherein,
Figure GDA0003012484460000093
maximum flood discharge flow of nth reservoir in period i, Vi nThe water storage capacity of the nth reservoir in the period of i, ZnDead level of the nth reservoir at time i, FZnFor checking flood level, Zend, of the nth reservoir at period inScheduling end of term water level, Ze, for the nth reservoirnThe water level is finally controlled for the nth reservoir,
Figure GDA0003012484460000094
the discharge quantity qm of the nth reservoir at the end of the period before the period inThe maximum amplitude of the permitted delivery flow of the nth reservoir in the adjacent time interval. The scheduling end-of-term water level constraint and the ex-warehouse flow amplitude constraint are used in the subsequent application of a POA algorithm, and the water storage capacity of the nth reservoir in each discrete i time period is used
Figure GDA0003012484460000095
Calculating the target value, and if any constraint condition is not satisfied during the calculation process, discarding the state value
Figure GDA0003012484460000096
And S3, calculating the average warehousing flow of the time periods based on the warehousing flow of the time periods in the nth reservoir scheduling period of the i time period.
Further, in step S3, the time-interval average warehousing traffic qenThe following formula is used to obtain:
Figure GDA0003012484460000097
t is the total number of periods within the scheduling period.
S4, judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge interval, and determining an initial solution of the discharge of the nth reservoir in each time interval according to the judgment result.
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, judging the time interval average warehousing flow qe of the nth reservoirnWhether or not to discharge the flow rate interval in time period
Figure GDA0003012484460000101
Within;
s42a, if the time interval average warehousing flow qe of the nth reservoirnInterval of time interval discharge
Figure GDA0003012484460000102
Within the time interval, the initial solution of the discharge of the nth reservoir in each time interval
Figure GDA0003012484460000103
Time-interval-based average warehousing flow qenA value of (d);
s42b, if the time interval average warehousing flow qe of the nth reservoirnInterval of discharge quantity not in time period
Figure GDA0003012484460000104
In the initial stage of the discharge of the nth reservoir in each periodBeginning to solve
Figure GDA0003012484460000105
Interval of time interval and discharge
Figure GDA0003012484460000106
Average warehousing flow qe of the period closest to the nth reservoirnThe value of (c).
S5, obtaining an initial dispatching line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time interval of the nth reservoir and the time interval warehousing flow, and determining the dispatching line of the nth reservoir by solving a total target value of a dispatching period based on the initial dispatching line and a water balance equation.
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, converting the starting water level of the nth reservoir into the water storage capacity of the 1 st time period of the reservoir by using the water level reservoir capacity curve, and then determining an initial dispatching line according to the water storage capacity of the 1 st time period of the nth reservoir, the initial solution of the reservoir discharge capacity of the nth reservoir in each time period and the time 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 GDA0003012484460000107
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 GDA0003012484460000108
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 GDA0003012484460000109
With constant value for other periodsUnder the condition, the water storage amount of the ith time interval of the nth reservoir in the initial scheduling line is changed every time, and the water storage amount in each discrete state is calculated
Figure GDA00030124844600001010
And comparing the target values in each discrete state to select the optimal water storage amount.
In the above step, hold
Figure GDA0003012484460000111
And
Figure GDA0003012484460000112
calculating the value of the 1 st moment of the nth reservoir after dispersion without change
Figure GDA0003012484460000113
Comparing the target values in each discrete state to select the optimal storage capacity of the reservoir
Figure GDA0003012484460000114
Make it
Figure GDA0003012484460000115
Then hold
Figure GDA0003012484460000116
And
Figure GDA0003012484460000117
unchanged, using the value of the discrete nth reservoir at time 2
Figure GDA0003012484460000118
Computing at each discrete state
Figure GDA0003012484460000119
Comparing the target values of each state to select the optimal storage capacity of the reservoir
Figure GDA00030124844600001110
Make it
Figure GDA00030124844600001111
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 quantity balance equation, and by combining the interval inflow between the nth reservoir in the period i and the flood control point, the total target value of the dispatching period is obtained:
Figure GDA00030124844600001112
wherein,
Figure GDA00030124844600001113
the discharge quantity of the nth reservoir in the period i,
Figure GDA00030124844600001114
and the inflow of the interval from the nth reservoir to the flood control point of the nth reservoir in the period i, wherein N is the number of the reservoirs, i is 1, 2, and T, N is 1, 2.
And 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 sequentially connecting the water storage amounts corresponding to the scheduling period total target value of each period to obtain the scheduling line of the nth 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 GDA00030124844600001115
zh、tiThe water storage capacity of the reservoir corresponding to the starting water level of the nth reservoir can be combined to obtain the second reservoirnVariation process of water storage capacity of individual reservoir
Figure GDA00030124844600001116
And drawing a curve by using the water storage capacity point corresponding to the total target value of the scheduling period of each time interval of the nth reservoir, wherein the curve is the scheduling line of the nth 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 GDA0003012484460000121
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
Figure GDA0003012484460000122
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.
And the time interval discharge rate interval calculation module 302 is used for solving the time interval discharge rate interval of the nth reservoir in the i time interval according to the water quantity balance equation and the constraint condition.
And the average warehousing flow calculation module 303 is configured to calculate the average warehousing flow in a time period based on the warehousing flow in the time period in the nth reservoir scheduling period in the time period i.
The judging module 304 is configured to judge whether the time-interval average warehousing flow rate of the nth reservoir is within the corresponding time-interval discharge rate interval.
And an initial solution determining module 305, configured to determine an initial solution of the discharge amount of the nth reservoir in each time interval according to the determination result.
And the dispatching line calculating module 306 is used for obtaining an initial dispatching line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time period of the nth reservoir and the time period warehousing flow, and determining the dispatching line of the nth reservoir by solving a total target value of a dispatching period based on the initial dispatching line and a water quantity balance equation.
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 407 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.
And S2, solving the time interval discharge rate interval of the nth reservoir in the i time interval according to the water quantity balance equation and the constraint condition.
And S3, calculating the average warehousing flow of the time periods based on the warehousing flow of the time periods in the nth reservoir scheduling period of the i time period.
S4, judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge interval, and determining an initial solution of the discharge of the nth reservoir in each time interval according to the judgment result.
S5, obtaining an initial dispatching line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time interval of the nth reservoir and the time interval warehousing flow, and determining the dispatching line of the nth reservoir by solving a total target value of a dispatching period based on the initial dispatching line and a water balance equation.
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 (7)

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 a time interval discharge interval of the nth reservoir in the i time interval according to a water quantity balance equation and a constraint condition;
the water quantity balance equation is as follows:
Figure FDA0003012484450000011
wherein,
Figure FDA0003012484450000012
the discharge quantity of the nth reservoir in the period i,
Figure FDA0003012484450000013
for forecasting the storage flow of the nth reservoir in the period ii nThe water storage capacity of the nth reservoir in the period i,
Figure FDA0003012484450000014
the water storage capacity of the nth reservoir at the end of the previous period of the period i, zh is a conversion coefficient, and tiIs the time length of the i time period and has the unit of second;
the constraint conditions comprise flood discharge capacity constraint and reservoir water storage capacity constraint;
the flood discharge capacity constraints are:
Figure FDA0003012484450000015
the reservoir water storage capacity constraint is as follows:
Zn≤Vi n≤FZn，
wherein,
Figure FDA0003012484450000016
maximum flood discharge flow of nth reservoir in period i, Vi nThe water storage capacity of the nth reservoir in the period of i, ZnDead level of the nth reservoir at time i, FZnChecking the flood level of the nth reservoir in the period i;
step S2 includes:
s21, converting the characteristic water level of the nth reservoir into corresponding water storage capacity by using the water level reservoir capacity curve; 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, carrying out the following discrete processing on the storage capacity of the nth reservoir in the period i meeting the reservoir storage capacity constraint:
Figure FDA0003012484450000017
wherein, S is discrete number of parts, S1, 2nThe amount of water stored, VF, corresponding to the dead water level of the nth reservoirnChecking the water storage amount corresponding to the flood level for the nth reservoir,
Figure FDA0003012484450000021
discrete water storage capacity of the nth reservoir in the period i;
s23, discrete water storage amount of the nth reservoir in a plurality of i periods
Figure FDA0003012484450000022
Substitution intoObtaining the discharge quantity of the nth reservoir in N i time periods in the water quantity balance equation
Figure FDA0003012484450000023
From N discharge flows
Figure FDA0003012484450000024
Selecting a maximum value that meets said flood discharge capacity constraint
Figure FDA0003012484450000025
And minimum value
Figure FDA0003012484450000026
Further determining the discharge interval of the nth reservoir in the period i
Figure FDA0003012484450000027
S3, calculating the average warehousing flow of the time periods based on the warehousing flow of the time periods in the nth reservoir scheduling period of the i time period;
s4, judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge interval, and determining an initial solution of the discharge of the nth reservoir in each time interval according to the judgment result;
s5, obtaining an initial dispatching line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time interval of the nth reservoir and the time interval warehousing flow, and determining the dispatching line of the nth reservoir by solving a total target value of a dispatching period based on the initial dispatching line and the water quantity balance equation.
2. The improved POA-based reservoir flood control scheduling method according to claim 1, wherein in step S3, the time-interval average warehousing flow rate is obtained by the following formula:
Figure FDA0003012484450000028
t is the total number of periods within the scheduling period.
3. The improved POA-based reservoir flood control scheduling method of claim 2, wherein the step S4 comprises:
s41, judging the time interval average warehousing flow qe of the nth reservoirnWhether or not to discharge the flow rate interval in time period
Figure FDA0003012484450000029
Within;
s42a, if the time interval average warehousing flow qe of the nth reservoirnInterval of time interval discharge
Figure FDA00030124844500000210
Within the time interval, the initial solution of the discharge of the nth reservoir in each time interval
Figure FDA00030124844500000211
Time-interval-based average warehousing flow qenA value of (d);
s42b, if the time interval average warehousing flow qe of the nth reservoirnInterval of discharge quantity not in time period
Figure FDA00030124844500000212
Within the time interval, the initial solution of the discharge of the nth reservoir in each time interval
Figure FDA00030124844500000213
Interval of time interval and discharge
Figure FDA00030124844500000214
Average warehousing flow qe of the period closest to the nth reservoirnThe value of (c).
4. The improved POA-based reservoir flood control scheduling method of claim 3, wherein the step S5 comprises:
s51, converting the starting water level of the nth reservoir into the water storage capacity of the 1 st time period of the reservoir by using a water level reservoir capacity curve, and then determining an initial scheduling line according to the water storage capacity of the 1 st time period of the nth reservoir, the initial solution of the reservoir discharge capacity of the nth reservoir in each time period and the time 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 FDA0003012484450000031
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 FDA0003012484450000032
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 quantity balance equation, and by combining the interval inflow between the nth reservoir in the period i and the flood control point, the total target value of the dispatching period is obtained:
Figure FDA0003012484450000033
wherein,
Figure FDA0003012484450000034
the method comprises the steps that inflow flows from the nth reservoir to a flood control point of the nth reservoir in an i period, N is the number of reservoirs, i is 1, 2, T, N is 1, 2, and N;
and 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 sequentially connecting the water storage amounts corresponding to the scheduling period total target value of each period to obtain the scheduling line of the nth reservoir.
5. 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;
the time interval discharge interval calculation module is used for solving a time interval discharge interval of the nth reservoir in the i time interval according to a water quantity balance equation and a constraint condition;
the water quantity balance equation is as follows:
Figure FDA0003012484450000035
wherein,
Figure FDA0003012484450000036
the discharge quantity of the nth reservoir in the period i,
Figure FDA0003012484450000037
for forecasting the storage flow of the nth reservoir in the period ii nThe water storage capacity of the nth reservoir in the period i,
Figure FDA0003012484450000038
the water storage capacity of the nth reservoir at the end of the previous period of the period i, zh is a conversion coefficient, and tiIs the time length of the i time period and has the unit of second;
the constraint conditions comprise flood discharge capacity constraint and reservoir water storage capacity constraint;
the flood discharge capacity constraints are:
Figure FDA0003012484450000041
the reservoir water storage capacity constraint is as follows:
Zn≤Vi n≤FZn，
wherein,
Figure FDA0003012484450000042
maximum flood discharge flow of nth reservoir in period i, Vi nThe water storage capacity of the nth reservoir in the period of i, ZnDead level of the nth reservoir at time i, FZnChecking the flood level of the nth reservoir in the period i;
the time interval leakage flow rate interval calculation module is specifically configured to:
converting the characteristic water level of the nth reservoir into corresponding water storage capacity by using a water level reservoir capacity curve; 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;
carrying out the following discrete processing on the water storage capacity of the nth reservoir in the period i meeting the reservoir water storage capacity constraint:
Figure FDA0003012484450000043
wherein, S is discrete number of parts, S1, 2nThe amount of water stored, VF, corresponding to the dead water level of the nth reservoirnChecking the water storage amount corresponding to the flood level for the nth reservoir,
Figure FDA0003012484450000044
discrete water storage capacity of the nth reservoir in the period i;
discrete water storage amount of nth reservoir in a plurality of i periods
Figure FDA0003012484450000045
Substituting the water quantity into the water quantity balance equation to obtain the discharge quantity of the nth reservoir in N i time periods
Figure FDA0003012484450000046
From N discharge flows
Figure FDA0003012484450000047
Selecting a maximum value that meets said flood discharge capacity constraint
Figure FDA0003012484450000048
And minimum value
Figure FDA0003012484450000049
Further determining the discharge interval of the nth reservoir in the period i
Figure FDA00030124844500000410
The average warehousing flow calculation module is used for calculating the average warehousing flow of the time periods based on the warehousing flow of the time periods in the nth reservoir scheduling period of the time period i;
the judgment module is used for judging whether the time-interval average warehousing flow of the nth reservoir is within the corresponding time-interval discharge flow interval or not;
the initial solution determining module is used for determining an initial solution of the discharge capacity of the nth reservoir in each time interval according to the judgment result;
and the scheduling line calculation module is used for obtaining an initial scheduling line of the nth reservoir according to the starting water level of the nth reservoir, the initial solution of the reservoir discharge of each time interval of the nth reservoir and the time interval warehousing flow, and determining the scheduling line of the nth reservoir by solving a total target value of the scheduling period based on the initial scheduling line and the water quantity balance equation.
6. An electronic device, comprising:
a processor;
a memory for storing instructions for the processor to control the steps of the method for improved POA-based reservoir flood control scheduling according to any one of claims 1-4.
7. 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 according to any of claims 1-4.
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 CN112561213A (en) 2021-03-26
CN112561213B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN114118554B (en) * 2021-11-18 2024-05-14 武汉大学 Dynamic control method and system for reservoir flood period running water level based on typhoon path
CN113934777B (en) * 2021-12-16 2022-03-04 长江水利委员会水文局 Method and system for quantifying influence of backwater jacking on water level change

Family Cites Families (2)

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

Also Published As

Publication number Publication date
CN112561213A (en) 2021-03-26

Similar Documents

Publication Publication Date Title
CN112561213B (en) Reservoir flood control scheduling method, system, equipment and medium based on improved POA
CN112633578B (en) Cascade reservoir group optimal scheduling method under influence of diversion and water diversion engineering
CN107045785B (en) A method of the short-term traffic flow forecast based on grey ELM neural network
CN108109076B (en) Method for analyzing risk of abandoned water in power generation dispatching of cascade hydropower station group by considering runoff forecasting
Azari et al. Multi-objective optimization of the reservoir system operation by using the hedging policy
CN110334879A (en) Power grid bus reactive load forecasting method and device
CN104636834B (en) A kind of improved joint probability plan model system optimization method
Chen et al. Applying genetic algorithm and neural network to the conjunctive use of surface and subsurface water
Chen et al. Flood control operation of reservoir group using Yin-Yang Firefly Algorithm
CN107300853A (en) A kind of subregion calculates the fluctuating wind speed time series input method of spatial coherence
CN114492233A (en) Basin water simulation method based on webGIS platform and considering comprehensive utilization requirements
CN117873686B (en) Hierarchical load balancing-based distributed hydrological model conflux parallel scheduling algorithm
Fang et al. Parallel improved DPSA algorithm for medium-term optimal scheduling of large-scale cascade hydropower plants
Zhu et al. Optimal operation of a parallel multireservoir system for flood control using a stagewise compensation method
Papageorgiou Certainty equivalent open-loop feedback control applied to multireservoir networks
CN109325610A (en) A kind of piecewise linearity dispatching method towards multi-water resources cascade pumping station water transfer system
CN116307709B (en) Comprehensive assessment method and system for flood control capacity of transformer substation based on information gain fusion
CN115983733B (en) Electric power market clear data processing method and device based on water level control
CN109886543B (en) Reservoir optimal scheduling method based on uniform design
Zhang et al. Daily reservoir inflow forecasting combining QPF into ANNs model
CN110989040B (en) Artificial intelligent lightning approach early warning method and system based on slice processing
CN109190160B (en) Matrixing simulation method of distributed hydrological model
CN112613654A (en) Comprehensive energy system flexibility evaluation method based on multi-type energy storage
Grigoriu et al. Probabilistic models for directionless wind speeds in hurricanes
CN113379134B (en) Intelligent scheduling method for adjusting reservoir for multiple years

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