Next Article in Journal
Simulation Research on the Heating Performance of the Combined System of Solar Energy and Heat-Source Tower Heat Pump in a Hot Summer and Cold Winter Area
Next Article in Special Issue
Harmonics Compensation by Using a Multi-Modular H-Bridge-Based Multilevel Converter
Previous Article in Journal
Grid Parity Analysis of China’s Centralized Photovoltaic Generation under Multiple Uncertainties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CPC-Based Minimizing of Balancing Compensators in Four-Wire Nonsinusoidal Asymmetrical Systems

by
Zbigniew Sołjan
*,
Grzegorz Hołdyński
and
Maciej Zajkowski
Faculty of Electrical Engineering, Bialystok University of Technology, ul. Wiejska 45D, 15-351 Białystok, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(7), 1815; https://doi.org/10.3390/en14071815
Submission received: 18 February 2021 / Revised: 15 March 2021 / Accepted: 22 March 2021 / Published: 24 March 2021
(This article belongs to the Special Issue Control Strategies Applied to Active Power Filters)

Abstract

:
The article presents the possibility of using the currents’ physical components (CPC) theory to generate the reference current of the active power filter (APF). The solution proposed by the authors is based on the cooperation of minimizing balancing compensators (MBC), which, due to their use in 4-wire systems, have been divided into two structures. The first compensator, which purpose is to minimize and balance the reactive current and the unbalanced current of the zero sequence, is built in the star system (STAR-MBC). The purpose of the second compensator, which operation occurs in the delta system (DELTA-MBC), is to minimize and balance the other two components, i.e., the unbalanced current of the negative sequence and the unbalanced current of the positive sequence. The two structures cooperating with each other significantly reduce the currents associated with the reactive elements, i.e., reactive current, and the unbalanced current. As mentioned, these currents are reduced but not compensated to zero or to the reference value. In order for the compensation and balancing to bring the preferable effect, an APF system should be included, which will cooperate with MBC compensators. This solution is presented in this publication. The control of the active part of the hybrid active power filter (HAPF), which was presented in the paper, consists of the reflection of the waveform of the nonsinusoidal active current. In this approach, no current shift in relation to voltage is obtained, but the waveforms of these quantities remain distorted. The reactive current is compensated and the unbalanced currents are balanced. The second definition of generating a reference current can also be used. Through this approach, the active current with a sinusoidal waveform is achieved. The second approach allows for an additional reduction of the three-phase RMS value of the load’s current. In both of these approaches, the active currents flowing through the lines will reflect the amplitude and phase asymmetry that is present in the supply voltage. The APF system will follow the changes in power or load conditions and generate the correct value for the reference current. The calculations presented in the article, as well as the current and voltage waveforms, were created as a result of the constructed mathematical models, which were used for theoretical illustrations. Calculations and waveforms were generated based on a script written in Matlab.

Graphical Abstract

1. Introduction

Low voltage electrical installations in industrial plants are four- or five-wire systems. Due to the very frequent unbalance load in the supply lines, there exists an imbalance of the current. In addition to the unbalanced current, there is also a reactive current associated with nonlinear elements and reactive elements [1,2].
If to the assumption of the existence of the current asymmetry, we also add the existence of the supply source imbalance, which is a real phenomenon in power networks, we obtain a three-phase four-wire system with an asymmetrical nonsinusoidal voltage source and an unbalanced (most often also nonlinear) load [3,4,5,6].
Another aspect that should be discussed is the possibility of application of the compensation devices in systems with currents of several or tens kA. In such systems, hybrid devices should be used [7,8]. They have reactive elements whose function is to initially compensate for high values of line currents. Hybrid systems also include switching devices widely known as active power filters (APF). Their main function is to compensate and balance the remaining currents of significantly lower values [8].
There are many solutions relating to the APF strategy control in the literature [9,10,11,12,13,14,15,16]. The most common approach is to use a mathematical description known as the instantaneous reactive power (IRP) p-q theory [12,13,17,18] or its expand [19,20]. This approach is based on the time domain and therefore is by definition convenient for generating reference currents because the APF strategy control is then significantly faster than methods based on the frequency approach. However, there are also significant disadvantages to this approach. The first is the generation of additional distortion in the form of the 3rd harmonic appearing in the current [1,21,22] and the lack of recapture in physical phenomena, especially in not including the asymmetry of the supply voltage. The lack of such a reflection contributes to the incorrect control of switching elements in the APF. Besides, the most popular time-domain method, which is the IRP p-q, the control is also done by a synchronous reference frame (SRF) [23,24,25] or an approach named simplified SRF (sSRF) [25]. This approach is also known as the d-q method or simplified d-q method. The most popular frequency domain approaches are discrete Fourier transform (DFT) [26,27] or fast Fourier transform (FFT) [28,29]. It is also possible to combine the d-q method and the Fourier transform (dqFourier) [30,31].
For several years now, an increasingly popular method of control is the one based on model predictive control (MPC) [32,33,34]. To improve the quality of energy, MPC controllers cooperate with APF systems, which are controlled by, for example, conservative power theory (CPT) [33]. At this point, we should also mention artificial neural networks (ANN) [35,36] or the wavelet transform [37].
The article presents a description of energy properties in 3-phase 4-wire systems. The considered voltage supply has amplitude and phase asymmetry. Additionally, it is also a nonsinusoidal voltage [3,38,39]. Three parameters related only to one physical quantity were not previously considered in the mathematical description of these types of systems. As it can be seen, it affects the voltage unbalance negatively affecting the load, causing an increase in the unbalanced current as a result of an additional component depending on the harmonic sequence. Determining the equivalent parameters of the load allows for the correct determination of the current components and their physical interpretation. The possibility of physical interpretation and the knowledge of the equivalent parameters allows determining the parameters of an ideal balancing compensator (a description of such a compensator for an asymmetric sinusoidal supply can be found in [40,41,42]). At the next stage, with the knowledge of the ideal parameters of the balancing compensators, it is possible to determine the parameters of the minimizing balancing compensator [1,6,43,44,45,46]. At the moment of writing this article, no one, except the authors, has determined the parameters of an ideal balancing compensator dedicated to 4-wire systems powered by asymmetrical nonsinusoidal voltage. If there is no description of the ideal compensator, the mathematical description and determination of MBC parameters are also impossible. This publication presents the methodology used in the design of minimizing balancing compensators. As it is known, the voltage source and the load are time-varying. Therefore, the article proposes an approach allowing the determination of the reference currents of the active part of the hybrid active power filter. When properly implemented, both approaches will allow the injection of current waveforms of the appropriate shape and value to the line and to the neutral wire. The first approach allows to obtain an active current identical to the supply voltage waveform, so this waveform will be a nonsinusoidal waveform. The second approach makes it possible to generate a reference current of an APF so that only the active current waveform of the fundamental harmonic is visible from the side of the voltage source. Presentation of the possibility of using the CPC Theory in designing both a minimizing balancing compensator and in generating the reference current of an active power filter in one publication is rare.
The energy description of three-phase four-wire systems supplied from a nonsinusoidal asymmetrical voltage source has been presented in Section 2 and Section 3 shows the theoretical illustration of a three-phase four-wire system presenting the correctness of the mathematical description proposed in Section 2. The Section 4 includes considerations and equations presenting the mathematical synthesis of the minimizing balancing compensator. The theoretical illustration of the solution proposed in Section 4 is shown in Section 5 and Section 6 discusses the possibility of generating reference currents injected into the supply lines and the neutral wire. The final Section 7 presents the values and waveforms that can be obtained by connecting an active part to the passive system in the form of an active power filter.

2. Description of Currents’ Physical Components (CPC) Theory in 3-Phase 4-Wire Systems with Asymmetrical Voltage Waveforms

The nonsinusoidal voltage of the source of the distribution system can be described as a three-phase vector, the elements of which are the voltages at the terminals R, S, T, namely:
e n = e R n e S n e T n ,
where: e R n , e S n , e T n are the nonsinusoidal temporal voltages values at the terminals R, S, and T of the voltage source.
This voltage, as asymmetrical and nonsinusoidal, may have symmetrical components for individual harmonics, i.e., the symmetrical component of the positive sequence, the symmetrical component of the negative sequence, and the symmetrical component of the zero sequence. Consequently, the asymmetric nonsinusoidal line voltage is given by the vector:
u n = u R n u S n u T n T ,
where: u R n , u S n , u T n are the nonsinusoidal temporal voltages values at the terminals R, S, and T of the load.
Which may have three components for each harmonic:
u n = u n p + u n n + u n z ,
where: u n p is the symmetrical component of the positive sequence of the supply voltage, u n n is the symmetrical component of the negative sequence of the supply voltage, u n z is the symmetrical component of the zero sequence of the supply voltage.
Figure 1 shows a three-phase four-wire system supplying an unbalanced linear time-invariant load (LTI).
Figure 1 shows the three-phase four-wire system supplying an unbalanced linear time-invariant load (LTI):
u t = u R t u S t u T t = 2 Re n N U R n U S n U T n e j n ω 1 t = 2 Re n N U n e j n ω 1 t ,
where: U R n , U S n , U T n are the nonsinusoidal complex RMS (CRMS) values of the voltages at the terminals R, S, and T of the load; U n is a three-phase vector of CRMS values of harmonics of phase voltages.
The vector of the currents of lines can be represented in the same way, namely:
i t = i R t i S t i T t = 2 Re n N I R n I S n I T n e j n ω 1 t = 2 Re n N I n e j n ω 1 t ,
where: I R n , I S n , I T n are the nonsinusoidal CRMS values of the currents in lines R, S, and T; I n is a three-phase vector of CRMS values of harmonics of line currents.
Considering the purpose of the article, only the final relationships necessary for understanding the next section will be presented in this section. Detailed information is presented in the authors’ publication [3].
Due to the currents’ physical components (CPC) theory, the current of an unbalanced LTI load supplied with an asymmetrical nonsinusoidal voltage [3] may have four components in total. The component of the current related to the permanent energy transmission from the source to the load is the active current i a , the waveform of which takes the form:
i a = G e u = n N i a n = 2 Re n N G e U n e j n ω 1 t = 2 Re n N G e U n p + U n n + U n z e j n ω 1 t = 2 Re n N G e 1 p U n p + 1 n U n n + 1 z U n z e j n ω 1 t ,
where: U n p , U n n , U n z are the CRMS values of symmetrical voltage components accomplished from the Fortescue transformation; 1 p , 1 n , 1 z are the unit vectors of symmetrical components described by the Fortescue transformation; U n p , U n n , U n z are the vectors of the CRMS values of symmetrical voltage components accomplished from the Fortescue transformation.
The three-phase RMS value of an active current is associated to the equivalent conductance of the entire system G e and to the three-phase RMS value of the supply voltage u :
i a = G e u = P u ,
The component of the current associated with the presence of the difference between the equivalent conductance of the entire system G e and the equivalent conductance of individual harmonics G b n is the scattered current i s , the waveform of which can be presented as follows:
i s = n N i s n = 2 Re n N G b n G e U n e j n ω 1 t = 2 Re n N G b n G e U n p + U n n + U n z e j n ω 1 t = 2 Re n N G b n G e 1 p U n p + 1 n U n n + 1 z U n z e j n ω 1 t
The three-phase RMS value of the scattered current i s is calculated from the relationship:
i s = n N G b n G e 2 u n 2 ,
The component of the current associated with the shift between the voltage waveform and the current waveform for individual harmonics and related to the equivalent parameter of the circuit in the form of equivalent susceptance B b n is the reactive current i r , and its waveform is described as follows:
i r = n N i r n = 2 Re n N j B b n U n e j n ω 1 t = 2 Re n N j B b n U n p + U n n + U n z e j n ω 1 t = 2 Re n N j B b n 1 p U n p + 1 n U n n + 1 z U n z e j n ω 1 t ,
The three-phase RMS value of the reactive current i r is equal:
i r = n N B b n 2 u n 2 ,
where: u n is the three-phase RMS value of the voltage of the respective harmonic.
In the system, as a result of imbalanced of a load, there is also a fourth component of the current, called the unbalance current i u . The waveform of the unbalanced current is resulted by subtracting from the waveform of the total current i three components in the form of active current i a , reactive current i r , and scattered current i s :
  i u = n N i u n = 2 Re n N I u n e j n ω 1 t = 2 Re n N I n I b n e j n ω 1 t   = 2 Re n N I n Y b n U n e j n ω 1 t = 2 Re n N I n G e + G b n G e + j B b n U n e j n ω 1 t           ,   = 2 Re n N I n I b n 1 p U n p + 1 n U n n + 1 z U n z e j n ω 1 t = 2 Re n N I n I b n U n p + U n n + U n z e j n ω 1 t
where I b n = Y b n U n is the vector of line currents consequent from the active and reactive powers of respective harmonics.
The three-phase RMS value of the unbalanced current i u is:
i u = i 2 i a 2 + i r 2 + i s 2 ,
Additionally, the unbalanced current can be divided into three symmetrical components of the appropriate order [3,4,5]. The first component of the unbalanced current i u is the unbalanced current of the positive sequence i u p , which is related to the voltage asymmetry dependent admittance Y d n , the generalized unbalanced admittance of the positive sequence Y n p , the generalized unbalanced admittance of the negative sequence Y n n , and the generalized unbalanced admittance of the zero sequence Y n z :
i u p = n N i u n p = 2 Re n N I u n p e j n ω 1 t = 2 Re n N 1 p Y d n U n p + Y n p U n p + Y n z U n n + Y n n U n z e j n ω 1 t ,
The three-phase RMS value of the unbalanced current of the positive sequence is equal:
i u p = 3 n N Y d n U n p + Y n p U n p + Y n z U n n + Y n n U n z 2 ,
The unbalanced current of the negative sequence i u n can be determined analogically, the waveform of which can be described by the relationship:
i u n = n N i u n n = 2 Re n N I u n n e j n ω 1 t = 2 Re n N 1 n Y d n U n n + Y n n U n p + Y n p U n n + Y n z U n z e j n ω 1 t ,
The three-phase RMS value of the unbalanced current of the negative sequence is equal:
i u n = 3 n N Y d n U n n + Y n n U n p + Y n p U n n + Y n z U n z 2 ,
The last component is the unbalanced current of the zero sequence i u z . The waveform of this current is:
i u z = n N i u n z = 2 Re n N I u n z e j n ω 1 t = 2 Re n N 1 z Y d n U n z + Y n z U n p + Y n n U n n + Y n p U n z e j n ω 1 t ,
The three-phase RMS value of the unbalanced current of the zero sequence is equal:
i u z = 3 n N Y d n U n z + Y n z U n p + Y n n U n n + Y n p U n z 2 ,
In relation to the fact that all components are mutually orthogonal [3], the three-phase RMS value of the load’s current, according to currents’ physical components theory, is:
i 2 = i a 2 + i s 2 + i r 2 + i u p 2 + i u n 2 + i u z 2 ,
The respective powers are obtained by multiplying the three-phase RMS value of the supply voltage and the three-phase RMS values of the individual components of the load current (7), (9), (11), (12), (17), and (19):
S = P 2 + Q 2 + D s 2 + D u p 2 + D u n 2 + D u z 2 ,
The power factor λ can be expressed in terms of the powers’ components or the currents’ components:
λ = P S = P P 2 + Q 2 + D s 2 + D u p 2 + D u n 2 + D u z 2 = i a i = i a i a 2 + i r 2 + i s 2 + i u p 2 + i u n 2 + i u z 2 ,
As can be seen, the power factor λ consists of six powers or currents which, apart from the active power P or the three-phase RMS value of the active current, are lower than its value. It should also be mentioned that in systems with nonsinusoidal waveforms, the power factor λ will not reach the value equal to 1, as long as there is the scattered current in the considered system.

