Doubly-fed virtual synchronous generator control method based on fast reactive power control and power angle compensation method

Abstract

This patent relates to the field of control technology of doubly-fed motors, and specifically to a method for fast reactive power control and power angle compensation control of a doubly-fed virtual synchronous machine. The present invention proposes a decoupling scheme that uses a fast reactive power control strategy to eliminate the influence of active power disturbance on reactive power, and uses a power angle compensation strategy to eliminate the influence of reactive reference value disturbance on active power. On the basis of designing the overall control strategy of the doubly-fed virtual synchronous machine, a small signal model of the doubly-fed virtual synchronous machine taking into account the reactive control link is constructed. By analyzing the influence of the PI controller parameters of the reactive link on the power angle stability, a fast reactive control is designed. The grid frequency disturbance and reactive reference value disturbance conditions are set for verification. The results show that the proposed control strategy can effectively avoid the mutual influence between active power and reactive power when the doubly-fed virtual synchronous machine is running at a large power angle, and realize the stable operation of the doubly-fed virtual synchronous machine under grid frequency disturbance and reactive reference value disturbance.

Description

Technical field:

[0001] This patent relates to the field of control technology of doubly-fed motors, and specifically to a doubly-fed virtual synchronous generator control method based on fast reactive power control and power angle compensation method. Background technology:

[0002] With the trend of building a cleaner and low-carbon energy system, renewable energy will gradually become the main energy form for the future. Compared with traditional power sources, renewable energy with power electronic converters as grid interface lacks frequency inertia, and the increasing penetration rate of renewable power will further deteriorate the frequency stability control of the grid.

[0003] Doubly-fed wind turbines have a large installed capacity, and enhancing their inertia and frequency regulation capabilities will be of great significance to the frequency safety of the power grid. Frequency droop control is introduced into the traditional control strategy of doubly-fed units to reflect the frequency regulation characteristics, but this method cannot achieve frequency inertia support. Regarding the virtual inertia control strategy, the active power reference value is calculated by extracting the grid frequency deviation and its change rate, reflecting the inertia output and frequency regulation characteristics.

[0004] The virtual synchronous machine control strategy of the doubly fed unit is also a solution. By introducing the motion equation of the synchronous generator into the doubly fed unit control strategy, its output reflects the frequency inertia. At this time, the doubly fed unit outputs in the form of a voltage source, which can achieve self-synchronization with the power grid and has flexible on-grid and off-grid operation modes. Under the virtual synchronous control strategy, the power control mode is Pf and QU, which depends on the two conditions that the transmission line is inductive and the power angle is small. Wind farms are often connected to the end of the power grid, and the short-circuit ratio is relatively small; in the distributed scenario of local access to the load center, the voltage level of the access grid is low, and the line has a resistive-inductive characteristic. In such a scenario, the transmission power coupling is serious, which is easy to cause synchronous frequency resonance and endanger the stable operation of the power grid.

[0005] Introducing virtual resistance is another solution. The voltage drop caused by active power passing through the resistance is compensated by reactive power to achieve power decoupling transmission. The actual grid resistance is difficult to measure, and this parameter will change with the maximum and minimum operation modes of the grid, making it difficult to achieve accurate decoupling.

[0006] Some people also use the method of adding virtual reactance to increase the impedance angle of the transmission model so that it can approximately reflect the reactance characteristics and achieve approximate power decoupling. This method has good robustness but limited accuracy.

[0007] Some literatures disclose an adaptive virtual impedance method, which linearizes the doubly-fed virtual synchronous machine model using a small signal at a specific operating point, establishes the relationship between different steady-state operating points and the compensation impedance, and realizes the adaptation of the virtual impedance at different operating points.

[0008] Some literatures disclose solutions for characterizing the degree of coupling between active power and reactive power through d-axis current, without paying attention to the grid resistance value. The model for characterizing the degree of power coupling through current needs to be adjusted continuously with the steady-state operating point.

[0009] Some literatures disclose that the active disturbance compensation caused by voltage and the reactive disturbance compensation caused by power angle are converted into inner loop current compensation in reactive control through a virtual current feedforward power decoupling strategy.

[0010] The above decoupling schemes have a certain degree of similarity: through different compensation methods, the transmission model can reflect the inductive characteristics as much as possible and eliminate power coupling.

