1. Introduction
In this context, this paper presents an analytical model of the VSC under unbalanced conditions for the direct computation of the non-characteristic third harmonics through rigorous mathematical derivations. The new model reveals the mechanism of the generation of third harmonics and the characteristics of third harmonics. The new model shows that a Thevenin circuit can represent the VSC for both the positive-sequence and negative-sequence third harmonics. In addition, there are coupling effects between the positive-sequence and negative-sequence third harmonics. The model also shows that the positive-sequence third harmonic is much larger than the negative-sequence third harmonic. The computation of the source and the impedance of the developed Thevenin model requires positive-sequence and negative-sequence fundamental frequency voltage and current on the ac side of the VSC, which can be obtained through load flow studies.
Power electronic devices have experienced a dramatic increase in modern power systems due to the initiative of renewable energies such as wind power generation and photovoltaic (PV) [
1], and other applications such as flexible AC transmission system (FACTS) and high voltage DC (HVDC) transmission [
2,
3]. In recent years, the voltage source converter (VSC) has become the most popular basis of power electronic devices due to its advantages of bidirectional power flow and independent control of active and reactive power [
4]. With the massive integration of these VSC devices, there is a concern about their harmonic impact on power systems [
5,
6,
7]. Under the normal condition, VSC devices only produce very-high-frequency harmonics, which is produced by the PWM switching, but such harmonics can generally be impeded by a passive filter in front of the VSC [
8,
9]. However, VSC devices can generate low-order non-characteristic third harmonics under the unbalanced condition, which is a power quality concern of power systems. Such a phenomenon occurs because the unbalanced fundamental frequency voltage and current on the ac side of the VSC can induce second harmonic ripples in the dc-link voltage of the VSC through the switching process [
9]. These second harmonic ripples in the dc-link voltage, in turn, will lead to a third harmonic voltage on the ac side of the VSC, again through the switching process. As a result, under the unbalanced condition, the VSC becomes a voltage source at the third harmonic. Therefore, it is essential to understand the mechanism of the production of the third harmonic voltage, especially in the distribution system, which can be unbalanced since single-phase loads are dominant [
10]. A computation method of the third harmonics under the unbalanced conditions is also needed to assess their impacts.
To date, there are numerous papers presenting the studies of the VSC under unbalanced conditions, but most of these works focus on the control strategies that can suppress the unbalanced condition-induced third harmonics [
11,
12,
13,
14]. These papers have shown that the unbalanced fundamental frequency voltage and current will lead to double-frequency power oscillation on the ac side of the VSC. As a result, the dc side experiences second harmonic voltage, which can lead to third harmonics on the ac side of the VSC. Based on this finding, the control schemes that can be used to suppress the second harmonic on the dc-link are designed in [
12,
13]. However, these papers do not provide an analytical algorithm to compute the second harmonic on the dc side or the third harmonic on the ac side. A few other works are proposed specifically for assessing VSC’s third harmonic emission under unbalanced conditions [
15,
16,
17]. In [
15], the switching function of the VSC and the harmonic interaction between the ac side and dc side of the VSC are analyzed to obtain the model that can be used to compute the third harmonic distortion. Nevertheless, the impact of the third harmonic on the switching function is not included in this method. Thus, the result could be inaccurate. In [
16], the oscillation between the dc side and ac side of the VSC is analyzed based on the ac/dc power interaction, which reveals the mechanism of the production of the third harmonic, but the final analytical method of the uncharacteristic third harmonics is not presented. The reference [
17] further provides a method for assessing the third harmonic distortions using the power interaction between the ac side and dc side of the VSC. However, the proposed model omits the effect of the control loop so that the model is not accurate. Moreover, this method is not able to fully show the differences between the positive-sequence and negative-sequence third harmonics and the coupling effect between these two kinds of third harmonics. References [
18,
19] propose a method that can be used to reveal the harmonic characteristics of the VSC by analyzing the control and AC/DC interaction, but the non-characteristic third harmonics are not analyzed. In summary, a full model that reveals the characteristics of third harmonics produced by VSC is still desired.
This paper presents an analytical model of the VSC under unbalanced conditions for the direct computation of the non-characteristic third harmonics through rigorous mathematical derivations. The development of the new model reveals the mechanism of the generation of third harmonics and their characteristics. The new model shows that a Thevenin circuit can represent the VSC at both the positive-sequence and negative-sequence third harmonics. In addition, there are coupling effects between the positive-sequence and negative-sequence third harmonics. The model also shows that the positive-sequence third harmonic is much larger than the negative-sequence third harmonic. The computation of the source and the impedance of the developed Thevenin model requires a positive-sequence and negative-sequence fundamental frequency voltage and current on the ac side of the VSC, which can be obtained through load flow studies.
The rest of this paper is organized as follows.
Section 2 explains the mechanism of generating the non-characteristic third harmonics by the VSC under unbalanced conditions.
Section 3 presents the development of the model of VSC under unbalanced conditions for a non-characteristic third harmonic. The proposed model is validated via time-domain simulations in
Section 4.
Section 5 shows some discussions on the proposed model.
2. Mechanism of the Production of Noncharacteristic Third Harmonics
A typical VSC system is illustrated in
Figure 1a, and the control strategy of this VSC is presented in
Figure 1b. As one can see, a two-loop control scheme is used, and the outer control loop is used to generate the reference current for the inner current control. The current is controlled using a PI regulator. A capacitor contains the dc-link voltage of the VSC. An inductor
Lf is connected at the front end of the VSC to filter out the switching harmonics.
Under the unbalanced operational condition, there are positive-sequence and negative-sequence fundamental frequency voltage and current on the AC side of the VSC, and the three-phase current can be expressed as
It is noted that the modulation signal is a result of the inner current control, so the modulation signal contains the same types of components as those in the current. Thus, the three-phase modulation signals can be expressed as
Since we focus on the study of the third harmonic, the high-frequency harmonics in the switching function can be omitted. Consequently, the modulation signal can be used to represent the switching function of the VSC. As a result, based on (1) and (2), the current on the DC side can be computed as
According to the above DC side current, the voltage on the DC-link capacitor can be written as
where
Vdc0 is the steady-state dc-link voltage. It can be seen that, due to the unbalanced condition, there are second harmonic ripples in the dc-link voltage. It is noteworthy that the expression in (4) is obtained on the basis that all dc-side harmonic currents pass through the dc-link capacitor. In practice, there will be a load (or source) on the other side of the dc-link that has an equivalent impedance
Zload, and the harmonic components in
Idc will be distributed between the load and dc-link capacitor (
Idc1 and
Idc2 in
Figure 1a), that is, all of the harmonic components in (4) will be multiplied by
Zload/(
Zload +
ZC). For simplification but not losing generality,
Zload is assumed to be very large so it can be ignored. The DC-link voltage in (4) is converted to AC voltage through the switching process, obtaining the expressions of the three-phase voltages, which are expressed in (5)~(7).
As one can see from (5)~(7), there are third harmonic components in the three-phase voltage of the VSC. In addition, it is found that the third harmonic voltages contain both positive-sequence or negative-sequence components, which are different from the conventional zero-sequence third harmonic. It can also be noted that the third harmonic is in the form of a voltage that is only affected by the positive-sequence and negative-sequence fundamental frequency components.
3. Modeling VSC under Unbalanced Condition
The above analysis has shown that an unbalanced condition can introduce a third harmonic voltage for the VSC. Such a third harmonic voltage will further introduce a third harmonic current in the interconnected grid. As a result of the current feedback control strategy, the third harmonic current would lead to an additional third harmonic voltage on the ac side of the VSC via the control section. As a result, the third harmonic voltage of the VSC shall be computed as shown in
Figure 2. As can be seen, there are two parts to be included. The first is the third harmonic voltage caused by the fundamental frequency voltage and current, and the second is the distribution of this third harmonic voltage. To compute the final voltage, it is necessary to include the third harmonic on the ac side of the VSC, and, as a result, the steady-state current and modulation signal can be expressed as (8) and (9).
Figure 2.
Illustration of modeling of the third harmonic.
Figure 2.
Illustration of modeling of the third harmonic.
The equations in (8) and (9) can be used to obtain the current on the DC side, which is given as (10):
Accordingly, the DC-link voltage can be computed as
The above dc-link voltage is converted to the ac-side voltage by the switching function, and the phase-A voltage is given as
and the phase-B and phase-C voltages can be computed similarly. As one can see, there will be numerous types of components in the ac voltage. Since we are only interested in the third harmonic components, the third harmonic components in (12) are extracted as in (13).
The expression in (13) shows that the AC-side voltage of the converter contains third harmonic components, which means that the converter acts as a source at the third harmonic. Specifically, we can see that the terms of the third harmonic voltage in (13) can be categorized into two types. The first type is the terms that are determined by fundamental frequency components (m±1 and I±1). The second type is the terms that are affected by external third harmonic components (m±3 and I±3). The first type of third harmonic voltage is an independent voltage source and the second type of third harmonic voltage exhibits an impedance nature. Therefore, (13) indicates that the VSC can be represented as a Thevenin circuit at the third harmonic. In the following, the detailed Thevenin circuit will be presented.
In (13), it is noted that the third harmonic components in the modulation signal can be related to the third harmonic current on the ac side of the VSC by (14):
It is noted that the outer control loop of the VSC generally has a very limited bandwidth (e.g., 100 Hz), so the reference current that is generated by the outer control loop can be considered as harmonic-free (i.e.,
Iref = 0). Substituting (14) into (13) leads to an equation that contains the third harmonic voltage and current, which gives the following model
where
The above final model gives the following key findings:
A Thevenin circuit can represent the VSC at both the positive-sequence and negative-sequence third harmonics, and the analytical expression of the Thevenin circuit can be established.
There is a coupling effect between the positive-sequence third harmonic and negative-sequence third harmonic. The coupling effect mainly originates from the interaction between positive-sequence fundamental frequency components and negative-sequence third harmonic components and vice versa.
The source of the Thevenin circuit is determined only by the positive-sequence and negative-sequence fundamental frequency components, whereas the impedances are determined by both the inner PI regulator and the fundamental frequency voltage and current.
According to the above results, the fundamental frequency components are needed to complete the developed model. For this purpose, the load flow methods regarding the unbalanced power systems with VSC devices are needed. At present, numerous works have been conducted on this aspect [
20,
21,
22,
23]. These methods can be used to compute the positive-sequence and negative-sequence fundamental frequency voltage and current. The procedure to obtain the final model of the VSC can be summarized as the flow chart in
Figure 3.
Compared with other existing models, the proposed model for third harmonic computation has the following advantages:
The mechanism of the production of the non-characteristic third harmonics are fully revealed so that the contributions of the different components to the non-characteristic third harmonics emission can be easily analyzed. Accordingly, the control scheme that can suppress the third harmonic emission can be designed.
The computation of the proposed model is very easy as the analytical expression of the proposed model is given. Compared with other methods, it is more straightforward for users.
The proposed model can be integrated into the existing harmonic power flow tools. As a result, the non-characteristic third harmonic can be computed along with other harmonic components, which can be easily adopted by the users.
4. Verification of the Proposed Model by Time-Domain Simulations
In this section, the developed model of VSC under unbalanced conditions is verified via time-domain simulations in the VSC system of
Figure 4 in MATLAB/Simulink. The VSC has the same control strategy as that in
Figure 1. The system parameters are listed in
Table 1 [
24,
25]. The unbalanced condition is created by adding a negative-sequence fundamental frequency voltage in the AC source (
Vs in
Figure 4). Therefore, the ac voltage source can be expressed as
For illustration, the unbalanced level of the system is defined as
The basic strategy to verify the proposed model is summarized as follows:
Implement the simulation system in MATLAB/Simulink, then set different unbalanced levels to run the simulations. The simulation results give the third harmonic distortion levels.
Based on the parameters of the simulation system and the unbalanced levels, the Thevenin circuit can be computed.
Once the Thevenin circuit is obtained, the third harmonic distortion levels in the system can be computed for different unbalanced levels.
The calculated third harmonic distortion levels are compared with the results obtained from the time-domain simulations in step (1). If these two types of results match well, then it verifies the correctness of the proposed model at third harmonics.
The developed third harmonic model is tested under an unbalanced level of 5%, 15%, 25%, and 35%. The simulated waveforms of the AC-side voltage and current of the VSC are presented in
Figure 5. The waveforms of the VSC’s voltage and current can be used to compute the third harmonic voltage and current using FFT, and the IHDs of the third harmonics for different unbalanced levels are shown in
Figure 6. It can be seen that a larger unbalanced level leads to larger third harmonics. When the unbalanced level is only 5%, the IHD of the positive-sequence third harmonic current is around 2%, while it increases to 11% when the unbalanced level rises to 35%. Additionally, one can also see that the positive-sequence third harmonic is much higher than the negative-sequence third harmonic. The reason is that, according to Equations (16)~(21), the positive-sequence third harmonic is more relevant with the positive-sequence fundamental frequency components, whereas the negative-sequence third harmonic is more relevant with the negative-sequence fundamental frequency components. As a result of the large positive-sequence third harmonic current, the coupling between the negative-sequence third harmonic and positive-sequence third harmonic current cannot be ignored.
It has been mentioned previously that, to obtain the developed model, the positive-sequence and negative-sequence fundamental frequency voltage and current are needed. In the following, the calculated positive-sequence and negative-sequence voltage and current are compared with the simulation results first, and then these values are used to compute the developed model.
4.1. Computation of the Fundamental Frequency Components
The positive-sequence and negative-sequence fundamental frequency voltage and current can be obtained through load flow studies. However, it is noted from the derived model that both the positive-sequence and negative-sequence fundamental frequency components in the modulation signal are needed to compute the developed model. Due to the control and switching process features, the positive-sequence and negative-sequence fundamental-frequency components in the modulation signal can be computed using Equations (24) and (25), respectively.
For the studied four cases, the calculated ac-side positive-sequence fundamental frequency voltage and negative-sequence fundamental frequency current and modulation signals regarding various unbalanced levels are compared in
Table 2 and
Table 3. It can be seen that the calculated fundamental frequency components under various unbalanced conditions match well with the simulated results. Based on the obtained positive-sequence voltage and the negative-sequence current, the corresponding modulation signal can be computed using (24) and (25).
4.2. Verification of the Proposed Model
According to the obtained fundamental frequency components above, the parameters of the developed model in (15) can be computed, and the parameters of this model for different unbalanced levels are listed in
Table 4. It can be seen that, since the positive-sequence third harmonic current is much larger than the negative-sequence third harmonic current, the coupling between the negative-sequence third harmonic voltage and the positive-sequence third harmonic current is more important. According to the calculated parameters, the third harmonic current can be computed using the following equation:
and the third harmonic voltage on the ac side of the VSC can be computed as
The comparison results between the calculated harmonic voltage and current on the ac side of the VSC for the four studied cases are listed in
Table 5 and
Table 6. One can see that both the calculated positive-sequence and negative-sequence third harmonic voltages match well with the simulated results, which verifies the correctness of the proposed model. In the calculations, the coupling between the positive-sequence and negative sequence third harmonics is included. This is very important for the computation for the negative-sequence third harmonics due to the large positive-sequence third harmonic current.
4.3. Impact of the Dc-Link Capacitor
According to the derivations in
Section 2 and
Section 3, the value of the DC-side capacitor plays an important role in the developed model. For the same unbalanced level, the larger the DC-side capacitor, the smaller the third harmonic. This can be verified by the simulation results in
Figure 7, which are based on an unbalanced level of 25%. As one can see, when the DC-side capacitor is increased to 100 μF, the third harmonic voltage and current are significantly decreased. The results can also be supported by the harmonic voltage of the dc side, which is shown in
Figure 8. It is found that the harmonic ripples (mainly second harmonics) in the dc-link voltage for a dc-link capacitor of 100 μF are much smaller than those for a dc-link capacitor of 50 μF. In summary, increasing the DC-link capacitor’s value can resist the generation of the uncharacteristic third harmonic.
5. Discussion
The above results have shown the correctness of the developed model, but it has to be pointed out that the new model can only be developed with the accurate information of the positive-sequence and negative-sequence fundamental frequency voltage and current. If these parameters cannot be provided, then the model cannot be accurately developed. It is noteworthy that the positive-sequence third harmonic distortion is much larger than the negative-sequence third harmonic. Therefore, for the evaluation of the harmonic impact, the level of the positive-sequence third harmonic is sufficient. Moreover, since the negative-sequence harmonic is quite small, its coupling effect to the positive-sequence voltage can be ignored as well. As a result, for the evaluation of the non-characteristic third harmonics, the model can be practically simplified as a Thevenin circuit, and the computation can be simplified accordingly.
The proposed model is able to show the effect of different components on the Thevenin circuit, and, thus, an optimal parameter design can be achieved to reduce the third harmonic distortion. Additionally, the proposed model can help to implement the control schemes for suppressing the third harmonic distortion caused by the VSC under unbalanced conditions.