摘要
针对传统同步相量算法受电网间谐波影响较大的问题,本文提出了一种基于新型改进矩阵束的电网同步相量测量算法.将电网信号构建成Hankel矩阵,对其进行奇异值分解并使用自适应定阶方法滤除其中的噪声干扰分量,通过矩阵运算得到准确的电压幅值、频率及相角.经过仿真验证,本算法对含有间谐波分量的电网信号测量结果优于流行的傅里叶变换测量算法.并根据此算法设计了一款便携式电网同步相量测量装置,可解决传统同步相量测量设备体积庞大、成本较高的问题,利用高精度卫星同步信号监测晶振实际频率,动态调整采样控制参数降低采样时间误差,将同步采样的电网波形参数与时间标签等无线传输至移动终端,完成电网信号幅值、频率及相角计算.实际测试表明,本装置具有较高的测量精度,其中电压幅值、频率和相角测量误差分别为0.038 5%、0.000 51 Hz和0.053 7°,满足《电力系统同步相量测量装置检测规范》(GB/T 26862—2011),具有较强的实际应用价值.
非线性、冲击性负荷和电力电子设备的不断增加,导致电网扰动问题越来越严重,直接影响公用电网可靠性和安全经济运
早期的同步相量测量主要使用过零点检测方法,该方法对硬件装置要求较低且易于实现,但受电网中的谐波及噪声影响,测量误差较
基于DFT的电网相量测量算法具有良好的抗谐波能力和较高的计算效率,但由于非整周期采样下的频谱泄露影响,此类算法在间谐波或带外干扰存在时性能较
在完成相量测量后,相应的硬件装置需要被部署以满足实际的工程需求.目前广泛使用的基于DFT的同步相量测量单元 (Phasor Measurement Unit, PMU
为解决上述问题,本文首先提出一种基于新型改进矩阵束的电网相量测量方法,将电网信号构建成Hankel矩阵,并进行奇异值分解,使用自适应定阶方法,滤除信号中的干扰信号,可显著降低间谐波对电网同步相量测量的影响,从而得到准确电网参数.对于硬件测量装置的设计,本文采用具有可移植性高、性能强等优点的Android移动终端实现电网相量测量中的数据处理,配合以高性能微控制单元(Microcontroller Unit, MCU)STM32F407为核心的处理电路实现电压数据的采集与发送,可节约硬件成本,减小装置体积,提升数据处理速度.
此外,为解决晶振频率可能出现偏
1 基于新型改进矩阵束的相量测量算法
1.1 同步相量定义
同步相量即以时间信号作为采样过程的基准,通过对采样数据计算而得到的相量数据.因此,电力系统电气信号模拟量的相量之间存在确定的相位关系.设模拟信号为:
(1) |
式中:Am为幅值;ω为角频率;θ为相角.
(2) |
如

图1 波形信号与同步相量的转化关系
Fig.1 Convention for synchrophasor representation
1.2 新型改进矩阵束同步相量测量算法
复杂电网信号可以看作由M个余弦分量和噪声信号组成:
(3) |
式中:M为余弦分量的个数;Ai为幅值;αi为衰减因子;ƒi为频率;θi为相位;n(t)为噪声信号.对
(4) |
式中:0 ≤ k ≤ N-1,N为信号y(t)的采样点数;代表采样信号的留数信息;zi=表示极点信息,其中为角频率.
通过采样数据y(kΔt)构造出Hankel矩阵Y:
(5) |
式中:L为矩阵束参数,一般选取N/4≤ L≤ N/3.
对矩阵Y进行奇异值分解:
(6) |
式中:U为(N-L)×(N-L)的正交矩阵;V为(L+1)×(L+1)正交矩阵;S为(N-L)×(L+1)对角阵,主对角元素σi为Y的奇异值.
对于理想的M阶信号,有σ1 ≥ σ2 ≥‧‧‧ ≥ σM ≥ 0,且其余奇异值均为零;对于含噪声信号,其余特征值可能因噪声的干扰变为非零值.对此,本文提出一种新型改进矩阵束算法,自适应判断奇异值矩阵阶数,消除因干扰变为非零值的奇异值,有效降低噪声对测量结果的影响.
依据奇异值在某点处会发生突变的现
(7) |
再定义相邻对数值之差:
(8) |
当相邻奇异值产生突变,则相邻奇异值的数量级会发生较大变化,也会取得相应的较大值,定义取最大值时的下标作为模型的阶数,即:
(9) |
将对应的奇异值σi全部置零,得到新的奇异值矩阵S',计算出频率成分数量为,即:
(10) |
此时取出V矩阵中前个右奇异向量构成(L+1)×矩阵V',删去V'的最后一行得到V1,删去V'的第一行得到V2,生成Y1和Y2:
(11) |
其中,Y1和Y2可分别表示为:
(12) |
(13) |
通过构造矩阵束,zi为的广义特征值,将求解zi的问题转化为求解广义特征值的问题,对上述极点信息的求取转化为求解
(14) |
式中:Y
获得zi与后,留数pi可利用最小二乘法的线性方程组计算得到:
(15) |
由
(16) |
综上,本文提出的新型改进矩阵束同步相量测量算法流程图如

图2 新型改进矩阵束同步相量测量算法流程
Fig.2 New improved matrix bundle synchronous vector measurement algorithm flowchart
1.3 算法仿真实验
为验证本文算法的优越性,利用MATLAB分别生成信噪比为60 dB的50 Hz标准电网信号、三次谐波、含有幅值大小为基波幅值10%的50.5 Hz、51 Hz、52 Hz以及55 Hz间谐波信号,对本文提出的新型改进矩阵束同步相量测量算法进行仿真测试,并与传统矩阵束算法、文献[
检验条件/Hz | 本文算法 | 传统矩阵束算法 | 文献[ | |||
---|---|---|---|---|---|---|
平均误差 | 最大误差 | 平均误差 | 最大误差 | 平均误差 | 最大误差 | |
50 |
2.2×1 |
2.66×1 |
2.53×1 |
2.83×1 |
3.04×1 |
2.34×1 |
50.5 |
3.69×1 |
2.05×1 |
4.69×1 |
2.51×1 |
8.15×1 |
5.48×1 |
51 |
2.56×1 |
3.83×1 |
1.70×1 |
1.10×1 |
1.80×1 |
1.05×1 |
52 |
4.75×1 |
1.16×1 |
2.60×1 |
2.20×1 |
3.19×1 |
2.01×1 |
55 |
1.08×1 |
1.94×1 |
8.60×1 |
5.56×1 |
4.17×1 |
1.60×1 |
由

图3 本文算法与文献[
Fig.3 Comparison of the measurement frequency results between this algorithm and the literature [10] algorithm

图4 本文算法与原始矩阵束算法测量频率结果对比
Fig.4 Comparison of the measurement frequency results between this algorithm and Traditional MP algorithm
为验证本文算法在噪声环境下测量精度影响,利用MATLAB分别生成含谐波和简谐波分量的不同电网模拟信号后进行仿真分析,如
检验条件 | 最大误差/Hz | 平均误差/Hz |
---|---|---|
2次谐波 |
2.47×1 |
1.10×1 |
6次谐波 |
2.37×1 |
1.17×1 |
10次谐波 |
2.79×1 |
1.03×1 |
14次谐波 |
2.27×1 |
1.18×1 |
18次谐波 |
2.53×1 |
1.13×1 |
含49.5 Hz间谐波 |
5.50×1 |
9.20×1 |
含55 Hz间谐波 |
1.94×1 |
1.08×1 |
含65 Hz间谐波 |
2.57×1 |
2.40×1 |
含75 Hz间谐波 |
2.68×1 |
1.28×1 |
含125 Hz间谐波 |
2.38×1 |
9.52×1 |
2 便携式电网同步相量测量装置设计
2.1 硬件单元设计
本文中便携式电网同步相量测量装置主要由稳压模块、电压调理模块、模数转换器AD7606、蓝牙模块HC-05、GPS/北斗模块ATK1218以及基于Android系统的移动终端等组成,系统硬件框图如

图5 装置硬件设计
Fig.5 Hardware circuit design

图6 便携式硬件装置实物图
Fig.6 Picture of portable hardware device
稳压模块将输入电压源信号转换为±9 V、5 V与3.3 V直流输出为硬件各模块供电, 使用芯片型号与转换示意如

图7 稳压模块示意图
Fig.7 Schematic diagram of voltage stabilizing module
电压采集电路包括微型精密电压互感器和低通滤波器两部分,220 V交流电压通过电压互感器转换为±5 V范围的交流电压,经低通滤波器有效抑制高频噪声.如
(17) |
式中:TR为电流型电压互感器一次绕组与二次绕组的匝数比,TR、R1与R2取值分别为1 、120 与1.5 .R3与C1构成的低通滤波器截止频率为:
(18) |
式中:R3与C1分别为820 、1 μF,截止频率为194 Hz,选用该截止频率可有效滤除高次谐波干扰.在此条件下,输出电压约为7.55P-P,在满足AD7606输入范围同时保持了较大电压,提高了测量结果的精度.

图8 电压调理电路
Fig.8 Voltage acquisition circuit
GPS接收器能够接收PPS信号,以此触发MCU的外部中断,实现同步采样.AD7606是16位SAR-ADC,可通过SPI接口实现与MCU之间的双向通信,采样率最高可达200 kHz,在保证蓝牙传输数据完整性且所采集电网信息可有效计算出电网参数的前提下,本文设定其采样率为1.2 kHz.蓝牙模块选用高性能主从一体蓝牙串口模块HC-05,波特率范围为4 800 bps~ 1 382 400 bps,本文选用波特率为115 200 bps,可保证传输的稳定性及足够的传输速率.通过PPS信号外部触发AD7606进行电压同步采样,以10 Hz的速率将采样数据通过HC-05发送至移动终端.
2.2 软件设计
2.2.1 同步采样控制方法
同步采样控制过程如
(19) |
式中:G为每个采样周期的时钟脉冲数;ƒs为采样率.若偏移后的频率为,则
(20) |
式中:TP为在X个PPS脉冲中实际测得的晶振脉冲数量.

图9 同步采样过程
Fig.9 Synchronized sampling process
选取最接近G的脉冲数G1与G2,且满足:
(21) |
将G1与G2轮流设置为定时器的定时控制间隔,动态调整AD7606的采样率,降低采样时间误差.
2.2.2 相量测量软件设计
本文在Windows操作系统中进行Android软件开发,Android系统开发环境主要包括JDK(Java Development Kit)、Android Studio和Android SDK(Software Development Kit)三部分.开发语言为Java编程语言,使用Redmi K20Pro(内存6 GB)作为软件的测试平台,处理器为高通骁龙855(主频2.84 GHz).
终端通过蓝牙接收MCU采集的电网数据与时间标签,利用新型改进矩阵束算法对数据进行计算处理,得到的电压幅值、频率及相角信息,以曲线形式直观展示测量结果.软件功能结构及使用界面如

图10 软件功能结构及使用界面
Fig.10 Software function structure and user interface
2.3 硬件误差分析与修正
电网同步测量系统中,硬件装置误差主要来源于ADC量化误差及失调误差.量化误差是一种原理性误差,由于芯片位数有限,量化时进行最低位数确定存在取舍问题;失调误差则是因为模拟输入电压与数字编码对应关系出现偏移,如

图11 ADC失调误差示意图
Fig.11 Schematic diagram of ADC misalignment erro
3 测试结果与分析
同步相量测量装置研制完成后,对其进行测试实验,如

图12 实验装置
Fig.12 Experimental equipment
分别使用本装置中的变间隔同步采样方式与固定间隔采样方式对信号发生器生成的标准50 Hz正弦波进行测量,补偿前后相角测量误差如

图13 固定间隔采样和变间隔采样的相角误差
Fig.13 Phase angle measurement results for fixed and variable interval sampling
为评估本文装置对实际电网信号幅值、频率和相角测量的准确性,使用传统的频率扰动记录
输入幅值/V | 测量均值/V | 电压幅值测量误差/% |
---|---|---|
46.66 | 46.652 3 | 0.010 9 |
52.50 | 52.526 9 | 0.038 5 |
58.33 | 58.344 2 | 0.020 3 |
64.16 | 64.182 4 | 0.031 9 |
70.00 | 70.013 0 | 0.018 6 |
输入频率/Hz | 测量均值/Hz | 频率测量误差/mHz |
---|---|---|
49.90 | 49.899 79 | 0.21 |
49.98 | 49.979 76 | 0.24 |
50.00 | 49.999 92 | 0.08 |
50.02 | 50.020 51 | 0.51 |
50.10 | 50.099 94 | 0.06 |
输入相角/(°) | 测量均值/(°) | 相角测量误差/(°) |
---|---|---|
10 | 9.978 6 | 0.021 4 |
60 | 59.965 1 | 0.034 9 |
90 | 90.001 6 | 0.001 6 |
180 | 179.966 0 | 0.034 0 |
210 | 209.946 3 | 0.053 7 |
电压幅值测试时使用高压放大器放大信号发生器电压,额定电压值为100/,即55.74 V,电压基准值为1.2倍的额定电压值,即70 V.由
为测量不同频率时装置的测量精度,调整信号发生器的输入频率,测量结果参见
为测量不同相角时装置的测量精度,输入频率固定为50 Hz,通过设置信号发生器中正弦信号的初相角改变输入电压信号的相角,相角测试结果如
使用本装置与传统设备分别对信号发生器产生的50 Hz标准交流电进行了频率测量,实验结果如

图14 50 Hz信号频率测量结果对比
Fig.14 Comparison of 50 Hz signal frequency measurement result
为测试实际环境下的本装置测量精度,使用本装置与FDR设备同时测量220 V实际电网信号,使用时间标签将二者的测量结果对齐,频率、相角测量结果分别如

图15 实际电网信号频率测量结果
Fig.15 Actual grid signal frequency measurement result

图16 实际电网信号相角测量结果
Fig.16 Actual grid signal phase angle measurement result
为测试本装置在实际运行中的可靠性,通过Android Profiler工具对移动终端的实时CPU占用率及内存使用情况进行了监测,结果如

图17 工作时移动终端的CPU占用率及内存使用情况
Fig.17 CPU usage and memory usage of the mobile terminal during work
4 结 论
本文针对传统同步相量测量算法受电网中间谐波影响较大的问题,提出了基于新型改进矩阵束的同步相量测量方法,经过仿真验证了本文算法对比传统算法的优越性,并以此设计了便携式电网同步相量测量装置,解决传统同步相量测量设备成本较高、布置较复杂、体积庞大的问题.之后将电网信号采集电路与移动终端相结合,降低了装置的硬件成本,经过与现有设备的对比测试,其对电网频率和相角测量误差分别小于0.000 51 Hz和0.053 7°,验证了本装置对电网频率及相角测量结果的准确性.本文所设计的装置实现了电网同步相量测量装置的低成本化、便携化以及轻量化,利于同步相量测量装置的普及,在电网相量参数的动态测量领域具有应用价值.
参考文献
ANGIONI A,LIPARI G,PAU M,et al.A low cost PMU to monitor distribution grids[C]//2017 IEEE International Workshop on Applied Measurements for Power Systems (AMPS).Liverpool,UK: IEEE,2017:1-6. [百度学术]
PEGORARO P A,BRADY K,CASTELLO P,et al.Line impedance estimation based on synchrophasor measurements for power distribution systems[J].IEEE Transactions on Instrumentation and Measurement,2019,68(4):1002-1013. [百度学术]
AHMAD T,SENROY N.Statistical characterization of PMU error for robust WAMS based analytics[J].IEEE Transactions on Power Systems,2020,35(2):920-928. [百度学术]
柳昀.适用于动态条件的同步相量测量算法的研究[D].重庆:重庆大学,2019. [百度学术]
LIU Y.Study on synchronous phasor measurement algorithm for dynamic conditions[D]. Chongqing: Chongqing University, 2019.(in Chinese) [百度学术]
刘洁波,黄纯,江亚群,等.基于强跟踪泰勒-卡尔曼滤波器的动态相量估计算法[J].电工技术学报,2018,33(2):433-441. [百度学术]
LIU J B,HUANG C,JIANG Y Q,et al.Dynamic phasor estimator based on strong tracking Taylor-Kalman filter[J].Transactions of China Electrotechnical Society,2018,33(2):433-441.(in Chinese) [百度学术]
曹磊,赵庆生,王旭平,等.基于改进泰勒加权最小二乘法的相量测量算法[J].电力系统保护与控制,2018,46(7):24-30. [百度学术]
CAO L,ZHAO Q S,WANG X P,et al.A phasor measurement algorithm based on improved Taylor weighted least squares[J].Power System Protection and Control,2018,46(7):24-30.(in Chinese) [百度学术]
张婕.同步相量测量装置的硬件设计与算法研究[D].太原:太原理工大学,2020. [百度学术]
ZHANG J. Hardware design and algorithm research of synchronous phasor measurement device[D]. Taiyuan:Taiyuan University of Technology,2020.(in Chinese) [百度学术]
金涛,陈毅阳,段小华,等.基于改进DFT的电力系统同步相量测量算法研究[J].电工技术学报,2017,32(17):1-10. [百度学术]
JIN T,CHEN Y Y,DUAN X H,et al.Research on synchronous phasor measurement algorithm of power system based on improved DFT[J]. Transactions of China Electrotechnical Society,2017, 32(17): 1-10.(in Chinese) [百度学术]
于华楠, 杜瑶, 郭树旭.基于压缩传感的高精度同步相量测量方法[J].吉林大学学报(工学版), 2018, 48(1): 312-318. [百度学术]
YU H N,DU Y,GUO S X.High-precision synchronous phasor measurement based on compressed sensing[J].Journal of Jilin University (Engineering and Technology Edition),2018,48(1):312-318.(in Chinese) [百度学术]
姚文轩.基于广域同步测量的电能质量参数检测关键技术研究[D].长沙:湖南大学,2017. [百度学术]
YAO W X. Research on key technologies of power quality parameter detection based on wide area synchronized measurement[D]. Changsha:Hunan University,2017.(in Chinese) [百度学术]
黄瑞,宋健,肖宇,等.基于改进矩阵束的动态同步相量测量算法[J].湖南大学学报(自然科学版),2022,49(2):135-142. [百度学术]
HUANG R,SONG J,XIAO Y,et al.Dynamic synchronous phasor estimation algorithm based on improved matrix pencil[J].Journal of Hunan University (Natural Sciences),2022,49(2):135-142.(in Chinese) [百度学术]
姚文轩,滕召胜,唐求,等.非同步采样下电力系统相量测量修正算法[J].仪器仪表学报,2013,34(7):79-85. [百度学术]
YAO W X,TENG Z S,TANG Q,et al.Power system phasor measurement and correction algorithm under asynchronous sampling[J].Chinese Journal of Scientific Instrument,2013, 34(7):79-85.(in Chinese) [百度学术]
刘晏池,王雪,刘佑达,等.基于分布式压缩感知的智能用电网络高分辨率谐波与间谐波测量[J].仪器仪表学报,2015, 36(10):2161-2166. [百度学术]
LIU Y C,WANG X,LIU Y D,et al.High-resolution harmonic and interharmonic measurement in smart power utilization network based on distributed compressed sensing[J].Chinese Journal of Scientific Instrument,2015,36(10): 2161-2166.(in Chinese) [百度学术]
段刚,严亚勤,谢晓冬,等.广域相量测量技术发展现状与展望[J].电力系统自动化,2015,39(1):73-80. [百度学术]
DUAN G,YAN Y Q,XIE X D,et al.Development status quo and tendency of wide area phasor measuring technology[J].Automation of Electric Power Systems,2015,39(1):73-80.(in Chinese) [百度学术]
尹宇翔.电能质量检测分析方法研究及其DSP实现[D].武汉:华中科技大学,2019. [百度学术]
YIN Y X.Research and the implementation on DSP of power quality detection and analysis methods[D].Wuhan:Huazhong University of Science and Technology,2019.(in Chinese) [百度学术]
滕召威. 基于嵌入式平台的电能质量在线监测系统研究 [D]. 北京: 华北电力大学, 2020. [百度学术]
TENG Z W. Research on power quality online monitoring system based on embedded platform[D]. Beijing: North China Electric Power University, 2020.(in Chinese) [百度学术]
刘志文,朱增坤,周治国,等.一种基于FPGA+ARM架构的μPMU实现[J].电测与仪表,2015,52(22):8-13. [百度学术]
LIU Z W,ZHU Z K,ZHOU Z G,et al.A design of μPMU based on the structure of FPGA+ARM[J].Electrical Measurement & Instrumentation,2015,52(22):8-13.(in Chinese) [百度学术]
WALLS F L,VIG J R.Fundamental limits on the frequency stabilities of crystal oscillators[J].IEEE Transactions on Ultrasonics,Ferroelectrics,and Frequency Control,1995,42(4):576-589. [百度学术]
ZHAN L W,ZHAO J Y,CULLISS J,et al.Universal Grid Analyzer design and development[C]//2015 IEEE Power & Energy Society General Meeting.Denver,CO,USA: IEEE,2015:1-5. [百度学术]
JIANG X C,ZHANG J M, HARDING B J, et al.Spoofing GPS receiver clock offset of phasor measurement units[J]. IEEE Transactions on Power Systems, 2013, 28(3):3253-3262. [百度学术]
KONONOV Y G, OSPISHCHEV M A, RUDNEV V Y. Current and voltage synchronous measurements using NI сRIO system with sigma-delta ADC NI 92хх and S.E.A. GPS modules [C]//18th International Workshop on ADC Modelling and Testing, 2014. [百度学术]
钱志鸿,刘丹.蓝牙技术数据传输综述[J].通信学报,2012,33(4):143-151. [百度学术]
QIAN Z H,LIU D.Survey on data transmission in Bluetooth technology[J].Journal on Communications,2012,33(4):143-151.(in Chinese) [百度学术]
刘思议,金涛,刘对.基于改进小波阈值去噪和RCRSV-MP算法的电力系统低频振荡模态辨识[J].电力自动化设备,2017,37(8):166-172. [百度学术]
LIU S Y,JIN T,LIU D.Power system low-frequency oscillation mode identification base on improved wavelet threshold de-noising and RCRSV-MP algorithm[J].Electric Power Automation Equipment,2017,37(8):166-172.(in Chinese) [百度学术]
ZHAN L W,LIU Y,YAO W X,et al.Utilization of chip-scale atomic clock for synchrophasor measurements[J]. IEEE Transactions on Power Delivery,2016,31(5): 2299-2300. [百度学术]