[0011] The doubly-fed virtual synchronous machine control method that controls power transmission with magnetizing potential introduces stator reactance into the power transmission model to reduce the power coupling problem caused by resistance. The larger stator reactance makes its operating power angle larger and has a wide range of variation, which brings new power coupling problems. Currently, there are few related studies. The doubly-fed virtual synchronous machine strategy that uses stator voltage as the target control voltage avoids the stator reactance appearing in the power transmission link, thereby avoiding the problem of large power angle operation of the virtual synchronous machine, but it cannot use the stator reactance to increase the reactance ratio in the transmission link. Summary of the invention:

[0012] This patent aims to solve the problem of high power angle operation caused by introducing stator reactance in the power transmission link, and proposes a fast reactive power control method by taking advantage of the different response time scales of the active loop and reactive loop; and proposes a power angle compensation control method to solve the problem that the fast reactive control method strengthens the power angle and active power disturbance during reactive disturbance. The technical solution is as follows:

[0013] The control method of the double-fed virtual synchronous generator based on the fast reactive power control and power angle compensation method, the power angle control expression is:

[0014]

[0015] Where, T Ω is the mechanical torque, T e is the electromagnetic torque, ω v ,ω 0 ,ω g are virtual synchronous angular frequency, power frequency angular frequency and stator voltage angular frequency respectively, D p is the damping coefficient, k f is the frequency modulation coefficient, J is the inertia coefficient, p is the differential operator; θ v ,θ g are the d-axis virtual synchronous phase angle and the stator voltage phase angle respectively, and σ is the power angle;

[0016] The reactive power control expression is:

[0017] E * =(k p +k i ∫)(Q ref +k u (U * -U)-Q),

[0018] In the formula, E * is the target value of the magnetizing potential, Q ref is the reactive power given value, k p ,k i is the reactive power control PI parameter, k u is the voltage regulation coefficient, U * is the grid voltage reference value;

[0019] The excitation voltage control expression is:

[0020]

[0021] In the formula, is the excitation voltage in the dq coordinate system, k pd ,k id ,k qd ,k iq The current inner loop corresponds to the PI control parameters, ω r is the rotor angular frequency, L m ,L r is the stator-rotor mutual inductance and rotor inductance; I rd , I rq is the rotor current; and the rotor current given value is:

[0022]

[0023] Optimal solution 1: reactive power control PI parameter k p ,k i Through fast reactive power control optimization,

[0024] choose

[0025] Among them, U s is the grid voltage, X s is the stator synchronous reactance, at this time k p The parameters can make the reactive power control system have a faster response speed and avoid k i The unstable state caused by the parameters; choose a certain k i Parameters are adjusted to meet the convergence speed requirements of the links corresponding to the negative real axis poles in the time domain, thus achieving static error-free tracking of reactive power.

[0026] Preferred solution 2: the power angle link includes a reactive power angle compensation function, which can avoid active power impact under reactive given disturbance;

[0027] The steady-state component in the reactive power reference value is filtered out through differential-integral operation, and the reactive power change is extracted:

[0028]

[0029] Let ΔP = Kk σ1 ΔE+KE 0 k σ2 Δσ=0, the relationship between power angle and reactive power is:

[0030]

[0031] Then the reactive power angle compensation function is constructed. Description of the drawings:

[0032] Figure 1 It is a schematic diagram of the overall control method of the double-fed virtual synchronous machine in the embodiment.

[0033] Figure 2 It is a schematic diagram of a small signal model of power angle control including an excitation link in an embodiment.

[0034] Figure 3 Schematic diagram of the root locus of the transfer function of the power angle to the grid frequency disturbance in the embodiment.

[0035] Figure 4 Schematic diagram of the relationship between the exciting potential and the active power changing with the power angle in the embodiment.

[0036] Figure 5 It is a schematic diagram of the reactive power angle compensation control method in the embodiment.

[0037] Figure 6 2 is a schematic diagram of the response principle of the double-fed virtual synchronous system under the grid frequency disturbance in the embodiment.

[0038] Figure 7 2 is a schematic diagram of the principle of transient change of virtual synchronous frequency in an embodiment.

[0039] Figure 8 It is a schematic diagram of the response principle of the double-fed virtual synchronous system under the reactive power set value disturbance in the embodiment. Specific implementation method:

[0040] Example:

[0041] As a rotating electromechanical device, the doubly-fed machine needs to meet frequency constraints for stable operation, namely:

[0042] n p ω Ω +ω r =ω,

