CN104993523A - Pumped storage power station characteristic accurate simulation method for optimized operation of wind power contained power grid system - Google Patents

Abstract

The invention relates to a pumped storage power station characteristic accurate simulation method for optimized operation of a wind power contained power grid system. The pumped storage power station characteristic accurate simulation method comprises the steps of: modeling a wind curtailment state of the wind power contained power grid system into a wind curtailment state linear function related to a binary variable; modeling relevant operating states of all pumped storage power stations in the system into a power station operating state linear function related to the binary variable on the basis of the binary variable; constraining the system through the combination of the two linear functions; and simulating output states of the pumped storage power stations when the system is in the wind curtailment state. The pumped storage power station characteristic accurate simulation method has the beneficial effects of conducting linear modeling on the situation that power generation by releasing water is forbidden by the pumped storage power stations when the system is in the wind curtailment state so that a model made by a power grid annual operating plan accurately considering pumped storage power station characteristics is still a mixed integer programming model, ensuring solving efficiency of the mode, being more inline with the actual power system dispatching condition and the actual operation conditions of the pumped storage power stations, and being capable of providing most intuitive judgment basis for dispatchers through calculation results.

Description

Make the hydroenergy storage station characteristic accurate analog method containing wind-powered electricity generation network system optimizing operation

Technical field

The present invention relates to the method that network optimization runs, particularly relate to a kind of hydroenergy storage station characteristic accurate analog method made containing wind-powered electricity generation network system optimizing operation.

Background technology

Because the stability of wind power generation capacity is not high, therefore the electricity of the actual online of wind-powered electricity generation has relatively large deviation with prediction electricity, this formulates to the annual power scheduling operating scheme containing wind-powered electricity generation electrical network comparatively accurately and brings difficulty, forms the bottleneck better implementing annual power scheduling scheme.

Run in the formulation of year scheme in the power system dispatching containing wind-powered electricity generation, for improving electric power system comprehensive utilization rate of energy source, solve Wind Power Development bottleneck problem, recent domestic research institution has carried out multinomial relevant subject study.And to utilize energy storage technology to solve the uncertainty of wind power output and anti-peaking problem be scheme the most intuitively.In existing energy storage technology, hydroenergy storage station technology is the most ripe, applies also extensive.Hydroenergy storage station is a kind of special hydroelectric station, is power supply and load, and has and start rapid, reliable advantage flexibly, can the change of effective tracking system load, the safe and stable operation of safeguards system.So far, wind-powered electricity generation-pumped storage cooperation scheme promotes operating scheme the most ripe in power grid wind receiving ability.Therefore, how setting up suitable and close to operation of power networks reality hydroenergy storage station model is the focus and difficulties studied in electric power system.Particularly large-scale wind power grid-connected after, how to make hydroenergy storage station model both meet wind-powered electricity generation management and running state, can ensure that again model computational efficiency is problem in the urgent need to address.

In existing model, document one " the wind-powered electricity generation annual plan formulating method based on time stimulatiom " (Automation of Electric Systems the 38th volume o. 11th the 13rd page) proposes a kind of wind-powered electricity generation annual plan method based on time stimulatiom, consider wind power output characteristic, part throttle characteristics, peak load regulation characteristic, electrical network send the factors such as ability, optimizing the power balance of the whole network containing wind-powered electricity generation by the period, establishing the Optimized model for studying the wind-powered electricity generation plan of provincial power network year.But, do not have in this Optimized model to consider the foundation to hydroenergy storage station model, cause peak load regulation network ability and actual electric network deviation to some extent in model, and then affect the formulation of final wind-powered electricity generation annual plan.Document two " Joint environmental andeconomic power dispatch considering wind power integration:Empirical analysis from Liaoning Province ofChina " (Renewable Energy the 52nd phase the 260th page) utilizes actual fired power generating unit, Hydropower Unit, wind power output and load data, constructs the electric power system economic environment integrated distribution model containing wind energy turbine set and hydroenergy storage station.But when it is to hydroenergy storage station modeling, do not consider the constraint of the capacity constrain of hydroenergy storage station, generated output and electric power constraint of drawing water.Too coarse hydroenergy storage station modeling causes optimum results can not run actual hydroenergy storage station provides effective guidance.Document three " mating capacity research with pumped storage containing windfarm system wind-powered electricity generation " (solar energy journal the 33rd volume the 6th phase the 1037th page), sets up wind-powered electricity generation-pumped storage cooperation Optimized model with the maximum target that turns to of social benefit.Document four " wind-powered electricity generation-pumped storage associating day operation Optimal Operation Model " (Automation of Electric Systems the 36th volume the 2nd phase the 36th page) with the maximizing the benefits of wind-powered electricity generation-pumped storage cooperation for target, consider start and stop restriction and the generating-pumping operation mode conversion restriction of pump-storage generator, characterize Unit Commitment Constraint and the constraint of hydroenergy storage station operating condition with quadratic constraints.It is set up the model of hydroenergy storage station, not realistic wind-powered electricity generation dispatch situation.In real system scheduling, if system occurs abandoning wind, so now, hydroenergy storage station can not carry out the generating that discharges water, and does not all carry out modeling for this constraint in existing scheduling model.

