Disclosure of Invention
The invention aims to overcome the defects in the prior art and designs a novel dynamic virtual current feedforward virtual synchronizer power second-order decoupling controller and a power decoupling control method thereof.
The invention can overcome the precision problem existing in the traditional VSG power decoupling, fully considers the second-order power angle variable quantity on the basis of the traditional virtual synchronous machine power first-order decoupling strategy, and realizes the coupling power compensation by utilizing a dynamic virtual current compensation mode based on the instantaneous reactive power theory.
In order to achieve the technical aim, the technical scheme provided by the invention is as follows:
the invention relates to a virtual synchronous machine (VSG) Power Second-Order Decoupling Controller (SOPDC), which is formed by sequentially connecting a given current calculation module (1), a reference current synthesis module (2) and a cross Decoupling module (3) in series, wherein the SOPDC is arranged at the rear stage of a traditional VSG Controller, input signals of the SOPDC are from a superior VSG Controller and a measuring unit, and output signals of the SOPDC are input into a PWM wave generation device of a VSG grid-connected inverter;
the given current calculation module (1) has 2 input ends and 2 output ends in total; the first and second input ends are respectively input with the induced electromotive force e output by the VSG controller and the voltage v of the grid-connected point provided by the measuring unitgThe first output end and the second output end of the signal are respectively connected with the first input end and the second input end of the reference current synthesis module (2);
the given current calculation module (1) consists of a transfer function module G(s) and a Park converter; the input of the transfer function G(s) is VSG induced electromotive force e and grid-connected point voltage vgIs output as a given current signal i*The calculation has been madeThe equation is shown in formula (1):
in the formula, s represents a complex variable, LeqIs the system equivalent inductance, ReqIs the system equivalent resistance, e is the induced electromotive force of the virtual synchronous machine, vgIs the grid-connected point voltage;
will give a current signal i
*The d-axis and q-axis outputs of the Park converter are respectively connected with a first output end and a second output end of the given current calculation module (1), and the output signals are given current signals i
*D-axis and q-axis components of
And
the reference current synthesis module (2) has 6 input ends and 2 output ends; the first input end and the second input end are respectively connected with the first output end and the second output end of the given current calculation module (1), and the third input end and the fourth input end are respectively input with a virtual current d-axis component signal delta i required by P-V, Q-delta decouplingd(P-V)And Δ id(Q-δ)The fifth input end and the sixth input end respectively input a virtual current q-axis component signal delta i required by P-V, Q-delta decouplingq(P-V)And Δ iq(Q-δ)(ii) a The first output end and the second output end are respectively connected with the first input end and the second input end of the cross decoupling module (3);
in the reference current synthesis module (2), a first input end inputs a signal
Virtual current d-axis component Δ i required for decoupling from P-V, Q- δ
d(P-V)、Δi
d(Q-δ)D-axis component signal i of the final value of the reference current is superposed and synthesized
d_refAnd output at a first output end of the reference current synthesis module (2) and input at a second input end
Virtual current q-axis components Δ i required for decoupling from P-V, Q- δ, respectively
q(P-V)、Δi
q(Q-δ)Superposing and synthesizing the q-axis component i of the final value of the reference current
q_refAnd is output at a second output end of the reference current synthesis module (2);
the cross decoupling module (3) has 6 input ends and 1 output end; the first and second input ends are respectively connected with the first and second output ends of the reference current synthesis module (2), and the input signals of the third and fourth input ends are respectively d-axis and q-axis components i of the output current of the virtual synchronous machinedAnd iqThe input signals of the fifth and sixth input terminals are respectively d-axis and q-axis components v of the grid-connected point voltagegdAnd vgqThe output end of the PWM wave generator is connected with a VSG grid-connected inverter;
the cross decoupling module (3) comprises two Gain modules Gain1 and Gain2 (the Gain is respectively omega L
eqAnd- ω L
eqWhere ω is the electrical angular velocity of VSG), two proportional-integral controllers PI1, PI2, a Park inverse transform module; input signal i of the first input terminal
d_refInput signal i to the third input terminal
dThe subtracted signals are input into a PI1 and an input signal i of a fourth input end
qThe output signal of PI1, the output signal of Gain2 and the input signal v of the fifth input end are sent to Gain2
gdD-axis component for generating induced electromotive force given value after superposition
Inputting a d-axis input end of a Park inverse transformation module; input signal i of the second input terminal
q_refAnd an input signal i to a fourth input terminal
qThe subtracted signals are input into a PI2 and an input signal i of a third input end
dThe output signal of PI2, the output signal of Gain1 and the input signal v of a sixth input end are sent into Gain1
gqQ-axis component for generating induced electromotive force given value after superposition
Inputting a q-axis input end of a Park inverse transformation module; park inverse transformationThe output of the module is used as the output end of the cross decoupling module (3), and the output signal is the set value e of the induced electromotive force
*;
The invention relates to a power decoupling control method of a virtual synchronizer power second-order decoupling controller, which comprises the following steps:
step 1: calculating the equivalent inductance L of the system according to the main circuit parameterseqEquivalent resistance Req:
Wherein, the equivalent inductance LeqIs calculated as in formula (2):
Leq=Ls+Lv+Lg (2)
in the formula, wherein LsFilter inductance of VSG, LvFor decoupling the virtual inductance, LgA network side line inductor;
equivalent resistance ReqIs calculated as shown in equation (4):
Req=Rs+Rv+Rg (3)
in the formula, wherein RsFilter resistance, R, for VSGvTo decouple the virtual resistance, RgA network side line resistor;
step 2: determining k of proportional-integral element in cross-decoupling module (3) according to conventional methodp,kiA parameter;
and step 3: based on the system equivalent inductance and resistance calculated in step 1, a given current i is calculated according to formula (1)*;
And 4, step 4: the dynamic virtual current required for power decoupling is calculated by using an instantaneous power theory in the following specific mode:
step 41: calculating the decoupling components delta i of the virtual current d axis and the q axis P-V needed by P-V decouplingd(P-V)And Δ iq(P-V)As shown in formula (4), wherein idAnd iqD-axis and q-axis components of the VSG steady-state stator current, respectively, and Δ V is the voltage variation caused by the disturbance, V0Is the VSG steady state operating point voltage;
step 42: calculating the required virtual current d-axis and Q-axis Q-delta decoupling components delta i for Q-delta decouplingd(Q-δ)And Δ iq(Q-δ)The method comprises the following steps:
step 421: calculating the decoupling components delta i of the d-axis and the Q-axis of the virtual current for compensating the Q-delta first-order coupling quantityd(Q-δ(1))And Δ iq(Q-δ(1))The calculation method is shown as formula (5), wherein idAnd iqD-axis and q-axis components of the VSG steady-state stator current, respectively, Δ δ is the power angle variation caused by the disturbance:
step 422: calculating decoupling components delta i of a d axis and a Q axis of virtual current for compensating Q-delta second-order coupling quantityd(Q-δ(2))And Δ iq(Q-δ(2))The calculation method is shown as formula (6), wherein idAnd iqD-and q-axis components, V, respectively, of VSG steady-state stator currentgTo the grid-connected point voltage amplitude, XeqIs a virtual inductor LeqThe corresponding reactance, Δ δ, is the power angle variation caused by the disturbance:
step 423: synthesizing the virtual current decoupling components obtained in the step 421 and the step 422 into a Q-delta virtual current decoupling component Δ id(Q-δ)And Δ iq(Q-δ)As shown in formula (7):
and 5: as shown in equation (8), a current i is given*D-axis component i ofd *And q-axis component iq *And the one obtained in step 4The P-V decoupling component and the Q-delta virtual current decoupling component are synthesized into a reference current final value id_refAnd iq_ref:
Step 6: the dq axis component i of the final value of the reference currentd_refAnd iq_refAn input cross decoupling module for outputting a given value e of an induced electromotive force*After the voltage is sent to a PWM generating device, a control signal is generated, and then the output of VSG is controlled;
the invention has the advantages that: the VSG power second-order decoupling controller based on dynamic virtual current feedforward is designed, the precision problem existing in the traditional VSG power decoupling is solved, the second-order power angle variable quantity is fully considered on the basis of the traditional first-order power decoupling strategy, and the coupling power compensation is realized by utilizing a dynamic virtual current compensation mode based on the instantaneous reactive power theory.
Detailed Description
The present invention will be described in further detail with reference to examples and drawings, but the present invention is not limited to these examples. The embodiment adopts a typical VSG grid-connected system as shown in FIG. 1, the topology is widely applied to research of VSG, and the effect of the invention can be tested by performing simulation analysis on the topology. The specific parameters of the main circuit are shown in table 1.
Table 1 main circuit parameters using VSG
Parameter(s)
|
Numerical value
|
Parameter(s)
|
Numerical value
|
Voltage V at DC sidedc |
42V
|
Network side inductor Lg |
4.5mH
|
Filter inductance Ls |
0.15mH
|
Network side resistor Rg |
0.045Ω
|
Filter capacitor Cs |
22uf
|
Effective value v of grid voltageg |
12V
|
Filter resistor Rs |
0.27Ω
|
|
|
Fig. 2 shows a control block diagram of the SOPDC proposed in the present patent in cooperation with a conventional VSG control method, wherein the VSG control shown in the first part of fig. 2 adopts a classical VSG structure, and its control parameters are listed in table 2:
TABLE 2 control parameters of VSG
Parameter(s)
|
Numerical value
|
Parameter(s)
|
Numerical value
|
Active droop coefficient DP |
0.2026
|
Reactive sag factor DQ |
117.88
|
Reactive sag factor DQ |
117.88
|
Coefficient of inertia without power K
|
74.0286
|
Virtual moment of inertia J
|
0.004
|
Reference angular velocity omegar |
314rad/s
|
Virtual resistance Rv |
0.005Ω
|
Virtual inductor Lv |
0.15mH |
In addition, according to the socdc connection structure, as shown in the second part of fig. 2, building the socdc controller, and according to the built VSG and the structure frame of the socdc, the steps of applying the socdc are as follows:
step 1: calculating the equivalent inductance L of the system by using the formula (2) and the formula (3) according to the main circuit parameterseqEquivalent resistance Req:
Leq=Ls+Lv+Lg=0.15mH+0.15mH+4.5mH=4.8mH
Req=Rs+Rv+Rg=0.27Ω+0.005Ω+0.045Ω=0.32Ω
Step 2: the parameters in the cross-decoupling module (3) are determined according to conventional methods, in this embodiment Kp=2,Ki=5;
And step 3: based on the system equivalent inductance and resistance calculated in step 1, given current i is calculated by using formula (1)*
And 4, step 4: the dynamic virtual current required by power decoupling is calculated by utilizing an instantaneous reactive power theory in the following specific mode:
step 41: calculating the d-axis and q-axis components Delta i of the virtual current required for P-V decoupling by using the formula (4)d(P-V)And Δ iq(P-V)
Step 42: calculating the required d-axis and Q-axis components Δ i of virtual current for Q-delta decouplingd(Q-δ)And Δ iq(Q-δ)The method comprises the following steps;
step 421: according to the equation (5), the d-axis and Q-axis components Δ i of the virtual current for compensating the Q- δ first-order coupling amount are calculatedd(Q-δ(1))And Δ iq(Q-δ(1));
Step 422:according to the equation (6), the d-axis and Q-axis components Δ i of the virtual current for compensating the Q- δ second-order coupling amount are calculatedd(Q-δ(2))And Δ iq(Q-δ(2));
Step 423: synthesizing the virtual currents obtained in the steps 321 and 322 into Δ i according to equation (7)d(Q-δ)And Δ iq(Q-δ)
And 5: according to equation (8), a given current i*Synthesized with the dummy current to a reference current final value id_refAnd iq_ref;
Step 6: the final value i of the reference currentd_refAnd iq_refAn input cross decoupling module for outputting a given value e of an induced electromotive force*And sending the control signal to a PWM generating device to generate a control signal so as to control the output of VSG;
in order to verify the superiority of the SOPDC provided by the present invention, this section sets the power decoupling effect of the VSG under two modes of common constant power control and Q-V droop control, comparing the following three control strategies under common inductive working condition: no decoupling strategy, a first-order decoupling strategy and a second-order decoupling strategy.
When VSG adopts constant power control, the simulation sets the initial state P of the systemset=0W,QsetAt 0.8s, the breaker is closed and the VSG is tied to the grid, 0 Var. Adjusting Q at 1.5ssetAdjusting P to 60Var, 6ssetTo 80W. Under the inductance resistance line parameters, active and reactive dynamic curves of VSG are respectively shown in fig. 3a and fig. 3b, both decoupling strategies have decoupling effects, the second-order decoupling strategy is higher in adjusting speed and smaller in reactive power impact compared with the first-order decoupling strategy, and the control effect of the SOPDC is obvious at the moment;
when VSG adopts Q-V control, setting initial time P in simulationset=0W,QsetWhen the value is 0Var, VSG grid connection is set at 0.8s, and P is set at 3ssetWhen the voltage is changed from 0 to 80W, Q is kept to be 0Var, and the active and reactive dynamic curves of VSG are respectively shown in fig. 4a and 4b under the resistance-inductance line parameter, after the SOPDC is adopted, the steady state deviation of the reactive power is obviously reduced compared with the other two control strategies, and the best decoupling effect among the three is achieved.
The embodiments described in this specification are merely illustrative of implementations of the inventive concept and the scope of the present invention should not be considered limited to the specific forms set forth in the embodiments but rather by the equivalents thereof as may occur to those skilled in the art upon consideration of the present inventive concept.