一种基于反馈迭代式频率跟踪的pmu动态数据处理方法
【技术领域】
[0001] 本发明属于电力系统及其自动化技术领域，具体设及一种基于反馈迭代式频率跟 踪的PMU动态数据处理方法。
【背景技术】
[0002] 数据采集单元的固定采样率为4曲Z，采样间隔A T=250 yS，对应于工作在额定 频率f"= 50化的电力系统，相当于每周波N=80点采样。当系统频率发生偏移时，原采样 信号不能构成一完整周波。通过软件测频获得当前频率化e，然后根据此频率对原采样点进 行二阶Lagrange插值拟合出对应于工频信号50化的等相位点，然后用运些等相位点作傅 里叶变换。W实现采样值的同步修正，保证DFT的计算为整周波数据，提高相量测量精度。 如图1所示。但是，当系统频率不稳定时，用于插值的频率如果存在偏差，便会使得插值数 据有所偏离，从而造成测量的误差。
[0003]同步相量算法模型BSEM度asic Sync虹0地asor Estimation Model)，低通滤波器 直接作用在DFT的实虚部上，如图2所示。该模型于IE邸标准C37. 118. r-2011中提出， 其主要思想为：在传统DFT算法的基础上，设置低通滤波器滤除带外频率W及频谱泄漏产 生的高频振荡分量，消除干扰量的影响。对应的公式为：
[000引其中，《。=23if。，f。为额定频率巧OHz) ;L为有限长单位冲击响应数字滤波器FIR (Finite Impulse Response)阶数；Wk为低通滤波器系数。
[0007] 离散傅里叶变换DFT (Discrete Fourier Transform)的平均化效应，要求测量相 量的时标选取为处理数据窗的窗中，使得平均化效应造成的量测误差达到最小。例如在工 频下，信号幅值发生振荡时（频率和相位为常量），有：
[0008]
稱
[0009] 可见计算所得相量的幅值为时间窗内信号幅值的平均值，该效应将影响PMU的动 态测量精度。当系统信号线性变化或者近似线性变化时，其时间窗内相量平均值与时间窗 中点时刻的实际相量相等，即相量时标选为时间窗中点时，由DFT平均化效应造成的误差 能被抵消掉，如图3所示；而信号非线性变化时，时间窗内相量平均值虽然无法与时间窗中 点时刻相量真实值一致，但将相量时标选在时间窗中点仍可大幅降低，使相量计算值最大 限度的与信号真值一致，如图4所示。
【发明内容】
[0010] 本发明的目的在于：针对现有技术存在的不足，提出了一种基于反馈迭代式频率 跟踪的PMU动态数据处理方法，该方法在测量系统出现频率偏移、幅值振荡、相位振荡等动 态条件下，仍然具有较高的测量精度。本发明提高了电力系统在非稳定状态下的测量数据 和控制状态的可靠性，有利于电网系统数字化高精度测控的安全有效。
[0011] 本发明采用的技术方案为：一种基于反馈迭代式频率跟踪的PMU动态数据处理方 法，在利用Lagrange插值计算等相位点实现电力系统频率跟踪计算的基础上，进一步规范 了插值所依据的频率参数，采用差分迭代的方法保证了频率值的周期稳定和精确对应，同 时，提出了依据此方法对插值方向和滤波器设计的规范要求，该方法包括如下步骤：
[0012] 步骤一，系统频率值由当前周波相量与上一周波相量（20ms前）的相角取差分得 至IJ，经过延时对应之后，作为频率跟踪处理环节的输入频率，使得系统的跟频补偿校正能够 紧随系统实时频率的变化而调整，W在随时间变化的动态环境中具有很高的测量精度；
[0013] 步骤二，基于DFT的平均化效应，要求测量相量的时标选取为处理数据窗的窗中， 由此对频率跟踪的Lagrange插值和数字滤波器的要求为时标点为中屯、，其两边的处理 数据点数完全相等。
[0014] 其中，步骤一中的频率值由当前周波相量与上一周波相量的相角取差分得到具体 为：频率值由当前时标点一前一后的两个相量的相位角作差分计算得到（保证计算频率 的时标位置对应），故频率的精确度取决于相量相位角度的精确计算，而相量的精确度又 取决于频率跟踪模块和滤波处理的效果。于是由此构成了一个反馈型的迭代处理方式，如 图5所示，初始频率即首次计算的缺省值选用工频50化，由于采用了频率跟踪算法，经过 Lagrange插值处理之后，所得相量为对应工频50化的等相位处理值，则每间隔20ms的计算 相量对应间隔一整周波信号；根据要求输出相量的频度为巧S=lOOHz，即每10ms间隔计 算一次相量，每两个计算间隔对应于一整个系统周波，因此则迭代的公式为：
[0015]
(1!'
[0017]其中，巧为相位角，化e为频率，为频率变化率，后方标注为时标位置描述，i表 at 示当前时标的计算点，i-l表示lOms之前一个计算点位，i-2表示20ms之前一个计算点位， 同理i+1表示未来下一次的计算点。
[001引即频率的计算值总要比相量晚一个点（10ms)，而频率变化率又要比频率晚一个点 (10ms)，此处的系统延迟至少为20ms，而频率跟踪采样的是上一周波即20ms的计算频率。 如图6所示，系统保留两个点的缓存，W最新的一组完整数据[i-2]作为输出，此时频率最 新点为[i-1]，相量最新点为[i]，然后W最新的频率[i-1]来作为Lagrange插值计算下一 个相量[i+1]，此二者从时标上而言正好相差一整个周波W至于不会造成系统相位的计算 偏离，再计算出下一组的频率和频率变化率数据，由此形成循环迭代。
[0019] 其中，步骤二中，要求依据测量相量的时标位置，其两边的处理数据点数完全相 等，具体为：插值时需要W时标选取点为中屯、，正反双向插值，W抵消插值数据位置造成的 时标偏移，如图7所示；且设计BSEM数字滤波器时，需要保证滤波器阶数L为偶数。由于滤 波器的卷积效应，每次计算相量时的实际处理数据长度为[N+kl]，一般情况下对整周波的 采样点数N都为偶数，故只有当滤波器阶数L也为偶数时，BSEM算法才能满足计算数据相 对时标点的中屯、对称，如图8所示。
[0020] 另外，该基于反馈迭代式频率跟踪的PMU动态数据处理方法，对系统的合理测量 范围需要有补充限定。具体为：实际测量系统中可能发生异常输入，需要设置系统的合理测 量范围作为口限。较为常见的一种明显异常状况为：系统无输入，即系统未有效接入时，由 环境噪声、电磁福射干扰等形成的微小信号或输入频率为零的直流信号，可能对频率跟踪 和滤波模块造成异常溢出的影响。因此需要限定正常的测量范围，如20~80化，超出此范 围的异常信号，使系统强制认定为50化工频或跳过跟频滤波处理模块的计算。
[0021] 本发明与现有技术相比的有益效果在于：本发明基于原有相量测量单元（PMU)的 算法模型BSEM，针对在发生频率偏移、幅值振荡、相位振荡等动态条件下，提出迭代反馈式 频率跟踪方法，同时进一步明确了数字滤波器的设计要求。在原有基础上，显著提高了PMU 在动态条件下的基波相量测量精度。
【附图说明】
[0022] 图1为连续采样点插值示意图；
[002引图2为同步相量算法模型BSEM图示；
[0024]图3为信号线性段时标位置窗首、窗中误差比对，其中，图3(a)为窗首时标，图 3(b)为窗中时标；
[0025]图4为信号非线性段时标位置窗首、窗中误差比对，其中，图4(a)为窗首时标，图 4(b)为窗中时标；
[0026] 图5为本发明一种基于反馈迭代式频率跟踪的PMU动态数据处理方法中，反馈迭 代式的频率跟踪算法处理流程图；
[0027] 图6为本发明一种基于反馈迭代式频率跟踪的PMU动态数据处理方法中，频率生 成和迭代方式示意图；