Next Article in Journal
Influence of Combustion Parameters on Fouling Composition after Wood Pellet Burning in a Lab-Scale Low-Power Boiler
Next Article in Special Issue
Adaptation of Powerline Communications-Based Smart Metering Deployments to the Requirements of Smart Grids
Previous Article in Journal
Parametric Optimization of Regenerative Organic Rankine Cycle System for Diesel Engine Based on Particle Swarm Optimization
Previous Article in Special Issue
A Novel Choice Procedure of Magnetic Component Values for Phase Shifted Full Bridge Converters with a Variable Dead-Time Control Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Application of Spectral Kurtosis to Separate Hybrid Power Quality Events

by
Juan José González De la Rosa
1,2,*,*,
José María Sierra-Fernández
1,2,
José Carlos Palomares-Salas
1,2,
Agustín Agüera-Pérez
1,2 and
Álvaro Jiménez Montero
1,2
1
Research Group Plan Andaluz Investigación Desarrollo Tecnológico e Innovación-Tecnologías de la Información y Comunicación-168 (PAIDI-TIC-168): Computational Instrumentation and Industrial Electronics (ICEI), E-11202 Algeciras, Spain
2
University of Cádiz, Area of Electronics, Polytechnic School of Algeciras, Av. Ramón Puyol, S/N,E-11202 Algeciras, Spain
*
Author to whom correspondence should be addressed.
*
Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +34-956-028-069; Fax: +34-956-028-001.
Energies 2015, 8(9), 9777-9793; https://doi.org/10.3390/en8099777
Submission received: 1 July 2015 / Revised: 18 August 2015 / Accepted: 26 August 2015 / Published: 9 September 2015
(This article belongs to the Special Issue Smart Metering)

Abstract

:
For the development of the future smart grid, the detection of power quality events is a key issue for the power system monitoring. Voltage sags, swells, harmonics (variations) and interruptions, which produce large losses for commercial and industrial consumers, are the main events to be considered due to the sensitivity of equipment to these electrical anomalies. The steady-state events are even more frequently accompanied by transients, the discrimination and localization being far more exigent and requiring advanced signal separating tools to be incorporated in the measurement equipment. This paper shows the event detection performance of the spectral kurtosis as a signal separating tool in the frequency domain. The disturbances under test are hybrid signals resulting from the coupling between amplitude defects and non-desired higher frequencies. Being a fourth-order spectrum, the kurtosis is confirmed as a noise-resistant tool that enhances impulsiveness, therefore characterizing the electrical anomalies. In the beginning of the analysis, the voltage sag is established as a reference; then, the disturbances (oscillatory transients and harmonics) are coupled at the starting and ending instants of the sag, resulting in complex hybrid events. The results show that the spectral kurtosis enhances the detection of both types of events (steady state and transients), which are outlined in a bump shape in the fourth-order frequency pattern and centered in the main carrier frequency. Indeed, while the oscillatory transients are associated with softer and lower-amplitude peaks, the harmonics correspond to crisper and higher ones. As these mixed electrical faults are very common in the actual power grid, the article postulates the higher-order spectra to be implemented in prospective online measurement instruments.

Graphical Abstract

1. Introduction

With the uncontrolled growth of society, the consequent energy demand forces drastically adapting the power grid topology. Examples of this change are electric and hybrid cars, which suppose important loads in the power system, and nowadays, they are widely introduced [1,2].
In order to protect the environment, big fossil fuel power plants are being replaced by small renewable energies generators, which suppose enormous implications for the environment and for the power grid. A big synchronous power generator, as can be found in an average power plant, delivers a perfect sinusoidal waveform. In addition, any wave distortion is compensated and converted into a perfect waveform, and it can correct some grid outbalances by generating or consuming reactive power.
At the same time, the power grid is turning into a meshed system. This configuration allows the distribution system to prevent faults, creating new alternate routes for the power flow. In addition, using the information of power generated and consumed in any node and the grid characteristics, it can be configured to minimize losses. However, this has implications when a disturbance occurs. In a mesh system, many lines are connected to any node. In this situation, when a disturbance takes place at a point, it will leave all of the lines connected to that node, expanding quickly through the power grid.
Evolution of the topology of a high voltage power grid could be seen in [3]; the effects of different topologies in a normal function are explained in [4]; and an idea of the complexity level reached by the power gird is given in [5]. All of these complex network schemes require changes in the control, as seen in [6], which implies more elements in the grid. Any additional element could introduce disturbances in the case of a fault.
Collaterally, the new advantages in electrical power generation and power distribution systems have led to a power quality (PQ) degradation. Equipment, such as transformers, motors or conductors, undergo overload and generate low PQ conditions, due to harmonics, insufficient voltage, oscillatory transients or other disturbances. Due to the shared responsibility of PQ, it is important to monitor the voltage and currents at the connection point and in the low, medium and high voltage grid.
One of the most important PQ disturbances is the sag or voltage dip, explained and regulated in [7]. It consists of a reduction of the amplitude value with a duration from a half cycle to a few minutes. Changes in the amplitude usually occur suddenly and carry coupled disturbances: oscillatory transients and temporal harmonic distortions are common during or after amplitude transitions. All of those disturbances are examined and explained carefully by IEEE in [8]. Because the amplitude change in the base waveform is usually higher, coupled disturbances are difficult to detect.
In this frame, this work postulates the spectral kurtosis (SK) as a computing tool to detect and characterize transient disturbances that are coupled to steady-state sags. The capabilities of higher-order statistics (HOS) to enhance impulsiveness (non-Gaussian behavior) have been formerly proven in some fields and areas of science and technology [9,10,11,12,13,14,15,16]. More precisely, in the field of the PQ characterization, the research team has developed intelligent methods to detect power quality events, using HOS in the time and frequency domains [17,18,19].
Regarding direct backgrounds of HOS applications in this field, in the time-domain, several notable works are worthy, e.g., Bollen et al. introduced new statistical features to PQ event detection [20]. In the same direction, Gu and Bollen [21] found relevant characteristics associated with PQ events in the time and frequency domains. The work by Ribeiro et al. [22] is also remarkable, which extracted new time-domain features based on cumulants. The same authors performed the classification of single and multiples disturbances using HOS in the time domain and Bayes’ theory-based techniques [23]. HOS techniques and estimators have also been implemented to specifically detect sags and swells [24]. The categorization of PQ anomalies had been formerly performed by Nezih and Ece in [25], where they proved that HOS and quadratic classifiers improve the second order-based methods. The most aligned work is the one developed by Liu et al. in [26], where the researchers extracted a rich set of higher-order features related to the analysis as different types of PQ events using spectral kurtosis.
The paper is structured as follows. Firstly, the SK procedure is explained in Section 2, then coupled disturbances are examined using the SK analysis in the Section 3. Finally, conclusions are drawn in Section 4.