3. Theoretical Illustration—Currents’ Components

For theoretical illustration, the system shown in Figure 2 has been prepared. This is a three-phase four-wire system supplied with an asymmetrical nonsinusoidal voltage. The load is unbalanced and is an LTI load type. In the voltage source, in addition to the fundamental harmonic with a frequency of 50 Hz, harmonics 3rd, 5th, and 7th are also presence.
The values of resistance as well as capacitive and inductive reactance for the fundamental harmonic, shown in Figure 2, are presented in Table 1.
Table 2 presents the values of the voltages harmonics and the values of the harmonics of the line currents of the load.
Figure 3 shows the three-phase voltage waveform at the terminals of the analyzed load, and the three-phase waveform of the line current is shown in Figure 4.
The equivalent parameters of the load presented in Figure 1 (in accordance with CPC Theory [3]) are listed in Table 3.
Extensive information why the generalized unbalanced admittance of the positive sequence is always 0 is presented in [3].
Based on the information included in Table 2, the three-phase RMS values of the supply voltages are 392.43 V, and the three-phase RMS values of the currents of the line are 294.09 A.
Table 4 presents the three-phase RMS values of currents’ components for appropriate harmonics described in CPC Theory (7), (9), (11), (15), (17), and (19). In Table 4 are also included the total three-phase values for specific currents’ components and the total value of the three-phase current (20).
Based on (7), (9), (11), (15), (17), and (19) and using the already known three-phase RMS values of appropriate components, the power factor λ (22) is 0.593.
As can be seen from the value of the power factor λ , it follows that the considered system should be subjected to compensation and balancing. The balancing compensation can be performed passively by using reactance elements in the form of chokes and/or capacitors. Compensation and balancing can also be carried out actively using an active power filter (APF). It is also possible to use a hybrid system in the form of a designed balancing compensator made of passive elements and an active part in the form of an APF. The third solution is a very often used compromise between the accuracy of compensation and symmetrization and the economic justification for improving the quality of electricity.

4. The Compensator Minimizes the Reactive Current and the Unbalanced Current on the Basis of the Currents’ Physical Components Theory in 3-p 4-w Systems

Due to the number of elements in one branch, the most advantageous condition is that such a branch has only one reactance element. Unfortunately, in the real conditions of the power system, if the compensator branch has a capacitive impedance, it can very often resonate with the system for higher harmonics. To prevent such a phenomenon, the compensator branch should be purely inductive or it should be connected in series with the capacitor, forming the LC branch [6,43,44,45,46].
The vector of the complex RMS values of nonsinusoidal voltages in the branches of balancing compensators takes the form:
U n = U RS n U ST n U TR n U R n U S n U T n T ,
where the nonsinusoidal phase-to-phase voltages are described in [46,47], and the nonsinusoidal phase voltages by Equation (4)
The vector of the branch currents of an ideal balancing compensator can be represented as:
j T = 2 Re n N j T n U n e j n ω 1 t ,
where the vector T n is:
T n = T RS n 0 0 0 0 0 0 T ST n 0 0 0 0 0 0 T TR n 0 0 0 0 0 0 T R n 0 0 0 0 0 0 T S n 0 0 0 0 0 0 T T n ,
The vector T n consists of two components relating to ideal balancing compensators. The elements T ST n , T RS n , T TR n refer to the value of the ideal balancing compensator (the index “n” means the order of harmonics) in the delta structure. So they are phase-to-phase susceptances. The elements of the T R n , T S n , T T n refer to the value of the ideal balancing compensator in the star structure. So they are phase-to-neutral susceptances.
For the minimizing balancing compensator (MBC), the reactive current, and the unbalanced currents, the vector of the branch currents can be represented as follows:
j D = 2 Re n N j D n U n e j n ω 1 t ,
where the vector D n is:
D n = D RS n 0 0 0 0 0 0 D ST n 0 0 0 0 0 0 D TR n 0 0 0 0 0 0 D R n 0 0 0 0 0 0 D S n 0 0 0 0 0 0 D T n ,
The vector D n , identical to the vector T n , consists of two compensators of different structures. The main difference between the elements in both vectors is the ideality of the constituent elements. The susceptances D RS n , D ST n , D TR n in the delta system are designed to minimize the value of unbalanced currents of the negative and positive sequences, but not for ideal compensation. The same is the case with susceptances D R n , D S n , D T n . The selection of the parameters of these susceptances is to contribute to the reduction of reactive current and the unbalanced current of the zero sequence. In order to be able to optimally minimize the reactive current and the unbalanced current, the minimum distance should be determined for these two vectors.
The distance of the j T and j D vectors can be defined as the RMS value of the difference of their currents, namely:
d = j T j D = Δ j ,
Two compensators with the susceptances D RS n , D ST n , D TR n and D R n , D S n , D T n , compensate and balance the unbalanced currents of the positive and negative sequences (the balancing compensator in delta structure) and the reactive current and the unbalanced current of the zero sequence (the balancing compensator in star structure), respectively. If the RMS value of the difference of hypothetical currents, ideal balancing compensators, and the minimizing balancing compensators built of at most two-elements LC, equal to:
Δ j = 2 Re n N T n D n U n e j n ω 1 t = 2 Re n N Δ J n e j n ω 1 t ,
is as small as possible, then the square of the RMS value of this current, otherwise the square of the distance of the currents j T and j D is equal to:
d 2 = Δ j 2 = Re n N Δ J n T Δ J n * = n N Δ J RS n 2 + Δ J ST n 2 + Δ J TR n 2 + Δ J R n 2 + Δ J S n 2 + Δ J T n 2 = n N T RS n D RS n 2 U RS n 2 + T ST n D ST n 2 U ST n 2 + T TR n D TR n 2 U TR n 2 + T R n D R n 2 U R n 2 + T S n D S n 2 U S n 2 + T T n D T n 2 U T n 2
The square of the currents distances j T and j D is therefore the sum of three nonnegative, mutually separate terms for each of the balancing compensator structures, the general form of which can be written as follows:
d = j T j D = n N Δ J n T Δ J n * = n N k K T k n D k n 2 U k n 2 ,
where k is a symbol in the set K = RS , ST , TR , R , S , T .
On the grounds of that the order of summing the harmonic orders n and the branch indexes k can be different, the distance d could be expressed as follows:
d = n N k K T k n D k n 2 U k n 2 = k K d k 2 ,
As can be seen, the values under the square root will be nonnegative numbers, so it is not possible for the values to eliminate.
The distance d is minimal when the individual distances d k for individual branches of the balancing compensators are also minimal, namely:
n N T k n D k n 2 U k n 2 = minimum ,
It is known that the RMS value of the fundamental voltage harmonic, whether phase-to-phase or phase, is much higher than the RMS value of any other harmonic [6,43,44,45,46]. Consequently, if in individual branches of one of the two compensators no resonance appears in the nearby of the harmonic frequency of the supply voltage, then the dominant component in the distance d of the j T and j D vectors is the difference in susceptance between the ideal and minimizing compensator in a given branch for the fundamental frequency, i.e.,: TXY1—DXY1 or TX1—DX1. This assumption determines the choice of susceptance for the fundamental frequency DXY1 or DX1, namely, their signs must be the same as those of TXY1 or TX1 susceptance.
It was assumed that if the susceptance of TXY1 or TX1 is negative, then the compensator with a delta structure between the TXY terminals should have a branch with a choke, because its susceptance is also negative, or if the compensator with a star structure between the phase terminal and the neutral wire TX should also have branch with a choke. The inductance of such a choke, regardless of whether it is connected between phases or phase and neutral wire, should minimize the value:
n N T k n + 1 n   ω 1 L k 2 U k n 2 = minimum ,
If the line-to-line susceptance TXY or the line susceptance TX is positive, the compensator should have a capacitor in the branch, but it should be connected in series with the choke, so that for higher order harmonics this branch has an inductive character and positive for the fundamental frequency:
d k 1 = ω 1 C k 1 ω 1 2   L k C k ,
This requires that the resonant frequency of such a branch be higher than the fundamental frequency, i.e.,:
ω r = 1 L k C k ,
The inductance and capacitance of such a branch should minimize the following expression:
n N T k n n ω 1 C k 1 n 2   ω 1 2   L k C k 2 U k n 2 = minimum ,
The inductance L k satisfying Equation (34) can be found from the condition:
d d L k n N T k n + 1 n ω 1   L k 2 U k n 2 = 0 ,
By fulfilling the expression (38), i.e., by calculating the inductance of a given branch, we can treat it as the optimal inductance of that branch, therefore:
L k , opt = n N U k n 2 n 2 ω 1 n N T k n U k n 2 n ,
The left side of the expression (37) is a function of two variables, i.e., inductance L k and capacitance C k . The inductance of L k is a continuously decreasing function, which means that there is no minimum for any finite value. Its choice depends on the person accomplishing the calculations [43,44,45]. After choosing the inductance, the capacitance C k can be calculated from formula (37). When looking for the minimum in expression (37), we have to fulfill the condition:
d d C k n N T k n n ω 1 C k 1 n 2   ω 1 2   L k C k 2 U k n 2 = 0 ,
The distance defined by the expression (40) for any branch has a minimum due to the optimal capacitance of the capacitor when:
n N n T k n U k n 2 1 n 2   ω 1 2   L k , opt C k , opt 2 n N n 2 ω 1 C k , opt U k n 2 1 n 2   ω 1 2   L k , opt C k , opt 3 = 0 ,
Considering the optimal capacity C k , opt occurring in the dependence (41), this equation is an implicit equation. Therefore, calculating this capacity requires numerical methods, especially in an iterative process. The iterative formula for determining the optimal capacity takes the form:
C k , i + 1 = n N n T k n U k n 2 1 n 2 ω 1 2 L k , opt C k , i 2 ω 1 n N n 2 U k n 2 1 n 2 ω 1 2 L k , opt C k , i 3 ,
After the ith iteration, it will concur to obtain several capacities, which in the final stage of calculation should cause to the optimal capacity.

