CN108418242A - A kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest - Google Patents
A kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest Download PDFInfo
- Publication number
- CN108418242A CN108418242A CN201810064146.0A CN201810064146A CN108418242A CN 108418242 A CN108418242 A CN 108418242A CN 201810064146 A CN201810064146 A CN 201810064146A CN 108418242 A CN108418242 A CN 108418242A
- Authority
- CN
- China
- Prior art keywords
- double
- fan motor
- motor machine
- fed fan
- indicate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000010355 oscillation Effects 0.000 claims abstract description 107
- 230000001052 transient effect Effects 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 31
- 230000001427 coherent effect Effects 0.000 claims abstract description 23
- 238000013016 damping Methods 0.000 claims abstract description 20
- 230000008569 process Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000006116 polymerization reaction Methods 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000011160 research Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000002776 aggregation Effects 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- H02J3/386—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest, the described method comprises the following steps：Step S10, the output power of double-fed fan motor machine is fitted to the linear combination of amplitude, phase, frequency and the exponential function of damping using prony algorithms；Step S20, the linear combination obtained according to step S10 determines the similar oscillation mode between double-fed fan motor machine using the theory of similarity；Step S30, people having the same aspiration and interest similarity quantizating index between double-fed fan motor machine is defined with similar oscillation mode, differentiates the Coherent Generator Group in wind power plant；Step S40, it obtains etc. the equivalent transient state built-in potential of check-ins using the active power of double-fed fan motor machine output as weight, and the equivalent transient state built-in potential busbar of the double-fed fan motor machine of same group is passed through into the equivalent transient state built-in potential busbar for the check-ins such as multiple no-load voltage ratio phase-shifting transformer is connected to respectively；Step S50, the equivalent parameters for calculating the check-ins such as each will wait the equivalent transient state built-in potential busbar of check-ins to be connected to equal check-ins by equivalent double-fed fan motor machine impedance.
Description
Technical field
The present invention relates to the wind power plant equivalent modeling technical fields in new energy power generation technology, and in particular to one kind being based on phase
Like the double-fed fan motor machine Dynamic Equivalence of the degree people having the same aspiration and interest.
Background technology
With the rapid development of modern electric power industry, extensive alternating current-direct current interconnected electric power system occurs, multimachine (thousands of
Machine), bulk power grid (thousands of circuits, thousands of busbares) feature it is more and more prominent, this makes planning and design of power system and operation
The calculating of mode becomes abnormal very complicated.On the one hand, it completes these in detail and calculates to need to take a substantial amount of time and space,
On the other hand, it is limited by simulation scale and hardware device, above-mentioned calculating is difficult to realize in practice.It is needed thus to original
System carries out certain and reasonably simplifies, to reduce system scale.In practice, to the dynamic property research of a large power system
General only most interested to some region, which is referred to as research system, and the region with this region farther out, in research only
Count and its influence to survey region, inside need not describe in detail, depression of order and simplification can be often done, to save the people of research
Power and material resources, the region which simplifies are known as external system.The external system for not needing detailed analysis is carried out equivalent simplification
It is very necessary.
The new energy that wind-powered electricity generation is developed as most development potentiality and scale, just receives the world from the 1970s
The extensive concern of various countries.With the fast development of wind-powered electricity generation, all there is fulminant growth in scale and number in wind power plant, because
This, the grid-connected influence to electric system in large-scale wind power field is also growing day by day.In research large-scale wind power access to electric system
When the influence of dynamic characteristic, the farm model of an energy accurate characterization wind power plant overall dynamics characteristic is needed.Wind power plant is usual
It is made of tens even up to a hundred Wind turbines, wherein the model of each double-fed wind turbine reaches more than ten ranks, using detailed
Model carries out modeling and simulating efficiency to wind power plant will be very low, it is therefore necessary to simplify polymerizeing to large-scale wind power field and build
Mould.
Currently, usually by Equivalent Network to higher voltage class in bulk power grid, and in low voltage grade
Wind power plant equivalence is that an active power is negative load, such to handle the dynamic characteristic for having neglected wind power plant.And with
The fast development of wind power technology, single-machine capacity have reached MW class, and wind-power electricity generation capacity is shared in power grid total power generating capacity
Ratio is increasing, and influence of the dynamic characteristic to operation of power networks characteristic is also increasingly notable, when carrying out dynamic equivalent, wind power plant
The characteristic in the area of concentration cannot be neglected simply again, therefore, be badly in need of the equivalent simplification method of further investigation wind power plant.
It is the important content of wind power plant dynamic equivalent to divide group and parameter aggregation.Currently, how to wind power plant people having the same aspiration and interest characteristic into
The pertinent literature that row differentiates is less, and pertinent literature proposes the method according to characteristic quantity coherence grouping, but such method has spy
Sign amount choose whether complete problem, this method with the presence or absence of universality be hinder its application bottleneck.Based on disturbed path
Similarity divide group to be to divide the common method of group to multi-computer system, often using the generator rotor angle of generator as observed quantity, and with synchronous electricity
Unlike machine, due to the presence of power electronic equipment in double-fed wind turbine so that its external dynamic characteristic and conventional synchronization hair
Motor difference is larger, therefore when carrying out people having the same aspiration and interest differentiation to wind turbine, for being compared using the disturbed path of which variable
There are no build consensus.Further, since double-fed wind turbine uses decoupling control, therefore inertia is not externally embodied, to the shadow of system
Ring only with send out power number characterize, rarely have the pertinent literature of the technique study for wind turbine parameter aggregation at present.
Therefore, it is badly in need of a kind of people having the same aspiration and interest suitable for wind power plant to differentiate and the Dynamic Equivalence of parameter aggregation, to meet containing advising greatly
The dynamic equivalent demand of the electric system of mould wind power plant.
Invention content
The purpose of the present invention is in view of the above shortcomings of the prior art, provide a kind of double-fed wind based on the similarity people having the same aspiration and interest
Motor dynamics equivalence method, to realize, the accurate of double-fed fan motor machine Coherent Generator Group divides and moves during wind power plant equivalent modeling
The reasonable polymerization of state parameter can be greatly reduced the scale of Wind turbines and retain the dynamic characteristic of Wind turbines.
The purpose of the present invention can be achieved through the following technical solutions：
A kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest, the described method comprises the following steps：
Step S10, the output power of double-fed fan motor machine is fitted to amplitude, phase, frequency and damping using prony algorithms
Exponential function linear combination；
Step S20, the linear combination obtained according to step S10 is determined similar between double-fed fan motor machine using the theory of similarity
Oscillation mode；
Step S30, people having the same aspiration and interest similarity quantizating index between double-fed fan motor machine is defined with similar oscillation mode, differentiates wind power plant
In Coherent Generator Group；
Step S40, the equivalent transient state built-in potential of check-ins is obtained etc. using the active power of double-fed fan motor machine output as weight, and
By the equivalent transient state built-in potential busbar of the double-fed fan motor machine of same group respectively by the check-ins such as multiple no-load voltage ratio phase-shifting transformer is connected to etc.
Imitate transient state built-in potential busbar；
Step S50, the equivalent parameters for calculating the check-ins such as each passes through the equivalent transient state of check-ins such as equivalent double-fed fan motor machine impedance general
The check-ins such as built-in potential busbar is connected to.
Further, the linear combination that step S10 is obtained is as follows：
Wherein,Indicate the approximation of X (n), bi=Aiexp(jθi), Zi=exp ((αi+j2πfi) Δ t), AiIt indicates
Amplitude, the θ of double-fed fan motor machine output poweriIndicate initial phase angle, the f of double-fed fan motor machine output poweriIndicate the output of double-fed fan motor machine
The frequency of oscillation of power, αiIndicate that sampling interval, p indicate arbitrary amplitude, phase, frequency for the decay factor of frequency of oscillation, Δ t
With the number of the linear combination of the exponential function of damping, J indicates the exponent number that prony algorithms decompose.
Further, the detailed process of step S20 is：The ENERGY E of double-fed fan motor machine i-th of oscillation mode of output poweri
It is obtained by the quadratic sum of each sampled point modulus value under the oscillation mode：
Ei=Ai exp(jθi)·exp(n·Δt·(αi+j2πfi)) (2)
Wherein, n=0,1 ..., J-1, J indicate the exponent number that prony algorithms decompose；
An energy threshold v is set, the numerical value of energy threshold v should be higher than that 99%, and concrete numerical value can be according to people having the same aspiration and interest precision
It is required that setting, all oscillation modes are sorted from high to low according to energy, and added up successively until gained energy is not less than total energy
AmountWith the product of energy threshold v：
After redundancy of the above method elimination before and after equivalent in oscillation mode, it can obtain comprising most information
Equivalence before l rank oscillation modes, l be much smaller than prony algorithms decompose exponent number J, simplify the searching process of similar oscillation mode；
According to frequency, damping, the principle of energy three Euclidean distance minimum in characteristic information, with a wherein double-fed fan motor
On the basis of machine output power curve l rank oscillation modes, found in other double-fed fan motor machine output power curve m rank oscillation modes
Similar oscillation mode, for i=1,2 ..., l, k=1,2 ..., m are calculated：
Wherein, fAiIndicate the frequency with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, fBkIndicate other
The frequency of double-fed fan motor machine k-th of oscillation mode of output power curve, αAiIndicate fAiDecay factor, αBkIndicate fBkDecaying
The factor, EAiIndicate the energy with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, EBkIndicate other double-fed fan motors
The energy of k-th of oscillation mode of machine output power curve calculates for oscillation mode i and obtains minimum dik, then double-fed wind is referred to
Output power of motor curve oscillation mode i is a pair of similar oscillation to other double-fed fan motor machine output power curve oscillation modes k
Pattern removes oscillation mode k in oscillation mode, is used for the searching process of the similar oscillation mode of lower a pair, finally obtains l to phase
Like oscillation mode.
Further, the detailed process of step S30 is：The similar oscillation of l ranks of each double-fed fan motor machine is obtained according to step S20
After pattern, the frequency similarity q (f between i-th of similar oscillation modei) and damping similarity q (αi) be respectively：
Wherein, fAiIndicate the frequency with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, fBiIndicate other
The frequency of double-fed fan motor machine i-th of oscillation mode of output power curve, αAiIndicate fAiDecay factor, αBiIndicate fBiDecaying
The factor considers the similarity weight of frequency and damping, therefore the similarity between i-th of similar oscillation mode is：
Wherein a, b indicate frequency and the similarity weight of damping respectively；
In oscillation mode the ratio of the energy of each merotype and gross energy as the merotype energy in similarity calculation
In account for the weights of whole oscillation mode energy：
The similarity of all similar oscillation modes is after consideration energy weight：
Wherein EiThe energy for indicating double-fed fan motor machine i-th of oscillation mode of output power calculates same between double-fed fan motor machine
After phase modulation likelihood metric index, the Coherent Generator Group in wind power plant is determined.
Further, the detailed process of step S40 is：First, by the stator of each double-fed fan motor machine, rotor voltage side
Journey indicates as follows：
Wherein, u, i, ψ indicate the voltage, electric current and magnetic linkage of double-fed fan motor machine respectively；R, L indicates double-fed fan motor machine respectively
Resistance and inductance；S, r in subscript indicate the stator of double-fed fan motor machine, rotor respectively；D, q in subscript indicate respectively d-axis,
Quadrature axis component；LmIndicate the mutual inductance of double-fed fan motor machine；P is differential operator；S indicates revolutional slip；
In formula,It enables：
Wherein, Xss＞ ＞ Rs, so the stator pressure drop of double-fed fan motor machine can ignore that and disregard, it is equivalent to define double-fed fan motor machine
Transient state built-in potential is as follows：
E'=Us-jXssIs (12)
Wherein, UsIndicate the set end voltage of double-fed fan motor machine, IsIndicate the stator current of double-fed fan motor machine, it is super that j represents phase
Preceding 90 degree of symbol；
Since double-fed fan motor machine is using decoupling control, inertia is not externally embodied, double-fed fan motor machine is to system
It influences to characterize only by the power of its injection, uses the polymerization using output power as weight, if double in Coherent Generator Group
The number for presenting wind turbine is N, waits the equivalent transient state built-in potential of check-insFor：
Wherein, piIndicate the active power of i-th of double-fed fan motor machine in Coherent Generator Group；The equivalent transient state built-in potential of equal check-ins
It is determined by double-fed fan motor machine all in Coherent Generator Group, according to the external characteristics of double-fed fan motor machine, the shared power of each double-fed fan motor machine
Their output power is depended on again；
Wherein, the no-load voltage ratio of the multiple no-load voltage ratio phase-shifting transformer is：
Indicate the equivalent transient state built-in potential of i-th of double-fed fan motor machine in Coherent Generator Group, and by the double-fed fan motor machine of same group
Equivalent transient state built-in potential busbar passes through the equivalent transient state built-in potential busbar for the check-ins such as multiple no-load voltage ratio phase-shifting transformer is connected to respectively.
Further, in step S50, the equivalent parameters of each equal check-ins is：
Wherein, N indicates that the number of double-fed fan motor machine in Coherent Generator Group, P indicate that the active power of double-fed fan motor machine, Q indicate
The reactive power of double-fed fan motor machine, XsIndicate the stator reactance of double-fed fan motor machine, XrIndicate the rotor reactance of double-fed fan motor machine, Hg
Indicate the rotor inertia time constant of perceptual generator in double-fed fan motor machine, HtIndicate that the rotor of wind energy conversion system in double-fed fan motor machine is used
Property time constant, KsIndicate the axis rigidity coefficient of wind energy conversion system in double-fed fan motor machine, the eq in subscript, which is respectively represented, waits check-ins ginseng
Number.
Compared with prior art, the present invention having the following advantages that and advantageous effect：
1, a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest of the invention, can be effectively to advising greatly
Mould double-fed fan motor machine carries out Coherent Generator Group division, and carries out Reasonable Parameters polymerization to people having the same aspiration and interest unit.This method can retain wind
The scale that Wind turbines are significantly reduced under the premise of motor group dynamic characteristic is closed to simplify perimeter to give top priority to what is the most important
The area to be studied of note.This method first with prony algorithms by Power Output for Wind Power Field be fitted to amplitude, phase, frequency and
The linear combination of the exponential function of damping determines similar oscillation mode followed by the theory of similarity, is determined with similar oscillation mode
People having the same aspiration and interest similarity quantizating index between adopted wind turbine differentiates the people having the same aspiration and interest unit in wind power plant；Then it is exported with double-fed fan motor machine active
Power is the equivalent transient state built-in potential for the check-ins such as weight obtains, and the double-fed fan motor machine transient state built-in potential busbar of same group is led to respectively
It crosses no-load voltage ratio phase-shifting transformer again and is connected to equivalent transient state built-in potential busbar；The equivalent parameters for finally calculating each equivalent unit, passes through
Equivalent double-fed fan motor machine impedance is by the check-ins such as equivalent transient state built-in potential is connected to.In this way, meeting the unit containing large-scale wind power
Bulk power grid high-precision equivalent simplification demand, the polymerization can more reflect the true dynamic characteristic of wind power plant, also be science, correct
The ground analysis grid-connected influence to bulk power grid dynamic characteristic in large-scale wind power field is laid a good foundation.
It is 2, of the invention that dynamic equivalent is carried out to wind power plant based on the double-fed fan motor machine Dynamic Equivalence of the similarity people having the same aspiration and interest,
The scale of Wind turbines can be greatly reduced and retain the power characteristic of Wind turbines, can more reflect the true of wind power plant after equivalence
Characteristic, to meet the bulk power grid high-precision equivalent simplification demand of the unit containing large-scale wind power.
Description of the drawings
Fig. 1 is the flow chart of double-fed fan motor machine Dynamic Equivalence of the embodiment of the present invention based on the similarity people having the same aspiration and interest.
The step of Fig. 2 is the check-ins such as construction of the embodiment of the present invention is schemed.
Specific implementation mode
Present invention will now be described in further detail with reference to the embodiments and the accompanying drawings, but embodiments of the present invention are unlimited
In this.
Embodiment：
Present embodiments provide a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest, the stream of the method
Journey figure is as shown in Figure 1, include the following steps：
Step S10, it is amplitude, phase, frequency by the characteristics of output power curve matching of double-fed fan motor machine using prony algorithms
The linear combination of rate and the exponential function of damping；Linear combination is as follows：
Wherein,Indicate the approximation of X (n), bi=Aiexp(jθi), Zi=exp ((αi+j2πfi) Δ t), AiIt indicates
Amplitude, the θ of double-fed fan motor machine output poweriIndicate initial phase angle, the f of double-fed fan motor machine output poweriIndicate the output of double-fed fan motor machine
The frequency of oscillation of power, αiIndicate that sampling interval, p indicate arbitrary amplitude, phase, frequency for the decay factor of frequency of oscillation, Δ t
With the number of the linear combination of the exponential function of damping, J indicates the exponent number that prony algorithms decompose；
Step S20, the linear combination obtained according to step S10 is determined similar between double-fed fan motor machine using the theory of similarity
Oscillation mode；Detailed process is：The ENERGY E of double-fed fan motor machine i-th of oscillation mode of output poweriBy respectively being adopted under the oscillation mode
The quadratic sum of sampling point modulus value obtains：
Ei=Aiexp(jθi)·exp(n·Δt·(αi+j2πfi)) (2)
Wherein, n=0,1 ..., J-1, J indicate the exponent number that prony algorithms decompose；
An energy threshold v is set, the numerical value of energy threshold v should be higher than that 99%, and concrete numerical value can be according to people having the same aspiration and interest precision
It is required that setting, will be set as 99.9% by the numerical value of energy threshold v herein, all oscillation modes sort from high to low according to energy,
And it is added up successively until gained energy is not less than gross energyWith the product of energy threshold v：
After redundancy of the above method elimination before and after equivalent in oscillation mode, it can obtain comprising most information
Equivalence before l rank oscillation modes, l be much smaller than prony algorithms decompose exponent number J, simplify the searching process of similar oscillation mode；
According to frequency, damping, the principle of energy three Euclidean distance minimum in characteristic information, with a wherein double-fed fan motor
On the basis of machine output power curve l rank oscillation modes, found in other double-fed fan motor machine output power curve m rank oscillation modes
Similar oscillation mode, for i=1,2 ..., l, k=1,2 ..., m are calculated：
Wherein, fAiIndicate the frequency with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, fBkIndicate other
The frequency of double-fed fan motor machine k-th of oscillation mode of output power curve, αAiIndicate fAiDecay factor, αBkIndicate fBkDecaying
The factor, EAiIndicate the energy with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, EBkIndicate other double-fed fan motors
The energy of k-th of oscillation mode of machine output power curve calculates for oscillation mode i and obtains minimum dik, then double-fed wind is referred to
Output power of motor curve oscillation mode i is a pair of similar oscillation to other double-fed fan motor machine output power curve oscillation modes k
Pattern removes oscillation mode k in oscillation mode, is used for the searching process of the similar oscillation mode of lower a pair, finally obtains l to phase
Like oscillation mode.
Step S30, people having the same aspiration and interest similarity quantizating index between double-fed fan motor machine is defined with similar oscillation mode, differentiates wind power plant
In Coherent Generator Group；Detailed process is：After the similar oscillation mode of l ranks for obtaining each double-fed fan motor machine according to step S20, i-th
Frequency similarity q (f between similar oscillation modei) and damping similarity q (αi) be respectively：
Wherein, fAiIndicate the frequency with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, fBiIndicate other
The frequency of double-fed fan motor machine i-th of oscillation mode of output power curve, αAiIndicate fAiDecay factor, αBiIndicate fBiDecaying
The factor considers the similarity weight of frequency and damping, therefore the similarity between i-th of similar oscillation mode is：
Wherein a, b indicate frequency and the similarity weight of damping respectively, think frequency and the similarity weight of damping herein
Equal, the similarity obtained between i-th of similar oscillation mode is：
q(λi)=0.5q (fi)+0.5q(αi)
In oscillation mode the ratio of the energy of each merotype and gross energy as the merotype energy in similarity calculation
In account for the weights of whole oscillation mode energy：
The similarity of all similar oscillation modes is after consideration energy weight：
Wherein EiThe energy for indicating double-fed fan motor machine i-th of oscillation mode of output power calculates same between double-fed fan motor machine
After phase modulation likelihood metric index, the Coherent Generator Group in wind power plant is determined.
Step S40, the equivalent transient state built-in potential of check-ins is obtained etc. using the active power of double-fed fan motor machine output as weight, and
By the equivalent transient state built-in potential busbar of the double-fed fan motor machine of same group respectively by the check-ins such as multiple no-load voltage ratio phase-shifting transformer is connected to etc.
Imitate transient state built-in potential busbar；Detailed process is：First, the stator of each double-fed fan motor machine, rotor voltage equation are indicated such as
Under：
Wherein, u, i, ψ indicate the voltage, electric current and magnetic linkage of double-fed fan motor machine respectively；R, L indicates double-fed fan motor machine respectively
Resistance and inductance；S, r in subscript indicate the stator of double-fed fan motor machine, rotor respectively；D, q in subscript indicate respectively d-axis,
Quadrature axis component；LmIndicate the mutual inductance of double-fed fan motor machine；P is differential operator；S indicates revolutional slip；
In formula,It enables：
Wherein, Xss＞ ＞ Rs, so the stator pressure drop of double-fed fan motor machine can ignore that and disregard, it is equivalent to define double-fed fan motor machine
Transient state built-in potential is as follows：
E'=Us-jXssIs (12)
Wherein, UsIndicate the set end voltage of double-fed fan motor machine, IsIndicate the stator current of double-fed fan motor machine, it is super that j represents phase
Preceding 90 degree of symbol；
Since double-fed fan motor machine is using decoupling control, inertia is not externally embodied, double-fed fan motor machine is to system
It influences to characterize only by the power of its injection, uses the polymerization using output power as weight, if double in Coherent Generator Group
The number for presenting wind turbine is N, waits the construction flow of check-ins as shown in Fig. 2, waiting the equivalent transient state built-in potential of check-insFor：
Wherein, piIndicate the active power of i-th of double-fed fan motor machine in Coherent Generator Group；The equivalent transient state built-in potential of equal check-ins
It is determined by double-fed fan motor machine all in Coherent Generator Group, according to the external characteristics of double-fed fan motor machine, the shared power of each double-fed fan motor machine
Their output power is depended on again；
Wherein, the no-load voltage ratio of the multiple no-load voltage ratio phase-shifting transformer is：
Indicate the equivalent transient state built-in potential of i-th of double-fed fan motor machine in Coherent Generator Group, and by the double-fed fan motor machine of same group
Equivalent transient state built-in potential busbar passes through the equivalent transient state built-in potential busbar for the check-ins such as multiple no-load voltage ratio phase-shifting transformer is connected to respectively；
Step S50, the equivalent parameters for calculating the check-ins such as each passes through the equivalent transient state of check-ins such as equivalent double-fed fan motor machine impedance general
The check-ins such as built-in potential busbar is connected to, wherein the equivalent parameters of each equal check-ins is：
Wherein, N indicates that the number of double-fed fan motor machine in Coherent Generator Group, P indicate that the active power of double-fed fan motor machine, Q indicate
The reactive power of double-fed fan motor machine, XsIndicate the stator reactance of double-fed fan motor machine, XrIndicate the rotor reactance of double-fed fan motor machine, Hg
Indicate the rotor inertia time constant of perceptual generator in double-fed fan motor machine, HtIndicate that the rotor of wind energy conversion system in double-fed fan motor machine is used
Property time constant, KsIndicate the axis rigidity coefficient of wind energy conversion system in double-fed fan motor machine, the eq in subscript, which is respectively represented, waits check-ins ginseng
Number.
The above, patent preferred embodiment only of the present invention, but the protection domain of patent of the present invention is not limited to
This, any one skilled in the art is in the range disclosed in patent of the present invention, according to the skill of patent of the present invention
Art scheme and its patent of invention design are subject to equivalent substitution or change, belong to the protection domain of patent of the present invention.
Claims (6)
1. a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest, which is characterized in that the method includes following
Step：
Step S10, the output power of double-fed fan motor machine is fitted to amplitude, phase, frequency and the finger of damping using prony algorithms
The linear combination of number function；
Step S20, the linear combination obtained according to step S10 determines the similar oscillation between double-fed fan motor machine using the theory of similarity
Pattern；
Step S30, people having the same aspiration and interest similarity quantizating index between double-fed fan motor machine is defined with similar oscillation mode, differentiated in wind power plant
Coherent Generator Group；
Step S40, the equivalent transient state built-in potential of check-ins is obtained etc. using the active power of double-fed fan motor machine output as weight, and will be same
The equivalent transient state built-in potential busbar of double-fed fan motor machine of group respectively by the way that multiple no-load voltage ratio phase-shifting transformer is connected to equal check-ins it is equivalent temporary
State built-in potential busbar；
Step S50, the equivalent parameters for calculating the check-ins such as each passes through electricity in the equivalent transient state of check-ins such as equivalent double-fed fan motor machine impedance general
The check-ins such as gesture busbar is connected to.
2. a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest according to claim 1, feature exist
In the linear combination that step S10 is obtained is as follows：
Wherein,Indicate the approximation of X (n), bi=Aiexp(jθi), Zi=exp ((αi+j2πfi) Δ t), AiIndicate double-fed
Amplitude, the θ of wind turbine output poweriIndicate initial phase angle, the f of double-fed fan motor machine output poweriIndicate double-fed fan motor machine output power
Frequency of oscillation, αiIndicate that sampling interval, p indicate arbitrary amplitude, phase, frequency and resistance for the decay factor of frequency of oscillation, Δ t
The number of the linear combination of the exponential function of Buddhist nun, J indicate the exponent number that prony algorithms decompose.
3. a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest according to claim 1, feature exist
In the detailed process of step S20 is：The ENERGY E of double-fed fan motor machine i-th of oscillation mode of output poweriBy under the oscillation mode
The quadratic sum of each sampled point modulus value obtains：
Ei=Aiexp(jθi)·exp(n·Δt·(αi+j2πfi)) (2)
Wherein, n=0,1 ..., J-1, J indicate the exponent number that prony algorithms decompose；
An energy threshold v is set, the numerical value of energy threshold v should be higher than that 99%, and concrete numerical value can be according to people having the same aspiration and interest required precision
Setting, all oscillation modes are sorted from high to low according to energy, and are added up successively until gained energy is not less than gross energyWith the product of energy threshold v：
Using the above method eliminate it is equivalent before and after after redundancy in oscillation mode, can obtain comprising most information etc.
L rank oscillation modes before value, l are much smaller than the exponent number J that prony algorithms decompose, simplify the searching process of similar oscillation mode；
It is defeated with a wherein double-fed fan motor machine according to frequency, damping, the principle of energy three Euclidean distance minimum in characteristic information
Go out on the basis of power curve l rank oscillation modes, is found in other double-fed fan motor machine output power curve m rank oscillation modes similar
Oscillation mode, for i=1,2 ..., l, k=1,2 ..., m are calculated：
Wherein, fAiIndicate the frequency with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, fBkIndicate other double-fed wind
The frequency of k-th of oscillation mode of output power of motor curve, αAiIndicate fAiDecay factor, αBkIndicate fBkDecay factor,
EAiIndicate the energy with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, EBkIndicate that other double-fed fan motor machines export
The energy of k-th of oscillation mode of power curve calculates for oscillation mode i and obtains minimum dik, then it is defeated to refer to double-fed fan motor machine
It is a pair of similar oscillation mode to go out power curve oscillation mode i to other double-fed fan motor machine output power curve oscillation modes k,
Oscillation mode k is removed in oscillation mode, is used for the searching process of the similar oscillation mode of lower a pair, is finally obtained l to similar oscillation
Pattern.
4. a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest according to claim 1, feature exist
In the detailed process of step S30 is：After the similar oscillation mode of l ranks for obtaining each double-fed fan motor machine according to step S20, i-th of phase
Like the frequency similarity q (f between oscillation modei) and damping similarity q (αi) be respectively：
Wherein, fAiIndicate the frequency with reference to double-fed fan motor machine i-th of oscillation mode of output power curve, fBiIndicate other double-fed wind
The frequency of i-th of oscillation mode of output power of motor curve, αAiIndicate fAiDecay factor, αBiIndicate fBiDecay factor, examine
Consider the similarity weight of frequency and damping, therefore the similarity between i-th of similar oscillation mode is：
Wherein a, b indicate frequency and the similarity weight of damping respectively；
The ratio of the energy of each merotype and gross energy is accounted for as the merotype energy in similarity calculation in oscillation mode
The weight of whole oscillation mode energy：
The similarity of all similar oscillation modes is after consideration energy weight：
Wherein EiIndicate the energy of double-fed fan motor machine i-th of oscillation mode of output power, the people having the same aspiration and interest is similar between calculating double-fed fan motor machine
After metrization index, the Coherent Generator Group in wind power plant is determined.
5. a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest according to claim 1, feature exist
In the detailed process of step S40 is：First, the stator of each double-fed fan motor machine, rotor voltage equation are indicated as follows：
Wherein, u, i, ψ indicate the voltage, electric current and magnetic linkage of double-fed fan motor machine respectively；R, L indicates the resistance of double-fed fan motor machine respectively
And inductance；S, r in subscript indicate the stator of double-fed fan motor machine, rotor respectively；D, q in subscript indicate d-axis, quadrature axis respectively
Component；LmIndicate the mutual inductance of double-fed fan motor machine；P is differential operator；S indicates revolutional slip；
In formula,It enables：
Wherein, Xss＞ ＞ Rs, so the stator pressure drop of double-fed fan motor machine can ignore that and disregard, define the equivalent transient state of double-fed fan motor machine
Built-in potential is as follows：
E'=Us-jXssIs (12)
Wherein, UsIndicate the set end voltage of double-fed fan motor machine, IsIndicate that the stator current of double-fed fan motor machine, j represent phase advanced 90
The symbol of degree；
Since double-fed fan motor machine is using decoupling control, externally inertia, influence of the double-fed fan motor machine to system are not embodied
It is characterized only by the power of its injection, the polymerization using output power as weight is used, if double-fed wind in Coherent Generator Group
The number of motor is N, waits the equivalent transient state built-in potential of check-insFor：
Wherein, piIndicate the active power of i-th of double-fed fan motor machine in Coherent Generator Group；The equivalent transient state built-in potential of equal check-ins is by same
Double-fed fan motor machine all in a group of planes is adjusted to determine that, according to the external characteristics of double-fed fan motor machine, weight takes shared by each double-fed fan motor machine
Certainly in their output power；
Wherein, the no-load voltage ratio of the multiple no-load voltage ratio phase-shifting transformer is：
Indicate the equivalent transient state built-in potential of i-th of double-fed fan motor machine in Coherent Generator Group, and the double-fed fan motor machine of same group is equivalent temporarily
State built-in potential busbar passes through the equivalent transient state built-in potential busbar for the check-ins such as multiple no-load voltage ratio phase-shifting transformer is connected to respectively.
6. a kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest according to claim 1, feature exist
In in step S50, the equivalent parameters of each equal check-ins is：
Wherein, N indicates that the number of double-fed fan motor machine in Coherent Generator Group, P indicate that the active power of double-fed fan motor machine, Q indicate double-fed
The reactive power of wind turbine, XsIndicate the stator reactance of double-fed fan motor machine, XrIndicate the rotor reactance of double-fed fan motor machine, HgIt indicates
The rotor inertia time constant of perceptual generator, H in double-fed fan motor machinetWhen indicating the rotor inertia of wind energy conversion system in double-fed fan motor machine
Between constant, KsIndicate the axis rigidity coefficient of wind energy conversion system in double-fed fan motor machine, the eq in subscript, which is respectively represented, waits check-ins parameter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810064146.0A CN108418242B (en) | 2018-01-23 | 2018-01-23 | Doubly-fed wind turbine dynamic equivalence method based on similarity coherence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810064146.0A CN108418242B (en) | 2018-01-23 | 2018-01-23 | Doubly-fed wind turbine dynamic equivalence method based on similarity coherence |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108418242A true CN108418242A (en) | 2018-08-17 |
| CN108418242B CN108418242B (en) | 2020-06-19 |
Family
ID=63126175
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810064146.0A Active CN108418242B (en) | 2018-01-23 | 2018-01-23 | Doubly-fed wind turbine dynamic equivalence method based on similarity coherence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108418242B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109657902A (en) * | 2018-10-30 | 2019-04-19 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Consider the double-fed blower Dynamic Equivalence of sub-synchronous oscillation characteristic |
| CN109787277A (en) * | 2018-12-28 | 2019-05-21 | 清华大学 | Power grid method is accessed based on the wind power plant direct current transportation that midfrequent AC bus collects |
| CN112928781A (en) * | 2021-03-23 | 2021-06-08 | 南方电网科学研究院有限责任公司 | Double-fed fan transient stability control method and system, computer equipment and medium |
| CN116402017A (en) * | 2023-06-02 | 2023-07-07 | 华北电力大学 | Simplified modeling method of doubly-fed fan for dynamic frequency analysis of power system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101882167A (en) * | 2009-11-25 | 2010-11-10 | 中国电力科学研究院 | Wind power station equivalent modeling method of large-scale wind power concentration access power grid |
| CN102624309A (en) * | 2012-03-29 | 2012-08-01 | 河海大学 | Variable-speed double-fed wind power station clustering equating method based on mechanical and electrical dynamic characteristics |
| CN103887815A (en) * | 2014-02-21 | 2014-06-25 | 华南理工大学 | Wind power plant parameter identification and dynamic equivalence method based on operation data |
-
2018
- 2018-01-23 CN CN201810064146.0A patent/CN108418242B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101882167A (en) * | 2009-11-25 | 2010-11-10 | 中国电力科学研究院 | Wind power station equivalent modeling method of large-scale wind power concentration access power grid |
| CN102624309A (en) * | 2012-03-29 | 2012-08-01 | 河海大学 | Variable-speed double-fed wind power station clustering equating method based on mechanical and electrical dynamic characteristics |
| CN103887815A (en) * | 2014-02-21 | 2014-06-25 | 华南理工大学 | Wind power plant parameter identification and dynamic equivalence method based on operation data |
Non-Patent Citations (3)
| Title |
|---|
| ZHENHUA XU等: "Aggregation of doubly-fed induction generator based wind farm considering storage device and coordinated control strategy", 《2011 INTERNATIONAL CONFERENCE ON ADVANCED POWER SYSTEM AUTOMATION AND PROTECTION》 * |
| 夏安俊，等: "双馈异步发电机风电场聚合模型研究", 《电网技术》 * |
| 朱林，等: "基于相似理论的电力系统动态等值精确度量化评价方法", 《电力科学与技术学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109657902A (en) * | 2018-10-30 | 2019-04-19 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Consider the double-fed blower Dynamic Equivalence of sub-synchronous oscillation characteristic |
| CN109657902B (en) * | 2018-10-30 | 2020-11-13 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Doubly-fed wind turbine dynamic equivalence method considering subsynchronous oscillation characteristics |
| CN109787277A (en) * | 2018-12-28 | 2019-05-21 | 清华大学 | Power grid method is accessed based on the wind power plant direct current transportation that midfrequent AC bus collects |
| CN109787277B (en) * | 2018-12-28 | 2021-07-16 | 清华大学 | Method for accessing wind power plant direct-current transmission into power grid based on medium-frequency alternating-current bus collection |
| CN112928781A (en) * | 2021-03-23 | 2021-06-08 | 南方电网科学研究院有限责任公司 | Double-fed fan transient stability control method and system, computer equipment and medium |
| CN116402017A (en) * | 2023-06-02 | 2023-07-07 | 华北电力大学 | Simplified modeling method of doubly-fed fan for dynamic frequency analysis of power system |
| CN116402017B (en) * | 2023-06-02 | 2023-08-22 | 华北电力大学 | Simplified modeling method of doubly-fed fan for dynamic frequency analysis of power system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108418242B (en) | 2020-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103887815B (en) | Based on wind energy turbine set parameter identification and the Dynamic Equivalence of service data | |
| CN101882167B (en) | Wind power station equivalent modeling method of large-scale wind power concentration access power grid | |
| CN108418242A (en) | A kind of double-fed fan motor machine Dynamic Equivalence based on the similarity people having the same aspiration and interest | |
| CN107086606A (en) | A kind of equivalent asynchronous motor load model parameters discrimination method of power distribution network synthesis | |
| CN112818491A (en) | Wind power plant aggregation equivalent modeling method based on principal component analysis and clustering algorithm | |
| CN103219725A (en) | Wind power plant equivalent modeling method based on real-time operation data | |
| Wang et al. | Dynamic equivalent modeling for wind farms with DFIGs using the artificial bee colony with K-means algorithm | |
| CN105678033B (en) | A kind of wind-powered electricity generation station equivalent modeling method suitable for electromagnetic transient simulation | |
| CN102801158A (en) | Method for calculating time-lag electric power system eigenvalue and discriminating stability based on Pade approximation | |
| CN109657902B (en) | Doubly-fed wind turbine dynamic equivalence method considering subsynchronous oscillation characteristics | |
| CN103793585B (en) | A kind of method for improving theoretical line loss caluclation integrality and accuracy | |
| CN103761349A (en) | Wind power plant equivalent modeling method based on probability coherence of wind power units | |
| Biserica et al. | Network partition for coordinated control in active distribution networks | |
| CN108649562A (en) | A kind of power system load modeling parameter identification method based on blue wolf algorithm | |
| CN109921426A (en) | Wind-electricity integration system probability load flow calculation method based on CV-KDE | |
| CN104617578B (en) | Method for acquiring available power transmission capability of power system with wind power plant | |
| CN108258725A (en) | A kind of double-fed fan motor machine Dynamic Equivalence based on the equivalent generator rotor angle people having the same aspiration and interest | |
| CN116738636A (en) | Multi-machine equivalent method of doubly-fed wind power station considering impedance characteristics and synchronization mechanism | |
| CN111130135A (en) | Power system inertia calculation method suitable for high-proportion new energy access | |
| CN107069808A (en) | A kind of new energy power station model error localization method and device | |
| CN108460228B (en) | Wind power plant equivalence method based on multi-objective optimization algorithm | |
| CN106410862A (en) | Wind power plant single machine equivalent method based on active recovery slope correction | |
| CN111737919B (en) | Direct-drive wind farm grouping method suitable for subsynchronous oscillation analysis | |
| CN109088436B (en) | VSC multi-scale transient modeling method | |
| CN112736927A (en) | Wind power plant uncertainty trend modeling method based on probability box theory |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |