CN108008246B - 基于时域自适应窗的线路故障测距方法及装置 - Google Patents

基于时域自适应窗的线路故障测距方法及装置 Download PDF

Info

Publication number
CN108008246B
CN108008246B CN201711175222.7A CN201711175222A CN108008246B CN 108008246 B CN108008246 B CN 108008246B CN 201711175222 A CN201711175222 A CN 201711175222A CN 108008246 B CN108008246 B CN 108008246B
Authority
CN
China
Prior art keywords
fault
time
traveling wave
terminal
intelligent monitoring
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
CN201711175222.7A
Other languages
English (en)
Other versions
CN108008246A (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.)
HUNAN XIANGNENG SMART ELECTRICAL EQUIPMENT CO Ltd
Original Assignee
HUNAN XIANGNENG SMART ELECTRICAL EQUIPMENT CO Ltd
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 HUNAN XIANGNENG SMART ELECTRICAL EQUIPMENT CO Ltd filed Critical HUNAN XIANGNENG SMART ELECTRICAL EQUIPMENT CO Ltd
Priority to CN201711175222.7A priority Critical patent/CN108008246B/zh
Publication of CN108008246A publication Critical patent/CN108008246A/zh
Application granted granted Critical
Publication of CN108008246B publication Critical patent/CN108008246B/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
    • 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/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • 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/088Aspects of digital computing
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage 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

基于时域自适应窗的线路故障测距方法及装置
技术领域
本发明涉及电力系统技术领域,尤其涉及一种基于时域自适应窗的线路故障测距方法及装置。
背景技术
智能电网的发展,构建更多的高压、长距离输电线路,电力行业非常重视输电线路故障测距难题。当输电线路故障发生时,为减少巡线时间,缩短停电时间,提高供电可靠性,迫切需要能及时准确地查找故障位置的适用方法。
对比电力系统传统的阻抗法,行波定位技术具有定位精度高、不受线路阻抗影响等优点。目前应用较多采用变电站内的行波监测装置,利用两端法或网络法定位故障点,但行波长距离的传输的衰减和畸变会给波头的捕捉带来影响,因此,为提高故障定位精确度,亟需一种新的定位方法。
发明内容
本发明提供一种基于时域自适应窗的线路故障测距方法及装置,用以解决现有技术中的测距方法对故障定位精确度不高的技术问题。
本发明一方面提供一种基于时域自适应窗的线路故障测距方法,输电线路上预先设置有监测点,每个监测点处的A相、B相和C相上分别安装有智能监测终端,包括:
步骤101,通过智能监测终端采集第一工频电流;
步骤102,对智能监测终端采集的所述第一工频电流进行分析,以判定所述智能监测终端所位于的线路是否发生故障,若是,转步骤103执行;
步骤103,通过位于故障相上的智能监测终端采集故障电流行波,以获得故障电流行波录波文件,所述故障相为线路发生故障的某相,所述故障电流行波录波文件包括故障电流行波录波波形;
步骤104,对所述故障电流行波录波文件进行分析,以获得故障电流行波录波波形的过零时间点,其中,相邻所述过零时间点为一个时域自适应窗;
步骤105,计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗;
步骤106,在有效行波阵时间窗内,以η为固定时间窗,以采样时间T为滑窗时间间隔,根据表达式
Figure BDA0001478090570000021
分别计算有效行波阵时间窗内的各个滑窗的单位行波能量,并以值最大的行波能量所对应的滑动时间窗的最早时间点为行波首波波头时间点,其中,△Tm为有效行波阵时间窗,k1∈{0,1,2,…,N-1},滑窗次数N=(△Tm-η)/T,
Figure BDA0001478090570000022
表示第k1个过零时间点处的时间,i(t1)表示t1时刻的电流值,
Figure BDA0001478090570000023
步骤107,根据计算出的行波首波波头时间点,基于双端定位进行拟合测距。
进一步的,步骤105具体包括:根据表达式
Figure BDA0001478090570000024
计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所述对应的时域自适应窗作为有效行波阵时间窗;
其中,k-1∈{0,1,2,…,n-1},n为过零时间点总个数,i(t)表示t时刻的电流值,t∈(Tk-1,Tk),Tk-1表示第K-1个过零时间点处的时间,Tk表示与第K-1个过零时间点相邻的第K个过零时间点处的时间,Tk>Tk-1,Tk-Tk-1为一个时域自适应窗。
进一步的,步骤102具体包括:
对所述智能监测终端采集的第一工频电流进行分析,判断第一工频电流是否在预设时间段内增加预设增量后减少到零,若是,则所述智能监测终端所位于的线路为故障相,转步骤103执行。
进一步的,在判定所述智能监测终端所位于的线路发生故障后,还包括步骤103’,识别所述故障相上的故障区域;
步骤103’具体包括:根据故障相上的所述智能监测终端采集的故障工频电流计算相位;
判断相邻两个智能监测终端采集的故障工频电流相位是否相反,若相反,以该两个智能监测终端中的一个智能监测终端为第一终端,另一个智能监测终端为第二终端,判断位于第一终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第一终端采集的故障工频电流相位相同,且位于第二终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第二终端采集的故障工频电流相位相同,其中,第一终端第一侧为第一终端的不与第二终端相邻的一侧;第二终端第一侧为第二终端的不与第一终端相邻的一侧;若判断结果为是,则第一终端与第二终端之间的区域为故障区域。
进一步的,步骤107具体包括:
将第一终端及位于第一终端第一侧的所有智能监测终端为节点归集为集合G1,将第二终端及位于第二终端第一侧的所有智能监测终端为节点归集为集合G2
在集合G1与集合G2的节点中,将集合G1中所包括的a个节点以距离故障点近远顺序依次与集合G2中所包括的b个节点以距离故障点近远顺序两两进行双端定位,以获得定位结果;其中,所述故障点位于故障区域中;
根据表达式Lab=-ΔTab×ν+2daF对所述定位结果进行拟合,选择故障位置重合次数最多的结果作为故障定位结果;其中:
Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
ΔTab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;
daF为集合G1的第a个节点到故障点F的距离;
ν为输电线路的行波波速度。
本发明另一方面提供一种基于时域自适应窗的线路故障测距装置,输电线路上预先设置有监测点,每个监测点处的A相、B相和C相上分别安装有智能监测终端,包括:
第一采集模块,用于通过智能监测终端采集第一工频电流;
第一判断模块,用于对智能监测终端采集的所述第一工频电流进行分析,以判定所述智能监测终端所位于的线路是否发生故障,若是,触发第二采集模块;
第二采集模块,用于通过位于故障相上的智能监测终端采集故障电流行波,以获得故障电流行波录波文件,所述故障相为线路发生故障的某相,所述故障电流行波录波文件包括故障电流行波录波波形;
过零时间点获取模块,用于对各所述故障电流行波录波文件进行分析,以获得故障电流行波录波波形的过零时间点,其中,相邻所述过零时间点为一个时域自适应窗;
有效行波阵时间窗获取模块,用于计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗;
第一计算模块,用于在有效行波阵时间窗内,以η为固定时间窗,以采样时间T为滑窗时间间隔,根据表达式
Figure BDA0001478090570000041
分别计算有效行波阵时间窗内的各个滑窗的单位行波能量,并以值最大的行波能量所对应的时间窗的最早时间点为行波首波波头时间点,其中,△Tm为有效行波阵时间窗,k1∈{0,1,2,…,N-1},滑窗次数N=(△Tm-η)/T,
Figure BDA0001478090570000042
表示第k1个过零时间点处的时间,i(t1)表示t1时刻的电流值,
Figure BDA0001478090570000043
第二计算模块,根据计算出的行波首波波头时间点,基于双端定位进行拟合测距。
进一步的,有效行波阵时间窗获取模块具体用于:
根据表达式
Figure BDA0001478090570000044
计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所述对应的时域自适应窗作为有效行波阵时间窗;
其中,k-1∈{0,1,2,…,n-1},n为过零时间点总个数,i(t)表示t时刻的电流值,t∈(Tk-1,Tk),Tk-1表示第K-1个过零时间点处的时间,Tk表示与第K-1个过零时间点相邻的第K个过零时间点处的时间,Tk>Tk-1,Tk-Tk-1为一个时域自适应窗。
进一步的,第一判断模块具体用于:
对所述智能监测终端采集的第一工频电流进行分析,判断第一工频电流是否在预设时间段内增加预设增量后减少到零,若是,则所述智能监测终端所位于的线路为故障相,触发第二采集模块。
进一步的,上述装置还包括第二判断模块,所述第二判断模块具体用于:
根据故障相上的所述智能监测终端采集的故障工频电流计算相位;
判断相邻两个智能监测终端采集的故障工频电流相位是否相反,若相反,以该两个智能监测终端中的一个智能监测终端为第一终端,另一个智能监测终端为第二终端,判断位于第一终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第一终端采集的故障工频电流相位相同,且位于第二终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第二终端采集的故障工频电流相位相同,其中,第一终端第一侧为第一终端的不与第二终端相邻的一侧;第二终端第一侧为第二终端的不与第一终端相邻的一侧;若判断结果为是,则第一终端与第二终端之间的区域为故障区域。
进一步的,所述第二计算模块具体用于:
将第一终端及位于第一终端第一侧的所有智能监测终端作为节点归集到集合G1中,将第二终端及位于第二终端第一侧的所有智能监测终端作为节点归集到集合G2中;
在集合G1与集合G2的节点中,将集合G1中所包括的a个节点以距离故障点近远顺序依次与集合G2中所包括的b个节点以距离故障点近远顺序两两进行双端定位,以获得定位结果;其中,所述故障点位于故障区域中;
根据表达式Lab=-ΔTab×ν+2daF对所述定位结果进行拟合,选择故障位置重合次数最多的结果作为故障定位结果;
其中:Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
ΔTab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;
daF为集合G1的第a个节点到故障点F的距离;
ν为输电线路的行波波速度。
本发明提供一种基于时域自适应窗的线路故障测距方法及装置,通过对智能监测终端采集的第一工频电流进行分析,以判定某相是否发生故障,然后对位于故障相上的智能监测终端采集的故障电流行波录波波形进行分析,以获得过零时间点,从而获得时域自适应窗,再计算单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗,最后计算获得有效行波阵时间窗内的行波能量,并以值最大的行波能量所对应的滑动时间窗的最早时间点为行波首波波头时间点,然后根据计算出的各行波首波波头时间点,基于双端定位进行拟合测距,从而确定故障定位结果,进一步提高定位精度。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中:
图1为本发明实施例提供的基于时域自适应窗的线路故障测距方法的一流程示意图;
图2为本发明实施例提供的输电线路C相发生接地故障示意图;
图3为本发明实施例提供的基于时域自适应窗的线路故障测距装置的一结构示意图。
在附图中,相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
实施例一
图1为本发明实施例提供的基于时域自适应窗的线路故障测距方法的一流程示意图;如图1所示,本实施例提供一种基于时域自适应窗的线路故障测距方法,输电线路上预先设置有监测点,每个监测点处的A相、B相和C相上分别安装有智能监测终端,包括步骤101-步骤107。
其中,步骤101,通过智能监测终端采集第一工频电流。在输电线路上每隔5公里至50公里选择一个监测点,每个监测点的A相、B相和C相上分别安装一台智能监测终端。每台智能监测终端实时采集所在线路的第一工频电流。此次第一工频电流即为工频电流,使用“第一”做限定是为了与后续出现的故障工频电流做区分,并未有特殊含义。
步骤102,对智能监测终端采集的第一工频电流进行分析,以判定智能监测终端所位于的线路是否发生故障,若是,转步骤103执行。
具体的,对各智能监测终端采集的实时第一工频电流进行分析,判断第一工频电流是否在预设时间段内增加预设增量后减少到零,若是,则智能监测终端所位于的线路为故障相,转步骤103执行。预设时间段和预设增量可分别根据实际情况进行设置,在此不做限定,一般的,预设增量为30%。
步骤103,通过位于故障相上的智能监测终端采集故障电流行波,以获得故障电流行波录波文件,故障相为线路发生故障的某相,故障电流行波录波文件包括故障电流行波录波波形。
具体的,故障相包括A相、B相和C相中的一相或多相,在确定故障相后,需要对故障相上的故障电流行波做采集,以便后续分析。
步骤104,对故障电流行波录波文件进行分析,以获得故障电流行波录波波形的过零时间点,其中,相邻过零时间点为一个时域自适应窗。
具体的,过零时间点即故障电流行波录波波形中,纵坐标为零时的点,相邻的过零时间点为一个时域自适应窗,即相邻的过零时间点的横坐标(横坐标为时间)之差(值大的时间减去值小的时间)的值为时域自适应窗的大小。
进一步的,可对上述获得的过零时间点进行进一步优化,若获得的过零时间点为T0、T1、T2、…、Tn(n+1为过零时间点总个数),设置过零时间点有效的边界条件为△Tk=Tk-Tk-1≥δ(δ的单位为微秒),其中,k-1∈{0,1,2,…,n-1},δ可根据实际情况进行设置,对获得的过零时间点采用上述边界条件进行筛选,不满足上述边界条件的过零时间点视为无效点,应删除,刷选完成之后未删除的过零时间点为有效过零时间点,相邻的有效过零时间点组成一个时域自适应窗,该时域自适应窗参与到后续计算处理中。
步骤105,计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量对应的时域自适应窗作为有效行波阵时间窗。
步骤106,在有效行波阵时间窗内,以η为固定时间窗,以采样时间T为滑窗时间间隔,根据表达式
Figure BDA0001478090570000081
分别计算有效行波阵时间窗内的各个滑窗的单位行波能量,并以值最大的行波能量所对应的滑动时间窗的最早时间点为行波首波波头时间点,其中,△Tm为有效行波阵时间窗,k1∈{0,1,2,…,N-1},滑窗次数N=(△Tm-η)/T,
Figure BDA0001478090570000082
表示第k1个过零时间点处的时间,i(t1)表示t1时刻的电流值,
Figure BDA0001478090570000083
滑动时间窗即
Figure BDA0001478090570000084
组成的时间区间。
步骤107,根据计算出的各行波首波波头时间点,基于双端定位进行拟合测距。
本实施例提供的基于时域自适应窗的线路故障测距方法,通过对智能监测终端采集的第一工频电流进行分析,以判定某相是否发生故障,然后对位于故障相上的智能监测终端采集的故障电流行波录波波形进行分析,以获得过零时间点,从而获得时域自适应窗,再计算单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗,最后计算获得有效行波阵时间窗内的行波能量,并以值最大的行波能量所对应的滑动时间窗的最早时间点为行波首波波头时间点,然后根据计算出的各行波首波波头时间点,基于双端定位进行拟合测距,从而确定故障定位结果,进一步提高定位精度。
在本发明另一个具体实施例中,在上述实施例的基础上,步骤105具体包括:根据表达式
Figure BDA0001478090570000085
计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量对应的时域自适应窗作为有效行波阵时间窗;
其中,k-1∈{0,1,2,…,n-1},n为过零时间点总个数,i(t)表示t时刻的电流值,t∈(Tk-1,Tk),Tk-1表示第K-1个过零时间点处的时间,Tk表示与第K-1个过零时间点相邻的第K个过零时间点处的时间,Tk>Tk-1,Tk-Tk-1为一个时域自适应窗。
进一步的,在判定智能监测终端所位于的线路发生故障后,还包括步骤103’,识别故障相上的故障区域;步骤103’具体包括:根据故障相上的智能监测终端采集的故障工频电流计算相位;
判断相邻两个智能监测终端采集的故障工频电流相位是否相反,若相反,以该两个智能监测终端中的一个智能监测终端为第一终端,另一个智能监测终端为第二终端,判断位于第一终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第一终端采集的故障工频电流相位相同,且位于第二终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第二终端采集的故障工频电流相位相同,其中,第一终端第一侧为第一终端的不与第二终端相邻的一侧;第二终端第一侧为第二终端的不与第一终端相邻的一侧;若判断结果为是,则第一终端与第二终端之间的区域为故障区域。
具体的,故障工频电流即在监测到故障发生后智能监测终端采集到的工频电流。
进一步的,步骤107具体包括:
将第一终端及位于第一终端第一侧的所有智能监测终端为节点归集为集合G1,将第二终端及位于第二终端第一侧的所有智能监测终端为节点归集为集合G2
在集合G1与集合G2的节点中,将集合G1中所包括的a个节点以距离故障点近远顺序依次与集合G2中所包括的b个节点以距离故障点近远顺序两两进行双端定位,以获得定位结果;其中,故障点位于故障区域中;
根据表达式Lab=-ΔTab×ν+2daF对定位结果进行拟合,选择故障位置重合次数最多的结果作为故障定位结果;其中:
Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
ΔTab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;
daF为集合G1的第a个节点到故障点F的距离;
ν为输电线路的行波波速度。
具体的,Lab=-ΔTab×ν+2daF,可转化为y=at+b;其中:
y=Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
t=ΔTab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;
b/2=为集合G1的第a个节点到故障点F的距离;
a=ν为输电线路的行波波速度。
以下列举具体实施例对上述过程进行详细说明。
步骤1:采集输电线路第一工频电流,或者在故障发生时采集故障相上的故障工频电流和故障电流行波。
在输电线路每隔5公里至50公里选择一个监测点,每个监测点的A相、B相和C相上分别安装一台智能监测终端。每台智能监测终端采集故障时段带时标的故障工频电流和故障电流行波。智能监测终端包括电源模块、采集处理模块、守时模块和通讯模块。
步骤2:识别输电线路故障相和故障区域。
故障相为故障发生的某相或某几相线路,故障区域为故障点所在的相邻两个监测点间的区域。识别故障相的依据为监测点对应的三相上的智能监测终端中,某台终端的工频电流突增30%后减少到零,则此智能监测终端所安装的线路为故障相。故障区域识别的依据为相邻二个监测点在故障后的故障电流相位相反。如图2所示,监测点1、监测点2、监测点3和监测点的A相、B相的故障后工频电流变化很小,则A相、B相为非故障相,C相电流突增66%,则C相为故障相。监测点1和监测点2的C相工频故障电流相位同相,监测点1和监测点2间的区域不是故障区域;监测点3和监测点4的C相工频故障电流相位同相,监测点3和监测点4间的区域不是故障区域;监测点2和监测点3的C相工频故障电流相位相反,监测点2和监测点3间的区域为故障区域。
步骤3:获取自适应时间窗。
实施步骤S3时,求解所有监测点的故障相故障电流行波波形的过零时间点为T0、T1、T2、…、Tn,过零时间点有效的边界条件为△Tk=Tk-Tk-1≥2微秒,否则过零点为无效点,应删除,形成时域自适应窗。如图3所示,故障行波电流波形中过零时间点为T0、T1、T2、…、T21,时间差均大于2微秒,均为有效过零时间点。
步骤4:能量法寻找有效行波阵时间窗。
在各个相邻的时域自适应窗中,利用表达式
Figure BDA0001478090570000111
计算单位时间的电流行波能量,求得电流行波单位能量为Q1、Q2、…、Qn,求取单位时间的最大电流行波能量为Qmax=max{Q1、Q2、…、Qn},Qmax对应的时间窗为有效行波阵时间窗△Tm。对某一个监测点故障相的故障电流行波单位能量运用上式求解,得到故障相的故障电流行波单位能量值,如表1所示,由表1可知,Q1=488为最大值,T0=162μs,T1=204μs,因此,有效行波阵时间窗为△T1=T1-T0=42μs。
表1
Figure BDA0001478090570000112
步骤5:能量法捕捉行波首波波头时间点。
如图3所示,在T0至T1时间窗内,以η=1μs为固定时间窗,以采样时间T=0.2μs为滑窗,滑窗次数为(42-1)/0.2=205,利用表达式
Figure BDA0001478090570000121
计算行波能量Q1’、Q2’、…、Q205’,其中Q1’=778,Q1’为最大值,因此,T0=162μs为行波首波波头时间点。
步骤6:基于双端定位的拟合测距。
如图2所示,确定的输电线路的故障区域在监测点2到监测点3之间,集合G1={监测点2,监测点1},集合G2={监测点3,监测点4}。利用双端定位法,监测点2与监测点3定位结果为d2f,监测点2与监测点4定位结果为d2f+α,α为定位误差;监测点1与监测点3定位结果为d1f,监测点1与监测点4定位结果为d1f+β,β为定位误差。
根据一元一次方程:y=at+b,设距离y的单位为km,时间t的单位为us,斜率a的单位为km/us,距离b的单位为km,在由距离y和时间t组成的直角坐标系中,对集合y={L23、L24、L13、L14}和集合t={(t2-t3)、(t2-t4)、(t1-t3)、(t1-t4)},分别进行绘点,选中拟合一元一次方程最多的点所在的直线,获得该直线与y轴的交点b,从而获得故障距离d=b/2,d即为故障点位置。
实施例二
本实施例为装置实施例,用于执行上述实施例中的方法。
如图3所示,本实施例提供一种基于时域自适应窗的线路故障测距装置,输电线路上预先设置有监测点,每个监测点处的A相、B相和C相上分别安装有智能监测终端,包括:
第一采集模块201,用于通过智能监测终端采集第一工频电流;
第一判断模块202,用于对智能监测终端采集的第一工频电流进行分析,以判定智能监测终端所位于的线路是否发生故障,若是,触发第二采集模块;
第二采集模块203,用于通过位于故障相上的智能监测终端采集故障电流行波,以获得故障电流行波录波文件,故障相为线路发生故障的某相,故障电流行波录波文件包括故障电流行波录波波形;
过零时间点获取模块204,用于对各故障电流行波录波文件进行分析,以获得故障电流行波录波波形的过零时间点,其中,相邻过零时间点为一个时域自适应窗;
有效行波阵时间窗获取模块205,用于计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗;
第一计算模块206,用于在有效行波阵时间窗内,以η为固定时间窗,以采样时间T为滑窗时间间隔,根据表达式
Figure BDA0001478090570000131
分别计算有效行波阵时间窗内的各个滑窗的单位行波能量,并以值最大的行波能量所对应的滑动时间窗的最早时间点为行波首波波头时间点,其中,△Tm为有效行波阵时间窗,k1∈{0,1,2,…,N-1},滑窗次数N=(△Tm-η)/T,
Figure BDA0001478090570000132
表示第k1个过零时间点处的时间,i(t1)表示t1时刻的电流值,
Figure BDA0001478090570000133
第二计算模块207,根据计算出的各行波首波波头时间点,基于双端定位进行拟合测距。
进一步的,有效行波阵时间窗获取模块205具体用于:
根据表达式
Figure BDA0001478090570000134
计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量对应的时域自适应窗作为有效行波阵时间窗;
其中,k-1∈{0,1,2,…,n-1},n为过零时间点总个数,i(t)表示t时刻的电流值,t∈(Tk-1,Tk),Tk-1表示第K-1个过零时间点处的时间,Tk表示与第K-1个过零时间点相邻的第K个过零时间点处的时间,Tk>Tk-1,Tk-Tk-1为一个时域自适应窗。
进一步的,第一判断模块202具体用于:
对智能监测终端采集的第一工频电流进行分析,判断第一工频电流是否在预设时间段内增加预设增量后减少到零,若是,则智能监测终端所位于的线路为故障相,触发第二采集模块。
进一步的,上述装置还包括第二判断模块(图中未示出),第二判断模块具体用于:
根据故障相上的智能监测终端采集的故障工频电流计算相位;
判断相邻两个智能监测终端采集的故障工频电流相位是否相反,若相反,以该两个智能监测终端中的一个智能监测终端为第一终端,另一个智能监测终端为第二终端,判断位于第一终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第一终端采集的故障工频电流相位相同,且位于第二终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第二终端采集的故障工频电流相位相同,其中,第一终端第一侧为第一终端的不与第二终端相邻的一侧;第二终端第一侧为第二终端的不与第一终端相邻的一侧;若判断结果为是,则第一终端与第二终端之间的区域为故障区域。
进一步的,第二计算模块207具体用于:
将第一终端及位于第一终端第一侧的所有智能监测终端作为节点归集到集合G1中,将第二终端及位于第二终端第一侧的所有智能监测终端作为节点归集到集合G2中;
在集合G1与集合G2的节点中,将集合G1中所包括的a个节点以距离故障点近远顺序依次与集合G2中所包括的b个节点以距离故障点近远顺序两两进行双端定位,以获得定位结果;其中,故障点位于故障区域中;
根据表达式Lab=-ΔTab×ν+2daF对定位结果进行拟合,选择故障位置重合次数最多的结果作为故障定位结果;
其中:Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
ΔTab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;
daF为集合G1的第a个节点到故障点F的距离;
ν为输电线路的行波波速度。
由于实施例二是与实施例一对应的装置实施例,具体可参见实施例一种的记载,在此不做赘述。
应该理解的是,本发明所公开的实施例不限于这里所公开的特定结构、处理步骤或材料,而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是,在此使用的术语仅用于描述特定实施例的目的,而并不意味着限制。
说明书中提到的“一个实施例”或“实施例”意指结合实施例描述的特定特征、结构或特性包括在本发明的至少一个实施例中。因此,说明书通篇各个地方出现的短语“一个实施例”或“实施例”并不一定均指同一个实施例。
虽然上述示例用于说明本发明在一个或多个应用中的原理,但对于本领域的技术人员来说,在不背离本发明的原理和思想的情况下,明显可以在形式上、用法及实施的细节上作各种修改而不用付出创造性劳动。因此,本发明由所附的权利要求书来限定。