2. The Spectral Kurtosis: Procedure of Analysis

The SK is a fourth-order spectrum that deals with the peakedness of the probability density function associated with a measurement data series. From a practical point of view, the traditional power spectrum does not reflect these impulsive features of the time series, which indeed constitute the key attributes of PQ events. For that reason, the SK pattern should reflect the transient characteristics associated with the frequencies under test [12].
The procedure of computation is explained as follows. The time domain register is chopped into realizations or segments overlapped by 50%. This means that the second half of a time realization is the first half of the next realization. Then, the FFT of each realization is calculated and stored in a matrix, where every new FFT is added in a new row. With the data presented this way, Equation (1) represents the estimator involved in the SK computation:
G ^ 2 , X N , M m = M M - 1 M + 1 i = 1 M X N i m 4 i = 1 M X N i m 2 2 - 2
where X is the FFT matrix, in which the upper index i means the row (realization) and the index m is associated with the column (frequency). M realizations with N points have been considered. The final SK vector has the same points as the FFT vector; equal to N / 2 . The more realizations are involved in the calculation, the more precise is this procedure. However, the computational load increases with the size and the number of realizations, so a correct size and number of points must be indicated for the realizations to achieve suitable resolution and precision, but with an acceptable computational burden; in this case, 1000 realizations of 0.1 s each, using a sampling frequency of 20,000 Hz for all analysis. This means an SK vector of 1000 points and a maximum frequency of 10,000 Hz.
As said formerly, 1000 realizations of 0.1 s are needed; however, each register has only a 1-s length. This problem is overcome using a sliding version for the SK estimator. The spectral data are accumulated for 100 s before the SK could be calculated.
Hereinafter, the results are presented, dividing them into groups corresponding to different types of electrical perturbations.

3. Coupled Disturbances

