电能的准确计量关系到供电、配电企业及用户等多方的经济利益,人们对电能计量准确度的要求越来越高[1],其中有功功率测量准确度直接关系到电能计量的准确度。随着电力技术和信息技术的发展,一方面,分布式发电备如风能、太阳能和非线性用电设备如汽车充电桩、电力机车等并入电网,导致电压电流信号产生畸变,电流大小呈现动态变化[2-3],使得输入的瞬时功率信号产生动态的波动;另一方面,受硬件条件的限制,非整周期采样不可避免的存在于电能表电压电流信号的采样中[4]。信号源的畸变以及动态波动和非整周期采样的限制影响了电能表有功功率的测量甚至电能计量的准确性,因此,对于电能表有功功率测量算法的研究一直是国内外学者研究的重点。
目前的功率测量算法主要是通过构造滤波器来实现的[5]。频域分析方面,文献[6-8]采用加窗,如矩形窗、Hamming窗、Blackman窗以及改进的窗函数进行快速傅立叶变换(FFT)分析,提高了频谱分析精度,所以更适用于频谱分析,但该方法会造成频谱泄露和衰减,同时针对畸变信号和动态信号,这些窗难以同时满足主瓣能量集中和旁瓣衰减快的要求,不适用于有功功率的测量。时域分析方面,准同步算法及其改进算法提高了测量的准确度,但缺点是硬件开销加大,不适用于电能表集成情况[9-10]。时域分析上点集合算法计算量小,易龙强等[11]对点集合算法进行了改进,提高了测量的准确度,但依然受到整周期采样的限制;俞集辉等[12]通过分析非整周期采样下的谐波功功率时域和频域计算方法,仿真了非整周期采样造成的畸变波形有功功率测量误差。
以上文献主要通过FFT算法分析研究了畸变信号时频特性、整周期采样下的功率、固定参量的电能测量算法及非整周期采样下的谐波功率测量误差,但对如何在非整周期采样下提高有功功率滤波效果和电能累计准确度,研究较为缺乏。
本文基于窗函数和电能表有功功率测量结构化模型构造了一种新的有功功率滤波器,仿真分析了非整周期采样下该滤波器的滤波效果,针对稳态和动态测试信号仿真计算了不同有功功率测量算法下电能累计的误差,以期为电能表的改进提供参考。
1 多项Sa函数自卷积窗构建 1.1 多项Sa函数多项Sa函数的一般形式是K个相继间隔为T的Sa函数之和,其时域表达式为
wSa(t)=K∑kαksinπT(t−kT)πT(t−kT) | (1) |
频谱函数可表示为
WSa(ω)={TK∑kαke−jωkT|ω|≤πT0|ω|>πT | (2) |
式中,αk为每项Sa函数的加权系数。从频谱函数可以看出,多项Sa函数的频率响应全部集中在[-π/T, π/T]区间内,当多项Sa函数的加权系数α0=1、其他加权系数为零时,则为单项Sa函数。
1.2 离散多项Sa函数窗将连续的多项Sa函数离散化,并用矩形窗截短,得到有限长的多项Sa函数(polynomial Sa,P-Sa)窗,形式如下
w(n)=K∑kαksinπT(n−kT)πT(n−kT) | (3) |
wP−Sa(n)=w(n)RN(n)={w(n),n=0,1,2,⋯,N−10,其他 | (4) |
式中,w(n)为离散的多项Sa函数,wP-Sa(n)为截短得到的窗函数,N是窗函数的长度,RN(n)是长度为N的矩形窗。由频域复卷积公式可以得到P-Sa窗函数的频域表达式
WP−Sa(ω)=12π∫π−πWSa(θ)WR(ω−θ)dθ=WSa(ω)∗WR(ω) | (5) |
式中,WSa(ω)为多项Sa函数的频域特性;WR(ω)为矩形窗的频域特性,且满足
根据卷积运算理论,由离散化的长度为N′的两项P-Sa窗卷积可以得到新窗的长度为2N′-1。在自卷积窗末尾补零可以得到长度为2N′的P-Sa自卷积窗为
w2P−Sa(n)=wP−Sa(n)∗wP−Sa(n) | (6) |
由卷积定理可得,wP-Sa2(n)窗的频域表达式为
W2P−Sa(ω)=[12π∫π−πWSa(θ)WR(ω−θ)dθ][12π∫π−πWSa(θ)WR(ω−θ)dθ]=[WSa(ω)∗WR(ω)]2 | (7) |
式(7)表明自卷积窗的窗函数在主瓣越窄的情况下,旁瓣衰减越快,卷积后的旁瓣衰减更优。
2 P-Sa自卷积窗频域特性分析经验表明,3项Sa函数窗在加权系数比为1:2:1时具有较好的旁瓣衰减特性,wP-Sa(n)的表达式为
wP−Sa(n)=1∑k=−1αksinπT(n−kT)πT(n−kT),n=0,1,2,⋯,N′−1,[α−1α0α1]=[0.250.50.25] | (8) |
式中,[α-1 α0 α1]为三项Sa函数窗的加权系数,N′为窗长度,且满足N′=2i;i为正整数;Sa函数的相继间隔时长T=N′/4。由N′=64的三项Sa函数窗构造P-Sa自卷积窗,构造后的窗长度N=128。对比P-Sa自卷积窗与矩形窗、Hamming以及Blackman窗的时域波形和频域响应如图 1所示。
![]() |
图 1 不同窗函数的时域波形和频域响应 Fig.1 Function waveforms and frequency responses of different window functions |
从图 1(b)中可以看出,矩形窗旁瓣衰减幅度最小,Hamming窗旁瓣衰减速度较慢,Blackman窗旁瓣衰减幅度较大,P-Sa自卷积窗的旁瓣衰减幅度和速度均为最大。
3 P-Sa滤波器有功功率测量算法 3.1 电能表有功功率测量结构化模型智能电能表的有功功率测量结构化模型如图 2所示。有功功率测量过程为:利用电压分压器和电流互感器(分流器)采集负载电压信号u(t)和电流信号i(t),将模/数(ADC)转化后采集到的离散电压通道和电流通道信号经数字乘法器相乘得到瞬时功率信号,再由低通滤波器(LPF)对瞬时功率信号滤波得到有功功率信号。有功功率信号作为电能表有功功率测量部分的输出,同时也是电能累计单元的输入,电能累计单元将有功功率信号进行数频转换后得到电能脉冲,再通过累计输出电能脉冲实现电能计量。
![]() |
图 2 智能电能表有功功率测量结构化模型 Fig.2 Structural model of active power measurements with a smart meter |
采集到的负载电压u(t)和负载电流i(t)信号的数学表达式为
u(t)=∞∑m=1Umsin(2πmf0t+φm) | (9) |
i(t)=v(t)∞∑m′=1Im′sin(2πm′f0t+fm′) | (10) |
式(9)、(10)中的电压、电流信号均包含基波和多次谐波;m、m′=1时表示基波信号参数,m、m′≠1时表示各次谐波信号参数;Um为基波电压和各次谐波电压的幅值;Im′为基波电流和各次谐波电流的幅值;f0为基波频率,取值50 Hz;φm为电压基波信号和各次谐波信号的初相位;v(t)为包络信号,v(t)=1时i(t)为稳态电流信号,v(t)为正弦包络、梯形包络[13-14]及OOK包络信号[3]时,i(t)为动态电流信号。其中正弦包络信号v1(t)表达式为
v1(t)=12+12mAsin(2πf1t) | (11) |
式中,f1为正弦调幅信号的频率;mA为调幅系数,由于v(t)≥0,故|mA|≤1。
符号函数sgn(t)可通过阶跃信号表示为
sgn(t)={1,t>0−1,t≤0 | (12) |
梯形包络信号v2(t)的表达式为
v2(t)=7∑i=0[αi+βi(t−ti)][1+sgn(t−ti)2][1−sgn(t−ti+1)2] | (13) |
式(13)中,sgn(t)为符号函数;αi, βi为不同时间段内线性函数的参数,满足0≤αi≤1;ti为包络幅度阶梯变化时间点。通过符号函数来截取不同时段内的线性函数αi+βi(t-ti),分段的线性函数组成了梯形包络函数v2(t)。文中后续仿真的梯形包络参数取值如表 1所示。
![]() |
下载CSV 表 1 梯形包络函数参数 Table 1 Trapezoidal envelope function parameters |
OOK包络信号v3(t)的表达式为
v3(t)={1,0≤t<Ton(电流信号通)0,Ton≤t<Ton+Toff(电流信号断) | (14) |
式中Ton为电流通的时长,Toff为电流断的时长。
负载信号经过AD转化采样后得到离散的电压和电流信号,表达式为
u(n)=M∑m=1Umsin(2πmf0n/fs+φm) | (15) |
i(n)=v(n)M∑m′=1Im′(n)sin(2πm′f0n/fs+ϕm′)=v(n)i′(n) | (16) |
式中,fs为采样频率,n为离散采样点数,M为采样得到的最高次谐波,ϕm′为电流基波信号和各次谐波信号的初相位。
3.2 瞬时功率测试信号建模瞬时功率测试信号pi(n)的表达式为
pi(n)=u(n)i′(n)v(n) | (17) |
式中,v(n)为包络信号,v(n)=1时pi(n)为稳态功率测试信号,v(n)为正弦包络、梯形包络及OOK包络信号时,pi(n)为动态功率测试信号;u(n)、i′(n)分别为稳态电压、稳态电流信号。本文仿真的稳态测试信号波形为IEC62052标准中规定的5种畸变波形和50 Hz工频正弦波形,5种畸变波形分别为平顶波、尖顶波、脉冲波、多个电流过零点波和多个电压过零点波。
3.3 P-Sa滤波器设计有功功率滤波器结构化模型如图 3所示,瞬时功率信号pi(n)经过低通滤波器输出po(n),po(n)即为有功功率。
![]() |
图 3 有功功率滤波结构化模型 Fig.3 Structural model of the active power filter |
若滤波器输入信号pi(n)的长度为n′,则经过低通滤波器后输出的有功功率信号po(n)为
po(n)=pi(n)h(N−1)+pi(n−1)h(N−2)+⋯+pi(n−N+1)h(0)=N−1∑k=0pi(n−k)h(N−1−k)po(n′+1)=pi(n′)h(N−2)+pi(n′−1)h(N−3)+⋯+pi(n′−N+2)h(0)po(n′+2)=pi(n′)h(N−3)+pi(n′−1)h(N−4)+⋯+pi(n′−N+3)h(0)⋮po(n′+N−1)=pi(n′)h(0) | (18) |
当输入信号pi(n)长度为n′时,输入信号的理论有功电能为
ei=∫T′0p(t)dt≈n′∑n=1pi(n)Ts | (19) |
式中,Ts为采样间隔时长,T′为输入瞬时功率信号的时长。经有功功率滤波器滤波之后,累加有功功率信号po(n)得到有功电能
eo=n′+N−1∑n=1po(n)Ts=[n′∑n=1(N−1∑k=0pi(n−k)h(N−1−k))+n′+N−1∑n=n′+1po(n)]Ts=[pi(1)h(N−1)+pi(2)h(N−1)+pi(1)h(N−2)+⋯+pi(N)h(N−1)+pi(N−1)h(N−2)+⋯+pi(1)h(0)+n′∑n=N+1(N−1∑k=0pi(n−k)h(N−1−k))+pi(n′)h(N−2)+pi(n′−1)h(N−3)+⋯+pi(n′−N)h(0)+⋯+pi(n′)h(0)]Ts | (20) |
由于滤波后对有功功率累计得到的电能与输入瞬时功率信号的理论有功电能相等,即ei=eo,故可得
n′∑n=1pi(n)Ts=n′∑n=1(N−1∑k=0pi(n−k)h(N−1−k))Ts+n′+N−1∑n=n′+1po(n)Ts=[h(N−1)n′∑n=1pi(n)+h(N−2)n′∑n=1pi(n)+⋯+h(0)n′∑n=1pi(n)]Ts=N−1∑k=0h(k)[n′∑n=1pi(n)Ts] | (21) |
由式(21)可知,滤波器系数需要满足
hP−Sa(n)=1∑n″ | (22) |
式(22)中,
3.2节中的稳态基波和5种畸变波形瞬时功率信号的归一化波形如图 4所示。仿真参数如表 2所示。
![]() |
图 4 基波及畸变波瞬时功率波形 Fig.4 Instantaneous power waveforms of fundamental and distortion waves |
![]() |
下载CSV 表 2 测试信号仿真条件 Table 2 Simulation conditions of test signals |
以平顶波为例,通过仿真得到动态瞬时功率测试信号归一化波形如图 5所示。
![]() |
图 5 平顶波动态瞬时功率测试信号波形 Fig.5 Quadriform wave dynamic instantaneous power test signal waveform |
图 4中的瞬时功率波形为稳态测试信号,图 5中的动态测试信号由正弦包络、梯形包络和OOK包络对稳态测试信号调制得到,稳态和动态测试信号作为后续误差仿真的输入信号。
4.2 有功功率滤波效果工频正弦波形及5种畸变波形的有功功率理论值表达式为
P = \sum\limits_{m = 1}^M {{U_m}{I_m}\cos {\theta _m}} | (23) |
式中,θm=φm-ϕm,非整周期采样时,瞬时功率经有功功率滤波器滤波后会产生纹波波动。P-Sa滤波器由长度为64的三项Sa函数窗自卷积后再乘以加权系数构造,构造的P-Sa滤波器长度N=128,响应时长为N/fs。对于时长大于2N/fs的瞬时功率信号,其有功功率Po(n)在N~2N点内的纹波波动εP大小可表示为
\begin{array}{l} {\varepsilon _{\rm{P}}} = \left\{ {\max \left[ {{P_{\rm{o}}}\left( N \right),{P_{\rm{o}}}\left( {N + 1} \right), \cdots ,{P_{\rm{o}}}\left( {N + N - } \right.} \right.} \right.\\ \left. {\left. 1 \right)} \right] - \min \left[ {{P_{\rm{o}}}\left( N \right),{P_{\rm{o}}}\left( {N + 1} \right), \cdots ,{P_{\rm{o}}}\left( {N + N - } \right.} \right.\\ \left. {\left. {\left. 1 \right)} \right]} \right\}/P \times 100\% \end{array} | (24) |
以平顶波为例,设不同滤波器仿真采样频率均为3 300 Hz,滤波器长度均为128,通过仿真得到其瞬时功率和经不同滤波器滤波后的有功功率波形如图 6所示。从图中可以发现,在非整周期采样下,矩形窗算法有功功率测量纹波波动最大,达到5.07%;Hamming窗和Blackman窗滤波后的纹波波动分别为0.42%和0.09%;而P-Sa自卷积滤波器滤波后波动最小,为0.02%,比矩形窗滤波效果提升252倍,比Hamming窗提升21倍,比Blackman窗提升4.5倍。
![]() |
图 6 平顶波瞬时功率及有功功率滤波形 Fig.6 Quadriform wave instantaneous power and active power waveform |
电能表通过累计有功功率来计量电能,电能累计误差表达式为
{\varepsilon _{\rm{e}}} = \frac{{{e_{\rm{o}}} - {e_{\rm{i}}}}}{{{e_{\rm{i}}}}} \times 100\% | (25) |
式中,ei为滤波前输入信号的理论有功电能,eo为采用不同的测试信号通过不同的有功功率算法滤波后的电能。累计误差如表 3~6所示。表 3为稳态测试信号,稳态仿真时长10.00 s;表 4为正弦包络动态测试信号,仿真时长2.10 s;表 5为梯形包络动态测试信号,仿真时长5.04 s;表 6为通断比20:20时的OOK动态测试信号,仿真时长2.00 s。
![]() |
下载CSV 表 3 稳态波形不同功率算法的电能计量误差 Table 3 Energy measurement errors for steady-state waveforms with different power algorithms |
![]() |
下载CSV 表 4 正弦包络不同功率算法的电能计量误差 Table 4 Energy measurement errors for sinusoidal envelopes with different power algorithms |
![]() |
下载CSV 表 5 梯形包络不同功率算法的电能计量误差 Table 5 Energy measurement errors for trapezoidal envelopes with different power algorithms |
![]() |
下载CSV 表 6 OOK包络不同功率算法的电能计量误差 Table 6 Energy measurement errors for OOK envelopes with different power algorithms |
从表 3~6可以发现,基于矩形窗、Hamming窗和Blackman窗3种算法测量的累计电能都小于实际电能,最大仿真误差达到-1.64×10-2,而本文算法的最大仿真误差不超过1.76×10-7,显著提高了有功电能计量的准确性。
4.2节和4.3节的仿真结果表明,P-Sa自卷积窗具有较好旁瓣衰减效果,有效抑制了非整周期采样下有功功率的纹波波动;同时该算法能有效测量稳态测试信号和动态测试信号的有功功率,电能测量更加准确。
5 结论提出一种P-Sa自卷积窗函数,并根据该窗构造了P-Sa滤波器,仿真结果表明:P-Sa自卷积窗具有很好的旁瓣衰减特性;P-Sa滤波器能有效抑制非整周期采样下有功功率的纹波波动;利用P-Sa滤波器测量有功功率能显著提高有功电能计量的准确性。本文研究结果为电能表非整周期采样下有功功率的测量以及电能的准确计量提供了改进的方向。
[1] |
CZARNECKI L S. Comments on active power flow and energy accounts in electrical systems with nonsinusoidal voltage and asymmetry[J]. IEEE Transactions on Power Delivery, 1996, 11(3): 1244-1250. DOI:10.1109/61.517478 |
[2] |
陆祖良, 王磊, 李敏. 对电能表动态测量功能评价的讨论[J]. 电测与仪表, 2010, 47(4): 1-4. LU Z L, WANG L, LI M. Discussion for evaluation of dynamic measurement function of electric energy meter[J]. Electric Measuring & Instrumentation, 2010, 47(4): 1-4. (in Chinese) DOI:10.3969/j.issn.1001-1390.2010.04.001 |
[3] |
WANG X W, CHEN J X, YUAN R M, et al. OOK power model based dynamic error testing for smart electricity meter[J]. Measurement Science & Technology, 2017, 28(2): 025015. |
[4] |
吴静, 赵伟. 非整周期采样下电功率测量的误差分析[J]. 清华大学学报(自然科学版), 2006, 46(7): 1197-1199. WU J, ZHAO W. Error analysis for electric power measurement with non-integer-period sampling[J]. Journal of Tsinghua University(Science and Technology), 2006, 46(7): 1197-1199. (in Chinese) DOI:10.3321/j.issn:1000-0054.2006.07.004 |
[5] |
肖勇, 赵伟, 罗睿希, 等. 数字化电能计量算法综述[J]. 电测与仪表, 2018, 55(7): 1-7. XIAO Y, ZHAO W, LUO R X, et al. Survey of digital energy measurement algorithms[J]. Electrical Measurement & Instrumentation, 2018, 55(7): 1-7. (in Chinese) DOI:10.3969/j.issn.1001-1390.2018.07.001 |
[6] |
WEN H, ZHANG J H, MENG Z, et al. Harmonic estimation using symmetrical interpolation FFT based on triangular self-convolution window[J]. IEEE Transactions on Industrial Informatics, 2017, 11(1): 16-26. |
[7] |
邱伟, 唐求, 唐璐, 等. 基于准同步序列重构的非稳态电力谐波分析[J]. 中国电机工程学报, 2018, 38(2): 456-464,676. QIU W, TANG Q, TANG L, et al. Power system harmonic analysis under non-stationary situations based on quasi-synchronous and sequence reconstruction[J]. Proceedings of the CSEE, 2018, 38(2): 456-464,676. (in Chinese) |
[8] |
高云鹏, 滕召胜, 卿柏元. 基于Kaiser窗双谱线插值FFT的谐波分析方法[J]. 仪器仪表学报, 2010, 31(2): 287-292. GAO Y P, TENG Z S, QING B Y. Harmonic analysis based on Kaiser window double spectrum line interpolation FFT[J]. Chinese Journal of Scientific Instrument, 2010, 31(2): 287-292. (in Chinese) |
[9] |
戴先中. 准同步采样及其在非正弦功率测量中的应用[J]. 仪器仪表学报, 1984(4): 55-61. DAI X Z. The quasisynchronous sampling and its application in the measurement of nonsinusoidal power[J]. Chinese Journal of Scientific Instrument, 1984(4): 55-61. (in Chinese) |
[10] |
王康, 滕召胜, 陈继斌, 等. 基于时域准同步的介损角快速测量算法[J]. 中国电机工程学报, 2015, 35(4): 873-881. WANG K, TENG Z S, CHEN J B, et al. A fast dielectric loss angle measurement algorithm based on the time-domain quasi-synchronous sampling technique[J]. Proceedings of the CSEE, 2015, 35(4): 873-881. (in Chinese) |
[11] |
易龙强, 宋慧娜, 陈金玲. 复化Newton-Cotes积分算法在电能计量中的应用[J]. 电测与仪表, 2010, 47(10): 15-20. YI L Q, SONG H N, CHEN J L. The application of compound Newton-Cotes integral algorithm to energy measurement[J]. Electrical Measurement & Instrumentation, 2010, 47(10): 15-20. (in Chinese) DOI:10.3969/j.issn.1001-1390.2010.10.004 |
[12] |
俞集辉, 刘金强, 郑可, 等. 非整周期采样谐波有功功率计量的时频域分析[J]. 计算机仿真, 2011, 28(1): 307-310. YU J H, LIU J Q, ZHENG K, et al. Analysis for harmonic active power measurement with non-integer period sampling in time-frequency domain[J]. Computer Simulation, 2011, 28(1): 307-310. (in Chinese) DOI:10.3969/j.issn.1006-9348.2011.01.078 |
[13] |
王学伟, 温丽丽, 袁瑞铭, 等. 智能电能表动态误差确定型测试激励信号的讨论[J]. 电测与仪表, 2014, 51(4): 1-7. WANG X W, WEN L L, YUAN R M, et al. Discussion of definitive test excitation signals for dynamic error of intelligent electric meter[J]. Electrical Measurement & Instrumentation, 2014, 51(4): 1-7. (in Chinese) DOI:10.3969/j.issn.1001-1390.2014.04.001 |
[14] |
张永旺, 朱孟, 王学伟, 等. 畸变波形动态测试信号模型及电能表动态误差分析[J]. 电测与仪表, 2018, 55(12): 92-99. ZHANG Y W, ZHU M, WANG X W, et al. Distortion waveform dynamic test signal model and dynamic error analysis of smart meter[J]. Electrical Measurement & Instrument, 2018, 55(12): 92-99. (in Chinese) DOI:10.3969/j.issn.1001-1390.2018.12.015 |