CN112952861B - Additional virtual double-PSS control method for active support type new energy unit - Google Patents

Additional virtual double-PSS control method for active support type new energy unit Download PDF

Info

Publication number
CN112952861B
CN112952861B CN202110331007.1A CN202110331007A CN112952861B CN 112952861 B CN112952861 B CN 112952861B CN 202110331007 A CN202110331007 A CN 202110331007A CN 112952861 B CN112952861 B CN 112952861B
Authority
CN
China
Prior art keywords
pss
virtual
active support
support type
angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN202110331007.1A
Other languages
Chinese (zh)
Other versions
CN112952861A (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.)
State Grid Liaoning Electric Power Co Ltd
Northeast Electric Power University
Original Assignee
Northeast Dianli University
State Grid Liaoning Electric Power Co Ltd
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 Northeast Dianli University, State Grid Liaoning Electric Power Co Ltd filed Critical Northeast Dianli University
Priority to CN202110331007.1A priority Critical patent/CN112952861B/en
Publication of CN112952861A publication Critical patent/CN112952861A/en
Application granted granted Critical
Publication of CN112952861B publication Critical patent/CN112952861B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses an additional virtual double PSS control method for an active support type new energy unit, and particularly relates to an additional virtual double PSS control method for an active support type new energy unit, which comprises the following steps: the method comprises the steps of establishing an active support type VSC infinite system linearization model, analyzing an oscillation mechanism of the active support type new energy VSC, designing an active support type VSC virtual double PSS, introducing a virtual double PSS link into a control loop on the basis of active support control of a new energy power grid, compensating a damping torque of a virtual excitation regulator, improving damping characteristics, effectively inhibiting power oscillation, improving system stability, and enabling the double PSS to exert the maximum advantages of the double PSS through parameter optimization of an excitation PSS and an angle PSS.

Description