Claims (10)

1.一种基于时域自适应窗的线路故障测距方法,输电线路上预先设置有监测点,每个监测点处的A相、B相和C相上分别安装有智能监测终端,其特征在于,包括:
步骤101,通过智能监测终端采集第一工频电流;
步骤102,对智能监测终端采集的所述第一工频电流进行分析,以判定所述智能监测终端所位于的线路是否发生故障,若是,转步骤103执行;
步骤103,通过位于故障相上的智能监测终端采集故障电流行波,以获得故障电流行波录波文件,所述故障相为线路发生故障的某相,所述故障电流行波录波文件包括故障电流行波录波波形;
步骤104,对所述故障电流行波录波文件进行分析,以获得故障电流行波录波波形的过零时间点,其中,相邻所述过零时间点为一个时域自适应窗;
步骤105,计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗;
步骤106,在有效行波阵时间窗内,以η为固定时间窗,以采样时间T为滑窗时间间隔,根据表达式
Figure FDA0002381356340000011
分别计算有效行波阵时间窗内的各个滑窗的单位行波能量,并以值最大的行波能量所对应的滑动时间窗的最早时间点为行波首波波头时间点,其中,k1∈{0,1,2,…,N-1},滑窗次数N=(△Tm-η)/T,△Tm为有效行波阵时间窗,
Figure FDA0002381356340000012
表示第k1个过零时间点处的时间,i(t1)表示t1时刻的电流值,
Figure FDA0002381356340000013
步骤107,根据计算出的行波首波波头时间点,基于双端定位进行拟合测距。
2.根据权利要求1所述的方法,其特征在于,步骤105具体包括:根据表达式
Figure FDA0002381356340000014
计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所述对应的时域自适应窗作为有效行波阵时间窗;
其中,k-1∈{0,1,2,…,n-1},n为过零时间点总个数,i(t)表示t时刻的电流值,t∈(Tk-1,Tk),Tk-1表示第K-1个过零时间点处的时间,Tk表示与第K-1个过零时间点相邻的第K个过零时间点处的时间,Tk>Tk-1,Tk-Tk-1为一个时域自适应窗。
3.根据权利要求1所述的方法,其特征在于,步骤102具体包括:
对所述智能监测终端采集的第一工频电流进行分析,判断第一工频电流是否在预设时间段内增加预设增量后减少到零,若是,则所述智能监测终端所位于的线路为故障相,转步骤103执行。
4.根据权利要求1-3中任一项所述的方法,其特征在于,在判定所述智能监测终端所位于的线路发生故障后,还包括步骤103’,识别所述故障相上的故障区域;
步骤103’具体包括:根据故障相上的所述智能监测终端采集的故障工频电流计算相位;
判断相邻两个智能监测终端采集的故障工频电流相位是否相反,若相反,以该两个智能监测终端中的一个智能监测终端为第一终端,另一个智能监测终端为第二终端,判断位于第一终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第一终端采集的故障工频电流相位相同,且位于第二终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第二终端采集的故障工频电流相位相同,其中,第一终端第一侧为第一终端的不与第二终端相邻的一侧;第二终端第一侧为第二终端的不与第一终端相邻的一侧;若判断结果为是,则第一终端与第二终端之间的区域为故障区域。
5.根据权利要求4所述的方法,其特征在于,步骤107具体包括:
将第一终端及位于第一终端第一侧的所有智能监测终端为节点归集为集合G1,将第二终端及位于第二终端第一侧的所有智能监测终端为节点归集为集合G2
在集合G1与集合G2的节点中,将集合G1中所包括的a个节点以距离故障点近远顺序依次与集合G2中所包括的b个节点以距离故障点近远顺序两两进行双端定位,以获得定位结果;其中,所述故障点位于故障区域中;
根据表达式Lab=-△Tab×v+2daF对所述定位结果进行拟合,选择故障位置重合次数最多的结果作为故障定位结果;
其中:
Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
△Tab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;
daF为集合G1的第a个节点到故障点F的距离;
v为输电线路的行波波速度。
6.一种基于时域自适应窗的线路故障测距装置,输电线路上预先设置有监测点,每个监测点处的A相、B相和C相上分别安装有智能监测终端,其特征在于,包括:
第一采集模块,用于通过智能监测终端采集第一工频电流;
第一判断模块,用于对智能监测终端采集的所述第一工频电流进行分析,以判定所述智能监测终端所位于的线路是否发生故障,若是,触发第二采集模块;
第二采集模块,用于通过位于故障相上的智能监测终端采集故障电流行波,以获得故障电流行波录波文件,所述故障相为线路发生故障的某相,所述故障电流行波录波文件包括故障电流行波录波波形;
过零时间点获取模块,用于对各所述故障电流行波录波文件进行分析,以获得故障电流行波录波波形的过零时间点,其中,相邻所述过零时间点为一个时域自适应窗;
有效行波阵时间窗获取模块,用于计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所对应的时域自适应窗作为有效行波阵时间窗;
第一计算模块,用于在有效行波阵时间窗内,以η为固定时间窗,以采样时间T为滑窗时间间隔,滑窗次数为(△Tm-η)/T,根据表达式
Figure FDA0002381356340000031
分别计算有效行波阵时间窗内的各个滑窗的单位行波能量并以值最大的行波能量所对应的时间窗的最早时间点为行波首波波头时间点,其中,△Tm为有效行波阵时间窗,k1∈{0,1,2,…,N-1},滑窗次数N=(△Tm-η)/T,
Figure FDA0002381356340000032
表示第k1个过零时间点处的时间,i(t1)表示t1时刻的电流值,
第二计算模块,用于根据计算出的行波首波波头时间点,基于双端定位进行拟合测距。
7.根据权利要求6所述的装置,其特征在于,有效行波阵时间窗获取模块具体用于:
根据表达式
Figure FDA0002381356340000033
计算获得单位时间的电流行波能量,并以值最大的单位时间的电流行波能量所述对应的时域自适应窗作为有效行波阵时间窗;
其中,k-1∈{0,1,2,…,n-1},n为过零时间点总个数,i(t)表示t时刻的电流值,t∈(Tk-1,Tk),Tk-1表示第K-1个过零时间点处的时间,Tk表示与第K-1个过零时间点相邻的第K个过零时间点处的时间,Tk>Tk-1,Tk-Tk-1为一个时域自适应窗。
8.根据权利要求6所述的装置,其特征在于,第一判断模块具体用于:
对所述智能监测终端采集的第一工频电流进行分析,判断第一工频电流是否在预设时间段内增加预设增量后减少到零,若是,则所述智能监测终端所位于的线路为故障相,触发第二采集模块。
9.根据权利要求6-8中任一项所述的装置,其特征在于,还包括第二判断模块,所述第二判断模块具体用于:
根据故障相上的所述智能监测终端采集的故障工频电流计算相位;
判断相邻两个智能监测终端采集的故障工频电流相位是否相反,若相反,以该两个智能监测终端中的一个智能监测终端为第一终端,另一个智能监测终端为第二终端,判断位于第一终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第一终端采集的故障工频电流相位相同,且位于第二终端第一侧的所有智能监测终端采集的故障工频电流相位是否与第二终端采集的故障工频电流相位相同,其中,第一终端第一侧为第一终端的不与第二终端相邻的一侧;第二终端第一侧为第二终端的不与第一终端相邻的一侧;若判断结果为是,则第一终端与第二终端之间的区域为故障区域。
10.根据权利要求9所述的装置,其特征在于,所述第二计算模块具体用于:
将第一终端及位于第一终端第一侧的所有智能监测终端作为节点归集到集合G1中,将第二终端及位于第二终端第一侧的所有智能监测终端作为节点归集到集合G2中;
在集合G1与集合G2的节点中,将集合G1中所包括的a个节点以距离故障点近远顺序依次与集合G2中所包括的b个节点以距离故障点近远顺序两两进行双端定位,以获得定位结果;其中,所述故障点位于故障区域中;
根据表达式Lab=-△Tab×v+2daF对所述定位结果进行拟合,选择故障位置重合次数最多的结果作为故障定位结果;
其中:Lab为集合G1的第a个节点到集合G2的第b个节点的距离;
△Tab为集合G1的第a个节点到集合G2的第b个节点的行波首波波头时间点的时间差;daF为集合G1的第a个节点到故障点F的距离;v为输电线路的行波波速度。
CN201711175222.7A 2017-11-22 2017-11-22 基于时域自适应窗的线路故障测距方法及装置 Active CN108008246B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711175222.7A CN108008246B (zh) 2017-11-22 2017-11-22 基于时域自适应窗的线路故障测距方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711175222.7A CN108008246B (zh) 2017-11-22 2017-11-22 基于时域自适应窗的线路故障测距方法及装置