[0043] Where np is the pole pair number, ω Ω is the rotor mechanical angular frequency, ω r is the rotor current angular frequency, ω is the stator rotating magnetic field angular frequency; it can be seen that the rotor magnetic field rotation speed is the same as the stator magnetic field rotation speed, and has similar characteristics to the synchronous generator; therefore, introducing the excitation potential widely used in synchronous generator control to construct a doubly fed machine model will be more conducive to the design of its virtual synchronous control method.

[0044] The stator and rotor of the doubly fed motor follow the convention of the motor. In the dq coordinate system, the direction of the exciting potential vector is the d-axis direction, and the magnetic flux equation is:

[0045]

[0046] In the formula, ψ s ,ψ r are the stator and rotor flux, I s ,I r are the stator and rotor currents, L s ,L r ,L m are the stator and rotor self-inductance and mutual inductance respectively; the voltage equation is:

[0047]

[0048] Where U s ,U r are the stator and rotor voltages respectively, and p is the differential operator;

[0049] Therefore, the magnetizing potential E is:

[0050]

[0051] Right now

[0052] In the synchronous rotating coordinate system, all variables are constant in steady state, and the differential components are neglected to obtain:

[0053]

[0054] Electromagnetic torque T e The equation is:

[0055] T e =n p L m (I sq I rd -I sd I rq ),

[0056] The equation of motion is:

[0057]

[0058] Where, T Ω is the mechanical torque, ω v ,ω 0 ,ω g are virtual synchronous angular frequency, power frequency angular frequency and stator voltage angular frequency respectively, k f is the frequency modulation coefficient, D p is the damping coefficient, J is the inertia coefficient; θ v ,θ g are the d-axis virtual synchronous phase angle and the stator voltage phase angle respectively, and σ is the power angle;

[0059] Ignoring the stator internal resistance, the stator power equation is:

[0060]

[0061] Where X s is the stator reactance;

[0062] The vector relationship between the magnetizing potential, stator voltage, and stator current is:

[0063] E=U s +jI s ,

[0064] Active and reactive power in the magnetizing potential E and stator voltage U s When the current is transmitted between the stator and the stator, the current passes through the stator reactance X s , the power expression is not affected by the impedance characteristics of the power grid; compared with the virtual synchronous machine control method that takes the stator voltage as the target control voltage, the double-fed virtual synchronous machine control method based on the exciting potential has better adaptability to typical weak power grids and resistive-inductive access systems such as distributed wind power;

[0065] It can be seen from the stator power equation that, under the condition of constant rated active and reactive power, the higher the transmission voltage and the smaller the stator reactance value, the smaller the power angle; when the power angle σ is small enough, sinσ can be approximated to σ, and cosσ can be approximated to 1. At this time, the active power P is approximately proportional to σ; similarly, the reactive power Q is approximately linearly related to E, and the control relationship between active power and reactive power is linearized and decoupled;

[0066] Based on the motion equation, the relationship from torque to power angle is established; based on the magnetizing potential expression, the reactive power closed-loop control expression is designed as follows:

[0067] E * =(k p +k i ∫)(Q ref +k u (U * -U)-Q),

[0068] In the formula, Q ref is the reactive power given value, k p ,k i is the reactive power control PI parameter, k u is the voltage regulation coefficient, U * is the grid voltage reference value;

[0069] According to the flux equation and voltage equation, the relationship between the excitation potential E and the rotor excitation voltage is established, and the excitation voltage controller is designed as follows:

[0070]

[0071] The overall control block diagram of the system is as follows: Figure 1 shown.

[0072] In actual application scenarios, since the voltage level of the doubly-fed generator connected to the grid is relatively low; under rated conditions, the stator current will produce a large voltage phase difference on the stator reactance, that is, the power angle σ is large, and the above approximate relationship will no longer be satisfied; when the power angle σ is large to a certain extent, the excitation potential E will play a dominant role in active power regulation, and the power angle σ will play a dominant role in reactive power regulation. When the doubly-fed virtual synchronous machine is in a large power angle operation state, the active and reactive power control is nonlinear and strongly coupled.

[0073] In order to solve the problem that the system transfer function is difficult to solve under nonlinear and strong coupling conditions, a control method of establishing a small signal model at a specific working point is adopted. Since the reactive current inner loop control responds quickly, the output of the outer loop E can be realized. * Fast tracking; when studying the stability of the outer loop control, the inner loop transfer function is approximately simplified to 1.

[0074] Under high power angle operation conditions, the influence of the change of the exciting potential on the active power cannot be ignored, E cannot be approximated as a constant, and the influence of the reactive link σ on E should be included in the power angle control link.

[0075] At the working point (E 0 ,σ 0 ), the small signal expressions of active and reactive power are:

[0076]

[0077] Where K = U s / X s ,k σ1 = sinσ 0 ,k σ2 =cosσ 0 ;

[0078] The equation of motion is:

[0079]

[0080] The stator power equation is:

[0081]

[0082] The transfer function of the magnetizing potential to the power angle in the reactive link is:

[0083]

[0084] Where k p ,k i (k i ≠0) is the reactive PI controller parameter, A=k p / k i ,B=(1+Kk σ2 k p ) / Kk σ2 k i ;because:

[0085]

[0086] It can be seen that k p The smaller the G Eσ (s) The more hysteresis characteristics are reflected, the larger the k i It can speed up the process of exponential approach to the final value in the time domain and improve the tracking speed; p Increase, G Eσ (s) gradually reflects the proportional characteristic; the steady-state operating point σ 0 increases, the hysteresis increases.

[0087] The transfer function of the reactive link excitation potential to the power angle small signal model G Eσ (s) is substituted into the active power expression to obtain the small signal model of power angle control including the reactive link, such as Figure 2 shown.

[0088] The feedback loop can be combined into:

[0089]

[0090] Power Angle to Grid Frequency Disturbance Δω g The closed-loop transfer function is:

[0091]

[0092] In the parameter k p ,k i During the increase process, the root locus of the transfer function of the power angle closed-loop system is as follows: Figure 3 As shown, it has the following characteristics:

[0093] In k p ,ki When it is 0, the magnetizing potential E is a constant and is not affected by σ; at this time, the power angle system is a second-order system, and the transfer function G σω (s) has a conjugate pole A,k′, whose position is determined by the moment of inertia J v and the damping coefficient D. In k i = 0, the power angle system is still a second-order system, k p The parameter will affect the pole location; with k p Increase, the poles A, A′ will move along k i = 0 The trajectory moves toward B, B′ and converges to B, B′;

[0094] In k i ≠0, the power angle system is a third-order system, with an additional negative real axis pole. p unchanged, k i When the value increases, the position of the conjugate pole changes perpendicularly to the trajectory curve in the figure. The symbol × in the figure represents k p =0, different k i The corresponding pole position. It can be seen that as k i Increase, the pole first approaches the imaginary axis and then moves away, and there is a possibility of instability; at k i unchanged, k p When k increases, the conjugate poles change along the trajectory curve in the figure. p When it is larger, the pole position approaches B, B′, and k i The effect on conjugate poles is small;

[0095] Negative real axis poles follow k p Increases closer to the imaginary axis, as k i Increase away from the imaginary axis;

[0096] Based on the above characteristics, we choose k p >>1 / Kk σ2 By making the conjugate poles tend to B, B′, a faster response speed can be obtained while avoiding k i The unstable state caused by the parameters. Choose a certain k i Parameters are adjusted to meet the convergence speed requirements of the links corresponding to the negative real axis poles in the time domain, while achieving static error-free tracking of reactive power.

[0097] According to the above analysis, optimize the reactive outer loop control parameter k p ,k i At this time G Eσ (s)≈E 0 k σ1 / k σ2 , reflecting the proportional characteristics. Eσ(s)ΔσSubstituting this into (18), we get ΔQ = 0. Therefore, when the power angle is disturbed, the reactive power Q = Q n Maintain stability and have strong resistance to power angle disturbance. The expression of magnetizing potential is:

[0098]

[0099] Active power expression under the current control method:

[0100]

[0101] The relationship between the magnetizing potential E and the active power P with the power angle σ is as follows: Figure 4 shown.

[0102] Curve L 1 is P(E,σ), and its projection curve L 2 Denotes P(σ), whose projection curve L 3 It is known that as σ increases, E and P both increase, and the power angle characteristic changes according to the tangent law. Em When the magnetizing potential reaches its maximum value E max , the power angle characteristics change into sinusoidal changes.

[0103] From the active power expression, it is known that the active power disturbance caused by the change of reactive power set value will be eliminated by adjusting the power angle σ. Therefore, the reactive power command generated by the change of reactive power set value or the deviation of grid voltage from the rated value will cause power angle disturbance. Due to the inertia characteristics of the virtual synchronous machine control method, the active power disturbance will be larger and last longer.

