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:
where
m is the harmonic order of the voltage and n is the harmonic order of the current.
and 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.
and are the root mean square (RMS) values for the harmonic voltage and harmonic current of orders m and n, respectively.
and 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:
where the first term,
, is the instantaneous active power, which is decomposed into three instantaneous components:
Here,
is the fundamental instantaneous active power, which is calculated as
This expression has two terms: the fundamental active power , and the oscillating component . The first term is equal to the average value of 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),
, is the harmonic instantaneous active power of order
h for
h ≠ 1, and is calculated as
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 . The sum is equal to the average value of and the sum of and 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 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),
, is the harmonic instantaneous active power, which is caused by the harmonic voltage of order
m and the harmonic current of order
n, where
m ≠
n. This term is calculated as
The average value of is zero.
The second term in (3),
, is the instantaneous reactive power, which is also decomposed into three instantaneous components:
where
is the fundamental instantaneous reactive power, which is calculated from (9).
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.
is harmonic instantaneous reactive power, which is caused by the harmonic voltage of order
m and the harmonic current of order
n, where
m ≠
n. This power is calculated from (11).
In a non-sinusoidal three-phase system, the instantaneous power and its components are expressed by the following equation:
For each of the phases (a, b and c), the instantaneous active power and the instantaneous reactive power 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
where
,
and
are line-to-neutral harmonic voltages of order
m and
,
and
are harmonic currents of order
n.
These voltages can be expressed in terms of the following symmetrical components:
Similarly, the symmetric current components are
where
and are the RMS values for the positive-sequence harmonic voltage and current of orders m and n, respectively.
and are the RMS values for the negative-sequence harmonic voltage and current of orders m and n, respectively.
and are the RMS values for the zero-sequence harmonic voltage and current of orders m and n, respectively.
, and are, respectively, the angles of positive-, negative- and zero-sequence of the harmonic voltage of order m.
, and 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 (
and
) become zero. The positive-sequence voltages for each phase (
a,
b and
c) are determined as follows:
where
.
Similarly, the phase-specific positive-, negative- and zero-sequence currents are
Under these conditions, the instantaneous power of the system is given by
where
Here,
is the instantaneous power caused by positive-sequence voltages and currents.
is the instantaneous power caused by positive-sequence voltages and negative-sequence currents.
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
where
is the angle between the positive-sequence harmonic voltage of order m and the positive-sequence harmonic current of order n.
is the angle between the positive-sequence harmonic voltage of order m and the negative-sequence harmonic current of order n.
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:
where
and
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 , , and
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:
Substituting these equalities into (24) and (25) produces the following expression:
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:
which can be decomposed into
,
,
and
as follows:
where
These instantaneous parameters are sinusoidal waveforms with zero average value and RMS values given by:
3.2. Phasor of Harmonic Apparent Power
For each harmonic voltage of order
m ≠ 1 and harmonic current of order
n ≠ 1, the phasor of harmonic apparent power
can be defined as:
where
, , and are the unit vectors associated with parameters , , and , respectively (see (35)–(38)).
and are the unit vectors associated with the positive-sequence active power and the positive-sequence reactive power (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:
For
m =
n = 1, the phasor of fundamental apparent power is obtained as:
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 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:
3.4. Application of between Two Buses in a System
Figure 3 shows two buses in a system (
i and
j) linked by a power line with impedance
.
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’:
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:
Substituting the voltage unbalance factors of order
m,
and
, into (44) produces the following expression:
where
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
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
as follows:
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 11 compares the following total apparent power values:
(positive-sequence fundamental apparent power), which is calculated from
and
;
. (fundamental apparent power), which is calculated from (45);
(total harmonic apparent power), which is calculated from (47); and
(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.