Additional virtual double-PSS control method for active support type new energy unit
Technical Field
The invention belongs to the technical field of electromechanical oscillation of a new-generation energy power system, and particularly relates to an additional virtual double-PSS control method for an active support type new energy unit.
Background
At present, the overuse of fossil energy aggravates environmental pollution, causes a series of problems such as greenhouse effect, and the like, so that the existing energy structure is urgently needed to be improved, and clean energy is vigorously developed to replace the fossil energy. The centralized access of large-scale new energy stations to the power grid is a typical characteristic of a new generation of power system. However, the high-proportion new energy grid-connected operation replaces the traditional unit which provides inertia support by a large rotating shaft, so that the rigidity inertia of the whole system is reduced, and the safe operation and the reliable supply of electric energy of the power system are influenced. Therefore, the active support type new energy unit can simulate the output characteristic of the traditional generator, provide flexible inertia for the system and improve the disturbance resistance of the system. Since the active-support-type new energy source unit simulates the electromechanical transient characteristics of the traditional synchronous generator, the power oscillation problem of the traditional synchronous generator is also introduced into the active-support-type new energy source unit.
Disclosure of Invention
The invention aims to provide an additional virtual double PSS control method for an active support type new energy unit, and solves the problem of power oscillation generated by grid connection of the active support type new energy unit in a new generation of power system.
The technical scheme adopted by the invention is that the control method of the additional virtual double PSS facing the active support type new energy unit is implemented according to the following steps:
step 1, establishing an active support type VSC infinite system linearization model;
step 2, designing an active support type VSC virtual double PSS;
and 3, optimizing the parameters of the active support type VSC virtual double PSS.
The invention is also characterized in that:
the specific process of the step 1 is as follows:
step 1.1, constructing an active support type VSC electromechanical transient model;
and 1.2, deducing a Heffron-Phillips model of the active support type VSC.
The step 1.1 comprises the following specific processes:
according to a second-order rotor motion equation, a first-order excitation loop equation and an excitation voltage regulation equation of the synchronous generator, constructing an active support type VSC electromechanical transient model:
Figure BDA0002996155300000021
in formula (1): p is m Is mechanical power, P e Is electromagnetic power, H is virtual inertia, D is damping coefficient, delta is power angle of the generator, delta omega is deviation between rated rotating speed and actual rotating speed, E q ' is a transient electromotive force, E q Is no-load electromotive force, E fd To force no-load electromotive force, E fd ' is the exciting electromotive force, deltaU, output from the automatic voltage regulator t As a deviation value of the inverter terminal voltage, K A For the gain of the automatic voltage regulator, T A Time constant, T, of the automatic voltage regulator d0 ' is the time constant of the field winding of the synchronous generator, omega 0 The angular frequency is a rated angular frequency, d omega/dt is the change rate of the angular frequency, and d delta/dt is the change rate of a power angle;
in the formula (1), the electromagnetic power P e And no-load electromotive force E q The expression is as follows:
Figure BDA0002996155300000031
in formula (2): u shape d 、U q Is infinite bus voltage U b D-axis component and q-axis component of (I) d 、I q For line current I t D-axis component and q-axis component, X' d For direct-axis transient reactance, X d Is a direct axis synchronous reactance.
The specific process of the step 1.2 is as follows:
the virtual impedance model for constructing the active support type VSC is as follows:
Figure BDA0002996155300000032
u in formula (3) td 、U tq For voltage U at VSC t D-axis component and q-axis component of (a);
the voltage equation of the power transmission line of the active support type VSC infinite system is as follows:
U t =jX t I t +U b (4)
in the formula (4), X t Is the line reactance, I t For line current, U b Is infinite bus voltage, U t Is the inverter terminal voltage;
obtained from equations (1) to (4):
Figure BDA0002996155300000033
linearized by equation (5) at the steady state operating point and simplified, resulting in:
Figure BDA0002996155300000041
wherein the coefficient K 1 ~K 6 Is a constant related to system structure, parameters, operating conditions;
substituting formula (6) into formula (1) to obtain a Heffron-Phillips model of the active support type VSC:
Figure BDA0002996155300000042
the active support type VSC virtual double PSS design in the step 2 comprises a virtual excitation PSS and a virtual angle PSS of electromagnetic torque of a synchronous motor, and the specific process is as follows:
constructing a virtual PSS according to a Heffron-Phillips model of the active support VSC;
the virtual excitation PSS design process comprises the following steps:
taking delta omega as a control signal, enabling the control signal to pass through a transfer function G' pss (s) post-feedback to the excitation circuit to produce additional electromagnetic torque Δ T' pss
Additional electromagnetic torque Delta T 'provided by virtual excitation PSS' pss Comprises the following steps:
Figure BDA0002996155300000043
the virtual angle PSS design process is as follows:
correction quantity delta introduced to power angle delta pss As additional control, where Δ δ pss =G″ pss (s). Δ ω such that the corrected work angle Δ δ = Δ δ pss +Δδ 0 An additional electromagnetic torque Δ T ″' is generated in the excitation circuit pss Additional electromagnetic torque Δ T ″, of virtual angle PSS pss Comprises the following steps:
Figure BDA0002996155300000051
the specific process of the step 3 is as follows: optimizing virtual excitation PSS parameters and virtual angle PSS parameters in the active support type VSC virtual double PSS;
optimizing virtual excitation PSS parameters:
let the phase-corrected value of the control signal Δ ω be Δ u pss I.e. Δ u pss =Δω*G′ pss (s),Δu pss To Δ T' pss Has a transfer function of F' pss (s), the formula is as follows:
Figure BDA0002996155300000052
let damping coefficient provided by virtual excitation PSS be D' pss I.e. by
D′ pss =G′ pss (s)F′ pss (s) (14)
Is prepared from F' pss (s)、G′ pss (s) is expressed in the form:
Figure BDA0002996155300000053
in the formula (I), the compound is shown in the specification,
Figure BDA0002996155300000054
is a transfer function F' pss The phase angle, γ ', of(s) is the transfer function G' pss (s) phase angle;
according to the phase compensation method, there are
Figure BDA0002996155300000055
In formula (II) T' pssd Additional damping torque, T ', provided for virtual excitation PSS' psss An additional synchronous torque provided for the virtual excitation PSS;
namely that
Figure BDA0002996155300000056
Setting a transfer function G' pss (s) is a proportional and lead-lag element, expressed as:
Figure BDA0002996155300000061
in the formula, T 1 ~T 4 All represent the time constant, K ', of the lead-lag link' 1 ,K′ 2 Setting a parameter of the virtual excitation PSS;
assuming that the virtual excitation PSS provides a positive damping torque, available from (16):
Figure BDA0002996155300000062
given T 1 、T 3 Obtaining T from the value of (19) 2 ,T 4 ,K′ 1 ,K′ 2 Obtaining a parametrically optimized transfer function G 'from equations (15) and (18)' pss (s);
And (3) optimizing parameters of the virtual angle PSS:
assuming that the virtual angle PSS provides a positive pure damping torque, Δ T ″ pss Correction quantity delta of power angle delta for electromagnetic torque provided by virtual angle PSS pss For outputting stable control signals of the virtual angle PSS, the electromagnetic torque provided by the virtual angle PSS is delta T ″ pss ,F″ pss (s) is from Δ δ pss To Δ T ″) pss The transfer function of (2) is then:
Figure BDA0002996155300000063
ΔT″ pss =F″ pss (s)Δδ pss (21)
order transfer function G ″) pss (s) phase compensating transfer function F ″) pss (s) phase, setting the transfer function G ″) pss (s) is a proportional and lead-lag element, expressed as:
Figure BDA0002996155300000064
in formula (II) T' 1 ~T′ 4 All represent the time constant, K ″, of the lead-lag link 1 ,K″ 2 Setting a parameter of a virtual angle PSS;
by the same virtual excitation PSS, the method can obtain
Figure BDA0002996155300000071
In the formula (I), the compound is shown in the specification,
Figure BDA0002996155300000072
is a transfer function F ″ pss The phase angle, γ ", of(s) is the transfer function G ″ pss (s) phase angle;
given T' 1 、T′ 3 Is obtained from the value of (23)' 2 ,T′ 4 ,K′ 1 ,K′ 2 The parameter-optimized transfer function G' is obtained from the formula (22) pss (s)。
The invention has the beneficial effects that:
the invention provides an additional virtual double PSS control method for an active support type new energy unit, and the additional virtual double PSS control method can provide an additional damping torque through a virtual double PSS control strategy when a system is interfered to generate power oscillation, so that the power oscillation is inhibited, and the disturbance resistance of a new generation energy power system is improved.
Drawings
FIG. 1 is a block diagram of an electromechanical transient model of an actively-supported new energy unit according to the present invention;
FIG. 2 is a diagram of an active support VSC infinite system of the present invention;
FIG. 3 is a Heffron-Phillips model of the active support VSC of the present invention;
FIG. 4 shows the range from Δ δ to Δ T of the present invention e A block diagram of the transfer function of (1);
FIG. 5 is an infinite electromagnetic torque vector diagram of the active support VSC of the present invention;
fig. 6 is a damping control block diagram of the additional virtual excitation PSS of the present invention;
FIG. 7 shows the equation Δ u according to the present invention pss To Δ T' pss A block diagram of the transfer function of (1);
fig. 8 is an electromagnetic torque vector diagram of the present invention incorporating a virtual excitation PSS;
FIG. 9 is a block diagram of damping control of the additional virtual angle PSS according to the present invention;
FIG. 10 is a plot of delta of the present invention pss To Δ T ″) pss A block diagram of the transfer function of (1);
FIG. 11 is an electromagnetic torque vector diagram incorporating the virtual angle PSS of the present invention;
FIG. 12 is a damping control block diagram of the additional dual PSS of the present invention;
FIG. 13 is a frequency response waveform of an actively-supported VSC infinite system of the present invention;
fig. 14 is a power response waveform of the actively-supported VSC infinite system of the present invention.
Detailed Description
The invention is described in detail below with reference to the drawings and the detailed description.
The invention relates to an additional virtual double PSS control method for an active support type new energy unit, which is implemented according to the following steps:
step 1, establishing an active support type VSC infinite system linearization model for analyzing the power oscillation problem of a power system; the specific process of the step 1 is as follows:
step 1.1, obtaining an active support type VSC electromechanical transient model according to a second-order rotor motion equation, a first-order excitation loop equation and an excitation voltage regulation equation of the synchronous generator, wherein a control block diagram of the model is shown in FIG. 1, and the formula is as follows:
Figure BDA0002996155300000081
in formula (1): p m As mechanical power, P e Is electromagnetic power, H is virtual inertia, D is damping coefficient, delta is power angle of the generator, delta omega is deviation between rated rotating speed and actual rotating speed, E q Is a transient electromotive force, E q Is no-load electromotive force, E fd To force no-load electromotive force, E fd ' is the exciting electromotive force, deltaU, output from the automatic voltage regulator t As a deviation value of the inverter terminal voltage, K A For the gain of the automatic voltage regulator, T A Is the time constant of the automatic voltage regulator, T d0 ' is the time constant, omega, of the field winding of the synchronous generator 0 The angular frequency is a rated angular frequency, d omega/dt is the change rate of the angular frequency, and d delta/dt is the change rate of a power angle;
in the formula (1), the electromagnetic power P e And no-load electromotive force E q The expression is as follows:
Figure BDA0002996155300000091
in the formula (2): u shape d 、U q Is infinite bus voltage U b D-axis component and q-axis component of (I) d 、I q For line current I t D-and q-axis components, X' d For direct-axis transient reactance, X d Is a direct axis synchronous reactance.
Step 1.2, constructing a virtual impedance model of the active support type VSC:
Figure BDA0002996155300000092
u in formula (3) td 、U tq For the voltage U at the VSC t D-axis component and q-axis component of (a);
the active support type VSC infinite system is shown in fig. 2, and the voltage equation of the transmission line is:
U t =jX t I t +U b (4)
in the formula (4), X t Is line reactance, I t For line current, U b Is infinite bus voltage, U t Is the inverter terminal voltage;
obtained from equations (1) to (4):
Figure BDA0002996155300000093
linearized by equation (5) at the steady state operating point and simplified, resulting in:
Figure BDA0002996155300000101
wherein the coefficient K 1 ~K 6 Is a constant related to system structure, parameters, operating conditions;
substituting formula (6) into formula (1) to obtain a Heffron-Phillips model of the active support VSC as shown in FIG. 3:
Figure BDA0002996155300000102
the state space expression of the active support type VSC infinite system is as follows:
Figure BDA0002996155300000103
step 2, designing an active support type VSC virtual double PSS in order to reduce the power oscillation problem caused by negative damping torque and improve the system robustness; designing the VSC virtual double PSS of the active support type comprises designing a virtual excitation PSS and a virtual angle PSS of electromagnetic torque of the synchronous motor,
FIG. 4 shows the range from Δ δ to Δ T e As can be found from fig. 4, the electromagnetic torque provided by the active support VSC is
Figure BDA0002996155300000111
An active support type VSC electromagnetic torque vector diagram can be drawn according to equation (9), as shown in fig. 5. Due to K in FIG. 4 4 The branch does not generate negative damping torque, and for the sake of simplicity of analysis, K is not considered in FIG. 5 4 The influence of the branch. It can be seen that under the environment of a quick-response and high-amplification-factor quick excitation system, the resultant torque delta T is e In the fourth quadrant of the vector diagram, damping torque Δ T d The negative number easily causes power oscillation problem and reduces system stability.
Therefore, the present invention provides an additional damping torque by adding the virtual excitation PSS and the virtual angle PSS, suppressing the power oscillation.
The specific process is as follows:
constructing a virtual PSS according to a Heffron-Phillips model of the active support type VSC;
the virtual excitation PSS design process comprises the following steps:
as shown in fig. 6, the control block diagram of the additional virtual excitation PSS is such that Δ ω is used as a control signal, and the control signal is subjected to a transfer function G' pss (s) post-feedback to the excitation loop to produce additional electromagnetic torque Δ T' pss
As can be seen from fig. 7, the additional electromagnetic torque Δ T 'provided by the virtual excitation PSS' pss Comprises the following steps:
Figure BDA0002996155300000112
electromagnetic torque DeltaT 'can be plotted according to equation (10)' pss Vector diagram, for simplicity of analysis, temporarily set G' pss (s) is a constant coefficient, and G 'is heavily considered in the parameter optimization link' pss Parameter(s), electromagnetic torque Δ T' pss The vector diagram is shown in fig. 8.
As can be seen from fig. 8, additional electromagnetic torque Δ T 'provided by virtual excitation PSS' pss In the first quadrant, a positive damping torque Δ T 'can be provided' d And the original negative damping torque of the system is reduced.
The virtual angle PSS design process is as follows:
when G 'due to limited damping torque provided by virtual excitation PSS' pss Too large a value of(s) may affect the system stability. Therefore, the angle virtual PSS is designed, so that the angle virtual PSS can provide additional damping torque on the basis of the excitation PSS.
A control block diagram of the additional virtual angle PSS is shown in fig. 9. Introducing correction quantity delta of power angle delta in system overall control framework pss As additional control, where Δ δ pss =G″ pss (s). Δ ω such that the corrected power angle Δ δ = Δ δ pss +Δδ 0 An additional electromagnetic torque Δ T ″' is generated in the excitation circuit pss As shown in FIG. 10, the additional electromagnetic torque Δ T ″' of the virtual angle PSS pss Comprises the following steps:
Figure BDA0002996155300000121
the electromagnetic torque Δ T ″, which can be plotted according to equation (11) pss Vector diagram, temporarily setting G' for simplifying analysis pss (s) is a constant coefficient, neglecting K 4 The influence of branch is considered in the parameter optimization link pss (s) parameter, electromagnetic Torque Δ T ″) pss The vector diagram is shown in fig. 11.
As can be seen from FIG. 11, the additional electromagnetic torque Δ T ″, provided by the virtual angle PSS pss In the first quadrant of the light source,can provide a positive damping torque DeltaT d And the original negative damping torque of the system is reduced.
In summary, the control block diagram of the virtual double PSS is shown in fig. 12, and the total electromagnetic torque Δ T provided by the virtual double PSS e * Is composed of
ΔT e * =ΔT e +ΔT pss +ΔT″ pss (12)
Through the additional control of the virtual double PSS, the damping torque of the system can be greatly improved, and the stability of the system is improved.
And 3, in order to improve the oscillation mode of the system, improve the controllability of the virtual double PSS and optimize the parameters of the active support type VSC virtual double PSS. The specific process of the step 3 is as follows:
optimizing virtual excitation PSS parameters and virtual angle PSS parameters in the active support type VSC virtual double PSS;
note that for simplicity of analysis, the transfer function G 'is temporarily set' pss (s)、G″ pss (s) is a constant coefficient, and in order for the virtual dual PSS to provide a positive pure damping torque, the virtual dual PSS parameters are optimized below.
Optimizing virtual excitation PSS parameters:
let the value after phase correction of the control signal Δ ω be Δ u pss Transfer function G 'is required to enable the virtual excitation PSS to provide a positive pure damping torque' pss The phase of(s) can be compensated for Δ u pss To Δ T' pss Of the transfer function, i.e. au pss =Δω*G′ pss (s),Δu pss To delta T' pss Is F' pss (s), the formula is as follows:
Figure BDA0002996155300000131
let damping coefficient provided by virtual excitation PSS be D' pss I.e. by
D′ pss =G′ pss (s)F′ pss (s) (14)
F' pss (s)、G′ pss (s) is expressed in the form:
Figure BDA0002996155300000132
in the formula (I), the compound is shown in the specification,
Figure BDA0002996155300000133
is a transfer function F' pss The phase angle, γ 'of(s) is the transfer function G' pss (s) phase angle;
according to the phase compensation method, there are
Figure BDA0002996155300000134
In formula (II) T' pssd Additional damping torque, T ', provided for virtual excitation PSS' psss An additional synchronous torque provided for the virtual excitation PSS;
namely, it is
Figure BDA0002996155300000135
Setting a transfer function G' pss (s) is a proportional and lead-lag element, expressed as:
Figure BDA0002996155300000136
in the formula, T 1 ~T 4 All represent the time constant, K ', of the lead-lag link' 1 ,K′ 2 Setting a parameter of the virtual excitation PSS;
assuming that the virtual excitation PSS provides a positive damping torque, it is obtained from (16):
Figure BDA0002996155300000141
given T 1 、T 3 The value of (a) is,obtaining T from equation (19) 2 ,T 4 ,K′ 1 ,K′ 2 Obtaining a parametrically optimized transfer function G 'from equations (15) and (18)' pss (s);
And (3) optimizing parameters of the virtual angle PSS:
assuming that the virtual angle PSS provides a positive pure damping torque, Δ T ″ pss Correction quantity delta of power angle delta for electromagnetic torque provided by virtual angle PSS pss For outputting a stable control signal for the virtual angle PSS, the electromagnetic torque provided by the virtual angle PSS is delta T ″ pss ,F″ pss (s) is from Δ δ pss To Δ T ″) pss The transfer function of (2) is then:
Figure BDA0002996155300000142
ΔT″ pss =F″ pss (s)Δδ pss (21)
in order for the virtual angle PSS to provide a positive pure damping torque, it is necessary to let the transfer function G ″' pss (s) phase compensating transfer function F ″) pss (s) phase, setting the transfer function G ″) pss (s) is a proportional and lead-lag element, expressed as:
Figure BDA0002996155300000143
in formula (II) T' 1 ~T′ 4 All represent the time constant, K ″, of the lead-lag link 1 ,K″ 2 Setting a parameter of a virtual angle PSS;
by the same virtual excitation PSS, the method can obtain
Figure BDA0002996155300000151
In the formula (I), the compound is shown in the specification,
Figure BDA0002996155300000152
to transmitFunction F ″) pss The phase angle, γ ", of(s) is the transfer function G ″ pss (s) phase angle;
given T' 1 、T′ 3 T 'is obtained from the value of (23)' 2 ,T′ 4 ,K′ 1 ,K′ 2 The parameter-optimized transfer function G' is obtained from the formula (22) pss (s)。
And finally, a single-machine infinite power grid simulation example is set up in DIgSILENT/PoweFactory simulation software, system power disturbance is set when t =200s, a system frequency response waveform is shown in fig. 10, and a converter output power waveform is shown in fig. 11. It can be obviously seen that the oscillation amplitude of the system frequency and power can be reduced by introducing the virtual excitation PSS, and the damping characteristic of the system can be further improved by introducing the virtual angle PSS.
Through the mode, the additional virtual double PSS control method for the active support type new energy unit is used for compensating the damping torque of the virtual excitation regulator and improving the damping characteristic by adding the virtual double PSS link into the control loop, so that power oscillation can be effectively inhibited, and the system stability is improved. The accuracy and the rationality of the method are verified through simulation.

