吴智利;赵庆生;陈惠英;韩肖清;郭贺宏
【摘 要】In consideration of unfavorable integration between tracking speed and measurement precision of conventional phasor measurement methods based on various constant-interval sampling, a ten-state linear Kalman-filter model for tracking power system RMS voltage, frequency, frequency variation rate as well as initial phase angle under the condition of low sampling frequency (4~8 sampling points per cycle) is presented. The algorithm is derived using the second-order Taylor expansion formula, and then the formulas of parameters above are established. According to the signal model, an appropriate time interval is selected to guarantee the response speed of algorithm to signal mutation. Sampling frequency, data window duration, and tracking precision to mutation signals are studied. Simulation results show that the algorithm has higer parameter estimation accuracy at low sampling frequency.%各种基于定间隔采样的传统相量测量算法在跟踪速度和测量精度上不能很好地统一，在低频采样情况(每周波4~8
个采样点)下提出一种基于线性卡尔曼滤波技术的10状态卡尔曼滤波模型，用于实时跟踪电力系统的电压有效值、频率、频率变化率和初相角，为系统提供精确的参数。结合二元泰勒展开公式，推导出各个参数的计算公式，并依据信号模型选取合适的时间区间保证算法对突变信号的响应速度。从采样频率、数据窗持续时间和算法对突变信号的跟踪精度方面对算法进行了讨论，结果表明使用该算法可以在低频采样下获得较高的参数估计精度。 【期刊名称】《电力系统保护与控制》
【年(卷),期】2014(000)015
【总页数】6页(P94-99)
【关键词】低频采样;卡尔曼滤波;二元泰勒展开;电力系统;实时相量测量
【作 者】吴智利;赵庆生;陈惠英;韩肖清;郭贺宏
【作者单位】太原理工大学电气与动力工程学院，山西 太原 030024;太原理工大学电气与动力工程学院，山西 太原 030024;太原理工大学电气与动力工程学院，山西 太原 030024;
太原理工大学电气与动力工程学院，山西 太原 030024;国网晋中供电公司，山西 晋中 030600
【正文语种】中 文
【中图分类】TM723;TM732
广域测量系统(WAMS)是一种有效的电网动态安全监控系统，相量测量算法的精度直接影响到电力系统稳定控制、故障分析及继电保护等高级应用的准确性[1]。常用的相量测量方法包括离散傅里叶变换(DFT)[2-4]、小波变换[5-6]、dq变换法[7-8]以及以最小二乘法[9-10]、卡尔曼滤波和最大似然估计为基础的线性或非线性的参数估计方法[11-12]。卡尔曼滤波是一种基于模型的最优估计方法，是一个不断预测、修正的递推过程。由于其在求解时不需要存储大量的观测数据，并且每次迭代运算量少，有利于实现高速运算，因此卡尔曼滤波被越来越多地应用于实时动态数据处理[13]。文献[14]使用线性卡尔曼滤波器的误差序增益和误差协方差阵以减小计算量。文献[15]探讨了复数型扩展卡尔曼滤波(ECKF)算法在电力系统畸变信号频率估计中的应用，但当算法收敛后，卡尔曼增益和误差协方差阵被限定在很小的数值上，影响了算法对突变信号的响应速度。文献[16]和[17]都使用两状态线性
纳米烟嘴卡尔曼技术计算电压的幅值和相位，但是无法直接算出频率，因此文献[17]使用锁频环(FLL)计算频率并与卡尔曼方法相配合以减小模型误差。
考虑到上述问题，本文在两状态卡尔曼滤波器模型的基础上引入频率和频率变化率参数(估计频率变化率的意义在于预测系统是否可能发生严重的频率骤降，并在事故发生之前提供合理的规避方案[18])，建立了10状态的卡尔曼滤波模型，并根据信号模型选取合适的时间区间提高算法对突变信号的响应速度。最后使用仿真算例验证算法有效性。
假定在采样过程中频率随着时间的推移线性变动，幅值在数据采样窗口内保持不变，电压信号可以表述为
式中：为电压峰值；为信号频率；为频率变化率，单位为；为信号初相角，单位为；为时间，单位为。
公式(1)可以写作如下形式：
使用多元泰勒展开公式将上述公式在和附近进行展开，展开后表达式如式(3)~式(5)所示：
其中定义状态变量为
参数为
应当注意的是，在式(3)中，仅仅只截取至泰勒展开式的一阶展开项，抛弃二阶及高阶展开项。如果按照预先选定的采样间隔对电压波形进行采样，则式(3)中的时间可以用第步的采样时间表示，式(3)可写作：
或者写成矩阵形式为
其中：为第步采样时刻的采样值；为第步采样时刻的测量矩阵；为第步采样时刻的状态矩阵；为第步采样时刻的随机噪声。
矩阵的取值取决于时间基准、采样频率及泰勒展开公式中保留的项数。可以预先进行离线计算再应用到之后的在线实时计算中。
假定系统的状态保持稳定(仅有每次采样时产生的随机误差)，那么相邻两次系统状态变量之间的关系表达式为
其中：为10×10的状态转移矩阵；为状态更新方程误差阵。式(7)和式(8)分别是电力系统电压信号模型的测量更新方程和状态更新方程，均以状态变量的形式表示。
综上所述，状态量及其内各个元素含义已知，根据卡尔曼递推公式，各个采样时刻的电压幅值、频率、频率变化率和电压相角的求取如下所示。在第个采样点处电压信号的峰值估计式为
同时，该采样点处的频率偏移估计式可以使用下述两个公式之一：
但是当或者的值偏小时，使用式(10)或者式(11)得到的估计值会有较大的误差，因此使用下述公式来得到频率偏移的精确值。
并使用式(10)和式(11)联合确定频率偏移量的正负：
同样出于计算精度的考虑，第次采样的频率变化率使用如下公式来计算：
初相角的计算可采用下面的公式：
2.1 滤波器模型
卡尔曼的状态更新方程和测量值更新方程需要写成如式(7)和式(8)的形式，在获得先验估计值以及先验估计方差后，卡尔曼滤波方程的一般递推步骤如下所示[19]。
导线测量法(1) 计算卡尔曼增益
(2) 计算后验估计值的方差
(3) 使用测量值校正先验估计值得到后验估计值
(4) 使用后验估计值推导下一次先验估计值及先验估计方差
当采样间隔固定时，式(16)、式(17)和式(20)是独立于测量值和估计值的，因此可以离线计算，只有式(18)和式(19)需要实时计算。图像识别
2.2 滤波器初始估计值
卡尔曼滤波需要使用一个先验估计值和对应的方差值作为滤波器的初始值进行后续的递推运算。比较简单的方法是对若干个测量值进行最小二乘估计求出状态量的初始估计值。使用下述方法，可以达到相同的精度。
其中：是由m个采样点对应的参数构成的m×n维观测方程矩阵；是m×1维当前测量值矩阵；是n×1先验估计向量。与状态量初值对应的误差阵为
2.3 滤波器迭代过程
观察观测方程系数矩阵式(5)及卡尔曼增益式(16)可知，随着时间取值增大，卡尔曼增益值逐渐降低，导致算法对信号突变的响应速度降低。为了避免这个问题，在迭代过程中，时间不可取值过大。经验证在式(5)中取值在0~0.04 s或0~0.06 s区间时，算法对突变信号的反应快速准确。以、采样频率200 Hz为例，在0.04 s内得到采样点称为一个数据窗，数据窗内依据卡尔曼公式进行迭代。迭代过程进入下一数据窗时，使用前一窗口末计算得到的最优估计值和方差作为本窗口的迭代初值，并强制将本窗口的时间参数置零，继续卡尔曼迭代过程。当采样间隔固定时，可根据采样时间预先计算得到窗口内每个采样点对应的参数值，无需重复计算。计算过程如图1中曲线所示。
图中表示电压瞬时值；为时间；表示滤波过程中的第个数据窗。H、K、P、Y与式(16)~式(20)中的含义相同。
2.4 算法适用范围
仿真恐龙制作算例使用的源信号频率偏移范围为-0.2 ~ -1 Hz、频率变化率为-0.05 ~ -1 Hz/s，信号有效
值为220 V，初相角为。源信号的采样频率为200 Hz，采样持续时间0.35 s。图2~图5显示了不同源信号下的参数估计结果。
根据仿真结果可得如下结论：
排油烟气防火止回阀
(1)信号的10状态模型适用于小幅度的频率偏移和频率变化率，在这种情况下，估算得到的频率值、频率变化率、电压有效值、相角值与真实值几乎相同。
(2)当频率偏移和频率变化率比较大时(=-1 Hz、=-1 Hz/s)，上述模型不再满足计算精度的要求，需要在泰勒展开式中保留更多的项以保证估计结果的精度。
(3) 经验证，在采样频率为200~400 Hz(源信号频率50 Hz)时，算法结果受采样频率的影响不大，频率变化率估计值对采样频率的变动最为敏感。当采样频率位于500~600 Hz时，频率变化率参数的绝对误差达到10-3级。随着采样频率的进一步提高，算法估计值将围绕真实值做大幅度波动。基于此，本算法更适用于200~400 Hz的低频率采样。
触摸调光ic(4) 算法性能同时受到源信号频率和频率变化率的影响，假定信号实际频率可用表示，则随时间变化，信号频率越接近工频则算法性能越好，否则越差。对频率和频率变化率参数分
别测试时可发现：当时，信号频率在工频信号±1 Hz范围内或当时，频率变化率在±2 Hz/s范围内时，算法仍可保证良好的精确性。
正常运行状态下电力系统母线电压是时刻变化的，母线电压的三个参量中，电压有效值和相位可能发生较大幅度的非正常变化，甚至是突变。算例使用上述信号波形对信号模型及其推导公式进行验证，求取频率、频率变化率、电压有效值以及电压信号的相角。信号初相角均为，采样频率为200 Hz，采样时间1 s；电压突变信号频率保持50 Hz不变，电压有效值由220 V突然跌落至190 V；相角突变信号频率保持50 Hz不变，初相角由突变至；以上两组信号突变的时间均位于0.2 s处；另一组为频率线性变化，信号频率在50 Hz的基础上以=0.05 Hz/s速度线性增加。计算结果如图6~图9所示，实线表示电压有效值突降时的估计结果曲线；点划线为相角突变时的估计结果曲线；虚线为频率线性增加时的估计结果曲线。算法的估计误差参见表1。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议