Volume 45 Issue 5
Oct.  2016
Article Contents

WANG Hong-bin, ZENG Xing-xing, ZHAO Hong, XIONG Xiao-fu, OUYANG Jin-xin, ZHANG You-qiang, CHEN Tao. Kalman Filter Algorithm Based on Backward Differentiation Formula of Delta Operator with Simulation[J]. Journal of University of Electronic Science and Technology of China, 2016, 45(5): 772-777. doi: 10.3969/j.issn.1001-0548.2016.05.011
Citation: WANG Hong-bin, ZENG Xing-xing, ZHAO Hong, XIONG Xiao-fu, OUYANG Jin-xin, ZHANG You-qiang, CHEN Tao. Kalman Filter Algorithm Based on Backward Differentiation Formula of Delta Operator with Simulation[J]. Journal of University of Electronic Science and Technology of China, 2016, 45(5): 772-777. doi: 10.3969/j.issn.1001-0548.2016.05.011

Kalman Filter Algorithm Based on Backward Differentiation Formula of Delta Operator with Simulation

doi: 10.3969/j.issn.1001-0548.2016.05.011
  • Received Date: 2015-06-11
  • Rev Recd Date: 2015-11-09
  • Publish Date: 2016-09-01
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(5)  / Tables(3)

Article Metrics

Article views(4661) PDF downloads(137) Cited by()

Related
Proportional views

Kalman Filter Algorithm Based on Backward Differentiation Formula of Delta Operator with Simulation

doi: 10.3969/j.issn.1001-0548.2016.05.011

Abstract: The relay protection based on the current principle is the most widely-used main-protection devices in power systems, its measurement circuit abnormity will directly lead to incorrect tripping or missoperation of protection equipment, which is the high-risk incident that power systems should be vigorously prevented. Polarity is an important characteristic in current measurement circuits of protection secondary systems, but its existing detection methods are still cumbersome and complex, and often bring wrong or leakage inspection events, which will seriously affect the safety and reliability of power systems. Therefore, it is necessary to research more simple and reliable polarity detection measures to ensure the correctness of polarity of current measurement circuits. So this paper defines the generalized ratio, and establishes a polarity diagnosis model of protection current measurement circuits, and analyzes the influence of comprehensive error of measurement circuits to the diagnosis of the polarity faults, and finally obtains the polarity fault criteria of current measurement circuits. Simulation analysis demonstrates the effectiveness of the proposed method.

