CN102496962B - Method for identifying and controlling wind power consumption capability of power system under peak load and frequency regulation constraints - Google Patents

Abstract

The invention discloses a method for identifying and controlling the wind power consumption capability of a power system under peak load and frequency regulation constraints, and belongs to the field of power system running and wind power grid-connection control. The method comprises the following steps of: in view of the peak load and frequency regulation constraints of the power system, optimizing the starting and stopping states of conventional machine sets, and identifying a three-dimensional controlled network for the wind power consumption of the power system the day before; according to the three-dimensional controlled network for the wind power consumption of the power system, judging whether predicted wind power can be completely consumed by the power system or not the day before; according to a judgment result indicating whether the predicted wind power can be completely consumed by the power system or not, determining an optimal control strategy for the predicted wind power consumption of the power system the day before; and controlling a wind generating set and the conventional machine sets the next day according to the optimal control strategy determined the day before. By the method, power system scheduling staff can be helped to quickly judge whether the predicted wind power can be completely consumed by the power system or not, the optimal control strategy for the predicted wind power consumption is provided, and practical significance and broad application prospect are presented.

Description

The lower electric power system of peak-frequency regulation constraint can dissolve identification and the control method of wind-powered electricity generation

Technical field

The invention belongs to power system operation and wind-electricity integration control field, particularly the lower electric power system of peak-frequency regulation constraint identification and the control method of wind-powered electricity generation of can dissolving.

Background technology

Since the eighties in last century, oil crisis, climate change, energy problem become international focus, and the clean energy resource take wind energy as representative is fast-developing, become important alternative energy source at a specified future date.Usually, wind-powered electricity generation is exerted oneself and is shown the distinguishing feature that is different from normal power supplies: 1) short time yardstick (hour level and hour level in) wind-powered electricity generation go out fluctuation be about installed capacity of wind-driven power ± 10% to ± 35%, went out fluctuation above ± 50% in 4～12 hours; The expansion of the district's scope that 2) covers along with wind energy turbine set, the complementarity between the wind energy turbine set unit strengthens, and wind energy turbine set integral body larger, that regional extent is disperseed goes out the fluctuation less; 3) the daily output characteristic of wind-powered electricity generation is opposite trend with part throttle characteristics more, i.e. anti-peak regulation characteristic.These characteristics are that the safe operation of electric power system has brought severe challenge with stable control, the method that science therefore is provided is to realize electric power system can dissolve identification and the control of wind-powered electricity generation, be main application of the present invention, it will deeply affect each function link such as operation, scheduling, control of electric power system.

The prerequisite of formulating the wind-electricity integration control strategy is can the dissolve identification of wind-powered electricity generation of electric power system, namely judges a few days ago (actual motion the previous day) prediction gained wind-powered electricity generation whether can fully be dissolved by electric power system (not abandoning wind).In the existing research, predict the gained wind-powered electricity generation for electric power system a few days ago, the basic skills of identification is: based on the classical scheduling model of electric power system (three public scheduling, economic dispatch, energy-saving distribution, low-carbon (LC) scheduling), whether checking exists feasible Unit Combination to exert oneself and can fully dissolve to guarantee a day wind-powered electricity generation.The method remains in deficiency:

1) start-up mode is determined, the reasonable conventional unit of closed portion can not excavate the electric power system ability of dissolving to greatest extent, the loss part wind-powered electricity generation of can dissolving;

2) ignore the frequency security constraint of electric power system, lack on the short time yardstick to the Security Checking of wind-electricity integration control strategy the frequency unstability risk that exists degree of depth peak regulation to bring;

3) identification result presents with a dimension indicator, can not characterize wind-powered electricity generation style characteristic (positive peak regulation, Heibei provincial opera peak, anti-peak regulation) and the mapping relations between the capacity of can dissolving, and can not control criterion for the spot dispatch operations staff provides;

The limiting behaviour that 4) can not provide electric power system can dissolve wind-powered electricity generation can not provide the visual electric power system wind-powered electricity generation " security domain " of dissolving for the spot dispatch operations staff.

5) calculation scale is large, operation time is long, can do nothing to help the spot dispatch operations staff and makes quick judgement.

And from the angle of the control strategy of wind-electricity integration, current institute Adopts measure only relates to the control of closing down of wind-powered electricity generation unit itself, does not cooperate the start-stop with conventional rack, is limiting electric power system dissolve controlled original paper and the controlled range of wind-powered electricity generation.Under the open state that conventional unit is determined, can't dissolve too much wind-powered electricity generation when exerting oneself when electric power system, the closed portion wind turbine consists of unique possible strategy, thereby has wasted the part wind power resources.

In sum, need one to overlap more electric power system and comprehensively can dissolve wind-powered electricity generation identification and control method, to take into account electric power system peak-frequency regulation ability, consider the wind-powered electricity generation power producing characteristics, the conventional unit start-up mode of optimal control is for the power system dispatching operations staff provides the instrument of Fast Identification with the control wind-powered electricity generation.

Summary of the invention

The objective of the invention is to overcome can the dissolve deficiency of wind-powered electricity generation identification and control method of existing electric power system, the identification and the control method that provide the lower electric power system of a kind of peak-frequency regulation constraint can dissolve wind-powered electricity generation, the present invention can help the power system dispatching operations staff a few days ago with regard to the dissolve three-dimensional controllable domain of wind-powered electricity generation of clear and definite electric power system, judge fast prediction gained wind-powered electricity generation and whether can fully be dissolved by electric power system, and the optimal control policy of the prediction gained wind-powered electricity generation of dissolving is provided.

The present invention proposes the lower electric power system of a kind of peak-frequency regulation constraint can dissolve identification and the control method of wind-powered electricity generation, it is characterized in that, comprise: 1) take into account electric power system peak-frequency regulation constraint, optimize conventional unit startup-shutdown state, at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago; 2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system; 3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago; 4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day;

1) takes into account the constraint of electric power system peak-frequency regulation; optimize conventional unit startup-shutdown state; (to reach the wind-powered electricity generation that belongs to arbitrarily this controllable domain is exerted oneself at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago; the purpose that electric power system all can fully be dissolved by rational control), specifically may further comprise the steps:

1-1) the dissolve three-dimensional controllable domain of wind-powered electricity generation of definition electric power system:

(I) the equivalent load rate η of the wind-powered electricity generation of dissolving, expression formula is as follows:

$\mathrm{\η}=\frac{P_{\mathrm{WIND}}^{\left(\mathrm{mean}\right)}}{P_{\mathrm{WIND}}^{\left(1\right)}}---\left(1\right)$

(II) peak interval of time of the wind-powered electricity generation of the dissolving poor δ that exerts oneself, expression formula is as follows:

$\mathrm{\δ}=\frac{P_{\mathrm{WIND}}^{\left(1\right)}-P_{\mathrm{WIND}}^{\left(2\right)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(2\right)$

(III) go out power rate γ the peak period of the wind-powered electricity generation of dissolving, expression formula is as follows:

$\mathrm{\γ}=\frac{P_{\mathrm{WIND}}^{\left(1\right)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(3\right)$

Wherein, t ⁽¹⁾Be power system load peak period; t ⁽²⁾Be the power system load low-valley interval;

For predict gained electric power system load peak period a few days ago;

Be the wind-powered electricity generation day part average output of being dissolved;

Wind-powered electricity generation is exerted oneself peak period in order to be dissolved;

The wind-powered electricity generation low-valley interval is exerted oneself in order to be dissolved;

1-2) according to the historical statistical data (comprising the historical statistical data that it is poor that wind-powered electricity generation equivalent load rate and peak interval of time are exerted oneself) of wind-powered electricity generation, preset dissolve the wherein value of two dimensions of the three-dimensional controllable domain of wind-powered electricity generation of electric power system, expression formula is as follows:

Ω＝{(η _h,δ _l)|η _h∈∏;δ _l∈Δ;h∈1L S ₁;l∈1L S ₂} (4)

In the formula (4),

∏ is the wind-powered electricity generation equivalent load rate of the being dissolved collection of setting, and expression formula is as follows:

∏={η _s|s∈1L S ₁}(5)

In the formula (5), η ₁Be wind-powered electricity generation equivalent load rate minimum in the historical statistics;

Be wind-powered electricity generation equivalent load rate maximum in the historical statistics; S ₁Number of elements for the wind-powered electricity generation equivalent load rate collection of dissolving set;

Δ is the wind-powered electricity generation peak interval of time of being dissolved of the setting difference set of exerting oneself, and expression formula is as follows:

Δ={δ _s|s∈1L S ₂}(6)

In the formula (6), δ ₁For wind-powered electricity generation peak interval of time minimum in the historical statistics exert oneself poor;

For wind-powered electricity generation peak interval of time maximum in the historical statistics exert oneself poor; S ₂Be the exert oneself number of elements of difference set of the wind-powered electricity generation peak interval of time of dissolving of setting;

1-3) according to the Ω that has set, selected identification variable specifically comprises:

The wind-powered electricity generation equivalent load rate η that dissolves that is setting _hWith the peak interval of time poor δ that exerts oneself _lCondition under, select to comprise that wind-powered electricity generation that dissolve peak period exerts oneself

With the conventional unit open state vector of electric power system D _{H, l}As the identification variable;

Wherein, D _{H, l}={ d _i| i=1L N}, if d _i=0, unit i keeps start; If d _i=1, then unit i is closed condition; N is conventional unit quantity;

1-4) according to the two-dimentional controllable domain of setting and selected identification variable, make up the dissolve Optimal Identification target function of wind-powered electricity generation controllable domain third dimension degree of electric power system, expression formula is as follows:

$\mathrm{\γ}({\mathrm{\η}}_h,{\mathrm{\δ}}_l)=\max \frac{P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\η}}_h,{\mathrm{\δ}}_l)}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(7\right)$

