CN109342888B - 一种基于重合闸动作分析的单端输电线路故障测距方法 - Google Patents

一种基于重合闸动作分析的单端输电线路故障测距方法 Download PDF

Info

Publication number
CN109342888B
CN109342888B CN201811407428.2A CN201811407428A CN109342888B CN 109342888 B CN109342888 B CN 109342888B CN 201811407428 A CN201811407428 A CN 201811407428A CN 109342888 B CN109342888 B CN 109342888B
Authority
CN
China
Prior art keywords
fault
state
reclosing
impedance
current
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811407428.2A
Other languages
English (en)
Other versions
CN109342888A (zh
Inventor
季亮
陶小杰
施佳斌
曹喆
倪泽宇
米阳
杨兴武
李振坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN201811407428.2A priority Critical patent/CN109342888B/zh
Publication of CN109342888A publication Critical patent/CN109342888A/zh
Application granted granted Critical
Publication of CN109342888B publication Critical patent/CN109342888B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/085Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Locating Faults (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

本发明公开了一种基于重合闸动作分析的单端输电线路故障测距方法。该方法根据电路叠加原理,分别分析故障状态和故障后自动重合闸状态下的附加状态网络,建立正序故障分量状态下本侧测量电压电流关于故障距离、过渡电阻及对侧系统阻抗的方程组。最后,利用信赖域反射算法求解故障状态及重合闸状态的方程组,获得精确的故障距离。本发明完全克服了传统单端法因可利用信息有限而无法同时消除故障距离、过渡电阻及对侧系统阻抗对测距精度影响的弊端。

Description

一种基于重合闸动作分析的单端输电线路故障测距方法
技术领域
本发明涉及电力系统继电保护技术领域,尤其是涉及一种基于重合闸动作分析的单端输电线路故障测距方法。
背景技术
快速、准确的故障定位技术对于修复故障,迅速恢复供电,提高电力系统可靠性和连续性有着至关重要的意义。根据测距原理,故障测距主要分为行波法及故障分析法。
行波法根据故障后行波到达测量点的时间差测算故障距离,原理上具有很高的精度。但波头及波速的检测存在一定难度,且设备成本过高。故障分析法利用系统故障前后的稳态电压电流等信息建立方程,求解故障距离。故障分析法根据信息来源,可分为单端法和双端法。双端故障分析法利用两边信息进行分析,可消除过渡电阻等因素的影响,拥有较高测距精度。此类方法的可靠性受线路两侧信号同步性影响,且经济性较差,同时,大部分双端法存在伪根判别的问题且无法在所有情况下都能判别出伪根。单端故障分析法通过分析单端工频电压与电流之间的关系方程,获取故障距离。传统的单端法测距精度不高,易受到负荷水平、实际故障距离、过渡电阻、对侧系统运行条件变化等因素的影响。究其根本原因,传统的单端法由于利用的测量信息有限,未知量个数大于方程个数,无法同时消除实际故障距离、过渡电阻及对侧系统运行条件变化的影响。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种为解决单端法可利用测量信息有限的问题,将分析多个时间层面(故障状态和重合闸状态)以提取更多有效信息,采用分布式线路参数,建立故障状态和单相重合闸状态下的单端电气量与故障距离、过渡电阻及对侧系统阻抗的方程组,通过信赖域反射算法求解方程组,获取精确的故障距离,应用方便,仅通过变电站中故障录波仪及线路参数即可应用于现有的输电系统的一种基于重合闸动作分析的单端输电线路故障测距方法。
本发明的目的可以通过以下技术方案来实现:
一种基于重合闸动作分析的单端输电线路故障测距方法,包括以下步骤:
步骤1:将线路发生永久性故障前后整个过程分为5个时间面状态;
步骤2:针对5个时间面状态中的两个时间面状态分别建模并分别建立包含未知故障距离值和其他相关未知数的关系方程组;
步骤3:对包含未知故障距离值和其他相关未知数的关系方程组进行求解,得出故障距离值;
步骤4:将故障距离值进行循环迭代,以满足误差的设定值以内并得出最终故障距离值。
进一步地,所述步骤1中的5个时间面状态包括:状态A:故障前状态,线路两侧断路器闭合、状态B:故障状态,线路发生故障,但线路两侧断路器仍然闭合、状态C:故障清除状态,保护动作后,线路两侧断路器断开、状态D:重合闸状态,线路一侧断路器先行闭合,另外一侧断路器保持断开和状态E:最终隔离状态,故障未能清除,保护动作,再次断开合闸侧的断路器。
进一步地,所述步骤2中的两个时间面状态为状态B:故障状态,线路发生故障,但线路两侧断路器仍然闭合和状态D:重合闸状态,线路一侧断路器先行闭合,另外一侧断路器保持断开。
进一步地,所述的状态B的关系方程组为:
Figure BDA0001877696670000021
Figure BDA0001877696670000022
式中,
Figure BDA0001877696670000023
Figure BDA0001877696670000024
分别为正序故障附加状态网络中本侧电压和电流,Vm和Im分别为本侧电压和电流,AL(1),BL(1),Dm(1),Bm(1),CL(1),DL(1)为本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL的正序分量,S=1/3[1 α α2]且α=ej120°,本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL和对侧状态网络模型参数矩阵AR,BR,Dn,Bn,CR,DR中:Bm=ALZm+BL,Dm=CLZm+DL,Bn=ARZn+BR,Dn=CRZn+DR,AL=DL=coshγx,BL=Zc·sinhγx,
Figure BDA0001877696670000025
AR=DR=cosh(γ(d-x)),BR=Zc·sinh(γ(d-x)),
Figure BDA00018776966700000312
Zm和Zn分别为本侧系统阻抗和对侧系统阻抗,Zc和γ分别为波阻抗和传播系数,d为输电线路总长度,x表示故障点到本侧测量点的距离,
Figure BDA0001877696670000031
Zf为过渡阻抗。
进一步地,所述的状态D的关系方程组为:
Figure BDA0001877696670000032
Figure BDA0001877696670000033
Figure BDA0001877696670000034
Figure BDA0001877696670000035
Figure BDA0001877696670000036
式中,Xij,Yij(i,j=1,2,3)分别为X,Y参数矩阵中的元素,矩阵
Figure BDA0001877696670000037
矩阵
Figure BDA0001877696670000038
Zms和Zmm分别为本侧系统的自阻抗和互阻抗,
Figure BDA0001877696670000039
Figure BDA00018776966700000310
分别为正序重合闸附加状态网络中本侧电压和电流,
Figure BDA00018776966700000311
为重合闸附加状态网络中的相电流。
进一步地,所述步骤3的包含未知故障距离值和其他相关未知数的关系方程组的求解方法采用信赖域反射求解方法。
进一步地,所述步骤4中的设定值为0.1%。
与现有技术相比,本发明具有以下优点:
(1)适用性强,本发明的线路模型为分布参数线路模型,不受分布电容电流的影响,适用于高压长距离输电线路。
(2)精确度高,本发明通过分别分析故障状态和故障后自动重合闸状态下的附加状态网络,建立正序故障分量状态下关于故障距离、过渡电阻及对侧系统阻抗的方程组,并利用信赖域反射算法计算出精确的故障位置。
(3)抗干扰性强,本发明完全克服了传统单端法因可利用信息有限而无法同时消除故障距离、过渡电阻及对侧系统阻抗对测距精度影响的弊端,同时,可通过信赖域反射算法量化这三种影响因素,根据EMTP-ATP的仿真验证结果,最终的测距精度很高,并且测距结果完全不受故障距离、过渡电阻及对侧系统阻抗的影响。
附图说明
图1为永久性故障时单相重合闸过程中的电流波形示意图,其中图1(a)为单相重合闸中本侧断路器先合闸情况的电流波形示意图,图1(b)为单相重合闸中对侧断路器先合闸情况的电流波形示意图;
图2为故障前网络示意图;
图3为故障附加状态网络(故障分量网络)示意图;
图4为故障后网络示意图;
图5为正序故障附加状态网络示意图;
图6为重合闸时附加状态网络示意图;
图7为重合闸时网络示意图;
图8为正序重合闸时附加状态网络示意图;
图9为输电系统模型示意图;
图10为本发明测距算法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
本发明的原理分析:
(1)重合闸动作过程分析
输电线路故障主要分为瞬时故障和永久性故障。瞬时性故障时,可通过自动重合闸重合断路器,迅速恢复供电,提高供电可靠性及连续性。永久性故障时,自动重合闸在重合失败后,将再次断开断路器,等待检修人员排除故障。
根据统计,在输电线路中,单相接地故障发生概率可高达80%以上。因此,本发明以单相接地故障和单相重合闸为对象,展开研究。在单相重合闸中,由于只断开故障相,原则上线路两测断路器不需要检同期而同时进行合闸。但这种做法使两边断路器同时再次受到故障的冲击。有文献提出了一种改进的单相重合闸方式,借用三相重合闸在永久性故障下,检无压后不再重合对侧断路器的原理,减少故障对一侧断路器的冲击。此外,现场中线路两侧由于设备机械反映时间、保护动作、保护返回时间不同等原因,也很难达到断路器同时合闸。因此,在改进的单相重合闸以及现场实际的单相重合闸动作情况下,断路器在重合闸的过程中,将会出现一侧故障相断路器闭合,另一侧故障相断路器断开的系统状态,在本发明中称为重合闸状态。
重合闸在重合过程中,线路两侧的断路器都可能先合闸。因此,将重合闸状态分为两种,a)本侧断路器先合闸b)对侧断路器先合闸。如图1所示表示在单相重合闸情况下,两种不同重合状态下的电流波形。在自动重合闸策略下,线路发生永久性故障后可分为5个时间面:
-状态A:故障前状态,线路两侧断路器闭合
-状态B:故障状态,线路发生故障,但线路两侧断路器仍然闭合
-状态C:故障清除状态,保护动作后,线路两侧断路器断开
-状态D:重合闸状态,线路一侧断路器先行闭合,另外一侧断路器保持断开
-状态E:最终隔离状态,永久性故障未能清除,保护动作,再次断开合闸侧的断路器。
在重合闸的5个状态中,状态B、D直接与故障距离、过渡电阻及对侧系统阻抗相关。传统单端故障测距方法一般只用到状态B或者A、B,此类方法无法提供足够的已知量去求解未知量。本发明提出来的基于重合闸动作过程的测距方法将同时利用状态B及状态D,建立多状态系统方程组,求解故障距离。本发明中,将状态B和D称为故障状态和重合闸状态。
(2)根据叠加电路原理
故障后网络(图4)可以表示为故障前网络(图2)叠加故障附加状态网络(图3)。显然,图4中的正序分量也可由图3和图2中的正序分量叠加得到。本部分将先利用图2和图3求出其中各自电压电流正序分量,再建立故障后网络正序分量方程,以简化模型。
图2、图3和图4中,根据输电线路分布参数的模型特性,可知:AL=DL=coshγx,BL=Zc·sinhγx,
Figure BDA0001877696670000061
AR=DR=cosh(γ(d-x)),BR=Zc·sinh(γ(d-x)),
Figure BDA0001877696670000062
其中,Zm和Zn分别为本侧系统阻抗和对侧系统阻抗,Zc和γ分别为波阻抗和传播系数,d为输电线路总长度,x表示故障点到本侧测量点的距离,Ef表示故障前故障点电压,图3中
Figure BDA0001877696670000063
分别表示故障后故障点左侧和右侧电流,图3及图4中变量的上标B表示状态B(故障状态),图中Zf为包含过渡电阻Rf的矩阵,Zf=RfK,发生A相短路故障时,
Figure BDA0001877696670000064
对于故障附加状态网络,如图3所示,有:
Figure BDA0001877696670000065
Figure BDA0001877696670000066
其中,Bm=ALZm+BL,Dm=CLZm+DL,Bn=ARZn+BR,Dn=CRZn+DR,Am=AL,Cm=CL,An=AR,Cn=CR,AL,BL,Dm,Bm,CL,DL,Am,Cm为本侧状态网络模型参数矩阵,AL,BL,Dm,Bm,CL,DL,An,Cn为对侧状态网络模型参数矩阵,
Figure BDA0001877696670000067
为故障后故障状态模型矩阵变量,
Figure BDA0001877696670000068
分别为故障后本侧和对侧故障状态网络模型电流;
由(1-1)和(1-2)可以推得:
Figure BDA0001877696670000069
Figure BDA00018776966700000610
在故障点对地支路,可推得故障后故障点电流
Figure BDA00018776966700000611
有如下的等式:
Figure BDA00018776966700000612
Figure BDA00018776966700000613
其中,
Figure BDA0001877696670000071
Ef可由故障前状态推得,有:
Ef=ALVm-BLIm (1-7)
Vm和Im分别为本侧电压和电流,由(1-3)-(1-7)推得故障后故障点电压
Figure BDA0001877696670000072
Figure BDA0001877696670000073
令S=1/3[1 α α2]且α=ej120°,则故障后故障点的正序电压
Figure BDA0001877696670000074
为:
Figure BDA0001877696670000075
式中,下标(1)表示正序分量。
图5即为正序故障附加状态网络,对于正序分量,其同样满足(1-1)-(1-4),故可推得正序故障附加状态网络中左侧故障电流
Figure BDA0001877696670000076
与正序故障附加状态网络中故障点电压
Figure BDA0001877696670000077
的关系:
Figure BDA0001877696670000078
式中,Dm(1)和Bm(1)为本侧状态网络模型参数矩阵Dm,Bm的正序分量,推得
Figure BDA0001877696670000079
的表达式后,正序故障附加状态网络中本侧电压和电流也可推得:
Figure BDA00018776966700000710
将(1-10)代入(1-11)便可得到正序故障附加状态的本侧电压电流,综上所述,如图5所示,最终确定的对于状态B的关系方程组为:
Figure BDA00018776966700000711
Figure BDA00018776966700000712
式中,
Figure BDA00018776966700000713
Figure BDA00018776966700000714
分别为正序故障附加状态网络中本侧电压和电流,Vm和Im分别为本侧电压和电流,AL(1),BL(1),Dm(1),Bm(1),CL(1),DL(1)为本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL的正序分量,S=1/3[1 α α2]且α=ej120°,本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL和对侧状态网络模型参数矩阵AR,BR,Dn,Bn,CR,DR中:Bm=ALZm+BL,Dm=CLZm+DL,Bn=ARZn+BR,Dn=CRZn+DR,AL=DL=coshγx,BL=Zc·sinhγx,
Figure BDA00018776966700000715
A=D=cosh(γ(d-x)),BR=Zc·sinh(γ(d-x)),
Figure BDA00018776966700000716
Zm和Zn分别为本侧系统阻抗和对侧系统阻抗,Zc和γ分别为波阻抗和传播系数,d为输电线路总长度,x表示故障点到本侧测量点的距离,
Figure BDA00018776966700000717
Zf为过渡阻抗。
(3)重合闸状态系统建模
重合闸状态根据线路两侧断路器先后重合的顺序,可分为本侧断路器先闭合和对侧断路器先闭合两种情况。两种情况的分析原理相同,以对侧断路器先进行重合闸为例,展开分析。图7表示重合闸状态网络图,根据叠加原理,重合闸状态(图7)可由故障后状态(图4)和重合闸附加状态(图6)叠加。在单相重合闸对侧先闭合情况下,重合闸状态中,对侧故障相已重合,本侧故障相断路器断开。此时,根据叠加原理,在重合闸附加状态网络中,可以用电流源来表示断开的故障相,其值为故障状态下相电流的负值:
Figure BDA0001877696670000081
式中,
Figure BDA0001877696670000082
为故障状态下相电流,
Figure BDA0001877696670000083
为重合闸附加状态网络中的相电流。
在故障附加状态网络中,故障点对地支路有:
Figure BDA0001877696670000084
Figure BDA0001877696670000085
Figure BDA0001877696670000086
分别为重合闸故障点电压和电流,与状态B的式(1-4)相同,状态D也有:
Figure BDA0001877696670000087
Figure BDA0001877696670000088
为重合闸故障点右侧电流,由式(2-2)和式(2-3)可得:
Figure BDA0001877696670000089
Figure BDA00018776966700000810
为重合闸故障点左侧电流,由附加状态网络图(即图6),可得到:
Figure BDA00018776966700000811
Figure BDA00018776966700000812
Figure BDA00018776966700000813
重合闸附加状态网络中本侧电压和电流,将式(2-4)代入式(2-5)得:
Figure BDA00018776966700000814
Figure BDA00018776966700000815
其中,
Figure BDA00018776966700000816
Figure BDA00018776966700000817
在图6中,本侧测量电压也可由电源侧表示,如式(2-10)所示:
Figure BDA00018776966700000818
将式(2-6)和式(2-7)代入式(2-10),可得:
Figure BDA0001877696670000091
其中,
Figure BDA0001877696670000092
分别为A,B,C三相的重合闸故障点左侧电流,Zmm为本侧系统互阻抗,E1,E2和E3为:
Figure BDA0001877696670000093
Figure BDA0001877696670000094
Figure BDA0001877696670000095
式中,Xij,Yij(i,j=1,2,3)分别为X,Y矩阵中的元素,Zms和Zmm分别为本侧系统的自阻抗和互阻抗,
式(2-4)可化成如下形式:
Figure BDA0001877696670000096
与状态B类似,重合闸故障点的正序电压
Figure BDA0001877696670000097
为:
Figure BDA0001877696670000098
式中,下标(1)表示正序分量。
为减少计算量,本发明提出了与故障状态类似的正序附加网络,如图8所示,综上所述,最终确定的对于状态D的关系方程组为:
Figure BDA0001877696670000099
Figure BDA0001877696670000101
Figure BDA0001877696670000102
Figure BDA0001877696670000103
Figure BDA0001877696670000104
式中,Xij,Yij(i,j=1,2,3)分别为X,Y参数矩阵中的元素,矩阵
Figure BDA0001877696670000105
矩阵
Figure BDA0001877696670000106
Zms和Zmm分别为本侧系统的自阻抗和互阻抗,
Figure BDA0001877696670000107
Figure BDA0001877696670000108
分别为正序重合闸附加状态网络中本侧电压和电流,
Figure BDA0001877696670000109
为重合闸附加状态网络中的相电流。
针对重合闸另外一种情况,即本侧先重合闸时,其状态B方程与对侧先重合闸相同。状态D方程与对侧先重合闸基本相似,在重合闸附加状态网络中统一以故障状态下的电流负值代表电流源,进行分析,因此本发明不再累述。
(4)通过故障状态和重合闸状态的分析,可得到4个复数方程式,即状态B和D的四个方程组,将其实部和虚部分离后,便转化为8个实数方程。其中,未知数为故障距离,过渡电阻,对侧系统正序阻抗和负序阻抗的实部和虚部,共6个未知数。使用信赖域反射算法(trust-region reflective)求解方程组,便能得高精度的故障距离值,同时能计算出过渡电阻和对侧系统阻抗。
实施例详细说明如下:
应用本发明的一种输电系统模型如图9所示。线路长100km,电压等级为500kV;线路正序参数:Zl1=0.02+j0.25Ω,Yl1=j4.7×10-6S。零序参数:Zl0=0.1+j0.76Ω,Yl0=j2.5×10-6S。系统的本侧阻抗为:Zm1=9.2+j52.2Ω,Zm0=6.7+j37.9Ω;系统的对侧阻抗设置为:Zn1=3.132+j39.189Ω,Zn0=41.271+j121.315Ω。
应用本发明方法的实施例具体步骤如下:
(1)故障状态系统建模。通过分析故障附加状态网络(图3)和正序故障附加状态网络(图5),建立本侧电压电流关于故障距离、过渡电阻及对侧系统阻抗的方程,即状态B的关系方程组:
Figure BDA0001877696670000111
Figure BDA0001877696670000112
式中,
Figure BDA0001877696670000113
Figure BDA0001877696670000114
分别为正序故障附加状态网络中本侧电压和电流,Vm和Im分别为本侧电压和电流,AL(1),BL(1),Dm(1),Bm(1),CL(1),DL(1)为本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL的正序分量,S=1/3[1 α α2]且α=ej120°,本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL和对侧状态网络模型参数矩阵AR,BR,Dn,Bn,CR,DR中:Bm=ALZm+BL,Dm=CLZm+DL,Bn=ARZn+BR,Dn=CRZn+DR,AL=DL=coshγx,BL=Zc·sinhγx,
Figure BDA0001877696670000115
AR=DR=cosh(γ(d-x)),BR=Zc·sinh(γ(d-x)),
Figure BDA0001877696670000116
Zm和Zn分别为本侧系统阻抗和对侧系统阻抗,Zc和γ分别为波阻抗和传播系数,d为输电线路总长度,x表示故障点到本侧测量点的距离,
Figure BDA0001877696670000117
Zf为过渡阻抗。
(2)重合闸状态系统建模。通过分析重合闸时附加状态网络(图6)和正序重合闸时附加状态网络(图8),同样建立本侧电压电流关于故障距离、过渡电阻及对侧系统阻抗的方程,即状态D的关系方程组:
Figure BDA0001877696670000118
Figure BDA0001877696670000119
Figure BDA0001877696670000121
Figure BDA0001877696670000122
Figure BDA0001877696670000123
式中,Xij,Yij(i,j=1,2,3)分别为X,Y矩阵中的元素,矩阵
Figure BDA0001877696670000124
矩阵
Figure BDA0001877696670000125
Zms和Zmm分别为本侧系统的自阻抗和互阻抗,
Figure BDA0001877696670000126
Figure BDA0001877696670000127
分别为正序重合闸附加状态网络中本侧电压和电流,
Figure BDA0001877696670000128
为重合闸附加状态网络中的相电流。
(3)处理并简化所得方程。根据故障状态和重合闸状态分析,得到了4个复数方程,将其实部虚部分开,得到8个实数方程组,其中未知数为故障距离,过渡电阻,系统对侧正序阻抗和负序阻抗的实部和虚部,共6个未知数,小于方程数,理论上可由迭代法求得:
f(x)=0
x=[x1,x2,···,x6]T
f(x)=[f1(x),f2(x),···,f8(x)]T
上式用2范数表示,可得:
Figure BDA0001877696670000129
因此,对于方程组的求解,转换成对于
Figure BDA00018776966700001210
的求解。
(4)采用信赖域反射求解方法(trust-region reflective)求解方程组,该求解方法求解方程组的具体步骤如下:
步骤1.在当前迭代点Xk的邻域内,建立一个信赖域反射求解方法的子问题,如下式所示:
Figure BDA0001877696670000131
其中,s=x-xk是尝试迭代步,hk是目标函数g(x)在Xk处的梯度,Bk是g(x)的Hessian矩阵,Δk>0是信赖域半径,||s||是s的2范数。
步骤2.应用共轭梯度求解方法求解,获得尝试迭代步s。
步骤3.令Ak=g(xk)-g(xk+sk),其中,Ak是第K次迭代的实际下降量,sk为信赖域子问题的解;令Pk=qk(0)-g(sk),其中,qk(0)为二次模型函数,Pk是第K次迭代的是预测下降量;定义
Figure BDA0001877696670000132
若Rk满足条件(Rk接近于1,例如Rk>0.75),则Xk+1=Xk+sk,否则减小信赖域半径Δk,并且跳转到步骤1。
步骤4.检测是否收敛,若收敛,则结束迭代并输出结果。
(5)本发明的具体流程如图10所示,为保证利用信赖域反射算法得到的迭代结果的准确性,本发明将第一次迭代得到的结果代入第二次迭代中,并比较两次迭代的结果,若两次迭代的结果相差0.1%以上,则用第二次迭代结果代入第三次迭代中,如此循环下去,直到两次迭代的结果的误差在0.1%以内,得到最终测距算法结果。
应用EMTP-ATP软件的仿真结果如表1所示:
表1本发明故障测距结果(过渡电阻为100欧姆)
实际故障距离(km) 计算故障距离(km) 误差(%)
10 9.97 0.03
20 19.99 0.01
30 29.75 0.25
40 39.63 0.37
50 49.58 0.42
60 59.67 0.33
70 69.53 0.47
80 79.66 0.34
90 89.53 0.47
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (4)