Claims (1)

1. The control method of the additional virtual double PSS facing the active support type new energy unit is characterized by comprising the following steps:
step 1, establishing an active support type VSC infinite system linearization model; the specific process is as follows:
step 1.1, constructing an active support type VSC electromechanical transient model; the specific process is as follows:
according to a second-order rotor motion equation, a first-order excitation loop equation and an excitation voltage regulation equation of the synchronous generator, constructing an active support type VSC electromechanical transient model:
Figure FDA0003780305060000011
in formula (1): p m As mechanical power, P e Is electromagnetic power, H is virtual inertia, D is damping coefficient, delta is power angle of the generator, delta omega is deviation between rated rotating speed and actual rotating speed, E q ' is a transient electromotive force, E q Is no-load electromotive force, E fd To force no-load electromotive force, E fd ' is the exciting electromotive force, delta U, output by the automatic voltage regulator t Is the value of the offset of the inverter terminal voltage, K A For the gain of the automatic voltage regulator, T A Is the time constant of the automatic voltage regulator, T d0 ' is the time constant, omega, of the field winding of the synchronous generator 0 D omega/dt is the change rate of the angular frequency, and d delta/dt is the change rate of the power angle;
in the formula (1), the electromagnetic power P e And no-load electromotive force E q The expression is as follows:
Figure FDA0003780305060000012
in formula (2): u shape d 、U q Is infinite bus voltage U b D-axis component and q-axis component of (I) d 、I q For line current I t D-axis component and q-axis component, X' d For direct-axis transient reactance, X d Is a direct-axis synchronous reactance;
step 1.2, deducing a Heffron-Phillips model of the active support type VSC; the specific process is as follows:
the virtual impedance model for constructing the active support type VSC is as follows:
Figure FDA0003780305060000021
u in formula (3) td 、U tq For voltage U at VSC t D-axis component and q-axis component of (a);
the voltage equation of the power transmission line of the active support type VSC infinite system is as follows:
U t =jX t I t +U b (4)
in the formula (4), X t Is the line reactance, I t For line current, U b Is infinite bus voltage, U t Is the inverter terminal voltage;
obtained from equations (1) to (4):
Figure FDA0003780305060000022
linearized by equation (5) at the steady state operating point and simplified to give:
Figure FDA0003780305060000023
wherein the coefficient K 1 ~K 6 Is a constant related to system structure, parameters, operating conditions;
bringing the formula (6) into the formula (1), and obtaining a Heffron-Phillips model of the active support type VSC, namely a linearization control model of an infinite system of the active support type VSC:
Figure FDA0003780305060000031
step 2, designing an active support type VSC virtual double PSS;
the design initiative supports the virtual two PSS of type VSC including the virtual excitation PSS and the virtual angle PSS of the electromagnetic torque of design synchronous machine, and the concrete process is:
constructing a virtual PSS according to a Heffron-Phillips model of the active support type VSC;
the virtual excitation PSS design process comprises the following steps:
Δ ω is used as a control signal, and the control signal is passed through a transfer function G' pss (s) post-feedback to the excitation loop to produce additional electromagnetic torque Δ T' pss
Additional electromagnetic torque delta T 'provided by virtual excitation PSS' pss Comprises the following steps:
Figure FDA0003780305060000032
the virtual angle PSS design process is as follows:
correction quantity delta introduced into power angle delta pss As additional control, where Δ δ pss =G″ pss (s). Δ ω such that the corrected power angle Δ δ = Δ δ pss +Δδ 0 An additional electromagnetic torque Δ T ″' is generated in the excitation circuit pss Additional electromagnetic torque Δ T ″' of virtual angle PSS pss Comprises the following steps:
Figure FDA0003780305060000033
step 3, optimizing the parameters of the active support type VSC virtual double PSS; the specific process is as follows: optimizing virtual excitation PSS parameters and virtual angle PSS parameters in the active support type VSC virtual double PSS;
optimizing virtual excitation PSS parameters:
let the phase-corrected value of the control signal Δ ω be Δ u pss I.e. Δ u pss =Δω*G′ pss (s),Δu pss To Δ T' pss Has a transfer function of F' pss (s), the formula is as follows:
Figure FDA0003780305060000041
let damping coefficient provided by virtual excitation PSS be D' pss I.e. by
D′ pss =G′ pss (s)F′ pss (s) (14)
F' pss (s)、G′ pss (s) is expressed in the form:
Figure FDA0003780305060000042
in the formula (I), the compound is shown in the specification,
Figure FDA0003780305060000043
is a transfer function F' pss The phase angle, γ ', of(s) is the transfer function G' pss (s) phase angle;
according to the phase compensation method, there are
Figure FDA0003780305060000044
In formula (II) T' pssd Additional damping torque, T ', provided for virtual excitation PSS' psss Additional synchronous torque provided for the virtual excitation PSS;
namely, it is
Figure FDA0003780305060000045
Setting a transfer function G' pss (s) is a proportional link and a lead-lag link, and is represented as:
Figure FDA0003780305060000046
in the formula, T 1 ~T 4 All represent the time constant, K 'of the lead-lag link' 1 ,K′ 2 Setting a parameter of the virtual excitation PSS;
assuming that the virtual excitation PSS provides a positive damping torque, available from (16):
Figure FDA0003780305060000051
given T 1 、T 3 Obtaining T from the value of (19) 2 ,T 4 ,K′ 1 ,K′ 2 Obtaining a parametrically optimized transfer function G 'from equations (15) and (18)' pss (s);
And (3) optimizing parameters of the virtual angle PSS:
assuming that the virtual angle PSS provides a positive pure damping torque, Δ T ″ pss Correction quantity delta of power angle delta for electromagnetic torque provided by virtual angle PSS pss For outputting stable control signals of the virtual angle PSS, the electromagnetic torque provided by the virtual angle PSS is delta T ″ pss ,F″ pss (s) is from Δ δ pss To Δ T pss The transfer function of (c) then has:
Figure FDA0003780305060000052
ΔT″ pss =F″ pss (s)Δδ pss (21)
order transfer function G ″) pss (s) phase compensating transfer function F ″) pss (s) phase, setting the transfer function G ″) pss (s) is a proportional link and a lead-lag link, and is represented as:
Figure FDA0003780305060000053
in formula (II) T' 1 ~T′ 4 All represent the time constant, K ″, of the lead-lag link 1 ,K″ 2 Setting a parameter of a virtual angle PSS;
virtual excitation PSS in the same way
Figure FDA0003780305060000061
In the formula (I), the compound is shown in the specification,
Figure FDA0003780305060000062
is a transfer function F ″) pss The phase angle, γ ", of(s) is the transfer function G ″ pss (s) phase angle;
given T' 1 、T′ 3 Is obtained from the value of (23)' 2 ,T′ 4 ,K′ 1 ,K′ 2 The parameter-optimized transfer function G' is obtained from the formula (22) pss (s)。
CN202110331007.1A 2021-03-29 2021-03-29 Additional virtual double-PSS control method for active support type new energy unit Expired - Fee Related CN112952861B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110331007.1A CN112952861B (en) 2021-03-29 2021-03-29 Additional virtual double-PSS control method for active support type new energy unit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110331007.1A CN112952861B (en) 2021-03-29 2021-03-29 Additional virtual double-PSS control method for active support type new energy unit

