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

本发明属于电力系统自动化技术领域，特别涉及了一种基于电弧暂态分量的配电网单相接地故障选线方法。

背景技术：

我国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)bti＝1,2,...,kmax-kmin+1

g2＝hl(i).kl(i)bsi＝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)bti＝1,2,...,kmax-kmin+1

g2＝hl(i).kl(i)bsi＝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-phaseseriesrlcload”模型，有功负荷分别为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实验结果：

上表的选线结果证明了基于电弧暂态分量的配电网单相接地故障选线方法能实现准确选线。

以上所术仅是本发明的一种实施方式，应当说明，所属领域的普通技术人员应当理解，依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。