1.一种基于重合闸动作分析的单端输电线路故障测距方法,其特征在于,包括以下步骤:
步骤1:将线路发生永久性故障前后整个过程分为5个时间面状态;
步骤2:针对5个时间面状态中的两个时间面状态分别建模并分别建立包含未知故障距离值和其他相关未知数的关系方程组;
步骤3:针对步骤2中的两个时间面状态,即故障状态和重合闸状态各自所建立的包含未知故障距离值和其他相关未知数的关系方程组进行求解,得出故障距离值;
步骤4:将故障距离值进行循环迭代,以满足误差的设定值以内并得出最终故障距离值;
所述步骤2中的两个时间面状态为状态B:故障状态,线路发生故障,但线路两侧断路器仍然闭合和状态D:重合闸状态,线路一侧断路器先行闭合,另外一侧断路器保持断开;
所述的状态B的关系方程组为:
Figure FDA0002619563100000011
Figure FDA0002619563100000012
式中,
Figure FDA0002619563100000013
Figure FDA0002619563100000014
分别为正序故障附加状态网络中本侧电压和电流,Vm和Im分别为本侧电压和电流,AL(1),BL(1),Dm(1),Bm(1),CL(1),DL(1)为本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL的正序分量,S=1/3[1 α α2]且α=ej120°,本侧状态网络模型参数矩阵AL,BL,Dm,Bm,CL,DL和对侧状态网络模型参数矩阵AR,BR,Dn,Bn,CR,DR中:Bm=ALZm+BL,Dm=CLZm+DL,Bn=ARZn+BR,Dn=CRZn+DR,AL=DL=coshγx,BL=Zc·sinhγx,
Figure FDA0002619563100000017
AR=DR=cosh(γ(d-x)),BR=Zc·sinh(γ(d-x)),
Figure FDA0002619563100000015
Zm和Zn分别为本侧系统阻抗和对侧系统阻抗,Zc和γ分别为波阻抗和传播系数,d为输电线路总长度,x表示故障点到本侧测量点的距离,
Figure FDA0002619563100000016
Zf为过渡阻抗;
所述的状态D的关系方程组为:
Figure FDA0002619563100000021
Figure FDA0002619563100000022
Figure FDA0002619563100000023
Figure FDA0002619563100000024
Figure FDA0002619563100000025
式中,Xij,Yij(i,j=1,2,3)分别为X,Y参数矩阵中的元素,矩阵
Figure FDA0002619563100000026
矩阵
Figure FDA0002619563100000027
Zms和Zmm分别为本侧系统的自阻抗和互阻抗,
Figure FDA0002619563100000028
Figure FDA0002619563100000029
分别为正序重合闸附加状态网络中本侧电压和电流,
Figure FDA00026195631000000210
为重合闸附加状态网络中的相电流。
2.根据权利要求1所述的一种基于重合闸动作分析的单端输电线路故障测距方法,其特征在于,所述步骤1中的5个时间面状态包括:状态A:故障前状态,线路两侧断路器闭合、状态B:故障状态,线路发生故障,但线路两侧断路器仍然闭合、状态C:故障清除状态,保护动作后,线路两侧断路器断开、状态D:重合闸状态,线路一侧断路器先行闭合,另外一侧断路器保持断开和状态E:最终隔离状态,故障未能清除,保护动作,再次断开合闸侧的断路器。
3.根据权利要求1所述的一种基于重合闸动作分析的单端输电线路故障测距方法,其特征在于,所述步骤3的包含未知故障距离值和其他相关未知数的关系方程组的求解方法采用信赖域反射求解方法。
4.根据权利要求1所述的一种基于重合闸动作分析的单端输电线路故障测距方法,其特征在于,所述步骤4中的设定值为0.1%。
CN201811407428.2A 2018-11-23 2018-11-23 一种基于重合闸动作分析的单端输电线路故障测距方法 Active CN109342888B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811407428.2A CN109342888B (zh) 2018-11-23 2018-11-23 一种基于重合闸动作分析的单端输电线路故障测距方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811407428.2A CN109342888B (zh) 2018-11-23 2018-11-23 一种基于重合闸动作分析的单端输电线路故障测距方法

