CN103913676B - 基于可变行波辨识时窗的输电线路单端故障定位方法 - Google Patents
基于可变行波辨识时窗的输电线路单端故障定位方法 Download PDFInfo
- Publication number
- CN103913676B CN103913676B CN201410142294.1A CN201410142294A CN103913676B CN 103913676 B CN103913676 B CN 103913676B CN 201410142294 A CN201410142294 A CN 201410142294A CN 103913676 B CN103913676 B CN 103913676B
- Authority
- CN
- China
- Prior art keywords
- fault
- reflection wave
- window
- wave
- ripple
- 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.)
- Expired - Fee Related
Links
Landscapes
- Locating Faults (AREA)
Abstract
本发明公开了一种基于可变行波辨识时窗的输电线路单端故障定位方法,该方法通过对端***阻抗的在线估计,构建改进的分布参数法,有效改善常规方法误差分布的不规律性,提升依据工频量定位误差构建组合定位算法的可行性;根据不同故障条件和误差规律计算可变的故障距离窗口,锁定故障行波反射波的对应时间区域,由此划定行波反射波辨识的可变时窗;基于该可变时窗,综合考虑测量本端和对端的母线接线方式、相邻母线反射波的影响等因素,利用行波法构建可变时窗内的反射波辨识方法及定位方法。仿真验证结果表明,该方法显著提高了组合定位方法的可靠水平,并具有较高的故障定位精度,具有一定工程应用价值。
Description
技术领域
本发明涉及故障定位领域,尤其涉及一种基于可变行波辨识时窗的输电线路单端故障定位方法。
背景技术
单端测距与双端测距相比,不受通讯条件、经济条件、管理权限等因素制约,其固有优点使单端测距拥有广阔的发展前景。但目前单端测距仍存在一些问题:工频量法由于理论限制,无法根本上消除对端助增电流的影响,实际应用误差相对较大;行波法理论上不受***运行方式、故障过渡电阻影响,测距精度高于工频量法,但行波法的关键在于行波折反射波的准确辨识,而工程现场行波波形复杂,波头辨识困难,无效的波头辨识将导致错误的定位结果。
目前从行波波头辨识的角度探讨测距精度的提升方法主要有以下几种:
1.利用极性关系识别反射波的性质,但该极性关系受限于现场母线类型。
2.针对“三一类”母线结构利用初始反极性行波的单端行波测距方法,该方法的应用范围具有限制性。
3.考虑到母线结构的综合电压行波极性和电流行波极性识别反射波的方法,但是该方法的前提条件是精确测量现场的电压行波。
4.利用单端阻抗法对行波法的多个结果进行筛选,基于阻抗法误差规律修正线路中点附近故障测距结果,在现场故障情况复杂的情况下,阻抗法测距精度难以遵循一定规律。
5.利用阻抗法测距结果及10%的阻抗法误差限定故障点反射波的识别范围,但测量端有两条出线时故障点反射波难以检测,且10%的误差范围不利于误差较小的故障情况有效筛除干扰波。
由此可见,组合故障定位方法具有较强理论性和工程应用可行性,但同时受限于母线接线方式和行波反射波限定范围的有效性。因此,如何结合实际误差情况灵活确定限定反射波的区域,并在接线方式多样化的情况下准确标识并选择限定范围内的定位波头,对于进一步提升组合定位方法的可行性具有重要意义。
发明内容
本发明的目的就是为了解决上述问题,提出了一种基于可变行波辨识时窗的输电线路单端故障定位方法,该方法可同时提高测距结果的精度和可靠性,以及对不同线路结构的通用性。
为了实现上述目的,本发明采用如下技术方案:
一种基于可变行波辨识时窗的输电线路单端故障定位方法,包括以下步骤:
(1)假设对端电源***:M端母线总出线数为m,N端母线总出线数为n,带变压器出线不计入母线总出线数;在线路MN之间一点F发生故障,提取故障前若干组工频电气量,利用最小二乘原理在线估计对端电源***的阻抗。
(2)根据在线估计的对端电源***阻抗,结合分布参数法,得到故障定位结果D,根据D与线路全长L的关系,得到窗口长度随故障位置变化的故障距离窗口。
(3)检测故障初始行波到达测量端的时刻对应采样点t0,将可变的故障距离窗口转换为时域上两个变长时窗:第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2。
(4)假设M端、N端的相邻母线分别为P、Q,相邻母线总出线数分别为p、q,相应线路长度分别为|MP|、|NQ|;判断是否存在相邻线路的干扰波,若存在,排除干扰波影响。
(5)根据线路结构和故障情况,采用不同的波头辨识策略,确定定位波头的类型是故障点反射波还是对端母线反射波;并确定故障点反射波到达测量端时刻对应的采样点t1或者对端母线反射波到达测量端时刻对应的采样点t2。
(6)根据定位波头的类型,构建故障定位方程,确定故障定位结果。
所述步骤(1)中在线估计对端电源***阻抗的方法为:
其中,为故障前N端第i组工频量,为故障前M端第i组工频量,为对端电源***阻抗的在线估计值;Zc为线路的波阻抗,γ为线路正序传播系数,L为线路全长,k为故障前提取的工频量组数。
所述步骤(2)中的故障距离窗口为:
X=[D-εL,D+εL]
其中,ε1<ε2
其中,L为故障线路全长,D为改进的分布参数法定位结果,ε为定位误差;μ为线路长度百分比,ε1、ε2分别为不同故障距离时的定位误差。
所述步骤(3)中的第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2分别为:
其中,t0为初始故障行波到达测量端的时刻对应的采样点,f为行波采样频率,v为行波波速,L为故障线路全长,D为改进的分布参数法定位结果,ε为定位误差。
所述步骤(4)中判断是否存在相邻线路的干扰波的方法为(下文所提同极性、反极性均是相对于初始故障行波的极性):
若相邻母线反射波出现在故障线路窗口T1或T2内,且该相邻母线反射波的极性与窗口内待辨识波头极性相同,则表明存在相邻线路的干扰波,将此相邻母线反射波定义为相邻线路干扰波头。
所述步骤(4)中排除干扰波影响的方法为:
确定受影响的故障线路窗口T1或T2,根据干扰线路长度,确定相邻线路反射波到达测量端母线的时刻对应的采样点,锁定与该采样点最近且极性与窗口内待辨识波头的极性相同的波头,将该波头对应的小波模极大值置零。相邻线路干扰波对故障线路对应窗口T1或T2的影响如下表所示:
所述步骤(5)中波头辨识策略为:
故障线路窗口T1和T2中仅有一个窗口存在明显反射波的接线方式,直接利用该窗口内对应极性的奇异性最强的波头作为定位波头:
对于m=2、n=1类型的线路,在T2内选择与初始故障行波极性相同的小波模极大值的绝对值最大的波头作为对端母线反射波。
对于m=2、n>2类型的线路,在T2内选择与初始故障行波反极性的小波模极大值的绝对值最大的波头作为对端母线反射波。
对于m>2、n=2类型的线路,在T1内选择与初始故障行波同极性的小波模极大值的绝对值最大的波头作为故障点反射波。
所述步骤(5)中波头辨识策略为:
故障线路窗口T1和T2中均存在明显反射波的接线方式,即m>2、n>2和m>2、n=1类型的线路,
若故障线路窗口T1和T2无交集,利用各自窗口内对应极性的奇异性最强的波头作为定位波头:T1内选择与初始故障行波同极性的小波模极大值的绝对值最大的波头作为故障点反射波,n>2结构在T2内选择与初始故障行波反极性的小波模极大值的绝对值最大的波头作为对端母线反射波,n=1结构在T2内选择与初始故障行波同极性的小波模极大值的绝对值最大的波头作为对端母线反射波;在T1、T2内所选的反射波头中取幅值最大者作为最终的定位波头。
若故障线路窗口T1和T2相交,即时,有以下分析:
对于m>2、n>2类型的接线方式,两端母线反射系数均为正值,故T1、T2内故障点反射波和对端母线反射波极性相反,辨识方法与两窗口无交集时相同。
对于m>2、n=1类型的接线方式,测量端母线反射系数为正值,对端母线反射系数为负值,则T1、T2内故障点反射波和对端母线反射波极性相同,判定故障线路窗口T1和T2合集范围内初始故障行波同极性的小波模极大值的绝对值最大的波头为对端母线反射波。
所述步骤(6)的具体方法为:
若定位波头为故障点反射波,测距结果为:
若定位波头为对端母线反射波,测距结果为:
若波头检测失败,测距结果为改进分布参数法定位结果x=D。
其中,x为本发明方法最终确定的故障点距测量端母线的距离,t0为故障初始行波到达测量端的时刻对应采样点,t1为故障点反射波到达测量端时刻对应的采样点,t2为对端母线反射波到达测量端时刻对应的采样点,f为行波采样频率,v为行波波速。
本发明的有益效果是:
本发明能够提高工频定位结果的准确性,并有效改善误差的稳定性,通过提取其误差分布规律,得到长度可变的行波辨识时窗,能够提升行波辨识的抗干扰能力。初始故障行波与后续反射波的相对极性取决于故障线路的母线接线方式,通过分析其特征规律,可实现不同接线方式下特定时域窗口内行波波头的准确辨识,以及定位波头的筛选。
仿真验证结果表明,本发明的方法不受母线接线结构、故障位置、过渡电阻、***运行方式等因素影响,精度、可靠性高,具有一定工程应用价值。
附图说明
图1为对端电源***故障示意图;
图2为工频量法确定的行波反射波的划定时窗示意图;
图3为利用戴维南等值电路表达对端电源***示意图;
图4为不同故障位置的改进分布参数法测距误差示意图;
图5为故障行波的折反射示意图;
图6为母线接线结构图;
图7(a)为m=2、n=1线路结构的定位波头图;
图7(b)为m=2、n>2线路结构的定位波头图;
图7(c)为m>2、n=2线路结构的定位波头图;
图7(d)为m>2、n>2线路结构的定位波头图;
图7(e)为两窗口不相交情况下m>2、n=1线路结构的定位波头图;
图7(f)为两窗口相交情况下m>2、n=1线路结构的定位波头图;
图8为本发明基于可变行波辨识时窗的输电线路单端故障定位方法流程图;
图9为本发明算例二158km处故障时模型三数据处理结果;
图10(a)为本发明算例四中不做数据预处理的结果;
图10(b)为本发明算例四中做数据预处理的结果;
具体实施方式:
下面结合附图与实施例对本发明做进一步说明:
行波法受限于现场强噪声、复杂操作干扰等因素,反射波辨识较困难。
图1所示为对端电源***图,线路MN于F点处发生故障,其故障行波如图2所示,辨识过程中发现故障行波中存在与真实反射波同极性的干扰波头1、2、3,分别对应28.7km、47.3km、118.3km,与真实故障距离99.6km相比误差很大,给真实反射波的有效辨识造成困难。考虑到工频量法测距依靠的是物理量之间的电气关系,具有误差相对较大但可靠性较高的特点,因此,利用工频量法定位结果限定故障距离范围,则利于对故障行波限定反射波的辨识窗口,可有效提升行波定位的抗干扰能力。
图2所示为工频量法确定的行波反射波的划定时窗,若工频量测距误差偏大,如窗口W3,则无法包含真实波头,将可能导致完全错误的测距结果;若窗口包含真实反射波但窗口长度偏大,如窗口W2,则无法有效消除真实反射波辨识过程中的干扰,并可能引入奇异性更强的干扰波头3;可见,包含真实反射波头且时窗长度合理的W1更利于消除干扰影响,给出高精度定位结果。
上述分析可得,反射波时窗划定需满足准确性和稳定性的要求,为此需提升工频量定位方法误差分布的稳定性,并准确提取其误差分布规律。
改进工频量分布参数法
分布参数法的主要缺点在于对端***阻抗不可知,随着对端***运行方式的改变而变化,常规方法是假定一个固定值作为对端***阻抗值用于计算。本发明通过对端***阻抗在线估计,对分布参数法进行改进,故称为“改进的分布参数法”。具体过程与常规方法的区别在于改进法利用的对端***阻抗值为在线估计值,以代替常规方法的固定值。
相比阻抗法而言,常规工频量分布参数法采用精确线路模型,定位精度明显高于集中参数模型测距算法,更符合时窗对工频量定位方法准确性的要求。但目前分布参数法均基于对端***阻抗为定值的假设,该假设在对端***阻抗角变化时,即使对于同一故障情况,故障定位精度将会大幅改变,因此,其误差分布具有明显不稳定性。
同时,对于单端测距方式而言,对端母线运行方式的改变具有未知性,无法实时更改计算参数。由此,若能实现对端***阻抗参数的在线辨识,真实反映对端***运行方式的变化,则可实时修正常规分布参数法中的对端***阻抗参数,提升定位误差分布稳定性。可见,改进分布参数法的关键问题在于对端***阻抗的在线估算方法。
利用戴维南等值电路表达对端***如图3所示,其中,ZN分别为对端等值电源和对端***阻抗。单端测距中对端***的运行方式未知,但其在线路故障时刻前后几个周波内发生变化的可能性几乎为零,可认为该时间段内ZN均为定值。而分布参数法测距利用的正是该时间段内故障期间的电气量,因此可由节点N处故障前几个周波内的工频电气量确定分布参数测距算法内的ZN值,实现对端***阻抗参数的实时修正。
电力***的正常运行过程始终伴随负荷波动、电压波动等小扰动的存在,因此线路正常运行时测量端的电压电流相量是实时变化的。故障录波装置依据行业标准,需记录***扰动开始前至少两个周波的实时模拟量数据,这为故障前多组工频稳态数据的获取提供了实际应用基础。提取故障前多组工频电气量,利用最小二乘原理估算本次故障期间对端***阻抗值:
其中,为故障前N端第i组工频量,为故障前M端第i组工频量,为对端电源***阻抗的在线估计值。
反射波变长时窗划定
长度合理的窗口需在满足准确性的前提下,尽可能缩小反射波辨识范围以排除干扰。若能跟踪不同故障情形下的误差规律,据此动态确定窗口长度,则可有效提高行波辨识的抗干扰能力。定位误差的规律性在改进分布参数法中得到了有效提升,在此基础上,需探索总结改进分布参数法定位精度的规律,同时需考虑不同故障条件对有效窗口长度的影响,由此划定反射波变长时窗范围。
在实时估算对端***阻抗后,利用改进分布参数法定位故障点的过程实质上是求解关于故障距离的一元方程,算法本身基本不受故障点过渡电阻值和故障类型的影响,故只需考虑故障点位置对定位误差的作用。针对一典型500kV对端电源***设计参数模型,线路全长范围内改进分布参数法的误差如图4所示。可见,改进分布参数法对于线路近测量端的大部分长度范围内的故障,定位误差小且范围集中,线路末端长度范围内的故障,定位误差随故障距离增加而加速变大。
改进分布参数法具有明显的故障位置分布规律,据此可得到窗口长度随故障位置变化的故障距离窗口:
X=[D-εL,D+εL](3)
其中,ε1<ε2(4)
式(3)中,L为故障线路全长,D为改进的分布参数法定位结果,ε为定位误差。式(4)中,μ为线路长度百分比,ε1、ε2为不同故障距离时的定位误差,μ、ε1、ε2三者取值需依据实际***情况与故障情况确定。
变长故障距离窗可转换为时域上两个变长时窗,分别包含第一个故障点反射波和第一个对端母线反射波到达测量端母线的时刻。转换如下:
式(5)中,t0为初始故障行波到达测量端的时刻对应的采样点,T1为第一个故障点反射波对应的采样点窗口,T2为第一个对端母线反射波对应的采样点窗口,f为行波采样频率,v为行波波速。
单端行波法关键在于对反射波的准确辨识。在故障窗口给定的条件下,反射波的辨识被限定在特定的较小的区间内,此时在给定窗口内的行波波头辨识成为单端测距的重要步骤。
如图5所示,输电线路M的F点于t0时刻发生故障。M端检测到的初始电流行波i0(t)、第一个故障点反射波i1(t)、第一个对端母线反射波i2(t)可表达为:
式(6)中,ρM、ρN、ρF分别为电流行波在母线M、N处和故障点F处的反射系数,γF为电流行波在F处的折射系数,τMF、τNF分别为电流行波在MF和NF线路的传播时间,ib、if分别为故障时刻故障点产生的沿线路和方向传播的故障电流行波。其中, Z1、Z2分别是不连续点前后线路波阻抗;显然(1+ρM)>0。
设故障点F存在过渡电阻值为RF,则阻抗不连续点后的等效特征阻抗为:
可见,Z2<Z1,故ρF>0,γF>0,那么M端检测到的电流行波的极性与母线MN处的反射系数ρM、ρN有关,因此,i1(t)、i2(t)分别与i0(t)的相对极性取决于故障线路两端母线的接线方式。
基于反射波极性与母线结构的分析,对于不同接线结构,有以下行波波头辨识方法(下文提到的同极性、反极性均是相对初始故障行波的极性而言),其中母线接线方式如图6所示,M端母线总出线数为m,N端母线总出线数为n,带变压器出线不计入母线总出线数。
1)对于时窗T1、T2仅其中一个窗口存在明显反射波的接线方式,直接利用该窗口内对应极性的奇异性最强的波头作为定位波头。
如图7(a)~7(c)所示,m=2、n=1,m=2、n>2类型线路在T2内分别选择同极性和反极性的绝对最大模极大值的波头作为对端母线反射波,m>2、n=2类型线路在T1内选择同极性的绝对模极大值最大的波头作为故障点反射波。
2)对于时窗T1、T2二者均存在明显反射波的接线方式,即m>2、n>2与m>2、n=1类型线路,若两窗口无交集,利用各自窗口内对应极性的奇异性最强的波头作为定位波头,如图7(d)-(e)所示。
若两窗口相交,即时,有以下分析:
m>2、n>2类型的接线方式,两端母线反射系数均为正值,故T1、T2内故障点反射波和对端母线反射波极性相反,可利用极性辨识波头。m>2,n>2类型的接线方式故障点反射波和对端母线反射波极性相反,所以两窗口相交不影响极性相反的两个波头的辨识。辨识方法与两窗口不相交时相同。
m>2、n=1类型的接线方式,测量端母线反射系数为正值,对端母线反射系数为负值,则T1、T2内故障点反射波和对端母线反射波极性相同,极性判别方法失效,需利用行波的奇异性实现波头辨识。行波传输过程中奇异性的衰减主要受传输长度和故障点过渡电阻影响,由于对端n=1的结构导致N处电流行波发生全反射,对端母线反射波的奇异性较强,然而近M端故障时,故障点反射波传输距离短于对端母线反射波,对端母线反射波的传播衰减更大,故不能直接比较二者的奇异性强弱,需进一步考虑过渡电阻值对行波衰减的影响。
现针对该线路结构利用ATP构建500kV输电***仿真模型,故障线路MN长300km,设置在T1、T2存在交集条件下,近M端在不同过渡电阻时故障,其中,波头1、2、3分别为初始行波、第一个故障点反射波和第一个对端母线反射波。
对端母线反射波的奇异性始终远强于故障点反射波,且二者差异随故障电阻值的增加而扩大。据此,判定T1、T2两窗口合集范围内与初始行波同极性的奇异性最强的波头为对端母线反射波,如图7(f)所示。
相邻母线反射波干扰时的数据预处理
如图6所示,M端、N端相邻母线分别为P、Q,母线总出线数分别为p、q,相应线路长度分别为|MP|、|NQ|。若相邻母线反射波出现在故障线路窗口T1、T2内,且与窗口内待辨识波头极性相同,则将此相邻母线反射波定义为相邻线路干扰波头。
相邻线路干扰波头的极性与相邻母线的接线方式有关,表1给出了不同线路结构下相邻线路干扰波对故障线路对应窗口内波头辨识的影响。
表1相邻线路干扰波对故障线路对应窗口的影响
获得变长时窗后,测量端结合表1检查是否存在满足产生干扰波头条件的相邻线路。若存在,有以下处理对策:确定受影响的时间窗口,根据干扰线路长度,锁定与对应的采样点最近且极性匹配的波头,将其对应的小波模极大值置零。极端情况下,干扰波头可能与待辨识波头重合,置零操作会湮没定位波头,但改进分布参数法定位窗口仍可保证测距结果的可靠性与精度。
本发明算法流程图如图8所示:
1)在线估计对端***阻抗。
2)改进的分布参数法给出定位结果D,根据D与线路全长L的关系,得到可变的故障距离窗口。
3)对故障后电流线模进行小波变换及模极大值计算,检测故障初始行波到达测量端的时刻对应采样点数t0,计算变长时窗T1和T2。
4)检测存在相邻线路干扰波的条件是否成立,若成立,排除干扰波影响。
5)根据线路结构和故障情况,采用不同的波头辨识策略,得到故障点反射波到达测量端的时刻对应采样点t1或对端母线反射波到达测量端的时刻对应采样点t2,确定构建测距方程的定位波头。
6)若定位波头为故障点反射波,测距结果为:
若定位波头为对端母线反射波,测距结果为:
若波头检测失败,测距结果为改进分布参数法定位结果x=D。
仿真验证及分析
利用ATP构建500kV输电线路仿真模型,如图10所示。线路MN长度为300km,故障时刻为0.2s,工频采样频率为20kHz。改变对端负载S2模拟对端***运行方式变化,故障前几个周波内改变测量端M上级线路在母线处的负载值S1,模拟上级线路小扰动,获取MN线路测量端的多组数据,据此估计对端***阻抗。
表2给出了对端***运行方式改变时,基于对端***阻抗估计的分布参数法对于不同故障距离、过渡电阻情况下的定位误差。
表2改进分布参数法定位误差/%
表2显示,改进的分布参数法对于对端***运行方式变化不敏感,定位结果准确性、可靠性均满足要求。在此基础上,验证本文算法的有效性。
可变时窗内的行波辨识及定位
为了体现算法对各种母线接线方式的适用性,本文分别建立5种母线结构的500kV输电线路仿真模型(模型一至五的接线方式分别为:m=2、n=1,m=2、n>2,m>2、n=1,m>2、n=2,m>2、n>2),如图6所示,故障线路MN长度为300km,M端为测量端,五种模型仅母线结构存在差异,MP、NQ长度均为300km。故障时刻为0.01s,工频采样频率为20kHz,采样0.2s;行波采样频率为1MHz,采样0.02s,行波线模波速为0.294km/μs。对线路MN分别设置近端故障、中点附近故障、远端高阻故障,且每种故障类型均对五种模型进行仿真。
算例1:线路MN近端60km处发生故障,过渡电阻值20Ω,行波中包含5%噪声。近端故障时,改进分布参数法定位较精确,令其误差ε=0.03。五种线路结构的测距结果如表3所示。
表3近区低阻故障时测距结果
算例2:线路中点附近故障。线路MN距M端158km处发生故障,过渡电阻值30Ω,行波中包含7%噪声。非末端故障时,改进分布参数法定位较精确,令其误差ε=0.05,那么T1、T2两窗口存在交集,模型一至五的测距结果如表4所示。
表4线路中区故障测距结果
算例3:线路MN末端275km处发生高过渡电阻故障,过渡电阻值100Ω,行波中包含7%噪声。末端故障时,改进分布参数法定位窗口偏大,令其误差ε=0.05。五种线路结构模型的测距结果见表5。
表5末端高阻故障测距结果
算例4:调整模型五中相邻线路长度:MP、NQ长度均为100km。线路MN距M端205km处发生故障,过渡电阻值20Ω,行波中包含7%噪声。两时窗内故障行波小波变换模极大值见图10(a)和图10(b)。
由图10(a)和图10(b)可看出,在此接线类型和故障条件下,时窗T1、T2内分别出现来自Q端母线和P端母线的干扰波头。不考虑相邻线路影响时,测距结果为199.893km,误差5.107km,利用本文方法对相邻线路影响进行处理后,测距结果为204.891km,误差0.109km。
综上所述,算法对于各种线路结构具有通用性,且测距精度较高。模型一、二测量端母线有两条出线,故障点反射波微弱,算法依靠对端母线反射波区间T2测距,测距结果较精确。模型四由于没有明显的对端母线反射波,仅依靠窗口T1测距,在高阻故障时,故障点反射波较弱,算法误差较其余模型略大,在0.2%以内。相邻母线反射波造成干扰时,数据预处理的操作可消除干扰影响,明显提高测距精度。
针对现场波形复杂、接线方式多样导致的行波波头难以准确辨识的问题,本文提出一种基于可变行波辨识时窗的输电线路单端故障定位方法。基于对端***阻抗在线估计的改进分布参数法能够提高工频定位结果的准确性,并有效改善误差的稳定性,通过提取其误差分布规律,得到长度可变的行波辨识时窗,能够提升行波辨识的抗干扰能力。初始故障行波与后续反射波的相对极性取决于故障线路的母线接线方式,通过分析其特征规律,可实现不同接线方式下特定时域窗口内行波波头的准确辨识,以及定位波头的筛选。
仿真验证结果表明,本文所提方法不受母线接线结构、故障位置、过渡电阻、***运行方式等因素影响,精度、可靠性高,具有一定工程应用价值。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,包括以下步骤:
(1)假设对端电源***:M端母线总出线数为m,N端母线总出线数为n,带变压器出线不计入母线总出线数;在线路MN之间一点F发生故障,提取故障前若干组工频电气量,利用最小二乘原理在线估计对端电源***的阻抗;
(2)根据在线估计的对端电源***阻抗,结合分布参数法,得到故障定位结果D,根据D与线路全长L的关系,得到窗口长度随故障位置变化的故障距离窗口;
(3)检测故障初始行波到达测量端的时刻对应采样点t0,将可变的故障距离窗口转换为时域上两个变长时窗:第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2;
(4)假设M端、N端的相邻母线分别为P、Q,相邻母线总出线数分别为p、q,相应线路长度分别为|MP|、|NQ|;判断是否存在相邻线路的干扰波,若存在,排除干扰波影响;
(5)根据线路结构和故障情况,采用不同的波头辨识策略,确定定位波头的类型是故障点反射波还是对端母线反射波;并确定故障点反射波到达测量端时刻对应的采样点t1或者对端母线反射波到达测量端时刻对应的采样点t2;
(6)根据定位波头的类型,构建故障定位方程,确定故障定位结果。
2.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(1)中在线估计对端电源***阻抗的方法为:
其中,为故障前N端第i组工频量,为故障前M端第i组工频量,为对端电源***阻抗的在线估计值;Zc为线路的波阻抗,γ为线路正序传播系数,L为线路全长,k为故障前提取的工频量组数。
3.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(2)中的故障距离窗口为:
X=[D-εL,D+εL]
其中,ε1<ε2
其中,L为故障线路全长,D为改进的分布参数法定位结果,ε为定位误差;μ为线路长度百分比,ε1、ε2分别为不同故障距离时的定位误差。
4.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(3)中的第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2分别为:
其中,t0为初始故障行波到达测量端的时刻对应的采样点,f为行波采样频率,v为行波波速,L为故障线路全长,D为改进的分布参数法定位结果,ε为定位误差。
5.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(4)中判断是否存在相邻线路的干扰波的方法为:
若相邻母线反射波出现在第一个故障点反射波对应的采样点窗口T1或第一个对端母线反射波对应的采样点窗口T2内,且该相邻母线反射波的极性与窗口内待辨识波头极性相同,则表明存在相邻线路的干扰波,将此相邻母线反射波定义为相邻线路干扰波头。
6.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(4)中排除干扰波影响的方法为:
确定受影响的第一个故障点反射波对应的采样点窗口T1或第一个对端母线反射波对应的采样点窗口T2,根据干扰线路长度,确定相邻线路反射波到达测量端母线的时刻对应的采样点,锁定与该采样点最近且极性与窗口内待辨识波头的极性相同的波头,将该波头对应的小波模极大值置零。
7.如权利要求6所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,相邻线路干扰波对第一个故障点反射波对应的采样点窗口T1或第一个对端母线反射波对应的采样点窗口T2的影响如下表所示:
8.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(5)中波头辨识策略为:
第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2中仅有一个窗口存在明显反射波的接线方式,直接利用该窗口内对应极性的奇异性最强的波头作为定位波头:
对于m=2、n=1类型的线路,在T2内选择与初始故障行波极性相同的小波模极大值的绝对值最大的波头作为对端母线反射波;
对于m=2、n>2类型的线路,在T2内选择与初始故障行波反极性的小波模极大值的绝对值最大的波头作为对端母线反射波;
对于m>2、n=2类型的线路,在T1内选择与初始故障行波同极性的小波模极大值的绝对值最大的波头作为故障点反射波。
9.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(5)中波头辨识策略为:
第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2中均存在明显反射波的接线方式,即m>2、n>2和m>2、n=1类型的线路,
若第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2无交集,利用各自窗口内对应极性的奇异性最强的波头作为定位波头:T1内选择与初始故障行波同极性的小波模极大值的绝对值最大的波头作为故障点反射波,n>2结构在T2内选择与初始故障行波反极性的小波模极大值的绝对值最大的波头作为对端母线反射波,n=1结构在T2内选择与初始故障行波同极性的小波模极大值的绝对值最大的波头作为对端母线反射波;在T1、T2内所选的反射波头中取幅值最大者作为最终的定位波头;
若第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2相交,即时,有以下分析:
对于m>2、n>2类型的接线方式,两端母线反射系数均为正值,故T1、T2内故障点反射波和对端母线反射波极性相反,辨识方法与两窗口无交集时相同;
对于m>2、n=1类型的接线方式,测量端母线反射系数为正值,对端母线反射系数为负值,则T1、T2内故障点反射波和对端母线反射波极性相同,判定第一个故障点反射波对应的采样点窗口T1和第一个对端母线反射波对应的采样点窗口T2合集范围内初始故障行波同极性的小波模极大值的绝对值最大的波头为对端母线反射波。
10.如权利要求1所述的一种基于可变行波辨识时窗的输电线路单端故障定位方法,其特征是,所述步骤(6)的具体方法为:
若定位波头为故障点反射波,测距结果为:
若定位波头为对端母线反射波,测距结果为:
若波头检测失败,测距结果为改进分布参数法定位结果x=D;
其中,x为本发明方法最终确定的故障点距测量端母线的距离,t0为故障初始行波到达测量端的时刻对应采样点,t1为故障点反射波到达测量端时刻对应的采样点,t2为对端母线反射波到达测量端时刻对应的采样点,f为行波采样频率,v为行波波速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410142294.1A CN103913676B (zh) | 2014-04-10 | 2014-04-10 | 基于可变行波辨识时窗的输电线路单端故障定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410142294.1A CN103913676B (zh) | 2014-04-10 | 2014-04-10 | 基于可变行波辨识时窗的输电线路单端故障定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103913676A CN103913676A (zh) | 2014-07-09 |
CN103913676B true CN103913676B (zh) | 2016-06-01 |
Family
ID=51039509
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410142294.1A Expired - Fee Related CN103913676B (zh) | 2014-04-10 | 2014-04-10 | 基于可变行波辨识时窗的输电线路单端故障定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103913676B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569744B (zh) * | 2014-11-26 | 2018-03-20 | 国家电网公司 | 一种适用于配电网线路的综合单端故障定位方法 |
CN104459470B (zh) * | 2014-12-11 | 2018-04-06 | 东北大学 | 一种适用于多类型母线结构的行波故障定位方法 |
CN105548803B (zh) * | 2015-12-04 | 2019-09-10 | 国网山东省电力公司电力科学研究院 | 基于加权最小二乘算法的特高压输电线路故障定位方法 |
CN106054023B (zh) * | 2016-06-30 | 2019-03-19 | 国网江西省电力公司检修分公司 | 一种输电线路单端测距中估计两侧***阻抗的方法 |
CN108594068B (zh) * | 2018-04-04 | 2020-09-08 | 南京南瑞继保电气有限公司 | 一种行波测距方法 |
CN109870628B (zh) * | 2018-08-31 | 2020-12-04 | 国网江苏省电力有限公司苏州供电分公司 | 一种多端柔性直流输电***故障线路识别方法 |
CN110895299A (zh) * | 2018-09-12 | 2020-03-20 | 中国石油化工股份有限公司 | 一种电力电缆测距识别方法 |
CN109470987B (zh) * | 2018-09-30 | 2021-04-09 | 昆明理工大学 | 一种基于区段匹配算法t接输电线路单端行波测距方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4499417A (en) * | 1981-10-27 | 1985-02-12 | General Electric Company P.L.C. | Determining location of faults in power transmission lines |
US4797805A (en) * | 1985-12-20 | 1989-01-10 | Asea Aktiebolag | Fault location in a power supply network |
CN101923139A (zh) * | 2010-04-19 | 2010-12-22 | 昆明理工大学 | 一种输电线路单端行波故障测距智能方法 |
CN102426323A (zh) * | 2011-10-31 | 2012-04-25 | 广东电网公司广州供电局 | 行波法定位中计算入射波与反射波时间差的方法及装置 |
-
2014
- 2014-04-10 CN CN201410142294.1A patent/CN103913676B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4499417A (en) * | 1981-10-27 | 1985-02-12 | General Electric Company P.L.C. | Determining location of faults in power transmission lines |
US4797805A (en) * | 1985-12-20 | 1989-01-10 | Asea Aktiebolag | Fault location in a power supply network |
CN101923139A (zh) * | 2010-04-19 | 2010-12-22 | 昆明理工大学 | 一种输电线路单端行波故障测距智能方法 |
CN102426323A (zh) * | 2011-10-31 | 2012-04-25 | 广东电网公司广州供电局 | 行波法定位中计算入射波与反射波时间差的方法及装置 |
Non-Patent Citations (1)
Title |
---|
行波法与阻抗法结合的综合单端故障测距新方法;卢继平 等;《电力***自动化》;20071210;第31卷(第23期);第65-69页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103913676A (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103913676B (zh) | 基于可变行波辨识时窗的输电线路单端故障定位方法 | |
CN102288869B (zh) | 一种输电线路单端行波故障测距方法 | |
CN108508320B (zh) | 基于谐波能量和波形畸变特征的弧光接地故障辨识方法 | |
CN103323741B (zh) | 一种针对强故障的基于故障电压初始行波幅值比较的d型线缆混合线路故障区段判别的方法 | |
CN104090214B (zh) | 一种电缆故障检测及老化分析方法 | |
CN103293449B (zh) | 一种消除煤矿高压电网单端行波故障定位死区的方法 | |
CN110007141B (zh) | 一种基于电压电流谐波相似度的谐振点检测方法 | |
CN106841913B (zh) | 配电线路故障测距方法 | |
CN105445624A (zh) | 结合小波变换及曲线拟合的电缆故障定位方法 | |
CN101907437A (zh) | 一种基于小波差分算法的电缆故障测距方法 | |
CN111766470B (zh) | 高压直流输电线路的故障定位方法、***及直流输电线路 | |
CN110514963B (zh) | 一种改进双端行波故障测距方法 | |
CN102087332A (zh) | 一种基于波速寻优的直流行波故障测距方法 | |
CN106771861A (zh) | 基于广域行波能量和时间差的复杂电网故障定位方法 | |
CN107632236A (zh) | 一种基于对端母线反射波识别的单出线输电线路单端故障测距方法 | |
CN103412240B (zh) | 一种不依赖波头辨识的同塔双回输电线路单端行波故障测距方法 | |
CN110247420B (zh) | 一种hvdc输电线路故障智能识别方法 | |
CN109190272A (zh) | 基于弹性波和机器学习的混凝土结构缺陷检测方法 | |
CN107179473B (zh) | 一种输电线路故障定位方法 | |
CN103543378A (zh) | 台片关系检测方法 | |
CN113884818B (zh) | 一种基于lstm的配电网故障行波到达时间精确估算方法 | |
CN105445614A (zh) | 一种基于小波分析的双端行波故障定位的方法和*** | |
CN103163428B (zh) | 一种提高单端行波测距可靠性的方法 | |
CN102313858B (zh) | 一种初始反极性方向行波的辨识方法 | |
CN105260521A (zh) | 一种基于层次分析法的负荷融合建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160601 Termination date: 20180410 |
|
CF01 | Termination of patent right due to non-payment of annual fee |