CN105914779A - Coordination control method for wind turbine generator set participating in automatic power generation control of power system - Google Patents

Abstract

The invention discloses a coordination control method for a wind turbine generator set participating in the automatic power generation control of a power system, and the method comprises the steps that a wind power plant transmits the ultrashort term predicted values of real-time power output and wind power output of the wind turbine generator set to a provincial dispatching center; the provincial dispatching center calculates the adjustment capacity and creep speed of the wind turbine generator set, capable of participating in the automatic power generation control, in an administrative region; a regional dispatching center employs model prediction to control and calculate the adjustment power of the wind turbine generator set and a conventional set; the provincial dispatching center distributes the adjustment power according to a proportion; and the wind power plant adjusts the wind active power output through a wind power plant AGC system. The method is advantageous in that the method effectively evaluates the adjustment capacity and performance of the wind turbine generator set participating in the automatic power generation control of the power system through the mutual coordination of the regional dispatching center, the provincial dispatching center and the wind power plant; the method gives full consideration to the operation characteristics and adjustment capabilities of various types of frequency modulation resources in the power system; the method achieves a purpose that the wind turbine generator set participates the automatic power generation control of the power system; and the method improves the stability and economics of a system.

Description

A kind of Wind turbines participates in the control method for coordinating of power system Automatic Generation Control

Technical field

The invention belongs to power system Automatic Generation Control field, be specifically related to a kind of Wind turbines and participate in electric power The control method for coordinating of system Automatic Generation Control.

Background technology

Power system frequency is the important indicator of the quality of power supply, that reflects between generating active power and load Equilibrium relation.In order to ensure the frequency stable of system, meet subscriber equipment, power plant equipment and power train Unite properly functioning needs, power system Automatic Generation Control (Automation Generation Control, AGC) it is widely used, exerts oneself adjusting generating set, meet the user's request being continually changing.

In recent years, under the dual promotion of energy crisis and environmental protection pressure, the new forms of energy such as wind-powered electricity generation, solar energy are sent out Power technology development is swift and violent, and the new forms of energy ratio in power system improves constantly.But, due to new energy such as wind-powered electricity generations Source generating has the feature such as undulatory property, randomness, many by generation of electricity by new energy in electric power system dispatching field at present Regard as " bearing " load to process, conventional rack provide standby for it and go out fluctuation compensating it, this will cause be Increasing the weight of of the tracking frequency rising of demand and the abrasion of frequency modulation unit, improves the cost of investment of power system and runs into This.Additionally, in the northern area of China, Wind turbines is big more than entrance in winter sends out period, but now by heat supply Demands etc. affect, the frequency regulation capacity shortage of power system, reduce power system further to intermittent renewable The receiving ability of the energy, causes " abandoning wind " phenomenon occur, causes huge economic loss and environmental protection pressure, no It is beneficial to national energy-saving reduce discharging and the realization of carbon emission reduction target.

Along with wind-power electricity generation controls the development of technology, Wind turbines accepts dispatching of power netwoks instruction, participates in electricity The ability of net management and running constantly promotes.Especially, wind-powered electricity generation participates in power system Automatic Generation Control firmly It is possibly realized on part: according to statistics, at present in the wind energy turbine set of China Ji Beidiqu, wind energy turbine set AGC is installed The ratio of system alreadys more than 70%, and threads at different levels had configured the dispatching patcher adapted already.Consider Due under the conditions of the limit wind that peak regulation difficulty causes, if wind power output and profile constraints still have more than needed, can be about Restraint and on the premise of its frequency modulation is exerted oneself more than zero, make Wind turbines participate in power system Automatic Generation Control, promoting While power system frequency modulation performance, it is achieved the additional issue of Wind turbines, promote stability and warp that system is run Ji property.

But, Wind turbines participates in the performance of power system Automatic Generation Control by the shadow by following factor Ring: the frequency modulation of the Wind turbines upper limit of exerting oneself is retrained by the wind power output upper limit, although wind power output at present is ultrashort Phase precision of prediction constantly promotes, but it still has certain error；It is different from conventional rack, the climbing of Wind turbines Speed not constant, it changes with the change of wind-powered electricity generation EIAJ；For realizing the purpose of wind-powered electricity generation additional issue, wind-powered electricity generation When participating in power system Automatic Generation Control, its frequency modulation is exerted oneself and should be consistently greater than zero.It is thus impossible to simply Wind turbines is made no exception with conventional rack, is directly brought in existing AGC system.