Publications (2)

Publication Number Publication Date
CN109342888A CN109342888A (zh) 2019-02-15
CN109342888B true CN109342888B (zh) 2020-10-23

Family

ID=65317416

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811407428.2A Active CN109342888B (zh) 2018-11-23 2018-11-23 一种基于重合闸动作分析的单端输电线路故障测距方法

Country Status (1)

Country Link
CN (1) CN109342888B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110082634B (zh) * 2019-04-30 2020-06-09 中国矿业大学 广域电流时间序列的配电网单相接地故障定位方法
CN112485601B (zh) * 2020-12-11 2023-08-25 国网四川省电力公司电力科学研究院 基于双端线路电气量信息的故障分析方法和系统
CN113866561B (zh) * 2021-09-29 2023-10-27 天津大学 适用于输电线路经过渡电阻单相接地故障的单端量测距方法
CN115932477B (zh) * 2022-12-28 2024-01-23 国网湖北省电力有限公司电力科学研究院 一种融合多源信息的架空输电线路故障原因诊断方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101877479A (zh) * 2009-12-10 2010-11-03 贵州电力试验研究院 输电线路单相自适应重合闸相位比较判别方法
CN103207351A (zh) * 2013-03-12 2013-07-17 西安工程大学 一种基于重合闸的输电线路故障定位方法
CN103630798A (zh) * 2013-09-12 2014-03-12 清华大学 输电线路单相接地故障单端测距方法
CN107589344A (zh) * 2017-09-15 2018-01-16 广西电网有限责任公司电力科学研究院 一种基于探测电流频率分析的低压电力线路故障定位方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8289665B2 (en) * 2009-07-17 2012-10-16 The Invention Science Fund I Llc Systems and methods for grounding power line sections to clear faults

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101877479A (zh) * 2009-12-10 2010-11-03 贵州电力试验研究院 输电线路单相自适应重合闸相位比较判别方法
CN103207351A (zh) * 2013-03-12 2013-07-17 西安工程大学 一种基于重合闸的输电线路故障定位方法
CN103630798A (zh) * 2013-09-12 2014-03-12 清华大学 输电线路单相接地故障单端测距方法
CN107589344A (zh) * 2017-09-15 2018-01-16 广西电网有限责任公司电力科学研究院 一种基于探测电流频率分析的低压电力线路故障定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
单相自适应重合闸新原理的研究;李博通;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20090415;C042-297 *

