CN108832658A - A kind of wind power penetration limit calculation method considering frequency constraint and wind-powered electricity generation frequency modulation - Google Patents

A kind of wind power penetration limit calculation method considering frequency constraint and wind-powered electricity generation frequency modulation Download PDF

Info

Publication number
CN108832658A
CN108832658A CN201810652509.2A CN201810652509A CN108832658A CN 108832658 A CN108832658 A CN 108832658A CN 201810652509 A CN201810652509 A CN 201810652509A CN 108832658 A CN108832658 A CN 108832658A
Authority
CN
China
Prior art keywords
wind
electricity generation
powered electricity
frequency
power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810652509.2A
Other languages
Chinese (zh)
Other versions
CN108832658B (en
Inventor
李世春
唐红艳
刘道兵
钟浩
舒征宇
雷小林
吕翔生
王扬
曹润杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN201810652509.2A priority Critical patent/CN108832658B/en
Publication of CN108832658A publication Critical patent/CN108832658A/en
Application granted granted Critical
Publication of CN108832658B publication Critical patent/CN108832658B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)
  • Wind Motors (AREA)

Abstract

A kind of wind power penetration limit calculation method considering frequency constraint and wind-powered electricity generation frequency modulation acquires the active increment of wind power plant primary frequency modulation response transfer function model and wind-powered electricity generation primary frequency modulation using weighted equivalent polymerization;And consider to be based on wind-powered electricity generation virtual inertia control strategy, calculating has solved the equivalent inertia time constant of electric system under different wind-powered electricity generation permeabilities;Pass through meter and load damping effect, the mains frequency response model of equivalent rotor vacillation equation is established to describe power system frequency response characteristic, and the calculation method of power grid wind penetration is proposed to constrain with power system frequency steady-state deviation boundary and frequency change rate boundary.The active increment that the response of true wind-powered electricity generation wind-powered electricity generation virtual inertia and primary frequency modulation response generate and power system frequency constraint are introduced into the calculating of wind power penetration limit by the present invention, the robustness for improving power grid under wind-electricity integration situation, to ensuring that power network safety operation and wind-electricity integration planning construction have great importance.

Description