In this section, an evolution to a coupled signal will be done. Firstly, the base response for a healthy signal will be set. Then, the sag will be examined alone, and finally, a sag with an oscillatory transient and harmonic distortion coupled is studied.
Firstly, and with the goal of getting a reference, the analysis of a healthy signal is performed (shown in Figure 1). The healthy condition is stated as a sinusoid of unitary amplitude. It will be the base amplitude assumed throughout the work.
With the frequency pattern characterization of this ideal signal, the anomalies associated with the disturbances are detected more easily. Synthetics are used in order to control all of the signal parameters, allowing in this way a correct comparison among them. In order to introduce real-life conditions, additive Gaussian noise with a standard deviation (σ) of 0.05 is added to all signals.
Figure 1. Spectral kurtosis (SK) analysis over a healthy signal. SK results corresponding to a 50-Hz sinusoid with additive Gaussian noise with (σ = 0.05), swept from 0 to 500 Hz.
Figure 1. Spectral kurtosis (SK) analysis over a healthy signal. SK results corresponding to a 50-Hz sinusoid with additive Gaussian noise with (σ = 0.05), swept from 0 to 500 Hz.
Energies 08 09777 g001
Signals have been synthesized with the same properties as described in the SK definition, i.e., 20,000 samples per second. That means a maximum frequency of 10,000 Hz for the SK analysis. However, the graph only shows up to 500 Hz. In the graph, the same response from 100 to 500 Hz can be seen. It continues equally for the rest of the frequencies. As no new information could be extracted from frequencies higher than 500 Hz, only frequencies under that value are shown, obtaining in this way a more detailed representation. Hereinafter, this is performed for all experiments, due to all frequencies involved in the considered disturbances being lower than 500 Hz.
In Figure 1, the same pattern is outlined for almost all of the frequencies. These are small variations around the zero, with a maximum deviation of 0.2. This value is related to the additive noise, which is rejected by the SK. In the lower frequency band of the graph, a −1-valued peak appears, centered at a frequency of 50 Hz. With a discrete spectral resolution of 10 Hz, the −1 value copes with the frequencies of 40 and 60 Hz. This zone is associated with the power signal frequency and indicates the feature associated with a constant amplitude signal.
Being a higher-order spectrum, the SK is postulated as a noise-resistant analysis tool. This property is tested studying different noise contamination levels. Table 1 shows the kurtosis values for the base frequency (50 Hz) and maximum variation around zero for the rest of the spectrum for Gaussian additive noise with σ from 0.05 to 0.80. If a 50-Hz kurtosis value is examining, no important change could be observed. As noise contamination increases, kurtosis takes values a litter higher, from −0.99998 to −0.9962. However, in this range, all could be considered −1, with an error lower than 1%. In relation to the maximum variation around zero, there are no relations with the noise contamination. Different values could be seen for different contamination levels, with a maximum value of 0.345210 and a minimum value of 0.206976. As a test of the randomness of these values, 1000 analyses with the σ = 0.05 condition have been performed. Each single analysis has returned a different maximum variation value around zero, the higher one being 0.417784 (only in two analyses) and the lower one 0.177371 (in some of them). The average value obtained in this experiment is 0.25.
Table 1. SK analysis results for a healthy signal with different noise contamination levels.
Table 1. SK analysis results for a healthy signal with different noise contamination levels.
Noise level (σ)50-Hz Peak valueMaximum variation around zero
0.05−0.9999860.250188
0.10−0.9999370.209871
0.15−0.9998660.251025
0.20−0.9997610.206976
0.25−0.9995990.234948
0.30−0.9994800.210391
0.35−0.9993110.259246
0.40−0.9990890.213317
0.50−0.9983680.263187
0.60−0.9978460.345210
0.70−0.9973000.216441
0.80−0.9962230.238456
The SK analysis associated with a non-defective state is examined. Now, the sag is introduced in the signal, and the differences in the SK pattern are analyzed, in order to establish the base analysis condition for isolated sags, so as to be able to study coupled disturbances. Signals are presented without Gaussian additive noise, with the objective of getting a cleaner representation. In Figure 2, the sag disturbance is represented, and Figure 3 depicts its associated SK analysis.
Figure 2. Sag. Time series of a sag with a depth of 30% and a duration of 10 cycles.
Figure 2. Sag. Time series of a sag with a depth of 30% and a duration of 10 cycles.
Energies 08 09777 g002
Figure 3. SK analysis over a sag. SK analysis shown from 0 to 500 Hz of a sag with Gaussian additive noise with σ = 0.05.
Figure 3. SK analysis over a sag. SK analysis shown from 0 to 500 Hz of a sag with Gaussian additive noise with σ = 0.05.
Energies 08 09777 g003
As said above, a sag is a reduction in the amplitude of the base waveform. In this experiment, the amplitude has been reduced by 30%, with crisp transitions.
In relation to the SK analysis, there are some changes. First, the vertical axis limits have been changed. The graph in Figure 1 shows the kurtosis values from −1.1 to 0.4. Now, it spreads from −15 to 190.
From 100 Hz onwards, the kurtosis response is constant, close to zero. If a zoom is made over this area, the same response is observed. However, with this vertical scale, it seems a constant value. The main change occurs in frequencies under 100 Hz. Two peaks appear, one of them before 50 Hz and the other one after this frequency component. This separation is occasioned by the same −1 peak at 50 Hz, seen in the previous analysis. As an important length is computed, even with the momentary change of amplitude, 50 Hz is still considered as a constant amplitude. That is represented with the −1 value at 50 Hz. Here, the - 1 kurtosis value could not be seen due to the vertical scale, but it shows exactly the same value as in the healthy signal. When the amplitude of a single frequency changes suddenly, the surrounding frequencies are affected, as well. That is the situation now. A kurtosis increase for frequencies around 50 Hz means an increase in the variability of these components. They do not reach important amplitudes, however, as they do not show any amplitude before; this means important amplitude variations.
Different sag conditions have been tested in the next experiment, keeping the noise contamination constant in σ = 0.05. Depths from 10 to 90% are settled, considering durations from 5 to 20 cycles. In each signal, kurtosis values for frequencies of 30, 50, 70, 100 and 150 Hz are taken. These frequencies are related to the first peak, the base frequency, the second peak, a near frequency in which the system is slightly altered and a far frequency in which the system is almost unaffected. All information related to those experiments is shown in the Table 2.
Table 2. SK analysis results for sag disturbance with different depths and lengths.
Table 2. SK analysis results for sag disturbance with different depths and lengths.
Sag deepSag lengthKurtosis
(%)(Cycles)30 Hz50 Hz70 Hz100 Hz150 Hz
1058794−0.99948420.0670.126
10109044-0.9992,643−0.0530.074
10208175−0.9992592−0.0470.084
20581,769−0.99933,7010.0510.127
201081,857−0.99935,095−0.0050.077
202077,130−0.99827,1540.1930.080
305173,911−0.99998,02910660.146
3010181,324−0.999102,50110990.088
3020178,341−0.99895,58311600.084
405256,598−0.999174,31126600.185
4010261,075−0.998172,65711270.242
4020260,819−0.996168,98522340.108
505319,816−0.999233,25061210.199
5010329,233−0.997237,40735340.208
5020323,410−0.995244,31052460.151
605361,178−0.998289,91811,3500.305
6010368,614−0.997288,90695360.237
6020366,564−0.994290,23711,3250.242
705392,615−0.998329,94516,6590.334
7010398,564−0.996326,26016,3001,121
7020396,867−0.993327,96217,3680.366
805415,655−0.998360,68729,4810.828
8010418,429−0.996357,97627,2130.890
8020418,173−0.992359,11426,5200.769
905431,010−0.998386,57338,4831488
9010433,368−0.996384,00242,9271236
9020433,365−0.992383,94140,5640.830
As indicated above, at 50 Hz, the peaks could be considered −1, even in the higher depth and length conditions. Values in this table are truncated, not rounded.
The durations of the disturbances have no effect on the kurtosis values, keeping the same value for each depth. In relation to the depth, a higher one involves a higher alteration of the waveform. SK detects frequency variability, showing higher values of kurtosis as sag depth increases.
In Table 2 and in Figure 3, it is seen that the 30-Hz frequency peak has a kurtosis value higher than the peak at 70 Hz. In a visual examination of all experiments, 100 Hz seems not to be altered by the sag disturbance. However, a detailed analysis reveals that the kurtosis value is increased when peaks reach important values. Nevertheless, the value associated with 150 Hz is hardly altered, even in the higher depth and length conditions.
Up to now, the base condition of a healthy signal and a sag disturbance alone has been intensely studied. Sag will be coupled to some disturbances. SK analysis will be studied in all of the situations, searching the characteristics of all disturbances introduced, with the goal of differentiating all of them. In the rest of the experiments, a depth of 30% and a length of 10 cycles have been selected.
Then, the oscillatory transient will be coupled to the sag. It consists of oscillations around the normal waveform value, with a crisp start and the exponential decay of its amplitude. The signal shown in Figure 4 is composed of the sag indicated before and the oscillatory transient with oscillations of 200 Hz, a duration of four cycles of the power waveform and an initial amplitude of 0.2. It has one oscillatory transient at the sag start and another at the sag end. As in previous graphs, this is shown without noise, in order to observe the disturbance better. In Figure 5, the SK analysis of the coupled signal is shown.
Figure 4. Sag with a coupled oscillatory transient. Time series of a sag with a depth of 30% and a duration of 10 cycles, with an oscillatory transient with a maximum amplitude of a 0.2, a duration of four cycles of the power signal and a frequency of 200 Hz.
Figure 4. Sag with a coupled oscillatory transient. Time series of a sag with a depth of 30% and a duration of 10 cycles, with an oscillatory transient with a maximum amplitude of a 0.2, a duration of four cycles of the power signal and a frequency of 200 Hz.
Energies 08 09777 g004
Figure 5. SK analysis over a sag with a coupled oscillatory transient. SK results shown from 0 to 500 Hz of a sag with an coupled oscillatory transient, with Gaussian additive noise with σ = 0.05.
Figure 5. SK analysis over a sag with a coupled oscillatory transient. SK results shown from 0 to 500 Hz of a sag with an coupled oscillatory transient, with Gaussian additive noise with σ = 0.05.
Energies 08 09777 g005
The main frequency introduced with the oscillatory transient is 200 Hz. As mentioned above, the sags do not alter the frequency over 150 Hz. At first glance, a new element could be seen. In the zone of the graph corresponding to lower frequencies, the same peaks as in the sag situation could be seen. Thus, the sags are detected the same as if this were the only disturbance. However, another peak appears around 200 Hz, which is the oscillatory frequency. This last peak is associated with the oscillatory transient.
As in the sag case, an oscillatory transient supposes an amplitude variation. However, on the contrary, it does not show a constant amplitude during all of the analyzed vectors, so it does not have a −1 value in the central frequency. The sudden start in an oscillatory transient creates an increment in the variability of oscillation frequency. In addition, the transient amplitude is continuously decreasing, which increases the variability of the oscillation frequency more. All of this affects the surrounding frequencies, increasing their variability, as well. The result of all of this is the peak shown in the graph, with a maximum value in the central frequency and a soft decreasing of the kurtosis value as the frequency moves away from 200 Hz, symmetrically in both directions.
Some oscillatory transient conditions are considered, keeping the same sag conditions as indicated before. The oscillatory transient amplitudes taken are from 0.15 to 0.3, and the lengths considered for each amplitude are from two to 10 power signal cycles. As symmetrical peak is returned; only the information of the half peak is taken, being similar on the other side. The kurtosis values for frequencies of 200, 210, 220 and 230 Hz have been taken, and the peak width is measured with the number of points around 200 Hz that have kurtosis values over three. This can be seen in Table 3.
Table 3. SK analysis results for a sag with an oscillatory transient coupled with different conditions.
Table 3. SK analysis results for a sag with an oscillatory transient coupled with different conditions.
TransientTransientKurtosis valuePeak
amplitudelength200 Hz210 Hz220 Hz230 Hzwidth
0.1524026365415190.4464
0.15428,87720,731648419586
0.15658,32140,09712,97924616
0.151080,05750,00010,98920506
0.20212,74310,8044,66518826
0.20462,29349,00712,88030777
0.20693,30471,55122,57634457
0.2010118,72591,33925,93942128
0.25232,12523,63910,084579411
0.25499,07282,61435,33410,22111
0.256140,299117,12045,02211,40412
0.2510150,465132,73856,85211,74913
0.30254,75539,57119,989949313
0.304148,846117,60853,01717,71713
0.306167,414152,66578,88723,71414
0.3010169,960166,17883,78223,32613
It is worth remarking about the following issues. Firstly, the length of two cycles reduces the kurtosis value in the oscillation frequency, for all amplitude situations. Simulation creates a sinusoidal wave with the desired length and applies an exponential decrease in the transient amplitude, which ends in zero. The shorter the transients were, the more the amplitude of the first cycle would be reduced. This reduces the initial change of the amplitude. Moreover, the variability is related to the change of amplitude during the transient, and as it has a lower extension, the variability is reduced, as well.
This trend continues as length increases. All kurtosis values are higher as the lengths of transients increase. This could be observed in any frequency for any amplitude considered in the table.
As indicated before, the frequencies selected in the table are related to the half peak. The alteration of SK is symmetrical, as can be seen in Figure 5. As the transient amplitude increases, the peak turns softer. This means a higher variation as the frequency moves away from the central one. For the 0.3 amplitude and a 10-cycle length, the kurtosis associated with 200 and 210 Hz is almost the same, and the next frequency shows half value of them. However, for the same length, taking an amplitude of 0.2, a change from 200 to 210 Hz supposes a loss of 23%, and in the next step, it losses 71%. Higher variations, as seen for the second one, mean a pointier peak. Lower kurtosis variations among near frequencies, as in the first one, create a softer peak, with a rounder shape.
In addition, the higher amplitude means a higher variability itself. At the starting moment, it creates a bigger step, and then, during the exponential decay, a bigger change is observed. Taking the same two situations, with an amplitude of 0.3 and a length of 10 cycles, the maximum value returned is 169.96. With the same length, and for an amplitude of 0.2, a peak with a maximum of 118.72 is originated.
The last column shows the width of the peaks, taking a threshold of three. An increase in the length multiplies the kurtosis value associated with each frequency. As frequencies near the end of each peak have very low values, multiplying them does not change the situation very much. This could be observed by the fact that almost the same length is indicated for each amplitude. On the other hand, as the transient amplitude increases, the kurtosis value is incremented, increasing at the same time the peak widths and turning them softer.
The last disturbance considered is the harmonic distortion at the sag start and end. The harmonic selected is the fourth one, which implies a frequency of 200 Hz. This one has been used in order to introduce the same frequency as in the previous test. In the oscillatory transient situation, oscillations starts with a high amplitude, and this decays in an exponential way. On the other hand, the harmonics take the same amplitude during all of the disturbed segment. In Figure 6, the signals with harmonic distortion are shown, and the SK analysis is shown in Figure 7 from 0 to 500 Hz.
The conditions for the sag are the same as in the previous experiments. The same frequency has been introduced; any difference observed in relation to the oscillatory transient are due to harmonic distortion. The oscillatory transient introduced in Figure 4 starts with a high alteration; however, the signal recovers its normal state very quickly. On the other hand, in Figure 6, changes in the signal are at the same level for all of the disturbed segment.
In Figure 7, the same response for frequencies under 100 Hz as in Figure 5 can be seen. This is related to the sag: as the same conditions have been considered, its response remains equal, as well. There is another peak around 200 Hz, in the same position as in the previous experiment. However, the peak’s shape is different. Now, it takes a shorter extension and reaches a higher kurtosis value, resulting in a pointier peak than in the previous situation.
Figure 6. Sag with harmonic distortion. Time series of a sag with a depth of 30% and a duration of 10 cycles with harmonic distortion with an amplitude of a 0.07 and a duration of four cycles of the power signal. The fourth harmonic has been selected (200 Hz).
Figure 6. Sag with harmonic distortion. Time series of a sag with a depth of 30% and a duration of 10 cycles with harmonic distortion with an amplitude of a 0.07 and a duration of four cycles of the power signal. The fourth harmonic has been selected (200 Hz).
Energies 08 09777 g006
Figure 7. Sag analysis with harmonic distortion. SK results shown from 0 to 500 Hz of a sag with harmonic distortion, with Gaussian additive noise with σ = 0.05.
Figure 7. Sag analysis with harmonic distortion. SK results shown from 0 to 500 Hz of a sag with harmonic distortion, with Gaussian additive noise with σ = 0.05.
Energies 08 09777 g007
As in the oscillatory transient experiment, some harmonic disturbance conditions will be tested, keeping the sag conditions and the harmonic distortion frequency equal. Harmonic distortion amplitudes have been taken from 0.07 to 0.20, using lengths from one to four cycles of the base waveform. The results of these experiments are shown in Table 4. As in the oscillatory transient situation, the kurtosis value related to each signal has been taken for frequencies of 200, 210, 220 and 230 Hz, and the peak widths have been measured as the number of points with a kurtosis value over three around 200 Hz.
In this experiment, lower extensions and amplitudes have been considered. Even with that, higher kurtosis values are obtained.
Examining Figure 6, the harmonic distortion peak looks pointier than the one related to the oscillatory transient, shown in Figure 4. Comparing the conditions shown in the figures, for the same length, four cycles of the base waveform, harmonic distortion with an amplitude 0.07 creates a three-point peak with a maximum value of 104, and an oscillatory transient, with an amplitude of 0.20, creates a seven-point peak with a maximum value of 62. The peak related to harmonic distortion shows a higher kurtosis value and a lower width than that related to the oscillatory transient.
If the same conditions are used for harmonic distortion and the oscillatory transient, for example an amplitude of 0.20 and a length of four cycles of the base waveform, different results are obtained. Harmonic distortion gives a peak with a maximum value of 261 and a width of nine, and the oscillatory transient returns a peak with a maximum value of 62 and a width of seven. Now, the peak related to harmonic distortion is wider; however, if the ratio between the maximum value and the width is shown, the harmonic distortion has 260/9 = 28.9 and the oscillatory transient 62/7 = 8.8. This means that the harmonic peak is taller in relation to its width than the one related to the oscillatory transient.
Table 4. SK analysis results for sag with harmonic distortion coupled with different conditions.
Table 4. SK analysis results for sag with harmonic distortion coupled with different conditions.
HarmonicsTransientKurtosis valuePeak
amplitudelength200 Hz210 Hz220 Hz230 Hzwidth
0.07110,048580430280.5495
0.07228,62418,55321120.1564
0.07358,13825,19620110.2004
0.074104,39534,22010690.2093
0.10132,18925,50511,37321506
0.10264,87344,32968940.5405
0.103116,76268,75666540.5715
0.104169,96472,61049260.6375
0.15184,98372,38641,95090149
0.152127,73489,24820,80419357
0.153174,153126,72923,23526746
0.154236,331129,79713,63413995
0.201149,340135,24084,17427,4548
0.202178,761133,98149,13471899
0.203207,419164,96345,99966928
0.204261,294166,90641,01650879
As indicated before, with the same conditions, the harmonic distortion generates kurtosis responses higher than the oscillatory transients. Taking the conditions of an amplitude of 0.20 and a duration of two cycles, the maximum kurtosis values are 12.7 for the oscillatory transient and 178.7 for the harmonic distortion. This difference is accentuated for a duration for four cycles and the same amplitude, reaching values of 62.3 for the oscillatory transient and 261.3 for the harmonic distortion. If the same maximum value is searched, an amplitude of 0.07 and a duration of three cycles must be used for harmonic distortion, obtaining a value of 58.13, and an amplitude of 0.15 with a duration of six cycles for the oscillatory transient, obtaining a value of 58.32.
As said before, the peaks associated with harmonic distortion are sharper than the ones related to the oscillatory transients, and now, this will be detected over the kurtosis evolution. First, the kurtosis change from the peak central frequency (200 Hz) to the next one (210 Hz) will be examined, for common conditions. These are amplitudes of 0.15 and 0.20 and durations of two and four base signal cycles. In these conditions, the change rates are from 9% to 30% for oscillatory transients and from 25% to 45% for harmonic distortion. If all simulations are examined, variation rates of 2% and 37% are obtained for oscillatory transients, and for harmonics distortion, these take values from 9% to 67%. Taking common conditions or taking all simulations, a higher variation could be observed for harmonic distortion. This higher change implies a crisper peak.
With the goal of considering a real-life situation, we have introduced Figure 8, where a sag is accompanied by a severe harmonic distortion (the same as in Figure 8) in a Gaussian additive noise floor with σ=0.5. In this case, the measurement bandwidth is wider, taking the whole profit of the data acquisition equipment, connected via a differential probe to the phase voltage A-B. The complete analysis is also shown, where the time domain variance (power variations), skewness (indicator of the waveform symmetry) and the kurtosis (peakedness of the time domain signal) are also presented. The peaks in the SK graph are the same as those in Figure 7. In fact, the amplitude transition is reflected as a spectral leakage around the harmonics. Indeed, the y-axis in the SK graphs accounts for variability, as the kurtosis is an estimator of the peakedness of the probability density function associated with the measurement time series.
Figure 8. SK analysis over a sag with harmonic distortion, with Gaussian additive noise with σ = 0.05 . If harmonics at 25 and 75 Hz are present and there is no sag, two “clean” spectral lines would appear. In the case of a hybrid event, the SK outlines the spectral leakage around the harmonic frequencies.
Figure 8. SK analysis over a sag with harmonic distortion, with Gaussian additive noise with σ = 0.05 . If harmonics at 25 and 75 Hz are present and there is no sag, two “clean” spectral lines would appear. In the case of a hybrid event, the SK outlines the spectral leakage around the harmonic frequencies.
Energies 08 09777 g008
Finally, and with the two-fold purpose of illustrating harmonic detection and noise rejection at the same time, two pure harmonics (at 50 and 75 Hz) have been contaminated by colored Gaussian noise (σ = 0.05, bandwidth: 250 to 350 Hz). The effect is the clear distinction of the two spectral lines along with the noise vanishing. In Figure 9, the spectral lines are outlined both in the second and in the fourth-order domain (SK); in the latter, noise is completely rejected, as expected by a higher-order spectrum.
Figure 9. Harmonics’ discrimination and noise rejection. Harmonics at 50 and 75 Hz have been mixed with colored Gaussian noise (σ = 0.05, bandwidth: 250 to 350 Hz).
Figure 9. Harmonics’ discrimination and noise rejection. Harmonics at 50 and 75 Hz have been mixed with colored Gaussian noise (σ = 0.05, bandwidth: 250 to 350 Hz).
Energies 08 09777 g009