Therefore, a kind of Mathematical Modeling Methods for running the hydroenergy storage station of year scheme containing the power system dispatching of wind-powered electricity generation is provided to be very urgent.

Summary of the invention

The present invention is to solve large-scale wind power grid-connected after, hydroenergy storage station model how is made not only to meet wind-powered electricity generation management and running state but also model computational efficiency can be ensured, the problem of the formulation of year scheme is run for the power system dispatching containing wind-powered electricity generation, and a kind of hydroenergy storage station characteristic accurate analog method made containing wind-powered electricity generation network system optimizing operation provided.

For achieving the above object, hydroenergy storage station characteristic accurate analog method of the present invention, comprise and described system is abandoned wind state be modeled as wind state of the abandoning linear function being relevant to binary variable, and based on this binary variable, the relevant operational state of hydroenergy storage stations all in system is modeled as the linear function of the described power station running status being relevant to described binary variable, linear function associating constrained system described in two, the system that simulates to be exerted oneself state abandoning hydroenergy storage station under wind state.

The described linear function table of wind state of abandoning is levied: the wind-powered electricity generation theoretical prediction power of t period system wind power is received with t period welding system difference for being relevant to binary variable A ^tlinear function, receive the size of wind power to judge currently to abandon wind state according to system.

The described mathematic(al) representation abandoning wind state linear function is:

$(A^t-1)\×P_w^{t,n*}+eps\≤P_w^{t,n*}-P_w^t\≤A^t\×P_w^{t,n*}---\left(1\right)$

In formula: when system occurs abandoning wind, A ^t=1; When system does not occur abandoning wind, A ^twhen=0; for non-negative continuous variable; for constant; Eps is positive dimensionless.

The linear function of described power station running status at least comprises output of power station state linear function, and the linear function table of this output of power station state is levied: t period pump storage plant generator power for being relevant to binary variable A ^tlinear function, when system occurs abandoning wind, exerting oneself of hydroenergy storage station is zero, and when system does not occur abandoning wind, hydroenergy storage station is exerted oneself between theoretical maximum generation power and the minimum generated output of theory.

The mathematic(al) representation of described system output of power station state linear function is:

$P_{PSG\min }^t\×(1-A^t)\≤P_{PSG}^t\≤P_{PSGmax}^t\×(1-A^t)---\left(2\right)$

In formula: for nonnegative variable, with for constant, represent hydroenergy storage station theoretical maximum generated output and minimum theoretical generated output respectively.

The linear function of the running status in described power station also wraps the linear function of power station extraction water state, and the linear function of this extraction water state characterizes: draw water state or the state of discharging water of t period is followed successively by binary variable any time, a kind of operate condition that described power station can only occur extraction or discharge water.

The mathematic(al) representation of the linear function of described power station t period extraction water state is:

$0\≤{\mathrm{PS}}_a^t+{\mathrm{PS}}_b^t\≤1---\left(3\right)$

In formula, represent that pumping operation is not carried out in power station, represent that pumping operation is carried out in power station; represent that the operation that discharges water is not carried out in power station, represent that the operation that discharges water is carried out in power station.

The linear function of described power station running status also comprises the linear function of described output of power station constraint and power station and to draw water the linear function of state; The linear function of described output of power station fluctuation status characterizes: when described power station discharges water and generates electricity, generated output is at the maximum generated output that discharges water of theory discharge water generated output minimum with theory between fluctuate arbitrarily, the mathematic(al) representation of this linear function is:

${\mathrm{PS}}_b^t\×P_{PSGmin}^t\≤P_{PSG}^t\≤{\mathrm{PS}}_b^t\×P_{PSGmax}^t---\left(4\right)$

Described power station draw water state linear function characterize: time drawing water in described power station, draw water unit once open, be just in full hair-like state, the mathematic(al) representation of this linear function is:

