北极星智能电网在线讯:摘要：针对目前多端线路故障定位中分支接点附近测距死区较大的问题,提出一种基于正序分量的精确故障定位新方法。该方法不以数值范围作为故障判别依据,而是采用双端方程有解和无解两种对立状态实现故障支路的判别,从而有效减小分支接点附近测距死区的范围。在求解过程中,同时实现故障支路判别与故障测距,无需事先判别故障支路。新方法定位及测距结果不受故障类型、过渡阻抗及系统运行方式等的影响。PSCAD仿真结果,验证了所提算法的正确性、有效性和高精度性。
关键词 : 故障定位; 正序基波分量; 测距死区; 虚拟故障点; 抽象双端线路; 判别矩阵;
（来源：微信公众号 电网技术 ID：dwjs1957 作者：刘瑞麟, 邰能灵, 范春菊, 黄文焘, 刘志远）
0 引言
传输线路故障导致的停电、电压不稳等事故通常会造成巨大的经济损失,且事故时间越长损失越大,因此故障发生之后,精准、及时的锁定故障所在位置,并迅速派出工作人员进行故障点的事故抢修,可有效保障负载端的可靠用电,减小损失[1-6]。
国内外对于多端线路故障定位问题的相关研究已经进行了数十年,目前主流的方法主要基于：行波、长线方程和高频脉冲等。行波法由于对波头的识别精度不高,结果通常不够精确[7-12],而高频脉冲在实际运行系统中应用效果也不理想[13-14]。基于长线方程的传统算法在分支接点附近存在一个较大的测距死区[15-18],而文献[19]依据故障点前后测距函数相位会发生突变的特性,将相位突变点视为故障点进行故障定位,文献[20]利用过渡阻抗的纯电阻性质实现故障定位,两者对减小测距死区的范围都有一定的效果,但均只适用于T型线路。
为解决上述问题,本文提出精确故障定位算法,本算法基于长线方程,利用双端正序基波分量找到远离盲区的虚拟故障点,再以所得虚拟故障点与邻近分接点构造新的双端问题,最后用各端两两对应求得的结果作为元素构造判别矩阵。通过所得矩阵特征量,判别出故障支路,同时给出精确测距结果。
1 精确故障定位原理
为实现高精度定位,消除架空平行线路间互感的影响,本文选用基于均匀传输线的长线方程作为基本测距原理,以正序基波分量作为计算电气量。与传统定位原理通常需要先进行故障支路判别,再通过简化线路实现定位不同,本文的精确故障定位原理可同时实现故障支路判别和精确测距。
1.1 双端定位原理基础
当传输线路上无支路连接时,如图1所示。
图1 双端线路简图
定义电流从各端母线向区域内汇集为正方向,U˙MU˙M、U˙NU˙N、I˙MI˙M、I˙NI˙N分别为M、N端母线处的正序电压和电流,LMLM、LNLN分别为故障点到M、N两端母线的距离,U˙U˙、I˙I˙为M端母线到故障点处之间任
一点的正序电压和电流,且该点到M端母线处的距离如图1中所示为LL,如图1中k处,则均匀传输线的长线方程为
U˙=U˙Mch(−γL)+I˙MZCsh(−γL)U˙=U˙Mch(−γL)+I˙MZCsh(−γL) (1)
I˙=I˙Mch(−γL)+U˙MZCsh(−γL)I˙=I˙Mch(−γL)+U˙MZCsh(−γL) (2)
式中：ZCZC为传输线路的正序特性阻抗;γγ为传输线路的正序传播常数。
将线路均匀分成DD段,从两端开始分别计算每小段端点的正序基波电压和电流值,可以得到故障点的位置
|ΔU˙F(ik)|=min{|ΔU˙F(i)||i=1,2,...,D}|ΔU˙F(ik)|=min{|ΔU˙F(i)||i=1,2,...,D} (3)
LF=(LM+LN)⋅ikDLF=(LM+LN)⋅ikD (4)
ΔU˙F(i)= (U˙M+ZCI˙M2e−γ(LM+LN)iD+U˙M−ZCI˙M2eγ(LM+LN)iD)− (U˙N+ZCI˙N2e−γ(LM+LN)D−iD+U˙N−ZCI˙N2eγ(LM+LN)D−iD)ΔU˙F(i)= (U˙M+ZCI˙M2e−γ(LM+LN)iD+U˙M−ZCI˙M2eγ(LM+LN)iD)− (U˙N+ZCI˙N2e−γ(LM+LN)D−iD+U˙N−ZCI˙N2eγ(LM+LN)D−iD)(5)
式中：ikik指与实际故障点最接近的小段端点序号;LFLF为计算所得故障点距离M端母线的距离。
1.2 多端线路精确故障定位原理
由于多端线路中,非故障支路会在分支接点处注入额外的电流,直接使用双端量得出的定位结果并不准确。因此传统方法通常采用先识别故障所在支路,再简化为双端线路后实现定位。
区别于传统原理,本文首先定义如下概念：
1)基准支路。作为双端线路中的一端,用于计算虚拟故障点f′Vf′V的所在位置,且计算抽象双端线路长度时,以本端母线处作为基准。
2)对照支路。与基准支路对应,作为求解虚拟故障点f′Vf′V的另一端。
3)虚拟故障点。从基准支路出发,与对照支路构成双端线路,直接按照求解双端线路的方法,求得参考点f′Vf′V,该参考点并非实际的故障所在位置。
4)抽象双端线路。以虚拟故障点和邻近分支接点作为两侧母线,人为构造的双端线路,线路长度可以为负数。
5)计算周期。任选两条支路,分别作为基准、对照支路,直至求得抽象双端线路对应方程解的过程。
6)旁路。一个计算周期中的其余支路。
多端线路中,任取两条支路分别作为基准支路、对照支路,其余为旁路,并根据双端原理求得
对应的虚拟故障点f′Vf′V。
根据传输线路长线方程有
U˙F=U˙Mch(−γL)+I˙MZCsh(−γL)= U˙M+ZCI˙M2e−γL+U˙M−ZCI˙M2eγL= U˙Me−γL+eγL2+ZCI˙Me−γL−eγL2U˙F=U˙Mch(−γL)+I˙MZCsh(−γL)= U˙M+ZCI˙M2e−γL+U˙M−ZCI˙M2eγL= U˙Me−γL+eγL2+ZCI˙Me−γL−eγL2 (6)
为不失一般性,令
I˙M=ex⋅ej θ I˙M=ex⋅ej θ (7)
代入式(6)并对xx求偏导
∂U˙F∂x=ZCe−γL−eγL2ex⋅ej θ =ZCsh(−γL)ex⋅ej θ ∂U˙F∂x=ZCe−γL−eγL2ex⋅ej θ =ZCsh(−γL)ex⋅ej θ (8)
式(8)显然满足复数域求导的柯西黎曼条件,因此：
dU˙Fdx=∂U˙F∂xdU˙Fdx=∂U˙F∂x (9)
再将双曲线泰勒展开：
sh(−γL)=−γL−(γL)33!−(γL)55!⋯sh(−γL)=−γL−(γL)33!−(γL)55!⋯ (10)
代入式(8)可得：
dU˙Fdx=ZC(−γL−(γL)33!−o[(γL)3])ex⋅ej θ dU˙Fdx=ZC(−γL−(γL)33!−o[(γL)3])ex⋅ej θ (11)
分别求解实部和虚部：
同时
其中：z=R+jX为传输线单位长度串联阻抗;y=G+jB为传输线单位长度并联导纳。
因此可得：
ZC(−γL−(γL)33!−o[(γL)3])= zy−−√(−zy−−√L−(zy−−√)36L3−o[(zy−−√L)3])≈ −zL−z2y6L3= −(R+jX)L−L36jB(R+jX)2= (13RXBL3−RL)+j[16(X2−R2)BL3−XL]ZC(−γL−(γL)33!−o[(γL)3])= zy(−zyL−(zy)36L3−o[(zyL)3])≈ −zL−z2y6L3= −(R+jX)L−L36jB(R+jX)2= (13RXBL3−RL)+j[16(X2−R2)BL3−XL](14)
即：
进而可知：
又由于系统阻抗ZM、ZN、ZS及传输线路阻抗ZC均是高感抗的(至少X>R),故电流滞后电压45°~ 90°,即电流相角θ在(-90°~-45°)范围内,因此：
由式(17)可知,电压(矢量)关于电流幅值(标量)的导数实部和虚部均为负,因此随着电流幅值增大,传输线路上的压降程度增大,电压沿线路所呈下行趋势的斜率变大,程度加强,进而虚拟故障点
f′Vf′V必在故障实际位置的远离分接点一侧。故该虚
拟故障点和邻近分接点可以作为两个虚拟端点构
造抽象双端线路。同时由于虚拟故障点f′Vf′V更加远
离传统算法中的测距死区,故有助于减小测距死区的范围。
在构造抽象双端线路时,由于每次选取不同支路依次定义为基准支路、对照支路,故会出现如下3种情况。
1)当以故障支路为基准支路时,可以构成存在有效解的抽象双端线路,因此求得精确的故障位置,如图2所示,其中横轴为三端输电线路,纵轴为沿线路分布的电压。
图2 故障位于基准支路时的原理
2)当以非故障支路为基准支路,且故障发生在对照支路上时,由于虚拟故障点比分接点距离基准支路侧母线更远,对应构造的抽象双端线路问题中,线路长度为负,电压呈现发散趋势,如图3所示,因此双端方程无解,可以-1标识。
图3 故障位于对照支路时的原理
3)当以非故障支路为基准支路,且故障发生在旁路上时,由于线路L1、L2上无电流分流或注入,计算所得虚拟故障点理论上即为分支接点,考虑误
差(包括二次侧设备测量误差),实际计算结果既可能在基准支路上,也可能在对照支路上,因此双端方程可能有解也可能无解,如图4所示。
图4 故障位于旁路时的原理
综上,将每个计算周期中计算出来的故障点位置信息作为矩阵元素,构造成判别矩阵,该矩阵行坐标对应基准支路,列坐标对应对照支路,以图5中三端线路故障位于L1支路为例,判别矩阵如下：
其中：D_Lmn(D_Lms)表示以M端L1支路为基准支路;N(S)端L2(L3)支路为对照支路时方程所得解,单位为km;*表示该项对应方程可能有解,也可能无解;-1表示对应方程无解。
图5 三端线路简图
上述判别矩阵具有如下特征：
1)故障所在支路选为基准支路时,对应的行内矩阵元素均为有效解,如式(18)中第1行;
2)故障所在支路选为对照支路时,对应的列内矩阵元素均为无解,标记为-1,如式(18)中第1列;
3)其它各行、列均不满足1)2)中所述特征;
4)故障所在支路对应行中的各个有效解为实际故障点到基准支路母线处的测距结果。
由上述判别矩阵的特征可知,有且仅有故障所在支路对应的行和列满足：对应的行均有解,而列均无解这一特征,因此利用该判别矩阵可同时实现故障支路判别与故障测距。
2 精确故障定位算法
2.1 三端故障定位算法
根据精确故障定位原理,三端线路算法分为如下4个步骤：
1)按顺序将各支路依次分别选为基准支路、对照支路和旁路,并求得对应的虚拟故障点f′Vf′V。见图2—4中的虚线。
2)以所得虚拟故障点和邻近分接点作为抽象双端线路的两侧虚拟端点,依据精确故障定位原理,求解对应双端方程,完成一个计算周期。见图2—3中的实线。
3)循环步骤1)、2)保证每次选取不同支路依次定义为基准支路、对照支路和旁路,并遍历所有可能的排列组合情况,即含有N条支路的线路,则循环(N2-N)次。
4)以上述(N2-N)次计算周期中计算出来的故障点位置信息,构造判别矩阵,如式(18)。
根据精确故障定位原理可知,此时矩阵中必然存在某列数值全为-1,而某行数据全为大于零的有效解,上述行、列同时对应的支路,即为故障所在支路,且上述行中各有效解即为计算所得的故障实际位置,对该行数据取算数平均值,得到测距结果。
2.2 多端故障定位算法
依据精确故障定位原理,求解判别矩阵中元素时,只计算拓扑结构上直接相连的对应支路。此处以故障在L1支路上为例：
此时,若多端线路在拓扑结构上近似对称,则式(19)中矩阵第4或5列有可能出现全无解,同时列所在支路对应行可能出现全有解的情况,进而无法判断是L1支路故障,还是L3或L5支路发生故障,为避免这种极端情况的发生,需将拓扑上不直接相连的支路,视为直接相连,并对算法进行相应修正。
各支路轮流被视为基准支路计算时,若出现基准支路与对照支路拓扑结构不直接相连的情况时,将对照支路带入计算的线路长度修正为：
修正后对照支路长度=基准支路邻近分接点到对照支路对应母线处的最短距离
如图6中,选L1为基准支路,L3为对照支路,则支路L3修正之后的长度为原L3支路与L2支路线路长度之和。
图6 四端线路简图
修正后,判别矩阵如下,满足本算法中的要求：
又例如,故障位于L2支路,则判别矩阵如式(21)所示,其余步骤同三端算法一致。
对于图7中含有3个分支接点,7条支路的五端线路,计算时若出现基准支路与对照支路拓扑结
构不直接相连的情况时,与上文一致,对照支路带入计算的线路长度修正为：
图7 五端线路简图
修正后对照支路长度=基准支路邻近分接点到对照支路对应母线处的最短距离
同时在计算时考虑本次计算循环中所选取的基准支路、对照支路在拓扑结构上相连的旁路在分支接点处汇入的电流,其余具体步骤与上文中所述一致。通过上述修正,可得其判别矩阵如下：
对于端数更多的线路,例如含有4个分支接点,9条支路的六端线路等等,当线路中总计存在N (N>4)条支路时,令其各条支路分别命名为T1~N,且假设故障位于T2支路上,则依据本文所述原理,可以得到其判别矩阵的通用表达形式如下：
3 算例仿真
为验证本文所提出的算法正确性,本文利用PSCAD软件进行仿真分析,采用全波傅里叶滤波,窗口时间20 ms,采样频率2400 Hz,每周波采48个点。
3.1 三端线路故障定位仿真验证
以图5所示线路中故障发生在L1支路,故障类型为A相短路接地,过渡电阻10 Ω,故障点距离分支接点O处30 km的事故为例,并选取L1为基准支路,L2为对照支路,其余参数见表1。此时依据长线方程计算所得正序基波电压沿上述两支路分布情况如图8所示,其中粗虚线、粗实线分别为算法第一、二步中电压沿线路的分布情况,细实线为辅助延长线,UPFCUPFC指正序基波电压。可见仿真结果所绘图8与图3中所示原理完全一致。
表1 三端线路参数
图8 线路沿线正序基波电压
由式(24)可知,仿真计算所得判别矩阵中,L1支路所在行均有解,列均无解,符合算法原理中所述特征,有效识别出故障位于L1支路,并对所在行各有效值求解算术平均值,测得故障距离分接点O处29.990 34 km,测距误差仅为0.016%。故算法的正确性在仿真中得到证实。
进一步检验不同故障支路、故障类型和过渡电阻下算法的有效性和高精度性,见表2。表3为多端线路参数。
3.2 多端线路故障定位仿真验证
表5表明当故障位于分支接点附近时,本文方法能够正确判别故障所在支路,同时保证较好的测距精度,故本算法对减小多端线路中测距死区的范围具有较好的效果。
本文方法中由于虚拟故障点必然比实际故障点更加远离分支接点,同时图4中所示原理在构造判别矩阵时已经考虑了互感器可能存在的测量误差,当二次侧设备测量误差较大时,本文方法会受一定影响,但是可以加强滤波,减小二次侧设备测
表2 三端线路在不同故障类型、过渡电阻下的仿真结果
表3 多端线路参数
量误差的影响。因此本文方法对于故障支路的判别较为准确,在测量过程中存在一定误差时,依然能够保持较好的测距精度。
表6列出了三端线路分支接点附近发生A相高阻接地故障时,本文方法与传统方法(以文献[20]为例)的比较结果,其中故障距离均指距离O点的长度,对比仿真均采用本文仿真时所用系统和参数,但因数据处理及其他因素影响,仿真和计算分析环境可能略有差异。从表6可以看出,传统方法在故障位于分支接点附近时失效,由于无法正确判别故障所在支路,进而无法正确实现故障测距,表6中以*号表示传统方法无法正确测距的情况。
4 结语
针对传统多端线路故障定位中普遍存在的测距死区问题,本文提出一种基于正序基波分量的多端线路精确故障定位算法。本算法基于正序基波分
表4 多端线路在不同故障类型、过渡电阻下的仿真结果
表5 多端线路分支接点附近不同故障类型、过渡 电阻下的仿真结果
量,采用双端方程有解和无解两种对立状态作为故障支路判别依据。大量仿真结果显示该算法不受故障类型、过渡阻抗、系统阻抗及系统运行方式等的影响,能保证较高的测距精度。此外本算法将多端
表6 三端线路中分支接点附近故障时测距结果对比
线路故障支路判别与故障测距融为一体,无需事先判别故障所在支路,同时对减小分支接点附近的测距死区范围也具有较好的效果。
作者简介：
刘瑞麟(1992),男,硕士研究生,研究方向为电力系统保护与控制、电网可视化等,E-mail：lrldyx16380@163.com;邰能灵(1972),男,教授,通信作者,博士生导师,研究方向为电力系统继电保护、高压直流输电、新能源等,E-mail：nltai@sjtu.edu.cn;
范春菊(1967),女,副教授,长期从事继电保护及电站综合自动化的研究及教学工作;
黄文焘(1989),男,博士,讲师,研究方向为电力系统保护与控制、微电网等;
刘志远(1992),女,硕士研究生,主要研究方向为电力系统稳定与控制。基金项目： 国家重点研发计划项目(2016YFB0900601); 国家自然科学基金项目(51377104);
参考文献
[1] 王波,江全元,顾伟,等.一种适于多分支输电线路的故障定位新方法[J].电网技术,2013,37(4):1152-1158. Wang Bo,Jiang Quanyuan,Gu Wei,et al.A novel fault location method suitable for transmission line with multi branches[J].Power System Technology,2013,37(4):1152-1158(in Chinese).
[2] 康小宁,索南加乐.基于参数识别的单端电气量频域法故障测距原理[J].中国电机工程学报,2005,25(2):22-27. Kang Xiaoning,Suonan Jiale.Frequency domain method of fault location based on parameter identification using one-terminal data[J].Proceedings of the CSEE,2005,25(2):22-27(in Chinese).
[3] 田羽,范春菊,龚震东.带同杆双回线的T型线路故障分支判定算法[J].电力系统自动化,2006,30(21):71-76. Tian Yu,Fan Chunju,Gong Zhendong.Faulted line ing method of T circuit of parallel transmission lines[J].Automation of Electric Power Systems,2006,30(21):71-76(in Chinese).
[4] Izykowski J,Rosolowski E,Saha M M,et al.A fault-location method for application with current differential relays of three-terminal lines[J].IEEE Trans on Power Delivery,2007,22(4):2099-2107.
[5] 陈旭,朱永利,赵雪松,等.考虑线路长度变化的T型线路行波测距[J].电网技术,2015,39(5):1438-1443. Chen Xu,Zhu Yongli,Zhao Xuesong,et al.Traveling wave fault location for T-Shaped transmission line considering change of line length[J].Power System Technology,2015,39(5):1438-1443(in Chinese).
[6] 刘鹏辉,黄纯.基于动态时间弯曲距离的小电流接地故障区段定位方法[J].电网技术,2016,40(3):952-957. Liu Penghui,Huang Chun.A fault section location method for small-current grounding fault in distribution network based on dynamic time warping distance[J].Power System Technology,2016,40(3):952-957(in Chinese).
[7] 范新桥,朱永利.基于双端行波原理的多端输电线路故障定位新方法[J].电网技术,2013,37(1):261-269. Fan Xinqiao,Zhu Yongli.A novel fault location scheme for multi-terminal transmission lines based on principle of double-ended traveling wave[J].Power System Technology,2013,37(1):261-269(in Chinese).
[8] 张永健,胥杰,孙嘉.基于静态小波变换的T型输电线路行波测距方法[J].电网技术,2012,36(6):84-88. Zhang Yongjian,Xu Jie,Sun Jia.Travelling wave fault location for three-terminal transmission line based on static wavelet transform[J].Power System Technology,2012,36(6):84-88(in Chinese).
[9] 贾惠彬,赵海锋,方强华,等.基于多端行波的配电网单相接地故障定位方法[J].电力系统自动化,2012,36(2):96-100. Jia Huibin,Zhao Haifeng,FangQianghua,et al.Single-phase earth fault location method for distribution network based on multi-terminal traveling wave[J].Automation of Electric Power Systems,2012,36(2):96-100(in Chinese).
[10] 张广斌,束洪春,于继来,等.不依赖双侧时钟同步的输电线双端行波测距[J].电工技术学报,2015,30(20):199-209. Zhang Guangbin,Shu Hongchun,Yu Jilai,et al.Double-ended travelling wave fault location independent[J].Transactions of China Electro Technical Society,2015,30(20):199-209(in Chinese).
[11] 姜宪国,李博通,张云柯,等.基于电压序量变化量的超高压混合线路故障测距方法[J].电网技术,2015,39(12):3578-3583. Jiang Xianguo,Li Botong,Zhang Yunke,et al.A fault location method for extra-high voltage mixed lines based on variation of sequence voltage[J].Power System Technology,2015,39(12): 3578-3583(in Chinese).
[12] 牛睿,梁军,贠志皓,等.基于故障距离区间的混联线路单端行波测距[J].电网技术,2015,39(1):156-163.Niu Rui,Liang Jun,Yun Zhihao,et al.A fault distance interval based single terminal traveling wave fault location method for hybrid transmission lines[J].Power System Technology,2015,39(1): 156-163(in Chinese).
[13] 雷勇,何晓,周聪聪.一种单端注入多端检测的配电网故障分支判定方法[J].电网技术,2016,40(10):3243-3249. Lei Yong,He Xiao,Zhou Congcong.A distribution network fault branch identification method using single-ended injecting and multi- terminal detecting[J].Power System Technology,2016,40(10): 3243-3249(in Chinese).
[14] 杨荦源,舒勤.基于线模行波突变的配电网全架空线路单相接地的多端检测定位新方法[J].电力系统保护与控制,2016,44(9):21-28. Yang Luoyuan,Shu Qin.A new multiterminal detecting location method of single-phase grounding to the overhead-line for distribution network based on travelling-wave mutation of line-mode[J].Power System Protection and Control,2016,44(9):21-28(in Chinese).
[15] 王波,周昱勇.基于PMU的多端传输线路故障定位新方法[J].电力系统保护与控制,2009,37(12):32-35.Wang Bo,Zhou Yuyong.A novel PMU-based fault location algorithm for multi-terminal transmission line[J].Power System Protection and Control,2009,37(12):32-35(in Chinese).
[16] Hong Yushen,Wen Junzhang,Yu Boduan,et al.A new PMU-based fault location algorithm for three-terminal[J].Advanced Materials Research,2013,2203(634):3925-3929.
[17] 陈旭,朱永利,高艳丰,等.基于故障分支快速辨识的T型高压输电线路故障定位新算法[J].电力系统自动化,2016,40(4):105-111. Chen Xu,Zhu Yongli,Gao Yanfeng,et al.A new fault location algorithm for high-voltage three-terminal transmission lines based on fault branch fast identification[J].Automation of Electric Power Systems,2016,40(4):105-111(in Chinese).
[18] 束洪春,高峰,陈学允,等.T型输电系统故障测距算法研究[J].中国电机工程学报,1998,18(6):28-33. ShuHongchun,Gao Feng,Chen Xueyun,et al.A study on accurate fault location algorithm of EHV T-connection to three terminals[J].Proceedings of the CSEE,1998,18(6):28-33(in Chinese).
[19] 林富洪,王增平,曾惠敏.基于测距函数相位特性的T型高压线路故障定位原理[J].中国电机工程学报,2011,31(13):107-113. Lin Fuhong,Wang Zengping,Zeng Huimin.A novel fault location algorithm based on phase acteristics of fault location function for three-terminal transmission lines[J].Proceedings of the CSEE,2011, 31(13):107-113(in Chinese).
[20] 施世鸿,何奔腾,张武军.T型高压输电线路故障测距[J].中国电机工程学报,2008,28(25):28-32. ShiShihong,He Benteng,Zhang Wujun.Fault location for HV three-terminal transmission lines[J].Proceedings of the CSEE,2008, 28(25):28-32(in Chinese).