Publications (2)

Publication Number Publication Date
CN108008246A CN108008246A (zh) 2018-05-08
CN108008246B true CN108008246B (zh) 2020-04-21

Family

ID=62053114

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711175222.7A Active CN108008246B (zh) 2017-11-22 2017-11-22 基于时域自适应窗的线路故障测距方法及装置

Country Status (1)

Country Link
CN (1) CN108008246B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113702767B (zh) * 2021-08-31 2023-01-10 中南大学 基于小波滑窗能量的孤岛直流微电网故障诊断方法
CN113917833A (zh) * 2021-09-18 2022-01-11 广西电网有限责任公司柳州供电局 一种对时装置及对时方法
CN113848428B (zh) * 2021-09-29 2022-06-14 华南理工大学 一种输电线路双端故障测距方法、系统、装置及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116112A (zh) * 2013-01-06 2013-05-22 广东电网公司电力科学研究院 一种同杆并架双回线路故障测距方法
CN104133980A (zh) * 2013-05-02 2014-11-05 弗兰克公司 自适应频域加窗
CN105242179A (zh) * 2015-09-25 2016-01-13 山东山大电力技术有限公司 一种阻抗法和行波法相结合的行波综合测距方法
WO2016188574A1 (en) * 2015-05-28 2016-12-01 Abb Technology Ltd Travelling wave pilot protection of a transmission line based on time synchronization
CN106771861A (zh) * 2016-12-08 2017-05-31 中国矿业大学 基于广域行波能量和时间差的复杂电网故障定位方法
CN107037316A (zh) * 2016-11-23 2017-08-11 国网湖北省电力公司检修公司 一种适应于输电线路行波测距的单端波形自动识别方法
CN107219440A (zh) * 2017-04-11 2017-09-29 中国矿业大学 单端辐射型配电网单相接地故障的定位方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116112A (zh) * 2013-01-06 2013-05-22 广东电网公司电力科学研究院 一种同杆并架双回线路故障测距方法
CN104133980A (zh) * 2013-05-02 2014-11-05 弗兰克公司 自适应频域加窗
WO2016188574A1 (en) * 2015-05-28 2016-12-01 Abb Technology Ltd Travelling wave pilot protection of a transmission line based on time synchronization
CN105242179A (zh) * 2015-09-25 2016-01-13 山东山大电力技术有限公司 一种阻抗法和行波法相结合的行波综合测距方法
CN107037316A (zh) * 2016-11-23 2017-08-11 国网湖北省电力公司检修公司 一种适应于输电线路行波测距的单端波形自动识别方法
CN106771861A (zh) * 2016-12-08 2017-05-31 中国矿业大学 基于广域行波能量和时间差的复杂电网故障定位方法
CN107219440A (zh) * 2017-04-11 2017-09-29 中国矿业大学 单端辐射型配电网单相接地故障的定位方法
CN107219440B (zh) * 2017-04-11 2018-11-27 中国矿业大学 单端辐射型配电网单相接地故障的定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CSC-HVDC输电线路单端行波自动故障定位方法;宋国兵 等;《高电压技术》;20140228;第40卷(第2期);第588-596页 *
配电网故障自动定位技术研究综述;唐金锐 等;《电力自动化设备》;20130531;第33卷(第5期);第7-13页 *

Also Published As

Publication number Publication date
CN108008246A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
CN104597376B (zh) 一种考虑实测波速的高压直流输电线路故障测距方法
CN102819004B (zh) 智能变电站数字化电能计量系统性能综合检测与分析平台
CN108008246B (zh) 基于时域自适应窗的线路故障测距方法及装置
CN102122823B (zh) 电力系统振荡扰动源定位的方法
CN105044524A (zh) 一种适用于智能变电站的避雷器监测方法及系统
CN109116289B (zh) 一种高压互感器群运行状态实时采集系统
CN104142421B (zh) 变电设备绝缘在线监测系统及其工作方法
CN102323494A (zh) 一种多谐波源辨识方法
CN112946424B (zh) 一种故障精确定位的方法及装置
CN108896879B (zh) 基于局部放电信号特征的诊断图谱相位开窗参数调整方法
CN109524972B (zh) 基于gso和svm算法的低频振荡参数估计方法
CN106019076A (zh) 一种高压直流输电线路故障测距方法
CN110907755A (zh) 一种输电线路在线监测故障识别方法
CN106646034A (zh) 一种避雷器结构健康在线监测分析系统
CN106405230A (zh) 频率测量方法和装置
CN110514963B (zh) 一种改进双端行波故障测距方法
CN104865467A (zh) 应用于主动配电网的谐波源定位和谐波责任划分方法
CN110247420B (zh) 一种hvdc输电线路故障智能识别方法
CN117192292B (zh) 一种雷击接地极线路故障测距方法及系统
CN113945806B (zh) 一种柔性直流输电线路单端故障测距方法及系统
CN103941149A (zh) 电气化铁道牵引网故障定位方法
CN109829652B (zh) 一种长时间尺度动态谐波责任划分方法
CN111398761A (zh) 局部放电信号采集装置和局部放电采集分析仪
CN102401856B (zh) 一种基于模糊聚类方法的工频电压信号实时频率测量方法
CN104865508A (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