基于零序电流包络线拟合斜率的配电网故障选线方法

㊀㊀
㊀㊀
收稿日期:2021-03-31;修回日期:2021-09-11
基金项目:陕西省工业科技攻关项目(2015G Y 049);西安市科技计划项目(G X Y D 13.10
)通信作者:范文璐(1997-),女,硕士研究生,主要从事配电网小电流接地故障定位研究;E -m a i l :302883468@q q
.c o m 第37卷第3期
电力科学与技术学报
V o l .37N o .3
2022年5月
J O U R N A LO FE I E C T R I CP O W E RS C I E N C EA N DT E C H N O L O G Y
M a y 2
022㊀
基于零序电流包络线拟合斜率的配电网故障选线方法
赵建文,范文璐,胡雨佳
(西安科技大学电气与控制工程学院,陕西西安710054
)摘㊀要:配电网发生高阻接地故障时,传统的选线方法可靠性不高㊂对健全线路与故障线路在高阻接地条件下暂态零序电流的差异性进行分析,发现健全线路与故障线路的暂态零序电流趋势不同,两者对应趋势的斜率正负相反㊂
因此提出了一种基于暂态零序电流包络线拟合斜率的故障选线方法㊂首先对暂态零序电流利用牛顿插值算法求取包络线;然后通过可变遗忘因子递推最小二乘法计算各条线路的拟合斜率,利用健全线路与故障线路包络线拟合斜率值正负相反的特点构建出选线判据;最后大量仿真结果表明:提出的选线方法适用于各种故障条件㊂
关㊀键㊀词:高阻接地故障;故障选线;包络线拟合斜率;牛顿插值法;可变遗忘因子递推最小二乘法
D O I :10.19781/j .i s s n .1673-9140.2022.03.014㊀㊀中图分类号:TM 76㊀㊀文章编号:1673-9140(2022)03-0117-09Ad i s t r i b u t i o nn e t w o r k f a u l t l i n e s e l e c t i o nm e t h o db a s e d o n
z e r o -s e q u e n c e c u r r e n t e n v e l o p e f i t t i n g s l o p
e Z H A OJ i a n w e n ,F A N W e n l u ,HU Y u j
i a (S c h o o l o fE l e c t r i c a l a n dC o n t r o l E n g i n e e r i n g ,X i a nU n i v e r s i t y o f S c i e n c e a n dT e c h n o l o g y
,X i a n710054,C h i n a )A b s t r a c t :W h e nah i g
h -r e s i s t a n c e g r o u n d f a u l t o c c u r s i n t h e d i s t r i b u t i o nn e t w o r k ,t h e t r a d i t i o n a l r o u t e s e l e c t i o nm e t h -o d i sn o t s u f f i c i e n t l y r e l i a b l e .T h e d i f f e r e n c e o f t r a n s i e n t z e r o -s e q u e n c e c u r r e n t u n d e r t h e c o n d i t i o no f h i g hr e s i s t a n c e g r o u n d i n g w h i l e b e t w e e n t h e n o r m a l o p e r a t e d a n d f a u l t l i n e s i s a n a l y z e d .T h em a j o r d i f f e r e n c e i s t h a t t h e s l o p
e s o
f t h e c o r r e s p o n d i n
g t r e n d s o f t
h o s e t w o t y p e s o f c u r r e n t a r e o p p o s
i t e .T h e r e f o r e ,a f a u l t l i n e s e l e c t i o nm e t h o db a s e do n t h e s l o p e o f t r a n s i e n t z e r o -s e q u e n c e c u r r e n t e n v e l o p e c u r v e f i t t i n g i s p r o p o s e d .F i r s t l y ,t h e e n v e l o p
e c u r v e o
f t h e t r a n s i e n t z e r o -s e q u e n c e c u r r e n t i s o b t a i n e d b y t h eN e w t o n i n t e r p o l a t i o n a l
g o r i t
h m ,a n d t h e n t h e f
i t t i n g s l o p e o f e a c h c u r v e i s c a l -c u l a t e db y t h e v a r i a b l e f o r g e t t i n g f a c t o r r e c u r s i v e l e a s t s q u a r em e t h o d ,a n d t h e r e s u l t i n gp o s i t i v e a n dn e g a t i v e c h a r a c -t e r i s t i c s o f t h e s l o p e v a l u e s f i t t e db y t h e e n v e l o p e c u r v e s o f t h e o p e r a t i n g a
n d t h e f a u l t l i n e c a nb e u t i l i z e d t o c o n s t r u c t t h e p r o p o s e d f a u l t l i n e s e l e c t i o n c r i t e r i o n .Al a r g en u m b e r o f s i m u l a t i o nr e s u l t s s h o wt h a t t h e r o u t e s e l e c t i o nm e t h o d p r o p o s e d i n t h i s p a p
e r i s s u i t a b l e
f o r v a r i o u s f a u l t c o n d i t i o n s .K e y w
o r d s :h i g hr e s i s t a n c e g r o u n df a u l t ;f a u l t l i n es e l e c t i o n ;s l o p eo f e n v e l o p ec u r v e ;N e w t o n i n t e r p o l a t i o n ;r e c u r s i v e l e a s t s q u a r em e t h o dw i t ha v a r i a b l e f o r g e t t i n g f
a c t o r Copyright ©博看网. All Rights Reserved.
电㊀㊀力㊀㊀科㊀㊀学㊀㊀与㊀㊀技㊀㊀术㊀㊀学㊀㊀报2022年5月
㊀㊀中国配电网多为小电流接地系统(中性点不接地系统或谐振接地系统)[1]㊂在谐振接地系统中,低阻接地故障的问题已基本解决,而高阻接地故障利用传统的保护方法难以可靠选出故障线路[2]㊂因此研究高阻接地故障选线成为小电流接地系统单相接地故障处理的关键㊂
在高阻接地情况下,故障信号微弱,检测难度大,故障选线较为困难㊂高阻接地故障占配电网故障总次数的10%以上,因此快速可靠地选出高阻接地故障对于电力网安全运行有重要意义㊂针对高阻接地,许多文献给出了有价值的研究,文献[3]提出利用零序电流波形畸变凹凸性的高阻接地故障检测方法,该方法具有较高的灵敏性,但在故障电流很小或噪声很大的情况下可能会失效;文献[4]计算出各条线路故障相的重心频率与相间相关系数作为选线特征量,这种方法在线路高阻接地故障时选线较为可靠,但在母线故障和接地电阻较小时可靠性有待提高;文献[5]提出一种基于电容电流极性及直流分量能量的选线方法,这种方法在小故障角情况下可以准确选线,但在高阻接地情况下,健全线路及故障线路电容电流的极性差异不明显,选线的准确性有待提高;文献[6]利用小波变换提取出故障突变特征量进行选线,这种方法在小故障角条件下有较高的准确率,但高阻条件下各频段的故障信号特征非常微弱,可能会出现误选的情况;文献[7]利用相关性原理选取小波相关系数最小的故障线路,但在高阻接地情况下零序电流振荡幅度很大,波形极其复杂,选线准确性有待提高;文献[8]通过判别零序功率方向来进行选线,这种方法在高阻接地时,零序电流㊁电压极性检测较困难㊂
高阻接地时线路暂态零序电流波形复杂,现有的选线方法存在误选的可能㊂本文提出一种基于暂态零序电流包络线拟合斜率的单相接地故障选线方法,这种方法不受零序电流波形复杂程度的影响㊂此方法引入零序电流信号包络线的概念,包络线反映的即为零序电流的变化趋势,利用牛顿插值法求取包络线,再通过可变遗忘因子递推最小二乘法计算得到包络线拟合斜率,从而构建选线判据㊂此方法对暂态零序电流包络线进行斜率拟合,使得其斜
率特征明显,减小故障特征的辨识难度㊂经仿真验证,该方法可以准确选出故障线路,且不受故障角及过渡电阻等条件的影响㊂
1㊀高阻接地故障零序电流特征分析配电网谐振接地系统单相接地故障等值电路如图1所示[9-11]㊂
i
u0（-1）
01
0J
0J0100（J-1）
C0J
i0f
u f
0f
0Lk
图1㊀配电网谐振接地系统单相高阻接地故障等值电路
F i g u r e1㊀S i n g l e-p h a s e g r o u n d f a u l t e q u i v a l e n t c i r c u i t o f r e s o n a n t g r o u n d i n g s y s t e mi nad i s t r i b u t i o nn e t w o r k
图1中共有J条馈线,u f=U m s i n(ω0t+φ)为故障点的虚拟电源,U m为故障相电压幅值,ω0为工频角频率,φ+π为故障初相角,u0f为母线零序电压,C0j(j=1,2, ,J)为第j条线路的对地零序电容,L k为消弧线圈零序等效电感,R为等效电阻,其等于线路等效电阻与3倍的接地点过渡电阻R f之和,i0j(i C0j)为第j条线路的零序电流(对地零序电容电流),i0f为故障点零序电流,i0L k为流经消弧线圈的零序电流㊂当开关K闭合时为谐振接地系统,打开时为中性点不接地系统㊂
根据图1可知,健全线路与故障线路暂态零序电流方向相反,由此可得流过健全线路零序电流i0j 及故障线路零序电流i0J如下[12]:
i0j=I0e-tτ+L k C0j(δ2A1-A1ω2f-
2δA2ωf)e-δt c o s(ωf t)+L k C0j㊃
(δ2A2-A2ω2f-2δA1ωf)e-δt s i n(ωf t)(1) i0J=-I0e-tτ-[A1+L k(C0Σ-C0J)(δ2A1-A1ω2f-2δA2ωf)]e-δt c o s(ωf t)-[A2+L k(C0Σ-C0J)(δ2A2-A2ω2f+2δA1ωf)]e-δt s i n(ωf t)
(2)
A1=-B s i nφ
A2=-δB s i nφ-ω0B c o sφ
ωf
ì
î
í(3)
811
Copyright©博看网. All Rights Reserved.
第37卷第3期
赵建文,等:基于零序电流包络线拟合斜率的配电网故障选线方法
B =-U m /Z (1-ω20
L k C 0Σ)(4)δ=1/2R C 0Σ
ωf =1/L k C 0Σ-(1/2R C 0Σ)
2
{
(5
)式(4)㊁(5)中㊀Z 为从故障点输入系统阻抗;ωf 为电流自振角频率;C 0ð为所有线路的对地零序电容之和㊂对式(1)㊁(2)求解关于时间t 的一阶导数,可得到健全线路及故障线路暂态零序电流的斜率表达
式,即
i ᶄ0j (t )=-1τ
I 0e -t τ+L k C 0j (
δ2A 1-A 1ω2f -2δA 2ωf )㊃[-δe -δt c o s (ωf t )-ωf e -δt s i n (ωf t )]+L k C 0j (
δ2A 2-A 2ω2
f -2δA 1ωf )㊃[-δe -δt s i n (ωf t )+ωf e -δt c o s (ωf t )](6
)i ᶄ0J (t )=1τ
I 0e -t τ-[A 1+L k (C 0Σ-C 0J )
㊃(δ2A 1-A 1ω2f -2δA 2ωf )][-δe -δt c o s (ωf t )-
ωf e -δt s i n (ωf t )]-[A 2+L k (C 0Σ-C 0J )
㊃(δ2A 2-A 2ω2
f +2
δA 1ωf )]㊃[-δe -δt s i n (ωf t )+ωf
e -δt c o s (ω
f t )](7
)根据式(6)㊁(7)可知,当发生高阻接地时,健全线路与故障线路暂态零序电流的变化趋势相反㊂谐振接地系统故障点处暂态等值电路如图2
所示㊂
u Σ图2㊀谐振接地系统单相接地故障点暂态等值电路F i g
u r e 2㊀T r a n s i e n t e q u i v a l e n t c i r c u i t o f s i n g l e -p h a s e g r o u n d f a u l t p o i n t i na r e s o n a n c e g r o u n d i n g s y
s t e m 当过渡电阻较大时,R 可近似等于3倍的故障
点过渡电阻R f ,此时流过线路的电流较小,电容的充放电速度较慢,主谐振频率相对较低
[13]
㊂
流过故障点的零序电流为
i 0f =-L k
R
(A 2ωf -δA 1)e -δt c o s (ωf t )-L k
R
(-A 1
ωf -δA 2)e -δt s i n (ωf t )(8
)㊀㊀由式(1)㊁(2)㊁(8)可知,发生高阻接地故障后暂态零序电流含有衰减的高频正弦分量㊂由于各种谐振状态的存在,导致暂态零序电流的波形复杂,并且
高阻接地情况下信号微弱,因此特征辨识难度大㊂现有对高阻接地故障的研究用现代信号处理的方法,通常采用幅值㊁极性㊁能量等作为选线特征量,容易受到复杂波形的影响,导致暂态零序电流信号幅值较小;且由于衰减直流分量的影响,存在两者零序电流极性相同的情况,健全线路与故障线路的能量㊁
极性等特征差异不明显[14]
,因此现有方法在高阻接
地故障条件下选线时可能会出现误选的情况㊂
根据式(6)㊁(7)
可知,当发生高阻接地时,健全线路与故障线路暂态零序电流的变化趋势相反,由于故障暂态过程零序电流的数值非常大,且故障发生后的暂态过程非常短,一般仅为2τ~3τ,所以暂态零序电流变化速率很大,其整体趋势对应的斜率数值很大,因此本文先对各条线路的暂态零序电流利用局部极大值的定义进行包络线的求取,包络线反映的即为暂态零序电流的趋势,再对包络线的斜率进行估计,将拟合出的包络线斜率作为选线特征量㊂在高阻接地情况下,此斜率特征显著,不受暂态零序电流波形复杂程度的影响㊂
发生单相接地故障后各条线路的暂态零序电流波形如图3所示㊂当发生高阻单相接地故障时,故障信号较微弱且复杂,但健全线路与故障线路都需要经历暂态衰减过程最终趋于平稳进入稳态过程,且由图3可知,健全线路与故障线路在暂态过程中将从2个不同的方向分别向稳态过程过渡,因此两者暂态零序电流对应趋势的斜率正负相反㊂同理,当过渡电阻较小时,健全线路与故障线路的趋势也
存在同样的特征㊂由上述分析可知发生低阻接地及高阻接地时都存在健全线路与故障线路包络线斜率正负相反的特征,因此可根据这一差异进行单相接地故障选线㊂
健全线路1
健全线路2
健全线路3故障线路
12.510.07.5
5.02.50.0-2.5-5.0-7.5
0.5零序电流/A
1.0 1.5
2.0 2.5
3.0 3.5
4.0 4.5
5.0
0.0时间/（10-3s ）
图3㊀各出线暂态零序电流
F i g
u r e 3㊀T r a n s i e n t z e r o -s e q u e n c e c u r r e n t o f e a c hb r a n c h 9
11Copyright ©博看网. All Rights Reserved.
电㊀㊀力㊀㊀科㊀㊀学㊀㊀与㊀㊀技㊀㊀术㊀㊀学㊀㊀报2022年5月
2㊀基于可变遗忘因子递推最小二乘法的包络线拟合
㊀㊀由于发生单相接地故障后,暂态零序电流的主谐振频率较高,暂态零序电流存在振荡现象,暂态零序电流中每个采样点对应的斜率都不相同,若直接用原始暂态零序电流信号来求取其趋势斜率,则会使得健全线路与故障线路的斜率数值较小,甚至会出现同正或同负的情况,特征差异不明显,因此不能直接利用暂态零序电流原始信号来求取斜率值㊂为了解决这一问题,此处先采用牛顿插值法对各条线路的暂态零序电流利用局部极大值的定义进行包络线的求取,这种算法实际计算时间短,计算量小,求出的插值多项式是唯一的,且不易出现过拟合的现象,得到的包络线与原始数据的趋势更接近,求出的包络线信号频率较低,不存在振荡现象㊂再对包络线利用可变遗忘因子递推最小二乘法进行拟合,此时求得的健全线路与故障线路暂态零序电流信号趋势的斜率正负一定相反,不会出现前述的问题㊂2.1㊀牛顿插值法包络线求取
为获取暂态零序电流的包络线需要先对暂态零序电流采样值进行筛选,即挑选出暂态零序电流的极大值i0m(m=1,2, ,M),其中m表示暂态零序电流极大值的个数,再通过牛顿插值法[15-17]求取暂态零序电流包络线㊂
根据牛顿插值法可得到关于采样时刻及采样值的均差,即
i[t1,t k]=i(t k)-i(t1)
t k-t1
(9)
i[t1,t2, ,t k]=
i[t2,t3, ,t k]-i[t1,t2, ,t k-1]
t k-t1
(10)其中,t1~t k为k个不同的采样时刻,式(9)为零序电流采样值i(t)关于t1㊁t k的一阶均差,式(10)为采样值i的k阶均差,i(t k)为暂态零序电流的极大值i0m㊂对极大值i0m进行筛选,求取每个极值点两侧相邻2个极值的均值,并对该均值设置一定的裕
度,按照中间极值点在该均值裕度范围内保留的原则进行筛选㊂
由式(10)可以写出零序电流包络线的插值公式,即
N(t)=i(t1)+i[t1,t2](t-t1)+ +
i[t1,t2, ,t k,t](t-t1) (t-t k-1)(t-t k)
(11)根据式(11)可求得暂态零序电流包络线㊂
2.2㊀可变遗忘因子递推最小二乘法的包络线拟合
最小二乘算法可以对参数进行估计,其中递推形式的最小二乘算法在有新的数据加入后,不需要对大矩阵进行求逆运算,只需要对原有的估计值进行修正,即可得到新的估计值㊂本文采用可变遗忘因子的递推最小二乘法对参数进行估计,这种算法在数据量较大时可以防止 数据饱和 现象发生,且算法精度较高[18-21]㊂
可变遗忘因子的递推最小二乘法目标函数为
T=ðL n=1λL-n[i n-t T n k]2(12)式中㊀L为总迭代次数;i n为零序电流原始信号值;t n(n=1,2, ,N)为采样时刻;k为待估计的斜率参数;λ为遗忘因子,0<λɤ1㊂根据式(12)可以得到可变遗忘因子的递推最小二乘法参数估计迭代公式为
k n=k n-1+Q n[i n-t T n k n-1](13)
P n=
1
λ
[I-Q n t T n]P n-1(14)
Q n=P n-1t n
λ+t T n P n-1t n
(15)
e n=i n-t n k T n(16)式(16)中㊀e n为误差㊂遗忘因子λ在取值范围内,λ的取值越小,数据遗忘的速度越快,λ越接近1,则算法精度高㊂为了平衡这2种性能,可变遗忘因子λ的取值为
λ[k]=m i nσqσv
ξ+|σe-σv|,λm a x
{}(17)式中㊀ξ为一个很小的常数,防止分母为零;σ2e为误差信号的功率;σ2v为系统噪声功率㊂
因此本文引入特征参量k来表示暂态零序电
021
Copyright©博看网. All Rights Reserved.
第37卷第3期
赵建文,等:基于零序电流包络线拟合斜率的配电网故障选线方法
流包络线趋势的斜率,并利用可变遗忘因子递推最小二乘法对暂态零序电流包络线进行拟合,并估计拟合出的直线斜率,此斜率k 即为包络线的拟合斜率㊂
3㊀基于暂态零序电流包络线拟合斜率的故障选线方法
㊀㊀根据前文分析可知,
健全线路与故障线路暂态零序电流波形复杂,但两者的变化趋势完全不同,其趋势对应的斜率正负相反㊂本文所用的包络线反映的即为暂态零序电流的趋势,因此健全线路与故障线路包络线的斜率正负相反㊂依据这个结论即可构造出适用于高阻接地故障的基于暂态零序电流包络线拟合斜率的小电流接地系统单相接地故障选线方法㊂
当小电流接地系统发生单相接地故障时,提取各条出线1/4周波暂态零序电流i 0j (
j =1,2, ,J )
,并求取其包络线,利用包络线进行拟合斜率的计算,从而消除暂态过程中零序电流振荡现象对拟合的影响,使得利用可变遗忘因子递推最小二乘法求出的健全线路与故障线路暂态零序电流包络线拟合斜率数值k j 更大,差别更显著更易被系统辨识㊂各条出线1/4周波时间窗内的暂态零序电流包络线拟合斜率为
k j n =k j n -1+Q j n i 0j n +i 0j (0)-t T
j
n k j n -1()(18
)式中㊀i 0j (0)为0时刻的零序电流值㊂由各条出线的拟合斜率构造出故障选线向量,即
K =[k 1,k 2, ,k J ]
(19
)㊀㊀由此可得到故障选线判据:
若故障选线向量K 中只有第j 条出线对应的暂态零序电流包络线拟合斜率k j 与其他出线的暂态零序电流包络线拟合斜率正负相反,则该条出线为故障线路;如果所有出线的暂态零序电流包络线拟合斜率k j 正负都相同,则母线发生单相接地故障㊂
故障选线方法的流程如图4所示㊂
求取暂态零序电流包络线
获取暂态零序电流Y
U 0>0.15U m
计算包络线拟合斜率
各条线路拟合斜率
是否存在数值正负相反
输出选线结果
Y 结束
开始
N
母线故障
N 图4㊀故障选线方法流程
F i g
u r e 4㊀T h e f l o wc h a r t o f f a u l t l i n e s e l e c t i o nm e t h o d 4㊀仿真分析及实验验证
4.1㊀仿真模型
利用仿真软件M a t l a b /S i m u l i n k 搭建如图5所
示的10k V 配电网系统仿真模型,共有4条馈出线路,长度分别为L 1=8k m ㊁L 2=12k m ㊁L 3=
16L 4所示㊂E L E s
1
3
图5㊀系统仿真模型
F i g
u r e 5㊀D i a g r a mo f t h e s i m u l a t i o n s y s t e m 表1㊀线路模型参数
T a b l e 1㊀P a r a m e t e r s o f t h e l i n em o d e l
架空线路
电阻/
(Ω/k m )
R 1
R 0
电感/(mH /k m )
L 1
L 0
电容/(μ
F /k m )C 1
C 0
L G J -240
0.17
0.23
1.21
5.48
0.061
0.038
4.2㊀仿真算例
在馈出线L 3距离母线1
0k m 处,设置单相接地故障,系统采用谐振接地过补偿方式运行,补偿度
1
21Copyright ©博看网. All Rights Reserved.
电㊀㊀力㊀㊀科㊀㊀学㊀㊀与㊀㊀技㊀㊀术㊀㊀学㊀㊀报2022年5月
为P =10%,故障点高阻过渡电阻为R g =1000Ω,故障初相角为φ=45ʎ㊂设置各采样点的仿真采样频率为10k H z
㊂各条线路的波形分别如图6~9所示㊂
零序电流包络线
零序电流采样值
0.500.250.00-0.25-0.50
0.5零序电流/A
1.0 1.5
2.0 2.5
3.0 3.5
4.0 4.5
5.0
0.0时间/（10-3s ）
-0.75
图6㊀线路L 1暂态零序电流及包络线
F i g
u r e 6㊀T r a n s i e n t z e r o s e q u e n c e c u r r e n t a n d e n v e l o p e c u r v e o fL 1
零序电流包络线
零序电流采样值
1.00.50.0-0.5-1.0
0.5零序电流/A
1.0 1.5
2.0 2.5
3.0 3.5
4.0 4.5
5.0
0.0时间/（10-3s ）
-1.5
图7㊀线路L 2暂态零序电流及包络线
F i g
u r e 7㊀T r a n s i e n t z e r o s e q u e n c e c u r r e n t a n d e n v e l o p e c u r v e o fL 2
零序电流包络线
零序电流采样值
3210.5零序电流/A
1.0 1.5
2.0 2.5
3.0 3.5
4.0 4.5
5.0
时间/（10-3s ）
图8㊀线路L 3暂态零序电流及包络线
F i g
u r e 8㊀T r a n s i e n t z e r o s e q u e n c e c u r r e n t a n d e n v e l o p e c u r v e o fL 3
零序电流包络线
零序电流采样值
10-10
0.5零序电流/A
1.0 1.5
2.0 2.5
3.0 3.5
4.0 4.5
5.0
时间/（10-3
s ）
图9㊀线路L 4暂态零序电流及包络线
F i g
u r e 9㊀T r a n s i e n t z e r o s e q u e n c e c u r r e n t a n d e n v e l o p e c u r v e o fL 4
图6~9中,虚线为各条出线的暂态零序电流波形,实线为用牛顿插值法求出的零序包络线波形㊂
利用可变遗忘因子递推最小二乘法计算各条馈
出线的包络线拟合斜率,并构造出故障选线向量
K =[46.10,144.64,-213.76,158.97],由选线向量可以看出,线路L 3的拟合斜率与其他线路的拟合斜率正负相反,则线路L 3发生了单相接地故障,
因此该方法可以准确选出高阻接地故障㊂4.3㊀方法适应性分析
为验证此方法能够在各种条件下都适用,尤其是经高阻接地时,此处设置不同的过渡电阻R g ㊁故障角φ㊁系统运行方式(中性点不接地㊁谐振接地欠补偿P =-20%㊁过补偿P =10%)㊁故障位置及加入不同信噪比的噪声进行仿真验证㊂L 3-15为故障点位于线路L 3距母线15k m 处,k j (j =1,2,3,4)为各条线路的暂态零序电流包络线拟合斜率㊂线路及母线故障选线结果如表2~6所示㊂
表2㊀不同过渡电阻故障选线结果
(φ=3
0ʎ,P =10%,L 3-10)T a b l e 2㊀F a u l t l i n e s e l e c t i o n r e s u l t s o f d i f f e r e n t t r a n s i t i o n r e s i s t a n c e s (φ=3
0ʎ,P =10%,L 3-10)R g /ΩK =[k 1,k 2,k 3,k 4]
选线结果㊀100
[604.56,1361.70,-2385.65,721.61]L 3500[46.78,176.11,-299.76,87.98]L 31000[69.09,98.58,-213.01,168.42]L 32000[41.61,43.15,-88.54,83.62]L 33000[36.98,33.04,-60.67,53.97]L 34000
[10.69,26.17,-37.10,39.73
]L 3
表3㊀不同故障角故障选线结果
(R g =1000Ω,P =10%,L 3-10)
T a b l e 3㊀F a u l t l i n e s e l e c t i o n r e s u l t s o f d i f f e r e n t f a u l t a n g l e s (R g =1000Ω,P =10%,L 3-10)
φ/(ʎ)K =[k 1,k 2,k 3,k 4]
选线结果0[-23.86,-25.82,56.87,-58.97]L 330[43.05,98.58,-139.92,104.42
]L 390
[45.30,150.24,-538.69,122.87]L 3
221Copyright ©博看网. All Rights Reserved.
第37卷第3期赵建文,等:基于零序电流包络线拟合斜率的配电网故障选线方法表4㊀不同运行方式故障选线结果
(φ=30ʎ,R g=1000Ω,L3-10)
T a b l e4㊀F a u l t l i n e s e l e c t i o n r e s u l t s o f d i f f e r e n t s y s t e m
o p e r a t i o nm o d e s(φ=30ʎ,R g=1000Ω,L3-10)
系统运行方式K=[k1,k2,k3,k4]选线结果
谐振接地过补偿[43.05,98.58,-139.92,104.42]L3
谐振接地欠补偿[40.45,90.36,-168.82,105.00]L3
不接地[16.23,42.11,-320.96,67.19]L3
表5㊀不同故障位置故障选线结果(φ=30ʎ,P=10%)
T a b l e5㊀F a u l t l i n e s e l e c t i o n r e s u l t s o f d i f f e r e n t
f a u l t l o c a t i o n s(φ=30ʎ,P=10%)
R g/Ω故障
位置
K=[k1,k2,k3,k4]选线结果
1000L3-10[43.05,98.58,-139.92,104.42]L3 1000L3-15[63.21,80.89,-106.66,139.52]L3 1000母线[50.83,74.78,110.08,83.44]母线表6㊀不同信噪比故障选线结果(φ=30ʎ,
R g=1000Ω,P=10%,L3-10)
T a b l e6㊀F a u l t l i n e s e l e c t i o n r e s u l t s o f d i f f e r e n t s i g n a l-t o-
n o i s e r a t i o s(φ=30ʎ,R g=1000Ω,P=10%,L3-10)
信噪比K=[k1,k2,k3,k4]选线结果
10[34.52,260.96,-154.61,7.72]L3
20[5.70,365.31,-162.18,97.48]L3
30[24.26,201.01,-158.94,52.74]L3
由表2~5可知,当线路发生单相接地故障时,健全线路与故障线路的包络线拟合斜率数值正负相反,且此数值的正负不受过渡电阻㊁故障角㊁运行方式及故障位置的影响,在最不利的故障条件下(故障角φ=0ʎ,高阻R g=4000Ω),健全线路与故障线路包络线拟合斜率正负仍相反㊂表6为加入不同信噪比噪声的仿真结果,由表6数据可看出,当原始信号受到噪声干扰时,选线结果不会受到影响㊂由此可见本方法可以准确选出故障线路㊂当母线发生单相高阻接地故障时,各条出线的包络线拟合斜率正负相同,由此可见本方法可以准确判别出母线故障㊂4.4㊀实验验证
将本文所提出的选线方法在具有4条馈出线的配电网模拟系统中进行实验验证,在不同接地电阻
条件下可准确选出故障线路㊂例如当模拟系统采用谐振接地过补偿10%,设置线路L1发生单相接地
故障,过渡电阻为1000Ω㊂通过计算得到的包络线拟合斜率向量为K=[-58.97,13.75,6.57,5.66],根据本文选线判据,故障线路L1与其余3条健全线路包络线拟合斜率正负相反,具有明显的故障特征㊂5㊀结语
本文提出了一种基于暂态零序电流包络线拟合斜率的单相高阻接地故障选线方法㊂
1)该方法利用发生单相高阻接地故障后,健全线路与故障线路暂态零序电流趋势不同,其对应的斜率正负相反的特点进行故障选线,在不同故障条件下选线的结果不受影响㊂
2)由于高阻接地条件下,故障特征不明显,因此本文采用包络线拟合的方法使得拟合斜率值更大,从而使得高阻接地条件下健全线路和故障线路的特征差异更显著,减小了故障特征的辨识难度㊂
参考文献:
[1]程路,陈乔夫.小电流接地系统单相接地选线技术综述
[J].电网技术,2009,33(18):219-224.
C H E N GL u,C H E N Q i a o f u.As u r v e y o n f a u l t y l i n e s e-l e c t i o nt e c h n o l o g y f o rs i n g l e-p h a s e g r o u n d e dt r a n s m i s-s i o n l i n e i ns m a l l c u r r e n tn e u t r a l g r o u n d e ds y s t e m[J].
P o w e r S y s t e m T e c h n o l o g y,2009,33(18):219-224. [2]方毅,薛永端,宋华茂,等.谐振接地系统高阻接地故障
暂态能量分析与选线[J].中国电机工程学报,2018,38
(19):5636-5645+5921.
F A N
G Y i,X U EY o n g d u a n,S O N G
H u a m a o,e t a l.T r a n-s i e n te n e r g y a n a l y s i s a n d f a u l t y f e e d e ri d e n t i f i c a t i o n m e t h o do f h i g h i m p e d a n c e f a u l t i n t h e r e s o n a n t g r o u n d-i n g s y s t e m[J].P r o c e e d i n g s o f t h eC S E E,2018,38(19): 5636-5645+5921.
[3]耿建昭,王宾,董新洲,等.中性点有效接地配电网高阻
接地故障特征分析及检测[J].电力系统自动化,2013, 37(16):85-91.
G E N GJ i a n z h a o,WA N G B i n,D O N G X i n z h o u,e t a l.A-
n a l y s i s a n dd e t e c t i o no f h i g h i m p e d a n c e g r o u n d i n g f a u l t
321
Copyright©博看网. All Rights Reserved.
电㊀㊀力㊀㊀科㊀㊀学㊀㊀与㊀㊀技㊀㊀术㊀㊀学㊀㊀报2022年5月
i nn e u t r a l p o i n te f f e c t i v e l yg r o u n d i n g d i s t r i b u t i o nn e t-w o r k[J].A u t o m a t i o no fE l e c t r i cP o w e rS y s t e m s,2013, 37(16):85-91.
[4]刘谋海,王媛媛,曾祥君,等.基于暂态相电流特征分析的故障选线新方法[J].电力系统及其自动化学报, 2017,29(1):30-36.
L I U M o u h a i,WA N G Y u a n y u a n,Z E N G X i a n g j u n,e t a l. N o v e lm e t h o do f f a u l t l i n e s e l e c t i o nb a s e do nc h a r a c t e r-i s t i c a n a l y s i so f t r a n s i e n t p h a s e-c u r r e n t[J].P r o c e e d i n g s o f t h eC S U-E P S A,2017,29(1):30-36.
[5]魏科文,张靖,何宇,等.基于VM D和相关性聚类的谐振接地系统单相接地故障选线[J].电力系统保护与控制,2021,49(22):105-113.
W E IK e w e n,Z HA N GJ i n g,H E Y u,e t a l.S i n g l e-p h a s e g r o u n d i n g f a u l tl i n es e l e c t i o ni nar e s o n a n t g r o u n d i n g s y s t e mb a s e do n VM D a n dc o r r e l a t i o nc l u s t e r i n g[J]. P o w e r s y s t e m p r o t e c t i o na n dc o n t r o l,2021,49(22): 105-113.
[6]王栋,邱志斌,魏巍,等.220k V某变电站G I S单相接地故障分析及处理措施[J].高压电器,2020,56(11):259-265.
WA N GD o n g,Q I UZ h i b i n,W E IW e i,e t a l.S i n g l e-p h a s e G r o u n d i n g F a u l tA n a l y s i s o f aG I S i na220k VS u b s t a-t i o n a n dT r e a t m e n tM e a s u r e s[J].H i g hV o l t a g eA p p a r a-t u s,2020,56(11):259-265.
[7]魏向向,温渤婴.基于2阶累加生成相关性的谐振接地系统故障选线方法[J].电网技术,2017,41(5):1674-1682.
W E IX i a n g x i a n g,W E NB o y i n g.An o v e l f a u l t l i n e d e t e c-t i o n m e t h o d b a s e do n2-o r d e ra c c u m u l a t e d g e n e r a t i n g o p e r a t i o nc o r r e l a t i o na n a l y s i s f o rr e s o n a n te a r t h e ds y s-t e m[J].P o w e rS y s t e m T e c h n o l o g y,2017,41(5):1674-1682.
[8]郑康霖,耿增玮,王红霞,等.基于同步量测大数据的配网单相接地故障检测及定位方法[J].电网与清洁能源,2020,36(9):50-56.
Z H E N G K a n g l i n,G E N G Z e n g w e i,WA N G H o n g x i a,e t a l.S i n g l e-p h a s e g r o u n d i n g f a u l td e t e c t i o na n dl o c a l i z a-t i o nb a s e do nb i g d a t a f r o ms y n c h r o n i z e d p h a s o rm e a s-u r e m e n t[J].P o w e rS y s t e m a n dC l e a nE n e r g y,2020,36 (9):50-56.
[9]杨帆,刘鑫星,沈煜,等.基于零序电流投影系数的小电
阻接地系统高阻接地故障保护[J].电网技术,2020,44 (3):1128-1133.
Y A N GF a n,L I U X i n x i n g,S H E N Y u,e t a l.H i g hr e s i s t-a n c e g r o u n d f a u l t p r o t e c t i o n o f l o wr e s i s t a n c e g r o u n d i n g s y s t e mb a s e do n z e r o s e q u e n c e c u r r e n t p r o j e c t i o n c o e f f i-c i e n t[J].P o w e rS y s t e m T e c h n o l o g y,2020,44(3):1128-1133.
[10]高俊青,李题印,胡晓琴,等.考虑频段筛选的配电网单
相接地选线方法研究[J].智慧电力,2020,48(2):92-97+118.
G A OJ u n q i n g,L IT i y i n,HU X i a o q i n,e t a l.L i n ed e t e c-
t i o nm e t h o d s f o r s i n g l e-p h a s e g r o u n d i n g f a u l t i nd i s t r i-
b u t i o nn e t w o r k
c o n s i
d
e r i n g
f r e q u e n c y b a n ds e l e c t i n g
[J].S m a r tP o w e r,2020,48(2):92-97+118. [11]P A U L P D.P h a s e-g r o u n df a u l tc u r r e n ta n a l y s i sa n d
p r o t e c t i o no f a h i g h-r e s i s t a n c e g r o u n d e d p o w e r s y s t e m
[J].I E E E T r a n s a c t i o n s o n I n d u s t r y A p p l i c a t i o n s, 2020,56(4):3306-3314.
[12]束洪春,龚振,田鑫萃,等.基于故障特征频带及形态谱
的单相接地故障选线[J].电网技术,2019,43(3): 1041-1053.
S HU H o n g c h u n,G O N GZ h e n,T I A NX i n c u i,e t a l.S i n-
g l el i n e-t o-g r o u n df a u l tl i n es e l e c t i o n b a s e d o nf a u l t
c h a r a c t e r i s t i c f r e q u e n c y b a n da n
d m o r p h o l o g i c a l s p
e c-
t r u m[J].P o w e r S y s t e m T e c h n o l o g y,2019,43(3): 1041-1053.
[13]王庭华,薛永端,赵雪霖,等.不同接地电阻下消弧线圈
接地系统接地故障暂态特征[J].电力系统及其自动化
学报,2019,31(9):43-51.
WA N G T i n g h u a,X U E Y o n g d u a n,Z HA O X u e l i n,e t
a l.T r a n s i e n tc h a r a c t e r i s t i c so fe a r t hf a u l t i na r cs u p-
p r e s s i o n c o i l g r o u n d e d s y s t e m w i t hd i f f e r e n t g r o u n d i n g r e s i s t a n c e s[J].P r o c e e d i n g so ft h eC S U-E P S A,2019, 31(9):43-51.
[14]李福志,郑卫宾,张文海.基于回路直流电阻测量的输
电线路单相接地故障离线故障定位[J].中国电力, 2021,54(2):140-146.
L IF u z h i,Z H E N G W e i b i n,Z HA N G W e n h a i.F a u l t p a t h d i r e c t-c u r r e n t r e s i s t a n c e b a s e d o f f-l i n e s i n g l e-p h a s e-t o-g r o u n df a u l tl o c a t i o n[J].E l e c t r i c P o w e r, 2021,54(2):140-146.
[15]Z HA N G Y,S HA N G C.C o m b i n i n g n e w t o ni n t e r p o l a-
421
Copyright©博看网. All Rights Reserved.
第37卷第3期赵建文,等:基于零序电流包络线拟合斜率的配电网故障选线方法
t i o n a n d d e e p l e a r n i n g f o r i m a g e c l a s s i f i c a t i o n[J].E l e c-t r o n i c sL e t t e r s,2015,51(1):40-42.
[16]宋智威,熊成林,黄路,等.基于牛顿插值的单相整流器
功率前馈无差拍控制[J].电网技术,2018,42(11): 3623-3629.
S O N GZ h i w e i,X I O N G C h e n g l i n,HU A N G L u,e ta l.
P o w e r f e e d b a c k-f o r w a r d a n d d e a d b e a t c o n t r o l o f s i n g l e-p h a s e r e c t i f i e r b a s e do nn e w t o n i n t e r p o l a t i o n[J].P o w-
e r S y s t e m T e c h n o l o g y,2018,42(11):3623-3629.
[17]盛立锃,曾喆昭,李莎.基于代数多项式模型的用电量
预测研究[J].电力科学与技术学报,2015,30(1):34-
40.
S H E N G L i z e n g,Z E N G Z h e z h a o,L IS h a.R e s e a r c ho n
e l e c t r i c i t y c o n s u m p t i o n p r e d i c t i o n w i t ht h ea l g e b r a i c
p o l y n o m i a lm o d e l[J].J o u r n a l o f E l e c t r i c P o w e r S c i e n c e
a n dT e c h n o l o g y,2015,30(1):34-40.
[18]B E Z A M,B O N G I O R N O M.A p p l i c a t i o no fr e c u r s i v e
l e a s t s q u a r e s a l g o r i t h m w i t hv a r i a b l e f o r g e t t i n g f a c t o r
f o r f r e q u e n c y c o m p o n e n t e s t i m a t i o n i na
g e n e r i c i n p u t
s i g n a l[J].I E E E T r a n s a c t i o n so nI n d u s t r y A p p l i c a-t i o n s,2014,50(2):1168-1176.
[19]尹良震,李奇,洪志湖,等.P E M F C发电系统F F R L S
在线辨识和实时最优温度广义预测控制方法[J].中国
电机工程学报,2017,37(11):3223-3235+3378.
Y I N L i a n g z h e n,L IQ i,HO N GZ h i h u,e t a l.F F R L So n-l i n ei d e n t i f i c a t i o n a n dr e a l-t i m eo p t i m a lt e m p e r a t u r e
g e n e r a l i z e d p r e d i c t i v e c o n t r o lm e t h o do fP E M F C p o w-
e r g e n e r a t i o ns y s t e m[J].P r o c e e d i n g so ft h e C S E E,
2017,37(11):3223-3235+3378.
[20]范忠,史明明,张宸宇,等.基于三点筛选法与偏最小
二乘法的系统谐波阻抗估计[J].电测与仪表,2020, 57(17):34-38.
F A NZ h o n g,S H IM i n g m i n g,Z HA N GC h e n y u,e t a l.
E s t i m a t i o no fs y s t e m h a r m o n i ci m p e d a n c e b a s e d o n
t h r e e-p o i n t s c r e e n i n g a n d p a r t i a l l e a s t s q u a r e s[J].E-l e c t r i c a l M e a s u r e m e n t&I n s t r u m e n t a t i o n,2020,57
(17):34-38.
[21]荀倩,王培良,李祖欣,等.基于递推最小二乘法的永磁
伺服系统参数辨识[J].电工技术学报,2016,31(17): 161-169.
X U N Q i a n,WA N G P e i l i a n g,L IZ u x i n,e ta l.P M S M
p a r a m e t e r s i d e n t i f i c a t i o n b a s e d o n r e c u r s i v e l e a s t s q u a r em e t h o d[J].T r a n s a c t i o n so fC h i n aE l e c t r o t e c h-n i c a l S o c i e t y,2016,31(17):161-169.
(上接第108页㊀C o n t i n u e d f r o m p a g e108) [12]刘幸蔚,李永丽.逆变型电源T接电网后对纵差保护
的影响及准入容量分析[J].电网技术,2016,40(5): 1595-1600.
L I U X i n g w e i,L IY o n g l i.E f f e c to f I I D Gc o n n e c t e dt o
g r i da s a t e e d l i n e o n l o n g i t u d i n a l d i f f e r e n t i a l p r o t e c t i o n
a n d m a x i m u m p e n e t r a t i o nc a p a c i t y[J].P o w e rS y s t e m
T e c h n o l o g y,2016,40(5):1595-1600.
[13]黄河,朱磊,高松,等.提升分布式电源消纳的配网快速
遍历重构方法[J].电力科学与技术学报,2019,34(3): 166-172.
HU A N G H e,Z HU L e i,G A OS o n g,e t a l.R e c o n f i g u r a-t i o n m e t h o do fd i s t r i b u t i o ns y s t e m f o r i n c r e a s i n g t h e p e n e t r a t i o no f d i s t r i b u t e d g e n e r a t i o n[J].J o u r n a l o fE-l e c t r i cP o w e rS c i e n c ea n d T e c h n o l o g y,2019,34(3): 166-172.[14]郭小强,邬伟扬,漆汉宏.电网电压畸变不平衡情况下
三相并网逆变器控制策略[J].中国电机工程学报, 2013,33(3):22-28+19.
G U O X i a o q i a n g,WU W e i y a n g,Q I H a n h o n g.C o n t r o l
s t r a t e g i e so ft h r e e-p h a s e P V g r i d-c o n n e c t e di n v e r t e r u n d e r d i s t o r t e da n du n b a l a n c e dv o l t a g e c o n d i t i o n s[J].
P r o c e e d i n g s o f t h eC S E E,2013,33(3):22-28+19.
[15]韩伟,刘素梅,肖超,等.计及发电单元故障特性差异的
光伏电站短路电流计算模型[J].电力系统保护与控
制,2021,49(20):129-135.
HA N W e i,L I US u m e i,X I A OC h a o,e t a l.S t e a d y s h o r t-
c i r c u i t c a l c u l a t i o nm o
d
e l
f o r a p h o t o v o l t a i c s t a t i o n c o n-
s i d e r i n g d i f f e r e n t g e n e r a t i o nu n i t s'f a u l t c h a r a c t e r i s t i c s
[J].P o w e rS y s t e m P r o t e c t i o na n d C o n t r o l,2021,49
(20):129-135.
521
Copyright©博看网. All Rights Reserved.。