CN102136036B - Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system - Google Patents

Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system Download PDF

Info

Publication number
CN102136036B
CN102136036B CN 201110070510 CN201110070510A CN102136036B CN 102136036 B CN102136036 B CN 102136036B CN 201110070510 CN201110070510 CN 201110070510 CN 201110070510 A CN201110070510 A CN 201110070510A CN 102136036 B CN102136036 B CN 102136036B
Authority
CN
China
Prior art keywords
double
unit
wind energy
energy turbine
turbine set
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.)
Expired - Fee Related
Application number
CN 201110070510
Other languages
Chinese (zh)
Other versions
CN102136036A (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN 201110070510 priority Critical patent/CN102136036B/en
Publication of CN102136036A publication Critical patent/CN102136036A/en
Application granted granted Critical
Publication of CN102136036B publication Critical patent/CN102136036B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a double-feed wind power station equivalent modeling method applied to analysis on the small signal stability of a power system. Point of interconnection frequency serves as a boundary signal, and the wind power station is equivalent to a power source controlled by the frequency of a power grid, so that a high-dimensional dynamic model of the double-feed wind power station is simplified into an algebraic constraint equation. Compared with the prior art, the invention simplifies the high-dimensional dynamic model of the wind power station into the simple algebraic constraint equation so as to greatly simplify the mathematical modeling of the wind power station. Simulation researches made under a plurality of conditions prove that the equivalent simplified model is effective for the analysis on characteristic values of the power system with the wind power station in an electromechanical mode.

Description

Be suitable for the double-fed fan motor station equivalent modeling method of Study of Power System Small Disturbance stability analysis
Technical field
The present invention relates to electrical engineering technology system, particularly relate to the modeling method of grid type double-feedback constant frequency wind generator group wind power station in the electric system.
Background technology
In recent years, China's Construction of Wind Power keeps high speed development situation always, and along with the rapid growth of wind-powered electricity generation installation scale, wind energy turbine set highlights day by day on the impact of electric system, wind energy turbine set modeling and the Electrical Power System Dynamic properties influence is become the study hotspot of this area.Because it is tens of to over one hundred that the unit number of Large Scale Wind Farm Integration can reach, when power system stability was analyzed, if all units all use detailed model, its computing cost was difficult to bear, so it is very necessary to set up the wind energy turbine set Equivalent Model.Document [1] has proposed to contain the continuous tide computation model of wind-powered electricity generation unit, can be used for static voltage stability research; Document [2] utilizes singular perturbation theory to set up the wind energy turbine set reduced-order model; Document [3-5] replaces wind energy turbine set with the combination of a typhoon turbine and a double feedback electric engine, is used for research wind energy turbine set dynamic behaviour and on the impact of electrical network damping characteristic; After document [6] hived off the wind energy turbine set unit, equivalence was multi-computer system, and simulation precision improves; The Equivalent Model that document [7] provides has been considered the wake effect of wind energy turbine set, has taken into account the blower fan position to injecting the impact of wind speed.
In stability of power system research, Equivalent Model in the past still has the following disadvantages: 1. put undue emphasis on wind speed and change the wind energy turbine set power response that causes, when Transient Stability Analysis, the time scale that head pendulum destabilization problems is concerned about is no more than 2 seconds, and the variation of discussion wind speed waveform and wind speed has little significance within the so short period.In addition, for the microvariations stable problem, no matter be stability of balance point or structural stability, all do not contact directly with the wind speed waveform.2. the curves that covets degree.When verification Equivalent Model validity, the curves degree varies is the evaluation index of science surely.For transient stability, should be take power limit or critical clearing time as evaluation index; For microvariations stability, frequency and the damping of testing dominant pattern in the presence of re-graduation.Therefore, on wind energy turbine set equivalent modeling problem, should according to problem character equivalent strategy be proposed towards the needs of particular problem, and the validity of scientific evaluation Equivalent Model.
[1] Wang Chengshan, Sun Wei, Wang Xinggang (WANG Chengshan, SUN Wei, WANG Xinggang).The Power System Total Transfer Capacity that contains Large Scale Wind Farm Integration calculates (Total transfer capability of power system including large-scale wind farm) [J]. Automation of Electric Systems (Automation of Electric Power Systems), 2007,31 (2): 17-21,31;
[2] M.G. Castro, J.M. Ferreira. A wind park reduced-order model using singular perturbations theory[J]. IEEE transactions on Energy Conversion, 1996, 11(4):735-741；
[3] L. M. Fernandez, F. Juradob, J. R. Saenz. Aggregated dynamic model for wind farms with doubly fed induction generator wind turbines [J]. Renewable Energy, 2008, 33(1):129-140；
[4] G. G. Miguel, M. P. Comech, J. Sallan, A. Llombart. Modelling wind farms for grid disturbance studies [J]. Renewable Energy, 2008, 33(1):2109-2121；
[5] Wang Chen, Shi Libao, Yao Liangzhong, Deng. (WANG Chen, SHI Libao, YAO Liangzhong, et al), the microvariations stability analysis (Small signal stability analysis of the large-scale wind farm with DFIGs) [J] of extensive Double-feed wind power field. Proceedings of the CSEE (Proceedings of the CSEE), 2010,30 (4): 63-70;
[6] Akhmatov V, Knudsen H. An aggregated model of a grid-connected, large-scale offshore wind farm for power stability investigations—importance of windmill mechanical system. Electric Power Energy Syst. 2002, 24: 709–717；
[7] L. M. Fernandez, C. A. Garcia, F. Jurado, J. R. Saenz. Aggregation of Doubly Fed Induction Generators Wind Turbines Under Different Incoming Wind Speeds[C]// IEEE Power Tech conference, St. Petersburg,Russia:2005。
Summary of the invention
Based on above-mentioned prior art, the present invention proposes a kind of double-fed fan motor station equivalent modeling method that is suitable for the Study of Power System Small Disturbance stability analysis, proposed to be incorporated into the power networks dot frequency as sideband signal, be a power source that is subjected to mains frequency control with the wind energy turbine set equivalence, realize simplifiedly the method for wind energy turbine set equivalent modeling.
The present invention proposes a kind of double-fed fan motor station equivalent modeling method that is suitable for the Study of Power System Small Disturbance stability analysis, the dot frequency that will be incorporated into the power networks is as sideband signal, be a power source that is subjected to mains frequency control with the wind energy turbine set equivalence, so that the double-fed fan motor field becomes the Algebraic Constraint equation by the higher-dimension Simplification of Dynamic Model, the method may further comprise the steps:
Step 1 is set up the microvariations power response model of double feedback electric engine under the sinusoidal excitation
If have a double-fed fan motor unit in the double-fed fan motor field, wherein and the mains frequency of site B be
Figure 656027DEST_PATH_IMAGE001
, voltage magnitude is
Figure 514261DEST_PATH_IMAGE002
If
Figure 798612DEST_PATH_IMAGE001
Little value sine-wave oscillation occurs, namely
Figure 742297DEST_PATH_IMAGE003
,
Figure 504717DEST_PATH_IMAGE004
Amplitude,
Figure 220388DEST_PATH_IMAGE005
It is oscillation frequency.Get
Figure 675640DEST_PATH_IMAGE006
,
Figure 106621DEST_PATH_IMAGE007
, the first squiggle of the meritorious output of unit;
Again order ,
Figure 239979DEST_PATH_IMAGE006
,
Figure 866133DEST_PATH_IMAGE007
, the second squiggle of the meritorious output of unit;
Wherein, the oscillation amplitude of the first squiggle is 60kW, and the oscillation amplitude of the second squiggle is 2 kW;
Step 2; Transfer function matrix with single unit
Figure 253252DEST_PATH_IMAGE009
Simplify processing, further comprising the steps of:
The first step: according to the analysis conclusion of the power response relevant treatment of double feedback electric engine under the sinusoidal excitation, ignore input signal
Figure 888633DEST_PATH_IMAGE010
Effect, namely at the transfer function matrix of wind energy turbine set dynamic subsystem
Figure 44807DEST_PATH_IMAGE009
In, ignore
Figure 373021DEST_PATH_IMAGE011
With
Figure 981856DEST_PATH_IMAGE012
, will
Figure 686507DEST_PATH_IMAGE009
Be reduced to
Figure 697189DEST_PATH_IMAGE013
:
Figure 196303DEST_PATH_IMAGE014
Wherein,
Figure 823593DEST_PATH_IMAGE015
Represent the input of the bus frequency that is incorporated into the power networks, take the meritorious transport function as output of blower fan;
Figure 269618DEST_PATH_IMAGE016
Expression is with the input of the bus frequency that is incorporated into the power networks, take the idle transport function as output of blower fan;
Figure 931544DEST_PATH_IMAGE011
Expression is take the busbar voltage that is incorporated into the power networks as input, take the meritorious transport function as output of blower fan;
Figure 70401DEST_PATH_IMAGE012
Expression is take the busbar voltage that is incorporated into the power networks as input, take the idle transport function as output of blower fan;
Transfer function matrix for conventional electrical network dynamic subsystem
Figure 184988DEST_PATH_IMAGE017
,
Figure 434703DEST_PATH_IMAGE010
Its output signal, correspondingly,
Figure 951135DEST_PATH_IMAGE017
Be reduced to
Figure 260894DEST_PATH_IMAGE018
, have:
Figure 862777DEST_PATH_IMAGE019
；
Second step: consider Steady-state characteristic under sinusoidal excitation, under the Stable State of Sine situation, the output quantity of wind-powered electricity generation unit
Figure 287122DEST_PATH_IMAGE021
With
Figure 767782DEST_PATH_IMAGE022
It all is sinusoidal phasor With
Figure 448479DEST_PATH_IMAGE024
Figure 142765DEST_PATH_IMAGE020
Effect be a complex gain matrix , wherein,
Figure 370801DEST_PATH_IMAGE026
,
Figure 297169DEST_PATH_IMAGE027
Electrical network subsystem after the simplification (
Figure 845962DEST_PATH_IMAGE018
) be two single-input single-output system (SISO system)s; The double-fed unit (
Figure 473951DEST_PATH_IMAGE020
) be single input two output systems;
Step 3, the simplification interface of double-fed fan motor field and electrical network is processed, and comprising:
If wind energy turbine set is installed n platform double-fed fan motor unit, the upstream of this n platform unit is common bus B, and when little value vibration appearred in conventional electrical network, for arbitrary unit, its input variable all was phasor
Figure 740984DEST_PATH_IMAGE028
In conjunction with the coordination of n platform unit, the power stage that wind energy turbine set is total With
Figure 874343DEST_PATH_IMAGE022
It is the phasor splicing of each double-fed unit power stage.The reciprocation expression formula of whole n machine wind energy turbine set and electrical network is:
Figure 133286DEST_PATH_IMAGE029
Based on the reciprocation of described wind energy turbine set and electrical network, realization contains the simplification analysis of the dynamo-electric pattern of electric system of double-fed fan motor field, and step is as follows:
The first step: the output of wind electric field that calculates the equilibrium point place by trend
Figure 153194DEST_PATH_IMAGE030
With , with power be With
Figure 374594DEST_PATH_IMAGE031
Constant power source access bus B;
Second step according to the mathematical model of separate unit double-fed fan motor unit, calculates after the linearization
Figure 147378DEST_PATH_IMAGE020
, and obtain the complex gain matrix
Figure 688081DEST_PATH_IMAGE032
The 3rd step is according to formula
Figure 597131DEST_PATH_IMAGE029
The dynamic behaviour equivalence of wind energy turbine set is served as reasons
Figure 197877DEST_PATH_IMAGE028
The controlled power source of control, controlled power source and constant power source are combined as Equivalent Model;
The 4th step embedded conventional electrical network with equivalent model, obtained the merit angle rocking curve of synchronous motor by the microvariations time-domain-simulation, and match obtains the leading electromechanical characteristics value of system through Prony again.Described mathematical model according to separate unit double-fed fan motor unit is calculated after the linearization
Figure 661219DEST_PATH_IMAGE020
, and obtain the complex gain matrix
Figure 536771DEST_PATH_IMAGE032
Step, have identical for the unit of model of the same race in the wind energy turbine set
Figure 34748DEST_PATH_IMAGE032
Battle array is only calculated a double-fed fan motor unit.
Compared with prior art, the invention enables wind energy turbine set to become a simple Algebraic Constraint equation by the higher-dimension Simplification of Dynamic Model, greatly simplified the Mathematical Modeling Methods to wind energy turbine set, and the simulation study under multiple situation confirms that Equivalent Simplification model in this paper is effective for the dynamo-electric mode characteristic values of electric system that analysis contains wind energy turbine set.
Description of drawings
Fig. 1 is grid connected wind power field synoptic diagram;
Fig. 2 is the reciprocation synoptic diagram of wind energy turbine set and conventional electrical network;
Fig. 3 is that grid type double-feedback wind-powered electricity generation unit structural drawing Fig. 4 is the meritorious response synoptic diagram under the sinusoidal signal excitation;
Fig. 5 is the simplification process synoptic diagram of single unit dynamic model;
Fig. 6 is the simplification schematic diagram of wind energy turbine set dynamic model;
Fig. 7 is the example system synoptic diagram of checking double-fed fan motor field Equivalent Model validity;
Fig. 8 is the merit angular response comparison diagram based on Equivalent Model and detailed model.
Embodiment
Under the microvariations meaning, the reciprocation of conventional electrical network and wind energy turbine set shows: the little change of sinusoidal pattern of the point voltage amplitude that is incorporated into the power networks is very little on the impact of wind energy turbine set, and sinusoidal little change of electric voltage frequency is larger on the wind energy turbine set impact.Based on this fact, the present invention proposes the dot frequency that to be incorporated into the power networks as sideband signal, be a power source that is subjected to mains frequency control with the wind energy turbine set equivalence, such processing has greatly been simplified the mathematical modeling to wind energy turbine set so that wind energy turbine set becomes a simple Algebraic Constraint equation by the higher-dimension Simplification of Dynamic Model.And the simulation study under multiple situation confirms that the Equivalent Simplification model that the present invention proposes is effective for the dynamo-electric mode characteristic values of electric system that analysis contains wind energy turbine set.
As shown in Figure 1, the double-fed fan motor field is at bus
Figure 71975DEST_PATH_IMAGE033
The place is incorporated into the power networks, and sideband signal is respectively
Figure 288192DEST_PATH_IMAGE010
,
Figure 701856DEST_PATH_IMAGE028
,
Figure 54340DEST_PATH_IMAGE021
With
Figure 262468DEST_PATH_IMAGE022
As shown in Figure 2, the reciprocation of expression wind energy turbine set and conventional electrical network, in this reciprocal process, for wind energy turbine set, the amplitude of the point voltage that is incorporated into the power networks and the variation of frequency will cause the meritorious and idle variation of wind energy turbine set, so
Figure 700402DEST_PATH_IMAGE010
With
Figure 917757DEST_PATH_IMAGE028
The input signal of wind energy turbine set,
Figure 124747DEST_PATH_IMAGE021
With
Figure 769355DEST_PATH_IMAGE022
It is its output signal; Otherwise just, the variation of injecting power will cause the variation of electrical network busbar voltage and frequency for electrical network, as seen for electrical network,
Figure 694586DEST_PATH_IMAGE021
With
Figure 918894DEST_PATH_IMAGE022
Become input signal,
Figure 777128DEST_PATH_IMAGE010
With
Figure 61479DEST_PATH_IMAGE028
Be output signal.
Two dynamic subsystem that feedback connects between electrical network and the wind energy turbine set.
Among Fig. 2
Figure 739585DEST_PATH_IMAGE017
Express conventional electrical network dynamic subsystem among Fig. 1;
Figure 767584DEST_PATH_IMAGE009
Express the wind energy turbine set dynamic subsystem.
Figure 480325DEST_PATH_IMAGE017
Output quantity be simultaneously
Figure 935577DEST_PATH_IMAGE009
Input,
Figure 366559DEST_PATH_IMAGE009
Output be simultaneously Input.Therefore, the reciprocation of conventional electrical network and wind energy turbine set is to connect by the feedback between subsystem to embody.According to the input, output signal dimension as can be known,
Figure 502846DEST_PATH_IMAGE017
With
Figure 129000DEST_PATH_IMAGE009
It all is the transfer function matrix of 2 * 2 dimensions.That is:
Figure 781698DEST_PATH_IMAGE034
（1）
Figure 151500DEST_PATH_IMAGE035
（2）
Each element is all expressed the dynamic perfromance of single-input single-output system in the transfer function matrix, as
Figure 838833DEST_PATH_IMAGE036
Expression is injected meritorious as input take the bus that is incorporated into the power networks, take the transport function of bus frequency as output, it has reflected that wind power injects effect and the result that electrical network is caused;
Figure 635888DEST_PATH_IMAGE015
Expression is take the bus frequency that is incorporated into the power networks as input, and take the meritorious transport function as output of wind-powered electricity generation, it has embodied the response characteristic that the bus frequency change that is incorporated into the power networks causes wind energy turbine set power.The physical significance of other transport functions can similar reasoning.
According to conventional electrical network and wind energy turbine set characteristics, transfer function matrix
Figure 775882DEST_PATH_IMAGE017
With
Figure 683795DEST_PATH_IMAGE009
In each element all be the high transport function of exponent number.In addition, above-mentioned reciprocation concerns, is equally applicable to the reciprocation of separate unit wind-powered electricity generation unit and electrical network, also is applicable in the interactive environment of separate unit wind-powered electricity generation unit and electrical network among the present invention.
The interactive simplification of double-fed fan motor field and electrical network is expressed, and comprises following processing:
1, the power response of double feedback electric engine under the sinusoidal excitation
When having weak (or negative) damping electromechanic oscillation mode in the synchronised grids, can cause low-frequency oscillation.In the low-frequency oscillation process, the voltage magnitude of associated bus and electric voltage frequency all show as the sine-wave oscillation of same frequency.As seen from Figure 1, the low-frequency oscillation behavior of conventional electrical network is equivalent to wind energy turbine set has been applied sinusoidal excitation signal; Wind energy turbine set then depends on the features of response of wind energy turbine set offset of sinusoidal excitation to the reaction result of conventional electrical network.Therefore, investigate the power response of double feedback electric engine under the sinusoidal excitation, and this " response equivalence " theoretical foundation as the wind energy turbine set equivalence.In addition, the excitation of wind energy turbine set has voltage and two signals of frequency, and the two response is had why not with also being analyzed.
If only have a double-fed fan motor unit in the wind energy turbine set, as shown in Figure 3.And the mains frequency of site B is
Figure 225635DEST_PATH_IMAGE001
, voltage magnitude is
Figure 193591DEST_PATH_IMAGE002
During the synchronised grids steady-state operation,
Figure 820881DEST_PATH_IMAGE001
With
Figure 532485DEST_PATH_IMAGE002
It all is constant.Now do following sinusoidal excitation microvariations l-G simulation test:
If
Figure 194411DEST_PATH_IMAGE001
Little value sine-wave oscillation occurs, namely
Figure 333268DEST_PATH_IMAGE003
,
Figure 182275DEST_PATH_IMAGE004
Amplitude,
Figure 697570DEST_PATH_IMAGE005
It is oscillation frequency.Get
Figure 214002DEST_PATH_IMAGE006
,
Figure 523761DEST_PATH_IMAGE007
, the waveform of the meritorious output of unit is shown in curve among Fig. 41.Again order
Figure 860064DEST_PATH_IMAGE008
, With Still get initial value, the waveform of the meritorious output of unit is shown in curve among Fig. 42.Curve data shows that the oscillation amplitude of curve 1 is 60kW, and the oscillation amplitude of curve 2 only is 2 kW.
Hence one can see that, and under sinusoidal signal excitation, double feedback electric engine is far longer than it to the response of voltage to the response of frequency, namely has the following fact: in the frequency range of being concerned about, as
Figure 765069DEST_PATH_IMAGE037
, have
Figure 854248DEST_PATH_IMAGE038
Be far longer than
Figure 445766DEST_PATH_IMAGE039
, equally also can verify:
Figure 405632DEST_PATH_IMAGE040
Be far longer than
Figure 588352DEST_PATH_IMAGE041
2, single unit
Figure 368089DEST_PATH_IMAGE009
Simplification
This section is with right Doing for two steps simplifies.According to the analysis conclusion of the power response relevant treatment of double feedback electric engine under the sinusoidal excitation, can ignore input signal
Figure 108829DEST_PATH_IMAGE010
Effect, namely at transfer function matrix
Figure 462450DEST_PATH_IMAGE009
In, ignore
Figure 729483DEST_PATH_IMAGE011
With
Figure 459542DEST_PATH_IMAGE012
So will
Figure 656649DEST_PATH_IMAGE009
Be reduced to
Figure 650013DEST_PATH_IMAGE020
, have:
Figure 935501DEST_PATH_IMAGE014
（3）
For
Figure 938092DEST_PATH_IMAGE017
,
Figure 992636DEST_PATH_IMAGE010
Its output signal, correspondingly,
Figure 156901DEST_PATH_IMAGE017
Be reduced to
Figure 929685DEST_PATH_IMAGE018
, have:
Figure 470387DEST_PATH_IMAGE019
（4）
Therefore, Fig. 2 is reduced to Fig. 5 (a), and this simplifies for the first step.
Second step is simplified emphasis and considered: the access of wind energy turbine set is how to affect the synchronous motor low frequency oscillation mode.Only need to consider
Figure 379438DEST_PATH_IMAGE020
Steady-state characteristic under sinusoidal excitation can satisfy the needs that electromechanic oscillation mode is analyzed.Under the Stable State of Sine situation, the output quantity of wind-powered electricity generation unit
Figure 980183DEST_PATH_IMAGE021
With
Figure 974684DEST_PATH_IMAGE022
It all is sinusoidal phasor
Figure 319078DEST_PATH_IMAGE023
With
Figure 551476DEST_PATH_IMAGE024
At this moment,
Figure 854281DEST_PATH_IMAGE020
Effect be a complex gain matrix , wherein, ,
Figure 571068DEST_PATH_IMAGE027
So Fig. 5 (a) further is reduced to Fig. 5 (b).Electrical network subsystem after the simplification (
Figure 44774DEST_PATH_IMAGE018
) be two single-input single-output system (SISO system)s; The double-fed unit ( ) be single input two output systems.
3, the simplification interface of double-fed fan motor field and electrical network calculates
If the wind energy turbine set among Fig. 1 is installed n platform double-fed fan motor unit, because the upstream of this n platform unit is common bus B, when little value vibration appearred in conventional electrical network, for arbitrary unit, its input variable all was phasor
Figure 700063DEST_PATH_IMAGE028
Consider simultaneously the coordination of n platform unit, the power stage that wind energy turbine set is total
Figure 907054DEST_PATH_IMAGE021
With It is the phasor splicing of each double-fed unit power stage.Therefore, the reciprocation of whole n machine wind energy turbine set and electrical network can use Fig. 6 (a) to express.Therefore, following formula is arranged:
Figure 476893DEST_PATH_IMAGE029
（5）
Based on Fig. 6 (a), can realize containing the simplification analysis of the dynamo-electric pattern of electric system of double-fed fan motor field, step is as follows: the output of wind electric field that 1. calculates the equilibrium point place by trend
Figure 232359DEST_PATH_IMAGE030
With
Figure 559435DEST_PATH_IMAGE031
, with power be
Figure 109365DEST_PATH_IMAGE030
With
Figure 521892DEST_PATH_IMAGE031
Constant power source access bus B; 2. according to the mathematical model of separate unit double-fed fan motor unit, calculate after the linearization
Figure 81049DEST_PATH_IMAGE020
, and obtain the complex gain matrix
Figure 262632DEST_PATH_IMAGE032
(unit for model of the same race in the wind energy turbine set has identical Battle array is calculated one and is got final product).3. according to formula (5), the dynamic behaviour equivalence of wind energy turbine set is served as reasons
Figure 883286DEST_PATH_IMAGE028
The controlled power source of control, controlled power source and constant power source are combined as Equivalent Model, such as Fig. 6 (b).4. equivalent model is embedded conventional electrical network, obtain the merit angle rocking curve of synchronous motor by the microvariations time-domain-simulation, match obtains the leading electromechanical characteristics value of system through Prony again.
Below for the sample calculation analysis of the above-mentioned modeling method of the present invention:
In order to verify double-fed fan motor in this paper field Equivalent Model validity, utilize Matlab/Simulink to build example system shown in Figure 7.SG1 and SG2 are conventional synchronous generators, and WG is wind energy turbine set, and wind energy turbine set has comprised the double-fed unit of 6 two kinds of models.
Sample calculation analysis step: 1. utilize joint 3 described methods, form the Equivalent Model shown in Fig. 6 (b), successively access on bus B5, B6 and the B7, at the B3 place voltage is set and falls disturbance, obtain the power angle of synchronous generator rocking curve; 2. use wind energy turbine set detailed model (simulate in detail each unit electromagnetism, electromechanical transient, do not do any polymerization, equivalence or simplify), do and 1. identical l-G simulation test.Fig. 8 (a) and (b) and (c) based on the simulation result of above-mentioned two kinds of models, curve is that the merit angle rocking curve of SG1 is (owing to only be concerned about Perturbation Analysis, so only intercept the follow-up small size oscillating curve of disturbance among the figure; SG2 and the SG1 people having the same aspiration and interest is not so the SG2 power-angle curve provides).Find out that from the variation tendency of power-angle curve Equivalent Model and former detailed model have high consistency.By to all curve Prony Fitting Analysis among the figure, can obtain the leading dynamo-electric mode characteristic values under each situation again, gather
In table 1.
Table 1 is based on the dominant characteristics value of Equivalent Model and detailed model
The wind energy turbine set position Equivalent Model Detailed model
B6 -0.9968±j12.777 -1.0128±j12.767
B5 -1.0074±j12.689 -1.0431±j12.656
B7 -0.9985±j12.660 -1.0618±j12.620
As seen from Table 1, based on the very approaching analysis result based on detailed model of the dynamo-electric pattern analysis result of Equivalent Model.The eigenwert imaginary part error that causes behind the wind energy turbine set equivalence is all less than 0.5%; Wind energy turbine set is when B6, B5 and B7, and the real part error that causes is respectively 1.78%, 3.45% and 5.96%.As seen, Equivalent Model in this paper is not only highly succinct but also have a desirable precision.

Claims (2)

1. double-fed fan motor station equivalent modeling method that is suitable for the Study of Power System Small Disturbance stability analysis, the dot frequency that will be incorporated into the power networks is as sideband signal, with the power source of wind energy turbine set equivalence for controlled by mains frequency, so that the double-fed fan motor field becomes the Algebraic Constraint equation by the higher-dimension Simplification of Dynamic Model, the method may further comprise the steps:
Step (1), set up the microvariations power response model of double feedback electric engine under the sinusoidal excitation:
If only have a double-fed fan motor unit in the double-fed fan motor field, wherein and the mains frequency of site B be ω B(t), voltage magnitude is V B(t);
If ω B(t) little value sine-wave oscillation, i.e. Δ ω occur B(t)=Asin (ω gT), A is amplitude, ω gIt is oscillation frequency; Get A=0.01p.u., ω g=2.0Hz, the first squiggle of the meritorious output of unit;
Make again Δ V B(t)=Asin (ω gT), A=0.01p.u., ω g=2.0Hz, the second squiggle of the meritorious output of unit;
Wherein, the oscillation amplitude of the first squiggle is 60kW, and the oscillation amplitude of the second squiggle is 2kW;
Step (2) is with the transfer function matrix G of single unit d(s) simplify processing, further comprising the steps of:
The first step: according to the power response model of double feedback electric engine under the sinusoidal excitation, ignore input signal Δ V BEffect, namely at the transfer function matrix G of wind energy turbine set dynamic subsystem d(s) in, ignore G PV(s) and G QV(s), with G d(s) be reduced to
Figure FDA00002067059400012
Wherein, G P ω(s) the be incorporated into the power networks input of bus frequency of representative is take the meritorious transport function for exporting of blower fan; G Q ω(s) expression is with the input of the bus frequency that is incorporated into the power networks, take the idle transport function as output of blower fan; G PV(s) expression is take the busbar voltage that is incorporated into the power networks as input, take the meritorious transport function as output of blower fan; G QV(s) expression is take the busbar voltage that is incorporated into the power networks as input, take the idle transport function as output of blower fan;
Transfer function matrix G for conventional electrical network dynamic subsystem s(s), Δ V BIts output signal, correspondingly, G s(s) be reduced to
Figure FDA00002067059400021
Have:
Figure FDA00002067059400022
Second step: consider
Figure FDA00002067059400023
Steady-state characteristic under sinusoidal excitation, under the Stable State of Sine situation, the output quantity Δ P of wind-powered electricity generation unit BWith Δ Q BIt all is sinusoidal phasor
Figure FDA00002067059400024
With
Figure FDA00002067059400025
Effect be a complex gain matrix K=e J φ| K|, wherein,
Figure FDA00002067059400026
Figure FDA00002067059400027
Electrical network subsystem after the simplification
Figure FDA00002067059400028
Be two single-input single-output system (SISO system)s; The double-fed unit
Figure FDA00002067059400029
Be single input two output systems;
Step (3), the simplification interface of double-fed fan motor field and electrical network is processed, and comprising:
If wind energy turbine set is installed n platform double-fed fan motor unit, the upstream of this n platform unit is common bus B, and when little value vibration appearred in conventional electrical network, for arbitrary unit, its input variable all was phasor Δ ω BIn conjunction with the coordination of n platform unit, the power stage Δ P that wind energy turbine set is total BWith Δ Q BIt is the phasor splicing of each double-fed unit power stage; The reciprocation expression formula of whole n machine wind energy turbine set and electrical network is:
Figure FDA000020670594000210
Based on the reciprocation of described wind energy turbine set and electrical network, realization contains the simplification analysis of the dynamo-electric pattern of electric system of double-fed fan motor field, and step is as follows:
The first step: the output of wind electric field P that calculates the equilibrium point place by trend B0And Q B0, be P with power B0And Q B0Constant power source access bus B;
Second step according to the mathematical model of separate unit double-fed fan motor unit, calculates after the linearization
Figure FDA00002067059400031
And obtain the complex gain matrix K;
The 3rd step is according to formula
Figure FDA00002067059400032
Be by Δ ω with the dynamic behaviour equivalence of wind energy turbine set BThe controlled power source of control, controlled power source and constant power source are combined as Equivalent Model;
The 4th step embedded conventional electrical network with equivalent model, obtained the merit angle rocking curve of synchronous motor by the microvariations time-domain-simulation, and match obtains the leading electromechanical characteristics value of system through Prony again.
2. the double-fed fan motor station equivalent modeling method that is suitable for the Study of Power System Small Disturbance stability analysis as claimed in claim 1 is characterized in that, described mathematical model according to separate unit double-fed fan motor unit is calculated after the linearization
Figure FDA00002067059400033
And obtain the step of complex gain matrix K, the unit of model of the same race has identical K battle array in wind energy turbine set, only calculates a double-fed fan motor unit.
CN 201110070510 2011-03-23 2011-03-23 Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system Expired - Fee Related CN102136036B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110070510 CN102136036B (en) 2011-03-23 2011-03-23 Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110070510 CN102136036B (en) 2011-03-23 2011-03-23 Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system

Publications (2)

Publication Number Publication Date
CN102136036A CN102136036A (en) 2011-07-27
CN102136036B true CN102136036B (en) 2013-01-16

Family

ID=44295820

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110070510 Expired - Fee Related CN102136036B (en) 2011-03-23 2011-03-23 Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system

Country Status (1)

Country Link
CN (1) CN102136036B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102682222B (en) * 2012-05-23 2014-03-12 甘肃省电力公司电力科学研究院 Continuous tide calculation method based on wind power fluctuation rule
CN103279641A (en) * 2013-04-24 2013-09-04 广东电网公司电力科学研究院 Method for conducting multi-modal dynamic equivalence on complex electrical power systems except for generators
CN103441522B (en) * 2013-07-04 2015-08-12 贵州大学 In automatic regulating wind power field blower fan gain merit export method
CN103475009B (en) * 2013-08-02 2016-08-17 李啸骢 A kind of Multi-index nonlinear control method for direct-drive permanent-magnetism wind turbine current transformer
CN103956767B (en) * 2014-02-21 2016-04-20 国家电网公司 A kind of wind farm grid-connected method for analyzing stability considering wake effect
CN103777525B (en) * 2014-02-26 2016-01-13 华北电力大学 The self defined interface of wind energy turbine set replicating machine and RTDS emulator
CN104167735B (en) * 2014-08-29 2016-02-10 河海大学 A kind of non-mechanism equivalent modeling method of wind energy turbine set and device
CN105140957B (en) * 2015-08-20 2017-10-10 华北电力大学 Electromechanic oscillation mode evaluation method based on wind power plant and photovoltaic plant polymerization model
CN105184418B (en) * 2015-09-23 2019-04-02 重庆大学 Based on the equivalent interconnected network optimal load flow calculation method of sensitivity and constraint
CN105375513B (en) * 2015-11-06 2017-08-25 国家电网公司 A kind of 110 kilovolts of wind power field automatic voltage controls based on real-time on-line equivalent
CN108616125A (en) * 2018-05-16 2018-10-02 天津大学 Contain asynchronous wind power plant Study of Power System Small Disturbance Method of Stability Analysis based on security domain
CN108964122B (en) * 2018-07-16 2021-01-15 中国电力科学研究院有限公司 Verification method and device for frequency response characteristic of wind turbine generator electrical model
CN111382525B (en) * 2020-03-24 2024-04-05 国网河北省电力有限公司沧州供电分公司 Low-frequency power transmission closed-loop test system and method
CN113270878B (en) * 2021-07-06 2022-08-30 华北电力大学 Resonance analysis method and system for wind power plant grid-connected system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101661530A (en) * 2009-09-30 2010-03-03 华北电力大学（保定） Method for acquiring steady-state equivalent wind speed and generated power in wind power station based on correlation analysis
CN101882167A (en) * 2009-11-25 2010-11-10 中国电力科学研究院 Wind power station equivalent modeling method of large-scale wind power concentration access power grid
CN101937483A (en) * 2010-09-09 2011-01-05 苏勋文 Dynamic equivalence modeling method for double-fed set wind power station based on action condition of pitch angle

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101661530A (en) * 2009-09-30 2010-03-03 华北电力大学（保定） Method for acquiring steady-state equivalent wind speed and generated power in wind power station based on correlation analysis
CN101882167A (en) * 2009-11-25 2010-11-10 中国电力科学研究院 Wind power station equivalent modeling method of large-scale wind power concentration access power grid
CN101937483A (en) * 2010-09-09 2011-01-05 苏勋文 Dynamic equivalence modeling method for double-fed set wind power station based on action condition of pitch angle

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
F. Michael Hughes et al..A Power System Stabilizer for DFIG-Based Wind Generation.《IEEE TRANSACTIONS ON POWER SYSTEMS》.2006,第21卷(第2期), *
J. Conroy et al..Aggregate modelling of wind farms containing.《IET Renewable Power Generation》.2009,第3卷(第1期),39–52. *
J.G. Slootweg et al..Aggregated Modelling of Wind Parks in Plower.《2003 IEEE Bologna PowerTech Conference》.2003, *

Also Published As

Publication number Publication date
CN102136036A (en) 2011-07-27

Similar Documents

Publication Publication Date Title
CN102136036B (en) Double-feed wind power station equivalent modeling method applied to analysis on small signal stability of power system
CN101937483B (en) Dynamic equivalence modeling method for double-fed set wind power station based on action condition of pitch angle
CN103761349B (en) A kind of wind energy turbine set equivalent modeling method based on Wind turbines probability same tone
CN102760191B (en) Double-fed machine set wind power station equivalent modeling system and method based on rotation speed grouping
CN103473478B (en) Power Network Transient Stability appraisal procedure based on energy function
CN102012956A (en) Wind farm equivalent method based on wind farm input wind speed and wind direction chance fluctuation
CN103558768B (en) A kind of based on the equivalent modeling method of wind speed distribution characteristics in wind energy turbine set
CN103777525B (en) The self defined interface of wind energy turbine set replicating machine and RTDS emulator
CN105703382A (en) Subsynchronous oscillation suppression method and system for wind power plant
CN103094920A (en) Equivalence method of direct-drive-type wind turbine generator wind power plant
CN107666155A (en) System Stochastic Stability Analysis method of providing multiple forms of energy to complement each other based on Markov model
CN102034008B (en) Variable-parameter capacitance compensation equivalent modeling method for wind power station of constant-speed wind power unit at fluctuating air speed
Gao et al. Impact of DFIG‐based wind farm integration on sub‐synchronous torsional interaction between HVDC and thermal generators
CN106208158A (en) The Inertia Matching method of the most micro-source paired running in microgrid
CN103248049A (en) Method of obtaining electrical power system voltage stability domain tangent plane containing DFIG (Doubly Fed Induction Generator) wind power plant
CN107769227B (en) A kind of wind-powered electricity generation station equivalent modeling method suitable for subsynchronous research
Eisa et al. Time domain study of a type-3 dfig wind turbine's dynamics: Q drop function effect and attraction vs control limits analysis
CN110210170B (en) Modeling method for equivalent small signal model of large-scale wind turbine group
Yao et al. RBF neural network based self-tuning PID pitch control strategy for wind power generation system
CN104200027A (en) Wind power plant equivalent modeling method and device with low voltage ride through characteristic considered
Yan et al. Transient modelling of doubly‐fed induction generator based wind turbine on full operation condition and rapid starting period based on low voltage ride‐through testing
Thomas et al. Event analysis and real‐time validation of doubly fed induction generator‐based wind energy system with grid reactive power exchange under sub‐synchronous and super‐synchronous modes
Ni et al. Cross‐gramian‐based dynamic equivalence of wind farms
Wang et al. Wind farm model with DFIG for small signal stability study
CN109657380A (en) A kind of double-fed fan motor field Dynamic Equivalence based on Extended Kalman filter

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130116

Termination date: 20210323