5. Theoretical Illustration—Symmetrization and Compensation of Currents’ Components and Minimalization of the Structure of Balancing Compensators

For the theoretical illustration of the minimizing balancing compensation, the system analyzed in Figure 2 was chosen. According to the information in Chapter 4, the minimization is primarily focused on incomplete balancing compensation of reactive current and the unbalanced current, as well as limiting the number of reactive elements to a maximum of two in one branch [6,43,44,45,46].
Based on Equation (36), the resonant frequency was chosen with the value ω r = 1.07 ω 1 . Table 5 presents the values of the reactance elements of the minimizing balancing compensators (MBC) with star (STAR-MBC) and delta (DELTA-MBC) structures, calculated on the basis of expressions (39) and (42). The required optimal capacity in each branch was estimated in the fifth iteration.
Figure 5 shows the analyzed system which connected the minimizing balancing compensators in both structures.
Table 6 presents the three-phase RMS values of the component currents, defined in the CPC theory, after the minimizing the balancing compensation shown in Figure 5.
As can be observed in Table 6, as a result of connecting the compensators MBC, the three-phase RMS currents for higher harmonics have increased. It is related to the purposeful inaccuracy of the inductive susceptances so that there is no resonance phenomenon between the MBC and the inductance of the power network. The increase in the three-phase RMS values of the component currents does not, however, result in the flow of a large amount of energy, because the amplitudes of the harmonic voltages are also relatively low in relation to the fundamental harmonic voltage.
Figure 6 shows the waveform of line currents after the minimizing the balancing compensation.
Due to the use of the minimizing balancing compensation, the power factor λ (22) increased to the value of 0.986.
This means that in systems with the scattered current associated with a change in conductance on the order of harmonics, it is not possible to obtain the power factor λ equal to unity. For this to become, such a system should use power electronic elements, known as an active power filter.

6. CPC Theory—The Possibility of Using the Active Power Filter to Generate the Reference Current

