Pages

Thursday, April 4, 2013

POWER SYSTEM HARMONICS IN DETAIL: Subject 2: Effect of FFT Window Size on Harmonic Spectrum

     Power system harmonics are integral parts of Power Quality concerns. Harmonics causes lots of harms on fragile loads in the system and they cause overloading of the power system elements such as transformers or cables. Therefore, they should be eliminated, however, in order to eliminate them totally ,at first we should correctly measure their magnitudes by our control system.
     
Our aim : Measuring the whole harmonic spectrum of a signal of any Power System variable.
Tools we have: We have measuring probes, interface electronic card and a microcontroller based control system. 

Methodology: At first we should determine frequency resolution for our harmonic spectrum. Generally 10 Hz resolution is enough for ordinary power system load currents and voltages as stated in the IEC 61000-4-7 recommendation for 50 Hz synchronous frequency. For 10 Hz resolution, the window size will be 1/10 = 0.1 sec. After knowing the window size, we can sample all the data included by the window and calculate the FFT(Fast Fourier Transform) for finding the harmonic spectrum.

     Effect of  FFT Window Size on Harmonic Spectrum

     Up to now it is seem to be well defined, however, this is not the actual case. If we have integer number of complete cycles for a harmonic in a measuring window, the result of FFT gives us the correct magnitudes. However, if the frequency of harmonic component is not the integer multiplication of fundemental frequency, then there won't be complete cycles in an FFT window. Therefore, the magnitude of the actual signal is shared by the nearby frequencies and most importantly high frequency components,which does not exist in reality, will appear at the harmonic spectrum. Let me try to explain what is the reason of this situation. FFT is a linear transform that transforms the time domain signals into the frequency domain. Actually for FFT calculations we never use the whole signal (it can not be practical since the signal is from -infinity to +infinity), instead of that we use a part of it and assume that the signal continues periodically as if it continues from the beginning of the window. Figure 1 illustrates what I try to explain.
Figure 1: Case when window size contains complete cycles of periodic signal

     As can be seen in figure 1 when we have complete cycles of periodic signal, the FFT behaves correctly to the signal. and if we have a 250 Hz signal at window size 100 msec, corresponding harmonic spectrum is as in figure 2.
Figure 2: Harmonic Spectrum of 250 Hz signal in 100 msec window(10 Hz resolution)

     As expected, 250 Hz component is appeared at harmonic order of 25. This is because we have 10 Hz resolution. Please draw your attention to that the signal magnitude is calculated correctly as 1 and there will be no high frequency components, their magnitudes are too small and they can be ignored. As a result, our FFT code calculates the magnitude of our signal which is sampled at 25 kHz correctly as 1, since the window contains 25 complete cycles. 
     Let me pass through the most impressive result. What if our frequency is slightly deviates to 255 Hz. Actually please consider that 255 Hz component is not the harmonic component of an ordinary load, this value can be an individual harmonic in addition to the characteristic harmonics. In other words, 255 Hz can not be generally appear as a characteristic 5th harmonic, this can only be happen if the system frequency  deviates to 51 Hz, which is an unlike situation 50 Hz deviates +-0.2Hz at max. We ,for this article, are dealing with the individual load harmonic of 255 Hz. This kind of harmonic components (inter-harmonics) can be encountered mostly in the loads of metal industry, especially for metal melting Induction Furnaces. If you want to reach this kind of load characteristics please contact me. Since we have 255 Hz signal, our FFT window contains 25 complete cycles and 1 half cycle, which means that our control system calculates FFT of a signal of 255 Hz which is distorted at every 100 msec. I try to explain this case in figure 3.
Figure 3: Case when window size does not contain complete cycles of periodic signal

     Figure 3 shows an exaggerated version of the case I try to mention but this is a good illustration. Actually in our case, the distortion happens once for every 25 cycle. Let us see what happens to frequency spectrum for this case.

Figure 4: Harmonic Spectrum of 255 Hz signal in 100 msec window(10 Hz resolution)

     Figure 4 shows the frequency spectrum when we have 255 Hz signal in 100 msec window. I can hear your criticisms like it is normal that the magnitude of the signal is reduced since our resolution is 10 Hz, thus the magnitude is shared by the nearby points. The criticism is totally correct, because in the IEC recommendation the harmonic components are calculated with geometric sum of nearby components for this kind of cases. However, the point I want to draw attention is that there are higher frequency components appearing in the harmonic spectrum which are actually does not exist. When you thoroughly look at the figure 4, you can see that the magnitude of 400 Hz component is about 3% of the individual harmonic, which is an important result. Unlike the harmonic spectrum in figure 2, this case creates the high frequency harmonic components which do not exist in reality. This situation is totally the result of wrong measurement and may cause limit exceeding for higher order harmonics, whose magnitude in reality do not exceed the limits.
     I tried to explain the effect of FFT window size on Harmonic Spectrum of a signal. Actually, the results I show are inevitable for inter-harmonics, since we can not form a window which contains the complete cycles of all harmonic components. However, if we increase our resolution, meaning size of the FFT window, we can significatly reduce the effect of this measurement defects. This time we should consider the storage and execution time concerns. Our control system can not calculate the FFT if our data is huge ,actually we can not even store them, and our process time increases significantly.


No comments:

Post a Comment