$\left\{\begin{array}{c}{P}_{PSP}^t={\mathrm{PS}}_a^t\×N_{PS}^t\×P_{PSPunit}^t\\ 0\≤N_{PS}^t\≤N_{PSmax}^t\end{array}\right.---\left(5\right),$ In formula, for nonnegative variable, represent that hydroenergy storage station is drawn water power; for positive integer variable, represent that the water pumper that hydroenergy storage station participates in optimizing is organized a performance number; for separate unit water pumper kludge capacity; for normal number, represent that hydroenergy storage station is drawn water unit head station number.

The linear function of described power station running status also comprises the linear function of the linear function of set optimization power constraint, the linear function of minimum start and stop time-constrain, the linear function of heat supply phase thermal power plant unit units limits and start and stop logic state constraint,

The linear function of described set optimization power constraint characterizes: region jth platform conventional power unit power output size for being relevant to binary variable linear function, the mathematic(al) representation of this linear function is:

$P_{j,\min }\·X_j^t\≤P_j^t\≤P_{j,max}\·X_j^t---\left(6\right)$

In formula: P _{j, max}, P _{j, min}for constant, be expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit; the binary variable of j platform unit in the running status of period t, represent that unit does not run, represent that unit runs;

The linear function of described minimum start and stop time-constrain characterizes: any time, unit j can only occur opening machine operate condition at period t or shutdown action state in one; The mathematic(al) representation of this linear function is:

$Y_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1---\left(7\right)$

Or $Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1---\left(8\right)$

In formula: for binary variable, represent that unit starts, represent that unit is not at starting state, represent that unit is shut down, represent that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; Dissimilar Unit Commitment machine time parameter is different;

The linear function of described heat supply phase thermal power plant unit units limits characterizes: the coupled relation that back pressure type thermal power plant unit and the electricity of bleeder thermal power plant unit within the heat supply phase are exerted oneself and heat is exerted oneself; The mathematic(al) representation of this linear function is:

$P_{j,BY}^t=C_j^b\·H_j^t---\left(9\right)$

$H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,\max }-H_j^t\·C_j^v---\left(10\right)$

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; be respectively jth platform unit thermal power plant unit coupled thermomechanics coefficient, represent the coupling coefficient of lower limit of exerting oneself, represent the coupling coefficient of the upper limit of exerting oneself; for t period load of heat; P _{j, max}, P _{j, min}for constant, be expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit, for unit output size of bleeding;

The linear function of described start and stop logic state constraint characterizes: operating states of the units, open machine state and running status meets logical constraint; The mathematic(al) representation of this linear function is:

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(11\right)$

The running status in described power station also comprises the following linear function of following nonbinary variable modeling: the linear function of the linear function of the linear function that the linear function of power station capacity constrain, described power station water balance retrain, unit climbing rate constraint, the linear function of interregional line transmission capacity-constrained, spinning reserve constraint, the linear function of region account load balancing constraints, the linear function of wind power constraint and the linear function of target function

The linear function of described hydroenergy storage station capacity constrain characterizes: in hydroenergy storage station, water yield size need meet the constraint of self storage capacity size, and when drawing water, the water yield can not exceed the maximum amount of water that storage capacity allows; When discharging water, the water yield can not be less than the least quantity that storage capacity allows; The mathematic(al) representation of this linear function is:

$W^t-W_{max}\≤P_{PSG}^t\·{\mathrm{\η}}_g-P_{PSP}^t\·{\mathrm{\η}}_p\≤W^t-W_{min}---\left(12\right)$

In formula, W _maxand W _minfor constant, represent maximum/minimum reservoir storage of this hydroenergy storage station; W ^tfor positive variable, represent the water yield of this hydroenergy storage station current time; η _gand η _pbe constant, the hydroenergy storage station water yield/electricity conversion coefficient when expression discharges water and draws water respectively.

The linear function of described hydroenergy storage station water balance constraint characterizes: the corresponding relation of storage capacity and its energy output in hydroenergy storage station; The mathematic(al) representation of this linear function is:

$\left\{\begin{array}{c}{P}_{PSG}^t\·{\mathrm{\η}}_g-P_{PSP}^t\·{\mathrm{\η}}_p=W^t-W_{end}^t\\ W^t=W_{end}^{t-1}\end{array}\right.---\left(13\right)$

W ¹＝W _InitialCap(14)

In formula, for positive variable, represent that this hydroenergy storage station stops electricity, the initial quantity of electricity of its value and subsequent time identical; W ¹for hydroenergy storage station first period storage capacity water yield size, its value should be the initial storage in power station, W _initialCapfor the hydroenergy storage station initial storage water yield;

The linear function of the linear function of described unit climbing rate constraint characterizes: what every platform unit unit interval can increase or reduce exerts oneself; The mathematic(al) representation of this linear function is:

$P_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}---\left(15\right)$

$P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}---\left(16\right)$

In formula: Δ P _{j, up}, Δ P _{j, down}be respectively the maximum upper creep speed of unit j and lower creep speed;

The linear function of described interregional line transmission capacity-constrained characterizes: between region, circuit allows the power of transmission can not exceed its physical restriction; The mathematic(al) representation of this linear function is:

$-L_{i,max}\≤L_i^t\≤L_{i,max}---\left(17\right)$

for the transmitted power of t period i-th transmission lines; And L _{i, max}with-L _{i, max}be respectively the i-th transmission lines transmission capacity bound; Setting current reference direction is: inflow region is positive direction, and outflow region is negative direction.So can positive and negative values be got, positive and negative, represent the direction of power delivery

