The present disclosure relates to inverters and, more specifically, relates to a synchronverter power control during unbalanced grid conditions.
Standard inverters are known to be having low inertia. If there are faults or sudden changes in load leading to transient periods, the standard inverters follow changes rapidly. As a result, grids taking support of the standard inverters may experience a worser condition. To overcome the deficiencies of standard inverters, Synchronverters (SVs) are introduced. The SVs are inverters with defined control strategy to behave as synchronous generators (SGs) to avoid traditional inverters' unwanted low inertia behaviors. The SVs may operate in a connected mode or in an islanded mode from a hosting power grid. However, generic SVs are not capable of providing safe and reliable operation during unbalanced faults or grid conditions. Consequently, the SVs may generate currents that may exceed their nominal values, inject active and reactive power with oscillations at twice the grid's fundamental frequency.
In one aspect of the present disclosure, a system for synchronverter power control during unbalanced grid conditions is disclosed. Methods and systems for a synchronverter power control during unbalanced grid conditions is disclosed. The system includes a synchronverter coupled with a power supply grid, a power reference generator, configured to receive a terminal voltage measurement vector vt and a current measurement vector i from the synchronverter, and generate an active power Pf and a reactive power Qf, a synchronverter control unit connected to the power reference generator and configured to process the active power Pf and the reactive power Qf and generate an electromotive force (EMF) vector e, and an active and reactive power control unit, connected between the synchronverter control unit and the synchronverter, configured to receive the electromotive force (EMF) vector e and the terminal voltage measurement vector vt, and regulate the current measurement vector i to eliminate power oscillations and current harmonics in the synchronverter during unbalanced grid conditions.
In another aspect, a method for synchronverter power control is disclosed. The method includes electrically coupling a synchronverter coupled with a power supply grid, generating an active power Pf and a reactive power Qf through a power reference generator based on a terminal voltage measurement vector vt and a current measurement vector i from the synchronverter, generating an electromotive force (EMF) vector e through a synchronverter control unit by processing an active power Pf and a reactive power Qf, and regulating the current measurement vector i to eliminate power oscillations and current harmonics in the synchronverter during unbalanced grid conditions through an active and reactive power control unit, based on the electromotive force (EMF) vector e and the terminal voltage measurement vector vt.
The foregoing as well as other features and advantages of the present disclosure will be more fully understood from the following description, examples, and claims.
A more complete appreciation of this disclosure and many of the attendant advantages thereof will be readily obtained as the same becomes better understood by reference to the following detailed description when considered in connection with the accompanying drawings, wherein:
It should be appreciated by those skilled in the art that any diagram herein represents conceptual views of illustrative systems embodying the principles of the present disclosure.
Reference will now be made in detail to specific embodiments or features, examples of which are illustrated in the accompanying drawings. A skilled artisan will appreciate that various alternate embodiments and forms may be prepared. Examples, therefore, given are only for illustration purposes without any intention to restrict the embodiments to a given set of examples. Specific functional aspects are provided merely to enable a person skilled in the art to perform the invention and should not be construed as limitations of the invention. Any method steps, and processes described herein are not to be construed as necessarily requiring their performance in the particular order discussed or illustrated, unless specifically identified as an order of performance. It is also to be understood that additional or alternative steps may be employed, unless otherwise indicated.
The use of the singular herein includes the plural (and vice versa) unless specifically stated otherwise. The use of the terms “include,” “includes”, “including,” “have,” “has,” or “having,” “comprise,” “comprises,” “comprising” or the like should be generally understood as open-ended and non-limiting unless specifically stated otherwise. It is understood that the order of steps or order for performing certain actions can be changed so long as the intended result is obtained. Moreover, two or more steps or actions may be conducted simultaneously. As used herein, the term “about” or “between” refers to a ±20% to ±10% variation from the nominal value unless otherwise indicated.
Embodiments of the present disclosure are directed to methods and systems for synchronverter power control during unbalanced grid conditions. The disclosure may equip the Synchronverters (SVs) with fault ride-through (FRT) units to avoid any adverse impact of unbalanced grid conditions and maintain a reliable SV operation. The FRT units are configured to eliminate active and reactive power oscillations problems and ensure current generation is within the defined inverter ratings. The FRT units maintain SVs' intrinsic features, guarantee a seamless transition during faults, and facilitate SVs integration in microgrids.
E=Mfif{dot over (θ)}; (1)
where Mf if is the virtual mutual field multiplied by the virtual field current. The Mf if is determined by following equation:
where Vt is vt magnitude, K is a reactive power regulating coefficient, Qset is a setpoint for a reactive power, Q is a produced reactive power by the SV, Dq is a voltage droop coefficient that regulates Q generation for a specific change in Vt, and E* is a voltage magnitude setpoint. The produced Q and Dq are provided by:
On the other hand, SV angular frequency {dot over (θ)} is regulated by the following equation:
where Pset is an active power setpoint, {dot over (θ)}* is the angular frequency setpoint, J is a virtual inertia, Te is an electromagnetic torque, and Dp is a frequency droop coefficient. Dq, Dp stimulates the SV to generate a certain amount of Te denoted as ΔTe and thus active power for each specific change in {dot over (θ)} denoted as Δ{dot over (θ)} as shown below as:
where θ is a balanced three-phase sinusoidal. The e vector is represented as:
e=E θ (8)
The SV control unit as illustrated in
where θt, Zs, and θs are a terminal voltage phase, a filter impedance magnitude, and an angle, respectively. From equation (9), it is apparent that the current depends on the voltage difference between SV and the measured Vta∠θt as well as the filter impedance Zs∠θs. Thus, the current may exceed the inverter's ratings and an interfacing filter during unbalanced conditions, which raises a requirement for limiting the current to avoid any damages.
According to instantaneous power theory, SVs generated P and Q are expressed as
P=vt·i
Q=vt⊥·i (10)
where the operator (.) may represent the dot product of vectors, and the subscript (⊥) may denote the orthogonal version of vt. Both vt and i may be redefined by applying symmetrical component theory as given by:
vt=vt+vt−
i=i++i− (11)
The superscripts (+) and (−) refer to positive and negative sequence components, respectively. By substituting (11) in (10), P and Q are rewritten as:
P=vt+·i++vt−·i−+vt+·i−+vt−·i+ (12)
Q=vt
Each of (12) and (13) include two terms, a constant term that is a result of vt and i interaction from a same sequence and an oscillatory term results from the v and i interaction from different sequences. The terms P and Q are defined as:
Ps=vt+·i++vt−·i− (14)
{tilde over (P)}=vt+·i−+vt−·i+ (15)
P=Ps+{tilde over (P)} (16)
Qs=vt⊥+·i++vt⊥−·i− (17)
{tilde over (Q)}=vt⊥+·i−+vt⊥−·i+ (18)
Q=Qs+{tilde over (Q)} (19)
where Ps is a constant term of P which is usually equal to the active power setpoint, i.e., Pset. The Qs is a constant term of Q, which is equal to the reactive power setpoint Qset. The oscillatory terms of P and Q are {tilde over (P)} and {tilde over (Q)}, respectively, having a zero-mean value. During unbalanced conditions, the {tilde over (P)} and {tilde over (Q)} are not null since the negative sequence components of vt and i exist. As a result, it may be required to mitigate the effect of unbalanced grid conditions by reducing/eliminating {tilde over (P)} and {tilde over (Q)}.
The problem of power oscillations and overcurrent generation is resolved for SVs as described henceforth. To significantly eliminate {tilde over (P)} and {tilde over (Q)}, an instantaneous active and reactive power control (IARC) technique may be used. If i aligns with the vt leads to P generation. Conversely, the Q generation is performed when the i aligns with the vt⊥. This concept is used to form the reference current that is essential to achieve IARC as following:
i*=ip*+iq* (20)
where i* is the three-phase reference current that performs the IARC function. ip* is a three-phase component of i that is responsible for P generation. Similarly, ig* is responsible for Q generation. ip*, and iq* are expressed as:
ip*=g*vt (21)
ig*=b*vt⊥ (22)
where g and b are the instantaneous conductance and susceptance, respectively. For SVs, g and b can be defined as:
g=Ps/|vt|2 (23)
b=Qs/|vt|2 (24)
|vt|2=|vt+|2+|vt−|2+2|vt+∥vt−|cos(2ωt+ϕ++ϕ−) (25)
where ϕ+ and ϕ− are the phases of vt+ and vt−, respectively. Based on the equations (20)-(24), the reference current that achieves IARC for SVs can be written as
From the equation (26), i* depends on Ps and Qs that are generated internally in the SV unit, which is explained in the methodology section. Also, it is well known that the IARC has a trade-off between the oscillations and the generated current quality due to the existence of a cosine term in the denominator. Allowing the complete existence of the cosine term gives the highest degree of eliminating the oscillations. However, the corresponding generated current includes higher-order harmonic components, mainly the third harmonic component. This third harmonic component may be avoided by eliminating the cosine term, which leads to an average active and reactive power control (AARC) methods that ensures the delivery of Ps and Qs equal to their corresponding setpoints, i.e., Pset and Qset, respectively. ip* and iq* for AARC are defined as:
ip*=G*vt (27)
iq*=B*vt⊥ (28)
where G and B are the average value of the conductance and susceptance, respectively. The G and B are constant during unbalanced grid conditions since they do not exhibit oscillations as their instantaneous counterparts g and b. The G and B are given by:
G=Ps/vtΣ2 (29)
B=Qs/vtΣ2 (30)
vtΣ=√{square root over (|vt+|2+|vt−|2)} (31)
With vtΣ, which is the collective RMS value of vt, the AARC reference current is expressed as:
The reference current provided in the equation (32) prevents the generation of higher-order components of the current; thus, resulting in a sinusoidal current generation. Also, the reference current in (32) is identical to (26) if the cosine term is removed.
The IARC and AARC may be concluded as a relation between the power oscillations that adversely affect the current quality. As described, the higher the oscillations are eliminated, the higher the current is distorted. Therefore, a comprehensive IARC (C-IARC) is described in the disclosure to eliminate the oscillations without affecting the current or its quality. The C-IARC as described herein is developed in an αβ frame. vt can be represented in the αβ frame as follows:
where vta, vtb, and vtc, are the abc quantities of vt, respectively. vα and vβ are the αβ frame quantities of vt, respectively. C-IARC involves the modification of the instantaneous power theory to avoid the drawbacks of IARC and AARC. The following expressions are express original instantaneous power theory and the modified instantaneous power theory, the active-power P and the reactive-power Q are defined by:
Equation (35) is the modified instantaneous power theory to estimate P and Q. Term {circumflex over (v)}a is a π/2 delayed version of vα while {circumflex over (v)}β is a π/2 advanced version of vβ. The aforementioned delay can be achieved by using a compensation filter with a transfer function defines as:
The quantity ωd represents the delayed quantity frequency and ωc is the filter cut-off frequency. ξ is the filter damping, which can be set as 1 to ensure good filter damping. The Fd(s) delays vα by π/2 and the same Fd(s) can advanced {circumflex over (v)}β by multiplying by (−1) as following:
{circumflex over (v)}α=vaFd(S) (37)
{circumflex over (v)}62=−vβFd(S) (38)
Consequently, the reference current αβ frame quantities {circumflex over (ι)}α* and {circumflex over (ι)}β* to perform the C-IARC can be derived as:
The {circumflex over (ι)}α(p)* and {circumflex over (ι)}β(p)* are responsible for generating Ps while {circumflex over (ι)}α(q)* and {circumflex over (ι)}β(q)* are on the other hand, responsible for generating Qs. Each sub-component of {circumflex over (ι)}α* and {circumflex over (ι)}β* based in (39) are formulated as
The reference currents described in equations (26), (32), and (40) achieve different FRT strategies. The C-IARC outperforms both the IARC and AARC. However, all of the aforementioned strategies have a common drawback; the unlimited current generation may lead to inverter damages. This damages may be prevented by having a current limiter. The current limiter may be designed by modifying Pset and Qset to a form that ensures active and reactive power delivery during normal conditions and limited current generation during unbalanced grid conditions. As described earlier the Ps and the Qs are equal to the Pset and the Qset during normal and abnormal conditions. During unbalanced conditions, SVs generate high currents to compensate for the voltage drop. This high current ensures the generated P and Q are equal to Pset and Qset. Therefore, the inverter is subjected to damages due to the high current generation. In conclusion, changing Pset and Qset to lower values may reduce the values of Ps and Qs lower, and thus the generated current may be lower since the Ps and Qs tracks the reference current of Equations (26), or (32), or (40), which all are affected by the value of Ps and Qs. The new active power setpoint Pf and reactive power setpoint Qf are derived as:
where Vd, Id, and Iq are the d axis vt component, and dq axis components of i, respectively. During faults, the quantity Vd is represented as:
Vd=(V+−V−) (42)
Also, Id, and Iq may be defined as a function of a maximum peak current Imax as:
Id=Imax cos(θi)
Iq=Imax sin(θi) (43)
where θi is the angle of a required power factor. This new set of power references ensures the current does not exceed the pre-defined limit Imax as demonstrated in the results sections.
The IARC, AARC, and C-IARC may be added in series to the generic SVs control unit. This configuration allows SVs to maintain their intrinsic features and guarantee a seamless transition without any switching to achieve a fault ride-through (FRT). As explained, IARC and AARC may be formulated and designed jointly.
Further, a maximum peak current reference input 312 is provided to a d axis vt component generator 314 that generates Id, and Iq based on Equation (43). The Id is provided to a multiplier 310 that multiplies Id with product of the Vd and 1.5 to generate Pf as provided in Equation (41). Further, Iq is multiplied with −1 by a multiplier 316 and provided to a multiplier 318. The multiplier 318 multiplies −Iq with the product of the Vd and 1.5 to generate Qf as provided in Equation (41).
The power references replacement ensures the required power delivery during normal conditions is identical if Pset and Qset are used. Also, Pf and Qf maintain the peak current during faults below a specific limit. Thus, Pf and Qf may lead to a safe current generation and a reliable operation during normal and unbalanced conditions.
For a traditional converter, IARC or AARC may be achieved by replacing Ps and Qs in equations (26) and (32) by the active power and reactive power setpoints that deliver a required amount of power. For the SVs, Pset and Qset are the power reference setpoints. However, it may not be possible to replace Ps and Qs by Pset and Qset since it leads to cancelation of the SV generic control unit. On the other hand, e may depend on Pset and Qset as described in equations (1)-(8). Thus, it may be possible to generate Ps and Qs from e by generating a current is that is obtained from e and vt difference divided by the filter impedance as shown in
Ps=vt+·is++vt−·is− (44)
Qs=vt⊥+·is++vt⊥−·is− (45)
As seen in (44) and (45), and also in
In an example implementation, the current is may be generated based on e and vt through an adder 352 that combines e and subtracts vt which is then provided to a component 354 to generate the current is. The current is is input to the sequence extractor 356 to generate is+ and is−. The is+ and is− are input to an internal power reference generator 364. Also using vt, the sequence extractor 358 generates vt+ and vt−, that are input to the orthogonal creator 360 and the orthogonal creator 362 to generate vt⊥+ and vt⊥−, respectively. The vt⊥+ and vt⊥− are provided to the internal power reference generator 364. The internal power reference generator 364 generates Ps and Qs as described in Equations (44) and (45). The Ps and Qs are input to an IARC and AARC reference current generator unit 372. Since the Ps and Qs are generic for both IARC and AARC, and switching between these FRT strategies is related to the cosine term, i* can be rewritten to switch flexibly between IARC and AARC as:
where D=|vt+|2+|vt−|2+2k|vt+∥vt−|cos(2ωt+ϕ++ϕ−) (49)
The parameter k in the denominator D determines whether i* achieves IARC or AARC. Setting the parameter k to 1, the denominator is equal to the one in Equation (26), and hence achieving an IARC method. On the other hand, setting k to 0 makes the denominator as in (32), which consequently leads to an AARC method. Thus, D is defined as:
The IARC and AARC reference current generator unit 372 generates i* based on Ps and Qs, vt, vt⊥ (received from the orthogonal creator 366), and D (received from component 370 based on vt and k input 368). For the control purpose design, creating D when k=1 is relatively complicated since it involves phase estimation. Thus, a more straightforward approach is followed by using D when k=1 as
D=vta2+vtb2+vtc2 (51)
where vta, vtb, and vtc, are the abc quantities of vt. The generated i* from IARC or IARC is tracked by i through subtractor 374 and proportional-resonant (PR) controller 376 along with feedforward vt to obtain through adder 378 the new EMF vector e* to create PWM signals through PWM 382.
The implementation of the proposed C-IARC is done similarly to the IARC and AARC implementation.
In an example implementation, the current is may be generated based on e and vt through an adder 452 that combines e and subtracts vt which is then provided to a component 454 to generate the current is. The current is is input to a sequence extractor 456 to generate is+ and is−. The is+ and is− are input to an internal power reference generator 464. Also using vt, the sequence extractor 458 generates vt+ and vt−, that are input to an orthogonal creator 460 and the orthogonal creator 462 to generate vt⊥+ and vt⊥−, respectively. The vt⊥+ and vt⊥− are provided to an internal power reference generator 464. The internal power reference generator 464 generates Ps and Qs as described in Equations (44) and (45). The Ps and Qs are input to a reference current generator unit 472.
The modifications made in
Simulation results are described henceforth. In an example, using a MATLAB/SIMULINK, a simulation model is built to simulate the SV topology in
The simulated SV may be subjected to an unbalanced grid condition at its terminal.
The generic SV simulation model described above is modified to include the IARC and AARC units to resolve a problem of the power oscillations demonstrated in
The results of the IARC and AARC show the advantages and disadvantages for both. C-IARC merges the advantages of IARC and AARC while discarding their disadvantages, that is, power oscillations and distorted current generation.
The disclosure provides a seamless FRT strategies. The addition of the FRT to the SV control unit does not affect the normal operation of the SV control unit. Thus, the SV intrinsic features are unaffected. The current limitation of the disclosure may be suitable for a generic SV as well as SV with the FRT strategies. The activation of the current limitation and any proposed FRTs strategies is a seamless transition process without requiring any switching. Thus, the SV stability may not be affected by any switching action. The IARC and AARC are flexibly coupled to achieve power oscillation reduction or current harmonic elimination. The C-IARC eliminates the power oscillations and current harmonics simultaneously. Further, there is no requirement to control the sequence components and add complexity to the SV control unit with the FRT strategies. Further, no hardware modification may be required for the current SV control unit.
It is understood that the examples, embodiments and teachings presented in this application are described merely for illustrative purposes. Any variations or modifications thereof are to be included within the scope of the present application as discussed.
Number | Name | Date | Kind |
---|---|---|---|
10615716 | Zhong | Apr 2020 | B2 |
20140067138 | Rodriguez Cortes | Mar 2014 | A1 |
20140152292 | Andresen | Jun 2014 | A1 |
20190190276 | Liu | Jun 2019 | A1 |
20190245458 | Wang | Aug 2019 | A1 |
20210249862 | Awal | Aug 2021 | A1 |
Number | Date | Country |
---|---|---|
108808704 | Jul 2021 | CN |
Entry |
---|
J. O. M. B. de Almeida, et al., “Three-phase photovoltaic inverters during unbalanced voltage sags: Comparison of control strategies and thermal stress analysis,” 2016 12th IEEE International Conference on Industry Applications, Curitiba, PR, Brazil, 2016, pp. 1-7, doi: 10.1109/INDUSCON.2016.7874564. (Year: 2016). |
M. Zhu, L. Hang, G. Li and X. Jiang, “Protected Control Method for Power Conversion Interface Under Unbalanced Operating Conditions in AC/DC Hybrid Distributed Grid,” in IEEE Transactions on Energy Conversion, vol. 31, No. 1, pp. 57-68, Mar. 2016, doi: 10.1109/TEC.2015.2480757. (Year: 2016). |
S. F. Zarei, H. Mokhtari, M. A. Ghasemi, S. Peyghami, P. Davari and F. Blaabjerg, “Control of Grid-Following Inverters Under Unbalanced Grid Conditions,” in IEEE Transactions on Energy Conversion, vol. 35, No. 1, pp. 184-192, Mar. 2020, doi: 10.1109/TEC.2019.2945699. (Year: 2020). |
Prasanna Piya “Grid fault ride-through capability of voltage-controlled inverters for distributed generation” May 5, 2017. |
Weiyi Zhang “Control of grid connected power converters with grid support functionalities” May 2017. |