CN102496962A - 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

Peak-frequency regulation constraint is can the dissolve identification and the control method of wind-powered electricity generation of electric power system down

Technical field

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

Background technology

Since the eighties in last century, oil crisis, climate change, energy problem become international focus, are that the clean energy resource of representative is fast-developing with the wind energy, 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 conventional power supply: 1) short time yardstick (hour level and hour level in) wind-powered electricity generation go out fluctuation be about the wind-powered electricity generation installed capacity ± 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 it is less relatively that wind energy turbine set integral body larger, that regional extent is disperseed goes out fluctuation; 3) the daily output characteristic of wind-powered electricity generation and part throttle characteristics are opposite trend more, promptly anti-peak regulation characteristic.These characteristics are that the safe operation of electric power system has brought severe challenge with stable control; Therefore provide scientific methods to realize electric power system can the dissolve identification and the control of wind-powered electricity generation; Be main application of the present invention, it will deeply influence 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, judges promptly that a few days ago (actual motion the previous day) forecasting institute gets wind-powered electricity generation whether can fully be dissolved by electric power system (not abandoning wind).In the existing research; To electric power system a few days ago forecasting institute get wind-powered electricity generation; 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.This method remains in deficiency:

1) start mode confirms, 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 to go up the safety of wind-electricity integration control strategy is checked the frequency unstability risk that exists degree of depth peak regulation to bring at short time yardstick (in 15 minutes);

3) identification result appears 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) calculating scale is big, 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, the current measure of taking 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 the dissolve controlled original paper and the controlled range of wind-powered electricity generation.Under the open state that conventional unit is confirmed, can't dissolve too much wind-powered electricity generation when exerting oneself when electric power system, closed portion wind-powered electricity generation machine consists of unique possible strategy, thereby has wasted part wind-powered electricity generation resource.

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 of optimal control mode is for the power system dispatching operations staff provides the instrument of quick 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 power system; The identification and the control method that provide the following 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 forecasting institute fast and get wind-powered electricity generation and whether can fully be dissolved, and provide the forecasting institute of dissolving to get the optimal control policy of wind-powered electricity generation by electric power system.

The present invention proposes a kind of peak-frequency regulation constraint can the dissolve identification and the control method of wind-powered electricity generation of electric power system down; 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 electric 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 that a few days ago forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system; 3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation; 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; (wind-powered electricity generation that belongs to this controllable domain is arbitrarily exerted oneself at the dissolve three-dimensional controllable domain of wind-powered electricity generation of identification electric power system a few days ago to reach; The purpose that electric power system all can be able to fully dissolve through 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 following:

$\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 following:

$\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 following:

$\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 forecasting institute a few days ago gets electric power system load peak period;

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, preestablish 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 following:

Ω＝{(η _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 following:

∏＝{η _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 following:

Δ＝{δ _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 setting of dissolving;

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=1 ... If N} is d _i=0, unit i keeps start; If d _i=1, then unit i is a closed condition; N is conventional unit quantity;

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

$\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-powered electricity generation that electric power system can be dissolved peak period is exerted oneself;

1-5) based on 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 following:

$\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 following:

$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 following:

$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)$

is the EIAJ of 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

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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)$

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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 ∈ 1 ... T;

(II) frequency fluctuation

expression formula that causes of electric power system short time peak period yardstick (in 15 minutes) wind-powered electricity generation fluctuation is following:

$\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 (in 15 minutes) causes is Δ P ⁽¹⁾, expression formula is following:

$\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 following:

$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, and expression formula is following:

$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_{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 pairing position of firm outputs such as heat supply, interconnector among the M, and other positions are 0;

goes out force vector peak period for the no wind-powered electricity generation time routine unit that is incorporated into the power networks;

(III) frequency fluctuation

expression formula that causes of electric power system low-valley interval short time yardstick (in 15 minutes) wind-powered electricity generation fluctuation is following:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{{\mathrm{<figure-callout\; id="0"\; label="f"\; filenames="HDA0000127997390000021.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 (in 15 minutes) causes is Δ P ⁽²⁾, expression formula is following:

${\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 following:

$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, and expression formula is following:

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

Conventional unit low-valley interval went out force vector when was incorporated into the power networks for no wind-powered electricity generation;

(IV) the electric power system wind-powered electricity generation capacity of can dissolving peak period is following for

expression formula under the fully dark peak regulation state:

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

In the formula (21);

is conventional machine group minimum output peak period of electric power system, and expression formula is following:

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

is conventional unit minimum output vector;

(V) the electric power system low-valley interval wind-powered electricity generation capacity of can dissolving is following for

expression formula under the fully dark peak regulation state:

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

In the formula (23);

is the minimum output of the conventional machine group low-valley interval of electric power system, and expression formula is following:

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

is conventional unit minimum output vector;

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

expression formula is following:

$P_{\mathrm{SYS}}^{(\max -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{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 that expression formula is following:

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

1-6) related based on operation states of electric power system variable and identification variable, confirm between the feasible region of control range and identification variable of operation states of electric power system variable, and then establish constraints controllable domain third dimension degree optimized recognition:

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

$\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 following:

${\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):

be not for having under the wind-electricity integration condition; But electric power system power supply capacity when conventional unit i fault is only arranged, expression formula is following:

$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)$

be not for having under the wind-electricity integration condition; But electric power system power supply capacity when conventional unit i and j fault are only arranged, expression formula is following:

$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;

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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)$

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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 following:

${\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 institute;

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 following:

${\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}_{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 institute;

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 following:

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

Wherein,

is 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 following:

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

Wherein, is 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 following:

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

Wherein,

is electric power system maximum power supply capacity peak period;

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

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

Wherein,

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

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

1-7) constraints of the optimized recognition target function of composite type (7) foundation and formula (27)～(39) foundation forms the optimized recognition 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) cycle repeats step 1-3), 1-4), 1-5), 1-6), 1-7), based on the two-dimentional controllable domain Ω that sets, accomplish optimized recognition to controllable domain third dimension degree, obtain the dissolve three-dimensional controllable domain Z of wind-powered electricity generation of electric power system, expression formula is following:

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 forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system, specifically comprise:

2-1) forecasting institute gets the fundamental characteristics of wind-powered electricity generation, specifically comprises:

(I) to get wind-powered electricity generation

expression formula of exerting oneself peak period following for forecasting institute:

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

(II) to get wind-powered electricity generation low-valley interval

expression formula of exerting oneself following for forecasting institute:

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

(III) to get wind-powered electricity generation wind-powered electricity generation average output

expression formula following for forecasting institute:

$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,

gets exerting oneself of wind-powered electricity generation day part for forecasting institute;

2-2) get the fundamental characteristics of wind-powered electricity generation based on forecasting institute, the measuring and calculating forecasting institute gets three dimension indicators of wind-powered electricity generation: η ⁽⁰⁾, δ ⁽⁰⁾, γ ⁽⁰⁾, specifically comprise:

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

${\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 following:

${\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 following:

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

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

(I), electric power system optimizes location prediction gained wind-powered electricity generation in dissolving the two-dimentional controllable domain Ω of wind-powered electricity generation (promptly in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system; Forecasting institute is got wind-powered electricity generation to navigate to and the minimum point

of its Euclidean distance), it is following to optimize the location model expression formula:

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

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

Find the solution this optimization location model (47) and obtain

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

$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, judge that the electric power system forecasting institute of whether can fully dissolving gets wind-powered electricity generation, the criterion expression formula is following:

Forecasting institute gets wind-powered electricity generation and exerts oneself if electric power system can be dissolved, then target r ₀=1; Otherwise, r then ₀=0;

3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation, specifically comprises:

3-1) confirm that forecasting institute is got the wind-powered electricity generation day part abandons the Optimization Model of wind control

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

$\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)}$

Optimization model (50) in,

