In the Introduction the negative (or detrimental) effects of supraharmonics have been briefly reviewed and may be synthesized as follows:
Such phenomena are now analyzed more in detail in the following.
2.1. Power losses
The basic rationale for power losses is related to resistive components in conductors, that increase in value roughly proportional to the square root of frequency, as a consequence of the skin effect. If conducting materials are ferromagnetic then the phenomenon is more evident (being the skin depth inversely proportional to the magnetic permeability). In addition, losses may occur also due to the magnetic properties of the material and hysteresis phenomena, as commonplace for transformers. When windings (and tightly packed conductors) are considered, a second effect affects the distribution of the current in the conductors cross section, the proximity effect, also approximately proportional to frequency.
An increase in conductor losses due to frequency is known and stigmatized in the [
25] (Figure 14A.4) when testing Instrument Transformers (ITs) for thermal losses at AC (50 or 60 Hz) rather than at DC: for large conductors with diameter
D ranging between 40 mm and 100 mm additional losses with respect to an equivalent DC current in terms of rms are about 150% at 50 mm diameter and 270% at 100 mm diameter for Alu conductors, and 180% and 330% for Copper conductors. The estimate is done using the Levasseur formula:
with
D the conductor diameter and
the skin depth in the conductor material.
Figure 12 - ref.
A more accurate evaluation of SH effect up to 100 kHz at frequencies higher than those characteristics of the mains and its harmonics [
26]. The work takes into account the formulas of the IEC 60287-1-1 standard [
27] providing separately the correction factors
for skin effect,
for proximity effect and
for the screen losses (so well applicable to MV cables).
with
containing the frequency dependence.
The IEC 60287-1-1 proposes then two formulations of
for the case of two-core cables (or two single-core cables) and three-core cables (or three single-core cables). The two cases are indicated by a subscript “2” and “3”.
with
containing the frequency dependence,
is the conductor diameter and
s is the separation of conductor axes (so equal to the sum of the diameter and the gap).
The two “Y” quantities are based on the calculation of the corresponding “X” factors, that in turn contain two “k” coefficients, related to the shape of the conductors: is largely equal to unity for round solid and stranded conductors for both copper and aluminum; is also unity for solid conductors and takes the value 0.8 for stranded conductors.
Relevant to note, this formula is stated as valid for
, and, as observed in [
26], this does not fit many practical cases, especially for SH frequencies.
A similar, but different, “
” quantity is thus proposed [
26], taken from [
28]:
where the two factors
and
are again expressed with different approximations at variable
, as done before for
, and
m is a factor for the cable geometry.
For the last quantity, , that accounts for the effect of the cable screen, the IEC 60287-1-1 proposes a simplified approach, taking the average of the two other factors, and .
To conclude, the three “Y” factors are put together to determine the “new” resistance
including dependence on frequency and geometry:
If the simplified expression for
is considered, then (
8) can be rewritten as:
For the accuracy of the published methods to predict additional losses, Topolski et al. [
29] provided a comprehensive overview and comparison of various formulas, unfortunately with the frequency interval limited to about 2 kHz. The reference formulation is the one achieved by means of Bessel functions and the selected formulations for the comparison are the Levasseur equation (
1) (that the Authors refer to a Schneider Electric application note) and the IEC 60287-1-1 formulation considered above.
The results for a single-strand and multi-strand cable show that the IEC 60287-1-1 formulation for the two skin and proximity effect terms is not so accurate, and underestimates losses for frequencies above about 1 kHz. However, when the terms of skin and proximity effect are taken together (so the sum
is evaluated), the correspondence with the theoretical reference method based on Bessel functions is almost perfect (up to
kHz. In [
26] the comparison was carried out against Finite Element Method (FEM) results, showing a very good correspondence (maximum error of about 6%), if also
is accounted for.
The relevance of supraharmonic components evaluated in [
26] over the 2 to 100 kHz interval is confirmed with an increase of total resistance by a factor of 6.5-7, proportional to the square root of the frequency ratio (=100/2=50). The minimum ratio with respect to the fundamental is about 3.5.
Looking for a suitable weight of SH amplitudes with the objective of limiting the losses to the same amount as caused by harmonic components, we observe a factor inversely proportional to the square root of frequency giving the already mentioned factor of 7 between 2 kHz and 100 kHz. Before proceeding the different standpoint of the IEEE Std. C57.110 [
30] that defines a quadratic function of the harmonic order as impacting on eddy current losses of transformers, whereas other stray losses have an exponent of 0.8 or less on the harmonic order (that agrees somewhat with the “square root of frequency” behavior). The standard, in its Annex C, reports some considerations on the effect of distortion at higher frequency (such as in the SH interval), observing that skin effect becomes relevant and deviations from the postulated “
” behavior are significant. A reduction of 10 % is then estimated already for the first 20 harmonics if skin effect is taken into account. This correction is, however, not sufficient in the SH Interval to match the
behavior.
In general, the total loss increase due to supraharmonics for the observed amplitude values is in the order of 10%-20% of the losses ascribed to harmonics (up to 2 kHz). The present harmonic limits have the objective of limiting losses to a value that requires an over-rating of cables of about 5% [
16], considering such amount of harmonic losses acceptable. This ratio is bounded to increase with the extensive use of static power conversion units with increasing switching frequency (renewables, EV chargers, etc.) shifting distortion to higher frequency, namely reducing the harmonic contribution in favor of increased SH components.
Proposing a SH limit curve originating from a weight starting from the harmonic limits for residential and industrial environments would imply additional losses equal to those allowed by present limits in the harmonic interval. Going down to a limit value for the single SH component is then more difficult if the number of SH components is not known a priori. A conservative standpoint may be assumed, by adopting a sufficiently large resolution bandwidth (e.g. the frequently adopted 2 kHz), dividing then the SH interval into 74 sub-intervals evenly distributed between 2 and 150 kHz with steps of 2 kHz. The risk is of course an excessive penalization, considering risky and non compliant situations that see an uneven distribution of components, as it is often the case.
2.4. Damage to MV cable terminations
Aging and failure of MV cable terminations is caused by localized heat possibly combined with worsened PD occurrence. Localized heat in particular has been observed in resistive stress-grading elements, where SHs affect the electric field distribution causing hot spots [
17]: the large intensity of the flowing current is due to capacitive current proportionally higher for increasing frequency.
A catastrophic failure (that represents a field evidence) occurred at Eagle Pass substation [
32] and in successive laboratory tests at ABB, Sweden, also reported in the same reference. The measured on-site voltage distortion was distributed among some characteristic SH components with significantly large amplitude, when compared to the
kV nominal value. The SH components were:
kHz with
kV amplitude,
kHz (the 3rd harmonic of the switching fundamental) with
kV amplitude and an atypical value of
kHz with
kV amplitude, due to network resonance. Relevant voltage distortion components around the
kHz resonance are in reality more than one, with the two adjacent to the resonance peak showing amplitude of 2 kV and
kV, leading to an overall rms value of
kV.
During lab tests it was determined that using the worst-case observed SH voltage levels, the heat caused at kHz was 20 times higher than at kHz. This is reasonable as the electric field was approximately 3 times larger (due to a different spatial distribution estimated by finite element method simulations) and the capacitive current approximately 10 times larger (following the frequency ratio).
The mitigation adopted was a different type of cable termination and addition of a filter to reduce the kHz component caused by the resonance.
An approach is proposed in [
17] starting from the determination of the heat
Q dissipated in the stress grading area and in the dielectric of the cable joint. This quantity
Q is proportional to the frequency and to the square of the voltage, as it is put in relationship to the electric field intensity. Any temperature increase in the cable joint is shown to be proportional to the excess heat caused by the supraharmonics. For this reason the
quantity is calculated separating and weighting the SH excess heat with that at the nominal frequency. The weighted
of SHs, each with amplitude
multiplied by frequency, is divided by the weighted
of the fundamental (i.e.
multiplied by the fundamental frequency
). It may be said that the weighting with frequency accounts for the capacitive nature of the current, increasing with frequency. For a single SH component of frequency
and amplitude
the resulting
Q value is:
When
approaches the value of 20, failure is very likely, based on the site and lab evidence of [
32], but subject to a significant uncertainty as the Eagle Pass cable was a unique case and the cable might have been defective, or particularly susceptible to thermal stress. Some safety margins are then proposed [
17] consisting of a multiplying factor
m, that takes values of 0.25, 0.5 or 1, identifying a green (
), yellow (
), orange (
) and red (
) area (the equal signs have been omitted and assigning them to the lower or upper interval does not represent an issue).
A further simplification shown in [
17] is to consider supraharmonics below and above 20 kHz, providing two straight limits rather than a more complex curve. Maybe this approach is oversimplifying since the capacitive current is proportional to frequency and between 2 and 20 kHz there is a decade of frequency, as well as 20 kHz and 150 kHz differ by a factor 7.5. In addition, the distribution of the electric field at the termination was quite variable depending on frequency as well. For both factors the variability is much larger than the trimming operated by
m, so that a more accurate frequency discrimination seems necessary.
A weighted sum, as commonly done for the weighted total harmonic distortion index, could be adopted to provide an equivalent representation of the heat stress caused by a complex pattern of SH components. In the end what is relevant is the heat stress, that may be assumed proportional to the frequency. Since the SH spectrum may be characterized by a mix of narrowband and broadband components, aggregation with a large resolution bandwidth (e.g. 2 kHz) is more robust and the result less exposed to changes in the lateral bands of switching groups and other modulation byproducts. The examples in [
17] indicate a risk of cable termination failure for single SH components in the order of 3 % to 8 %, that means that for a SH spectrum populated by
n components of not-too-different frequency the limit must be reduced by
and with different frequencies for the mentioned components, each must receive a limit that is proportional to
.
2.5. Interference to equipment
This and the next subsection distinguish the two cases of interference to equipment in general and to PLC systems, in particular.
In general, due to the so called “secondary emissions”, amplification of normal emissions is expected, as caused by the presence of other loads connected to the same grid: the reason is that the feeding impedance may be altered by the presence of e.g. input filters and capacitive components, with consequential amplification of current emissions at some frequencies in the SH interval.
This was tested experimentally by combining in pairs an electric vehicle, a LED lamp and a TV set in [
33]. The amount of measured emissions from the EV were in the order of 5mA–10mA below 10 kHz and 5 mA located at an isolated component at 55 kHz. Theoretically this was demonstrated in [
20,
34] considering a grid feeding a set of identical EVs located at different distances from the PCC on the same or two different branches of the LV grid. The EV was modeled with an experimentally determined input impedance, and the mutual effect between EVs caused amplification of current emissions as soon as they were connected to not-to-far feeding points (in the order of 20m–50m.
For comparison the EN 50627 [
35] reports various cases of measured emissions from a wide range of sources, including portable tools, electrical appliances, and inverters (such as active infeed inverters, AICs). Low-power devices (such as LED lamps, tools and various power supplies feeding e.g. a fibre switch, a computer, a surveillance camera to cite the most prominent examples) are in any case able to cause emissions in the order of 90dB
V–100dB
V in the frequency interval 9kHz–70kHz when measured with a standardized Line Impedance Stabilization Network (LISN): the equivalent current emission level may be obtained considering a LISN impedance variable between 5
and 15
, leading to a range of amplitudes of about 5mA–20mA for the largest emissions. This is in agreement with [
33] and provides an indication of worst-case emissions from low-power devices.
Inverters interfacing renewables and energy storage devices are responsible in general for higher emissions, located as narrowband components at the switching frequency harmonics, usually with a fundamental located in the 10kHz–20kHz range (the examples of the EN 50627 [
35] report fundamentals of 16 and 18 kHz and the three PV inverters of Figure 20 of the EN 50669 [
36] have fundamental of 15, 16 and 20 kHz). Except for dramatic cases, as reported for Belgium in the EN 50627, usual emission levels are about 20 dB higher than the previous ones, thus implying as well ten times higher current emissions and significant chances of interference.
This type of emissions has been observed in normal domestic LV grids, as electrical appliances, EVs and PV inverters are commonplace in the residential grids. The maximum values of emissions were in the order of 0.3V–1V up to 20 kHz, decreasing then gradually by no more than an order of magnitude.
A complex scheme of interference has been observed also for residual current devices (RCDs) [
37] with different responses depending on models. Besides undesirable tripping by superposition of SH components, what is most relevant is the possible RCD desensitization tripping then at higher fundamental (50 or 60 Hz) current values than required for electrical safety. Tests performed in [
37] at significant, but not unlikely, injected SH current components ( 800 mA between 1 and 5 kHz) confirmed a desensitization of 2–2.5 for two different RCDs, tripping then at residual current values of 58 mA and 76 mA at 50 Hz with significant implications in terms of electrical safety. Such current values under the usual assumption of a 5
grid impedance (the LISN value), correspond to 4 V of SH voltage distortion. By extrapolation of the two SH disturbance levels ( 800 mA and 1300 mA) for the tested RCD model (“A-30”), the minimum SH current causing a deviation above the desired 30 mA tripping threshold is about 600 mA, that corresponds to about 3 V or 130 dB
V at low frequency (approximately over 2kHz–10kHz). About the assumed 5
value of grid impedance, it represents for sure an overestimation if compared for example with the impedance curves shown in Figure 18 of the IEC 62578 [
31], lying in the 0.5
–4
interval. Such an overestimation provides a margin when calculating voltage distortion from an observed current, but not the opposite, and we have seen that current, rather than voltage, disturbance is often relevant.
The EN 61000-4-19 [
38] has the objective of establishing relevant and reliable immunity test levels for the differential-mode disturbance in the SH interval. Unfortunately, this standard has not been applied yet, nor enforced by inclusion in any product or generic EMC standard. The first three test profiles (Level 1, Level 2 and Level 3) have amplitude of
, 3 and 12 V set up to 9 kHz and the decreasing linearly in a log scale by a factor of 5 reaching values of
,
and
V at 95 kHz, where they are then held constant up to 150 kHz.
A straightforward comparison to the emission levels of PLC devices will be done at the end of this section. It is clear that in light of aggregation from different sources, as well as some margin of variability due to e.g. grid resonances and slightly different products characteristics, Level 1 does not provide any confidence that a tested equipment will be immune to SH disturbance. Level 2 provides a safety margin of 10dB–15dB with respect to the observed values, barely sufficient to cover the said aggregation. In fact, 10 sources disturbing the same frequency band within about 50 m would aggregate with a factor of 3.2 minimum, as resulting from a
assumption. This was demonstrated for harmonic and SH emissions of EVs in [
39]. Level 3 thus seems highly advisable for all grids, noting that there is no big distinction any longer between residential, light industrial and industrial environments for what regards disturbance in the SH interval.
Some cases of highly disturbing inverters should instead be mitigated at the source, providing the necessary decoupling, e.g. by means of series reactors.
2.6. Interference to PLC
Besides some episodes of interference caused by PLC systems (e.g. to lighting, as cited in sec. 5.2.1 of the EN 50627), PLC technology is in general the victim of various forms of emissions. Emissions in the SH frequency interval are not in general disciplined by applicable limits and even for those products with an emission standard (e.g. lighting and electrical appliances with the EN 55014-1 [
40]) there are episodes of non complying devices. It is observed that the limits of the EN 55014-1 in the first 9kHz–50kHz range are quite high, allowing 110 dB
V of quasi-peak amplitude.
Emissions measurements in the SH interval carried out before 2010 seemed to lead to the conclusion that interference to PLC was unlikely, as observed in [
41]. However, for the successive years various episodes of interference to PLC devices have been described in the EN 50627 [
35], mainly based on documented investigations carried out in Sweden by the local authority. The interfering levels at various frequencies are confirmed of being in excess of 100 dB
V and interference was suppressed with the application of mitigations (namely EMC filters), bringing emissions to a safe lower level. These two concomitant amplitude values for selected components (with and without interference) for the same system under the same conditions are extremely useful to define a reference line of interference below the lowest of such interfering cases.
A study of the various apparatuses connected to a microgrid (PV inverters, turbine, pump and battery charger) [
42] has shown that the PRIME PLC system (an Orthogonal Frequency-Division Multiplexing, OFDM, Narrow-Band device) is significantly disturbed by narrowband components originating from the battery charger and falling into the PLC operating band, between approximately 40 and 80 kHz: the relevant disturbing components are 104 dB
V at 48 kHz and
dB
V at 72 kHz. It is curious that two narrowband components can disturb a system that in principle should be able to hop on different frequencies and tolerate momentary interference on some of them (if not to exclude such frequencies from the list of the OFDM ones).
The various emissions cases discussed so far are collectively represented in
Figure 1, showing also the cases of levels with confirmed interference to PLC devices discussed above.
2.7. SH transfer efficiency between MV and LV levels
The negative effects discussed above can be seen as peculiar of the MV or LV level when the involved physical phenomena and devices are considered with a closer look: interference to equipment and to PLC systems is likely to occur at the LV level, whereas cable termination were said characteristic of MV grids, as well as aging of insulating materials especially at a pre-existing high fundamental voltage stress. Conversely dielectric aging may be seen equally occurring at LV and MV levels, as filters are used indistinctly in both cases.
What turns out to be relevant to liaise commented SH distortion levels and to derive comprehensive acceptance limits is the quantification of the transfer between the MV and the LV sides of distribution transformers at the SH frequencies. The SH frequency interval is such that the relevant coupling between the primary and secondary windings is a mix of inductive terms (designed magnetic coupling between windings operating at the fundamental frequency and above it) and capacitive terms (mutual stray capacitance).
A few publications [
43,
44] provide some measured curves of the transfer ratio between the LV and MV sides and between phases. Usually such transfer ratio is expressed by a multiplication coefficient of the nominal transfer ratio at the fundamental: it is thus expected around unity at low frequency and then decreasing at high frequency, except in the case of resonances.
The behavior of a small 100 kVA, 20/
kV (50:1 nominal ratio), delta-wye connected, transformer was studied in [
43], by measuring the transfer ratios from LV to MV first (test named “P1”, transfer ratios
,
,
), with
k indicating LV to MV and “11” used for same winding at LV and MV respectively, “12” for
finire finire) and then applying the test signal on the MV side, to measure the opposite MV-to-LV transfer capability (test “P2”, transfer ratios
,
,
), with
h indicating MV to LV and the numeric indexes having the same meaning). What may be concluded is:
the transfer between LV and MV sides undergoes a significant resonance that is not visible from the MV side; apart from this (where the transfer ratio peaks to almost 10x the 50:1 nominal ratio), the values are low between about 0.2 and 0.01;
the usual transfer behavior is that the phase L1 on the LV side influences the corresponding phase L1 on the MV side and similarly the phase L2, but not phase L3; this for a delta-wye transformer as customary used for MV to LV distribution;
the transformer has a symmetric behavior for which the self transfer ratios (each LV phase to the corresponding one on the MV side) is the same (Figure 8);
the transfer from MV to LV for the same phase is more effective and does not show significant variations vs. frequency, being at around unity (ranging between 0.64 and 2.25 up to 80 kHz).
The more recent study in [
44] has tested a 1 MVA, 20/
kV, delta-wye connected transformer for both current and voltage transfer characteristics. The frequency interval is unfortunately limited to below 8 kHz, thus limiting the generality of the results that occur around the transformer resonance at some kHz (as identified in [
43]) with the risk of being pessimistic. The reported results may be synthesized as:
regarding the current LV-to-MV transfer the results for L1 to L1 show an almost unity transfer ratio between kHz and kHz with a slight amplification (30%) at some components;
voltage, instead (always for L1 to L1 from LV to MV) is attenuated by a factor of 3 to 10 in the same frequency range;
for the MV to LV transfer the results provided by [
43] are not fully confirmed, having found a slightly larger variation (more persistently around a factor of 2 to 3); what is relevant is that in case of an unloaded transformer (not magnetized) the behavior is quite different and variable;
last, also the LV-to-LV transfer occurring between the secondary windings of two different transformers through the MV grid was studied and the observed transfer ratio is more than unity (e.g. 2 to 3) at several frequency points, whereas some attenuation should be in general expected; this is a relevant result for what regards the propagation of interference within the same LV grid but on different feeders and parts of the grid.