CN108572303A - 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法 - Google Patents
一种基于随机矩阵和Hausdorff距离的配电网故障选线方法 Download PDFInfo
- Publication number
- CN108572303A CN108572303A CN201810376560.5A CN201810376560A CN108572303A CN 108572303 A CN108572303 A CN 108572303A CN 201810376560 A CN201810376560 A CN 201810376560A CN 108572303 A CN108572303 A CN 108572303A
- Authority
- CN
- China
- Prior art keywords
- matrix
- hausdorff
- feeder
- feeder line
- line
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 102
- 238000010187 selection method Methods 0.000 title claims abstract description 21
- 229940050561 matrix product Drugs 0.000 claims abstract description 15
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 238000013519 translation Methods 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 abstract description 8
- 230000005611 electricity Effects 0.000 abstract 1
- 238000010606 normalization Methods 0.000 abstract 1
- 230000007935 neutral effect Effects 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000001629 suppression Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification 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/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
-
- 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
-
- 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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
-
- 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
本发明公开一种基于随机矩阵和Hausdorff距离的配电网故障选线方法,选取故障前后馈线的三相电流采样值,通过分块、平移处理,加入高斯白噪声,生成状态数据矩阵,采用随机矩阵奇异值等价变换得到乘积矩阵,归一化得到标准矩阵积,获得特征值向量,对它们进行概率统计,将概率P<10%的特征值当作离群点过滤掉,采用Hausdorff距离算法,计算某馈线与其它馈线的特征值向量之间的Hausdorff距离,去掉最大值,加以平均得到该馈线的平均Hausdorff距离,如果此平均距离大于阈值,则判断该馈线故障;如果各馈线的平均Hausdorff距离都小于阈值,则判断所连接的母线故障。本发明可准确判断故障馈线和母线,不依赖于配电网模型,不受故障位置、过渡电阻、初相角、线路类型的影响,有较好的实用性。
Description
技术领域
本发明属于配电网故障选线的技术领域,具体涉及一种基于随机矩阵和Hausdorff距离的配电网故障选线方法。
背景技术
配电网线路较短,主要采用中线点不直接接地方式,单相接地故障是配电网最常见的故障,占小电流接地***故障的80%以上。国内外学者深入研究配电网故障选线方法,现有技术中配电网故障选线存在以下方法:
利用馈线暂态零序电流的衰减直流分量、工频分量和高频分量,结合不同频段的线路阻抗特性,提出基于高低频段暂态信号的相关分析选线法;将暂态零序电流进行小波包分解,小波分解后的暂态零序电流在非故障馈线中变化趋势相同,而与故障馈线的变化趋势相反;若故障发生在母线上,则各馈线波形特征变化趋势相同。通过暂态零序能量在故障后各馈线上的大小和方向差异判断故障发生在哪一条馈线或母线上。利用D-S证据理论,提出了一种利用多种故障信息、连续多周期融合的配电网故障选线方法,充分利用多方面、不同时段的故障信息,将多判据多周期信息融合选线的问题转化为证据推理问题,用D-S证据理论组合证据后,得到各条线路的综合基本信度分配值,基本信度最大的为故障线路。利用“S注人法”向故障后的配电网人为输入某一信号,然后跟踪检测该信号的流通通路,由于只有故障线路的故障相在此通路中,所以该信号流通的线路为故障线路,从而实现故障选线的功能;根据故障后位移电压大小的不同,从消弧线圈电压互感器或故障相电压互感器注入可变频率恒流信号,利用零序信号电压计算接地电阻和阻尼率大小,进行故障选线。
以上方法大多依赖于模型的建立,都是从模型角度出发,但是对故障数据本身的分析不足。目前缺少将随机矩阵理论应用于配电网故障选线的方法。随机矩阵理论是一种数据处理方法,起源于量子物理的发展和研究。
基于随机矩阵理论的数据处理方法关注于数据本身的特点而非模型,可非常自然地处理海量数据,分析过程方法依赖于数据变换,对模型的依赖性低,更加适用于大规模复杂网络。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于随机矩阵和Hausdorff距离的配电网故障选线方法,以解决现有配电网故障选线依赖于***模型且少有将随机矩阵理论应用于配电网故障选线问题。
为了达到上述目的,本发明采取的技术方案是:
提供一种基于随机矩阵和Hausdorff距离的配电网故障选线方法,其
包括:
获取配电网馈线三相电流的采样数据,形成原始数据矩阵X0,通过分块、平移处理,形成矩阵Xc,将矩阵Xc和高斯白噪声矩阵叠加,形成状态数据矩阵X;
对状态数据矩阵X进行奇异值等价变换,得到奇异值等价矩阵Xu,定义若干个奇异值等价矩阵Xu的乘积矩阵为对它进行归一化处理,得到标准矩阵积Z;
计算标准矩阵积Z的特征值向量;
统计各特征值的分布,将概率P<10%的特征值当作离群点,将概率P<10%的特征值过滤掉,形成过滤后的特征值向量;
针对某馈线,计算该馈线与其它m-1条馈线的特征值向量之间的Hausdorff距离Hi-j,去掉其中的最大距离,计算剩余Hausdorff距离的平均值Hi.avg:
其中,Hi-j max为Hi-j的最大值,i、j代表不同馈线的特征值向量,m为馈线条数;
根据各条馈线正常运行时的历史采样数据,计算得到正常运行的各条馈线之间的Hausdorff距离Hi.his,计算得到故障选线阈值Hset:
其中,kH为可靠系数,其取值为1.3~1.5,Hi.his为正常运行时各条馈线之间的Hausdorff距离;
构造馈线故障选线判据:
在m条馈线中,如果馈线i的平均Hausdorff距离Hi.avg>Hset,则馈线i为故障馈线;
构造母线故障的判据:
如果与某母线连接的m条馈线中平均Hausdorff距离均满足Hi.avg<Hset
(i=1,2,…,m),则判断该母线故障。
优选地,在一个时间窗中,对于配电网中N个采样时刻的A、B、C三相电流采样数据,形成原始数据矩阵X0:
其中,矩阵元素的下标第1个数字的1、2、3分别代表A、B、C相电流;
由于原始数据矩阵X0的行数相对于列数较小,对其行列元素通过分块、平移调整,来获取较优的行列比;
保持矩阵X0行数不变,将X0的N列元素按顺序拆分成l块,则X0被分割成l个子矩阵,即:
其中,分别为3行、N/l列的子矩阵;
将子矩阵依次平移到的下方,扩展形成矩阵Xc:
矩阵Xc为3×l行、N/l列的矩阵,其行列比c=3×l2/N;
将矩阵Xc和高斯白噪声矩阵mag·W叠加,形成状态数据矩阵X:
X=Xc+mag·W
式中,mag为噪声幅值,其大小根据故障发生后非故障馈线采样数据波动程度确定。矩阵W的元素服从标准正态分布,维数与矩阵Xc相同。
优选地,奇异值等价矩阵Xu为:
其中,Xu为n×n的方阵,n=3×l,U是Haar酉矩阵,H表示共轭转置;
乘积矩阵为为:
其中,L为奇异值等价矩阵Xu的最大数量;
标准矩阵积Z为:
其中,σ表示标准差;
计算该标准矩阵积Z的特征值向量;
每条馈线在复平面中形成1个n维的特征值向量Λ=(λ1,λ2,…,λn),λ1,λ2,…,λn为标准矩阵积Z的n个复数特征值;
对于有m条馈线的配电网,分别计算m条馈线的状态数据矩阵X的标准矩阵积Z的特征值向量,在复平面中形成m个n维的特征值向量Λ1,Λ2,…,Λm。
优选地,离群点过滤的方法为:
将单环定理中的外圆沿半径方向均分为T个圆环,统计各圆环中特征值的概率Pi(i=1,2,…,T);
定义概率Pi<10%的圆环中的特征值作为离群点;
将所述离群点从特征值中过滤掉,保留概率Pi≥10%的圆环中的特征值。
优选地,计算各馈线的特征值向量之间的Hausdorff距离:
假设有A、B两组向量(向量元素个数可不等):
A={a1,…,ap}
B={b1,…,bq}
向量ai(i=1,2,…,p)与向量B中之间的最小距离定义为:
其中,d(ai,bi)=||ai-bi||,||·||表示2-范数;
对于向量A中所有的元素,满足上式的最大值为从集合A到集合B的Hausdorff单向距离,即:
同理,向量到向量A的Hausdorff单向距离为:
取从向量A到向量B、从向量B到向量A两者之间单向Hausdorff距离较大值,作为向量A与向量B之间Hausdorff距离:
H(B,A)=max[h(A,B),h(B,A)]
优选地,计算m条馈线的特征值向量的平均Hausdorff距离Hi.avg:
首先求取某条馈线i与其它m-1条馈线的特征值向量之间的m-1个Hausdorff距离;由于与故障馈线相关的Hausdorff距离较大,为了消除与故障馈线之间的Hausdorff距离对非故障馈线的平均Hausdorff距离的影响,先去掉馈线i的m-1个平均Hausdorff距离中的最大值然后对其余m-2个Hausdorff距离求平均值,得到馈线i与其它馈线的特征值向量之间的平均Hausdorff距离Hi.avg:
优选地,利用各馈线的历史数据,计算在各条馈线正常运行时它们之间的Hausdorff距离Hi.his中的最大值,并考虑一定的裕度,获得故障选线阈值如下:
由于电缆线路和架空线路的电抗、电纳和对地电容参数差别较大,可靠系数kH取1.3~1.5。
优选地,给出配电网故障选线判据:
配电网馈线故障选线判据:
如果m条馈线中存在馈线i的平均Hausdorff距离Hi.avg>Hset,则判断馈线i为故障馈线;
当母线故障时,与其连接的m条馈线的标准矩阵积的特征值分布接近,它们相互之间的Hausdorff距离较小,所有Hi.avg都小于Hset,由此产生针对母线故障的故障判据;
配电网母线故障判据:
如果与某母线连接的m个馈线中平均Hausdorff距离都满足Hi.avg<Hset(i=1,2,…,m),则判断该母线故障。
本发明提供的基于随机矩阵和Hausdorff距离的配电网故障选线方法,具有以下有益的效果:
1)本发明利用随机矩阵分析方法,只利用故障数据,不需要针对各种运行方式分别建立模型,不受中性点接地运行方式的影响,适用于中性点不接地和中性点经消弧线圈接地运行方式。
2)本方法的故障选线结果准确,不受故障位置、过渡电阻、故障初相角等影响,具有很高的可靠性。
3)本方法可应用于架空线和电缆线混合线路故障情况下的故障选线。
附图说明
图1为配电网故障选线方法流程。
图2为单环定理示意图。
图3为故障馈线和非故障馈线数据特征值分布示意图。
图4为5条馈线的配电网结构示意图。
图5为5条馈线特征值的分布示意图。
图6为5条馈线两两之间的Hausdorff距离针状图。
图7为不同故障情形下各馈线的平均Hausdorff距离。
图8为不同过渡电阻下平均Hausdorff距离(θ=45°、Lf=50%)。
图9为不同故障初相角下平均Hausdorff距离(Rg=10Ω、Lf=95%)。
图10为不同故障位置下平均Hausdorff距离(Rg=10Ω、θ=45°)。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
根据本申请的一个实施例,参考图1,提供一种基于随机矩阵和Hausdorff距离的配电网故障选线方法,其包括:
获取配电网馈线三相电流的采样数据,形成状态数据矩阵X,其具体步骤如下:
在一个时间窗中,对于配电网中N个采样时刻A、B、C三相电流的采样数据,可形成原始数据矩阵X0:
其中矩阵元素的下标第1个数字的1、2、3分别代表A、B、C相电流。
由于原始数据矩阵X0的行数相对于列数较小,对其行列元素通过分块和平移调整,来获取较优的行列比。
保持矩阵X0行数不变,将X0的N列元素按顺序拆分成l块,则X0被分割成l个子矩阵,即:
其中,分别为3行、N/l列的子矩阵。
将子矩阵依次平移到的下方,扩展形成矩阵Xc:
矩阵Xc为3×l行、N/l列的矩阵。其行列比c=3×l2/N。
将矩阵Xc和高斯白噪声矩阵叠加,形成状态数据矩阵X:
X=Xc+mag·W
式中,mag为噪声幅值,其大小根据故障发生后非故障馈线采样数据波动程度确定。矩阵W的元素服从标准正态分布,维数与纯净信号的状态数据矩阵Xc相同。
对状态数据矩阵X进行奇异值等价变换,得到状态数据矩阵X的奇异值等价矩阵Xu,定义若干个所述奇异值等价矩阵Xu的乘积矩阵为对它归一化,得到标准矩阵积Z,其详细步骤如下:
变换状态数据矩阵X,得到其奇异值等价矩阵Xu:
其中,Xu为n×n的方阵,n=3×l,U是Haar酉矩阵,H代表共轭转置,且有
L个状态数据矩阵Xi(i=1,2,…L),经过上述变换,每个都可得到对应的奇异值等价矩阵Xu,i(i=1,2,…L),本发明中L=1。
定义若干个奇异值等价矩阵Xu的乘积矩阵为
对乘积矩阵进行归一化处理,得到标准矩阵积Z,
其中,σ表示标准差。
计算得到标准矩阵积Z的特征值向量和特征值的分布;
计算该标准矩阵积Z的特征值向量。每条馈线在复平面中形成1个n维的特征值向量Λ=(λ1,λ2,…,λn),λ1,λ2,…,λn为标准矩阵积Z的n个复数特征值。
对于有m条馈线的配电网,分别计算m条馈线的状态数据矩阵X的标准矩阵积Z的特征值向量,在复平面中形成m个n维的特征值向量Λ1,Λ2,…,Λm
参考图3,当配电网正常运行时,状态数据矩阵X计算得到的标准矩阵积Z,标准矩阵积Z里的元素服从独立同分布的高斯随机变量,且均值为0,方差为1,其特征值大致分布在一个外径r1=1、内径r2=(1-c)L/2的圆环内。
当配电网故障时,状态数据矩阵X计算得到的标准矩阵积Z,标准矩阵积Z元素的随机性被破坏,特征值在复平面上的分布会“收缩”到半径为r2=(1-c)L/2的内圆中。
根据特征值的分布,采用概率统计的离群点过滤处理方法,对离群点特征值进行过滤,其具体步骤如下:
参考图2,由于离群点对Hausdorff距离的计算结果有影响,故需过滤去除离群点。利用概率统计的离群点过滤方法对离群点进行过滤。将单环定理中外圆沿半径方向到圆心均分为T个圆环,每个圆环中特征值出现的概率为Pi(i=1,2,…,T)。例如,随机生成一个服从高斯随机分布的120×300非厄米特矩阵,通过计算,可得到一个维数为120的特征值向量。令T=10,则该矩阵特征值在每个圆环中出现的概率如表1所示。
表1特征值的分布范围
从表1可看出,大部分特征值分布在(0.7,1]之间,占到92.5%,概率P<10%的特征值出现的概率较小,将P<10%的特征值作为离群点,将其从特征值向量中过滤掉,只保留特征值出现概率P≥10%的圆环中的特征值,供后续的Hausdorff距离求解使用。
利用Hausdorff算法,计算各馈线的特征值向量之间的Hausdorff距离。
Hausdorff距离计算的方法为:
假设有A、B两组向量(向量元素个数可不等):
A={a1,…,ap}
B={b1,…,bq}
向量ai(i=1,2,…,p)与向量B中之间的最小距离定义为:
其中,d(ai,bi)=||ai-bi||,||·||表示2-范数。
对于向量A中所有的元素,满足上式的最大值为从集合A到集合B的Hausdorff单向距离,即:
同理,向量B到向量A的Hausdorff单向距离为:
取从向量A到向量B、从向量B到向量A两者之间单向Hausdorff距离较大值,作为向量A与向量B之间Hausdorff距离:
H(B,A)=max[h(A,B),h(B,A)]
首先利用上述Hausdorff距离的计算方法,求取各馈线i与其它m-1条馈线的特征值向量之间的m-1个Hausdorff距离;
然后计算m条馈线的特征值向量的平均Hausdorff距离Hi.avg,方法如下:由于与故障馈线相关的Hausdorff距离较大,为了消除与故障馈线之间的Hausdorff距离对非故障馈线的平均Hausdorff距离的影响,先去掉馈线i的m-1个平均Hausdorff距离中的最大值再对其余m-2个Hausdorff距离求平均值,得到馈线i与其它馈线的特征值向量之间的平均Hausdorff距离Hi.avg:
利用各馈线的历史数据,计算在各条馈线正常运行时它们之间的Hausdorff距离Hi.his中的最大值,并考虑一定的裕度,获得故障选线阈值如下:
由于电缆线路和架空线路的电抗、电纳和对地电容参数差别较大,可靠系数kH取1.3~1.5。
给出配电网故障选线判据:
配电网馈线故障选线判据:
如果m条馈线中存在馈线i的平均Hausdorff距离Hi.avg>Hset,则判断馈线i为故障馈线。
当母线故障时,与其连接的m条馈线的标准矩阵积的特征值分布接近,它们相互之间的Hausdorff距离较小,所有Hi.avg都小于Hset,由此产生针对母线故障的故障判据。
配电网母线故障判据:
如果与某母线连接的m个馈线中平均Hausdorff距离都满足Hi.avg<Hset(i=1,2,…,m),则判断该母线故障。
实施例
根据本申请的一个实施例,参考图4,利用MATLAB/Simulink进行配电网的搭建和仿真。
在MATLAB/Simulink中搭建5条配电网,配电变压器低压侧有5条馈线,馈线l1、l2为架空线路模型,l3、l4为电缆线路模型,l5为缆-线混合模型(75%架空线+25%电缆)。其中,架空线正序参数R1=0.096Ω/km,L1=1.22mH/km,C1=0.011uF/km;零序参数R0=0.23Ω/km,L0=3.66mH/km,C0=0.007uF/km;电缆线正序参数R11=0.11Ω/km,L11=0.52mH/km,C11=0.29uF/km;零序参数R00=0.34Ω/km,L00=1.54mH/km,C00=0.19uF/km。
在消弧线圈接地***仿真时,采用过补偿度方式运行,过补偿度为8%,消弧线圈的电感为0.4379H、电阻为4.127Ω。状态变量选用三相电流幅值,算例中判断阈值Hset=0.24。
参考图5,馈线l1发生A相接地故障。由5条馈线三相电流幅值(对于互感器不完全星形接线方式可选取两相电流幅值)采样数据分别构成各自的状态数据矩阵,计算得到各自的特征值分布。
由图5可知,馈线l1的特征值分布在内圆中,偏离正常分布,与非故障馈线l2~l4特征值之间的距离较远,而馈线l2~l4的特征值基本上分布在圆环中,它们之间的距离较近。
参考图6,与馈线l1相关的Hausdorff距离数值都较大,其它馈线的距离数值都较小,通过故障选线判据,可得到馈线l1发生故障。
根据本申请的一个实施例,验证本发明对各种故障情形的普适性,在变压器中性点不接地***中,在馈线l150%处分别设置A相接地故障、AB相间故障、AB两相接地故障、ABC三相故障。平均Hausdorff距离计算结果如图7所示。由图可知,馈线l1的平均Hausdorff距离在各种故障情况下均大于Hset,判断馈线l1为故障馈线。
由于中性点不接地***中单相接地故障特征不明显,针对不同过渡电阻、故障角、故障位置做了大量的试验,如图8-图10是不同情况下各馈线的平均Hausdorff距离,其中Rg为接地过渡电阻,θ为故障角,Lf为故障点距母线的距离百分数。
表2为选取8种典型故障情形下的故障选线结果,包括架空线故障、电缆线故障、混合线中架空线部分、电缆线故障及母线故障。
表2中性点不接地***单相接地故障选线结果
试验结果表明本发明方法在不同线路类型、过渡电阻、故障角和故障距离条件下都能准确选线,在图8-图10中,非故障馈线的平均Hausdorff距离在不同故障条件下差别较小,而故障馈线的平均Hausdorff距离在高阻、零故障角、线路始端故障条件下最小,但仍然大于Hset。在极端故障条件下也能正确选线,对各种故障状态及不同类型线路,本发明方法具有较好的适应性。
在谐振接地***中,做单相接地故障实验。分别在馈线l1~l5设置高阻接地故障(Rg=1000Ω、θ=0°、Lf=5%),选线结果如表3所示。由表3可知,在谐振接地***中,极端故障条件下,本发明方法仍能准确选线,同时也验证了故障选线阈值(Hset=0.24)计算方法的合理性。
表3谐振接地***单相接地下的故障选线结果
本发明不依赖于***模型,充分分析了故障数据,基于随机矩阵和Hausdorff距离的配电网故障选线方法,利用配电网采集到的一个时间窗的电流幅值数据构建一个适用于故障选线的高维随机矩阵,分析该随机矩阵的特征值分布,基于Hausdorff距离获得各馈线的随机矩阵特征值之间的距离,计算故障选线阀值,构造选线判据,识别故障馈线和母线。本发明不依赖于***模型,实用性强;且不受故障位置、过渡电阻、故障初相角等影响,具有很高的可靠性;同时还适用于架空线路和电缆线路故障选线,具有较强的适用性。
虽然结合附图对发明的具体实施方式进行了详细地描述,但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内,本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。
Claims (8)
1.一种基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,包括:
获取配电网馈线三相电流的采样数据,形成原始数据矩阵X0,通过分块、平移处理,形成矩阵Xc,将矩阵Xc和高斯白噪声矩阵叠加,形成状态数据矩阵X;
对状态数据矩阵X进行奇异值等价变换,得到奇异值等价矩阵Xu,定义若干个所述奇异值等价矩阵Xu的乘积矩阵为对它进行归一化处理,得到标准矩阵积Z;
计算标准矩阵积Z的特征值向量;
统计各特征值的分布,将概率P<10%的特征值当作离群点,将概率P<10%的特征值过滤掉,形成过滤后的特征值向量;
针对某馈线,计算该馈线与其它m-1条馈线的特征值向量之间的Hausdorff距离Hi-j,去掉其中的最大距离,计算剩余Hausdorff距离的平均值Hi.avg:
其中,Hi-j max为Hi-j的最大值,i、j代表不同馈线的特征值向量,m为馈线条数;
根据各条馈线正常运行时的历史采样数据,计算得到正常运行的各条馈线之间的Hausdorff距离Hi.his,并计算得到故障选线阈值Hset:
其中,kH为可靠系数,其取值为1.3~1.5,Hi.his为正常运行时各条馈线之间的Hausdorff距离;
构造馈线故障选线判据:
在m条馈线中,如果馈线i的平均Hausdorff距离Hi.avg>Hset,则馈线i为故障馈线;
构造母线故障的判据:
若与某母线连接的m条馈线中平均Hausdorff距离均满足Hi.avg<Hset(i=1,2,…,m),则判断该母线故障。
2.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,在一个时间窗中,对于配电网中N个采样时刻的A、B、C三相电流采样数据,形成原始数据矩阵X0:
其中,矩阵元素的下标第1个数字的1、2、3分别代表A、B、C相电流;
由于原始数据矩阵X0的行数相对于列数较小,对其行列元素通过分块、平移调整,来获取较优的行列比;
保持矩阵X0行数不变,将X0的N列元素按顺序拆分成l块,则X0被分割成l个子矩阵,即:
其中,分别为3行、N/l列的子矩阵;
将子矩阵依次平移到的下方,扩展形成矩阵Xc:
矩阵Xc为3×l行、N/l列的矩阵,其行列比c=3×l2/N;
将矩阵Xc和高斯白噪声矩阵mag·W叠加,形成状态数据矩阵X:
X=Xc+mag·W
式中,mag为噪声幅值,其大小根据故障发生后非故障馈线采样数据波动程度确定。矩阵W的元素服从标准正态分布,维数与矩阵Xc相同。
3.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,所述奇异值等价矩阵Xu为:
其中,Xu为n×n的方阵,n=3×l,U是Haar酉矩阵,H表示共轭转置;
所述乘积矩阵为为:
其中,L为奇异值等价矩阵Xu的最大数量;
所述标准矩阵积Z为:
其中,σ表示标准差;
计算该标准矩阵积Z的特征值向量;
每条馈线在复平面中形成1个n维的特征值向量Λ=(λ1,λ2,…,λn),λ1,λ2,…,λn为标准矩阵积Z的n个复数特征值;
对于有m条馈线的配电网,分别计算m条馈线的状态数据矩阵X的标准矩阵积Z的特征值向量,在复平面中形成m个n维的特征值向量Λ1,Λ2,…,Λm。
4.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,所述离群点过滤的方法为:
将单环定理中的外圆沿半径方向均分为T个圆环,统计各圆环中特征值的概率Pi(i=1,2,…,T);
定义概率Pi<10%的圆环中的特征值作为离群点;
将所述离群点从特征值中过滤掉,保留概率Pi≥10%的圆环中的特征值。
5.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,计算各馈线的特征值向量之间的Hausdorff距离:
假设有A、B两组向量(向量元素个数可不等):
A={a1,…,ap}
B={b1,…,bq}
向量ai(i=1,2,…,p)与向量B中之间的最小距离定义为:
其中,d(ai,bi)=||ai-bi||,||·||表示2-范数;
对于向量A中所有的元素,满足上式的最大值为从集合A到集合B的Hausdorff单向距离,即:
同理,向量到向量A的Hausdorff单向距离为:
取从向量A到向量B、从向量B到向量A两者之间单向Hausdorff距离较大值,作为向量A与向量B之间Hausdorff距离:
H(B,A)=max[h(A,B),h(B,A)]。
6.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,计算m条馈线的特征值向量的平均Hausdorff距离Hi.avg:
首先求取某条馈线i与其它m-1条馈线的特征值向量之间的m-1个Hausdorff距离;由于与故障馈线相关的Hausdorff距离较大,为了消除与故障馈线之间的Hausdorff距离对非故障馈线的平均Hausdorff距离的影响,先去掉馈线i的m-1个平均Hausdorff距离中的最大值然后对其余m-2个Hausdorff距离求平均值,得到馈线i与其它馈线的特征值向量之间的平均Hausdorff距离Hi.avg:
7.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,利用各馈线的历史数据,计算在各条馈线正常运行时它们之间的Hausdorff距离Hi.his中的最大值,并考虑一定的裕度,获得故障选线阈值如下:
由于电缆线路和架空线路的电抗、电纳和对地电容参数差别较大,可靠系数kH取1.3~1.5。
8.根据权利要求1所述的基于随机矩阵和Hausdorff距离的配电网故障选线方法,其特征在于,给出配电网故障选线判据:
配电网馈线故障选线判据:
如果m条馈线中存在馈线i的平均Hausdorff距离Hi.avg>Hset,则判断馈线i为故障馈线;
当母线故障时,与其连接的m条馈线的标准矩阵积的特征值分布接近,它们相互之间的Hausdorff距离较小,所有Hi.avg都小于Hset,由此产生针对母线故障的故障判据;
配电网母线故障判据:
如果与某母线连接的m个馈线中平均Hausdorff距离都满足Hi.avg<Hset(i=1,2,…,m),则判断该母线故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810376560.5A CN108572303B (zh) | 2018-04-25 | 2018-04-25 | 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810376560.5A CN108572303B (zh) | 2018-04-25 | 2018-04-25 | 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108572303A true CN108572303A (zh) | 2018-09-25 |
CN108572303B CN108572303B (zh) | 2019-08-27 |
Family
ID=63575211
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810376560.5A Active CN108572303B (zh) | 2018-04-25 | 2018-04-25 | 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108572303B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975661A (zh) * | 2019-04-22 | 2019-07-05 | 西南交通大学 | 一种基于斯皮尔曼相关系数的输电线路故障检测方法 |
CN110137911A (zh) * | 2019-06-03 | 2019-08-16 | 西南交通大学 | 一种基于Dijkstra算法的阀组最后断路器识别方法 |
CN110532329A (zh) * | 2019-09-02 | 2019-12-03 | 智慧谷(厦门)物联科技有限公司 | 一种基于区块链技术的智能手环数据处理与共享方法 |
CN110579661A (zh) * | 2019-09-19 | 2019-12-17 | 中国电力科学研究院有限公司 | 一种检测超导电缆的失超的方法和*** |
CN110988604A (zh) * | 2019-12-30 | 2020-04-10 | 武汉理工大学 | 一种配电网单相接地故障选相方法 |
CN111445108A (zh) * | 2020-03-03 | 2020-07-24 | 国网江苏省电力有限公司电力科学研究院 | 数据驱动的配电网线变关系诊断方法、装置及*** |
CN112034390A (zh) * | 2020-09-23 | 2020-12-04 | 山东科汇电力自动化股份有限公司 | 一种基于Hausdorff距离的低压配电***漏电流故障监测方法 |
CN112083286A (zh) * | 2020-09-07 | 2020-12-15 | 昆明理工大学 | 一种基于零序电流波形相似度的风电场集电线路单相接地故障选线方法 |
CN112305374A (zh) * | 2020-10-22 | 2021-02-02 | 西安工程大学 | 一种配电网单相接地故障选线方法 |
CN112345858A (zh) * | 2020-11-03 | 2021-02-09 | 西华大学 | 一种计及量测篡改攻击引发虚假故障的电网故障诊断方法 |
CN112557812A (zh) * | 2020-11-24 | 2021-03-26 | 山东理工大学 | 基于Hausdorff距离的小电流接地故障定位方法及*** |
CN115469260A (zh) * | 2022-09-29 | 2022-12-13 | 武汉格蓝若智能技术有限公司 | 一种基于Hausdorff的电流互感器异常识别方法、*** |
-
2018
- 2018-04-25 CN CN201810376560.5A patent/CN108572303B/zh active Active
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975661A (zh) * | 2019-04-22 | 2019-07-05 | 西南交通大学 | 一种基于斯皮尔曼相关系数的输电线路故障检测方法 |
CN110137911B (zh) * | 2019-06-03 | 2021-01-29 | 西南交通大学 | 一种基于Dijkstra算法的阀组最后断路器识别方法 |
CN110137911A (zh) * | 2019-06-03 | 2019-08-16 | 西南交通大学 | 一种基于Dijkstra算法的阀组最后断路器识别方法 |
CN110532329A (zh) * | 2019-09-02 | 2019-12-03 | 智慧谷(厦门)物联科技有限公司 | 一种基于区块链技术的智能手环数据处理与共享方法 |
CN110579661A (zh) * | 2019-09-19 | 2019-12-17 | 中国电力科学研究院有限公司 | 一种检测超导电缆的失超的方法和*** |
CN110579661B (zh) * | 2019-09-19 | 2022-12-13 | 中国电力科学研究院有限公司 | 一种检测超导电缆的失超的方法和*** |
CN110988604A (zh) * | 2019-12-30 | 2020-04-10 | 武汉理工大学 | 一种配电网单相接地故障选相方法 |
CN111445108A (zh) * | 2020-03-03 | 2020-07-24 | 国网江苏省电力有限公司电力科学研究院 | 数据驱动的配电网线变关系诊断方法、装置及*** |
CN112083286A (zh) * | 2020-09-07 | 2020-12-15 | 昆明理工大学 | 一种基于零序电流波形相似度的风电场集电线路单相接地故障选线方法 |
CN112083286B (zh) * | 2020-09-07 | 2021-09-28 | 昆明理工大学 | 一种风电场集电线路单相接地故障选线方法 |
CN112034390A (zh) * | 2020-09-23 | 2020-12-04 | 山东科汇电力自动化股份有限公司 | 一种基于Hausdorff距离的低压配电***漏电流故障监测方法 |
CN112034390B (zh) * | 2020-09-23 | 2023-02-07 | 山东科汇电力自动化股份有限公司 | 一种基于Hausdorff距离的低压配电***漏电流故障监测方法 |
CN112305374A (zh) * | 2020-10-22 | 2021-02-02 | 西安工程大学 | 一种配电网单相接地故障选线方法 |
CN112305374B (zh) * | 2020-10-22 | 2024-05-24 | 西安工程大学 | 一种配电网单相接地故障选线方法 |
CN112345858A (zh) * | 2020-11-03 | 2021-02-09 | 西华大学 | 一种计及量测篡改攻击引发虚假故障的电网故障诊断方法 |
CN112345858B (zh) * | 2020-11-03 | 2021-06-01 | 西华大学 | 一种计及量测篡改攻击引发虚假故障的电网故障诊断方法 |
CN112557812A (zh) * | 2020-11-24 | 2021-03-26 | 山东理工大学 | 基于Hausdorff距离的小电流接地故障定位方法及*** |
CN115469260A (zh) * | 2022-09-29 | 2022-12-13 | 武汉格蓝若智能技术有限公司 | 一种基于Hausdorff的电流互感器异常识别方法、*** |
CN115469260B (zh) * | 2022-09-29 | 2024-05-17 | 武汉格蓝若智能技术股份有限公司 | 一种基于Hausdorff的电流互感器异常识别方法、*** |
Also Published As
Publication number | Publication date |
---|---|
CN108572303B (zh) | 2019-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108572303A (zh) | 一种基于随机矩阵和Hausdorff距离的配电网故障选线方法 | |
CN103713237B (zh) | 一种电力***输电线路短路故障诊断方法 | |
CN103424669B (zh) | 一种利用故障馈线零序电流矩阵主成分分析第一主成分的选线方法 | |
CN103344875B (zh) | 谐振接地***单相接地故障分类选线方法 | |
CN102135558B (zh) | 过电压数据采集分层识别***及类型分层模式识别方法 | |
CN102854437A (zh) | 应用时频原子分解理论的小电流接地***故障选线方法 | |
CN103954884B (zh) | 一种基于形态谱分析和形态学峰谷检测可拓融合的配网故障选线方法 | |
CN108663602A (zh) | 柔性直流配电网单极故障选线与区段定位方法及*** | |
CN106199342B (zh) | 一种配电网单相接地故障选线方法 | |
CN108510005A (zh) | 一种基于卷积神经网络的配电网高阻接地故障识别方法 | |
CN107167702A (zh) | 一种配电网馈线故障类型识别方法及装置 | |
CN101924354A (zh) | 利用s变换能量样本属性的配电网故障的人工神经网络选线方法 | |
CN103941163A (zh) | 利用模糊k均值聚类的谐振接地***故障选线方法 | |
CN104614638A (zh) | 小电流***接地选线方法 | |
CN103675605A (zh) | 基于故障信号暂态相关分析的小电流接地故障选线方法 | |
CN108132422A (zh) | 一种用于容性无功补偿装置故障的监测***及方法 | |
CN106505536A (zh) | 一种基于分布参数模型下利用电流突变特性的特高压直流线路保护方法 | |
Li et al. | Application of signal processing and analysis in detecting single line‐to‐ground (SLG) fault location in high‐impedance grounded distribution network | |
CN104466922B (zh) | 一种电流极性比较式快速母线保护方法 | |
CN108599114A (zh) | 一种高压交直流混联***交流线路暂态方向保护方法 | |
CN103163430A (zh) | 一种利用复小波和ann结合的谐振接地***故障选线方法 | |
CN103267927A (zh) | 一种利用工频分量小波系数直线拟合检测的小电流接地***故障选线方法 | |
CN110244122A (zh) | 电力***谐振过电压检测与控制方法 | |
CN109596945B (zh) | 一种基于相关系数向量相似程度的配电网故障选线新方法 | |
CN107045093A (zh) | 基于快速s变换的小电流单相接地故障选线方法 |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240620 Address after: 230000 b-1018, Woye Garden commercial office building, 81 Ganquan Road, Shushan District, Hefei City, Anhui Province Patentee after: HEFEI WISDOM DRAGON MACHINERY DESIGN Co.,Ltd. Country or region after: China Address before: 610031 No. two, section 111, ring road, Chengdu, Sichuan, China Patentee before: SOUTHWEST JIAOTONG University Country or region before: China |
|
TR01 | Transfer of patent right |