CN106875049A - A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect - Google Patents
A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect Download PDFInfo
- Publication number
- CN106875049A CN106875049A CN201710071010.8A CN201710071010A CN106875049A CN 106875049 A CN106875049 A CN 106875049A CN 201710071010 A CN201710071010 A CN 201710071010A CN 106875049 A CN106875049 A CN 106875049A
- Authority
- CN
- China
- Prior art keywords
- unit
- dynamic economic
- max
- point effect
- periods
- 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 50
- 230000000694 effects Effects 0.000 title claims abstract description 41
- 230000003068 static effect Effects 0.000 claims abstract description 12
- 238000005457 optimization Methods 0.000 claims abstract description 7
- 230000033001 locomotion Effects 0.000 claims description 9
- 238000009987 spinning Methods 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 6
- 238000010977 unit operation Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 5
- 238000012804 iterative process Methods 0.000 claims description 5
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000005587 bubbling Effects 0.000 claims description 3
- 230000009194 climbing Effects 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000004870 electrical engineering 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
- 230000010349 pulsation Effects 0.000 description 1
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)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of meter and the Electrical Power System Dynamic economic load dispatching method of valve point effect, comprise the following steps:(1) the dynamic economic dispatch model of meter and valve point effect, is set up, the dynamic economic dispatch object function for considering valve point effect is given;(2) the Electrical Power System Dynamic economic load dispatching constraints of meter and valve point effect, is set up;(3), dynamic economic dispatch object function is simplified by the method for piece-wise linearization, the problem of MILP is converted into, then is solved with the MILP solvers in CPLEX, so as to obtain initial approximation solution;(4), by resetting exerting oneself bound and realizing the decoupling between adjacent moment for each period unit, dynamic economic dispatch problem is converted into static economy scheduling problem;(5) rapid solving, is carried out to the unit output of each period using improved dimension prompt drop method, final optimization solution is obtained.
Description
Technical field
The present invention relates to electrical engineering field, more particularly, to a kind of meter and valve point of power-system short-term traffic control
The dynamic economic dispatch method of effect.
Background technology
Electrical Power System Dynamic economic load dispatching is in the case where system operation constraints is met, by each generating set of Optimized Operation
Exert oneself, make system all operation periods overall running costs it is minimum, be one of key issue for facing in Operation of Electric Systems.
In systems in practice, the hot candied phenomenon occurred when being opened suddenly due to steam turbine intake valve can make unit consumption characteristic produce valve point
Effect, can be superimposed a pulsation effect on the consumption characteristic curve of unit, that is, produce valve point effect.Practical operating experiences table
Bright, ignoring valve point effect can be significantly affected the operation result of economic load dispatching.
Mathematically, it is considered to the Electrical Power System Dynamic economic load dispatching of valve point effect be a typical higher-dimension, dynamic it is non-
Convex, nonlinear optimal problem, solves more difficult.Analysis current research, mainly has two class methods to can be used to solve this problem.
The first kind is meta-heuristic optimized algorithm, and meta-heuristic algorithm and its hybrid algorithm are required for setting the parameter related to algorithm,
The setting of these parameters can produce influence to the performance of algorithm.Equations of The Second Kind is optimization algorithm, due to considering valve point effect
The non-linear and non-convex characteristic of dynamic economic dispatch problem, classical optimization algorithm is not appropriate for solving this problem.
The content of the invention
It is an object of the invention to provide a kind of meter and the Electrical Power System Dynamic economic load dispatching method of valve point effect, valve is being considered
On the basis of the static economy scheduling optimum result characteristic of point effect, by the way that constraint will be coupled between the dynamic economic dispatch period
(Climing constant) is decoupled, and introduces the concept of class singular point, finally using improved dimension prompt drop method (Modified
Dimensional Steepest Decline Method, MDSD) and MILP (Mixed Integer
Linear Programming, MILP) method solves to the problem.
To achieve the above object, the present invention uses following technical proposals:
A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect, comprises the following steps:
(1) the dynamic economic dispatch model of meter and valve point effect, is set up, the dynamic economic dispatch for considering valve point effect is given
Object function;
(2) power balance of, the Electrical Power System Dynamic economic load dispatching constraints of foundation meter and valve point effect, including system
Constraint and the constraint of system spinning reserve;
(3), dynamic economic dispatch object function is simplified by the method for piece-wise linearization, MIXED INTEGER is converted into linear
The problem of planning, then solved with the MILP solvers in CPLEX, so as to obtain initial approximation solution;
(4), on the basis of the initial approximation solution obtained in step (3), by resetting exerting oneself for each period unit
Bound realizes the decoupling between adjacent moment, and dynamic economic dispatch problem is converted into static economy scheduling problem;
(5), for the static economy scheduling problem after decoupling in step (4), using improved dimension prompt drop method to each
The unit output of period carries out rapid solving, obtains final optimization solution.
Further, in step (1) foundation meter and the dynamic economic dispatch model of valve point effect, it is considered to valve point effect
Dynamic economic dispatch object function is:
In formula:TC is total generation cost;CiIt is i-th group of cost function of t periods of generator;NGIt is generating set number;
NTIt is the period sum of dispatching cycle;ai, bi, ciIt is fuel cost coefficient;PI, tIt is to export i-th group of t period of generator
Active power;Pi minIt is the minimum load of generating set i;ei, fiIt is valve point effect coefficient.
Further, in step (2), the power balance of system is constrained to:
In formula:Pt lossIt is the system losses of t periods, is tried to achieve by B Y-factor method Ys;DtIt is the system loading of t periods, Bij
It is the coefficient of B matrixes,
The power output constraint of generator:
In formula, Pi maxIt is the EIAJ of i-th group of generator;
Generator Climing constant:
In formula, URi,DRiIt is the upper and lower climbing rate of fired power generating unit;
2) system spinning reserve constraint:
In formula, SRi,tIt is spinning reserve that the generating set i t periods provide;SSRt reqBe t period systems rotation it is standby
Use demand.
Further, in step (3), the dynamic economic dispatch object function of generating set i t periods is segmented linearisation
For:
Will generating set dynamic economic dispatch object function Ci(PI, t) it is expressed as being segmented the summation of line segment:Parameter KiRepresent
Generating set i from minimum unit output to the segments of maximum unit output, if each sinusoidal cycles is averagely divided into M
Section, then KiFor:
In formula, bM, iAnd cM, iThe slope and intercept of each segmentation line segment are represented respectively:
Wherein, the bound of each segmentationWithThe corresponding function representation of point on abscissa:
Setting parameter UM, i, tIt is 0,1 variable, represents in the t periods, whether generating set i's exerts oneself at m sections of linear segmented
On;The then power output P of unitI, tWith exert oneself P of the generating set i ts in the m sections of linear segmentedM, i, tShould meet with
Lower constraints:
3) each sinusoidal cycles is averaged and is divided into two sections (M=2), and the Optimized model for finally obtaining directly passes through MILP
Solved, obtained an initial approximation solution, convergence during solution can be set to the convergence gap of MILP less than one to
Fixed value.
Further, the bound equation of exerting oneself of unit is as follows in step (4):
In formula,WithRespectively new units limits of generating set i t periods;
It is determined that a certain unit of period t exert oneself new bound when, it is necessary to the period of using t-1 and period t+1 this
Unit is exerted oneself, and the unit output of t-1 periods needs the new power generating value that is obtained with dimension prompt drop method, and the unit of t+1 periods
Exert oneself, with the initial power generating value obtained in the first step.
Further, solving static economy scheduling problem using improved dimension prompt drop method can be divided into two steps:
1) first step needs to obtain unit output initial solution;
Improved dimension prompt drop method needs to use the singular point of unit, with P_SP 0I, jTo represent that i-th unit operation exists
On its j-th singular point:
P_SP is used hereinI, t, jRepresent that generating set i is the P_SP on its j-th singular point in exerting oneself for period tI, t, jBag
IncludeAnd in P_SP 0 between the twoI, j, NSp, tRepresent the sum of t period singular points;With DI, t, jCarry out table
Show the rate of descent of improved dimension prompt drop method, it is the t periods, and generating set i is by+1 singular point output drop of jth to
Every megawatt of cost when being exerted oneself at j singular point;Use LenghtI, t, jRepresent+1 singular point of jth of t periods, generating set i
Exert oneself and the difference exerted oneself at j-th singular point at place;
LengthI, t, j=P_SPI, t, j+1-P_SPI, t, j (20)
Initial solution need not strictly observe power-balance constraint, can be using bubbling method to DI, t, jIt is ascending to be ranked up,Represent the D after sequenceI, t, jOrder is order,AsSorted accordingly,
In first time iteration, every unit operation is exerted oneself on lower limit at it, per iteration once, according to the order of order in each unit
Exert oneself lower limit and on the basis of plus one sectionThe difference EER of power-balancetCalculated by following formula:
In iterative process, in EERtIn the case of zero, when it is closest to zero, iteration stopping;Then will be in iteration
During useThe exerting oneself on lower limit of every unit is added to respectively according to machine group # i, is obtained just
Every power generating value of unit in beginning solution, wherein parameter I0+1 are set to the value of final iterations;
2) second step is to adjust the initial solution that obtains in the first step, records the result after all adjustment, so as to remake into
One successive step:It is variableness to select I0 iteration and its state of first three vertical motion and the state of rear four vertical motions,
Further adjustment is done for each case, exactly selects wherein one unit as lax unit to meet power-balance constraint, it is right
All in the t periods meet the constraint of unit output boundSituation sum up the costs, choosing
Take the feasible solution that totle drilling cost can be made minimum.
The beneficial effects of the invention are as follows,
The present invention is adjusted on the basis of the static economy scheduling optimum result characteristic for considering valve point effect by by dynamic economy
Climing constant decoupling between spending the period, and the concept of class singular point is introduced, finally utilize improved dimension prompt drop method and mix
Integral linear programming method is solved to the problem, with obtain globally optimal solution in terms of feasibility and validity, preferably
Solve the problems, such as non-linear and non-convex characteristic dynamic economic dispatch.
Brief description of the drawings
Fig. 1 is flow chart of the present invention.
Specific embodiment
As shown in figure 1, a kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect, comprises the following steps:
(1) the dynamic economic dispatch model of meter and valve point effect, is set up, the dynamic economic dispatch for considering valve point effect is given
Object function, it is considered to which the dynamic economic dispatch object function of valve point effect is:
In formula:TC is total generation cost;CiIt is i-th group of cost function of t periods of generator;NGIt is generating set number;
NTIt is the period sum of dispatching cycle;ai, bi, ciIt is fuel cost coefficient;PI, tIt is to export i-th group of t period of generator
Active power;Pi minIt is the minimum load of generating set i;ei, fiIt is valve point effect coefficient.
(2) the Electrical Power System Dynamic economic load dispatching constraints of meter and valve point effect, the power balance constraint of system, are set up
For:
In formula:Pt lossIt is the system losses of t periods, is tried to achieve by B Y-factor method Ys;DtIt is the system loading of t periods, Bij
It is the coefficient of B matrixes,
The power output constraint of generator:
In formula, Pi maxIt is the EIAJ of i-th group of generator;
Generator Climing constant:
In formula, URi,DRiIt is the upper and lower climbing rate of fired power generating unit;
System spinning reserve is constrained:
In formula, SRi,tIt is spinning reserve that the generating set i t periods provide;SSRt reqBe t period systems rotation it is standby
Use demand.
(3), dynamic economic dispatch object function is simplified by the method for piece-wise linearization, MIXED INTEGER is converted into linear
The problem of planning, then solved with the MILP solvers in CPLEX, so as to obtain initial approximation solution.
The dynamic economic dispatch object function of generating set i t periods is segmented and linearly turns to:
Will generating set dynamic economic dispatch object function Ci(RI, t) it is expressed as being segmented the summation of line segment:Parameter KiRepresent
Generating set i from minimum unit output to the segments of maximum unit output, if each sinusoidal cycles is averagely divided into M
Section, then KiFor:
In formula, bM, iAnd cM, iThe slope and intercept of each segmentation line segment are represented respectively:
Wherein, the bound of each segmentationWithThe corresponding function representation of point on abscissa:
Setting parameter UM, i, tIt is 0,1 variable, represents in the t periods, whether generating set i's exerts oneself at m sections of linear segmented
On;PM, i, tExpression generating set i exerting oneself in the m sections of linear segmented in t, sometime, generates electricity to certain
Machine, in all of linear segmented, the generator output value in only one segmentation is greater than 0, in remaining all segmentation
Generator output is all 0, but the power generating value of the generator in which segmentation is determined by optimization more than 0.Then unit is defeated
Go out power PI, tWith exert oneself P of the generating set i ts in the m sections of linear segmentedM, i, tFollowing constraints should be met:
Each sinusoidal cycles is averaged and is divided into two sections (M=2), and the Optimized model for finally obtaining directly is entered by MILP
Row is solved, and obtains an initial approximation solution, and the convergence gap that convergence during solution can be set to MILP gives less than one
Value.
(4), on the basis of the initial approximation solution obtained in step (3), by resetting exerting oneself for each period unit
Bound realizes the decoupling between adjacent moment, and dynamic economic dispatch problem is converted into static economy scheduling problem;Step (4)
The bound equation of exerting oneself of middle unit is as follows:
In formula,WithRespectively new units limits of generating set i t periods;
It is determined that a certain unit of period t exert oneself new bound when, it is necessary to the period of using t-1 and period t+1 this
Unit is exerted oneself, and the unit output of t-1 periods needs the new power generating value that is obtained with dimension prompt drop method, and the unit of t+1 periods
Exert oneself, with the initial power generating value obtained in the first step.
(5), for the static economy scheduling problem after decoupling in step (4), using improved dimension prompt drop method to each
The unit output of period carries out rapid solving, obtains final optimization solution.Static warp is solved using improved dimension prompt drop method
Ji scheduling problem can be divided into two steps:
1) first step needs to obtain unit output initial solution;
Improved dimension prompt drop method needs to use the singular point of unit, with P_SP 0I, jTo represent that i-th unit operation exists
On its j-th singular point:
P_SP is used hereinI, t, jRepresent that generating set i is the P_SP on its j-th singular point in exerting oneself for period tI, t, jBag
IncludeAnd in P_SP 0 between the twoI, j, NSp, tRepresent the sum of t period singular points;With DI, t, jCarry out table
Show the rate of descent of improved dimension prompt drop method, it is the t periods, and generating set i is by+1 singular point output drop of jth to
Every megawatt of cost when being exerted oneself at j singular point;Use LengthI, t, jRepresent+1 singular point of jth of t periods, generating set i
Exert oneself and the difference exerted oneself at j-th singular point at place;
LengthI, t, j=P_SPI, t, j+1-P_SPI, t, j (20)
Initial solution need not strictly observe power-balance constraint, can be using bubbling method to DI, t, jIt is ascending to be ranked up,Represent the D after sequenceI, t, jOrder is order,AsSorted accordingly,
In first time iteration, every unit operation is exerted oneself on lower limit at it, per iteration once, according to the order of order in each unit
Exert oneself lower limit and on the basis of plus one sectionThe difference EER of power-balancetCalculated by following formula:
In iterative process, in EERtIn the case of zero, when it is closest to zero, iteration stopping;Then will be in iteration
During useThe exerting oneself on lower limit of every unit is added to respectively according to machine group # i, obtains initial
Every power generating value of unit in solution, wherein parameter I0+1 are set to the value of final iterations, obtain the value of parameter I0;
2) second step is to adjust the initial solution that obtains in the first step, records the result after all adjustment, so as to remake into
One successive step, makes it strictly observe power-balance constraint.Finally, the feasible solution that totle drilling cost can be made minimum is chosen.By the first small step
In iterative process be referred to as initial rise process, occur in that uphill process, then vertical motion be exactly open state, if I≤
I0, vertical motion is the state opened in the I times iteration during initial rise, if conversely, I > I0, are to close
State.On the basis of ensureing that every unit is all operated on singular point, adjust the I0 times during initial rise and change
The instead of state of the several uphill process after preceding several uphill process to the I0 times iteration.In general, state is changed
The quantity of vertical motion is more, and totle drilling cost increase is also more.Therefore, it is same by the quantity set of state change in uphill process
When change state no more than 4.Select I0 iteration and its state and the shape of rear four vertical motions of first three vertical motion
State is variableness, it will have 163 kinds of situations.For each case do further adjustment, exactly select wherein one unit by for
Lax unit meets power-balance constraint.All for the t periods meet the constraint of unit output boundSituation sum up the costs, selection can make the minimum feasible solution of totle drilling cost.
Although above-mentioned be described with reference to accompanying drawing to specific embodiment of the invention, not to present invention protection model
The limitation enclosed, one of ordinary skill in the art should be understood that on the basis of technical scheme those skilled in the art are not
Need the various modifications made by paying creative work or deformation still within protection scope of the present invention.
Claims (6)
1. it is a kind of count and valve point effect Electrical Power System Dynamic economic load dispatching method, it is characterised in that comprise the following steps:
(1) the dynamic economic dispatch model of meter and valve point effect, is set up, the dynamic economic dispatch target for considering valve point effect is given
Function;
(2) the Electrical Power System Dynamic economic load dispatching constraints of meter and valve point effect, is set up, including the power balance of system is constrained
And the constraint of system spinning reserve;
(3), dynamic economic dispatch object function is simplified by the method for piece-wise linearization, MILP is converted into
Problem, then solved with the MILP solvers in CPLEX, so as to obtain initial approximation solution;
(4), on the basis of the initial approximation solution obtained in step (3), by resetting exerting oneself up and down for each period unit
Limit to realize the decoupling between adjacent moment, dynamic economic dispatch problem is converted into static economy scheduling problem;
(5), for the static economy scheduling problem after decoupling in step (4), using improved dimension prompt drop method to each period
Unit output carry out rapid solving, obtain final optimization solution.
2. it is as claimed in claim 1 it is a kind of count and valve point effect Electrical Power System Dynamic economic load dispatching method, it is characterised in that
The dynamic economic dispatch model of foundation meter and valve point effect in step (1), it is considered to the dynamic economic dispatch target of valve point effect
Function is:
In formula:TC is total generation cost;CiIt is i-th group of cost function of t periods of generator;NGIt is generating set number;NTIt is
One period sum of dispatching cycle;ai, bi, ciIt is fuel cost coefficient;PI, tIt is having of exporting i-th group of t period of generator
Work(power;It is the minimum load of generating set i;ei, fiIt is valve point effect coefficient.
3. it is as claimed in claim 1 it is a kind of count and valve point effect Electrical Power System Dynamic economic load dispatching method, it is characterised in that
In step (2), the power balance of system is constrained to:
In formula:Pt lossIt is the system losses of t periods, is tried to achieve by B Y-factor method Ys;DtIt is the system loading of t periods, BijIt is B
The coefficient of matrix,
The power output constraint of generator:
In formula, Pi maxIt is the EIAJ of i-th group of generator;
Generator Climing constant:
In formula, URi,DRiIt is the upper and lower climbing rate of fired power generating unit;
2) system spinning reserve constraint:
In formula, SRi,tIt is spinning reserve that the generating set i t periods provide;SSRt reqThe spinning reserve for being t period systems need to
Ask.
4. it is as claimed in claim 1 it is a kind of count and valve point effect Electrical Power System Dynamic economic load dispatching method, it is characterised in that
In step (3), the dynamic economic dispatch object function of generating set i t periods is segmented and linearly turns to:
Will generating set dynamic economic dispatch object function Ci(PI, t) it is expressed as being segmented the summation of line segment:Parameter KiRepresent and generate electricity
Unit i from minimum unit output to the segments of maximum unit output, if each sinusoidal cycles is averagely divided into M sections,
Then KiFor:
In formula, bM, iAnd cM, iThe slope and intercept of each segmentation line segment are represented respectively:
Wherein, the bound of each segmentationWithThe corresponding function representation of point on abscissa:
Setting parameter UM, i, tIt is 0,1 variable, represents in the t periods, whether generating set i's exerts oneself in the m sections of linear segmented;Then
The power output P of unitI, tWith exert oneself P of the generating set i ts in the m sections of linear segmentedM, i, tShould meet it is following about
Beam condition:
3) each sinusoidal cycles is averaged and is divided into two sections (M=2), and the Optimized model for finally obtaining directly is carried out by MILP
Solve, obtain an initial approximation solution, the convergence gap that convergence during solution can be set to MILP is given less than one
Value.
5. it is as claimed in claim 1 it is a kind of count and valve point effect Electrical Power System Dynamic economic load dispatching method, it is characterised in that
The bound equation of exerting oneself of unit is as follows in step (4):
In formula,WithRespectively new units limits of generating set i t periods;
It is determined that a certain unit of period t exert oneself new bound when, it is necessary to the period of using t-1 and period t+1 this unit
Exert oneself, the unit output of t-1 periods needs the new power generating value obtained with dimension prompt drop method, and the unit output of t+1 periods
Then with the initial power generating value obtained in the first step.
6. it is as claimed in claim 5 it is a kind of count and valve point effect Electrical Power System Dynamic economic load dispatching method, it is characterised in that
Solving static economy scheduling problem using improved dimension prompt drop method can be divided into two steps:
1) first step needs to obtain unit output initial solution;
Improved dimension prompt drop method needs to use the singular point of unit, with P_SP0I, jTo represent i-th unit operation in its
On j singular point:
P_SP is used hereinI, t, jRepresent that generating set i is the P_SP on its j-th singular point in exerting oneself for period tI, t, jIncluding And in P_SP0 between the twoI, j, NSp, tRepresent the sum of t period singular points;With DI, t, jTo represent improvement
Dimension prompt drop method rate of descent, it is the t periods, and generating set i is strange to j-th by+1 singular point output drop of jth
Every megawatt of cost when being exerted oneself at dissimilarity;Use LengthI, t, jGo out at+1 singular point of jth for representing t periods, generating set i
The difference exerted oneself at power and j-th singular point;
LengthI, t, j=P_SPI, t, j+1-P_SPI, t, j (20)
Initial solution need not strictly observe power-balance constraint, can be using bubbling method to DI, t, jIt is ascending to be ranked up,
Represent the D after sequenceI, t, jOrder is order,AsSorted accordingly, for the first time
In iteration, every unit operation is exerted oneself on lower limit at it, per iteration once, according to order the exerting oneself in each unit of order
Lower limit adds one section with the basis ofThe difference EER of power-balancetCalculated by following formula:
In iterative process, in EERtIn the case of zero, when it is closest to zero, iteration stopping;Then will in an iterative process
UseThe exerting oneself on lower limit of every unit is added to respectively according to machine group # i, in obtaining initial solution
Every power generating value of unit, wherein parameter I0+1 are set to the value of final iterations;
2) second step is to adjust the initial solution obtained in the first step, records the result after all adjustment, so as to remake further
Adjustment:It is variableness to select I0 iteration and its state of first three vertical motion and the state of rear four vertical motions, for
Each case does further adjustment, exactly selects wherein one unit as lax unit to meet power-balance constraint, during for t
The all of section meet the constraint of unit output boundSituation sum up the costs, selection can make
The minimum feasible solution of totle drilling cost.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710071010.8A CN106875049A (en) | 2017-02-09 | 2017-02-09 | A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710071010.8A CN106875049A (en) | 2017-02-09 | 2017-02-09 | A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106875049A true CN106875049A (en) | 2017-06-20 |
Family
ID=59166532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710071010.8A Pending CN106875049A (en) | 2017-02-09 | 2017-02-09 | A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106875049A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108321792A (en) * | 2018-01-15 | 2018-07-24 | 国网山东省电力公司莱芜供电公司 | The more scene Cooperative Optimization Algorithms of Electrical Power System Dynamic economic load dispatching of meter and wind-powered electricity generation |
CN112396232A (en) * | 2020-11-19 | 2021-02-23 | 燕山大学 | Economic dispatching method and system for electric power system with valve point effect |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003134666A (en) * | 2001-10-16 | 2003-05-09 | Takeo Kawamura | Optimum power flow calculating system by nonlinear programming method |
CN105322571A (en) * | 2015-11-30 | 2016-02-10 | 华南理工大学 | Two-stage mixed integer programming based electric power system dispatching method |
CN106295884A (en) * | 2016-08-11 | 2017-01-04 | 中国电力科学研究院 | The Electrical Power System Dynamic economic load dispatching method of valve point effect is considered for rapid solving |
-
2017
- 2017-02-09 CN CN201710071010.8A patent/CN106875049A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003134666A (en) * | 2001-10-16 | 2003-05-09 | Takeo Kawamura | Optimum power flow calculating system by nonlinear programming method |
CN105322571A (en) * | 2015-11-30 | 2016-02-10 | 华南理工大学 | Two-stage mixed integer programming based electric power system dispatching method |
CN106295884A (en) * | 2016-08-11 | 2017-01-04 | 中国电力科学研究院 | The Electrical Power System Dynamic economic load dispatching method of valve point effect is considered for rapid solving |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108321792A (en) * | 2018-01-15 | 2018-07-24 | 国网山东省电力公司莱芜供电公司 | The more scene Cooperative Optimization Algorithms of Electrical Power System Dynamic economic load dispatching of meter and wind-powered electricity generation |
CN108321792B (en) * | 2018-01-15 | 2020-09-04 | 国网山东省电力公司莱芜供电公司 | Wind power-related dynamic economic dispatching multi-scene collaborative optimization algorithm for power system |
CN112396232A (en) * | 2020-11-19 | 2021-02-23 | 燕山大学 | Economic dispatching method and system for electric power system with valve point effect |
CN112396232B (en) * | 2020-11-19 | 2022-03-08 | 燕山大学 | Economic dispatching method and system for electric power system with valve point effect |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108879759B (en) | Harmonic analysis and treatment method for grid connection of double-fed wind generating set | |
CN107979111A (en) | A kind of energy management method for micro-grid based on the optimization of two benches robust | |
CN108846507A (en) | Electric-gas coupled system based on MIXED INTEGER Second-order cone programming economic load dispatching method a few days ago | |
CN105005940A (en) | Correlation-considered GEPOPF calculation method | |
CN105449713A (en) | Active power distribution network soft normally open point (SNOP) planning method considering distributed generator characteristics | |
CN104573875B (en) | A kind of method of the power generating facilities and power grids optimization planning of low-carbon | |
CN103425878B (en) | Power system Quasi dynamic trend and grid operation situation quick calculation method | |
CN104584357A (en) | Power distribution system loss reduction with distributed energy resource control | |
CN103904641B (en) | The micro-electrical network intelligent power generation of isolated island control method based on correlated equilibrium intensified learning | |
CN106505633A (en) | A kind of honourable access capacity determines method and device | |
CN107767086A (en) | New energy station output lower limit rolling amendment method based on generated power forecasting | |
CN104734175A (en) | Method for intelligently modifying wind speed and power curve of wind turbine generator | |
CN106779313A (en) | Multiple target distributed power source addressing constant volume method based on mixed integer programming | |
CN106875049A (en) | A kind of Electrical Power System Dynamic economic load dispatching method of meter and valve point effect | |
CN106712031A (en) | Sequence-robustness optimizing dispatching system and dispatching method of active power distribution network considering uncertainty | |
CN113988714A (en) | Dynamic planning method, device and medium for park integrated energy system considering multiple uncertainties | |
CN102790389A (en) | Method for calculating probability distribution of transient stability margin | |
CN115481856A (en) | Comprehensive energy system multi-scale scheduling method and system considering comprehensive demand response | |
Almasoudi et al. | Nonlinear coordination strategy between renewable energy sources and fuel cells for frequency regulation of hybrid power systems | |
CN107294101A (en) | A kind of multiple target unit built-up pattern and method for solving based on security domain target and constraint | |
CN106684931A (en) | Active control method and system of new-energy power plant | |
CN106208162A (en) | A kind of MINLP model Electrical Power System Dynamic dispatching method | |
CN117578601B (en) | Low-carbon power system source load interactive scheduling method and device | |
CN103559387B (en) | Fractal dimension steepest descent calculation method for valve point effect power grid economic dispatch problem | |
CN110048428B (en) | Power system probability load flow calculation method based on probability conservation principle |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170620 |