The linear function of the linear function of described spinning reserve constraint characterizes: in order to realize the balance of system active power, system should have certain reserve capacity.The reserve capacity of system refers to that, in system peak load situation, the power available capacity of system is greater than the part of generation load.The mathematic(al) representation of this linear function is:

$-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,\max }\·X_j^t\≤-\underset{n=1}{\overset{N}{\Σ}}{P}_l^{t,n}-\mathrm{Pr}e---\left(18\right)$

In formula: then represent the electric load of n region t period, Pre is that positive rotation is for subsequent use;

The linear function of described region account load balancing constraints characterizes: electric power system generated output and station service power load should Real-time Balancings, and the mathematic(al) representation of this linear function is:

$P_{all}^{t,n}+P_w^{t,n}+P_{PSG}^t-P_{PSP}^t+L_i^t=P_l^{t,n}---\left(19\right)$

In formula: for fired power generating unit gross capability, for the wind power that electrical network is received.

The linear function of described wind power constraint characterizes: the wind power that electrical network is received should be less than its predicted power, and the mathematic(al) representation of this linear function is:

$0\≤P_w^{t,n}\≤P_w^{t,n^{*}}---\left(20\right)$

In formula: wind-powered electricity generation predicted power size

Described target function characterizes: electrical network Preferred Acceptance wind-powered electricity generation electric power; The mathematic(al) representation of this function is:

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_w^{t,n}---\left(21\right).$

Useful fruit of the present invention is: 1) the present invention can carry out linear modelling for hydroenergy storage station generating of forbidding when system occurs abandoning wind discharging water, the formulation model making the electrical network year operational plan accurately considering hydroenergy storage station characteristic is still mixed-integer programming model, ensure that the solution efficiency of model; 2) the more realistic power system dispatching situation of method of the present invention and hydroenergy storage station actual operating mode, result of calculation can provide basis for estimation the most intuitively for dispatcher.

Accompanying drawing explanation

Fig. 1 is that certain regional water 3 subregion load in non-leap year is exerted oneself sequence.

Fig. 2 is certain annual wind-powered electricity generation sequence in regional water 3 subregions in non-leap year.

Fig. 3 is that wind state and pumped storage running state analysis figure are abandoned in certain region.

Embodiment

Below in conjunction with an embodiment, the present invention will be further described.

The first step, wind state is abandoned to system and carries out linear modelling:

$(A^t-1)\×P_w^{t,n*}+eps\≤P_w^{t,n*}-P_w^t\≤A^t\×P_w^{t,n*}---\left(1\right)$

In formula, A ^tfor binary variable, represent that t period system abandons wind state variable, when its value is 1, expression system occurs abandoning wind; When its value is 0, expression system does not occur abandoning wind, and namely wind-powered electricity generation is all received by system this moment. for nonnegative variable, represent that t period system receives wind power; for constant, represent t period wind-powered electricity generation theoretical prediction watt level; Eps is positive dimensionless.When system occurs abandoning wind, namely theoretical wind power output is greater than system receiving wind power output owing to being subject to the constraint of formula (1), now A ^t=0; When system does not occur abandoning wind, wind power output and theory is namely received to exert oneself equal, due to the constraint of formula (1), now A ^t=1.Especially, current time wind-powered electricity generation theory is exerted oneself when being 0, namely time, regulation A ^t=0.Therefore, this linear restriction can receive wind-powered electricity generation size according to the t period, and simple and clear system of judging abandons wind state at current time.

Second step, should abandon wind state model, carry out modeling to the running status of hydroenergy storage stations all in system based on above-mentioned system.

1) when system occurs abandoning wind, hydroenergy storage station goes out force modeling

$P_{PSG\min }^t\×(1-A^t)\≤P_{PSG}^t\≤P_{PSGmax}^t\×(1-A^t)---\left(2\right)$

In formula, for nonnegative variable, represent pump storage plant generator power; with for constant, represent hydroenergy storage station theoretical maximum generated output and minimum theoretical generated output respectively.As shown in Equation 1, when system occurs abandoning wind, A ^t=1, now the exerting oneself of hydroenergy storage station when system does not occur abandoning wind, A ^t=1, now hydroenergy storage station exert oneself between its theoretical minimax is exerted oneself run.

Formula (1-2) combine constraint when making system occur abandoning wind, hydroenergy storage station forbids the generating that discharges water, and this model conforms to hydroenergy storage station actual motion state.All the other constraintss of hydroenergy storage station are consistent with conventional model, only do simple introduction here.

2) state and units limits are taken out/discharged water to hydroenergy storage station

$0\≤{\mathrm{PS}}_a^t+{\mathrm{PS}}_b^t\≤1---\left(3\right)$

${\mathrm{PS}}_b^t\×P_{PSG\min }^t\≤P_{PSG}^t\≤{\mathrm{PS}}_b^t\×P_{PSGmax}^t---\left(4\right)$

