CN104167762A - Combination optimization method of unit with chemical energy storage system and intermittent power source - Google Patents

Abstract

The invention provides a combination optimization method of a unit with a chemical energy storage system and an intermittent power source. The intermittent power source comprises a wind power generator unit and a photovoltaic power generator unit, and the chemical energy storage system, the intermittent power source and a thermal generator unit form an electrical power system jointly. Under the condition that the chemical energy storage system and the intermittent power source are assessed, a unit combination and output plan of all time frames of the coming day is made. A robust optimization method is adopted, the forecast error and the chemical energy storage output feature of the power of the intermittent power source are fully considered, the access capability of the intermittent power source is improved, and operation safety of the system is fully guaranteed.

Description

A kind of Unit Combination optimization method containing chemical energy storage system and fitful power

Technical field

The present invention relates to a kind of optimization method, be specifically related to a kind of Unit Combination optimization method containing chemical energy storage system and fitful power.

Background technology

Unit Combination method be take system cost minimum as target, according to short-term load forecasting, adjusts in real time unit output, meets system safety constraint, and the fitful power of fully dissolving, to realize energy-conservation Unit Combination.Unit Combination had obtained successful Application in Real-Time Scheduling field abroad in recent years, almost can complete all analytical works of electrical production and scheduling with Unit Combination related tool, and had also become at home research and application focus.

Chemical energy storage comprises lead-acid battery, lithium ion battery, flow battery, sodium-sulphur battery etc.Flow battery has the potentiality of extensive energy storage, but current most popular or lead-acid battery.The application of chemical energy storage technology in electric power system is still in the starting stage, but along with the proportion day of fitful power in electrical network more increases, having come into operation into trend of the times of chemical energy storage, how in operation plan, making full use of the access capability that chemical energy storage system improves fitful power is the hot issue of research.

Summary of the invention

In order to overcome above-mentioned the deficiencies in the prior art, the invention provides a kind of Unit Combination optimization method containing chemical energy storage system and fitful power, in the situation that considering chemical energy storage and fitful power access, Unit Combination and the plan of exerting oneself of the establishment day part of following a day.Adopt robust Optimal methods, take into full account predicated error and the chemical energy storage power producing characteristics of fitful power power, improve the access capability of fitful power, and the fail safe of fully assurance system operation.

In order to realize foregoing invention object, the present invention takes following technical scheme:

The invention provides a kind of Unit Combination optimization method containing chemical energy storage system and fitful power, described fitful power comprises wind turbine generator and photovoltaic generation unit, and chemical energy storage system, fitful power and thermal power generation unit form electric power system jointly; Said method comprising the steps of:

Step 1: set up the Unit Combination Mathematical Modeling of considering chemical energy storage characteristic and fitful power receiving ability;

Step 2: initialization clearance-type power supply actual power, and iterations k=0 is set;

Step 3: reset iterations, make k=k+1, before calculating, under all clearance-type power supply actual power scenes, take the target function minimum value α exerting oneself as variable, and calculate and under this condition of exerting oneself, take the target function maximum β that clearance-type power supply actual power is variable;

Step 4: if the minimum value α of target function maximum β and target function difference be greater than convergence precision ε, return to step 3; Otherwise finish, the Unit Combination result now obtaining is final optimization pass result.

In described step 1, the target function of Unit Combination Mathematical Modeling is as follows:

$\begin{array}{c}\underset{P_{i,h,}{I}_{i,h,}{P}_{f,m,h}{P}_{\mathrm{stor},s,h}}{\min }\underset{P_{f,m,h}^{\mathrm{actual}}}{\sup }f=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}\underset{h=1}{\overset{H}{\mathrm{\Σ}}}[F_{\mathrm{ci}}(P_{i,h},I_{i,h})+{\mathrm{SU}}_{i,h}(I_{i,h},I_{i,h-1}...)]+\underset{m=1}{\overset{W}{\mathrm{\Σ}}}\underset{h=1}{\overset{H}{\mathrm{\Σ}}}(P_{f,m,h}\×F_{\mathrm{wind}})\\ +N\×\underset{m=1}{\overset{W}{\mathrm{\Σ}}}\underset{h=1}{\overset{H}{\mathrm{\Σ}}}(P_{f,m,h}^{\mathrm{actual}}-P_{f,m,h})\end{array}---\left(1\right)$

Wherein, i is fired power generating unit sequence number, and h is Unit Combination period sequence number, and m is clearance-type power supply unit sequence number, and s is chemical energy storage system sequence number, and f is target function, I _i,hfor the state of fired power generating unit i in the h period, 0 for shutting down, and 1 is start; P _i,hfor fired power generating unit i is at h value of exerting oneself of period, P _{f, m, h}for the power planning value of fitful power unit m in the h period, for the power actual value of fitful power unit m in the h period, P _{stor, s, h}for the performance number of chemical energy storage system s in the h period;

for fired power generating unit fuel and start expense, NG is fired power generating unit sum, and H is Unit Combination hop count when total, F _ci(P _i,h, I _i,h) be the fuel cost function of fired power generating unit i, SU _i,h(I _i,h, I _{i, h-1}) be that fired power generating unit i is at the start cost function of period h;

for fitful power power purchase expense, W is fitful power unit sum, F _windfor fitful power rate for incorporation into the power network;

for system is abandoned fitful power penalty function, N is penalty function coefficient.

The constraints that described target function is corresponding comprises and discharges and recharges power constraint, power capacity constraint, power-electric energy equality constraint, chemical energy storage system climbing constraint, the constraint of chemical energy storage system reserve, plans last period electric energy constraint, fired power generating unit units limits, system loading Constraints of Equilibrium, fired power generating unit climbing constraint, the constraint of fitful power power planning value and fired power generating unit startup-shutdown and retrain.

The constraint representation that discharges and recharges about chemical energy storage system is:

${-P}_{\mathrm{stor},s,h}^{\mathrm{cha},\max }\≤P_{\mathrm{stor},s,h}\≤P_{\mathrm{stor},s,h}^{\mathrm{discha},\max }(s=1,...,S;h=1,...,H)---\left(2\right)$

Wherein, for the maximum charge power of chemical energy storage system s in the h period, for the maximum discharge power of chemical energy storage system s in the h period, S is chemical energy storage system sum;

