CN111030174B - Grid-connected inverter VSG mode and current source mode undisturbed switching control method - Google Patents

Grid-connected inverter VSG mode and current source mode undisturbed switching control method Download PDF

Info

Publication number
CN111030174B
CN111030174B CN201911308959.0A CN201911308959A CN111030174B CN 111030174 B CN111030174 B CN 111030174B CN 201911308959 A CN201911308959 A CN 201911308959A CN 111030174 B CN111030174 B CN 111030174B
Authority
CN
China
Prior art keywords
inverter
vsm
current
voltage
mode
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.)
Active
Application number
CN201911308959.0A
Other languages
Chinese (zh)
Other versions
CN111030174A (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201911308959.0A priority Critical patent/CN111030174B/en
Publication of CN111030174A publication Critical patent/CN111030174A/en
Application granted granted Critical
Publication of CN111030174B publication Critical patent/CN111030174B/en
Active 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/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers

Abstract

The invention discloses a method for controlling undisturbed switching between a grid-connected inverter VSG mode and a current source mode, which is applied toThe inverter topology structure using the control method comprises a direct current voltage source and a direct current side filter capacitor C in The three-phase inverter, the LC filter and the power grid; the DC voltage source and the DC side filter capacitor C in The grid-connected inverter is connected with a direct-current voltage source in series, and the output of the grid-connected inverter is filtered by an LC filter and then passes through a power grid impedance L g And accessing to a power grid. When the inverter works in a current source mode, tracking current source mode output by a VSG mode algorithm; when the inverter works in the VSG mode, the current source mode algorithm tracks the voltage source mode output, so that the double-mode undisturbed switching of the inverter is ensured, and the grid-connected stability of the inverter is improved.

Description

Grid-connected inverter VSG mode and current source mode undisturbed switching control method
Technical Field
The invention belongs to the technical field of distributed power generation and power electronics, and particularly relates to a grid-connected inverter VSG mode and current source mode undisturbed switching control method.
Background
With the worsening of the environment and the worsening of the energy crisis, new energy plays an increasingly important role in economic development and social progress. Compared with the traditional synchronous generator, the power accessed into the power grid through the current source control type inverter has the advantages of high control speed, capability of better tracking the maximum power point, good economic benefit and the like, but also has the defects of insufficient inertia and damping. When the new energy power generation permeability is improved, the instability problems such as harmonic resonance and the like are easily caused.
In order to ensure the safety and the friendliness of the electric power transmitted to the power grid by the new energy power generation equipment, many researchers at home and abroad propose a VSG (virtual synchronous generator) control strategy. The existing research shows that VSG control can simulate the inertia and damping characteristics of the traditional synchronous generator, and has better power grid friendliness and safety than current source control in a weak power grid state. However, the grid state changes from moment to moment, and it is desirable that the new energy power generation equipment can select a corresponding control mode according to the grid state, so as to deliver safe electric power to the grid to serve thousands of households.
Therefore, the switching between the inverter voltage control type strategy and the current control type strategy is a problem to be researched. Regarding the problem, the existing research schemes at home and abroad mainly focus on the undisturbed switching aspect of grid-connected and off-grid inverters.
At present, for a disturbance-free switching scheme of multiple control modes of an inverter, a plurality of academic papers have been analyzed and proposed to solve the problem, for example:
1. the subject is a research on a microgrid on-line/off-line smooth switching control strategy based on an energy storage generalized control algorithm, the Chinese Motor engineering project, and an article on 2840-2853 of No. 10 in 2019. The paper provides a generalized control algorithm suitable for an energy storage grid-connected inverter, and microgrid on-grid/off-grid smooth switching is achieved. However, the method aims at the grid-connected and grid-disconnected scene, and the inverter is disconnected when working in the voltage source mode, so that the method is not suitable for the scene that the inverter is always connected to the grid.
2. The invention discloses a grid-connected and grid-disconnected switching device and system which are disclosed in the patent document of China (publication number CN 108899935B) in 2019, 10 and 18.A grid-connected and grid-disconnected switching device and system are provided by the invention, wherein the grid-connected and grid-disconnected switching device comprises a grid-connected point switch, a bidirectional tide inverter, an island detection device and an energy storage device, but the grid-connected and grid-disconnected switching device needs the energy storage device and a communication device, so that the complexity of the system is increased, and the coordination control is not facilitated.
3. In the invention patent document of china (publication No. CN 110021959A) in "a grid-connected inverter dual-mode control method based on short-circuit ratio under weak grid" published in 2019, 7, 16, the invention provides a switching method of droop mode control and current source control, but the voltage source control method of the invention is droop control, and does not relate to a switching method of VSG control and current source control. Compared with droop control, VSG control has inertia and damping effects, and can better simulate the characteristics of a synchronous generator.
In view of the above documents, the dual-mode undisturbed switching of the existing inverter has the following disadvantages:
1. the existing research on the undisturbed switching of the inverter mainly focuses on grid-connected and off-grid switching, and the grid-connected and off-grid undisturbed switching mode is not suitable for the application scene of grid-connected power generation all the time;
2. the existing grid-connected inverter undisturbed switching mode mainly focuses on undisturbed switching between current source control and droop control, and VSG control has more damping and inertia than droop control, so that the voltage and frequency of a synchronous generator supporting grid can be better simulated.
Disclosure of Invention
The invention provides a grid-connected inverter VSG mode and current source mode undisturbed switching control method, an inverter applying the control method can complete switching from VSG control to current source control or from current source control to VSG control, and the output voltage and current of the inverter are smooth and undisturbed in the switching process.
The object of the invention is thus achieved. The invention provides a method for controlling undisturbed switching of a VSG mode and a current source mode of a grid-connected inverter, which can realize switching of inverter control modes under different power grid states. Inverter control includes voltage source control and current source control. When the inverter operates in a voltage source mode, the current source control algorithm tracks the voltage source control output; when the inverter is operating in a current source mode, the voltage source control algorithm tracks the current source output, thereby ensuring undisturbed switching between the voltage source and current source modes. The modulation signal passes through the dsp and the switching device to obtain the inverter side voltage, and passes through the filter inductor L of the LC filter f And filter capacitor C of LC filter f Through the network impedance L g And accessing to a power grid.
In particular to a grid-connected inverterThe inverter topology structure applying the control method comprises a direct current voltage source V dc DC side filter capacitor C in Grid-connected inverter, LC filter and grid impedance L g And Grid; the DC voltage source V dc And a DC side filter capacitor C in In parallel, the grid-connected inverter and the DC side filter capacitor C in In parallel connection, the output of the inverter is filtered by an LC filter and then passes through a power grid impedance L g Accessing to a power grid;
the control method comprises the following steps:
step 1, sampling output phase voltage U of inverter oa ,U ob ,U oc Sampling inverter side filter inductor phase current I of inverter inva ,I invb ,I invc Sampling of the output phase current I of the inverter oa ,I ob ,I oc And setting an inverter operation mode and setting: 0 is the current source mode, 1 is the VSG mode, all voltage source modes in the text represent VSG mode;
step 2, firstly, setting the operation mode of the inverter to be 0, and sampling the output phase voltage U of the three-phase inverter according to the step 1 to obtain oa ,U ob ,U oc Obtaining q-axis voltage V of output phase voltage of the three-phase inverter through a conversion formula from three-phase voltage to a two-phase rotating coordinate system csm_q And d-axis voltage V csm_d Then obtaining the angular frequency omega of the current source through a single synchronous coordinate system phase-locking formula csm Wherein the d axis is an active axis and the q axis is a reactive axis;
the transformation formula from the three-phase voltage to the two-phase rotating coordinate system is as follows:
Figure BDA0002323979560000041
Figure BDA0002323979560000042
the phase locking formula of the single synchronous coordinate system is as follows:
ω csm =V csm_q ×(K spll_p +K spll_i /s)+100π
wherein K spll_p Is the spll PLL proportioner coefficient, K spll_i The coefficient of a phase-locked loop integral regulator of a single synchronous coordinate system is obtained, s is a Laplace operator, and theta' is the output voltage angle of the three-phase inverter in the previous period;
step 3, according to the angular frequency omega of the current source obtained in the step 2 csm Obtaining the output voltage angle theta of the three-phase inverter through a mode selection integral formula 1;
the mode selection integral formula 1 is:
Figure BDA0002323979560000043
step 4, according to the inverter output voltage angle theta obtained in the step 3, the inverter side filter inductance phase current I obtained in the step 1 through the current source mode inverter side phase current rotating coordinate transformation equation 1 inva ,I invb ,I invc Coordinate transformation is carried out to obtain d-axis feedback I of filter inductance current at the inverter side in the current source mode of the inverter csm_dfeedback And current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback
The current source mode inverter side phase current rotation coordinate transformation equation 1 is as follows:
Figure BDA0002323979560000051
Figure BDA0002323979560000052
step 5, according to the filtered inductive current d-axis feedback I obtained in the step 4, on the inverter side of the inverter in the current source mode csm_dfeedback Current source mode inverter side filter inductance current q axis feedback I csm_qfeedback And externally given electricityStream source d-axis reference I csm_dref Current source q-axis reference I csm_qref Obtaining the current source mode modulation d-axis reference e through a current source mode current loop equation csm_dref And current source mode modulation q-axis reference e csm_qref
The current loop equation of the current source mode is as follows:
Figure BDA0002323979560000053
Figure BDA0002323979560000054
wherein, K csm_ip Is the current source mode current loop proportional regulator coefficient, K csm_ii Current source mode current loop integral regulator coefficients;
step 6, modulating d-axis reference e according to the current source mode obtained in the step 5 csm_dref And current source mode modulation q-axis reference e csm_qref Selecting equation 1 through a modulation mode to obtain a modulation wave e under a dq coordinate system d And e q Then modulating the wave to make the inverter pass through the grid impedance L g Outputting electric energy to a power grid;
the modulation mode selection equation 1 is:
e d =e csm_dref
e q =e csm_qref
and 7, converting the three-phase inverter output phase voltage U obtained by sampling in the step 1 according to a conversion formula 1 from the output three-phase voltage to a two-phase static coordinate system oa ,U ob ,U oc Converting to obtain two-phase static coordinate system component U of inverter output phase voltage ,U Converting the inverter output phase current I obtained by sampling in the step 1 according to the conversion formula 1 from the output three-phase current to the two-phase static coordinate system oa ,I ob ,I oc Converting to obtain the two-phase static coordinate system component I of the inverter output phase current ,I (ii) a Then, the inverse is obtained by the power calculation formula 1The inverter outputs active power P and the inverter outputs reactive power Q; then the VSG angular frequency omega is obtained through a VSG reactive power control formula 1 vsm Obtaining d-axis voltage amplitude reference U through VSG active power control formula 1 dref
The conversion formula 1 from the output three-phase voltage to the two-phase static coordinate system is as follows:
Figure BDA0002323979560000061
Figure BDA0002323979560000062
the conversion formula 1 from the output three-phase current to the two-phase static coordinate system is as follows:
Figure BDA0002323979560000063
Figure BDA0002323979560000064
the power calculation formula 1 is:
P=U I +U I
Q=U I -U I
the VSG active power control formula 1 is:
Figure BDA0002323979560000065
wherein, P ref To rated active power, omega 0 At the fundamental angular frequency, ω csm The current source angular frequency is adopted, m is an active droop coefficient, J is a virtual moment of inertia, and D is a virtual damping coefficient;
the VSG reactive power control formula 1 is as follows:
Figure BDA0002323979560000066
wherein: q ref Rated reactive power, V grms Is the effective value of the fundamental wave of the power grid, n is the reactive droop coefficient, I d_ref_vsm "is a given value of a period on d-axis current of a current loop in VSG mode, and I is set when a current source mode is started d_ref_vsm "has an initialization value of 0;
step 8, according to the filtered inductive current d-axis feedback I obtained in the step 4 at the side of the inverter in the current source mode of the inverter csm_dfeedback And current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback Obtaining the d-axis feedback I of the filter inductance current at the inverter side in the VSG mode through a voltage source mode inverter side phase current rotation coordinate transformation equation 1 vsm_dfeedback And VSG mode inverter side filter inductor current q-axis feedback I vsm_qfeedback And (3) converting the inverter output phase voltage U sampled in the step (1) through a voltage source mode inverter side phase voltage rotating coordinate transformation equation 1 oa ,U ob ,U oc Coordinate transformation is carried out to obtain d-axis feedback U of voltage source mode output voltage of the inverter vsm_dfeedback And inverter voltage source mode output voltage q-axis feedback U vsm_qfeedback
The voltage source mode inverter side phase current rotation coordinate transformation equation 1 is as follows:
I vsm_dfeedback =I csm_dfeedback
I vsm_qfeedback =I csm_qfeedback
the voltage source mode inverter side phase voltage rotating coordinate transformation equation 1 is as follows:
Figure BDA0002323979560000071
Figure BDA0002323979560000072
step 9, d-axis feedback U is carried out according to the output voltage of the inverter voltage source mode obtained in the step 8 vsm_dfeedback Inverter voltage source mode output voltage q-axis feedback U vsm_qfeedback And U obtained in step 7 drefv Obtaining d-axis reference I of voltage source mode current loop through voltage source mode voltage loop equation 1 vsm_dref And voltage source mode current loop q-axis reference I vsm_qref
Equation 1 for the voltage source mode voltage ring is:
I vsm_dref =(U drefv -U vsm_dfeedback )×(K vsm_vp +K vsm_vi /s)
I vsm_qref =(0-U vsm_qfeedback )×(K vsm_vp +K vsm_vi /s)
wherein K vsm_vp Is the VSG mode voltage loop proportional regulator coefficient, K vsm_vi Is the VSG mode voltage loop integral regulator coefficient;
step 10, according to the d-axis reference I of the voltage source mode current loop obtained in the step 9 vsm_dref Voltage source mode current loop q-axis reference I vsm_qref And d-axis feedback I of inverter VSG mode inverter side filter inductor current obtained in step 8 vsm_dfeedback And VSG mode inverter side filter inductor current q-axis feedback I vsm_qfeedback Obtaining the voltage source mode modulation d-axis reference e through the voltage source mode current loop equation 1 vsm_drefv And voltage source mode modulating q-axis reference e vsm_qrefv
The voltage source mode current loop equation 1 is:
e vsm_dref =(I vsm_dref -I vsm_dfeedback )×K vsm_ip +e csm_dref
e vsm_qref =(I vsm_qref -I vsm_qfeedback )×K vsm_ip +e csm_qref
wherein K vsm_ip Is the VSG mode current loop proportioner coefficient;
step 11, switching into a VSG mode, namely setting the operation mode as 1, and recording the operation of an inverter in the mode as a new period;
sampling output phase voltage U of new cycle of inverter oa ',U ob ',U oc ', sampling the inverter side filter inductance phase current I of the new cycle of the inverter inva ',I invb ',I invc ', sampling the output phase current I of the new cycle of the inverter oa ',I ob ',I oc ';
Step 12, according to the conversion formula 2 from the output three-phase voltage to the two-phase static coordinate system, the output phase voltage U of the new period of the inverter sampled in the step 11 is obtained oa ',U ob ',U oc ' conversion to obtain two-phase stationary coordinate system component U of output phase voltage in new cycle of inverter ',U '; converting the output three-phase current to a two-phase static coordinate system according to a conversion formula 2, and outputting the output phase current I of the inverter in the new period sampled in the step 11 oa ',I ob ',I oc ' conversion to obtain the output phase current two-phase stationary coordinate system component I of the new cycle of the inverter ',I '; secondly, obtaining the output active power P 'of the inverter in the new period and the output reactive power Q' of the inverter in the new period through a power calculation formula 2; then, a VSG angular frequency omega of a new period is obtained through a VSG reactive power control formula 2 vsm ', VSG active power control formula 2 obtains d-axis voltage amplitude reference U of a new period dref ';
The conversion formula 2 from the output three-phase voltage to the two-phase static coordinate system is as follows:
Figure BDA0002323979560000091
Figure BDA0002323979560000092
the conversion formula 2 from the output three-phase current to the two-phase stationary coordinate system is as follows:
Figure BDA0002323979560000093
Figure BDA0002323979560000094
the power calculation formula 2 is:
P=U 'I '+U 'I '
Q'=U 'I '-U 'I '
the active power control formula 2 is:
Figure BDA0002323979560000095
the reactive power control formula 2 is:
Figure BDA0002323979560000096
step 13, obtaining a new periodic angular frequency ω of the voltage source according to the step 12 csm ', obtaining an output voltage angle theta' of a new period of the inverter through a mode selection integral formula 2;
the mode selection integral equation 2 is:
Figure BDA0002323979560000097
step 14, according to the output voltage angle θ' of the inverter in the new period obtained in step 13, the inverter side filter inductor phase current I of the inverter in the new period obtained in step 11 is sampled by the voltage source mode inverter side phase current rotating coordinate transformation equation 2 inva ',I invb ',I invc ' coordinate transformation is carried out to obtain the feedback I of the side filter inductive current d axis of the voltage source mode inverter of a new period of the inverter vsm_dfeedback ' and voltage source module of new cycleFormula dc-to-ac converter side filter inductive current q axle feedback I vsm_qfeedback ' the output phase voltage U of the new cycle of the inverter sampled in the step 11 is converted by the voltage source mode inverter side phase voltage rotating coordinate transformation equation 2 oa ',U ob ',U oc ' coordinate transformation is carried out to obtain d-axis feedback U of output voltage of the voltage source mode inverter in a new period of the inverter vsm_dfeedback ' and new period of q-axis feedback U of output voltage of voltage source mode inverter vsm_qfeedback ';
The voltage source mode inverter side phase current rotation coordinate transformation equation 2 is as follows:
Figure BDA0002323979560000101
Figure BDA0002323979560000102
the voltage source mode inverter side phase voltage rotating coordinate transformation equation 2 is as follows:
Figure BDA0002323979560000103
Figure BDA0002323979560000104
step 15, according to the filter inductance current d-axis feedback I of the voltage source mode inverter side of the new cycle of the inverter obtained in the step 14 vsm_dfeedback ' with new period of the voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ', the current source mode inverter side filter inductance current d-axis feedback I of the new cycle of the inverter is obtained through the current source mode inverter side phase current rotation coordinate transformation equation 2 csm_dfeedback ' and new periodic current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback ';
The current source mode inverter side phase current rotation coordinate transformation equation 2 is:
I csm_dfeedback '=I vsm_dfeedback '
I csm_qfeedback '=I vsm_qfeedback '
step 16, according to the d-axis feedback U of the output voltage of the voltage source mode inverter of the new period of the inverter obtained in the step 14 vsm_dfeedback ' and new period of q-axis feedback U of output voltage of voltage source mode inverter vsm_qfeedback ' and step 12. The new cycle of U drefv ' obtaining a new cycle of d-axis reference I of the voltage source mode current loop through the voltage source mode voltage loop equation 2 vsm_dref ' and new periodic Voltage Source mode Current Loop q-axis reference I vsm_qref ';
The voltage source mode voltage loop equation 2 is:
I vsm_dref '=(U drefv '-U vsm_dfeedback ')×(K vsm_vp +K vsm_vi /s)
I vsm_qref '=(0-U vsm_qfeedback ')×(K vsm_vp +K vsm_vi /s)
step 17, referring to I according to the d-axis of the new periodic voltage source mode current loop obtained in step 16 vsm_dref ', new periodic voltage source mode current loop q-axis reference I vsm_qref ', new cycle voltage source mode inverter side filter inductor current d-axis feedback I of new cycle obtained in step 14 vsm_dfeedback ' with new period of the voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ' and current source mode modulation d-axis reference e obtained in step 5 csm_dref Current source mode modulated q-axis reference e csm_qref Obtaining a new period of voltage source mode modulation d-axis reference e through a voltage source mode current loop equation 2 vsm_dref ' modulating the q-axis reference e with the new periodic voltage source pattern vsm_qref ';
The voltage source mode current loop equation 2 is:
e vsm_dref '=(I vsm_dref '-I vsm_dfeedback ')×K vsm_ip +e csm_dref
e vsm_qref '=(I vsm_qref '-I vsm_qfeedback ')×K vsm_ip +e csm_qref
step 18, obtaining a new period of voltage source mode modulation d-axis reference e according to step 17 vsm_dref ' modulating the q-axis reference e with the new periodic voltage source pattern vsm_qrefv ' obtaining a modulation wave e of a new period in the dq coordinate system by selecting equation 2 through the modulation mode d ' and e q ', then modulating the wave to pass the inverter through the grid impedance L g Outputting electric energy to a power grid;
the modulation mode selection equation 2 is:
e d '=e vsm_dref '
e q '=e vsm_qref '
step 19, according to the new period obtained in step 15, the d-axis feedback I of the current source mode inverter side filter inductor current csm_dfeedback ', new period current source mode inverter side filter inductance current q axis feedback I csm_qfeedback ', externally given current source d-axis reference I csm_dref Current source q-axis reference I csm_qref And current source mode modulated d-axis reference e obtained in step 5 csm_dref And current source mode modulation q-axis reference e csm_qref Obtaining a new cycle of current source mode modulation d-axis reference e through a current source mode current loop equation 2 csm_dref ' modulating the q-axis reference e with a new periodic current source pattern csm_qref ';
The current loop equation 2 for the current source mode is:
Figure BDA0002323979560000121
preferably, the three-phase inverter output voltage angle θ ″ of the previous cycle in step 2 is initialized to 0 in the first cycle of the operation of the present control method.
Compared with the prior art, the invention has the beneficial effects that:
1. the method for controlling the undisturbed switching between the VSG mode and the current source mode of the grid-connected inverter has the advantages of simple structure, economy and effectiveness, and no need of adding other equipment;
2. the undisturbed switching realized by the invention can transmit green electric energy friendly to the power grid, and the environment is protected and the sustainable development is realized;
3. the VSG control of the invention has inertia and damping, can better simulate the performance of the synchronous generator, and plays a role in supporting the voltage and frequency of a power grid.
Drawings
Fig. 1 is a grid-connected topological diagram of a grid-connected inverter VSG mode and current source mode undisturbed switching control method.
Fig. 2 is a control block diagram of a method for controlling undisturbed switching between a grid-connected inverter VSG mode and a current source mode according to the invention.
Fig. 3 is a voltage waveform diagram of the VSG mode and current source mode undisturbed switching control grid connection when the short circuit ratio is 10.
Fig. 4 is a waveform diagram of the VSG mode and current source mode undisturbed switching control grid-connected current when the short-circuit ratio is 10.
Fig. 5 is a voltage waveform diagram of the VSG mode and current source mode undisturbed switching control grid when the short circuit ratio is 5.
Fig. 6 is a waveform diagram of the VSG mode and current source mode undisturbed switching control grid-connected current when the short-circuit ratio is 5.
Detailed Description
The present embodiment will be described in detail below with reference to the accompanying drawings.
Fig. 1 is a topological diagram of a method for controlling the undisturbed switching between a grid-connected inverter VSG mode and a current source mode. As can be seen from the figure, the topology includes a DC voltage source V dc DC side filter capacitor C in Grid-connected inverter, LC filter and grid impedance L g And a Grid; the DC voltage source V dc And a DC side filter capacitor C in In parallel, the grid-connected inverter and the DC side filter capacitor C in In parallel, the output of the inverter passes through an LC filterFiltered through the grid impedance L g And accessing to a power grid. In addition, in FIG. 1, L f Is the filter inductance of an LC filter, C f Is the filter capacitance of the LC filter, L g Is the grid impedance.
The specific parameters are as follows: the rated output line voltage of the inverter is 380V/50Hz, and the filter capacitor C at the direct current side in Filter inductance L of =15mf, lc filter f =0.56mh, filter capacitance C of lc filter f =270uF, rated inverter capacity of 100kVA, and grid impedance L when the short-circuit ratio is 5 g =0.56mH; when the short circuit ratio is 10, the network impedance L g =0.56mH。
Fig. 2 is a control block diagram of the method for controlling the switching between the grid-connected inverter in the VSG mode and the current source mode without disturbance. As can be seen from the figure, the steps of the invention are as follows:
step 1, sampling inverter output phase voltage U oa ,U ob ,U oc Sampling inverter side filter inductor phase current I inva ,I invb ,I invc Sampling inverter output phase current I oa ,I ob ,I oc Setting an inverter operation mode and setting: the voltage source mode is represented by the voltage source mode, i.e., VSG mode, where 0 is the current source mode and 1 is the VSG mode.
Step 2, firstly setting the operation mode of the inverter to be 0, and obtaining the output phase voltage U of the three-phase inverter according to the sampling in the step 1 oa ,U ob ,U oc Obtaining the q-axis voltage V of the output phase voltage of the three-phase inverter through a conversion formula from the three-phase voltage to the two-phase rotating coordinate system csm_q And d-axis voltage V csm_d Then obtaining the angular frequency omega of the current source through a single synchronous coordinate system phase-locking formula csm Wherein the d axis is an active axis and the q axis is a reactive axis;
the transformation formula from the three-phase voltage to the two-phase rotating coordinate system is as follows:
Figure BDA0002323979560000141
Figure BDA0002323979560000142
the phase locking formula of the single synchronous coordinate system is as follows:
ω csm =V csm_q ×(K spll_p +K spll_i /s)+100π
wherein K spll_p Is the spll PLL proportioner coefficient, K spll_i The method is characterized in that the method is a single synchronous coordinate system phase-locked loop integral regulator coefficient, s is a Laplace operator, and theta' is the output voltage angle of the three-phase inverter in the previous period. The three-phase inverter output voltage angle θ ″ of the previous cycle is initialized to 0 in the first cycle of the operation of the control method.
Step 3, according to the angular frequency omega of the current source obtained in the step 2 csm And obtaining the output voltage angle theta of the three-phase inverter through the mode selection integral formula 1.
The mode selection integral formula 1 is:
Figure BDA0002323979560000151
step 4, according to the inverter output voltage angle theta obtained in the step 3, the inverter side filter inductance phase current I obtained in the step 1 through the current source mode inverter side phase current rotating coordinate transformation equation 1 inva ,I invb ,I invc Coordinate transformation is carried out to obtain d-axis feedback I of filter inductance current at the inverter side in the current source mode of the inverter csm_dfeedback And current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback
The current source mode inverter side phase current rotation coordinate transformation equation 1 is as follows:
Figure BDA0002323979560000152
Figure BDA0002323979560000153
step 5, according to the filtered inductive current d-axis feedback I obtained in the step 4, on the inverter side of the inverter in the current source mode csm_dfeedback Current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback And an externally given current source d-axis reference I csm_dref Current source q-axis reference I csm_qref Obtaining the current source mode modulation d-axis reference e through a current source mode current loop equation csm_dref And current source mode modulation q-axis reference e csm_qref
The current loop equation of the current source mode is as follows:
Figure BDA0002323979560000161
Figure BDA0002323979560000162
wherein, K csm_ip Is the current source mode current loop proportional regulator coefficient, K csm_ii The current source mode current loop integrates the regulator coefficients. In this embodiment, K csm_ip =1.8,K csm_ii =5000。
Step 6, modulating the d-axis reference e according to the current source mode obtained in the step 5 csm_dref And current source mode modulation q-axis reference e csm_qref Selecting equation 1 through a modulation mode to obtain a modulation wave e under a dq coordinate system d And e q Then modulating the wave to make the inverter pass through the line inductance L g And outputting the electric energy to a power grid.
The modulation mode selection equation 1 is:
e d =e csm_dref
e q =e csm_qref
step 7, converting the inverter output phase voltage U obtained by sampling in the step 1 according to the conversion formula 1 from the output three-phase voltage to the two-phase static coordinate system oa ,U ob ,U oc Converting to obtain two-phase static coordinate system component U of inverter output phase voltage ,U And (3) converting the three-phase inverter output phase current I obtained by sampling in the step (1) according to a conversion formula 1 for outputting the three-phase current to a two-phase static coordinate system oa ,I ob ,I oc Converting to obtain the two-phase static coordinate system component I of the inverter output phase current ,I (ii) a Secondly, obtaining the output active power P of the inverter and the output reactive power Q of the inverter through a power calculation formula 1; then the VSG angular frequency omega is obtained through a VSG reactive power control formula 1 vsm Obtaining d-axis voltage amplitude reference U through VSG active power control formula 1 dref
The conversion formula 1 from the output three-phase voltage to the two-phase static coordinate system is as follows:
Figure BDA0002323979560000171
Figure BDA0002323979560000172
the conversion formula 1 from the output three-phase current to the two-phase static coordinate system is as follows:
Figure BDA0002323979560000173
Figure BDA0002323979560000174
the power calculation formula 1 is:
P=U I +U I
Q=U I -U I
the VSG active power control formula 1 is:
Figure BDA0002323979560000175
wherein, P ref To rated active power, omega 0 At the fundamental angular frequency, ω csm The current source angular frequency is shown, m is an active droop coefficient, J is a virtual moment of inertia, and D is a virtual damping coefficient. In this example P ref =100000,ω 0 =100×π,m=3.1415×10 -5 ,J=5,D=2。
The VSG reactive power control formula 1 is as follows:
Figure BDA0002323979560000176
wherein: q ref Rated reactive power, V grms Is the effective value of the fundamental wave of the power grid, n is the reactive droop coefficient, I d_ref_vsm "is a given value of a period on d-axis current of a current loop in VSG mode, and I is set when a current source mode is started d_ref_vsm "has an initialization value of 0. In this example Q ref =0,V grms =220,n=7×10 -5
Step 8, according to the filtered inductor current d-axis feedback I obtained in the step 4 on the inverter side of the current source mode of the inverter csm_dfeedback And current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback Obtaining the d-axis feedback I of the filter inductance current at the inverter side in the VSG mode through a voltage source mode inverter side phase current rotation coordinate transformation equation 1 vsm_dfeedback And VSG mode inverter side filter inductor current q-axis feedback I vsm_qfeedback (ii) a Converting equation 1 through the side phase voltage rotating coordinate of the voltage source mode inverter to obtain inverter output voltage U of the inverter sampled in step 1 oa ,U ob ,U oc Coordinate transformation is carried out to obtain d-axis feedback U of voltage source mode output voltage of the inverter vsm_dfeedback And inverter voltage source mode output voltage q-axis feedback U vsm_qfeedback
The voltage source mode inverter side phase current rotation coordinate transformation equation 1 is as follows:
I vsm_dfeedback =I csm_dfeedback
I vsm_qfeedback =I csm_qfeedback
the voltage source mode inverter side phase voltage rotating coordinate transformation equation 1 is as follows:
Figure BDA0002323979560000181
Figure BDA0002323979560000182
step 9, d-axis feedback U is output according to the voltage source mode output voltage of the inverter obtained in the step 8 vsm_dfeedback Inverter voltage source mode output voltage q-axis feedback U vsm_qfeedback And U obtained in step 7 drefv Obtaining d-axis reference I of voltage source mode current loop through voltage source mode voltage loop equation 1 vsm_dref And voltage source mode current loop q-axis reference I vsm_qref
Equation 1 for the voltage source mode voltage ring is:
I vsm_dref =(U drefv -U vsm_dfeedback )×(K vsm_vp +K vsm_vi /s)
I vsm_qref =(0-U vsm_qfeedback )×(K vsm_vp +K vsm_vi /s)
wherein K vsm_vp Is the VSG mode voltage loop proportional regulator coefficient, K vsm_vi Is the VSG mode voltage loop integral regulator coefficient. In this example K vsm_vp =1,K vsm_vi =400。
Step 10, referring to I according to d axis of voltage source mode current loop obtained in step 9 vsm_dref Voltage source mode current loop q-axis reference I vsm_qref And d-axis feedback I of inverter VSG mode inverter side filter inductor current obtained in step 8 vsm_dfeedback And VSG mode inverter side filter inductor current q-axis feedback I vsm_qfeedback With a voltage source mode current loop equation 1,obtaining a voltage source mode modulated d-axis reference e vsm_drefv And voltage source mode modulating q-axis reference e vsm_qrefv
The voltage source mode current loop equation 1 is:
e vsm_dref =(I vsm_dref -I vsm_dfeedback )×K vsm_ip +e csm_dref
e vsm_qref =(I vsm_qref -I vsm_qfeedback )×K vsm_ip +e csm_qref
wherein K vsm_ip Is the VSG mode current loop proportioner coefficient.
And step 11, finishing the current source mode operation, switching into the VSG mode, namely setting the operation mode as 1, and recording the operation of the inverter in the mode as a new period.
Sampling output phase voltage U of new cycle of inverter oa ',U ob ',U oc ', sampling the inverter side filter inductance phase current I of the new cycle of the inverter inva ',I invb ',I invc ', sampling the output phase current I of the new cycle of the inverter oa ',I ob ',I oc '。
Step 12, according to the conversion formula 2 from the output three-phase voltage to the two-phase static coordinate system, the output phase voltage U of the new period of the inverter sampled in the step 11 is obtained oa ',U ob ',U oc ' conversion to obtain two-phase stationary coordinate system component U of output phase voltage in new cycle of inverter ',U '; converting the output three-phase current to a two-phase static coordinate system according to a conversion formula 2, and outputting the output phase current I of the inverter in the new period sampled in the step 11 oa ',I ob ',I oc ' conversion to obtain the output phase current two-phase stationary coordinate system component I of the new cycle of the inverter ',I '; secondly, obtaining the output active power P 'of the inverter in the new period and the output reactive power Q' of the inverter in the new period through a power calculation formula 2; then, the VSG angular frequency omega of a new period is obtained through a VSG reactive power control formula 2 vsm ', VSG active power control formula 2 obtains d-axis voltage amplitude reference U of new period dref '。
The formula 2 for converting the output three-phase voltage into the two-phase static coordinate system is as follows:
Figure BDA0002323979560000201
Figure BDA0002323979560000202
the conversion formula 2 from the output three-phase current to the two-phase stationary coordinate system is as follows:
Figure BDA0002323979560000203
Figure BDA0002323979560000204
the power calculation formula 2 is:
P'=U 'I '+U 'I '
Q'=U 'I '-U 'I '
the active power control formula 2 is:
Figure BDA0002323979560000205
the reactive power control formula 2 is:
Figure BDA0002323979560000206
step 13, obtaining the angular frequency omega of the voltage source in the new period according to the step 12 csm 'the output voltage angle theta' of the new cycle of the inverter is obtained through the mode selection integral formula 2.
The mode selection integral formula 2 is:
Figure BDA0002323979560000207
step 14, according to the output voltage angle θ' of the inverter in the new period obtained in step 13, the inverter side filter inductor phase current I of the inverter in the new period obtained in step 11 is sampled by the voltage source mode inverter side phase current rotating coordinate transformation equation 2 inva ',I invb ',I invc ' coordinate transformation is carried out to obtain the feedback I of the side filter inductive current d axis of the voltage source mode inverter of a new period of the inverter vsm_dfeedback ' with new period of the voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ' the output phase voltage U of the new cycle of the inverter sampled in the step 11 is converted into the output phase voltage U of the new cycle of the inverter by the voltage source mode inverter side phase voltage rotating coordinate transformation equation 2 oa ',U ob ',U oc ' coordinate transformation is carried out to obtain the new period of d-axis feedback U of the filter inductive current at the side of the voltage source mode inverter vsm_dfeedback ' and new period of the q-axis feedback U of the filter inductor current on the side of the voltage source mode inverter vsm_qfeedback '。
The voltage source mode inverter side phase current rotation coordinate transformation equation 2 is as follows:
Figure BDA0002323979560000211
Figure BDA0002323979560000212
the voltage source mode inverter side phase voltage rotation coordinate transformation equation 2 is as follows:
Figure BDA0002323979560000213
Figure BDA0002323979560000214
step 15, according to the filtered inductive current d-axis feedback I of the inverter side of the voltage source mode of the new period of the inverter obtained in the step 14 vsm_dfeedback ' with new period of the voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ', the current source mode inverter side filter inductance current d-axis feedback I of the new cycle of the inverter is obtained through the current source mode inverter side phase current rotation coordinate transformation equation 2 csm_dfeedback ' and new periodic current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback ';
The current source mode inverter side phase current rotation coordinate transformation equation 2 is:
I csm_dfeedback '=I vsm_dfeedback '
I csm_qfeedback '=I vsm_qfeedback '
step 16, according to the d-axis feedback U of the output voltage of the voltage source mode inverter of the new period of the inverter obtained in the step 14 vsm_dfeedback ' New period voltage source mode inverter output voltage q-axis feedback U vsm_qfeedback ' and step 12 the new periodic d-axis voltage magnitude reference U drefv ' obtaining a new cycle of d-axis reference I of the voltage source mode current loop through the voltage source mode voltage loop equation 2 vsm_dref ' and new periodic voltage source mode current loop q-axis reference I vsm_qref '。
The voltage source mode voltage loop equation 2 is:
I vsm_dref '=(U drefv '-U vsm_dfeedback ')×(K vsm_vp +K vsm_vi /s)
I vsm_qref '=(0-U vsm_qfeedback ')×(K vsm_vp +K vsm_vi /s)
step 17, referring to I according to the d-axis of the new periodic voltage source mode current loop obtained in step 16 vsm_dref ', new periodic voltage source mode current loop q-axis referenceI vsm_qref ', new cycle voltage source mode inverter side filter inductor current d-axis feedback I obtained in step 14 vsm_dfeedback ', new period of filter inductor current q-axis feedback I at side of voltage source mode inverter vsm_qfeedback ' and current source mode modulation d-axis reference e obtained in step 5 csm_dref Current source mode modulation q-axis reference e csm_qref Obtaining a new period of voltage source mode modulation d-axis reference e through a voltage source mode current loop equation 2 vsm_dref ' modulating the q-axis reference e with the new periodic voltage source pattern vsm_qref '。
The voltage source mode current loop equation 2 is:
e vsm_dref '=(I vsm_dref '-I vsm_dfeedback ')×K vsm_ip +e csm_dref
e vsm_qref '=(I vsm_qref '-I vsm_qfeedback ')×K vsm_ip +e csm_qref
step 18, obtaining a new period of voltage source mode modulation d-axis reference e according to step 17 vsm_dref ' modulating the q-axis reference e with the new periodic voltage source pattern vsm_qrefv ' obtaining a modulation wave e of a new period in the dq coordinate system by selecting equation 2 through the modulation mode d ' and e q ' then modulating the wave to make the inverter pass through the network impedance L g And outputting the electric energy to a power grid.
The modulation mode selection equation 2 is:
e d '=e vsm_dref '
e q '=e vsm_qref '
step 19, according to the new period obtained in step 15, the d-axis feedback I of the current source mode inverter side filter inductor current csm_dfeedback ', new period current source mode inverter side filter inductance current q axis feedback I csm_qfeedback ', externally given current source d-axis reference I csm_dref Current source q-axis reference I csm_qref And current source mode modulated d-axis reference e obtained in step 5 csm_dref And current source mode modulated q-axis referencee csm_qref Obtaining a new cycle of current source mode modulation d-axis reference e through a current source mode current loop equation 2 csm_dref ' modulating the q-axis reference e with a new periodic current source pattern csm_qref '。
The current loop equation 2 for the current source mode is:
Figure BDA0002323979560000231
fig. 3 and 4 are a grid-connected voltage change waveform diagram and a grid-connected current waveform diagram respectively when the inverter with VSG control and current source control undisturbed switching control is connected to a power grid with short-circuit ratios of 10 and 5, the grid-connected current waveform diagram is that the grid-connected current ratio is 10 when 0-0.5s, the working mode signal input is 0, the grid-connected inverter works in the current source mode, the grid-connected current ratio is changed to 5 when 0.5s, at this time, the operation mode switching is required, the working mode signal input is 1 when 0.5-1s, and the grid-connected inverter works in the VSG mode; the grid short-circuit ratio is 10 when 1-1.5s, the working mode signal input is 0, the grid-connected inverter works in the current source mode, the grid short-circuit ratio is 5 when 1.5-2s, the working mode signal input is 1, the grid-connected inverter works in the VSG mode, and the operation is finished when 2 s. It can be seen that the voltage and current waveforms are undisturbed at the three switching moments of 0.5s, 1s and 1.5s, which illustrates the effectiveness of the scheme provided by the invention.
Fig. 5 and fig. 6 are a grid-connected voltage variation waveform diagram and a grid-connected current waveform diagram respectively when the inverter with VSG control and current source control undisturbed switching control is connected to a power grid with short-circuit ratios of 8 and 3, the grid-connected short-circuit ratio is 8 at 0-0.5s, the working mode signal input is 0, the grid-connected inverter works in the current source mode, the grid-connected short-circuit ratio is 3 at 0.5s, at this time, the operation mode switching is required, the working mode signal input is 1 at 0.5-1s, and the grid-connected inverter works in the VSG mode; the grid short-circuit ratio is 8 when 1-1.5s, the working mode signal input is 0, the grid-connected inverter works in the current source mode, the grid short-circuit ratio is changed into 3 when 1.5-2s, the working mode signal input is 1, the grid-connected inverter works in the VSG mode, and the operation is finished when 2 s. It can be seen that the voltage and current waveforms are undisturbed at the three switching moments of 0.5s, 1s and 1.5s, which illustrates the effectiveness of the scheme provided by the invention.
It can be seen from the figure that, under the weak grid, by adopting the VSG mode and current source mode undisturbed switching method provided by the patent, grid-connected voltage and current can be switched undisturbed well.

Claims (2)

1. A control method for undisturbed switching between VSG mode and current source mode of a virtual synchronous machine of a grid-connected inverter is disclosed, wherein the topological structure of the inverter applying the control method comprises a direct-current voltage source V dc DC side filter capacitor C in Grid-connected inverter, LC filter and grid impedance L g And a Grid; the DC voltage source V dc And a DC side filter capacitor C in Connected in parallel, the grid-connected inverter and a DC side filter capacitor C in In parallel connection, the output of the inverter is filtered by an LC filter and then passes through a power grid impedance L g Accessing a power grid;
the control method is characterized by comprising the following steps:
step 1, sampling output phase voltage U of an inverter oa ,U ob ,U oc Sampling inverter side filter inductor phase current I of inverter inva ,I invb ,I invc Sampling the output phase current I of the inverter oa ,I ob ,I oc Setting an inverter operation mode and setting: 0 is a current source mode, 1 is a VSG mode, and all voltage source modes represent the VSG mode;
step 2, firstly setting the operation mode of the inverter to be 0, and obtaining the output phase voltage U of the three-phase inverter according to the sampling in the step 1 oa ,U ob ,U oc Obtaining q-axis voltage V of output phase voltage of the three-phase inverter through a conversion formula from three-phase voltage to a two-phase rotating coordinate system csm_q And d-axis voltage V csm_d Then obtaining the angular frequency omega of the current source through a single synchronous coordinate system phase-locking formula csm Wherein the d axis is an active axis and the q axis is a reactive axis;
the transformation formula from the three-phase voltage to the two-phase rotating coordinate system is as follows:
Figure FDA0003784820060000011
Figure FDA0003784820060000012
the phase locking formula of the single synchronous coordinate system is as follows:
ω csm =V csm_q ×(K spll_p +K spll_i /s)+100π
wherein K spll_p For the spll PLL proportioner coefficient, K spll_i The coefficient of a phase-locked loop integral regulator of a single synchronous coordinate system is obtained, s is a Laplace operator, and theta' is the output voltage angle of the three-phase inverter in the previous period;
step 3, according to the angular frequency omega of the current source obtained in the step 2 csm Obtaining the output voltage angle theta of the three-phase inverter through a mode selection integral formula 1;
the mode selection integral formula 1 is:
Figure FDA0003784820060000021
step 4, according to the inverter output voltage angle theta obtained in the step 3, the inverter side filter inductance phase current I obtained in the step 1 through the current source mode inverter side phase current rotating coordinate transformation equation 1 inva ,I invb ,I invc Coordinate transformation is carried out to obtain the d-axis feedback I of the filter inductance current at the inverter side of the inverter in the current source mode of the inverter csm_dfeedback And current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback
The current source mode inverter side phase current rotation coordinate transformation equation 1 is as follows:
Figure FDA0003784820060000022
Figure FDA0003784820060000023
step 5, according to the filtered inductive current d-axis feedback I obtained in the step 4, on the inverter side of the inverter in the current source mode csm_dfeedback Current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback And an externally given current source d-axis reference I csm_dref Current source q-axis reference I csm_qref Obtaining the current source mode modulation d-axis reference e through a current source mode current loop equation csm_dref And current source mode modulation q-axis reference e csm_qref
The current loop equation of the current source mode is as follows:
Figure FDA0003784820060000031
Figure FDA0003784820060000032
wherein, K csm_ip Is the current source mode current loop proportional regulator coefficient, K csm_ii Current source mode current loop integral regulator coefficients;
step 6, modulating d-axis reference e according to the current source mode obtained in the step 5 csm_dref And current source mode modulation q-axis reference e csm_qref Selecting equation 1 through a modulation mode to obtain a modulation wave e under a dq coordinate system d And e q Then modulating the wave to make the inverter pass through the network impedance L g Outputting electric energy to a power grid;
the modulation mode selection equation 1 is:
e d =e csm_dref
e q =e csm_qref
step 7, converting the three-phase voltage to a two-phase static coordinate system according to the output three-phase voltageFormula 1 is that the three-phase inverter output phase voltage U obtained by sampling in the step 1 oa ,U ob ,U oc Converting to obtain two-phase static coordinate system component U of inverter output phase voltage ,U Converting the inverter output phase current I obtained by sampling in the step 1 according to the conversion formula 1 from the output three-phase current to the two-phase static coordinate system oa ,I ob ,I oc Converting to obtain the two-phase static coordinate system component I of the inverter output phase current ,I (ii) a Secondly, obtaining the output active power P of the inverter and the output reactive power Q of the inverter through a power calculation formula 1; then obtaining the VSG angular frequency omega through a VSG reactive power control formula 1 vsm Obtaining d-axis voltage amplitude reference U through VSG active power control formula 1 dref
The conversion formula 1 from the output three-phase voltage to the two-phase static coordinate system is as follows:
Figure FDA0003784820060000033
Figure FDA0003784820060000034
the conversion formula 1 from the output three-phase current to the two-phase static coordinate system is as follows:
Figure FDA0003784820060000041
Figure FDA0003784820060000042
the power calculation formula 1 is:
P=U I +U I
Q=U I -U I
the VSG active power control formula 1 is:
Figure FDA0003784820060000043
wherein, P ref To rated active power, omega 0 At the fundamental angular frequency, ω csm The current source angular frequency is adopted, m is an active droop coefficient, J is a virtual moment of inertia, and D is a virtual damping coefficient;
the VSG reactive power control formula 1 is as follows:
Figure FDA0003784820060000044
wherein: q ref Rated reactive power, V grms Is the effective value of the fundamental wave of the power grid, n is the reactive droop coefficient, I d_ref_vsm "is a given value of a period on d-axis current of a current loop in VSG mode, and I is set when a current source mode is started d_ref_vsm "has an initialization value of 0;
step 8, according to the filtered inductive current d-axis feedback I obtained in the step 4 at the side of the inverter in the current source mode of the inverter csm_dfeedback And current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback Obtaining the d-axis feedback I of the filter inductance current at the inverter side in the VSG mode through a voltage source mode inverter side phase current rotation coordinate transformation equation 1 vsm_dfeedback And VSG mode inverter side filter inductor current q-axis feedback I vsm_qfeedback And (3) converting the inverter output phase voltage U sampled in the step (1) through a voltage source mode inverter side phase voltage rotating coordinate transformation equation 1 oa ,U ob ,U oc Coordinate transformation is carried out to obtain d-axis feedback U of voltage source mode output voltage of the inverter vsm_dfeedback And inverter voltage source mode output voltage q-axis feedback U vsm_qfeedback
The voltage source mode inverter side phase current rotation coordinate transformation equation 1 is as follows:
I vsm_dfeedback =I csm_dfeedback
I vsm_qfeedback =I csm_qfeedback
the voltage source mode inverter side phase voltage rotating coordinate transformation equation 1 is as follows:
Figure FDA0003784820060000051
Figure FDA0003784820060000052
step 9, d-axis feedback U is output according to the voltage source mode output voltage of the inverter obtained in the step 8 vsm_dfeedback Inverter voltage source mode output voltage q-axis feedback U vsm_qfeedback And U obtained in step 7 drefv Obtaining d-axis reference I of voltage source mode current loop through voltage source mode voltage loop equation 1 vsm_dref And voltage source mode current loop q-axis reference I vsm_qref
Equation 1 for the voltage source mode voltage ring is:
I vsm_dref =(U drefv -U vsm_dfeedback )×(K vsm_vp +K vsm_vi /s)
I vsm_qref =(0-U vsm_qfeedback )×(K vsm_vp +K vsm_vi /s)
wherein K vsm_vp Is the VSG mode voltage loop proportional regulator coefficient, K vsm_vi Is the VSG mode voltage loop integral regulator coefficient;
step 10, according to the d-axis reference I of the voltage source mode current loop obtained in the step 9 vsm_dref Voltage source mode current loop q-axis reference I vsm_qref And d-axis feedback I of inverter VSG mode inverter side filter inductor current obtained in step 8 vsm_dfeedback And VSG mode inverter side filter inductor current q-axis feedback I vsm_qfeedback Obtaining a voltage source mode modulation d-axis reference e through a voltage source mode current loop equation 1 vsm_drefv And voltage source modeModulating q-axis reference e vsm_qrefv
The voltage source mode current loop equation 1 is:
e vsm_dref =(I vsm_dref -I vsm_dfeedback )×K vsm_ip +e csm_dref
e vsm_qref =(I vsm_qref -I vsm_qfeedback )×K vsm_ip +e csm_qref
wherein K vsm_ip Is the VSG mode current loop proportioner coefficient;
step 11, switching into a VSG mode, namely setting the operation mode as 1, and recording the operation of an inverter in the mode as a new period;
sampling output phase voltage U of new cycle of inverter oa ',U ob ',U oc ', sampling the inverter side filter inductance phase current I of the new cycle of the inverter inva ',I invb ',I invc ', sampling the output phase current I of the new cycle of the inverter oa ′,I ob ′,I oc ′;
Step 12, according to the conversion formula 2 from the output three-phase voltage to the two-phase static coordinate system, the output phase voltage U of the new cycle of the inverter sampled in the step 11 is obtained oa ',U ob ',U oc ' converting to obtain the output phase voltage two-phase stationary coordinate system component U of the new cycle of the inverter ',U '; converting the output three-phase current to a two-phase static coordinate system according to a conversion formula 2, and outputting the output phase current I of the inverter in the new period sampled in the step 11 oa ',I ob ',I oc ' conversion to obtain the output phase current two-phase stationary coordinate system component I of the new cycle of the inverter ',I '; secondly, obtaining the output active power P 'of the inverter in a new period and the output reactive power Q' of the inverter in the new period through a power calculation formula 2; then, the VSG angular frequency omega of a new period is obtained through a VSG reactive power control formula 2 vsm ', VSG active power control formula 2 obtains d-axis voltage amplitude reference U of new period dref ';
The conversion formula 2 from the output three-phase voltage to the two-phase static coordinate system is as follows:
Figure FDA0003784820060000061
Figure FDA0003784820060000062
the conversion formula 2 from the output three-phase current to the two-phase stationary coordinate system is as follows:
Figure FDA0003784820060000071
Figure FDA0003784820060000072
the power calculation formula 2 is:
P′=U ′I ′+U ′I
Q'=U 'I '-U 'I '
the active power control formula 2 is:
Figure FDA0003784820060000073
the reactive power control formula 2 is:
Figure FDA0003784820060000074
step 13, obtaining the angular frequency omega of the voltage source in the new period according to the step 12 csm ', obtaining an output voltage angle theta' of a new period of the inverter through a mode selection integral formula 2;
the mode selection integral equation 2 is:
Figure FDA0003784820060000075
step 14, according to the output voltage angle θ' of the inverter in the new period obtained in step 13, the inverter side filter inductor phase current I of the inverter in the new period obtained in step 11 is sampled by the voltage source mode inverter side phase current rotating coordinate transformation equation 2 inva ',I invb ',I invc ' coordinate transformation is carried out to obtain the feedback I of the side filter inductive current d axis of the voltage source mode inverter of a new period of the inverter vsm_dfeedback ' with new period of the voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ' the output phase voltage U of the new cycle of the inverter sampled in the step 11 is converted into the output phase voltage U of the new cycle of the inverter by the voltage source mode inverter side phase voltage rotating coordinate transformation equation 2 oa ',U ob ',U oc ' coordinate transformation is carried out to obtain d-axis feedback U of output voltage of the voltage source mode inverter in a new period of the inverter vsm_dfeedback ' and new period of q-axis feedback U of output voltage of voltage source mode inverter vsm_qfeedback ';
The voltage source mode inverter side phase current rotation coordinate transformation equation 2 is as follows:
Figure FDA0003784820060000081
Figure FDA0003784820060000082
the voltage source mode inverter side phase voltage rotating coordinate transformation equation 2 is as follows:
Figure FDA0003784820060000083
Figure FDA0003784820060000084
step 15, according to the filter inductance current d-axis feedback I of the voltage source mode inverter side of the new cycle of the inverter obtained in the step 14 vsm_dfeedback ' and new cycle voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ' obtaining the filter inductance current d-axis feedback I at the side of the current source mode inverter of a new period of the inverter through a current source mode inverter side phase current rotation coordinate transformation equation 2 csm_dfeedback ' and new periodic current source mode inverter side filter inductor current q-axis feedback I csm_qfeedback ';
The current source mode inverter side phase current rotation coordinate transformation equation 2 is:
I csm_dfeedback ′=I vsm_dfeedback
I csm_qfeedback ′=I vsm_qfeedback
step 16, according to the d-axis feedback U of the output voltage of the voltage source mode inverter of the new period of the inverter obtained in the step 14 vsm_dfeedback ' and new period of q-axis feedback U of output voltage of voltage source mode inverter vsm_qfeedback ' and step 12. The new cycle of U drefv ' obtaining a new cycle of d-axis reference I of the voltage source mode current loop through the voltage source mode voltage loop equation 2 vsm_dref ' and new periodic Voltage Source mode Current Loop q-axis reference I vsm_qref ';
The voltage source mode voltage loop equation 2 is:
I vsm_dref '=(U drefv '-U vsm_dfeedback ')×(K vsm_vp +K vsm_vi /s)
I vsm_qref '=(0-U vsm_qfeedback ')×(K vsm_vp +K vsm_vi /s)
step 17, referring to I according to the d-axis of the new periodic voltage source mode current loop obtained in step 16 vsm_dref ', new cycle voltage source mode electricityFlow ring q-axis reference I vsm_qref ', new cycle voltage source mode inverter side filter inductor current d-axis feedback I for the new cycle obtained in step 14 vsm_dfeedback ' and new cycle voltage source mode inverter side filter inductor current q-axis feedback I vsm_qfeedback ' and current source mode modulation d-axis reference e obtained in step 5 csm_dref Current source mode modulated q-axis reference e csm_qref Obtaining a new period of voltage source mode modulation d-axis reference e through a voltage source mode current loop equation 2 vsm_dref ' modulating the q-axis reference e with the new periodic voltage source pattern vsm_qref ';
The voltage source mode current loop equation 2 is:
e vsm_dref '=(I vsm_dref '-I vsm_dfeedback ')×K vsm_ip +e csm_dref
e vsm_qref '=(I vsm_qref '-I vsm_qfeedback ')×K vsm_ip +e csm_qref
step 18, obtaining a new period of voltage source mode modulation d-axis reference e according to step 17 vsm_dref ' modulating the q-axis reference e with the new periodic voltage source pattern vsm_qrefv ' obtaining a modulation wave e of a new period in the dq coordinate system by selecting equation 2 through the modulation mode d ' and e q ', then modulating the wave to pass the inverter through the grid impedance L g Outputting electric energy to a power grid;
the modulation mode selection equation 2 is:
e d ′=e vsm_dref
e q ′=e vsm_qref
step 19, according to the current source mode inverter side filter inductor current d-axis feedback I of the new period obtained in the step 15 csm_dfeedback ', new period current source mode inverter side filter inductance current q axis feedback I csm_qfeedback ', externally given current source d-axis reference I csm_dref Current source q-axis reference I csm_qref And current source mode modulation d-axis reference e obtained in step 5 csm_dref And a current sourceMode-modulated q-axis reference e csm_qref Obtaining a new cycle of current source mode modulation d-axis reference e through a current source mode current loop equation 2 csm_dref ' modulating the q-axis reference e with a new periodic current source pattern csm_qref ';
The current loop equation 2 for the current source mode is:
Figure FDA0003784820060000101
2. the method according to claim 1, wherein the three-phase inverter output voltage angle θ ″ of the previous cycle in step 2 is initialized to 0 in the first cycle of the operation of the control method.
CN201911308959.0A 2019-12-18 2019-12-18 Grid-connected inverter VSG mode and current source mode undisturbed switching control method Active CN111030174B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911308959.0A CN111030174B (en) 2019-12-18 2019-12-18 Grid-connected inverter VSG mode and current source mode undisturbed switching control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911308959.0A CN111030174B (en) 2019-12-18 2019-12-18 Grid-connected inverter VSG mode and current source mode undisturbed switching control method