WANG Hong-bin, ZENG Xing-xing, ZHAO Hong, XIONG Xiao-fu, OUYANG Jin-xin, ZHANG You-qiang, CHEN Tao. Kalman Filter Algorithm Based on Backward Differentiation Formula of Delta Operator with Simulation[J]. Journal of University of Electronic Science and Technology of China, 2016, 45(5): 772-777. doi: 10.3969/j.issn.1001-0548.2016.05.011
Citation: WANG Hong-bin, ZENG Xing-xing, ZHAO Hong, XIONG Xiao-fu, OUYANG Jin-xin, ZHANG You-qiang, CHEN Tao. Kalman Filter Algorithm Based on Backward Differentiation Formula of Delta Operator with Simulation[J]. Journal of University of Electronic Science and Technology of China, 2016, 45(5): 772-777. doi: 10.3969/j.issn.1001-0548.2016.05.011
  • 近年来,随着智能电网技术和数字化变电站技术的不断发展,基于新型智能电子设备、网络通信以及多测点信息综合比较判断的广域保护受到了国内外学者的广泛认可,但同时也增大了电力系统继电保护错误动作的风险。因此,有必要加强保护设备投运之前的检测工作,避免保护装置因错检漏检而造成重大损失。

    电流互感器(CT)作为继电保护测量回路的重要组成,其极性的连接直接影响了保护的正确动作,因此现有规程均要求CT在投运前进行离线校验和带负荷检测以确保极性的正确。目前,离线检测方法主要包括直流法、交流法和仪表法[1]。直流法将直流电源和毫安表按相同的正负极性分别接于CT一、二次侧,根据毫安表的指针位置来判断极性。交流法则需短接相同极性端口的CT一、二次线圈,并在二次侧接交流电压,通过测量两个端口的电压与一次侧电压的关系来判断极性。仪表法利用互感器校验仪直接检测CT极性,若其指示器没有指示,则其极性为正,否则为负。

    离线检测具有较高的精度,但无法避免接线等造成的极性错误,因此带负荷检测仍是确保CT极性正确的必要环节[2]。常用的极性带负荷检测方法主要通过电流电压间相角与有功无功间相角的对比、测量差动保护的差动电流以及比较线路同名相电流相位[3-4]。然而,无论是离线检测还是带负荷检测,现有方法均需要复杂接线和操作,不仅工作效率较低,而且容易因为人为疏忽导致漏检或错检情况出现[5]。文献[6-7]介绍了两起变压器差动保护误动的案例,文献[8]介绍了母线差动保护误动的案例,其原因均是由于电流互感器极性错误,且都造成了较为严重的停运事故。

    为此,有必要研究更加简单有效的极性检测方法以确保继电保护的可靠工作,避免测量回路极性错误引起保护装置错误动作,以及可能造成的区域性大面积停电事故。本文研究分析了继电保护电流测量回路极性错误的影响,通过测量回路广义变比的定义,建立了基于广义变比的极性辨识模型,从而构建了继电保护电流测量回路极性检测新方法。该方法可避免复杂接线和人工操作,无需添加硬件,具有精度高、易于实现的优点,是现有互感器极性检测方法的有效补充。

  • 目前,继电保护系统广泛使用的电流互感器分为电磁型和电子型,二者并存于智能电网建设初期,如图 1所示。图 1a为微机保护电流测量回路,一次侧大电流经电流互感器变换为二次侧小电流,再通过电缆、电流变换器、模拟低通滤波器(ALF)、采样保持电路(S/H)、多路模拟开关(MUX)和A/D转换电路转换为计算机可识别的电流采样值[9]图 1b为数字化保护电流测量回路,罗氏线圈通过电磁感应关系将一次系统电流变换为二次采样电流之后,经积分放大、低通滤波、相位补偿、A/D转换、编码、电光转换和合并单元转换为保护CPU所采用的数字量[10]

    以上各环节与测量回路极性直接相关的主要是电磁互感器的一、二次绕组的接法和罗氏线圈的具体绕向[9-10]。另外,测量回路各环节维护、更换以及合并单元软件程序配置或版本升级都有可能导致极性错误。一般情况下极性故障都会使保护装置误动,而小电流支路极性错误却不易被发现,一旦系统故障就会诱发该隐藏故障爆发。

    以差动保护为例,若接入一个节点(如母线、变压器单元等)的各支路CT极性正确,则在正常运行及差动保护区外故障时母线保护差动电流应小于其动作电流,有:

    式中,${i_{\rm{d}}}$为差动电流;${i_{2j}}$为接入母线第j条支路的二次电流;${i_{{\rm{op}}}}$为母差保护动作电流。

    根据文献[11],若接入差动保护回路第k条支路电流测量回路极性接反,即第k条支路二次电流${i_{2k}}$反向,则式(1)变为:

    式中,${i'_{\rm{d}}}$为第k条支路回路极性接反时的差动电流。

    通常情况下,由于$ - 2{i_{2k}}$的影响${i'_{\rm{d}}}$将大于${i_{{\rm{op}}}}$,保护装置动作,进而导致母线各支路被切除,致使接入母线所有支路用户停电。事实上,母线本身并没发生故障,而是其保护测量回路极性错误所导致的误跳闸。而当第k条支路为小电流支路时,系统正常运行下,此时$ - 2{i_{2k}}$非常小,基本不会改变差动电流的大小,保护装置不误动;当系统发生故障时,此时$ - 2{i_{2k}}$将远大于${i_{\rm{d}}}$,从而导致保护装置错误动作,以致给电网造成巨大损失。

  • 定义电流测量回路广义变比为流进CT一次侧电流与保护CPU所得二次电流的比值,用${N_{\rm{g}}}$表示,即图 1所示电流测量回路首端一次系统电流${i_1}$与末端保护CPU所得电流采样值${i_2}$之比,有:

    通常用额定变比N表征CT变换电流的关系,定义为一次额定电流与二次额定电流的比值。如对于二次额定电流为5 A的电磁型电流互感器,当一次额定电流为1 000 A时,选择的电磁型电流互感器的额定变比为N=1 000/5。由此可得广义变比标幺值,即为广义变比有名值${N_{\rm{g}}}$与其额定变比N的比值,表示为${n_{\rm{g}}}$,有:

    约定正常情况下,电流测量回路极性正确时,${N_{\rm{g}}}$的符号为正,即${n_{\rm{g}}}$为正;当电流测量回路的极性接反时,${i_2}$反向、符号变化,此时${N_{\rm{g}}}$的符号变为负,即${n_{\rm{g}}}$为负。因此,电流测量回路极性异常的监视可转化为对测量回路广义变比正负符号的辨识。

  • 根据广义变比的定义,欲通过求解其辨识值检测电流测量回路极性异常情况,必先获取系统一、二次电流值。然而,工程实际中只能采集到二次电流值,因此就要解决在一次电流未知的情况下求解广义变比的问题。

    实际系统中,每套保护装置特别是电流差动原理的保护,其接入电流反应了一次系统的连接关系。因此,根据一次系统拓扑关系,利用基尔霍夫电流定律,即在任一瞬时,流出(或流入)某元件电流的代数和恒等于零,即:

    式中,${i_{1i}}$为流过元件第i条支路的一次电流。在电力一次系统中,无论是连接在母线上的各支路电流还是任意划定的多个元件组成的广义节点,始终满足上述定律。在此约定流入某广义节点电流为正,则流出为负。如图 2所示,${i_{11}}{i_{12}} \cdots {i_{1i}} \cdots {i_{1n}}$为流过广义节点G的某相各一次电流,经电流测量回路转换为对应保护二次电流${i_{21}}{i_{22}} \cdots {i_{2i}} \cdots {i_{2n}}$。根据广义变比定义得到:

    式中,${N_{{\rm{g}}i}}$为第i条支路电流测量回路的广义变比。将式(6)代入式(5),等式两边同时除以额定变比N,则任意采样时段有:

    式中,$n_{{\rm{g}}i}^{(j)}$为在采样时刻j下第i条支路电流测量回路广义变比的标幺值,即为在采样时刻j下第i条支路电流测量回路广义变比的有名值$N_{{\rm{g}}i}^{(j)}$与其额定变比N的比值。

    设对测量回路的监视时间段为${t_1}$~${t_m}$,其中第j个采样时刻第i条支路的保护CPU计算值为$i_{2i}^{(j)}$,每组计算值均应满足式(7),则m(m>n)个时刻的计算值可列出如下方程组:

    式(8)的解为m个采样时刻不同广义变比的集合。然而,由于电流互感器角差一般较小[2],各采样时刻广义变比的差别很小,利用奇异值分解法的迭代求解,可获得广义变比的平均辨识值ngi。因此,由式(8)可得:

    将上式表示为矩阵形式,有:

    式中,

    由此,根据奇异值理论[12-13]可将A(ARm×n)分解为:

    式中,SD分别是$m \times m$阶和$n \times n$阶正交矩阵;$V = {\rm{diag}}({\lambda _1},{\lambda _2}, \cdots ,{\lambda _n})$是对角矩阵,其对角元素为A的奇异值,并按降序排列。

    计算中使采样组数m远大于电流测量回路数n,由式(5)、式(6)可知A的列向量是线性相关的[14],且当m足够大时,可使${\rm{rank}}(A) = n - 1$,此时超定齐次方程组式(10)只有一个基础解,为排除零解,加入约束条件$||{n_{\rm{g}}}|| = 1$。将式(11)代入式(10),可得如下转换关系:

    其中,单位正交矩阵DT可逆,上式可化简为:

    式中,${V^{\rm{T}}}V = {\rm{diag}}(\lambda _1^2\;\lambda _2^2\; \cdots \;\lambda _{n - 1}^2\;0)$;$D = {[{d_1}\;{d_2}\; \cdots \;{d_n}]^{\rm{T}}}$。

    由此可得DT,故可解得:

    式中,dn为正交矩阵D的最后一行,即ATA的最小特征值所对应的特征向量。

  • 考虑到继电保护电流测量回路各环节误差对其极性故障诊断的影响,可视为其二次电流采样值中引入一个综合误差项,有:

    式中,${i'_2}$为当测量回路无误差时的二次电流;e为综合误差。将式(15)代入式(4)可得:

    当电流测量回路极性正确时,一次电流与二次电流之比为额定变比,即${i_1}/{i'_2} = N$。由式(15)、式(16)可得:

    当电流测量回路极性接反时,${i'_2}$反向、符号变化,则${i_1}/{i'_2} = - N$。同理可得:

    正常运行状态下,电流测量回路综合误差绝对值$|\varepsilon |$应小于10%[15],由式(18)得极性故障时广义变比的取值范围${n_{\rm{g}}} \in ( - 1.111, - 0.909)$。为使方法可靠性更高,设定可靠系数${r_n} = 1.1$,得极性故障判据为:

    根据极性故障判据的取值范围,由2.2节所得电流测量回路abc相各支路广义变比辨识值,即可得如图 3所示的继电保护电流测量回路极性在线诊断流程。

  • 图 4所示为采用PSCAD/EMTDC搭建的母线差动保护电流测量回路极性故障检测仿真系统。

    图 4中各电流互感器额定变比均为N=1 200/5,单相电源MN相同,其额定电压为220 kV、频率为50 Hz,其负载分别为Z1=2 500 Ω+j942 Ω,Z2= 1000Ω+j157 Ω和Z3=150 Ω+j31.4 Ω。

    设定图 4仿真系统中各支路电流测量回路综合误差为零且极性正确,正常运行状态下,根据图 2采集图 4中母线保护装置各二次电流${i_{2i}}$,采样时段为0.1~0.2 s,如图 5所示。

  • 假定图 4仿真系统CT1和CT2的相位误差为0.225°且无幅值误差。分别设定测量回路所有支路极性正常、仅CT1或CT2所在测量回路支路极性接反以及CT1和CT2所在测量回路支路极性同时接反,根据2.2节方法计算测量回路各支路广义变比,如表 1所示。

    广义 变比 各支路极性均正确 仅CT1所在 支路极性接反 仅CT2所在支路极性接反 CT1、CT2所在 支路极性同时接反
    ng1 0.856 5 -0.856 5 0.856 5 -0.856 5
    ng2 0.922 0 0.922 0 -0.922 0 -0.922 0
    ng3 1.029 5 1.029 5 1.029 5 1.029 5
    ng4 1.000 1 1.000 1 1.000 1 1.000 1
    ng5 1.000 0 1.000 0 1.000 0 1.000 0

    表 1可知,当各测量回路支路极性均正确时,各支路广义变比符号全为正;仅CT1或CT2所在测量回路支路极性接反时,CT1或CT2所在支路广义变比符号为负,其他支路广义变比符号为正;当CT1、CT2所在测量回路支路极性同时接反时,CT1、CT2所在支路广义变比同时变负,其他支路广义变比符号为正。由此,验证了可根据电流测量回路广义变比的正负符号辨识其极性故障的结论。

  • 系统实际运行条件下,电流测量回路各环节误差不可回避,为验证测量回路误差对极性故障诊断的影响,在图 4的仿真系统中设定CT1所在支路电流测量回路极性接反且存在如表 2所示的误差情况,其中ε为幅值误差、θ为相位误差(超前为正、滞后为负)。计算CT1所在支路广义变比ng1,如表 2所示,CT1所在支路测量回路极性接反时在不同误差下的广义变比辨识值。

    q/(°) e/%
    -10 -5 0 +5 +10
    1 -1.224 3 -1.159 7 -1.101 8 -1.049 4 -1.001 6
    0.45 -1.153 1 -1.092 4 -1.037 4 -0.988 2 -0.943 2
    0 -1.099 1 -1.041 1 -0.998 3 -0.944 2 -0.899 4
    -0.45 -1.041 7 -0.986 7 -0.937 3 -0.892 6 -0.852 3
    -1 -1.019 0 -0.963 7 -0.915 7 -0.871 7 -0.828 3

    表 2可知,电流测量回路误差对其极性故障的诊断影响非常显著。当电流测量回路极性接反时,其广义变比的辨识值始终为负;随着测量回路幅值误差的增大,广义变比辨识值增大;随着测量回路相位误差的增大,广义变比辨识值减小。另外,综合误差数值越大,其广义变比辨识值偏离理想变比-1越大,这与实际运行经验相吻合。测量回路综合误差在规程规定范围之内时[3],即$\varepsilon \in ( - 10{\rm{\% }},{\rm{ + }}10{\rm{\% }}),\theta \in ( - 1^\circ ,{\rm{ + }}1^\circ )$,在本文研究所设定的可靠系数1.1下其广义变比辨识值满足极性故障判据。由此可见,通过广义变比辨识值取值范围的判断可以排除测量回路综合误差对极性故障检测的影响,从而验证了本文方法的有效性。

    广义变比相对负载灵敏度定义为单位负载量变化所导致的变比值变化为:

    设定图 4仿真系统仅CT1所在支路极性接反,其幅值和相位误差分别为-10%与1°,逐一改变CT1所在支路负载大小,分别为0.5Z1、1.5Z1、1.9Z1、2Z1以及2.5Z1,并统计仿真系统所有支路二次电流采样值最大值与最小值的平均倍数f,并判断能否诊断其极性故障。表 3为本文研究的极性故障诊断方法灵敏度仿真结果。

    负载 f 广义变比 能否诊断极性故障
    ng1 ng2 ng3 ng4 ng5
    0.5Z1 6.8 -1.199 1.105 0.951 0.999 1.000
    Z1 13.3 -1.224 1.093 0.971 1.000 1.000
    1.5Z1 19.9 -1.187 1.055 0.984 1.000 1.000
    1.9Z1 25.3 -1.223 1.049 0.985 1.000 1.000
    2Z1 26.7 -1.233 1.049 0.985 1.000 1.002 不能
    2.5Z1 33.4 -1323 1.050 0.982 1.001 1.000 不能

    根据表 3的仿真结果可知:即便保护系统角差和幅值误差处于极端情况下,系统各二次采样电流最大值与最小值的平均倍数f高达25时,本文方法都能正确检测出极性故障情况,具有较高的灵敏度。

  • 目前,电力系统中二次系统测量回路极性错误问题难以避免。本文提出了继电保护测量回路的广义变比概念,建立了继电保护电流测量回路极性诊断的数学模型,并基于奇异值分解法提出了电流测量回路极性故障的判据,建立了基于广义变比的继电保护测量回路极性检测方法,通过仿真算例验证了该方法的有效性。该方法仅通过软件算法即可实现,无需额外接线或添加硬件设备,可以作为现有极性检测方法的补充或正确性参照,是一种简单有效的继电保护电流测量回路极性检测新方法。

Reference (15)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return