Summary of the invention
The object of the invention is to, a kind of year-end level of multi-year regulating storage reservoir for timed position prediction method and system are provided, solution procedure is simple, considered carry-over storage later scheduling later stage benefit at the year end, can improve the accuracy of Forecasting Methodology, compare with the method such as BP neural network, the present invention more is applicable to the long reservoir of Streamflow Data series.
For solving the problems of the technologies described above, the present invention adopts following technical scheme: a kind of year-end level of multi-year regulating storage reservoir for timed position prediction method may further comprise the steps:
S1 sets up year-end level of multi-year regulating storage reservoir for timed position prediction model, and it is stored in the model bank;
S2, predictive server is the Optimized Operation target to the maximum according to the step gross energy, calculates and selects corresponding year-end level of multi-year regulating storage reservoir for timed position prediction model;
S3, predictive server is called the model solution algorithm in the algorithms library, finds the solution year-end level of multi-year regulating storage reservoir for timed position prediction model, obtains the position rule of falling into water that disappears the year end.
In the aforesaid year-end level of multi-year regulating storage reservoir for timed position prediction method, year-end level of multi-year regulating storage reservoir for timed position prediction model comprises multiple goal coupling prediction model and the Successive Regression forecast model that excavates based on mathematical statistics.
In the aforesaid year-end level of multi-year regulating storage reservoir for timed position prediction method,
(1) objective function of multiple goal coupling prediction model is: E=Max (E
d+ E
s), in the formula:
E is step reservoir generating gross energy, E
dFor step reservoir is worked as annual electricity generating capacity, E
sBe the accumulation of energy at the year end of carry-over storage,
Be objective function Ob:
In the formula: n is whole step reservoir number, and m is carry-over storage number in the step reservoir, and T is the period sum, K
iBe the comprehensive power factor of the reservoir that is numbered i, Q
I, tBe the reservoir that the is numbered i generating flow at period t, H
I, tBe the reservoir that the is numbered i generating net head at period t, A
jFor coefficient is A
j=K
j/ 3600, V
NjBe the carry-over storage that the is numbered j storage capacity at the year end,
For the carry-over storage that is numbered j holds the average productive head that stays storage capacity at the year end, if this carry-over storage downstream still has reservoir, then
Should be lower reservoir and hold the average productive head sum of staying storage capacity;
(2) constraint condition of multiple goal coupling prediction model comprises:
A water balance constraint: V (i, t+1)=V (i, t)+(Q
I(i, t)-Q
O(i, t)) * Δ t,
In the formula: V (i, t), V (i, t+1) represent respectively first, the last storage capacity of i reservoir t period; Q
I(i, t), Q
O(i, t) represents respectively warehouse-in and the outbound flow of i reservoir t period;
B flow equilibrium constraint: Q
I(i+1, t)=Q
O(i, t)+q (i, t),
In the formula: q (i, t) expression i and the local inflow of i+1 reservoir t period;
C water storage level constraint: Z
Min(i, t)≤Z (i, t)≤Z
Max(i, t),
In the formula: Z
Max(i, t), Z
Min(i, t) represents that respectively i reservoir t period allows the upper and lower limit of water level;
D vent flow constraint: Q
Omin(i, t)≤Q
O(i, t)≤Q
Omax(i, t),
In the formula: Q
Omax(i, t), Q
Omin(i, t) represents respectively the upper and lower limit of i reservoir vent flow;
E output of power station constraint: N
Min(i, t)≤N (i, t)≤N
Max(i, t),
In the formula: N
Max(i, t), N
Min(i, t) represents that respectively i reservoir t period allows the upper and lower limit of exerting oneself.
In the aforesaid year-end level of multi-year regulating storage reservoir for timed position prediction method, application progressively optimize-approaches one by one dynamic programming (DPSA-POA) hybrid algorithm and finds the solution the multiple goal coupling prediction model.
In the aforesaid year-end level of multi-year regulating storage reservoir for timed position prediction method, the finding the solution of Successive Regression forecast model of excavating based on mathematical statistics is based on the long serial combined optimization adjusting result of calculation of Cascade Reservoirs, use the stepwise regression analysis method, set up the nonlinear function of year-end level of multi-year regulating storage reservoir for timed position correlation factor, and it is carried out reasonableness test.
In the aforesaid year-end level of multi-year regulating storage reservoir for timed position prediction method, stepwise regression method is when fitting function, according to the level of confidence of drafting, utilize the F check, operation is introduced or kicked out of to each factor, thereby automatically select optimum factor and make up fitting function, have larger elasticity and fitting precision.
Realize a kind of year-end level of multi-year regulating storage reservoir for timed position prediction system of preceding method, comprise forecast model storehouse, predictive server and algorithms library, predictive server is provided with the Model Selection module and the database that predicts the outcome, the Model Selection module is connected with forecast model storehouse, algorithms library respectively, algorithms library is connected with the database that predicts the outcome, wherein
The forecast model storehouse is used for storage year-end level of multi-year regulating storage reservoir for timed position prediction model;
The Model Selection module is used for predictive server and is the Optimized Operation target to the maximum according to the step gross energy, selects corresponding year-end level of multi-year regulating storage reservoir for timed position prediction model;
Algorithms library is used for storage and finds the solution year-end level of multi-year regulating storage reservoir for timed position prediction model;
The database that predicts the outcome is used for the position rule of falling into water that disappears at year end that the storage solving model obtains.
In the aforesaid year-end level of multi-year regulating storage reservoir for timed position prediction system, also be provided with long serial combined optimization on the predictive server and regulate computing module, long serial combined optimization is regulated computing module and is connected with algorithms library, the database that predicts the outcome respectively, be used for the long serial Streamflow Data of watershed and calculate, obtain long serial combined optimization and regulate result of calculation.
Compared with prior art, the present invention utilizes disappear when falling into water the position year end of multiple goal coupling prediction model research carry-over storage, can process preferably the generated energy of reservoir and the relation between the reservoir accumulation of energy contradiction target, this method has been taken into account reservoir generated energy then and the energy that stored in the reservoir in follow-up several years.But owing to the unknown to the water situation, the multi-objective predictive model does not fully take into account a few years from now on water to the impact of step reservoir scheduling, therefore, the present invention adopts the method for statistical law to set up year-end level of multi-year regulating storage reservoir for timed position prediction model, the Successive Regression forecast model that namely excavates based on mathematical statistics.Regulation o f reservoir operation is with concrete summarizes and the expression of historical summary to the reasonable moving law research of reservoir, can be used for instructing the following operation of reservoir.Based on the thinking of obtaining scheduling graph, consider to have storehouse group's cooperation of carry-over storage participation, can set up really qualitative mathematics model of multi-reservoir combined dispatching, then grow series and regulate calculating, obtain the long-term optimal operation plan of multi-reservoir, wherein also comprise the reasonable change rule of year-end level of multi-year regulating storage reservoir for timed position, it is carried out analytical calculation, then can draw control (prediction) rule (or function) of the drowning position that disappears the year end.Solution procedure of the present invention is simple, and has considered carry-over storage later scheduling later stage benefit at the year end, so that the accuracy of Forecasting Methodology has improved 5 ~ 8%.Because the methods such as the Forecasting Methodology that the present invention relates to and BP neural network are compared, and do not need long serial runoff sample data training network and then draft network parameter, therefore comparatively speaking, the present invention is more practical to the long reservoir of runoff data system.
Embodiment
Embodiments of the invention: the position prediction method of falling into water that disappears at the end of a kind of carry-over storage (Hong Jiadu of Wujiang River Basin, two reservoirs in Goupitan) as shown in Figure 1, may further comprise the steps:
S1 sets up year-end level of multi-year regulating storage reservoir for timed position prediction model, and it is stored in the model bank;
S2, predictive server is the Optimized Operation target to the maximum according to the step gross energy, calculates and selects corresponding year-end level of multi-year regulating storage reservoir for timed position prediction model;
S3, predictive server is called the model solution algorithm in the algorithms library, finds the solution year-end level of multi-year regulating storage reservoir for timed position prediction model, obtains the position rule of falling into water that disappears the year end.
Year-end level of multi-year regulating storage reservoir for timed position prediction model comprises multiple goal coupling prediction model and the Successive Regression forecast model that excavates based on mathematical statistics.
1, multiple goal coupling prediction model:
(1) objective function of multiple goal coupling prediction model is: E=Max (E
d+ E
s), in the formula:
E is step reservoir generating gross energy, E
dFor step reservoir is worked as annual electricity generating capacity, E
sBe the accumulation of energy at the year end of carry-over storage,
Be objective function Ob:
In the formula: n is whole step reservoir number, and m is carry-over storage number in the step reservoir, and T is the period sum, K
iBe the comprehensive power factor of the reservoir that is numbered i, Q
I, tBe the reservoir that the is numbered i generating flow at period t, H
I, tBe the reservoir that the is numbered i generating net head at period t, A
jFor coefficient is A
j=K
j/ 3600, V
NjBe the carry-over storage that the is numbered j storage capacity at the year end,
For the carry-over storage that is numbered j holds the average productive head that stays storage capacity at the year end, if this carry-over storage downstream still has reservoir, then
Should be lower reservoir and hold the average productive head sum of staying storage capacity;
(2) constraint condition of multiple goal coupling prediction model comprises:
A water balance constraint: V (i, t+1)=V (i, t)+(Q
I(i, t)-Q
O(i, t)) * Δ t,
In the formula: V (i, t), V (i, t+1) represent respectively first, the last storage capacity of i reservoir t period; Q
I(i, t), Q
O(i, t) represents respectively warehouse-in and the outbound flow of i reservoir t period;
B flow equilibrium constraint: Q
I(i+1, t)=Q
O(i, t)+q (i, t),
In the formula: q (i, t) expression i and the local inflow of i+1 reservoir t period;
C water storage level constraint: Z
Min(i, t)≤Z (i, t)≤Z
Max(i, t),
In the formula: Z
Max(i, t), Z
Min(i, t) represents that respectively i reservoir t period allows the upper and lower limit of water level;
D vent flow constraint: Q
Omin(i, t)≤Q
O(i, t)≤Q
Omax(i, t),
In the formula: Q
Omax(i, t), Q
Omin(i, t) represents respectively the upper and lower limit of i reservoir vent flow;
E output of power station constraint: N
Min(i, t)≤N (i, t)≤N
Max(i, t),
In the formula: N
Max(i, t), N
Min(i, t) represents that respectively i reservoir t period allows the upper and lower limit of exerting oneself.
Application progressively optimize-approaches one by one dynamic programming (DPSA-POA) hybrid algorithm and finds the solution the multiple goal coupling prediction model.Its solving result and analysis on its rationality are as follows:
Consider affect disappear the year end the fall into water factor of position of reservoir and mainly contain two: one is water level at the beginning of the reservoir, and another is exactly year to put a runoff in storage.At first, get identical warehouse-in runoff, analyze different beginning of the year water level to the impact of position of falling into water that disappears the year end of Hong Jiadu, goupitan reservoir; Secondly, get water level at the beginning of the identical reservoir, the warehouse-in runoff of analyzing different frequency is on the reservoir impact of position of falling into water that disappears the year end.In calculating, be to reduce the phase mutual interference between Hong Jiadu, two carry-over storages in Goupitan, cross reservoir in research flood man and disappear the year end when falling into water the position rule, the water level at the whole story of goupitan reservoir is decided to be 620m; Simultaneously, disappear at the end of the research goupitan reservoir when falling into water the position rule, Jiang Hongjia crosses reservoir water level at the whole story and is decided to be 1124m.
(1) impact of warehouse-in runoff
For analyzing the different runoff process of frequencies of putting in storage to the impact of year-end level of multi-year regulating storage reservoir for timed position, get flood man cross reservoir, Goupitan at the beginning of water level be respectively 1100m, 608m, use the DPSA algorithm above-mentioned Optimal Operation Model found the solution calculating, the result as table 1,2 and Fig. 2,3 shown in.
Station different frequency warehouse-in runoff process result of calculation table (unit: water level/m, energy/hundred million kWh) crosses in table 1 flood man
Frequency |
5% |
10% |
20% |
30% |
40% |
50% |
60% |
70% |
80% |
90% |
95% |
The year end water level |
1125 |
1120 |
1118 |
1113 |
1112 |
1110 |
1104 |
1098 |
1095 |
1091 |
1088 |
Energy |
381.3 |
358.7 |
355.4 |
334 |
325.6 |
300.3 |
290.2 |
273.5 |
266 |
209.7 |
198.5 |
Table 2 Goupitan different frequency warehouse-in runoff process result of calculation table (unit: water level/m, energy/hundred million kWh)
Frequency |
5% |
10% |
20% |
30% |
40% |
50% |
60% |
70% |
80% |
90% |
95% |
The year end water level |
620 |
619 |
618 |
617 |
616 |
615 |
614 |
612 |
608 |
602 |
600 |
Energy |
360 |
338.9 |
335.4 |
319.5 |
313.4 |
293.9 |
282.1 |
269.5 |
258.8 |
200.6 |
169.4 |
As shown in Table 1, flood man crossed reservoir and disappeared the year end and to fall into water the position generally about 1120m the high flow year, was not less than 1110m, played the effect of " holding rich "; Normal flow year generally about 1110m, is not less than water level at the beginning of the year; Low flow year generally is not higher than about 1100m, is not less than level of dead water 1076m, plays the effect of " mending withered ".
(2) impact analysis of water level at the beginning of
Discrete water level at the beginning of two carry-over storages, it is 50% corresponding 1992 warehouse-in runoff process that the warehouse-in runoff is got frequency, use DPSA and calculate, result of calculation as table 3,4 and Fig. 4,5 shown in.
The position result of calculation (P=50%) (unit: water level/m, energy/hundred million kWh) of falling into water that disappears at the end of the water level the different beginning of the years is crossed by table 3 flood man
The beginning of the year water level |
1080 |
1082 |
1084 |
1086 |
1088 |
1090 |
1092 |
1094 |
1096 |
1098 |
1100 |
The year end water level |
1085 |
1085 |
1086 |
1087 |
1089 |
1091 |
1094 |
1094 |
1099 |
1110 |
1111 |
Energy |
298.7 |
299.9 |
301.3 |
302.5 |
304 |
305.4 |
307.1 |
308.6 |
310.1 |
311.7 |
313.9 |
The beginning of the year water level |
1102 |
1104 |
1106 |
1108 |
1110 |
1115 |
1120 |
1125 |
1130 |
1135 |
1140 |
The year end water level |
1112 |
1113 |
1116 |
1117 |
1118 |
1121 |
1128 |
1130 |
1131 |
1136 |
1139 |
Energy |
316.6 |
319.1 |
321.1 |
323.3 |
325.6 |
331.3 |
337.4 |
343.5 |
349 |
356.6 |
364.2 |
The position result of calculation (P=50%) of falling into water disappears at the end of the water level at the beginning of the difference of table 4 Goupitan
(unit: water level/m, energy/hundred million kWh)
Under warehouse-in runoff and other equal conditions, from table 3 and Fig. 4,5 as can be known: if water level was lower than 1110m at the beginning of reservoir crossed in flood man, then its susceptibility that falls into water position influence to disappearing the year end is not clearly, and along with raising of water level at the beginning of the reservoir crossed by flood man, the position institute of falling into water of disappearing the year end is influenced larger.This mainly is because what adopt when calculating is normal flow year, the warehouse-in runoff belongs to not rich not withered situation, owing to being subject to the restriction that the step assurance is exerted oneself and minimum load retrains, if the beginning of the year, water level was lower, the water consumption rate of its unit of electrical energy is just larger, and corresponding water consumption is just larger, and if the beginning of the year water level higher, corresponding water consumption is less, and the water yield of saving just can increase its productive head again.Therefore the scheduling end of term, for having the leading reservoir of emphasizing to save performance in the step reservoir, begins have higher water level in schedule periods, for can be got back to high water stage and established certain basis.The goupitan reservoir gross energy is along with constantly the raising of water level at the beginning of the year as shown in Table 4, presents the trend of increase.
(3) Hong Jiadu, the goupitan reservoir position prediction result that falls into water that disappears the year end
Use above-mentioned model, can draw Hong Jiadu, the different warehouse-in of goupitan reservoir runoff and the position prediction result that falls into water that disappears at the end of the water level the different beginning of the years, shown in table 5,6.
The reservoir position prediction table (unit: m) that falls into water that disappears the year end crosses in table 5 flood man
Annotate: P is warehouse-in flow frequency (%).
Position prediction table (the unit: m) that falls into water disappears at the end of table 6 goupitan reservoir
Annotate: P is warehouse-in flow frequency (%).
The position multiple goal coupling of falling into water for disappearing at the end of Hong Jiadu, the Goupitan in the his-and-hers watches 5,6 predicts the outcome and carries out analysis on its rationality, and analog computation is carried out in this research application 1951.5 ~ 2007.4 totally 56 older series materials.At first according to the water level at the beginning of the year of the water frequency of calculating the time and Hong Jiadu, goupitan reservoir, fall into water according to determining in the table 5,6 that Hong Jiadu, goupitan reservoir disappear the year end thus; Use the maximum model of step reservoir generated energy again, carry out that optimal operation of cascade reservoirs calculates in year, for many years average generated energy, assurance exert oneself and generate electricity fraction and design load compares to step power station, and be as shown in table 7.
The long serial examination table of table 7 year-end level of multi-year regulating storage reservoir for timed position multiple goal coupling prediction model
As can be seen from the above table, the step after model calculates for many years average generated energy is 292.45 hundred million kWh, and comparing design load 289.7 hundred million kWh has increased by 0.9%, and the assurance of calculating is exerted oneself has increased by 8.37% for 2417.2MW compares design load 2230.4MW.Therefore, the drowning position rule that disappears at the end of the Hong Jiadu under this rule, the goupitan reservoir is rationally, reliably, can use in production practices.
2, the Successive Regression forecast model that excavates based on mathematical statistics:
It is based on the long serial combined optimization of Cascade Reservoirs and regulates result of calculation, uses the stepwise regression analysis method, sets up the nonlinear function of year-end level of multi-year regulating storage reservoir for timed position correlation factor, and it is carried out reasonableness test.
Stepwise regression method according to the level of confidence of drafting, utilizes the F check when fitting function, operation is introduced or kicked out of to each factor, makes up fitting function thereby automatically select optimum factor.The relevant initial regression vectors in position of falling into water of disappearing the year end is crossed by flood man to be had: Hong Jiadu then reservoir inflow and following several years reservoir inflow, Hong Jia cross water level at the beginning of the year, introduction cross reach then following several years outbound flow, Hong Jiadu-east wind-introduction crosses the following several years runoff reach of reservoir; The fall into water initial regression vectors of position of disappearing at goupitan reservoir year end has: goupitan reservoir is the water level and the outbound flow of following several years warehouse-in runoff (comprising the reservoir runoff reach) and Da Hua water at the beginning of the year of water level, the following several years runoff reach of Wu Jiangdu-goupitan reservoir, lower reservoir (Silin) at the beginning of reservoir inflow and following several years reservoir inflow, the goupitan reservoir then.The position influence factor of falling into water of disappearing the year end is crossed by flood man: flood man cross the beginning of the year water level and square, First Year reservoir inflow square, Second Year reservoir inflow, the 3rd year reservoir inflow square; The position influence factor of falling into water of disappearing at Goupitan year end is: at the beginning of the goupitan reservoir water level and square, Goupitan First Year reservoir inflow, Goupitan Second Year reservoir inflow, Goupitan the 4th year and the 5th year reservoir inflow.
When using the maximum model of step generated energy to grow series to calculate, for reduce as far as possible reservoir the whole story water level on the impact of regression model, this returns and calculates the adjusting calculating achievement that adopted 1952 ~ 2002 years and analyze.Choosing respectively level of confidence is 0.1 and 0.001, sets up the Successive Regression model, and result of calculation is as shown in table 8.Result in the associative list, and consider that the unitarity of upstream and downstream reservoir, this forecast model set up that unified to adopt level of confidence be 0.1, obtain Hong Jiadu and Goupitan Successive Regression forecast model respectively suc as formula 1, shown in the formula 2.
The different degree of confidence regression effects of table 8 relatively
Z
Year disappears=7144.898-11.06Z
1+ 5.0359 * 10
-3Z
1 2+ 5.5068 * 10
-4Q
1 2-6.2918 * 10
-3Q
2+ 1.2885 * 10
-4Q
3 2(formula 1)
Z
Year disappears=4320.053-11.96Z
1+ 0.01Z
1 2+ 0.1Q
1-0.04Q
2-0.06Q
4+ 0.03Q
5(formula 2)
Wherein: Z
Year disappearsBe the position of falling into water that disappears at year end of reservoir; Z
1Be water level at the beginning of this reservoir; Q
1, Q
2, Q
3, Q
4, Q
5For forecasting the water yield and forecast in rear 5 years then, this reservoir comes the water yield.
The forecast model solving result of Hong Jiadu, Goupitan is seen respectively Fig. 6,7.Flood man is crossed, two kinds of model predictions of goupitan reservoir disappear falls into water the position comparison diagram respectively shown in Fig. 8,9.
When the concrete use procedure of model, have the following suggestion can be for reference:
(1) because Goupitan year-end level of multi-year regulating storage reservoir for timed position influence many factors shows fully and statistical analysis method is very difficult, therefore can preferentially select the multi-objective predictive model;
(2) can find from Fig. 8, the predicting the outcome of two kinds of models that carry-over storage crosses in flood man is more or less the same, but for the low flow year, Statistical Prediction Model often predicted value is lower, cause like this step generated energy to increase, the adjusting function of considering the downstream goupitan reservoir is better, and therefore when water was more withered, the drowning position that disappears the year end is crossed by flood man can pay the utmost attention to Statistical Prediction Model;
(3) for the high flow year of Hong Jiadu, middle water year, two kinds of model prediction results are then comparatively close.In the actual use procedure, the drowning position that disappears the year end of Hong Jiadu, two carry-over storages in Goupitan can be predicted respectively by two kinds of models, then according to the actual conditions of producing, controls in conjunction with yardman's operating experience.
Realize a kind of year-end level of multi-year regulating storage reservoir for timed position prediction system of preceding method, as shown in figure 10, comprise forecast model storehouse 1, predictive server 4 and algorithms library 3, predictive server 4 is provided with Model Selection module 2 and the database 5 that predicts the outcome, Model Selection module 2 is connected with forecast model storehouse 1, algorithms library 3 respectively, algorithms library 3 is connected with the database 5 that predicts the outcome, wherein
Forecast model storehouse 1 is used for storage year-end level of multi-year regulating storage reservoir for timed position prediction model;
Model Selection module 2 is used for predictive server 4 and is the Optimized Operation target to the maximum according to the step gross energy, selects corresponding year-end level of multi-year regulating storage reservoir for timed position prediction model;
Algorithms library 3 is used for storage and finds the solution year-end level of multi-year regulating storage reservoir for timed position prediction model;
The database 5 that predicts the outcome is used for the position rule of falling into water that disappears at year end that the storage solving model obtains.
Also be provided with long serial combined optimization on the predictive server 4 and regulate computing module 6, long serial combined optimization is regulated computing module 6 and is connected with algorithms library 3, the database 5 that predicts the outcome respectively, be used for the long serial Streamflow Data of watershed and calculate, obtain long serial combined optimization and regulate result of calculation.