The implication of this target function (7) is: the wind-powered electricity generation equivalent load rate η that is setting _hWith the peak interval of time poor δ that exerts oneself _lCondition under, the maximum wind that electric power system can be dissolved peak period is exerted oneself;

1-5) based on the two-dimentional controllable domain of setting and selected identification variable, set up the related of operation states of electric power system variable (characterizing the parameter of power supply reliability, frequency modulation fail safe and peak regulation fail safe in the power system operation process) and identification variable, specifically comprise:

(I) mathematic expectaion (LOLE) of electric power system second order power failure hourage, expression formula is as follows:

$\mathrm{LOLE}(D_{h,l},{\mathrm{\η}}_h,P_{\mathrm{WIND}}^{\left(1\right)})=\underset{i}{\mathrm{\Σ}}\left\{\right[q_i\·\underset{k\≠i}{\mathrm{\Π}}(1-q_k)]\·t\left(P_{h,l}^{(\mathrm{loss}-i)}\right)\}---\left(8\right)$

$+\underset{i}{\mathrm{\Σ}}\underset{j}{\mathrm{\Σ}}\left\{\right[q_i\·q_j\·\underset{h\≠i,j}{\mathrm{\Π}}(1-q_h)]\·t\left(P_{h,l}^{(\mathrm{loss}-i,j)}\right)\}$

In the formula (8):

But be the average power supply capacity of electric power system when conventional unit i fault is only arranged, expression formula is as follows:

$P_{h,l}^{(\mathrm{loss}-i)}={\mathrm{\η}}_h{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\η}}_h,{\mathrm{\δ}}_l)+\underset{u=1}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }---\left(9\right)$

But be the average power supply capacity of electric power system when conventional unit i and j fault are only arranged, expression formula is as follows:

$P_{h,l}^{(\mathrm{loss}-i,j)}={\mathrm{\η}}_h{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\η}}_h,{\mathrm{\δ}}_l)+\underset{u=1}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }-P_{\mathrm{GEN},j}^{\max }---\left(10\right)$

Maximum output for conventional unit u;

For the conventional unit d that closes _k=1, have

q _iFailure rate for conventional unit i;

For the conventional unit d that closes _k=1, have

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{h,l}^{(\mathrm{loss}-i)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\Σ}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{h,l}^{(\mathrm{loss}-i)})+1]---\left(11\right)$

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{h,l}^{(\mathrm{loss}-i,j)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\Σ}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{h,l}^{(\mathrm{loss}-i,j)})+1]---\left(12\right)$

P _{LOAD, t}Be the load of electric power system period t, t ∈ 1L T;

(II) frequency fluctuation that causes of electric power system short time peak period yardstick wind-powered electricity generation fluctuation

Expression formula is as follows:

${\mathrm{\Δf}}_{h,l}^{\left(1\right)}=\frac{f_0{\mathrm{\ΔP}}^{\left(1\right)}}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(13\right)$

In the formula (13):

The power supply vacancy that the fluctuation of wind-powered electricity generation short time peak period yardstick causes is Δ P ⁽¹⁾, expression formula is as follows:

${\mathrm{\ΔP}}^{\left(1\right)}={\mathrm{\λ}}^{\left(1\right)}{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\η}}_h,{\mathrm{\δ}}_l)-P_R^{\left(1\right)}---\left(14\right)$

K is power system load-frequency effect adjustment factor; f ₀Be initial power system frequency; λ ⁽¹⁾For wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio of peak period;

Be electric power system peak period reserve capacity, expression formula is as follows:

$P_R^{\left(1\right)}=P_{\mathrm{GEN}}^{(\max -1)}-E\·{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(15\right)$

Can exert oneself for electric power system maximum peak period, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\max -1)}=P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\η}}_h,{\mathrm{\δ}}_l)+[E-D-M]\·{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\·{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(16\right)$

E is that dimension 1 * N and element are 1 vector; M is the vector of dimension 1 * N, is 1 with the corresponding position of the firm outputs such as heat supply, interconnector among the M, and other positions are 0;

Conventional unit goes out force vector peak period when being incorporated into the power networks without wind-powered electricity generation;

(III) frequency fluctuation that causes of electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation

Expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{{\mathrm{<figure-callout\; id="0"\; label="f"\; filenames="HDA0000127997390000011.png,HDA0000127997390000012.png"\; state="\{\{state\}\}">f</figure-callout>}}_0{\mathrm{\&Delta;P}}^{\left(2\right)}}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(2\right)}}}---\left(17\right)$

In the formula (17):

The power supply vacancy that the fluctuation of low-valley interval wind-powered electricity generation short time yardstick causes is Δ P ⁽²⁾, expression formula is as follows:

${\mathrm{\&Delta;P}}^{\left(2\right)}={\mathrm{\&lambda;}}^{\left(2\right)}{P}_{\mathrm{WIND}}^{\left(2\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(2\right)}---\left(18\right)$

λ ⁽²⁾For wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio of low-valley interval;

Be electric power system low-valley interval load;

Be low-valley interval electric power system reserve capacity, expression formula is as follows:

$P_R^{\left(2\right)}=P_{\mathrm{GEN}}^{(\max -2)}-E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(19\right)$

Can exert oneself for electric power system low-valley interval maximum, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\max -2)}=P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}+[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(20\right)$

Conventional unit low-valley interval goes out force vector when being incorporated into the power networks without wind-powered electricity generation;

(IV) the electric power system wind-powered electricity generation capacity of can dissolving peak period is under the fully dark peak regulation state

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(\max -1)}=E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T-P_{\mathrm{GEN}}^{(\min -1)}---\left(21\right)$

In the formula (21),

Be conventional machine group minimum output peak period of electric power system, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\min -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\min }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(22\right)$

Be conventional unit minimum output vector;

(V) the electric power system low-valley interval wind-powered electricity generation capacity of can dissolving is under the fully dark peak regulation state

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(\max -2)}=E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T-P_{\mathrm{GEN}}^{(\min -2)}---\left(23\right)$

In the formula (23),

Be the minimum output of the conventional machine group low-valley interval of electric power system, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\min -2)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\min }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(24\right)$

Be conventional unit minimum output vector;

(VI) electric power system maximum power supply capacity peak period is

Expression formula is as follows:

$P_{\mathrm{SYS}}^{(\max -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T+P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)---\left(25\right)$

(VII) the maximum power supply capacity of electric power system low-valley interval is

Expression formula is as follows:

$P_{\mathrm{SYS}}^{(\max -2)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T+P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}---\left(26\right)$

1-6) related based on operation states of electric power system variable and identification variable determined between the feasible region of the control range of operation states of electric power system variable and identification variable, and then establishes the constraints to controllable domain third dimension degree Optimal Identification:

(I) power system power supply reliability constraint, expression formula is as follows:

$\mathrm{LOLE}(D_{h,l},{\mathrm{\&eta;}}_h,P_{\mathrm{WIND}}^{\left(1\right)})-{\mathrm{LOLE}}_{\mathrm{BASE}}\&le;0---\left(27\right)$

The implication of this formula (27) is: the mathematic expectaion of electric power system second order power failure hourage during wind-electricity integration The mathematic expectaion LOLE that is not higher than electric power system benchmark second order power failure hourage _BASE

In the formula (27), LOLE _BASEBe the mathematic expectaion of electric power system benchmark second order power failure hourage, expression formula is as follows:

${\mathrm{LOLE}}_{\mathrm{BASE}}=\underset{i}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;\underset{k\&NotEqual;i}{\mathrm{\&Pi;}}(1-q_k)]\&CenterDot;t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}\right)\}+\underset{i}{\mathrm{\&Sigma;}}\underset{j}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;q_j\&CenterDot;\underset{h\&NotEqual;i,j}{\mathrm{\&Pi;}}(1-q_h)]\&CenterDot;t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}\right)\}---\left(28\right)$

In the formula (28):

For not having under the wind-electricity integration condition, but when conventional unit i fault is only arranged the electric power system power supply capacity, expression formula is as follows:

$P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}=\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }---\left(29\right)$

For not having under the wind-electricity integration condition, but when conventional unit i and j fault are only arranged the electric power system power supply capacity, expression formula is as follows:

$P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}=\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }-P_{\mathrm{GEN},j}^{\max }---\left(30\right)$

q _iStatistics failure rate for conventional unit i;

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{\mathrm{BASE}}^{(\mathrm{loss}-i)})+1]---\left(31\right)$

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)})+1---\left(32\right)$

(II) electric power system frequency modulation peak period constraint, expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(1\right)}=\frac{f_0[{\mathrm{\&lambda;}}^{\left(1\right)}{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(1\right)}]}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(1\right)}}}\&le;{\mathrm{\&Delta;f}}_{\max }---\left(33\right)$

In the formula (32), Δ f _MaxBe the receptible peak frequency fluctuation of electric power system;

The frequency fluctuation that causes for electric power system short time peak period yardstick wind-powered electricity generation fluctuation;

(III) electric power system low-valley interval frequency modulation constraint, expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{f_0[{\mathrm{\&lambda;}}^{\left(2\right)}(P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}})-P_R^{\left(2\right)}]}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(2\right)}}}\&le;{\mathrm{\&Delta;f}}_{\max }---\left(34\right)$

