CN112448399A - Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance - Google Patents

Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance Download PDF

Info

Publication number
CN112448399A
CN112448399A CN201910810798.9A CN201910810798A CN112448399A CN 112448399 A CN112448399 A CN 112448399A CN 201910810798 A CN201910810798 A CN 201910810798A CN 112448399 A CN112448399 A CN 112448399A
Authority
CN
China
Prior art keywords
rotor
inductance
doubly
control
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910810798.9A
Other languages
Chinese (zh)
Other versions
CN112448399B (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201910810798.9A priority Critical patent/CN112448399B/en
Publication of CN112448399A publication Critical patent/CN112448399A/en
Application granted granted Critical
Publication of CN112448399B publication Critical patent/CN112448399B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Abstract

The invention discloses a method for inhibiting doubly-fed wind power plant subsynchronous oscillation based on analog inductance, which comprises the following steps: firstly, establishing a system model of a doubly-fed wind power plant which is connected with the grid through a series compensation capacitor; then applying an analog inductance control strategy to the inner loop control of a double-fed fan rotor side controller to control the output voltage of the rotor; then, carrying out self-adaptive on-line setting on the simulated inductance value by utilizing a PSO algorithm, and determining the optimal value of the simulated inductance under a specific working condition; and finally, applying an analog inductance control strategy and a rotor side converter control strategy to the control of the rotor side of the doubly-fed fan to complete the suppression of the subsynchronous oscillation of the system. The invention improves the electrical damping in the subsynchronous frequency band, improves the stability of the system and can meet the subsynchronous oscillation suppression requirement under multiple working conditions.

Description

Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance
Technical Field
The invention belongs to the technical field of stability control of wind driven generators, and particularly relates to a method for inhibiting doubly-fed wind power plant subsynchronous oscillation based on analog inductance.
Background
From the viewpoint of sustainable development, fossil energy, which is currently relied on by mankind, is eventually exhausted, and in order to solve the current energy and environmental problems, the development and utilization of clean energy to replace the conventional energy is an inevitable trend of human social development. The wind energy is an ideal clean energy which is inexhaustible, environment-friendly and pollution-free. The wind power generation technology is mature, can be developed in a large scale, has good commercial prospect, and is the main form of the current new energy power generation. When wind energy resources are developed and utilized, different types of wind driven generators are used, and in the wind power plant which is put into operation at present, the doubly-fed induction wind driven generator is the most widely applied type. Due to the fact that resources and load centers are distributed in a reverse mode, large-capacity and long-distance wind power delivery is imperative. The series compensation capacitor technology is widely applied, but the problem of subsynchronous oscillation of a wind power system can be induced, and the safe and stable operation of a large-scale wind power base and a delivery system is influenced.
Due to uniqueness of a wind power plant, many control methods for obtaining ideal effects in thermal power plant subsynchronous oscillation suppression may be difficult to obtain the ideal effects in the wind power plant. In actual operation, different operation conditions exist, so that the self-adaption problem of the designed controller must be considered. At present, researches on a subsynchronous oscillation mechanism of a doubly-fed wind power plant are not completely thorough, and measures for inhibiting subsynchronous oscillation are relatively single. Therefore, further research on the suppression measures of the doubly-fed wind power plant subsynchronous oscillation is needed, and the method has important significance for stable grid connection of the wind power plant.
Disclosure of Invention
The invention aims to provide a method for inhibiting the sub-synchronous oscillation of a doubly-fed wind power plant based on analog inductance, which can inhibit the sub-synchronous oscillation caused by series compensation capacitance and improve the stability of a doubly-fed wind power plant through a series compensation capacitance grid-connected system.
The technical solution for realizing the purpose of the invention is as follows: a doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance comprises the following steps:
step 1, establishing a system model of a doubly-fed wind power plant which is connected with the grid through a series compensation capacitor;
step 2, applying an analog inductance control strategy to inner loop control of a double-fed fan rotor side controller to control the output voltage of the rotor;
step 3, carrying out self-adaptive on-line setting on the simulated inductance value by utilizing a PSO algorithm, and determining the optimal value of the simulated inductance under a specific working condition;
and 4, applying an analog inductance control strategy and a rotor side converter control strategy to the control of the rotor side of the doubly-fed wind turbine to complete the suppression of the subsynchronous oscillation of the system.
Further, the establishing of the system model of the doubly-fed wind farm through the series compensation capacitor grid connection in the step 1 is as follows:
step 1.1, obtaining the corresponding relation among the rotating speed, the output power and the wind speed of the rotor of the doubly-fed wind turbine according to a wind energy capture model;
step 1.2, writing out an electromagnetic equation of the system in a dq coordinate system as follows:
step 1.2.1, the stator and the rotor in the doubly-fed induction generator adopt the motor convention, the positive direction is the current flowing direction, and the rotating direction is the same as the positive direction of the electromagnetic torque; according to the coordinate transformation principle, current, voltage and magnetic flux variables in the stator and the rotor are transformed from a three-phase static coordinate to a dq coordinate system, and the electromagnetic relation in the stator and the rotor is obtained as follows:
Figure BDA0002184983430000021
Figure BDA0002184983430000022
Figure BDA0002184983430000023
wherein ω is1For synchronizing the angular velocities of rotation, ωs＝ω1-ωrIs the slip angular velocity, omegarIs the rotor rotational angular velocity; u. ofds、uqsIs the stator dq-axis voltage, udr、uqrIs the rotor dq axis voltage; i.e. ids、iqsIn order to obtain the stator dq-axis current,idr、iqris the rotor dq axis current; psids、ψqsFor stator dq axis flux linkage, psidr、ψqrA rotor dq axis flux linkage; rsIs stator resistance, RrIs the rotor resistance; l islsFor stator leakage inductance, LlrFor rotor leakage inductance, LmFor coaxial mutual inductance of stator and rotor, Ls、LrEquivalent self-inductance of stator and rotor windings;
step 1.2.2, controlling a target of a double-fed fan rotor side converter by tracking the maximum wind energy, keeping the frequency of generated energy constant and controlling reactive output; fixing the stator flux linkage direction on a d axis to realize active and reactive decoupling; and (3) adopting PI control to regulate the rotor voltage, and writing a relation column expressed by a control block diagram into a differential equation:
Figure BDA0002184983430000031
Figure BDA0002184983430000032
Figure BDA0002184983430000033
Figure BDA0002184983430000034
Figure BDA0002184983430000035
Figure BDA0002184983430000036
Figure BDA0002184983430000037
Figure BDA0002184983430000038
wherein Teref、QsrefReference value for electromagnetic torque and stator reactive power, KTe、KQsProportional coefficient, K, for the external power loop PI control of the rotor-side converteriq、KidProportional coefficient for PI control of current loop in the rotor side converter; t isTe、TQsIntegral coefficient, T, for rotor side converter outer power loop PI controliq、TidIntegral coefficient of current loop PI control in the rotor side converter; i.e. iqrref、idrrefRotor dq axis current reference, x, for rotor side converter outer loop PI control output1～x4Four intermediate variables were introduced.
Further, the step 2 of applying the analog inductance control strategy to the inner loop control of the double-fed fan rotor side converter to control the rotor output voltage specifically comprises the following steps:
step 2.1, detecting external wind speed to obtain a torque reference value corresponding to the wind speed, determining a reactive reference value of a stator according to a control target, establishing a relevant state equation of a simultaneous system, and calculating a steady-state expected value of a current dq axis component of a rotor of the doubly-fed fan by combining a corresponding control strategy
Figure BDA0002184983430000039
Step 2.2, after the actual value of the rotor current dq axis component is subjected to difference gain with the steady-state expected value, the difference gain is used as an additional signal to act on the input of the rotor voltage, and the formula is as follows:
Figure BDA00021849834300000310
and 2.3, enabling the rotor voltage output by the rotor side converter to meet the relation of an equation (13) under the action of the additional signal of the equation (12):
Figure BDA0002184983430000041
wherein G isPNamely the analog inductance.
Further, the step 3 of utilizing the PSO algorithm to perform adaptive online tuning on the analog inductance value, and determining the optimal value of the analog inductance under the specific working condition, specifically as follows:
step 3.1, writing a system state equation in a row, linearizing the system state equation at a balance point to obtain a small signal equation:
dΔX/dt＝AΔX (14)
in the formula, X is the state variable of the system, and the number of X is 20; a is a Jacobian matrix determined by system parameters and balance point states;
step 3.2, calculating the subsynchronous oscillation frequency f of the systemnThe formula is as follows:
Figure BDA0002184983430000042
in the formula (f)1For mains frequency, XCFor series compensation of capacitive reactance of capacitor, X∑Is the sum of effective inductive reactance in the system;
corresponding to a power oscillation frequency of fSSR＝f1-fnBased on the characteristic value, determining a corresponding subsynchronous oscillation mode in the characteristic value of the Jacobian matrix A;
step 3.3 at GPWhen the value is changed, two groups of characteristic value real parts are changed, wherein one group corresponds to a subsynchronous mode, and the corresponding characteristic value is sigma1±jω1(ii) a The other group corresponds to low-frequency mode, and the corresponding characteristic value is sigma2±jω2The PSO algorithm is adopted to carry out self-adaptive setting on the analog inductance value, the primary objective of the algorithm is to inhibit subsynchronous oscillation, and therefore the damping ratio corresponding to the subsynchronous modal characteristic value is taken as a fitness function, as shown in the formula (16):
Figure BDA0002184983430000043
step 3.4, limiting the real part of the characteristic value corresponding to the low-frequency mode in a certain range to obtain a constraint of the algorithm:
σ2≤σ2max (17)
in the formula sigma2maxA constant less than zero;
and 3.5, updating the position and the search speed of the ith particle according to the formula (18):
Figure BDA0002184983430000044
in the formula, k is the current iteration number; w is the inertial weight; c. C1And c2Is a learning factor; r is1And r2A random number in the range of 0 to 1;
Figure BDA0002184983430000051
is the current best position of the ith particle,
Figure BDA0002184983430000052
is the current global optimum position;
3.6, establishing a dynamic adjustment model of the inertial weight, and searching by using the fixed weight at the initial stage of the algorithm; and after the set number of iterations, varying w according to equation (19):
Figure BDA0002184983430000053
in the formula, wmaxAnd wminIs a fixed value, fiFor the current particle adaptation value, favAverage of the adapted values for all particles, fbestAdapting the value for the global optimum particle;
step 3.7, simulating the inductor GPThe optimum range of values and particle velocity are constrained:
Figure BDA0002184983430000054
wherein G isPminAnd GPmaxTo a fixed value selected as required, vmaxAnd vminThe upper and lower limits of the particle search speed can be selected according to typical values.
Further, the step 4 of applying the analog inductance control strategy and the rotor-side converter control strategy to the control of the rotor side of the doubly-fed wind turbine to complete the suppression of the subsynchronous oscillation of the system is as follows:
the simulated inductance G obtained by PSO algorithm settingPThe value is applied to the control of the converter, the control of the analog inductance of the double-fed fan is realized, and the suppression of the subsynchronous oscillation of the system is completed.
Compared with the prior art, the invention has the remarkable advantages that: (1) the method inhibits subsynchronous oscillation caused by series compensation capacitance, and improves the stability of a doubly-fed wind power plant through a series compensation capacitance grid-connected system; (2) the invention has few parameters needing to be set, and avoids complex parameter setting; (3) the PI control of the double-fed fan rotor side converter is improved, the basic control structure of the RSC converter is not changed, and the requirements of current fan manufacturers are met; (4) the method can suppress subsynchronous oscillation under multiple working conditions.
Drawings
Fig. 1 is a flow chart of a method for suppressing the subsynchronous oscillation of a doubly-fed wind power plant based on an analog inductor according to the present invention.
FIG. 2 is a schematic structural diagram of a doubly-fed wind power plant grid-connected model through series compensation capacitors.
Fig. 3 is a control block diagram of the double-fed wind turbine rotor side converter in the invention.
Fig. 4 is a control block diagram of the doubly-fed wind turbine rotor side converter based on the analog inductance in the invention.
Fig. 5 is a schematic diagram of the principle of analog inductance adaptive parameter tuning using the PSO algorithm in the present invention.
Fig. 6 is an oscillation curve diagram of the electromagnetic torque of the doubly-fed wind turbine when the wind speed is kept fixed and the series compensation degree is changed under the conventional PI control scheme in the embodiment of the invention.
Fig. 7 is an oscillation curve of the electromagnetic torque of the doubly-fed wind turbine under three different working conditions when the control strategy based on the simulated inductance is adopted in the embodiment of the present invention, wherein (a) is the wind speed of 7m/s and the series compensation degree of 70%; (b) the wind speed is 7m/s, and the series compensation degree is 80 percent; (c) the wind speed is 8m/s, and the degree of serial compensation is 80%.
Detailed Description
As shown in fig. 1, the invention provides a method for suppressing the sub-synchronous oscillation of a doubly-fed wind power plant based on analog inductance, which comprises the following steps:
step 1, establishing a system model of a doubly-fed wind power plant connected through a series compensation capacitor, which comprises the following specific steps:
step 1.1, obtaining the corresponding relation among the rotating speed, the output power and the wind speed of the rotor of the doubly-fed wind turbine according to the wind energy capture model. As shown in FIG. 2, a 100MW wind farm is formed by 50 fans of 2MW capacity per unit, where X istg、XTIs a net side smoothing reactor and a transformer reactor, RL、XLAnd XCThe resistance, reactance and series compensation capacitance reactance of the power transmission line are obtained; in order to fully utilize the wind speed, the unit operates under the maximum wind energy tracking strategy, and the corresponding relation of the rotor speed, the output power and the wind speed of the doubly-fed wind turbine is shown in table 1, wherein V iswIs the wind speed, ωrAs the rotor speed, PwFor wind turbine output power, TwOutputting torque for the wind turbine;
TABLE 1 reference table for rotor speed and wind turbine output power
Vw(m/s) 7 8 9 10 11 12
ωr(p.u) 0.75 0.85 0.95 1.05 1.15 1.25
Pw(p.u) 0.32 0.49 0.69 0.95 1.25 1.6
Tw＝Pw/ωr 0.43 0.58 0.73 0.90 1.09 1.28
Step 1.2, writing out an electromagnetic equation of the system in a dq coordinate system as follows:
step 1.2.1, the stator and the rotor in the doubly-fed induction generator all adopt the motor convention, the positive direction is the current flowing direction, and the rotating direction is the same as the positive direction of the electromagnetic torque. According to the coordinate transformation principle, variables such as current, voltage, magnetic flux and the like in the stator and the rotor are transformed from a three-phase static coordinate to a dq coordinate system, and the electromagnetic relation in the stator and the rotor is obtained as follows:
Figure BDA0002184983430000071
Figure BDA0002184983430000072
Figure BDA0002184983430000073
wherein ω is1For synchronizing the angular velocities of rotation, ωs＝ω1-ωrIs the slip angular velocity, omegarIs the rotor rotational angular velocity; u. ofds、uqsIs the stator dq-axis voltage, udr、uqrIs the rotor dq axis voltage; i.e. ids、iqsFor stator dq axis current, idr、iqrIs the rotor dq axis current; psids、ψqsFor stator dq axis flux linkage, psidr、ψqrA rotor dq axis flux linkage; rsIs stator resistance, RrIs the rotor resistance; l islsFor stator leakage inductance, LlrFor rotor leakage inductance, LmFor coaxial mutual inductance of stator and rotor, Ls、LrThe stator and rotor windings are equivalent to self inductance.
Step 1.2.2, controlling a target of a double-fed fan rotor side converter by tracking the maximum wind energy, keeping the frequency of generated energy constant and controlling reactive output; the stator flux linkage direction is fixed on the d axis, so that active and reactive decoupling can be realized; the rotor voltage is regulated using PI control, and the control block diagram of the rotor-side converter is shown in fig. 3. Wherein Teref、QsrefReference value for electromagnetic torque and stator reactive power, KTe、KQsProportional coefficient, K, for the external power loop PI control of the rotor-side converteriq、KidProportional coefficient for PI control of current loop in the rotor side converter; t isTe、TQsIntegral coefficient, T, for rotor side converter outer power loop PI controliq、TidIntegral coefficient of current loop PI control in the rotor side converter; i.e. iqrref、idrrefAnd a rotor dq axis current reference value output by the rotor side converter outer ring PI control.
The relation column expressed by the control block diagram is written into a differential equation as shown in formulas (4) to (11):
Figure BDA0002184983430000074
Figure BDA0002184983430000075
Figure BDA0002184983430000076
Figure BDA0002184983430000077
Figure BDA0002184983430000081
Figure BDA0002184983430000082
Figure BDA0002184983430000083
Figure BDA0002184983430000084
wherein x is1～x4Four intermediate variables were introduced.
Step 2, applying the analog inductance control strategy to the inner loop control of the rotor side converter of the doubly-fed wind turbine to control the output voltage of the rotor, and combining with the graph 4, the method specifically comprises the following steps:
step 2.1, detecting external wind speed, finding a torque reference value corresponding to the wind speed by referring to the table 1, determining a reactive reference value of a stator according to a specific control target, establishing a relevant state equation of a simultaneous system, and calculating a steady-state expected value of a current dq axis component of a rotor of the doubly-fed wind turbine by combining a corresponding control strategy
Figure BDA0002184983430000085
Step 2.2, taking the difference between the actual value of the rotor current dq axis component and the steady-state expected value to obtain a gain, and then using the gain as an additional signal to act on the input of the rotor voltage, wherein the formula is as follows:
Figure BDA0002184983430000086
and 2.3, enabling the rotor voltage output by the rotor side converter to meet the relation of an equation (13) under the action of the additional signal of the equation (12):
Figure BDA0002184983430000087
wherein G isPI.e. the analog inductor, the value principle of which is introduced in step 3.
Step 3, carrying out self-adaptive online setting on the simulated inductance value by utilizing a PSO algorithm, and determining the optimal value of the simulated inductance, wherein the method specifically comprises the following steps:
step 3.1, writing a system state equation in a row, linearizing the system state equation at a balance point to obtain a small signal equation:
dΔX/dt＝AΔX (14)
in the formula, X is the state variable of the system, and the number of X is 20; a is a Jacobian matrix determined by system parameters and balance point states;
step 3.2, determining the subsynchronous oscillation frequency f of the system according to the formula (15)nThe power oscillation frequency corresponding thereto is fSSR＝f1-fnBased on the above, determining the corresponding subsynchronous oscillation mode in the eigenvalue of the jacobian matrix a, the formula is:
Figure BDA0002184983430000091
in the formula (f)1For mains frequency, XCFor series compensation of capacitive reactance of capacitor, X∑Is the sum of effective inductive reactance in the system;
step 3.3 at GPWhen the value is changed, two groups of characteristic value real parts are changed, one group corresponds to a subsynchronous mode, the other group corresponds to a low-frequency mode, and the characteristic value corresponding to the subsynchronous mode is used as sigma1±jω1The characteristic value corresponding to the low-frequency mode is represented by σ2±jω2The method is characterized in that a Particle Swarm Optimization (PSO) algorithm is adopted to perform self-adaptive setting on an analog inductance value, the primary goal of the algorithm is to inhibit subsynchronous oscillation, and therefore a damping ratio corresponding to a subsynchronous modal characteristic value is taken as a fitness function, as shown in a formula (16):
Figure BDA0002184983430000092
step 3.4, GPThe improper value selection can introduce low-frequency oscillation, so that the real part of the characteristic value corresponding to the low-frequency mode needs to be limited in a certain range, and one constraint of the algorithm is obtained:
σ2≤σ2max (17)
in the formula sigma2maxA constant less than zero;
and 3.5, updating the position and the search speed of the ith particle according to the formula (18):
Figure BDA0002184983430000093
in the formula, k is the current iteration number; w is the inertial weight; c. C1And c2Is a learning factor; r is1And r2A random number in the range of 0 to 1;
Figure BDA0002184983430000094
is the current (iterated to kth) best position for the ith particle,
Figure BDA0002184983430000095
is the current global optimum position;
and 3.6, in the standard PSO algorithm, the inertia weight w is constant, and the algorithm is easy to fall into a local optimal solution under the condition. Therefore, a dynamic adjustment model of the inertial weight is established, and in the initial stage of the algorithm, the global search capability is enhanced by searching with fixed weight; and after the set number of iterations, varying w according to equation (19):
Figure BDA0002184983430000096
in the formula, wmaxAnd wminIs a fixed value, fiFor the current particle adaptation value, favAverage of the adapted values for all particles, fbestAdapting the value for the global optimum particle;
under the model, the larger the fitness of the particles is, the smaller the corresponding inertial weight is, and the stronger the local searching capability is, so that a new optimal position can be found;
step 3.7, in order to reduce the operation amount, the analog inductor G needs to be simulatedPThe optimizing range of the value is restricted; in order to avoid too large a particle search speed, the particle speed needs to be constrained, as shown in equation (20):
Figure BDA0002184983430000101
wherein G isPminAnd GPmaxTo a fixed value selected as required, vmaxAnd vminThe upper and lower limits of the particle search speed can be selected according to typical values. A simplified block diagram of the rotor side converter control strategy with the addition of the PSO optimization algorithm is shown in fig. 5.
The principle of analog inductance parameter tuning can be summarized as follows: detecting the current running state of the fan system at a fixed time point, and reading in system data; the PSO algorithm program optimizes parameters according to the current system state and adjusts the corresponding GPA value; and detecting the system working condition again after a fixed time period, and repeating the process.
And 4, applying the control strategies determined in the steps II and III of the analog inductance control strategy and the rotor side converter control strategy to the control of the rotor side of the double-fed fan to complete the suppression of the subsynchronous oscillation of the system.
The simulated inductor G obtained by PSO algorithm setting in the third stepPThe values are applied to the block diagram shown in fig. 4, and the control mode of the rotor converter shown in fig. 3 is replaced, so that the analog inductance control of the doubly-fed wind turbine is realized.
The present invention will be described in detail with reference to examples.
Examples
Under a Matlab platform, a system shown in fig. 2 is built, a wind power plant in the system is boosted by a transformer and then is connected to an infinite system through a series compensation line, wherein the 100MW wind power plant is formed by combining 50 double-fed wind power generators of 2MW, and specific system parameters are shown in tables 2 to 4:
TABLE 2 Induction Generator parameters
Figure BDA0002184983430000102
Figure BDA0002184983430000111
TABLE 3 Transmission line and shafting parameters
Name (R) Parameter(s) Name (R) Parameter(s)
Transformation ratio of transformer 690V/161kV Inertia of fan 4.29s
Reference capacity 100MVA Inertia of generator 0.90s
Line resistor 0.02pu Generator damping 0.00pu
Line reactance 0.50pu Shafting damping system 0.00pu
Reactance of transformer 0.14pu Wind turbine damping 1.50pu
System for controlling a power supplyImpedance XS 0.06pu Shafting stiffness system 0.15pu/rad
TABLE 4 RSC controller parameters
Name (R) Parameter(s) Name (R) Parameter(s)
TTe 0.05 Kiq 0.0001
TQs 0.025 Kid 0.0001
Tiq 0.005 KTe 0.0001
Tid 0.0025 KQs 0.0001
Fig. 6 is a simulation result of the conventional PI control strategy, and the operation conditions are as follows: the wind speed is 7m/s, the initial series compensation degree of the line is 20%, and the line is in steady-state operation. When the series compensation degrees are changed to 40%, 62.5% and 65% respectively at the time t being 0.5s, an oscillation curve of the fan electromagnetic torque Te is obtained. According to the simulation result, when the series compensation degree is 40%, the oscillation caused by disturbance can be automatically subsided; when the series compensation degree is 62.5%, the electromagnetic torque becomes constant-amplitude oscillation, the stability of the system is poor, and the system is in a critical stable state; when the series compensation degree is further increased to 65%, the stability of the system continues to be deteriorated, the subsynchronous oscillation component is rapidly dispersed, and the system is unstable. The above results indicate that the larger the degree of crosstalk, the worse the stability of the system, and the more likely subsynchronous oscillation occurs.
The method for restraining the doubly-fed wind power plant subsynchronous oscillation based on the simulated inductance is adopted for time domain simulation, and the algorithm parameters are set as follows: the particle swarm is 200; the number of iterations is 3000; inertial weight wmax＝1.2，wmin0.4; learning factor c1＝c22; particle search Range [0,50 ]](ii) a The particle search speed range is [ -5,5 [)]. Selecting three typical working conditions, and obtaining the optimal G of the three typical working conditions by using a PSO algorithm respectivelyPThe values and time domain simulation results are shown in fig. 7, wherein (a) is an oscillation curve of the electromagnetic torque when the series compensation degree is 70% at a wind speed of 7 m/s; (b) the system is still stable due to the oscillation curve of the electromagnetic torque at the wind speed of 7m/s and the series compensation degree of 80 percent; (c) the wind speed is 8m/s, and the electromagnetic torque oscillation curve when the series compensation degree is 80 percent, and the figure shows that the method for inhibiting the sub-synchronous oscillation of the doubly-fed wind power plant based on the simulated inductance can effectively inhibit the sub-synchronous oscillation under different working conditions, and can not bring other adverse effects to the stability of the system.
According to the analysis, the method for inhibiting the sub-synchronous oscillation of the doubly-fed wind power plant based on the simulated inductance has good self-adaptability when the wind speed and the series compensation degree change in a large range, and can effectively inhibit the sub-synchronous oscillation phenomenon caused by the series compensation grid connection of the doubly-fed wind power plant.

Claims (5)

1. A doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance is characterized by comprising the following steps:
step 1, establishing a system model of a doubly-fed wind power plant which is connected with the grid through a series compensation capacitor;
step 2, applying an analog inductance control strategy to the inner ring control of the rotor side converter of the double-fed fan to control the output voltage of the rotor;
step 3, carrying out self-adaptive on-line setting on the simulated inductance value by utilizing a PSO algorithm, and determining the optimal value of the simulated inductance under a specific working condition;
and 4, applying an analog inductance control strategy and a rotor side converter control strategy to the control of the rotor side of the doubly-fed wind turbine to complete the suppression of the subsynchronous oscillation of the system.
2. The method for suppressing the subsynchronous oscillation of the doubly-fed wind farm based on the analog inductance according to claim 1, wherein the step 1 of establishing the system model of the doubly-fed wind farm through the series compensation capacitance grid connection specifically comprises the following steps:
step 1.1, obtaining the corresponding relation among the rotating speed, the output power and the wind speed of the rotor of the doubly-fed wind turbine according to a wind energy capture model;
step 1.2, writing out an electromagnetic equation of the system in a dq coordinate system as follows:
step 1.2.1, the stator and the rotor in the doubly-fed induction generator adopt the motor convention, the positive direction is the current flowing direction, and the rotating direction is the same as the positive direction of the electromagnetic torque; according to the coordinate transformation principle, current, voltage and magnetic flux variables in the stator and the rotor are transformed from a three-phase static coordinate to a dq coordinate system, and the electromagnetic relation in the stator and the rotor is obtained as follows:
Figure FDA0002184983420000011
Figure FDA0002184983420000012
Figure FDA0002184983420000013
wherein ω is1For synchronizing the angular velocities of rotation, ωs＝ω1-ωrIs the slip angular velocity, omegarIs the rotor rotational angular velocity; u. ofds、uqsIs the stator dq-axis voltage, udr、uqrIs the rotor dq axis voltage; i.e. ids、iqsFor stator dq axis current, idr、iqrIs the rotor dq axis current; psids、ψqsFor stator dq axis flux linkage, psidr、ψqrA rotor dq axis flux linkage; rsIs stator resistance, RrIs the rotor resistance; l islsFor stator leakage inductance, LlrFor rotor leakage inductance, LmFor coaxial mutual inductance of stator and rotor, Ls、LrEquivalent self-inductance of stator and rotor windings;
step 1.2.2, controlling a target of a double-fed fan rotor side converter by tracking the maximum wind energy, keeping the frequency of generated energy constant and controlling reactive output; fixing the stator flux linkage direction on a d axis to realize active and reactive decoupling; and (3) regulating the rotor voltage by adopting PI control to obtain a differential equation:
Figure FDA0002184983420000021
Figure FDA0002184983420000022
Figure FDA0002184983420000023
Figure FDA0002184983420000024
Figure FDA0002184983420000025
Figure FDA0002184983420000026
Figure FDA0002184983420000027
Figure FDA0002184983420000028
wherein Teref、QsrefReference value for electromagnetic torque and stator reactive power, KTe、KQsProportional coefficient, K, for the external power loop PI control of the rotor-side converteriq、KidProportional coefficient for PI control of current loop in the rotor side converter; t isTe、TQsIntegral coefficient, T, for rotor side converter outer power loop PI controliq、TidIntegral coefficient of current loop PI control in the rotor side converter; i.e. iqrref、idrrefRotor dq axis current reference, x, for rotor side converter outer loop PI control output1～x4Four intermediate variables were introduced.
3. The method for suppressing the subsynchronous oscillation of the doubly-fed wind power plant based on the analog inductance according to claim 1, wherein the analog inductance control strategy in the step 2 is applied to the inner loop control of the rotor-side converter of the doubly-fed wind turbine to control the output voltage of the rotor, and specifically comprises the following steps:
step 2.1, detecting external wind speed to obtain a torque reference value corresponding to the wind speed, determining a reactive reference value of a stator according to a control target, establishing a relevant state equation of a simultaneous system, and calculating a steady-state expected value of a current dq axis component of a rotor of the doubly-fed fan by combining a corresponding control strategy
Figure FDA0002184983420000031
Step 2.2, after the actual value of the rotor current dq axis component is subjected to difference gain with the steady-state expected value, the difference gain is used as an additional signal to act on the input of the rotor voltage, and the formula is as follows:
Figure FDA0002184983420000032
and 2.3, enabling the rotor voltage output by the rotor side converter to meet the relation of an equation (13) under the action of the additional signal of the equation (12):
Figure FDA0002184983420000033
wherein G isPTo simulate an inductance.
4. The method for suppressing the subsynchronous oscillation of the doubly-fed wind power plant based on the simulated inductance according to claim 1, wherein the step 3 is to perform adaptive online setting on the simulated inductance value by using a PSO algorithm to determine the optimal value of the simulated inductance under a specific working condition, and specifically comprises the following steps:
step 3.1, writing a system state equation in a row, linearizing the system state equation at a balance point to obtain a small signal equation:
dΔX/dt＝AΔX (14)
in the formula, X is the state variable of the system, and the number of X is 20; a is a Jacobian matrix determined by system parameters and balance point states;
step 3.2,Calculating the subsynchronous oscillation frequency f of the systemnThe formula is as follows:
Figure FDA0002184983420000034
in the formula (f)1For mains frequency, XCFor series compensation of capacitive reactance of capacitor, X∑Is the sum of effective inductive reactance in the system;
corresponding to a power oscillation frequency of fSSR＝f1-fnBased on the characteristic value, determining a corresponding subsynchronous oscillation mode in the characteristic value of the Jacobian matrix A;
step 3.3 at GPWhen the value is changed, two groups of characteristic value real parts are changed, wherein one group corresponds to a subsynchronous mode, and the corresponding characteristic value is sigma1±jω1(ii) a The other group corresponds to low-frequency mode, and the corresponding characteristic value is sigma2±jω2The PSO algorithm is adopted to carry out self-adaptive setting on the analog inductance value, the primary objective of the algorithm is to inhibit subsynchronous oscillation, and therefore the damping ratio corresponding to the subsynchronous modal characteristic value is taken as a fitness function, as shown in the formula (16):
Figure FDA0002184983420000041
step 3.4, limiting the real part of the characteristic value corresponding to the low-frequency mode in a certain range to obtain a constraint of the algorithm:
σ2≤σ2max (17)
in the formula sigma2maxA constant less than zero;
and 3.5, updating the position and the search speed of the ith particle according to the formula (18):
Figure FDA0002184983420000042
in the formula, k is the current iteration number; w is the inertial weight; c. C1And c2Is a learning factor; r is1And r2A random number in the range of 0 to 1;
Figure FDA0002184983420000043
is the current best position of the ith particle,
Figure FDA0002184983420000044
is the current global optimum position;
3.6, establishing a dynamic adjustment model of the inertial weight, and searching by using the fixed weight at the initial stage of the algorithm; and after the set number of iterations, varying w according to equation (19):
Figure FDA0002184983420000045
in the formula, wmaxAnd wminIs a fixed value, fiFor the current particle adaptation value, favAverage of the adapted values for all particles, fbestAdapting the value for the global optimum particle;
step 3.7, simulating the inductor GPThe optimum range of values and particle velocity are constrained:
Figure FDA0002184983420000046
wherein G isPminAnd GPmaxTo a fixed value selected as required, vmaxAnd vminThe upper and lower limits of the particle search speed can be selected according to typical values.
5. The method for suppressing the subsynchronous oscillation of the doubly-fed wind power plant based on the analog inductance of claim 1, wherein the analog inductance control strategy and the rotor-side converter control strategy are applied to the control of the rotor side of the doubly-fed wind turbine in the step 4, so as to suppress the subsynchronous oscillation of the system, specifically as follows: analog inductor obtained by PSO algorithm settingGPThe value is applied to the control of the converter, the control of the analog inductance of the double-fed fan is realized, and the suppression of the subsynchronous oscillation of the system is completed.
CN201910810798.9A 2019-08-29 2019-08-29 Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance Active CN112448399B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910810798.9A CN112448399B (en) 2019-08-29 2019-08-29 Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910810798.9A CN112448399B (en) 2019-08-29 2019-08-29 Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance

Publications (2)

Publication Number Publication Date
CN112448399A true CN112448399A (en) 2021-03-05
CN112448399B CN112448399B (en) 2022-09-16

Family

ID=74742196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910810798.9A Active CN112448399B (en) 2019-08-29 2019-08-29 Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance

Country Status (1)

Country Link
CN (1) CN112448399B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112952863A (en) * 2021-04-19 2021-06-11 华北电力大学 Doubly-fed system switching type oscillation analysis method based on phase diagram
CN113193589A (en) * 2021-03-31 2021-07-30 哈尔滨工业大学 DFIG wind power plant subsynchronous oscillation suppression method based on digital twin simulation
CN113820571A (en) * 2021-09-03 2021-12-21 国网辽宁省电力有限公司沈阳供电公司 Wind power plant cable insulation online monitoring method and device

Citations (3)

* Cited by examiner, † Cited by third party
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
CN108631338A (en) * 2018-06-21 2018-10-09 国网宁夏电力有限公司电力科学研究院 A method of for inhibiting the grid-connected sub-synchronous oscillation in double-fed fan motor field
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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
CN108631338A (en) * 2018-06-21 2018-10-09 国网宁夏电力有限公司电力科学研究院 A method of for inhibiting the grid-connected sub-synchronous oscillation in double-fed fan motor field
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

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113193589A (en) * 2021-03-31 2021-07-30 哈尔滨工业大学 DFIG wind power plant subsynchronous oscillation suppression method based on digital twin simulation
CN112952863A (en) * 2021-04-19 2021-06-11 华北电力大学 Doubly-fed system switching type oscillation analysis method based on phase diagram
CN113820571A (en) * 2021-09-03 2021-12-21 国网辽宁省电力有限公司沈阳供电公司 Wind power plant cable insulation online monitoring method and device
CN113820571B (en) * 2021-09-03 2024-01-26 国网辽宁省电力有限公司沈阳供电公司 Wind power plant cable insulation on-line monitoring method and device