Publications (2)

Publication Number Publication Date
CN111030174A CN111030174A (en) 2020-04-17
CN111030174B true CN111030174B (en) 2022-10-04

Family

ID=70210414

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911308959.0A Active CN111030174B (en) 2019-12-18 2019-12-18 Grid-connected inverter VSG mode and current source mode undisturbed switching control method

Country Status (1)

Country Link
CN (1) CN111030174B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI746097B (en) 2020-07-31 2021-11-11 台達電子工業股份有限公司 Power conversion device with dual-mode control
CN113595137B (en) * 2021-07-20 2023-04-07 科华数据股份有限公司 Direct current source grid-connected power generation mode switching method and device and direct current source system
CN113765140A (en) * 2021-08-25 2021-12-07 国电南瑞科技股份有限公司 Operation mode control method, device and system suitable for current source type photovoltaic grid-connected inverter
CN114189007B (en) * 2021-11-15 2024-04-12 阳光电源股份有限公司 Vehicle-mounted charger working mode switching control method and device and vehicle-mounted charger
CN117277421B (en) * 2023-11-17 2024-02-02 合肥工业大学 Multi-model control method for grid-connected inverter

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108400603A (en) * 2018-01-19 2018-08-14 南京邮电大学 A kind of micro-capacitance sensor that idle delay based on VSG is adjusted is without frequency difference method for handover control
CN110098637A (en) * 2019-05-23 2019-08-06 合肥工业大学 Cascaded H-bridges photovoltaic combining inverter double mode smooth sliding control method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2634909B1 (en) * 2012-03-02 2017-02-15 ABB Research Ltd. Method for controlling a grid-connected boost-buck full-bridge current-source inverter cascade for photovoltaic applications and device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108400603A (en) * 2018-01-19 2018-08-14 南京邮电大学 A kind of micro-capacitance sensor that idle delay based on VSG is adjusted is without frequency difference method for handover control
CN110098637A (en) * 2019-05-23 2019-08-06 合肥工业大学 Cascaded H-bridges photovoltaic combining inverter double mode smooth sliding control method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
微电网逆变器VCS模式与CCS模式的切换技术;梁建钢 等;《电网技术》;20140430;第830-837页 *

