CN116544964B - Impedance optimization method of wind power generation system - Google Patents
Impedance optimization method of wind power generation system Download PDFInfo
- Publication number
- CN116544964B CN116544964B CN202310503296.8A CN202310503296A CN116544964B CN 116544964 B CN116544964 B CN 116544964B CN 202310503296 A CN202310503296 A CN 202310503296A CN 116544964 B CN116544964 B CN 116544964B
- Authority
- CN
- China
- Prior art keywords
- current
- impedance
- axis component
- expression
- net
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000010248 power generation Methods 0.000 title claims abstract description 15
- 239000003990 capacitor Substances 0.000 claims abstract description 17
- 230000014509 gene expression Effects 0.000 claims description 39
- 238000013178 mathematical model Methods 0.000 claims description 12
- 241000238633 Odonata Species 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000002955 isolation Methods 0.000 claims description 3
- 239000000178 monomer Substances 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/22—Current control, e.g. using a current control loop
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P27/00—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
- H02P27/04—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
- H02P27/06—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
- H02P27/08—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation
- H02P27/085—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation wherein the PWM mode is adapted on the running conditions of the motor, e.g. the switching frequency
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P27/00—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
- H02P27/04—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
- H02P27/06—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
- H02P27/08—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation
- H02P27/12—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation pulsing by guiding the flux vector, current vector or voltage vector on a circle or a closed curve, e.g. for direct torque control
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Power Engineering (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Linguistics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Water Supply & Treatment (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The invention discloses an impedance optimization method of a wind power generation system, which belongs to the field of power system stability, wherein a doubly-fed fan main circuit comprises the following components from left to right: the invention discloses an impedance optimization method for a doubly-fed fan system, which comprises a doubly-fed fan, a rotor side transformer T, a rotor side filter inductance Lr, a machine side converter MSC, a converter capacitor, a grid side converter GSC, a grid side filter inductance Lg, a grid side capacitor and a power grid, wherein the impedance optimization method is based on a control link of a virtual resistor of a doubly-fed fan rotor side control link, a control link of a doubly-fed grid side converter additional band-pass filter and an optimization debugging method of the doubly-fed grid side converter. The invention provides the impedance optimization scheme and the corresponding optimization and debugging method for the small wind working condition, which can pointedly solve the problem of subsynchronous oscillation under the small wind working condition, realize the impedance optimization of the doubly-fed fan, and have good technical effect and high practical value.
Description
Technical Field
The invention belongs to the field of power system stability, and particularly relates to an impedance optimization method of a wind power generation system.
Background
Under the goals of carbon peak and carbon neutralization, the large-scale development and utilization of new energy represented by wind energy and solar energy are quickened, the low-carbon transformation of energy power is facilitated, the power grid-connected installation 4757 ten thousand kW is newly increased by 2021, and the power grid-connected installation is increased by 16.6% in a same ratio. By the end of 2021, the wind power had accumulated installed capacity of 3.28 million kW, with an onshore wind turbine accumulated installed capacity of 3.02 million kW. High-power, high-voltage and long-distance conveying patterns are formed in the aspect of wind power. Wind power is sent out through a series compensation capacitor and is widely applied to practical engineering. However, the double-fed fan system has the problem of subsynchronous oscillation, and the subsynchronous oscillation of the system can be effectively restrained by carrying out impedance optimization on the double-fed fan. Therefore, the impedance optimization is carried out on the double-fed fan, and the suppression of the subsynchronous oscillation of the double-fed fan system is of great significance.
Li Penghan et al propose that the rotor side control link of the doubly-fed wind turbine utilizes the control strategy of the sliding mode variable structure to inhibit subsynchronous oscillation caused by the grid-connected subsynchronous control interaction of the doubly-fed wind turbine in the literature of feedback linearization sliding mode variable structure inhibition of subsynchronous control interaction of the doubly-fed wind turbine; su Tianyu et al in the literature of the phenomenon of grid-connected double-fed wind power plant subsynchronous/supersynchronous mixed oscillation and damping control scheme, aiming at subsynchronous/supersynchronous oscillation caused by series compensation grid connection of a double-fed wind power plant, a damper is added to control a rotor side of the double-fed wind power plant to inhibit the subsynchronous/supersynchronous oscillation; shan Bihan et al propose to use an additional damping control strategy in the rotor side control link to achieve subsynchronous resonance suppression in the literature "stator side analogue resistance based doubly fed wind farm subsynchronous oscillation suppression strategy study".
Under the working condition of small wind, namely the working condition of minimum active power (for example, the active power is the minimum power, 20% rated power, 50% rated power and the like), the impedance of the doubly-fed fan system is difficult to meet the technical requirements of the impedance characteristics of the new energy station. At present, a plurality of schemes are proposed for solving the problem of impedance optimization of the doubly-fed fans, but most of impedance optimization schemes based on rated power working conditions have poor optimization effect under small wind working conditions, are not suitable for impedance optimization of the doubly-fed fans under the small wind conditions, and are not designed for impedance optimization under the small wind working conditions; in addition, the existing scheme is not provided with a virtual resistor at the machine side, and a band-pass filter impedance optimization scheme is controlled and added at the network side.
Disclosure of Invention
The invention aims to provide an impedance optimization method of a wind power generation system, which can be used for carrying out impedance optimization on a doubly-fed fan under a small wind working condition and solving the problem of secondary super-synchronous oscillation under the small wind working condition.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows: an impedance optimization method of a wind power generation system, a main circuit of a doubly-fed wind power generation system comprises the following steps from left to right: double-fed fan, rotor side transformer T and rotor side filter inductance L r Side current transformer (MSC), current transformer capacitor, network side current transformer (GSC), network side filter inductance L g Network side capacitance and power grid.
Rotor side inductor current i r_abc The d-axis component i of the inductance current of the rotor side is obtained by collecting the inductance current from the rotor side line and transforming the coordinates dr And rotor side inductor current q-axis component i qr ；
Net side inductance current i g_abc Is acquired from a network side line and obtained after coordinate transformationNet side inductance current d-axis component i gd And net side inductor current q-axis component i gq ；
Network side capacitance voltage u g_abc The d-axis component u of the net-side capacitance voltage is obtained after the coordinate transformation is acquired from the net-side line gd And net side capacitor voltage q-axis component u gq ；
Rotor side phase angle theta r Frequency omega r Obtained via photoelectric encoder, net side phase angle theta g Frequency omega g Obtained via a Phase Locked Loop (PLL).
The modulation mode of MSC and GSC is space vector modulation (SVPWM), q-axis component control mode of MSC is torque outer ring-current inner ring double closed-loop control, d-axis component control mode is voltage outer ring-current inner ring double closed-loop control; the q-axis component control mode of the GSC is current closed-loop control, and the d-axis component control mode is voltage outer ring-current inner ring double closed-loop control.
Wherein the reference current i qr_ref From rotor-side reference torque T e_ref After adjustment, the regulator expression is: k omega g K is the active current coefficient, ω, calculated from the torque g Is the network side frequency; reference current i dr_ref From the q-axis component u of the net-side capacitor voltage gq After adjustment, the regulator expression is: 1/X m ，X m Excitation reactance at 50 Hz; reference current i d_ref By the capacitance voltage u of the current transformer dc Reference voltage u of converter capacitor dc_ref Obtained after passing through a third PI regulator.
The DFIG rotor side introduces virtual impedance, and the d-axis component i of the rotor side inductance current acquired in the rotor side line dr And rotor side inductor current q-axis component i qr And an additional signal is obtained through virtual impedance for the input component, and is added to the output link of the first PI regulator to play a role in impedance optimization. Wherein the virtual impedance comprises a resistance, an inductance or a capacitance, and the transfer function Z(s) of the virtual impedance controller is expressed as follows:
wherein R represents resistance, L represents inductance, and C represents capacitance.
Meanwhile, a band-pass filter is added in a network side control link, and a network side capacitance voltage d-axis component u acquired in a network side line gd And net side capacitor voltage q-axis component u gq The input component is input into the second PI regulator as negative feedback after passing through the band-pass filter, and plays roles in harmonic filtering and impedance optimization. The bandpass filter transfer function expression is:
wherein m is gain, k is damping coefficient, 0.707, omega 0 Center angular frequency=2pi f, f being the center frequency.
The impedance optimization debugging method of the fan comprises the following steps:
step S1, debugging parameters of band-pass filtering gain m and center frequency f by utilizing a dragonfly algorithm;
step S2, judging whether the system can stably run, if so, executing step S3, and if not, returning to execute step S1;
step S3, judging whether the impedance phase of the system in the frequency range of 0-H1 meets the requirement, if so, executing step S4, and if not, returning to execute step S1;
s4, debugging virtual impedance parameters including virtual resistance, capacitance or inductance by utilizing a dragonfly algorithm;
step S5, judging whether the system can stably run, if so, executing step S6, and if not, returning to execute step S4;
and S6, judging whether the impedance phase of the system in all frequency bands meets the requirement, if so, ending the optimization debugging, and if not, returning to the step S4.
The dragonfly algorithm in the optimization debugging method comprises the following steps:
step 1, setting related parameters, and randomly initializing a monomer position vector and a step length vector in an algorithm;
step 2, judging whether the current iteration times t is smaller than the maximum iteration times M, if yes, carrying out step 3; if not, outputting a parameter debugging result, and ending the debugging;
step 3, updating weight and field radius, calculating individual fitness, updating separation weight s, and adding Ji Quan weight a, cohesion weight c, food factor f and enemy factor e;
step 4, judging whether a field exists, namely whether a target solution exists in the field, if so, updating the step length vector and the position vector by using mathematical model expressions of the step length vector and the position vector, and then carrying out step 5; if not, directly carrying out the step 5;
the mathematical model expression of the step length vector is as follows:
ΔX t+1 ＝(sS i +aA i +cC i +fF i +eE i )+ωΔX t
wherein ω is an inertial weight; ΔX t Step length vector when the iteration times t are; s is S i For the ith individual isolation, the expression is:
x is the current individual position; xj is the j-th adjacent individual position; n is the number of adjacent individuals of the population;
A i aligned for the ith individual, the expression is:
V j the speed of the jth adjacent individual in the population;
C i cohesive for the ith individual, the expression is:
F i for the i-th individual food source, the expression is:
F i ＝X + -X
X + is the current position of the food;
E i for the i-th individual enemy location, the expression is:
E i ＝X - -X
X - is the current location of the enemy;
the mathematical model expression of the position vector is:
X t+1 ＝X t +ΔX t+1
step 5, updating the step length vector by using the mathematical model expression of the step length vector in the step 4, and returning to the step 2 after updating the position vector by using Levy flight; the Levy flight update position vector expression is:
X t+1 ＝X t +Levy(d)×X t
wherein d is the dimension of the position vector;
wherein r is 1 、r 2 Is [0,1]2 random numbers in (a); beta is a constant; the calculation formula of sigma is:
wherein: Γ (x) = (x-1) +.! .
By adopting the technical scheme, the invention has the following technical effects:
the invention provides the impedance optimization scheme and the corresponding optimization and debugging method for the small wind working condition, which can pointedly solve the problem of subsynchronous oscillation under the small wind working condition, realize the impedance optimization of the doubly-fed fan, and have good technical effect and high practical value.
Drawings
FIG. 1 is a schematic diagram of a system circuit and control scheme of the present invention;
FIG. 2 is a rotor side converter current loop control diagram of the present invention;
FIG. 3 is a diagram of the current loop control of the grid-side current transformer of the present invention;
FIG. 4 is a flow chart of an optimization debugging method of the present invention;
FIG. 5 is a diagram of the full band impedance phase requirement of the present invention;
fig. 6 is a comparison of the impedance curves before and after optimization.
Detailed Description
The technical scheme of the invention is further described in detail below with reference to the accompanying drawings:
as shown in fig. 1, an impedance optimization method of a wind power generation system, a main circuit of a doubly-fed wind power generation system includes, from left to right: double-fed fan, rotor side transformer T and rotor side filter inductance L r Side current transformer (MSC), current transformer capacitor, network side current transformer (GSC), network side filter inductance L g Network side capacitance and power grid. In the figure, K: an active current coefficient calculated from the torque; x is X m : excitation reactance at 50 Hz; k (K) 1r : leakage reactance coefficient of the motor; l (L) s : self-inductance of the stator winding; l (L) m The stator and rotor windings are mutually inductive.
Rotor side inductor current i r_abc The d-axis component i of the inductance current of the rotor side is obtained by collecting the inductance current from the rotor side line and transforming the coordinates dr And rotor side inductor current q-axis component i qr ；
Net side inductance current i g_abc The d-axis component i of the net-side inductance current is obtained after the coordinate transformation is acquired from the net-side line gd And net side inductor current q-axis component i gq ；
Network side capacitance voltage u g_abc The d-axis component u of the net-side capacitance voltage is obtained after the coordinate transformation is acquired from the net-side line gd And net side capacitor voltage q-axis component u gq ；
Rotor side phase angle theta r Frequency omega r Obtained via photoelectric encoder, net side phase angle theta g Frequency omega g Obtained via a Phase Locked Loop (PLL).
The modulation mode of MSC and GSC is space vector modulation (SVPWM), q-axis component control mode of MSC is torque outer ring-current inner ring double closed-loop control, d-axis component control mode is voltage outer ring-current inner ring double closed-loop control; the q-axis component control mode of the GSC is current closed-loop control, and the d-axis component control mode is voltage outer ring-current inner ring double closed-loop control.
Wherein the reference current i qr_ref From rotor-side reference torque T e_ref After adjustment, the regulator expression is: k omega g K is the active current coefficient, ω, calculated from the torque g Is the network side frequency; reference current i dr_ref From the q-axis component u of the net-side capacitor voltage gq After adjustment, the regulator expression is: 1/X m ，X m Excitation reactance at 50 Hz; reference current i d_ref By the capacitance voltage u of the current transformer dc Reference voltage u of converter capacitor dc_ref Obtained after passing through a third PI regulator.
The current loop control process of MSC and GSC is:
the MSC current loop control process comprises the following steps: as shown in fig. 2, the reference current i qr_ref /i dr_ref Q/d axis component i of inductor current on rotor side qr /i dr After difference, the output signal is adjusted by PI, and i qr /i dr The signals obtained by the virtual impedances are added to obtain a modulated wave. Wherein the phase angle of the modulated wave is theta g -θ r ，θ g For the net side phase angle, θ r The rotor side phase angle; the virtual impedance includes a resistance, an inductance, or a capacitance, and the transfer function Z(s) of the virtual impedance controller is expressed as:
wherein R represents resistance, L represents inductance, and C represents capacitance.
The GSC current loop control process comprises the following steps: as shown in FIG. 3, the reference current is 0/i d_ref Q/d axis component i of inductance current at network side gq /i gd After the difference is made, the filtered signal fed back is subtracted, and the PI regulation is carried outObtain output signal, and obtain q/d axis component u of network side capacitor voltage gq /u gd Subtracting the output signal thereof to obtain a modulated wave. Wherein the phase angle of the modulated wave is theta g ，θ g Is the network side phase angle; the negative feedback signal is the d-axis component u of the net-side capacitance voltage collected from the net-side line gd And net side capacitor voltage q-axis component u gq The band-pass filter transfer function expression is as follows:
wherein m is gain, k is damping coefficient, 0.707, omega 0 Center angular frequency=2pi f, f being the center frequency.
As shown in fig. 4, the impedance optimization and debugging method of the fan of the present invention comprises:
step S1, debugging parameters of band-pass filtering gain m and center frequency f by utilizing a dragonfly algorithm;
step S2, judging whether the system can stably run, if so, executing step S3, and if not, returning to execute step S1;
step S3, judging whether the impedance phase of the system in the frequency range of 0-H1 meets the requirement, if yes, executing step S4, and if not, returning to execute step S1, wherein the value of H1 is 100, the value range of the full frequency range is 0-2000, and the full frequency range impedance phase requirement is shown in figure 5;
s4, debugging virtual impedance parameters including virtual resistance, capacitance or inductance by utilizing a dragonfly algorithm;
step S5, judging whether the system can stably run, if so, executing step S6, and if not, returning to execute step S4;
and S6, judging whether the impedance phase of the system in all frequency bands meets the requirement, if so, ending the optimization debugging, and if not, returning to the step S4.
The method for debugging parameters by using the dragonfly algorithm comprises the following steps of:
step 1, setting related parameters, and randomly initializing a monomer position vector and a step length vector in an algorithm;
step 2, judging whether the current iteration times t is smaller than the maximum iteration times M, if yes, carrying out step 3; if not, outputting a parameter debugging result, and ending the debugging;
step 3, updating weight and field radius, calculating individual fitness, updating separation weight s, and adding Ji Quan weight a, cohesion weight c, food factor f and enemy factor e;
step 4, judging whether a field exists, namely whether a target solution exists in the field, if so, updating the step length vector and the position vector by using mathematical model expressions of the step length vector and the position vector, and then carrying out step 5; if not, directly carrying out the step 5;
the mathematical model expression of the step length vector is as follows:
ΔX t+1 ＝(sS i +aA i +cC i +fF i +eE i )+ωΔX t
wherein ω is an inertial weight; ΔX t Step length vector when the iteration times t are; s is S i For the ith individual isolation, the expression is:
x is the current individual position; xj is the j-th adjacent individual position; n is the number of adjacent individuals of the population;
A i aligned for the ith individual, the expression is:
V j the speed of the jth adjacent individual in the population;
C i cohesive for the ith individual, the expression is:
F i for the i-th individual food source, the expression is:
F i ＝X + -X
X + is the current position of the food;
E i for the i-th individual enemy location, the expression is:
E i ＝X - -X
X - is the current location of the enemy;
the mathematical model expression of the position vector is:
X t+1 ＝X t +ΔX t+1
step 5, updating the step length vector by using the mathematical model expression of the step length vector in the step 4, and returning to the step 2 after updating the position vector by using Levy flight; the Levy flight update position vector expression is:
X t+1 ＝X t +Levy(d)×X t
wherein d is the dimension of the position vector;
wherein r is 1 、r 2 Is [0,1]2 random numbers in (a); beta is a constant; the calculation formula of sigma is:
wherein: Γ (x) = (x-1) +.! .
The comparison result is shown in fig. 6, the blue line in fig. 6 is the result before optimization, the purple line is the result after optimization, and the impedance phase requirement of 0-100Hz is set in the red dotted line, so that the impedance curve after optimization through the scheme is more in the red dotted line, and the impedance phase requirement is more satisfied.
Claims (6)
1. The method comprises the following steps ofThe impedance optimization method of the wind power generation system comprises the following steps of from left to right: double-fed fan, rotor side transformer T, rotor side filter inductance Lr, machine side converter MSC, converter electric capacity, net side converter GSC, net side filter inductance Lg, net side electric capacity and electric wire netting, its characterized in that: the rotor side of the doubly-fed wind turbine introduces virtual impedance, and the d-axis component i of the rotor side current is collected in the rotor side line dr And rotor-side current q-axis component i qr The additional signal is obtained through virtual impedance for the input component and added to the output link of the first PI regulator, so that the function of impedance optimization is achieved; meanwhile, a band-pass filter is added in a network side control link, and a network side voltage d-axis component u acquired in a network side line gd And net side voltage q-axis component u gq The input component is input into a second PI regulator as negative feedback after passing through a band-pass filter; the optimization debugging method comprises the following steps:
step S1, debugging parameters of band-pass filtering gain m and center frequency f by utilizing a dragonfly algorithm;
step S2, judging whether the system can stably run, if so, executing step S3, and if not, returning to execute step S1;
step S3, judging whether the impedance phase of the system in the frequency range of 0-H1 meets the requirement, if so, executing step S4, and if not, returning to execute step S1;
s4, debugging virtual impedance parameters including virtual resistance, capacitance or inductance by utilizing a dragonfly algorithm;
step S5, judging whether the system can stably run, if so, executing step S6, and if not, returning to execute step S4;
and S6, judging whether the impedance phase of the system in all frequency bands meets the requirement, if so, ending the optimization debugging, and if not, returning to the step S4.
2. A method of optimizing the impedance of a wind power generation system according to claim 1, wherein: rotor side inductor current i r_abc The d-axis component i of the inductance current of the rotor side is obtained by collecting the inductance current from the rotor side line and transforming the coordinates dr And a rotorSide inductor current q-axis component i qr The method comprises the steps of carrying out a first treatment on the surface of the Net side inductance current i g_abc The d-axis component i of the net-side inductance current is obtained after the coordinate transformation is acquired from the net-side line gd And net side inductor current q-axis component i gq The method comprises the steps of carrying out a first treatment on the surface of the Network side capacitance voltage u g_abc The d-axis component u of the net-side capacitance voltage is obtained after the coordinate transformation is acquired from the net-side line gd And net side capacitor voltage q-axis component u gq The method comprises the steps of carrying out a first treatment on the surface of the Rotor side phase angle theta r Frequency omega r Obtained via photoelectric encoder, net side phase angle theta g Frequency omega g Obtained via a phase locked loop.
3. A method of optimizing the impedance of a wind power generation system according to claim 1, wherein: the virtual impedance includes a resistance, an inductance, or a capacitance, and the transfer function Z(s) of the virtual impedance controller is expressed as:
wherein R represents resistance, L represents inductance and C represents capacitance;
the bandpass filter transfer function expression is:
wherein m is gain, k is damping coefficient, 0.707, omega 0 Center angular frequency=2pi f, f being the center frequency.
4. A method of optimizing the impedance of a wind power generation system according to claim 1, wherein: the modulation mode of MSC and GSC is space vector modulation, q-axis component control mode of MSC is torque outer ring-current inner ring double closed-loop control, d-axis component control mode is voltage outer ring-current inner ring double closed-loop control; the q-axis component control mode of the GSC is current closed-loop control, and the d-axis component control mode is voltage outer ring-current inner ring double closed-loop control.
5. A method of optimizing the impedance of a wind power generation system according to claim 1, wherein: reference current i qr_ref From rotor-side reference torque T e_ref After adjustment, the regulator expression is: k omega g K is the active current coefficient, ω, calculated from the torque g Is the network side frequency; reference current i dr_ref From the q-axis component u of the net-side capacitor voltage gq After adjustment, the regulator expression is: 1/X m ，X m Excitation reactance at 50 Hz; reference current i d_ref By the capacitance voltage u of the current transformer dc Reference voltage u of converter capacitor dc_ref Obtained after passing through a third PI regulator.
6. A method of optimizing the impedance of a wind power generation system according to claim 1, comprising the steps of: the steps of applying dragonfly algorithm debugging parameters are as follows:
step 1, setting related parameters, and randomly initializing a monomer position vector and a step length vector in an algorithm;
step 2, judging whether the current iteration times t is smaller than the maximum iteration times M, if yes, performing a step S3; if not, outputting a parameter debugging result, and ending the debugging;
step 3, updating weight and field radius, calculating individual fitness, updating separation weight s, and adding Ji Quan weight a, cohesion weight c, food factor f and enemy factor e;
step 4, judging whether a field exists, namely whether a target solution exists in the field, if so, updating the step length vector and the position vector by using mathematical model expressions of the step length vector and the position vector, and then carrying out step 5; if not, directly carrying out the step 5;
the mathematical model expression of the step length vector is as follows:
ΔX t+1 ＝(sS i +aA i +cC i +fF i +eE i )+ωΔX t
wherein ω is an inertial weight; deltaX t Step length vector when the iteration times t are; s is S i For the ith individual isolation, the expression is:
x is the current individual position; xj is the j-th adjacent individual position; n is the number of adjacent individuals of the population;
A i aligned for the ith individual, the expression is:
V j the speed of the jth adjacent individual in the population;
C i cohesive for the ith individual, the expression is:
F i for the i-th individual food source, the expression is:
F i ＝X + -X
X + is the current position of the food;
E i for the i-th individual enemy location, the expression is:
E i ＝X - -X
X - is the current location of the enemy;
the mathematical model expression of the position vector is:
X t+1 ＝X t +ΔX t+1
step 5, updating the step length vector by using the mathematical model expression of the step length vector in the step 4, and returning to the step 2 after updating the position vector by using Levy flight; the Levy flight update position vector expression is:
X t+1 ＝X t +Levy(d)×X t
wherein d is the dimension of the position vector;
wherein r is 1 、r 2 Is [0,1]2 random numbers in (a); beta is a constant; the calculation formula of sigma is:
wherein: Γ (x) = (x-1) +.! .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310503296.8A CN116544964B (en) | 2023-05-06 | 2023-05-06 | Impedance optimization method of wind power generation system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310503296.8A CN116544964B (en) | 2023-05-06 | 2023-05-06 | Impedance optimization method of wind power generation system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116544964A CN116544964A (en) | 2023-08-04 |
| CN116544964B true CN116544964B (en) | 2023-12-01 |
Family
ID=87444727
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310503296.8A Active CN116544964B (en) | 2023-05-06 | 2023-05-06 | Impedance optimization method of wind power generation system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116544964B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107017646A (en) * | 2017-05-25 | 2017-08-04 | 东南大学 | The double-fed blower fan sub-synchronous oscillation suppression method controlled based on virtual impedance |
| CN109120001A (en) * | 2018-09-29 | 2019-01-01 | 华中科技大学 | Double-fed fan motor field grid-connected system sub-synchronous oscillation suppression method based on virtual resistance |
| CN109449958A (en) * | 2018-11-29 | 2019-03-08 | 云南电网有限责任公司 | A kind of double-fed blower grid-connected system method for analyzing stability |
| CN109698502A (en) * | 2018-12-26 | 2019-04-30 | 电子科技大学 | A kind of multi-inverter micro-capacitance sensor harmonic resonance suppressing method of adaptive virtual impedance remodeling |
| CN111245013A (en) * | 2020-02-27 | 2020-06-05 | 华北电力大学 | DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling |
| CN112653184A (en) * | 2020-12-15 | 2021-04-13 | 上海交通大学 | Method, device, terminal and medium for identifying black box model of wind power generation equipment |
| CN112671009A (en) * | 2021-01-13 | 2021-04-16 | 天津大学 | Doubly-fed fan subsynchronous oscillation suppression method with additional damping controller |
| CN113328458A (en) * | 2020-12-15 | 2021-08-31 | 新疆金风科技股份有限公司 | Feedforward control method and device for grid-connected inverter and computer readable storage medium |
| CN115473277A (en) * | 2022-10-19 | 2022-12-13 | 华北电力科学研究院有限责任公司 | Impedance shaping method and device for near-power-frequency doubly-fed wind turbine generator |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2554954B (en) * | 2016-10-17 | 2018-11-21 | Zhong Qingchang | Operating doubly-fed induction generators as virtual synchronous generators |
-
2023
- 2023-05-06 CN CN202310503296.8A patent/CN116544964B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107017646A (en) * | 2017-05-25 | 2017-08-04 | 东南大学 | The double-fed blower fan sub-synchronous oscillation suppression method controlled based on virtual impedance |
| CN109120001A (en) * | 2018-09-29 | 2019-01-01 | 华中科技大学 | Double-fed fan motor field grid-connected system sub-synchronous oscillation suppression method based on virtual resistance |
| CN109449958A (en) * | 2018-11-29 | 2019-03-08 | 云南电网有限责任公司 | A kind of double-fed blower grid-connected system method for analyzing stability |
| CN109698502A (en) * | 2018-12-26 | 2019-04-30 | 电子科技大学 | A kind of multi-inverter micro-capacitance sensor harmonic resonance suppressing method of adaptive virtual impedance remodeling |
| CN111245013A (en) * | 2020-02-27 | 2020-06-05 | 华北电力大学 | DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling |
| CN112653184A (en) * | 2020-12-15 | 2021-04-13 | 上海交通大学 | Method, device, terminal and medium for identifying black box model of wind power generation equipment |
| CN113328458A (en) * | 2020-12-15 | 2021-08-31 | 新疆金风科技股份有限公司 | Feedforward control method and device for grid-connected inverter and computer readable storage medium |
| CN112671009A (en) * | 2021-01-13 | 2021-04-16 | 天津大学 | Doubly-fed fan subsynchronous oscillation suppression method with additional damping controller |
| CN115473277A (en) * | 2022-10-19 | 2022-12-13 | 华北电力科学研究院有限责任公司 | Impedance shaping method and device for near-power-frequency doubly-fed wind turbine generator |
Non-Patent Citations (4)
| Title |
|---|
| 基于定子侧模拟电阻的双馈风电场次同步振荡抑制策略研究;单碧涵 等;《 电力系统保护与控制》(第5期);第10-16页 * |
| 弱电网下双馈风电机组电网电压扰动补偿控制策略;张学广 等;《电力系统自动化》(第6期);第146-154页 * |
| 微电网变换器可编程虚拟无源网络原理及分析;郭小强 等;《电力电子技术》(第10期);第30-32页 * |
| 蜻蜓算法研究综述;池建华 等;《 自动化与信息工程》(第3期);第7-14页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116544964A (en) | 2023-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109120001B (en) | Doubly-fed wind power plant grid-connected system subsynchronous oscillation suppression method based on virtual resistor | |
| CN108631338B (en) | Method for inhibiting grid-connected subsynchronous oscillation of doubly-fed wind power plant | |
| CN110797883B (en) | Wind power plant flexible direct grid-connected system subsynchronous oscillation suppression method based on impedance method | |
| CN107863783B (en) | Doubly-fed wind turbine virtual synchronous control method | |
| CN104362668B (en) | The control method of double-fed wind power generator during a kind of Voltage unbalance/harmonic distortion | |
| CN112448398B (en) | Stator side analog resistance-based doubly-fed wind power plant subsynchronous oscillation suppression method | |
| CN109473997A (en) | A kind of double-fed fan motor field sub-synchronous oscillation suppression method based on source net Collaborative Control | |
| Kerrouche et al. | Fractional-order sliding mode control for D-STATCOM connected wind farm based DFIG under voltage unbalanced | |
| CN105790269B (en) | Suppress the method and device of subsynchronous resonance by double-fed fan stator side converter | |
| CN111130142B (en) | Virtual capacitance control method for double-fed wind power generation grid-connected system | |
| CN111969649B (en) | Control method and system for improving power transmission limit of double-fed fan in weak grid | |
| CN105790270A (en) | Method and device for suppressing subsynchronous resonance through doubly-fed fan rotor side converter | |
| CN116014820B (en) | Asymmetric fault low-voltage ride-through control method and system based on super capacitor | |
| CN111130136A (en) | Subsynchronous oscillation suppression method based on additional virtual impedance control | |
| CN111245013B (en) | DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling | |
| CN112260302A (en) | Energy storage PCS multi-machine parallel resonance active suppression method of virtual RC branch circuit | |
| CN116073398A (en) | Additional damping control system for improving synchronous stability of doubly-fed wind turbine generator | |
| CN116054187A (en) | Subsynchronous oscillation suppression method for wind power generated through soft direct grid connection | |
| CN117543732B (en) | Active power control method and device for virtual synchronous generator of network-structured converter | |
| CN116544964B (en) | Impedance optimization method of wind power generation system | |
| CN116995693A (en) | Method for controlling energy storage virtual inertia of network super capacitor | |
| CN113937789B (en) | Voltage source type double-fed fan feedforward damping control method based on fractional order filtering | |
| CN113193574B (en) | Subsynchronous oscillation suppression method and device based on impedance characteristics of ports of doubly-fed wind turbine | |
| CN114844041A (en) | Feedforward control method for enhancing system stability and filtering power grid higher-order background harmonic | |
| CN113394824A (en) | MMC alternating current active power rapid regulation and control method without direct current side disturbance |
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 |