煤矿中一种实用的电网谐波监测方法

刘 霞

[文章摘要] 本文介绍了一种实用的电网谐波监测方法，这种方法同时采用了锁相环和快速抗混叠傅立叶变换。采用同步锁相环控制，避免了频谱泄漏的产生，使用FAFT算法分析采样数据，不仅速度快精度高，而且抗混叠。

[关 键 词] 谐波监测 锁相环 频谱泄漏 FAFT 抗混叠 1. 引言

在煤矿的电力系统中，谐波对各种煤矿电力设备、通讯设备及线路都会产生影响，严重时造成设备损坏和煤矿电力系统事故。为了保证电能质量，限制煤矿电网谐波，需要对电网中的谐波进行实时监测。

2.谐波测量方法简述

由于谐波具有固有的非线性、随机性、分布性、非平稳性和影响因素的复杂性等，因此难以对谐波进行准确测量，为此，许多学者对谐波测量问题进行广泛研究。现有谐波检测法按原理可分为：模拟滤波器、基于Fryze传统功率定义的谐波检测法、基于瞬时无功功率理论的谐波检测法、基于傅立叶变换的谐波检测法、基于神经网络的谐波检测法、自适应谐波检测法和基于小波分析的谐波检测法。 其中基于傅立叶变换的谐波测量是当今应用最多也是最广泛的一种方法，它由离散傅立叶变换过渡到快速傅立叶变换的基本原理构成，使用此方法测量谐波，精度较高、功能较多、使用方便。其缺点是需要一定时间的电流值，且需进行两次变换，计算量大、计算时间长，且这种方法未能解决极其重要的频谱泄漏、频谱混叠等重要问题。 3. 问题的解决

3.1频谱泄漏原因及解决方案的讨论

所谓频谱泄漏就是信号谱中各谱线之间相互影响，使测量值偏离实际值，同时在真实谱线两侧其它基波整数倍频率点上出现一些幅值较小的假谱。简单的讲，造成频谱泄漏的原因是因为采样频

率与信号频率的不同步，造成周期采样信号的相位在始端和终端不连续。

FFT变换存在着局限性，采集的点数必须是整周波或者整周波的倍数，因为FFT从整体上看是在信号整周波求信号加权平均值，如果由于采集频率与信号频率不同步，即没有进行整周波采集，势必会造成采集值的积分与实际整周波积分值之间存在偏差，使测量结果偏离实际值，同时在真实谱线两侧其它基波整数倍频率点上出现一些幅值较小的假谱，产生所谓的频谱泄漏问题。 减少频谱泄漏问题的方法有三种：

1）利用加窗插值算法对快速傅立叶算法进行修正，可选用不同的窗函数（如矩形窗、海宁窗、布莱克曼窗、布莱克曼——哈里斯窗）的插值算法，在实际测量过程中，选用矩形窗和海宁窗插值算法可满足测量精度的要求，但实时性较差，要采集多周期的数据，精度随着采样点及采样频率的增加而增加。

2）修正理想采样频率法，这种方法的主要思想是对每个采样点进行修正，得到理想采样频率下的采样值。该方法的计算量不大，并不需要添加任何硬件，实时性比上一种好，适合在线测量，但只能减少50％的泄漏。

3）利用锁相环(PLL)使信号频率和采样频率同步，锁相环把取自采样系统的电压信号的频率和相位与锁相环输出的同步反馈信号进行相位比较，当失步时，锁相环的相位比较器输出与二者相位差和频率差相关的电压，经滤波后控制并改变压控振荡

器的频率，直到输入频率与反馈信号频率同步为止，一旦锁定，便将跟踪输入信号频率变化，保持二者的频率同步，输出的同步信号去控制对信号的采样和加窗函数，此种方法实时性较好，同时可以解决频谱泄漏难题。 3.2 频谱混叠解决方案的讨论

由采样定理可知，信号的采样频率fs只有在2kjt1TTdt中的积分计算将复频谱系Akx(t)eT0分段进行。设Δ为采样时间间隔，则周期TN将

积分区间[0,T]以步长h2划分成N/2连续的长度相同的区间段，对其中任一区间段[tn1,tn1],t=1,3,5,…,N-1(tnn),记其内的信号

x(t)为xn(t)，则区间段内的积分计算可写成：

Akntn1xn(t)ej2ktTdt （2）

不小于信号的截至频率fc的两倍的条件下，即

fs2fc时，才不会产生频谱混叠的现象。

1）为了减少频谱混叠误差，可考虑采取增加周期内采样点数的办法，但这同样会大大增加计算量和硬件的负担。在实际的谐波测量中，通常采用在信号的输入通道上加前置低通抗混叠滤波器的办法对信号进行预处理。但这种做法在消除了频谱混叠误差的同时，又引入了另一方面的误差，原因是实际的低通滤波器的幅频特性不可能是理想的矩形，其相频特性也不可能是理想的直线，信号通过后，虽然滤去了导致频谱混叠的高频分量，但同时也使信号中有待进行分析计算的低频分量的幅值和相位产生了失真，从而给频谱的分析计算结果带来了相应的误差，所以加前置抗混叠滤波器也并非是抗频谱混叠误差的良策。

2）我国学者陈祥训提出了一种新的快速高精度的频谱分析方法FAFT测量谐波，可省去抗滤波器及有关部分，还可降低采样速率及对A/D精度的要求，使硬件成本降低，结构简化。 3.3 FAFT原理

电压、信号x(t)在区间[0,T]上可用傅立叶级数表示为：xtj2kTt

kAke2k则A1TjkTtT0x(t)edtAkejk

AkcoskjAksin k (1) 其中，Ak为k次谐波幅值，－k为其相位。

tn1则Ak就等于这N/2个区间段的积分的和再除以周期T，即：A1kTAkn （3） n将xn(t)在区间[tn1,tn1]内、在ttn处展开为泰勒级数，因为采样间隔Δ很小，因而可取泰勒展开式的前三项来近似地表示xn(t)，即：

xnta0a1(ttn)a2(ttn)2 （4） 式中，a0,a1,a2分别为泰勒展开式前三项的系数。 令

aj2kT，

可

得： Atn1attn1tattn1atkntan_10edttan11(tn)edtta2(ttn)2edt n1（5） 当k0时，有

A{8k22aatnsin2k222at0Ten2kknN8ka2esinN

8ka2k3at2k2T2eatncosN4a2TensinNj8k22aat2k2at1TencosNn2k34kaT1esinN}/(2k)

(6) 将上式中泰勒系数a0,a1,a2用该区间段内xn(t)的离散采样序列值xn(tn1)、xn(tn)、xn(tn1)的

函数来近似表示，则有:

a0xn(tn);

'xn(tn)xn(tn1)xn(tn1)a1;

1!2这样综合以上两种情况，FAFT复频谱系数实用计

算公式可写成以下形式：

Ak2ukAevenkRkAoddkWNkZkF （11）

''xn(tn)xn(tn1)2xn(tn)xn(tn1)a2 (7)

式中，Ak22!2即将x'''n(tn)和xn(tn)在该区间内用一阶差商和二

阶中心差商来近似表示，将TN，

2aj2kT ，WjNeN，经计算整理得：

NNA212rk21k2ukx2rWNRkx2r1W2rkNZkF (8) r0r0式中： ZZ'kkTukjvk;Fx0xN; RkR'kTNR22cos2kNN2ksin2kN;uku'kTN4kN8k223cosNksin4kN;v'1N2kvkT2k12k1cos4kN4kN4ksinN;xrr0,1,2,...N,为信号xn(t)的离散采样序列

的第r点的采样值，即xrxrTN

以上即为k0时的FAFT复频谱系数的实用计算

公式。

下面讨论k0时的情形。 令k0，则有

Atn1tn1tn1knta0dta1(ttn)dttan12(ttn)2dt(9) n_1tn1经计算整理可得： N21N21A02u0x2rR0x2r1Z0F (10) r0r0其中 u013N , R403N , v00 ， Zu100jv03N

even为信号偶次序列的离散傅立叶变N21换，即A2evenkx2rWrkN

r0Aoddk为信号奇次序列的离散傅立叶变换，即

NA 21xW2rkoddk2r1N；

r 0 u1403N , R03N ,v00 ;

Z1 0u0jv03N , Fx0xN; R2kN2kNk22cosN2ksinkNuNk8k223cos4kNN4kksinNv1N2k2k11cos4kN4k2kN4ksinNZkukjvk,k0

由于Akevenk,Aoddk和WN的周期递推性质

AmNevenk2Aevenk

AkmNodd2Aoddk

k1 W

kmNWk , Wm2NNNWkNN 可得FAFT复频谱系数计算公式的递推性质如下： A2ukRkkmNkmNAevenkmNAoddkWNZkmNFAkm1N2u2km1

2NAevenkRkm2NAoddkWk1NZkm1 (12) 2NF

从式（11）和（12）不难看出，由于只与k和N有关的常系数uk,Rk和Zk没有周期性重复的递推性质，因而用FAFT复频谱计算公式及其递推公式算出的复频谱系数不会出现采用FFT算法时发生的周期性重复现象，从而能够较好地抗混叠。而且利用信号的N+1个采样点就可以分析计算出N/2次以上的复频谱系数。

3.4 锁相倍频电路与FAFT算法的结合

根据国家标准关于谐波测量的规定，需测量谐波次数范围为从第2～63次。这样，根据FAFT频谱算法，每个电网周波采样64点即可满足测量要求。因此，可采用锁相倍频电路将输入信号进行N倍频，并与输入信号严格同步，产生N倍频同步触发信号，控制采样及保持电路进行A/D转换。可以解决频谱泄漏难题。 参考文献：

图1 锁相倍频电路

锁相倍频电路如图1所示。首先对电网信号uit进行取样，然后对它做整形处理（过零比较器），获得与基波信号频率一致的方波信号，将它进行锁相倍频，获得输出频率为64倍基波频率的整周期同步采样脉冲信号，用它去控制采样保持电路。 根据以上分析可知：

对于电压信号：AkUkejkUkcoskjUksink 对于电流信号：AkIkejkIkcoskjIksink 根据FAFT可得出煤矿电力系统电压各次谐波幅值 Uk、相角-k和电流各次谐波幅值Ik、相角-k。4. 结论

本文采用FAFT分析方法和锁相倍频技术检测煤矿电力系统谐波，不仅可以快速地利用信号的N+1个采样点就分析计算出N/2次以上的谐波幅值和相位，还可以达到避免频谱泄漏和抗混叠的作

uit整形电 路 锁相环 4046 计数器 64分频 采样保持用。

1.李圣清，朱英浩. 电网谐波检测方法的综述. 高电压技术. 2004, 30(3)：39－42

2.张伏生. 电力系统谐波分析的高精度FFT算法. 中国电机工程学报. 1999, 19(3)：62－65 3.潘文. 基于加窗插值FFT的电力谐波测量理论窗函数的研究. 仪器仪表学报. 1994，1：50－54 4.江道灼,马进. 锁相环在电力系统现场测控装置中的应用[J],继电器，2000，28(8):43－45