4. Conclusions

This paper shows that the SK has good noise resistance and is a robust spectral analysis tool, returning a similar response with different noise contamination levels (uniform and Gaussian noise produce the same results). As a variability analysis tool, associating the peakedness of the probability density function, and in the frame of PQ analysis, different sag depths have returned different kurtosis maximum values, but always with the same shape, showing the SK normal value at 150 Hz. The sag length has no effect on the kurtosis response, showing similar values for each sag depth.
In the coupled disturbances study, sag conditions have been set constant in order to obtain similar test conditions. A frequency of 200 Hz has been introduced with oscillatory transients and harmonic distortion, returning a peak centered at the same position. Peaks related to these show a maximum width of 140 Hz, which means a range from 130 to 270 Hz. In the considered sag conditions, the maximum kurtosis value is 0.8 at 130 Hz, with an average kurtosis value of 0.05. That confirms that there is no interference from the sag SK response for the oscillatory transient and harmonic disturbance SK response.
Oscillatory transients show an exponential amplitude decay, so the alteration of the base wave is lower than in harmonic distortion, because this second situation keeps its amplitude constant. In addition, the oscillatory transient has one hard transition, at its start, because of its end amplitude being zero. On the other hand, harmonic disturbance have two hard transitions, one at its start and one at its end. All of these create higher values in harmonic disturbance situations, with even higher impact on the central frequency, for the same amplitude and length conditions. In addition, crisper peaks are obtained in harmonic distortion, due to a higher impact over the central frequency than in the oscillatory transient situation, in which adjacent frequencies are more affected. Taking all of this into consideration, when a symmetrical peak appears in the SK response, wide and soft ones are related to oscillatory transients and narrow and crisp ones to harmonic distortion.
As a final remark, the performance of the SK over hybrid (mixed) PQ events is promising, and it has been proven over a set of synthetics. The next step will focus the implementation on an autonomous measurement system, which requires adapting the measurement requirements to norms or the International Electrotechnical Commission (IEC). Compliance with the standards [27,28] that define the conditions for data acquisition in harmonic measurements for power supply systems implies the use of a basic acquisition window corresponding to a time period of Δt=200 ms. Sampling frequency and the rest of the parameters have to be adapted to these requirements.

