CN107329044A - 一种基于电弧暂态分量的配电网单相接地故障选线方法 - Google Patents

一种基于电弧暂态分量的配电网单相接地故障选线方法 Download PDF

Info

Publication number
CN107329044A
CN107329044A CN201710531764.7A CN201710531764A CN107329044A CN 107329044 A CN107329044 A CN 107329044A CN 201710531764 A CN201710531764 A CN 201710531764A CN 107329044 A CN107329044 A CN 107329044A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
msubsup
criterion
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
Application number
CN201710531764.7A
Other languages
English (en)
Other versions
CN107329044B (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.)
State Grid Corp of China SGCC
Xuzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Xuzhou Power Supply Co of State Grid Jiangsu Electric Power 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 State Grid Corp of China SGCC, Xuzhou Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201710531764.7A priority Critical patent/CN107329044B/zh
Publication of CN107329044A publication Critical patent/CN107329044A/zh
Application granted granted Critical
Publication of CN107329044B publication Critical patent/CN107329044B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Emergency Protection Circuit Devices (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种基于电弧暂态分量的配电网单相接地故障选线方法，当发生接地故障时，求出故障相和故障时刻，选取故障前后十分之一工频周期的电流暂态信号作为采样数据，进行深度滤波和分解，得到IMF1分量，根据IMF1分量得出一阶差分极性和瞬间相位进行故障选线的判据，综合一阶差分极性判据和瞬间相位判据进行故障选线。本发明通过深度滤波和CEEMD分解抑制电弧接地故障过程中伴随的宽带白噪声和孤立脉冲噪声，根据消除干扰后的特征分量进行故障选线，有效的提高了故障选线的准确性和可靠性。

Description

一种基于电弧暂态分量的配电网单相接地故障选线方法
技术领域
本发明属于电力系统自动化技术领域，特别涉及了一种基于电弧暂态分量的配电网单相接地故障选线方法。
背景技术
我国35kV及以下的电压等级的小电流接地系统中为减小接地电流，往往采用经消弧线圈接地方式。接地电流的减小是通过消弧线圈补偿电感电流来实现的。由于小电流接地系统发生单相接地故障时故障特征比较复杂，故障信号稳态量往往很小，再加上电网本身及负荷的复杂性，使得实际中很难提取小电流谐振接地系统中的故障信号，因此其故障选线尤其困难。而且准确性往往不能得到很好的保证。
为此，大量电力工作者也提出了一些利用接地故障后的零序分量特征实现的选线方法。传统的小电流接地系统中的选线方法，如零序电流无功功率方向法、极性比较法、群体比幅比相发，这类方法在发生电弧性接地故障，尤其是间歇性电弧急切地故障时其准确性会受到较大影响，同时不太适合谐振接地系统中。而且当电网发生高阻等稳态接地故障时，由于稳态分量很小或几乎没有时，该类方法也很有可能出现误选的情况。因此，利用暂态量来进行故障选线是近几年研究的主要方向，暂态量在数值上往往远大于稳态量，因此可以在一定程度上提高选线的准确性。而如何有效利用暂态分量则成为该类方法应用的难点和重点。
发明内容
为了解决上述背景技术提出的技术问题，本发明旨在提供一种通过改变电弧暂态分量权重来实现配电网单相接地故障选线的方法，其克服了现有技术中故障选线方法准确性差、可靠性低等缺陷，该方法通过深度滤波和CEEMD分解抑制电弧接地故障过程中伴随的宽带白噪声和孤立脉冲噪声，根据消除干扰后的特征分量进行故障选线，提高故障选线的准确性和可靠性。
为了实现上述技术目的，本发明的技术方案为：
当检测到配电网发生单相接地故障时，执行以下步骤进行选线：
a、求取发生单相接地故障时的故障相和故障时刻；
b、选取故障时刻前后十分之一工频周期的电流暂态信号作为采样数据，并对采样数据依次进行深度滤波处理和CEEMD分解，得到各馈线的IMF1分量；
c、根据各馈线的IMF1分量，得出一阶差分极性和瞬间相位进行故障选线的判据；
d、综合一阶差分极性判据和瞬间相位判据进行故障选线。
优选的是，所述步骤a中，故障相和故障时刻的求取过程为，
a1、令
a2、计算
上式中，分别为三相相电动势，中幅值最小所对应的相就是故障相；
a3、分别取零序电压采样值和故障相电压采样值代入式中，分别得到满足上式的零序电压采样点对应的时刻T1和满足上式的故障相电压采样点对应的时刻T2，将这两个时刻取平均得到故障时刻：
上式中，u(j)为电压的第j个采样值，则T1、T2分别为零序电压、故障相电压采样值带入上式后对应的时刻，M为设定的阈值，N为一个周波采样点数。
优选的是，所述步骤b中，深度滤波处理的公式为：
上式中，Zn1、Zn2为深度滤波处理后的数据，f(n)为定义在{0,1,...,N-1}上的离散函数，f(n)是输入的采样数据，g(m)是定义在{0,1,...,M-1}上的离散函数，g(m)是结构元素集合，N为采样长度，n为采样点数，N1是结构元素类型数，M为大于零的常数，N>>M，gi(m)是g(m)中的第i种结构元素，N(n)和N(-n)为一对大小相等、符号相反的白噪声。
优选的是，所述结构元素集合g(m)的确定过程如下：
(1)确定使用的2种结构元素：三角形结构元素和半圆形结构元素，它们的单位结构元素分别记为BT和BS：
(2)设零序电流采样信号的极大值序列为PE＝{PEi|i＝1,2,...,NPE},NPE为极大值序列个数，局部极小值序列为NE＝{NEi|i＝1,2,...,NNE},NNE为极小值序列个数。定义局部极大值间隔和局部极小值间隔：
Dp＝{dp|dpi＝PEi+1-PEi,i＝1,2,...,NPE-1}
DN＝{dN|dNi＝NEi+1-NEi,i＝1,2,...,NNE-1}
由上式可得结构元素的长度尺度的最值如下：
Klmin＝min([(min(dpi)-1)/2),[(min(dNi)-1)/2))
Klmax＝max(((max(dpi)-1)/2],((max(dNi)-1)/2])
上式中，[·)为向上取整计算，(·]为向下取整计算。
由上式可得结构元素的长度序列Kl：
Kl＝{Klmin,Klmin+1,...,Klmax-1,Klmax}
(3)确定结构元素的高度，根据采样的暂态零序电流的局部极值的幅值大小来确定结构元素高度，设极大值序列PE和极小值序列NE的最值分别为 ppmax、ppmin、pnmax、pnmin。定义信号的局部极值的高度值是He：
He＝max(ppmax-ppmin,pnmax-pnmin)
因此，结构元素高度序列可定义为：
Hl＝{α·[He/(Kmax-Kmin+1)+(j-1)·He/(Kmax-Kmin+1)]}
其中，j＝1,2,...,Kmax-Kmin+1，α为小于1的正数，是高度比例系数，本实施例取0.05；
(4)得到各种尺度的结构元素
G1＝Hl(i).Kl(i)BT i＝1,2,...,Kmax-Kmin+1
G2＝Hl(i).Kl(i)BS i＝1,2,...,Kmax-Kmin+1
G1、G2为经由上式构造得到的三角形结构元素和半圆形结构元素集合，即g(m)＝(G1，G2)。
优选的是，所述步骤b中，CEEMD分解的具体过程如下：
(1)在深度滤波后的信号里成对地添加符号相反、大小相同的白噪声；
(2)将加入白噪声后的信号进行EMD分解，获得IMF分量；
(3)重复添加不同的白噪声，并进行EMD分解，对获得的所有IMF分量取平均值，得到最后的分解结果。
优选的是，所述步骤c中，
所述一阶差分极性判据为，当某一馈线所对应的IMF1分量的一阶差分极性唯一，且与其余馈线对应的IMF1分量的一阶差分极性均不同，则该馈线为可疑故障线路；
所述瞬时相位判据为，当某一馈线所对应的IMF1分量的瞬时相位与其余馈线对应的IMF1分量的瞬时相位的差值的绝对值均大于等于则该馈线为可疑故障线路。
优选的是，所述步骤d的具体过程为：
d1、分配一阶差分极性判据和瞬间相位判据的权值；
d2、根据两种判据的权值计算综合判据，Ck值最大的馈线即为故障线路。
8、根据权利要求7所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤d1采用的分配权重方案为：
设置暂态零序电流幅值的阈值为IZ；当暂态零序电流幅值I≥IZ，则设定瞬时相位判据的权值为0.9，一阶差分极性判据的权值为0.1；当暂态零序电流幅值I＜IZ，则设定一阶差分极性判据和瞬间相位判据的权值均为0.5。
优选的是，所述步骤d2中，Ck计算公式为：
上式中，Ck表示可疑程度综合值，sak表示第k条馈线的瞬时相位判据的可疑程度，sdk表示第k条馈线的一阶差分极性判据的可疑程度，Nm为馈线数。
优选的是，所述sak和sdk的取值方案为：
当第k条馈线满足瞬时相位判据，则sak取0.85，若第k条馈线不满足瞬时相位判据，则sak取0.15；当第k条馈线满足一阶差分极性判据，则sdk取0.72，若第k条馈线不满足一阶差分极性判据，则sdk取0.18。
本发明有益效果：
(1)本发明可以有效地抑制电力系统故障过程伴随的宽带白噪声和孤立脉冲噪声对单相接地选线地影响，能有效提取可用于电弧接地选线的暂态特征分量；
(2)本发明综合利用特征分量的一阶差分极性和瞬时相位综合判断故障线路，对提高电弧接地选线的准确性具有很好的应用价值，可有效地提高配电网供电的安全性和可靠性。
(3)本发明无需增加新的设备，只需对原装置进行算法改进即可，因此具有经济性强，便于大规模的推广应用的特点。
附图说明
图1为本发明提供的基于电弧暂态分量变权重配电网单相接地故障选线流程图；图2为本发明实施例所提供的配电系统仿真模型结构图；
图3为本发明实施例所提供的所有线路的IMF1的hht变换波形图；
图4是本发明实施例所提供的仿真环境下L3CEEMD分解波形；
图5是本发明实施例所提供的仿真环境下L1的CEEMD分解波形。
具体实施方式
为了加深对本发明的理解，下面结合附图，对本发明的具体实施方式作进一步的详细说明。
本实施例采用根据检测到的零序电压判断是否发生单相接地故障。
当时，判断发生了单相接地故障，其中，为监测到的零序电压的幅值，U0zd为零序电压的整定值。
当检测到配电网发生单相接地故障时，执行下列选线步骤。
a、若发生了故障，求取故障相和故障时刻。
首先，令
然后，计算
上式中，分别为三相相电动势；
中幅值最小所对应的相就是故障相；
分别取零序电压和故障相电压代入式中，分别得到两个时刻T1和T2，将这两个时刻取平均得到故障时刻：
上式中，u(j)为电压的第j个采样值，则T1、T2分别为零序电压、故障相电压的采样值代入上式后对应的时刻，M为设定的阈值，N为采样长度。
b、选取故障时刻前后十分之一工频周期的电流暂态信号作为采样数据，并对采样数据依次进行深度滤波处理和CEEMD分解，得到各馈线的IMF1分量。
所述深度滤波处理的方法如下：
上式中，Zn1、Zn2为深度滤波处理后的数据，f(n)为定义在{0,1,...,N-1}上的离散函数，f(n)是输入的采样数据，g(m)是定义在{0,1,...,M-1}上的离散函数，g(m)是结构元素集合，N为采样长度，n为采样点数，M为大于零的常数，N>>M，gi(m)是 g(m)中的第i种结构元素，N(n)和N(-n)为一对大小相等、符号相反的白噪声。
下文是关于结构元素集合g(m)的确定。
首先确定使用的2种结构元素：三角形结构元素和半圆形结构元素，它们的单位结构元素分别记为BT和BS：
设零序电流采样信号的极大值序列为PE＝{PEi|i＝1,2,...,NPE},NPE为极大值序列个数，局部极小值序列为NE＝{NEi|i＝1,2,...,NNE},NNE为极小值序列个数。定义局部极大值间隔和局部极小值间隔：
Dp＝{dp|dpi＝PEi+1-PEi,i＝1,2,...,NPE-1}
DN＝{dN|dNi＝NEi+1-NEi,i＝1,2,...,NNE-1}
由上式可得结构元素的长度尺度的最值如下：
Klmin＝min([(min(dpi)-1)/2),[(min(dNi)-1)/2))
Klmax＝max(((max(dpi)-1)/2],((max(dNi)-1)/2])
上式中，[·)为向上取整计算，(·]为向下取整计算。
由上式可得结构元素的长度序列Kl：
Kl＝{Klmin,Klmin+1,...,Klmax-1,Klmax}
然后就是确定结构元素的高度，根据采样的暂态零序电流的局部极值的幅值大小来确定结构元素高度，设极大值序列PE和极小值序列NE的最值分别为 ppmax、ppmin、pnmax、pnmin。定义信号的局部极值的高度值是He：
He＝max(ppmax-ppmin,pnmax-pnmin)
因此，结构元素高度序列可定义为：
Hl＝{α·[He/(Kmax-Kmin+1)+(j-1)·He/(Kmax-Kmin+1)]}
其中，j＝1,2,...,Kmax-Kmin+1，α为小于1的正数，是高度比例系数，本实施例取0.05。
于是可以得到各种尺度的结构元素
G1＝Hl(i).Kl(i)BT i＝1,2,...,Kmax-Kmin+1
G2＝Hl(i).Kl(i)BS i＝1,2,...,Kmax-Kmin+1
G1、G2为经由上式构造得到的三角形结构元素和半圆形结构元素集合，即g(m)＝(G1，G2)。
在深度滤波后，进行CEEMD分解。在信号里成对地添加符号相反、大小相同的白噪声；将加入白噪声后的信号进行EMD分解，获得IMF分量；重复添加不同的白噪声，并进行EMD分解，对获得的所有IMF分量取平均值，得到最后的分解结果。
c、根据各馈线的IMF1分量，得出一阶差分极性和瞬间相位进行故障选线的判据。
一阶差分极性判据为，当某一馈线所对应的IMF1分量的一阶差分极性唯一，且与其余馈线对应的IMF1分量的一阶差分极性均不同，则该馈线为可疑故障线路。瞬时相位判据为，当某一馈线所对应的IMF1分量的瞬时相位与其余馈线对应的IMF1分量的瞬时相位的差值的绝对值均大于等于则该馈线为可疑故障线路。
d、综合一阶差分极性判据和瞬间相位判据进行故障选线。
分配一阶差分极性判据和瞬间相位判据的权值；
设置暂态零序电流幅值的阈值IL，当暂态零序电流幅值I≥IL，则设定瞬时相位判据的权值为0.9，一阶差分极性判据的权值为0.1，当暂态零序电流幅值 I＜IZ，则设定两种判据的权值均为0.5；
根据两种判据的权值计算综合判据，Ck值最大的馈线即为故障线路：
上式中，sak表示第k条馈线的瞬时相位判据的可疑程度，当第k条馈线满足瞬时相位判据，则sak取0.85，若第k条馈线不满足瞬时相位判据，则sak取0.15， sdk表示第k条馈线的一阶差分极性判据的可疑程度，当第k条馈线满足一阶差分极性判据，则sdk取0.72，若第k条馈线不满足一阶差分极性判据，则sdk取0.28； Nm为馈线数。
在Matlab/Simulink仿真环境下，供配电系统原理图的仿真模型如图2所示。该系统采用中性点经消弧线圈过补偿方式，过补偿度为5％，经计算消弧线圈电感L＝0.8697H，R＝30Ω。在MATLAB仿真中，线路的正序参数为：R1＝0.01273 Ω/Km，L1＝9.337×10-4H/Km,C1＝1.274×10-8F/Km；零序参数为R0＝0.3863Ω /Km,L0＝4.1264×10-3H/Km，C0＝7.751×10-9F/Km。模型中线路长度分别为103Km、 175Km、151Km。电源采用的是“Three-phasesource”模型，电源功率为10000MV.A, 电压是110Kv，输出电压为11Kv。变压器额定容量为Sn＝20MV.A,短路电压 Us％＝10.5,短路损耗ΔPs＝135Kw，空载损耗ΔP0＝22Kw，空载电流I0％＝0.8，变比Kt＝110/10，高低压绕组均为Y型联结。线路负荷均采用“Three-phaseSeries RLC Load”模型，有功负荷分别为1MW、0.2MW、2MW，采样频率设为100Khz,初始故障时间先设为0.04S。选取故障时刻前后十分之一工频周期的电流暂态信号作为采样数据，并对采样数据依次进行深度滤波处理和CEEMD分解，得到各馈线的 IMF1分量，图3为所有线路的IMF1的hht变换波形图。图4、图5为过渡电阻 Rf＝0，故障时间t＝0.04s时的故障线路L3和一条非故障线路L1的部分CEEMD 分解结果波形图。从图4、图5可以看出，非故障线路与故障线路在故障后的1/4 周期相位相差180度。
表1实验结果：
上表的选线结果证明了基于电弧暂态分量的配电网单相接地故障选线方法能实现准确选线。
以上所术仅是本发明的一种实施方式，应当说明，所属领域的普通技术人员应当理解，依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：
当检测到配电网发生单相接地故障时，执行以下步骤进行选线：
a、求取发生单相接地故障时的故障相和故障时刻；
b、选取故障时刻前后十分之一工频周期的电流暂态信号作为采样数据，并对采样数据依次进行深度滤波处理和CEEMD分解，得到各馈线的IMF1分量；
c、根据各馈线的IMF1分量，得出一阶差分极性和瞬间相位进行故障选线的判据；
d、综合一阶差分极性判据和瞬间相位判据进行故障选线。
2.根据权利要求1所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤a中，故障相和故障时刻的求取过程为，
a1、令
a2、计算
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <msub> <mover> <mi>E</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> <mo>+</mo> <msubsup> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>b</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <msub> <mover> <mi>E</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>b</mi> </msub> <mo>+</mo> <msubsup> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>c</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <msub> <mover> <mi>E</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>c</mi> </msub> <mo>+</mo> <msubsup> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced>
上式中，分别为三相相电动势，中幅值最小所对应的相就是故障相；
a3、分别取零序电压采样值和故障相电压采样值代入式中，分别得到满足上式的零序电压采样点对应的时刻T1和满足上式的故障相电压采样点对应的时刻T2，将这两个时刻取平均得到故障时刻：
<mrow> <mi>T</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> </mrow> <mn>2</mn> </mfrac> </mrow>
上式中，u(j)为电压的第j个采样值，则T1、T2分别为零序电压、故障相电压采样值带入上式后对应的时刻，M为设定的阈值，N为一个周波采样点数。
3.根据权利要求1所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤b中，深度滤波处理的公式为：
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mn>1</mn> </msub> </munderover> <msubsup> <mi>Z</mi> <mi>n</mi> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>Z</mi> <mi>n</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Z</mi> <mrow> <mi>n</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mn>1</mn> </msub> </munderover> <msubsup> <mi>Z</mi> <mi>n</mi> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>Z</mi> <mi>n</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>Z</mi> <mi>n</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mi>min</mi> <mo>{</mo> <mi>max</mi> <mo>{</mo> <mi>max</mi> <mo>{</mo> <mi>min</mi> <mo>&amp;lsqb;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>+</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>}</mo> <mo>,</mo> <mi>m</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>Z</mi> <mi>n</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msubsup> <mo>=</mo> <mi>max</mi> <mo>{</mo> <mi>min</mi> <mo>{</mo> <mi>min</mi> <mo>{</mo> <mi>max</mi> <mo>&amp;lsqb;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>+</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>}</mo> <mo>,</mo> <mi>m</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced>
上式中，Zn1、Zn2为深度滤波处理后的数据，f(n)为定义在{0,1,...,N-1}上的离散函数，f(n)是输入的采样数据，g(m)是定义在{0,1,...,M-1}上的离散函数，g(m)是结构元素集合，N为采样长度，n为采样点数，N1是结构元素类型数，M为大于零的常数，N>>M，gi(m)是g(m)中的第i种结构元素，N(n)和N(-n)为一对大小相等、符号相反的白噪声。
4.根据权利要求3所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述结构元素集合g(m)的确定过程如下：
(1)确定使用的2种结构元素：三角形结构元素和半圆形结构元素，它们的单位结构元素分别记为BT和BS：
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>B</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mo>|</mo> <mi>i</mi> <mo>|</mo> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>B</mi> <mi>S</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>i</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> </mfenced>
(2)设零序电流采样信号的极大值序列为PE＝{PEi|i＝1,2,...,NPE},NPE为极大值序列个数，局部极小值序列为NE＝{NEi|i＝1,2,...,NNE},NNE为极小值序列个数。定义局部极大值间隔和局部极小值间隔：
Dp＝{dp|dpi＝PEi+1-PEi,i＝1,2,...,NPE-1}
DN＝{dN|dNi＝NEi+1-NEi,i＝1,2,...,NNE-1}
由上式可得结构元素的长度尺度的最值如下：
Klmin＝min([(min(dpi)-1)/2),[(min(dNi)-1)/2))
Klmax＝max(((max(dpi)-1)/2],((max(dNi)-1)/2])
上式中，[·)为向上取整计算，(·]为向下取整计算。
由上式可得结构元素的长度序列Kl：
Kl＝{Klmin,Klmin+1,...,Klmax-1,Klmax}
(3)确定结构元素的高度，根据采样的暂态零序电流的局部极值的幅值大小来确定结构元素高度，设极大值序列PE和极小值序列NE的最值分别为ppmax、ppmin、pnmax、pnmin。定义信号的局部极值的高度值是He：
He＝max(ppmax-ppmin,pnmax-pnmin)
因此，结构元素高度序列可定义为：
Hl＝{α·[He/(Kmax-Kmin+1)+(j-1)·He/(Kmax-Kmin+1)]}
其中，j＝1,2,...,Kmax-Kmin+1，α为小于1的正数，是高度比例系数，本实施例取0.05；
(4)得到各种尺度的结构元素
G1＝Hl(i).Kl(i)BT i＝1,2,...,Kmax-Kmin+1
G2＝Hl(i).Kl(i)BS i＝1,2,...,Kmax-Kmin+1
G1、G2为经由上式构造得到的三角形结构元素和半圆形结构元素集合，即g(m)＝(G1，G2)。
5.根据权利要求1所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤b中，CEEMD分解的具体过程如下：
(1)在深度滤波后的信号里成对地添加符号相反、大小相同的白噪声；
(2)将加入白噪声后的信号进行EMD分解，获得IMF分量；
(3)重复添加不同的白噪声，并进行EMD分解，对获得的所有IMF分量取平均值，得到最后的分解结果。
6.根据权利要求1所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤c中，
一阶差分极性判据为，当某一馈线所对应的IMF1分量的一阶差分极性唯一，且与其余馈线对应的IMF1分量的一阶差分极性均不同，则该馈线为可疑故障线路；
瞬时相位判据为，当某一馈线所对应的IMF1分量的瞬时相位与其余馈线对应的IMF1分量的瞬时相位的差值的绝对值均大于等于则该馈线为可疑故障线路。
7.根据权利要求1所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤d的具体过程为：
d1、分配一阶差分极性判据和瞬间相位判据的权值；
d2、根据两种判据的权值计算综合判据，Ck值最大的馈线即为故障线路。
8.根据权利要求7所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤d1采用的分配权重方案为：
设置暂态零序电流幅值的阈值为IZ；当暂态零序电流幅值I≥IZ，则设定瞬时相位判据的权值为0.9，一阶差分极性判据的权值为0.1；当暂态零序电流幅值I＜IZ，则设定一阶差分极性判据和瞬间相位判据的权值均为0.5。
9.根据权利要求7所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述步骤d2中，Ck计算公式为：
<mrow> <msub> <mi>C</mi> <mi>k</mi> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <mrow> <mo>(</mo> <mn>0.5</mn> <msub> <mi>s</mi> <mrow> <mi>a</mi> <mi>k</mi> </mrow> </msub> <mo>+</mo> <mn>0.5</mn> <msub> <mi>s</mi> <mrow> <mi>d</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>I</mi> <mo>&lt;</mo> <msub> <mi>I</mi> <mi>Z</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <mrow> <mo>(</mo> <mn>0.9</mn> <msub> <mi>s</mi> <mrow> <mi>a</mi> <mi>k</mi> </mrow> </msub> <mo>+</mo> <mn>0.1</mn> <msub> <mi>s</mi> <mrow> <mi>d</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>I</mi> <mo>&amp;GreaterEqual;</mo> <msub> <mi>I</mi> <mi>Z</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
上式中，Ck表示可疑程度综合值，sak表示第k条馈线的瞬时相位判据的可疑程度，sdk表示第k条馈线的一阶差分极性判据的可疑程度，Nm为馈线数。
10.根据权利要求9所述的基于电弧暂态分量的配电网单相接地故障选线方法，其特征在于：所述sak和sdk的取值方案为：
当第k条馈线满足瞬时相位判据，则sak取0.85，若第k条馈线不满足瞬时相位判据，则sak取0.15；当第k条馈线满足一阶差分极性判据，则sdk取0.72，若第k条馈线不满足一阶差分极性判据，则sdk取0.18。
CN201710531764.7A 2017-06-30 2017-06-30 一种基于电弧暂态分量的配电网单相接地故障选线方法 Active CN107329044B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710531764.7A CN107329044B (zh) 2017-06-30 2017-06-30 一种基于电弧暂态分量的配电网单相接地故障选线方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710531764.7A CN107329044B (zh) 2017-06-30 2017-06-30 一种基于电弧暂态分量的配电网单相接地故障选线方法

Publications (2)

Publication Number Publication Date
CN107329044A true CN107329044A (zh) 2017-11-07
CN107329044B CN107329044B (zh) 2020-08-04

Family

ID=60197770

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710531764.7A Active CN107329044B (zh) 2017-06-30 2017-06-30 一种基于电弧暂态分量的配电网单相接地故障选线方法

Country Status (1)

Country Link
CN (1) CN107329044B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108254657A (zh) * 2018-03-28 2018-07-06 山东大学 基于暂态能量分析的配电网小电流接地故障区段定位方法
CN108919046A (zh) * 2018-05-15 2018-11-30 国网山东省电力公司济南供电公司 一种配电网接地故障试停线路序位决策方法及系统
CN110501631A (zh) * 2019-08-19 2019-11-26 重庆大学 一种在线间歇故障检测与诊断方法
CN111257687A (zh) * 2020-01-19 2020-06-09 太原理工大学 一种基于改进的meemd的单相接地故障选线方法
CN111337767A (zh) * 2020-02-25 2020-06-26 河北工业大学 一种谐振波减速器故障分析方法
CN111650476A (zh) * 2020-07-06 2020-09-11 国网江苏省电力有限公司沛县供电分公司 基于采样值法的配电网单相电弧性接地故障选线方法
CN111948568A (zh) * 2019-04-30 2020-11-17 国家电网有限公司 基于瞬时相位一致的电网线路故障处理方法
CN112255493A (zh) * 2020-09-07 2021-01-22 西安理工大学 一种多判据综合投票配电网故障选线方法
CN115184729A (zh) * 2022-07-20 2022-10-14 国网辽宁省电力有限公司电力科学研究院 一种自动检测并处理间歇性电弧接地故障的模拟系统和方法
CN117434406A (zh) * 2023-12-20 2024-01-23 天津航空机电有限公司 一种基于互补集合经验模态分解的电弧故障检测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551433A (zh) * 2009-05-05 2009-10-07 昆明理工大学 一种利用hht检测的配电网馈出线路故障选线方法
US20120086459A1 (en) * 2010-10-12 2012-04-12 Howard University Apparatus and Method for Fault Detection and Location Determination
CN103344875A (zh) * 2013-07-02 2013-10-09 福州大学 谐振接地系统单相接地故障分类选线方法
CN104237731A (zh) * 2014-09-25 2014-12-24 福州大学 基于eemd与能量法的谐振接地配电网单相接地故障选线方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551433A (zh) * 2009-05-05 2009-10-07 昆明理工大学 一种利用hht检测的配电网馈出线路故障选线方法
US20120086459A1 (en) * 2010-10-12 2012-04-12 Howard University Apparatus and Method for Fault Detection and Location Determination
CN103344875A (zh) * 2013-07-02 2013-10-09 福州大学 谐振接地系统单相接地故障分类选线方法
CN104237731A (zh) * 2014-09-25 2014-12-24 福州大学 基于eemd与能量法的谐振接地配电网单相接地故障选线方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
田书等: "基于ＨＨＴ 的谐振接地配电网暂态量故障选线研究", 《测控技术》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108254657A (zh) * 2018-03-28 2018-07-06 山东大学 基于暂态能量分析的配电网小电流接地故障区段定位方法
CN108254657B (zh) * 2018-03-28 2019-11-26 山东大学 基于暂态能量分析的配电网小电流接地故障区段定位方法
CN108919046A (zh) * 2018-05-15 2018-11-30 国网山东省电力公司济南供电公司 一种配电网接地故障试停线路序位决策方法及系统
CN111948568A (zh) * 2019-04-30 2020-11-17 国家电网有限公司 基于瞬时相位一致的电网线路故障处理方法
CN110501631B (zh) * 2019-08-19 2021-10-29 重庆大学 一种在线间歇故障检测与诊断方法
CN110501631A (zh) * 2019-08-19 2019-11-26 重庆大学 一种在线间歇故障检测与诊断方法
CN111257687A (zh) * 2020-01-19 2020-06-09 太原理工大学 一种基于改进的meemd的单相接地故障选线方法
CN111257687B (zh) * 2020-01-19 2022-03-22 太原理工大学 一种基于改进的meemd的单相接地故障选线方法
CN111337767A (zh) * 2020-02-25 2020-06-26 河北工业大学 一种谐振波减速器故障分析方法
CN111650476A (zh) * 2020-07-06 2020-09-11 国网江苏省电力有限公司沛县供电分公司 基于采样值法的配电网单相电弧性接地故障选线方法
CN112255493A (zh) * 2020-09-07 2021-01-22 西安理工大学 一种多判据综合投票配电网故障选线方法
CN115184729A (zh) * 2022-07-20 2022-10-14 国网辽宁省电力有限公司电力科学研究院 一种自动检测并处理间歇性电弧接地故障的模拟系统和方法
CN117434406A (zh) * 2023-12-20 2024-01-23 天津航空机电有限公司 一种基于互补集合经验模态分解的电弧故障检测方法
CN117434406B (zh) * 2023-12-20 2024-04-09 天津航空机电有限公司 一种基于互补集合经验模态分解的电弧故障检测方法

Also Published As

Publication number Publication date
CN107329044B (zh) 2020-08-04

Similar Documents

Publication Publication Date Title
CN107329044A (zh) 一种基于电弧暂态分量的配电网单相接地故障选线方法
CN103018627B (zh) 非有效接地系统故障类型自适应接地选线方法
CN106872852B (zh) 基于零序电压的配电网单相断线接地复故障类型诊断方法
CN103675605B (zh) 一种基于故障信号暂态相关分析的配电网小电流接地故障选线方法
CN109307824B (zh) 一种基于聚类的配电网单相接地故障区段定位方法
CN103346542B (zh) 基于分布参数的高压直流输电线路高阻接地故障识别方法
CN101814731A (zh) 一种应用非工频暂态分量的配电网故障选线方法
CN103207352A (zh) 利用选线阻抗幅值特性实现配电网单相接地故障选线方法
CN103226176A (zh) 一种配电网单相接地故障选线方法
CN107765139A (zh) 一种高准确率的谐振接地系统单相接地故障选线方法
Biswas et al. A dual-time transform assisted intelligent relaying scheme for the STATCOM-compensated transmission line connecting wind farm
CN101924354A (zh) 利用s变换能量样本属性的配电网故障的人工神经网络选线方法
CN110488155A (zh) 一种应用于柔性接地系统的单相接地故障选线方法
CN103207354A (zh) 基于选线系数最大原理的配电网单相接地故障选线方法
CN107167726A (zh) 一种断路器内部击穿电弧建模方法
Yuan et al. Study on fault line detection methods based on multi‐feature fusion in distribution systems
CN105606955A (zh) 一种基于数值微分与经验模态分解的故障线路判别方法
CN105445618B (zh) 一种小电流接地系统故障选线方法及装置
CN112234579B (zh) 大型水轮发电机的注入式定子接地保护方法及系统
CN110687390A (zh) 一种基于可拓理论的小电流融合选线方法
Zadeh et al. A compensation scheme for CVT transient effects using artificial neural network
CN113552441B (zh) 一种单相接地故障检测方法及装置
CN110208655B (zh) 一种消弧线圈接地方式下的配电网单相接地故障模型分析方法
CN112649694B (zh) 一种小电流接地系统单相接地故障的判定方法
CN112234578A (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