$\left\{\begin{array}{c}{P}_{PSP}^t={\mathrm{PS}}_a^t\×N_{PS}^t\×P_{PSPunit}^t\\ 0\≤N_{PS}^t\≤N_{PS\max }^t\end{array}\right.---\left(5\right)$

In formula (3), for binary variable, represent that state is taken out/discharged water to hydroenergy storage station respectively, be that pumping operation is not carried out in 0 expression power station, be that pumping operation is carried out in 1 expression power station; be that the operation that discharges water is not carried out in 0 expression power station, be that pumping operation is carried out in 1 expression power station, at any time, can only there is a kind of operate condition in hydroenergy storage station.

Formula (4) represents that its power can be maximum in its theory when hydroenergy storage station discharges water generating with minimum fluctuate arbitrarily between generated output.

In formula (5), for nonnegative variable, represent that hydroenergy storage station is drawn water power; for positive integer variable, represent that the water pumper that hydroenergy storage station participates in optimizing is organized a performance number; for separate unit water pumper kludge capacity; for normal number, represent 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 open, just must be in full hair-like state.

3) hydroenergy storage station capacity constrain

$W^t-W_{max}\≤P_{PSG}^t\·{\mathrm{\η}}_g-P_{PSP}^t\·{\mathrm{\η}}_p\≤W^t-W_{\min }---\left(6\right)$

In formula (6), W _maxand W _minfor constant, represent maximum/minimum reservoir storage of this hydroenergy storage station; W ^tfor positive variable, represent the water yield of this hydroenergy storage station current time; η _gand η _pbe constant, the hydroenergy storage station water yield/electricity conversion coefficient when expression discharges water and draws water respectively.

4) hydroenergy storage station water balance constraint

$\{\begin{array}{c}{P}_{PSG}^t\·{\mathrm{\η}}_g-P_{PSP}^t\·{\mathrm{\η}}_p=W^t-W_{end}^t\\ W^t=W_{end}^{t-1}\end{array}---\left(7\right)$

W ¹＝W _InitialCap(8)

In formula (7), for positive variable, represent that this hydroenergy storage station stops electricity, its value is identical with the initial quantity of electricity of subsequent time, and formula (8) expression composes an initial reservoir capacity value to hydroenergy storage station.

Other constraint is such as: the constraint of set optimization power constraint, minimum start and stop time-constrain, heat supply phase thermal power plant unit units limits, start and stop logic state, the constraint of unit climbing rate, the constraint of interregional line transmission capacity-constrained, spinning reserve, the constraint of region account load balancing constraints, wind power are substantially identical with prior art with target function, specific as follows:

5) set optimization power constraint is in order to characterize certain region jth platform conventional power unit t power output size, for being positively correlated with binary variable linear function, its mathematic(al) representation is:

$P_{j,\min }\·X_j^t\≤P_j^t\≤P_{j,max}\·X_j^t---\left(9\right)$

In formula: P _{j, max}, P _{j, min}for constant, be expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit; the binary variable of j platform unit in the running status of period t, represent that unit does not run, represent that unit runs;

6) minimum start and stop time-constrain, in order to characterize any time, unit j can only occur opening machine operate condition at period t or shutdown action state in one; The mathematic(al) representation of its linear function is

$Y_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1---\left(10\right)$

Or $Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1---\left(11\right)$

In formula: for binary variable, represent that unit starts, represent that unit is not at starting state, represent that unit is shut down, represent that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; Dissimilar Unit Commitment machine time parameter is different;

7) heat supply phase thermal power plant unit units limits, in order to characterize the coupled relation that back pressure type thermal power plant unit and the electricity of bleeder thermal power plant unit within the heat supply phase are exerted oneself and heat is exerted oneself, linear function of its constraint is characterized by:

$P_{j,BY}^t=C_j^b\·H_j^t---\left(12\right)$

$H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,max}-H_j^t\·C_j^v---\left(13\right)$

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; for thermal power plant unit coupled thermomechanics coefficient, represent the coupling coefficient of lower limit of exerting oneself, represent the coupling coefficient of the upper limit of exerting oneself; for t period load of heat; P _{j, max}, P _{j, min}for constant, be expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit, for unit output size of bleeding.

8) start and stop logic state constraint, in order to characterizing operating states of the units, open the logical constraint that machine state and stopped status should meet, represent with following mathematic(al) representation:

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(14\right)$

9) unit climbing rate constraint, in order to exerting oneself of characterizing that every platform unit unit interval can increase or reduce, its mathematics reaches formula and is:

$P_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}---\left(15\right)$

$P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}---\left(16\right)$

15, Δ P in 16 formula _{j, up}, Δ P _{j, down}be respectively the maximum upper creep speed of unit j and lower creep speed;

10) interregional line transmission capacity-constrained, allow the power of transmission can not exceed its physical restriction in order to characterize circuit between region, its mathematics reaches formula and is:

$-L_{i,max}\≤L_i^t\≤L_{i,max}---\left(17\right)$