Summary of the invention

It is contemplated that at least solve one of above-mentioned technical problem.

To this end, it is an object of the present invention to propose a kind of Wind turbines to participate in power system automatic generation control The control method for coordinating of system.

To achieve these goals, embodiment of the invention discloses that a kind of Wind turbines participates in power system certainly The control method for coordinating of dynamic Generation Control, comprises the following steps: S1: wind energy turbine set is uploaded to provincial control centre Exerting oneself in real time and the ultra-short term predictive value of wind power output of Wind turbines；S2: described provincial control centre according to Profile constraints information and wind energy turbine set unit performance information, calculate the local wind that may participate in Automatic Generation Control Group of motors pondage and creep speed；S3: subdispatch center is according to district control deviation, it is considered to wind Group of motors and conventional rack participate in characteristic and the constraint of power system Automatic Generation Control, use model prediction control System calculates the regulation power of Wind turbines and conventional rack；S4: described provincial control centre is according to each wind energy turbine set Pondage pro-rata regulation power；And S5: wind energy turbine set is adjusted by wind energy turbine set AGC system Wind-powered electricity generation is meritorious exerts oneself, it is achieved the response to regulation power instruction.

Wind turbines according to embodiments of the present invention participates in the coordination controlling party of power system Automatic Generation Control Method, takes into full account that Wind turbines and conventional rack participate in characteristic and the constraint of power system Automatic Generation Control, While promoting power system frequency modulation performance, it is achieved the additional issue of Wind turbines, promote Operation of Electric Systems Stability and economy.

It addition, Wind turbines according to the above embodiment of the present invention participates in the association of power system Automatic Generation Control Control method, it is also possible to there is following additional technical characteristic:

Further, in step sl, the P that exerts oneself in real time of described Wind turbines_W,i, i=1,2 ..., N, wind The ultra-short term predictive value P that electricity is exerted oneself_F,i(t), i=1,2 ..., N, t=1,2 ..., P, wherein, N is in provincial scheduling The wind energy turbine set quantity of the heart local installation AGC system, P is that described subdispatch center performs model prediction The prediction time domain length controlled.

Further, step S2 farther includes:

S201: calculate the local Wind turbines regulation that may participate in Automatic Generation Control by below equation and hold Amount:

$C_{W,k}^{up}\left(t\right)=\min (S_k-\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{W,i}-\underset{i=1}{\overset{N_{G,k}}{\Σ}}{P}_{G,i},\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{F,i}(t)-\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{W,i}),k=1,2,\mathrm{...},J$

$C_W^{up}\left(t\right)=\underset{k=1}{\overset{J}{\Σ}}{C}_{W,k}^{up}\left(t\right)$

$C_W^{down}\left(t\right)=\underset{i=1}{\overset{N}{\Σ}}(P_{W,i}-P_{C,i})$

Wherein,S_kBe respectively under section k is adjusted up capacity and section capacity-constrained, N_W,kWith N_G,kIt is respectively the wind-powered electricity generation number of fields under section k and conventional power plant number, P_G,iExerting oneself in real time for conventional power unit, J is Section sum, P_C,iHaving under the conditions of wind energy turbine set limit wind during for not considering Wind turbines participation Automatic Generation Control Merit power instruction；And

S202: calculate the local creep speed that may participate in Automatic Generation Control by below equation:

$R_W^{up}\left(t\right)=min\left(C_W^{up}\right(t)/T,\underset{i=1}{\overset{N}{\Σ}}{R}_{W,i}^{up}(t\left)\right)$

$R_W^{down}\left(t\right)=min\left(C_W^{down}\right(t)/T,\underset{i=1}{\overset{N}{\Σ}}{R}_{W,i}^{down}(t\left)\right)$

Wherein,It is respectively the upwards creep speed of each wind energy turbine set and downward creep speed, Described upwards creep speed and described downward creep speed all with the ultra-short term predictive value P of wind power output_F,i(t) phase Close.

