CN115184718B - 一种基于时域分析的三芯电缆故障定位方法 - Google Patents
一种基于时域分析的三芯电缆故障定位方法 Download PDFInfo
- Publication number
- CN115184718B CN115184718B CN202210414945.2A CN202210414945A CN115184718B CN 115184718 B CN115184718 B CN 115184718B CN 202210414945 A CN202210414945 A CN 202210414945A CN 115184718 B CN115184718 B CN 115184718B
- Authority
- CN
- China
- Prior art keywords
- phase
- fault
- shielding layer
- metal shielding
- virtual
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 239000002184 metal Substances 0.000 claims description 71
- 239000013598 vector Substances 0.000 claims description 27
- 239000004020 conductor Substances 0.000 claims description 7
- 239000013256 coordination polymer Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 240000005373 Panax quinquefolius Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/083—Locating faults in cables, transmission lines, or networks according to type of conductors in cables, e.g. underground
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
一种基于时域分析的三芯电缆故障定位方法,首先,利用故障发生后电缆线路首端和末端测量的线芯电流、对地电压和接地线电流时域信号,计算虚拟故障点处三相线芯和屏蔽层之间的电压、三相线芯流向屏蔽层的电流;然后,估计虚拟故障点处三相线芯和屏蔽层间的虚拟故障电阻,计算虚拟故障点处三相线芯和屏蔽层之间的电压差在一段时间内的积分,确定虚拟故障点在线路全长范围内变动时积分最小值对应的虚拟故障距离;最后,通过比较虚拟故障距离处计算得到的三相虚拟故障电阻在一段时间内的均值确定故障相,该故障相对应的最小虚拟故障距离即为实际的故障距离。仿真表明该方法能准确辨识不同故障状况下的故障相和故障精确位置,具有较高的工程应用价值。
Description
技术领域
本发明涉及电网设备故障排查技术领域,具体涉及一种基于时域分析的三芯电缆故障定位方法。
背景技术
随着城市化进程的加快,电缆在城网供电中所占的份量也越来越重,已逐步取代了架空供电线路,被广泛应用于配电网。电缆线路容易因外力、环境等因素导致故障,而电缆线路大多敷设于地下,故障查找困难。为了提高运维效率、提升配网电缆故障感知水平、保证配网安全稳定运行,进行三线电缆故障定位势在必行。
目前,仅有针对中、高压单芯电缆线路的故障定位方法,而尚未由针对三芯电缆结构提出有效的故障定位方法。单芯电缆故障定位方法主要分为两类:行波法和阻抗法。行波法主要利用波头到达时刻与故障距离的关系实现定位;阻抗法则利用线路一端或者两端稳态电压、电流相量,通过求解包含故障距离、线路参数的方程(组)实现定位。文献《基于全波形信息的混联线路单端行波定位方法》(邓丰、李欣然、曾祥君,电工技术学报)构建了不同故障位置和故障状况下的行波时频能量谱矩阵,通过比较不同位置的全波波形能量相关性确定电缆-架空混联线路的故障位置。文献《利用故障特征频带和TT变换的电缆单端行波测距》(束洪春、田鑫萃、董俊、杨毅,中国电机工程学报)采用TT变换精确标定两个行波波头到达时刻,利用单端法实现单芯电缆故障精确测距。文献《计及金属护层结构的电缆单端故障测距方法》(唐进、张姝、林圣、何正友,中国电机工程学报)通过迭代正常相线芯和护层电流故障稳态分量在故障点前后的差值以确定单芯的故障精确位置。
不同于单芯电缆,三芯电缆结构复杂,三相导芯和金属屏蔽层之间的存在强电磁耦合,现有研究所提出的故障定位方法无法适用这种复杂情况。此外,由于三芯电缆线路距离短,行波波头辨识困难,研究行波法的意义不大;阻抗法虽然不需要辨识行波波头,但实际稳态故障信号有时无法获取,研究阻抗法的价值不高。因此,基于故障时域等效模型,利用故障时域信号实现三芯电缆的故障定位方法具有重要价值和意义。
发明内容
为了解决现有关于电缆故障定位研究存在的不足,提出一种基于时域分析的三芯电缆故障定位方法。
一种基于时域分析的三芯电缆故障定位方法,包括如下步骤:
步骤1,在三芯电缆发生故障后,测量线路首端三相导芯电流、导芯对地电压、接地线电流,以及线路末端三相导芯电流、导芯对地电压、屏蔽层对地电压,获取线路阻抗和接地线电阻;
步骤2,对于所有采集得到的电流、电压,截取有效时间窗[t1,t2] 内的信号;其中,t1和t2分别为有效时间窗的起始、终止时刻;
步骤3,令虚拟故障点与线路首端的距离x的初始值为0;
步骤4,将线路首端测量的电流、电压信号带入下式计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层对地电压信号uLFAG(t,x)、 uLFBG(t,x)、uLFCG(t,x)和uLFSG(t,x):
其中,uLA(t)、uLB(t)、uLC(t)和iLS(t)分别为故障发生后线路首端 A相、B相、C相导芯对地电压和接地线电流信号,iLP(t)和iLS(t)分别为故障发生后线路首端P相导芯的电流和接地线电流信号,P∈{A, B,C},RAP、RBP、RCP、RSP分别为A相、B相、C相导芯、金属屏蔽层与P相导芯之间的互电阻,RAS、RBS、RCS分别为A相、B相、C 相导芯与金属屏蔽层之间的互电阻,RSS为金属屏蔽层的自电阻,RS为接地线电阻,LAP、LBP、LCP、LSP分别为A相、B相、C相导芯、金属屏蔽层与P相导芯之间的互电感,LAS、LBS、LCS分别为A相、B 相、C相导芯与金属屏蔽层之间的互电感,LSS为金属屏蔽层的自电感;
步骤5,根据uLFAG(t,x)、uLFBG(t,x)、uLFCG(t,x)和uLFSG(t,x)计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号 uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x);
步骤6,将线路末端测量的电流、电压信号带入下式计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层对地电压信号uRFAG(t,x)、 uRFBG(t,x)、uRFCG(t,x)和uRFSG(t,x):
其中,uRA(t)、uRB(t)、uRC(t)和uS(t)分别为故障发生后线路末端 A相、B相、C相导芯对地电压和金属屏蔽层对地电压,iRP(t)为故障发生后线路末端P相导芯电流,l为线路总长;
步骤7,根据uRFAG(t,x)、uRFBG(t,x)、uRFCG(t,x)和uRFSG(t,x)计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号 uRFAS(t,x)、uRFBS(t,x)、uRFCS(t,x);
步骤8,根据下式,计算虚拟故障点x处,A相、B相、C相流向金属屏蔽层的电流信号iFAS(t,x)、iFBS(t,x)、iFCS(t,x):
步骤9,根据电压信号uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x)和电流信号iFAS(t,x)、iFBS(t,x)、iFCS(t,x)计算出虚拟故障点处A相、B相、 C相导芯与金属屏蔽层之间的虚拟故障电阻RFAS(t,x)、RFBS(t,x)、RFCS (t,x);
步骤10,将步骤5和步骤7的计算结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗[t1,t2]内的积分DA(x)、DB(x)、DC(x):
步骤11,如果虚拟故障距离x小于等于线路总长l,则令x=x+1m,重复执行步骤4至步骤10;否则,得到如下向量:
其中,RFPS(t2)表示当有效时间窗终止时刻为t2时,在不同虚拟故障距离下,P相导芯和金属屏蔽层之间的虚拟故障电阻向量,DP(t2) 表示当有效时间窗终止时刻为t2时,在不同虚拟故障距离下,P相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗内的积分向量, Nx为线路上设置的虚拟故障点的总个数,函数Fz()用于取一个实数的整数部分;
步骤12,根据步骤11的结果,确定向量RFPS(t2)的平均值和向量DP(t2)的最小值,分别记为RFPS_mean(t2)和DP_min(t2);
步骤13,如果有效时间窗的终止时刻t2小于等于有效时间窗起始时刻t1与数据窗平移总长度Tt的和,则令t1=t1+δ,t2=t2+δ,δ表示每次时间窗的平移量,重复执行步骤3至步骤12;否则,得到如下向量:
其中,RFPS表示在不同的有效时间窗终止时刻和虚拟故障距离下, P相导芯和金属屏蔽层之间的虚拟故障电阻平均值向量,DP表示不同的有效时间窗终止时刻和虚拟故障距离下,P相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗内的积分最小值向量,Nt为设置的不同有效时间窗终止时刻的总个数;
步骤14,根据步骤13的结果,确定向量RFPS的平均值,记为 RFPS_mean;
步骤15,利用如下判据辨识故障相,估计故障电阻:
其中,函数min()用于确定的RFAS_mean、RFBS_mean和RFCS_mean三个数中的最小值以及对应的索引(A,B,C);
步骤16,利用下式计算实际故障精确距离xF:
其中,k为一个连续变量且k∈Z+。
进一步地,步骤2和步骤13中,有效时间窗的起止时刻和数据窗平移总长度按照下式来确定:
其中,t0为行波首波头到达时刻,可由小波模极大值法得到;Tt为数据窗平移总长度。
进一步地,步骤5中,将步骤4的结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号:
其中,uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x)分别为利用线路首端测量的电流、电压信号计算得到的虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号。
进一步地,步骤7中,将步骤6的结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号:
其中,uRFAS(t,x)、uRFBS(t,x)、uRFCS(t,x)分别为利用线路末端测量的电流、电压信号计算得到的虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号。
进一步地,步骤9中,将步骤5和步骤8的计算结果带入下式,估计虚拟故障点处A相、B相、C相导芯与金属屏蔽层之间的虚拟故障电阻RFAS(t,x)、RFBS(t,x)、RFCS(t,x):
进一步地,步骤13中,每次时间窗的平移量δ的值由下式确定:
其中,fs为信号的采样频率。
本发明达到的有益效果为:
本方法建立了故障发生后三芯电缆的等效时域模型,通过测量线路首末端导芯电流、对地电压及接地线电流时域信号,利用虚拟故障点法构建了故障相辨识、故障电阻估算以及故障精确定位判据。本方法无需相模变换、无需获取稳态相量,只需要提取故障后一段时间内的电压、电流信号,在实现故障精确定位的同时还能够可靠辨识故障相、估算故障电阻,具有较高的工程实践意义。
附图说明
图1为本发明实施例中的三芯电缆故障定位方法流程示意图。
图2为本发明实施例中的故障后三芯电缆时域等效模型示意图。
具体实施方式
下面结合说明书附图对本发明的技术方案做进一步的详细说明。
结合附图,对本发明的技术方案进行详细说明。如图2所示,故障发生后三芯电缆模型可看作由A相、B相、C相导芯和金属屏蔽层等效自阻抗和互阻抗构成的集总参数电路。其中,L端为线路首端、 R端为线路末端,l为线路全长,x为虚拟故障距离。iLA、iLB、iLC和 iLS分别为流过线路首端A相、B相、C相导芯和接地线的电流信号, uLA、uLB、uLC分别为线路首端A相、B相和C相导芯的对地电压信号,iRA、iRB、iRC分别为流过线路末端A相、B相、C相导芯的电流信号,uRA、uRB、uRC、uS分别为线路末端A相、B相、C相导芯、金属屏蔽层的对地电压信号,RS为接地线电阻,ZAA、ZAB、ZAC、ZAS分别为A相导芯的自阻抗、A相和B相导芯的互阻抗、A相和C相导芯的互阻抗、A相导芯和金属屏蔽层的互阻抗,ZBB、ZBC、ZBS分别为B相导芯自阻抗、B相和C相导芯的互阻抗、B相和金属屏蔽层的互阻抗,ZCC和ZCS分别为C相导芯自阻抗、C相导芯和金属屏蔽层的互阻抗,ZLA、ZLB、ZLC为A相、B相和C相等效负荷。根据电路分析,三芯电缆故障测距的具体步骤如下。
步骤1,在三芯电缆发生故障后,测量线路首端三相导芯电流、导芯对地电压、接地线电流,以及线路末端三相导芯电流、导芯对地电压、屏蔽层对地电压,获取线路阻抗和接地线电阻参数。
步骤2,对于所有采集得到的电流、电压,截取有效时间窗[t1,t2] 内的信号;其中,t1和t2分别为有效时间窗的起始、终止时刻,按照下式确定:
其中,t0为暂态行波首波头到达时刻,该值可由小波模极大值法得到;Tt表示数据窗平移总长度。
步骤3,令虚拟故障点与线路首端的距离x的初始值为0。
步骤4,将线路首端测量的电流、电压信号带入下式计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层对地电压信号uLFAG(t,x)、 uLFBG(t,x)、uLFCG(t,x)和uLFSG(t,x):
其中,uLA(t)、uLB(t)、uLC(t)和iLS(t)分别为故障发生后线路首端 A相、B相、C相导芯对地电压和接地线电流信号,iLP(t)和iLS(t)分别为故障发生后线路首端P相(P∈{A,B,C})导芯的电流和接地线电流信号,RAP、RBP、RCP、RSP分别为A相、B相、C相导芯、金属屏蔽层与P相导芯之间的互电阻,RAS、RBS、RCS分别为A相、B相、 C相导芯与金属屏蔽层之间的互电阻,RSS为金属屏蔽层的自电阻, RS为接地线电阻,LAP、LBP、LCP、LSP分别为A相、B相、C相导芯、金属屏蔽层与P相导芯之间的互电感,LAS、LBS、LCS分别为A相、B 相、C相导芯与金属屏蔽层之间的互电感,LSS为金属屏蔽层的自电感。
步骤5,将步骤4的结果带入下式,计算虚拟故障点处A相、B 相、C相导芯和金属屏蔽层之间的电压信号:
其中,uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x)分别为利用线路首端测量的电流、电压信号计算得到的虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号。
步骤6,将线路末端测量的电流、电压信号带入下式计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层对地电压信号uRFAG(t,x)、 uRFBG(t,x)、uRFCG(t,x)和uRFSG(t,x):
其中,uRA(t)、uRB(t)、uRC(t)和uS(t)分别为故障发生后线路末端 A相、B相、C相导芯对地电压和金属屏蔽层对地电压,iRP(t)为故障发生后线路末端P相(P∈{A,B,C})导芯电流,l为线路总长。
步骤7,将步骤6的结果带入下式,计算虚拟故障点处A相、B 相、C相导芯和金属屏蔽层之间的电压信号:
其中,uRFAS(t,x)、uRFBS(t,x)、uRFCS(t,x)分别为利用线路末端测量的电流、电压信号计算得到的虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号。
步骤8,根据下式,计算虚拟故障点x处,A相、B相、C相流向金属屏蔽层的电流信号iFAS(t,x)、iFBS(t,x)、iFCS(t,x):
其中,iLA(t)、iLB(t)、iLC(t)为线路首端测量的A相、B相、C相导芯电流,iRA(t)、iRB(t)、iRC(t)分别为线路末端测量的A、B、C相导芯电流。
步骤9,将步骤5和步骤8的计算结果带入下式,估计虚拟故障点处A相、B相、C相导芯与金属屏蔽层之间的虚拟故障电阻RFAS(t, x)、RFBS(t,x)、RFCS(t,x):
步骤10,将步骤5和步骤7的计算结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗[t1,t2]内的积分DA(x)、DB(x)、DC(x):
步骤11,如果虚拟故障距离x小于等于线路总长l,则令x=x+1m,重复执行步骤4至步骤10;否则,得到如下向量:
其中,RFPS(t2)表示当有效时间窗终止时刻为t2时,在不同虚拟故障距离下,P(P∈{A,B,C})相导芯和金属屏蔽层之间的虚拟故障电阻向量,DP(t2)表示当有效时间窗终止时刻为t2时,在不同虚拟故障距离下,P相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗内的积分向量,Nx为线路上设置的虚拟故障点的总个数,函数Fz() 用于取一个实数的整数部分。
步骤12,根据步骤11的结果,确定向量RFPS(t2)的平均值和向量DP(t2)的最小值,分别记为RFPS_mean(t2)和DP_min(t2)。
步骤13,如果有效时间窗的终止时刻t2小于等于有效时间窗起始时刻t1与数据窗平移总长度Tt的和,则令t1=t1+δ,t2=t2+δ(δ表示每次时间窗的平移量,其值为采样频率的倒数,即采样间隔),重复执行步骤3至步骤12;否则,得到如下向量:
其中,RFPS表示在不同的有效时间窗终止时刻和虚拟故障距离下, P(P∈{A,B,C})相导芯和金属屏蔽层之间的虚拟故障电阻平均值向量,DP表示不同的有效时间窗终止时刻和虚拟故障距离下,P相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗内的积分最小值向量,Nt为设置的不同有效时间窗终止时刻的总个数,执行步骤13。
步骤14,根据步骤13的结果,确定向量RFPS的平均值,记为 RFPS_mean;
步骤15,利用如下判据辨识故障相,估计故障电阻:
其中,函数min()用于确定的RFAS_mean、RFBS_mean和RFCS_mean三个数中的最小值以及对应的索引(A,B,C)、
步骤16,利用下式计算实际故障精确距离xF:
其中,k为一个连续变量且k∈Z+。
为了验证本方法的正确性,在PSCAD上搭建一个10kV单端辐射型三芯电缆配电网故障仿真模型。三芯电缆采用依频特性相位模型。出线电缆1、2、3总长分别为:1.5km、2km和1km。在第三条电缆上模拟不同情况的故障,所得到的故障测距结果如表1所示。从表中可以看出:所提出的方法能够辨识故障电阻,实现不同故障情况下三芯电缆的故障精确定位,最大和最小定位误差分别为25m和2m。
表1不同情况下三芯电缆故障定位结果
/>
以上所述仅为本发明的较佳实施方式,本发明的保护范围并不以上述实施方式为限,但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化,皆应纳入权利要求书中记载的保护范围内。
Claims (6)
1.一种基于时域分析的三芯电缆故障定位方法,其特征在于:包括如下步骤:
步骤1,在三芯电缆发生故障后,测量线路首端三相导芯电流、导芯对地电压、接地线电流,以及线路末端三相导芯电流、导芯对地电压、屏蔽层对地电压,获取线路阻抗和接地线电阻;
步骤2,对于所有采集得到的电流、电压,截取有效时间窗[t1,t2]内的信号;其中,t1和t2分别为有效时间窗的起始、终止时刻;
步骤3,令虚拟故障点与线路首端的距离x的初始值为0;
步骤4,将线路首端测量的电流、电压信号带入下式计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层对地电压信号uLFAG(t,x)、uLFBG(t,x)、uLFCG(t,x)和uLFSG(t,x):
其中,uLA(t)、uLB(t)、uLC(t)和iLS(t)分别为故障发生后线路首端A相、B相、C相导芯对地电压和接地线电流信号,iLP(t)和iLS(t)分别为故障发生后线路首端P相导芯的电流和接地线电流信号,P∈{A,B,C},RAP、RBP、RCP、RSP分别为A相、B相、C相导芯、金属屏蔽层与P相导芯之间的互电阻,RAS、RBS、RCS分别为A相、B相、C相导芯与金属屏蔽层之间的互电阻,RSS为金属屏蔽层的自电阻,RS为接地线电阻,LAP、LBP、LCP、LSP分别为A相、B相、C相导芯、金属屏蔽层与P相导芯之间的互电感,LAS、LBS、LCS分别为A相、B相、C相导芯与金属屏蔽层之间的互电感,LSS为金属屏蔽层的自电感;
步骤5,根据uLFAG(t,x)、uLFBG(t,x)、uLFCG(t,x)和uLFSG(t,x)计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x);
步骤6,将线路末端测量的电流、电压信号带入下式计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层对地电压信号uRFAG(t,x)、uRFBG(t,x)、uRFCG(t,x)和uRFSG(t,x):
其中,uRA(t)、uRB(t)、uRC(t)和uS(t)分别为故障发生后线路末端A相、B相、C相导芯对地电压和金属屏蔽层对地电压,iRP(t)为故障发生后线路末端P相导芯电流,l为线路总长;
步骤7,根据uRFAG(t,x)、uRFBG(t,x)、uRFCG(t,x)和uRFSG(t,x)计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号uRFAS(t,x)、uRFBS(t,x)、uRFCS(t,x);
步骤8,根据下式,计算虚拟故障点x处,A相、B相、C相流向金属屏蔽层的电流信号iFAS(t,x)、iFBS(t,x)、iFCS(t,x):
步骤9,根据电压信号uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x)和电流信号iFAS(t,x)、iFBS(t,x)、iFCS(t,x)计算出虚拟故障点处A相、B相、C相导芯与金属屏蔽层之间的虚拟故障电阻RFAS(t,x)、RFBS(t,x)、RFCS(t,x);
步骤10,将步骤5和步骤7的计算结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗[t1,t2]内的积分DA(x)、DB(x)、DC(x):
步骤11,如果虚拟故障距离x小于等于线路总长l,则令x=x+1m,重复执行步骤4至步骤10;否则,得到如下向量:
其中,RFPS(t2)表示当有效时间窗终止时刻为t2时,在不同虚拟故障距离下,P相导芯和金属屏蔽层之间的虚拟故障电阻向量,DP(t2)表示当有效时间窗终止时刻为t2时,在不同虚拟故障距离下,P相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗内的积分向量,Nx为线路上设置的虚拟故障点的总个数,函数Fz()用于取一个实数的整数部分;
步骤12,根据步骤11的结果,确定向量RFPS(t2)的平均值和向量DP(t2)的最小值,分别记为RFPS_mean(t2)和DP_min(t2);
步骤13,如果有效时间窗的终止时刻t2小于等于有效时间窗起始时刻t1与数据窗平移总长度Tt的和,则令t1=t1+δ,t2=t2+δ,δ表示每次时间窗的平移量,重复执行步骤3至步骤12;否则,得到如下向量:
其中,RFPS表示在不同的有效时间窗终止时刻和虚拟故障距离下,P相导芯和金属屏蔽层之间的虚拟故障电阻平均值向量,DP表示不同的有效时间窗终止时刻和虚拟故障距离下,P相导芯和金属屏蔽层之间电压的绝对差值在有效时间窗内的积分最小值向量,Nt为设置的不同有效时间窗终止时刻的总个数;
步骤14,根据步骤13的结果,确定向量RFPS的平均值,记为RFPS_mean;
步骤15,利用如下判据辨识故障相,估计故障电阻:
其中,函数min()用于确定的RFAS_mean、RFBS_mean和RFCS_mean三个数中的最小值以及对应的索引(A,B,C);
步骤16,利用下式计算实际故障精确距离xF:
其中,k为一个连续变量且k∈Z+。
2.根据权利要求1所述的一种基于时域分析的三芯电缆故障定位方法,其特征在于:步骤2和步骤13中,有效时间窗的起止时刻和数据窗平移总长度按照下式来确定:
其中,t0为行波首波头到达时刻,可由小波模极大值法得到;Tt为数据窗平移总长度。
3.根据权利要求1所述的一种基于时域分析的三芯电缆故障定位方法,其特征在于:步骤5中,将步骤4的结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号:
其中,uLFAS(t,x)、uLFBS(t,x)、uLFCS(t,x)分别为利用线路首端测量的电流、电压信号计算得到的虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号。
4.根据权利要求1所述的一种基于时域分析的三芯电缆故障定位方法,其特征在于:步骤7中,将步骤6的结果带入下式,计算虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号:
其中,uRFAS(t,x)、uRFBS(t,x)、uRFCS(t,x)分别为利用线路末端测量的电流、电压信号计算得到的虚拟故障点处A相、B相、C相导芯和金属屏蔽层之间的电压信号。
5.根据权利要求1所述的一种基于时域分析的三芯电缆故障定位方法,其特征在于:步骤9中,将步骤5和步骤8的计算结果带入下式,估计虚拟故障点处A相、B相、C相导芯与金属屏蔽层之间的虚拟故障电阻RFAS(t,x)、RFBS(t,x)、RFCS(t,x):
6.根据权利要求1所述的一种基于时域分析的三芯电缆故障定位方法,其特征在于:步骤13中,每次时间窗的平移量δ的值由下式确定:
其中,fs为信号的采样频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210414945.2A CN115184718B (zh) | 2022-04-20 | 2022-04-20 | 一种基于时域分析的三芯电缆故障定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210414945.2A CN115184718B (zh) | 2022-04-20 | 2022-04-20 | 一种基于时域分析的三芯电缆故障定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115184718A CN115184718A (zh) | 2022-10-14 |
CN115184718B true CN115184718B (zh) | 2023-11-14 |
Family
ID=83512176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210414945.2A Active CN115184718B (zh) | 2022-04-20 | 2022-04-20 | 一种基于时域分析的三芯电缆故障定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115184718B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117092453B (zh) * | 2023-10-19 | 2023-12-19 | 华南理工大学 | 三芯电缆的故障定位方法、装置、设备及存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1997034160A1 (en) * | 1996-03-11 | 1997-09-18 | Electric Power Research Institute, Inc. | Locating a fault in an electrical power cable |
JP2011017637A (ja) * | 2009-07-09 | 2011-01-27 | Kansai Electric Power Co Inc:The | ケーブル地絡箇所検出方法およびケーブル地絡箇所検出装置 |
CN106443329A (zh) * | 2016-09-07 | 2017-02-22 | 南京南瑞继保电气有限公司 | 一种直流接地极线路故障检测***与方法 |
CN106443351A (zh) * | 2016-11-02 | 2017-02-22 | 南京工程学院 | 一种船舶多芯电缆的故障检测***和检测方法 |
CN108693446A (zh) * | 2018-05-25 | 2018-10-23 | 中国矿业大学 | 一种非同步采样电网暂态行波模量时间差的故障定位方法 |
CN112083271A (zh) * | 2020-08-18 | 2020-12-15 | 昆明理工大学 | 一种基于护层电流行波时频复合分析的10kV电缆在线测距方法 |
CN113093050A (zh) * | 2021-03-31 | 2021-07-09 | 中国矿业大学 | 电缆接地线电流时频特征的电缆早期故障辨识方法及*** |
CN113655342A (zh) * | 2021-09-23 | 2021-11-16 | 国网江苏省电力有限公司连云港供电分公司 | 一种基于多导体耦合模型的三芯电缆单相故障定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3171185B1 (de) * | 2015-11-20 | 2018-06-27 | Siemens Aktiengesellschaft | Verfahren und einrichtung zum bestimmen des fehlerortes im falle eines fehlers auf einer elektrischen leitung |
KR102180748B1 (ko) * | 2019-08-21 | 2020-11-19 | 연세대학교 산학협력단 | 멀티코어 케이블 결함 분류 시스템 및 방법 |
-
2022
- 2022-04-20 CN CN202210414945.2A patent/CN115184718B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1997034160A1 (en) * | 1996-03-11 | 1997-09-18 | Electric Power Research Institute, Inc. | Locating a fault in an electrical power cable |
JP2011017637A (ja) * | 2009-07-09 | 2011-01-27 | Kansai Electric Power Co Inc:The | ケーブル地絡箇所検出方法およびケーブル地絡箇所検出装置 |
CN106443329A (zh) * | 2016-09-07 | 2017-02-22 | 南京南瑞继保电气有限公司 | 一种直流接地极线路故障检测***与方法 |
CN106443351A (zh) * | 2016-11-02 | 2017-02-22 | 南京工程学院 | 一种船舶多芯电缆的故障检测***和检测方法 |
CN108693446A (zh) * | 2018-05-25 | 2018-10-23 | 中国矿业大学 | 一种非同步采样电网暂态行波模量时间差的故障定位方法 |
CN112083271A (zh) * | 2020-08-18 | 2020-12-15 | 昆明理工大学 | 一种基于护层电流行波时频复合分析的10kV电缆在线测距方法 |
CN113093050A (zh) * | 2021-03-31 | 2021-07-09 | 中国矿业大学 | 电缆接地线电流时频特征的电缆早期故障辨识方法及*** |
CN113655342A (zh) * | 2021-09-23 | 2021-11-16 | 国网江苏省电力有限公司连云港供电分公司 | 一种基于多导体耦合模型的三芯电缆单相故障定位方法 |
Non-Patent Citations (11)
Title |
---|
Comparison of Location Accuracy between Frequency Domain Reflectometry and Line Resonance Analysis for Power Cables;Lu Chen等;2020 IEEE International Conference on High Voltage Engineering and Application (ICHVE);全文 * |
Permanent Fault Location for Three-Core Cable Using Multiconductor Coupling Model;Zehua Tang等;IEEE Transactions on Power Delivery;全文 * |
噪声时域反射法实现多根电缆故障的同时检测;马培国等;激光与光电子学进展;全文 * |
基于TDR多芯信号电缆故障测试装置的研制;吕瀚;吴雄伟;关伟智;吴林斌;;大地测量与地球动力学(S1);全文 * |
基于双层阻抗模型的三相单芯电缆自恢复故障定位;张姝;林圣;唐进;何正友;;电工技术学报(17);全文 * |
基于多源暂稳态信息的电缆型配电网故障辨识与定位研究;彭楠;中国博士学位论文全文数据库工程科技Ⅱ辑;全文 * |
基于时域反射的多芯线缆检测装置设计;汪兴海;计算机与数字工程;全文 * |
基于暂态特征模量分析的 配网三芯铠装电缆故障感知与测距;彭楠等;中国电机工程学报;全文 * |
矿山电网电弧性接地故障区段定位与辨识研究;程成;中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑;全文 * |
考虑金属护套和铠装结构的电缆故障测距方法;宋勤;王媛媛;;电气技术(10);全文 * |
计及金属护层结构的电缆单端故障测距方法;唐进;张姝;林圣;何正友;;中国电机工程学报(06);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115184718A (zh) | 2022-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Suonan et al. | A novel fault-location method for HVDC transmission lines | |
CN103792465B (zh) | 一种基于零序电压的配电网单相接地故障测距的方法 | |
CN102175954B (zh) | 一种线路相间故障单端测距方法 | |
CN110082636B (zh) | 一种电力电缆故障定位方法及*** | |
CN109633367B (zh) | 基于故障前后电压、电流变化的输电线路故障定位方法 | |
CN103941147B (zh) | 利用暂态主频分量的配网电缆单相接地故障测距方法 | |
CN102135571B (zh) | 超高压/特高压多回路输电线路零序阻抗抗干扰测量方法 | |
CN111308264B (zh) | 一种基于余弦相似度的配电网单相接地故障区段定位方法 | |
Qi et al. | Calculation of interference voltage on the nearby underground metal pipeline due to the grounding fault on overhead transmission lines | |
CN109100605B (zh) | 利用故障边界条件的高压电缆单相接地故障的单端定位方法 | |
CN111638423B (zh) | 一种电力电缆护套层和铠装层接地故障的定位方法 | |
Zhang et al. | A new fault-location algorithm for extra-high-voltage mixed lines based on phase characteristics of the hyperbolic tangent function | |
CN115184718B (zh) | 一种基于时域分析的三芯电缆故障定位方法 | |
CN113655342A (zh) | 一种基于多导体耦合模型的三芯电缆单相故障定位方法 | |
Cao et al. | Asynchronous fault location scheme based on voltage distribution for three-terminal transmission lines | |
Zhang et al. | Voltage-sag-profiles-based fault location in high-speed railway distribution system | |
CN112526290A (zh) | 一种基于广域行波侧后模拟的复杂电网接地故障定位方法 | |
Qin et al. | A novel pilot protection scheme for HVDC lines based on waveform matching | |
CN103424627B (zh) | 双端测量平行电网线路零序阻抗的方法 | |
Zhang et al. | Fault location method for three-core cable using amplitude ratio of shield-grounding wire currents | |
CN112904233B (zh) | 基于地线电流无功分量的接地故障定位方法和装置 | |
CN114859168A (zh) | 一种基于双端非同步测量的三芯电缆故障测距方法 | |
CN113805012B (zh) | 适用于中性点经小电阻接地电缆配电网故障区段辨识方法 | |
Guo et al. | A method of 20 kV cable line fault location based on sheath grounding current | |
CN103424628A (zh) | 测量平行电网线路正序阻抗的方法 |
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 |