for the transmitted power of t period i-th transmission lines, and L _{i, max}with-L _{i, max}be respectively the i-th transmission lines transmission capacity bound; Setting current reference direction is: inflow region is positive direction, and outflow region is negative direction.So can positive and negative values be got, positive and negative, represent the direction of power delivery.

11) spinning reserve constraint, in order to characterize the balance realizing system active power, should have certain reserve capacity.The reserve capacity of system refers to that, in system peak load situation, the power available capacity of system is greater than the part of generation load.Its mathematics reaches formula:

$-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,\max }\·X_j^t\≤-\underset{n=1}{\overset{N}{\Σ}}{P}_l^{t,n}-\mathrm{Pr}e---\left(18\right)$

In formula: then represent the electric load of n region t period, Pre is that positive rotation is for subsequent use.

12) region account load balancing constraints, in order to characterize electric power system generated output and station service power load should Real-time Balancing.Its mathematics reaches formula:

$P_{all}^{t,n}+P_w^{t,n}+P_{PSG}^t-P_{PSP}^t+L_i^t=P_l^{t,n}---\left(19\right)$

In formula: for fired power generating unit gross capability, for the wind power that electrical network is received.

13) wind power constraint, the wind power received in order to characterize electrical network should be less than its theoretical power (horse-power).Its mathematics reaches formula:

$0\≤P_w^{t,n}\≤P_w^{t,n^{*}}---\left(20\right)$

In formula: wind-powered electricity generation theoretical power (horse-power) size.

14) target function, in order to characterize electrical network Preferred Acceptance wind-powered electricity generation electric power.Its mathematics reaches formula:

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_w^{t,n}---\left(21\right)$

In formula: then represent the electric load of n region t period; Pre is that positive rotation is for subsequent use; Discontinuity surface when T is total optimization; N is all numbers of partitions of electrical network.

The emulation embodiment of a checking

Simulate by the network system operation conditions of above-mentioned the inventive method to a certain region, the network system in region is divided into subregion 1, subregion 2, region, 3 three, subregion.3 the annual wind-powered electricity generation sequences of subregional forcasted years, loads exert oneself sequence as shown in Figure 1 and Figure 2, and simulation time step-length is 1 hour, and tracking unit distribution situation is in Table 1-4.Region 1 is 1800MW to the transmission capacity limits in region 3; Subregion 3 is 1500MW with the transmission capacity limits of subregion 2.The positive reserve capacity of system is 660MW.

Solidifying gas formula unit tables of data in table 1 real system

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 table 2 real system

Installed capacity/(MW)	Minimum load/(MW)	Number of units	Heat load/(MW)	C bValue	C vValue	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 table 3 real system

Unit capacity/(MW)	Minimum load/(MW)	Number of units	Heat load/(MW)	C bValue	Region
						64	38	1	350	0.16	1

Pump-storage generator tables of data in table 4 real system

Different case optimum results is as shown in table 5.

Table 5 each case result of calculation comparative analysis table

Analyzed from table 5, when formulating the power system dispatching operation year scheme containing wind-powered electricity generation, need consider the modeling to hydroenergy storage station, improve obviously the whole network wind-powered electricity generation ability of receiving, wind-powered electricity generation is abandoned wind rate and is reduced 5.75%, region many receivings wind-powered electricity generation 538639MW.Because this method carries out linear modelling to hydroenergy storage station, so little on impact computing time of management and running year scheme, increase only 0.3 minute.

As seen from Figure 3: when 1) hydroenergy storage station is drawn water, as long as the unit that draws water is opened, the power that will overfill is run; 2) when system occurs abandoning wind, hydroenergy storage station forbids the constraint generated electricity that discharges water, the more realistic running situation of result of calculation.

Claims (10)

1. make the hydroenergy storage station characteristic accurate analog method containing wind-powered electricity generation network system optimizing operation, it is characterized in that: comprise and described system is abandoned wind state be modeled as wind state of the abandoning linear function being relevant to binary variable, and based on this binary variable, the relevant operational state of hydroenergy storage stations all in system is modeled as the linear function of the described power station running status being relevant to described binary variable, linear function associating constrained system described in two, the system that simulates to be exerted oneself state abandoning hydroenergy storage station under wind state.

2. hydroenergy storage station characteristic accurate analog method according to claim 1 is characterized in that: described in abandon the linear function table of wind state and levy: the wind-powered electricity generation theoretical prediction power of t period system wind power is received with t period system difference for being relevant to binary variable A ^tlinear function, receive the size of wind power to judge currently to abandon wind state according to system.

3. hydroenergy storage station characteristic accurate analog method according to claim 2 is characterized in that: described in abandon wind state linear function mathematic(al) representation be:

$(A^t-1)\×P_w^{t,n*}+eps\≤P_w^{t,n*}-P_w^t\≤A^t\×P_w^{t,n*}---\left(1\right)$