The active power filter (APF) as a compensator for reactive power (reactive current) and unbalanced power (unbalanced current) has to be equipped with a measurement and identification system for the load’s parameters. Currents and voltages can be measured at the spot, where it will be operated without feedback, or where the feedback loop will be required. In the first case, the system will be investigated as a stable system, while in the second case, problems with stability may appear. The problem of stability must therefore be included in the APF control algorithm.
Stability is not investigated in this article. Essentially, the identification of the parameters of an unbalanced three-phase load, which is performed by measuring the currents and voltages at its terminals, is significant for compensation. The identification of the load has to be made on the basis of the knowledge of the current and voltage waveforms, and more precisely, the sequences of their discrete values, provided by the waveform sampling system and the analog-to-digital converter. When the sampled sequences are known, the complex RMS values of the supply voltage, load currents, and the current in the neutral conductor can be calculated, i.e.,: U R , U S , U T , I R , I S , I T and I N , which will be the resultant of the line currents. The most common method for this is the discrete Fourier transform (DFT) together with the fast Fourier transform (FFT) algorithm [9,16,48]. It is also possible to abbreviate the computational time by using the algorithm described in [49].
Before launching the capability of generating a reference current in an active power filter, the appropriate APF system should be chosen first. Among many available solutions [9,10,12,13,14,15,16,17,18,50,51,52], four-wire systems [10,15] should be used. The use of four-wire systems has the advantage over three-wire systems that in addition to injecting reference currents to compensate for the line currents iR (t), iS (t), and iT (t), it is also possible to inject a reference current to compensate for the current in the neutral wire iN (t). As it is known, the current in the neutral wire with non-linear loads or with unbalanced linear loads can reach very high values.
To obtain the reference current i RST r e f t of an active power filter or a hybrid power filter (the active part of this filter), we use the system of equations:
i RST r e f t = i R r e f t = 2 Re n = 1 N i rR n t + n = 1 N i sR n t + n = 1 N i uR n p t + n = 1 N i uR n n t + n = 1 N i uR n z t i S r e f t = 2 Re n = 1 N i rS n t + n = 1 N i sS n t + n = 1 N i uS n p t + n = 1 N i uS n n t + n = 1 N i uS n z t   i T r e f t = 2 Re n = 1 N i rT n t + n = 1 N i sT n t + n = 1 N i uT n p t + n = 1 N i uT n n t + n = 1 N i uT n z t  
In this way, despite the significant number of calculations, it is possible to generate a reference current which, when injected into the line, guarantees the active current waveform seen from the power source. As can be seen, the active current waveform will be identical to the supply voltage waveform, i.e., it will be a nonsinusoidal waveform with the same asymmetry and nonsinusoidal character as the voltage waveform. Additionally, the system of Equation (43) should be supplemented with the equation of the reference current injected into the neutral wire, namely:
i N r e f t = 2 Re n = 1 N i aZ n t + n = 1 N i rZ n t + n = 1 N i sZ n t + n = 1 N i u n z t ,
From Equation (44), it can be seen that not only a nonsinusoidal unbalanced current of the zero sequence flows in the neutral wire, but also the zero sequences of nonsinusoidal active, reactive, and scattered currents. It is caused only by the asymmetry of the supply voltage. At the moment when the supply voltage will be symmetrical, the components related to the zero sequence current of the respective components will disappear from Equation (44).
Figure 7 shows an example of APF connection to the LTI load.
As has already been mentioned in Chapter 5, due to the use of the minimizing balancing compensation, the power factor (22) was equal to 0.986. Using the system of Equations (43) and (44), we can obtain the power factor equal to 1.
It should be noted, however, that the definition of this power factor consists of the sum of the active powers P and the sum of apparent powers S. If we define the power factor as the value of the active power of the fundamental harmonic P 1 or the value of the active current of the fundamental harmonic i a 1 , then the power factor λ 1 :
λ 1 = P 1 P = i a 1 i a ,
It will still be lower than unity and it will be 0.979.
To obtain the active current of the fundamental harmonic, the relation (43) should be modified to the form:
i RST r e f t = i R r e f t = 2 Re n = 2 N i aR n t + n = 1 N i rR n t + n = 1 N i sR n t + n = 1 N i uR n p t + n = 1 N i uR n n t + n = 1 N i uR n z t i S r e f t = 2 Re n = 2 N i aS n t + n = 1 N i rS n t + n = 1 N i sS n t + n = 1 N i uS n p t + n = 1 N i uS n n t + n = 1 N i uS n z t     i T r e f t = 2 Re n = 2 N i aT n t + n = 1 N i rT n t + n = 1 N i sT n t + n = 1 N i uT n p t + n = 1 N i uT n n t + n = 1 N i uT n z t  
Compensation of the current in the neutral wire will still be achieved with the use of the reference current described by Equation (44), because the zero sequence active current is associated with each harmonic, including the fundamental harmonic.
By using Equations (44) and (46), it is possible to receive the power factor λ equal to unity and associated with the fundamental harmonic of the active current. The waveform of the active current of the fundamental harmonic will have the same asymmetry as the fundamental harmonic of the voltage waveform.

7. Theoretical and Simulation Illustration—The Analysis of the System with an Active Power Filter Connected

The summary of the theoretical analysis is the determination of the reference currents in the X-X’ and Y-Y’ cross-sections and the depiction of the waveforms of the line current that can be obtained on the basis of the relationship (43) and (46) in the Z-Z’ cross-section.
Figure 8 does not present the load considered earlier, only the minimizing balancing compensators (DELTA-MBC and STAR-MBC) and the connected energy power filter (APF) with marked cross-sections that will be considered.
Basing the APF control algorithm based on the relationship (43), the waveform of the three-phase reference current is received, shown in Figure 9.
Based on the reference current generated in this way (Figure 9) and the relationship (43) in the cross-section Z-Z’, a nonsinusoidal active current waveform will be obtained as shown in Figure 10.
Due to the use of the passive compensation in the form of the minimizing balancing compensation, it was capable to reduce the RMS value of the currents in the respective lines (Table 7). The application of the active filtration and the injection of the reference currents to the supply lines will contribute to a further reduction in the RMS value of the currents of the line.
The values presented in Table 7 have been calculated based on an algorithm written in the form of a script in the Matlab software. The algorithm created in Matlab was based on the relationships described in Section 2, Section 4, and Section 6.
Based on Equation (44), the reference current injected into the neutral wire (cross-section Y-Y’) has the waveform shown in Figure 11.
Due to the use of passive compensation in the form of the minimizing balancing compensation, it was possible to reduce the RMS value of the current in the neutral wire (Table 7). The use of active filtration and the injection of the reference current (44) to the neutral wire will contribute to the complete (in a theoretical example) compensation of the current (Table 7).
Based on the APF control algorithm ensuing from Equation (46), the waveform of the three-phase reference current is obtained, as shown in Figure 12.
Based on the reference current generated in this way (Figure 12) and the relation (46) in the Z-Z’ cross-section, a sinusoidal active current waveform is obtained (Figure 13).
The application of the active filtration expressed in relation (46) and the injection of the reference currents to the supply lines will contribute to a further reduction of the RMS values of the currents of the lines (Table 7).
As can be seen from Figure 13, the waveform of the active currents in the respective lines is charged with the difference in amplitude and phase shift, as it responds to the values of amplitude and phase shift in the supply voltage of the fundamental harmonic. When using CPC theory (approach for the description of systems with asymmetrical voltage waveforms) as a method for generating reference currents, it should be remembered that without symmetrization of the voltage supply, the active currents (version with sinusoidal or nonsinusoidal waveforms), will reflect the amplitude and phase asymmetry present in this supply voltage.

8. Conclusions

The article presents an unbalanced linear time-invariant load supplied with asymmetrical nonsinusoidal voltage in four-wire systems. It can be almost perfectly balanced and its reactive power can be compensated also almost perfectly. The use of the minimizing balancing compensation (the passive part of the hybrid filter) increases the power factor to a value close to unity. The application of an active power filter results in additional benefits of the compensation of the reactive, scattered, and unbalanced currents that cannot be compensated using passive elements in the form of the minimizing balancing compensator.
The solution proposed on the basis of Equation (43) in generating reference currents results in obtaining only the nonsinusoidal active current. In such an approach, the reactive power and the unbalanced power do not exist, because the active current waveform is identical to the nonsinusoidal supply voltage waveform (there is no shift between voltage and current waveform). The second approach, which is described with the relation (46), focuses only on obtaining a sinusoidal active current reflecting the voltage waveform of the fundamental harmonic. Additionally, the second approach improves the stability and energy efficiency of the supply network.
The control methods of APF proposed in the article clearly show the possibility of including the asymmetry of the supply voltage, which is so often omitted by other authors. The phenomenon of voltage asymmetry is present to a greater or lower rank in every supply system. As it has been shown in the article, accepting this phenomenon allows the correct mathematical description of LTI loads and, at the same time, contributes to a more precise control of generation of the reference current of the APF. This is because previously no one considers the three coefficients of asymmetry, which together can correctly describe the admittance dependent on voltage asymmetry.
The authors of the article, for some time, have been focusing on the correct description of energy properties, which are closely related to physical phenomena. Due to this approach, a correct mathematical description is obtained, which is a further stage to use for theoretical and then real improvement of the power supply conditions in single- or three-phase systems.

Author Contributions