Acknowledgments

This work is funded and supported by the Spanish Ministry of Economy and Competitiveness in the frame of the Statal Plan for Excellency in Scientific and Technical Research, via Research Project T E C 2013 - 47316 - C 3 - 2 - P (SCEMS-AD-TED-PQR). In addition, the authors would like to thank the Andalusian Government for funding Research Unit PAIDI-TIC-168.

Author Contributions

Juan José González de la Rosa and José María Sierra-Fernández have arranged the electronic instruments and acquired the data series of electrical perturbations. Agustín Agüera Pérez, José Carlos Palomares Salas and Álvaro Jiménez-Montero have analyzed the time series from the measurements (feature extraction from the SK graphs).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Esmaili, M.; Rajabi, M. Optimal charging of plug-in electric vehicles observing power grid constraints. IET Gener. Transm. Distrib. 2014, 8, 583–590. [Google Scholar] [CrossRef]
  2. Clement-Nyns, K.; Haesen, E.; Driesen, J. The impact of Charging plug-in hybrid electric vehicles on a residential distribution grid. IEEE Trans. Power Syst. 2010, 25, 371–380. [Google Scholar] [CrossRef] [Green Version]
  3. Buzna, L.; Issacharoff, L.; Helbing, D. The evolution of the topology of high-voltage electricity networks. Int. J. Crit. Infrastruct. 2009, 5, 72–85. [Google Scholar] [CrossRef]
  4. Rohden, M.; Sorge, A.; Witthaut, D.; Timme, M. Impact of network topology on synchrony of oscillatory power grids. Chaos 2014, 24, 1–19. [Google Scholar] [CrossRef] [PubMed]
  5. Pagani, G.; Aiello, M. Power grid complex network evolutions for the smart grid. Phys. A: Stat. Mech. Appl. 2014, 396, 248–266. [Google Scholar] [CrossRef]
  6. Blaabjerg, F.; Teodorescu, R.; Liserre, M.; Timbus, A. Overview of control and grid synchronization for distributed power generation systems. IEEE Trans. Ind. Electron. 2006, 53, 1398–1409. [Google Scholar] [CrossRef]
  7. International Electrotechnical Commission. IEC 61000-2-8-2002: Part 2-8: Environment—Voltage Dips and Short Interruptions on Public Electric Power Supply Systems with Statistical Measurement Results; International Electrotechnical Commission (IEC): Geneva, Switzerland, 2008. [Google Scholar]
  8. IEEE. IEEE Recommended Practice for Monitoring Electric Power Quality; IEEE Std 1159-2009 (Revision of IEEE Std 1159-1995); IEEE: Piscataway, NJ, USA, 2009; pp. c1–c81. [Google Scholar]
  9. Ribeiro, M.; Marques, C.; Duque, C.; Cerqueira, A.; Pereira, J. Detection of disturbances in voltage signals for power quality analysis using HOS. Eurasip J. Adv. Signal Process. 2007, 2007, 177–177. [Google Scholar] [CrossRef]
  10. Gonzalez de la Rosa, J.; Lloret, I.; Puntonet, C.; Piotrkowski, R.; Moreno, A. Higher-order spectra measurement techniques of termite emissions. A characterization framework. Measurement 2008, 41, 105–118. [Google Scholar]
  11. Mansour, A.; Mesleh, R.; Aggoune, E.H. Blind estimation of statistical properties of non-stationary random variables. Eurasip J. Adv. Signal Process. 2014, 2014, 21. [Google Scholar] [CrossRef]
  12. Agueera-Perez, A.; Palomares-Salas, J.C.; Gonzalez de la Rosa, J.J.; Sierra-Fernandez, J.M.; Ayora-Sedeno, D.; Moreno-Munoz, A. Characterization of electrical sags and swells using higher-order statistical estimators. Measurement 2011, 44, 1453–1460. [Google Scholar] [CrossRef]
  13. Gonzalez de la Rosa, J.J.; Agueera-Perez, A.; Palomares-Salas, J.C.; Sierra-Fernandez, J.M.; Moreno-Munoz, A. A novel virtual instrument for power quality surveillance based in higher-order statistics and case-based reasoning. Measurement 2012, 45, 1824–1835. [Google Scholar] [CrossRef]
  14. Gonzalez de la Rosa, J.J.; Agueera-Perez, A.; Palomares-Salas, J.C.; Moreno-Munoz, A. Higher-order statistics: Discussion and interpretation. Measurement 2013, 46, 2816–2827. [Google Scholar] [CrossRef]
  15. Yin, J.; Wang, W.; Man, Z.; Khoo, S. Statistical modeling of gear vibration signals and its application to detecting and diagnosing gear faults. Inf. Sci. 2014, 259, 295–303. [Google Scholar] [CrossRef]
  16. Miyazaki, R.; Saruwatari, H.; Nakamura, S.; Shikano, K.; Kondo, K.; Blanchette, J.; Bouchard, M. Musical-noise-free blind speech extraction integrating microphone array and iterative spectral subtraction. Signal Process. 2014, 102, 226–239. [Google Scholar] [CrossRef]
  17. Salas, J.C.P.; De la Rosa, J.J.G.; Sierra-Fernández, J.M.; Agüera-Pérez, A. HOS network-based classification of power quality events via regression algorithms. Series: Advanced signal processing techniques and telecommunications network infrastructures for Smart Grid analysis, monitoring and management. EURASIP J. Adv. Signal Process 2015, 49, 1–11. [Google Scholar]
  18. De la Rosa, J.J.G.; Sierra-Fernández, J.M.; Agüera-Pérez, A.; Salas, J.C.P.; Moreno-Muñoz, A. An application of the spectral kurtosis to characterize power quality events. Electr. Power Energy Syst. 2013, 49, 386–398. [Google Scholar] [CrossRef]
  19. De la Rosa, J.J.G.; Salas, J.C.P.; Agüera-Pérez, A.; Sierra-Fernández, J.M.; Moreno-Muñoz, A. A novel virtual instrument for power quality surveillance based in higher-order statistics and case-based reasoning. Measurement 2011, 44, 1453–1460. [Google Scholar] [CrossRef]
  20. Bollen, M.H.J.; Gu, I.Y.H.; Axelberg, P.G.V.; Styvaktakis, E. Classification of underlying causes of power quality disturbances: Deterministic versus statistical methods. EURASIP J. Adv. Signal Process. 2007, 2007, 1–17. [Google Scholar] [CrossRef]
  21. Gu, Y.H.; Bollen, M.H.J. Time-frequency and time-scale domain analysis of voltage disturbances. IEEE Trans. Power Deliv. 2000, 15, 1279–1283. [Google Scholar] [CrossRef]
  22. Ribeiro, M.V.; Marques, C.A.G.; Duque, C.A.; Cerqueira, A.S.; Pereira, J.L.R. Detection of disturbances in voltage signals for power quality analysis using HOS. EURASIP J. Adv. Signal Process. 2007, 2007, 1–13. [Google Scholar] [CrossRef]
  23. Ribeiro, M.V.; Pereira, J.L.R. Classification of single and multiple disturbances in electric signals. EURASIP J. Adv. Signal Process. 2007, 2007, 1–18. [Google Scholar] [CrossRef]
  24. Agüera-Pérez, A.; Palomares-Salas, J.C.; de la Rosa, J.J.G.; Sierra-Fernández, J.M.; Ayora-Sedeño, D.; Moreno-Muñoz, A. Characterization of electrical sags and swells using higher-order statistical estimators. Measurement 2011, 44, 1453–1460. [Google Scholar] [CrossRef]
  25. Gerek, O.N.; Ece, D.G. Power-quality event analysis using higher order cumulants and quadratic classifiers. IEEE Trans. Power Deliv. 2006, 21, 883–889. [Google Scholar] [CrossRef]
  26. Liu, Z.; Zhang, Q.; Han, Z.; Chen, G. A new classification method for transient power quality combining spectral kurtosis with neural network. Neurocomputing 2014, 125, 95–101. [Google Scholar] [CrossRef]
  27. International Electrotechnical Commission. IEC 61000-4-7: Testing and Measurement Techniques. General Guide on Harmonics and Interharmonics Measurements and Instrumentation, for Power Supply Systems and Equipment Connected Thereto; International Electrotechnical Commission (IEC): Geneva, Switzerland, 2002. [Google Scholar]
  28. International Electrotechnical Commission. IEC 61000-4-30: Testing and Measurement Techniques. Power Quality Measurement Methods; International Electrotechnical Commission (IEC): Geneva, Switzerland, 2008. [Google Scholar]

