CN108537449B - Reservoir coordination scheduling strategy obtaining method considering river flood demand - Google Patents

Reservoir coordination scheduling strategy obtaining method considering river flood demand Download PDF

Info

Publication number
CN108537449B
CN108537449B CN201810326686.1A CN201810326686A CN108537449B CN 108537449 B CN108537449 B CN 108537449B CN 201810326686 A CN201810326686 A CN 201810326686A CN 108537449 B CN108537449 B CN 108537449B
Authority
CN
China
Prior art keywords
reservoir
flow
demand
flood
water level
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
CN201810326686.1A
Other languages
Chinese (zh)
Other versions
CN108537449A (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.)
Changjiang Institute of Survey Planning Design and Research Co Ltd
Original Assignee
Changjiang Institute of Survey Planning Design and Research 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 Changjiang Institute of Survey Planning Design and Research Co Ltd filed Critical Changjiang Institute of Survey Planning Design and Research Co Ltd
Priority to CN201810326686.1A priority Critical patent/CN108537449B/en
Publication of CN108537449A publication Critical patent/CN108537449A/en
Application granted granted Critical
Publication of CN108537449B publication Critical patent/CN108537449B/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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • G06Q10/06375Prediction of business process outcome or impact based on a proposed change
    • 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

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Control Of Non-Electrical Variables (AREA)

Abstract

The invention provides a reservoir coordination scheduling strategy obtaining method considering river flood demand, which comprises the following steps: step one, acquiring a forecast value of reservoir warehousing runoff of a current day, an expected control value of reservoir water level at the end of the current day and a power generation command regulating output process of a hydropower station of the current day; secondly, considering the hierarchical flood demand of the downstream river channel, and constructing a reservoir coordination scheduling model considering the flood demand of the river channel; thirdly, integrating constraint conditions in the coordination scheduling model on the basis of the objective function, and solving the established reservoir coordination scheduling model by using a dynamic programming method; and fourthly, providing key indexes from the aspects of reservoir dispatching safety and downstream river flood safety, constructing a safety evaluation index system suitable for coordinating the characteristics of the dispatching model, and comprehensively evaluating and feeding back the safety of the obtained dispatching strategy. Compared with the prior art, the downstream river flood demand is not considered in detail, and the method can provide a technical support scheme for promoting river flood safety management.

Description

