Next Article in Journal
Practical Analysis and Design of a Battery Management System for a Grid-Connected DC Microgrid for the Reduction of the Tariff Cost and Battery Life Maximization
Previous Article in Journal
Control Scheme and Power Electronics Architecture for a Wirelessly Distributed and Enabled Battery Energy Storage System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Formulation of the Phasors of Apparent Harmonic Power: Application to Non-Sinusoidal Three-Phase Power Systems

1
Departamento de Ingeniería Eléctrica, Universitat Politècnica de València, Plaza Ferrandiz y Carbonell s/n, 03801 Alcoy, Spain
2
Departamento de Ingeniería Mecánica y Materiales, Universitat Politècnica de València, Plaza Ferrandiz y Carbonell s/n, 03801 Alcoy, Spain
*
Author to whom correspondence should be addressed.
Energies 2018, 11(7), 1888; https://doi.org/10.3390/en11071888
Submission received: 27 June 2018 / Revised: 12 July 2018 / Accepted: 16 July 2018 / Published: 19 July 2018
(This article belongs to the Section F: Electrical Engineering)

Abstract

:
In this work, the expression of the phasor of apparent power of harmonic distortion is formulated in the time domain. Applying this phasor along with the phasor of apparent unbalance power allows us to obtain a new set of phasors that include all of the inefficient power components appearing in the transfer of energy in non-linear and unbalanced systems. In this manner, a new model of inefficient power in electrical systems is developed. For each voltage harmonic of order ‘m’ and current harmonic of order ‘n’, a phasor of harmonic apparent power is obtained. Accuracy in the determination of the total apparent power of a system depends on the number of harmonics considered. Each phasor of apparent harmonic power is formed from six mutually orthogonal parameters or components that are calculated from the harmonic voltages at the nodes of the network and the circulating harmonic currents. To demonstrate the validity of the proposed formulation, a four-wire non-linear system formed by two nodes is assessed.

1. Introduction

In addition to useful power for the transfer of energy (or active power), the power generated by a three-phase electrical system includes non-useful components that must be taken into account in the analysis of energy transfer [1]. Phase differences between voltage and current result in reactive power [2], and an unbalanced three-phase electrical system produces a so-called power of unbalance [3]. In addition, such electrical systems are designed to work with voltages and sinusoidal currents, and the use of non-linear loads and/or variants over time distorts the voltage and current waveforms, resulting in the creation of harmonic components. The analysis of the generation and propagation of harmonic components through an electrical system is called harmonic power flow analysis [4].
These types of non-useful power can be critically detrimental to an electrical system through their effects on the operation of protection and measuring equipment, which can result in unexpected opening and erroneous measurement [5,6]. In addition, increased currents can enhance power losses in transmission lines [7] and in motors can cause winding overheating; loss of thermal insulation, copper and iron; reduction in the motor torque; decrease in performance, etc. [8]. In generators, non-useful components can hinder automatic synchronization with the network, and in neutral conductors they produce dangerous overheating because they add harmonic currents of the third-order (or integer multiples thereof). The physical significance of non-useful power components has been widely discussed in the scientific literature [9,10,11,12].
Non-linear loads are increasingly used in general applications, including switch mode power supplies, motors and variable speed drives, photocopiers, personal computers, laser printers, fax machines, battery chargers, etc. There is also currently a high amount of government support for the development of power generation using photovoltaic technology. Solar plants are considered to be an important source of harmonics as a result of the power electronics technology used in the electricity generation process [13,14,15]. The non-linear loads and solar energy sourcing has therefore increased the overall harmonic power present in electrical systems.
Harmonic power has also been widely discussed in the literature [16,17,18], with most recent studies focusing on the development of active filters to mitigate the effects of harmonic components [19,20,21,22,23,24]. Other research focuses on the development of techniques to assess and differentiate the harmonic contribution of customers and networks [25,26]. In these contributions, the harmonic power of the system is not evaluated phasorially. Currently, electrical power expressions developed in IEEE Std. 1459–2010 [10], obtained in the time domain, are accepted as valid, although approaches in the frequency domain are also possible [27]. Under the time-domain standard, harmonic current and voltage distortion rates (THDi and THDv, respectively) are established as proportions of the fundamental components and are expressed as percentages of the respective components.
In 2016, the authors of this paper defined the unbalance power phasor [28] and developed an equivalent circuit for the analysis of unbalance power in three-wire electrical systems [29]. The unbalance power and apparent power phasors of the positive sequence can be combined to produce the apparent power phasor of an unbalanced linear electrical system; this phasor has a modulus equivalent to that obtained by Buchholz [30] and IEEE Std. 1459–2010 [10].
This paper extends the method used in the formulation of the unbalance power phasor [28] to non-linear systems for calculating the harmonic distortion power. This application to non-linear power systems is not immediate. In this case, there are voltage and current harmonics of the same order, as well as of a different order. This gives rise to different harmonic powers, depending on the harmonic order of voltage and current that is considered. In this work, these powers are formulated phasorially.
The main advantage of the proposed method is that it allows analysis by phasor components. This allows us to obtain a new set of phasors that encompasses all balanced, unbalanced, sinusoidal or non-sinusoidal inefficient power components appearing in the transfer of energy in an electrical system in an approach that maintains all of the advantages of the phasor formulation.
The difference compared to the methods listed in the literature, mainly the IEEE Std. 1459–2010 [9], is that it only determines the harmonic powers in absolute values. This value matches with the module of the phasor proposed in this work. Moreover, another added advantage, presented by the phasor of apparent harmonic power, is that it allows us to know the contribution of each load to each of the mentioned harmonic powers at any point in the network.
The proposed procedure is not intended to replace the procedure established in IEEE Std. 1459-2010 but instead represents a new approach to understanding the inefficient power in an electrical system from a time-domain-based perspective.
The rest of this paper is organized as follows. In Section 2, the production of instantaneous power in non-linear systems is reviewed. In Section 3, the phasor of apparent harmonic power is formulated for a three-phase system with balanced voltages and non-linear loads. In Section 4, the apparent harmonic power phasor for a three-phase system with unbalanced voltages and non-linear loads is analyzed. In Section 5, the total apparent power of an electrical system and its modulus are analyzed. To aid in the understanding of the proposed calculation method and its application, Section 6 examines a case involving a system with unbalanced voltages and non-linear loads.

2. Review of Instantaneous Power in a Non-Sinusoidal System

The voltage and instantaneous current in a non-sinusoidal single-phase system are given by the following equations:
v ( t ) = V 0 + 2   m V m sin ( m ω t + α m )
i ( t ) = I 0 + 2   n I n sin ( n ω t + β n )
where
  • m is the harmonic order of the voltage and n is the harmonic order of the current.
  • V 0 and I 0 are, respectively, the voltage and current constant terms. In alternating-current power systems, these values are generally very small and, therefore, will not be considered further in this work.
  • V m and I n are the root mean square (RMS) values for the harmonic voltage and harmonic current of orders m and n, respectively.
  • α m and β n are the phase angle harmonic voltage and phase angle harmonic current of orders m and n, respectively.
The instantaneous power is given by the following equation:
p ( t ) = v ( t )   i ( t ) =   p a +   p q
where the first term, p a , is the instantaneous active power, which is decomposed into three instantaneous components:
p a = p a 1 + p a h + p a ( m , n )
Here, p a 1 is the fundamental instantaneous active power, which is calculated as
p a 1 = V 1 I 1   cos θ 1 [ 1 cos ( 2 ω t + 2 α 1 ) ]
This expression has two terms: the fundamental active power P 1 = V 1 I 1   cos θ 1 , and the oscillating component P 1 cos ( 2 ω t + 2 α 1 ) . The first term is equal to the average value of p a 1 and determines the active power caused by the fundamental harmonics of voltage and current. The second term, which has an average value of zero, is always present when net energy is transferred to the load; however, it does not cause power loss in the conductors.
The second term in (4), p a h , is the harmonic instantaneous active power of order h for h ≠ 1, and is calculated as
p a h = h 1 V h I h   cos θ h [ 1 cos ( 2 h ω t + 2 α h ) ]  
This term is caused by harmonic voltages of order h and harmonic currents of order h, that is, for harmonics of equal voltage and current, and has a value of zero for a linear electrical system.
The first term in (6) represents the sum of the harmonic active power components P h = V h I h   cos θ h . The sum h 1 P h is equal to the average value of p a h and the sum of P 1 and h 1 P h is the total active power of a non-sinusoidal electrical system. The second term of (6) represents the sum of oscillating harmonic components of amplitude P h and a frequency of 2h. The average value of the second term is zero, and the term does not correspond to power loss in the conductors.
The third term in (4), p a ( m , n ) , is the harmonic instantaneous active power, which is caused by the harmonic voltage of order m and the harmonic current of order n, where mn. This term is calculated as
p a ( m , n ) = m n m , n = 1 V m I n   cos θ n m [ cos ( m n ) ω t cos ( ( m + n ) ω t + 2 α m ) ]  
The average value of p a ( m , n ) is zero.
The second term in (3), p q , is the instantaneous reactive power, which is also decomposed into three instantaneous components:
p q = p q 1 + p q h + p q ( m , n )
where
  • p q 1 is the fundamental instantaneous reactive power, which is calculated from (9).
  • p q h is harmonic instantaneous reactive power of order h ≠ 1, which is calculated from (10). This power component is caused by the harmonic voltage and harmonic current of orders h and h, respectively.
  • p q ( m , n ) is harmonic instantaneous reactive power, which is caused by the harmonic voltage of order m and the harmonic current of order n, where mn. This power is calculated from (11).
    p q 1 = V 1 I 1   sin θ 1 sin ( 2 ω t + 2 α 1 )
    p q h = h 1 V h I h   sin θ h sin ( 2 h ω t + 2 α h )
    p q ( m , n ) = m n m , n = 1 V m I n   sin θ n m [ sin ( m n ) ω t + sin ( ( m + n ) ω t + 2 α m ) ]
In a non-sinusoidal three-phase system, the instantaneous power and its components are expressed by the following equation:
p ( t ) = z = a , b , c v z i z =   z = a , b , c p a z + z = a , b , c p q z
For each of the phases (a, b and c), the instantaneous active power p a z and the instantaneous reactive power p q z are calculated from (4)–(7) and (8)–(11), respectively.

3. Harmonic Power in a Non-Linear Three-Phase Power System with Non-Sinusoidal Balanced Voltages

Figure 1 shows a non-linear load that is connected to an electrical system with non-sinusoidal balanced voltages.
where
v a =   m = 1 v a m i a =   n = 1 i a n
v b =   m = 1 v b m i b =   n = 1 i b n
v c =   m = 1 v c m i c =   n = 1 i c n .
where v a m , v b m and v c m are line-to-neutral harmonic voltages of order m and i a n , i b n and i c n are harmonic currents of order n.
These voltages can be expressed in terms of the following symmetrical components:
v m + = 2   m = 1 V m + sin ( m ω t + α m + )
v m = 2   m = 1 V m sin ( m ω t + α m )
v m 0 = 2   m = 1 V m 0 sin ( m ω t + α m 0 )
Similarly, the symmetric current components are
i n + = 2   n = 1 I n + sin ( n ω t + β n + )
i n = 2   n = 1 I n sin ( n ω t + β n )
i n 0 = 2   n = 1 I n 0 sin ( n ω t + β n 0 )
where
  • V m + and I n + are the RMS values for the positive-sequence harmonic voltage and current of orders m and n, respectively.
  • V m and I n are the RMS values for the negative-sequence harmonic voltage and current of orders m and n, respectively.
  • V m 0 and I n 0 are the RMS values for the zero-sequence harmonic voltage and current of orders m and n, respectively.
  • α m + , α m and α m 0   are, respectively, the angles of positive-, negative- and zero-sequence of the harmonic voltage of order m.
  • β n + , β n and β n 0 are, respectively, the angles of positive-, negative- and zero-sequence of the harmonic current of order n.