Also Published As

Publication number Publication date
CN109342888A (zh) 2019-02-15

Similar Documents

Publication Publication Date Title
CN109342888B (zh) 一种基于重合闸动作分析的单端输电线路故障测距方法
Das et al. Transmission line fault detection and location using wide area measurements
Suonan et al. A novel fault-location method for HVDC transmission lines
Dobakhshari et al. A novel method for fault location of transmission lines by wide-area voltage measurements considering measurement errors
Liao et al. Online optimal transmission line parameter estimation for relaying applications
Azizi et al. Locating faults on untransposed, meshed transmission networks using a limited number of synchrophasor measurements
Gao et al. Design and evaluation of a directional algorithm for transmission-line protection based on positive-sequence fault components
CN104898021B (zh) 一种基于k‑means聚类分析的配电网故障选线方法
CN104391229A (zh) 一种基于s变换的输电线路故障快速选相方法
CN102129012B (zh) 一种利用形态峰谷检测的配电网故障选线方法
CN103543376A (zh) 用于小电流接地系统故障选线的径向基神经网络方法
CN103018632A (zh) 基于费歇信息的小电流接地系统单相接地故障选线方法
Cao et al. Asynchronous fault location scheme based on voltage distribution for three-terminal transmission lines
CN107817420A (zh) 非全程同杆并架双回输电线路非同步数据故障测距方法
CN104237677A (zh) 220kV及以上电网单回输电线路零序参数在线辨识方法
Taheri et al. A fault‐location algorithm for parallel line based on the long short‐term memory model using the distributed parameter line model
Lin et al. A signal-superimposed technique for fault location in transmission lines through IED measurements considering communication service failure
Han et al. Fault location on a mixed overhead and underground transmission feeder using a multiple-zone quadrilateral impedance relay and a double-ended travelling wave fault locator
Ghaedi et al. Modified WLS three-phase state estimation formulation for fault analysis considering measurement and parameter errors
Galvez et al. Fault location in meshed and active power distribution networks
CN112782532A (zh) 基于断路器合闸产生行波信号的配电网故障测距方法
CN109001589B (zh) 基于非全相运行故障录波数据的线路参数时域计算方法
CN110879332A (zh) 一种适用于小电流接地系统的单相接地故障选相方法
CN112363009B (zh) 一种用于同塔线路接地故障的单端故障测距方法及系统
George et al. Traveling wave based autoreclosure scheme for multi-terminal lines

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant