CN109188193B - 一种基于特征频带收敛性闵氏距离的配电网故障选线方法 - Google Patents
一种基于特征频带收敛性闵氏距离的配电网故障选线方法 Download PDFInfo
- Publication number
- CN109188193B CN109188193B CN201810834672.0A CN201810834672A CN109188193B CN 109188193 B CN109188193 B CN 109188193B CN 201810834672 A CN201810834672 A CN 201810834672A CN 109188193 B CN109188193 B CN 109188193B
- Authority
- CN
- China
- Prior art keywords
- fault
- frequency band
- feeder
- characteristic frequency
- distance
- 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
Images
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/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Locating Faults (AREA)
Abstract
本发明涉及一种基于特征频带收敛性闵氏距离的配电网故障选线方法,属于电力系统继电保护技术领域。当配电网发生单相接地故障时,采集各馈线的零序电流数据,在短时窗内对零序电流采用小波变换计算出各馈线在各尺度下的总能量,根据能量和最大原则确定特征频带,计算各馈线间特征频带的收敛性闵氏距离,形成特征频带的故障距离矩阵,利用特征频带的故障距离矩阵计算出各馈线之间的综合故障距离,根据馈线间的综合故障距离得到各馈线发生单相接地故障的相对概率,根据发生故障相对概率的大小进行故障选线。理论分析和大量仿真表明本发明效果良好。
Description
技术领域
本发明涉及一种基于特征频带收敛性闵氏距离的配电网故障选线方法,属于电力系统继电保护技术领域。
背景技术
配电网小电流接地系统发生单相接地故障时会使非故障相对地电压升高,电压升高对配电网的绝缘产生破坏,进而发展为相间或多点接地故障,引起系统过电压,因此必须准确、快速找到故障线路并及时使故障线路与系统隔离。为找到故障线路,传统的方法有以下几种:(1)人工拉路,人工逐个线路拉闸以确定故障线路,该方法会造成正常供电线路瞬时停电,降低供电的可靠性;(2)利用信号注入确定故障线路,需要增加额外的设备装置,增加成本;(3)利用行波确定线路,但对于缆线混合线路,因线路参数不同,在二者接点处,行波会发生衰减甚至被淹没,无法判断故障线路。且目前配电网单相接地故障选线的研究对象模型都是直配线路,并无分支,但配电网多含有分支线路。若在自动重合闸后,还无法排除故障,停电时间将会延长,产生的操作过电压和谐振过电压会对配电网络造成冲击影响,因此准确、可靠的选出故障馈线极为重要。
发明内容
本发明要解决的技术问题是针对配电网发生单相接地故障选线现有的技术问题,本发明提供一种原理简单、各种故障情况下均能正确选线的基于特征频带收敛性闵氏距离的配电网故障选线方法。
本发明采用的技术方案是:一种基于特征频带收敛性闵氏距离的配电网故障选线方法,包括如下步骤:当配电网发生单相接地故障时,采集各馈线的零序电流数据,在短时窗内对零序电流采用小波变换计算出各馈线在各个频带下的总能量,根据能量和最大原则确定特征频带,计算各馈线间特征频带的收敛性闵氏距离,形成特征频带的故障距离矩阵,利用特征频带的故障距离矩阵计算出各馈线之间的综合故障距离,根据馈线间的综合故障距离得到各馈线发生单相接地故障的相对概率,根据发生故障相对概率的大小进行故障选线。
具体步骤如下:
(1)母线零序电压瞬时值大于设定阈值时,判断为发生了单相接地故障;
(2)若判定系统发生单相接地故障,采集各馈线零序电流,采用小波对故障零序电流作多分辨率分析,按照式(7)计算分解后各频带信号对应的能量和;
信号xi(n)在尺度j下的能量和Eij为
式(7)中,Dj(k)为信号xi(n)在尺度j下的小波重构系数,k=1,2,...,n,n为采样点数;
(3)除去工频量所在的频带,能量集中的频带包含了暂态电容电流的特征,根据能量和最大原则来选择特征频带来获取暂态电流的特征,则所有信号在尺度j的总能量Ej为
式(8)中,h表示共有h个信号;
式中,xi(k)表示馈线i的零序电流特征频带,xj(k)表示馈线j的零序电流特征频带,i,j=1,2,…,m,m为馈线数,表示馈线i与馈线j之间在特征频带下p次幂的收敛性闵氏距离,,考虑到保留小数点后两位数字,ε取0.005,矩阵D中对角线元素均为0,非对角线元素dij表示馈线i与馈线j之间特征频带下的收敛性闵氏距离;
(5)根据特征频带的故障距离矩阵求取综合故障距离矩阵d,定义馈线i与其他馈线的综合故障距离的平均作为馈线i的综合故障距离di,即
d=[d1 d2 … dm]
(6)根据综合故障距离求取每条馈线相对于其他馈线发生故障的可能性,定义馈线i发生故障的概率是本馈线与其他馈线的距离比上所有馈线综合故障距离的和,形成馈线i发生故障的相对概率Pi,如式(10)所示;
P=[P1 P2 … Pm]
(7)选出发生单相接地故障相对概率最大的三个按照大小排序分别为Pa,Pb,Pc,当满足Pa>Pb+Pc时,相对概率最大的Pa对应的馈线为故障馈线,否则判断为母线故障。
本发明的有益效果是:
1、本方法选用的基于特征频带收敛性闵式距离算法具有距离算法的代表性与统一性,更能体现信号之间的相似程度;
2、本方法利用了故障后1/4周期内各馈线零序电流数据,此时消弧线圈还未对系统进行补偿,因此该方法避免了消弧线圈补偿作用的影响;
3、本方法计算的收敛性闵氏距离是通过不断增加式(1)中p的值,来确定各馈线之间最终的闵氏距离,该闵氏距离具有较强的收敛性;
4、本方法使用的最终判据是基于各馈线发生单相接地故障的相对概率,在各种故障情况下均拉大了健全线路与故障线路的差别,提高了选线的准确性。
附图说明
图1为本发明用于仿真的小电流接地系统结构图;其中,G为无限大电源,T为主变压器,Tz为Z字型变压器,L为消弧线圈,R为消弧线圈的阻尼电阻;
图2为本发明的故障选线流程图;
图3为例1各尺度下的总能量;
图4为例1零序电流特征频带波形;
图5为例2各尺度下的总能量;
图6为例2零序电流特征频带波形;
图7为例3各尺度下的总能量;
图8为例3零序电流特征频带波形;
图9为例4各尺度下的总能量;
图10为例4零序电流特征频带波形。
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步说明。
实施例1:一种基于特征频带收敛性闵氏距离的配电网故障选线方法,包括如下步骤:当配电网发生单相接地故障时,采集各馈线的零序电流数据,在短时窗内对零序电流采用小波变换计算出各馈线在各个频带下的总能量,根据能量和最大原则确定特征频带,计算各馈线间特征频带的收敛性闵氏距离,形成特征频带的故障距离矩阵,利用特征频带的故障距离矩阵计算出各馈线之间的综合故障距离,根据馈线间的综合故障距离得到各馈线发生单相接地故障的相对概率,根据发生故障相对概率的大小进行故障选线。
具体步骤如下:
(1)母线零序电压瞬时值大于设定阈值时,判断为发生了单相接地故障,如果瞬时值大于相电压的15%,则判断发生了单相接地故障;
(2)若判定系统发生单相接地故障,采集各馈线零序电流,采用小波对故障零序电流作多分辨率分析,按照式(7)计算分解后各频带信号对应的能量和;
信号xi(n)在尺度j下的能量和Eij为
式(7)中,Dj(k)为信号xi(n)在尺度j下的小波重构系数,k=1,2,...,n,n为采样点数;
(3)除去工频量所在的频带,能量集中的频带包含了暂态电容电流的特征,根据能量和最大原则来选择特征频带来获取暂态电流的特征,则所有信号在尺度j的总能量Ej为
式(8)中,h表示共有h个信号;
式中,xi(k)表示馈线i的零序电流特征频带,xj(k)表示馈线j的零序电流特征频带,i,j=1,2,…,m,m为馈线数,表示馈线i与馈线j之间在特征频带下p次幂的收敛性闵氏距离,考虑到保留小数点后两位数字,ε取0.005,矩阵D中对角线元素均为0,非对角线元素dij表示馈线i与馈线j之间特征频带下的收敛性闵氏距离;
(5)根据特征频带的故障距离矩阵求取综合故障距离矩阵d,定义馈线i与其他馈线的综合故障距离的平均作为馈线i的综合故障距离di,即
d=[d1 d2 … dm]
(6)根据综合故障距离求取每条馈线相对于其他馈线发生故障的可能性,定义馈线i发生故障的概率是本馈线与其他馈线的距离比上所有馈线综合故障距离的和,形成馈线i发生故障的相对概率Pi,如式(10)所示;
P=[P1 P2 … Pm]
(7)选出发生单相接地故障相对概率最大的三个按照大小排序分别为Pa,Pb,Pc,当满足Pa>Pb+Pc时,相对概率最大的Pa对应的馈线为故障馈线,否则判断为母线故障。
本发明的原理是:
配电网发生单相接地故障后,采集各馈线的零序电流数据,在短时窗内对零序电流采用小波变换计算出各馈线在各频带的总能量,根据能量和最大原则确定特征频带,计算各馈线间特征频带的收敛性闵氏距离,形成特征频带的故障距离矩阵,利用特征频带的故障距离矩阵计算出各馈线之间的综合故障距离,根据馈线间的综合故障距离得到各馈线发生故障的相对概率,通过选出发生单相接地故障相对概率最大的三个按照大小排序分别为Pa,Pb,Pc,当满足Pa>Pb+Pc时,相对概率最大的Pa对应的馈线为故障馈线,否则判断为母线故障。
举例说明:下面结合具体的例子,对本发明进行详细说明。
例1:如图1所示110kV/35kV配电网单相接地故障仿真模型,有4条架空馈线且都含分支线路。Z字型变压器中性点通过消弧线圈串联电阻接地。设置馈线L1的分支线路1距离主干线2公里处发生金属性单相接地故障,故障角为90°,采样率为5kHz。本例子的基于特征频带收敛性闵氏距离的配电网单相接地故障选线方法的具体步骤如下:
谐振接地系统发生单相接地故障后,对各馈线故障后5ms短时窗零序电流采用小波变换,如图3所示,尺度5能量和最大,根据能量和最大原则,确定特征频带为156.25Hz至312.5Hz,其零序电流的特征频带如图4所示。
计算各馈线间特征频带的收敛性闵氏距离,形成配电网特征频带故障距离矩阵D。
用特征频带故障距离矩阵求取每条馈线相对于其他馈线的综合故障距离,形成配电网综合故障距离矩阵d;
d=[58.54 25.67 25.44 29.69]
用综合故障距离矩阵求取每条馈线相对于其他馈线发生故障的可能性,即各馈线发生单相接地故障的相对概率,形成相对概率矩阵P;
P=[0.43 0.18 0.18 0.21]
选出发生单相接地相对概率最大的前三个按照大小排序分别为Pa=0.43,Pb=0.21,Pc=0.18,满足Pa>Pb+Pc,相对概率最大的Pa对应的馈线L1为故障馈线,由此判定故障线路为L1,选线结果正确。
例2:110kV/35kV谐振接地系统配电网单相接地故障仿真模型与例1相同,馈线L4的分支线路2距离主干线3km处发生单相接地故障,接地电阻200Ω,故障角60°,采样频率为5kHz。
按例1相同的方法,如图5所示,尺度6的能量和最大,确定特征频带为78.125Hz至156.25Hz,其零序电流的特征频带如图6所示,计算各馈线间的特征频带故障距离矩阵D、综合故障距离矩阵d及相对概率矩阵P;
d=[5.36 5.38 5.41 14.07]
P=[0.18 0.18 0.18 0.46]
选出发生单相接地相对概率最大的前三个按照大小排序分别为Pa=0.46,Pb=0.18,Pc=0.18,满足Pa>Pb+Pc,相对概率最大的Pa对应的馈线L4为故障馈线,由此判定故障线路为L4,选线结果正确。
例3:110kV/35kV谐振接地系统配电网单相接地故障仿真模型与例1相同,馈线L2距离母线3km处发生单相接地故障,接地电阻600Ω,故障角30°,采样频率为5kHz。
按例1相同的方法,如图7所示,尺度6的能量和最大,确定特征频带为78.125Hz至156.25Hz,其零序电流的特征频带如图8所示,计算各馈线间的特征频带故障距离矩阵D、综合故障距离矩阵d及相对概率矩阵P;
d=[2.26 6.10 2.25 2.78]
P=[0.18 0.47 0.17 0.18]
选出发生单相接地相对概率最大的前三个按照大小排序分别为Pa=0.47,Pb=0.18,Pc=0.18,满足Pa>Pb+Pc,相对概率最大的Pa对应的馈线L2为故障馈线,由此判定故障线路为L2,选线结果正确。
例4:110kV/35kV谐振接地系统配电网单相接地故障仿真模型与例1相同,母线发生单相接地故障,接地电阻150Ω,故障角60°,采样频率为5kHz。
按例1相同的方法,如图9所示,尺度6的能量和最大,确定特征频带为78.125Hz至156.25Hz,其零序电流的特征频带如图10所示,计算各馈线间的特征频带故障距离矩阵D、综合故障距离矩阵d及相对概率矩阵P;
d=[3.91 3.92 4.09 4.50]
P=[0.24 0.24 0.25 0.27]
选出发生单相接地相对概率最大的前三个按照大小排序分别为Pa=0.27,Pb=0.25,Pc=0.24,不满足Pa>Pb+Pc,由此判定母线故障,选线结果正确。
以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (1)
1.一种基于特征频带收敛性闵氏距离的配电网故障选线方法,其特征在于:包括如下步骤:当配电网发生单相接地故障时,采集各馈线的零序电流数据,在短时窗内对零序电流采用小波变换计算出各馈线在各个频带下的总能量,根据能量和最大原则确定特征频带,计算各馈线间特征频带的收敛性闵氏距离,形成特征频带的故障距离矩阵,利用特征频带的故障距离矩阵计算出各馈线之间的综合故障距离,根据馈线间的综合故障距离得到各馈线发生单相接地故障的相对概率,根据发生故障相对概率的大小进行故障选线;
具体步骤如下:
(1)母线零序电压瞬时值大于设定阈值时,判断为发生了单相接地故障;
(2)若判定系统发生单相接地故障,采集各馈线零序电流,采用小波对故障零序电流作多分辨率分析,按照式(7)计算分解后各频带信号对应的能量和;
信号xi(n)在尺度j下的能量和Eij为
式(7)中,Dj(k)为信号xi(n)在尺度j下的小波重构系数,k=1,2,...,n,n为采样点数;
(3)除去工频量所在的频带,能量集中的频带包含了暂态电容电流的特征,根据能量和最大原则来选择特征频带来获取暂态电流的特征,则所有信号在尺度j的总能量Ej为
式(8)中,h表示共有h个信号;
式中,xi(k)表示馈线i的零序电流特征频带,xj(k)表示馈线j的零序电流特征频带,i,j=1,2,…,m,m为馈线数,
(5)根据特征频带的故障距离矩阵求取综合故障距离矩阵d,定义馈线i与其他馈线的综合故障距离的平均作为馈线i的综合故障距离di,即
d=[d1 d2 … dm]
(6)根据综合故障距离求取每条馈线相对于其他馈线发生故障的可能性,定义馈线i发生故障的概率是本馈线与其他馈线的距离比上所有馈线综合故障距离的和,形成馈线i发生故障的相对概率Pi,如式(10)所示;
P=[P1 P2 … Pm]
(7)选出发生单相接地故障相对概率最大的三个按照大小排序分别为Pa,Pb,Pc,当满足Pa>Pb+Pc时,相对概率最大的Pa对应的馈线为故障馈线,否则判断为母线故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810834672.0A CN109188193B (zh) | 2018-07-26 | 2018-07-26 | 一种基于特征频带收敛性闵氏距离的配电网故障选线方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810834672.0A CN109188193B (zh) | 2018-07-26 | 2018-07-26 | 一种基于特征频带收敛性闵氏距离的配电网故障选线方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109188193A CN109188193A (zh) | 2019-01-11 |
CN109188193B true CN109188193B (zh) | 2020-09-25 |
Family
ID=64937613
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810834672.0A Active CN109188193B (zh) | 2018-07-26 | 2018-07-26 | 一种基于特征频带收敛性闵氏距离的配电网故障选线方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109188193B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109856506B (zh) * | 2019-01-17 | 2021-07-27 | 陕西科技大学 | 基于邻点差异法的单相接地故障区域定位方法 |
CN109709448A (zh) * | 2019-03-06 | 2019-05-03 | 南京工程学院 | 一种基于同步挤压小波变换的配电网单相高阻接地故障选线方法 |
CN111337791A (zh) * | 2020-03-25 | 2020-06-26 | 国网河南省电力公司电力科学研究院 | 一种基于梯度提升树算法的配电网单相接地故障选线方法 |
CN113640621B (zh) * | 2021-08-27 | 2023-09-15 | 国网四川省电力公司成都供电公司 | 一种小电流接地故障线路全自动检测方法及检测系统 |
CN114814593B (zh) * | 2022-04-29 | 2024-06-11 | 哈尔滨工业大学(威海) | 一种两步检测策略的电池组多故障快速诊断方法 |
CN117909909B (zh) * | 2024-03-19 | 2024-06-14 | 青岛鼎信通讯股份有限公司 | 一种配电网弧光接地故障识别方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000357168A (ja) * | 1999-06-14 | 2000-12-26 | Nippon Telegr & Teleph Corp <Ntt> | 多次元空間検索における最小包囲領域作成方法、最小包囲球符号化方法、多次元空間データ構造、多次元空間データ更新方法および多次元空間データ探索方法と前記方法を実施するプログラムを記録した記録媒体および多次元空間データ構造を記録した記録媒体 |
CN101162838A (zh) * | 2007-11-29 | 2008-04-16 | 昆明理工大学 | 一种利用小波包分解和相关分析的小电流接地系统故障选线方法 |
CN103257304A (zh) * | 2013-04-10 | 2013-08-21 | 昆明理工大学 | 一种利用零序电流特征频带内cwt系数rms值的ann故障选线方法 |
CN107064759A (zh) * | 2017-04-25 | 2017-08-18 | 国网安徽省电力公司电力科学研究院 | 一种特高压设备绝缘缺陷类型判断方法和系统 |
-
2018
- 2018-07-26 CN CN201810834672.0A patent/CN109188193B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000357168A (ja) * | 1999-06-14 | 2000-12-26 | Nippon Telegr & Teleph Corp <Ntt> | 多次元空間検索における最小包囲領域作成方法、最小包囲球符号化方法、多次元空間データ構造、多次元空間データ更新方法および多次元空間データ探索方法と前記方法を実施するプログラムを記録した記録媒体および多次元空間データ構造を記録した記録媒体 |
CN101162838A (zh) * | 2007-11-29 | 2008-04-16 | 昆明理工大学 | 一种利用小波包分解和相关分析的小电流接地系统故障选线方法 |
CN103257304A (zh) * | 2013-04-10 | 2013-08-21 | 昆明理工大学 | 一种利用零序电流特征频带内cwt系数rms值的ann故障选线方法 |
CN107064759A (zh) * | 2017-04-25 | 2017-08-18 | 国网安徽省电力公司电力科学研究院 | 一种特高压设备绝缘缺陷类型判断方法和系统 |
Non-Patent Citations (1)
Title |
---|
改进的基于闵氏距离的软子空间聚类算法;支晓斌等;《西安邮电大学学报》;20151130;第20卷(第6期);56-60 * |
Also Published As
Publication number | Publication date |
---|---|
CN109188193A (zh) | 2019-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109188193B (zh) | 一种基于特征频带收敛性闵氏距离的配电网故障选线方法 | |
Wu et al. | Ultra-high-speed directional protection of transmission lines using mathematical morphology | |
CN103424669B (zh) | 一种利用故障馈线零序电流矩阵主成分分析第一主成分的选线方法 | |
Nanayakkara et al. | Fault location in extra long HVDC transmission lines using continuous wavelet transform | |
CN109655713B (zh) | 一种单相接地故障定位方法和系统 | |
CN108562820B (zh) | 小电流接地系统单相接地故障选线方法及装置 | |
CN109444657B (zh) | 一种配电网高阻接地故障区段定位方法 | |
CN103163417B (zh) | 一种基于短时窗高低频暂态能量比值的虚幻接地识别方法 | |
CN103474981B (zh) | 一种基于零序电流多阶差分变换方向的配网单相接地保护方法 | |
CN103197202A (zh) | 一种基于三相突变电流分量特征频带内小波系数相关分析的配网故障选线方法 | |
CN112557950B (zh) | 基于矩阵相似性的配电网谐振接地系统故障选线方法 | |
CN112069691A (zh) | 一种基于避雷器动作电流频带分析的配电网直击雷与感应雷辨识方法 | |
CN103823158B (zh) | 采用不变矩的谐振接地系统故障选线方法 | |
CN112255493A (zh) | 一种多判据综合投票配电网故障选线方法 | |
CN110736900A (zh) | 一种基于方向行波的配电线路单相接地故障选线方法 | |
CN109596935A (zh) | 一种结合磁场能衰减特征的输电线路故障行波法 | |
CN108957234A (zh) | 基于多量测点和故障纵向特征的配电网接地选线方法及系统 | |
CN112083278A (zh) | 一种基于站端故障电流宽频检测的配电网直击雷与感应雷辨识方法 | |
Grebovic et al. | Lightning outage performance of power distribution line located in Mountain Lovćen area | |
CN112162173A (zh) | 一种基于故障电流频带分布差异的配电网雷击与非雷击故障辨识方法 | |
Kezunovic et al. | Improved transmission line fault location using automated correlation of big data from lightning strikes and fault-induced traveling waves | |
Borghetti et al. | A statistical approach for estimating the correlation between lightning and faults in power distribution systems | |
CN115575771A (zh) | 一种适用于多回输电线路雷电绕击与反击的识别方法 | |
Dong et al. | Optimizing solution of fault location using single terminal quantities | |
CN112162170B (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 |