As the voltages are balanced, the negative and zero-sequence voltages ( v m and v m 0 ) become zero. The positive-sequence voltages for each phase (a, b and c) are determined as follows:
  v a m + = v m +       v b m + = a 2 v m +       v c m + = a   v m +
where a = 1   e j 120 .
Similarly, the phase-specific positive-, negative- and zero-sequence currents are
i a n + = i n +       i b n + = a 2 i n +       i c n + = a   i n +
i a n = i n       i b n = a   i n       i c n = a 2   i n
i a n 0 = i b n 0 = i c n 0 = i n 0
Under these conditions, the instantaneous power of the system is given by
p ( t ) = p + + p + + p 0 +
where
p + = m = 1 n = 1 z = a , b , c V z m + I z n +
p + = m = 1 n = 1 z = a , b , c V z m + I z n
p 0 + = m = 1 n = 1 z = a , b , c V z m + I z n 0
Here,
  • p + is the instantaneous power caused by positive-sequence voltages and currents.
  • p + is the instantaneous power caused by positive-sequence voltages and negative-sequence currents.
  • p 0 + is the instantaneous power caused by positive-sequence voltages and zero-sequence currents.
Expressions of instantaneous power for voltage and current harmonics of order m and n, respectively, can then be obtained as
p + ( m , n ) = 3   V m + I n + cos ( ( m n ) ω t + θ n + m + )
p ( m , n ) + = 3   V m + I n cos ( ( m + n ) ω t + 2 α m + θ n m + )
p 0 ( m , n ) + = 3   V m + I n 0 cos ( ( m + n ) ω t + 2 α m + θ n 0 m + )
where
  • θ n + m + is the angle between the positive-sequence harmonic voltage of order m and the positive-sequence harmonic current of order n.
  • θ n m + is the angle between the positive-sequence harmonic voltage of order m and the negative-sequence harmonic current of order n.
  • θ n 0 m + is the angle between the positive-sequence harmonic voltage of order m and the zero-sequence harmonic current of order n.
Equation (23) can be decomposed into the following two terms:
P n + m + = 3   V m + I n + cos θ n + m +
Q n + m + = 3   V m + I n + sin θ n + m +
where P n + m + and Q n + m + are, respectively, the positive-sequence active power and the positive-sequence reactive power for a harmonic voltage of order m and a harmonic current of order n.

3.1. Harmonic Parameters A n m , B n m , C n m and D n m

A three-phase four-wire star-connected system with a single-phase load in the A-phase with balances and sinusoidal voltages will satisfy the following conditions:
I n = I n 0 =   I n + =   1 3 I a n
V a m = V a m + = V m +
α m + = α a m = α a m +
θ n + m + = θ n m + = θ n 0 m + = θ a n a m = θ a n a m +
Substituting these equalities into (24) and (25) produces the following expression:
p ( m , n ) + + p 0 ( m , n ) + = 2   V a m + I a n cos ( ( m + n ) ω t + 2 α a m + θ a n a m + )
If we proceed in this manner in connecting the load to Phase B and then to Phase C, applying the superposition theorem results in the following expression:
p ( m , n ) + + p 0 ( m , n ) + = 2 z = a , b , c V z m + I z n cos ( ( m + n ) ω t + 2 α z m + θ z n z m + )
which can be decomposed into A n m ( t ) , B n m ( t ) , C n m ( t ) and D n m ( t ) as follows:
p ( m , n ) + + p 0 ( m , n ) + = A n m ( t ) + B n m ( t ) + C n m ( t ) + D n m ( t )
where
A n m ( t ) = 2 z = a , b , c V z m + I z n cos θ z n z m + cos 2 α z m + cos ( m + n ) ω t
B n m ( t ) = 2 z = a , b , c V z m + I z n cos θ z n z m + sin 2 α z m + cos ( m + n ) ω t
C n m ( t ) = 2 z = a , b , c V z m + I z n sin θ z n z m + sin 2 α z m + cos ( m + n ) ω t
D n m ( t ) = 2 z = a , b , c V z m + I z n sin θ z n z m + cos 2 α z m + sin ( m + n ) ω t
These instantaneous parameters are sinusoidal waveforms with zero average value and RMS values given by:
A n m = 2 z = a , b , c V z m + I z n cos θ z n z m + cos 2 α z m +
B n m = 2 z = a , b , c V z m + I z n cos θ z n z m + sin 2 α z m +
C n m = 2 z = a , b , c V z m + I z n sin θ z n z m + sin 2 α z m +
D n m = 2 z = a , b , c V z m + I z n sin θ z n z m + cos 2 α z m +

3.2. Phasor of Harmonic Apparent Power S n m

For each harmonic voltage of order m ≠ 1 and harmonic current of order n ≠ 1, the phasor of harmonic apparent power S n m can be defined as:
S n m = u n ( A ) m   A n m + u n ( B ) m   B n m + u n ( C ) m   C n m + u n ( D ) m   D n m + u n ( X ) m   P n + m + + u n ( Y ) m   Q n + m +
where
  • u n ( A ) m , u n ( B ) m , u n ( C ) m and u n ( D ) m are the unit vectors associated with parameters A n m , B n m , C n m and D n m , respectively (see (35)–(38)).
  • u n ( X ) m and u n ( Y ) m are the unit vectors associated with the positive-sequence active power P n + m + and the positive-sequence reactive power Q n + m + (see (26) and (27)).
As these unit vectors are all mutually orthogonal, the moduli of harmonic apparent power for a harmonic voltage of order m and harmonic current of order n is given by:
| S n m | = A n m 2 + B n m 2 + C n m 2 + D n m 2 + P n + m + 2 + Q n + m + 2
For m = n = 1, the phasor of fundamental apparent power is obtained as:
  S 1 1 = u 1 ( A ) 1   A 1 1 + u 1 ( B ) 1   B 1 1 + u 1 ( C ) 1   C 1 1 + u 1 ( D ) 1   D 1 1 + u 1 ( X ) 1   P 1 + 1 + + u 1 ( Y ) 1   Q 1 + 1 +