The paper was written by Z.S.; conceptualization, Z.S., G.H. and M.Z.; methodology, Z.S.; software, Z.S.; validation, Z.S., G.H. and M.Z.; formal analysis, Z.S.; investigation, Z.S.; resources, Z.S.; data curation, Z.S.; writing—original draft preparation, Z.S.; writing—review and editing, Z.S. and M.Z.; visualization, Z.S.; supervision, G.H. and M.Z.; project administration, Z.S.; funding acquisition, Z.S., M.Z. and G.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Science and Higher Education in Warsaw, Poland at the Białystok University of Technology in Białystok, Poland under research subsidies No. WZ/WE-IA/2/2020, and No. WZ/WE-IA/3/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Czarnecki, L.S. Powers in Electric Systems with Nonsinusoidal Current and Voltage Waveforms; Oficyna Wydawnicza Politechniki Warszawskiej: Warszawa, Poland, 2005. (In Polish) [Google Scholar]
  2. Hanzelka, Z. Jakość Dostawy Energii Elektrycznej: Zaburzenia Wartości Skutecznej Napięcia; Oficyna Wydawnicza AGH: Kraków, Poland, 2013. [Google Scholar]
  3. Sołjan, Z.; Hołdyński, G.; Zajkowski, M. CPC-based current′s description in three-phase four-wire systems with asymmetrical and nonsinusoidal waveforms. Bull. Pol. Acad. Sci. Tech. Sci. 2020, 68, 1127–1134. [Google Scholar] [CrossRef]
  4. Czarnecki, L.S.; Haley, P.M. Currents′ Physical Components (CPC) in Four-Wire Systems with Nonsinusoidal Symmetrical Voltage. Przeg. Elektr. 2015, 91, 48–53. [Google Scholar]
  5. Czarnecki, L.S.; Haley, P.M. Power Properties of Four-Wire Systems at Nonsinusoidal Supply Voltage. IEEE Trans. Power Deliv. 2015, 31, 513–521. [Google Scholar] [CrossRef]
  6. Czarnecki, L.S. CPC-based reactive balancing of linear loads in four-wire supply systems with nonsinusoidal voltage. Przeg. Elektr. 2019, 95, 1–8. [Google Scholar]
  7. Czarnecki, L.S. Currents′ Physical Components (CPC)—based Power Theory. W review, Part I: Power properties of electrical Circuits and Systems. Przeg. Elektr. 2019, 95, 1–11. [Google Scholar]
  8. Czarnecki, L.S. Currents′ Physical Components (CPC)—based Power Theory. The review, Part II: Filters and reactive, switching and hybrid compensators. Przeg. Elektr. 2020, 96, 1–11. [Google Scholar]
  9. Karbasforooshan, M.S.; Monfared, M. An adaptive recursive discrete fourier transform technique for the reference current generation of single-phase shunt active power filters. In Proceedings of the 7th Power Electronics, Drive Systems & Technologies Conference (PEDSTC), Tehran, Iran, 16–18 February 2016; pp. 253–259. [Google Scholar]
  10. Quinn, C.; Mohan, N. Active filtering of harmonic currents in three-phase, four-wire systems with three-phase and single-phase nonlinear loads. In Proceedings of the IEEE 7th Applied Power Electronics Conference, Boston, MA, USA, 23–27 February 1992; pp. 829–836. [Google Scholar]
  11. Czarnecki, L.S.; Almousa, M. Adaptive Balancing of Three-Phase Loads at Four-Wire Supply with Reactive Compensators and Nonsinusoidal Voltage. In Proceedings of the 2020 IEEE Texas Power and Energy Conference (TPEC), College Station, TX, USA, 6–7 February 2020; pp. 1–6. [Google Scholar] [CrossRef]
  12. Patrascu, A.; Popescu, M. Comparative active current calculation by p-q and CPC Theories. Ann. Univ. Craiova 2011, 35, 25–30. [Google Scholar]
  13. Bitoleanu, A.; Popescu, M.; Dobriceanu, M.; Nastasoiu, F. Current decomposition methods on p-q and CPC Theories for active filtering reasons. WSEAS Trans. Circ. Syst. 2008, 7, 869–878. [Google Scholar]
  14. Akagi, H.; Isozaki, K. A Hybrid Active Filter for a Three-Phase 12-Pulse Diode Rectifier Used as the Front End of a Medi-um-Voltage Motor Drive. IEEE Trans. Power Electron. 2012, 27, 69–77. [Google Scholar]
  15. Wang, Y.; Wong, M.; Lam, C. Historical review of parallel hybrid active power filter for power quality improvement. In Proceedings of the IEEE Region 10 Humanitarian Technology Conference, Cebu, Philippines, 9–12 December 2015; pp. 1–6. [Google Scholar]
  16. Nicolae, P.M.; Popa, D.L. Improvement of the controller algorithms for active power filters operating under distorted and unbalanced grid voltages. In Proceedings of the 14th International Conference On Industrial Informatics (INDIN), Futuroscope-Poitiers, France, 8–21 July 2016; pp. 276–281. [Google Scholar]
  17. Abdullah, A.; Biswal, G.R.; Roy, A.K. Modeling and control of hybrid power filter using p-q theory. In Proceedings of the IEEE 1st International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES), Delhi, India, 4–6 July 2016; pp. 1–6. [Google Scholar]
  18. Akagi, H. New trends in active filters for power conditioning. IEEE Trans. Ind. Appl. 1996, 32, 1312–1322. [Google Scholar] [CrossRef] [Green Version]
  19. Kim, H.; Blaabjerg, F.; Bak-Jensen, B.; Choi, J. Instantaneous power compensation in three-phase systems by using p-q-r theory. IEEE Trans. Power Electron. 2002, 17, 701–710. [Google Scholar]
  20. Kim, K.; Blaabjerg, F.; Bak-Jensen, B. Spectral analysis of instantaneous powers in single-phase and three-phase systems with use of p-q-r theory. IEEE Trans. Power Electron. 2002, 17, 711–720. [Google Scholar]
  21. Czarnecki, L.S. Effect of supply voltage asymmetry on IRP p-q-based switching compensator control. IET Power Electron. 2010, 3, 11–17. [Google Scholar]
  22. Czarnecki, L.S. Effect of supply voltage harmonics on IRP-based switching compensator control. IEEE Trans. Power Electron. 2009, 24, 483–488. [Google Scholar]
  23. Monfared, M.; Golestan, S.; Guerrero, J.M. A New Synchronous Reference Frame-Based Method for Single-Phase Shunt Active Power Filters. J. Power Electron. 2013, 13, 692–700. [Google Scholar] [CrossRef] [Green Version]
  24. Jain, N.; Gupta, A. Comparison between two compensation current control methods of shunt active power filter. Int. J. Eng. Res. Gen. Sci. 2014, 2, 603–615. [Google Scholar]
  25. Hoon, Y.; Radzi, M.A.M.; Hassan, M.K.; Mailah, N.F.; Wahab, N.I. A Simplified Synchronous Reference Frame for Indirect Current Controlled Three-level Inverter-based Shunt Active Power Filters. J. Power Electron. 2016, 16, 1964–1980. [Google Scholar] [CrossRef] [Green Version]
  26. Mattavelli, P.; Marafao, F. Repetitive-Based Control for Selective Harmonic Compensation in Active Power Filters. IEEE Trans. Ind. Electron. 2004, 51, 1018–1024. [Google Scholar] [CrossRef]
  27. Maza-Ortega, J.M.; Rosendo-Macias, J.A.; Gomez-Exposito, A.; Ceballos-Mannozzi, S.; Barragan-Villarejo, M. Reference Current Computation for Active Power Filters by Running DFT Techniques. IEEE Trans. Power Deliv. 2010, 25, 1986–1995. [Google Scholar] [CrossRef]
  28. Vodyakho, O.; Mi, C.C. Three-Level Inverter-Based Shunt Active Power Filter in Three-Phase Three-Wire and Four-Wire Systems. IEEE Trans. Power Electron. 2009, 24, 1350–1363. [Google Scholar] [CrossRef]
  29. Chen, K.F.; Mei, S.L. Composite Interpolated Fast Fourier Transform with the Hanning Window. IEEE Trans. Instrum. Meas. 2010, 59, 1571–1579. [Google Scholar] [CrossRef]
  30. Areerak, K.-L.; Areerak, K.-N. The comparison study of harmonic detection methods for shunt active power filters. World Acad. Sci. Eng. Technol. 2010, 4, 204–209. [Google Scholar]
  31. Sujitjorn, S.; Kulworawanichpong, T.; Areerak, K.-L. The DQ Axis With Fourier (DQF) Method for Harmonic Identification. IEEE Trans. Power Deliv. 2006, 22, 737–739. [Google Scholar] [CrossRef]
  32. Tarisciotti, L.; Formentini, A.; Gaeta, A.; Degano, M.; Zanchetta, P.; Rabbeni, R.; Pucci, M. Model Predictive Control for Shunt Active Filters With Fixed Switching Frequency. IEEE Trans. Ind. Appl. 2017, 53, 296–304. [Google Scholar] [CrossRef]
  33. Taher, S.A.; Alaee, M.H.; Arani, Z.D. Model predictive control of PV-based shunt active power filter in single phase low voltage grid using conservative power theory. In Proceedings of the 2017 8th Power Electronics, Drive Systems & Technologies Conference (PEDSTC), Mashhad, Iran, 7–14 January 2017; pp. 253–258. [Google Scholar] [CrossRef]
  34. Nademi, H.; Burgos, R.; Soghomonian, Z. Power Quality Characteristics of a Multilevel Current Source with Optimal Predictive Scheme From More-Electric-Aircraft Perspective. IEEE Trans. Veh. Technol. 2018, 67, 160–170. [Google Scholar] [CrossRef]
  35. Saribulut, L.; Teke, A.; Tümay, M. Artificial neural network-based discrete-fuzzy logic controlled active power filter. IET Power Electron. 2014, 7, 1536–1546. [Google Scholar] [CrossRef]
  36. Radzi, M.A.M.; Rahim, N.A. Neural Network and Bandless Hysteresis Approach to Control Switched Capacitor Active Power Filter for Reduction of Harmonics. IEEE Trans. Ind. Electron. 2009, 56, 1477–1484. [Google Scholar] [CrossRef]
  37. Firouzjah, K.G.; Sheikholeslami, A.; Karami-Mollaei, M.R.; Khaleghi, M. A new harmonics detection method for shunt active filter based on wavelet transform. J. Appl. Sci. Res. 2008, 4, 1561–1568. [Google Scholar]
  38. Soljan, Z.; Bialostocka, P. Part 1: Decomposition of the Load′s Current in Three-phase Four-wire Systems Supplied with Asymmetrical Sinusoidal Voltage Following Currents′ Physical Components Theory. Autom. Elektr. Zakłócenia 2020, 11, 34–42. [Google Scholar] [CrossRef]
  39. Soljan, Z.; Holdynski, G.; Zajkowski, M. Decomposition of the load′s current supplied from a sinusoidal and asymmetrical voltage source in accordance with the Currents’ Physical Components (CPC) Theory. In Proceedings of the 2019 20th International Scientific Conference on Electric Power Engineering (EPE), Kouty nad Desnou, Czech Republic, 15–17 May 2019; Institute of Electrical and Electronics Engineers (IEEE): Picataway, NJ, USA, 2019; pp. 1–6. [Google Scholar]
  40. Soljan, Z.; Bialostocka, P. Part 2: Reactive Compensation of the Linear Load in Three-Phase Four-Wire Systems Supplied with Asymmetrical Sinusoidal Voltage Following Currents’ Physical Components Theory. Autom. Elektr. Zakłócenia 2020, 11, 52–64. [Google Scholar] [CrossRef]
  41. Sołjan, Z.; Hołdyński, G.; Zajkowski, M. Balancing reactive compensation at three-phase four-wire systems with a sinusoidal and asymmetrical voltage source. Bull. Pol. Acad. Sci. Tech. Sci. 2020, 68, 71–79. [Google Scholar] [CrossRef]
  42. Holdynski, G.; Soljan, Z.; Zajkowski, M. CPC-based reactive compensation at three-phase four-wire systems powered from asymmetric sinusoidal voltage source. In Proceedings of the 2019 20th International Scientific Conference on Electric Power Engineering (EPE), Kouty nad Desnou, Czech Republic, 15–17 May 2019; Institute of Electrical and Electronics Engineers (IEEE): Picataway, NJ, USA, 2019; pp. 1–6. [Google Scholar]
  43. Czarnecki, L.S. Minimisation of distortion power of nonsinusoidal source applied to linear loads. Proc. IEE 1981, 128, 208–210. [Google Scholar]
  44. Czarnecki, L.S. Minimisation of unbalanced and reactive currents in three-phase asymmetrical circuits with nonsinusoidal voltage. IEE Proc. B 1992, 139, 347–354. [Google Scholar]
  45. Czarnecki, L.S. Minimization of reactive power under nonsinusoidal conditions. IEE Trans. Instrum. Meas. 1987, 36, 18–22. [Google Scholar]
  46. Czarnecki, L.S.; Bhattarai, P.D. CPC-based reactive compensation of linear loads supplied with asymmetrical nonsinusoidal voltage. Prz. Elektrotechniczny 2019, 1, 3–10. [Google Scholar]
  47. Czarnecki, L.S.; Bhattarai, P. Power and compensation in circuits with distorted and asymmetrical current and voltage waveforms. Vol. 8, Power and compensation of circuits with asymmetrical supply voltage. Aut. Elektr. Zakł. 2015, 6, 8–18. (In Polish) [Google Scholar] [CrossRef]
  48. Buła, D.; Grabowski, D.; Lewandowski, M.; Maciążek, M.; Piwowar, A. Software Solution for Modeling, Sizing, and Allocation of Active Power Filters in Distribution Networks. Energies 2020, 14, 133. [Google Scholar] [CrossRef]
  49. Czarnecki, L.S. Two algoritkms of the fundamental harmonic complex RMS value calculation. Arch. Fur Elektr. 1992, 73, 163–168. [Google Scholar]
  50. Patrascu, A.; Popescu, M.; Suru, V.; Patrascu, A. CPC theory implementation for active filtering and its limits. In Proceedings of the 2012 International Conference on Applied and Theoretical Electricity (ICATE), Craiova, Romania, 25–27 October 2012; pp. 1–6. [Google Scholar] [CrossRef]
  51. Rahmani, S.; Hamadi, A.; Al-Haddad, K. A Lyapunov-Function-Based control for a three-phase shunt hybrid active filter. IEEE Trans. Ind. Electron. 2012, 59, 1418–1429. [Google Scholar]
  52. Asiminoael, L.; Blaabjerg, F.; Hansen, S. Detection is key—Harmonic detection methods for active power filter applications. IEEE Ind. Appl. Mag. 2007, 13, 22–33. [Google Scholar]