In the formula (34), Δ f _MaxBe the receptible peak frequency fluctuation of electric power system;

The frequency fluctuation that causes for electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation;

(IV) electric power system peak regulation peak period constraint, expression formula is as follows:

$P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)\&le;P_{\mathrm{WIND}}^{(\max -1)}---\left(35\right)$

Wherein,

Be the wind-powered electricity generation capacity of can dissolving peak period of electric power system under the fully dark peak regulation state;

(V) electric power system low-valley interval peak regulation constraint, expression formula is as follows:

$P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\&le;P_{\mathrm{WIND}}^{(\max -2)}---\left(36\right)$

Wherein,

Be the wind-powered electricity generation capacity of can dissolving of electric power system low-valley interval under the fully dark peak regulation state;

(VI) electric power system power supply capacity peak period constraint, expression formula is as follows:

$P_{\mathrm{LOAD},t^{\left(1\right)}}\&le;P_{\mathrm{SYS}}^{(\max -1)}---\left(37\right)$

Wherein,

Be electric power system maximum power supply capacity peak period;

(VII) electric power system low-valley interval power supply capacity constraint, expression formula is as follows:

$P_{\mathrm{LOAD},t^{\left(2\right)}}\&le;P_{\mathrm{SYS}}^{(\max -2)}---\left(38\right)$

Wherein,

Be the maximum power supply capacity of electric power system low-valley interval;

(VIII) conventional unit startup-shutdown constraint, expression formula is as follows:

1-7) constraints of the Optimal Identification target function of composite type (7) foundation and formula (27)～(39) foundation forms the Optimal Identification model, finds the solution this model, the equivalent load rate η that is setting _hWith the peak interval of time poor δ that exerts oneself _lUnder the condition, obtain the dissolve controllable domain third dimension degree γ (η of wind-powered electricity generation of electric power system _h, δ _l);

1-8) be cycled to repeat step 1-3), 1-4), 1-5), 1-6), 1-7), the two-dimentional controllable domain Ω based on setting finishes the Optimal Identification to controllable domain third dimension degree, obtains the dissolve three-dimensional controllable domain Z of wind-powered electricity generation of electric power system, expression formula is as follows:

Z={(η _h,δ _l,γ _h,l)|(h,l)∈Ω} (40)

2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge that a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system, specifically comprises:

2-1) the fundamental characteristics of prediction gained wind-powered electricity generation specifically comprises:

(I) prediction gained wind-powered electricity generation is exerted oneself peak period Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-1)}=P_{\mathrm{WIND},t^{\left(1\right)}}^{\left(0\right)}---\left(41\right)$

(II) prediction gained wind-powered electricity generation low-valley interval is exerted oneself

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-2)}=P_{\mathrm{WIND},t^{\left(2\right)}}^{\left(0\right)}---\left(42\right)$

(III) prediction gained wind-powered electricity generation wind-powered electricity generation average output

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-\mathrm{mean})}=\frac{1}{T}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{P}_{\mathrm{WIND},t}^{\left(0\right)}---\left(43\right)$

Wherein,

Be exerting oneself of prediction gained wind-powered electricity generation day part;

2-2) based on the fundamental characteristics of prediction gained wind-powered electricity generation, three dimension indicators of gained wind-powered electricity generation: η is predicted in measuring and calculating ⁽⁰⁾, δ ⁽⁰⁾, γ ⁽⁰⁾, specifically comprise:

(I) wind-powered electricity generation equivalent load rate η ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&eta;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-\mathrm{mean})}}{P_{\mathrm{WIND}}^{(0-1)}}---\left(44\right)$

(II) the peak interval of time wind-powered electricity generation poor δ that exerts oneself ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&delta;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-1)}-P_{\mathrm{WIND}}^{(0-2)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(45\right)$

(III) peak period, wind-powered electricity generation went out power rate γ ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&gamma;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-1)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(46\right)$

2-3) based on three dimension indicators of prediction gained wind-powered electricity generation and the electric power system of the identification gained three-dimensional controllable domain of wind-powered electricity generation of dissolving, judge and predict that whether the gained wind-powered electricity generation can fully be dissolved by electric power system, specifically comprises:

(I) optimum position prediction gained wind-powered electricity generation (namely in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system, navigates to point with its Euclidean distance minimum with prediction gained wind-powered electricity generation in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system

), optimum position model tormulation formula is as follows:

min (η ⁽⁰⁾-η _h) ²+(δ ⁽⁰⁾-δ _l) ² (47)

s.t. (η _h,δ _l)∈Ω

Finding the solution this optimum position model (47) obtains

(II) mapped location in the three-dimensional controllable domain Z of wind-powered electricity generation is dissolved in electric power system, expression formula is as follows:

$z^{*}=({\mathrm{\&eta;}}_{h^{*}},{\mathrm{\&delta;}}_{l^{*}},{\mathrm{\&gamma;}}_{h^{*},l^{*}})---\left(48\right)$

(III) according to the result of mapped location among the three-dimensional controllable domain Z of wind-powered electricity generation of dissolving in electric power system, to judge whether electric power system can fully dissolve to predict the gained wind-powered electricity generation, the criterion expression formula is as follows:

Prediction gained wind-powered electricity generation is exerted oneself if electric power system can be dissolved, then target r ₀=1; Otherwise, r then ₀=0;

3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago, specifically comprises:

3-1) determine prediction gained wind-powered electricity generation day part is abandoned wind control

Optimized model (close minimum total wind-powered electricity generation and exert oneself to guarantee the whole day T period), expression formula is as follows:

$\min \underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}$

$s.t.r\left({\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}\right)=1---\left(50\right)$

$0\&le;{\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}\&le;P_{\mathrm{WIND},t}^{\left(0\right)}$

In the Optimized model (50), For the wind-powered electricity generation that day part is closed is exerted oneself; By step 2) described method obtains For closing

The target amount whether rear electric power system can fully dissolve to the residue wind-powered electricity generation;

3-2) determine the control of closing down of the conventional unit of electric power system, specifically comprise:

(I) establish control variables, specifically comprise:

The conventional unit open state vector of N platform D ₀Wherein, D ₀={ d _{0, i}| i=1L N}, if d _{0, i}=0, conventional unit i start; If d _{0, i}=1, then conventional unit i closes;

(II) target function (so that the total capacity of the unit of closing is maximum) of the conventional unit optimal control of structure, expression formula is as follows:

$\mathrm{<figure-callout\; id="0"\; label="max\; D"\; filenames="HDA0000127997390000011.png,HDA0000127997390000012.png"\; state="\{\{state\}\}">max</figure-callout>}{D}_{}{}_0\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T---\left(51\right)$

(III) integrated objective function (51) forms conventional unit optimal control Optimized model with the constraints of formula (27)～(39) foundation, finds the solution this Optimized model, obtains the controlled quentity controlled variable D of conventional unit ₀

4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day, specifically comprise:

4-1) control output of wind electric field specifically comprises:

In the actual motion of next day, at period t, dispatch command is assigned to wind energy turbine set in the power system dispatching center, closes

Effective output;

4-2) control conventional unit output, specifically comprise:

In the actual motion of next day, dispatch command is assigned by the conventional unit of mind-set power plant in the power system dispatching: if d _{0, i}=1, then assign out code to conventional unit; If d _{0, i}=0, then assign start-up command to conventional unit.

Technical characterstic of the present invention and beneficial effect:

The present invention has jumped out existing electric power system wind-powered electricity generation identification and the control method constraint in flow scheme design and theoretical method aspect of can dissolving, the lower electric power system of a cover peak-frequency regulation constraint can dissolve identification and the control method of wind-powered electricity generation have been set up, complete electric power system peak regulation, the fm capacity taken into account, take into full account the wind-powered electricity generation power producing characteristics, the conventional unit start-up mode of scientific optimization is for the power system dispatching operations staff provides the instrument of a cover Fast Identification with the control wind-powered electricity generation.The present invention can help the power system dispatching operations staff a few days ago with regard to the dissolve three-dimensional controllable domain of wind-powered electricity generation of clear and definite electric power system, whether judge fast prediction gained wind-powered electricity generation can fully be dissolved by electric power system, and the optimal control policy of the prediction gained wind-powered electricity generation of dissolving is provided, each function links such as the operation of electric power system, scheduling, control are had important practical significance and good application prospect.

Description of drawings

Fig. 1 is the dissolve three-dimensional controllable domain of wind-powered electricity generation of embodiments of the invention electric power system;

Fig. 2 is that embodiments of the invention prediction gained wind-powered electricity generation is exerted oneself;

Fig. 3 is that the wind-powered electricity generation that embodiments of the invention wind energy turbine set day part is closed is exerted oneself;

Embodiment

Below in conjunction with drawings and Examples, be elaborated as follows to can dissolve identification and the control method of wind-powered electricity generation of the lower electric power system of peak frequency modulation constraint.The invention provides the lower electric power system of a kind of peak-frequency regulation constraint can dissolve identification and the control method of wind-powered electricity generation, it is characterized in that, comprise: 1) take into account electric power system peak-frequency regulation constraint, optimize conventional unit startup-shutdown state, at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago; 2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system; 3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago; 4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day;

1) takes into account the constraint of electric power system peak-frequency regulation; optimize conventional unit startup-shutdown state; (to reach the wind-powered electricity generation that belongs to arbitrarily this controllable domain is exerted oneself at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago; the purpose that electric power system all can fully be dissolved by rational control), specifically may further comprise the steps:

1-1) the dissolve three-dimensional controllable domain of wind-powered electricity generation of definition electric power system:

(I) the equivalent load rate η of the wind-powered electricity generation of dissolving, expression formula is as follows:

$\mathrm{\&eta;}=\frac{P_{\mathrm{WIND}}^{\left(\mathrm{mean}\right)}}{P_{\mathrm{WIND}}^{\left(1\right)}}---\left(1\right)$

(II) peak interval of time of the wind-powered electricity generation of the dissolving poor δ that exerts oneself, expression formula is as follows:

$\mathrm{\&delta;}=\frac{P_{\mathrm{WIND}}^{\left(1\right)}-P_{\mathrm{WIND}}^{\left(2\right)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(2\right)$

(III) go out power rate γ the peak period of the wind-powered electricity generation of dissolving, expression formula is as follows:

$\mathrm{\&gamma;}=\frac{P_{\mathrm{WIND}}^{\left(1\right)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(3\right)$

Wherein, t ⁽¹⁾Be power system load peak period; t ⁽²⁾Be the power system load low-valley interval;

For predict gained electric power system load peak period a few days ago;

Be the wind-powered electricity generation day part average output of being dissolved;

Wind-powered electricity generation is exerted oneself peak period in order to be dissolved;

The wind-powered electricity generation low-valley interval is exerted oneself in order to be dissolved;

1-2) according to the historical statistical data (comprising the historical statistical data that it is poor that wind-powered electricity generation equivalent load rate and peak interval of time are exerted oneself) of wind-powered electricity generation, preset dissolve the wherein value of two dimensions of the three-dimensional controllable domain of wind-powered electricity generation of electric power system, expression formula is as follows:

Ω＝{(η _h,δ _l)|η _h∈∏;δl∈Δ;h∈1L S ₁;l∈1L S ₂} (4)

In the formula (4),

∏ is the wind-powered electricity generation equivalent load rate of the being dissolved collection of setting, and expression formula is as follows:

∏={η _s|s∈1L S ₁} (5)

In the formula (5), η ₁Be wind-powered electricity generation equivalent load rate minimum in the historical statistics;

Be wind-powered electricity generation equivalent load rate maximum in the historical statistics; S ₁Number of elements for the wind-powered electricity generation equivalent load rate collection of dissolving set;

Δ is the wind-powered electricity generation peak interval of time of being dissolved of the setting difference set of exerting oneself, and expression formula is as follows:

Δ={δ _s|s∈1L S ₂} (6)

In the formula (6), δ ₁For wind-powered electricity generation peak interval of time minimum in the historical statistics exert oneself poor;

For wind-powered electricity generation peak interval of time maximum in the historical statistics exert oneself poor; S ₂Be the exert oneself number of elements of difference set of the wind-powered electricity generation peak interval of time of dissolving of setting;

1-3) according to the Ω that has set, selected identification variable specifically comprises:

The wind-powered electricity generation equivalent load rate η that dissolves that is setting _hWith the peak interval of time poor δ that exerts oneself _lCondition under, select to comprise that wind-powered electricity generation that dissolve peak period exerts oneself

With the conventional unit open state vector of electric power system D _{H, l}As the identification variable;

Wherein, D _{H, l}={ d _i| i=1L N}, if d _i=0, unit i keeps start; If d _i=1, then unit i is closed condition; N is conventional unit quantity;

1-4) according to the two-dimentional controllable domain of setting and selected identification variable, make up the dissolve Optimal Identification target function of wind-powered electricity generation controllable domain third dimension degree of electric power system, expression formula is as follows:

$\mathrm{\&gamma;}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)=\max \frac{P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(7\right)$

The implication of this target function (7) is: the wind-powered electricity generation equivalent load rate η that is setting _hWith the peak interval of time poor δ that exerts oneself _lCondition under, the maximum wind that electric power system can be dissolved peak period is exerted oneself;

1-5) based on the two-dimentional controllable domain of setting and selected identification variable, set up the related of operation states of electric power system variable (characterizing the parameter of power supply reliability, frequency modulation fail safe and peak regulation fail safe in the power system operation process) and identification variable, specifically comprise:

(I) mathematic expectaion (LOLE) of electric power system second order power failure hourage, expression formula is as follows:

$\mathrm{LOLE}(D_{h,l},{\mathrm{\&eta;}}_h,P_{\mathrm{WIND}}^{\left(1\right)})=\underset{i}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;\underset{k\&NotEqual;i}{\mathrm{\&Pi;}}(1-q_k)]\&CenterDot;t\left(P_{h,l}^{(\mathrm{loss}-i)}\right)\}---\left(8\right)$

$+\underset{i}{\mathrm{\&Sigma;}}\underset{j}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;q_j\&CenterDot;\underset{h\&NotEqual;i,j}{\mathrm{\&Pi;}}(1-q_h)]\&CenterDot;t\left(P_{h,l}^{(\mathrm{loss}-i,j)}\right)\}$

In the formula (8):

But be the average power supply capacity of electric power system when conventional unit i fault is only arranged, expression formula is as follows:

$P_{h,l}^{(\mathrm{loss}-i)}={\mathrm{\&eta;}}_h{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)+\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }---\left(9\right)$

But be the average power supply capacity of electric power system when conventional unit i and j fault are only arranged, expression formula is as follows:

$P_{h,l}^{(\mathrm{loss}-i,j)}={\mathrm{\&eta;}}_h{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)+\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }-P_{\mathrm{GEN},j}^{\max }---\left(10\right)$

Maximum output for conventional unit u;

For the conventional unit d that closes _k=1, have

q _iFailure rate for conventional unit i;

For the conventional unit d that closes _k=1, have

For power system load greater than The accumulation hourage, expression formula is as follows:

$t\left(P_{h,l}^{(\mathrm{loss}-i)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{h,l}^{(\mathrm{loss}-i)})+1]---\left(11\right)$

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{h,l}^{(\mathrm{loss}-i,j)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{h,l}^{(\mathrm{loss}-i,j)})+1]---\left(12\right)$

P _{LOAD, t}Be the load of electric power system period t, t ∈ 1L T;

(II) frequency fluctuation that causes of electric power system short time peak period yardstick wind-powered electricity generation fluctuation Expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(1\right)}=\frac{{\mathrm{<figure-callout\; id="0"\; label="f"\; filenames="HDA0000127997390000011.png,HDA0000127997390000012.png"\; state="\{\{state\}\}">f</figure-callout>}}_0{\mathrm{\&Delta;P}}^{\left(1\right)}}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(13\right)$

In the formula (13):

The power supply vacancy that the fluctuation of wind-powered electricity generation short time peak period yardstick causes is Δ P ⁽¹⁾, expression formula is as follows:

${\mathrm{\&Delta;P}}^{\left(1\right)}={\mathrm{\&lambda;}}^{\left(1\right)}{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(1\right)}---\left(14\right)$

K is power system load-frequency effect adjustment factor; f ₀Be initial power system frequency; λ ⁽¹⁾For wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio of peak period;

Be electric power system peak period reserve capacity, expression formula is as follows:

$P_R^{\left(1\right)}=P_{\mathrm{GEN}}^{(\max -1)}-E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(15\right)$

Can exert oneself for electric power system maximum peak period, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\max -1)}=P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)+[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(16\right)$

E is that dimension 1 * N and element are 1 vector; M is the vector of dimension 1 * N, is 1 with the corresponding position of the firm outputs such as heat supply, interconnector among the M, and other positions are 0;

Conventional unit goes out force vector peak period when being incorporated into the power networks without wind-powered electricity generation;

(III) frequency fluctuation that causes of electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation

Expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{{\mathrm{<figure-callout\; id="0"\; label="f"\; filenames="HDA0000127997390000011.png,HDA0000127997390000012.png"\; state="\{\{state\}\}">f</figure-callout>}}_0{\mathrm{\&Delta;P}}^{\left(2\right)}}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(2\right)}}}---\left(17\right)$

In the formula (17):

The power supply vacancy that the fluctuation of low-valley interval wind-powered electricity generation short time yardstick causes is Δ P ⁽²⁾, expression formula is as follows:

${\mathrm{\&Delta;P}}^{\left(2\right)}={\mathrm{\&lambda;}}^{\left(2\right)}{P}_{\mathrm{WIND}}^{\left(2\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(2\right)}---\left(18\right)$

λ ⁽²⁾For wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio of low-valley interval;

Be electric power system low-valley interval load;

Be low-valley interval electric power system reserve capacity, expression formula is as follows:

$P_R^{\left(2\right)}=P_{\mathrm{GEN}}^{(\max -2)}-E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(19\right)$

Can exert oneself for electric power system low-valley interval maximum, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\max -2)}=P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}+[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(20\right)$

Conventional unit low-valley interval goes out force vector when being incorporated into the power networks without wind-powered electricity generation;

(IV) the electric power system wind-powered electricity generation capacity of can dissolving peak period is under the fully dark peak regulation state

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(\max -1)}=E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T-P_{\mathrm{GEN}}^{(\min -1)}---\left(21\right)$

In the formula (21), Be conventional machine group minimum output peak period of electric power system, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\min -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\min }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(22\right)$

Be conventional unit minimum output vector;

(V) the electric power system low-valley interval wind-powered electricity generation capacity of can dissolving is under the fully dark peak regulation state

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(\max -2)}=E\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T-P_{\mathrm{GEN}}^{(\min -2)}---\left(23\right)$

In the formula (23),

Be the minimum output of the conventional machine group low-valley interval of electric power system, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\min -2)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\min }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(24\right)$

Be conventional unit minimum output vector;

(VI) electric power system maximum power supply capacity peak period is Expression formula is as follows:

$P_{\mathrm{SYS}}^{(\max -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\right]}^T+P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)---\left(25\right)$

(VII) the maximum power supply capacity of electric power system low-valley interval is

Expression formula is as follows:

$P_{\mathrm{SYS}}^{(\max -2)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{\mathrm{GEN},t^{\left(2\right)}}^{\mathrm{typical}}\right]}^T+P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}---\left(26\right)$

1-6) related based on operation states of electric power system variable and identification variable determined between the feasible region of the control range of operation states of electric power system variable and identification variable, and then establishes the constraints to controllable domain third dimension degree Optimal Identification:

(I) power system power supply reliability constraint, expression formula is as follows:

$\mathrm{LOLE}(D_{h,l},{\mathrm{\&eta;}}_h,P_{\mathrm{WIND}}^{\left(1\right)})-{\mathrm{LOLE}}_{\mathrm{BASE}}\&le;0---\left(27\right)$

The implication of this formula (27) is: the mathematic expectaion of electric power system second order power failure hourage during wind-electricity integration The mathematic expectaion LOLE that is not higher than electric power system benchmark second order power failure hourage _BASE

In the formula (27), LOLE _BASEBe the mathematic expectaion of electric power system benchmark second order power failure hourage, expression formula is as follows:

${\mathrm{LOLE}}_{\mathrm{BASE}}=\underset{i}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;\underset{k\&NotEqual;i}{\mathrm{\&Pi;}}(1-q_k)]\&CenterDot;t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}\right)\}+\underset{i}{\mathrm{\&Sigma;}}\underset{j}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;q_j\&CenterDot;\underset{h\&NotEqual;i,j}{\mathrm{\&Pi;}}(1-q_h)]\&CenterDot;t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}\right)\}---\left(28\right)$

In the formula (28):

For not having under the wind-electricity integration condition, but when conventional unit i fault is only arranged the electric power system power supply capacity, expression formula is as follows:

$P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}=\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }---\left(29\right)$

For not having under the wind-electricity integration condition, but when conventional unit i and j fault are only arranged the electric power system power supply capacity, expression formula is as follows:

$P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}=\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }-P_{\mathrm{GEN},j}^{\max }---\left(30\right)$

q _iStatistics failure rate for conventional unit i;

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{\mathrm{BASE}}^{(\mathrm{loss}-i)})+1]---\left(31\right)$

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)})+1---\left(32\right)$

(II) electric power system frequency modulation peak period constraint, expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(1\right)}=\frac{f_0[{\mathrm{\&lambda;}}^{\left(1\right)}{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(1\right)}]}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(1\right)}}}\&le;{\mathrm{\&Delta;f}}_{\max }---\left(33\right)$

In the formula (32), Δ f _MaxBe the receptible peak frequency fluctuation of electric power system;

The frequency fluctuation that causes for electric power system short time peak period yardstick wind-powered electricity generation fluctuation;

(III) electric power system low-valley interval frequency modulation constraint, expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{f_0[{\mathrm{\&lambda;}}^{\left(2\right)}(P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}})-P_R^{\left(2\right)}]}{{\mathrm{KP}}_{\mathrm{LOAD},t^{\left(2\right)}}}\&le;{\mathrm{\&Delta;f}}_{\max }---\left(34\right)$

In the formula (34), Δ f _MaxBe the receptible peak frequency fluctuation of electric power system; The frequency fluctuation that causes for electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation;

(IV) electric power system peak regulation peak period constraint, expression formula is as follows:

$P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)\&le;P_{\mathrm{WIND}}^{(\max -1)}---\left(35\right)$

Wherein,

Be the wind-powered electricity generation capacity of can dissolving peak period of electric power system under the fully dark peak regulation state;

(V) electric power system low-valley interval peak regulation constraint, expression formula is as follows:

$P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{\mathrm{GEN},t^{\left(1\right)}}^{\mathrm{typical}}\&le;P_{\mathrm{WIND}}^{(\max -2)}---\left(36\right)$

Wherein,

Be the wind-powered electricity generation capacity of can dissolving of electric power system low-valley interval under the fully dark peak regulation state;

(VI) electric power system power supply capacity peak period constraint, expression formula is as follows:

$P_{\mathrm{LOAD},t^{\left(1\right)}}\&le;P_{\mathrm{SYS}}^{(\max -1)}---\left(37\right)$

Wherein,

Be electric power system maximum power supply capacity peak period;

(VII) electric power system low-valley interval power supply capacity constraint, expression formula is as follows:

$P_{\mathrm{LOAD},t^{\left(2\right)}}\&le;P_{\mathrm{SYS}}^{(\max -2)}---\left(38\right)$

Wherein,

Be the maximum power supply capacity of electric power system low-valley interval;

(VIII) conventional unit startup-shutdown constraint, expression formula is as follows:

1-7) constraints of the Optimal Identification target function of composite type (7) foundation and formula (27)～(39) foundation forms the Optimal Identification model, finds the solution this model, the equivalent load rate η that is setting _hWith the peak interval of time poor δ that exerts oneself _lUnder the condition, obtain the dissolve controllable domain third dimension degree γ (η of wind-powered electricity generation of electric power system _h, δ _l);

1-8) be cycled to repeat step 1-3), 1-4), 1-5), 1-6), 1-7), the two-dimentional controllable domain Ω based on setting finishes the Optimal Identification to controllable domain third dimension degree, obtains the dissolve three-dimensional controllable domain Z of wind-powered electricity generation of electric power system, expression formula is as follows:

Z={(η _h,δ _l,γ _h,l)|(h,l)∈Ω} (40)

2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge that a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system, specifically comprises:

2-1) the fundamental characteristics of prediction gained wind-powered electricity generation specifically comprises:

(I) prediction gained wind-powered electricity generation is exerted oneself peak period

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-1)}=P_{\mathrm{WIND},t^{\left(1\right)}}^{\left(0\right)}---\left(41\right)$

(II) prediction gained wind-powered electricity generation low-valley interval is exerted oneself

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-2)}=P_{\mathrm{WIND},t^{\left(2\right)}}^{\left(0\right)}---\left(42\right)$

(III) prediction gained wind-powered electricity generation wind-powered electricity generation average output

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-\mathrm{mean})}=\frac{1}{T}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{P}_{\mathrm{WIND},t}^{\left(0\right)}---\left(43\right)$

Wherein,

Be exerting oneself of prediction gained wind-powered electricity generation day part;

2-2) based on the fundamental characteristics of prediction gained wind-powered electricity generation, three dimension indicators of gained wind-powered electricity generation: η is predicted in measuring and calculating ⁽⁰⁾, δ ⁽⁰⁾, γ ⁽⁰⁾, specifically comprise:

(I) wind-powered electricity generation equivalent load rate η ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&eta;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-\mathrm{mean})}}{P_{\mathrm{WIND}}^{(0-1)}}---\left(44\right)$

(II) the peak interval of time wind-powered electricity generation poor δ that exerts oneself ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&delta;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-1)}-P_{\mathrm{WIND}}^{(0-2)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(45\right)$

(III) peak period, wind-powered electricity generation went out power rate γ ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&gamma;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-1)}}{P_{\mathrm{LOAD},t^{\left(1\right)}}}---\left(46\right)$

2-3) based on three dimension indicators of prediction gained wind-powered electricity generation and the electric power system of the identification gained three-dimensional controllable domain of wind-powered electricity generation of dissolving, judge and predict that whether the gained wind-powered electricity generation can fully be dissolved by electric power system, specifically comprises:

(I) optimum position prediction gained wind-powered electricity generation (in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system, navigates to point with its Euclidean distance minimum with prediction gained wind-powered electricity generation in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system

), optimum position model tormulation formula is as follows:

min (η ⁽⁰⁾-η _h) ²+(δ ⁽⁰⁾-δ _l) ² (47)

s.t. (η _h,δ _l)∈Ω

Finding the solution this optimum position model (47) obtains

(II) mapped location in the three-dimensional controllable domain Z of wind-powered electricity generation is dissolved in electric power system, expression formula is as follows:

$z^{*}=({\mathrm{\&eta;}}_{h^{*}},{\mathrm{\&delta;}}_{l^{*}},{\mathrm{\&gamma;}}_{h^{*},l^{*}})---\left(48\right)$

(III) according to the result of mapped location among the three-dimensional controllable domain Z of wind-powered electricity generation of dissolving in electric power system, to judge whether electric power system can fully dissolve to predict the gained wind-powered electricity generation, the criterion expression formula is as follows:

Prediction gained wind-powered electricity generation is exerted oneself if electric power system can be dissolved, then target r ₀=1; Otherwise, r then ₀=0;

3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago, specifically comprises:

3-1) determine prediction gained wind-powered electricity generation day part is abandoned wind control Optimized model (close minimum total wind-powered electricity generation and exert oneself to guarantee the whole day T period), expression formula is as follows:

$\min \underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}$

$s.t.r\left({\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}\right)=1---\left(50\right)$

$0\&le;{\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}\&le;P_{\mathrm{WIND},t}^{\left(0\right)}$

In the Optimized model (50),

For the wind-powered electricity generation that day part is closed is exerted oneself; By step 2) described method obtains

For closing

