CN103558512B - 一种基于矩阵运算的配电网10kV馈线故障定位方法 - Google Patents
一种基于矩阵运算的配电网10kV馈线故障定位方法 Download PDFInfo
- Publication number
- CN103558512B CN103558512B CN201310579589.0A CN201310579589A CN103558512B CN 103558512 B CN103558512 B CN 103558512B CN 201310579589 A CN201310579589 A CN 201310579589A CN 103558512 B CN103558512 B CN 103558512B
- Authority
- CN
- China
- Prior art keywords
- matrix
- prime
- fault
- node
- feeder
- 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)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
本发明是一种基于矩阵运算的配电网10kV馈线故障定位方法,属于配电网自动化技术领域。所要解决的技术问题是配电网中压10kV馈线的故障定位问题,本发明方法是先根据配电网馈线拓扑结构形成配电网初始供电路径矩阵H,为解决故障重复定位对H进行修正得到修正路径矩阵Hˊ,通过馈线终端装置(FTU)获得故障信息,再根据故障信息对修正路径矩阵Hˊ进行运算处理,形成初始判断矩阵Pˊ,最终得到故障判断矩阵P,根据P中的相关元素可以有效地判断出故障区域的父节点,从而得到故障区域。基于本发明方法实现配电网馈线故障定位,判断原理简单、有效,无需进行繁琐的矩阵相乘运算,并且可以很好的解决馈线末端故障定位的问题。
Description
技术领域
本发明属于配电网自动化领域,涉及一种基于矩阵运算的配电网10kV馈线故障定位新方法。
背景技术
配电自动化系统的一个重要功能就是馈线自动化,而馈线自动化的主要功能在于准确定位故障,迅速隔离故障区域并恢复非故障区域供电。目前,现有的故障定位原理和方法主要可分为三类:一类是以在线路端点处测量故障距离为目的的故障测距法;一类是故障发生后通过向系统注入特定信号实现寻迹的信号注入法;还有一类是根据馈线各开关处安装的馈线终端装置(FeederTerminalUnit,FTU)所提供的故障信息(一般是过流信息)确定故障区段的故障区段定位法。第一类从理论上可行,但由于配电网结构复杂,混合线路接头处、线路分支处和负荷处均为波阻抗不连续点,行波在波阻抗不连续点的折射和反射造成线路一端测得的行波波形将特别复杂,很难正确识别出故障点的反射波,使测距实现困难;第二类具有测距结果不受负载参数变化影响的优点,但由于其取地模网络做为故障定位信息依据,不能解决只存在线模分量的相间短路故障的定位问题;第三类利用线路FTU实现线路故障信息的检测,可以自动定位出故障所在区段,在配电自动化系统中最为常用。
基于故障信息的故障区段定位方法有多种,其中应用较为广泛的是利用矩阵算法实现故障的分段定位,其基本思路是将网络描述矩阵和故障信息矩阵相乘,再进行规格化处理,根据最后所得的判别矩阵判断故障区域。但对于N个节点的配电网,网络描述矩阵是一个N×N方阵,需占用较多存储空间,有的算法甚至需要繁杂的矩阵相乘运算,运算量大,处理时间长;对于馈线末端故障,缺乏有效定位,或者定位过程比较繁琐。本发明定义了馈线的供电路径矩阵,并在该矩阵的基础上进行运算,与已有的矩阵算法相比,矩阵规模更小且矩阵运算更简单,并能有效地定位馈线末端故障,判断原理简单、有效,可满足实时性要求。
发明内容
本发明所要解决的技术问题是一种基于矩阵运算的配电网中压10kV馈线故障定位方法。该方法根据配电网馈线拓扑结构形成配电网初始供电路径矩阵H,为解决故障重复定位对H进行修正得到修正路径矩阵H',通过FTU获得故障信息,再根据故障信息对修正路径矩阵H'进行运算处理,形成初始判断矩阵P',最终得到故障判断矩阵P,根据P中的相关元素可以有效地判断出故障区域的父节点,从而得到故障区域。基于本发明方法实现配电网馈线故障定位,判断原理简单、有效,无需进行繁琐的矩阵相乘运算,并且可以很好的解决馈线末端故障定位的问题。
本发明的技术解决方案如下:
一种基于矩阵运算的配电网10kV馈线故障定位新方法,其特征在于,包括如下步骤:
步骤1:对于一个辐射状结构的包含一个电源、M个负荷及N个开关的10kV馈线,将电源、各负荷及各开关均当做节点,将相邻节点间的各支路当做边,则该馈线构成一个无向连通图;采用深度优先搜索法对连通图中的开关节点进行搜索,即先从电源节点出发搜索其一个邻接开关节点,从该节点搜索下一个子节点,如此不断地向纵深前进直到不能再前进(到达负荷节点)时,才从当前节点返回到上一级节点,沿另一方向又继续前进,这种方法保证了的搜索树是从树根开始一枝一枝逐渐形成的;当搜索完整个馈线后,按搜索到各开关节点的顺序对开关节点编号;
步骤2:由电源节点向一个负荷节点供电所经历的开关节点的有序集合构成一个供电路径,一条含M个负荷的馈线有M个供电路径,故障发生前,建立由M个供电路径和N个开关节点所构成的M行N列的配电网初始路径矩阵H:
其中,若路径i包含节点j,则矩阵H第i行j列元素hij=1,否则hij=0;
步骤3:若初始路径矩阵H中第j列存在元素hij=1,则令第j列中所有满足k>i的元素hkj=0,如此对各列进行处理后得到修正路径矩阵H':
H'每列只会有一个元素的值为1,使公共节点只在一条供电路径中出现,保证故障不重复定位;
步骤4:馈线某一处发生相间短路故障时,由各开关节点处的FTU检测各节点处电流值,该值远大于正常电流值时,则有故障电流流过,若检测到节点j有故障电流流过,则保持修正路径矩阵H'的第j列元素不变,否则第j列中所有元素置0,由此形成初始判断矩阵P':
若初始判断矩阵P'中的元素p′ij=1,说明节点j有故障电流流过;
步骤5:若矩阵P'中第j列存在元素p′ij=1,则将矩阵P'中满足k<j的所有第k列元素均置0,经处理后得到故障判断矩阵P,P矩阵中只会有一个元素值为1:
在故障判断矩阵P中,值为1的元素所对应的列编号即为故障节点编号,以该故障节点为父节点的最小供电区域为故障区段。
有益效果:
(1)基于矩阵运算的配电网馈线故障定位方法判断原理简单、有效,实时性好。
(2)对于一个辐射状结构的包含一个电源、M个负荷及N个开关的10kV馈线,传统的网络描述矩阵是一个(N+M+1)×(N+M+1)方阵,采用本发明方法只需构成M×N矩阵,同时无需复杂的矩阵相乘运算,减少了存储空间,运算量小,计算速度快。
(3)修正后的供电路径矩阵成为稀疏矩阵,解决了故障重复定位的问题(即多条供电路径中的公共节点都显示有故障存在),减轻了运算负担,当定位对象是复杂配电网时尤为明显。
(4)该方法不仅可以对一般故障定位,而且可以对馈线末端故障定位,无需假定虚拟节点和正方向,存储空间也有所减少。
(5)在配电网拓扑结构不发生变化的情况下,馈线供电路径矩阵H和修正路径矩阵H'不会发生变化,完全可以根据馈线拓扑结构事先设定相应的路径矩阵,在故障发生时只需直接进入故障判定的阶段,可节约运算时间,提高定位速度。
附图说明
图1基于矩阵运算的配电网馈线故障定位方法的实现流程框图;
图2一个典型的配电网馈线简化网络拓扑模型。
具体实施方式
以下将结合附图和具体实施例对本发明做进一步详细说明:
实施例1:
以下为利用矩阵运算对配电网馈线一般区域故障进行定位的实施步骤。
(1)采用深度优先搜索法对附图2中的开关节点进行编号,从电源点开始出发,先搜索到与之邻接的第一个开关子节点,编号为1,与节点1相邻接的有两个开关子节点,访问其中一个开关子节点并编号为2,继续深度优先搜索,访问到节点2的子节点并编号为3,直至访问到负荷节点为止,得到供电路径1由节点集合[1,2,3]构成;供电路径1搜索完后,退回到节点1,访问节点1的另一个子节点并编号为4,……,如此继续访问搜索,直至图中所有和电源节点有路径相通的节点均已被搜索到,得到供电路径2由节点[1,4,5,6,7,8]构成,供电路径3由节点[1,4,5,9,10,11]构成,整个馈线是一个3路径11节点的模型,各开关节点的标号如附图2所示;
(2)若路径i包含节点j,则矩阵H第i行j列元素hij=1,否则hij=0。例如,路径2包含节点4,因此矩阵中第2行4列元素h24=1;路径1不包含节点5,因此矩阵中第1行5列元素h15=0。由此得到一个3行11列的配电网馈线供电路径矩阵H:
(3)通过附图2和配电网供电路径矩阵H可知,节点1是路径1、2、3的公共节点,节点4、5是路径2和路径3的公共节点。当故障发生在以公共节点(如节点1)为父节点的最小配电区域中时,就会出现故障重复定位的问题,即多条路径都显示有故障存在,需要对供电路径矩阵H作出修正,对路径中的公共节点做处理,节点1在第1列存在元素h11=1,h21=1,h31=1,令h21=0,h31=0。其他节点以此类推,对配电网供电路径矩阵H进行处理,从而得到修正路径矩阵H':
(4)假设故障发生于节点4和5之间,各节点处的FTU均正常工作,可以有效检测到故障电流,则FTU检测到节点1有故障电流流过,修正路径矩阵H'中第1列中元素不变;检测到节点5没有故障电流流过,将修正路径矩阵H'中第5列中所有元素置0,各节点检测结束后得到初始判断矩阵P':
(5)为了得到要寻找的故障节点,需对初始判断矩阵P'进一步处理,在初始判断矩阵P'中,第4列中存在元素p′ij=1,则矩阵中第k(k<4)列所有元素均置0,得到故障判断矩阵P:
在故障判断矩阵P中,值为1的元素为p24,元素对应第4列,因此故障区域为以节点4为父节点的区段,反映到附图2中即为节点4和节点5之间的区段。
实施例2:
以下为利用矩阵运算对配电网馈线末端故障进行定位的实施步骤。
(1)假设附图2中配电网中的故障位于节点3之后的馈线末端,各节点处的FTU均正常工作。定位过程中形成配电网馈线供电路径矩阵H和修正路径矩阵H'的步骤与实施例1相同,在此不加赘述。
(2)FTU检测到节点1、2、3有故障电流流过,修正路径矩阵H'中第1、2、3列中元素不变;其他节点没有故障电流流过,修正路径H'中其他列中所有元素置0,得到初始判断矩阵P':
由初始判断矩阵P'判断故障位于路径1中。
(3)在初始判断矩阵P'中第1、2、3列中存在元素p′ij=1,需对矩阵进行处理:矩阵中第k(k<3)列所有元素均置0,由此得到故障判断矩阵P:
在故障判断矩阵P中,值为1的元素为p13,元素对应第3列,因此故障区域为以节点3为父节点的末端区段。
Claims (1)
1.一种基于矩阵运算的配电网10kV馈线故障定位新方法,其特征在于,包括如下步骤:
步骤1:对于一个辐射状结构的包含一个电源、M个负荷及N个开关的10kV馈线,将电源、各负荷及各开关均当做节点,将相邻节点间的各支路当做边,则该馈线构成一个无向连通图,采用深度优先搜索法对连通图中的开关节点进行搜索,并按各开关节点搜索到的顺序对开关节点进行编号;
步骤2:由电源节点向一个负荷节点供电所经历的开关节点的有序集合构成一个供电路径,一条含M个负荷的馈线有M个供电路径,故障发生前,建立由M个供电路径和N个开关节点所构成的M行N列的配电网初始路径矩阵H:
其中,若路径i包含节点j,则矩阵H第i行j列元素hij=1,否则hij=0;
步骤3:若初始路径矩阵H中第j列存在元素hij=1,则令第j列中所有满足k>i的元素hkj=0,如此对各列进行处理后得到修正路径矩阵H',H'每列只会有一个元素的值为1:
步骤4:馈线某一处发生相间短路故障时,由各开关节点处的馈线终端装置(FeederTerminalUnit,FTU)检测故障电流,若检测到节点j有故障电流流过,则保持修正路径矩阵H'的第j列元素不变,否则第j列中所有元素置0,由此形成初始判断矩阵P':
步骤5:若矩阵P'中第j列存在元素p′ij=1,则将矩阵P'中满足k<j的所有第k列元素均置0,经处理后得到故障判断矩阵P,P矩阵中只会有一个元素值为1:
在故障判断矩阵P中,值为1的元素所对应的列编号即为故障节点编号,以该故障节点为父节点的最小供电区域为故障区段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310579589.0A CN103558512B (zh) | 2013-11-19 | 2013-11-19 | 一种基于矩阵运算的配电网10kV馈线故障定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310579589.0A CN103558512B (zh) | 2013-11-19 | 2013-11-19 | 一种基于矩阵运算的配电网10kV馈线故障定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103558512A CN103558512A (zh) | 2014-02-05 |
CN103558512B true CN103558512B (zh) | 2016-02-10 |
Family
ID=50012814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310579589.0A Expired - Fee Related CN103558512B (zh) | 2013-11-19 | 2013-11-19 | 一种基于矩阵运算的配电网10kV馈线故障定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103558512B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104820163A (zh) * | 2015-04-13 | 2015-08-05 | 同济大学 | 一种基于ftu的配电网故障定位方法 |
CN107024638A (zh) * | 2016-02-02 | 2017-08-08 | 天津理工大学 | 一种基于改进矩阵算法的配电网故障定位方法及装置 |
CN105977932B (zh) * | 2016-06-14 | 2018-07-31 | 华北电力大学 | 一种配电网区域保护关联域在线计算系统及其方法 |
CN106124929A (zh) * | 2016-06-30 | 2016-11-16 | 湖南大学 | 一种配电网物理故障与信息故障辨识方法 |
CN106226657B (zh) * | 2016-09-07 | 2021-04-23 | 上海市政工程设计研究总院(集团)有限公司 | 一种用于电网故障诊断的停电区域识别方法 |
CN106443360B (zh) * | 2016-11-11 | 2019-11-12 | 国网新疆电力公司巴州供电公司 | 基于关联矩阵的配电网行波故障定位方法 |
CN106849023B (zh) * | 2017-03-09 | 2019-04-02 | 西安科技大学 | 基于行波信息的广域后备保护关联域识别方法 |
CN107271852B (zh) * | 2017-07-18 | 2019-08-20 | 广东双新电气科技有限公司 | 基于电压暂降信息的复杂配电网故障定位方法 |
CN108206522B (zh) * | 2018-01-10 | 2021-03-30 | 云南电网有限责任公司电力科学研究院 | 一种电力设备状态监控方法及系统 |
CN109254225A (zh) * | 2018-08-13 | 2019-01-22 | 中国南方电网有限责任公司超高压输电公司南宁监控中心 | 一种电网故障检测和故障线路识别方法 |
CN109342887B (zh) * | 2018-11-22 | 2021-11-09 | 广州供电局有限公司 | 基于相量测量装置的配电网故障统一定位方法及系统 |
CN110609205A (zh) * | 2019-08-15 | 2019-12-24 | 国网江苏省电力有限公司南京供电分公司 | 一种基于综合信息矩阵的配电网故障定位的方法 |
CN110596528A (zh) * | 2019-08-23 | 2019-12-20 | 国网山东省电力公司寿光市供电公司 | 一种基于统一矩阵算法的配电网故障定位方法及系统 |
CN110794263B (zh) * | 2019-12-17 | 2021-11-12 | 国网新疆电力有限公司昌吉供电公司 | 一种含分布式电源的配电网线路故障区段定位方法 |
CN111474442B (zh) * | 2020-03-30 | 2022-05-24 | 国网山东省电力公司德州供电公司 | 一种配电网单相接地故障定位方法及系统 |
CN111766475A (zh) * | 2020-07-14 | 2020-10-13 | 国网江苏省电力有限公司泗洪县供电分公司 | 一种断线故障定位矩阵算法 |
CN112698150B (zh) * | 2020-12-11 | 2022-05-24 | 河北工业大学 | 基于配电变压器监测终端的配电网行波故障定位方法 |
CN112766639B (zh) * | 2020-12-28 | 2023-03-28 | 上海电力大学 | 一种配网多故障抢修任务快速分配方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1112589C (zh) * | 1999-11-20 | 2003-06-25 | 深圳市中兴通讯股份有限公司 | 对配电网故障进行定位的方法 |
CN101188359B (zh) * | 2007-11-16 | 2011-05-25 | 湖南大学 | 基于元件的故障传递特性的配电网可靠性评估方法 |
CN101216530A (zh) * | 2007-12-29 | 2008-07-09 | 湖南大学 | 一种电子电路测试与故障诊断参数识别优化方法 |
CN101718834B (zh) * | 2009-12-01 | 2012-07-04 | 湖南大学 | 一种模拟电路故障传播特性的分析方法 |
JP5161930B2 (ja) * | 2010-07-12 | 2013-03-13 | 三菱電機株式会社 | 事故点標定システム |
CN102185315B (zh) * | 2011-05-16 | 2013-10-09 | 广东电网公司深圳供电局 | 基于EMS的110kV地区电网网络重构方法 |
-
2013
- 2013-11-19 CN CN201310579589.0A patent/CN103558512B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103558512A (zh) | 2014-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103558512B (zh) | 一种基于矩阵运算的配电网10kV馈线故障定位方法 | |
CN101968525B (zh) | 一种仿真计算与实时监测相结合的配电网故障定位方法 | |
CN103576048B (zh) | 一种用于电压暂降源定位的可能故障线路集提取方法 | |
CN104820163A (zh) | 一种基于ftu的配电网故障定位方法 | |
CN104133156B (zh) | 一种基于故障距离区间的混联线路单端行波故障测距方法 | |
CN104635114A (zh) | 一种电能质量扰动源定位系统及定位方法 | |
CN103454559A (zh) | 一种配电网单相接地故障区段定位方法及定位装置 | |
CN103645423A (zh) | 一种输电网络的故障定位方法和系统 | |
CN109474023B (zh) | 智能配电网区段实时更新方法、系统、存储介质及终端 | |
CN103267926A (zh) | 基于差分进化算法故障特征匹配的含dg配电网故障测距 | |
CN108693446B (zh) | 一种非同步采样电网暂态行波模量时间差的故障定位方法 | |
CN104599193A (zh) | 一种基于规则库的配电网单相接地故障定位方法 | |
CN107219501A (zh) | 一种基于gis局部放电的精确定位方法 | |
CN105004969A (zh) | 架空线配电网单相接地故障的定位方法 | |
CN102508097B (zh) | 线路故障分段检测装置及检测方法 | |
CN105203951A (zh) | 一级马尔科夫模型开关磁阻电机系统可靠性定量评估方法 | |
CN104750878A (zh) | 一种基于混合搜索策略的拓扑故障诊断方法 | |
CN105425105B (zh) | 电网故障线段综合定位方法及系统 | |
CN106600452A (zh) | 基于时间分析矩阵及聚类分析的配电网行波故障定位方法 | |
CN109738760A (zh) | 一种融合多种配电终端数据的配网短路故障定位方法 | |
CN107561405B (zh) | 一种基于非接触传感装置的故障选线定位方法及系统 | |
CN107229004A (zh) | 一种多分支输电线路故障的定位方法 | |
CN105929305A (zh) | 一种非全程混压双回线故障区段识别与精确测距方法 | |
CN109828185B (zh) | 针对含屋顶光伏发电系统的配电网故障定位方法 | |
Tashakkori et al. | Optimal placement of synchronized voltage traveling wave sensors in a radial distribution network |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160210 Termination date: 20181119 |