CN106168994B - A kind of subsynchronous frequency equivalent modeling method of wind power plant - Google Patents

A kind of subsynchronous frequency equivalent modeling method of wind power plant Download PDF

Info

Publication number
CN106168994B
CN106168994B CN201610515963.4A CN201610515963A CN106168994B CN 106168994 B CN106168994 B CN 106168994B CN 201610515963 A CN201610515963 A CN 201610515963A CN 106168994 B CN106168994 B CN 106168994B
Authority
CN
China
Prior art keywords
equivalent
wind
generator
power plant
wind power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610515963.4A
Other languages
Chinese (zh)
Other versions
CN106168994A (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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201610515963.4A priority Critical patent/CN106168994B/en
Publication of CN106168994A publication Critical patent/CN106168994A/en
Application granted granted Critical
Publication of CN106168994B publication Critical patent/CN106168994B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Eletrric Generators (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention belongs to power system securities to stablize field more particularly to a kind of subsynchronous frequency equivalent modeling method of wind power plant.It is characterized in that, seeking the equivalent impedance of single wind generator the method includes establishing the equivalent circuit of single wind generator, wind power plant cluster electric network is established according to wind farm network topological structure;Seek the subsynchronous frequency characteristic of wind power plant cluster electric network;K-means clustering algorithm is carried out to the wind-driven generator in wind power plant and divides group, the subsynchronous frequency for then carrying out wind power plant primary network is equivalent;Carry out the subsynchronous frequency Equivalent Model amendment of wind power plant.The present invention can be used in the wind power plant equivalent modeling of sub-synchronous oscillation analysis, entire wind power plant is characterized with one or more wind-driven generator, and there is similar subsynchronous frequency characteristic with wind power plant, it is the basis analyzed wind power plant sub-synchronous oscillation mechanism and further take effective braking measure.

Description

A kind of subsynchronous frequency equivalent modeling method of wind power plant
Technical field
The invention belongs to power system securities to stablize field more particularly to a kind of subsynchronous frequency equivalent modeling side of wind power plant Method.
Background technique
THE WIND ENERGY RESOURCES IN CHINA is mainly contained in western part, and a large amount of electric loads are located at eastern region, this makes " Xi electricity Dong Send " principles and policies of the strategy as China's electric power development.Since series capacitor compensation technology can significantly improve ac transmission Power delivery capabilities, so that it is worldwide widely used.In the big of the wind power bases such as China Inner Mongol, Xinjiang, Gansu For type wind power plant since local load is lower, this wind power plant, which mends ac transmission through string, will become a kind of universal transmission of electricity of China's power grid Mode.Under certain condition since Wind turbines and series compensation interact, wind power plant sub-synchronous oscillation problem can be caused, Seriously threaten the safe and stable operation of electric system.
Wind power plant is the Wind turbines cluster for converting wind energy into electric energy, and string mends ac transmission can cause under certain condition Sub-synchronous oscillation seriously threatens the safety of wind power plant.Sub-synchronous oscillation repeatedly occurs both at home and abroad causes wind farm device to damage thing Therefore cause serious economic loss.Wind power plant equivalent modeling is carried out, simulation analysis sub-synchronous oscillation phenomenon is analysis wind power plant Sub-synchronous oscillation mechanism and the basis for further taking effective braking measure have great economic and social benefit.
Currently used wind-powered electricity generation station equivalent modeling method is for transient stability analysis of power system or Low Frequency Oscillation Analysis , it is not suitable for sub-synchronous oscillation analysis.Up to the present, there are no a set of wind power plant for sub-synchronous oscillation analysis is equivalent Modeling method.Entire wind power plant often is simulated with a Fans in research, whether entire wind power plant can be characterized for the blower Subsynchronous frequency characteristic do not have still studies have reported that.In consideration of it, the invention patent proposes a kind of subsynchronous frequency of Large Scale Wind Farm Integration Equivalent modeling method, to lay the foundation for wind power plant sub-synchronous oscillation simulation analysis.
Summary of the invention
To solve the above-mentioned problems, the technical scheme adopted by the invention is as follows it is as follows:
A kind of subsynchronous frequency equivalent modeling method of wind power plant, which is characterized in that the described method comprises the following steps:
Step 1, the equivalent circuit for establishing single wind generator, seek the equivalent impedance of single wind generator, according to wind Electric field network topological structure establishes wind power plant cluster electric network;If the wind power plant cluster electric network of trunk line type distribution then will be each The equivalent impedance of wind-driven generator is sequentially ingressed into system in trunk line type;If the wind power plant cluster electric network of emanant distribution is then By the equivalent impedance of each wind-driven generator in emanant access system successively in parallel;
Step 2, the subsynchronous frequency characteristic for seeking wind power plant cluster electric network;
Step 3 divides group to the wind-driven generator progress k-means clustering algorithm in wind power plant, then carries out wind power plant one The subsynchronous frequency of secondary network is equivalent;
Step 4 carries out the subsynchronous frequency Equivalent Model amendment of wind power plant.
The step 1 establishes the equivalent circuit of single wind generator, seeks the tool of the equivalent impedance of single wind generator Body process is
When wind-driven generator is through serial supplementary line access system, it is f that frequency is flowed through in generator unit statorerSubsynchronous electric current When, it is 2 π f that angular frequency is generated in air gaperRotating mmf, wind-driven generator is entered in terms of stator terminal, wherein ωrFor wind The rotor frequency of power generator, natural resonance frequency
Series compensation degrees k and the relationship of serial supplementary line parameter are
Natural resonance frequency ωerAlso it can be acquired by following formula
The equivalent impedance of single wind generator
Wherein, revolutional slipXrFor wind power generator rotor reactance, XmFor mutual induction reactance, XsFor wind-driven generator Stator reactance, XTFor transformer reactance, ReqFor the equivalent resistance of single wind generator, xeqFor the equivalent of single wind generator Reactance;RrFor the rotor resistance of wind-driven generator, LrFor the rotor leakage inductance of wind-driven generator, RsFor the stator electricity of wind-driven generator Resistance, LsFor the stator leakage inductance of wind-driven generator, L is the inductance of serial supplementary line, C is the capacitor of serial supplementary line, R is serial supplementary line Resistance, LTFor transformer equivalent inductance, j is empty unit.
The detailed process for seeking the subsynchronous frequency characteristic of wind power plant cluster electric network in the step 2 is
Using series connection frequency sweep method, the system of research is simulated with positive sequence net, non-generator to be studied is used in network Subtranient reactance is equivalent, and generator to be studied using its asynchronous generator model it is equivalent, each single wind generator etc. A unit current source is placed in the generator amature side imitated in circuit, then changes the unitary current source frequency, measures the unit Current source both end voltage, the voltage value are the equivalent impedance looked over from generator amature to be studied to system side, according to etc. Imitate impedance ZeqThe real imaginary part of expression formula draws equivalent resistance, equivalent reactance curve varying with frequency;
Equivalent impedance expression is in the wind power plant cluster electric network of trunk line type distribution
Equivalent impedance expression is in the wind power plant cluster electric network of emanant distribution
Zeq=Zeq_k+(Xd1//Xd2//…//Xdk-1//Xdk+1//…//Xdn//Xs)=Req+jXeq
Wherein X "diFor the subtranient reactance of the i-th typhoon power generator, XsIt " is the reactance of Infinite bus system, Zeq_kFor wait grind The equivalent impedance for the single wind generator studied carefully, n are the number of units of wind-driven generator in wind power plant.
Carrying out k-means clustering algorithm to the wind-driven generator in wind power plant in the step 3 divides the detailed process of group to be
Using the subsynchronous frequency characteristic situation of change of wind-driven generator as group's judgment criteria is divided, group's index d will be dividedijIt is identical Wind-driven generator is divided into a group;When there are n typhoon power generators in wind power plant, n* (n-1)/2 distance can be obtained;
dijTo divide group's index, the subsynchronous frequency characteristic of wind-driven generator is characterized;
Wherein, Chebyshev's distance d1ij=max | Xi(ω)-Xj(ω) |, observation be single wind generator time frequency Rate characteristic;Euclidean distanceWhat is observed is the similarity degree between each typhoon power generator;Xi (ω) is subsynchronous frequency reactance value of the i-th typhoon power generator in series compensation degrees k and when generator speed is ω;Xj(ω) is the Subsynchronous frequency reactance value of the j typhoon power generator in series compensation degrees k and when generator speed is ω;Q is Chebyshev's distance d1ij Weighted value, take 50%.
The detailed process of the subsynchronous frequency equivalence of progress wind power plant primary network is in the step 3
Step 301, wind-powered electricity generation field parameters are equivalent;
1. generator parameter
xsIt is the reactance of generator unit stator, xrIt is the reactance of generator amature, xs_eqIt is the equivalent reactance of generator unit stator, xr_eqIt is the equivalent reactance of generator amature respectively;xmIt is the excitation reactance of generator, xm_eqIt is the equivalent excitation electricity of generator It is anti-;rsIt is the resistance of generator amature, rrIt is the resistance of generator amature, rs_eqIt is the equivalent resistance of generator unit stator, rr_eqIt is The equivalent resistance of generator amature;M is the number of units of this group of wind-driven generators;
2. transformer parameter
STFor transformer capacity, ST_eqFor Equivalent Transformer capacity;ZTFor transformer impedance, ZT_eqFor Equivalent Transformer resistance It is anti-;
Step 302, cluster electric network parameter predigesting
Each wind-driven generator in identical a group before the voltage at the grid-connected place of wind-driven generator after equivalence is equivalent to equivalence The average value of the weighting of grid-connected place's voltage, the weight of each wind-driven generator are their own output powers;
A. the wind power plant cluster electric network of trunk line type distribution
Points of common connection PCC with for kth typhoon power generator impedance on voltage UkBetween voltage difference be U is grid entry point voltage, ZiIt is the impedance of i-th unit connecting cable, PjIt is that jth typhoon power generator is defeated Power out;
It will 1 to the m wind-driven generator in the equivalent preceding same group be equivalent obtains if equivalent wind-driven generator 1 is , the average voltage before merging after the weighting of m typhoon power generator is poorPiIt is that the i-th typhoon power generator is defeated Power out, Zeq_1For the equivalent impedance of wind-driven generator 1 after equivalent merging, wind-driven generator 1 and public coupling after equivalent merging Voltage difference between chalazaΔUiFor the voltage drop in the impedance of i-th equivalent wind-driven generator;
Equivalent front and backObtain the equivalent impedance of serial supplementary line Wherein LeqFor the equivalent inductance of serial supplementary line, CeqFor the equivalent capacity of serial supplementary line, ReqLine is mended for string The equivalent resistance on road;
B. the wind power plant cluster electric network of emanant distribution
If the power factor of wind-driven generator is approximately 1, UkThe difference of voltage between bus PCC isU It is the voltage value of wind-driven generator grid entry point, ZiIt is the impedance value of the cable of the i-th typhoon power generator connection, ZkIt is kth platform etc. Imitate the impedance of wind-driven generator, PkIt is the active power of the equivalent wind-driven generator of kth platform;
If equivalent wind-driven generator 1 for will it is equivalent before the same group 1 Dao m wind-driven generator equivalence acquisition, The average voltage of the weighting of this m typhoon power generator is poor before equivalentZeq_1It is sent out for wind-force after equivalent merging The equivalent impedance of motor 1, voltage difference between equivalent rear wind-driven generator 1 and common coupling node
Equivalent front and backObtain the equivalent impedance of serial supplementary line
The subsynchronous modified detailed process of frequency Equivalent Model of wind power plant is carried out in the step 4 is
Step 401 obtains subsynchronous frequency characteristic after wind power plant equivalence
Single wind generator equivalent impedance Z ' after equivalenceeqExpression formula are as follows:
Wherein Rr_eqFor the equivalent rotor resistance of wind-driven generator, Xr_eqFor the equivalent rotor reactance of wind-driven generator, Xm_eq It is equivalent mutual induction reactance, Rs_eqFor the equivalent stator resistance of wind-driven generator, Xs_eqFor the equivalent stator reactance of wind-driven generator, XT_eqFor transformer equivalent reactance, LeqFor the equivalent inductance of serial supplementary line, CeqFor the equivalent capacity of serial supplementary line, ReqIt is mended for string The equivalent resistance of route;
Z′eq=R 'eq+jX′eq
Wherein R 'eqFor wind power plant equivalent resistance, the X ' after equivalenceeqFor the wind power plant equivalent reactance after equivalence;
The error analysis of the equivalent front and back of step 402 wind power plant
The error of the equivalent front and back of wind power plant
Beneficial effect
The present invention can be used in the wind power plant equivalent modeling of sub-synchronous oscillation analysis, with one or more wind-driven generator table Levy entire wind power plant, and there is similar subsynchronous frequency characteristic with wind power plant, be analysis wind power plant sub-synchronous oscillation mechanism and Further take the basis of effective braking measure.
Detailed description of the invention
The subsynchronous frequency equivalent modeling method flow chart of Fig. 1 wind power plant;
Fig. 2 wind power plant trunk line type distribution schematic diagram;
The emanant distribution schematic diagram of Fig. 3 wind power plant;
Fig. 4 wind-driven generator string mends circuit diagram;
Fig. 5 wind-driven generator equivalent impedance schematic diagram;
Fig. 6 trunk line type is distributed equivalent impedance schematic diagram;
The emanant distribution equivalent impedance schematic diagram of Fig. 7;
The wind power plant equivalent model schematic diagram of tri- machine of Fig. 8 characterization;
The equivalent subsynchronous frequency impedance figure in front and back of Fig. 9 wind power plant.
Specific embodiment
With reference to the accompanying drawing, it elaborates to the present invention.Using certain three typhoon power generator as example, if this three typhoons power Generator in emanant access system, wherein No. 1 and No. 2 be same model wind-driven generator, generator, transformer and The relevant parameter of route such as following table, generator amature reactance therein, stator reactance, transformer reactance are to join under synchronizing frequency Number, when calculating, need to convert.
Form wind power plant cluster electric network, in the wind power plant cluster electric network of emanant distribution each wind-driven generator etc. Imitate impedance Zeq1、Zeq2、Zeq3In in emanant access system successively in parallel.
Using series connection frequency sweep method, the system of research is simulated with positive sequence net, non-generator to be studied is used in network Subtranient reactance is equivalent, and generator to be studied is equivalent using its asynchronous generator model, the power generation in above-mentioned equivalent circuit A unit current source is placed in machine rotor side, then changes its frequency, measures its both end voltage, as turns from generator to be studied The equivalent impedance that son is looked over to system side draws equivalent resistance, equivalent reactance curve varying with frequency according to imaginary part in fact.
The subsynchronous frequency equivalent modeling method flow chart of Fig. 1 wind power plant.
Step 1, the equivalent circuit for establishing single wind generator, seek the equivalent impedance of single wind generator, according to wind Electric field network topological structure establishes wind power plant cluster electric network;If the wind power plant cluster electric network of trunk line type distribution then will be each The equivalent impedance of wind-driven generator is sequentially ingressed into system in trunk line type;Fig. 2 is wind power plant trunk line type distribution schematic diagram;If putting The wind power plant cluster electric network of formula distribution is penetrated then by the equivalent impedance of each wind-driven generator in emanant access system successively in parallel In system;Fig. 3 is the emanant distribution schematic diagram of wind power plant;
Step 2, the subsynchronous frequency characteristic for seeking wind power plant cluster electric network;
Step 3 divides group to the wind-driven generator progress k-means clustering algorithm in wind power plant, then carries out wind power plant one The subsynchronous frequency of secondary network is equivalent;
Step 4 carries out the subsynchronous frequency Equivalent Model amendment of wind power plant.
Fig. 4 is that wind-driven generator string mends circuit diagram;The step 1 establishes the equivalent circuit of single wind generator, The detailed process for seeking the equivalent impedance of single wind generator is
When wind-driven generator is through serial supplementary line access system, it is f that frequency is flowed through in generator unit statorerSubsynchronous electric current When, it is 2 π f that angular frequency is generated in air gaperRotating mmf, wind-driven generator is entered in terms of stator terminal, wherein ωrFor wind The rotor frequency of power generator, natural resonance frequency
Series compensation degrees k and the relationship of serial supplementary line parameter are
Natural resonance frequency ωerAlso it can be acquired by following formula
The equivalent impedance of single wind generator
Wherein, revolutional slipXrFor wind power generator rotor reactance, XmFor mutual induction reactance, XsFor wind-driven generator Stator reactance, XTFor transformer reactance, ReqFor the equivalent resistance of single wind generator, XeqFor the equivalent of single wind generator Reactance;RrFor the rotor resistance of wind-driven generator, LrFor the rotor leakage inductance of wind-driven generator, RsFor the stator electricity of wind-driven generator Resistance, LsFor the stator leakage inductance of wind-driven generator, L is the inductance of serial supplementary line, C is the capacitor of serial supplementary line, R is serial supplementary line Resistance, LTFor transformer equivalent inductance, j is empty unit.Fig. 5 is wind-driven generator equivalent impedance schematic diagram;
The detailed process for seeking the subsynchronous frequency characteristic of wind power plant cluster electric network in the step 2 is
Using series connection frequency sweep method, the system of research is simulated with positive sequence net, non-generator to be studied is used in network Subtranient reactance is equivalent, and generator to be studied using its asynchronous generator model it is equivalent, each single wind generator etc. A unit current source is placed in the generator amature side imitated in circuit, then changes the unitary current source frequency, measures the unit Current source both end voltage, the voltage value are the equivalent impedance looked over from generator amature to be studied to system side, according to etc. Imitate impedance ZeqThe real imaginary part of expression formula draws equivalent resistance, equivalent reactance curve varying with frequency;
Equivalent impedance expression is in the wind power plant cluster electric network of trunk line type distribution
Fig. 6 is that trunk line type is distributed equivalent impedance schematic diagram;
Equivalent impedance expression is in the wind power plant cluster electric network of emanant distribution
Zeq=Zeq_k+(Xd1//Xd2//…//Xdk-1//Xdk+1//…//Xdn//Xs)=Req+jXeq
Wherein X "diFor the subtranient reactance of the i-th typhoon power generator, XsIt " is the reactance of Infinite bus system, Zeq_kFor wait grind The equivalent impedance for the single wind generator studied carefully, n are the number of units of wind-driven generator in wind power plant, and Fig. 7 is that emanant distribution is equivalent Impedance schematic diagram.
Carrying out k-means clustering algorithm to the wind-driven generator in wind power plant in the step 3 divides the detailed process of group to be
Using the subsynchronous frequency characteristic situation of change of wind-driven generator as group's judgment criteria is divided, group's index d will be dividedijIt is identical Wind-driven generator is divided into a group;When there are n typhoon power generators in wind power plant, n* (n-1)/2 distance can be obtained;
dijTo divide group's index, the subsynchronous frequency characteristic of wind-driven generator is characterized;
Wherein, Chebyshev's distance d1ij=max | Xi(ω)-Xj(ω) |, observation be single wind generator time frequency Rate characteristic;Euclidean distanceWhat is observed is the similarity degree between each typhoon power generator;Xi (ω) is subsynchronous frequency reactance value of the i-th typhoon power generator in series compensation degrees k and when generator speed is ω;Xj(ω) is the Subsynchronous frequency reactance value of the j typhoon power generator in series compensation degrees k and when generator speed is ω;Q is Chebyshev's distance d1ij Weighted value, take 50%.
The detailed process of the subsynchronous frequency equivalence of progress wind power plant primary network is in the step 3
Step 301, wind-powered electricity generation field parameters are equivalent;
1. generator parameter
xsIt is the reactance of generator unit stator, xrIt is the reactance of generator amature, xs_eqIt is the equivalent reactance of generator unit stator, xr_eqIt is the equivalent reactance of generator amature respectively;xmIt is the excitation reactance of generator, xm_eqIt is the equivalent excitation electricity of generator It is anti-;rsIt is the resistance of generator amature, rrIt is the resistance of generator amature, rs_eqIt is the equivalent resistance of generator unit stator, rr_eqIt is The equivalent resistance of generator amature;M is the number of units of this group of wind-driven generators;
2. transformer parameter
STFor transformer capacity, ST_eqFor Equivalent Transformer capacity;ZTFor transformer impedance, ZT_eqFor Equivalent Transformer resistance It is anti-;
Step 302, cluster electric network parameter predigesting
Each wind-driven generator in identical a group before the voltage at the grid-connected place of wind-driven generator after equivalence is equivalent to equivalence The average value of the weighting of grid-connected place's voltage, the weight of each wind-driven generator are their own output powers;
A. the wind power plant cluster electric network of trunk line type distribution
Points of common connection PCC with for kth typhoon power generator impedance on voltage UkBetween voltage difference be U is grid entry point voltage, ZiIt is the impedance of i-th unit connecting cable, PjIt is that jth typhoon power generator is defeated Power out;
It will 1 to the m wind-driven generator in the equivalent preceding same group be equivalent obtains if equivalent wind-driven generator 1 is , the average voltage before merging after the weighting of m typhoon power generator is poorPiIt is that the i-th typhoon power generator is defeated Power out, Zeq_1For the equivalent impedance of wind-driven generator 1 after equivalent merging, wind-driven generator 1 and public coupling after equivalent merging Voltage difference between chalazaΔUiFor the voltage drop in the impedance of i-th equivalent wind-driven generator;
Equivalent front and backObtain the equivalent impedance of serial supplementary line Wherein LeqFor the equivalent inductance of serial supplementary line, CeqFor the equivalent capacity of serial supplementary line, ReqLine is mended for string The equivalent resistance on road;
B. the wind power plant cluster electric network of emanant distribution
If the power factor of wind-driven generator is approximately 1, UkThe difference of voltage between bus PCC isU It is the voltage value of wind-driven generator grid entry point, ZiIt is the impedance value of the cable of the i-th typhoon power generator connection,kIt is kth platform etc. Imitate the impedance of wind-driven generator, PkIt is the active power of the equivalent wind-driven generator of kth platform;
If equivalent wind-driven generator 1 for will it is equivalent before the same group 1 Dao m wind-driven generator equivalence acquisition, The average voltage of the weighting of this m typhoon power generator is poor before equivalentZeq_1It is sent out for wind-force after equivalent merging The equivalent impedance of motor 1, voltage difference between equivalent rear wind-driven generator 1 and common coupling node
Equivalent front and backObtain the equivalent impedance of serial supplementary line
The subsynchronous modified detailed process of frequency Equivalent Model of wind power plant is carried out in the step 4 is
Step 401 obtains subsynchronous frequency characteristic after wind power plant equivalence
Single wind generator equivalent impedance Z ' after equivalenceeqExpression formula are as follows:
Wherein Rr_eqFor the equivalent rotor resistance of wind-driven generator, Xr_eqFor the equivalent rotor reactance of wind-driven generator, Xm_eq It is equivalent mutual induction reactance, Rs_eqFor the equivalent stator resistance of wind-driven generator, Xs_eqFor the equivalent stator reactance of wind-driven generator, XT_eqFor transformer equivalent reactance, LeqFor the equivalent inductance of serial supplementary line, CeqFor the equivalent capacity of serial supplementary line, ReqIt is mended for string The equivalent resistance of route;
Z′eq=R 'eq+jX′eq
Wherein R 'eqFor wind power plant equivalent resistance, the X ' after equivalenceeqFor the wind power plant equivalent reactance after equivalence;
The error analysis of the equivalent front and back of step 402 wind power plant
The error of the equivalent front and back of wind power plant
K-means clustering algorithm is carried out to wind-driven generator before equivalent and divides group, by the subsynchronous frequency characteristic of wind-driven generator Situation of change as group's judgment criteria is divided will there is the wind-driven generator of similar subsynchronous frequency characteristic to be divided into a group.1,2, No. 3 Wind-driven generator is divided into same a group.Fig. 8 is the wind power plant equivalent model schematic diagram of three machines characterization.
Wind-powered electricity generation field parameters equivalence and cluster electric network parameter predigesting are carried out, subsynchronous frequency characteristic after wind power plant equivalence is obtained Curve.The equivalent front and back contrast schematic diagram of wind power plant is shown in Fig. 9.

Claims (5)

1. a kind of subsynchronous frequency equivalent modeling method of wind power plant, which is characterized in that the described method comprises the following steps:
Step 1, the equivalent circuit for establishing single wind generator, seek the equivalent impedance of single wind generator, according to wind power plant Network topology structure establishes wind power plant cluster electric network;If the wind power plant cluster electric network of trunk line type distribution is then by each wind-force The equivalent impedance of generator is sequentially ingressed into system in trunk line type;If the wind power plant cluster electric network of emanant distribution then will be each The equivalent impedance of a wind-driven generator is in emanant access system successively in parallel;
Step 2, the subsynchronous frequency characteristic for seeking wind power plant cluster electric network;
Step 3 divides group to the wind-driven generator progress k-means clustering algorithm in wind power plant, then carries out the primary net of wind power plant The subsynchronous frequency of network is equivalent;
Step 4 carries out the subsynchronous frequency Equivalent Model amendment of wind power plant;
The step 1 establishes the equivalent circuit of single wind generator, seeks the specific mistake of the equivalent impedance of single wind generator Cheng Wei
When wind-driven generator is through serial supplementary line access system, it is f that frequency is flowed through in generator unit statorerSubsynchronous electric current when, It is 2 π f that it generates angular frequency in air gaperRotating mmf, wind-driven generator is entered in terms of stator terminal, wherein ωrFor wind-force hair The rotor frequency of motor, natural resonance frequency
Series compensation degrees k and the relationship of serial supplementary line parameter are
Natural resonance frequency ωerAlso it can be acquired by following formula
The equivalent impedance of single wind generator
Wherein, revolutional slipXrFor wind power generator rotor reactance, XmFor mutual induction reactance, XsFor aerogenerator stator Reactance, XTFor transformer reactance, ReqFor the equivalent resistance of single wind generator, XeqFor the equivalent electricity of single wind generator It is anti-;RrFor the rotor resistance of wind-driven generator, LrFor the rotor leakage inductance of wind-driven generator, RsFor the stator resistance of wind-driven generator, LsFor the stator leakage inductance of wind-driven generator, L is the inductance of serial supplementary line, C is the capacitor of serial supplementary line, the electricity that R is serial supplementary line Resistance, LTFor transformer equivalent inductance, j is empty unit.
2. the subsynchronous frequency equivalent modeling method of wind power plant according to claim 1, which is characterized in that in the step 2 The detailed process for seeking the subsynchronous frequency characteristic of wind power plant cluster electric network is
Using series connection frequency sweep method, the system of research is simulated with positive sequence net, non-generator to be studied is used secondary temporary in network State reactance is equivalent, and generator to be studied is equivalent using its asynchronous generator model, in the equivalent electricity of each single wind generator A unit current source is placed in generator amature side in road, then changes the unitary current source frequency, measures the unitary current Source both end voltage, which is the equivalent impedance looked over from generator amature to be studied to system side, according to equivalent resistance Anti- ZeqThe real imaginary part of expression formula draws equivalent resistance, equivalent reactance curve varying with frequency;
Equivalent impedance expression is in the wind power plant cluster electric network of trunk line type distribution
Equivalent impedance expression is in the wind power plant cluster electric network of emanant distribution
Zeq=Zeq_k+(X″d1//X″d2//...//X″dk-1//X″dk+1//...//X″dn//X″s)=Req+jXeq
Wherein X "diFor the subtranient reactance of the i-th typhoon power generator, X "sFor the reactance of Infinite bus system, Zeq_kIt is to be studied The equivalent impedance of single wind generator, n are the number of units of wind-driven generator in wind power plant.
3. the subsynchronous frequency equivalent modeling method of wind power plant according to claim 1, which is characterized in that in the step 3 Carrying out k-means clustering algorithm to the wind-driven generator in wind power plant divides the detailed process of group to be
Using the subsynchronous frequency characteristic situation of change of wind-driven generator as group's judgment criteria is divided, group's index d will be dividedijIdentical wind-force Generator is divided into a group;When there are n typhoon power generators in wind power plant, n × (n-1)/2 distance can be obtained;
dijTo divide group's index, the subsynchronous frequency characteristic of wind-driven generator is characterized;
Wherein, Chebyshev's distance d1ij=max | Xi(ω)-Xj(ω) |, observation be single wind generator secondary frequencies it is special Property;Euclidean distanceWhat is observed is the similarity degree between each typhoon power generator;Xi(ω) Subsynchronous frequency reactance value for the i-th typhoon power generator in series compensation degrees k and when generator speed is ω;Xj(ω) is jth platform Subsynchronous frequency reactance value of the wind-driven generator in series compensation degrees k and when generator speed is ω;Q is Chebyshev's distance d1ij's Weighted value takes 50%.
4. the subsynchronous frequency equivalent modeling method of wind power plant according to claim 1, which is characterized in that in the step 3 The detailed process of subsynchronous frequency equivalence for carrying out wind power plant primary network is
Step 301, wind-powered electricity generation field parameters are equivalent;
Generator parameter includes:
xsIt is the reactance of generator unit stator, xrIt is the reactance of generator amature, xs_eqIt is the equivalent reactance of generator unit stator, xr_eqPoint It is not the equivalent reactance of generator amature;xmIt is the excitation reactance of generator, xm_eqIt is the equivalent excitation reactance of generator;rsIt is The resistance of generator amature, rrIt is the resistance of generator amature, rs_eqIt is the equivalent resistance of generator unit stator, rr_eqIt is generator The equivalent resistance of rotor;M is the number of units of this group of wind-driven generators;
Transformer parameter includes:
STFor transformer capacity, ST_eqFor Equivalent Transformer capacity;ZTFor transformer impedance, ZT_eqFor Equivalent Transformer impedance;
Step 302, cluster electric network parameter predigesting
Each wind-driven generator is grid-connected in identical a group before the voltage at the grid-connected place of wind-driven generator after equivalence is equivalent to equivalence Locate the average value of the weighting of voltage, the weight of each wind-driven generator is their own output power;
The wind power plant cluster electric network of situation a. trunk line type distribution
Points of common connection PCC with for kth typhoon power generator impedance on voltage UkBetween voltage difference beU is grid entry point voltage, ZiIt is the impedance of i-th unit connecting cable, PjIt is jth typhoon power hair The power of motor output;
If equivalent wind-driven generator 1 be will 1 to the m wind-driven generator in the equivalent preceding same group is equivalent obtains, conjunction And the average voltage after the weighting of preceding m typhoon power generator is poorPiIt is that the i-th typhoon power generator exports Power, Zeq_1For the equivalent impedance of wind-driven generator 1 after equivalent merging, wind-driven generator 1 and point of common coupling after equivalent merging Between voltage differenceΔUiFor the voltage drop in the impedance of i-th equivalent wind-driven generator;
Equivalent front and backObtain the equivalent impedance of serial supplementary line Wherein LeqFor the equivalent inductance of serial supplementary line, CeqFor the equivalent capacity of serial supplementary line, ReqLine is mended for string The equivalent resistance on road;
The wind power plant cluster electric network of the emanant distribution of situation b.
If the power factor of wind-driven generator is approximately 1, UkThe difference of voltage between bus PCC isU is wind The voltage value of power generator grid entry point, ZiIt is the impedance value of the cable of the i-th typhoon power generator connection, ZkIt is the equivalent wind of kth platform The impedance of power generator, PkIt is the active power of the equivalent wind-driven generator of kth platform;
If equivalent wind-driven generator 1 for will it is equivalent before the same group 1 Dao m wind-driven generator equivalence acquisition, it is equivalent The average voltage of the weighting of this m typhoon power generator is poor beforeZeq_1For wind-driven generator after equivalent merging 1 equivalent impedance, voltage difference between equivalent rear wind-driven generator 1 and common coupling node
Equivalent front and backObtain the equivalent impedance of serial supplementary line
5. the subsynchronous frequency equivalent modeling method of wind power plant according to claim 1, which is characterized in that in the step 4 Carrying out the subsynchronous modified detailed process of frequency Equivalent Model of wind power plant is
Step 401 obtains subsynchronous frequency characteristic after wind power plant equivalence
Single wind generator equivalent impedance Z ' after equivalenceeqExpression formula are as follows:
Wherein Rr_eqFor the equivalent rotor resistance of wind-driven generator, Xr_eqFor the equivalent rotor reactance of wind-driven generator, Xm_eqIt is Imitate mutual induction reactance, Rs_eqFor the equivalent stator resistance of wind-driven generator, Xs_eqFor the equivalent stator reactance of wind-driven generator, XT_eqFor Transformer equivalent reactance, LeqFor the equivalent inductance of serial supplementary line, CeqFor the equivalent capacity of serial supplementary line, ReqFor serial supplementary line Equivalent resistance;
Z′eq=R 'eq+jX′eq
Wherein R 'eqFor wind power plant equivalent resistance, the X ' after equivalenceeqFor the wind power plant equivalent reactance after equivalence;
The error analysis of the equivalent front and back of step 402 wind power plant
The error of the equivalent front and back of wind power plant
CN201610515963.4A 2016-07-01 2016-07-01 A kind of subsynchronous frequency equivalent modeling method of wind power plant Active CN106168994B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610515963.4A CN106168994B (en) 2016-07-01 2016-07-01 A kind of subsynchronous frequency equivalent modeling method of wind power plant

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610515963.4A CN106168994B (en) 2016-07-01 2016-07-01 A kind of subsynchronous frequency equivalent modeling method of wind power plant

Publications (2)

Publication Number Publication Date
CN106168994A CN106168994A (en) 2016-11-30
CN106168994B true CN106168994B (en) 2019-03-01

Family

ID=58064741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610515963.4A Active CN106168994B (en) 2016-07-01 2016-07-01 A kind of subsynchronous frequency equivalent modeling method of wind power plant

Country Status (1)

Country Link
CN (1) CN106168994B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106991229B (en) * 2017-03-30 2020-08-21 华北电力大学 Wind power plant equivalent modeling method for complex topology
CN108134399B (en) * 2018-01-12 2020-10-13 清华大学 Method and device for optimizing full working condition of network side subsynchronous damping controller
CN108154315B (en) * 2018-01-19 2021-11-05 国电南瑞科技股份有限公司 Grid-connected converter subsynchronous oscillation risk analysis method considering influence of phase-locked loop
CN109787281A (en) * 2019-01-21 2019-05-21 华北电力大学 Large-scale double-fed fan motor play synchronized oscillation emulation modelling method
CN113357101B (en) * 2021-07-05 2022-08-26 湖南工程学院 Method for identifying electromagnetic space-time coupling load of wind driven generator

Citations (4)

* Cited by examiner, † Cited by third party
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
CN102760191A (en) * 2012-06-26 2012-10-31 西安交通大学 Double-fed machine set wind power station equivalent modeling system and method based on rotation speed grouping
CN103870923A (en) * 2014-03-03 2014-06-18 华北电力大学 Information entropy condensation type hierarchical clustering algorithm-based wind power plant cluster aggregation method
CN105703396A (en) * 2016-03-21 2016-06-22 广西电网有限责任公司电力科学研究院 Multi-turbine characterization equivalent modeling method of double-fed wind power plant

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102682185B (en) * 2011-03-10 2014-07-09 华锐风电科技(集团)股份有限公司 Single wind turbine wind power prediction method

Patent Citations (4)

* Cited by examiner, † Cited by third party
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
CN102760191A (en) * 2012-06-26 2012-10-31 西安交通大学 Double-fed machine set wind power station equivalent modeling system and method based on rotation speed grouping
CN103870923A (en) * 2014-03-03 2014-06-18 华北电力大学 Information entropy condensation type hierarchical clustering algorithm-based wind power plant cluster aggregation method
CN105703396A (en) * 2016-03-21 2016-06-22 广西电网有限责任公司电力科学研究院 Multi-turbine characterization equivalent modeling method of double-fed wind power plant

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
双馈风电场并网系统的次同步振荡特性及抑制策略研究;赵斌;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20160115;C042-84
基于聚类算法的风电场动态等值问题研究;王聪;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20121015;C042-875
风电场动态等值模型的多机表征方法;米增强 等;《电工技术学报》;20100526;第25卷(第5期);162-169

Also Published As

Publication number Publication date
CN106168994A (en) 2016-11-30

Similar Documents

Publication Publication Date Title
CN106168994B (en) A kind of subsynchronous frequency equivalent modeling method of wind power plant
CN103701148B (en) Starting control method for accessing VSC-MTDC system of large wind farm
Jamal et al. Performance evaluation of wind turbine with doubly-fed induction generator
CN105703396B (en) A kind of double-fed wind generating field multimachine characterizes equivalent modeling method
CN109787239B (en) Alternating iteration load flow calculation method and system of load flow model and multi-frequency alternating current series-parallel system
CN111509714B (en) Impedance modeling-based offshore wind turbine group grid-connected resonance stability judging method
Wang et al. Dynamic equivalent modeling for wind farms with DFIGs using the artificial bee colony with K-means algorithm
CN106058937A (en) Generalized load modeling method for distribution network containing direct-drive permanent magnet wind farm
CN109672221A (en) A kind of directly driven wind-powered field Dynamic Equivalence for sub-synchronous oscillation analysis
Al-Bayati et al. Aggregated models of wind farms: Current methods and future trends
CN109657902A (en) Consider the double-fed blower Dynamic Equivalence of sub-synchronous oscillation characteristic
Wu et al. Equivalent modeling method for regional decentralized photovoltaic clusters based on cluster analysis
CN108879723B (en) Method for determining type and main participation area of internal resonance of offshore wind power plant
CN114513004A (en) New energy station equivalence method based on improved k-means algorithm and application
CN109787281A (en) Large-scale double-fed fan motor play synchronized oscillation emulation modelling method
Zhu et al. High-frequency resonance of DFIG-based wind generation under weak power network
CN103259269B (en) Wind power plant reactive power optimization configuration method based on double-fed induction fan reactive power generating capacity
CN107370183B (en) The subsynchronous stable region control method and controller of double-fed air-blower control parameter
CN115659601A (en) Parameter identification method and device for single-machine equivalent grid-connected model of double-fed wind power plant
CN103886178A (en) Parameter equating method for aggregation models of wind power plants with permanent-magnetic direct-drive units
Chouket et al. Wind Turbine PI Controller's Optimization Using PSO Algorithm
CN107221958B (en) A kind of wind power generating set backup power system capacity determining methods
Sun et al. A preliminary study on the Thevenin equivalent parameters of DFIG wind farms
Jiang et al. Analysis of the combined impact of small-scale wind generators and electric vehicles on future power networks
Gao et al. Improved extended kalman filter based dynamic equivalent method of DFIG wind farm cluster

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant