Indices to Study the Electrical Power Signals in Active and Passive Distribution Lines: a Combined Analysis with Empirical Mode Decomposition

The broad diffusion of renewable energy-based technologies has introduced several open issues in the design and operation of smart grids (SGs) when distributed generators (DGs) inject a large amount of power into the grid. In this paper, a theoretical investigation on active and reactive power data is performed for one active line characterized by several photovoltaic (PV) plants with a great amount of injectable power and two passive lines, one of them having a small peak power PV plant and the other one having no PV power. The frequencies calculated via the empirical mode decomposition (EMD) method based on the Hilbert-Huang transform (HHT) are compared to the ones obtained via the fast Fourier transform (FFT) and the wavelet transform (WT), showing a wider spectrum of significant modes mainly due to the non-periodical behavior of the power signals. The results obtained according to the HHT-EMD analysis are corroborated by the calculation of three new indices that are computed starting from the electrical signal itself and not from the Hilbert spectrum. These indices give the quantitative deviation from the periodicity and the coherence degree of the power signals, which typically deviate from the stationary regime and have a nonlinear behavior in terms of amplitude and phase. This information allows to extract intrinsic features of power lines belonging to SGs and this is useful for their optimal operation and planning.


Introduction
The massive use of non-linear loads and electronic-based equipment in residential and industrial plants greatly affects the operation and management of distribution lines.These problems are rapidly increasing after the deregulation of the electrical power energy sector.Therefore, the real-time and off-line monitoring of electrical power systems has become a very important task.Several techniques have been used to extract relevant features from measured quantities in order to monitor, diagnose and model the power systems [1], especially for the distribution system.Smart distribution technologies provide new capabilities to automatically locate and isolate faults using automated feeder switches and re-closers, dynamically optimize voltage and reactive power levels, and monitor asset health to effectively guide the maintenance and replacement of equipment [2].As a result, the electrical grid becomes more flexible, reliable, reconfigurable, efficient, secure, sustainable, intelligent and observable at a wide range of voltage levels.Therefore, SGs are investigated and deployed to change the way electrical power systems are planned, designed, managed, and monitored [3,4].Nevertheless, the SGs require massive information and the ability to manage big data, i.e., "high volume, high velocity, and/or high variety information assets that require new forms of processing to enable enhanced decision making, insight discovery and process optimization" [5].Although a considerable amount of effort was spent to compact the electrical information [6][7][8][9][10][11], an analysis of the current literature reveals that the compression of signal waveforms from power systems is far from being as mature as for speech, image, and video compression [4].
The several electrical signals of an SG can be, at the same time, linear and nonlinear, stationary and non-stationary, periodic and non-periodic, because they belong to lines with very different loads (static, rotating, power electronics-based, and so on) and generators (deterministic-based, stochastic-based).For example, for SGs with high penetration of DGs, a significant amount of conventional generation can be replaced by distributed PV resources with the result of the lack of the reactive power [12].The reduced system inertia is another consequence of utilizing higher amounts of PV generators [13,14], which declines the overall power system stability [15], such as during transient periods.The power disturbances can be measured in several modes depending on line parameters (such as power fluxes, length of the line, low/medium/high voltage network, passive or active line, etc.).When unexpected fluctuations appear, they introduce anomalies in the correct operation of the system.For example, while a short circuit manifests itself as a high-frequency component [16], a load variation gives rise to a low-frequency component [17].These problems are even more important for active lines, i.e., lines that can either absorb or feed the active power, than for passive lines that can only absorb active power.
Roughly speaking, the mathematical tools for the electrical signal analysis by means of compression techniques can be grouped in two categories: lossless techniques [18][19][20][21] and lossy algorithms [6,10,[22][23][24][25].Table I of [4] reports a comparison among them, highlighting the compression ratio and the distortion values.Nevertheless, it is not easy, even for an expert, to decide what is the best method to use in any specific case.Sometimes, the electrical signals are stationary and periodic, then standard tools such as a discrete Fourier transform and its variants and wavelets are the best choice.Otherwise, more powerful mathematical tools have to be used.
In this paper, we propose the application of the EMD based on the calculation of the HHT [26] able to extract the main features of a non-linear, non-stationary, non-periodic electrical signal.Recently, the EMD and its variants were used as tools to study the main features of several kinds of signals including electrical ones [27][28][29][30][31][32].In this respect, the HHT-EMD analysis presents several advantages with respect to the FFT one that is able to give a correct quantitative description only for linear, stationary and periodic systems.The HHT-EMD analysis has been used on signals collected by the distribution system operator and related to several distribution lines in three different scenarios: (i) line without grid-connected PV plant; (ii) line with a very low PV power; (iii) line with a great amount of PV power (as monitored in [33]).In order to evaluate the effectiveness of the HHT on electrical power signals we compare HHT results to those of the FFT analysis.The analysis proposed in [34] for seasonal periods is now extended to the whole period of the data sets covering about nine months.
In addition, in this work the HHT-EMD investigation is supported by the definition of three indices able to extract the main features of electrical power signals.Indeed, to the best of our knowledge, in the literature the most representative index introduced to study the behavior of signals, including the electrical ones, is the index of stationarity.However, this index is expressed either in terms of the Hilbert spectrum and the Hilbert marginal spectrum [35] or from different distance measures between time-frequency representations to detect abrupt changes in non-stationary noisy signals [36].Instead, the indices introduced in this study are quite different.Firstly, they are defined starting from the electrical signals themselves in terms of statistical averages.Secondly, they do not give any direct information about the non-linear behavior of the electrical signals but they provide direct information about the degree of coherence, the time-correlation and the periodicity, whereas the degree of the non-linearity is obtained indirectly.The three proposed indices are not quantitatively related to the Energies 2016, 9, 211 3 of 18 degree of stationarity.The information contained in these indices completes the HHT-EMD analysis, while the comparison between the calculated values of these indices is useful to choose the proper algorithm (for example, either the FFT or HHT-EMD technique) for the signal under investigation.
The paper is organized as follows: Section 2 recalls the formulation of the EMD tool based on the HHT and introduces the proposed indices useful to characterize quantitatively the electrical signal.In Section 3 the system being tested is described and the measured power data are presented.Section 4 discusses the results obtained by the proposed approach giving special emphasis to the calculation of the indices applied to the measured power data for operating lines belonging to an SG.Finally, Section 5 concludes the paper.

Theoretical Framework
In this Section, the well-known EMD based on the HHT formalism is introduced, recalling the main definitions known from the literature.Secondly, three indices to quantitatively describe the electrical power signal features and to allow a more definite classification of power lines are proposed and briefly discussed.

HHT and EMD
The theoretical investigation is based on the well-known HHT formalism consisting of the calculation of the Hilbert transform Y(t) = H{X(t)} of a time-dependent signal X(t): Yptq " where P denotes the principal value.The Hilbert complex extension of X(t) gives the analytical signal associated to X(t), Z ptq " X ptq `j Y ptq " X ptq `j H tXptqu " a ptq e j θptq , where j is the imaginary unit and where the amplitude a(t) and the phase θ ptq are expressed by a ptq " " X 2 ptq `Y2 ptq ‰ 1 2 and θ ptq " arctan " Yptq Xptq ı , respectively.By applying the stationary phase method to the Fourier transform F(ω) of Z(t) yields ω ptq " dθptq dt , where ω(t) is the instantaneous angular frequency [37].According to the EMD tool [26], the signal X(t) is decomposed into a finite collection of individual characteristic intrinsic modes named Intrinsic Mode Functions (IMFs), viz.X ptq " N ř l"1 c l ptq `rN ptq, where c l (t) is a time-dependent signal with variable amplitude and r N (t) is the residue resulting from the sifting process.For each IMF coefficient c l (t) characterized by an angular frequency ω l (t) it is possible to write the complex representation: with a l (t) the amplitude of z l (t) and θ l ptq " ş ω l ptq dt, l " 1, 2, ...N.The original amplitude X(t) of the electrical signal expanded in a Fourier representation via the IMF components can be written in the form: where Re denotes the real part, f l (t)= ω l (t)/2π is the instantaneous frequency of the l-th IMF and where the residue r N is neglected for instance when it is a monotonic function from which no IMFs can be extracted.
In our framework, the electrical signal X(t) is represented either by the active power P(t) or by the reactive power Q(t).Therefore, the Hilbert transform expressed by Equation ( 1) and the Fourier representation given in Equation ( 3) are specifically expressed in terms either of P(t) or of Q(t) in place of the generic variable X(t).Equation (3) represents the Hilbert amplitude spectrum H(ω,t) or, more simply, the well-known Hilbert spectrum.

The Three Indices for an Electrical Signal
In this subsection, a formulation of characteristic indices of an electrical signal is proposed.The proposed indices are defined by taking into account the main features of the electrical signal without the need to know the corresponding Hilbert spectrum as is often done, for example, for computing the stationary index [35].In particular, the three indices refer to the deviation from a purely periodical behavior and to the coherence degree of the electrical signal giving only indirect information about the degree of non-linearity.The computation of these indices is important because it allows the extraction of general features of the power signal depending on the concentration of PV plants and on the type of users of a given power line.
In order to give a quantitative analysis of the periodicity and coherence of the electrical signal, three indices for classifying the active and reactive powers signal expressing, (1) the deviation from periodicity (DP), (2) the coherence index (CI), (3) the phase coherence (PC), have been defined.The DP and CI are time-dependent, while PC is a real and positive number.
where < . . .> denotes the time average over the whole period.Equation (4) measures the deviation of the power signal from a purely periodic function during the total period T. The angular frequency ω of the cosine function is the frequency of the first harmonic of the FFT of the power signal that can be regarded as the most representative frequency of the time-dependent signal.Concerning (2), let us consider two generic signals, S 1 (t) and S 2 (t) (with S 1 (t), S 2 (t) = P(t), Q(t)).The index CI expresses the time coherence degree and it can be written as a time-correlation function given by the convolution between the two signals [38], namely: In order to compute the indices DP and CI (see Equations ( 4) and ( 5)), the electrical signals have to be normalized to their maximum values.
The PC index computes the coherence degree between two different electrical signals (one of them can be a sinusoid) [39].This index is a number ranging between 0 and 1 and it is expressed in terms of the FFT of the phases of the two electrical signals as follows: where Φ 1 and Φ 2 are the phases summed over the time instants i.
The numerical computation of the indices given in Equations ( 4)-( 6) is presented and exploited in Section 4 for the power line datasets introduced in Section 3.

Description of the System Being Tested and Measured Power Data
In this Section, the system being tested is described in detail and the measured active and reactive power data are presented.The system is an SG located in the south of Italy and constituted by several passive and active lines.These lines feed both residential and commercial users.Each line absorbs a peak mean power in the range 50/350 kW over a length which varies from 250 to 470 m.In the first place, the analysis was performed on seven lines with different amounts of PV power.However, we have been able to collect the lines with similar features.Finally, at the end of the analysis, for the aims of this paper, we have chosen three representative lines, one active and two passive, with important load variations during the day and the week.The line L 1 is passive with a peak power of about 100 kW and no PV power.The line L 2 is passive, it has a very low amount of PV power (1 grid-connected PV plant, 6 kWp) and it has a peak-absorbed power of 50 kW.The line L 3 is an active line with a great amount of PV power (over 100 kWp) and with a maximum peak of about 80 kW of absorbed power.Power measurements have a sampling period of 10 minutes and have been captured between the end of September 2013 and the end of June 2014 for a total of 278 days with 144 samples per day (278 ˆ144 + 1 = 40033 events recorded), corresponding approximately to 9 months.In particular, the line L 1 has a great amount of residential users, the line L 2 has a prevalence of commercial users, and the line L 3 has several residential users also fed with grid-connected PV plants.
It is important to note that in the region where the power data were captured (the south of Italy) the relevant months for the customers, especially the ones fed with grid-connected PV plants, are May and June that can be considered for all practical purposes summer months.First, the duration of the day is longer especially during June, and second, the share of air conditioning in demand can be comparable or even higher with respect to that of the central summer months like July and August.For these reasons, we expect that the present analysis covering a period of 9 months is realistic to extract the described general trend that could be obtained for a 12 months period.
In Figure 1a-c, the active powers for L 1 , L 2 , and L 3 for the total sampling period are depicted.The corresponding reactive powers are displayed in Figure 1d-f.To highlight the main information, the samples of both the active and reactive powers contained in six hours have been grouped in only one sample equal to the mean value for a sake of a clear figure, even if the whole analysis is based on a 10 min sampling time.Due to the high DG penetration, the active line L 3 exhibits negative values of the active absorbed power, and this means that the line injects energy into the grid.It is worth noting that the negative values are concentrated in the first part, for the month of September, and especially towards the end, for the months of May-June, just when the daylight hours of the day increase; in these months L 1 and L 2 show a different behavior, characterized by narrower oscillations.During the autumn and winter months, the line L 3 becomes a globally passive line, even if there are some time-limited negative values.Furthermore, this means that the three lines are globally inductive and L 2 shows an increasing trend towards the end, probably due to the intensive use of air conditioners by commercial users.

Description of the System Being Tested and Measured Power Data
In this Section, the system being tested is described in detail and the measured active and reactive power data are presented.The system is an SG located in the south of Italy and constituted by several passive and active lines.These lines feed both residential and commercial users.Each line absorbs a peak mean power in the range 50/350 kW over a length which varies from 250 to 470 m.In the first place, the analysis was performed on seven lines with different amounts of PV power.However, we have been able to collect the lines with similar features.Finally, at the end of the analysis, for the aims of this paper, we have chosen three representative lines, one active and two passive, with important load variations during the day and the week.The line L1 is passive with a peak power of about 100 kW and no PV power.The line L2 is passive, it has a very low amount of PV power (1 grid-connected PV plant, 6 kWp) and it has a peak-absorbed power of 50 kW.The line L3 is an active line with a great amount of PV power (over 100 kWp) and with a maximum peak of about 80 kW of absorbed power.Power measurements have a sampling period of 10 minutes and have been captured between the end of September 2013 and the end of June 2014 for a total of 278 days with 144 samples per day (278 × 144 + 1 = 40033 events recorded), corresponding approximately to 9 months.In particular, the line L1 has a great amount of residential users, the line L2 has a prevalence of commercial users, and the line L3 has several residential users also fed with grid-connected PV plants.
It is important to note that in the region where the power data were captured (the south of Italy) the relevant months for the customers, especially the ones fed with grid-connected PV plants, are May and June that can be considered for all practical purposes summer months.First, the duration of the day is longer especially during June, and second, the share of air conditioning in demand can be comparable or even higher with respect to that of the central summer months like July and August.For these reasons, we expect that the present analysis covering a period of 9 months is realistic to extract the described general trend that could be obtained for a 12 months period.
In Figure 1a-c, the active powers for L1, L2, and L3 for the total sampling period are depicted.The corresponding reactive powers are displayed in Figure 1d-f.To highlight the main information, the samples of both the active and reactive powers contained in six hours have been grouped in only one sample equal to the mean value for a sake of a clear figure, even if the whole analysis is based on a 10 min sampling time.Due to the high DG penetration, the active line L3 exhibits negative values of the active absorbed power, and this means that the line injects energy into the grid.It is worth noting that the negative values are concentrated in the first part, for the month of September, and especially towards the end, for the months of May-June, just when the daylight hours of the day increase; in these months L1 and L2 show a different behavior, characterized by narrower oscillations.During the autumn and winter months, the line L3 becomes a globally passive line, even if there are some time-limited negative values.Furthermore, this means that the three lines are globally inductive and L2 shows an increasing trend towards the end, probably due to the intensive use of air conditioners by commercial users.As shown in Figure 1, the waveform of the active and reactive powers appears non-periodic and noisy.Moreover, it is also evident the deformed wave-profile with pointed crests and troughs expression of a highly non-linear oscillation.

Results and Discussion
In this Section, we discuss the main results found by applying the HHT-EMD and the indices presented in Section 2 to the power lines datasets.For a complete description of the frequency spectrum associated to all the power lines, a comparison between the HHT, the FFT and the WT has been performed.It has been found that the HHT-EMD analysis takes into account of non-periodical and non-linear features of an electrical signal and it allows a wider spectrum of modes to be found that could not be extracted by limiting the investigation to the FFT and WT techniques.The computation of the indices adds important information to the HHT-EMD for a more complete understanding of the power lines under study.

FFT Results
In Figure 2, the FFT spectra of the lines for both active and reactive powers are shown.For the three power lines studied, four main peaks are associated with the different spectral amplitudes of the main harmonics H1, H2, H3 and H4 both for the active and the reactive powers.The peaks are observable at the frequencies 11.50 × 10 −6 Hz, 23.10 × 10 −6 Hz, 34.70 × 10 −6 Hz and 46.30 × 10 −6 Hz corresponding to 24 h, 12 h, 8 h and 6 h, respectively.We can observe that L2 exhibits an additional peak at a frequency of 1.66 × 10 −6 Hz (1 week) related to another operation mode.We attribute the nature of this mode to the commercial loads, which are massively present in L2.We also note that the PV line L3 shows a similar frequency behavior for the active and reactive powers (see Figure 2c,f) with a main first harmonic and decreasing amplitude of the higher harmonics.As shown in Figure 1, the waveform of the active and reactive powers appears non-periodic and noisy.Moreover, it is also evident the deformed wave-profile with pointed crests and troughs expression of a highly non-linear oscillation.

Results and Discussion
In this Section, we discuss the main results found by applying the HHT-EMD and the indices presented in Section 2 to the power lines datasets.For a complete description of the frequency spectrum associated to all the power lines, a comparison between the HHT, the FFT and the WT has been performed.It has been found that the HHT-EMD analysis takes into account of non-periodical and non-linear features of an electrical signal and it allows a wider spectrum of modes to be found that could not be extracted by limiting the investigation to the FFT and WT techniques.The computation of the indices adds important information to the HHT-EMD for a more complete understanding of the power lines under study.

FFT Results
In Figure 2, the FFT spectra of the lines for both active and reactive powers are shown.For the three power lines studied, four main peaks are associated with the different spectral amplitudes of the main harmonics H 1 , H 2 , H 3 and H 4 both for the active and the reactive powers.The peaks are observable at the frequencies 11.50 ˆ10 ´6 Hz, 23.10 ˆ10 ´6 Hz, 34.70 ˆ10 ´6 Hz and 46.30 ˆ10 ´6 Hz corresponding to 24 h, 12 h, 8 h and 6 h, respectively.We can observe that L 2 exhibits an additional peak at a frequency of 1.66 ˆ10 ´6 Hz (1 week) related to another operation mode.We attribute the nature of this mode to the commercial loads, which are massively present in L 2 .We also note that the PV line L 3 shows a similar frequency behavior for the active and reactive powers (see Figure 2c,f) with a main first harmonic and decreasing amplitude of the higher harmonics.

Wavelet Results
In order to make a quantitative comparison with the FFT results on power signals and to check whether taking into account also the non-stationary behavior can lead to different results, we have also used the complex Morlet wavelet mother as prototype of WT.For technical details about the basic formalism, see [25].Note that this choice is not restrictive because other kinds of WT existing in literature would lead to very similar results.The calculated time-frequency behavior for both active and reactive powers is illustrated in Figure 3.
The frequencies extracted from the time-frequency domain are 11.55 × 10 −6 Hz, 23.10 × 10 −6 Hz, 34.70 × 10 −6 Hz for the three lines.Looking at Figure 3 the two FFT frequencies corresponding to the harmonics H1 and H2 (24 h and 12 h, respectively) are found also according to the WT-based method both for active and reactive powers for lines L1 and L3, while for line L2 only the one corresponding to 12 h was computed.Instead, the frequencies corresponding to the H3 and H4 harmonics (8 h and 6 h, respectively) lack in the WT analysis for the line L3, but the frequency corresponding to H3 compares well to the one found with WT analysis in lines L1 (small value) and L2.These results suggest that the non-stationary behavior included within the WT analysis does not lead to essential changes of the frequencies corresponding to the main H1 and H2 harmonics.On the other hand, a limited additional frequency present in the FFT spectra for the line L2 corresponding to a week's activity was also found.
Nevertheless, despite its solid basis, the above preliminary analysis performed according to the FFT and to the WT techniques does not allow the full characterization of the different lines-for instance, in terms of their non-linear behavior that is in turn strictly connected with the deviation from periodicity for an electrical power signal.For this reason, as a first further step the HHT analysis combined with the EMD has been applied to the three lines and the obtained results have been compared to the ones derived by means of the FFT and WT analysis.However, since the HHT-EMD tool did not completely quantitatively characterize the features of the power signal, the

Wavelet Results
In order to make a quantitative comparison with the FFT results on power signals and to check whether taking into account also the non-stationary behavior can lead to different results, we have also used the complex Morlet wavelet mother as prototype of WT.For technical details about the basic formalism, see [25].Note that this choice is not restrictive because other kinds of WT existing in literature would lead to very similar results.The calculated time-frequency behavior for both active and reactive powers is illustrated in Figure 3.
The frequencies extracted from the time-frequency domain are 11.55 ˆ10 ´6 Hz, 23.10 ˆ10 ´6 Hz, 34.70 ˆ10 ´6 Hz for the three lines.Looking at Figure 3 the two FFT frequencies corresponding to the harmonics H 1 and H 2 (24 h and 12 h, respectively) are found also according to the WT-based method both for active and reactive powers for lines L 1 and L 3 , while for line L 2 only the one corresponding to 12 h was computed.Instead, the frequencies corresponding to the H 3 and H 4 harmonics (8 h and 6 h, respectively) lack in the WT analysis for the line L 3 , but the frequency corresponding to H 3 compares well to the one found with WT analysis in lines L 1 (small value) and L 2 .These results suggest that the non-stationary behavior included within the WT analysis does not lead to essential changes of the frequencies corresponding to the main H 1 and H 2 harmonics.On the other hand, a limited additional frequency present in the FFT spectra for the line L 2 corresponding to a week's activity was also found.
Nevertheless, despite its solid basis, the above preliminary analysis performed according to the FFT and to the WT techniques does not allow the full characterization of the different lines-for instance, in terms of their non-linear behavior that is in turn strictly connected with the deviation from periodicity for an electrical power signal.For this reason, as a first further step the HHT analysis combined with the EMD has been applied to the three lines and the obtained results have been compared to the ones derived by means of the FFT and WT analysis.However, since the HHT-EMD tool did not completely quantitatively characterize the features of the power signal, the new suggested indices introduced in Section 2 have been computed.Note that the calculation of the three indices expressed in Equations ( 4)-( 6) is the main step of our investigation.
Energies 2016, 9, 211 8 of 17 new suggested indices introduced in Section 2 have been computed.Note that the calculation of the three indices expressed in Equations ( 4)-( 6) is the main step of our investigation.

HHT-EMD Results Compared to FFT and Wavelet Results: Calculation of Frequencies of IMFs
According to the EMD analysis based on the HHT technique, the electrical signal is decomposed into 14 IMFs for both the active and the reactive powers of the three lines.As an example, Figure 4 shows the time behavior of the 14 IMFs of the P3 active power for the whole sampling period, labeled as IMFi with I = 1, 2, … 14 in decreasing frequency order.The corresponding amplitudes of the IMFs related to the other active powers and reactive powers have a similar trend.
Looking at Figure 4, the number of oscillations of the IMFs decreases as the IMF order increases.This depends on the procedure used for the iterative extractions of each IMF.In fact, each step of EMD provides a new signal equal to the subtraction between the original one and the mean value of the envelope curves of its local maximum and minimum values.At the end of the process, it retains only the residual non-oscillatory term rN (not shown).
Looking at the amplitudes of the IMFs it can be seen that there are no changes in the oscillation frequency, confirming that for each IMF amplitude there are no superpositions of different monocomponents.The IMFs 5-7 (denoted with arrows in Figure 4) have frequencies comparable to those corresponding to the FFT harmonics and also to some of the WT frequencies.The amplitudes averaged over the whole period have zero mean values, as should be expected.

HHT-EMD Results Compared to FFT and Wavelet Results: Calculation of Frequencies of IMFs
According to the EMD analysis based on the HHT technique, the electrical signal is decomposed into 14 IMFs for both the active and the reactive powers of the three lines.As an example, Figure 4 shows the time behavior of the 14 IMFs of the P 3 active power for the whole sampling period, labeled as IMF i with I = 1, 2, . . .14 in decreasing frequency order.The corresponding amplitudes of the IMFs related to the other active powers and reactive powers have a similar trend.
Looking at Figure 4, the number of oscillations of the IMFs decreases as the IMF order increases.This depends on the procedure used for the iterative extractions of each IMF.In fact, each step of EMD provides a new signal equal to the subtraction between the original one and the mean value of the envelope curves of its local maximum and minimum values.At the end of the process, it retains only the residual non-oscillatory term r N (not shown).
Looking at the amplitudes of the IMFs it can be seen that there are no changes in the oscillation frequency, confirming that for each IMF amplitude there are no superpositions of different monocomponents.The IMFs 5-7 (denoted with arrows in Figure 4) have frequencies comparable to those corresponding to the FFT harmonics and also to some of the WT frequencies.The amplitudes averaged over the whole period have zero mean values, as should be expected.The frequencies of the IMFs according to the average procedure described below are summarized in Table 1.As a general remark, it is evident that there is a greater number of frequencies calculated via the HHT-EMD than the ones determined by means of the FFT and WT analysis.This is not surprising because it is strictly connected mainly to the non-periodical behavior of the power signals in all the lines with different degrees of deviation from periodicity (see Section 4 for details) and to a lesser extent also to the non-stationary behavior.The deviation from a purely periodical behavior is in turn related to the non-linear behavior of power signals in DLs of SGs.In particular, the HHT approach combined with the EMD procedure allows the extraction of additional modes compared to the FFT.Finally, just as EMD is based on the envelope curves of local maximum and minimum values, IMFs with lower frequencies can be extracted before IMFs with higher frequencies.For example, The frequencies of the IMFs according to the average procedure described below are summarized in Table 1.As a general remark, it is evident that there is a greater number of frequencies calculated via the HHT-EMD than the ones determined by means of the FFT and WT analysis.This is not surprising because it is strictly connected mainly to the non-periodical behavior of the power signals in all the lines with different degrees of deviation from periodicity (see Section 4 for details) and to a lesser extent also to the non-stationary behavior.The deviation from a purely periodical behavior is in turn related to the non-linear behavior of power signals in DLs of SGs.In particular, the HHT approach combined with the EMD procedure allows the extraction of additional modes compared to the FFT.Finally, just as EMD is based on the envelope curves of local maximum and minimum values, IMFs with lower frequencies can be extracted before IMFs with higher frequencies.For example, looking at the first row of Table 1 for P 1 , it can be observed that IMF 10 has 0.08 ˆ10 ´6 Hz, whereas the IMF 11 and IMF 12 have higher frequencies, equal to 0.13 ˆ10 ´6 Hz and 0.12 ˆ10 ´6 Hz, respectively.Note that the frequencies shown in Table 1 are not properly the instantaneous frequencies expressed in Equation ( 3) that by definition depend on time as is typically shown in the so-called Hilbert spectrum.As a matter of fact, in order to compare the obtained frequencies with the corresponding ones computed via FFT and WT techniques, we have performed an average of the instantaneous IMF frequencies over the total sampling period.This procedure recalls the one employed for the calculation of the Hilbert average marginal spectrum obtained from the Hilbert spectrum by integrating it over the total sampling period and by dividing the result of this integration by the total sampling period.Of course, this analogy is no longer valid when the instantaneous frequency of the tone does not oscillate in a very narrow band so there could be points assigned to a different frequency in the Hilbert marginal spectrum.However, we would also like to stress that the Hilbert marginal spectrum has relevant meaning in the analysis of the HHT frequencies, is strictly connected to the HHT analysis and also used for extracting the general features of an electrical power signal.The advantage of this procedure is that it is possible to directly compare these mean single frequency values with the ones obtained with the other two techniques described above.On the other hand, this approach does not change the discussion since we have found that the instantaneous IMF frequencies extracted from the EMD technique do not appreciably vary with time despite the daily variations of loads in the different power lines.The weak variation as a function of time also characterizes the frequencies calculated by means of the WT analysis (see Figure 3).
In order to clarify better the following discussion, we distinguish for every line three different frequency regions of the spectrum.In Figure 5a, the three frequency regions are shown by plotting the frequencies of the IMFs related to the three active powers.In addition to the central region, covering the frequency range of the typical harmonics found in the FFT spectrum, there are other two regions, the upper and the lower one.The lower region includes low frequencies with values up to about 5 ˆ10 ´6 Hz and, apart from one exception, the found frequencies are not present in the FFT spectrum.The central region has frequencies included in the range between about 10 ˆ10 ´6 Hz and about 50 ˆ10 ´6 Hz and corresponds more or less to the ones of the FFT spectrum.In the upper region there are frequencies having values above 50 ˆ10 ´6 Hz up to approximately 450 ˆ10 ´6 Hz; these frequencies are not present in the FFT spectrum.This wide range of frequencies computed via the HHT calculation is the combined result of the deviation from a purely periodical behavior, in turn strictly related to non-linearity effects that cannot be taken into account by the FFT analysis.For a detailed quantitative discussion on this point see next Section 4, where the calculated indices presented in Section 2 are discussed.
We now quantitatively discuss the frequencies related to the active powers.For the L 1 passive line and the L 3 active line the HHT frequencies differ appreciably especially in the upper and central regions (see also Table 1).This different trend is not surprising because of the strong difference in PV plants feeding the two lines which leads to different values of the active powers with P 3 assuming also negative values (see Figure 1c).Moreover, frequencies referred to the active powers of the L 2 passive line have overall different values with respect to those of P 1 and P 3 .
More specifically, it is useful to make a quantitative comparison among the frequencies belonging to the central region with the corresponding frequencies of the typical FFT harmonics for P 1 , P 2 and P 3 active powers.The comparison between the HHT-EMD frequencies and the FFT harmonics is illustrated in Figure 5b.For the first passive line L 1 , the two first most representative frequencies 8.50 and 17.28 ˆ10 ´6 Hz (green circles) strongly differ from the frequencies of the H 1 (11.50 ˆ10 ´6 Hz) and H 2 (23.10 ˆ10 ´6 Hz) harmonics, respectively determined via the FFT (black squares), while the third one at 32.65 ˆ10 ´6 Hz has a much closer value to the frequency (34.70 ˆ10 ´6 Hz) of the H 3 harmonic.On the other hand, there is no trace of the one corresponding to the 6 h peak as found according to the FFT.By contrast, the HHT-EMD frequencies corresponding to those of the H1, H2 and H4 (46.30 × 10 −6 Hz) harmonics for the second passive line L2 (red up triangles) are rather close to the FFT ones, while the frequency related to the H3 harmonic is not present in the HHT calculation.The frequencies obtained for the L3 active line are slightly different with respect to the FFT ones especially those corresponding to the H1 and H2 harmonics (blue down triangles).Also for the L3 active line there is the lack of the mode corresponding to the frequency of the H3 harmonic and a frequency at 88.64 × 10 −6 Hz not found according to the FFT calculation (it would be the equivalent of the frequency of the H5 harmonic).This frequency corresponds to 3 h.The same conclusions can be drawn if the HHT-EMD frequencies are compared to the WT frequencies, which are superimposable to the ones of the FFT.Special attention deserves also the analysis of other characteristic frequencies belonging to the IMFs in the lower region of the spectrum.In particular, in addition to several frequencies not present in the FFT spectrum, the HHT analysis also confirms, for the L2 passive line, the frequency (1.29 × 10 −6 Hz) corresponding to 1 week found via the FFT analysis, which is typical of commercial loads.Concerning the frequencies of the IMFs above 50 × 10 −6 Hz in the upper region, these would correspond approximately for every line to 2 h, 1 h and about ½ h, respectively.The above-mentioned additional features typical of the HHT spectrum and related to a combination of non-linear effects, deviation from a periodical behavior and non-stationary behavior were not extracted from the FFT and WT analysis.
To check further the findings obtained for the P(t) signal we have processed the same datasets for the corresponding reactive powers Q(t) illustrated in Table 1 where the Q2 reactive power has 13 IMFs frequencies and the Q3 has 4 additional IMFs in the lower frequency region (not shown).The frequencies of the IMFs are also plotted in Figure 6a for Q(t) and subdivided into the three regions as done in Figure 5a.In Figure 6b the FFT frequencies of the main FFT harmonics related to the reactive powers are compared to the HHT corresponding ones in the central region of the spectrum.On the other hand, there is no trace of the one corresponding to the 6 h peak as found according to the FFT.By contrast, the HHT-EMD frequencies corresponding to those of the H 1 , H 2 and H 4 (46.30ˆ10 ´6 Hz) harmonics for the second passive line L 2 (red up triangles) are rather close to the FFT ones, while the frequency related to the H 3 harmonic is not present in the HHT calculation.The frequencies obtained for the L 3 active line are slightly different with respect to the FFT ones especially those corresponding to the H 1 and H 2 harmonics (blue down triangles).Also for the L 3 active line there is the lack of the mode corresponding to the frequency of the H 3 harmonic and a frequency at 88.64 ˆ10 ´6 Hz not found according to the FFT calculation (it would be the equivalent of the frequency of the H 5 harmonic).This frequency corresponds to 3 h.The same conclusions can be drawn if the HHT-EMD frequencies are compared to the WT frequencies, which are superimposable to the ones of the FFT.Special attention deserves also the analysis of other characteristic frequencies belonging to the IMFs in the lower region of the spectrum.In particular, in addition to several frequencies not present in the FFT spectrum, the HHT analysis also confirms, for the L 2 passive line, the frequency (1.29 ˆ10 ´6 Hz) corresponding to 1 week found via the FFT analysis, which is typical of commercial loads.Concerning the frequencies of the IMFs above 50 ˆ10 ´6 Hz in the upper region, these would correspond approximately for every line to 2 h, 1 h and about ½ h, respectively.The above-mentioned additional features typical of the HHT spectrum and related to a combination of non-linear effects, deviation from a periodical behavior and non-stationary behavior were not extracted from the FFT and WT analysis.
To check further the findings obtained for the P(t) signal we have processed the same datasets for the corresponding reactive powers Q(t) illustrated in Table 1 where the Q 2 reactive power has 13 IMFs frequencies and the Q 3 has 4 additional IMFs in the lower frequency region (not shown).The frequencies of the IMFs are also plotted in Figure 6a for Q(t) and subdivided into the three regions as done in Figure 5a.In Figure 6b the FFT frequencies of the main FFT harmonics related to the reactive powers are compared to the HHT corresponding ones in the central region of the spectrum.Looking at Table 1 and Figure 6a, the frequencies of the IMFs corresponding to the Q1 and Q2 reactive powers are rather close each other due to the similar behavior of the signal (see Figure 1d,e), apart from some discrepancies which result from the difference between PV plants characterizing the L1 and the L2 passive lines.Instead, due to the different behavior of Q3 (see Figure 1f), most of the calculated frequencies of the IMFs are different with respect to those of L1 and L2.Unlike the FFT analysis that gives very similar results, the comparison of the frequencies related to the active and the reactive powers, leads to different results.Indeed, there are significant differences between the frequencies related to the active powers and the corresponding ones associated to the reactive powers especially in the central region of the IMF frequencies.For example, the frequencies corresponding to those of the H2 and H3 harmonics typical of the FFT spectra were not found (see Figure 6b).The frequencies 12.08 × 10 −6 Hz and 11.91 × 10 −6 Hz rather close to that of the H1 harmonic (11.50 × 10 −6 Hz) were found only for Q1 and Q2 belonging to the lines L1 and L2, respectively but not for Q3 belonging to the line L3.Even though the HHT calculated frequencies 51.55 × 10 −6 Hz, 51.72 × 10 −6 Hz and 54.83 × 10 −6 Hz for Q1, Q2 and Q3 are higher than the frequency of the typical FFT H4 harmonic (46.30 × 10 −6 Hz), they can be considered their corresponding HHT frequencies.The found discrepancies can be ascribed to the more accurate analysis given by HHT.
In fact, the HHT computation is able to highlight the local properties of a typical power signal related to the effects of non-linearity that in turn give rise to its not-fully periodical character.In this respect, we believe that the deviation from periodicity of the power signal (for a quantitative analysis see Section 4) could lead to the lack of some of the frequencies in the HHT spectrum and especially the ones corresponding to the higher FFT harmonics.
We also performed an analysis of the data for single year's quarters based on the HHT and some results are different (here not shown).For more details about this seasonal analysis see [34].

The Indices: Numerical Evaluation
In this subsection, the results of the numerical computation of the indices expressed by Equations ( 4)-( 6) are presented.In Figure 7 the DP index defined in Equation ( 4) and calculated for Looking at Table 1 and Figure 6a, the frequencies of the IMFs corresponding to the Q 1 and Q 2 reactive powers are rather close each other due to the similar behavior of the signal (see Figure 1d,e), apart from some discrepancies which result from the difference between PV plants characterizing the L 1 and the L 2 passive lines.Instead, due to the different behavior of Q 3 (see Figure 1f), most of the calculated frequencies of the IMFs are different with respect to those of L 1 and L 2 .Unlike the FFT analysis that gives very similar results, the comparison of the frequencies related to the active and the reactive powers, leads to different results.Indeed, there are significant differences between the frequencies related to the active powers and the corresponding ones associated to the reactive powers especially in the central region of the IMF frequencies.For example, the frequencies corresponding to those of the H 2 and H 3 harmonics typical of the FFT spectra were not found (see Figure 6b).The frequencies 12.08 ˆ10 ´6 Hz and 11.91 ˆ10 ´6 Hz rather close to that of the H 1 harmonic (11.50 ˆ10 ´6 Hz) were found only for Q 1 and Q 2 belonging to the lines L 1 and L 2 , respectively but not for Q 3 belonging to the line L 3 .Even though the HHT calculated frequencies 51.55 ˆ10 ´6 Hz, 51.72 ˆ10 ´6 Hz and 54.83 ˆ10 ´6 Hz for Q 1 , Q 2 and Q 3 are higher than the frequency of the typical FFT H 4 harmonic (46.30ˆ10 ´6 Hz), they can be considered their corresponding HHT frequencies.The found discrepancies can be ascribed to the more accurate analysis given by HHT.
In fact, the HHT computation is able to highlight the local properties of a typical power signal related to the effects of non-linearity that in turn give rise to its not-fully periodical character.In this respect, we believe that the deviation from periodicity of the power signal (for a quantitative analysis see Section 4) could lead to the lack of some of the frequencies in the HHT spectrum and especially the ones corresponding to the higher FFT harmonics.
We also performed an analysis of the data for single year's quarters based on the HHT and some results are different (here not shown).For more details about this seasonal analysis see [34].

The Indices: Numerical Evaluation
In this subsection, the results of the numerical computation of the indices expressed by Equations ( 4)-( 6) are presented.In Figure 7 the DP index defined in Equation ( 4) and calculated for the P 1 , P 2 and P 3 active powers is displayed, both for the total sampling period of 278 days in panels (a)-(c).In addition, a zoom of the signal for a temporal window ranging between the 140th and the 150th day (approximately located in the middle) is reported in panels (e)-(g) of Figure 6.As a reference, the coherence degree Coherence cos-cos as a function of time is displayed for a couple of cosine functions (see panel (d) and (h)).As expected, according to the definition of time-correlation functions [38], the general behavior is the decrease of the DP index amplitude with increasing time as illustrated in panels (a)-(d) marking a general reduction of the correlation of a given power signal with a cosine (or sine) function.From the analysis of the oscillation amplitudes, it turns out that the strongest deviation from the purely periodical behavior that is typical of a cosine (or sine) function having indeed the highest coherence degree characterizes the P 2 active power.This can be inferred by looking at the maximum oscillation amplitude of DP which for the P 2 active power is less pronounced (panel (f)) with respect to the ones involving the P 1 and P 3 active powers (panels (e) and (g)).In particular, in the shown temporal window of 10 days, the amplitude variation is about 750 for DP P 2 , about 1700 for DP P 1 and more than 4000 for DP P 3 .The oscillation amplitudes of DP P 1 and DP P 3 are much closer to the one characterizing the Coherence cos´cos (panel (h)) between two cosine functions thus marking their higher periodicity degree.This behavior is very similar also by focusing on other temporal windows.The conclusions drawn for the DP index computed for the active powers remain valid also for the corresponding reactive powers.Therefore, the quantitative analysis on the DP index indicates that the L 1 passive line and the L 3 active line are characterized by a signal deviating to a lesser extent from the purely periodical behavior.First, from an overall analysis looking at Figure 7 it is important to note that the periodicity degree of the power signal does not strictly depend on the injected PV power, but it is strongly related to the loads and to the operating conditions.In particular, in this study we observe that the lines characterized by power signals with a smaller deviation from a purely periodical behavior are the ones feeding residential users, namely the L 1 passive line and the L 3 active line, one with no PV power and the other with a high PV power.However, from a more specific and deeper analysis it can be seen, by comparing panels (e), (g) and (h) of Figure 7, that the L 3 active line having the largest injected PV power is the one with the smaller deviation from a purely periodical behavior because it is characterized by the largest oscillation amplitudes.Therefore, the computation of the DP index could be helpful for highlighting the future advantages of planning lines having a high percentage of PV plants.Indeed, we can conclude that the more the concentration of PV plants is higher, the more the power signals tends towards periodical behavior, and this is potentially advantageous for the future operating conditions in DLs of SGs.
one with the smaller deviation from a purely periodical behavior because it is characterized by the largest oscillation amplitudes.Therefore, the computation of the DP index could be helpful for highlighting the future advantages of planning lines having a high percentage of PV plants.Indeed, we can conclude that the more the concentration of PV plants is higher, the more the power signals tends towards periodical behavior, and this is potentially advantageous for the future operating conditions in DLs of SGs.In Figure 8 the CI between a given couple of active powers defined in Equation ( 5) is shown in normalized units.According to its definition, resulting from the time-correlation function and similarly to the DP index also the amplitude of the CI reduces with increasing time.In panels (d)-(f) of Figure 7 the CI is shown in a temporal window ranging between the 140th and the 150th day.It can be noted that the CI trend as a function of time is strongly altered with an important amplitude reduction when the coherence is computed between the P 1 and P 2 and between the P 2 and P 3 active powers.Moreover, in this case there is also a reduced CI amplitude.Instead, CI P 1 -P 3 exhibits regular oscillations marking a stronger coherence between the P 1 and P 3 active powers.The CI for the reactive powers shows very similar behavior with respect to that of the active powers, as occurred for the DP index.This can be understood by taking into account that the L 1 and L 3 lines feed residential users, while the L 2 line feeds commercial users.By means of this analysis we are able to quantify the coherence degree of two given power lines in different temporal windows and this can be done for any given available power line.As expected, the coherence degree resulting from the two residential lines, one active and the other passive, is higher than that computed between a residential line and a commercial line (one active and one passive or two passive lines).In Figure 8 the CI between a given couple of active powers defined in Equation ( 5) is shown in normalized units.According to its definition, resulting from the time-correlation function and similarly to the DP index also the amplitude of the CI reduces with increasing time.In panels (d)-(f) of Figure 7 the CI is shown in a temporal window ranging between the 140th and the 150th day.It can be noted that the CI trend as a function of time is strongly altered with an important amplitude reduction when the coherence is computed between the P1 and P2 and between the P2 and P3 active powers.Moreover, in this case there is also a reduced CI amplitude.Instead, P -P 13 CI exhibits regular In order to validate the previous results the PC defined in Equation ( 6) has also been calculated.The results of this calculation are displayed in Figure 9.The phase coherence between the P 1 , P 2 and P 3 active powers and the cosine function (empty circles) is higher for the two residential lines, being about 0.68 for P 1 -cos and P 3 -cos, while the one related to the commercial line (P 2 -cos) is about 0.54.oscillations marking a stronger coherence between the P1 and P3 active powers.The CI for the reactive powers shows very similar behavior with respect to that of the active powers, as occurred for the DP index.This can be understood by taking into account that the L1 and L3 lines feed residential users, while the L2 line feeds commercial users.By means of this analysis we are able to quantify the coherence degree of two given power lines in different temporal windows and this can be done for any given available power line.As expected, the coherence degree resulting from the two residential lines, one active and the other passive, is higher than that computed between a residential line and a commercial line (one active and one passive or two passive lines).
In order to validate the previous results the PC defined in Equation ( 6) has also been calculated.The results of this calculation are displayed in Figure 9.The phase coherence between the P1, P2 and P3 active powers and the cosine function (empty circles) is higher for the two residential lines, being about 0.68 for P1-cos and P3-cos, while the one related to the commercial line (P2-cos) is about 0.54.These results show that the phase coherence with a co-sinusoidal (or sinusoidal) electrical signal is higher for the residential lines confirming the previous predictions.Moreover, the phase coherence between P1 and P3 (full circles) is equal to 0.97, while the ones between P1 and P2 and P2 and P3 active powers are about 0.90.According to this analysis, it can be established the PC for each given line and it is confirmed that the degree of phase coherence is higher between the two residential lines.Therefore, this index is useful to characterize the type of users of the line and it is not influenced by the quantity of installed PV power.

Conclusions
In summary, it has been shown that the HHT-EMD analysis applied to three different power lines belonging to an SG allows the extraction of several features that were not found according to the FFT and WT investigations, and the quantitative differences between the results obtained by These results show that the phase coherence with a co-sinusoidal (or sinusoidal) electrical signal is higher for the residential lines confirming the previous predictions.Moreover, the phase coherence between P 1 and P 3 (full circles) is equal to 0.97, while the ones between P 1 and P 2 and P 2 and P 3 active powers are about 0.90.According to this analysis, it can be established the PC for each given line and it is confirmed that the degree of phase coherence is higher between the two residential lines.Therefore, this index is useful to characterize the type of users of the line and it is not influenced by the quantity of installed PV power.

Conclusions
In summary, it has been shown that the HHT-EMD analysis applied to three different power lines belonging to an SG allows the extraction of several features that were not found according to the FFT and WT investigations, and the quantitative differences between the results obtained by means of the FFT, WT and HHT techniques have been discussed.The frequencies spectrum calculated via HHT-EMD is wider with respect to the ones found according to the FFT and WT analysis and this is mainly due to the non-periodical behavior of the electrical signals of power lines belonging to an SG in turn connected to their strong deviation from linearity.
An analysis giving important quantitative information about the features of power lines was carried out according to the calculation of new indices that are able to better characterize the electrical signal in terms of deviation from a purely periodical behavior and coherence properties.This is the key result of this study and completes the findings obtained by means of the previous analysis based on the usual mathematical tools.However, the indices proposed are defined according to a statistical approach starting from the electrical signals, and their information strongly differs from that given by the index of stationarity discussed in the literature.
It has been found that the higher the concentration of PV plants, the more the power signal tends towards a periodical behavior, and this can be an important advantage for planning the future operating conditions in DLs of SGs.Moreover, the power trend associated with the two residential lines has a higher periodicity degree than that of the commercial line.Finally, by computing the time and phase coherence indices it was also possible to establish in a quantitative way the coherence degree of two given power lines.Specifically, we have quantified the time and phase coherence degree correlating the lines having either more residential users or more commercial users.It was found that, in general, the time and phase coherence degree is higher between the two lines having more residential users.The proposed approach has helped to highlight the intrinsic features and the behavior of the power lines not evidenced by the study of the time signal according to the HHT-EMD analysis or by the simple application of the FFT and WT techniques.This approach, mainly focusing on the calculation of the proposed indices, can be usefully employed for the optimal operation and planning of SGs in the presence of strongly nonlinear loads and intermittent ones and it can be useful for the characterization of the lines with respect to the connected loads.

Figure 2 .
Figure 2. FFT spectra.(a) FFT of the P1 active power; (b) FFT of the P2 active power; (c) FFT of the P3 active power; (d) FFT of the Q1 reactive power; (e) FFT of the Q2 reactive power and (f) FFT of the Q3 reactive power.

Figure 2 .
Figure 2. FFT spectra.(a) FFT of the P 1 active power; (b) FFT of the P 2 active power; (c) FFT of the P 3 active power; (d) FFT of the Q 1 reactive power; (e) FFT of the Q 2 reactive power and (f) FFT of the Q 3 reactive power.

Figure 3 .
Figure 3. Time-frequency domain plot computed according to the WT-based method.(a) WT of the P1 active power; (b) WT of the P2 active power; (c) WT of the P3 active power; (d) WT of the Q1 reactive power; (e) WT of the Q2 reactive power and (f) WT of the Q3 reactive power.

Figure 3 .
Figure 3. Time-frequency domain plot computed according to the WT-based method.(a) WT of the P 1 active power; (b) WT of the P 2 active power; (c) WT of the P 3 active power; (d) WT of the Q 1 reactive power; (e) WT of the Q 2 reactive power and (f) WT of the Q 3 reactive power.

Figure 4 .
Figure 4. Amplitude of the 14 IMFs as a function of time related to the P3 active power.An arrow indicates the amplitudes of the IMFs having frequencies close to the ones of the harmonics obtained via the FFT analysis (see the text for a discussion).

Figure 4 .
Figure 4. Amplitude of the 14 IMFs as a function of time related to the P 3 active power.An arrow indicates the amplitudes of the IMFs having frequencies close to the ones of the harmonics obtained via the FFT analysis (see the text for a discussion).

Figure 5 .
Figure 5. (a) HHT-EMD frequencies for the active power P of L1, L2 and L3 subdivided into the three frequency regions (upper, central and lower) of the HHT spectrum.Green circles: HHT frequencies for P1.Red up triangles: HHT frequencies for P2.Blue down triangles: HHT-EMD frequencies for P3; (b) main FFT frequencies compared to the HHT-EMD frequencies of the central region for the active power of L1, L2 and L3.Black squares: FFT harmonics H1, H2, H3 and H4.The meaning of the other symbols is the same as in panel (a).

Figure 5 .
Figure 5. (a) HHT-EMD frequencies for the active power P of L 1 , L 2 and L 3 subdivided into the three frequency regions (upper, central and lower) of the HHT spectrum.Green circles: HHT frequencies for P 1 .Red up triangles: HHT frequencies for P 2 .Blue down triangles: HHT-EMD frequencies for P 3 ; (b) main FFT frequencies compared to the HHT-EMD frequencies of the central region for the active power of L 1 , L 2 and L 3 .Black squares: FFT harmonics H 1 , H 2 , H 3 and H 4 .The meaning of the other symbols is the same as in panel (a).

Figure 6 .
Figure 6.(a) HHT-EMD frequencies for the reactive power of L1, L2 and L3 subdivided into the three frequency regions (upper, central and lower) of the HHT spectrum.Green squares: HHT frequencies for Q1.Red diamonds: HHT-EMD frequencies for Q2.Blue up triangles: HHT frequencies for Q3.(b) Main FFT frequencies compared to the HHT-EMD frequencies of the central region for the reactive power of L1, L2 and L3.Black circles: FFT harmonics H1, H2, H3 and H4.The meaning of the other symbols is the same as in panel (a).

Figure 6 .
Figure 6.(a) HHT-EMD frequencies for the reactive power of L 1 , L 2 and L 3 subdivided into the three frequency regions (upper, central and lower) of the HHT spectrum.Green squares: HHT frequencies for Q 1 .Red diamonds: HHT-EMD frequencies for Q 2 .Blue up triangles: HHT frequencies for Q 3 .(b) Main FFT frequencies compared to the HHT-EMD frequencies of the central region for the reactive power of L 1 , L 2 and L 3 .Black circles: FFT harmonics H 1 , H 2 , H 3 and H 4 .The meaning of the other symbols is the same as in panel (a).

Figure 7 .
Figure 7. DP index for the three active powers and for the cosine function.(a) DP P 1 ; (b) DP P 2 ; (c) DP P3 ; (e) DP P1 in a temporal window; (f) DP P2 in a temporal window; (g) DP P 3 in a temporal window; (d) Coherence cos´cos and (h) Coherence cos´cos in a temporal window.

Figure 7 . 3 P
Figure 7. DP index for the three active powers and for the cosine function.(a)

Figure 8 .
Figure 8. CI between couples of active powers.(a) PP 12 CI  ; (b) PP 13 CI  ; (c) PP 23 CI  ; (d) PP 12 CI  in a temporal window; (e) PP 13 CI  in a temporal window; and (f) PP 23 CI  in a temporal window.

Figure 9 .
Figure 9. PC between couples of active powers (full circles) and between a given active power and the cosine function (empty circles).

Figure 9 .
Figure 9. PC between couples of active powers (full circles) and between a given active power and the cosine function (empty circles).

Table 1 .
Frequencies of the 14 IMFs for L 1 , L 2 and L 3 lines.