for all time off wind power output; according to step 2) the method

is off

after the remaining wind power systems are subject to the full amount of consumptive;

3-2) confirm 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 (so that the total capacity of the unit of closing is maximum) of the conventional unit optimal control of structure, expression formula is following:

$\mathrm{<figure-callout\; id="0"\; label="max\; D"\; filenames="HDA0000127997390000021.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 Optimization Model with the constraints of formula (27)～(39) foundation, finds the solution this Optimization 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, specifically comprise in next day:

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 the effective output of

;

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; Set up a cover peak-frequency regulation constraint can dissolve identification and the control method of wind-powered electricity generation of electric power system down; 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 of scientific optimization mode is for the power system dispatching operations staff provides the instrument of the quick identification of a cover 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 forecasting institute fast gets wind-powered electricity generation and can fully be dissolved by electric power system; And the optimal control policy that provides the forecasting institute of dissolving to get wind-powered electricity generation, each function links such as the operation of electric power system, scheduling, control are had important practical significance and good prospects for application.

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 exerts oneself for the embodiments of the invention forecasting institute gets wind-powered electricity generation;

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

Embodiment

Below in conjunction with accompanying drawing and embodiment, to peak frequency modulation constraint down can the dissolve identification and the control method of wind-powered electricity generation of electric power system be elaborated as follows.The invention provides a kind of peak-frequency regulation constraint can the dissolve identification and the control method of wind-powered electricity generation of electric power system down; 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 electric 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 that a few days ago forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system; 3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation; 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; (wind-powered electricity generation that belongs to this controllable domain is arbitrarily exerted oneself at the dissolve three-dimensional controllable domain of wind-powered electricity generation of identification electric power system a few days ago to reach; The purpose that electric power system all can be able to fully dissolve through 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 following:

$\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 following:

$\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 following:

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

Wherein, t (1) is power system load peak period; T (2) is the power system load low-valley interval; gets electric power system load peak period for forecasting institute a few days ago; is 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, preestablish 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 following:

Ω＝{(η _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 following:

∏＝{η _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 following:

Δ＝{δ ₂|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 setting of dissolving;

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=1 ... If N} is d _i=0, unit i keeps start; If d _i=1, then unit i is a closed condition; N is conventional unit quantity;

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

$\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-powered electricity generation that electric power system can be dissolved peak period is exerted oneself;

1-5) based on 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 following:

$\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 following:

$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 following:

$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)$

is the EIAJ of 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

is the accumulation hourage of power system load greater than , and expression formula is following:

$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)$

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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)$

PLOAD, t are the load of electric power system period t, t ∈ 1 ... T;

(II) frequency fluctuation

expression formula that causes of electric power system short time peak period yardstick (in 15 minutes) wind-powered electricity generation fluctuation is following:

${\mathrm{\&Delta;f}}_{h,l}^{\left(1\right)}=\frac{{\mathrm{<figure-callout\; id="0"\; label="f"\; filenames="HDA0000127997390000021.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 (in 15 minutes) causes is Δ P ⁽¹⁾, expression formula is following:

${\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 following:

$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, and expression formula is following:

$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_{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 pairing position of firm outputs such as heat supply, interconnector among the M, and other positions are 0;

goes out force vector peak period for the no wind-powered electricity generation time routine unit that is incorporated into the power networks;

(III) frequency fluctuation

expression formula that causes of electric power system low-valley interval short time yardstick (in 15 minutes) wind-powered electricity generation fluctuation is following:

${\mathrm{\&Delta;f}}_{h,l}^{\left(2\right)}=\frac{{\mathrm{<figure-callout\; id="0"\; label="f"\; filenames="HDA0000127997390000021.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 (in 15 minutes) causes is Δ P ⁽²⁾, expression formula is following:

${\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 following:

$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, and expression formula is following:

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

Conventional unit low-valley interval went out force vector when

was incorporated into the power networks for no wind-powered electricity generation;

(IV) the electric power system wind-powered electricity generation capacity of can dissolving peak period is following for

expression formula under the fully dark peak regulation state:

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

In the formula (21);

is conventional machine group minimum output peak period of electric power system, and expression formula is following:

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

is conventional unit minimum output vector;

(V) the electric power system low-valley interval wind-powered electricity generation capacity of can dissolving is following for

expression formula under the fully dark peak regulation state:

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

In the formula (23);

is the minimum output of the conventional machine group low-valley interval of electric power system, and expression formula is following:

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

is conventional unit minimum output vector;

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

expression formula is following:

$P_{\mathrm{SYS}}^{(\max -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{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 that

expression formula is following:

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

1-6) related based on operation states of electric power system variable and identification variable, confirm between the feasible region of control range and identification variable of operation states of electric power system variable, and then establish constraints controllable domain third dimension degree optimized recognition:

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

$\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 following:

${\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):

be not for having under the wind-electricity integration condition; But electric power system power supply capacity when conventional unit i fault is only arranged, expression formula is following:

$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)$

be not for having under the wind-electricity integration condition; But electric power system power supply capacity when conventional unit i and j fault are only arranged, expression formula is following:

$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;

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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)$

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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 following:

${\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 institute;

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 following:

${\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}_{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 institute; 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 following:

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

Wherein,

is 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 following:

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

Wherein,

is 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 following:

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

Wherein,

is electric power system maximum power supply capacity peak period;

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

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

Wherein, is the maximum power supply capacity of electric power system low-valley interval;

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

1-7) constraints of the optimized recognition target function of composite type (7) foundation and formula (27)～(39) foundation forms the optimized recognition 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) cycle repeats step 1-3), 1-4), 1-5), 1-6), 1-7), based on the two-dimentional controllable domain Ω that sets, accomplish optimized recognition to controllable domain third dimension degree, obtain the dissolve three-dimensional controllable domain Z of wind-powered electricity generation of electric power system, expression formula is following:

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 forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system, specifically comprise:

2-1) forecasting institute gets the fundamental characteristics of wind-powered electricity generation, specifically comprises:

(I) to get wind-powered electricity generation

expression formula of exerting oneself peak period following for forecasting institute:

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

(II) to get wind-powered electricity generation low-valley interval expression formula of exerting oneself following for forecasting institute:

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

(III) to get wind-powered electricity generation wind-powered electricity generation average output

expression formula following for forecasting institute:

$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,

gets exerting oneself of wind-powered electricity generation day part for forecasting institute;

2-2) get the fundamental characteristics of wind-powered electricity generation based on forecasting institute, the measuring and calculating forecasting institute gets three dimension indicators of wind-powered electricity generation: η ⁽⁰⁾, δ ⁽⁰⁾, γ ⁽⁰⁾, specifically comprise:

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

${\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 following:

${\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 following:

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

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

(I) in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system, optimize location prediction gained wind-powered electricity generation (in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system, forecasting institute is got wind-powered electricity generation navigates to that to optimize the location model expression formula following with the minimum point

of its Euclidean distance:

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

(47)

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

Find the solution this optimization location model (47) and obtain

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

$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, judge that the electric power system forecasting institute of whether can fully dissolving gets wind-powered electricity generation, the criterion expression formula is following:

Forecasting institute gets wind-powered electricity generation and exerts oneself if electric power system can be dissolved, then target r ₀=1; Otherwise, r then ₀=0;

3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation, specifically comprises:

3-1) confirm that forecasting institute is got the wind-powered electricity generation day part abandons the Optimization Model of wind control

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

$\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)}$

Optimization model (50) in,

for all time off wind power output; according to step 2) the method

is off

After the remaining wind power systems are subject to the full amount of consumptive;

3-2) confirm 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 (so that the total capacity of the unit of closing is maximum) of the conventional unit optimal control of structure, expression formula is following:

$\mathrm{<figure-callout\; id="0"\; label="max\; D"\; filenames="HDA0000127997390000021.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 Optimization Model with the constraints of formula (27)～(39) foundation, finds the solution this Optimization 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, specifically comprise in next day:

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 the effective output of ;

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:

With certain provincial area is that example is set forth peak-frequency regulation constraint proposed by the invention can the dissolve identification and the control method of wind-powered electricity generation of electric power system down.This embodiment based on power planning decision support electric power system GOPT5.0 (by the power system planning software kit of Electric Motor Engineering and Applied Electronic Technology Department, Qinghua exploitation; Set up the unified planning Optimization Model of power industry sustainable development; Seek science, economy, continuable electric power development scheme); Under the condition of not considering wind-electricity integration, produce 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, preestablishes 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}；

Confirm 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, EIAJ, 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 electric power system a few days ago;

Adopt method provided by the invention, as shown in Figure 1 at the dissolve controllable domain of wind-powered electricity generation of identification electric power system a few days ago.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 is exerted oneself the poor correct time that is less than 1 (the logarithm value is less than 0) and peak interval of time; Electric power system has the higher wind-powered electricity generation ability of dissolving, and 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 that a few days ago forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system;

Shown in Figure 2 is that the wind-powered electricity generation of prediction a few days ago 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 judged and can't fully be dissolved by electric power system;

3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation;

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 1 (closing), and unit 11 start-up command are 0 (unlatching);

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 wind instruction.

Table 1

Above-described specific embodiment is merely explanation realization effect of the present invention, not in order to restriction the present invention.Modification, conversion and the improvement of any unsubstantiality of being done within all basic ideas and frameworks in method proposed by the invention all should be included within protection scope of the present invention.

Claims (1)

1. peak-frequency regulation constraint can the dissolve identification and the control method of wind-powered electricity generation of electric power system down; 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 electric 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 that a few days ago forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system; 3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation; 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,, specifically may further comprise the steps at the dissolve three-dimensional controllable domain of wind-powered electricity generation of identification electric power system a few days ago:

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 following:

$\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 following:

$\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 following:

$\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 forecasting institute a few days ago gets electric power system load peak period;

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, preestablish 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 following:

Ω＝{(η _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 following:

∏＝{η _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 following:

Δ＝{δ _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 setting of dissolving;

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=1 ... If N} is d _i=0, unit i keeps start; If d _i=1, then unit i is a closed condition; N is conventional unit quantity;

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

$\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-powered electricity generation that electric power system can be dissolved peak period is exerted oneself;

1-5) based on 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 following:

$\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 following:

$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 following:

$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)$

is the EIAJ of 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

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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)$

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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 expression formula that causes of electric power system short time peak period yardstick wind-powered electricity generation fluctuation is following:

${\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 following:

${\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 following:

$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, and expression formula is following:

$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_{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 pairing position of firm outputs such as heat supply, interconnector among the M, and other positions are 0;

goes out force vector peak period for the no wind-powered electricity generation time routine unit that is incorporated into the power networks;

(III) frequency fluctuation

expression formula that causes of electric power system low-valley interval short time yardstick wind-powered electricity generation fluctuation is following:

${\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 (in 15 minutes) causes is Δ P ⁽²⁾, expression formula is following:

${\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 following:

$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, and expression formula is following:

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

Conventional unit low-valley interval went out force vector when

was incorporated into the power networks for no wind-powered electricity generation;

(IV) the electric power system wind-powered electricity generation capacity of can dissolving peak period is following for

expression formula under the fully dark peak regulation state:

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

In the formula (21);

is conventional machine group minimum output peak period of electric power system, and expression formula is following:

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

is conventional unit minimum output vector;

(V) the electric power system low-valley interval wind-powered electricity generation capacity of can dissolving is following for

expression formula under the fully dark peak regulation state:

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

In the formula (23);

is the minimum output of the conventional machine group low-valley interval of electric power system, and expression formula is following:

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

is conventional unit minimum output vector;

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

expression formula is following:

$P_{\mathrm{SYS}}^{(\max -1)}=[E-D-M]\&CenterDot;{\left[P_{\mathrm{GEN}}^{\max }\right]}^T+M\&CenterDot;{\left[P_{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 that

expression formula is following:

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

1-6) related based on operation states of electric power system variable and identification variable, confirm between the feasible region of control range and identification variable of operation states of electric power system variable, and then establish constraints controllable domain third dimension degree optimized recognition:

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

$\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 following:

${\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):

be not for having under the wind-electricity integration condition; But electric power system power supply capacity when conventional unit i fault is only arranged, expression formula is following:

$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)$

be not for having under the wind-electricity integration condition; But electric power system power supply capacity when conventional unit i and j fault are only arranged, expression formula is following:

$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;

is the accumulation hourage of power system load greater than

, and expression formula is following:

$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)$

is the accumulation hourage of power system load greater than , and expression formula is following:

$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 following:

${\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 institute;

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 following:

${\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}_{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 institute;

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 following:

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

Wherein,

is 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 following:

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

Wherein,

is 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 following:

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

Wherein,

is electric power system maximum power supply capacity peak period;

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

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

Wherein, is the maximum power supply capacity of electric power system low-valley interval;

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

1-7) constraints of the optimized recognition target function of composite type (7) foundation and formula (27)～(39) foundation forms the optimized recognition 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) cycle repeats step 1-3), 1-4), 1-5), 1-6), 1-7), based on the two-dimentional controllable domain Ω that sets, accomplish optimized recognition to controllable domain third dimension degree, obtain the dissolve three-dimensional controllable domain Z of wind-powered electricity generation of electric power system, expression formula is following:

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 forecasting institute gets wind-powered electricity generation and whether can fully be dissolved by electric power system, specifically comprise:

2-1) forecasting institute gets the fundamental characteristics of wind-powered electricity generation, specifically comprises:

(I) to get wind-powered electricity generation

expression formula of exerting oneself peak period following for forecasting institute:

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

(II) to get wind-powered electricity generation low-valley interval

expression formula of exerting oneself following for forecasting institute:

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

(III) to get wind-powered electricity generation wind-powered electricity generation average output

expression formula following for forecasting institute:

$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,

gets exerting oneself of wind-powered electricity generation day part for forecasting institute;

2-2) get the fundamental characteristics of wind-powered electricity generation based on forecasting institute, the measuring and calculating forecasting institute gets three dimension indicators of wind-powered electricity generation: η ⁽⁰⁾, δ ⁽⁰⁾, γ ⁽⁰⁾, specifically comprise:

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

${\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 following:

${\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 following:

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

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

(I) in the two-dimentional controllable domain Ω of wind-powered electricity generation is dissolved in electric power system, optimize location prediction gained wind-powered electricity generation, it is following to optimize the location model expression formula:

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

(47)

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

Find the solution this optimization location model (47) and obtain

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

$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, judge that the electric power system forecasting institute of whether can fully dissolving gets wind-powered electricity generation, the criterion expression formula is following:

Forecasting institute gets wind-powered electricity generation and exerts oneself if electric power system can be dissolved, then target r ₀=1; Otherwise, r then ₀=0;

3) get the judged result whether wind-powered electricity generation can fully be dissolved by electric power system according to forecasting institute, confirming that a few days ago the electric power system forecasting institute of dissolving gets the optimal control policy of wind-powered electricity generation, specifically comprises:

3-1) confirm forecasting institute is got the Optimization Model that the wind-powered electricity generation day part is abandoned wind control

, expression formula is following:

$\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)}$

Optimization model (50) in,

time off of each wind power output; in step 2) obtained by the method described in

is off

After the wind power system is available for the remaining full the subject of consumptive amount;

3-2) confirm 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 following:

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

(III) integrated objective function (51) forms conventional unit optimal control Optimization Model with the constraints of formula (27)～(39) foundation, finds the solution this Optimization 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, specifically comprise in next day:

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 the effective output of

;

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