Also Published As

Publication number Publication date
CN112448399B (en) 2022-09-16

Similar Documents

Publication Publication Date Title
CN112448399B (en) Doubly-fed wind power plant subsynchronous oscillation suppression method based on analog inductance
CN109347097B (en) Doubly-fed wind power system subsynchronous oscillation suppression strategy based on improved particle swarm optimization algorithm
CN112448398B (en) Stator side analog resistance-based doubly-fed wind power plant subsynchronous oscillation suppression method
CN109861251B (en) Double-fed fan comprehensive control method for micro-grid transient-steady-state frequency optimization
Mohammadpour et al. SSR analysis of a DFIG-based wind farm interfaced with a gate-controlled series capacitor
CN106786677B (en) A kind of interconnected electric power system distributed dynamic matrix frequency control method
CN108011364B (en) A method of analysis DFIG kinetic characteristics and Electrical Power System Dynamic reciprocal effect
CN106059422B (en) A kind of fuzzy control method inhibited for double-fed fan motor play synchronized oscillation
Pulgar Painemal Wind farm model for power system stability analysis
CN110350587B (en) Doubly-fed fan subsynchronous oscillation suppression method and system based on converter control
CN108270241A (en) The control method of wind turbine gird-connected inverter virtual synchronous generator
CN106786673B (en) The suppressing method and device of double-fed blower compensated transmission system subsynchronous resonance
CN115313419B (en) Direct-drive wind power system control method and system based on network-building type converter
CN105762816B (en) A kind of method for suppressing subsynchronous interaction
CN109787274B (en) Virtual synchronous control method and rotor side frequency converter controller
CN108539760A (en) A kind of double-fed induction Wind turbines frequency modulation PID control method based on group's grey wolf optimization algorithm
CN116094035A (en) New energy unit inertia simulation control method based on super capacitor energy storage
CN106294959B (en) The modeling and simulating method that model reference self-adapting control is combined with double-fed blower reduced-order model
CN111478365B (en) Optimization method and system for control parameters of virtual synchronizer of direct-drive wind turbine generator
CN110417047B (en) Method for analyzing SSCI damping characteristics of doubly-fed fan based on complex torque coefficient
CN109755968A (en) A kind of neural network guaranteed cost virtual synchronous control method of double-fed fan motor unit
CN114759575A (en) Virtual synchronous double-fed fan subsynchronous oscillation suppression method and system
Yang et al. Frequency control scheme with dynamic droop characteristics of a DFIG for mitigating the frequency fluctuations
CN116488211A (en) VSG improved parameter self-adaption method for single-phase photovoltaic energy storage
CN115455687A (en) Wind power plant dynamic aggregation modeling method based on virtual synchronous wind generating set

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Baohua

Inventor after: Shan Bihan

Inventor after: Jiang Haifeng

Inventor after: Liu Yang

Inventor before: Shan Bihan

Inventor before: Wang Baohua

Inventor before: Jiang Haifeng

Inventor before: Liu Yang

GR01 Patent grant
GR01 Patent grant