Publications (2)

Publication Number Publication Date
CN112952861A CN112952861A (en) 2021-06-11
CN112952861B true CN112952861B (en) 2022-11-22

Family

ID=76227051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110331007.1A Expired - Fee Related CN112952861B (en) 2021-03-29 2021-03-29 Additional virtual double-PSS control method for active support type new energy unit

Country Status (1)

Country Link
CN (1) CN112952861B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105449699A (en) * 2016-01-11 2016-03-30 东北电力大学 Nonlinear fractional order auto disturbance rejection damping control method of doubly fed induction generators
CN107346892A (en) * 2017-04-24 2017-11-14 国家电网公司 PSS design methods based on the input hierarchical control of wide area multi signal
CN109149665A (en) * 2017-06-28 2019-01-04 中国电力科学研究院 Multi-rate simulating method and system for flexible direct current AC network associative simulation
CN109683037A (en) * 2018-12-08 2019-04-26 国网辽宁省电力有限公司电力科学研究院 A kind of power system stabilizer, PSS test method based on CoCo-80

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104659802B (en) * 2015-03-11 2016-10-05 云南电网有限责任公司电网规划研究中心 A kind of VSC-HVDC alternating voltage-frequency coordination control method improving AC system transient stability
CN106712055B (en) * 2015-11-12 2019-07-19 中国电力科学研究院 It is a kind of with the low power system stabilizer, PSS configuration method encouraging limitation function and mutually coordinating
CN106451421B (en) * 2016-09-27 2019-01-04 南方电网科学研究院有限责任公司 Modeling and ring-opening method for electromagnetic ring network-containing alternating current-direct current power grid
CN108011390B (en) * 2017-12-11 2019-10-08 华中科技大学 A kind of back-to-back flexible HVDC transmission system and bicyclic additional frequency control method
CN108199391B (en) * 2018-01-10 2018-12-07 云南电网有限责任公司电力科学研究院 A kind of control method for coordinating of generator UEL and PSS
WO2020173549A1 (en) * 2019-02-26 2020-09-03 Abb Power Grids Switzerland Ag Statcom arrangement comprising energy storage
CN110661273A (en) * 2019-08-20 2020-01-07 中国电力科学研究院有限公司 Damping control method and system for inhibiting low-frequency oscillation and subsynchronous oscillation
CN110970893A (en) * 2019-11-23 2020-04-07 国网辽宁省电力有限公司电力科学研究院 Method for adaptively optimizing and regulating network inertia level based on active support type VSC

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105449699A (en) * 2016-01-11 2016-03-30 东北电力大学 Nonlinear fractional order auto disturbance rejection damping control method of doubly fed induction generators
CN107346892A (en) * 2017-04-24 2017-11-14 国家电网公司 PSS design methods based on the input hierarchical control of wide area multi signal
CN109149665A (en) * 2017-06-28 2019-01-04 中国电力科学研究院 Multi-rate simulating method and system for flexible direct current AC network associative simulation
CN109683037A (en) * 2018-12-08 2019-04-26 国网辽宁省电力有限公司电力科学研究院 A kind of power system stabilizer, PSS test method based on CoCo-80