Reservoir coordination scheduling strategy obtaining method considering river flood demand
Technical Field
The invention relates to the technical field of hydraulic engineering management and hydropower energy scheduling, in particular to a reservoir coordination scheduling strategy acquisition method considering river flood demand.
Background
The economic situation of China gradually enters the normality, the annual electricity consumption increase trend of the society is slow, but the peak-valley difference is increased year by year. With the continuous expansion of the access scale of new energy such as wind power, photoelectricity and the like, the demand degree and the dependence degree of a power grid on power resources with the characteristic of peak shaving capacity are higher. The hydropower is used as a green energy source, has the characteristics of quick response, low operation cost and the like, is an ideal peak regulation power supply, and has great significance for fully exerting the excellent regulation performance of the hydropower to relieve peak regulation of a power grid and environmental stress.
When a hydropower station executes tasks such as power grid peak regulation, frequency regulation and the like, the generated power flow rapidly rises and falls, the unstable flow of the reservoir discharged from the reservoir easily causes large fluctuation of the section flow of a downstream river channel, and if residential life or production operation points are distributed along the river channel, the flood situation of the river channel is further aggravated. In order to relieve or even eliminate adverse effects of peak-shaving unsteady flow of the hydropower station on flood of a downstream river channel, a reservoir with a certain adjusting performance is usually built downstream of the hydropower station. Because the peak-shaving resources of the power grid are increasingly scarce, the reservoir with poor regulation performance sometimes also plays a role in peak-shaving and frequency-modulation of the power grid. Generally, the adjustment reservoir capacity of a built reservoir is limited, so that the lifting speed of the reservoir water level is high, and in order to avoid the dam overflow or emptying risk of the reservoir, a reservoir gate must be frequently used for adjustment. Therefore, on the premise of ensuring the operation and dispatching safety of the hydropower station reservoir engineering, how to coordinate the peak regulation capacity and the flow stabilization function is one of the problems to be solved urgently in theoretical research and actual dispatching work at present.
Patent document CN102817335A discloses a method and system for optimal scheduling of combined flood control of cascade reservoir groups. The invention can dynamically judge the level of flood and realize the hierarchical scheduling on the basis, and simultaneously, selects a proper optimized scheduling target to realize the requirements of protecting downstream protection objects and dam safety and flood reclamation according to the comprehensive information of the regulation performance of the reservoir, the scheduling period, the reliability of forecast incoming water and the like in the hierarchical scheduling.
Patent document CN105869065A discloses a scheduling method for coordinating flood control risk and power generation benefit in flood season of hydropower stations. The hydropower station flood season double-target two-stage optimization scheduling model is established by considering hydrologic forecast information and uncertainty influence thereof, the hydrologic forecast information is fully utilized, flood season scheduling is carried out according to a model solving result, and maximization of comprehensive benefits of power generation and flood control is achieved.
Although the two inventions have differences in main target problems, key technical methods and the like, when the flood demand of the downstream river channel is considered, the maximum flow or the maximum water level of the control station is used as a key index for measuring whether the flood situation is safe or not, but the flood index of the downstream river channel is not taken into the invention authority.
Patent document CN106327022A discloses a method and a device for stabilizing fluctuation of power generation flow of a cascade hydropower station. On the basis of obtaining the predicted flow rate in a warehouse, a daily load prediction curve, an optimization target and a peak-to-charge ratio, under the condition of meeting the constraints of hydropower station operation, water level amplitude variation and flow fluctuation, the maximization of peak regulation benefit is realized, and meanwhile, the flow pulsation of peak regulation output is stabilized by applying a sliding average filtering method, so that the stability of flow rate out of the warehouse is improved. The method takes flood indexes of a downstream river channel into consideration in a refining manner, but does not take supporting measures such as river channel flood discharge early warning implementation, gate operation pressure relief and the like into the invention permission, and if the measures are taken into consideration, the application is richer.
Disclosure of Invention
The invention aims to provide a reservoir coordination scheduling strategy obtaining method considering river channel flood demand, aiming at the defects of the prior art, and providing a technical support scheme for improving outlet flow water flow form and promoting river channel flood safety management.
The invention provides a reservoir coordination scheduling strategy acquisition method considering river flood demand, which is characterized by comprising the following steps of:
step one, acquiring a forecast value of reservoir warehousing runoff of a current day, an expected control value of reservoir water level at the end of the current day and a power generation command regulating output process of a hydropower station of the current day;
secondly, considering the hierarchical flood demand of the downstream river channel, and constructing a reservoir coordination scheduling model considering the flood demand of the river channel;
thirdly, integrating and coordinating all constraint conditions in the reservoir coordination scheduling model on the basis of the objective function;
fourthly, solving the established reservoir coordination scheduling model by using a dynamic programming method;
fifthly, providing key indexes from two aspects of reservoir dispatching safety and downstream river flood safety, constructing a safety evaluation index system suitable for coordinating the characteristics of a dispatching model, and comprehensively evaluating and feeding back the safety of the dispatching strategy obtained by calculation in the fourth step; and if the safety evaluation index requirement is not met, returning to the second step, adjusting the parameters of the reservoir coordination scheduling model and recalculating until the safety of the scheduling strategy obtained by calculation meets the safety evaluation index requirement.
In the above technical solution, the objective function is:
Figure RE-GDA0001735091050000041
in the formula: fabs (. cndot.) as a function of absolute value; qout, t is the delivery flow of the reservoir in time t; qout, t-1 is the delivery flow of the reservoir in the time period t-1; t is a scheduling period; and lambdat is a river flood early warning demand coefficient of the reservoir at the time t.
The hierarchical acquisition mode of the demand coefficient lambdat is as follows:
Figure RE-GDA0001735091050000042
in the formula: SET { indices } is the selected index SET; fi (-) is an early warning distinguishing mode corresponding to the ith level basic flow; λ i (i ═ 1,2, …, n) is a set constant; the concrete expression mode of the SET { indices } and fi (-) is determined by combining the actual characteristics of engineering.
In the above technical solution, the constraint condition includes:
a. water balance constraint
Vt+1=Vt+(It-Qpower,t-Qab,t)Δt (3)
In the formula: vt, Vt +1 are the initial and final storage capacity of the reservoir in the period t respectively; it, Qpower, t, Qab, t are respectively the warehousing flow, the power generation flow and the overflow and abandon flow of the reservoir in the period t; Δ t is the period length;
b. water level restraint
Zt,min≤Zt≤Zt,max (4)
In the formula: zt, min, Zt, max are respectively the water level, the lowest water level and the highest water level allowed by the reservoir in the time period t.
c. Let-down flow restriction
qt,min≤Qpower,t+Qab,t≤qt,max (5)
In the formula: qt, min, qt and max are respectively a lower discharge flow limit and an upper discharge flow limit of the reservoir in the time period t;
d. initial and final water level control
Z0=ZFirst stage,ZT=ZPowder (6)
In the formula: z0 and ZT are the initial and final water levels of the dispatching period of the reservoir respectively; zFirst stage,Z PowderThe water level values are respectively given at the beginning and the end of the day.
In the above technical solution, the fourth step includes the following steps:
s1: dividing stages and determining stage variables, dividing a reservoir dispatching date (one day) into 96 stages according to a 15-minute time interval, and regarding the stages as a multi-stage decision problem, wherein T represents the stage variables (T is 1-T); t is the facing time period, and T + 1-T is the remaining time period;
s2: determining state variable, selecting reservoir capacity as state variable, and recording Vt-1The initial reservoir capacity at the t-th time period, VtThe storage capacity of the reservoir at the end of the t period;
s3: determining decision variables and selecting reservoir discharge flow Qout,t=Qpower,t+Qab,tIs a decision variable;
s4: determining a state transition equation, namely a water quantity balance equation, an equation (3);
s5: determining stage indexes, namely river flood warning demand coefficient lambda t and reservoir delivery flow amplitude delta Q of each stageoutThe product of (a);
s6: determining an optimum function, i.e. ft *(Vt-1) Denotes an initial storage capacity V from the t-th staget-1Starting from the sum of the optimal water abandoning amplitude in the T-th time period, the reverse time sequence recurrence equation of the dynamic programming can be obtained as:
Figure RE-GDA0001735091050000061
In the formula: vt-1 m1Represents a bank capacity value V when the initial discrete point of the t-th period is m1t m2Denotes a storage capacity value at which the discrete point at the end of the t-th period (beginning of t +1 period) is taken as M2, and M1 is 0,1, …, M; m2 ═ 0,1, …, M; m is the discrete point number of the storage capacity; f. oft *(Vt-1 m1) Denotes the initial storage capacity V from the t-th periodt-1 m1Starting to sum of the optimal water abandoning amplitude in the T-th time period; f. of* t+1(Vt m2) Represents that the initial storage capacity is V from the t +1 th periodt m2Starting to sum of the optimal water abandoning amplitude in the T-th time period; omega t is a decision variable and represents the warehouse-out flow QtStorage capacity Vt-1A decision set meeting various constraints of the reservoir of the hydropower station at a given time;
s7: calculating reverse time sequence recursion, according to the formula (7), carrying out forward time sequence recursion from the last stage (stage T) to the first stage, and solving the reverse time sequence process with the minimum target function of the reservoir at each period in the whole dispatching period under the condition of meeting related constraint conditions;
s8: and (4) performing sequential recursive calculation to obtain the optimal strategy (Q) of the ex-warehouse flow process corresponding to the optimal valuetAnd corresponding reservoir capacity optimal state point value { V }t}。
In the above technical solution, the obtaining manner of the demand coefficient is as follows:
when the overflow flow variation is within the range of the flood demand of the river channel, early warning information does not need to be issued, the demand coefficient lambda is 1, and the dispatching objective is to pursue the stability of the flow process of leaving the reservoir;
when the scheduling strategy breaks through the early warning requirement of the overflow flow amplitude, starting an early warning program, and taking the overflow flow amplitude as the requirement coefficient;
the larger the demand coefficient is, the larger the overflow and abandon flow amplitude is, the more adverse to the downstream river flood situation is, so the demand for implementing early warning is higher, the pressure born by corresponding scheduling work is increased, and the pressure is avoided when an optimization decision is made.
In the above technical solution, the key indexes in the fifth step include: the number of times of early warning water level upper limit crossing of reservoir water level, the number of times of early warning water level lower limit crossing of reservoir water level, daily variation of reservoir water level, the number of times of early warning flood of river channel and the operation pressure of reservoir gate.
In the above technical solution, the safety evaluation index system in the fifth step is shown in the following table:
Figure RE-GDA0001735091050000071
Figure RE-GDA0001735091050000081
table 1 is a simple table of safety evaluation indexes of reservoir scheduling strategies
Wherein, the priority level decreases from top to bottom.
The invention provides a reservoir coordination scheduling strategy obtaining method considering river flood demand. The method aims at the practical engineering problem that the operation safety of a dam and the downstream river flood demand are difficult to coordinate in reservoir scheduling, provides a processing mode taking the river flood warning demand coefficient as a key technology according to the practical scheduling warning work requirement, constructs a reservoir coordination scheduling model taking the river flood demand into consideration, and provides a detailed process for solving the model by using a dynamic programming method. On the basis, according to the characteristics of the scheduling strategy obtained by the method, a safety evaluation index system is constructed from two aspects of reservoir scheduling safety and downstream river flood safety, and the purposes of comprehensive evaluation and decision feedback on the scheduling strategy safety can be realized. More importantly, the method takes reservoir dispatching regulations, hydropower station power generation plan approval processes and the like as calculation bases, so that the method is more mature and reliable in technology and wider in application prospect.
Drawings
FIG. 1 is a general design flow diagram of the present invention;
FIG. 2 is a schematic diagram of an optimization process for solving the established reservoir coordination scheduling model by using a dynamic programming method;
FIG. 3 is a process of output, warehousing flow and power generation flow of the first unfavorable condition of the reservoir A in the embodiment;
FIG. 4 illustrates a first adverse operating condition of the reservoir A;
FIG. 5 shows the first unfavorable condition of the reservoir A, namely the flow of leaving the reservoir and the flow of overflowing and discarding;
FIG. 6 is a schematic view of load regulation under the unfavorable condition II of the reservoir A in the embodiment;
FIG. 7 illustrates the flow of the reservoir A entering and exiting the reservoir A under the unfavorable operating condition II in the embodiment;
FIG. 8 is a process of generating flow and discarding flow of the second unfavorable condition of the reservoir A in the embodiment;
FIG. 9 illustrates a second embodiment of the reservoir level process for unfavorable conditions of the reservoir A;
FIG. 10 shows the output, warehousing flow and power generation flow of the third unfavorable condition of the reservoir A in the embodiment;
FIG. 11 shows the discharge and discharge flow of the third unfavorable condition of the reservoir A according to the embodiment;
fig. 12 shows the third reservoir level process of the unfavorable condition of the reservoir a in the embodiment.
Detailed Description
The invention will be further described in detail with reference to the following drawings and specific examples, which are not intended to limit the invention, but are for clear understanding.
As shown in fig. 1, the invention provides a reservoir coordination scheduling strategy obtaining method considering river flood demand.
Taking a certain reservoir A in Sichuan of China as an example, a coordinated dispatching strategy for considering river flood demand is obtained. And respectively selecting three unfavorable typical working conditions appearing in daily scheduling of the reservoir A for method application. The specific implementation mode is as follows:
step 1, acquiring data (power generation regulating command, forecast water and interval inflow shown in figure 1) of a current day warehousing runoff forecast value, a current day initial water level and last reservoir water level expected control value, a current day hydropower station power generation regulating command output process and the like. Wherein the water level control is shown in table 2.
Unfavorable operating conditions Initial water level control (m) End water level control (m) Remarks for note
Working condition one 1012.58 1012.23 Adverse natural factors
Working condition two 1012.87 1013.39 Adverse load factor
Three working conditions 1013.01 1013.42 Adverse combined factors
TABLE 2 control of end and beginning water level for three adverse conditions
And 2, considering the hierarchical flood demand of the downstream river channel, and constructing a reservoir A coordination scheduling model considering the flood demand of the river channel. The invention takes reservoir dispatching requirements into account, including reservoir water level out-of-limit, reservoir water level variation, power station peak regulation operation, and river flood demand, including reservoir outlet flow variation of the reservoir and river flood early warning mechanism.
The determination of the demand coefficient λ t in the objective function is based on the following: the local flood control and drought resistance department issues a guidance document related to the flood discharge early warning of the reservoir A, and indicates that basic flow and increase and decrease flow change data are adopted as early warning indexes in view of weak regulating capacity of the reservoir A and frequent opening and closing of a flood discharge gate. The method specifically comprises the following steps: when the basic flow is less than 1000m3(s) the flow of the opening and closing flood discharge gate is increased and decreased by more than 400m3S; base flow rate equal to or greater than 1000m3The/s is less than 2000m3S, the flow of the opening and closing flood discharge gate is increased and decreased by more than 600m3S; the basic flow is equal to or more than 2000 and less than 4000m3S, the flow of the open and close flood discharge gate is increased and decreased by more than 1000m3S; the basic flow rate is equal to or greater than 4000m3S, the flow of the opening and closing flood discharge gate is increased and decreased by more than 1500m3And/s, early warning is required to be implemented, and inspection along the river is organized.
Thus, the objective function of the coordinated scheduling model is:
Figure RE-GDA0001735091050000111
wherein the content of the first and second substances,
Figure RE-GDA0001735091050000112
the variables in the formula are characterized as defined above.
Step 3, integrating all constraint conditions in the coordinated scheduling model on the basis of the objective function, wherein the constraint conditions comprise the following steps: water balance constraint, water level constraint, downward discharge flow constraint and initial and final water level control.
And 4, solving by using the dynamic programming method to obtain the scheduling strategy of each working condition of the reservoir A. The optimizing process is shown in fig. 2 (in the figure, the abscissa is the number of time periods, and the ordinate is the reservoir capacity). The specific steps include 8 steps S1-S8, and formula (7) is a key formula in the stage optimization process.
The scheduling strategy of the unfavorable working condition one (unfavorable natural factors) is shown in fig. 3-5, and comprises an output process, an ex-warehouse flow process, a power generation flow process, an overflow and abandon flow process and a reservoir water level change process.
The scheduling strategy of the adverse working condition two (the load factor is adverse) is shown in fig. 6-9, and comprises an output process, an ex-warehouse flow process, a power generation flow process, an overflow and abandon flow process and a warehouse water level change process.
The scheduling strategy of the unfavorable working condition three (the comprehensive factor is unfavorable) is shown in fig. 10-12, and comprises an output process, an ex-warehouse flow process, a power generation flow process, an overflow and abandon flow process and a reservoir water level change process.
And 5, constructing a safety evaluation index system from the two aspects of reservoir scheduling safety and downstream river flood safety, and comprehensively evaluating and feeding back the safety of the obtained three adverse working condition scheduling strategies. The indexes are quantified according to the following steps:
index one: and early warning the upper limit times of the water level of the reservoir. The normal water level of the reservoir A is 1015.0m, a safety margin of 0.5m is reserved on the basis to set an early warning water level upper limit (namely 1014.5m), and when the reservoir water level exceeds the water level, early warning is needed to avoid the dam overflowing risk of the reservoir.
Index two: and early warning the lower limit times of the water level when the reservoir water level is over. The dead water level of the reservoir A is 1010.0m, a safety margin of 1.0m is reserved on the basis, an early warning water level lower limit (namely 1011.0m) is set, when the reservoir water level is lower than the water level, early warning is needed, and the risk of reservoir emptying is avoided.
Index three: the daily amplitude of reservoir water level. According to the regulation of reservoir A dispatching regulations, the daily variation of the reservoir water level of the reservoir A is not more than 3.0m/d (daily variation safety threshold), and when the daily variation exceeds the threshold, early warning is needed to avoid the unstable risk of the bank slope.
The index is four: considering the flood safety of the river channel at the downstream of the reservoir A, according to the regulation of the flood discharge early warning guide file of the reservoir A, when the basic flow and the increase and decrease flow values exceed the early warning threshold values, early warning is needed, and the risk of life and property loss at the downstream of the reservoir is avoided.
Index five: reservoir gate operating pressure. According to the dispatching characteristics of the reservoir A, the operation pressure of the reservoir gate is quantified in the following mode (see table 3 in detail).
Figure RE-GDA0001735091050000121
Figure RE-GDA0001735091050000131
Table 3 shows the quantitative manner of the operating pressure of the reservoir gate in the tung tree forest
And (3) analyzing the application effect:
the application effect of the invention is analyzed, and various safety evaluation indexes are shown in table 4. It can be seen that the present invention has a significant adverse impact on the improvement compared to the actual scheduling regime. Three unfavorable typical conditions were further analyzed in detail:
Figure RE-GDA0001735091050000132
table 1 shows the effect of the present invention
(1) Unfavorable working condition 1
Fig. 3 shows the daily warehousing flow and the whole plant output process, and the actual operation and the reservoir a power generation flow process of the method. The water regime log shows that the current day is moderate to heavy rain and local rainstorm, so the warehousing flow is large (the average warehousing flow is 2560m 3/s). It can be seen that: the output process of the hydropower station in the same day has high fluctuation, the current generation process lines of the two working conditions are basically consistent, and certain saw-tooth shapes are formed along with the non-uniformity of the output process.
Fig. 4 shows the actual operation condition of the day and the reservoir a water level process after the scheduling by the method of the invention. It can be seen that:
the water level of the reservoir is greatly influenced by the water, and the running range of the water level of the reservoir in the actual working condition is 1011.9-1014.4 m. Analyzing the change characteristic of reservoir water level in the day, wherein the reservoir water level slowly falls in the time period of 0: 00-08: 45, the reservoir water level rises obviously because the hydropower station does not generate power in the time period of 9: 00-9: 45, the reservoir water level keeps operating at a higher level all the time, and the reservoir water level (1014.3-1014.4 m) approaches the upper limit of early warning water level of 1014.5m in the time period of 21: 00-22: 45.
And secondly, under the scheduling of the method, the reservoir water level operation range is 1011.1-1012.5 m. And in the period of 0: 00-08: 45, the trend characteristics of the reservoir water level are consistent with those of the actual working conditions (both decrease), the reservoir water level still slowly decreases in the later period, reaches the lowest value 1011.1m of the reservoir water level in the current day at 18:45, then slowly rises, and keeps the same level with the reservoir water level in the actual working conditions at the end of the scheduling period in the current day.
The power generation flow process of the method is basically consistent with the actual working condition, but the reservoir water level processes of the power generation flow process and the actual working condition are obviously different. Thus, the overflow flow process of both needs to be further analyzed, see FIG. 5 for details. It can be seen that: in the time period of 9: 00-9: 45, when the current generation amount of the unit is 0, the method increases the overflow flow, which is the reason why the characteristics of reservoir water level change after the scheduling of the method are opposite to the actual working condition.
Figure RE-GDA0001735091050000141
Figure RE-GDA0001735091050000151
Table 5 shows statistical analysis of scheduling results of the present invention and actual conditions (unfavorable conditions one)
Table 5 is a statistical indicator value of the scheduling process, wherein the coefficient of variation is equal to the standard deviation divided by the mean value, and in statistics, the larger the value is, the larger the dispersion of the characterization data sequence is, and the larger the fluctuation of the characterization process is used here. It can be seen that:
the fluctuation of the power generation flow process of the invention is not greatly different from the fluctuation of the power generation flow process of the actual working condition, and is basically consistent with the fluctuation of the power generation flow process of the same day.
Compared with the actual working condition, the invention obviously improves the stability of the warehouse-out flow process by adjusting the overflow and abandon flow; when the product is put in storage, the coefficient of variation in the flow process is 0.08, and after the product is regulated by a reservoir, the coefficient of variation in the flow process is reduced to 0.01 when the product is taken out of the storage.
The analysis shows that on the basis of ensuring that the power generation task is smoothly completed on the same day, compared with the actual working condition, the method has the advantages that the regulation effect of the reservoir A is better exerted by calling the reservoir gate, the water flow form of the reservoir outlet flow is effectively improved, and the flood safety of the downstream river channel is improved.
(2) Unfavorable working condition two
Fig. 6 shows the load command process on the same day, which reaches 14 times of load commands in total, and the load rate is 36%. Fig. 7 shows the process of the warehouse-in flow and the warehouse-out flow at the same day, compared with the actual working condition, the stability of the process of the warehouse-out flow of the reservoir is improved to a certain extent.
The distribution process of the current generation flow and the overflow flow is further analyzed, as shown in fig. 8. In the aspect of the power generation flow process, the shape of the process line of the invention is similar to that of the process line of the actual working condition, and the goodness of fit is higher; in the aspect of the overflow flow process, the two operation modes have large fluctuation, and comparison analysis needs to be performed by combining statistical data.
In terms of reservoir scheduling safety, fig. 9 shows the reservoir water level change process of the tung forest reservoir according to the invention and the actual working conditions. In actual operation, in a period of 14: 00-15: 45, the reservoir level (1014.48-1014.54 m) approaches or even breaks through the upper limit of the early warning water level, and the scheduling work is tense; for the method provided by the invention, the reservoir water level does not break through the upper and lower limits of the early warning water level, and the pressure of the scheduling work of the reservoir in the same day is released to a certain extent.
Figure RE-GDA0001735091050000161
Table 6 shows statistical analysis of scheduling results of the present invention and actual conditions (unfavorable conditions two)
Table 6 shows the statistical indicator values of two scheduling procedures, which can be seen as follows:
the fluctuation of the invention and the fluctuation of the actual working condition power generation flow process are both larger, and the fluctuation of the invention and the daily output process are basically consistent.
Compared with the actual working condition, the invention reduces the fluctuation of the overflow and abandon flow process by using the gate operation, thereby reducing the fluctuation of the ex-warehouse flow.
Compared with the fluctuation (the coefficient of variation is 0.16) of the flow process during warehousing, the coefficient of variation of the flow process during ex-warehouse is reduced to 0.09 after the dispatching by the method, and the regulation and storage effect of the reservoir on the warehousing runoff is exerted.
The analysis shows that under the condition that the load changes temporarily in the current day process, the actual operation conditions have different degrees of adverse effects on the reservoir and the downstream river channel.
(3) Unfavorable working condition three
Table 7 shows statistics of data in the scheduling process of the present invention and the actual operating mode on the same day, fig. 10 and 11 show the power generation flow, the overflow flow, and the delivery flow processes of the two operating modes, and fig. 12 shows the reservoir water level change process of the two operating modes. And (4) analyzing by combining the charts:
Figure RE-GDA0001735091050000171
table 7 shows statistical analysis of scheduling results of the present invention and actual conditions (three unfavorable conditions)
The comparison of the scheduling results of the present invention and the actual conditions shown in table 7 shows that:
average daily warehouse entry flow rate is 4178m3And/s (the incoming water level is large, the process stability is high), the hydropower station unit is stopped in a period of 2: 00-6: 45, and the load rate is 58%.
Secondly, the scheduling strategy established by the invention ensures the stability of the flow out of the reservoir, the flood early warning times of the river channel are reduced from 5 times to 0 time, and the flood safety of the river channel is greatly improved.
Thirdly, in actual operation, in the time periods of 7: 00-7: 45 (reservoir level 1014.48m) and 22: 00-22: 45 (reservoir level 1014.57m), the reservoir level approaches or even breaks through the upper limit of the early warning level; the invention relieves the tension situation to a certain extent and provides space for reasonably controlling the reservoir water level.
The above analysis shows that the present invention plays a positive role in the day.
Details not described in this specification are within the skill of the art that are well known to those skilled in the art.

Claims (6)

1. A reservoir coordination scheduling strategy obtaining method considering river flood demand is characterized by comprising the following steps:
step one, acquiring a forecast value of reservoir warehousing runoff of a current day, an expected control value of reservoir water level at the end of the current day and a power generation command regulating output process of a hydropower station of the current day;
secondly, considering the hierarchical flood demand of the downstream river channel, and constructing a reservoir coordination scheduling model considering the flood demand of the river channel;
thirdly, integrating constraint conditions in a coordinated scheduling model on the basis of the objective function;
fourthly, solving the established reservoir coordination scheduling model by using a dynamic programming method;
fifthly, providing key indexes from two aspects of reservoir dispatching safety and downstream river flood safety, constructing a safety evaluation index system suitable for coordinating the characteristics of a dispatching model, and comprehensively evaluating and feeding back the safety of the dispatching strategy obtained by the fourth step; if the safety evaluation index requirements are not met, returning to the second step, adjusting the model parameters and recalculating;
the objective function is as follows:
Figure FDA0003233728940000011
in the formula: fabs (. cndot.) as a function of absolute value; qout,tIs waterThe flow of the warehouse out of the warehouse in the time period t; qout,t-1The flow rate of the reservoir in the time period t-1 is taken out; t is a scheduling period; lambda [ alpha ]tThe river flood early warning demand coefficient of the reservoir in the time period t is obtained; t is 1 to T;
wherein the demand coefficient lambdatThe grading acquisition mode is as follows:
Figure FDA0003233728940000021
in the formula: SET { indices } is the selected index SET; fi (-) is an early warning distinguishing mode corresponding to the ith level basic flow; λ i is a set constant, where i ═ 1,2, …, n; the concrete expression mode of the SET { indices } and fi (-) is determined by combining the actual characteristics of engineering.
2. The method for obtaining the reservoir coordination scheduling strategy considering the river flood demand according to claim 1, wherein the constraint conditions comprise:
a. water balance constraint
Vt+1=Vt+(It-Qpower,t-Qab,t)Δt (3)
In the formula: vt、Vt+1The initial storage capacity and the final storage capacity of the reservoir at the time t are respectively; i ist、Qpower,t、Qab,tRespectively the warehousing flow, the power generation flow and the overflow and abandon flow of the reservoir at the time t; Δ t is the period length;
b. water level restraint
Zt,min≤Zt≤Zt,max (4)
In the formula: zt、Zt,min、Zt,maxRespectively the water level, the lowest water level and the highest water level allowed by the reservoir in the time period t;
c. let-down flow restriction
qt,min≤Qpower,t+Qab,t≤qt,max (5)
In the formula: q. q.st,min、qt,maxRespectively is the lower discharge flow limit and the upper discharge flow limit of the reservoir in the time period t;
d. initial and final water level control
Z0=ZFirst stage,ZT=ZPowder (6)
In the formula: z0 and ZT are the initial and final water levels of the dispatching period of the reservoir respectively; zFirst stage,ZPowderThe water level values are respectively given at the beginning and the end of the day.
3. The method for obtaining the reservoir coordination scheduling strategy considering the river flood demand according to claim 1, wherein the fourth step comprises the following steps:
s1: dividing stages and determining stage variables, dividing a reservoir dispatching period into 96 stages according to a time interval of 15 minutes, regarding the stages as a multi-stage decision problem, and representing the stage variables by t; t is the facing time period, and T + 1-T is the remaining time period;
s2: determining state variable, selecting reservoir capacity as state variable, and recording Vt-1The initial reservoir capacity at the t-th time period, VtThe storage capacity of the reservoir at the end of the t period;
s3: determining decision variables and selecting reservoir discharge flow Qout,t=Qpower,t+Qab,tIs a decision variable; qpower,t、Qab,tRespectively representing the generating flow and the overflowing and abandoning flow of the reservoir at the time t;
s4: determining a state transition equation, namely a water quantity balance equation, an equation (3);
s5: determining stage indexes, namely river flood warning demand coefficient lambda t and reservoir delivery flow amplitude delta Q of each stageoutThe product of (a);
s6: determining an optimum function, i.e. ft *(Vt-1) Denotes an initial storage capacity V from the t-th staget-1Starting from the sum of the optimal water abandoning amplitude in the T-th time period, the obtained dynamically-planned reverse time sequence recursion equation is as follows:
Figure FDA0003233728940000041
in the formula: vt-1 m1Represents a bank capacity value V when the initial discrete point of the t-th period is m1t m2Denotes a storage capacity value at which the discrete point at the end of the t-th period (beginning of t +1 period) is taken as M2, and M1 is 0,1, …, M; m2 ═ 0,1, …, M; m is the discrete point number of the storage capacity; f. oft *(Vt-1 m1) Denotes the initial storage capacity V from the t-th periodt-1 m1Starting to sum of the optimal water abandoning amplitude in the T-th time period; f. of* t+1(Vt m2) Represents that the initial storage capacity is V from the t +1 th periodt m2Starting to sum of the optimal water abandoning amplitude in the T-th time period; f. of* T+1(VT m2) Represents that the initial storage capacity is V from the T +1 th periodT m2Starting to sum of the optimal water abandoning amplitude in the T-th time period; omega t is a decision variable and represents the warehouse-out flow QtStorage capacity Vt-1A decision set meeting various constraints of the reservoir of the hydropower station at a given time;
s7: performing reverse time sequence recursion calculation, carrying out forward time sequence recursion from the last stage to the first stage according to the formula (7), and solving a reverse time sequence process with the minimum target function of the reservoir in each time period in the whole dispatching period under the condition of meeting related constraint conditions;
s8: and (4) performing sequential recursive calculation to obtain the optimal strategy (Q) of the ex-warehouse flow process corresponding to the optimal valuetAnd corresponding reservoir capacity optimal state point value { V }t}。
4. The method for acquiring the reservoir coordination scheduling strategy considering the river flood demand according to claim 1, wherein the acquiring mode of the demand coefficient is as follows:
when the overflow flow variation is within the range of the flood demand of the river channel, early warning information does not need to be issued, the demand coefficient lambda is 1, and the dispatching objective is to pursue the stability of the flow process of leaving the reservoir;
when the scheduling strategy breaks through the early warning requirement of the overflow flow amplitude, starting an early warning program, and taking the overflow flow amplitude as the requirement coefficient;
the larger the demand coefficient is, the larger the overflow and abandon flow amplitude is, the more adverse to the downstream river flood situation is, so the demand for implementing early warning is higher, the pressure born by corresponding scheduling work is increased, and the pressure is avoided when an optimization decision is made.
5. The method for acquiring the reservoir coordination scheduling strategy considering the river flood demand according to claim 1, wherein the key indexes in the fifth step comprise: the number of times of early warning water level upper limit crossing of reservoir water level, the number of times of early warning water level lower limit crossing of reservoir water level, daily variation of reservoir water level, the number of times of early warning flood of river channel and the operation pressure of reservoir gate.
6. The method for acquiring the reservoir coordination scheduling strategy considering the river flood demand according to claim 5, wherein the safety evaluation index system in the fifth step is shown in the following table:
Figure FDA0003233728940000051
wherein, the priority level decreases from top to bottom.
CN201810326686.1A 2018-04-12 2018-04-12 Reservoir coordination scheduling strategy obtaining method considering river flood demand Active CN108537449B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810326686.1A CN108537449B (en) 2018-04-12 2018-04-12 Reservoir coordination scheduling strategy obtaining method considering river flood demand

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810326686.1A CN108537449B (en) 2018-04-12 2018-04-12 Reservoir coordination scheduling strategy obtaining method considering river flood demand

