Virtual synchronous generator active power control method and device based on virtual inertia self-adaption and power angle correction

Abstract

The invention discloses a virtual synchronous generator active power control method and device based on virtual inertia self-adaption and power angle correction, through designing a virtual inertia self-adaption method, frequency overshoot during active power sudden change is reduced, and through designing a power angle correction strategy, frequency overshoot during active power sudden change is reduced. The contradictory problem that the active power of a conventional virtual synchronous generator has oscillation dynamically and has power deviation in a steady state is solved. On one hand, a virtual inertia self-adaption method is designed according to the dynamic change rule of active power and frequency, on the other hand, according to the angular frequency output by the virtual synchronous generator, power angle correction is achieved after adjustment of the proportional controller, so that the virtual synchronous generator has smaller frequency overshoot when the active power suddenly changes, and the virtual inertia self-adaption method is more suitable for the virtual synchronous generator. When the frequency of the power grid changes suddenly, the oscillation and overshoot of the active power are smaller, and there is no active power deviation in a steady state.

Description

Technical Field

[0001] The present invention relates to the technical field of grid connection, and in particular to a method and a device for controlling active power of a virtual synchronous generator based on virtual inertia self-adaptation and power angle correction. Background Art

[0002] In recent years, with the increasing penetration rate of new energy grid-connected, grid-connected converters have received more and more attention in order to improve the operational stability of power systems. Virtual synchronous generator technology is a key technology for realizing grid-connected converter operation, and it has also received widespread attention and application. However, conventional virtual synchronous generator control technology has two main disadvantages. One is that the output angular frequency of the virtual synchronous generator will produce a large overshoot when the active power changes. The other is that when the virtual damping is small, the active power dynamics have large oscillations and overshoots. When the virtual damping is large, although the oscillations and overshoots of the active power can be suppressed, the steady state will produce active power deviations. Although there are existing literatures that have studied this problem, few literatures solve the above two problems at the same time, and conventional methods often have the disadvantages of requiring the use of differential operations and low stability.

[0003] The literature "Yang Yun, Mei Fei, Zhang Chenyu, et al. Cooperative adaptive control strategy of moment of inertia and damping coefficient of virtual synchronous generator [J]. Electric Power Automation Equipment, 2019, 39(03): 125-131." studied a virtual synchronous generator active power control strategy based on cooperative adaptation of moment of inertia and damping coefficient, which reduced power oscillation and overshoot, but did not solve the problem of deviation in steady-state active power.

[0004] The literature "Lan Zheng, Long Yang, Zeng Jinhui, et al. Virtual synchronous generator transient power oscillation suppression strategy considering overshoot [J]. Automation of Electric Power Systems, 2022, 46(11): 131-141." studied a virtual synchronous motor transient power oscillation suppression strategy, but this method requires the use of differential operations, which are prone to amplify sampling noise, resulting in reduced system stability.

[0005] The literature "Liu Qin'e, Ren Dongfeng, Kang Yiqun, et al. Virtual synchronous machine active power control method based on damping power adjustment [J]. Smart Power, 2024, 52(03): 63-70." studied a virtual synchronous machine active power control method based on damping power adjustment, using a proportional integral controller to dynamically adjust the damping power, solving the problems of active power dynamic oscillation and steady-state static error, and avoiding the use of differential operations. However, this method does not consider the angular frequency overshoot problem when the active power changes suddenly.

[0006] From the above analysis, it can be seen that although the existing literature has studied the improved active power control strategy of virtual synchronous generators, the existing methods either require differential operations or fail to simultaneously solve the problems of angular frequency overshoot, active power dynamic oscillation and steady-state static error. Summary of the invention

[0007] In order to overcome the problems existing in the active power control of conventional virtual synchronous generators, the present invention discloses a method and a device for controlling the active power of a virtual synchronous generator based on virtual inertia adaptation and power angle correction.

[0008] In order to achieve the above object, the technical solution adopted by the present invention is:

[0009] A virtual synchronous generator active power control method based on virtual inertia adaptation and power angle correction, characterized in that it comprises the following steps:

[0010] Step 1: Difference the set grid rated angular frequency ω0 and the angular frequency ω output by the virtual synchronous generator control strategy, and multiply it by the droop coefficient k ω The first power deviation ΔP1 is obtained, which satisfies:

[0011] ΔP1=(ω0-ω)*k ω

[0012] Among them, the rated angular frequency ω0 of the power grid is 100πrad / s;

[0013] Step 2: Set the virtual synchronous generator active power reference value P ref Add the first power deviation ΔP1 obtained in step 1 and subtract the active power P output by the virtual synchronous generator e , and obtain the second power deviation ΔP2, which satisfies:

[0014] ΔP2=P ref +ΔP1-P e

[0015] Among them, the active power P e It is calculated based on the power angle output by the virtual synchronous generator at the previous moment;

[0016] Step 3: Integrate the second power deviation ΔP2 obtained in step 2 according to the following formula to obtain the first frequency deviation Δω1, which satisfies:

[0017]

[0018] Where, s represents the Laplace operator, J represents the adaptive virtual inertia of the virtual synchronous generator;

[0019] Step 4: Add the first frequency deviation Δω1 obtained in step 3 to the set grid rated angular frequency ω0 to obtain the angular frequency ω output by the virtual synchronous generator, satisfying:

[0020] ω=ω0+Δω1

[0021] Step 5: Subtract the actual angular frequency ω of the power grid from the angular frequency ω of the virtual synchronous generator output obtained in step 4 g After that, the integration operation is performed to obtain the initial power angle δ0 output by the virtual synchronous generator, which satisfies:

[0022]

[0023] Step 6: Multiply the angular frequency ω of the virtual synchronous generator output obtained in step 4 by the proportional coefficient k p Then the power angle correction value Δδ is obtained, and the initial power angle δ0 obtained in step 5 is added to obtain the final power angle δ output by the virtual synchronous generator, which satisfies:

[0024] δ=δ0+Δδ

[0025] Δδ=k p ω

[0026] Among them, k p is the proportional coefficient, and the power angle δ is used to calculate the active power P output by the virtual synchronous generator at the next moment. e .

[0027] Furthermore, the active power P output by the virtual synchronous generator in step 2 is e The calculation formula satisfies:

[0028] P e =Kδ'

[0029] in, U is the grid voltage amplitude, E is the amplitude of the virtual synchronous generator output voltage, X is the grid equivalent line impedance, and δ' is the power angle of the virtual synchronous generator output at the previous moment.

[0030] Furthermore, the calculation formula of the adaptive virtual inertia J of the virtual synchronous generator in step 3 satisfies:

[0031]

[0032] Where J0 is the initial virtual inertia of the virtual synchronous generator, k J is the adaptive proportional coefficient of virtual inertia, e is the natural base, equal to 2.71828, sgn() is the sign function, ΔP3 is the third power deviation, Δω2 is the second angular frequency deviation, and its calculation method satisfies:

[0033] ΔP3=Pe -P e '

[0034] Δω2=ω-ω'

[0035] Among them, P e is the active power output by the virtual synchronous generator, P e ' is P after delay time t e , ω is the angular frequency of the virtual synchronous generator output, and ω' is ω after the delay time t.

[0036] A virtual synchronous generator active power control device based on virtual inertia self-adaptation and power angle correction, comprising:

[0037] The first power deviation calculation module is used to make a difference between the set grid rated angular frequency ω0 and the angular frequency ω output by the virtual synchronous generator control strategy, and multiply it by the droop coefficient k ω The first power deviation ΔP1 is obtained, which satisfies:

[0038] ΔP1=(ω0-ω)*k ω

[0039] Among them, the rated angular frequency ω0 of the power grid is 100πrad / s;

[0040] The second power deviation calculation module is used to set the virtual synchronous generator active power reference value P ref Add the first power deviation ΔP1 calculated by the first power deviation calculation module, and subtract the active power P output by the virtual synchronous generator. e , and obtain the second power deviation ΔP2, which satisfies:

[0041] ΔP2=P ref +ΔP1-P e

[0042] Among them, the active power P e It is calculated based on the power angle output by the virtual synchronous generator at the previous moment;

[0043] The first frequency deviation calculation module is used to integrate the second power deviation ΔP2 calculated by the second power deviation calculation module according to the following formula to obtain a first frequency deviation Δω1, which satisfies:

[0044]

[0045] Where, s represents the Laplace operator, J represents the adaptive virtual inertia of the virtual synchronous generator;

[0046] The angular frequency calculation module is used to add the first frequency deviation Δω1 calculated by the first frequency deviation calculation module to the set grid rated angular frequency ω0 to obtain the angular frequency ω output by the virtual synchronous generator, satisfying:

[0047] ω=ω0+Δω1

[0048] The initial power angle calculation module is used to subtract the actual angular frequency ω of the power grid from the angular frequency ω of the virtual synchronous generator output calculated by the angular frequency calculation module. g After that, the integration operation is performed to obtain the initial power angle δ0 output by the virtual synchronous generator, which satisfies:

[0049]

[0050] The final power angle is the angular frequency ω of the virtual synchronous generator output calculated by the angular frequency calculation module multiplied by the proportional coefficient k p Then the power angle correction Δδ is obtained, and the initial power angle δ0 calculated by the initial power angle calculation module is added to obtain the final power angle δ output by the virtual synchronous generator, which satisfies:

[0051] δ=δ0+Δδ

[0052] Δδ=k p ω

[0053] Among them, k p is the proportional coefficient, and the power angle δ is used to calculate the active power P output by the virtual synchronous generator at the next moment. e .

[0054] Furthermore, the active power P output by the virtual synchronous generator e The calculation formula satisfies:

[0055] P e =Kδ'

[0056] in, U is the grid voltage amplitude, E is the amplitude of the virtual synchronous generator output voltage, X is the grid equivalent line impedance, and δ' is the power angle of the virtual synchronous generator output at the previous moment.

[0057] Furthermore, the calculation formula of the adaptive virtual inertia J of the virtual synchronous generator satisfies:

[0058]

[0059] Where J0 is the initial virtual inertia of the virtual synchronous generator, k Jis the adaptive proportional coefficient of virtual inertia, e is the natural base, equal to 2.71828, sgn() is the sign function, ΔP3 is the third power deviation, Δω2 is the second angular frequency deviation, and its calculation method satisfies:

[0060] ΔP3=P e -P e '

[0061] Δω2=ω-ω'

[0062] Among them, P e is the active power output by the virtual synchronous generator, P e ' is P after delay time t e , ω is the angular frequency of the virtual synchronous generator output, and ω' is ω after the delay time t.

[0063] Compared with the prior art, the present invention has the following beneficial effects:

[0064] Compared with conventional methods, the method of the present invention simultaneously solves the problems of angular frequency overshoot, active power dynamic oscillation and steady-state deviation suppression, and the method of the present invention does not need to introduce differential operations, so it is simpler to implement. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 This is the block diagram of active power control of conventional virtual synchronous generator;

[0066] Figure 2 It is a flow chart of a virtual synchronous generator active power control method based on virtual inertia adaptation and power angle correction proposed by the present invention;

[0067] Figure 3 The active power comparison simulation results of the conventional method and the method of the present invention are shown;

[0068] Figure 4 The figure shows the angular frequency comparison simulation results of the conventional method and the method of the present invention. DETAILED DESCRIPTION

[0069] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0070] like Figure 2 As shown, an embodiment of the present invention provides a virtual synchronous generator active power control method based on virtual inertia adaptation and power angle correction, comprising the following steps:

[0071] Step 1: Difference the set grid rated angular frequency ω0 and the angular frequency ω output by the virtual synchronous generator control strategy, and multiply it by the droop coefficient k ω The first power deviation ΔP1 is obtained, which satisfies:

[0072] ΔP1=(ω0-ω)*k ω

[0073] Among them, the rated angular frequency ω0 of the power grid is 100πrad / s;

[0074] Step 2: Set the virtual synchronous generator active power reference value P ref Add the first power deviation ΔP1 obtained in step 1 and subtract the active power P output by the virtual synchronous generator e , and obtain the second power deviation ΔP2, which satisfies:

[0075] ΔP2=P ref +ΔP1-P e

[0076] Among them, the active power P e It is calculated based on the power angle of the virtual synchronous generator output at the previous moment, satisfying:

[0077] P e =Kδ'

[0078] in, U is the grid voltage amplitude, E is the amplitude of the virtual synchronous generator output voltage, X is the grid equivalent line impedance, and δ' is the power angle of the virtual synchronous generator output at the previous moment.

[0079] Step 3: Integrate the second power deviation ΔP2 obtained in step 2 according to the following formula to obtain the first frequency deviation Δω1, which satisfies:

[0080]

[0081] Where s represents the Laplace operator and J represents the adaptive virtual inertia of the virtual synchronous generator, satisfying:

[0082]

[0083] Where J0 is the initial virtual inertia of the virtual synchronous generator, k J is the adaptive proportional coefficient of virtual inertia, e is the natural base, equal to 2.71828, sgn() is the sign function, ΔP3 is the third power deviation, Δω2 is the second angular frequency deviation, and its calculation method satisfies:

[0084] ΔP3=P e -Pe '

[0085] Δω2=ω-ω'

[0086] Among them, P e is the active power output by the virtual synchronous generator, P e ' is P after delay time t e , ω is the angular frequency of the virtual synchronous generator output, and ω' is ω after the delay time t.

[0087] Step 4: Add the first frequency deviation Δω1 obtained in step 3 to the set grid rated angular frequency ω0 to obtain the angular frequency ω output by the virtual synchronous generator, satisfying:

[0088] ω=ω0+Δω1

[0089] Step 5: Subtract the actual angular frequency ω of the power grid from the angular frequency ω of the virtual synchronous generator output obtained in step 4 g After that, the integration operation is performed to obtain the initial power angle δ0 output by the virtual synchronous generator, which satisfies:

[0090]

[0091] Step 6: Multiply the angular frequency ω of the virtual synchronous generator output obtained in step 4 by the proportional coefficient k p Then the power angle correction value Δδ is obtained, and the initial power angle δ0 obtained in step 5 is added to obtain the final power angle δ output by the virtual synchronous generator, which satisfies:

[0092] δ=δ0+Δδ

[0093] Δδ=k p ω

[0094] Among them, k p is the proportional coefficient, and the power angle δ is used to calculate the active power P output by the virtual synchronous generator at the next moment. e .

[0095] The embodiment of the present invention further provides a virtual synchronous generator active power control device based on virtual inertia adaptation and power angle correction, comprising:

[0096] The first power deviation calculation module is used to make a difference between the set grid rated angular frequency ω0 and the angular frequency ω output by the virtual synchronous generator control strategy, and multiply it by the droop coefficient k ω The first power deviation ΔP1 is obtained, which satisfies:

[0097] ΔP1=(ω0-ω)*k ω

[0098] Among them, the rated angular frequency ω0 of the power grid is 100πrad / s;

[0099] The second power deviation calculation module is used to set the virtual synchronous generator active power reference value P ref Add the first power deviation ΔP1 calculated by the first power deviation calculation module, and subtract the active power P output by the virtual synchronous generator. e , and obtain the second power deviation ΔP2, which satisfies:

[0100] ΔP2=P ref +ΔP1-P e

[0101] Among them, the active power P e It is calculated based on the power angle output by the virtual synchronous generator at the previous moment;

[0102] The first frequency deviation calculation module is used to integrate the second power deviation ΔP2 calculated by the second power deviation calculation module according to the following formula to obtain a first frequency deviation Δω1, which satisfies:

[0103]

[0104] Where, s represents the Laplace operator, J represents the adaptive virtual inertia of the virtual synchronous generator;

[0105] The angular frequency calculation module is used to add the first frequency deviation Δω1 calculated by the first frequency deviation calculation module to the set grid rated angular frequency ω0 to obtain the angular frequency ω output by the virtual synchronous generator, satisfying:

[0106] ω=ω0+Δω1

[0107] The initial power angle calculation module is used to subtract the actual angular frequency ω of the power grid from the angular frequency ω of the virtual synchronous generator output calculated by the angular frequency calculation module. g After that, the integration operation is performed to obtain the initial power angle δ0 output by the virtual synchronous generator, which satisfies:

[0108]

[0109] The final power angle is the angular frequency ω of the virtual synchronous generator output calculated by the angular frequency calculation module multiplied by the proportional coefficient k p Then the power angle correction Δδ is obtained, and the initial power angle δ0 calculated by the initial power angle calculation module is added to obtain the final power angle δ output by the virtual synchronous generator, which satisfies:

[0110] δ=δ0+Δδ

[0111] Δδ=k p ω

[0112] Among them, k p is the proportional coefficient, and the power angle δ is used to calculate the active power P output by the virtual synchronous generator at the next moment. e .

[0113] In order to verify the effectiveness of the proposed method, the active power control strategy of the virtual synchronous generator of the conventional grid-connected converter is compared with that of the conventional grid-connected converter. Figure 1 During the simulation, the active power reference value P ref 5kW before 4s, P at 4s ref The actual angular frequency of the power grid increases from 5kW to 10kW at 7s. g It suddenly decreases from 100πrad / s (50Hz) to 99.8πrad / s (49.9Hz). The virtual inertia J is 1kg / m 2 , the grid rated angular frequency ω0 is 100πrad / s, the droop coefficient k ω The voltage peak of the grid phase is 311V, and the peak voltage of the virtual synchronous generator output is 311. The grid line resistance is 2Ω, and the line inductance is 3mH.

