CN109038613A - A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds - Google Patents
A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds Download PDFInfo
- Publication number
- CN109038613A CN109038613A CN201810719653.3A CN201810719653A CN109038613A CN 109038613 A CN109038613 A CN 109038613A CN 201810719653 A CN201810719653 A CN 201810719653A CN 109038613 A CN109038613 A CN 109038613A
- Authority
- CN
- China
- Prior art keywords
- wind
- frequency modulation
- load
- inertia
- primary frequency
- 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.)
- Pending
Links
- 230000005611 electricity Effects 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 18
- 230000004044 response Effects 0.000 claims abstract description 38
- 230000005540 biological transmission Effects 0.000 claims abstract description 17
- 230000008859 change Effects 0.000 claims abstract description 11
- 238000011217 control strategy Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000006116 polymerization reaction Methods 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 238000007665 sagging Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims 1
- 240000002853 Nelumbo nucifera Species 0.000 claims 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims 1
- 230000009467 reduction Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Abstract
The present invention provides a kind of meter and the adaptive low frequency deloading methods of wind-powered electricity generation virtual inertia/primary frequency modulation response.Wind turbines apply under typical virtual inertia control policy condition, the Analytical Solution virtual inertia of wind-powered electricity generation and the equivalent inertia time constant (H of electric system∑).According to required real-time H∑, in conjunction with initial system frequency change rate, calculate system accurately total active power shortage.Simultaneously when Wind turbines apply typical primary frequency modulation control strategy, the active increment of wind power plant primary frequency modulation response transmission function and primary frequency modulation response is solved, the power shortage operated during low-frequency load shedding by wheel is calculated in conjunction with the virtual inertia dynamic of wind-powered electricity generation of solution, and optimizes off-load amount accordingly.Sample calculation analysis the result shows that, because fully considering wind-powered electricity generation virtual inertia and primary frequency modulation effect, electric system to the frequency modulation containing wind-powered electricity generation, mentioned method can more acurrate computing system power vacancy, each coherence low-frequency load shedding of more precise manipulation, frequency retrieval dynamic process and stable state more meet truth.
Description
Technical field
The present invention relates to technical field of electric power, and in particular to a kind of meter and wind-powered electricity generation virtual inertia/primary frequency modulation response
Adaptive low frequency deloading method.
Background technique
When system unit failure or load increase suddenly, power shortage is generated, system frequency reduces.Especially work as system
When large disturbances occur, it cannot only inhibit the continuous decrease of frequency by the frequency adjustment effect of itself, it may under serious situation
Lead to system crash and large area blackout.Last line of defense of the low-frequency load shedding as system stable operation, should first cut off
Sub-load prevents uncontrollable chain reaction to restore frequency stabilization as early as possible, causes a wide range of, prolonged power failure.
In recent years, installed capacity of wind-driven power sustainable growth, the permeability of wind-powered electricity generation is higher and higher in electric system.And speed-changing draught fan
Rotor speed and mains frequency decouple so that Wind turbines do not have inertial response and auxiliary fm capacity, increase power grid frequency
Rate adjusts pressure.Therefore, more factors are being had to take into account that the research of low-frequency load shedding under the new situation, wind-powered electricity generation can be efficiently used
Unit participates in frequency modulation and is of great importance for system frequency stabilization.
During the low-frequency load shedding after considering that wind-powered electricity generation participates in frequency modulation, for being applied with the wind turbine of virtual inertia control
Group, equivalent inertia time constant are generally difficult to parse.For now, researcher is usually to be set to a definite value,
There is certain error with actual value, affects the accuracy of acquired power shortage.If equivalent inertia time constant can be acquired at any time
Between the dynamic value that changes, be of great importance to the correction update of power shortage during off-load.In addition, currently for Wind turbines
Low-frequency load shedding research seldom consider primary frequency modulation influence, and wind-powered electricity generation participate in primary frequency modulation have great role to frequency stabilization.
Therefore, the low-frequency load shedding after virtual inertia control is combined with primary frequency modulation response is studied under high wind-powered electricity generation permeability for frequency
Stable operation is of great importance.
Summary of the invention
Using large-scale wind power access electric system as object, for the aforementioned drawback of the low-frequency load shedding scheme prior art, originally
Invention provides a kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds, can more acurrate calculating system
System power shortage, each coherence low-frequency load shedding of more precise manipulation, frequency retrieval dynamic process and stable state more meet truth.
The technical scheme adopted by the invention is as follows:
A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds, Wind turbines apply allusion quotation
Under the conditions of type virtual inertia control strategy, the Analytical Solution virtual inertia of wind-powered electricity generation and the equivalent inertia time constant of electric system
(H∑), according to required real-time H∑, in conjunction with initial system frequency change rate, calculate system accurately total active power shortage.Together
When Wind turbines apply typical primary frequency modulation control strategy, solved wind power plant primary frequency modulation response transmission function and primary
The active increment of frequency modulation response calculates the power operated during low-frequency load shedding by wheel in conjunction with the virtual inertia dynamic of wind-powered electricity generation of solution
Vacancy, and optimize off-load amount accordingly.
A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds, comprising the following steps:
Step 1: to containing be applied with virtual inertia control Wind turbines the equivalent inertia time constant of electric system into
Row analytical Calculation；
Step 2: according to the equivalent inertia time constant of electric system of the virtual inertia containing wind-powered electricity generation acquired, calculating electric system
Active power shortage；
Step 3: being asked according to the primary frequency modulation transfer function model of separate unit Wind turbines using weighted equivalent polymerization
Obtain the active increment of wind power plant primary frequency modulation；
Step 4: doing off-load meter according to the initial active power shortage of electric system and draw；
Step 5: according to the significance level of load, determining off-load place；
Step 6: according to the equivalent inertia time constant of accurate electric system and the active increment of primary frequency modulation, by wheel
Secondary ground dynamic optimization off-load amount.
In the step 1, according to the definition of inertia time constant, the H of separate unit blowerequAre as follows:
In formula, ωnom, ωs0, ωr0, HDFIG,Tf,KdfRespectively rated angular velocity, system initial synchronisation angular speed, initially
Rotor velocity, double-fed blower built-in inertia time constant, time constant filter, inertia control gain.KpT,KiTFor speed control
Device parameter processed.
It is a unit by wind power plant equivalence when the wind power plant DFIG blower of platform containing p, then wind power plant polymerization inertia time is normal
Number HeqWFiAre as follows:
Wherein,For the mean speed of i-th of wind power plant inner blower, CHi=ωs0HDFIGi/ω2 nomi。
Power plant in system is divided into two class of wind power plant of conventional power unit, the control containing virtual inertia, and station quantity is distinguished
For a, b；
The then equivalent inertia time constant H of the step 1∑Expression formula are as follows:
In formula: H0, HeqWFi, Δ H is respectively the inertia time constant of conventional power plant, the control wind power plant containing virtual inertia
System equivalent inertia time constant increment when equivalent inertia time constant and consideration wind-powered electricity generation virtual inertia response, SCONi, SeqWFi
The respectively rated capacity of conventional power unit controls wind power plant containing virtual inertia.
In the step 2, the expression formula of electric system active power shortage are as follows:
Wherein, H∑For the equivalent inertia time constant of system, system frequency change rate when df/dt is disturbance.
In the step 3, primary frequency modulation transfer function model assists control strategy using the primary frequency modulation of revolving speed control
When, solve the transmission function h that wind power plant primary frequency modulation responds equivalent polymerization modelmWF(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) are as follows:
Wherein b0, a0, a1, a2, a3For hwt(s) transfer-function coefficient, KpfFor sagging control gain；
After polymerization, dynamic response model transmission function h of the wind power plant based on revolving speed primary frequency modulation control systemmWF(s) are as follows:
Wherein b0G, a0G, a1G, a2G, a3GRespectively transmission function hmWF(s) every equivalent parameters.
It is hereby achieved that the active support that primary frequency modulation provides:
Through reverse drawing Laplace transform, time-domain expression is obtained:
L-1[△PmWF(s)]=△ PmWF(t)
ΔPmWFIt (s) is the active increment of wind-powered electricity generation primary frequency modulation.
In the step 4, off-load meter is drawn as follows:
Basic wheel is divided into 4 wheels during off-load, and movement threshold value is set as: 49.2Hz, 49.0Hz, 48.8Hz, 48.6Hz, being pressed
According to actual electric network off-load mode, the original off-load amount of every round is successively set as the 30% of power shortage, 25%, 25%, 20%.Separately
Outside, the special wheel of a wheel is also set up, movement return frequency is set as 49.6Hz, time delay of action 15s.The criterion whether off-load terminates
Are as follows:
f′k<0&f≤fi
Wherein, f is system frequency at this time, fiI-th for setting takes turns frequency node, f 'kFor the frequency at kth wheel off-load moment
Change rate.In the step 5, excision sequence is determined according to the importance of load, utilizes " electric system automatic low-frequency load shedding skill
Art provides " can obtain load the comprehensive weight of significance level expression formula are as follows:
λ in formulai, βijThe whole network importance value of respectively i-th load bus and the jth class of i-th of load bus
The local significance level of load.
In the step 6, during executing off-load scheme, virtual inertia response/primary tune of Wind turbines is utilized
Frequency control carries out real-time update to the calculated value of power shortage, and power shortage corrected value is as follows before next round off-load is implemented:
In above formula, HΣiIt (t) is the equivalent inertia time constant of every wheel for following time-varying, Δ PmWF.iFor wind before the i-th wheel off-load
Machine primary frequency modulation controls the active increment provided system.Pd.oldTo correct preceding power shortage calculated, fa, fbRespectively
I takes turns the system frequency before and after off-load.
A kind of meter of the present invention and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response, advantage are:
Because fully considering wind-powered electricity generation virtual inertia and primary frequency modulation effect, the electric system to the frequency modulation containing wind-powered electricity generation, mentioned method can be more acurrate
Each coherence low-frequency load shedding of computing system power vacancy, more precise manipulation, frequency retrieval dynamic process and stable state more meet true feelings
Condition.
Detailed description of the invention
Present invention will be further explained below with reference to the attached drawings and examples:
Fig. 1 is flow chart of the invention.
Fig. 2 is the analogue system of the embodiment of the present invention.
Fig. 3 is the frequency variation diagram for implementing different off-load schemes after meter and wind-powered electricity generation virtual inertia response/primary frequency modulation control.
Fig. 4 is virtual inertia function and effect figure in comparison off-load strategy.
Specific embodiment
A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds, Wind turbines apply allusion quotation
Under the conditions of type virtual inertia control strategy, the Analytical Solution virtual inertia of wind-powered electricity generation and the equivalent inertia time constant of electric system
(H∑), according to required real-time H∑, in conjunction with initial system frequency change rate, calculate system accurately total active power shortage.Together
When Wind turbines apply typical primary frequency modulation control strategy, solved wind power plant primary frequency modulation response transmission function and primary
The active increment of frequency modulation response calculates the power operated during low-frequency load shedding by wheel in conjunction with the virtual inertia dynamic of wind-powered electricity generation of solution
Vacancy, and optimize off-load amount accordingly.
Specific step is as follows:
Step 1: to containing be applied with virtual inertia control Wind turbines the equivalent inertia time constant of electric system into
Row analytical Calculation；
Step 2: according to the equivalent inertia time constant of electric system of the virtual inertia containing wind-powered electricity generation acquired, calculating electric system
Active power shortage；
Step 3: being asked according to the primary frequency modulation transfer function model of separate unit Wind turbines using weighted equivalent polymerization
Obtain the active increment of wind power plant primary frequency modulation；
Step 4: doing off-load meter according to the initial active power shortage of electric system and draw；
Step 5: according to the significance level of load, determining off-load place；
Step 6: according to the equivalent inertia time constant of accurate electric system and the active increment of primary frequency modulation, by wheel
Secondary ground dynamic optimization off-load amount.
In the step 1, according to the definition of inertia time constant, the H of separate unit blowerequAre as follows:
In formula, ωnom,ωs0,ωr0,HDFIG,Tf,KdfRespectively rated angular velocity, system initial synchronisation angular speed, initially
Rotor velocity, double-fed blower built-in inertia time constant, time constant filter, inertia control gain.KpT,KiTFor speed control
Device parameter processed.S is Laplace operator.By the formula, the inertia time constant of wind power plant separate unit blower can be obtained.
It is a unit by wind power plant equivalence when the wind power plant DFIG blower of platform containing p, then wind power plant polymerization inertia time is normal
Number HeqWFiAre as follows:
Wherein,For the mean speed of i-th of wind power plant inner blower, HDFIGIt is inherently used for double-fed blower
Property time constant, if CHi=ωs0HDFIGi/ω2 nomi,CHiFor a constant.By the formula, can obtain entire wind power plant etc.
Imitate inertia time constant.
Power plant in system is divided into two class of wind power plant of conventional power unit, the control containing virtual inertia, and station quantity is distinguished
For a, b, i.e. conventional power plant has a, and the wind power plant of the control containing virtual inertia has b；
The then equivalent inertia time constant H of the step 1∑Expression formula are as follows:
In formula: H0,HeqWFi,ΔH,HCONiThe respectively inertia time constant of conventional power plant controls wind containing virtual inertia
The equivalent inertia time constant of electric field considers that system equivalent inertia time constant increment when the response of wind-powered electricity generation virtual inertia, separate unit are normal
Advise the inertia time constant of generator, SCONi,SeqWFiThe respectively rated capacity of conventional power unit controls wind-powered electricity generation containing virtual inertia
?.
By step 1, the electric system containing wind field can be obtained with the equivalent inertia time constant of time-varying, being difficult to parse
Equivalent inertia time constant parsed come out, lay the foundation to acquire accurate power shortage.
In the step 2, the expression formula of electric system active power shortage are as follows:
Wherein, H∑For the equivalent inertia time constant of system, system frequency change rate when df/dt is disturbance.
By step 2, the accurate power shortage of system is acquired, it is next compared to the case where inertia time constant is constant is set
It says, this method is due to being that the equivalent inertia time constant of actual measurement can effectively prevent cutting so that the data of power shortage are more accurate
Or owe the case where cutting.
In the step 3, primary frequency modulation transfer function model assists control strategy using the primary frequency modulation of revolving speed control
When, solve the transmission function h that wind power plant primary frequency modulation responds equivalent polymerization modelmWF(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) are as follows:
Wherein b0, a0, a1, a2, a3For hwt(s) transfer-function coefficient, KpfFor sagging control gain.
After polymerization, dynamic response model transmission function h of the wind power plant based on revolving speed primary frequency modulation control systemmWF(s) are as follows:
Wherein b0G, a0G, a1G, a2G, a3GRespectively transmission function hmWF(s) every equivalent parameters.
It is hereby achieved that the active support that primary frequency modulation provides:
Through reverse drawing Laplace transform, time-domain expression is obtained:
L-1[△PmWF(s)]=△ PmWF(t)
Wherein, Δ PmWF(s),ΔPmWF(t) be the active increment of wind-powered electricity generation primary frequency modulation frequency domain and frequency domain presentation form.
It can be the active increment that system provides in the hope of primary frequency modulation when system disturbs by step 3,
There is the support of the increment, it will so that off-load amount is less.
In the step 4, off-load meter is drawn as follows:
Basic wheel is divided into 4 wheels during off-load, and movement threshold value is set as: 49.2Hz, 49.0Hz, 48.8Hz, 48.6Hz, being pressed
According to actual electric network off-load mode, the original off-load amount of every round is successively set as the 30% of power shortage, 25%, 25%, 20%.Separately
Outside, the special wheel of a wheel is also set up, movement return frequency is set as 49.6Hz, time delay of action 15s.The criterion whether off-load terminates
Are as follows:
f′k<0&f≤fi
Wherein, f is system frequency at this time, fiI-th for setting takes turns frequency node, f 'kFor the frequency at kth wheel off-load moment
Change rate.
By step 4, the set cutting load amount of frequency movement threshold value and beginning provided with system makes off-load work
Planned progress.Off-load ending-criterion is set, so that system will not occur to cut situation.
In the step 5, excision sequence is determined according to the importance of load, utilizes " electric system automatic low-frequency load shedding
Technical stipulation " can obtain load the comprehensive weight of significance level expression formula are as follows:
λ in formulai, βijThe whole network importance value of respectively i-th load bus and the jth class of i-th of load bus
The local significance level of load.
Step 5 provides excision sequence for off-load process, first cuts off unessential load.As China's electric load is divided into one
Grade, second level, three stage loads.During Load Shedding According To Frequency, 3 stage loads should be preferentially cut off.
In the step 6, during executing off-load scheme, virtual inertia response/primary tune of Wind turbines is utilized
Frequency control carries out real-time update to the calculated value of power shortage, and power shortage corrected value is as follows before next round off-load is implemented:
In above formula, HΣiIt (t) is the equivalent inertia time constant of every wheel for following time-varying, Δ PmWF.iFor wind before the i-th wheel off-load
Machine primary frequency modulation controls the active increment provided system.Pd.oldTo correct preceding power shortage calculated, fa, fbRespectively
I takes turns the system frequency before and after off-load.
Step 6 all optimizes original off-load amount using virtual inertia/primary frequency modulation before every wheel off-load, so that
While meeting system frequency stability range, total off-load amount is less, to achieve the purpose that less cutting load.
Implement example:
Under Matlab/simulink environment, the analogue system of Fig. 4 is established, two regions pass through two connection in system
Winding thread connection, region 1 include a Hydropower Unit G2 and a wind power plant, and region 2 includes two fired power generating units G3 and G4, load
L1, L2, C1, C2 respectively at two Area Interfaces buses access, load L3 as disturbance load, by L3 access and cut off come
Simulate the frequency accident of the analogue system power shortage.
Wherein simulation parameter is as follows: 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, excitation
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.
Generator parameter (G2, G3, G4): Sn=900MVA, Un=20kV, Xd=1.8, Xq=1.7, Xa=0.2, Xd'=
0.3, Xq'=0.55, Xd"=0.25, Xq"=0.25, Ra=0.0025, Td0'=8.0, Tq0'=0.4, Td0"=0.03, Tq0″
=0.05, H=6.5 (G2), H=6.175 (G3, G4)
Transformer parameter (T1, T2, T3, T4): Sn=900MVA, Un1/Un2=20KV/230KV, Rt+jXt=0+
J0.15pu transmission line parameter (on the basis of 100MVA, 230kV):
RL=0.0001pu/km, XL=0.001pu/km, BC=0.00175pu/km
Load data: PL1=800MW, QL=100MVAR, QC1=-187MVAR, QC2=-200MVAR, PL2=800MW,
QL=100MVAR, QC1=-187MVAR, QC2=-350MVAR additional load PL3=900MW
In verification process, emulation project includes: (1) meter and wind-powered electricity generation virtual inertia response/primary frequency modulation control situation, with this
Scheme of the invention, Traditional Method, half adjustment procedure carry out excision load to the electric system under disturbance.The present invention program is in each of off-load
Wheel carries out real-time update optimization to the value of power shortage before system next one off-load is implemented according to step 6, obtains scheme one
Shown in curve.Traditional Method calculates to obtain system initial power vacancy by step 2, then according to step 4,5 excision loads, the side of obtaining
Curve shown in case two.The off-load amount optimization of half adjustment procedure first round more than Traditional Method, obtains curve shown in scheme three.
(2) meter and wind-powered electricity generation virtual inertia response/primary frequency modulation control action after off-load curve, only consider primary frequency modulation control action after
Off-load curve.
Fig. 3 is shown when considering wind-powered electricity generation virtual inertia response/primary frequency modulation control, with the present invention program, tradition
Scheme, half adjustment procedure cut off the change curve of system frequency after system loading, and the simulation process of corresponding diagram 3 leads to every numerical value
Table 1 is crossed to be compared.
The off-load amount and steady frequency of each off-load scheme of table 1
According to simulation result it can be found that: 1) from 1 → scheme of scheme 2, the scheme 1 in the case where system water averagely restores stable situation
That is the off-load amount of the present invention program is minimum, is 688.6MW, and the off-load amount of scheme 2 is 717.6MW, cuts off than scheme more than 1
29MW, an and important indicator for measuring low-frequency load shedding scheme superiority and inferiority is exactly off-load amount is guaranteeing that system frequency restores stable
In the case of off-load amount should lack as far as possible.Thus illustrate, meter and wind-powered electricity generation virtual inertia/primary frequency modulation response under high wind-powered electricity generation permeability
Adaptive low frequency deloading method off-load amount is less, and off-load amount can be effectively reduced by demonstrating the present invention program.2) from 2 → scheme of scheme
3, scheme 2 cuts off load using Traditional Method, and Traditional Method is that pre-set amount is cut off at frequency node, does not account for blower
Inertia support and the active increment of primary frequency modulation, therefore off-load amount is maximum, is 717.6MW, scheme 3 is half adjustment procedure, in off-load the
One wheel determines off-load amount according to the amount of frequency change rate, improves only the first run, can be seen that from figure, basic after the second wheel
It is consistent with traditional scheme.The off-load amount of Comprehensive Correlation, the present invention program is optimal.
Fig. 4, which is shown, demonstrates influence of the virtual inertia to off-load process, the emulation of corresponding diagram 4 using the present invention program
Every numerical value is compared by journey by table 2.
Whether table 2 considers off-load amount and steady frequency after virtual inertia
According to simulation result it can be found that, it is contemplated that blower virtual inertia effect can provide a system to effective inertia branch
Support, slows down frequency decrease speed, improves frequency descending depth.It is found in the comparison of off-load amount, does not consider virtual inertia
When system cut off 700.4MW, than having cut off 11.8MW after considering virtual inertia more.Thus, considering the virtual used of blower
Property recovery of the response to frequency, it is all helpful to system off-load amount.
Claims (8)
1. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response, it is characterised in that: wind turbine
Group applies under typical virtual inertia control policy condition, the Analytical Solution virtual inertia of wind-powered electricity generation and the equivalent inertia time of electric system
Constant (H∑), according to required real-time H∑, in conjunction with initial system frequency change rate, calculating system, accurately total active power is lacked
Volume；Simultaneously when Wind turbines apply typical primary frequency modulation control strategy, wind power plant primary frequency modulation transmission function and one have been solved
The active increment of secondary frequency modulation response calculates low-frequency load shedding in the process by wheel operation in conjunction with the virtual inertia dynamic of wind-powered electricity generation of solution
Power shortage, and optimize off-load amount accordingly.
2. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response, it is characterised in that including with
Lower step:
Step 1: the equivalent inertia time constant of electric system containing the Wind turbines for being applied with virtual inertia control is solved
Analysis calculates；
Step 2: according to the equivalent inertia time constant of electric system of the virtual inertia containing wind-powered electricity generation acquired, it is active to calculate electric system
Power shortage；
Step 3: wind is acquired using weighted equivalent polymerization according to the primary frequency modulation transfer function model of separate unit Wind turbines
The active increment of electric field primary frequency modulation；
Step 4: doing off-load meter according to the initial active power shortage of electric system and draw；
Step 5: according to the significance level of load, determining off-load place；
Step 6: according to the equivalent inertia time constant of accurate electric system and the active increment of primary frequency modulation, by round
Dynamic optimization off-load amount.
3. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response according to claim 2,
It is characterized by: in the step 1, according to the definition of inertia time constant, the H of separate unit blowerequAre as follows:
In formula, ωnom, ωs0, ωr0, HDFIG,Tf,KdfRespectively rated angular velocity, system initial synchronisation angular speed, initial rotor
Angular speed, double-fed blower built-in inertia time constant, time constant filter, inertia control gain；KpT,KiTFor speed control
Parameter；
It is a unit by wind power plant equivalence, then wind power plant polymerize inertia time constant when the wind power plant DFIG blower of platform containing p
HeqWFiAre as follows:
Wherein,For the mean speed of i-th of wind power plant inner blower, CHi=ωs0HDFIGi/ω2 nomi；
Power plant in system is divided into two class of wind power plant of conventional power unit, the control containing virtual inertia, and station quantity is respectively a,
b；
The then equivalent inertia time constant H of the step 1∑Expression formula are as follows:
In formula: H0, HeqWFi, Δ H is respectively the inertia time constant of conventional power plant, controls the equivalent of wind power plant containing virtual inertia
System equivalent inertia time constant increment when inertia time constant and consideration wind-powered electricity generation virtual inertia response, SCONi, SeqWFiRespectively
For the rated capacity of conventional power unit, wind power plant is controlled containing virtual inertia.
4. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response according to claim 2,
It is characterized by: in the step 2, the expression formula of electric system active power shortage are as follows:
Wherein, H∑For the equivalent inertia time constant of system, system frequency change rate when df/dt is disturbance.
5. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response according to claim 2,
It is characterized by: in the step 3, primary frequency modulation transfer function 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) are as follows:
Wherein b0, a0, a1, a2, a3For hwt(s) transfer-function coefficient, KpfFor sagging control gain；
After polymerization, dynamic response model transmission function h of the wind power plant based on revolving speed primary frequency modulation control systemmWF(s) are as follows:
Wherein b0G, a0G, a1G, a2G, a3GRespectively transmission function hmWF(s) every equivalent parameters；
It is hereby achieved that the active support that primary frequency modulation provides:
Through reverse drawing Laplace transform, time-domain expression is obtained:
L-1[△PmWF(s)]=△ PmWF(t)
Wherein, Δ PmWFIt (s) is the active increment of wind-powered electricity generation primary frequency modulation.
6. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response according to claim 2,
It is characterized by: off-load meter is drawn as follows in the step 4:
Basic wheel is divided into 4 wheels during off-load, and movement threshold value is set as: 49.2Hz, 49.0Hz, 48.8Hz, 48.6Hz, according to reality
The original off-load amount of every round is successively set as the 30% of power shortage, 25%, 25%, 20% by border power grid off-load mode；In addition,
The special wheel of a wheel is also set up, movement return frequency is set as 49.6Hz, time delay of action 15s；The criterion whether off-load terminates are as follows:
f′k<0&f≤fi
Wherein, f is system frequency at this time, fiI-th for setting takes turns frequency node, f 'kChange for the frequency at kth wheel off-load moment
Rate.
7. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response according to claim 2,
It is characterized by: determining excision sequence in the step 5 according to the importance of load, utilizing " electric system automatic low-frequency Reduction of Students' Study Load
Lotus technical stipulation " can obtain load the comprehensive weight of significance level expression formula are as follows:
λ in formulai, βijThe whole network importance value of respectively i-th load bus and the jth type load of i-th of load bus
Local significance level.
8. a kind of meter and the adaptive low frequency deloading method of wind-powered electricity generation virtual inertia/primary frequency modulation response according to claim 2,
It is characterized by: in the step 6, during executing off-load scheme, using Wind turbines virtual inertia response/it is primary
Frequency modulation control carries out real-time update to the calculated value of power shortage, and power shortage corrected value is as follows before next round off-load is implemented:
In above formula, HΣiIt (t) is the equivalent inertia time constant of every wheel for following time-varying, Δ PmWF.iFor blower one before the i-th wheel off-load
The active increment that secondary frequency modulation control provides system；Pd.oldTo correct preceding power shortage calculated, fa, fbRespectively i-th wheel
System frequency before and after off-load.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810719653.3A CN109038613A (en) | 2018-07-02 | 2018-07-02 | A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810719653.3A CN109038613A (en) | 2018-07-02 | 2018-07-02 | A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109038613A true CN109038613A (en) | 2018-12-18 |
Family
ID=65522230
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810719653.3A Pending CN109038613A (en) | 2018-07-02 | 2018-07-02 | A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109038613A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109728590A (en) * | 2018-12-28 | 2019-05-07 | 南京工业大学 | The self-adaptation control method of Wind turbines participation primary frequency modulation |
| CN109787244A (en) * | 2019-03-18 | 2019-05-21 | 重庆大学 | A kind of method that determining power distribution network low pressure accelerates off-load total amount |
| CN110518632A (en) * | 2019-08-05 | 2019-11-29 | 三峡大学 | A kind of wind farm grid-connected quantitative calculation method that power grid inertia is weakened |
| CN110594093A (en) * | 2019-08-29 | 2019-12-20 | 湖北工业大学 | Double-fed fan inertia control method based on second-order frequency differential of power system |
| CN110649596A (en) * | 2019-09-05 | 2020-01-03 | 大连理工大学 | Frequency full-response analytic model considering system initial state |
| CN111900742A (en) * | 2020-07-02 | 2020-11-06 | 武汉大学 | Frequency modulation method of wind storage system based on double-layer cooperative control |
| CN113007039A (en) * | 2021-03-18 | 2021-06-22 | 南方电网科学研究院有限责任公司 | Wind power plant inertia response testing method and device |
| CN113746134A (en) * | 2021-08-20 | 2021-12-03 | 华北电力大学（保定） | Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102394497A (en) * | 2011-10-18 | 2012-03-28 | 重庆大学 | Microgrid coordination control system and method based on state transition |
| CN105162164A (en) * | 2015-09-22 | 2015-12-16 | 河海大学 | Method of building low-order dynamic frequency response model with wind power integration system |
| CN106786620A (en) * | 2016-12-29 | 2017-05-31 | 北京四方继保自动化股份有限公司 | Emergency load control method under island network malfunction |
| CN106849088A (en) * | 2017-02-17 | 2017-06-13 | 三峡大学 | It is a kind of based on pitch control wind-powered electricity generation it is active/frequency coupling electrical power system response computational methods |
| CN106910142A (en) * | 2017-02-17 | 2017-06-30 | 三峡大学 | A kind of power system frequency characteristic computing method containing the active frequency coupling of wind-powered electricity generation |
| CN107482652A (en) * | 2017-08-22 | 2017-12-15 | 山东大学 | A kind of power system UFLS Implementation of Virtual Experiment and system |
-
2018
- 2018-07-02 CN CN201810719653.3A patent/CN109038613A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102394497A (en) * | 2011-10-18 | 2012-03-28 | 重庆大学 | Microgrid coordination control system and method based on state transition |
| CN105162164A (en) * | 2015-09-22 | 2015-12-16 | 河海大学 | Method of building low-order dynamic frequency response model with wind power integration system |
| CN106786620A (en) * | 2016-12-29 | 2017-05-31 | 北京四方继保自动化股份有限公司 | Emergency load control method under island network malfunction |
| CN106849088A (en) * | 2017-02-17 | 2017-06-13 | 三峡大学 | It is a kind of based on pitch control wind-powered electricity generation it is active/frequency coupling electrical power system response computational methods |
| CN106910142A (en) * | 2017-02-17 | 2017-06-30 | 三峡大学 | A kind of power system frequency characteristic computing method containing the active frequency coupling of wind-powered electricity generation |
| CN107482652A (en) * | 2017-08-22 | 2017-12-15 | 山东大学 | A kind of power system UFLS Implementation of Virtual Experiment and system |
Non-Patent Citations (3)
| Title |
|---|
| 李世春 等: "基于转速控制的双馈风电机组一次调频辅助控制系统建模", 《中国电机工程学报》 * |
| 李世春 等: "风电场等效虚拟惯性时间常数计算", 《电力系统自动化》 * |
| 李顺 等: "高风电渗透率下的自适应低频减载策略研究", 《电网技术》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109728590B (en) * | 2018-12-28 | 2022-08-19 | 南京工业大学 | Self-adaptive control method for wind turbine generator to participate in primary frequency modulation |
| CN109728590A (en) * | 2018-12-28 | 2019-05-07 | 南京工业大学 | The self-adaptation control method of Wind turbines participation primary frequency modulation |
| CN109787244A (en) * | 2019-03-18 | 2019-05-21 | 重庆大学 | A kind of method that determining power distribution network low pressure accelerates off-load total amount |
| CN109787244B (en) * | 2019-03-18 | 2021-07-30 | 重庆大学 | Method for determining total low-voltage accelerated load shedding amount of power distribution network |
| 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 |
| CN110594093A (en) * | 2019-08-29 | 2019-12-20 | 湖北工业大学 | Double-fed fan inertia control method based on second-order frequency differential of power system |
| CN110649596A (en) * | 2019-09-05 | 2020-01-03 | 大连理工大学 | Frequency full-response analytic model considering system initial state |
| CN110649596B (en) * | 2019-09-05 | 2022-09-06 | 大连理工大学 | Frequency full-response analytic model considering system initial state |
| CN111900742A (en) * | 2020-07-02 | 2020-11-06 | 武汉大学 | Frequency modulation method of wind storage system based on double-layer cooperative control |
| CN113007039A (en) * | 2021-03-18 | 2021-06-22 | 南方电网科学研究院有限责任公司 | Wind power plant inertia response testing method and device |
| CN113007039B (en) * | 2021-03-18 | 2022-06-14 | 南方电网科学研究院有限责任公司 | Wind power plant inertia response testing method and device |
| CN113746134A (en) * | 2021-08-20 | 2021-12-03 | 华北电力大学（保定） | Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109038613A (en) | A kind of adaptive low frequency deloading method counted and wind-powered electricity generation virtual inertia/primary frequency modulation responds | |
| CN109066770B (en) | Control method and device for accessing wind power to flexible direct current power transmission system | |
| Sarasúa et al. | Dynamic response and governor tuning of a long penstock pumped-storage hydropower plant equipped with a pump-turbine and a doubly fed induction generator | |
| US9822765B2 (en) | Wind power plant control system | |
| CN105162164B (en) | A kind of method for the low order dynamic frequency response model for establishing the system containing wind-electricity integration | |
| CN108832658A (en) | A kind of wind power penetration limit calculation method considering frequency constraint and wind-powered electricity generation frequency modulation | |
| CN108365633A (en) | A kind of doubly-fed wind turbine hypervelocity virtual inertia control method of Control of decreasing load variable element | |
| CN108695857B (en) | Automatic voltage control method, device and system for wind power plant | |
| CN104333037A (en) | Cooperative control method for participating in frequency modulation and pressure regulation of power system by wind storage cluster | |
| CN106910142A (en) | A kind of power system frequency characteristic computing method containing the active frequency coupling of wind-powered electricity generation | |
| CN110829487A (en) | Dynamic frequency prediction method for power system | |
| CN108631359A (en) | Wind power plant real-time simulation modeling and method | |
| CN102594244A (en) | Joint control method of primary frequency modulation for doubly-fed wind power generation set | |
| CN107895955A (en) | A kind of cooperative control method of wind-powered electricity generation compensation hydraulic turbine water hammer effect | |
| CN102340140A (en) | Method for automatically fast calculating stability limit of large-scale interconnected power grid and stability control strategy | |
| CN108711868A (en) | It is a kind of meter and islet operation voltage security GA for reactive power optimization planing method | |
| CN104485670B (en) | The control method of voltage sensitivity industrial load time-varying damping characteristic in island network | |
| CN105186511B (en) | Battery energy storage system participates in electric grid secondary frequency modulation control method | |
| Spichartz et al. | New stability concept for primary controlled variable speed wind turbines considering wind fluctuations and power smoothing | |
| CN107346889B (en) | Load reduction optimization model construction method considering primary and secondary frequency modulation and minimum frequency deviation | |
| CN102904266A (en) | Method for determining inactive compensation capacity network adaptability of wind power plant | |
| CN107979112B (en) | Fan control method, system, terminal and readable storage medium | |
| CN111130122B (en) | Online monitoring method for reactive power control capability of wind power plant | |
| CN117117901A (en) | Frequency control method of offshore wind power flexible-direct system | |
| CN115800296B (en) | Voltage frequency collaborative supporting method for open sea wind power through VSC-MTDC grid-connected system |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181218 |
|
| RJ01 | Rejection of invention patent application after publication |