In formula: when system occurs abandoning wind, A ^t=1; When system does not occur abandoning wind, A ^twhen=0; for non-negative continuous variable; for constant; Eps is positive dimensionless.

4. hydroenergy storage station characteristic accurate analog method according to claim 1 is characterized in that: the linear function of described power station running status at least comprises output of power station state linear function, and the linear function table of this output of power station state is levied: t period pump storage plant generator power for being relevant to binary variable A ^tlinear function, when system occurs abandoning wind, exerting oneself of hydroenergy storage station is zero, and when system does not occur abandoning wind, hydroenergy storage station is exerted oneself between theoretical maximum generation power and the minimum generated output of theory.

5. hydroenergy storage station characteristic accurate analog method according to claim 4 is characterized in that: the mathematic(al) representation of described system output of power station state linear function is:

$P_{PSG\min }^t\×(1-A^t)\≤P_{PSG}^t\≤P_{PSGmax}^t\×(1-A^t)---\left(2\right)$

In formula: for nonnegative variable, with for constant, represent hydroenergy storage station theoretical maximum generated output and minimum theoretical generated output respectively.

6. hydroenergy storage station characteristic accurate analog method according to claim 5 is characterized in that: the linear function of the running status in described power station also wraps the linear function of power station extraction water state, and the linear function of this extraction water state characterizes: draw water state or the state of discharging water of t period is followed successively by binary variable any time, a kind of operate condition that described power station can only occur extraction or discharge water.

7. hydroenergy storage station characteristic accurate analog method according to claim 6 is characterized in that: the mathematic(al) representation of the linear function of described power station t period extraction water state is:

$0\≤{\mathrm{PS}}_a^t+{\mathrm{PS}}_b^t\≤1---\left(3\right)$

In formula, represent that pumping operation is not carried out in power station, represent that pumping operation is carried out in power station; represent that the operation that discharges water is not carried out in power station, represent that the operation that discharges water is carried out in power station.

8. hydroenergy storage station characteristic accurate analog method according to claim 7 is characterized in that: the linear function of described power station running status also comprises the linear function of described output of power station constraint and power station and to draw water the linear function of state;

The linear function of described output of power station fluctuation status characterizes: when described power station discharges water and generates electricity, generated output is at the maximum generated output that discharges water of theory discharge water generated output minimum with theory between fluctuate arbitrarily, the mathematic(al) representation of this linear function is:

${\mathrm{PS}}_b^t\×P_{PSG\min }^t\≤P_{PSG}^t\≤{\mathrm{PS}}_b^t\×P_{PSGmax}^t---\left(4\right)$

Described power station draw water state linear function characterize: time drawing water in described power station, draw water unit once open, be just in full hair-like state, the mathematic(al) representation of this linear function is:

$\left\{\begin{array}{c}{P}_{PSP}^t={\mathrm{PS}}_a^t\×N_{PS}^t\×P_{PSPunit}^t\\ 0\≤N_{PS}^t\≤N_{PSmax}^t\end{array}\right.---\left(5\right)$

In formula, for nonnegative variable, represent that hydroenergy storage station is drawn water power; for positive integer variable, represent that the water pumper that hydroenergy storage station participates in optimizing is organized a performance number; for separate unit water pumper kludge capacity; for normal number, represent that hydroenergy storage station is drawn water unit head station number.

9. hydroenergy storage station characteristic accurate analog method according to claim 7 is characterized in that: the linear function of described power station running status also comprises the linear function of the linear function of set optimization power constraint, the linear function of minimum start and stop time-constrain, the linear function of heat supply phase thermal power plant unit units limits and start and stop logic state constraint

The linear function of described set optimization power constraint characterizes: certain region jth platform conventional power unit power output size for being relevant to binary variable linear function, the mathematic(al) representation of this linear function is:

$P_{j,\min }\·X_j^t\≤P_j^t\≤P_{j,max}\·X_j^t---\left(6\right)$

In formula: P _{j, max}, P _{j, min}for constant, be expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit; the binary variable of j platform unit in the running status of period t, represent that unit does not run, represent that unit runs;

The linear function of described minimum start and stop time-constrain characterizes: any time, unit j can only occur opening machine operate condition at period t or shutdown action state in one; The mathematic(al) representation of this linear function is:

$Y_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1---\left(7\right)$

Or $Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1---\left(8\right)$

In formula: for binary variable, represent that unit starts, represent that unit is not at starting state, represent that unit is shut down, represent that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; Dissimilar Unit Commitment machine time parameter is different;

The linear function of described heat supply phase thermal power plant unit units limits characterizes: the coupled relation that back pressure type thermal power plant unit and the electricity of bleeder thermal power plant unit within the heat supply phase are exerted oneself and heat is exerted oneself; The mathematic(al) representation of this linear function is:

$P_{j,BY}^t=C_j^b\·H_j^t---\left(9\right)$

$H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,max}-H_j^t\·C_j^v---\left(10\right)$

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; be respectively jth platform unit thermal power plant unit coupled thermomechanics coefficient, represent the coupling coefficient of lower limit of exerting oneself, represent the coupling coefficient of the upper limit of exerting oneself; for t period load of heat; P _{j, max}, P _{j, min}for constant, be expressed as the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit;