Further, step S3 farther includes: S301: according to current frequency deviation f and interconnection Exchange power deviation Δ P_tieZoning control deviation ACE:ACE=B Δ f+ Δ P_tie, wherein, B is system frequency Rate deviation factor；S302: near Operation of Electric Systems point, region difference control process is carried out linearisation Processing, dynamic mathematical models based on conventional rack Yu Wind turbines, by power system region difference control mistake The state space equation of journey dynamic mathematical models is written as form:Y=CX, Wherein, X represents state variable vector, and U represents input variable vector, and W represents disturbance quantity vector, Y generation Table output vector, A represents the state matrix of this system, and B represents the input matrix of this system, and R represents Perturbation matrix, C represents output matrix；State space equation is carried out sliding-model control, obtains power system The discrete form of region difference control process mathematical model:

X (k+1)=A_dx(k)+B_du(k)+R_dW(k)

In above formula, the expression formula of each matrix is as follows:

$A_d=e^{{\mathrm{AT}}_s},B_d=\underset{0}{\overset{T_s}{\∫}}{e}^{At}dt\·B,R_d=\underset{0}{\overset{T_s}{\∫}}{e}^{At}dt\·R$

Wherein, T_sSampling time for discrete system；S303: determine Controlling object function such as following formula:

$\begin{array}{c}\underset{U}{\min }J(U,x(k\left)\right)=\underset{i=1}{\overset{P}{\Σ}}{\[x_G(k+i)\]}^T{Q}_G{x}_G(k+i)+\underset{i=1}{\overset{P}{\Σ}}{\[x_W(k+i)\]}^T{Q}_W{x}_W(k+i)\\ +\underset{i=0}{\overset{M-1}{\Σ}}{\[u_G(k+i)\]}^T{R}_G{u}_G(k+i)+\underset{i=0}{\overset{M-1}{\Σ}}{\[u_W(k+i)\]}^T{R}_W{u}_W(k+i)\end{array}$

Wherein, Q_G、Q_W、R_G、R_WIt is respectively conventional rack quantity of state, Wind turbines quantity of state, tradition Unit allocation amount and the weight matrix of Wind turbines controlled quentity controlled variable, M is that the control time domain of Model Predictive Control is long Degree；S304: determine that constraints is as follows:

X (k+i+1 | k)=A_dx(k+i|k)+B_du(k+i|k)+R_dW(k+i|k),i≥0

Represent, the linear prediction model of power system region difference control process；

u_Gmin≤u_G(k+i|k)≤u_Gmax, i=0 ..., M-1

u_Wmin≤u_W(k+i|k)≤u_Wmax, i=0 ..., M-1

Representing, conventional rack controlled quentity controlled variable and Wind turbines controlled quentity controlled variable are all controlled the constraint in interval；

x_Gmin≤x_G(k+i|k)≤x_Gmax, i=1 ..., P

x_Wmin≤x_W(k+i|k)≤x_Wmax, i=1 ..., P

Representing, conventional rack quantity of state is all retrained by its pondage with Wind turbines state；

Δx_Gmin≤Δx_G(k+i|k)≤Δx_Gmax, i=1 ..., P

Δx_Wmin≤Δx_W(k+i|k)≤Δx_Wmax, i=1 ..., P

Representing, conventional rack quantity of state increment is all retrained by its creep speed with Wind turbines quantity of state increment. S305: solve according to active constraint set algorithm, calculates described conventional rack regulation power and described wind Group of motors regulation power.

Further, further comprising the steps of after step S305: S306: described conventional rack is adjusted Joint power and described Wind turbines regulation power carry out feedback compensation and rolling optimization；S307: by described tradition Unit regulation power and described Wind turbines regulation power are handed down to described provincial control centre, by described provincial Regulation power is allocated by control centre, and is issued to each wind energy turbine set.

Further, in step s 4, described provincial control centre leads to according to the pondage of each wind energy turbine set Cross below equation pro-rata regulation power A_w:

$\left\{\begin{array}{c}{A}_{W,i}\left(t\right)=\min (A_W\·\frac{C_{W,i}^{up}\left(t\right)}{C_W^{up}\left(t\right)},R_{W,i}^{up}(t)\·T),i=1,2,\mathrm{...},N,A_w\≥0\\ A_{W,i}\left(t\right)=\min (A_W\·\frac{C_{W,i}^{down}\left(t\right)}{C_W^{down}\left(t\right)},R_{W,i}^{down}(t)\·T),i=1,2,\mathrm{...},N,A_w\≤0\end{array}\right..$

The additional aspect of the present invention and advantage will part be given in the following description, and part will be retouched from following Become obvious in stating, or recognized by the practice of the present invention.

Accompanying drawing explanation

Above-mentioned and/or the additional aspect of the present invention and advantage are from combining the accompanying drawings below description to embodiment Will be apparent from easy to understand, wherein:

Fig. 1 is that the Wind turbines of one embodiment of the invention participates in the coordination control of power system Automatic Generation Control The flow chart of method processed；

Fig. 2 is the physics of the Wind turbines participation power system Automatic Generation Control of one embodiment of the invention Organization Chart；

Fig. 3 is that the Wind turbines of one embodiment of the invention participates in the coordination control of power system Automatic Generation Control The theory diagram of system strategy.

Detailed description of the invention

Embodiments of the invention are described below in detail, and the example of described embodiment is shown in the drawings, wherein certainly Begin to same or similar label eventually represent same or similar element or there is the unit of same or like function Part.The embodiment described below with reference to accompanying drawing is exemplary, is only used for explaining the present invention, and can not It is interpreted as limitation of the present invention.

In describing the invention, it is to be understood that term " " center ", " longitudinally ", " laterally ", " on ", D score, "front", "rear", "left", "right", " vertically ", " level ", " top ", " end ", " interior ", " outward " etc. Orientation or the position relationship of instruction are based on orientation shown in the drawings or position relationship, are for only for ease of description The present invention and simplification describe rather than indicate or imply that the device of indication or element must have specific side Position, with specific azimuth configuration and operation, be therefore not considered as limiting the invention.Additionally, term " first ", " second " are only used for describing purpose, and it is not intended that indicate or hint relative importance.

In describing the invention, it should be noted that unless otherwise clearly defined and limited, term " peace Dress ", should be interpreted broadly " being connected ", " connection ", for example, it may be fix connection, it is also possible to be removable Unload connection, or be integrally connected；Can be to be mechanically connected, it is also possible to be electrical connection；Can be to be joined directly together, Can also be indirectly connected to by intermediary, can be the connection of two element internals.General for this area For logical technical staff, above-mentioned term concrete meaning in the present invention can be understood with concrete condition.

With reference to explained below and accompanying drawing, it will be clear that these and other aspects of embodiments of the invention.At this In a little descriptions and accompanying drawing, specifically disclose some particular implementation in embodiments of the invention, represent Implement some modes of the principle of embodiments of the invention, but it is to be understood that the model of embodiments of the invention Enclose not limited.On the contrary, embodiments of the invention include falling into the spirit of attached claims and interior All changes, amendment and equivalent in the range of culvert.

Describe Wind turbines according to embodiments of the present invention below in conjunction with accompanying drawing and participate in power system automatic generation The control method for coordinating controlled.

Refer to Fig. 1-3, a kind of Wind turbines participates in the control method for coordinating of power system Automatic Generation Control, Comprise the following steps:

S1: wind energy turbine set uploads exerting oneself in real time and the ultra-short term of wind power output of Wind turbines to provincial control centre Predictive value.

Specifically, before each AGC instruction cycle T starts, each wind energy turbine set calling module 1 and module The P that exerts oneself in real time of 2 monitoring Wind turbines_W,i(MW), i=1,2 ..., N, calculate the ultra-short term prediction of wind power output Value P_F,i(t), i=1,2 ..., N, t=1,2 ..., P, wherein, N is the local installation in provincial control centre AGC The wind energy turbine set quantity of system, P is the prediction time domain length that subdispatch center performs Model Predictive Control.Adjust Above-mentioned information is reported to provincial control centre by module 3.

S2: provincial control centre, according to profile constraints information and wind energy turbine set unit performance information, calculates local May participate in Wind turbines pondage and the creep speed of Automatic Generation Control.

Specifically, after provincial control centre receives the information of its local all wind energy turbine set, calling module 4, read profile constraints information and wind energy turbine set unit performance information, calling module 5 calculates local may participate in To Wind turbines pondage and the creep speed of Automatic Generation Control, concrete calculation procedure is as follows:

S201: calculate local Wind turbines and participate in the pondage of Automatic Generation Control.According to Wind turbines Exert oneself in real time, the prediction of the ultra-short term of wind power output and profile constraints, calculate local Wind turbines and participate in Automatic Generation Control be adjusted up capacityPondage downwards

$C_{W,k}^{up}\left(t\right)=\min (S_k-\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{W,i}-\underset{i=1}{\overset{N_{G,k}}{\Σ}}{P}_{G,i},\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{F,i}(t)-\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{W,i}),k=1,2,\mathrm{...},J$

$C_W^{up}\left(t\right)=\underset{k=1}{\overset{J}{\Σ}}{C}_{W,k}^{up}\left(t\right)$

$C_W^{down}\left(t\right)=\underset{i=1}{\overset{N}{\Σ}}(P_{W,i}-P_{C,i})$

Wherein,S_kBe respectively under section k is adjusted up capacity (MW) and section capacity-constrained (MW), N_W,kAnd N_G,kIt is respectively the wind-powered electricity generation number of fields under section k and conventional power plant number, P_G,iReal for conventional power unit Time exert oneself (MW), J be section sum, P_C,iFor do not consider Wind turbines participate in Automatic Generation Control time wind-powered electricity generation Active power instruction (MW) under the conditions of the limit wind of field.

S202: calculate local Wind turbines and participate in the creep speed of Automatic Generation Control.According to local The Wind turbines pondage and the wind energy turbine set unit performance that may participate in Automatic Generation Control calculate area under one's jurisdiction endogenous wind Electric field participates in the upwards creep speed of Automatic Generation ControlUpwards creep speed

$R_W^{up}\left(t\right)=min\left(C_W^{up}\right(t)/T,\underset{i=1}{\overset{N}{\Σ}}{R}_{W,i}^{up}(t\left)\right)$

$R_W^{down}\left(t\right)=min\left(C_W^{down}\right(t)/T,\underset{i=1}{\overset{N}{\Σ}}{R}_{W,i}^{down}(t\left)\right)$

Wherein,It is respectively the upwards creep speed of each wind energy turbine set and downward creep speed, Its all with the ultra-short term predictive value P of wind power output_F,iT () is correlated with.

S3: subdispatch center is according to district control deviation, it is considered to Wind turbines participates in electric power with conventional rack The characteristic of system Automatic Generation Control and constraint, use Model Predictive Control to calculate Wind turbines and conventional rack Regulation power.Step S3 comprises the following steps:

S301: zoning control deviation ACE.Merit is exchanged with interconnection according to current frequency deviation f Rate deviation delta P_tieZoning control deviation ACE:

ACE=B Δ f+ Δ P_tie

Wherein, B is system frequency deviation coefficient (MW/Hz).

S302: calculate the forecast model of power system region difference control process.Attached at Operation of Electric Systems point Closely, region difference control process is carried out linearization process, dynamic number based on conventional rack Yu Wind turbines Learn model, the state space equation of power system region difference control process dynamics mathematical model is written as Form:

$\stackrel{\·}{X}=AX+BU+RW$

Y=CX

Wherein, X represents state variable vector, and U represents input variable vector, and W represents disturbance quantity vector, Y represents output vector, and A represents the state matrix of this system, and B represents the input matrix of this system, R generation Table perturbation matrix, C represents output matrix.

Owing to the dynamic mathematical models of power system region difference control process are continuous system model, and model The forecast model used in PREDICTIVE CONTROL is discrete system model, needs state space equation is carried out discretization Process, obtain the discrete form of power system region difference control process mathematical model:

X (k+1)=A_dx(k)+B_du(k)+R_dW(k)

In above formula, the expression formula of each matrix is as follows:

$A_d=e^{{\mathrm{AT}}_s},B_d=\underset{0}{\overset{T_s}{\∫}}{e}^{At}dt\·B,R_d=\underset{0}{\overset{T_s}{\∫}}{e}^{At}dt\·R$

Wherein, T_sSampling time for discrete system.

S303: determine control target.

$\begin{array}{c}\underset{U}{\min }J(U,x(k\left)\right)=\underset{i=1}{\overset{P}{\Σ}}{\[x_G(k+i)\]}^T{Q}_G{x}_G(k+i)+\underset{i=1}{\overset{P}{\Σ}}{\[x_W(k+i)\]}^T{Q}_W{x}_W(k+i)\\ +\underset{i=0}{\overset{M-1}{\Σ}}{\[u_G(k+i)\]}^T{R}_G{u}_G(k+i)+\underset{i=0}{\overset{M-1}{\Σ}}{\[u_W(k+i)\]}^T{R}_W{u}_W(k+i)\end{array}$

Wherein, Q_G、Q_W、R_G、R_WIt is respectively conventional rack quantity of state, Wind turbines quantity of state, tradition Unit allocation amount and the weight matrix of Wind turbines controlled quentity controlled variable, M is that the control time domain of Model Predictive Control is long Degree.

S304: determine constraints.

X (k+i+1 | k)=A_dx(k+i|k)+B_du(k+i|k)+R_dW(k+i|k),i≥0

Represent, the linear prediction model of power system region difference control process.

u_Gmin≤u_G(k+i|k)≤u_Gmax, i=0 ..., M-1

u_Wmin≤u_W(k+i|k)≤u_Wmax, i=0 ..., M-1

Representing, conventional rack controlled quentity controlled variable and Wind turbines controlled quentity controlled variable are all controlled the constraint in interval.

x_Gmin≤x_G(k+i|k)≤x_Gmax, i=1 ..., P

x_Wmin≤x_W(k+i|k)≤x_Wmax, i=1 ..., P

Representing, conventional rack quantity of state is all retrained by its pondage with Wind turbines state.

Δx_Gmin≤Δx_G(k+i|k)≤Δx_Gmax, i=1 ..., P

Δx_Wmin≤Δx_W(k+i|k)≤Δx_Wmax, i=1 ..., P

Representing, conventional rack quantity of state increment is all retrained by its creep speed with Wind turbines quantity of state increment.

S305: ask for conventional rack regulation power and regulate power with Wind turbines.The above-mentioned mathematical model set up Actual is a quadratic programming problem, it is possible to use action constaint set scheduling algorithm solves.Right for realizing The dynamic tracking control of target desired value, also needs to carry out feedback compensation, rolling optimization.After having solved, should Calculated conventional rack regulation power is issued with Wind turbines regulation power, wherein Wind turbines Regulation power should be issued to provincial control centre, provincial control centre be allocated regulation power, and under Send to each wind energy turbine set.

S4: provincial control centre pro rata distributes according to the pondage of each wind energy turbine set and regulates power.

Specifically, after provincial control centre receives the Wind turbines regulation power instruction issued, calling module 8, pro rata distribute according to the pondage of each wind energy turbine set and regulate power A_w, it is issued to each wind energy turbine set:

$\left\{\begin{array}{c}{A}_{W,i}\left(t\right)=\min (A_W\·\frac{C_{W,i}^{up}\left(t\right)}{C_W^{up}\left(t\right)},R_{W,i}^{up}(t)\·T),i=1,2,\mathrm{...},N,A_w\≥0\\ A_{W,i}\left(t\right)=\min (A_W\·\frac{C_{W,i}^{down}\left(t\right)}{C_W^{down}\left(t\right)},R_{W,i}^{down}(t)\·T),i=1,2,\mathrm{...},N,A_w\≤0\end{array}\right.$

Especially, in order to ensure to regulate the implementation effect of power instruction, reward the wind energy turbine set that implementation effect is good, It is considered as each wind energy turbine set response condition to upper period regulation power instruction, for being adjusted up power instruction, Can preferentially distribute to the wind energy turbine set good to upper period regulation power instruction response, for regulating downwards power Instruction, can be proportionately distributed to the wind energy turbine set of all participation Automatic Generation Control.

S5: wind energy turbine set is exerted oneself by wind energy turbine set AGC system adjustment wind-powered electricity generation is meritorious, it is achieved refer to regulation power The response of order.

It addition, the Wind turbines of the embodiment of the present invention participates in the coordination controlling party of power system Automatic Generation Control Other composition of method and effect are the most all known, in order to reduce redundancy, Do not repeat.

In the description of this specification, reference term " embodiment ", " some embodiments ", " example ", The description of " concrete example " or " some examples " etc. means to combine this embodiment or example describes specific features, Structure, material or feature are contained at least one embodiment or the example of the present invention.In this manual, The schematic representation of above-mentioned term is not necessarily referring to identical embodiment or example.And, the tool of description Body characteristics, structure, material or feature can be with properly in any one or more embodiments or example Mode combine.

Although an embodiment of the present invention has been shown and described, it will be understood by those skilled in the art that: These embodiments can be carried out in the case of without departing from the principle of the present invention and objective multiple change, amendment, Replacing and modification, the scope of the present invention is limited by claim and equivalent thereof.

Claims (6)

1. Wind turbines participates in a control method for coordinating for power system Automatic Generation Control, and its feature exists In, comprise the following steps:

S1: wind energy turbine set uploads exerting oneself in real time and the ultra-short term of wind power output of Wind turbines to provincial control centre Predictive value；

S2: described provincial control centre, according to profile constraints information and wind energy turbine set unit performance information, calculates linchpin Wind turbines pondage and the creep speed of Automatic Generation Control is may participate in district；

S3: subdispatch center is according to district control deviation, it is considered to Wind turbines participates in electric power with conventional rack The characteristic of system Automatic Generation Control and constraint, use Model Predictive Control to calculate Wind turbines and conventional rack Regulation power；

S4: described provincial control centre pro rata distributes according to the pondage of each wind energy turbine set and regulates power；With And

S5: wind energy turbine set is exerted oneself by wind energy turbine set AGC system adjustment wind-powered electricity generation is meritorious, it is achieved refer to regulation power The response of order.

Wind turbines the most according to claim 1 participates in the coordination control of power system Automatic Generation Control Method processed, it is characterised in that in step sl, the P that exerts oneself in real time of described Wind turbines_W,i, i=1,2 ..., N, The ultra-short term predictive value P of wind power output_F,i(t), i=1,2 ..., N, t=1,2 ..., P, wherein, N is provincial scheduling The wind energy turbine set quantity of the local installation in center AGC system, P is that described subdispatch center execution model is pre- The prediction time domain length of observing and controlling.

Wind turbines the most according to claim 2 participates in the coordination control of power system Automatic Generation Control Method processed, it is characterised in that step S2 farther includes:

S201: calculate the local Wind turbines regulation that may participate in Automatic Generation Control by below equation and hold Amount:

$C_{W,k}^{up}\left(t\right)=min(S_k-\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{W,i}-\underset{i=1}{\overset{N_{G,k}}{\Σ}}{P}_{G,i},\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{F,i}(t)-\underset{i=1}{\overset{N_{W,k}}{\Σ}}{P}_{W,i}),k=1,2,...,J$

$C_W^{up}\left(t\right)=\underset{k=1}{\overset{J}{\Σ}}{C}_{W,k}^{up}\left(t\right)$

$C_W^{down}\left(t\right)=\underset{i=1}{\overset{N}{\Σ}}(P_{W,i}-P_{C,i})$

Wherein,S_kBe respectively under section k is adjusted up capacity and section capacity-constrained, N_W,kWith N_G,kIt is respectively the wind-powered electricity generation number of fields under section k and conventional power plant number, P_G,iExerting oneself in real time for conventional power unit, J is Section sum, P_C,iHaving under the conditions of wind energy turbine set limit wind during for not considering Wind turbines participation Automatic Generation Control Merit power instruction；And

S202: calculate the local creep speed that may participate in Automatic Generation Control by below equation:

$R_W^{up}\left(t\right)=min\left(C_W^{up}\right(t)/T,\underset{i=1}{\overset{N}{\Σ}}{R}_{W,i}^{up}(t\left)\right)$

$R_W^{down}\left(t\right)=min\left(C_W^{down}\right(t)/T,\underset{i=1}{\overset{N}{\Σ}}{R}_{W,i}^{down}(t\left)\right)$

Wherein,It is respectively the upwards creep speed of each wind energy turbine set and downward creep speed, Described upwards creep speed and described downward creep speed all with the ultra-short term predictive value P of wind power output_F,i(t) phase Close.

Wind turbines the most according to claim 1 participates in the coordination control of power system Automatic Generation Control Method processed, it is characterised in that step S3 farther includes:

S301: according to current frequency deviation f and Tie line Power deviation delta P_tieZoning controls partially Difference ACE:

ACE=B Δ f+ Δ P_tie

Wherein, B is system frequency deviation coefficient；

S302: near Operation of Electric Systems point, region difference control process is carried out linearization process, base In the dynamic mathematical models of conventional rack Yu Wind turbines, by power system region difference control process dynamics number The state space equation learning model is written as form:

$\stackrel{\·}{X}=AX+BU+RW$

Y=CX

Wherein, X represents state variable vector, and U represents input variable vector, and W represents disturbance quantity vector, Y represents output vector, and A represents the state matrix of this system, and B represents the input matrix of this system, R generation Table perturbation matrix, C represents output matrix；

State space equation is carried out sliding-model control, obtains power system region difference control process mathematical modulo The discrete form of type:

X (k+1)=A_dx(k)+B_du(k)+R_dW(k)

In above formula, the expression formula of each matrix is as follows:

$A_d=e^{{\mathrm{AT}}_s},B_d=\underset{0}{\overset{T_s}{\∫}}{e}^{At}dt\·B,R_d=\underset{0}{\overset{T_s}{\∫}}{e}^{At}dt\·R$

Wherein, T_sSampling time for discrete system；

S303: determine Controlling object function such as following formula:

$\begin{array}{c}\underset{U}{min}J(U,x(k\left)\right)=\underset{i=1}{\overset{p}{\Σ}}{\[x_G(k+i)\]}^T{Q}_G{x}_G(k+i)+\underset{i=1}{\overset{p}{\Σ}}{\[x_W(k+i)\]}^T{Q}_W{x}_W(k+i)\\ +\underset{i=0}{\overset{M-1}{\Σ}}{\[u_G(k+i)\]}^T{R}_G{u}_G(k+i)+\underset{i=0}{\overset{M-1}{\Σ}}{\[u_W(k+i)\]}^T{R}_W{u}_W(k+i)\end{array}$

Wherein, Q_G、Q_W、R_G、R_WIt is respectively conventional rack quantity of state, Wind turbines quantity of state, tradition Unit allocation amount and the weight matrix of Wind turbines controlled quentity controlled variable, M is that the control time domain of Model Predictive Control is long Degree；

S304: determine that constraints is as follows:

X (k+i+1 | k)=A_dx(k+i|k)+B_du(k+i|k)+R_dW(k+i|k),i≥0

Represent, the linear prediction model of power system region difference control process；

u_Gmin≤u_G(k+i|k)≤u_Gmax, i=0 ..., M-1

u_Wmin≤u_W(k+i|k)≤u_Wmax, i=0 ..., M-1

Representing, conventional rack controlled quentity controlled variable and Wind turbines controlled quentity controlled variable are all controlled the constraint in interval；

x_Gmin≤x_G(k+i|k)≤x_Gmax, i=1 ..., P

x_Wmin≤x_W(k+i|k)≤x_Wmax, i=1 ..., P

Representing, conventional rack quantity of state is all retrained by its pondage with Wind turbines state；

Δx_Gmin≤Δx_G(k+i|k)≤Δx_Gmax, i=1 ..., P

Δx_Wmin≤Δx_W(k+i|k)≤Δx_Wmax, i=1 ..., P

Representing, conventional rack quantity of state increment is all retrained by its creep speed with Wind turbines quantity of state increment.

S305: solve according to active constraint set algorithm, calculates described conventional rack regulation power and institute State Wind turbines regulation power.

Wind turbines the most according to claim 4 participates in the coordination control of power system Automatic Generation Control Method processed, it is characterised in that further comprising the steps of after step S305:

S306: to described conventional rack regulation power and described Wind turbines regulation power carry out feedback compensation and Rolling optimization；

S307: described conventional rack regulation power and described Wind turbines regulation power are handed down to described provincial Control centre, is allocated regulation power by described provincial control centre, and is issued to each wind energy turbine set.

Wind turbines the most according to claim 5 participates in the coordination control of power system Automatic Generation Control Method processed, it is characterised in that in step s 4, described provincial control centre is according to the regulation of each wind energy turbine set Capacity is by below equation pro-rata regulation power A_w:

$\left\{\begin{array}{c}{A}_{W,i}\left(t\right)=\min (A_W\·\frac{C_{W,i}^{up}\left(t\right)}{C_W^{up}\left(t\right)},R_{W,i}^{up}(t)\·T),i=1,2,...,N,A_w\≥0\\ A_{W,i}\left(t\right)=\min (A_W\·\frac{C_{W,i}^{down}\left(t\right)}{C_W^{down}\left(t\right)},R_{W,i}^{down}(t)\·T),i=1,2,...,N,A_w\≤0\end{array}\right..$