This expression is the same as that defined in [22]. In this phasor, the positive-sequence active and reactive power components, as well as the unbalanced power caused by unbalanced currents, are all included.

3.3. Application of S n m on a Bus with Loads Connected in Parallel

Figure 2 shows two non-linear three-phase loads in a star configuration in which they are connected to bus ‘i’ of an electric power system.
Under these conditions, for a harmonic voltage of order m and harmonic current of order n, the phasor of harmonic apparent power in bus ‘i’ is given by the sum of the individual phasors of each load connected to the bus:
S ( b u s   i ) n m =   S ( l o a d   1 ) n m + S ( l o a d   2 ) n m

3.4. Application of S n m between Two Buses in a System

Figure 3 shows two buses in a system (i and j) linked by a power line with impedance Z i j .
Under these conditions, for a harmonic voltage of order m and harmonic current of order n the phasor of harmonic apparent power in the power line between the two buses is determined by subtracting the phasor in bus ‘j’ from that in bus ‘i’:
S ( b u s   l i n e ) n m =   S ( b u s   i ) n m S ( b u s   j ) n m

4. Harmonic Power in a Non-Linear Three-Phase Power System with Non-Sinusoidal Unbalanced Voltages

In Section 3, we assumed that the voltages were non-sinusoidal but balanced in each of the phases. In a real three-phase non-linear system, the voltages are normally both non-sinusoidal and unbalanced; in this case, the expression of the apparent power module defined by Buchholz in Reference [24] can be extended for voltage and current harmonics of order m and n, respectively, to obtain:
| S T n m | = 3   ( V m + 2 + V m 2 + V m 0 2 )   ( I n + 2 + I n 2 + I n 0 2 )
Substituting the voltage unbalance factors of order m, δ m = V m / V m + and δ m 0 = V m 0 / V m + , into (44) produces the following expression:
| S T n m | = 1 + F u m 2   9   V m + 2     ( I n + 2 + I n 2 + I n 0 2 )
where F u m = δ m 2 + δ m 0 2 and is referred to as the global harmonic voltage unbalance factor of order m. The expression under the second square root in (45) corresponds to the phasor magnitude apparent harmonic power defined in (39); this suggests that the modulus of the apparent power | S n m | for a harmonic voltage of order m and a harmonic current of order n at any point of the system must be corrected and multiplied by 1 + F u m 2 as follows:
| S T n m | = | S n m |   1 + F u m 2

5. Total Harmonic Apparent Power and Total Apparent Power

In any node of a non-linear system, the number of harmonics of voltage and current that are considered will determine the number of phasors of harmonic apparent power. As these phasors are all mutually orthogonal, they do not sum arithmetically. In any node, the module of total harmonic apparent power S T H is calculated from either of the following two expressions:
| S T H | = m 1 n 1 | S T n m | 2
| S T H | = m 1 n 1 | S n m | 2   ( 1 + F u m 2 )
These values coincide with those obtained in IEEE Std. 1459-2010 [9] and by Buchholz [24]. Equation (47) or (48) can be alternatively expressed as a function of the fundamental harmonic of voltage and current (m = n = 1) in either of the two following manners:
| S T H | = H   | S T 1 1 |
| S T H | = H   | S 1 1 |   ( 1 + F u 1 2 )
where
H =   T H D I 2 + T H D V 2 + T H D I 2   T H D V 2
The values T H D I and T H D V are, respectively, the rates of harmonic distortion of the currents and voltages of the three-phase system. Symmetrical components are determined by:
T H D I =   n 1 ( I n + 2 + I n 2 + I n 0 2 ) I 1 + 2 + I 1 2 + I 10 2
T H D V =   m 1 ( V m + 2 + V m 2 + V m 0 2 ) V 1 + 2 + V 1 2 + V 10 2
In terms of phase components, they are determined by:
T H D I =   n 1 ( I a n 2 + I b n 2 + I c n 2 ) I a 1 2 + I b 1 2 + I c 1 2
T H D V =   m 1 ( V a m 2 + V b m 2 + V c m 2 ) V a 1 2 + V b 1 2 + V c 1 2
At any point in the system, the module of total apparent power S T is given by:
| S T | = | S T 1 1 | 2 + | S T H | 2

6. Practical Application

In this section, a practical case study for checking all of the concepts discussed in the previous sections is developed. Figure 4 shows a four-wire, two-bus electrical system with two unbalanced three-phase non-linear loads and an unbalanced three-phase linear load. All of the loads are modelled at a constant impedance. To simulate loads 2 and 3, we use single-phase full-wave rectifiers (see Figure 5). The values of Ra, Rb and Rc for loads 2 and 3 are listed in Table 1.
Using the ‘PSPICE’ analysis software (version 9.2), the harmonic line-to-neutral voltages in the buses (1 and 2) and harmonic currents circulating in the loads are obtained. These magnitudes are displayed in Table 2, Table 3 and Table 4. The harmonic currents supplied by the generator and the circulating currents along the line between buses 1 and 2 are easily deduced from the load currents. Here, we consider only harmonics 1, 3, 5 and 7 of voltage and current.
Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10 list the results obtained using the equations formulated in the preceding sections, where:
  • the values of P n + m + and Q n + m + are calculated from (26), (27) and those of A n m , B n m , C n m and D n m are calculated from (35)–(38). The ‘line 1–2’ results in Table 10 are calculated from the difference between the corresponding node 1 and 2 values.
  • the values of S n m are calculated from (40).
Table 11 compares the following total apparent power values: S 1 + 1 + (positive-sequence fundamental apparent power), which is calculated from P 1 + 1 + and Q 1 + 1 + ; S T 1 1 . (fundamental apparent power), which is calculated from (45); S T H (total harmonic apparent power), which is calculated from (47); and S T (total apparent power), which is calculated from (55). These values are calculated from (48) and are the same as those obtained using IEEE Std. 1459–2010.

7. Conclusions

In a power system, there are several inefficient power components within the energy balance, and reducing or eliminating them requires that they be properly quantified. Three-phase sinusoidal systems are exclusively affected by fundamental voltage and current components, with imbalances in their magnitudes at fundamental frequencies leading to the appearance of unbalance power. This power component was analyzed by the authors in a previous work [22] through the application of the unbalance power phasor. In non-sinusoidal systems, the voltages and currents contain harmonic components other than the fundamental component that generate a harmonic power component that is also inefficient.
In this study, the authors extended the method they applied in analyzing unbalance power to the analysis of harmonic power by formulating the phasor of apparent harmonic power for the voltage harmonic of order ‘m’ and current harmonic of order ‘n’. The number of phasors depends on the number of voltage and current harmonics considered, and phasors of voltage and current harmonics of the same and different orders are considered. Each phasor of apparent harmonic power is formed from six parameters or components that are mutually orthogonal; these components are calculated from the harmonic voltages at the nodes of the network and from the circulating harmonic currents. For a node with several loads connected in parallel, the components of each resultant phasor are calculated as the arithmetic sum of the individual components of the loads. For nodes connected in serial, the components are calculated through arithmetic subtraction of the components. An expression of the modules of the total harmonic power and total apparent power at any point was also proposed. The results of these formulations were shown to coincide with those obtained using IEEE Std. 1459–2010.
The proposed set of phasors of harmonic power can be used in conjunction with the phasor of unbalance power to analyze all of the inefficient power components appearing in the transfer of energy in an electrical system, regardless of whether they are balanced or unbalanced or sinusoidal or non-sinusoidal. To validate the applicability of the proposed expressions and improve our understanding of them, we undertook practical case studies of four-wire systems with two nodes, unbalanced voltages and loads, and non-linear loads.

Author Contributions

Conceptualization, P.A.B., R.M. and J.M.D.; Methodology, P.A.B., R.M.-M., J.M.D. and R.M.; Validation, P.A.B., R.M.-M. and R.M.; Formal Analysis, P.A.B. and R.M.-M.; Investigation, P.A.B., R.M.-M., J.M.D. and R.M.; Resources, P.A.B. and R.M.-M.; Data Curation, P.A.B. and R.M.-M.; Writing—Original Draft Preparation, J.M.D. and R.M.; Writing—Review & Editing, R.M.-M., J.M.D. and R.M.; Visualization, P.A.B. and R.M.-M.; Supervision, J.M.D., R.M. and M.J.R.; Project Administration, J.M.D.

Funding