The target amount whether rear electric power system can fully dissolve to the residue wind-powered electricity generation;

3-2) determine the control of closing down of the conventional unit of electric power system, specifically comprise:

(I) establish control variables, specifically comprise:

The conventional unit open state vector of N platform D ₀Wherein, D ₀={ d _{0, i}| i=1L N}, if d _{0, i}=0, conventional unit i start; If d _{0, i}=1, then conventional unit i closes;

(II) target function (so that the total capacity of the unit of closing is maximum) of the conventional unit optimal control of structure, expression formula is as follows:

$\mathrm{<figure-callout\; id="0"\; label="max\; D"\; filenames="HDA0000127997390000011.png,HDA0000127997390000012.png"\; state="\{\{state\}\}">max</figure-callout>}{D}_{}{}_0\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T---\left(51\right)$

(III) integrated objective function (51) forms conventional unit optimal control Optimized model with the constraints of formula (27)～(39) foundation, finds the solution this Optimized model, obtains the controlled quentity controlled variable D of conventional unit ₀

4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day, specifically comprise:

4-1) control output of wind electric field specifically comprises:

In the actual motion of next day, at period t, dispatch command is assigned to wind energy turbine set in the power system dispatching center, closes

Effective output;

4-2) control conventional unit output, specifically comprise:

In the actual motion of next day, dispatch command is assigned by the conventional unit of mind-set power plant in the power system dispatching: if d _{0, i}=1, then assign out code to conventional unit; If d _{0, i}=0, then assign start-up command to conventional unit.

Embodiment:

Set forth the lower electric power system of peak-frequency regulation constraint proposed by the invention can dissolve identification and the control method of wind-powered electricity generation as an example of certain provincial area example.This embodiment is based on the power system planning software kit of power planning decision support electric power system GOPT5.0(by Electric Motor Engineering and Applied Electronic Technology Department, Qinghua's exploitation, set up the unified planning Optimized model of power industry sustainable development, seek science, economy, continuable electric power development scheme), under the condition of not considering wind-electricity integration, carry out production simulation, obtain certain operation day load prediction data, wind-powered electricity generation prediction data and conventional unit start plan a few days ago.

The setup parameter value according to the historical statistical data (comprising the historical statistical data that it is poor that wind-powered electricity generation equivalent load rate and peak interval of time are exerted oneself) of this area's wind-powered electricity generation, presets dissolve the wherein value of two dimensions of the three-dimensional controllable domain of wind-powered electricity generation of electric power system, specifically comprises:

Set the set of wind-powered electricity generation equivalent load rate collection ∏:

η＝{0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0}；

Set the peak interval of time wind-powered electricity generation difference set of exerting oneself:

Δ＝{-0.2 -0.18 -0.16 -0.14 -0.12 -0.10 0.7 0.8 0.9 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 0.2}；

Determine according to this area's electric power system conventional operation data: round-the-clock hop count T=24; Conventional unit quantity N=134; Power system load-frequency effect adjustment factor K=5.0; Initial power system frequency f ₀=50Hz; Wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio lambda of peak period ⁽¹⁾=10%; Wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio lambda of peak period ⁽²⁾=20%; The receptible peak frequency fluctuation Δ f of electric power system institute _Max=0.2Hz; Electric power system t peak period ⁽¹⁾=12; Electric power system load peak period

Electric power system low-valley interval t ⁽²⁾=7; Electric power system low-valley interval load Conventional unit failure rate, maximum output, minimum output, whether the AGC unit, whether the heat supply unit, whether can to close unit information as shown in table 1.

1) takes into account electric power system peak-frequency regulation constraint, optimize conventional unit startup-shutdown state, at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago;

Adopt method provided by the invention, Identification of Power System a few days ago dissolve wind-powered electricity generation controllable domain as shown in Figure 1.That the three-dimensional coordinate of Fig. 1 is respectively wind-powered electricity generation equivalent load rate, peak interval of time is exerted oneself is poor, peak period the wind-powered electricity generation ability of dissolving, wherein wind-powered electricity generation equivalent load rate is the result who gets denary logarithm.Can see, wind-powered electricity generation equivalent load rate less than 1(logarithm value less than 0) and the peak interval of time poor timing that is of exerting oneself, electric power system has the higher wind-powered electricity generation ability of dissolving, its physical significance is: wind-powered electricity generation power producing characteristics and power system load characteristic similarity degree are higher, and the ability of dissolving of electric power system is stronger;

2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system;

Shown in Figure 2 is that a few days ago prediction wind-powered electricity generation that GOPT5.0 provides is exerted oneself, and abscissa is the period, and ordinate is that wind-powered electricity generation is exerted oneself, and for example the 1st period predicted that it was 6100 megawatts that wind-powered electricity generation is exerted oneself.Adopt method provided by the invention, this wind-powered electricity generation is exerted oneself and is determined and can't be fully dissolved by electric power system;

3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago;

Adopt method provided by the invention, the wind-powered electricity generation that the wind energy turbine set day part is closed is exerted oneself as shown in Figure 3, and abscissa is the period, and ordinate is to abandon air quantity, and for example to abandon wind be 1210 megawatts the 1st period; The conventional rack control strategy is shown in " control command " in the table 1, and for example unit 1 start-up command is that 1(closes), unit 11 start-up command are that 0(opens);

4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day;

In the next day actual motion, " control command " is that 1 conventional unit is assigned shutdown command in the his-and-hers watches 1, closes such unit; " control command " is that 0 conventional unit is assigned start-up command in the his-and-hers watches 1, opens such unit; Assign to wind energy turbine set at day part by Fig. 3 and to abandon the wind instruction.

Table 1

Above-described specific embodiment is only for explanation realization effect of the present invention, not in order to limit the present invention.All in method proposed by the invention basic ideas and framework within modification, conversion and the improvement of any unsubstantiality of doing, all should be included within protection scope of the present invention.

Claims (1)

1. the lower electric power system of peak-frequency regulation constraint identification and the control method of wind-powered electricity generation of can dissolving, it is characterized in that, comprise: 1) take into account electric power system peak-frequency regulation constraint, optimize conventional unit startup-shutdown state, at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago; 2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system; 3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago; 4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day;

1) take into account electric power system peak-frequency regulation constraint, optimize conventional unit startup-shutdown state, at the dissolve three-dimensional controllable domain of wind-powered electricity generation of Identification of Power System a few days ago, specifically may further comprise the steps:

1-1) the dissolve three-dimensional controllable domain of wind-powered electricity generation of definition electric power system:

(I) the equivalent load rate η of the wind-powered electricity generation of dissolving, expression formula is as follows:

$\mathrm{\&eta;}=\frac{P_{\mathrm{WIND}}^{\left(\mathrm{mean}\right)}}{P_{\mathrm{WIND}}^{\left(1\right)}}---\left(1\right)$

(II) peak interval of time of the wind-powered electricity generation of the dissolving poor δ that exerts oneself, expression formula is as follows:

$\mathrm{\&delta;}=\frac{P_{\mathrm{WIND}}^{\left(1\right)}-P_{\mathrm{WIND}}^{\left(2\right)}}{P_{{\mathrm{LOAD},t}^{\left(1\right)}}}---\left(2\right)$

(III) go out power rate γ the peak period of the wind-powered electricity generation of dissolving, expression formula is as follows:

$\mathrm{\&gamma;}=\frac{P_{\mathrm{WIND}}^{\left(1\right)}}{P_{{\mathrm{LOAD},t}^{\left(1\right)}}}---\left(3\right)$

Wherein, t ⁽¹⁾Be power system load peak period; t ⁽²⁾Be the power system load low-valley interval;

For predict gained electric power system load peak period a few days ago;

Be the wind-powered electricity generation day part average output of being dissolved;

Wind-powered electricity generation is exerted oneself peak period in order to be dissolved; The wind-powered electricity generation low-valley interval is exerted oneself in order to be dissolved;

1-2) according to the historical statistical data of wind-powered electricity generation, preset dissolve the wherein value of two dimensions of the three-dimensional controllable domain of wind-powered electricity generation of electric power system, expression formula is as follows:

Ω＝{(η _h,δ _l)|η _h∈∏;δ _l∈Δ;h∈1…S ₁;l∈1…S ₂}(4)

In the formula (4),

∏ is the wind-powered electricity generation equivalent load rate of the being dissolved collection of setting, and expression formula is as follows:

∏={η _s|s∈1…S ₁}(5)

In the formula (5), η ₁Be wind-powered electricity generation equivalent load rate minimum in the historical statistics; Be wind-powered electricity generation equivalent load rate maximum in the historical statistics; S ₁Number of elements for the wind-powered electricity generation equivalent load rate collection of dissolving set;

Δ is the wind-powered electricity generation peak interval of time of being dissolved of the setting difference set of exerting oneself, and expression formula is as follows:

Δ={δ _s|s∈1…S ₂}(6)

In the formula (6), δ ₁For wind-powered electricity generation peak interval of time minimum in the historical statistics exert oneself poor;

For wind-powered electricity generation peak interval of time maximum in the historical statistics exert oneself poor; S ₂Be the exert oneself number of elements of difference set of the wind-powered electricity generation peak interval of time of dissolving of setting;

1-3) according to the Ω that has set, selected identification variable specifically comprises:

At the wind-powered electricity generation equivalent load rate η h that dissolves of institute that sets and the peak interval of time poor δ that exerts oneself _lCondition under, select to comprise that wind-powered electricity generation that dissolve peak period exerts oneself

With the conventional unit open state vector of electric power system D _{H, l}As the identification variable;

Wherein, D _{H, l}={ d _i| i=1 ... if N} is d _i=0, unit i keeps start; If d _i=1, then unit i is closed condition; N is conventional unit quantity;

1-4) according to the two-dimentional controllable domain of setting and selected identification variable, make up the dissolve Optimal Identification target function of wind-powered electricity generation controllable domain third dimension degree of electric power system, expression formula is as follows:

$\mathrm{\&gamma;}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)=\max \frac{P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)}{P_{{\mathrm{LOAD},t}^{\left(1\right)}}}---\left(7\right)$

The implication of this target function (7) is: the wind-powered electricity generation equivalent load rate η that is setting _hWith the peak interval of time poor δ that exerts oneself _lCondition under, the maximum wind that electric power system can be dissolved peak period is exerted oneself;

1-5) based on the two-dimentional controllable domain of setting and selected identification variable, set up the related of operation states of electric power system variable and identification variable, specifically comprise:

(I) mathematic expectaion of electric power system second order power failure hourage, expression formula is as follows:

$\mathrm{LOLE}(D_{h,l},{\mathrm{\&eta;}}_h,P_{\mathrm{WIND}}^{\left(1\right)})=\underset{i}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;\underset{k\&NotEqual;i}{\mathrm{\&Pi;}}(1-q_k)]\&CenterDot;t\left(P_{h,l}^{(\mathrm{loss}-i)}\right)\}$ $\left(8\right)$

$+\underset{i}{\mathrm{\&Sigma;}}\underset{j}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;q_j\&CenterDot;\underset{h\&NotEqual;i,j}{\mathrm{\&Pi;}}(1-q_h)]\&CenterDot;t\left(P_{h,l}^{(\mathrm{loss}-i,j)}\right)\}$

In the formula (8):

But be the average power supply capacity of electric power system when conventional unit i fault is only arranged, expression formula is as follows:

$P_{h,l}^{(\mathrm{loss}-i)}={\mathrm{\&eta;}}_h{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)+\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }---\left(9\right)$

But be the average power supply capacity of electric power system when conventional unit i and j fault are only arranged, expression formula is as follows:

$P_{h,l}^{(\mathrm{loss}-i,j)}={\mathrm{\&eta;}}_h{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)+\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }-P_{\mathrm{GEN},j}^{\max }---\left(10\right)$

Maximum output for conventional unit u;

For the conventional unit d that closes _k=1, have

q _iStatistics failure rate for conventional unit i;

For the conventional unit d that closes _k=1, have

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{h,l}^{(\mathrm{loss}-i)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{h,l}^{(\mathrm{loss}-i)})+1]---\left(11\right)$

For power system load greater than The accumulation hourage, expression formula is as follows:

$t\left(P_{h,l}^{(\mathrm{loss}-i,j)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{h,l}^{(\mathrm{loss}-i,j)})+1]---\left(12\right)$

P _{LOAD, t}Be the load of electric power system period t, t ∈ 1 ... T;

(II) frequency fluctuation that causes of electric power system short time peak period yardstick wind-powered electricity generation fluctuation Expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(1\right)}=\frac{f_0{\mathrm{\&Delta;P}}^{\left(1\right)}}{{\mathrm{KP}}_{{\mathrm{LOAD},t}^{\left(1\right)}}}---\left(13\right)$

In the formula (13):

The power supply vacancy that the fluctuation of wind-powered electricity generation short time peak period yardstick causes is Δ P ⁽¹⁾, expression formula is as follows:

${\mathrm{\&Delta;P}}^{\left(1\right)}={\mathrm{\&lambda;}}^{\left(1\right)}{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(1\right)}---\left(14\right)$

K is power system load-frequency effect adjustment factor; f ₀Be initial power system frequency; λ ⁽¹⁾For wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio of peak period;

Be electric power system peak period reserve capacity, expression formula is as follows:

$P_R^{\left(1\right)}=P_{\mathrm{GEN}}^{(\max -1)}-E\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(15\right)$

Can exert oneself for electric power system maximum peak period, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\max -1)}=P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)+[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(16\right)$

E is that dimension 1 * N and element are 1 vector; M is the vector of dimension 1 * N, is 1 with heat supply, the corresponding position of interconnector firm output among the M, and other positions are 0;

Conventional unit goes out force vector peak period when being incorporated into the power networks without wind-powered electricity generation;

(III) frequency fluctuation that causes of electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation

Expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{f_0{\mathrm{\&Delta;P}}^{\left(2\right)}}{{\mathrm{KP}}_{{\mathrm{LOAD},t}^{\left(2\right)}}}---\left(17\right)$

In the formula (17):

The power supply vacancy that the fluctuation of low-valley interval wind-powered electricity generation short time yardstick causes is Δ P ⁽²⁾, expression formula is as follows:

${\mathrm{\&Delta;P}}^{\left(2\right)}={\mathrm{\&lambda;}}^{\left(2\right)}{P}_{\mathrm{WIND}}^{\left(2\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(2\right)}---\left(18\right)$

λ ⁽²⁾For wind-powered electricity generation is exerted oneself in the maximum fluctuation ratio of low-valley interval;

Be electric power system low-valley interval load; Be low-valley interval electric power system reserve capacity, expression formula is as follows:

$P_R^{\left(2\right)}=P_{\mathrm{GEN}}^{(\max -2)}-E\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(19\right)$

Can exert oneself for electric power system low-valley interval maximum, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\max -2)}=P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}+[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(20\right)$

Conventional unit low-valley interval goes out force vector when being incorporated into the power networks without wind-powered electricity generation;

(IV) the electric power system wind-powered electricity generation capacity of can dissolving peak period is under the fully dark peak regulation state

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(\max -1)}=\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}\right]}^T-P_{\mathrm{GEN}}^{(\min -1)}---\left(21\right)$

In the formula (21),

Be conventional machine group minimum output peak period of electric power system, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\min -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\min }\right]}^T+M\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}\right]}^T---\left(22\right)$

Be conventional unit minimum output vector;

(V) the electric power system low-valley interval wind-powered electricity generation capacity of can dissolving is under the fully dark peak regulation state Expression formula is as follows:

$P_{\mathrm{WIND}}^{(\max -2)}=E\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(2\right)}}^{\mathrm{typical}}\right]}^T-P_{\mathrm{GEN}}^{(\min -2)}---\left(23\right)$

In the formula (23),

Be the minimum output of the conventional machine group low-valley interval of electric power system, expression formula is as follows:

$P_{\mathrm{GEN}}^{(\min -2)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\min }\right]}^T+M\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(2\right)}}^{\mathrm{typical}}\right]}^T---\left(24\right)$

Be conventional unit minimum output vector;

(VI) electric power system maximum power supply capacity peak period is

Expression formula is as follows:

$P_{\mathrm{SYS}}^{(\max -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}\right]}^T+P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)---\left(25\right)$

(VII) the maximum power supply capacity of electric power system low-valley interval is

Expression formula is as follows:

$P_{\mathrm{SYS}}^{(\max -2)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{{\mathrm{GEN},t}^{\left(2\right)}}^{\mathrm{typical}}\right]}^T+P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}---\left(26\right)$

1-6) related based on operation states of electric power system variable and identification variable determined between the feasible region of the control range of operation states of electric power system variable and identification variable, and then establishes the constraints to controllable domain third dimension degree Optimal Identification:

(I) power system power supply reliability constraint, expression formula is as follows:

$\mathrm{LOLE}(D_{h,l},{\mathrm{\&eta;}}_h,P_{\mathrm{WIND}}^{\left(1\right)})-{\mathrm{LOLE}}_{\mathrm{BASE}}\&le;0---\left(27\right)$

The implication of this formula (27) is: the mathematic expectaion of electric power system second order power failure hourage during wind-electricity integration

The mathematic expectaion LOLE that is not higher than electric power system benchmark second order power failure hourage _BASE

In the formula (27), LOLE _BASEBe the mathematic expectaion of electric power system benchmark second order power failure hourage, expression formula is as follows:

${\mathrm{LOLE}}_{\mathrm{BASE}}=\underset{i}{\mathrm{\&Sigma;}}\left\{\right[q_i\&CenterDot;\underset{k\&NotEqual;i}{\mathrm{\&Pi;}}(1-q_k)]\&CenterDot;t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}\right)\}+\underset{i}{\mathrm{\&Sigma;}}\underset{j}{\mathrm{\&Sigma;}}\left\{\right[q_i\text{\&CenterDot;}{q}_j\&CenterDot;\underset{h\&NotEqual;i,j}{\mathrm{\&Sigma;}}(1-q_h)]\&CenterDot;t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}\right)\}---\left(28\right)$

In the formula (28):

For not having under the wind-electricity integration condition, but when conventional unit i fault is only arranged the electric power system power supply capacity, expression formula is as follows:

$P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}=\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }---\left(29\right)$

For not having under the wind-electricity integration condition, but when conventional unit i and j fault are only arranged the electric power system power supply capacity, expression formula is as follows:

$P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}=\underset{u=1}{\overset{N}{\mathrm{\&Sigma;}}}{P}_{\mathrm{GEN},u}^{\max }-P_{\mathrm{GEN},i}^{\max }-P_{\mathrm{GEN},j}^{\max }---\left(30\right)$

q _iStatistics failure rate for conventional unit i;

For power system load greater than The accumulation hourage, expression formula is as follows:

$t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{\mathrm{BASE}}^{(\mathrm{loss}-i)})+1]---\left(31\right)$