[0104] For reactive power reference value Q * In order to solve the active disturbance problem caused by the disturbance, a reactive power angle compensation control method is proposed, which aims to synchronously adjust the current power angle through power angle compensation when the reactive disturbance occurs, so as to offset the active disturbance caused by the reactive change. Figure 5 shown.

[0105] make:

[0106]

[0107] In the formula, ΔP = 0, and the relationship between power angle and reactive power is:

[0108] Let ΔP = Kk σ1 ΔE+KE 0 k σ2 Δσ=0, the relationship between power angle and reactive power is:

[0109]

[0110] The steady-state component in the reactive power reference value is filtered out through differential-integral operation, and its variation is extracted:

[0111]

[0112] G 1 (s) replaces the differential link and uses the virtual angular frequency ω v The integral link between the reactive power reference value Q * The compensation relationship between the power angle change Δσ. * When stable, the compensation amount is zero and does not affect the frequency inertia of the double-fed virtual synchronous machine; * When it changes, the power angle is adjusted linearly according to its change amount to achieve active power stability under the disturbance of reactive reference value.

[0113] In order to verify the effectiveness of the proposed control method, a grid-connected double-fed virtual synchronous wind turbine system model was built in Matlab / Simulink. Two working conditions, grid frequency disturbance and reactive power reference disturbance, were set for verification.

[0114] Operation characteristics of the double-fed virtual synchronous system under grid frequency disturbance. The double-fed wind power system is set to the rated state initially, and the fast reactive power control method is adopted; at 1s, the grid frequency f g Reduce 0.5Hz, and recover at 3s; at 5s f g The frequency increases by 0.5 Hz and recovers at 7 seconds. The changes of various variables in the double-fed wind power system are as follows: Figure 6 shown.

[0115] At the power angle σ, the exciting potential E, and the stator active power P s and stator reactive power Q s A control is set in the waveform to verify the influence of different reactive power control PI parameters on the dynamic characteristics. 1 The PI parameters corresponding to the fast reactive power control method are the same as k in (15). p >>1 / Kk σ2 Consistent; waveform p 2 Corresponding to k in (15) p <<1 / Kk σ2 And k i Smaller PI parameters. Other variables are responses under the fast reactive power control method.

[0116] Depend on Figure 6 It can be seen that at 1s grid frequency f g When the virtual synchronization frequency f v Inertial tracking grid frequency f g ; The accumulation of transient positive frequency difference increases the power angle σ, and the stator active power P s Increase. Grid frequency fg The deviation of the power frequency generates a frequency modulation active power instruction, making the virtual synchronous frequency f v With the grid frequency f g After synchronization, the power angle σ is still greater than the rated value, and the double-fed unit continues to output transient excess active power. The wind wheel angular kinetic energy provides inertia and frequency modulation power reserve, and the rotor speed r mech Reduce; mech The influence on the slip power is greater, and the active power P of the grid-side converter gsc The change is obvious with the decrease of speed. Under the current disturbance, the power P of the doubly fed unit g The peak value is -4.3MW, which is 1.26 times the rated value. The dynamic change mechanism is the same at 3s, 5s and 7s.

[0117] Compare the fast reactive power control method under each p 1 Waveform and non-control method of each p 2 The waveform shows that under the fast reactive power control method, the exciting potential E responds quickly to the change of the power angle σ, and then the stator reactive power Q s Maintain stability and achieve decoupling control of active power disturbances; when the power angle changes, the tangent power angle characteristic has a faster active response speed than the sine power angle characteristic, so the stator active power P s It can respond to grid frequency disturbances more quickly; the transient change in power angle is small, which is conducive to the stability of the power angle of the double-fed virtual synchronous machine system.

[0118] Power grid\synchronous frequency waveform f g \f v The details of area B are as follows Figure 7 As shown. After the grid frequency jumps, after 0.1s, the virtual synchronous frequency f v Tracking to the grid frequency f g , which is related to the virtual inertia, damping and reactive power control PI parameters.

[0119] Operation characteristics of the double-fed virtual synchronous system under reactive power reference disturbance: The double-fed wind power system is set to the rated state initially, and the fast reactive power control method is adopted. At 1s, the reactive power reference Q ref It changes from 0 to 0.8 MVar and recovers at 3 s; at 5 s, it changes from 0 to -0.8 MVar and recovers at 7 s. Figure 8 shown.

[0120] At the power angle σ, the exciting potential E, and the stator active power P s and stator reactive power Q s A control is set in the waveform of to verify the effect of the power angle compensation control method. on Corresponding to enable the control method; c offThe corresponding control method is disabled. The other variables are the responses when the power angle compensation control method is enabled.

[0121] Depend on Figure 8 It can be seen that at 1s, due to the change of reactive power reference value, reactive power reference value Q * The fast reactive power control adjusts the exciting potential E so that the stator reactive power Q s Fast Track Reference Q * , which will make the exciting potential E and stator active power P s Large transient fluctuations occur. As the active power returns to balance, the power angle σ and the magnetizing potential E reach a new steady-state operating point. The dynamic change mechanism at 3s, 5s and 7s is the same.

[0122] Compare the c under the reactive power angle compensation control method on waveform and deactivate each c of this control method off The waveform shows that: while the fast reactive power control method adjusts the exciting potential E, the reactive power angle compensation control method can synchronously adjust the power angle to offset the impact of the change in exciting potential E on the active power, and achieve the decoupling of active power from reactive power disturbance. The power angle compensation control method can effectively avoid excessive exciting potential transients and active power disturbances caused by reactive reference value disturbances, allowing the double-fed virtual synchronous system to quickly reach a new steady-state operating point.

[0123] In order to improve the voltage support capability, voltage differential reactive control is set up, the formula is as follows:

[0124] E * =(k p +k i ∫)(Q ref +k u (U * -U)-Q);

[0125] When the grid voltage is disturbed, the response of the doubly-fed virtual synchronous machine system to the voltage regulation reactive power command is similar to its response to the reactive power set value disturbance under this condition.

[0126] The fast reactive power control and power angle compensation control method proposed in this patent utilizes the different response time scales of active power and reactive power under the double-fed virtual synchronous machine control method to achieve reactive power stability under active power disturbance and active power stability under reactive disturbance, eliminating power coupling under the high power angle operation state of the double-fed virtual synchronous machine. Simulations under various disturbance conditions have proven the effectiveness of this control method.

Claims

1. A double-fed virtual synchronous generator control method based on fast reactive power control and power angle compensation method, characterized in that: The power angle control expression is: Where, T Ω is the mechanical torque, T e is the electromagnetic torque, ω v ,ω0,ω g are virtual synchronous angular frequency, power frequency angular frequency and stator voltage angular frequency respectively, D p is the damping coefficient, k f is the frequency modulation coefficient, J is the inertia coefficient, p is the differential operator; θ v ,θ g are the d-axis virtual synchronous phase angle and the stator voltage phase angle respectively, and σ is the power angle; The reactive power control expression is: E * =(k p +k i ∫)(Q ref +k u (U * -U)-Q), In the formula, E * is the target value of the excitation potential, Q ref is the reactive power given value, k p ,k i is the reactive power control PI parameter, k u is the voltage regulation coefficient, U * is the grid voltage reference value; The excitation voltage control expression is: In the formula, is the excitation voltage in the dq coordinate system, k pd ,k id ,k qd ,k iq The current inner loop corresponds to the PI control parameters, ω r is the rotor angular frequency, L m ,L r is the stator-rotor mutual inductance and rotor inductance; I rd , I rq is the rotor current; and the rotor current given value is:

2. The control method of a double-fed virtual synchronous generator based on a fast reactive power control and power angle compensation method according to claim 1 is characterized in that: Reactive power control PI parameter k p ,k i Through fast reactive power control optimization, choose Among them, U s is the grid voltage, X s is the stator synchronous reactance, at this time k p The parameters can make the reactive power control system have a faster response speed and avoid k i The unstable state caused by the parameters; choose a certain k i Parameters are adjusted to meet the convergence speed requirements of the links corresponding to the negative real axis poles in the time domain, thus achieving static error-free tracking of reactive power.

3. The control method of a double-fed virtual synchronous generator based on a fast reactive power control and power angle compensation method according to claim 1, characterized in that: The power angle link includes a reactive power angle compensation function, which can avoid active power impact under reactive given disturbance; The steady-state component in the reactive power reference value is filtered out through differential-integral operation, and the reactive power change is extracted: Let ΔP = Kk σ1 ΔE+KE0k σ2 Δσ=0, the relationship between power angle and reactive power is: Then the reactive power angle compensation function is constructed.

Images