This work is supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and the European Regional Development Fund (ERDF) under Grant ENE2015-64087-C2-2R.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Emanuel, A.E. On the definition of power factor and apparent power in unbalanced polyphase circuits with sinusoidal voltage and currents. IEEE Trans. Power Deliv. 1993, 8, 841–852. [Google Scholar] [CrossRef]
  2. Jeon, S.J. Definitions of apparent power and power factor in a power system having transmission lines with unequal resistances. IEEE Trans. Power Deliv. 2005, 20, 1806–1811. [Google Scholar] [CrossRef]
  3. Jouane, A.V.; Banerjee, B. Assessment of voltage unbalance. IEEE Trans. Power Deliv. 2001, 16, 782–790. [Google Scholar] [CrossRef]
  4. Semlyen, A.; Acha, E.; Arrillaga, J. Newton-Type algorithms for the harmonic phasor analysis of nonlinear power circuits in periodical steady state with special reference to magnetic nonlinearities. IEEE Trans. Power Deliv. 1998, 3, 1090–1098. [Google Scholar] [CrossRef]
  5. Kersting, W.H. Causes and effects of unbalanced voltages serving an induction motor. IEEE Trans. Ind. Appl. 2001, 37, 165–170. [Google Scholar] [CrossRef]
  6. Wannous, K.; Toman, P. Evaluation of Harmonics Impact on Digital Relays. Energies 2018, 11, 893. [Google Scholar] [CrossRef]
  7. Viswanadha, G.K.; Bijwe, P.R. Efficient reconfiguration of balanced and unbalanced distribution systems for loss minimization. IET Gener. Transm. Distrib. 2008, 2, 7–12. [Google Scholar] [CrossRef]
  8. Pillay, P.; Manyage, M. Loss of life in induction machines operating with unbalanced supplies. IEEE Trans. Energy Convers. 2006, 21, 813–822. [Google Scholar] [CrossRef] [Green Version]
  9. Angarita, M.L.; Ramos, G.A. Power calculations in nonlinear and unbalanced conditions according to IEEE Std. 1459–2010. In Proceedings of the Power Electronics and Power Quality Applications (PEPQA), Bogota, Colombia, 6–7 July 2013; pp. 1–7. [Google Scholar] [CrossRef]
  10. IEEE Standards Association. IEEE Standard Definitions for the Measurement of Electric Power Quantities under Sinusoidal, Nonsinusoidal, Balanced, or Unbalanced Conditions, IEEE Std. 1459–2010; IEEE Standards Association: New York, NY, USA, 19 March 2010; pp. 1–50. [Google Scholar]
  11. Willems, J.L. Reflections on apparent power and power factor in non-sinusoidal and polyphase situations. IEEE Trans. Power Deliv. 2004, 19, 835–840. [Google Scholar] [CrossRef]
  12. Emanuel, A.E. Apparent powers definitions for three-phase systems. IEEE Trans. Power Deliv. 1999, 14, 767–772. [Google Scholar] [CrossRef]
  13. Chicco, G.; Schlabbach, J.; Spertino, F. Characterisation and assessment of the harmonic emission of grid-connected photovoltaic systems. In Proceedings of the 2005 IEEE Russia Power Tech, St. Petersburg, Russia, 27–30 June 2005; pp. 1–7. [Google Scholar] [CrossRef]
  14. Tan, Y.T.; Kirschen, D.S.; Jenkins, N. A model of PV generation suitable for stability analysis. IEEE Trans. Energy Convers. 2004, 19, 748–755. [Google Scholar] [CrossRef]
  15. Orchi, T.F.; Mahmud, M.A.; Oo, A.M.T. Generalized dynamical modeling of multiple photovoltaic units in a grid-connected system for analyzing dynamic interactions. Energies 2018, 11, 296. [Google Scholar] [CrossRef]
  16. Jayatunga, U.; Perera, S.; Ciufo, P.; Agalgaonkar, A.P. Deterministic methodologies for the quantification of voltage unbalance propagation in radial and interconnected networks. IET Gener. Transm. Distrib. 2015, 9, 1069–1076. [Google Scholar] [CrossRef] [Green Version]
  17. Bucci, G.; Ciancetta, F.; Fiorucci, E.; Ometto, A. Survey about classical and innovative definitions of the power quantities under nonsinusoidal conditions. Int. J. Emerg. Electr. Power Syst. 2017, 18, 1–16. [Google Scholar] [CrossRef]
  18. Barkas, D.A.; Psomopoulos, C.S.; Ioannidis, G.C.; Kaminaris, S.D.; Malatestas, P. Experimental and theoretical investigation of harmonic distortion in high voltage 3-phase transformers. In Proceedings of the MedPower 2014, Athens, Greece, 2–5 November 2014; pp. 1–8. [Google Scholar] [CrossRef]
  19. Marcos, B.K.; Cursino, B.J.; Lima, A.M.N. Shaping control strategies for active power filters. IET Power Electron. 2018, 11, 175–181. [Google Scholar] [CrossRef]
  20. Limongi, L.R.; Bradaschia, F.; de Oliveira Lima, C.H.; Cavalcanti, M.C. Reactive Power and Current Harmonic Control Using a Dual Hybrid Power Filter for Unbalanced Non-Linear Loads. Energies 2018, 11, 1392. [Google Scholar] [CrossRef]
  21. Xu, J.; Gu, X.; Liang, C.; Bai, Z.; Kubis, A. Harmonic suppression analysis of a harmonic filtering distribution transformer with integrated inductors based on field–circuit coupling simulation. IET Gener. Transm. Distrib. 2018, 12, 615–623. [Google Scholar] [CrossRef]
  22. Abdel Aleem, S.H.E.; Ibrahim, A.M.; Zobaa, A.F. Harmonic assessment-based adjusted current total harmonic distortion. J. Eng. 2016, 2016, 64–72. [Google Scholar] [CrossRef] [Green Version]
  23. Badr, M.; Maarouf, M.; Basyouni, M.M.; Ahmed, S.A. Reducing harmonic distortion and correcting power factor in distribution systems. In Proceedings of the 22nd International Conference and Exhibition on Electricity Distribution (CIRED 2013), Stockholm, Sweden, 10–13 June 2013; pp. 1–4. [Google Scholar] [CrossRef]
  24. Carvajal, W.; Ordoñez, G.; Moreno, A.L.; Duarte, C.A. Simulation of electric systems with non-linear and time-variant loads. Ingeniare Rev. Chilena Ing. 2011, 19, 76–92. [Google Scholar] [CrossRef]
  25. Sezgin, E.; Göl, M.; Salor, Ö. State-estimation-based determination of harmonic current contributions of iron and steel plants supplied from PCC. IEEE Trans. Ind. Appl. 2016, 52, 2654–2663. [Google Scholar] [CrossRef]
  26. Zebardast, A.; Mokhtari, H. Technique for online tracking of a utility harmonic impedance using by synchronising the measured samples. IET Gener. Transm. Distrib. 2016, 10, 1240–1247. [Google Scholar] [CrossRef]
  27. Ahmed, E.E.; Xu, W. Assessment of harmonic distortion level considering the interaction between distributed three-phase harmonic sources and power grid. IET Gener. Transm. Distrib. 2007, 1, 506–515. [Google Scholar] [CrossRef]
  28. Diez, J.M.; Blasco, P.A.; Montoya, R. Formulation of phasor unbalance power: Application to sinusoidal power systems. IET Gener. Transm. Distrib. 2016, 10, 4178–4186. [Google Scholar] [CrossRef]
  29. Montoya-Mira, R.; Diez, J.M.; Blasco, P.A.; Montoya, R. Equivalent circuit and calculation of unbalanced power in three-wire three-phase linear networks. IET Gener. Transm. Distrib. 2018, 12, 1466–1473. [Google Scholar] [CrossRef]
  30. Buchholz, F. Die drehstrom-scheinleistung bei ungleichmassiger belastung der drei zweige. Licht und Kraft 1922, 2, 9–11. [Google Scholar]
Figure 1. Three-phase electrical system with non-linear load.
Figure 1. Three-phase electrical system with non-linear load.
Energies 11 01888 g001
Figure 2. Two parallel non-linear loads connected to a system bus.
Figure 2. Two parallel non-linear loads connected to a system bus.
Energies 11 01888 g002
Figure 3. Two system buses linked by a power line.
Figure 3. Two system buses linked by a power line.
Energies 11 01888 g003
Figure 4. Three-phase electrical system with non-linear loads.
Figure 4. Three-phase electrical system with non-linear loads.
Energies 11 01888 g004
Figure 5. Single-phase full-wave rectifier (loads 2 and 3).
Figure 5. Single-phase full-wave rectifier (loads 2 and 3).
Energies 11 01888 g005
Table 1. Load 2 and 3 values (see Figure 5).
Table 1. Load 2 and 3 values (see Figure 5).
ResistanceLoad 2Load 3
Ra (Ω)75
Rb (Ω)1215
Rc (Ω)2025
Table 2. Harmonic line-to-neutral voltage.
Table 2. Harmonic line-to-neutral voltage.
BusOrderVam (V)Vbm (V)Vcm (V)
mModulusAngleModulusAngleModulusAngle
Bus 11227.5−0.086223.0−120.1217.6119.9
30.261−137.40.188−127.10.14−121.1
50.196−137.00.143−11.270.099114.0
70.174−137.60.126108.20.085−5.423
Bus 21226.4−0.109222.2−120.1217.0119.9
30.704−131.90.521−122.40.376−115.9
50.555−131.50.411−5.4500.278120.2
70.507−133.00.375112.00.248−1.351
Table 3. Harmonic currents in each phase and positive-sequence harmonic currents.
Table 3. Harmonic currents in each phase and positive-sequence harmonic currents.
LoadOrderIan (A)Ibn (A)Icn (A)In+ (A)
nModulusAngleModulusAngleModulusAngleModulusAngle
L1145.59−3.6839.1−118.551.51121.445.36−0.271
30.032−143.20.013−124.80.01−98.120.008−144.7
50.036−137.20.025−8.2180.025103.40.005−144.9
70.033−138.80.208102.10.024−3.6520.088−136.7
L219.174−13.126.168−137.24.64798.816.652−16.25
32.8503.8951.92812.581.47220.610.4367.065
51.679−9.5471.113117.20.823−116.60.266−41.31
71.183−19.320.78−133.00.573113.00.842−14.59
L317.435−15.475.363−139.54.19796.695.656−18.68
32.2977.581.67815.971.33723.770.3128.613
51.342−7.510.954119.10.733−114.60.192−43.64
70.946−18.040.667−131.80.509114.30.704−13.12
Table 4. Harmonic line-to-neutral voltage in sequence components and voltages unbalance factors.
Table 4. Harmonic line-to-neutral voltage in sequence components and voltages unbalance factors.
BusOrderVm+ (V)Vm− (V)Vm0 (V)δmδm0
mModulusAngleModulusAngleModulusAngle
Bus 11222.7−0.102.87233.212.852−32.80.0130.013
30.035−133.20.042−177.50.195−130.21.1855.525
50.030−174.90.145−132.60.028−121.04.7780.910
70.128−133.00.026−121.30.028−176.30.2020.218
Bus 21221.87−0.102.71233.212.725−33.820.0120.012
30.095−124.40.113−174.10.53−125.11.1925.607
50.088−171.10.413−126.90.079−114.64.6920.894
70.376−128.80.075−114.50.08−173.00.2010.214
Table 5. Results in Bus 1.
Table 5. Results in Bus 1.
Harmonic A n m (VA) B n m (VA) C n m (VA) D n m (VA) P n + m + (W) Q n + m + (VAr) F u m S n m (VA)
mn
11−2115.82−2638.24−138.238−1463.14138,155.52534.160.01838,423.27
13−2057.60−583.3531460.40188.710490.788−64.7490.0182643.823
15−1240.00−70.671−820.267−79.911223.934207.2230.0181521.736
17−261.988−52.893−84.365−128.994960.705287.2990.0181049.309
310.297−0.0800.4230.263−3.843−4.7025.65035.010
33−0.123−0.2370.2330.211−0.061−0.0505.6502.409
350.094−0.1320.147−0.0910.000−0.0485.6501.387
370.003−0.0140.043−0.018−0.071−0.1435.6500.956
510.3300.3060.0270.224−5.158−0.8174.86426.051
530.2890.009−0.2020.037−0.0680.0034.8641.793
550.168−0.0110.1140.006−0.028−0.0314.8641.032
570.032−0.0120.0160.018−0.127−0.0514.8640.711
711.072−0.2911.5350.947−13.864−17.0450.29723.027
73−0.449−0.8540.8490.757−0.219−0.1810.2971.584
750.339−0.4770.530−0.3300.000−0.1750.2970.912
770.0110.0490.156−0.067−0.255−0.5170.2970.629
Table 6. Results in Bus 2.
Table 6. Results in Bus 2.
Harmonic A n m (VA) B n m (VA) C n m (VA) D n m (VA) P n + m + (W) Q n + m + (VAr) F u m S n m (VA)
mn
11−2052.78742.814−72.357−181.1557821.8422430.9560.0178480.335
13−2058.52−587.6481456.722195.425493.339−67.6990.0172644.456
15−1246.02−61.539−825.547−69.403225.739204.5710.0171528.464
17−242.657−91.012−51.341−146.629999.69245.9240.0171073.082
31−0.299−0.3430.611−0.54−1.02−3.3375.73221.022
33−0.552−0.3870.8350.27−0.142−0.1575.7326.555
350.184−0.3670.359−0.3310.018−0.1295.7323.789
37−0.006−0.0050.096−0.085−0.153−0.4115.7322.66
510.666−0.5180.1320.166−2.916−1.4394.77716.426
530.822−0.12−0.5620.236−0.198−0.0044.7775.122
550.484−0.0780.3360.001−0.076−0.0944.7772.961
570.072−0.0570.0470.061−0.377−0.1584.7772.078
71−1.096−1.7352.602−1.667−5.063−12.9050.29314.954
73−1.804−2.0922.9591.712−0.611−0.580.2934.663
750.86−1.4491.505−1.1490.031−0.5150.2932.695
77−0.042−0.0550.426−0.274−0.733−1.5810.2931.892
Table 7. Results in load 1.
Table 7. Results in load 1.
Harmonic A n m (VA) B n m (VA) C n m (VA) D n m (VA) P n + m + (W) Q n + m + (VAr) F u m S n m (VA)
mn
11−55.554−3384.39−65.46−1281.0530,304.5993.0230.01830,525.2
138.8375.75−1.922−6.779−4.4133.1370.01813.791
1510.703−9.1268.391−10.302−2.6251.8530.01819.617
17−18.46438.519−32.77918.204−42.70740.3170.01881.743
310.3830.1170.1710.375−3.275−3.5275.6527.813
3300.0020−0.0010.00105.650.013
3500.00200.002005.650.018
370.008−0.00500.0010.0090.0015.650.074
510.0750.455−0.0010.178−4.12−0.3884.86420.696
53−0.001−0.00100.0010.00104.8640.009
55−0.0010.001−0.0010.001004.8640.013
570.0040.0050.004−0.0030.006−0.0054.8640.055
711.3860.4180.6241.357−11.819−12.7860.29718.293
73−0.0010.006−0.001−0.0030.0030.0010.2970.008
750.0010.00800.0070.00200.2970.012
770.028−0.0170.0010.0030.0340.0020.2970.049
Table 8. Results in load 2.
Table 8. Results in load 2.
Harmonic A n m (VA) B n m (VA) C n m (VA) D n m (VA) P n + m + (W) Q n + m + (VAr) F u m S n m (VA)
mn
11−1200.93419.477−38.06−101.7124253.0731231.7140.0174608.92
13−1133.33−294.699779.61870.706288.081−36.2320.0171438.41
15−685.076−42.472−440.704−48.682133.342116.7640.017836.27
17−142.086−51.436−28.3−84.158542.853140.2980.017587.532
31−0.174−0.2030.359−0.31−0.586−1.7935.73211.425
33−0.298−0.2020.4620.131−0.082−0.0935.7323.566
350.1−0.1980.202−0.1770.009−0.0755.7322.073
37−0.004−0.0030.056−0.049−0.081−0.2255.7321.456
510.392−0.2980.0740.096−1.591−0.7474.7778.927
530.448−0.075−0.2950.141−0.115−0.0044.7772.786
550.266−0.040.1810.006−0.045−0.0544.7771.62
570.042−0.0330.0270.035−0.204−0.0894.7771.138
71−0.632−1.0231.523−0.95−2.872−6.9220.2938.127
73−0.983−1.1041.6530.874−0.353−0.3420.2932.536
750.468−0.7820.8438−0.60930.0132−0.29970.2931.47
77−0.025−0.0340.2487−0.156−0.3889−0.86560.2931.04
Table 9. Results in load 3.
Table 9. Results in load 3.
Harmonic A n m (VA) B n m (VA) C n m (VA) D n m (VA) P n + m + (W) Q n + m + (VAr) F u m S n m (VA)
mn
11−851.848323.337−34.297−79.4443568.7681199.24120.0173875.11
13−925.192−292.949677.104124.719205.258−31.4670.0171208.05
15−560.943−19.067−384.842−20.72092.396587.80680.017692.781
17−100.571−39.576−23.041−62.471456.837105.6250.017485.841
31−0.125−0.1400.252−0.230−0.4337−1.54435.7329.606
33−0.254−0.1860.3730.139−0.0603−0.06475.7322.995
350.084−0.1700.157−0.1540.0088−0.05365.7321.717
37−0.002−0.0010.040−0.037−0.0724−0.18625.7321.204
510.274−0.2200.0580.070−1.3247−0.69174.7777.506
530.375−0.045−0.2670.096−0.0824−0.00044.7772.340
550.218−0.0380.156−0.005−0.0308−0.04024.7771.342
570.029−0.0240.0200.026−0.1726−0.074.7770.941
71−0.095−0.0670.161−0.233−0.2011−1.24140.2936.531
73−0.230−0.0620.3180.015−0.0413−0.0560.2932.036
750.050−0.1330.109−0.1400.0116−0.04100.2931.168
770.0020.0020.024−0.035−0.0399−0.15150.2930.819
Table 10. Results in line 1–2 (Bus 1–Bus 2).
Table 10. Results in line 1–2 (Bus 1–Bus 2).
Harmonic A n m (VA) B n m (VA) C n m (VA) D n m (VA) P n + m + (W) Q n + m + (VAr) F u m S n m (VA)
mn
11−7.4673.324−0.44−0.92829.05110.1850.25232.861
13−7.882−1.455.5690.0581.862−0.1880.25210.247
15−4.683−0.007−3.115−0.2060.820.7990.2525.923
17−0.865−0.398−0.246−0.5673.7221.0580.2524.158
310.1880.15−0.3340.4240.4522.1625.73913.318
330.3930.128−0.564−0.0370.0810.1075.7394.153
35−0.0950.232−0.2080.239−0.0190.085.7392.4
37−0.002−0.002−0.050.0640.0730.2685.7391.685
51−0.4080.366−0.107−0.1171.8781.0094.74210.699
53−0.5260.1280.353−0.1990.1290.0074.7423.336
55−0.3130.065−0.2230.0040.0480.0634.7421.928
57−0.0430.04−0.035−0.0390.2430.1124.7421.354
710.7621.025−1.6711.2593.0188.6470.2919.875
731.331.209−2.082−0.9290.3890.3980.2913.079
75−0.5250.962−0.9710.813−0.0340.3390.2911.78
770.0220.024−0.2680.2030.4441.0620.2911.25
Table 11. Summary of apparent power components.
Table 11. Summary of apparent power components.
ID S 1 + 1 +   ( VA ) S T 1 1   ( VA ) S T H   ( VA ) S T   ( VA )
Bus 138,239.5538,423.2683328.83138,567.196
Bus 28190.8948480.3353339.879114.319
Load 130,304.73730,525.19195.77530,525.341
Load 24427.8394608.9181820.7844955.54
Load 33764.8753875.1091521.0954162.956
Line 1–230.78532.86126.8442.429

Share and Cite

MDPI and ACS Style

Blasco, P.A.; Montoya-Mira, R.; Diez, J.M.; Montoya, R.; Reig, M.J. Formulation of the Phasors of Apparent Harmonic Power: Application to Non-Sinusoidal Three-Phase Power Systems. Energies 2018, 11, 1888. https://doi.org/10.3390/en11071888

AMA Style

Blasco PA, Montoya-Mira R, Diez JM, Montoya R, Reig MJ. Formulation of the Phasors of Apparent Harmonic Power: Application to Non-Sinusoidal Three-Phase Power Systems. Energies. 2018; 11(7):1888. https://doi.org/10.3390/en11071888

Chicago/Turabian Style

Blasco, Pedro A., Rafael Montoya-Mira, José M. Diez, Rafael Montoya, and Miguel J. Reig. 2018. "Formulation of the Phasors of Apparent Harmonic Power: Application to Non-Sinusoidal Three-Phase Power Systems" Energies 11, no. 7: 1888. https://doi.org/10.3390/en11071888

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