CN106815656B - Method for acquiring cascade reservoir energy storage dispatching diagram - Google Patents
Method for acquiring cascade reservoir energy storage dispatching diagram Download PDFInfo
- Publication number
- CN106815656B CN106815656B CN201611266466.1A CN201611266466A CN106815656B CN 106815656 B CN106815656 B CN 106815656B CN 201611266466 A CN201611266466 A CN 201611266466A CN 106815656 B CN106815656 B CN 106815656B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- time interval
- energy storage
- initial
- dispatching
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 98
- 238000004146 energy storage Methods 0.000 title claims abstract description 81
- 238000010586 diagram Methods 0.000 title claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 90
- 230000008569 process Effects 0.000 claims abstract description 41
- 230000008859 change Effects 0.000 claims description 15
- 230000003247 decreasing effect Effects 0.000 claims description 12
- 238000010187 selection method Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000010248 power generation Methods 0.000 abstract description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 26
- 238000011144 upstream manufacturing Methods 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000005457 optimization Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 230000003321 amplification Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 241000283707 Capra Species 0.000 description 2
- 241000283086 Equidae Species 0.000 description 2
- 241001494479 Pecora Species 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
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
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a step reservoir energy storage dispatching diagram obtaining method, which belongs to the field of hydropower energy optimized operation and power system power generation optimized dispatching and comprises an upper basic dispatching line and a lower basic dispatching line obtaining step, a typical runoff obtaining step corresponding to the upper basic dispatching line and the lower basic dispatching line and an increased and reduced output line obtaining step. The method for acquiring the cascade reservoir energy storage dispatching diagram overcomes the problems that the determination of the time interval initial discrimination coefficient is difficult and the calculation of the typical runoff process is difficult in the reverse time sequence recursion calculation process, and has better rationality and effectiveness.
Description
Technical Field
The invention belongs to the field of optimal operation of hydroelectric energy and optimal scheduling of power generation of a power system, and particularly relates to a discrimination coefficient-based backward calculation method for a cascade reservoir energy storage scheduling graph.
Background
The discrimination coefficient method is widely applied to the cascade reservoir combined dispatching due to the clear physical significance, particularly, the method is mature in research and application in the aspect of combined application with a cascade energy storage dispatching diagram, can not only fully play the upstream and downstream compensation functions of the cascade reservoir, but also can effectively avoid some defects when the discrimination coefficient method is singly applied. However, the existing research mainly focuses on the aspect of manufacturing the stepped energy storage dispatching diagram based on various optimization algorithms, such as a genetic algorithm, a particle swarm algorithm, a step-by-step optimization algorithm and the like, and for a discriminant coefficient method, the existing results are only used in the simulation dispatching after the energy storage dispatching diagram is obtained, and are not unified in the manufacturing and application of the energy storage dispatching diagram. And the research results of obtaining the cascade energy storage dispatching graph by adopting typical year hydrological runoff series through reverse time recursive calculation based on a discriminant coefficient method are less. Although the method for obtaining the cascade energy storage dispatching graph through the optimization algorithm is simple and convenient, the optimization algorithm is usually not tightly combined with the actual background of the problem, certain physical significance is lacked, and the reliability of the optimization result is not high.
Different from a single reservoir dispatching diagram, the factors considered in the step energy storage dispatching diagram manufacturing process are numerous, the drawing process is more difficult and complicated, and particularly in the process of carrying out inverse time sequence recursion calculation by adopting typical annual runoff series, the two difficulties are as follows:
(1) the reverse time sequence recursion calculation of each time interval is to combine the last state of the system time interval with a discrimination coefficient and time interval output to calculate the initial state of the system time interval, the discrimination coefficient used in the calculation is a value corresponding to the initial state of the time interval, but the initial state of the system is unknown when the reverse time sequence recursion calculation is started, so that the discrimination coefficient (size relationship) corresponding to the initial state of the time interval is difficult to obtain at the moment, and a mature method for solving the problem does not exist at present.
(2) In the process of calculating the typical runoff corresponding to the upper basic dispatching line and the lower basic dispatching line, the initial judgment coefficient of each time period of the system is unknown, and the incoming flow size of the cascade system is the amount to be calculated, so that the calculation of enlarging and reducing the output line of the cascade energy storage dispatching graph is extremely difficult. At present, no effective method for solving the problem exists, and therefore, the application of the method for calculating the cascade energy storage dispatching graph by adopting a reverse time sequence recursion mode is greatly limited.
Patent document CN102080366A discloses a method for making a cascade reservoir joint dispatching graph. The method comprises the steps of taking a long-series historical runoff process as input, taking ten days as a calculation time interval, ensuring output according to the design of each hydropower station, and making each reservoir dispatching diagram according to a single-reservoir dispatching diagram making method step by step from upstream to downstream according to an equal-output method. The invention takes the output indicated by the cascade joint dispatching diagram as coordination, can play the role of compensation dispatching of the cascade reservoir to a certain extent and improve the utilization of water energy. However, the method for making the cascade reservoir united dispatching graph disclosed by the document has the following defects:
(1) the invention sets up each reservoir dispatching diagram according to the method for making the single-reservoir dispatching diagram from upstream to downstream step by the equal output method, considers the water quantity relation between the upstream and the downstream of the cascade reservoir, but fails to fully consider the water head relation, and the energy stored in the cascade hydroelectric system is related to the water quantity and the water head of each power station in the system at the same time, so the reservoir dispatching diagram obtained by the invention can not fully exert the whole energy benefit of the cascade system.
(2) When the watershed cascade hydroelectric system performs combined dispatching, the electric power system usually only gives out the total output of the system, and when the dispatching graph obtained by the invention is used for carrying out actual dispatching operation of the cascade reservoirs, each reservoir can determine one output according to the dispatching graph of the reservoir, so that the problem that the total output of the given system of the power grid is not equal to the sum of the outputs of all cascade power stations exists, and the difficulty in distributing the outputs among the cascade reservoirs is easily caused.
Disclosure of Invention
Aiming at the defects or the improvement requirements of the prior art, the invention provides a reverse calculation method based on a discrimination coefficient to obtain a cascade reservoir energy storage scheduling graph, and aims to provide a total flow for calculating the cascade reservoir energy storage scheduling graph based on a water storage and supply discrimination coefficient and a reverse time sequence recursion mode.
In order to achieve the aim, the invention provides a discrimination coefficient-based cascade reservoir energy storage dispatching diagram back-pushing calculation method, which comprises an upper basic dispatching line and a lower basic dispatching line obtaining step, a typical runoff obtaining step corresponding to the upper basic dispatching line and the lower basic dispatching line and an increased output-reduced line obtaining step,
firstly, according to the typical year selection method in the acquisition of the single-base dispatch graph, selecting N-year typical hydrological runoff series from long series historical runoff data,
then, from the last period at the end of the scheduling period, a judgment coefficient method is used to ensure that the output mode is used to obtain a step total energy storage change curve corresponding to each typical year runoff process through reverse time recursive calculation,
and finally, taking an upper envelope of each typical annual step total energy storage change curve as an upper basic dispatching line and taking the following envelopes as lower basic dispatching lines, wherein in the reverse time sequence recursion calculation, the method for determining the time interval initial discrimination coefficient comprises the following steps:
s1: calculating all possible permutation combinations of the size sorting of the discrimination coefficients of the banks according to the number of the cascade reservoirs, wherein one permutation combination corresponds to one reservoir number combination,
s2: assuming a time interval initial reservoir number combination, then carrying out reverse time sequence recursion calculation on the time interval initial reservoir number combination to obtain a time interval initial state of the reservoir system, finally calculating each reservoir discrimination coefficient corresponding to the time interval initial state of the reservoir system, and sequencing the reservoir numbers by the each reservoir discrimination coefficient to obtain another time interval initial reservoir number combination,
s3: if the initial reservoir number combination in the other time interval is consistent with the assumed initial reservoir number combination in the time interval, the calculation of the current time interval is finished, the size of each reservoir discrimination coefficient corresponding to the assumed initial reservoir number combination in the time interval is taken as the initial discrimination coefficient in the current time interval,
otherwise, the number combination of the initial reservoirs in another time interval is assumed, and the calculation is carried out again until the number combination of the initial reservoirs in the assumed time interval is consistent with the number combination of the initial reservoirs in another time interval obtained by calculation.
Further, in the typical runoff obtaining step corresponding to the upper and lower basic dispatching lines,
firstly, supposing that the total inflow of the cascade reservoir system in the current time period is TQ, distributing the total inflow of the TQ to each power station according to the annual average inflow proportion of each reservoir in the cascade reservoir system to obtain the supposed typical runoff corresponding to each reservoir in the cascade reservoir system,
then, the initial time interval discrimination coefficient determining method in the upper and lower basic dispatching line obtaining step is adopted to obtain the initial time interval state by the end time interval state, the total energy storage E of the cascade corresponding to the initial time interval state is calculated, the total energy storage E' of the cascade corresponding to the initial time interval state is read from the upper and lower basic dispatching lines,
finally, comparing the total energy storage E with the total energy storage E ', if the total energy storage E and the total energy storage E' are equal, finishing the calculation of the current time interval, taking the assumed typical runoff corresponding to each reservoir in the cascade reservoir system as the typical runoff corresponding to each reservoir in the actual cascade reservoir system,
otherwise, the assumed system total incoming flow TQ is updated with the difference E-E ', and the calculation is performed again until E ═ E',
and repeating the calculation in a sequential mode until the typical runoff corresponding to each time period of each reservoir in the cascade reservoir system is obtained.
Furthermore, in the step of obtaining the increased and decreased output lines, the set increased and decreased output values are used as output values in the equal output calculation, the reverse time recursion calculation is carried out based on a discriminant coefficient method, and corresponding increased and decreased output lines are obtained,
in the reverse time-series recursive computation, the time interval initial discrimination coefficient determination method is the same as the time interval initial discrimination coefficient determination method in the upper and lower basic scheduling line acquisition steps.
In general, compared with the prior art, the above technical solution contemplated by the present invention can achieve the following beneficial effects:
(1) the invention provides a discrimination coefficient-based reverse time recursive calculation method for a cascade reservoir energy storage dispatching graph. Compared with a combined dispatching diagram obtained directly through an optimization algorithm, the method provided by the invention is based on the reservoir water storage and supply discrimination coefficient and the typical annual hydrological runoff series, so that the obtained cascade energy storage dispatching diagram has better physical background significance and higher reliability, and more importantly, the method well solves the problems of time period initial discrimination coefficient determination and typical runoff process deduction in the reverse time sequence recursion calculation process of the energy storage dispatching diagram, so that the technology is more mature, and the application prospect is wider.
(2) The invention takes the cascade reservoir of the Li Xian river basin in China as an example, the cascade reservoir energy storage dispatching diagram is obtained by the method, and long-series simulation calculation is carried out, and the calculation result shows that the method is superior to the traditional combined dispatching method in the aspects of total generated energy and guaranteed output under the premise that the guarantee rate is not reduced, the increase in the generated energy is 0.26%, the increase in the guaranteed output is particularly obvious and can reach 6.8%. Therefore, the simulation result well verifies the rationality and effectiveness of the method.
Drawings
FIG. 1 is a flow chart of upper and lower basic scheduling line calculation and determination of time interval initial discrimination coefficients;
fig. 2 is a diagram of a typical runoff calculation process corresponding to an upper basic dispatching line and a lower basic dispatching line;
FIG. 3 is a general flow chart of a step reservoir energy storage dispatching diagram backward calculation method based on the discrimination coefficient according to the embodiment of the invention;
FIG. 4 is a trend evolution diagram of the generated energy and the guarantee rate related to the step reservoir energy storage dispatch graph backward calculation method based on the discrimination coefficient according to the embodiment of the present invention along with the change of the guaranteed output;
FIG. 5 is a LiXianjiang cascade reservoir energy storage dispatching graph which meets the requirement of guarantee rate and has maximum power generation quantity and is related to a cascade reservoir energy storage dispatching graph backward calculation method based on a discrimination coefficient of the embodiment of the invention;
FIG. 6 is a diagram of a process of the annual average water level of a cliff goat mountain reservoir related to a method for calculating the backward push of a cascade reservoir energy storage dispatching diagram based on a discrimination coefficient according to the embodiment of the invention;
fig. 7 is a process diagram of the multi-year average water level of the rockfill reservoir, which relates to the discrimination coefficient-based method for calculating the inverse of the cascade reservoir energy storage dispatching diagram according to the embodiment of the invention;
fig. 8 is a process diagram of the mean water level of the dragon horse reservoir over many years, which relates to a discrimination coefficient-based step reservoir energy storage dispatching diagram back-stepping calculation method.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
The invention obtains the cascade reservoir energy storage dispatching diagram by selecting the typical annual runoff series and performing reverse time recursive calculation, and well solves the problems of time interval initial discrimination coefficient determination and typical runoff process deduction in the reverse time recursive calculation process of the energy storage dispatching diagram. The method mainly comprises the steps of calculating an upper basic dispatching line and a lower basic dispatching line, and drawing parts for calculating and increasing the reduced output line in the typical incoming flow process corresponding to the upper basic dispatching line and the lower basic dispatching line. The method specifically comprises the following steps:
(1) upper and lower basic dispatch line calculation
Firstly, selecting a typical hydrological runoff series from long-series historical runoff data according to a typical year selection method when a single-base dispatch graph is drawn; then, from the last period at the end of the scheduling period, combining a discriminant coefficient method to obtain a step total energy storage change curve corresponding to each typical annual runoff process through reverse time recursive calculation in an equal output (guaranteed output) mode; and finally, taking an upper envelope of each typical annual step total energy storage change curve as an upper basic dispatching line and taking a lower envelope as a lower basic dispatching line.
For the aforementioned problem that the time interval initial discrimination coefficient cannot be effectively determined in the reverse time-series recursive calculation, the invention provides a solution by taking a three-base cascade system as an example. Specifically, assuming that the cascade reservoir is numbered (0,1,2) in sequence from upstream to downstream, if the reservoirs are numbered according to the magnitude of the discrimination coefficients of the first time interval, there are six combinations of (0,1,2), (0,2,1), (1,0,2), (1,2,0), (2,0,1) and (2,1,0), and in the reverse time-series recursion calculation, it is the magnitude sequence of the discrimination coefficients of the first time interval, not the values of the discrimination coefficients themselves, that is, all that is required is. Therefore, the initial state of the cascade reservoir system can be calculated in a trial mode in a traversing combination mode. That is, firstly, the reservoir number combination corresponding to the initial state of one time interval is assumed, then reverse time order recursion calculation is carried out on the combination to obtain the initial state of the system time interval, finally, the discrimination coefficients of all the reservoirs corresponding to the initial state are calculated, and the reservoir numbers are sequenced according to the discrimination coefficients, so that the initial reservoir number combination of another time interval is obtained. If the calculated combination is consistent with the assumed combination, the calculation of the current time interval is finished, otherwise, another combination is assumed, and the calculation is carried out again until the assumed combination is equal to the calculated combination.
The general flow of the calculation of the upper and lower basic scheduling lines is shown in fig. 1, fig. 1 is a flow chart of the calculation of the upper and lower basic scheduling lines and the determination of the initial time interval coefficient thereof, and the calculation flow of the solution of the initial time interval coefficient determination problem is shown in the left part of fig. 1.
Specifically, in fig. 1, N-year typical runoff data starts to be input, and a back-thrust calculation is performed from the T-th period by using the ith (i ═ 1,2, …, N) typical runoff data: assuming that the current time interval is the t-th time interval, and assuming a magnitude sequence of discrimination coefficients K of each primary step library in the t-th time interval, wherein the magnitude sequence of the K value is combined with the time interval incoming flow, the t-th time interval end water level and the guaranteed output value, iterative calculation is carried out to obtain the time interval primary library capacity of each step library, so that another K value magnitude sequence can be calculated, and if the calculated K value magnitude relation is consistent with the assumption, the calculation in the current time interval is finished to obtain the total energy storage of the primary step system in the current time interval; making t equal to t-1, and entering calculation of the next time interval until t equal to 1, and finishing the reverse calculation aiming at the ith typical year to obtain a step total energy storage (initial time interval) change process; similarly, similar inverse calculation is carried out on other typical annual runoff data to obtain different step total energy storage change curves, and finally, an upper envelope line and a lower envelope line are taken to obtain an upper basic dispatching line and a lower basic dispatching line.
(2) Typical incoming flow process estimation corresponding to upper and lower basic scheduling lines
Before the increase and decrease output lines are calculated, a typical incoming flow process corresponding to the upper and lower basic scheduling lines needs to be determined first to serve as a typical incoming flow process in the calculation of the increase and decrease output lines. In the typical incoming flow process calculation, not only the size of the incoming flow is a to-be-solved quantity, but also the initial discrimination coefficient of each time period of the cascade system is an unknown quantity, so that the operation is difficult to realize.
The present invention provides the following solutions:
firstly, assuming that the total flow of the cascade system in the current time period is TQ, and distributing the total flow to each power station through the average incoming flow proportion of each cascade library for many years; then, by combining a discrimination coefficient processing method used in the calculation of the basic dispatching line, the initial state of the system time period is calculated according to the end state of the system time period, the total energy storage E of the cascade corresponding to the initial state of the time period is calculated, and the total energy storage E' of the other cascade corresponding to the initial state of the time period is read from the basic dispatching line; and finally, comparing the total energy storage E with the total energy storage E ', if the total energy storage E and the total energy storage E' are equal, finishing the calculation of the current time interval, updating the assumed total system incoming flow TQ by the difference E-E ', and calculating again until the E is equal to the E'.
The calculation flow of the method is shown in fig. 2, fig. 2 is a typical runoff calculation process diagram corresponding to upper and lower basic scheduling lines, and the calculation flow for solving the typical runoff calculation process corresponding to the upper and lower basic scheduling lines is shown in the left part of fig. 2.
Specifically, in fig. 2, the total energy storage variation curve of the cascade system corresponding to the upper and lower basic lines starts to be input, and the back-stepping calculation is performed from the T-th time period: supposing that the current time period is t and the total inflow rate of the steps in the current time period is TQ, distributing the TQ to each reservoir of the steps according to the inflow distribution proportion of each reservoir of the steps obtained by long series historical runoff data statistics; on the basis, a size sequence of discrimination coefficients K of each primary step library in the t-th time period is assumed, the size sequence of the K values is combined with the allocated time period incoming flow, the water level at the end of the t-th time period and the guaranteed output value, the time period primary library capacity of each primary step library is obtained through iterative calculation, so that another K value size sequence can be calculated, and if the calculated K value size relationship is consistent with the assumption, the time period primary library capacity of each primary step library can be obtained, so that a total energy storage E of the primary step system in the current time period can be obtained; in addition, another total energy storage E ' of the cascade system at the beginning of the current time interval can be found out from a known basic dispatching line, if E is equal to E ', the calculation of the current time interval is finished, t is equal to t-1, the calculation of the next time interval is carried out, otherwise, the total inflow TQ of the cascade at the current time interval is assumed again, and the calculation is carried out again until E is equal to E '; and finally, when t is 1, all time interval calculation is finished, and a typical incoming flow process of each time interval of each reservoir of the steps corresponding to the upper basic dispatching line or the lower basic dispatching line can be obtained.
(3) Increase and decrease of output line calculation
After typical incoming flow processes corresponding to the upper basic dispatching line and the lower basic dispatching line are obtained, the set increasing and decreasing output values are respectively used as output values in equal output calculation, reverse time sequence recursion calculation is carried out by combining a discriminant coefficient method, and the corresponding increasing and decreasing output lines can be obtained, wherein a time interval initial discriminant coefficient processing method in the calculation process is consistent with the time for calculating the upper basic dispatching line and the lower basic dispatching line. The total flow of the cascade reservoir energy storage dispatching graph reverse time sequence recursion calculation method based on the reservoir water storage and supply discrimination coefficient is shown in the attached figure 3. Fig. 3 is a general flow chart of a step reservoir energy storage dispatching diagram back-stepping calculation method based on a discrimination coefficient according to an embodiment of the present invention, which includes three parts of drawing of upper and lower basic dispatching lines, typical runoff calculation corresponding to the basic dispatching lines, and drawing of increasing and decreasing output lines.
Specifically, in fig. 3, two typical runoff processes (corresponding to the upper and lower basic dispatching lines, respectively) corresponding to each time period of each library of the cascade system obtained by the method described in fig. 2 are calculated by reverse thrust from the T-th time period: assuming that the current time interval is the t-th time interval, and assuming a magnitude sequence of discrimination coefficients K of each base of the primary cascade in the t-th time interval, wherein the magnitude sequence of the K value is combined with the time interval incoming flow, the t-th time interval end water level and a corresponding increased output (or decreased output) value, and the time interval primary base capacity of each base of the cascade can be obtained through iterative calculation, so that another magnitude sequence of the K value can be calculated, if the magnitude relation of the calculated K value is consistent with the assumption, the calculation in the time interval is finished, and the total energy storage of the primary cascade system in the current time interval is obtained; and (3) making t equal to t-1, calculating the next time interval until t equal to 1, finishing the reverse calculation, and obtaining a step total energy storage (initial time interval) change process, wherein the process is a corresponding increase or decrease output line.
The cascade energy storage dispatching graph is manufactured by using a reverse time order recursion calculation method, long series simulation dispatching calculation is carried out, and the result is compared and analyzed with the traditional combined dispatching method so as to show the effect achieved by the method.
The Li Xian Jiangjiang river basin is located in the Yunnan province of China, seven hydropower stations are planned and built on the main stream, cliff sheep mountains, rock canes, New plain villages, dragon horses, Jupu ferres, Golan beaches and native karivers are sequentially arranged, the cliff sheep mountains, the rock canes and the dragon horses have adjusting performance, and the three adjusting type hydropower stations are selected as research objects in the embodiment. The normal water storage levels of the three hydropower stations are 835m, 756m and 639m respectively, the dead water levels are 818m, 740m and 605m respectively, and the design power generation guarantee rates are all 95%. Flood limit water levels of the cliff goat reservoir and the rock threshold reservoir are 818m and 740m respectively, and the flood limit water level of the dragon horse reservoir is 639 m.
The implementation steps of the invention are as follows:
the method comprises the following steps: selecting typical year runoff series from long series historical runoff data to ensure that output is taken as output of each time period, and calculating initial energy storage of each time period of the cascade system in a reverse time sequence from the last time period of a scheduling period by combining a discriminant coefficient method so as to obtain an energy storage change curve of the cascade system. The upper envelope and the lower envelope of each typical annual energy storage change curve are respectively taken as the upper basic dispatching line and the lower basic dispatching line of the energy storage dispatching graph, and the specific calculation flow is shown in figure 1.
Step two: and (3) calculating a corresponding typical incoming flow process by taking the obtained upper and lower basic dispatching lines as a system energy storage evolution curve and combining a judgment coefficient method and a guaranteed output value, wherein the specific flow is shown in fig. 2.
Step three: according to the obtained typical incoming flow process and the corresponding increasing and decreasing output values, each increasing and decreasing output process line is calculated in a reverse time sequence from the last time period of the scheduling period by combining a discriminant coefficient method, the process is basically consistent with the reverse calculation process when the upper and lower basic scheduling lines are manufactured in the step one, and the specific process is shown in the rightmost part of the figure 3.
The results obtained according to the process of the invention are as follows:
when the output change coefficients of the energy storage dispatching diagram are 1.2,1.1,1,0.9 and 0.8, a process curve of the gradient total power generation and the change of the guarantee rate along with the change of the guarantee output can be obtained, as shown in fig. 4. Since the power generation guarantee rate of the Lixianjiang cascade reservoir is not lower than 95%, as can be seen from FIG. 4, the total guaranteed output of the cascade is 115100kW, the annual average power generation amount of the cascade is 23.56 hundred million kWh, and the corresponding cascade energy storage scheduling graph is shown in FIG. 5.
In order to fully embody the implementation effect of the present invention, the above result is compared and analyzed with the result of the conventional joint scheduling. The traditional combined dispatching adopts a single-warehouse dispatching diagram of each cascade warehouse to carry out dispatching operation from top to bottom, and the power generation guarantee rate of each power station is not lower than 95%. The guaranteed output of the obtained three storehouses of the sheepson mountain, the rock threshold and the dragon horse is 23200kW,23406kW and 61200kW respectively, the sum of the guaranteed output is 107806kW, and the average annual power generation amount of the cascade is 23.500 hundred million kWh. Therefore, compared with the traditional combined dispatching, the energy storage dispatching diagram obtained in the invention has the advantages that the output is ensured to be increased by 7294kW, the amplification is 6.8%, and the total stepped power generation is increased by 0.060 hundred million kWh, and the amplification is 0.26%.
Therefore, under the same power generation guarantee rate constraint, the obtained cascade energy storage dispatching graph is superior to the traditional combined dispatching method in both the aspects of power generation and guaranteed output, and although the amplification of the cascade energy storage dispatching graph is small in the aspect of power generation, the amplification of the cascade energy storage dispatching graph in the aspect of guaranteed output is very obvious. Therefore, the energy storage dispatching diagram obtained by the invention can improve the power supply stability and reliability of the cascade power generation system to a great extent, thereby improving the power supply quality and having remarkable economic and social benefits.
In addition, according to the obtained step energy storage dispatching diagram, the processes of the average water levels of the three reservoirs of the cliff mountain, the rock threshold and the dragon horse can be obtained for many years as shown in fig. 6, 7 and 8 respectively. It can be seen from the three figures that the upstream reservoir often discharges water first to generate electricity, so that the water quantity of the upstream reservoir fully utilizes the high water head of the downstream reservoir to generate electricity, and the downstream reservoir discharges water backwards to generate electricity, so as to maintain high water level operation, and fully improve the electricity generation benefit of the water quantity of the upstream reservoir. Therefore, the water level process of the simulated scheduling of the cascade banks is consistent with the actual scheduling operation principle, and the effectiveness and the reasonability of the method provided by the invention are further verified.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (3)
1. A method for acquiring a cascade reservoir energy storage dispatching diagram is characterized by comprising an upper basic dispatching line and a lower basic dispatching line acquiring step, a typical runoff acquiring step corresponding to the upper basic dispatching line and the lower basic dispatching line and an increased output-reduced line acquiring step,
firstly, according to the typical year selection method in the acquisition of the single-base dispatch graph, selecting N-year typical hydrological runoff series from long series historical runoff data,
then, from the last period at the end of the scheduling period, a judgment coefficient method is used to ensure that the output mode is used to obtain a step total energy storage change curve corresponding to each typical year runoff process through reverse time recursive calculation,
and finally, taking an upper envelope of each typical annual step total energy storage change curve as an upper basic dispatching line and taking the following envelopes as lower basic dispatching lines, wherein in the reverse time sequence recursion calculation, the method for determining the time interval initial discrimination coefficient comprises the following steps:
s1: calculating all possible permutation combinations of the size sorting of the discrimination coefficients of the banks according to the number of the cascade reservoirs, wherein one permutation combination corresponds to one reservoir number combination,
s2: assuming a time interval initial reservoir number combination, then carrying out reverse time sequence recursion calculation on the time interval initial reservoir number combination to obtain a time interval initial state of the reservoir system, finally calculating each reservoir discrimination coefficient corresponding to the time interval initial state of the reservoir system, and sequencing the reservoir numbers by the each reservoir discrimination coefficient to obtain another time interval initial reservoir number combination,
s3: if the initial reservoir number combination in the other time interval is consistent with the assumed initial reservoir number combination in the time interval, the calculation of the current time interval is finished, the size of each reservoir discrimination coefficient corresponding to the assumed initial reservoir number combination in the time interval is taken as the initial discrimination coefficient in the current time interval,
otherwise, the number combination of the initial reservoirs in another time interval is assumed, and the calculation is carried out again until the number combination of the initial reservoirs in the assumed time interval is consistent with the number combination of the initial reservoirs in another time interval obtained by calculation.
2. The method for obtaining the cascade reservoir energy storage dispatching diagram as claimed in claim 1, wherein in the step of obtaining the typical runoff corresponding to the upper and lower basic dispatching lines,
firstly, supposing that the total inflow of the cascade reservoir system in the current time period is TQ, distributing the total inflow of the TQ to each power station according to the annual average inflow proportion of each reservoir in the cascade reservoir system to obtain the supposed typical runoff corresponding to each reservoir in the cascade reservoir system,
then, the initial time interval discrimination coefficient determining method in the upper and lower basic dispatching line obtaining step is adopted to obtain the initial time interval state by the end time interval state, the total energy storage E of the cascade corresponding to the initial time interval state is calculated, the total energy storage E' of the cascade corresponding to the initial time interval state is read from the upper and lower basic dispatching lines,
finally, comparing the total energy storage E with the total energy storage E ', if the total energy storage E and the total energy storage E' are equal, finishing the calculation of the current time interval, taking the assumed typical runoff corresponding to each reservoir in the cascade reservoir system as the typical runoff corresponding to each reservoir in the actual cascade reservoir system,
otherwise, the assumed system total incoming flow TQ is updated with the difference E-E ', and the calculation is performed again until E ═ E',
and repeating the calculation in this way until the typical runoff corresponding to each time period of each reservoir in the cascade reservoir system is obtained.
3. The method for obtaining the cascade reservoir energy storage dispatching map as claimed in claim 2, wherein in the step of obtaining the increased and decreased output lines,
taking the set increasing and decreasing output values as output values in equal output calculation, carrying out reverse time order recursion calculation based on a discriminant coefficient method to obtain corresponding increasing and decreasing output lines,
in the reverse time-series recursive computation, the time interval initial discrimination coefficient determination method is the same as the time interval initial discrimination coefficient determination method in the upper and lower basic scheduling line acquisition steps.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611266466.1A CN106815656B (en) | 2016-12-31 | 2016-12-31 | Method for acquiring cascade reservoir energy storage dispatching diagram |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611266466.1A CN106815656B (en) | 2016-12-31 | 2016-12-31 | Method for acquiring cascade reservoir energy storage dispatching diagram |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106815656A CN106815656A (en) | 2017-06-09 |
| CN106815656B true CN106815656B (en) | 2020-07-10 |
Family
ID=59110767
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611266466.1A Expired - Fee Related CN106815656B (en) | 2016-12-31 | 2016-12-31 | Method for acquiring cascade reservoir energy storage dispatching diagram |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106815656B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107808041B (en) * | 2017-10-16 | 2020-11-17 | 南京南瑞继保电气有限公司 | Method for optimally drawing cascade reservoir hydraulic connection diagram |
| CN108181938B (en) * | 2017-12-07 | 2020-12-18 | 南京南瑞继保电气有限公司 | Water level adjusting method for reversely dissolving waste water of hydropower station |
| CN110348692B (en) * | 2019-06-13 | 2022-03-18 | 华中科技大学 | Large-scale series-parallel reservoir group multi-target energy storage scheduling graph calculation method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102080366A (en) * | 2011-01-01 | 2011-06-01 | 国网电力科学研究院 | Method for drawing joint scheduling graph of step reservoir |
| 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 |
| CN106087847A (en) * | 2016-06-16 | 2016-11-09 | 长江勘测规划设计研究有限责任公司 | Water-supply reservoir scheduling graph dynamic control method based on Runoff Forecast |
-
2016
- 2016-12-31 CN CN201611266466.1A patent/CN106815656B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102080366A (en) * | 2011-01-01 | 2011-06-01 | 国网电力科学研究院 | Method for drawing joint scheduling graph of step reservoir |
| 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 |
| CN106087847A (en) * | 2016-06-16 | 2016-11-09 | 长江勘测规划设计研究有限责任公司 | Water-supply reservoir scheduling graph dynamic control method based on Runoff Forecast |
Non-Patent Citations (2)
| Title |
|---|
| Instructional Mutation Ant Colony Algorithm in Application of Reservoir Operation Chart Optimization;Shan Yu等;《2011 Fourth International Symposium on Knowledge Acquisition and Modeling》;20120123;第1-4页 * |
| 水库调度图优化研究;张铭等;《武汉大学学报（工学版）》;20040630;第37卷(第3期);第5-7页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106815656A (en) | 2017-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106682810B (en) | Long-term operation method of cross-basin cascade hydropower station group under dynamic production of giant hydropower station | |
| CN112633578B (en) | Cascade reservoir group optimal scheduling method under influence of diversion and water diversion engineering | |
| CN105719091B (en) | A kind of parallel Multiobjective Optimal Operation method of Hydropower Stations | |
| CN106355292B (en) | Cascade reservoir optimal scheduling method and system based on quantum particle swarm optimization | |
| Zhang et al. | Optimal operation of large-scale cascaded hydropower systems in the upper reaches of the Yangtze River, China | |
| CN106815656B (en) | Method for acquiring cascade reservoir energy storage dispatching diagram | |
| CN107248015B (en) | Reservoir real-time water storage scheduling method based on ensemble prediction | |
| CN110163420B (en) | Multi-target ecological scheduling method and system based on decomposition culture evolution algorithm | |
| CN103580020B (en) | A kind of based on NSGA-II and Look-ahead containing wind energy turbine set power system multiobjective Dynamic Optimization dispatching method | |
| CN103971174A (en) | Hydropower station group optimized dispatching method based on improved quantum-behaved particle swarm algorithm | |
| CN112036632B (en) | Optimal scheduling method based on cascade reservoir ecological power generation multi-target medium and long term random scheduling model | |
| CN106355285B (en) | Power load prediction method based on parameter correction | |
| CN110222938A (en) | A kind of Hydropower Stations head relation cooperative optimization method and system | |
| CN104182806B (en) | A kind of GROUP OF HYDROPOWER STATIONS Optimization Scheduling based on orthogonal dimensionality reduction searching algorithm | |
| CN104504455B (en) | A kind of lower GROUP OF HYDROPOWER STATIONS Long-term Optimal Dispatch method of step accumulation of energy control | |
| CN109472463B (en) | Multi-power-grid peak regulation scheduling method and system for hydropower station group in extra large basin | |
| CN108537370A (en) | Especially big basin water station group Optimization Scheduling based on hybrid intelligent dimension-reduction algorithm | |
| CN104166887A (en) | Orthogonal discrete differential dynamic programming method for cascade hydropower station group joint optimization scheduling | |
| CN108564231B (en) | Agent optimization dimension reduction method for large-scale hydropower station reservoir group joint scheduling | |
| CN108683174B (en) | Network loss optimization method based on multi-time scale demand response model | |
| CN109635999B (en) | Hydropower station scheduling method and system based on particle swarm-bacterial foraging | |
| CN104408531B (en) | A kind of uniform dynamic programming method of multidimensional multistage complicated decision-making problems | |
| CN111697607A (en) | Multi-terminal flexible direct-current transmission receiving-end power grid access method and system | |
| CN104182909A (en) | Multi-core parallel successive approximation method of hydropower system optimal scheduling | |
| CN105207233B (en) | Based on the idle work optimization method that Metropolis Hastings are combined with PSO |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200710 Termination date: 20211231 |