For power system load greater than

The accumulation hourage, expression formula is as follows:

$t\left(P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)}\right)=\frac{1}{2}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}[\mathrm{sign}(P_{\mathrm{LOAD},t}-P_{\mathrm{BASE}}^{(\mathrm{loss}-i,j)})+1]---\left(32\right)$

(II) electric power system frequency modulation peak period constraint, expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(1\right)}=\frac{f_0[{\mathrm{\&lambda;}}^{\left(1\right)}{P}_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-P_R^{\left(1\right)}}{{\mathrm{KP}}_{{\mathrm{LOAD},t}^{\left(1\right)}}}\&le;{\mathrm{\&Delta;f}}_{\max }---\left(33\right)$

In the formula (32), Δ f _MaxBe the receptible peak frequency fluctuation of electric power system;

The frequency fluctuation that causes for electric power system short time peak period yardstick wind-powered electricity generation fluctuation;

(III) electric power system low-valley interval frequency modulation constraint, expression formula is as follows:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{f_0[{\mathrm{\&lambda;}}^{\left(2\right)}(P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_l{P}_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}})-P_R^{\left(2\right)}]}{{\mathrm{KP}}_{{\mathrm{LOAD},t}^{\left(2\right)}}}\&le;{\mathrm{\&Delta;f}}_{\max }---\left(34\right)$

In the formula (34), Δ f _MaxBe the receptible peak frequency fluctuation of electric power system;

The frequency fluctuation that causes for electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation;

(IV) electric power system peak regulation peak period constraint, expression formula is as follows:

$P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)\&le;P_{\mathrm{WIND}}^{(\max -1)}---\left(35\right)$

Wherein, Be the wind-powered electricity generation capacity of can dissolving peak period of electric power system under the fully dark peak regulation state;

(V) electric power system low-valley interval peak regulation constraint, expression formula is as follows:

$P_{\mathrm{WIND}}^{\left(1\right)}({\mathrm{\&eta;}}_h,{\mathrm{\&delta;}}_l)-{\mathrm{\&delta;}}_{\text{l}}{P}_{{\mathrm{GEN},t}^{\left(1\right)}}^{\mathrm{typical}}\&le;P_{\mathrm{WIND}}^{(\max -2)}---\left(36\right)$

Wherein,

Be the wind-powered electricity generation capacity of can dissolving of electric power system low-valley interval under the fully dark peak regulation state;

(VI) electric power system power supply capacity peak period constraint, expression formula is as follows:

$P_{{\mathrm{LOAD},t}^{\left(1\right)}}\&le;P_{\mathrm{SYS}}^{(\max -1)}---\left(37\right)$

Wherein,

Be electric power system maximum power supply capacity peak period;

(VII) electric power system low-valley interval power supply capacity constraint, expression formula is as follows:

$P_{{\mathrm{LOAD},t}^{\left(2\right)}}\&le;P_{\mathrm{SYS}}^{(\max -2)}---\left(38\right)$

Wherein,

Be the maximum power supply capacity of electric power system low-valley interval;

(VIII) conventional unit startup-shutdown constraint, expression formula is as follows:

1-7) constraints set up of the Optimal Identification target function set up of composite type (7) and formula (27) ~ (39) forms the Optimal Identification model, finds the solution this model, at the equivalent load rate η h of setting and the peak interval of time poor δ that exerts oneself _lUnder the condition, obtain the dissolve controllable domain third dimension degree γ (η of wind-powered electricity generation of electric power system _h, δ _l);

1-8) be cycled to repeat step 1-3), 1-4), 1-5), 1-6), 1-7), the two-dimentional controllable domain Ω based on setting finishes the Optimal Identification to controllable domain third dimension degree, obtains the dissolve three-dimensional controllable domain Z of wind-powered electricity generation of electric power system, expression formula is as follows:

Z={(η _h，δ _l,γ _h,l)|(h,l)∈Ω}(40)

2) according to the dissolve three-dimensional controllable domain of wind-powered electricity generation of electric power system, judge that a few days ago whether prediction gained wind-powered electricity generation can fully be dissolved by electric power system, specifically comprises:

2-1) the fundamental characteristics of prediction gained wind-powered electricity generation specifically comprises:

(I) prediction gained wind-powered electricity generation is exerted oneself peak period

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-1)}=P_{{\mathrm{WIND},t}^{\left(1\right)}}^{\left(0\right)}---\left(41\right)$

(II) prediction gained wind-powered electricity generation low-valley interval is exerted oneself

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-2)}=P_{{\mathrm{WIND},t}^{\left(2\right)}}^{\left(0\right)}---\left(42\right)$

(III) prediction gained wind-powered electricity generation wind-powered electricity generation average output

Expression formula is as follows:

$P_{\mathrm{WIND}}^{(0-\mathrm{mean})}=\frac{1}{T}\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{P}_{\mathrm{WIND},t}^{\left(0\right)}---\left(43\right)$

Wherein,

Be exerting oneself of prediction gained wind-powered electricity generation day part;

2-2) based on the fundamental characteristics of prediction gained wind-powered electricity generation, three dimension indicators of gained wind-powered electricity generation: η is predicted in measuring and calculating ⁽⁰⁾, δ ⁽⁰⁾, γ ⁽⁰⁾, specifically comprise:

(I) wind-powered electricity generation equivalent load rate η ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&eta;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-\mathrm{mean})}}{P_{\mathrm{WIND}}^{(0-1)}}---\left(44\right)$

(II) the peak interval of time wind-powered electricity generation poor δ that exerts oneself ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&delta;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-1)}-P_{\mathrm{WIND}}^{(0-2)}}{P_{{\mathrm{LOAD},t}^{\left(1\right)}}}---\left(45\right)$

(III) peak period, wind-powered electricity generation went out power rate γ ⁽⁰⁾, expression formula is as follows:

${\mathrm{\&gamma;}}^{\left(0\right)}=\frac{P_{\mathrm{WIND}}^{(0-1)}}{P_{{\mathrm{LOAD},t}^{\left(1\right)}}}---\left(46\right)$

2-3) based on three dimension indicators of prediction gained wind-powered electricity generation and the electric power system of the identification gained three-dimensional controllable domain of wind-powered electricity generation of dissolving, judge and predict that whether the gained wind-powered electricity generation can fully be dissolved by electric power system, specifically comprises:

(I) optimum position prediction gained wind-powered electricity generation in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system, optimum position model tormulation formula is as follows:

min(η ⁽⁰⁾-η _h) ²+(δ ⁽⁰⁾-δ _l) ²(47)

s.t.(η _h,δ _l)∈Ω

Finding the solution this optimum position model (47) obtains

(II) mapped location in the three-dimensional controllable domain Z of wind-powered electricity generation is dissolved in electric power system, expression formula is as follows:

$z^{*}=({\mathrm{\&eta;}}_{h^{*}},{\mathrm{\&delta;}}_{l^{*}},{\mathrm{\&gamma;}}_{h^{*},l^{*}})---\left(48\right)$

(III) according to the result of mapped location among the three-dimensional controllable domain Z of wind-powered electricity generation of dissolving in electric power system, to judge whether electric power system can fully dissolve to predict the gained wind-powered electricity generation, the criterion expression formula is as follows:

Prediction gained wind-powered electricity generation is exerted oneself if electric power system can be dissolved, then target r ₀=1; Otherwise, r then ₀=0;

3) judged result that whether can fully be dissolved by electric power system according to prediction gained wind-powered electricity generation is being determined the dissolve optimal control policy of prediction gained wind-powered electricity generation of electric power system a few days ago, specifically comprises:

3-1) determine prediction gained wind-powered electricity generation day part is abandoned wind control

Optimized model, expression formula is as follows:

$\min \underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}$

$s.t.r\left({\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}\right)=1---\left(50\right)$

$0\&le;{\mathrm{\&Delta;P}}_{\mathrm{WIND},t}^{\left(0\right)}\&le;P_{\mathrm{WIND},t}^{\left(0\right)}$

In the Optimized model (50),

For the wind-powered electricity generation that day part is closed is exerted oneself; By step 2) described method obtains

For closing

The target amount whether rear electric power system can fully dissolve to the residue wind-powered electricity generation;

3-2) determine the control of closing down of the conventional unit of electric power system, specifically comprise:

(I) establish control variables, specifically comprise:

The conventional unit open state vector of N platform D ₀Wherein, D ₀={ d _{0, i}| i=1 ... if N} is d _{0, i}=0, conventional unit i start; If d _{0, i}=1, then conventional unit i closes;

(II) target function of the conventional unit optimal control of structure, so that the total capacity of the unit of closing is maximum, expression formula is as follows:

$\max D_0\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T---\left(51\right)$

(III) integrated objective function (51) forms conventional unit optimal control Optimized model with the constraints of formula (27) ~ (39) foundation, finds the solution this Optimized model, obtains the controlled quentity controlled variable D of conventional unit ₀

4) according to determined optimal control policy a few days ago, control wind-powered electricity generation unit and conventional unit in next day, specifically comprise:

4-1) control output of wind electric field specifically comprises:

In the actual motion of next day, at period t, dispatch command is assigned to wind energy turbine set in the power system dispatching center, closes

Effective output;

4-2) control conventional unit output, specifically comprise:

In the actual motion of next day, dispatch command is assigned by the conventional unit of mind-set power plant in the power system dispatching: if d _{0, i}=1, then assign out code to conventional unit; If d _{0, i}=0, then assign start-up command to conventional unit.

Images