The linear function of described start and stop logic state constraint characterizes: operating states of the units, open machine state and running status meets logical constraint; The mathematic(al) representation of this linear function is:

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(11\right)$

10. hydroenergy storage station characteristic accurate analog method according to claim 7 is characterized in that: the running status in described power station also comprises the following linear function of following nonbinary variable modeling: the linear function of the linear function of the linear function that the linear function of power station capacity constrain, described power station water balance retrain, unit climbing rate constraint, the linear function of interregional line transmission capacity-constrained, spinning reserve constraint, the linear function of region account load balancing constraints, the linear function of wind power constraint and the linear function of target function

The linear function of described hydroenergy storage station capacity constrain characterizes: in hydroenergy storage station, water yield size need meet the constraint of self storage capacity size, and when drawing water, the water yield can not exceed the maximum amount of water that storage capacity allows; When discharging water, the water yield can not be less than the least quantity that storage capacity allows; The mathematic(al) representation of this linear function is:

$W^t-W_{max}\≤P_{PSG}^t\·{\mathrm{\η}}_g-P_{PSP}^t\·{\mathrm{\η}}_p\≤W^t-W_{min}---\left(12\right)$

In formula, W _maxand W _minfor constant, represent maximum/minimum reservoir storage of this hydroenergy storage station; W ^tfor positive variable, represent the water yield of this hydroenergy storage station current time; η _gand η _pbe constant, the hydroenergy storage station water yield/electricity conversion coefficient when expression discharges water and draws water respectively.

The linear function of described hydroenergy storage station water balance constraint characterizes: the corresponding relation of storage capacity and its energy output in hydroenergy storage station; The mathematic(al) representation of this linear function is:

$\left\{\begin{array}{c}{P}_{PSG}^t\·{\mathrm{\η}}_g-P_{PSP}^t\·{\mathrm{\η}}_p=W^t-W_{end}^t\\ W^t=W_{end}^{t-1}\end{array}\right.---\left(13\right)$

W ¹＝W _InitialCap(14)

In formula, for positive variable, represent that this hydroenergy storage station stops electricity, the initial quantity of electricity of its value and subsequent time identical; W ¹for hydroenergy storage station first period storage capacity water yield size, its value should be the initial storage in power station, W _initialCapfor the hydroenergy storage station initial storage water yield;

The linear function of the linear function of described unit climbing rate constraint characterizes: what every platform unit unit interval can increase or reduce exerts oneself; The mathematic(al) representation of this linear function is:

$P_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}---\left(15\right)$

$P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}---\left(16\right)$ In formula: Δ P _{j, up}, Δ P _{j, down}be respectively the maximum upper creep speed of unit j and lower creep speed;

The linear function of described interregional line transmission capacity-constrained characterizes: between region, circuit allows the power of transmission can not exceed its physical restriction; The mathematic(al) representation of this linear function is:

$-L_{i,max}\≤L_i^t\≤L_{i,max}---\left(17\right)$

for the transmitted power of t period i-th transmission lines; And L _{i, max}with-L _{i, max}be respectively the i-th transmission lines transmission capacity bound; Setting current reference direction is: inflow region is positive direction, and outflow region is negative direction.So can positive and negative values be got, positive and negative, represent the direction of power delivery

The linear function of the linear function of described spinning reserve constraint characterizes: in order to realize the balance of system active power, system should have certain reserve capacity.The reserve capacity of system refers to that, in system peak load situation, the power available capacity of system is greater than the part of generation load.; The mathematic(al) representation of this linear function is:

$-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,max}\·X_j^t\≤-\underset{n=1}{\overset{N}{\Σ}}{P}_l^{t,n}-\mathrm{Pr}e---\left(18\right)$

In formula: then represent the electric load of n region t period, Pre is that positive rotation is for subsequent use;

The linear function of described region account load balancing constraints characterizes: electric power system generated output and station service power load should Real-time Balancings; The mathematic(al) representation of this linear function is:

$P_{all}^{t,n}+P_w^{t,n}+P_{PSG}^t-P_{PSP}^t+L_i^t=P_l^{t,n}---\left(19\right)$

In formula: for fired power generating unit gross capability, for the wind power that electrical network is received

The linear function of described wind power constraint characterizes: the wind power that electrical network is received should be less than its predicted power; The mathematic(al) representation of this linear function is:

$0\≤P_w^{t,n}\≤P_w^{t,n^{*}}---\left(20\right)$

In formula: wind-powered electricity generation predicted power size

Described target function characterizes: electrical network Preferred Acceptance wind-powered electricity generation electric power; The mathematic(al) representation of this function is:

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_w^{t,n}---\left(21\right).$