[0114] Figure 3 The active power comparison simulation results of the conventional method and the method proposed by the present invention are given. Figure 4 The angular frequency comparison simulation results of the conventional method and the method proposed in the present invention are given. During the simulation, the virtual damping coefficient D of the conventional method is 0 and 10 Ws / rad respectively, and the proportional coefficient k of the method proposed in the present invention is p is 0.1, k J is 1, and the delay time t is 0.1s.

[0115] Depend on Figure 3 It can be seen that when the virtual damping coefficient D is zero, the active power and frequency output by the conventional virtual synchronous generator show obvious oscillation and overshoot. When the virtual damping coefficient D increases to 10, the oscillation is significantly suppressed. Figure 3 It can be seen that the increase of virtual damping coefficient will cause the grid frequency to deviate (after 7s), and the steady-state deviation of active power will also increase significantly. It can be seen that the conventional virtual synchronous generator active power control strategy cannot achieve the dynamic oscillation suppression of active power and the elimination of steady-state static error by adjusting the virtual damping coefficient. At the same time, when the active power suddenly changes at 4s, the angular frequency overshoot of the conventional method is also large, such as Figure 4 shown.

[0116] At the same time, by Figure 3 It can be seen that when the method proposed in the present invention is adopted, the active power output by the virtual synchronous generator has neither dynamic oscillation and overshoot nor steady-state active power deviation. Figure 4It can be seen that when the active power suddenly changes at 4s, the angular frequency overshoot of the method proposed in the present invention is also significantly reduced, which shows that the method proposed in the present invention can simultaneously achieve angular frequency overshoot suppression, active power dynamic oscillation suppression and steady-state static error elimination. This proves the effectiveness of the method proposed in the present invention.

[0117] The present invention is different from the conventional method in that the present invention introduces an adaptive virtual inertia J to reduce the angular frequency overshoot when the active power suddenly changes. At the same time, the present invention uses the angular frequency output by the virtual synchronous generator to correct the power angle, so that the active power control strategy of the virtual synchronous generator of the grid-type converter has no active power oscillation and overshoot in dynamic state and no active power static difference in steady state.

[0118] The above description is only a specific implementation mode of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by technicians in this technical field within the technical scope disclosed by the present invention should be covered by the protection scope of the present invention.

Claims

1. A virtual synchronous generator active power control method based on virtual inertia adaptation and power angle correction, characterized in that: The steps include: Step 1: Difference the set grid rated angular frequency ω0 and the angular frequency ω output by the virtual synchronous generator control strategy, and multiply it by the droop coefficient k ω The first power deviation ΔP1 is obtained, which satisfies: ΔP1=(ω0-ω)*k ω Among them, the rated angular frequency ω0 of the power grid is 100πrad / s; Step 2: Set the virtual synchronous generator active power reference value P ref Add the first power deviation ΔP1 obtained in step 1 and subtract the active power P output by the virtual synchronous generator e , and obtain the second power deviation ΔP2, which satisfies: ΔP2=P ref +ΔP1-P e Among them, the active power P e It is calculated based on the power angle output by the virtual synchronous generator at the previous moment; Step 3: Integrate the second power deviation ΔP2 obtained in step 2 according to the following formula to obtain the first frequency deviation Δω1, which satisfies: Where, s represents the Laplace operator, J represents the adaptive virtual inertia of the virtual synchronous generator; Step 4: Add the first frequency deviation Δω1 obtained in step 3 to the set grid rated angular frequency ω0 to obtain the angular frequency ω output by the virtual synchronous generator, satisfying: ω=ω0+Δω1 Step 5: Subtract the actual angular frequency ω of the power grid from the angular frequency ω of the virtual synchronous generator output obtained in step 4 g After that, the integration operation is performed to obtain the initial power angle δ0 output by the virtual synchronous generator, which satisfies: Step 6: Multiply the angular frequency ω of the virtual synchronous generator output obtained in step 4 by the proportional coefficient k p Then the power angle correction value Δδ is obtained, and the initial power angle δ0 obtained in step 5 is added to obtain the final power angle δ output by the virtual synchronous generator, which satisfies: δ=δ0+Δδ Δδ=k p oh Among them, k p is the proportional coefficient, and the power angle δ is used to calculate the active power P output by the virtual synchronous generator at the next moment. e .

2. The method for controlling active power of a virtual synchronous generator based on virtual inertia adaptation and power angle correction according to claim 1 is characterized in that: The active power P output by the virtual synchronous generator in step 2 e The calculation formula satisfies: P e =Kδ' in, U is the grid voltage amplitude, E is the amplitude of the virtual synchronous generator output voltage, X is the grid equivalent line impedance, and δ' is the power angle of the virtual synchronous generator output at the previous moment.

3. The method for controlling active power of a virtual synchronous generator based on virtual inertia adaptation and power angle correction according to claim 1, characterized in that: The calculation formula of the adaptive virtual inertia J of the virtual synchronous generator in step 3 satisfies: Where J0 is the initial virtual inertia of the virtual synchronous generator, k J is the adaptive proportional coefficient of virtual inertia, e is the natural base, equal to 2.71828, sgn() is the sign function, ΔP3 is the third power deviation, Δω2 is the second angular frequency deviation, and its calculation method satisfies: ΔP3=P e -P e ' Δω2=ω-ω' Among them, P e is the active power output of the virtual synchronous generator, P′ e is P after delay time t e , ω is the angular frequency of the virtual synchronous generator output, and ω' is ω after the delay time t.

4. A virtual synchronous generator active power control device based on virtual inertia adaptation and power angle correction, characterized in that: include: The first power deviation calculation module is used to make a difference between the set grid rated angular frequency ω0 and the angular frequency ω output by the virtual synchronous generator control strategy, and multiply it by the droop coefficient k ω The first power deviation ΔP1 is obtained, which satisfies: ΔP1=(ω0-ω)*k ω Among them, the rated angular frequency ω0 of the power grid is 100πrad / s; The second power deviation calculation module is used to set the virtual synchronous generator active power reference value P ref Add the first power deviation ΔP1 calculated by the first power deviation calculation module, and subtract the active power P output by the virtual synchronous generator. e , and obtain the second power deviation ΔP2, which satisfies: ΔP2=P ref +ΔP1-P e Among them, the active power P e It is calculated based on the power angle output by the virtual synchronous generator at the previous moment; The first frequency deviation calculation module is used to integrate the second power deviation ΔP2 calculated by the second power deviation calculation module according to the following formula to obtain a first frequency deviation Δω1, which satisfies: Where, s represents the Laplace operator, J represents the adaptive virtual inertia of the virtual synchronous generator; The angular frequency calculation module is used to add the first frequency deviation Δω1 calculated by the first frequency deviation calculation module to the set grid rated angular frequency ω0 to obtain the angular frequency ω output by the virtual synchronous generator, satisfying: ω=ω0+Δω1 The initial power angle calculation module is used to subtract the actual angular frequency ω of the power grid from the angular frequency ω of the virtual synchronous generator output calculated by the angular frequency calculation module. g After that, the integration operation is performed to obtain the initial power angle δ0 output by the virtual synchronous generator, which satisfies: The final power angle is the angular frequency ω of the virtual synchronous generator output calculated by the angular frequency calculation module multiplied by the proportional coefficient k p Then the power angle correction Δδ is obtained, and the initial power angle δ0 calculated by the initial power angle calculation module is added to obtain the final power angle δ output by the virtual synchronous generator, which satisfies: δ=δ0+Δδ Δδ=k p oh Among them, k p is the proportional coefficient, and the power angle δ is used to calculate the active power P output by the virtual synchronous generator at the next moment. e .

5. The virtual synchronous generator active power control device based on virtual inertia adaptation and power angle correction according to claim 4 is characterized in that: The active power P output by the virtual synchronous generator e The calculation formula satisfies: P e =Kδ' in, U is the grid voltage amplitude, E is the amplitude of the virtual synchronous generator output voltage, X is the grid equivalent line impedance, and δ' is the power angle of the virtual synchronous generator output at the previous moment.

6. The virtual synchronous generator active power control device based on virtual inertia adaptation and power angle correction according to claim 4 is characterized in that: The calculation formula of the adaptive virtual inertia J of the virtual synchronous generator satisfies: Where J0 is the initial virtual inertia of the virtual synchronous generator, k J is the adaptive proportional coefficient of virtual inertia, e is the natural base, equal to 2.71828, sgn() is the sign function, ΔP3 is the third power deviation, Δω2 is the second angular frequency deviation, and its calculation method satisfies: ΔP3=P e -P e ' Δω2=ω-ω' Among them, P e is the active power output of the virtual synchronous generator, P′ e is P after delay time t e , ω is the angular frequency of the virtual synchronous generator output, and ω' is ω after the delay time t.

Images