Also Published As

Publication number Publication date
CN111030174A (en) 2020-04-17

Similar Documents

Publication Publication Date Title
CN111030174B (en) Grid-connected inverter VSG mode and current source mode undisturbed switching control method
CN106849184B (en) Asynchronous Wind turbines virtual synchronous grid-connected control method and system
CN108233415B (en) Control method for two-stage photovoltaic inverter virtual synchronous generator
CN111030139B (en) Series compensation power grid resonance suppression method based on virtual synchronous generator
CN112290561B (en) LC resonance suppression method of virtual synchronous machine based on virtual parallel resistor
CN106159999B (en) Bavin storage mixing independent micro-grid dynamic frequency stable strategy based on virtual synchronous generator
CN108847669B (en) Multi-synchronous rotation coordinate system-based multifunctional grid-connected inverter harmonic treatment method
CN101702583B (en) Method for controlling direct-drive wind power generation convertor
CN107196344B (en) Self-synchronizing virtual synchronous inverter grid-connected controller and method with local load based on SPF-PLL
CN112234654B (en) LC resonance suppression method of virtual synchronous machine based on capacitance voltage differential feedback
CN110021959B (en) Grid-connected inverter dual-mode control method based on short-circuit ratio under weak power grid
CN107257141B (en) Utilize the self synchronous three-phase grid-connected converter control method of DC capacitor dynamic implement
CN107069828B (en) Virtual synchronous generator self-synchronization control method based on real-time phase difference adjustment
CN110838731A (en) Three-phase four-bridge arm photovoltaic grid-connected inverter and multi-target cooperative fault control method
CN113991755A (en) Self-synchronizing voltage source control method for new energy power generation unit
CN113765140A (en) Operation mode control method, device and system suitable for current source type photovoltaic grid-connected inverter
WO2022022202A1 (en) Inductance current differential feedback-based virtual synchronous machine strong power grid stability control method
Ge et al. Inverter control based on virtual impedance under unbalanced load
CN107623338A (en) The independent excitation control method of three-phase four-arm virtual synchronous generator
CN113629763B (en) Current control method and system for medium-voltage direct-hanging energy storage converter under non-ideal power grid
CN113937789A (en) Voltage source type double-fed fan feedforward damping control method based on fractional order filtering
Raju et al. Enhanced decoupled power control of wind turbine driven DFIG using DVR under unbalanced grid voltage
Singh et al. Power control in centralized distributed AC load for wind energy system
Pang et al. Voltage Harmonics Optimization for Weak Grid-Tied Doubly-Fed Induction Generator With the Capability of Suppressing Current Harmonics
CN110557037A (en) direct current control method of single-phase PWM rectifier

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