Power capacity constraint representation about chemical energy storage system is:

η _stor,s,min×C _stor,s,max≤C _stor,s,h≤η _stor,s,max×C _stor,s,max(s＝1,…,S；h＝1,…,H) (3)

Wherein, η _{stor, s, min}and η _{stor, s, max}be respectively electric energy lower limit proportionality coefficient and the power upper limit proportionality coefficient of chemical energy storage system s, C _{stor, s, max}for the power capacity of chemical energy storage system s, C _{stor, s, h}for the storage of electrical energy value of chemical energy storage system s in the h period;

Power-electric energy equality constraint about chemical energy storage system is expressed as (4):

$\left\{\begin{array}{c}{C}_{\mathrm{stor},s,h-1}-{\mathrm{\η}}_s^{\mathrm{discha}}\×P_{\mathrm{stor},s,h}\×1_H(P_{\mathrm{stor},s,h}>0)\\ C_{\mathrm{stor},s,h-1}-{\mathrm{\η}}_s^{\mathrm{cha}}\×P_{\mathrm{stor},s,h}\×1_H(P_{\mathrm{stor},s,h}\≤0)\end{array}\right.=C_{\mathrm{stor},s,h}(s=1,...,S;h=1,...,H)---\left(4\right)$

Wherein, C _{stor, s, h-1}for the storage of electrical energy value of chemical energy storage system s in the h-1 period, for the discharge energy efficiency factor of chemical energy storage system s, for the rechargeable energy efficiency factor of chemical energy storage system s, P _{stor, s, h}for the performance number of chemical energy storage system s in the h period, 1 _hit is 1 chronomere;

Chemical energy storage system climbing constraint representation is:

$\left\{\begin{array}{c}{P}_{\mathrm{stor},s,h}-P_{\mathrm{stor},s,h-1}\≤{\mathrm{US}}_s\\ P_{\mathrm{stor},s,h-1}-P_{\mathrm{stor},s,h}\≤{\mathrm{DS}}_s\end{array}\right.(s=1,...,S;h=1,...,H)---\left(5\right)$

Wherein, P _{stor, s, h-1}chemical energy storage system s is at the performance number of h-1 period, US _sand DS _sbe respectively creep speed and the landslide speed of chemical energy storage system s;

Chemical energy storage system reserve constraint representation is:

$\left\{\begin{array}{c}{R}_{\mathrm{stor},\mathrm{up},s,h}=\min \{P_{\mathrm{stor},s,h}^{\mathrm{discha},\max }-P_{\mathrm{stor},s,h},\frac{C_{\mathrm{stor},s,h}-{\mathrm{\η}}_{\mathrm{stor},s,\min }\×C_{\mathrm{stor},s,\max }}{1_H}\}\\ R_{\mathrm{stor},\mathrm{down},s,h}=\min \{P_{\mathrm{stor},s,h}+P_{\mathrm{stor},s,h}^{\mathrm{cha},\max },\frac{{\mathrm{\η}}_{\mathrm{stor},s,\max }\×C_{\mathrm{stor},s,\max }-C_{\mathrm{stor},s,h}}{1_H}\end{array}\right.(s=1,...,S;h=1,...,H)---\left(6\right)$

Wherein, R _{stor, up, s, h}and R _{stor, down, s, h}for rise reserve capacity and the downward reserve capacity that chemical energy storage system s provided in the h period, C _{stor, s, h}for the storage of electrical energy value of chemical energy storage system s in the h period;

Requirement according to chemical energy storage system to reserve capacity, the constraint of chemical energy storage system reserve is specifically expressed as:

$\left\{\begin{array}{c}\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}(P_{i,\max }-P_{i,h})\×I_{i,h}+R_{\mathrm{stor},\mathrm{up},s,h}\≥R_{\mathrm{up},h}\\ \underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}(P_{i,h}-P_{i,\min })\×I_{i,h}+R_{\mathrm{stor},\mathrm{down},s,h}\≥R_{\mathrm{down},h}\end{array}\right.(h=\mathrm{1,2},...,H)---\left(7\right)$

Wherein, P _{i, max}and P _{i, min}be respectively the fired power generating unit i upper and lower bound of exerting oneself, R _{up, h}and R _{down, h}be respectively the rise reserve capacity of chemical energy storage system in the h period and require and lower reserve capacity requirement;

Electric energy constraint representation of last period of plan about chemical energy storage system is:

$C_{\mathrm{stor},s,h}\≥C_{\mathrm{stor},s,\max }\×{\mathrm{\η}}_{\mathrm{stor},s,\min }^{\mathrm{cap}}(s=1,...,S;h=H)---\left(8\right)$

Wherein, for the last period energy storage of chemical energy storage system system s electric energy proportionality coefficient.

Fired power generating unit units limits is expressed as:

P _i,min×I _i,h≤P _i,h≤P _i,max×I _i,h(i＝1,…,NG；h＝1,…,H) (9)

Wherein, P _{i, max}and P _{i, min}be respectively the fired power generating unit i upper and lower bound of exerting oneself;

System loading Constraints of Equilibrium is expressed as:

$\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}{P}_{i,h}\×I_{i,h}+\underset{m=1}{\overset{W}{\mathrm{\Σ}}}{P}_{f,m,h}+\underset{s=1}{\overset{S}{\mathrm{\Σ}}}{P}_{\mathrm{stor},s,h}=P_{D,h}---\left(10\right)$

Wherein, P _d,hfor the load value of system in the h period;

Fired power generating unit climbing constraint representation is:

$\left\{\begin{array}{c}{P}_{i,h}-P_{i,h-1}\≤[1-I_{i,h}(1-I_{i,h-1})]{\mathrm{UR}}_i+I_{i,h}(1-I_{i,h-1})P_{i,\min }\\ P_{i,h-1}-P_{i,h}\≤[1-I_{i,h-1}(1-I_{i,h})]{\mathrm{DR}}_i+I_{i,h-1}(1-I_{i,h})P_{i,\min }\end{array}\right.(i=1,...,\mathrm{NG};h=1,...,H)---\left(11\right)$

Wherein, P _{i, h-1}for fired power generating unit i is at h-1 value of exerting oneself of period, I _{i, h-1}for the state of fired power generating unit i in the h-1 period, 0 for shutting down, and 1 is start; UR _iand DR _ibe respectively creep speed and the landslide speed of fired power generating unit i;

Fitful power power planning value constraint representation is:

$\left\{\begin{array}{c}0\≤P_{f,m,h}\≤P_{f,m,h}^{\mathrm{forecast}}\\ P_{f,m,h}^{\mathrm{forecast}}+P_{\mathrm{given}}\≥P_{f,m,h}^{\mathrm{actual}}\\ P_{f,m}^{\mathrm{cap}}\≥P_{f,m,h}^{\mathrm{actual}}\≥0\end{array},\≥P_{f,m,h}^{\mathrm{forecast}}-P_{\mathrm{given}}\right.(m=1,...,W;h=1,...,H)---\left(12\right)$

Wherein, for the power prediction value of fitful power unit m in the h period, P _givenfor the fitful power power of the assembling unit predicated error upper limit, the installed capacity of fitful power unit m;

Fired power generating unit startup-shutdown constraint representation is:

$\left\{\begin{array}{c}[X_{i,h-1}^{\mathrm{on}}-T_i^{\mathrm{on}}]\×[I_{i,h-1}-I_{i,h}]\≥0\\ [X_{i,h-1}^{\mathrm{off}}-T_i^{\mathrm{off}}]\×[I_{i,h-1}-I_{i,h}]\≥0\end{array}\right.(i=1,...,\mathrm{NG};h=1,...,H)---\left(13\right)$

Wherein, with be respectively fired power generating unit i hop count and hop count while having shut down when the start of h-1 period; with hop count while being respectively minimum when start hop count of fired power generating unit i and minimum shutdown.

In described step 2, initialization clearance-type power supply actual power,

$P_{f,m,h}^{\mathrm{actual}}=P_{f,m,h}^{\mathrm{forecast}}(m=\mathrm{1,2},...,W;h=\mathrm{1,2},...,H)---\left(14\right)$

Wherein, for the power actual value of fitful power unit m in the h period, for the power prediction value of fitful power unit m in the h period, W is fitful power unit sum, and H is Unit Combination hop count when total, and iterations k=0 is set.

Described step 3 comprises the following steps:

Step 3-1: the Unit Commitment and exerting oneself as the target function minimum value α of variable of take under all clearance-type power supply actual power scenes is expressed as:

$\begin{array}{c}\underset{P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k}{\min }\mathrm{\α}\\ s.t.f(P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k,P_{f,m,h}^{\mathrm{actual},j})\≤\mathrm{\α},j=\mathrm{0,1,2},...,k-1\end{array}---\left(15\right)$

Wherein, be that the fired power generating unit i of the k time iteration is in h value of exerting oneself of period; be the fired power generating unit i of the k time iteration at the state of h period, 0 for shutting down, and 1 is start; be the power planning value of the fitful power unit m of the k time iteration in the h period; be the performance number of the chemical energy storage system s of the k time iteration in the h period; be the power actual value of the fitful power unit m of the j time iteration in the h period;

Step 3-2: the target function maximum β that the clearance-type power supply actual power of take under this condition of exerting oneself is variable is expressed as:

$\mathrm{\β}=\underset{P_{f,m,h}^{\mathrm{actual},k}}{\max }f(P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k,P_{f,m,h}^{\mathrm{actual},k})---\left(16\right)$

Wherein, be the power actual value of the fitful power unit m of the k time iteration in the h period.

Compared with prior art, beneficial effect of the present invention is:

1. considered chemical energy storage characteristic, comprised and discharge and recharge power upper limit, efficiency for charge-discharge etc., realized the optimization that energy storage is exerted oneself and calculated;

2. in target function, increase fitful power and abandon penalty function, at utmost improved the ability that system is received fitful power;

Robust optimization to initial point require lowly, do not need to carry out analog simulation through enumerating a large amount of scenes, fast convergence rate, is to solve containing a large amount of chemical energy storages and the efficient accurate method of the lower Unit Combination of fitful power access.

Accompanying drawing explanation

Fig. 1 contains the Unit Combination optimization method flow chart of chemical energy storage system and fitful power in the embodiment of the present invention;

Fig. 2 is that in the embodiment of the present invention, chemical energy storage system is filled (putting) electrical power upper limit schematic diagram.

Embodiment

Below in conjunction with accompanying drawing, the present invention is described in further detail.

The present invention adopts Shandong side optimized algorithm, takes into full account the fluctuation of chemical energy storage characteristic and fitful power, take and guarantees that systematically fail safe is as prerequisite, and at utmost raising system is received the ability of fitful power, realizes the Efficient Solution of Unit Combination.

As Fig. 1, the invention provides a kind of Unit Combination optimization method containing chemical energy storage system and fitful power, described fitful power comprises wind turbine generator and photovoltaic generation unit, and chemical energy storage system, fitful power and thermal power generation unit form electric power system jointly; Said method comprising the steps of:

Step 1: set up the Unit Combination Mathematical Modeling of considering chemical energy storage characteristic and fitful power receiving ability;

Step 2: initialization clearance-type power supply actual power, and iterations k=0 is set;

Step 3: reset iterations, make k=k+1, before calculating, under all clearance-type power supply actual power scenes, take the target function minimum value α exerting oneself as variable, and calculate and under this condition of exerting oneself, take the target function maximum β that clearance-type power supply actual power is variable;

Step 4: if the minimum value α of target function maximum β and target function difference be greater than convergence precision ε, return to step 3; Otherwise finish, the Unit Combination result now obtaining is final optimization pass result.

In described step 1, the target function of Unit Combination Mathematical Modeling is as follows:

$\begin{array}{c}\underset{P_{i,h,}{I}_{i,h,}{P}_{f,m,h}{P}_{\mathrm{stor},s,h}}{\min }\underset{P_{f,m,h}^{\mathrm{actual}}}{\sup }f=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}\underset{h=1}{\overset{H}{\mathrm{\Σ}}}[F_{\mathrm{ci}}(P_{i,h},I_{i,h})+{\mathrm{SU}}_{i,h}(I_{i,h},I_{i,h-1}...)]+\underset{m=1}{\overset{W}{\mathrm{\Σ}}}\underset{h=1}{\overset{H}{\mathrm{\Σ}}}(P_{f,m,h}\×F_{\mathrm{wind}})\\ +N\×\underset{m=1}{\overset{W}{\mathrm{\Σ}}}\underset{h=1}{\overset{H}{\mathrm{\Σ}}}(P_{f,m,h}^{\mathrm{actual}}-P_{f,m,h})\end{array}---\left(1\right)$

Wherein, i is fired power generating unit sequence number, and h is Unit Combination period sequence number, and m is clearance-type power supply unit sequence number, and s is chemical energy storage system sequence number, and f is target function, I _i,hfor the state of fired power generating unit i in the h period, 0 for shutting down, and 1 is start; P _i,hfor fired power generating unit i is at h value of exerting oneself of period, P _{f, m, h}for the power planning value of fitful power unit m in the h period, for the power actual value of fitful power unit m in the h period, P _{stor, s, h}for the performance number of chemical energy storage system s in the h period;

for fired power generating unit fuel and start expense, NG is fired power generating unit sum, and H is Unit Combination hop count when total, F _ci(P _i,h, I _i,h) be the fuel cost function of fired power generating unit i, SU _i,h(I _i,h, I _{i, h-1}) be that fired power generating unit i is at the start cost function of period h;

for fitful power power purchase expense, W is fitful power unit sum, F _windfor fitful power rate for incorporation into the power network;

for system is abandoned fitful power penalty function, N is penalty function coefficient.

The constraints that described target function is corresponding comprises and discharges and recharges power constraint, power capacity constraint, power-electric energy equality constraint, chemical energy storage system climbing constraint, the constraint of chemical energy storage system reserve, plans last period electric energy constraint, fired power generating unit units limits, system loading Constraints of Equilibrium, fired power generating unit climbing constraint, the constraint of fitful power power planning value and fired power generating unit startup-shutdown and retrain.

Chemical energy storage system discharges and recharges the control that is subject to its inverter that connects, and affected by the factors such as chemical energy storage system type, chemical energy storage system power value.The power characteristic that discharges and recharges of chemical energy storage system is given, and available segment broken line represents.In model, curve abscissa is corresponding to state-of-charge (the State of Charge of energy-storage system, SOC), also be the percentage that chemical energy storage system storage electric energy accounts for power capacity, ordinate is corresponding to maximum charge or the discharge power value of chemical energy storage system, and therefore the constraint representation that discharges and recharges about chemical energy storage system is:

${-P}_{\mathrm{stor},s,h}^{\mathrm{cha},\max }\≤P_{\mathrm{stor},s,h}\≤P_{\mathrm{stor},s,h}^{\mathrm{discha},\max }(s=1,...,S;h=1,...,H)---\left(2\right)$

Wherein, for the maximum charge power of chemical energy storage system s in the h period, for the maximum discharge power of chemical energy storage system s in the h period, S is chemical energy storage system sum;

For each chemical energy storage system, the power characteristic that discharges and recharges of broken line that can be by accompanying drawing 2 represents, wherein SOC is taken as fixed value, is 10%, 20% ... 100%.

Under actual conditions, chemical energy storage system discharge and recharge power characteristic difference, need the power characteristic parameter of respectively given chemical energy storage system charging, electric discharge with ratio between the maximum charge and discharge performance number of its representative energy-storage system in each state-of-charge situation and the specified charge and discharge performance number of chemical energy storage system.For the chemical energy storage system being formed by dissimilar chemical cell, discharge and recharge the stackable formation of power characteristic, corresponding charging, discharge power characteristic curve parameter are respectively with $\mathrm{\Σ}{P0}_{\mathrm{discha},\mathrm{pra}}^{\mathrm{stor}},...,\mathrm{\Σ}{P10}_{\mathrm{discha},\mathrm{pra}}^{\mathrm{stor}}.$

The life-span of chemical energy storage system with its discharge and recharge number of times, to discharge and recharge the degree of depth relevant, in order to extend the useful life of energy-storage system as far as possible, power capacity scope that generally can given chemical energy storage system operation, about the power capacity constraint representation of chemical energy storage system is:

η _stor,s,min×C _stor,s,max≤C _stor,s,h≤η _stor,s,max×C _stor,s,max(s＝1,…,S；h＝1,…,H)(3)

Wherein, η _{stor, s, min}and η _{stor, s, max}be respectively electric energy lower limit proportionality coefficient and the power upper limit proportionality coefficient of chemical energy storage system s, C _{stor, s, max}for the power capacity of chemical energy storage system s, C _{stor, s, h}for the storage of electrical energy value of chemical energy storage system s in the h period;

Self there is certain electric energy loss in chemical energy storage system, this can use the rechargeable energy efficiency factor of chemical energy storage system with discharge energy efficiency factor represent.By energy-storage system from A% * C _{stor, s, max}be charged to B% * C _{stor, s, max}required electric energy is C ₁, now have B% * C _{stor, s, max}-A% * C _{stor, s, max}< C ₁, reduce the charging loss that part is chemical energy storage system, now have in like manner, if by chemical energy storage system power from B% * C _{stor, s, max}discharge into A% * C _{stor, s, max}, release electric energy is C ₂, now have B% * C _{stor, s, max}-A% * C _{stor, s, max}> C ₂, reduce the discharge loss that part is chemical energy storage system, now have power-electric energy equality constraint about chemical energy storage system is expressed as (4):

$\left\{\begin{array}{c}{C}_{\mathrm{stor},s,h-1}-{\mathrm{\&eta;}}_s^{\mathrm{discha}}\&times;P_{\mathrm{stor},s,h}\&times;1_H(P_{\mathrm{stor},s,h}>0)\\ C_{\mathrm{stor},s,h-1}-{\mathrm{\&eta;}}_s^{\mathrm{cha}}\&times;P_{\mathrm{stor},s,h}\&times;1_H(P_{\mathrm{stor},s,h}\&le;0)\end{array}\right.=C_{\mathrm{stor},s,h}(s=1,...,S;h=1,...,H)---\left(4\right)$

Wherein, C _{stor, s, h-1}for the storage of electrical energy value of chemical energy storage system s in the h-1 period, for the discharge energy efficiency factor of chemical energy storage system s, for the rechargeable energy efficiency factor of chemical energy storage system s, P _{stor, s, h}for the performance number of chemical energy storage system s in the h period, 1 _hit is 1 chronomere;

The similar fired power generating unit of power climbing constraint of chemical energy storage system, difference is that its performance number can be for negative, and energy-storage system may operate in charged state and occurs as system loading, and chemical energy storage system climbing constraint representation is:

$\left\{\begin{array}{c}{P}_{\mathrm{stor},s,h}-P_{\mathrm{stor},s,h-1}\&le;{\mathrm{US}}_s\\ P_{\mathrm{stor},s,h-1}-P_{\mathrm{stor},s,h}\&le;{\mathrm{DS}}_s\end{array}\right.(s=1,...,S;h=1,...,H)---\left(5\right)$

Wherein, P _{stor, s, h-1}chemical energy storage system s is at the performance number of h-1 period, US _sand DS _sbe respectively creep speed and the landslide speed of chemical energy storage system s;

The full discharge condition of middle cell of chemical energy storage system is Millisecond to the fringe time of fully charged state.Energy-storage system consists of a large amount of battery cells, consider inverter that energy-storage system is supporting and control the time delay of communicator, energy-storage system can not surpass 1～2s from completely putting (filling) electricity condition to the fringe time that completely fills (putting) electricity condition yet.Therefore,, for chemical energy storage system, its power creep speed and landslide speed are much larger than fired power generating unit.

The reserve capacity that chemical energy storage system can provide is subject to its maximum in the different periods to discharge and recharge performance number on the one hand to be affected, be subject on the other hand chemical energy storage system in the impact of different period institute's storage of electrical energy values, the two has determined the reserve capacity value that chemical energy storage system can provide in the different periods jointly, and chemical energy storage system reserve constraint representation is:

$\left\{\begin{array}{c}{R}_{\mathrm{stor},\mathrm{up},s,h}=\min \{P_{\mathrm{stor},s,h}^{\mathrm{discha},\max }-P_{\mathrm{stor},s,h},\frac{C_{\mathrm{stor},s,h}-{\mathrm{\&eta;}}_{\mathrm{stor},s,\min }\&times;C_{\mathrm{stor},s,\max }}{1_H}\}\\ R_{\mathrm{stor},\mathrm{down},s,h}=\min \{P_{\mathrm{stor},s,h}+P_{\mathrm{stor},s,h}^{\mathrm{cha},\max },\frac{{\mathrm{\&eta;}}_{\mathrm{stor},s,\max }\&times;C_{\mathrm{stor},s,\max }-C_{\mathrm{stor},s,h}}{1_H}\end{array}\right.(s=1,...,S;h=1,...,H)---\left(6\right)$

Wherein, R _{stor, up, s, h}and R _{stor, down, s, h}for rise reserve capacity and the downward reserve capacity that chemical energy storage system s provided in the h period, C _{stor, s, h}for the storage of electrical energy value of chemical energy storage system s in the h period;

Requirement according to chemical energy storage system to reserve capacity, the constraint of chemical energy storage system reserve is specifically expressed as:

$\left\{\begin{array}{c}\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}(P_{i,\max }-P_{i,h})\&times;I_{i,h}+R_{\mathrm{stor},\mathrm{up},s,h}\&GreaterEqual;R_{\mathrm{up},h}\\ \underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}(P_{i,h}-P_{i,\min })\&times;I_{i,h}+R_{\mathrm{stor},\mathrm{down},s,h}\&GreaterEqual;R_{\mathrm{down},h}\end{array}\right.(h=\mathrm{1,2},...,H)---\left(7\right)$

Wherein, P _{i, max}and P _{i, min}be respectively the fired power generating unit i upper and lower bound of exerting oneself, R _{up, h}and R _{down, h}be respectively the rise reserve capacity of chemical energy storage system in the h period and require and lower reserve capacity requirement;

In some cases, need to set up the electric energy constraint of chemical energy storage system in the last period of planning cycle, require to execute after the plan of discharging and recharging, chemical energy storage system has certain electrical power storage, about the electric energy constraint representation of last period of plan of chemical energy storage system, is:

$C_{\mathrm{stor},s,h}\&GreaterEqual;C_{\mathrm{stor},s,\max }\&times;{\mathrm{\&eta;}}_{\mathrm{stor},s,\min }^{\mathrm{cap}}(s=1,...,S;h=H)---\left(8\right)$

Wherein, for the last period energy storage of chemical energy storage system system s electric energy proportionality coefficient.

In the fired power generating unit of open state, exerting oneself should be between its minimum and maximum.Fired power generating unit units limits is expressed as:

P _i,min×I _i,h≤P _i,h≤P _i,max×I _i,h(i＝1,…,NG；h＝1,…,H) (9)

Wherein, P _{i, max}and P _{i, min}be respectively the fired power generating unit i upper and lower bound of exerting oneself;

In model, the power supply and demand of system is balance constantly, ignores the network loss value of system in model, is also that fired power generating unit, fitful power unit planned value sum should equate with the load value of corresponding period of system.System loading Constraints of Equilibrium is expressed as:

$\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{P}_{i,h}\&times;I_{i,h}+\underset{m=1}{\overset{W}{\mathrm{\&Sigma;}}}{P}_{f,m,h}+\underset{s=1}{\overset{S}{\mathrm{\&Sigma;}}}{P}_{\mathrm{stor},s,h}=P_{D,h}---\left(10\right)$

Wherein, P _d,hfor the load value of system in the h period;

Fired power generating unit climbing constraint representation is:

$\left\{\begin{array}{c}{P}_{i,h}-P_{i,h-1}\&le;[1-I_{i,h}(1-I_{i,h-1})]{\mathrm{UR}}_i+I_{i,h}(1-I_{i,h-1})P_{i,\min }\\ P_{i,h-1}-P_{i,h}\&le;[1-I_{i,h-1}(1-I_{i,h})]{\mathrm{DR}}_i+I_{i,h-1}(1-I_{i,h})P_{i,\min }\end{array}\right.(i=1,...,\mathrm{NG};h=1,...,H)---\left(11\right)$

Wherein, P _{i, h-1}for fired power generating unit i is at h-1 value of exerting oneself of period, I _{i, h-1}for the state of fired power generating unit i in the h-1 period, 0 for shutting down, and 1 is start; UR _iand DR _ibe respectively creep speed and the landslide speed of fired power generating unit i;

Fitful power power planning value should be less than or equal to its predicted value.In addition, when the power deviation higher or lower than its planned value appears in fitful power, its minimum value can not be less than 0MW, the highlyest can not surpass its installed capacity, and fitful power power planning value constraint representation is:

$\left\{\begin{array}{c}0\&le;P_{f,m,h}\&le;P_{f,m,h}^{\mathrm{forecast}}\\ P_{f,m,h}^{\mathrm{forecast}}+P_{\mathrm{given}}\&GreaterEqual;P_{f,m,h}^{\mathrm{actual}}\\ P_{f,m}^{\mathrm{cap}}\&GreaterEqual;P_{f,m,h}^{\mathrm{actual}}\&GreaterEqual;0\end{array},\&GreaterEqual;P_{f,m,h}^{\mathrm{forecast}}-P_{\mathrm{given}}\right.(m=1,...,W;h=1,...,H)---\left(12\right)$

Wherein, for the power prediction value of fitful power unit m in the h period, P _givenfor the fitful power power of the assembling unit predicated error upper limit, the installed capacity of fitful power unit m;

Once its start of fired power generating unit or shutdown, in order to guarantee security reliability and the economy of its operation, hop count requirement when having corresponding minimum start or shutting down, fired power generating unit startup-shutdown constraint representation is:

$\left\{\begin{array}{c}[X_{i,h-1}^{\mathrm{on}}-T_i^{\mathrm{on}}]\&times;[I_{i,h-1}-I_{i,h}]\&GreaterEqual;0\\ [X_{i,h-1}^{\mathrm{off}}-T_i^{\mathrm{off}}]\&times;[I_{i,h-1}-I_{i,h}]\&GreaterEqual;0\end{array}\right.(i=1,...,\mathrm{NG};h=1,...,H)---\left(13\right)$

Wherein, with be respectively fired power generating unit i hop count and hop count while having shut down when the start of h-1 period; with hop count while being respectively minimum when start hop count of fired power generating unit i and minimum shutdown.

In described step 2, initialization clearance-type power supply actual power,

$P_{f,m,h}^{\mathrm{actual}}=P_{f,m,h}^{\mathrm{forecast}}(m=\mathrm{1,2},...,W;h=\mathrm{1,2},...,H)---\left(14\right)$

Wherein, for the power actual value of fitful power unit m in the h period, for the power prediction value of fitful power unit m in the h period, W is fitful power unit sum, and H is Unit Combination hop count when total, and iterations k=0 is set.

Described step 3 comprises the following steps:

Step 3-1: the Unit Commitment and exerting oneself as the target function minimum value α of variable of take under all clearance-type power supply actual power scenes is expressed as:

$\begin{array}{c}\underset{P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k}{\min }\mathrm{\&alpha;}\\ s.t.f(P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k,P_{f,m,h}^{\mathrm{actual},j})\&le;\mathrm{\&alpha;},j=\mathrm{0,1,2},...,k-1\end{array}---\left(15\right)$

Wherein, be that the fired power generating unit i of the k time iteration is in h value of exerting oneself of period; be the fired power generating unit i of the k time iteration at the state of h period, 0 for shutting down, and 1 is start; be the power planning value of the fitful power unit m of the k time iteration in the h period; be the performance number of the chemical energy storage system s of the k time iteration in the h period; be the power actual value of the fitful power unit m of the j time iteration in the h period;

Step 3-2: the target function maximum β that the clearance-type power supply actual power of take under this condition of exerting oneself is variable is expressed as:

$\mathrm{\&beta;}=\underset{P_{f,m,h}^{\mathrm{actual},k}}{\max }f(P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k,P_{f,m,h}^{\mathrm{actual},k})---\left(16\right)$

Wherein, be the power actual value of the fitful power unit m of the k time iteration in the h period.

Finally should be noted that: above embodiment is only in order to illustrate that technical scheme of the present invention is not intended to limit; those of ordinary skill in the field still can modify or be equal to replacement the specific embodiment of the present invention with reference to above-described embodiment; these do not depart from any modification of spirit and scope of the invention or are equal to replacement, within the claim protection range of the present invention all awaiting the reply in application.

Claims (7)

1. containing a Unit Combination optimization method for chemical energy storage system and fitful power, described fitful power comprises wind turbine generator and photovoltaic generation unit, and chemical energy storage system, fitful power and thermal power generation unit form electric power system jointly; It is characterized in that: said method comprising the steps of:

Step 1: set up the Unit Combination Mathematical Modeling of considering chemical energy storage characteristic and fitful power receiving ability;

Step 2: initialization clearance-type power supply actual power, and iterations k=0 is set;

Step 3: reset iterations, make k=k+1, before calculating, under all clearance-type power supply actual power scenes, take the target function minimum value α exerting oneself as variable, and calculate and under this condition of exerting oneself, take the target function maximum β that clearance-type power supply actual power is variable;

Step 4: if the minimum value α of target function maximum β and target function difference be greater than convergence precision ε, return to step 3; Otherwise finish, the Unit Combination result now obtaining is final optimization pass result.

2. the Unit Combination optimization method containing chemical energy storage system and fitful power according to claim 1, is characterized in that: in described step 1, the target function of Unit Combination Mathematical Modeling is as follows:

$\begin{array}{c}\underset{P_{i,h,}{I}_{i,h,}{P}_{f,m,h}{P}_{\mathrm{stor},s,h}}{\min }\underset{P_{f,m,h}^{\mathrm{actual}}}{\sup }f=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}\underset{h=1}{\overset{H}{\mathrm{\&Sigma;}}}[F_{\mathrm{ci}}(P_{i,h},I_{i,h})+{\mathrm{SU}}_{i,h}(I_{i,h},I_{i,h-1}...)]+\underset{m=1}{\overset{W}{\mathrm{\&Sigma;}}}\underset{h=1}{\overset{H}{\mathrm{\&Sigma;}}}(P_{f,m,h}\&times;F_{\mathrm{wind}})\\ +N\&times;\underset{m=1}{\overset{W}{\mathrm{\&Sigma;}}}\underset{h=1}{\overset{H}{\mathrm{\&Sigma;}}}(P_{f,m,h}^{\mathrm{actual}}-P_{f,m,h})\end{array}---\left(1\right)$

Wherein, i is fired power generating unit sequence number, and h is Unit Combination period sequence number, and m is clearance-type power supply unit sequence number, and s is chemical energy storage system sequence number, and f is target function, I _i,hfor the state of fired power generating unit i in the h period, 0 for shutting down, and 1 is start; P _i,hfor fired power generating unit i is at h value of exerting oneself of period, P _{f, m, h}for the power planning value of fitful power unit m in the h period, for the power actual value of fitful power unit m in the h period, P _{stor, s, h}for the performance number of chemical energy storage system s in the h period;

for fired power generating unit fuel and start expense, NG is fired power generating unit sum, and H is Unit Combination hop count when total, F _ci(P _i,h, I _i,h) be the fuel cost function of fired power generating unit i, SU _i,h(I _i,h, I _{i, h-1}) be that fired power generating unit i is at the start cost function of period h;

for fitful power power purchase expense, W is fitful power unit sum, F _windfor fitful power rate for incorporation into the power network;

for system is abandoned fitful power penalty function, N is penalty function coefficient.

3. the Unit Combination optimization method containing chemical energy storage system and fitful power according to claim 2, is characterized in that: the constraints that described target function is corresponding comprise discharge and recharge power constraint, power capacity constraint, power-electric energy equality constraint, chemical energy storage system climbing constraint, the constraint of chemical energy storage system reserve, plan last period electric energy constraint, fired power generating unit units limits, system loading Constraints of Equilibrium, fired power generating unit climbing constraint, fitful power power planning value retrains and fired power generating unit startup-shutdown retrains.

4. the Unit Combination optimization method containing chemical energy storage system and fitful power according to claim 3, is characterized in that: the constraint representation that discharges and recharges about chemical energy storage system is:

${-P}_{\mathrm{stor},s,h}^{\mathrm{cha},\max }\&le;P_{\mathrm{stor},s,h}\&le;P_{\mathrm{stor},s,h}^{\mathrm{discha},\max }(s=1,...,S;h=1,...,H)---\left(2\right)$

Wherein, for the maximum charge power of chemical energy storage system s in the h period, for the maximum discharge power of chemical energy storage system s in the h period, S is chemical energy storage system sum;

Power capacity constraint representation about chemical energy storage system is:

η _stor,s,min×C _stor,s,max≤C _stor,s,h≤η _stor,s,max×C _stor,s,max(s＝1,…,S；h＝1,…,H) (3)

Wherein, η _{stor, s, min}and η _{stor, s, max}be respectively electric energy lower limit proportionality coefficient and the power upper limit proportionality coefficient of chemical energy storage system s, C _{stor, s, max}for the power capacity of chemical energy storage system s, C _{stor, s, h}for the storage of electrical energy value of chemical energy storage system s in the h period;

Power-electric energy equality constraint about chemical energy storage system is expressed as (4):

$\left\{\begin{array}{c}{C}_{\mathrm{stor},s,h-1}-{\mathrm{\&eta;}}_s^{\mathrm{discha}}\&times;P_{\mathrm{stor},s,h}\&times;1_H(P_{\mathrm{stor},s,h}>0)\\ C_{\mathrm{stor},s,h-1}-{\mathrm{\&eta;}}_s^{\mathrm{cha}}\&times;P_{\mathrm{stor},s,h}\&times;1_H(P_{\mathrm{stor},s,h}\&le;0)\end{array}\right.=C_{\mathrm{stor},s,h}(s=1,...,S;h=1,...,H)---\left(4\right)$

Wherein, C _{stor, s, h-1}for the storage of electrical energy value of chemical energy storage system s in the h-1 period, for the discharge energy efficiency factor of chemical energy storage system s, for the rechargeable energy efficiency factor of chemical energy storage system s, P _{stor, s, h}for the performance number of chemical energy storage system s in the h period, 1 _hit is 1 chronomere;

Chemical energy storage system climbing constraint representation is:

$\left\{\begin{array}{c}{P}_{\mathrm{stor},s,h}-P_{\mathrm{stor},s,h-1}\&le;{\mathrm{US}}_s\\ P_{\mathrm{stor},s,h-1}-P_{\mathrm{stor},s,h}\&le;{\mathrm{DS}}_s\end{array}\right.(s=1,...,S;h=1,...,H)---\left(5\right)$

Wherein, P _{stor, s, h-1}chemical energy storage system s is at the performance number of h-1 period, US _sand DS _sbe respectively creep speed and the landslide speed of chemical energy storage system s;

Chemical energy storage system reserve constraint representation is:

$\left\{\begin{array}{c}{R}_{\mathrm{stor},\mathrm{up},s,h}=\min \{P_{\mathrm{stor},s,h}^{\mathrm{discha},\max }-P_{\mathrm{stor},s,h},\frac{C_{\mathrm{stor},s,h}-{\mathrm{\&eta;}}_{\mathrm{stor},s,\min }\&times;C_{\mathrm{stor},s,\max }}{1_H}\}\\ R_{\mathrm{stor},\mathrm{down},s,h}=\min \{P_{\mathrm{stor},s,h}+P_{\mathrm{stor},s,h}^{\mathrm{cha},\max },\frac{{\mathrm{\&eta;}}_{\mathrm{stor},s,\max }\&times;C_{\mathrm{stor},s,\max }-C_{\mathrm{stor},s,h}}{1_H}\end{array}\right.(s=1,...,S;h=1,...,H)---\left(6\right)$

Wherein, R _{stor, up, s, h}and R _{stor, down, s, h}for rise reserve capacity and the downward reserve capacity that chemical energy storage system s provided in the h period, C _{stor, s, h}for the storage of electrical energy value of chemical energy storage system s in the h period;

Requirement according to chemical energy storage system to reserve capacity, the constraint of chemical energy storage system reserve is specifically expressed as:

$\left\{\begin{array}{c}\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}(P_{i,\max }-P_{i,h})\&times;I_{i,h}+R_{\mathrm{stor},\mathrm{up},s,h}\&GreaterEqual;R_{\mathrm{up},h}\\ \underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}(P_{i,h}-P_{i,\min })\&times;I_{i,h}+R_{\mathrm{stor},\mathrm{down},s,h}\&GreaterEqual;R_{\mathrm{down},h}\end{array}\right.(h=\mathrm{1,2},...,H)---\left(7\right)$

Wherein, P _{i, max}and P _{i, min}be respectively the fired power generating unit i upper and lower bound of exerting oneself, R _{up, h}and R _{down, h}be respectively the rise reserve capacity of chemical energy storage system in the h period and require and lower reserve capacity requirement;

Electric energy constraint representation of last period of plan about chemical energy storage system is:

$C_{\mathrm{stor},s,h}\&GreaterEqual;C_{\mathrm{stor},s,\max }\&times;{\mathrm{\&eta;}}_{\mathrm{stor},s,\min }^{\mathrm{cap}}(s=1,...,S;h=H)---\left(8\right)$

Wherein, for the last period energy storage of chemical energy storage system system s electric energy proportionality coefficient.

5. the Unit Combination optimization method containing chemical energy storage system and fitful power according to claim 3, is characterized in that: fired power generating unit units limits is expressed as:

P _i,min×I _i,h≤P _i,h≤P _i,max×I _i,h(i＝1,…,NG；h＝1,…,H) (9)

Wherein, P _{i, max}and P _{i, min}be respectively the fired power generating unit i upper and lower bound of exerting oneself;

System loading Constraints of Equilibrium is expressed as:

$\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{P}_{i,h}\&times;I_{i,h}+\underset{m=1}{\overset{W}{\mathrm{\&Sigma;}}}{P}_{f,m,h}+\underset{s=1}{\overset{S}{\mathrm{\&Sigma;}}}{P}_{\mathrm{stor},s,h}=P_{D,h}---\left(10\right)$

Wherein, P _d,hfor the load value of system in the h period;

Fired power generating unit climbing constraint representation is:

$\left\{\begin{array}{c}{P}_{i,h}-P_{i,h-1}\&le;[1-I_{i,h}(1-I_{i,h-1})]{\mathrm{UR}}_i+I_{i,h}(1-I_{i,h-1})P_{i,\min }\\ P_{i,h-1}-P_{i,h}\&le;[1-I_{i,h-1}(1-I_{i,h})]{\mathrm{DR}}_i+I_{i,h-1}(1-I_{i,h})P_{i,\min }\end{array}\right.(i=1,...,\mathrm{NG};h=1,...,H)---\left(11\right)$

Wherein, P _{i, h-1}for fired power generating unit i is at h-1 value of exerting oneself of period, I _{i, h-1}for the state of fired power generating unit i in the h-1 period, 0 for shutting down, and 1 is start; UR _iand DR _ibe respectively creep speed and the landslide speed of fired power generating unit i;

Fitful power power planning value constraint representation is:

$\left\{\begin{array}{c}0\&le;P_{f,m,h}\&le;P_{f,m,h}^{\mathrm{forecast}}\\ P_{f,m,h}^{\mathrm{forecast}}+P_{\mathrm{given}}\&GreaterEqual;P_{f,m,h}^{\mathrm{actual}}\\ P_{f,m}^{\mathrm{cap}}\&GreaterEqual;P_{f,m,h}^{\mathrm{actual}}\&GreaterEqual;0\end{array},\&GreaterEqual;P_{f,m,h}^{\mathrm{forecast}}-P_{\mathrm{given}}\right.(m=1,...,W;h=1,...,H)---\left(12\right)$

Wherein, for the power prediction value of fitful power unit m in the h period, P _givenfor the fitful power power of the assembling unit predicated error upper limit, the installed capacity of fitful power unit m;

Fired power generating unit startup-shutdown constraint representation is:

$\left\{\begin{array}{c}[X_{i,h-1}^{\mathrm{on}}-T_i^{\mathrm{on}}]\&times;[I_{i,h-1}-I_{i,h}]\&GreaterEqual;0\\ [X_{i,h-1}^{\mathrm{off}}-T_i^{\mathrm{off}}]\&times;[I_{i,h-1}-I_{i,h}]\&GreaterEqual;0\end{array}\right.(i=1,...,\mathrm{NG};h=1,...,H)---\left(13\right)$

Wherein, with be respectively fired power generating unit i hop count and hop count while having shut down when the start of h-1 period; with hop count while being respectively minimum when start hop count of fired power generating unit i and minimum shutdown.

6. the Unit Combination optimization method containing chemical energy storage system and fitful power according to claim 1, is characterized in that: in described step 2, and initialization clearance-type power supply actual power,

$P_{f,m,h}^{\mathrm{actual}}=P_{f,m,h}^{\mathrm{forecast}}(m=\mathrm{1,2},...,W;h=\mathrm{1,2},...,H)---\left(14\right)$

Wherein, for the power actual value of fitful power unit m in the h period, for the power prediction value of fitful power unit m in the h period, W is fitful power unit sum, and H is Unit Combination hop count when total, and iterations k=0 is set.

7. the Unit Combination optimization method containing chemical energy storage system and fitful power according to claim 1, is characterized in that: described step 3 comprises the following steps:

Step 3-1: the Unit Commitment and exerting oneself as the target function minimum value α of variable of take under all clearance-type power supply actual power scenes is expressed as:

$\begin{array}{c}\underset{P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k}{\min }\mathrm{\&alpha;}\\ s.t.f(P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k,P_{f,m,h}^{\mathrm{actual},j})\&le;\mathrm{\&alpha;},j=\mathrm{0,1,2},...,k-1\end{array}---\left(15\right)$

Wherein, be that the fired power generating unit i of the k time iteration is in h value of exerting oneself of period; be the fired power generating unit i of the k time iteration at the state of h period, 0 for shutting down, and 1 is start; be the power planning value of the fitful power unit m of the k time iteration in the h period; be the performance number of the chemical energy storage system s of the k time iteration in the h period; be the power actual value of the fitful power unit m of the j time iteration in the h period;

Step 3-2: the target function maximum β that the clearance-type power supply actual power of take under this condition of exerting oneself is variable is expressed as:

$\mathrm{\&beta;}=\underset{P_{f,m,h}^{\mathrm{actual},k}}{\max }f(P_{i,h}^k,I_{i,h}^k,P_{f,m,h}^k,P_{\mathrm{stor},s,h}^k,P_{f,m,h}^{\mathrm{actual},k})---\left(16\right)$

Wherein, be the power actual value of the fitful power unit m of the k time iteration in the h period.