CN114897419A - Cascade power station load distribution method based on margin method - Google Patents
Cascade power station load distribution method based on margin method Download PDFInfo
- Publication number
- CN114897419A CN114897419A CN202210618163.0A CN202210618163A CN114897419A CN 114897419 A CN114897419 A CN 114897419A CN 202210618163 A CN202210618163 A CN 202210618163A CN 114897419 A CN114897419 A CN 114897419A
- Authority
- CN
- China
- Prior art keywords
- output
- power station
- distribution
- time period
- power
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 229910052739 hydrogen Inorganic materials 0.000 claims 1
- 239000001257 hydrogen Substances 0.000 claims 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000007115 recruitment Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a load distribution method of a cascade power station based on an allowance method, which comprises the following steps: step S1: acquiring initial distribution results of output of each power station in each time period; step S2: finding out a time period when the total output per hour is equal to the sum of the lower limits of the outputs of all the power stations, and distributing all the power stations according to the lowest output; finding out a time period when the total output per hour is equal to the sum of the upper output limits of all the power stations, and distributing all the power stations according to the highest output; step S3: calculating the output values violating the upper output limit and the lower output limit of each power station in each time period, and constructing a violating constraint matrix; step S4: distributing the output force of each power station in all time periods within the time periods, so that all elements violating the constraint matrix in the distributed result are 0; step S5: and distributing the output of each power station in each time period, so that the daily average output of each distributed power station is the same as the daily average output in the initial distribution result, and obtaining a final distribution result.
Description
Technical Field
The invention relates to the field of hydroelectric power energy systems, in particular to a load distribution method for a cascade power station based on an allowance method.
Background
The current cascade hydropower station makes up the limitation of a single reservoir through combined dispatching and compensation adjustment, and more effectively improves the hydropower benefit. In order to make the cascade hydropower stations better perform the compensation and adjustment function among all the hydropower stations, research on the load distribution technology of the hydropower stations needs to be carried out.
The existing method for distributing loads of cascade power stations realizes multiple adjustments of the cascade output process by limiting the maximum daily electric quantity deviation of each hydropower station, considers the constraint of minimum water abandon risk and maximum cascade energy storage, and introduces an electric quantity deviation rate index, but the method has incomplete constraint conditions, the determination of the deviation index needs to be determined according to the historical operating condition, the adjustment process is complicated, and the time for actual operation is long. In addition, a calculation method realized by nesting and coupling daily total electric quantity distribution and an operation scheme is established, the effect of the empirical coefficient on distribution of daily total electric quantity among stations is replaced, and the calculation is realized by optimization calculation, but the load distribution method carries out ideal assumption on output in a peak-valley period, the change of the output of the power station from time period to time period is not considered, and the obtained result is slightly different from the actual operation of the power station. Moreover, by providing a water abandoning correction strategy based on peak shaving response and a load balancing reduction strategy in a low-ebb period, the check analysis and the output adjustment of the power station in a fixed dispatching mode are realized, but load distribution is required to be carried out according to the capacity classification of the cascade power station during application, and the practical use is inconvenient.
The load distribution method of the cascade power stations does not relate to a method for distributing output power among the cascade power stations by time periods. The method is beneficial to providing beneficial reference for the formulation of the dispatching scheme of the cascade hydropower station.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a load distribution method of a cascade hydropower station based on an allowance method, which provides reference for a scheduling scheme system of the cascade hydropower station, thereby better exerting the comprehensive benefit of cascade shop recruitment.
The invention is realized by the following technical scheme:
a load distribution method of a cascade power station based on an allowance method comprises the following steps:
step S1: acquiring initial distribution results of output of each power station in each time period;
step S2: finding out time periods when the total output per hour is equal to the sum of the lower limits of the outputs of all the power stations in the initial distribution result, and distributing all the power stations according to the lowest output in the time periods; finding out time periods when the total output per hour is equal to the sum of the upper output limits of all the power stations, and distributing each power station according to the highest output in the time periods to obtain secondary distribution results of the output of each power station in each time period;
step S3: calculating the output values violating the upper output limit and the lower output limit of each power station in each time period, and constructing a violating constraint matrix;
step S4: distributing the output of each power station in all time intervals in the secondary distribution result within the time intervals, so that all elements violating the constraint matrix in the distributed result are 0, and obtaining a tertiary distribution result;
step S5: and distributing the output of each power station in each time period, so that the daily average output of each distributed power station is the same as the daily average output in the initial distribution result, and obtaining a final distribution result.
In the technical scheme, when the small total output in the initial distribution result is equal to the output lower limit or the output upper limit of each power station, the output of each power station is ensured to be not lower than the output lower limit or exceed the output upper limit according to the output lower limit or the output upper limit of each power station; redistributing the output values of the power stations in a time period which is lower than the lower output limit or higher than the upper output limit, so that the output values which violate the constraint do not appear in the distributed results; and further distributing the average daily output value of each power station to be equal to the average daily output value in the initial distribution result.
Further, the method for calculating the hourly output distribution data of each power station initiated in step S1 includes: and multiplying the distribution weight of the output of each reservoir input from the outside by the total output of the step per hour.
Further, the step S3 specifically includes:
step S31: calculating the output value violating the upper limit and the lower limit of the output of each power station in each time period, and specifically comprising the following steps:
in the formula:
N w，i，t the output value violating the upper limit and the lower limit of the output of the ith hydropower station in the tth period is i-1 … n, and T-1 … T, and the unit is ten thousand kW; n is a radical of i，t The output of the ith hydropower station in the t time period is ten thousand kW; n is a radical of down，i，t The lower limit of the output of the ith hydropower station in the t period is ten thousand kW; n is a radical of up，i，t The output upper limit of the ith hydropower station in the t time period is ten thousand kW;
s32: the violating constraint matrix is formed by the violating output values of the upper limit and the lower limit of the output of each power station in each time period, and specifically comprises the following steps:
in the formula: n is a radical of w，all Is a violation of the constraint matrix.
Further, the step S4 specifically includes the following steps:
step S41: obtaining the output distribution result of each power station in each time interval in the secondary distribution result as a vector to be operated, and calculating the element sum (N) of violation constraint matrix of all power stations in each time interval w，1，t +N w，2，t +…+N w，n，t ) (ii) a If the element sum is not equal to 0, go to step S42; if the element sum is equal to 0, go to step S44;
step S42: adjusting the output value of each power station in each time interval to ensure that the output value of each power station in each time interval meets the output constraint:
calculating the sum of the output values which are greater than the upper output limit:
in the formula: n is a radical of wk The sum of the output values is greater than the upper output limit; k is a power station serial number which is greater than or equal to the upper limit of the output, and K is 1, 2. Alpha is a serial number set of the power station which is greater than or equal to the upper limit of the output; calculating the sum of the output values smaller than the lower output limit:
in the formula: n is a radical of wq Is the sum of the output values less than the lower output limit; q is a power station serial number smaller than or equal to the lower limit of output, and Q is 1, 2. Beta is a serial number set of the power station which is less than or equal to the lower limit of the output;
calculating the sum of the remaining force values violating the force constraints:
N wl ＝N wk +N wq
in the formula: n is a radical of wl The sum of the remaining output values violating the output constraint;
adjusting the power plant output for this period:
step S43: calculating the sum of the elements based on the adjusted power plant output result, and if the value is not equal to 0, skipping to step S32 again; if this value is equal to 0, go to step S34;
step S44: if t < 24, t is t +1, and the process proceeds to step S31.
Further, in step S44, if t is 24, the process ends.
Further, if the element sum is equal to 0 in step S41, the flow proceeds to step S44.
Further, the step S5 specifically includes:
step S51: calculating the daily average output of each power station according to the initial distribution structure, and subtracting the daily average output of each power station from the daily average output of the result of S4 to obtain a difference vector:
[N d，1 ，N d，2 ，…，N d，n ]＝[N o，1 ，N o，2 ，…，N o，n ]-[N a，1 ，N a，2 ，…，N a，n ]
in the formula: n is a radical of d，i The difference value of the original daily average output of the ith hydropower station and the daily average output of the result of S4 is obtained; n is a radical of o，i The average output of the ith hydropower station in the original day; n is a radical of a，i The daily average output for the result of S4 for the ith hydropower station;
step S52: and finding out a time period in which the absolute values of the differences between the output of all the power stations and the respective output upper limit and output lower limit are greater than a preset buffer value in the three-time distribution result:
the output vector of the power station in the t-th period is N i，t ，N 2，t ，…，N n，t ]If the tth period is an allocation period, the following is satisfied:
in the formula: delta is a preset buffer value, and the unit is ten thousand kW;
step S53: equally dividing the difference vector by the number of the distribution time interval, and adding the equally divided vector into the power station output vector of the distribution time interval:
[N 1，z ，N 2，z ，…，N 3，z ]＝[N 1，z ，N 2，z ，…，N 3，z ]+[N d，1 ，N d，2 ，…，N d，n ]z∈γ
in the formula: [ N ] 1，z ，N 2，z ，…，N n，z ]The power station output vector of the distribution time interval is obtained, and z is the serial number of the distribution time interval; gamma is a sequence number set of the distribution time interval;
step S54: calculating whether the output of the power station in the distribution time period after distribution exceeds the upper limit and the lower limit; if yes, reducing the buffer value, and turning to step S52; if not, the allocation is ended.
Compared with the prior art, the invention has the beneficial effects that:
the invention provides a load distribution method for a cascade power station based on an allowance method, which can distribute loads to the cascade power station on the premise that the percentage value of the output distribution result of each hydropower station to the cascade output is as close as possible to a given weight, and the distribution result meets the following three constraint conditions:
constraint 1. the daily average output of each station in the distribution result is equal to the daily average output obtained by the given distribution weight;
constraint 2. the step force of each time interval in the distribution result is equal to the given step force of each time interval;
and 3, constraining the output of each power station in the distribution result to be less than the constraint of an upper limit and a lower limit.
The method not only considers various constraint conditions, but also has high efficiency and convenient use, and the obtained result meets the actual engineering requirements, thereby providing a beneficial reference for the formulation of the cascade hydropower station scheduling scheme.
Drawings
FIG. 1 is a schematic diagram of a method according to an embodiment of the invention;
FIG. 2 is a schematic diagram of an allocation of elemental sums to negative values according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of the distribution of elemental sums to positive values according to an embodiment of the present invention.
Detailed Description
The technical solutions of the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings, and it is to be understood that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without any inventive step, are within the scope of the present invention.
As shown in fig. 1, the present embodiment provides a method for allocating a load of a cascaded power station based on an allowance method, which specifically includes the following steps:
s1: obtaining distribution weights of all power stations according to the percentage of installed capacity of all power stations in the total capacity, obtaining step total output per hour through the existing program, and inputting the distribution weights and the step total output per hour into the method; multiplying the input distribution weight of the output of each reservoir and the total output of the output of each hourly step to obtain an initial distribution result of the output of each power station in each time period, wherein in the embodiment, the number of the hydropower stations is 6, and the output upper limit and the output lower limit of each power station are shown in table 1:
TABLE 2 Upper and lower limits of output for each station (ten thousand kW)
| Power station numbering | Upper limit of output | Lower limit of output |
| 1 | 154.22 | 103.41 |
| 2 | 280.02 | 187.77 |
| 3 | 72.01 | 48.29 |
| 4 | 200.04 | 134.14 |
| 5 | 42.36 | 28.4 |
| 6 | 203.85 | 136.69 |
S2: finding out a time period in which the total output per hour is equal to the sum of the lower limits of the outputs of all the power stations (the length of each time period is 1 hour, 0-24 hours in a day are sequentially divided into 24 time periods by taking one hour as the time length, and the time periods are named as 1-24 respectively), wherein in the time periods (the time periods 1-7 and the time periods 22-24 in the embodiment), all the power stations are distributed according to the lowest output; the time periods during which the total output per hour equals the sum of the upper limits of the outputs of all the stations are found, and in these time periods (time periods 19-21 in this example), the stations are assigned the highest outputs. After the distribution, a secondary distribution result is obtained, and the secondary distribution result is shown in table 1:
TABLE 1 Secondary distribution results (WankW)
| Time interval sequence | 1 electric station | 2 power station | 3 power station | 4 power station | 5 electric station | 6 power station | Sum of step forces |
| 1 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 2 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 3 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 4 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 5 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 6 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 7 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 8 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 9 | 133.05 | 241.58 | 62.13 | 172.58 | 36.54 | 175.86 | 821.74 |
| 10 | 105.68 | 191.90 | 49.35 | 137.09 | 29.03 | 139.69 | 652.74 |
| 11 | 71.33 | 129.52 | 33.31 | 92.53 | 19.59 | 94.29 | 440.57 |
| 12 | 113.16 | 205.48 | 52.84 | 146.79 | 31.08 | 149.58 | 698.93 |
| 13 | 111.23 | 201.97 | 51.94 | 144.28 | 30.55 | 147.02 | 686.99 |
| 14 | 93.28 | 169.38 | 43.56 | 121.00 | 25.62 | 123.30 | 576.13 |
| 15 | 57.81 | 104.97 | 26.99 | 74.99 | 15.88 | 76.41 | 357.05 |
| 16 | 105.25 | 191.10 | 49.14 | 136.52 | 28.91 | 139.12 | 650.04 |
| 17 | 101.93 | 185.09 | 47.60 | 132.22 | 28.00 | 134.74 | 629.58 |
| 18 | 113.66 | 206.39 | 53.07 | 147.44 | 31.22 | 150.24 | 702.02 |
| 19 | 154.22 | 280.02 | 72.01 | 200.04 | 42.36 | 203.85 | 952.50 |
| 20 | 154.22 | 280.02 | 72.01 | 200.04 | 42.36 | 203.85 | 952.50 |
| 21 | 154.22 | 280.02 | 72.01 | 200.04 | 42.36 | 203.85 | 952.50 |
| 22 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 23 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 24 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| Mean output of the sun | 108.61 | 197.20 | 50.71 | 140.88 | 29.83 | 143.56 | 670.79 |
S3: calculating values of the output force of each power station in each time period in the secondary distribution result, which is lower than the lower output force limit or higher than the upper output force limit, and constructing a violation constraint matrix based on the calculation result, specifically:
s31: calculating the output value violating the upper output limit and the lower output limit of each power station in each time period:
in the formula:
N w，i，t the output value violating the upper limit and the lower limit of the output of the ith hydropower station in the tth period is i-1 … n, and T-1 … T, and the unit is ten thousand kW;
N i，t for the departure of the ith hydropower station in the t-th periodForce, in units of ten thousand kW; n is a radical of down，i，t The lower limit of the output of the ith hydropower station in the t period is ten thousand kW; n is a radical of up，i，t The output upper limit of the ith hydropower station in the t period is ten thousand kW.
S32: and forming a violation constraint matrix based on the output values violating the output upper limit and the output lower limit of each power station in each time period:
in the formula: n is a radical of w，all Is a violation of the constraint matrix. The values of the elements in the violation constraint matrix in this embodiment are shown in table 3:
table 3 violates the unit of value of each element in the constraint matrix: thousands kw
| Time interval sequence | Power station 1 | Electric station 2 | Power station 3 | Power station 4 | Power station 5 | Electric power station 6 |
| 1 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 2 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 3 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 4 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 5 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 6 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 7 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 8 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 9 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 10 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 11 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 12 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 13 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 14 | 0，00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 15 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 16 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 17 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 18 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 19 | 4.22 | 20.02 | 0.01 | 0.00 | 1.86 | 0.00 |
| 20 | 4.22 | 20.02 | 0.01 | 0.00 | 1.86 | 0.00 |
| 21 | 4.22 | 20.02 | 0.01 | 0.00 | 1.86 | 0.00 |
| 22 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 23 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| 24 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
S4: as shown in fig. 2 and 3, the power output of each power station in all time periods in the secondary distribution data set is distributed in the time period, so that all elements violating the constraint matrix in the distributed result are 0, and a tertiary distribution result is obtained, which specifically includes the following steps:
s41, obtaining the output distribution result of each power station in each time interval in the secondary distribution result as the vector to be operated, and calculating the element sum (N) of violation constraint matrix of all power stations in each time interval w，1，t +N w，2，t +…+N w，n，t ) (ii) a If the element sum is not equal to 0, go to step S42; if the element sum is equal to 0, go to step S44;
s42, adjusting the output value of each power station in each time interval to enable the output value of each power station in each time interval to meet the output constraint:
calculating the sum of the output values which are greater than the upper output limit:
in the formula: n is a radical of wk The sum of the output values is greater than the upper output limit; k is a power station serial number which is greater than or equal to the upper limit of the output, and K is 1, 2. Alpha is a serial number set of the power station which is greater than or equal to the upper limit of the output; calculating the sum of the output values smaller than the lower output limit:
in the formula: n is a radical of wq Is the sum of the output values less than the lower output limit; q is a power station serial number smaller than or equal to the lower limit of output, and Q is 1, 2. Beta is a serial number set of the power station which is less than or equal to the lower limit of the output;
calculating the sum of the remaining force values violating the force constraints:
N wl ＝N wk +N wq
in the formula: n is a radical of wl The sum of the remaining output values violating the output constraint;
adjusting the power plant output for this period:
s43, calculating element sum based on the adjusted power station output result, and if the value is not equal to 0, skipping to the step S32 again; if this value is equal to 0, go to step S34;
s44, if t < 24, if t is t +1, then proceed to step S31; and if t is 24, ending.
The results of the three assignments are detailed in table 4:
table 4 units of triplicate assignments: ten thousand kW
| Time interval sequence | Power station 1 | Electric station 2 | Power station 3 | Power station 4 | Power station 5 | Power station 6 | Sum of step forces |
| 1 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 2 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 3 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 4 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 5 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 6 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 7 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 8 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 9 | 133.05 | 241.58 | 62.13 | 172.58 | 36.54 | 175.86 | 821.74 |
| 10 | 105.68 | 191.90 | 49.35 | 137.09 | 29.03 | 139.69 | 652.74 |
| 11 | 71.33 | 129.52 | 33.31 | 92.53 | 19.59 | 94.29 | 440.57 |
| 12 | 113.16 | 205.48 | 52.84 | 146.79 | 31.08 | 149.58 | 698.93 |
| 13 | 111.23 | 201.97 | 51.94 | 144.28 | 30.55 | 147.02 | 686.99 |
| 14 | 93.28 | 169.38 | 43.56 | 121.00 | 25.62 | 123.30 | 576.13 |
| 15 | 57.81 | 104.97 | 26.99 | 74.99 | 15.88 | 76.41 | 357.05 |
| 16 | 105.25 | 191.10 | 49.14 | 136.52 | 28.91 | 139.12 | 650.04 |
| 17 | 101.93 | 185.09 | 47.60 | 132.22 | 28.00 | 134.74 | 629.58 |
| 18 | 113.66 | 206.39 | 53.07 | 147.44 | 31.22 | 150.24 | 702.02 |
| 19 | 150.00 | 260.00 | 72.00 | 210.00 | 40.50 | 220.00 | 952.50 |
| 20 | 150.00 | 260.00 | 72.00 | 210.00 | 40.50 | 220.00 | 952.50 |
| 21 | 150.00 | 260.00 | 72.00 | 210.00 | 40.50 | 220.00 | 952.50 |
| 22 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 23 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| 24 | 103.41 | 187.77 | 48.29 | 134.14 | 28.40 | 136.69 | 638.70 |
| Mean output of the sun | 108.08 | 194.70 | 50.71 | 142.12 | 29.60 | 145.58 | 670.79 |
S5, distributing the output of each power station in a time period to make the daily average output of each distributed power station the same as the daily average output in the initial distribution result, and specifically comprising the following steps:
s51, calculating the daily average output of each power station according to the initial distribution structure, and subtracting the daily average output of each power station of the result of S4 from the daily average output of each power station to obtain a difference vector:
[N d，1 ，N d，2 ，…，N d，n ]＝[N o，1 ，N o，2 ，…，N o，n ]-[N a，1 ，N a，2 ，…，N a，n ]
in the formula:
N d，i the difference value of the original daily average output of the ith hydropower station and the daily average output of the third distribution result is obtained; n is a radical of o，i The average output of the ith hydropower station in the original day; n is a radical of a，i The daily average output in the results of the three allocations of the ith hydropower station is obtained.
S52, finding out the time period in which the absolute values of the differences between the output of all the power stations and the upper limit and the lower limit of the output of each power station are greater than the preset buffer value from the three distribution results:
the output vector of the power station in the t-th time period is N 1，t ，N 2，t ，…，N n，t ]If the tth period is an allocation period, the following is satisfied:
in the formula: delta is a preset buffer value, and the unit is ten thousand kW;
s53, dividing the difference vector equally by the number of the distribution time interval, and adding the average vector into the power station output vector of the distribution time interval:
[N 1，z ，N 2，z ，…，N 3，z ]＝[N 1，z ，N 2，z ，…，N 3，z ]+[N d，1 ，N d，2 ，…，N d，n ]z∈γ
in the formula: [ N ] 1，z ，N 2，z ，…，N n，z ]The power station output vector of the distribution time interval is obtained, and z is the serial number of the distribution time interval; gamma is a sequence number set of the distribution time interval;
s54, calculating whether the output of the power station in the distributed time period exceeds the upper limit and the lower limit or not; if yes, reducing the buffer value, and turning to step S42; if not, ending the distribution to obtain a final distribution result, and the final distribution result is shown in table 5:
TABLE 5 Final assignment results
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions deviate from the technical solutions of the embodiments of the present invention.
Claims (8)
1. A load distribution method of a cascade power station based on an allowance method is characterized by comprising the following steps:
step S1: acquiring initial distribution results of output of each power station in each time period;
step S2: finding out time periods when the total output per hour is equal to the sum of the lower limits of the outputs of all the power stations in the initial distribution result, and distributing all the power stations according to the lowest output in the time periods; finding out time periods when the total output per hour is equal to the sum of the upper output limits of all the power stations, and distributing each power station according to the highest output in the time periods to obtain secondary distribution results of the output of each power station in each time period;
step S3: calculating the output values violating the upper output limit and the lower output limit of each power station in each time period, and constructing a violating constraint matrix;
step S4: distributing the output of each power station in all time intervals in the secondary distribution result within the time intervals, so that all elements violating the constraint matrix in the distributed result are 0, and obtaining a tertiary distribution result;
step S5: and distributing the output of each power station in each time period, so that the daily average output of each distributed power station is the same as the daily average output in the initial distribution result, and obtaining a final distribution result.
2. The cascade power station load distribution method based on the margin method as claimed in claim 1, wherein the initial calculation method of the hourly output distribution data of each power station in step S1 is as follows: and multiplying the distribution weight of the output of each reservoir input from the outside by the total output of the step per hour.
3. The method for load distribution of a stepladder power station based on the margin method as claimed in claim 1, wherein said step S3 specifically comprises:
step S31: calculating the output value violating the upper limit and the lower limit of the output of each power station in each time period, and specifically comprising the following steps:
in the formula:
N w，i，t the output value violating the upper limit and the lower limit of the output of the ith hydropower station in the tth period is i-1 … n, and T-1 … T, and the unit is ten thousand kW; n is a radical of hydrogen i，t The output of the ith hydropower station in the t time period is ten thousand kW; n is a radical of down，i，t The lower limit of the output of the ith hydropower station in the t period is ten thousand kW; n is a radical of up，i，t The output upper limit of the ith hydropower station in the t time period is ten thousand kW;
step S32: the violating constraint matrix is formed by the violating output values of the upper limit and the lower limit of the output of each power station in each time period, and specifically comprises the following steps:
in the formula: n is a radical of w，all Is a violation of the constraint matrix.
4. The cascade power station load distribution method based on the margin method as claimed in claim 3, wherein the step S4 comprises the following steps:
step S41: acquiring output distribution results of all power stations in each time period in the secondary distribution results as vectors to be operated, and calculating the sum of elements of violation constraint matrixes of all power stations in each time period; if the element sum is not equal to 0, go to step S42;
step S42: adjusting the output value of each power station in each time period to enable the output value of each power station in each time period to meet the output constraint;
step S43: calculating the sum of the elements based on the adjusted power plant output result, and if the value is not equal to 0, skipping to step S42 again; if this value is equal to 0, go to step S44;
step S44: if t < 24, t is t +1, and the process proceeds to step S41.
5. The method for load distribution to cascaded power stations based on the margin method as claimed in claim 4, wherein in step S44, if t is 24, it is ended.
6. The method of claim 4 wherein if the element sum is equal to 0 in step S41, the method proceeds to step S44.
7. The cascade power station load distribution method based on the margin method as claimed in claim 1, wherein the step S5 is specifically as follows:
step S51: calculating the daily average output of each power station according to the initial distribution structure, and subtracting the daily average output of each power station of the result of S4 from the daily average output of each power station to obtain a difference vector;
step S52: finding out a time period when the absolute values of the differences between the output of all the power stations and the respective output upper limit and lower limit are larger than a preset buffer value in the third distribution result;
step S53: the difference vector is divided equally according to the number of the distribution time intervals, and the divided vector is added into the power station output vector of the distribution time intervals;
step S54: calculating whether the output of the power station in the distribution time period after distribution exceeds the upper limit and the lower limit; if yes, reducing the buffer value, and turning to step S42; if not, the allocation is ended.
8. The method of claim 7 wherein the buffer value is 5-20 ten thousand kw.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210618163.0A CN114897419A (en) | 2022-06-01 | 2022-06-01 | Cascade power station load distribution method based on margin method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210618163.0A CN114897419A (en) | 2022-06-01 | 2022-06-01 | Cascade power station load distribution method based on margin method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114897419A true CN114897419A (en) | 2022-08-12 |
Family
ID=82726311
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210618163.0A Pending CN114897419A (en) | 2022-06-01 | 2022-06-01 | Cascade power station load distribution method based on margin method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114897419A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1255387A1 (en) * | 2001-05-04 | 2002-11-06 | Siemens Information and Communication Networks S.p.A. | Estimation of a plurality of channels in a diversity receiver |
| CN104636830A (en) * | 2015-02-12 | 2015-05-20 | 华中科技大学 | Water power and thermal power generation real-time load adjusting method of provincial grid under inflow change |
| CN106655280A (en) * | 2016-11-26 | 2017-05-10 | 大连理工大学 | Electric quantity control-based stair-step hydropower short-term peak-regulating model and solving method |
| CN108470249A (en) * | 2018-03-16 | 2018-08-31 | 大连理工大学 | A kind of Hydropower Stations short-term electricity generation dispatching method of coupling clustering and decision tree |
| CN110400232A (en) * | 2019-06-28 | 2019-11-01 | 云南电网有限责任公司 | Consider the monthly trading program quantity division method of GROUP OF HYDROPOWER STATIONS of power grid profile constraints |
| CN110472826A (en) * | 2019-07-09 | 2019-11-19 | 大连理工大学 | A kind of step power station load variations real-time adaptive method considering daily electricity deviation |
-
2022
- 2022-06-01 CN CN202210618163.0A patent/CN114897419A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1255387A1 (en) * | 2001-05-04 | 2002-11-06 | Siemens Information and Communication Networks S.p.A. | Estimation of a plurality of channels in a diversity receiver |
| CN104636830A (en) * | 2015-02-12 | 2015-05-20 | 华中科技大学 | Water power and thermal power generation real-time load adjusting method of provincial grid under inflow change |
| CN106655280A (en) * | 2016-11-26 | 2017-05-10 | 大连理工大学 | Electric quantity control-based stair-step hydropower short-term peak-regulating model and solving method |
| CN108470249A (en) * | 2018-03-16 | 2018-08-31 | 大连理工大学 | A kind of Hydropower Stations short-term electricity generation dispatching method of coupling clustering and decision tree |
| CN110400232A (en) * | 2019-06-28 | 2019-11-01 | 云南电网有限责任公司 | Consider the monthly trading program quantity division method of GROUP OF HYDROPOWER STATIONS of power grid profile constraints |
| CN110472826A (en) * | 2019-07-09 | 2019-11-19 | 大连理工大学 | A kind of step power station load variations real-time adaptive method considering daily electricity deviation |
Non-Patent Citations (2)
| Title |
|---|
| 吴志远;黄显峰;李昌平;刘志佳;颜山凯;: "基于分段粒子群算法的梯级水库多目标优化调度模型研究", 水资源与水工程学报, no. 03, 15 June 2020 (2020-06-15), pages 145 - 154 * |
| 杨莉: "梯级水电站联合发电优化调度方法研究与应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, 15 February 2012 (2012-02-15), pages 037 - 135 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106169108B (en) | Active power distribution network short-term active power optimization method containing battery energy storage system | |
| CN109492861B (en) | Method for decomposing medium-term electricity quantity trading plan of cascade hydropower station group | |
| CN109256799B (en) | New energy power system optimal scheduling method based on sample entropy | |
| CN110994694B (en) | Micro-grid source-charge-storage coordination optimization scheduling method considering differentiated demand response | |
| CN109980636B (en) | Wind, water and fire coordinated optimization scheduling method based on improved Benders decomposition method | |
| CN109840636B (en) | Newton method-based power system random rolling scheduling method | |
| CN112116150A (en) | Method for regulating heat accumulating type electric heating power market by load aggregators | |
| CN114421459A (en) | Cluster division evaluation method and system for large-scale grid connection of distributed power supply | |
| CN107679670A (en) | Sale of electricity optimization method commercially available from a kind of industrial park sale of electricity based on TC CVAR | |
| CN108365637A (en) | A kind of power transmission plan and draw water energy storing and electricity generating plan optimization method and system | |
| CN107341567B (en) | Storage capacity replacement calculation method for cascade reservoir group | |
| CN110766210B (en) | Short-term optimized scheduling method and system for cascade reservoir group | |
| CN114819406A (en) | Water-light complementary scheduling large system decomposition coordination optimization method | |
| CN116418001A (en) | Reservoir group multi-energy complementary scheduling method and system for coping with uncertainty of new energy | |
| CN112150190B (en) | Demand response complementary electricity price system and method for high-component flexible load | |
| CN114897419A (en) | Cascade power station load distribution method based on margin method | |
| CN112651177A (en) | Power distribution network flexible resource allocation method and system considering flexible service cost | |
| CN108539799B (en) | method and device for scheduling wind power in power grid | |
| CN116454944A (en) | Energy storage device optimal configuration method and system based on random production simulation | |
| CN105868841A (en) | Wind-water-fire combined dispatching method based on wind power priority grid connection | |
| CN110969294B (en) | Virtual power plant sectional output plan determination and energy storage system configuration method and device | |
| CN114475331A (en) | Electric automobile aggregation optimization scheduling method considering load rebound | |
| CN113506185A (en) | Power generation optimization scheduling method and device for cascade hydropower station and computer equipment | |
| CN105550944B (en) | Method for determining guaranteed output and annual average generated energy of perennial regulation reservoir | |
| CN111835005A (en) | Method and device for determining power grid transmission capacity of new energy base and related equipment |
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 |