Share and Cite

MDPI and ACS Style

De la Rosa, J.J.G.; Sierra-Fernández, J.M.; Palomares-Salas, J.C.; Agüera-Pérez, A.; Montero, Á.J. An Application of Spectral Kurtosis to Separate Hybrid Power Quality Events. Energies 2015, 8, 9777-9793. https://doi.org/10.3390/en8099777

AMA Style

De la Rosa JJG, Sierra-Fernández JM, Palomares-Salas JC, Agüera-Pérez A, Montero ÁJ. An Application of Spectral Kurtosis to Separate Hybrid Power Quality Events. Energies. 2015; 8(9):9777-9793. https://doi.org/10.3390/en8099777

Chicago/Turabian Style

De la Rosa, Juan José González, José María Sierra-Fernández, José Carlos Palomares-Salas, Agustín Agüera-Pérez, and Álvaro Jiménez Montero. 2015. "An Application of Spectral Kurtosis to Separate Hybrid Power Quality Events" Energies 8, no. 9: 9777-9793. https://doi.org/10.3390/en8099777

APA Style

De la Rosa, J. J. G., Sierra-Fernández, J. M., Palomares-Salas, J. C., Agüera-Pérez, A., & Montero, Á. J. (2015). An Application of Spectral Kurtosis to Separate Hybrid Power Quality Events. Energies, 8(9), 9777-9793. https://doi.org/10.3390/en8099777

Article Metrics

Back to TopTop