Summary of the invention
Regarding the issue above, the present invention provides the wind-powered electricity generation virtual inertia response of a kind of meter and Frequency Index constraint
Consider with the wind power penetration limit calculation method of primary frequency modulation response using large-scale wind power access electric system as object
The response of wind-powered electricity generation virtual inertia and primary frequency modulation response effect, wind-powered electricity generation is calculated in terms of Frequency Index constraint penetrates power pole
Limit.
The technical scheme adopted by the invention is that a kind of wind power penetration limit for considering frequency constraint and wind-powered electricity generation frequency modulation
Calculation method specifically includes following steps,
Step 1:Using the dynamic equivalent parameter aggregation method of weighting, respectively to wind turbine group and thermal motor group, water
Motor group
Group carries out parameter equivalent calculation, seeks the parameter K of its equivalent unitG;
Step 2:The coupling for considering the response of wind-powered electricity generation virtual inertia, the system calculated under different wind-powered electricity generation permeabilities are equivalent used
The property time
Constant H∑;
Step 3:According to the dynamic response model transmission function h of the revolving speed primary frequency modulation control system of separate unit Wind turbineswt
(s), it adopts
With weighting dynamic equivalent parameter aggregation method, calculates and solve wind power plant based on the dynamic of revolving speed primary frequency modulation control system
State response
Model transfer function hmWF(s) and the active increment Delta P of wind-powered electricity generation primary frequency modulationmWF(s);
Step 4:Establish steam turbine-governor Equivalent Model hmT(s) and the hydraulic turbine-governor Equivalent Model hmH(s), table
Sign is conventional
The adjustment effect of unit;
Step 5:Meter and load damping effect are established and contain H∑And hmWF(s) the equivalent rotor vacillation equation of electric system
Power grid
Frequency response models describe power system frequency response characteristic;
Step 6:It is constraint with power system frequency steady-state deviation boundary, finds out wind power penetration limit α1;
Step 7:It is constraint with power system frequency change rate boundary, finds out wind power penetration limit α2;
Step 8:Compare wind power penetration limit α1With wind power penetration limit α2, take α1And α2Smaller value is as electricity
Net
Wind power penetration limit.
Further, in step 1, if machine group includes N platform unit set G={ 1,2 ... j ... N }, and assume with group's unit
Rotor velocity is mutually all ω, then the parameter of equivalent unit is:
Wherein, subscript j indicates jth platform unit in a group of planes, is thermoelectricity, water power or Wind turbines;Sj、KjRespectively unit
The equivalent inertia time constant of capacity and unit.
Further, in step 2, if power plant includes the wind power plant and conventional power plant of the control containing virtual inertia in system,
Wind power plant group of planes quantity is a, and conventional power plant group of planes quantity is b, ifThe equivalent inertia time constant H of system∑It can table
It is shown as:
Wherein HeqWFi、SeqWFiThe respectively equivalent inertia time constant of wind power plant, rated capacity of the control containing virtual inertia,
HCONi、SCONiRespectively conventional power plant inertia time constant, rated capacity, KdfIncreasing is controlled for the virtual inertia of wind power plant inner blower
Benefit;H0Indicate the inertia size of conventional power unit, Δ H indicates the inertia size of wind-powered electricity generation virtual inertia response.
Further, in step 3, the dynamic response model transmission function of the primary frequency modulation control system of separate unit Wind turbines
hwt(s) it is:
Wherein b0, a0, a1, a2, a3For transmission function hwt(s) coefficient, KpfFor sagging control gain;Using the polymerization of formula (1)
Equivalent method, dynamic response model transmission function h of the wind power plant based on revolving speed primary frequency modulation control systemmWF(s) it is:
Wherein b0G, a0G, a1G, a2G, a3GRespectively transmission function hmWF(s) every equivalent parameters, if
Thus Δ ω is obtainedsWith Δ PmWFRelationship:
Further, in step 4, the mechanical output increment of equivalent steam turbine-governor is:
The mechanical output increment of the equivalent hydraulic turbine-governor is represented by:
Wherein RTG、RHG、TRHG、FHPG、TwGRespectively steam turbine difference coefficient, hydraulic turbine difference coefficient, reheater time are normal
Number, high pressure turbine stage power accounting, the equivalent polymerization parameter of water hammer effect coefficient.
Further, in step 5, separate unit generating set j is with rated capacity SjFor benchmark capacity, after polymerization with equal check-ins
Capacity SGIt waves equation for the equal check-ins of benchmark capacity and is:
Wherein Hj、Dj、PmjAnd PejThe respectively inertia time constant of motor j, damped coefficient, original machine power and electromagnetic work
Rate, each parameter value are with its own rated capacity SjFor the perunit value of base value;Unit in system is divided into two classes, conventional machine
The Wind turbines of group and the control of virtual inertia containing wind-powered electricity generation;Meter and load damping effect and wind power penetration limit α contain wind-powered electricity generation
Virtual inertia and the mains frequency response model of wind-powered electricity generation primary frequency modulation response effect are:
Wherein Δ PLFor system power vacancy, Δ PmWF、ΔPmT、ΔPmHThe active increasing of conventional Wind turbines after respectively disturbing
Amount, the active increment of fired power generating unit, the active increment of Hydropower Unit, D are system equivalent load damped coefficient.
Further, in step 6, when Operation of Electric Systems tends to stable state, d Δ ω/dt=0, corresponding frequency-domain result has Δ
ω (s) .s=0, convolution (9) can obtain:
Wherein, δ is steady frequency deviation boundary, the anti-wind power penetration limit for solving the deviation constraint containing steady frequency
α1。
Further, in step 7, power shortage generate, frequency disturbance initial stage, frequency departure very little, can approximation take Δ
ω=0, convolution (9) can obtain:
Wherein, η is frequency change rate boundary, obtains the wind power penetration limit α constrained containing frequency change rate2。
Further, in step 8, convolution (10) and (11) can obtain the response of wind-powered electricity generation virtual inertia and primary frequency modulation response
Wind power penetration limit α:
α=min (α1,α2) (12)
Beneficial effects of the present invention:
1) active increment and electric power that the present invention generates the response of true wind-powered electricity generation wind-powered electricity generation virtual inertia and primary frequency modulation response
System frequency constraint is introduced into the calculating of wind power penetration limit, improves the robustness of power grid under wind-electricity integration situation;
2) present invention objectively conducts in-depth analysis to the threshold value of the stochastic volatility of wind power output, improves wind power plant
The adaptability of programme;
3) a kind of wind power penetration limit calculation method for considering frequency constraint and wind-powered electricity generation frequency modulation proposed by the present invention can
To be widely used in the simulation analysis of the electric system containing wind-powered electricity generation, it is ensured that power network safety operation and wind-electricity integration planning construction have
Great importance.
Specific embodiment
As shown in Figure 1, a kind of wind power penetration limit calculation method for considering frequency constraint and wind-powered electricity generation frequency modulation, specific to wrap
Include following steps,
Step 1:The dynamic equivalent parameter aggregation method based on weighting of introducing, it is conventional to a typhoon motor group and b platform respectively
Machine group carries out parameter equivalent calculation;If learning that some group of planes (thermoelectricity, water power or wind-powered electricity generation) includes N platform by the identification of people having the same aspiration and interest method
Unit set G={ 1,2 ... j ... N }, and assume with group's unit rotor angular speed to be mutually all ω, then the parameter of equivalent unit is:
Wherein, subscript j, G are respectively jth platform unit and equal check-ins, S in a group of planesj、KjThe respectively capacity and unit of unit
Equivalent inertia time constant;
Step 2:Using wind-powered electricity generation virtual inertia when, by power plant in system be divided into containing virtual inertia control wind power plant and
Two class of conventional power plant (unit), and station quantity is respectively a, b, ifThe then system etc. under different wind-powered electricity generation permeabilities
Imitate inertia time constant H∑It is represented by:
Wherein HeqWFi、SeqWFiThe respectively equivalent inertia time constant of wind power plant, rated capacity of the control containing virtual inertia,
HCONi、SCONiRespectively conventional power plant inertia time constant, rated capacity, KdfIncreasing is controlled for the virtual inertia of wind power plant inner blower
Benefit;H0Indicate the inertia size of conventional power unit, Δ H indicates the inertia size of wind-powered electricity generation virtual inertia response, no wind-powered electricity generation virtual inertia
Δ H=0 when effect;
Step 3:In wind power plant primary frequency modulation response polymerization model, using the primary frequency modulation auxiliary control plan of revolving speed control
When slightly, the transmission function h that wind power plant primary frequency modulation responds equivalent polymerization model is solvedmWF(s) and wind power plant mechanical output increment Delta
PmWF(s):The dynamic response model transmission function h of the primary frequency modulation control system of separate unit Wind turbineswt(s) it is:
Wherein b0, a0, a1, a2, a3For transmission function hwt(s) coefficient, KpfFor sagging control gain;Using the equivalence of formula (1)
Method, dynamic response model transmission function h of the wind power plant based on revolving speed primary frequency modulation control systemmWF(s) it is:
Wherein b0G, a0G, a1G, a2G, a3GRespectively transmission function hmWF(s) every equivalent parameters, if
Thus Δ ω is obtainedsWith Δ PmWFRelationship:
Step 4:Establish steam turbine-governor Equivalent Model hmT(s) and the hydraulic turbine-governor Equivalent Model hmH(s), table
Levy the adjustment effect of conventional power unit;For steam turbine-governor model, using weighted equivalent method, by fired power generating unit in system etc.
Value is a machine, and the mechanical output increment of equivalent steam turbine-governor is:
Similarly, the mechanical output increment of the equivalent hydraulic turbine-governor is represented by:
Wherein RTG、RHG、TRHG、FHPG、TwGRespectively steam turbine difference coefficient, hydraulic turbine difference coefficient, reheater time are normal
Number, high pressure turbine stage power accounting, the equivalent polymerization parameter of water hammer effect coefficient;
Step 5:To separate unit generating set j, with rated capacity SjFor benchmark capacity, after polymerization with equal check-ins capacity SGFor
The equal check-ins of reference capacity wave equation and are:
Wherein Hj、Dj、PmjAnd PejThe respectively inertia time constant of motor j, damped coefficient, original machine power and electromagnetic work
Rate, each parameter value are with its own rated capacity SjFor the perunit value of base value;The unit in system is still divided into two classes, often
Advise the Wind turbines of unit and the control of virtual inertia containing wind-powered electricity generation;Meter and load damping effect and wind power penetration limit α, build
Vertical virtual inertia containing wind-powered electricity generation and the mains frequency response model of wind-powered electricity generation primary frequency modulation response effect are:
Wherein Δ PLFor system power vacancy, Δ PmWF、ΔPmT、ΔPmHThe active increasing of conventional Wind turbines after respectively disturbing
Amount, the active increment of fired power generating unit, the active increment of Hydropower Unit, D are system equivalent load damped coefficient;
Step 6:It is constraint with power system mesomeric state frequency departure, when Operation of Electric Systems tends to stable state, d Δ ω/dt=
0, corresponding frequency-domain result has Δ ω (s) .s=0, and convolution (9) can obtain:
Wherein, δ is steady frequency deviation boundary, counter can solve the wind power penetration limit of the deviation constraint containing steady frequency
α1;
Step 7:It is constraint with frequency change rate, in power shortage generation, frequency disturbance initial stage, frequency departure very little can
Approximation takes Δ ω=0, and convolution (9) can obtain:
Wherein, η is frequency change rate boundary, can obtain the wind power penetration limit α constrained containing frequency change rate2;
Step 8:It takes with the smaller of the penetration under steady frequency deviation constraint and frequency change rate boundary constraint
It is worth the penetration as the electric system under Frequency Index constraint, convolution (10) and (11) can obtain wind-powered electricity generation virtual inertia
The wind power penetration limit α of response and primary frequency modulation response:
α=min (α1,α2) (12)
A kind of embodiment penetrates power pole by the wind-powered electricity generation that simulation example verifies meter and virtual inertia and primary frequency modulation response
The accuracy of limit analytical Calculation result and different Kpf、KdfInfluence to wind power penetration limit.
As shown in Fig. 2, establishing analogue system under Matlab/simulink environment, two regions pass through two in system
Interconnection connection, region 1 are provided with a Hydropower Unit G2, a wind power plant and transformer T1, transformer T2, and region 2 is arranged
There are fired power generating unit G3, fired power generating unit G4 and transformer T3, transformer T4, load L1, L2, C1, C2 are respectively in two Area Interfaces
It is accessed at bus, load L3 is accessed and cut off to simulate the frequency of the analogue system power shortage as disturbance load, by L3
Accident.The response of wind-powered electricity generation virtual inertia and primary frequency modulation response, 8 formula of verification step are based on to the Wind turbines of wind power plant in Fig. 2
(12) it falls into a trap and virtual inertia/accuracy and difference of the wind power penetration limit analytical Calculation result of primary frequency modulation response
Kpf、KdfInfluence to wind power penetration limit.
Wherein simulation parameter is as follows:
With capacity 100MVA, on the basis of voltage 230kV, double-fed fan parameter:Voltage rating Vn=575V, rated power Pn
=1.5MW, stator resistance Rs=0.023pu, stator inductance Ls=0.18pu, rotor resistance Rr=0.016pu, inductor rotor Lr
=0.16pu, magnetizing inductance Lm=2.9pu, built-in inertia time constant HDFIG=5.29s, speed control integral coefficient Ki=
0.6.Rated angular velocity ωnom=157.08rad/s, rated wind speed VwN=11.7m/s, current transformer timeconstantτ=0.02s.
With capacity 100MVA, on the basis of voltage 230kV, Hydropower Unit G2, fired power generating unit G3, fired power generating unit G4 parameter:Volume
Constant volume Sn=900MVA, voltage rating Un=20kV, d axle inductance Xd=1.8pu, q axle inductance Xq=1.7pu, Xa=0.2pu, d
Axis transient inductance Xd'=0.3pu, q axis transient inductance Xq'=0.55pu, d axis time transient inductance Xd"=0.25pu, q axis time transient state
Inductance Xq"=0.25pu, Ra=0.0025pu, d axis transient state time Td0'=8.0s, q axis transient state time Tq0'=0.4s, d axis
Transient state time Td0"=0.03s, q axis time transient state time Tq0"=0.05s, Hydropower Unit G2 inertia time constant H=6.5s, fire
Motor group G3, fired power generating unit G4 inertia time constant H=6.175s.
With capacity 100MVA, on the basis of voltage 230kV, transformer T1, transformer T2, transformer T3, transformer T4 parameter:
Rated capacity Sn=900MVA, voltage change ratio Un1/Un2=20Kv/230kV, impedance Rt+jXt=0+j0.15pu.
With capacity 100MVA, on the basis of voltage 230kV, transmission line parameter:Resistance RL=0.0001pu/km, inductance XL=
0.001pu/km, conductance BC=0.00175pu/km.
With capacity 100MVA, on the basis of voltage 230kV, load data:L1 active-power PL1=800MW, L1 reactive power
QL1=100MVAR, C1 reactive power QC1=-187MVAR, C2 reactive power QC2=-200MVAR, L2 active-power PL2=
800MW, L2 reactive power QL2=100MVAR, additional load active-power PL3=160MW.
In verification process, emulation project includes:(a) initial wind speed Vw=10m/s, Kdf=0, different K are setpf, different
Primary frequency modulation responds under function and effect, calculates wind power penetration limit by step 8 formula (12), and set wind with calculated result
Electric grid connection capacity is emulated by the frequency response under power shortage, verifies steady-state deviation in turn, whether frequency change rate approaches
Defined steady frequency deviation constraint value and frequency change rate binding occurrence;(b) initial wind speed VwDifferent K are arranged in=10m/sdf, no
Under same wind-powered electricity generation virtual inertia response function and effect, passes through step 8 formula (12) and calculate wind power penetration limit.Subsequent process with
Emulation project (a) is identical;(c) initial wind speed V is setwDifferent K are arranged in=10m/spfAnd Kdf, reflect different wind-powered electricity generation virtual inertias
Coupled in common function and effect are responded with primary frequency modulation, according to (a) emulation Project Process verifying.With reference to national standard and Europe Entso-e
Technical stipulation, is arranged steady frequency deviation δ=± 0.2Hz, and frequency change rate boundary value is η=± 0.005Hz/s.
As shown in figure 3, obtaining Kdf=0, KpfFrequency departure curve when=0,1 times, 2 times, more every frequency response
Value, as shown in table 1, lists 4 kinds of situations in table 1.According to simulation result it can be found that:1) 1 → situation of comparable situation 3, power grid
Steady frequency deviation is very close defined | δ |=0.2Hz, and be limited in the restrained boundary.Thus illustrate that frequency refers to
Under mark constraint, the correctness of the calculation method of the wind power penetration limit of the response of primary frequency modulation containing wind-powered electricity generation;2) from 1 → feelings of situation
Condition 3, the sagging control gain of primary frequency modulation is bigger, and wind power penetration limit is higher, KpfWhen by 0 → 2 times, wind-powered electricity generation penetrates power pole
Limiting scaling up is about 2%, and steady frequency deviation is still limited near binding side dividing value at the same time.This illustrates wind
Electric primary frequency modulation action intensity influences significantly, if not considering wind-powered electricity generation primary frequency modulation function and effect wind power penetration limit
(Kpf=0) result will be made relatively conservative, limit the grid-connected scale of wind-powered electricity generation;3) comparable situation 3 and situation 4, the two is to electricity
Wind turbines are applied with primary frequency modulation control in net, and operation of power networks condition is identical.But situation 3 is calculated according to step 8 formula (12)
α, and meter and the effect of wind-powered electricity generation primary frequency modulation during calculating α of situation 4 make to calculate two and different penetrate power pole
Limit, respectively 22.57% and 20.55%, difference 2%.It can be seen that in the practical electricity for applying the auxiliary control of wind-powered electricity generation primary frequency modulation
In Force system, if not including the active response model of primary frequency modulation in frequency response computation model, the α value of calculating is significantly less than normal
(deviateing true value), tend to be conservative in instructing wind-powered electricity generation planning construction.
Difference K when 1 sudden increase in load of tablepfFrequency response Indexes Comparison
As shown in figure 4, obtaining Kpf=0, KdfFrequency change rate curve when=0,1 times, 2 times, the emulation of corresponding diagram 4
Journey, more every frequency response values list 4 kinds of situations as shown in table 2 in table 2.According to simulation result it can be found that:1) from
1 → situation of situation 3, frequency change rate are very close defined | η |=0.005Hz/s, and be limited in the restrained boundary.
Thus illustrate Frequency Index constraint under, virtual inertia containing wind-powered electricity generation response wind power penetration limit calculation method it is correct
Property;2) 1 → situation of comparable situation 3, wind power penetration limit scaling up are about 3%.This illustrates that wind-powered electricity generation virtual inertia is rung
Action intensity is answered to influence wind power penetration limit significant, i.e., with KdfIncrease, the wind power penetration limit of system
It is bigger;3) comparable situation 3 and situation 4, penetration difference calculate two kind in the case of similar with emulation project (a)
For 22.84% and 20.13%, difference about 2.5%.Therefore, in the practical power systems for applying the control of wind-powered electricity generation virtual inertia, if
The Δ H increment of reflection wind-powered electricity generation virtual inertia response is not included in frequency response computation model, the α value of calculating is significant less than normal (inclined
From true value), it will also tend to be conservative in instructing wind-powered electricity generation planning construction.
Difference K when 2 sudden increase in load of tabledfUnder comparison result
As shown in Figure 5, Figure 6, K when sudden increase in loadpf=1 difference KdfSystem frequency deviation response and system frequency under coupling
The emulation of rate change rate response, more every frequency response index list 4 kinds of situations as shown in table 3 in table 3.According to emulation
As a result known to:1) from 1 → situation of situation 3, frequency change rate maximum value, steady frequency deviation is very close to defined | δ |=
0.2Hz, | η |=0.005Hz/s, and be limited in the restrained boundary.Thus illustrate to calculate under Frequency Index constraint
Wind power penetration limit α it is more accurate;2) 1 → situation of situation 3 in comparison sheet 1, table 2 and table 3, can must comprehensively consider used
Wind power penetration limit under property and primary frequency modulation coupling penetrates power pole than containing only the wind-powered electricity generation under a frequency modulation control
Limit increases;3) comparable situation 3 and situation 4, it is similar with emulation project (a), power pole is penetrated in the case of calculating two kinds
Limit is respectively 23.83% and 21.24%, difference 2.6%.Thus, if in wind power penetration limit computation model not
Meter and the control of wind-powered electricity generation virtual inertia and primary frequency modulation control action, the α value acquired are (deviateing true value) less than normal.
K when 3 sudden increase in load of tablepfAnd KdfComparison result under coupling