Also Published As

Publication number Publication date
CN112952861A (en) 2021-06-11

Similar Documents

Publication Publication Date Title
Sguarezi Filho et al. Model-based predictive control applied to the doubly-fed induction generator direct power control
CN107968591B (en) Grid-connected inverter virtual inertia power decoupling control method based on common point voltage prediction
CN111277001A (en) Fan grid-connected control method based on virtual synchronous generator parameter adaptive control
CN104079228B (en) The implicit trapezoidal rule electromagnetic transient modeling method with damping of double fed asynchronous machine
CN110611321B (en) Virtual power system stabilizer design method for compensating negative damping characteristic of virtual synchronous machine
CN111987956B (en) Torque ripple suppression method for direct-drive wind turbine generator
CN106533289B (en) A kind of non-linear voltage control method and system
CN110365051A (en) A kind of virtual synchronous motor control method of adaptive instruction filtering inverting
CN111049178A (en) Method for analyzing stability control of direct-drive permanent magnet wind turbine generator through VSC-HVDC grid connection
CN115882762A (en) Frequency optimization control method of grid-connected wind power system
CN103208817B (en) Second-order slip form-based method for controlling doubly-fed wind generator (DFIG)
Tsourakis et al. A power system stabilizer for variable-speed wind generators
CN116632866B (en) Hybrid energy storage self-adaptive inertia VSG control method for liquid flow super-capacity lithium battery
Dong et al. Voltage and reactive power control of front-end speed controlled wind turbine via H∞ strategy
CN112952861B (en) Additional virtual double-PSS control method for active support type new energy unit
CN110460098B (en) Wind turbine dual-mass-block shafting stability control method based on virtual mass block
CN112350346A (en) Method for establishing transient stability simulation model of power system
Yang et al. Adaptive passivity-based control of a TCSC for the power system damping improvement of a PMSG based offshore wind farm
Abulizi et al. Research of current control strategies for doubly-fed wind power generation system
CN109301867B (en) Virtual synchronous motor control method for simulating diesel generator set
CN115459301A (en) Virtual power system stabilizer design method for improving VSG stability
Phan et al. An improved control method for DFIG-based wind system supplying unbalanced stand-alone loads
CN112103986B (en) Wind power generation double-fed unit and synchronous grid-connected control method and device thereof
CN113612238A (en) Improved damping torque method for analyzing thermal power generating unit shafting small interference oscillation caused by wind power integration
Baesmat et al. Design of pole placement controllers for doubly-fed induction generators in the complex domain

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20221122

CF01 Termination of patent right due to non-payment of annual fee