The hydroenergy storage station characteristic for making network system containing wind-powered electricity generation optimize operation is accurately simulated
Method
Technical field
The present invention relates to network optimization operation method, more particularly to it is a kind of make network system containing wind-powered electricity generation optimize run take out
The accurate analogy method of water storage station characteristic.
Background technology
Because the stability of wind power generation capacity is not high, thus the actual online of wind-powered electricity generation electricity with prediction electricity have it is larger partially
Difference, this accurately formulates to the annual power scheduling operating scheme containing wind-powered electricity generation power network and brings difficulty, is formed and preferably implements year
Spend the bottleneck of power scheduling scheme.
In the formulation of the electric power system dispatching operation year scheme containing wind-powered electricity generation, to improve power system comprehensive utilization of energy effect
Rate, solves Wind Power Development bottleneck problem, and recent domestic research institution has carried out multinomial related subject study.And utilize energy storage
The uncertain and anti-peaking problem that technology solves wind power output is the most intuitively scheme.In existing energy storage technology, draw water storage
Energy power station technology is the most ripe, using also the most extensive.Hydroenergy storage station is a kind of special power station, be both power supply and
Load, and with starting rapid, reliable flexible advantage, can effective tracking system load change, the peace of safeguards system
Full stable operation.So far, wind-powered electricity generation-water-storage cooperation scheme is the most ripe in lifting power grid wind receiving ability
Operating scheme.Therefore, how to set up suitable and close to operation of power networks reality hydroenergy storage station model is electric power
The focus and difficulties studied in system.Particularly large-scale wind power it is grid-connected after, hydroenergy storage station model was both accorded with
Wind-powered electricity generation management and running state is closed, can guarantee that model computational efficiency is problem in the urgent need to address again.
In existing model, document one《Wind-powered electricity generation annual plan formulating method based on time stimulatiom》(power system is automatic
Change the o. 11th page 13 of volume 38) a kind of wind-powered electricity generation annual plan method based on time stimulatiom is proposed, consider wind power output
The factors such as characteristic, part throttle characteristics, peak load regulation characteristic, power network submitting ability, power balance of the whole network containing wind-powered electricity generation is optimized by the period,
Establish the Optimized model for studying the wind-powered electricity generation plan of provincial power network year.But, do not accounted for drawing water in the Optimized model
The foundation of storage station's model, causes in model peak load regulation network ability and actual electric network deviation, and then influence final wind-powered electricity generation
The formulation of annual plan.Document two《Joint environmental and economic power dispatch
considering wind power integration:Empirical analysis from Liaoning Province
of China》(the 52nd phases of Renewable Energy page 260) using actual fired power generating unit, Hydropower Unit, wind power output with
And load data, construct the power system economic environment integrated distribution model containing wind power plant and hydroenergy storage station.But its is right
When hydroenergy storage station is modeled, storage capacity constraint, generated output constraint and the electric power constraint of drawing water of hydroenergy storage station are not considered.
Excessively coarse hydroenergy storage station modeling causes optimum results effective to the operation offer of actual hydroenergy storage station can not instruct.
Document three《Wind-powered electricity generation containing windfarm system matches capacity research with water-storage》(solar energy journal the 6th phase the 1037th of volume 33
Page), target is turned to social benefit maximum and sets up wind-powered electricity generation-water-storage cooperation Optimized model.Document four《Wind-powered electricity generation-draw water
Day operation Optimal Operation Model is combined in accumulation of energy》(the Automation of Electric Systems phase page 36 of volume 36 the 2nd) is with wind-powered electricity generation-water-storage
The maximizing the benefits of cooperation is target, it is considered to the start and stop limitation of pump-storage generator and generating-pumping operation mode conversion limit
System, characterizes Unit Commitment Constraint and hydroenergy storage station operating condition is constrained with quadratic constraints.Its mould to hydroenergy storage station
Type is set up, and does not meet actual wind-powered electricity generation dispatch situation.In real system scheduling, if system occurs abandoning wind, then now, take out
Water storage station can not carry out the generating that discharges water, and be not modeled for this constraint in existing scheduling model.
Therefore it provides a kind of mathematics of hydroenergy storage station for the electric power system dispatching operation year scheme containing wind-powered electricity generation is built
Mould method is very urgent.
The content of the invention
The present invention is to solve after large-scale wind power is grid-connected, how to make hydroenergy storage station model both meet wind-powered electricity generation and dispatch
Running status can guarantee that model computational efficiency again, and the formulation for the electric power system dispatching operation year scheme containing wind-powered electricity generation is asked
Topic, and a kind of accurate analogy method of hydroenergy storage station characteristic for making network system containing wind-powered electricity generation optimize operation for providing.
For achieving the above object, the accurate analogy method of hydroenergy storage station characteristic of the invention, including by the system
What system abandoned that wind state is modeled as being relevant to binary variable abandons wind state linear function, and based on the binary variable by system
The relevant operational state of all hydroenergy storage stations is modeled as being relevant to the power station running status of the binary variable
Linear function, linear function joint constrained system, simulates system hydroenergy storage station in the case where wind state is abandoned and exerts oneself shape described in two
State.
It is described abandon that the linear function table of wind state levies be:The wind-powered electricity generation theoretical prediction power of t period systemsConnect with the t periods
System receives wind powerDifference to be relevant to binary variable AtLinear function, according to system receive wind power it is big
Small judgement currently abandons wind state.
It is described abandon wind state linear function mathematic(al) representation be:
In formula:When system occurs abandoning wind, At=1;When system does not occur abandoning wind, AtWhen=0;It is non-negative continuous variable;It is constant;Eps is positive dimensionless.
The linear function of the power station running status at least includes output of power station state linear function, the output of power station state
Linear function characterize be:T period pump storage plant generator powerTo be relevant to binary variable AtLinear function,
When wind occurs abandoning in system, exerting oneself for hydroenergy storage station is zero, and when wind does not occur abandoning in system, hydroenergy storage station is exerted oneself
Between theoretical maximum generated output and theoretical minimum generated output.
The mathematic(al) representation of system's output of power station state linear function is:
In formula:It is nonnegative variable,WithIt is constant, hydroenergy storage station theoretical maximum is represented respectively
Generated output and minimum theoretical generated output.
The linear function of the running status in the power station also wraps the linear function of power station extraction water state, the extraction water state
Linear function characterize be:The t periods draw water state or the state of discharging water is followed successively by binary variableAny time,
Extraction or a kind of operating state for discharging water can only occur in the power station.
The mathematic(al) representation of the linear function of the power station t period extractions water state is:
In formula,Represent that power station does not carry out pumping operation,Represent that power station carries out pumping operation;Represent that power station does not carry out the operation that discharges water,Represent that power station carries out the operation that discharges water.
The linear function of the power station running status also includes that the linear function of output of power station constraint and power station are drawn water
The linear function of state;What the linear function of the output of power station fluctuation status was characterized is:The power station discharges water when generating electricity, and generates electricity
Power discharges water generated output in theoretical maximumDischarged water generated output with theoretical minimumBetween arbitrarily fluctuate, should
The mathematic(al) representation of linear function is:
The power station draw water state linear function characterize be:When the power station is drawn water, the unit that draws water once is opened, just
In full hair-like state, the mathematic(al) representation of the linear function is:
In formula,
It is nonnegative variable, represents that hydroenergy storage station is drawn water power;It is positive integer variable, represents that hydroenergy storage station is participated in excellent
The unit number of units of drawing water changed;It is separate unit water pumper kludge capacity;It is normal number, represents hydroenergy storage station
Draw water unit head station number.
The linear function of the power station running status also linear function including set optimization power constraint, minimum start and stop
The constraint of the linear function of time-constrain, the linear function of heat supply phase thermal power plant unit units limits and start and stop logic state it is linear
Function,
What the linear function of the set optimization power constraint was characterized is:Region jth platform conventional power unit power output's
Size is to be relevant to binary variableLinear function, the mathematic(al) representation of the linear function is:
In formula:Pj,max, Pj,minIt is constant, is expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit;It is j platforms
Unit the running status of period t binary variable,Represent that unit does not run,Represent that unit is transported
OK;
What the linear function of the minimum start and stop time-constrain was characterized is:Any time, unit j can only go out in period t
Now open machine operating stateOr shutdown action stateIn one kind;The mathematic(al) representation of the linear function is:
Or
In formula:It is binary variable,Represent that unit starts,Represent that unit is not starting shape
State,Represent that unit is shut down,Represent unit not in stopped status;konFor unit minimum opens the machine time;koffFor
Unit minimum downtime;Different types of Unit Commitment machine time parameter is different;
What the linear function of the heat supply phase thermal power plant unit units limits was characterized is:Back pressure type thermal power plant unit and bleeder are supplied
The coupled relation that electricity of the heat engine group within the heat supply phase is exerted oneself and heat is exerted oneself;The mathematic(al) representation of the linear function is:
In formula:For back pressure type thermal power plant unit electricity is exerted oneself size,For bleeder thermal power plant unit electricity is exerted oneself size;Respectively jth platform unit thermal power plant unit coupled thermomechanics coefficient,The coefficient of coup of lower limit of exerting oneself is represented,Represent
The coefficient of coup of the upper limit of exerting oneself;It is t period load of heat;Pj,max, Pj,minIt is constant, is expressed as jth platform unit
The upper limit of exerting oneself and lower limit of exerting oneself,It is pumping unit output size;
What the linear function of start and stop logic state constraint was characterized is:Operating states of the units, open machine state and operation
The logical constraint of state symbol;The mathematic(al) representation of the linear function is:
The running status in the power station also following linear function including following nonbinary variable modeling:Power station storage capacity is about
The linear function of beam, the linear function of power station water balance constraint, linear function, the region top-stitching of the constraint of unit climbing rate
The linear function of road transmission capacity constraint, linear function, the linear function of region account load balancing constraints, the wind of spinning reserve constraint
The linear function of electrical power constraint and the linear function of object function,
What the linear function of hydroenergy storage station storage capacity constraint was characterized is:Water size need to expire in hydroenergy storage station
The constraint of sufficient itself storage capacity size, when drawing water, water is no more than the maximum amount of water that storage capacity is allowed;When discharging water, water can not
Less than the least quantity that storage capacity is allowed;The mathematic(al) representation of the linear function is:
In formula, WmaxAnd WminIt is constant, represents the hydroenergy storage station maximum/minimum reservoir storage;WtIt is positive variable, represents
The water at the hydroenergy storage station current time;ηgAnd ηpConstant is, hydroenergy storage station water when expression discharges water and draws water respectively
Amount/electricity conversion coefficient.
What the linear function of hydroenergy storage station water balance constraint was characterized is:In hydroenergy storage station storage capacity and its
The corresponding relation of generated energy;The mathematic(al) representation of the linear function is:
W1=WInitialCap (14)
In formula,It is positive variable, represents that the hydroenergy storage station terminates the initial quantity of electricity of electricity, its value and subsequent timeIt is identical;W1It is first period storage capacity water size of hydroenergy storage station, its value should be the initial storage in power station,
WInitialCapIt is hydroenergy storage station initial storage water;
What the linear function of the linear function of unit climbing rate constraint was characterized is:Every unit unit interval can increase
Or reduce exert oneself;The mathematic(al) representation of the linear function is:
In formula:ΔPj,up、ΔPj,downCreep speed and lower creep speed respectively in the maximum of unit j;
What the linear function of the interregional line transmission capacity-constrained was characterized is:Circuit allows the work(of transmission between region
Rate is no more than its physical limit;The mathematic(al) representation of the linear function is:
It is the transmitted power of the transmission lines of t periods i-th;And Li,maxWith-Li,maxThe transmission of respectively the i-th transmission lines is held
Amount bound;Setting electric current reference direction is:Inflow region is positive direction, and outflow region is negative direction.SoCan take just
Negative value, it is positive and negative, represent the direction of power transmission
What the linear function of the linear function of spinning reserve constraint was characterized is:In order to realize the flat of system active power
Weighing apparatus, system should have certain spare capacity.The spare capacity of system refer in the case of system peak load, system can electricity consumption
Part of the source capacity more than generation load.The mathematic(al) representation of the linear function is:
In formula:The electric load of n region t periods is then represented, Pre is that positive rotation is standby;
What the linear function of the region account load balancing constraints was characterized is:Power system generated output and station service power load
Should Real-time Balancing, the mathematic(al) representation of the linear function is:
In formula:It is fired power generating unit gross capability,For the wind power that power network is received.
What the linear function of wind power constraint was characterized is:The wind power that power network is received should be less than its pre- measurement of power
Rate, the mathematic(al) representation of the linear function is:
In formula:Wind-powered electricity generation predicts watt level
What the object function was characterized is:Power network Preferred Acceptance wind-powered electricity generation electric power;The mathematic(al) representation of the function is:
Beneficial fruit of the invention is:1) present invention can forbid putting for hydroenergy storage station when wind occurs abandoning in system
Water generates electricity and carries out linear modelling so that the formulation model of the accurate power network for considering hydroenergy storage station characteristic year operational plan is still
It is mixed-integer programming model, it is ensured that the solution efficiency of model;2) method of the present invention more conforms to practical power systems tune
Degree situation and hydroenergy storage station actual operating mode, result of calculation can provide most intuitively basis for estimation for dispatcher.
Brief description of the drawings
Fig. 1 is that 3 subregion loads of certain regional water non-leap year are exerted oneself sequence.
Fig. 2 is the annual wind-powered electricity generation sequence in 3 subregions of certain regional water non-leap year.
Fig. 3 is that wind state and water-storage running state analysis figure are abandoned in certain region.
Specific embodiment
With reference to a specific embodiment, the present invention will be further described.
The first step, wind state is abandoned to system carries out linear modelling:
In formula, AtIt is binary variable, represents that t period systems abandon wind state variable, when its value is 1, represents that system occurs
Abandon wind;When its value is 0, expression system does not occur abandoning wind, i.e., wind-powered electricity generation is all received by system this moment.It is nonnegative variable, represents t
Period system receives wind power;It is constant, represents t period wind-powered electricity generation theoretical prediction watt levels;Eps is positive infinitely small
Amount.When wind occurs abandoning in system, i.e., theoretical wind power output receives wind power output more than systemDue to by formula
(1) constraint, now At=0;When wind does not occur abandoning in system, that is, wind power output and theory is received to exert oneself equal, due to formula (1)
Constraint, now At=1.Particularly, current time wind-powered electricity generation theory exert oneself for 0 when, i.e.,When, it is stipulated that At=0.Therefore,
The linear restriction can receive wind-powered electricity generation size according to the t periods, and simple and clear judges that system abandons wind state at current time.
Second step, should abandon wind state model, to the operation shape of all hydroenergy storage stations in system based on above-mentioned system
State is modeled.
1) hydroenergy storage station goes out force modeling when wind occurs abandoning in system
In formula,It is nonnegative variable, represents pump storage plant generator power;WithIt is constant, point
Biao Shi not hydroenergy storage station theoretical maximum generated output and minimum theoretical generated output.As shown in Equation 1, when system occurs abandoning wind
When, At=1, now hydroenergy storage station exert oneselfWhen there is not abandoning wind in system, At=1, now water-storage electricity
Exert oneself and run between its theoretical maximum minimum load in station.
When the joint of formula (1-2) constrains the system of causing and occurs abandoning wind, hydroenergy storage station forbids generatings that discharge water, the model and
Hydroenergy storage station actual motion state is consistent.Remaining constraints of hydroenergy storage station is consistent with conventional model, only does here
It is simple to introduce.
2) state and units limits are taken out/discharged water to hydroenergy storage station
In formula (3),It is binary variable, represents that state is taken out/discharged water to hydroenergy storage station respectively,It is 0
Represent that power station does not carry out pumping operation,For 1 expression power station carries out pumping operation;For 0 expression power station is not entered
The capable operation that discharges water,For 1 expression power station carries out pumping operation, at any time, hydroenergy storage station can only occur one
Plant operating state.
Formula (4) represents that its power can be in its theoretical maximum when hydroenergy storage station discharges water generatingAnd minimumArbitrarily fluctuated between generated output.
In formula (5),It is nonnegative variable, represents that hydroenergy storage station is drawn water power;It is positive integer variable, table
Show that hydroenergy storage station participates in the unit number of units of drawing water of optimization;It is separate unit water pumper kludge capacity;For normal
Number, represents that hydroenergy storage station is drawn water unit head station number.When formula (5) represents that hydroenergy storage station is drawn water, the unit that draws water once is opened
Open, just necessarily be in full hair-like state.
3) hydroenergy storage station storage capacity constraint
In formula (6), WmaxAnd WminIt is constant, represents the hydroenergy storage station maximum/minimum reservoir storage;WtIt is positive variable,
Represent the water at hydroenergy storage station current time;ηgAnd ηpConstant is, water-storage electricity when expression discharges water and draws water respectively
Stand water/electricity conversion coefficient.
4) hydroenergy storage station water balance constraint
W1=WInitialCap (8)
In formula (7),It is positive variable, represents that the hydroenergy storage station terminates electricity, its value is initial with subsequent time
Electricity is identical, and formula (8) is represented to hydroenergy storage station and assigns an initial reservoir capacity value.
Other constraints are for example:Set optimization power constraint, minimum start and stop time-constrain, heat supply phase thermal power plant unit are exerted oneself
Constraint, the constraint of start and stop logic state, the constraint of unit climbing rate, the constraint of interregional line transmission capacity-constrained, spinning reserve, area
Domain account load balancing constraints, wind power constraint and object function and prior art are essentially identical, specific as follows:
5) set optimization power constraint is used to characterize certain region jth platform conventional power unit t power outputSize,
To be positively correlated with binary variableLinear function, its mathematic(al) representation is:
In formula:Pj,max, Pj,minIt is constant, is expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit;It is j platforms
Unit the running status of period t binary variable,Represent that unit does not run,Represent that unit is transported
OK;
6) minimum start and stop time-constrain, is used to characterize any time, and unit j can only occur opening and motor-driven makees shape in period t
StateOr shutdown action stateIn one kind;The mathematic(al) representation of its linear function is
Or
In formula:It is binary variable,Represent that unit starts,Represent that unit is not starting shape
State,Represent that unit is shut down,Represent unit not in stopped status;konFor unit minimum opens the machine time;koffFor
Unit minimum downtime;Different types of Unit Commitment machine time parameter is different;
7) heat supply phase thermal power plant unit units limits, are used to characterize back pressure type thermal power plant unit and bleeder thermal power plant unit in heat supply
The coupled relation that electricity in phase is exerted oneself and heat is exerted oneself, the linear function of its constraint is characterized as:
In formula:For back pressure type thermal power plant unit electricity is exerted oneself size,For bleeder thermal power plant unit electricity is exerted oneself size;It is thermal power plant unit coupled thermomechanics coefficient,The coefficient of coup of lower limit of exerting oneself is represented,Represent the coupling of the upper limit of exerting oneself
Coefficient;It is t period load of heat;Pj,max, Pj,minIt is constant, is expressed as exerting oneself for jth platform unit and the upper limit and exerts oneself
Lower limit,It is pumping unit output size.
8) start and stop logic state constraint, is used to characterize operating states of the units, opens machine state and stopped status should meet
Logical constraint, is represented with following mathematic(al) representation:
9) unit climbing rate constraint, is used to characterize that every unit unit interval can increase or decrease exerts oneself, and its mathematics reaches
Formula is:
15th, Δ P in 16 formulaj,up、ΔPj,downCreep speed and lower creep speed respectively in the maximum of unit j;
10) interregional line transmission capacity-constrained, circuit allows the power for transmitting no more than between being used to characterize region
Its physical limit, its mathematics is up to formula:
It is the transmitted power of the transmission lines of t periods i-th, and Li,maxWith-Li,maxThe transmission of respectively the i-th transmission lines is held
Amount bound;Setting electric current reference direction is:Inflow region is positive direction, and outflow region is negative direction.SoCan take just
Negative value, it is positive and negative, represent the direction of power transmission.
11) spinning reserve constraint, is used to characterize the balance for realizing system active power, should there is certain spare capacity.System
The spare capacity of system refers to the part of the power available capacity more than generation load of system in the case of system peak load.Its
Mathematics is up to formula:
In formula:The electric load of n region t periods is then represented, Pre is that positive rotation is standby.
12) region account load balancing constraints, being used to characterize power system generated output should put down in real time with station service power load
Weighing apparatus.Its mathematics is up to formula:
In formula:It is fired power generating unit gross capability,For the wind power that power network is received.
13) wind power constraint, the wind power for being used to characterize power network receiving should be less than its theoretical power (horse-power).Its mathematics reaches
Formula is:
In formula:Wind-powered electricity generation theoretical power (horse-power) size.
14) object function, is used to characterize power network Preferred Acceptance wind-powered electricity generation electric power.Its mathematics is up to formula:
In formula:Then represent the electric load of n region t periods;Pre is that positive rotation is standby;T is total optimization time
Section;N is all numbers of partitions of power network.
One emulation embodiment of checking
Some region of network system operation conditions is simulated with above-mentioned the inventive method, by the network system in region
It is divided into subregion 1, subregion 2, the region of subregion 3 three.3 subregional horizontal year wind-powered electricity generation sequences, load are exerted oneself sequence
As shown in Figure 1 and Figure 2, simulation time step-length is 1 hour to row, and system adjusts unit distribution situation to be shown in Table 1-4.Biography of the region 1 to region 3
Defeated capacity limitation is 1800MW;Subregion 3 is 1500MW with the transmission capacity limits of subregion 2.The positive spare capacity of system is
660MW。
Gas formula unit tables of data is coagulated in the real system of table 1
Installed capacity/(MW) |
Minimum load/(MW) |
Number of units |
Region |
300 |
156 |
2 |
1 |
600 |
312 |
1 |
1 |
660 |
343 |
1 |
1 |
600 |
312 |
2 |
2 |
Bleeder unit tables of data in the real system of table 2
Installed capacity/(MW) |
Minimum load/(MW) |
Number of units |
Thermic load/(MW) |
CbValue |
CvValue |
Region |
100 |
52 |
10 |
100 |
0.72 |
0.95 |
1 |
200 |
102 |
13 |
200 |
0.56 |
0.95 |
1 |
300 |
156 |
18 |
300 |
0.65 |
0.98 |
1 |
330 |
155 |
2 |
330 |
0.60 |
0.93 |
1 |
660 |
243 |
1 |
660 |
0.60 |
0.96 |
1 |
200 |
102 |
2 |
200 |
0.56 |
0.92 |
2 |
200 |
102 |
1 |
200 |
0.56 |
0.95 |
3 |
600 |
282 |
1 |
600 |
0.52 |
0.95 |
3 |
Back pressure type unit tables of data in the real system of table 3
Unit capacity/(MW) |
Minimum load/(MW) |
Number of units |
Thermic load/(MW) |
CbValue |
Region |
64 |
38 |
1 |
350 |
0.16 |
1 |
Pump-storage generator tables of data in the real system of table 4
Different case optimum results are as shown in table 5.
Each case result of calculation comparative analysis table of table 5
Analyzed from table 5, when the operation year scheme of the electric power system dispatching containing wind-powered electricity generation is formulated, need to considered to water-storage
The modeling in power station, receives ability to improve the whole network wind-powered electricity generation obvious, and wind-powered electricity generation abandons wind rate and reduces 5.75%, receives wind-powered electricity generation more in the region
538639MW.Because this method is to carry out linear modelling to hydroenergy storage station, so calculating the time to management and running year scheme
Influence is little, increase only 0.3 minute.
As seen from Figure 3:1) when hydroenergy storage station is drawn water, as long as drawing water unit unlatching, the power that will overfill operation;2)
When wind occurs abandoning in system, hydroenergy storage station is forbidden discharging water the constraint of generating, and result of calculation more conforms to real system operation
Situation.