A kind of wind power penetration limit calculation method considering frequency constraint and wind-powered electricity generation frequency modulation
Technical field
The invention belongs to wind-electricity integration control fields, and in particular to a kind of wind-powered electricity generation for considering frequency constraint and wind-powered electricity generation frequency modulation is worn Saturating power limit calculation method.
Background technique
As wind-powered electricity generation permeability is continuously improved, variable-speed wind-power unit substitutes more and more conventional power units and occupies in power grid More high specific gravity weakens the decline of electric system equivalent inertia, Primary frequency control ability, causes system response to deteriorate, limit Receiving ability of the power grid to wind-powered electricity generation.Wind power penetration limit refers to electric system institute's energy in the case where meeting various operation constraint conditions The wind power plant maximum installed capacity of receiving accounts for the ratio of system peak load.Wind-powered electricity generation virtual inertia control technology and primary frequency modulation are auxiliary Helping control technology is to improve system frequency to fall effective means with wind power integration capacity.Therefore, in speed-changing draught fan virtual inertia How control technology and the application of primary frequency modulation control technology popularity under the new situation, calculate wind-powered electricity generation to accurate quantitative analysis and penetrate function The rate limit, for ensuring that power network safety operation and wind-electricity integration planning construction have great importance.
About wind power penetration limit calculation method there are many.When not considering wind-powered electricity generation inertia/primary frequency modulation, Traditional wind power penetration limit calculation method is the penetration for assuming first that wind-powered electricity generation, and it is imitative then to carry out repeatedly dynamic True experiment, and assumed value is constantly corrected, until meeting restriction condition, this method computational efficiency and precision are low.Scholar is again A kind of wind power penetration limit calculation method based on mathematical optimization is proposed, to restrict the factor of wind power integration as constraint item Part optimizes calculating using wind power integration maximum capacity as optimization aim, solves wind power penetration limit.In wind-powered electricity generation The new situations that permeability continues to increase, speed-changing draught fan virtual inertia control technology and primary frequency modulation control technology popularity are applied Under, calculating the wind power penetration limit controlled containing virtual inertia and primary frequency modulation is necessary.According to not considering Wind-powered electricity generation inertia/primary frequency modulation method calculates the wind power penetration limit of the control containing virtual inertia/primary frequency modulation, it will makes As a result relatively conservative, the grid-connected scale of wind-powered electricity generation is limited, the wind power penetration limit not being suitable under the new situation calculates.From mesh From the point of view of the preceding constraint condition that research is calculated wind power penetration limit both at home and abroad, frequency steady state constraint and frequency are not considered temporarily Transient state constraint combines to calculate penetration.
Therefore, the wind-powered electricity generation of the response of wind-powered electricity generation virtual inertia and primary frequency modulation response of research meter and Frequency Index constraint penetrates function Rate Limit Calculation Method.
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 HAnd 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 systemIt 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 (α12) (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.
Detailed description of the invention
Present invention will be further explained below with reference to the attached drawings and examples.
Fig. 1 is the flow chart for considering the wind power penetration limit calculation method of frequency constraint and wind-powered electricity generation frequency modulation.
Fig. 2 is the analogue system schematic diagram of the embodiment of the present invention.
K when Fig. 3 is sudden increase in loaddf=0 difference KpfUnder system frequency deviation response curve.
K when Fig. 4 is sudden increase in loadpf=0 difference KdfUnder system frequency change rate response curve.
K when Fig. 5 is sudden increase in loadpf=1 times of difference KdfSystem frequency deviation response curve under coupling.
K when Fig. 6 is sudden increase in loadpf=1 times of difference KdfSystem frequency change rate response curve under coupling.
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 HIt 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 (α12) (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

Claims (9)

1. a kind of wind power penetration limit calculation method for considering frequency constraint and wind-powered electricity generation frequency modulation, which is characterized in that specific packet Include following steps,
Step 1:Using the dynamic equivalent parameter aggregation method of weighting, respectively to wind turbine group and thermal motor group, hydroelectric machine 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, calculates the equivalent inertial time of system under different wind-powered electricity generation permeabilities Between 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), Using weighting dynamic equivalent parameter aggregation method, calculates and solve dynamic response of the wind power plant based on revolving speed primary frequency modulation control system 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), characterization is conventional The adjustment effect of unit;
Step 5:Meter and load damping effect are established and contain HAnd hmWF(s) power grid of the equivalent rotor vacillation equation of electric system 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 α2Wind of the smaller value as power grid Electric penetration.
2. the wind power penetration limit calculation method according to claim 1 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, in step 1, if machine group includes N platform unit set G={ 1,2 ... j ... N }, and assumes with group's unit rotor angle speed Degree 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、KjThe respectively capacity of unit With the equivalent inertia time constant of unit.
3. the wind power penetration limit calculation method according to claim 2 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, in step 2, if power plant includes the wind power plant and conventional power plant of the control containing virtual inertia, wind power plant machine in system Group's quantity is a, and conventional power plant group of planes quantity is b, ifThe equivalent inertia time constant H of systemIt 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, KdfGain is controlled for the virtual inertia of wind power plant inner blower;H0 Indicate the inertia size of conventional power unit, Δ H indicates the inertia size of wind-powered electricity generation virtual inertia response.
4. the wind power penetration limit calculation method according to claim 2 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, in step 3, 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;It is equivalent using the polymerization 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, ifThus Obtain Δ ωsWith Δ PmWFRelationship:
5. the wind power penetration limit calculation method according to claim 2 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, 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 constant, High pressure turbine stage power accounting, the equivalent polymerization parameter of water hammer effect coefficient.
6. the wind power penetration limit calculation method according to claim 2 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, in step 5, separate unit generating set j is 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 power, Each parameter value is with its own rated capacity SjFor the perunit value of base value;
Unit in system is divided into two classes, the Wind turbines of conventional power unit and the control of virtual inertia containing wind-powered electricity generation;Meter and load resistance The mains frequency response of Buddhist nun's effect and wind power penetration limit α, virtual inertia containing wind-powered electricity generation and wind-powered electricity generation primary frequency modulation response effect Model is:
Wherein Δ PLFor system power vacancy, Δ PmWF、ΔPmT、ΔPmHThe active increment of routine Wind turbines, fire after respectively disturbing The active increment of motor group, the active increment of Hydropower Unit, D are system equivalent load damped coefficient.
7. the wind power penetration limit calculation method according to claim 6 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, 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 frequency1
8. the wind power penetration limit calculation method according to claim 7 for considering frequency constraint and wind-powered electricity generation frequency modulation, Be characterized in that, in step 7, power shortage generate, frequency disturbance initial stage, frequency departure very little, can approximation take Δ ω=0, tie Box-like (9) can obtain:
Wherein, η is frequency change rate boundary, obtains the wind power penetration limit α constrained containing frequency change rate2
9. the wind power penetration limit calculation method according to claim 8 for considering frequency constraint and wind-powered electricity generation frequency modulation, It is characterized in that, in step 8, convolution (10) and (11) can obtain the wind-powered electricity generation that the response of wind-powered electricity generation virtual inertia is responded with primary frequency modulation and penetrate Power limit α:
α=min (α12) (12)。
CN201810652509.2A 2018-06-22 2018-06-22 Wind power penetration power limit calculation method considering frequency constraint and wind power frequency modulation Active CN108832658B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810652509.2A CN108832658B (en) 2018-06-22 2018-06-22 Wind power penetration power limit calculation method considering frequency constraint and wind power frequency modulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810652509.2A CN108832658B (en) 2018-06-22 2018-06-22 Wind power penetration power limit calculation method considering frequency constraint and wind power frequency modulation