Figure 1. Linear time-invariant load (LTI) load supplied from four-wire line.
Figure 1. Linear time-invariant load (LTI) load supplied from four-wire line.
Energies 14 01815 g001
Figure 2. Three-phase four-wire system with LTI load.
Figure 2. Three-phase four-wire system with LTI load.
Energies 14 01815 g002
Figure 3. The three-phase voltage waveforms at the terminals of the load.
Figure 3. The three-phase voltage waveforms at the terminals of the load.
Energies 14 01815 g003
Figure 4. The three-phase current waveforms of the load.
Figure 4. The three-phase current waveforms of the load.
Energies 14 01815 g004
Figure 5. The analyzed system with connected minimizing balancing compensators (MBC).
Figure 5. The analyzed system with connected minimizing balancing compensators (MBC).
Energies 14 01815 g005
Figure 6. The three-phase current waveforms of the load after minimizing balancing compensation.
Figure 6. The three-phase current waveforms of the load after minimizing balancing compensation.
Energies 14 01815 g006
Figure 7. The analyzed system with connected the minimizing balancing compensators (MBC) and the active power filter (APF).
Figure 7. The analyzed system with connected the minimizing balancing compensators (MBC) and the active power filter (APF).
Energies 14 01815 g007
Figure 8. The analyzed theoretical illustration which connected the minimizing balancing compensators (MBC) and the active power filter (APF).
Figure 8. The analyzed theoretical illustration which connected the minimizing balancing compensators (MBC) and the active power filter (APF).
Energies 14 01815 g008
Figure 9. The waveform of three-phase referent current in cross-section X-X’—Equation (43).
Figure 9. The waveform of three-phase referent current in cross-section X-X’—Equation (43).
Energies 14 01815 g009
Figure 10. The waveform of three-phase nonsinusoidal active current in cross-section Z-Z’.
Figure 10. The waveform of three-phase nonsinusoidal active current in cross-section Z-Z’.
Energies 14 01815 g010
Figure 11. The waveform of the reference current in cross-section Y-Y’.
Figure 11. The waveform of the reference current in cross-section Y-Y’.
Energies 14 01815 g011
Figure 12. The waveform of three-phase referent current in cross-section X-X’—Equation (42).
Figure 12. The waveform of three-phase referent current in cross-section X-X’—Equation (42).
Energies 14 01815 g012
Figure 13. The waveform of three-phase sinusoidal active current in cross-section Z-Z’.
Figure 13. The waveform of three-phase sinusoidal active current in cross-section Z-Z’.
Energies 14 01815 g013
Table 1. List of phase impedances with division of resistance, inductive, and capacitive reactance values.
Table 1. List of phase impedances with division of resistance, inductive, and capacitive reactance values.
LineResistance R [Ω]Inductive Reactance XL [Ω]Capacitive reactance XC [Ω]
R11-
S21.52
TR1 = 0.5, R2 = 51-
Table 2. List of the values of the phase voltages and the values of the line currents.
Table 2. List of the values of the phase voltages and the values of the line currents.
Harmonic OrderPhase RPhase SPhase T
U n   [ V ]
1st 230 e j 5 215 e j 110 220 e j 99
3rd 20 e j 0 33 e j 10 30 e j 20
5th 40 e j 6 20 e j 114 27 e j 133
7th 19 e j 7 10 e j 100 30 e j 124
I n [A]
1st 162.64 e j 40 104.47 e j 169 220 e j 45.9
3rd 6.33 e j 71.6 15.04 e j 0.9 12.36 e j 71.9
5th 7.85 e j 84.7 9.64 e j 109.6 7.99 e j 175
7th 2.69 e j 74.9 4.91 e j 103 7.61 e j 89.9
Table 3. List of the values of the individuals admittances.
Table 3. List of the values of the individuals admittances.
Parameter [S]Harmonic Order
1st3rd5th7th
G e 0.445
G b n 0.454 0.318 0.152 0.180
B b n 0.572 0.205 0.171 0.133
Y d n 4.475 10 3 e j 178.8 0.056 e j 151.6 0.099 e j 16.8 0.067 e j 27.1
Y n p ----
Y n n 0.107 e j 142.8 0.158 e j 171.8 0.159 e j 161.5 0.157 e j 156.7
Y n z 0.193 e j 45.3 0.091 e j 97.1 0.116 e j 119.9 0.125 e j 127.5
Table 4. List of the values of the currents’ components.
Table 4. List of the values of the currents’ components.
Current [A]Harmonic OrderTotal Value [A]
1st3rd5th7th
i a 170.75 21.73 23.22 16.40 174.46
i s 3.81 6.21 15.30 9.75 19.55
i r 219.82 10.04 8.94 4.91 220.29
i u p 18.55 6.73 4.68 0.83 20.30
i u n 28.08 5.20 2.24 3.86 28.91
i u z 76.42 2.23 6.98 2.29 76.81
i 292.86 20.47 14.78 9.45 294.09
Table 5. List of the values of the reactive elements for the minimizing balancing compensators.
Table 5. List of the values of the reactive elements for the minimizing balancing compensators.
StructureRST
STAR-MBC C [ μ F ] L [ mH ] C [ μ F ] L [ mH ] C [ μ F ] L [ mH ]
95.21 92.75 203.91 43.38 402.29 21.93
RSSTTR
DELTA-MBC C [ μ F ] L [ mH ] C [ μ F ] L [ mH ] C [ μ F ] L [ mH ]
53.96 163.83 21.2 28.88 307.93
Table 6. List of the values of the three-phase currents’ components after the minimizing balancing compensation.
Table 6. List of the values of the three-phase currents’ components after the minimizing balancing compensation.
Current [A]Harmonic OrderTotal Value [A]
1st3rd5th7th
i a 170.75 21.73 23.22 16.40 174.46
i s 3.81 6.21 15.30 9.75 19.55
i r 3.56 12.11 10.96 6.45 17.92
i u p 0.82 7.23 3.54 0.95 8.14
i u n 2.04 5.70 2.04 4.11 7.59
i u z 2.89 2.43 7.23 2.29 8.47
i 170.87 27.35 31.02 20.71 177.02
Table 7. List of the RMS value of the currents of the lines before compensation and after using the passive and active compensation.
Table 7. List of the RMS value of the currents of the lines before compensation and after using the passive and active compensation.
PhaseRMS Value [A]
BEFOREMBCAPF (43)APF (46)
R162.97105.93104.95102.25
S106.1099.3798.3495.58
T220.62101.19100.6797.80
i 294.09177.02174.46170.75
N180.7562.9600
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sołjan, Z.; Hołdyński, G.; Zajkowski, M. CPC-Based Minimizing of Balancing Compensators in Four-Wire Nonsinusoidal Asymmetrical Systems. Energies 2021, 14, 1815. https://doi.org/10.3390/en14071815

AMA Style

Sołjan Z, Hołdyński G, Zajkowski M. CPC-Based Minimizing of Balancing Compensators in Four-Wire Nonsinusoidal Asymmetrical Systems. Energies. 2021; 14(7):1815. https://doi.org/10.3390/en14071815

Chicago/Turabian Style

Sołjan, Zbigniew, Grzegorz Hołdyński, and Maciej Zajkowski. 2021. "CPC-Based Minimizing of Balancing Compensators in Four-Wire Nonsinusoidal Asymmetrical Systems" Energies 14, no. 7: 1815. https://doi.org/10.3390/en14071815

APA Style

Sołjan, Z., Hołdyński, G., & Zajkowski, M. (2021). CPC-Based Minimizing of Balancing Compensators in Four-Wire Nonsinusoidal Asymmetrical Systems. Energies, 14(7), 1815. https://doi.org/10.3390/en14071815

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop