高精度FFT用于测量电力系统间谐波以及电压波动与闪变研究.pdfVIP

高精度FFT用于测量电力系统间谐波以及电压波动与闪变 张伏生 刘学 0049) (西安交通大学电力工程系，西安，71 1．引言 电力系统中谐波问题已被广泛的关注，通常的谐波一般指频率为工频整数倍的谐波成分，现行国家公 用电网谐波标准也只对这类谐波规定了限值和测试方法，而对于间谐波，目前国内缺乏测量手段，相关 的文献资料也很少，实际上间谐波及其影响广泛的存在于电力系统中。 问谐波是指非整数倍基波频率的谐波，这类谐波可以是离散频谱或连续频谱。根据傅立叶分解理论，周期 性的非正弦量只能分解出整数次的谐波。实际上许多非线性负载是波动的，其电流的幅值、相位或波形是 变化的。例如工业电弧炉、晶闸管整流供电的轧机是快速变化的冲击负荷，其电气量(电压或电流)的变化在 几毫秒或几十毫秒内就能观察到。在这种情况下，对于工频，“周期性”的前提已不存在，因而用傅立叶理论 分析的结果不符合或不完全符合实际。 目前工业上普遍采用由电容、电感和电阻构成的无源滤波器来吸收谐波和提高用电功率因数，这些滤波 器只能做成对某几次主要谐波(通常为整数次谐波，例如3,5，7，I1次)有明显的滤波效果，并没有考虑问谐波的影 响；而且，为了避免系统在某个频次下发生喈振，往往将滤波器的调谐频率偏移一定的范围，如果系统 在这个偏移的调谐频率点附近有较大的间谐波成分，该滤波装置将会通过较大的间谐波电流，给装置本身 以及系统带来较大的危害。所以，有必要对问谐波进行深入分析，从而对治理电网谐波污染提供更科学 的决策。 2．FFT计算谐波介绍 电力系统的谐波分析，通常都是通过快速傅立叶变换(rFr)实现的。利用FFvr算法可以实现整数次谐波的 检测。但是对于非整数次谐波，利用FFT算法得到的频谱分布却比较复杂： 用FFT算法，可以得到对应【o=21Tn／T，(n=0，l，…N一1)的离散频谱(其中T为采样周期)。 设频率为∞．．=2订n／T的信号，(七为0～N一1的任意数)，利用FFT算法得到的频谱分布为： Xo(n)=A。N (1) siIn聆(一(n．If-)zk万)rc如如拙(爿。复振幅) 若七为整数，则： xo(n)=华：硪 (2) 。X如)=0；n*k 可以看出，如果为整数，即被检测谐波为整数次谐波时，利用FFT算法得到的频谱分布为一条谱线， 利用这条谱线相关的参数就可以求出整数次谐波的频率、幅值和相角。因此利用FFT算法可以实现整数次 谐波的检测。但是对于非整数次谐波，即不是整数时，利用FFT算法得到的频谱分布比较复杂而且不能 计算出非整数次频谱。 F丌存在栅栏效应和泄漏现象，得出的信号幅值和相位不准，尤其是相位误差很大，无法满足准确的谐波测 量要求。本系统采用高精度FFT算法对FFT进行了改进，具体算法可参考文献[1】．实现了谐波(包括 整数次以及非整数次)的准确测量。 3．采用改进牛顿法提高计算速度 高精度FFT算法能精确测量出电力系统中的谐波和问谐波。为了提高计算速度，需要采取一些计算技 巧。 179 经统计，高精度FFT的计算量主要耗费在FFT上。为了尽量少用到FFT，我们将采样得到的电压电 流序列同时进行处理，即将电压信号作为复序列的实部，电流信号作为复序列的虚部，在得到计算结果 后再利用虚实奇偶对称公式将需要的结果分离出来心1。 (这里给出N等于1024时的公式)： · ∥ril．咖，置g=(a『i1．抽魍g-a[1024-i1．imag)／2； U[il．real=(a『i1．real+a[1024-il．real)12； (3) ／i11．imag=(a[1024一il．reM-ari]．real)12； ，fil．real=(a『i1．imag+a[1024-i]．imag)／2； 这样可以直接对二通道采样数据进行处理，节省了计算量。 高精度FFT算法巾在求解高阶方程时，方法有很多。在计算机中求解较为容易，可以调用现成的 函数，但