Publications (2)

Publication Number Publication Date
CN108832658A true CN108832658A (en) 2018-11-16
CN108832658B CN108832658B (en) 2021-06-04

Family

ID=64138044

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810652509.2A Active CN108832658B (en) 2018-06-22 2018-06-22 Wind power penetration power limit calculation method considering frequency constraint and wind power frequency modulation

Country Status (1)

Country Link
CN (1) CN108832658B (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109599865A (en) * 2018-12-13 2019-04-09 国网新疆电力有限公司 A kind of method of new energy accounting in determining electric system
CN110061521A (en) * 2019-04-30 2019-07-26 东南大学 A kind of maximum wind permeability fast evaluation method considering frequency accumulation effect
CN110175933A (en) * 2019-05-29 2019-08-27 哈尔滨工业大学 Direct-driving type wind power plant Dynamic Equivalence suitable for frequency modulation control
CN110445191A (en) * 2019-08-02 2019-11-12 浙江大学 With the electric system generator group frequency modulation non-firm power distribution method of support vector machines
CN110518632A (en) * 2019-08-05 2019-11-29 三峡大学 A kind of wind farm grid-connected quantitative calculation method that power grid inertia is weakened
CN110750882A (en) * 2019-09-29 2020-02-04 武汉大学 Wind power ratio limit value analytical calculation method considering frequency constraint
CN110829487A (en) * 2019-10-22 2020-02-21 电子科技大学 Dynamic frequency prediction method for power system
CN111864813A (en) * 2020-06-23 2020-10-30 国网辽宁省电力有限公司电力科学研究院 Wind/thermal power combined frequency control method based on virtual weight coefficient
CN111934328A (en) * 2020-07-03 2020-11-13 国家电网有限公司 Energy storage self-adaptive control method for coping with wind power limit scene
CN112436504A (en) * 2019-08-26 2021-03-02 南京理工大学 Wind power penetration power limit analysis method considering system frequency secondary drop
CN112531742A (en) * 2020-12-07 2021-03-19 江苏方天电力技术有限公司 Method for keeping low-inertia power grid frequency stable based on deep peak shaving state
CN112821419A (en) * 2021-01-21 2021-05-18 中国三峡新能源(集团)股份有限公司 Power grid frequency modulation limit analysis method and system
CN112838621A (en) * 2021-01-22 2021-05-25 上海交通大学 Electric power system frequency modulation capacity realization method considering new energy growth
CN112886586A (en) * 2021-03-10 2021-06-01 东北电力大学 Method for optimizing equivalent model power recovery control parameters of doubly-fed wind power plant
CN112994013A (en) * 2021-02-08 2021-06-18 四川大学 Multi-source power system day-ahead optimization scheduling method considering frequency constraints
CN113065237A (en) * 2021-03-19 2021-07-02 四川英杰电气股份有限公司 Method for automatically setting frequency modulation boundary and radio frequency power supply
CN113098030A (en) * 2021-04-07 2021-07-09 王巧 Frequency control optimization method for wind power plant participating in primary frequency modulation
CN113659639A (en) * 2021-08-13 2021-11-16 云南电网有限责任公司电力科学研究院 Wind power plant inertia response power distribution method considering rotating speed constraint
CN113746134A (en) * 2021-08-20 2021-12-03 华北电力大学(保定) Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit
CN116706944A (en) * 2023-08-07 2023-09-05 昆明理工大学 Method for determining new energy permeability under frequency safety constraint
CN113746134B (en) * 2021-08-20 2024-06-28 华北电力大学(保定) Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104167765A (en) * 2014-07-11 2014-11-26 海南电网公司 Admitting ability distribution-based maximum wind power installed capacity calculation method
CN105281358A (en) * 2014-05-30 2016-01-27 清华大学 Wind power limit grid-connected capacity calculating method under constraint of frequency modulation and peak-load regulation adequacy
CN105610159A (en) * 2016-03-02 2016-05-25 东北电力大学 Method for analyzing wind farm-provided primary frequency modulation reserve capacity required for power system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105281358A (en) * 2014-05-30 2016-01-27 清华大学 Wind power limit grid-connected capacity calculating method under constraint of frequency modulation and peak-load regulation adequacy
CN104167765A (en) * 2014-07-11 2014-11-26 海南电网公司 Admitting ability distribution-based maximum wind power installed capacity calculation method
CN105610159A (en) * 2016-03-02 2016-05-25 东北电力大学 Method for analyzing wind farm-provided primary frequency modulation reserve capacity required for power system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUAN FU等: "Integrated wind turbine controller with virtual inertia and primary frequency responses for grid dynamic frequency support", 《IET RENEWABLE POWER GENERATION》 *
付媛等: "变速风电机组的惯性与一次调频特性分析及综合控制", 《中国电机工程学报》 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109599865A (en) * 2018-12-13 2019-04-09 国网新疆电力有限公司 A kind of method of new energy accounting in determining electric system
CN109599865B (en) * 2018-12-13 2022-06-24 国网新疆电力有限公司 Method for determining new energy ratio in power system
CN110061521A (en) * 2019-04-30 2019-07-26 东南大学 A kind of maximum wind permeability fast evaluation method considering frequency accumulation effect
CN110175933B (en) * 2019-05-29 2021-06-22 哈尔滨工业大学 Direct-drive wind power plant dynamic equivalence method suitable for frequency modulation control
CN110175933A (en) * 2019-05-29 2019-08-27 哈尔滨工业大学 Direct-driving type wind power plant Dynamic Equivalence suitable for frequency modulation control
CN110445191A (en) * 2019-08-02 2019-11-12 浙江大学 With the electric system generator group frequency modulation non-firm power distribution method of support vector machines
CN110518632A (en) * 2019-08-05 2019-11-29 三峡大学 A kind of wind farm grid-connected quantitative calculation method that power grid inertia is weakened
CN110518632B (en) * 2019-08-05 2022-10-21 三峡大学 Quantitative calculation method for weakening power grid inertia by wind power plant grid connection
CN112436504A (en) * 2019-08-26 2021-03-02 南京理工大学 Wind power penetration power limit analysis method considering system frequency secondary drop
CN110750882A (en) * 2019-09-29 2020-02-04 武汉大学 Wind power ratio limit value analytical calculation method considering frequency constraint
CN110750882B (en) * 2019-09-29 2022-09-02 武汉大学 Wind power ratio limit value analytical calculation method considering frequency constraint
CN110829487A (en) * 2019-10-22 2020-02-21 电子科技大学 Dynamic frequency prediction method for power system
CN111864813A (en) * 2020-06-23 2020-10-30 国网辽宁省电力有限公司电力科学研究院 Wind/thermal power combined frequency control method based on virtual weight coefficient
CN111864813B (en) * 2020-06-23 2022-11-11 国网辽宁省电力有限公司电力科学研究院 Wind/thermal power combined frequency control method based on virtual weight coefficient
CN111934328B (en) * 2020-07-03 2022-05-06 国家电网有限公司 Energy storage self-adaptive control method for coping with wind power limit scene
CN111934328A (en) * 2020-07-03 2020-11-13 国家电网有限公司 Energy storage self-adaptive control method for coping with wind power limit scene
CN112531742A (en) * 2020-12-07 2021-03-19 江苏方天电力技术有限公司 Method for keeping low-inertia power grid frequency stable based on deep peak shaving state
CN112821419A (en) * 2021-01-21 2021-05-18 中国三峡新能源(集团)股份有限公司 Power grid frequency modulation limit analysis method and system
CN112838621A (en) * 2021-01-22 2021-05-25 上海交通大学 Electric power system frequency modulation capacity realization method considering new energy growth
CN112994013A (en) * 2021-02-08 2021-06-18 四川大学 Multi-source power system day-ahead optimization scheduling method considering frequency constraints
CN112886586A (en) * 2021-03-10 2021-06-01 东北电力大学 Method for optimizing equivalent model power recovery control parameters of doubly-fed wind power plant
CN113065237B (en) * 2021-03-19 2022-11-08 四川英杰电气股份有限公司 Method for automatically setting frequency modulation boundary and radio frequency power supply
CN113065237A (en) * 2021-03-19 2021-07-02 四川英杰电气股份有限公司 Method for automatically setting frequency modulation boundary and radio frequency power supply
CN113098030A (en) * 2021-04-07 2021-07-09 王巧 Frequency control optimization method for wind power plant participating in primary frequency modulation
CN113098030B (en) * 2021-04-07 2023-06-16 王巧 Frequency control optimization method for wind power plant participating in primary frequency modulation
CN113659639A (en) * 2021-08-13 2021-11-16 云南电网有限责任公司电力科学研究院 Wind power plant inertia response power distribution method considering rotating speed constraint
CN113659639B (en) * 2021-08-13 2023-11-21 云南电网有限责任公司电力科学研究院 Wind power plant inertia response power distribution method considering rotation speed constraint
CN113746134A (en) * 2021-08-20 2021-12-03 华北电力大学(保定) Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit
CN113746134B (en) * 2021-08-20 2024-06-28 华北电力大学(保定) Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit
CN116706944A (en) * 2023-08-07 2023-09-05 昆明理工大学 Method for determining new energy permeability under frequency safety constraint
CN116706944B (en) * 2023-08-07 2023-10-27 昆明理工大学 Method for determining new energy permeability under frequency safety constraint

Also Published As

Publication number Publication date
CN108832658B (en) 2021-06-04

Similar Documents

Publication Publication Date Title
CN108832658A (en) A kind of wind power penetration limit calculation method considering frequency constraint and wind-powered electricity generation frequency modulation
CN106910142B (en) Method for calculating frequency characteristic of power system with wind power active power-frequency coupling effect
CN106849088B (en) A kind of mains frequency characteristic computing method based on the control of wind-powered electricity generation pitch primary frequency modulation
CN109560573B (en) Method and device for optimizing frequency controller parameters of variable-speed wind turbine generator
CN103558768B (en) A kind of based on the equivalent modeling method of wind speed distribution characteristics in wind energy turbine set
CN102738809A (en) Optimized control method for wind power field reactive power compensation capacity considering wind power distribution rule
CN106786759B (en) The Wind turbines primary frequency modulation method mutually coordinated with Synchronous generator
CN106712058A (en) Coordinated control method for wind power plant of double-fed fan to participate in one-time frequency modulation of power system
CN102594244A (en) Joint control method of primary frequency modulation for doubly-fed wind power generation set
CN104283210A (en) Two-stage distributed simulation method for wind electricity collection area based on distributed load flow calculation
CN109038613A (en) A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds
CN110206686A (en) A kind of adaptive maximum power tracking and controlling method for wind power generating set
CN107895955A (en) A kind of cooperative control method of wind-powered electricity generation compensation hydraulic turbine water hammer effect
Margaris et al. Impact of wind power in autonomous power systems—power fluctuations—modelling and control issues
CN102769299B (en) Wind farm voltage control method based on voltage operating state
CN105186511A (en) Method for controlling participation of battery energy storage system in power grid secondary frequency modulation
CN102034008B (en) Variable-parameter capacitance compensation equivalent modeling method for wind power station of constant-speed wind power unit at fluctuating air speed
CN102904266A (en) Method for determining inactive compensation capacity network adaptability of wind power plant
CN102769292B (en) Method for enhancing voltage safety margin of wind farm
CN106712032A (en) Optimal power flow model construction method considering active power voltage regulation capacity of wind turbine generator set
Sanchez et al. Dynamic model of wind energy conversion systems with variable speed synchronous generator and full-size power converter for large-scale power system stability studies
CN113131526A (en) Static stability control method for wind-fire bundling system with virtual inertia control
Xueguang et al. Application of models of the wind energy conversion system to wind power dynamic analysis
CN104767205B (en) Method for establishing automatic generation control system of electric power system based on wind power plant access
CN104037815B (en) Eliminate the photovoltaic plant active power autocontrol method that electrical grid transmission power is out-of-limit

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20181116

Assignee: Hubei Yispatiotemporal Unmanned Aerial Vehicle Technology Co.,Ltd.

Assignor: CHINA THREE GORGES University

Contract record no.: X2023980044218

Denomination of invention: Calculation method for wind power penetration power limit considering frequency constraints and wind power frequency regulation

Granted publication date: 20210604

License type: Common License

Record date: 20231024