Publications (2)

Publication Number Publication Date
CN108537449A CN108537449A (en) 2018-09-14
CN108537449B true CN108537449B (en) 2021-11-12

Family

ID=63480173

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810326686.1A Active CN108537449B (en) 2018-04-12 2018-04-12 Reservoir coordination scheduling strategy obtaining method considering river flood demand

Country Status (1)

Country Link
CN (1) CN108537449B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109167367B (en) * 2018-09-19 2021-12-24 长江勘测规划设计研究有限责任公司 Dynamic hydropower station load adjustment method based on unilateral matching mechanism
JP7436291B2 (en) 2020-05-28 2024-02-21 株式会社日立パワーソリューションズ Time series optimization program and time series optimization system
CN111898253B (en) * 2020-07-15 2022-04-15 武汉大学 Reservoir dispatching and downstream river ecological environment protection cooperation value evaluation method
CN112199631B (en) * 2020-09-30 2023-07-21 昆明电力交易中心有限责任公司 Method, device and equipment for declarative electric quantity associated control of step power station
CN113506010A (en) * 2021-07-21 2021-10-15 华电云南发电有限公司 Large-scale watershed hydropower station flood discharge gate digital centralized regulation and control method
CN115712800B (en) * 2022-11-24 2023-07-28 国能大渡河流域水电开发有限公司 Reservoir water level fluctuation treatment method
CN116258278A (en) * 2023-05-10 2023-06-13 青岛研博数据信息技术有限公司 Method, system and equipment for deducing constructed water level reservoir capacity
CN117252406B (en) * 2023-11-20 2024-02-02 长江水利委员会长江科学院 Water replenishing scheduling method, device and medium facing downstream water taking requirement

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101303588A (en) * 2008-06-04 2008-11-12 东北大学 Method and system for automatically making steel-smelting continuous casting furnace sub batch plan
CN101599143A (en) * 2009-07-20 2009-12-09 武汉大学 A kind of dispatching zone of reservoir determinacy Optimization Dispatching is determined method and equipment thereof
CN101705671A (en) * 2009-11-19 2010-05-12 武汉大学 Yellow River upstream cascade hydroelectric station operation design and optimized dispatching method as well as equipment
CN102867275A (en) * 2012-08-14 2013-01-09 贵州乌江水电开发有限责任公司 Medium-term and long-term combined power generation optimal scheduling method and system in cascade reservoir group
CN103580020A (en) * 2013-03-07 2014-02-12 长沙理工大学 Method for multi-target dynamic optimal dispatching of power system with wind farm based on NSGA-II and Look-ahead
CN103810104A (en) * 2014-03-04 2014-05-21 中国人民解放军63863部队 Method and system for optimizing software test case
CN104181900A (en) * 2014-09-04 2014-12-03 徐雪松 Layered dynamic regulation method for multiple energy media
CN105956714A (en) * 2016-05-21 2016-09-21 华能澜沧江水电股份有限公司 Novel group searching method for optimal scheduling of cascade reservoir groups
CN107491635A (en) * 2017-07-26 2017-12-19 华北电力大学 A kind of step reservoir water sand joint optimal operation method based on nested dimension-reduction algorithm
CN107590556A (en) * 2017-08-31 2018-01-16 华北电力大学 A kind of flood control by reservoir regulation Optimization Scheduling implemented according to back water storage effect

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7627491B2 (en) * 2003-01-07 2009-12-01 Swiss Reinsurance Company Method for evaluating flood plain risks
US20120253746A1 (en) * 2011-03-28 2012-10-04 International Business Machines Corporation Autonomous data quality evaluation

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101303588A (en) * 2008-06-04 2008-11-12 东北大学 Method and system for automatically making steel-smelting continuous casting furnace sub batch plan
CN101599143A (en) * 2009-07-20 2009-12-09 武汉大学 A kind of dispatching zone of reservoir determinacy Optimization Dispatching is determined method and equipment thereof
CN101705671A (en) * 2009-11-19 2010-05-12 武汉大学 Yellow River upstream cascade hydroelectric station operation design and optimized dispatching method as well as equipment
CN102867275A (en) * 2012-08-14 2013-01-09 贵州乌江水电开发有限责任公司 Medium-term and long-term combined power generation optimal scheduling method and system in cascade reservoir group
CN103580020A (en) * 2013-03-07 2014-02-12 长沙理工大学 Method for multi-target dynamic optimal dispatching of power system with wind farm based on NSGA-II and Look-ahead
CN103810104A (en) * 2014-03-04 2014-05-21 中国人民解放军63863部队 Method and system for optimizing software test case
CN104181900A (en) * 2014-09-04 2014-12-03 徐雪松 Layered dynamic regulation method for multiple energy media
CN105956714A (en) * 2016-05-21 2016-09-21 华能澜沧江水电股份有限公司 Novel group searching method for optimal scheduling of cascade reservoir groups
CN107491635A (en) * 2017-07-26 2017-12-19 华北电力大学 A kind of step reservoir water sand joint optimal operation method based on nested dimension-reduction algorithm
CN107590556A (en) * 2017-08-31 2018-01-16 华北电力大学 A kind of flood control by reservoir regulation Optimization Scheduling implemented according to back water storage effect

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
水利水电工程综合效益数学模型研究;陈建明;《河海大学学报》;19920131;1-6页 *
水库调度与防洪度汛技术发展现状与展望;靳少波;《:水库调度与防洪度汛技术发展现状与展望》;20170228;全文 *

Also Published As

Publication number Publication date
CN108537449A (en) 2018-09-14

Similar Documents

Publication Publication Date Title
CN108537449B (en) Reservoir coordination scheduling strategy obtaining method considering river flood demand
CN107818385B (en) Method for predicting real-time operation trend of cascade hydropower station group
CN110472824B (en) Cascade hydropower station short-term multi-objective optimization scheduling method considering peak shaving requirements
CN101705671B (en) Yellow River upstream cascade hydroelectric station operation design and optimized dispatching method as well as equipment
CN112149984B (en) Reservoir flood regulation multidimensional uncertainty risk analysis method based on Bayesian network
CN110071530B (en) Wind power climbing coordination scheduling method of energy storage-containing power system based on LSTM
CN110472825B (en) Multistage scheduling mechanism coordinated step hydropower station real-time scheduling water abandoning reduction method
CN108710972B (en) Reservoir flood control scheduling method based on POA algorithm
CN112884232A (en) Reservoir ecological scheduling method considering inflow forecasting uncertainty difference of rich and low water periods
CN108090671B (en) Cross-basin load distribution method for multi-stage inter-domain-intra-domain feedback iteration
CN113112060A (en) Lifting pump station scheduling method and device based on soft flow measurement and computer equipment
CN107862408B (en) Minimum early warning coordinated rolling optimization method for water abandonment of hydraulic power plant
CN114722709A (en) Cascade reservoir group optimal scheduling method and system giving consideration to generated energy and minimum output
CN112215464A (en) Prediction balance scheduling system for blast furnace gas under multiple working conditions
CN106570786B (en) Reservoir adaptive scheduling method based on scheduling rule time-varying characteristics
CN102663517A (en) Annual maximum load prediction method based on engineering consultation industry expansion and temperature reduction model
CN116703134B (en) Multi-target scheduling method and system for large cross-river basin water diversion reservoir
CN106468879A (en) A kind of liquid level flow nonlinear area control method
CN114021902B (en) Reservoir dispatching method for dynamic planning dimension reduction based on dynamic rope collection and discrete mechanism
CN116011733B (en) Multi-scale cooperative control intelligent scheduling method and system for cascade hydropower station group
CN112418539A (en) Basin hydropower station cluster comprehensive linkage flood prevention method and system based on big data
CN115907402B (en) Method and system for evaluating joint guaranteed output of cascade hydropower station
CN107591847B (en) Method for adjusting Automatic Gain Control (AGC) of hydroelectric generating set by using variable parameter mode
CN115169825B (en) Method for determining drought alarm water levels of reservoirs with different adjusting capacities
CN112365126B (en) Large-scale water diversion and regulation project scheduling method applied to reservoir containing online storage and regulation

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