基于VMD的小电流接地故障融合选线方法

孙孔明, 李玉敦, 范荣奇, 薛永端, 庞清乐, 徐宪泽

孙孔明,李玉敦,范荣奇,等. 基于VMD的小电流接地故障融合选线方法[J]. 工矿自动化,2023,49(3):115-123. DOI: 10.13272/j.issn.1671-251x.2022080082
引用本文: 孙孔明,李玉敦,范荣奇,等. 基于VMD的小电流接地故障融合选线方法[J]. 工矿自动化,2023,49(3):115-123. DOI: 10.13272/j.issn.1671-251x.2022080082
SUN Kongming, LI Yudun, FAN Rongqi, et al. Fusion fault line selection method of small current grounding fault based on VMD[J]. Journal of Mine Automation,2023,49(3):115-123. DOI: 10.13272/j.issn.1671-251x.2022080082
Citation: SUN Kongming, LI Yudun, FAN Rongqi, et al. Fusion fault line selection method of small current grounding fault based on VMD[J]. Journal of Mine Automation,2023,49(3):115-123. DOI: 10.13272/j.issn.1671-251x.2022080082

基于VMD的小电流接地故障融合选线方法

基金项目: 国网山东省电力公司科技项目(5206002000TZ)。
详细信息
    作者简介:

    孙孔明(1989—),男,山东潍坊人,博士,主要从事继电保护技术、配电网优化与控制工作,E-mail: skming2008@163.com

  • 中图分类号: TD608

Fusion fault line selection method of small current grounding fault based on VMD

  • 摘要: 针对目前煤矿配电网故障选线方法在相关故障特征不明显时存在故障选线失效,基于单一模态分量和单一故障特征的故障选线方法的选线准确度较低等问题,提出了一种基于变分模态分解(VMD)的小电流接地故障融合选线方法。利用VMD将母线中各出线的故障零序电流分解为多个模态分量,根据模态分量的故障特征确定VMD层数,并选取故障特征明显的模态分量作为故障选线的有效模态分量;分别计算各出线故障零序电流有效模态分量的暂态能量和波形相似度;根据各出线有效模态分量的暂态能量占比和波形相似度占比,构建基于暂态能量的故障选线判据和基于波形相似度的故障选线判据,并将2种故障选线判据融合,形成基于VMD的故障融合选线算法。利用电磁暂态仿真软件ATP/EMTP搭建煤矿配电网模型,在不同接地故障电阻、故障初相角和故障位置的单相接地故障场景下,对所提出的故障融合选线方法进行验证,结果表明:在配电网发生各种单相接地故障时,基于VMD的小电流接地故障融合选线方法不受故障位置的影响,较能量法和相关性聚类法的故障选线正确率分别提高了17%和50%,且不受故障类型影响,可应用于小电流接地故障选线。
    Abstract: At present, the fault line selection method of coal mine distribution network has the problem of fault line selection failure when the relevant fault features are not obvious. The fault line selection method based on single modal component and single fault feature has low accuracy. In order to solve the above problems, a method of small current grounding fusion fault line selection based on variational mode decomposition (VMD) is proposed. VMD is used to decompose the fault zero-sequence current of each outgoing line in the bus into multiple modal components. The layer number of the VMD is determined according to the fault features of the modal components. The modal components with obvious fault features are selected as effective modal components for fault line selection. The transient energy and the waveform similarity of the effective modal component of the zero sequence current of each outgoing line fault are calculated respectively. According to the proportion of transient energy and the proportion of waveform similarity of effective modal components of each outgoing line, a fault line selection criterion based on transient energy and a fault line selection criterion based on waveform similarity are constructed. The two fault line selection criteria are fused to form a fault fusion line selection algorithm based on VMD. A coal mine distribution network model is built by using the electromagnetic transient simulation software ATP/EMTP. The proposed fusion fault line selection method is verified under single-phase-to-ground fault scenarios with different ground fault resistances, fault initial phase angles and fault locations. The results show that when various single-phase ground faults occur in the distribution network, the fusion line selection method of small-current grounding fault based on VMD is not affected by the fault location. The fault line selection accuracy is respectively improved by 17% and 50% compared with the energy method and the correlation clustering method. The fusion line selection method is not affected by the fault type, and can be applied to the small current grounding fault line selection.
  • 随着煤矿电气化和现代化水平的不断提高,先进的综采技术不断推广,煤矿供电容量不断增加,供电距离持续增长,导致煤矿配电网发生故障的几率不断提高[1]。在配电线路故障中,80%以上为单相接地故障[2]。煤矿配电网普遍采用中性点经消弧线圈接地方式,当发生单相接地故障时,受故障信号微弱、故障类型复杂和运行方式多变等因素影响,故障特征提取尤为困难,因此,单相接地故障选线一直是配电网领域的难题[3]。单相接地故障易引发设备绝缘破坏和电弧,导致故障蔓延,甚至引发安全事故。因此,研究快速而准确的故障选线技术对建设煤矿智能电网具有重要意义[4]

    国内外专家学者提出了多种故障选线方法,主要包括稳态信息法[5-6]、暂态信息法[7-8]、谐波法[9]、信号注入法[10]、行波法[11-12]和智能化方法[13-14]。在谐振接地系统中,消弧线圈一般采用过补偿运行方式,故障线路和非故障线路的稳态信息类似,因此,稳态信息法不适用于煤矿供电系统。暂态信息法仅需提取故障初始阶段的故障信息,可实现快速故障选线,但当某些暂态故障特征不明显时,故障选线会失效。谐波法通过检测故障电流中的谐波信息进行故障选线,但当谐波含量较低时无法实现故障选线。信号注入法需要沿线检测注入信号,不适用于煤矿供电系统。由于行波波头检测困难,行波法的故障选线准确度较低。智能化方法可以实现多种故障特征的融合故障选线,且具有容错能力,但该方法对供电系统的适应性差,且故障选线速度较慢。为了提高煤矿配电网故障选线速度和准确度,国内外学者提出了多种故障选线新方法。文献[1]提出了一种基于小波包的矿井供电系统单相接地故障选线方法,该方法通过比较线路特征频带小波包分解系数模最大值极性确定故障线路,可以实现快速故障选线,但在有些故障场景下极性特征不明显,无法实现正确的故障选线。文献[15]针对矿井电网附加直流源系统发生单相接地故障时故障线路暂态零序电流的直流分量大于非故障支路的特征,采用最小二乘矩阵束算法提取故障暂态零序电流中的直流分量,将各支路的直流分量组成矩阵,再利用Romanovsky准则进行判别,但直流分量受故障初相角影响较大,故障选线准确度较低。文献[16]对采集到的零序电流信号进行集合经验模态分解(Ensemble Empirical Mode Decomposition,EEMD),提取分解得到的高频分量进行小波阈值去噪,对去噪后的高频分量进行Hilbert变换并提取瞬时幅值,根据能量占比最大原则判断故障线路,但在各出线的电气参数差别较大时,故障线路幅值特征可能不明显,故障选线困难。针对上述2种故障选线方法存在的问题,文献[17]融合去噪后暂态零序电流的幅值和极性实现故障选线,克服了基于幅值和极性单一特征故障选线的局限性,提高了故障选线准确度,但存在经验模态分解(Empirical Mode Decomposition,EMD)固有的模态混叠现象,导致故障特征不明显,降低了故障选线准确度。为了克服EMD的缺点,有些专家学者将变分模态分解(Variational Mode Decomposition,VMD)应用于故障选线。文献[18]利用VMD对零序电流进行分解,分时段计算各分量能量,利用神经网络实现故障选线,但在能量特征不明显时,故障选线失效。文献[19]采用VMD分离暂态零序电流中的工频分量和噪声信号,获得暂态非工频分量,通过比较非工频分量的相关性实现故障选线,但在相关性故障特征不明显时故障选线失效。

    针对上述问题,本文提出了基于VMD的小电流接地故障融合选线方法。在分析小电流接地故障暂态特征的基础上,利用VMD提取故障零序电流的高频振荡分量,选择合适的模态分量,构建基于能量的故障选线判据和基于波形相似度的故障选线判据,并将2种故障选线判据融合,形成基于VMD的融合故障选线算法,通过仿真分析验证所提出融合选线方法的有效性。

    小电流接地系统包括中性点不接地系统和中性点经消弧线圈接地系统(也称谐振接地系统)。中性点不接地系统的故障选线较容易,本文不作分析,重点分析谐振接地系统的单相接地故障暂态特征。发生单相接地故障时谐振接地系统的零序等效电路如图1所示,其中C为系统三相对地分布电容总和,L,R分别为系统在零序电路中的等效电感和电阻; Lp为消弧线圈的等效电感,uf为故障点f的虚拟电源,ifiCiL分别为故障零序电流、故障电容电流和消弧线圈电感电流[7]

    图  1  单相接地故障零序等效电路
    Figure  1.  Zero sequence equivalent circuit of single-phase ground fault

    谐振接地系统发生单相接地故障时暂态等值电路为三阶电路,为便于求解,将其分为低阻故障和高阻故障进行降阶。

    发生低阻接地故障时,电容充放电速度较快,主谐振频率较大,消弧线圈的等效阻抗远大于系统三相对地分布电容总和的阻抗,因此,忽略消弧线圈对故障零序电流的影响,将三阶等值电路等效为等效电感L和系统三相对地分布电容总和C构成的二阶串联电路。设故障点虚拟电源${u_{\rm{f}}} = U\sin \left( {\omega t + \varphi } \right)$U为系统相电压峰值,ω为系统基准频率,t为时间,φ为故障初相角,则故障零序电流if由工频交流分量和高频振荡分量组成。

    $$ \begin{gathered} {i_{\rm{f}}} = U\omega C \bigg(\frac{{{1 {}}}}{{LC\omega {\omega _{{\rm{f}}1}}}}\sin \varphi \sin ({\omega _{{\rm{f}}1}}t )+ \frac{{{\gamma _1}}}{{{\omega _{{\rm{f}}1}}}}\cos \varphi \sin ({\omega _{{\rm{f1}}}}t )- \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \cos \varphi \cos( {\omega _{{\rm{f}}1}}t )\bigg)\exp ( - {\gamma _1}t) + U\omega C\cos \left( {\omega t + \varphi } \right) \\ \end{gathered} $$ (1)

    式中:ωf1为主谐振频率,${\omega _{{\rm{f1}}}} = \sqrt {\dfrac{1}{{LC}} - {{\left( { \dfrac{R}{{2L}} } \right)}^2}}$γ1为主谐振分量的衰减系数,γ1=R/(2L)。

    发生高阻接地故障时,电容充放电速度较慢,主谐振频率较低,消弧线圈不能忽略,等效电阻R远大于等效电感L,因此,忽略等效电感L对故障零序电流的影响,将三阶等值电路等效为消弧线圈Lp和系统三相对地分布电容总和C构成的二阶并联电路。当$R < \dfrac{1}{2}\omega {L_{\text{p}}}$时,为欠阻尼情况,故障零序电流if由工频交流分量和直流衰减分量组成。

    $$ \begin{gathered} {i_{\rm{f}}} = {i_{\rm{L}}} + {i_{\rm{C}}} = \left( {1 + {L_{\text{p}}}Cp_1^2} \right){A_1}\exp \left( {{p_1}t} \right) + \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {1 + {L_{\text{p}}}Cp_2^2} \right){A_2}\exp \left( {{p_2}t} \right) + \left( {1 - {\omega ^2}} \right)B\sin \left( {\omega t + \phi } \right) \\ \end{gathered} $$ (2)

    式中:p1p2为2个不相等负实根;A1A2分别为消弧线圈电感电流中2个直流衰减分量的幅值;B为消弧线圈电感电流中工频交流分量的幅值; $ \phi $为消弧线圈电感电流中工频交流分量的初相位[20]

    $R > \dfrac{1}{2}\omega {L_{\text{p}}}$时,为过阻尼情况,故障零序电流if由工频交流分量和高频振荡分量组成。

    $$ \begin{gathered} {i_{\rm{f}}} = {i_{\rm{L}}} + {i_{\rm{C}}} = \exp ( - {\gamma _2}t)({A_3}\cos ({\omega _{{\rm{f2}}}}t) + {A_4}\sin ({\omega _{{\rm{f2}}}}t)) \times \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} (1 + {L_{\text{p}}}C\gamma _2^2) + B\sin (\omega t + \phi )(1 - {L_{\text{p}}}C{\omega ^2}) - {L_{\text{p}}}C \times \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 2{\gamma _2}\exp ( - {\gamma _2}t)( - {\omega _{{\rm{f2}}}}{A_3}\sin ({\omega _{{\rm{f2}}}}t) + {\omega _{{\rm{f2}}}}{A_4}\cos ({\omega _{{\rm{f2}}}}t)) + \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {L_{\text{p}}}C \exp ( - {\gamma _2}t)(- \omega _{{\rm{f2}}}^2{A_3}\cos ({\omega _{{\rm{f2}}}}t) - \omega _{{\rm{f2}}}^2{A_4}\sin ({\omega _{{\rm{f2}}}}t)) \\ \end{gathered} $$ (3)

    式中:γ2为主谐振分量的衰减系数,${\gamma _2} = \dfrac{1}{{2RC}}$A3A4分别为消弧线圈电感电流中2个谐振分量的幅值[20]ωf2为主谐振频率,${\omega _{{\rm{f2}}}} = \sqrt {\dfrac{1}{{{L_{\text{p}}}C}} - {{\left( { \dfrac{1}{{2RC}} } \right)}^2}}$

    由上述分析可知,故障零序电流由直流衰减分量、工频交流分量和高频振荡分量组成,实际的故障零序电流中可能还存在高频干扰。直流衰减分量受故障初相角影响较大,工频交流分量受消弧线圈电流影响较大,均不能作为故障特征进行故障选线。而故障线路和健全线路的高频振荡分量在幅值和极性上存在较大差异,可作为暂态故障特征进行故障选线。为了消除直流衰减分量、工频交流分量和高频干扰对暂态故障特征的影响,采用VMD提取故障零序电流中的高频振荡分量作为有效暂态故障特征,以提高故障选线的准确度。

    采用文献[21]的煤矿10 kV配电网模型,利用电磁暂态仿真软件ATP/EMTP搭建仿真模型,如图2所示,其中S1S2S3为架空线路,长度分别为1.2,2.2,2 km,单位长度正序电阻为0.096 Ω/km,单位长度正序电感为1.22 mH/km,单位长度正序电容为0.011 μF/km,单位长度零序电阻为0.23 Ω/km,单位长度零序电感为3.66 mH/km,单位长度零序电容为0.007 μF/km。S4S11为电缆线路,长度分别为1.5,1.8,5,1.3,2,3,8,1.5 km,单位长度正序电阻为0.11 Ω/km,单位长度正序电感为0.52 mH/km,单位长度正序电容为0.29 μF/km,单位长度零序电阻为0.34 Ω/km,单位长度零序电感为1.54 mH/km,单位长度零序电容为0.19 μF/km。消弧线圈电感Lp=0.560 5 H,串联电阻RL=10 Ω。D1−D7为负载,负载功率分别为0.81,2,0.68,0.315,0.68,0.315,1.25 MW。f1,f2,f3为故障点位置。

    图  2  煤矿配电网模型
    Figure  2.  Coal mine distribution network model

    VMD是一种自适应的非平稳信号处理方法,其基本原理是根据预设模态分解个数K自适应地匹配各模态的中心频率和带宽,将信号分解为K个不同频段的本征模态函数(Intrinsic Mode Function,IMF)分量。VMD能够克服EMD的模态混叠问题,精确地分解出故障零序电流信号中的高频振荡分量,因此,采用VMD提取故障零序电流的故障特征[18-19]。VMD算法主要包括变分问题的构造和求解。

    变分问题构造为如下有约束的变分模型。

    $$ \left\{ \begin{gathered} \mathop {\min }\limits_{\left\{ {{i_k}} \right\},\left\{ {{\omega _k}} \right\}} \left( {\sum\limits_{k = 1}^K {\left\| {\left. {{\frac{\partial }{{{\partial t}}}}\left[ {\left( {\delta + \frac{{\rm{j}}}{{{\text{π}} t}}} \right) * {i_k}} \right]\exp \left( { - {\rm{j}}{\omega _k}t} \right)} \right\|_2^2} \right.} } \right) \\ {\text{s}}{\text{. t}}{\text{. }}\sum\limits_{k = 1}^K {{i_k}} = {i_{\rm{f}}} \\ \end{gathered} \right. $$ (4)

    式中:ik 为故障零序电流的第k个模态分量,k=1, 2, …, Kωk ik的中心频率;δ为冲激函数;*为卷积运算。

    为了将有约束的变分问题转换为无约束的变分问题,引入拉格朗日乘法算子λ和二次惩罚因子α,得到增广拉格朗日表达式

    $$ \begin{split} l\left( {\left\{ {{i_k}} \right\},\left\{ {{\omega _k}} \right\},\lambda } \right) = & \alpha \sum\limits_{k = 1}^K {\left\| {{\frac{\partial }{{{\partial t}}} }\left[ {\left( {\delta + \frac{{\rm{j}}}{{{\text{π}} t}}} \right) * {i_k}} \right]\exp \left( { - {\rm{j}}{\omega _k}t} \right)} \right\|} _2^2 + \\& \left\| {{i_{\rm{f}}} - \sum\limits_{k = 1}^K {{i_k}} } \right\|_2^2 + \left\langle {\lambda ,{i_{\rm{f}}} - \sum\limits_{k = 1}^K {{i_k}} } \right\rangle \\[-18pt] \end{split} $$ (5)

    式中〈 〉为内积运算。

    采用乘法算子交替方向法并结合傅里叶变换求解式(5),得到该表达式的“鞍点”及变分模型的最优解。各模态分量、中心频率和拉格朗日算子的更新表达式为

    $$ \hat i_k^{n + 1}(\omega ) = \dfrac{{{{\hat i}_{\rm{f}}}(\omega ) - \displaystyle \sum\limits_{w = 1,w \ne k}^K {\hat i_w^n(\omega ) + \frac{{{{\hat \lambda }^n}(\omega )}}{2}} }}{{1 + 2\alpha {{(\omega - \omega _k^n)}^2}}} $$ (6)
    $$ \omega _k^{n + 1} = \dfrac{{{{ \displaystyle \int_0^\infty {\omega \left| {\hat i_k^{n + 1}(\omega )} \right|} }^2}{\rm{d}}\omega }}{{{{ \displaystyle \int_0^\infty {\left| {\hat i_k^{n + 1}(\omega )} \right|} }^2}{\rm{d}}\omega }} $$ (7)
    $$ {\hat \lambda ^{n + 1}}(\omega ) = {\hat \lambda ^n}(\omega ) + \tau ({\hat i_{\rm{f}}}(\omega ) - \sum\limits_{k = 1}^K {\hat i_k^{n + 1}(\omega )} ) $$ (8)

    式中:$\hat i_k^{n + 1}(\omega )$为第k个模态分量ik第(n+1)次迭代值的傅里叶变换;${\hat i_{\rm{f}}}(\omega )$if的傅里叶变换;$\hat i_w^n(\omega )$为第w个模态分量iwn次迭代值的傅里叶变换;${\hat \lambda ^n}(\omega )$λn次迭代值的傅里叶变换;$\omega _k^{n + 1}$ωk第(n+1)次迭代值;τ为噪声容限参数。

    VMD算法的主要求解步骤如下:

    1) 初始化第k个模态分量ik的第1次迭代值{$i_k^1$}、{$\omega _k^1$}、λ第1次迭代值${\lambda ^1}$和最大迭代次数N,并令n=1。

    2) 更新ik、ωkλ,并令n=n+1。

    3) 若满足式(9)收敛条件,则迭代完成,输出最终的ikωk,若不满足收敛条件,则继续迭代。

    $$ \sum\limits_{k = 1}^K {\frac{{\left\| {\hat i_k^{n + 1}(\omega ) - \hat i_k^n(\omega )} \right\|_2^2}}{{\left\| {\hat i_k^n(\omega )} \right\|_2^2}}} < \varepsilon $$ (9)

    式中ε为迭代结束阈值,一般设为非常小的正数。

    在基于VMD的小电流接地故障选线方法中,通常取模态分量个数K=3,并利用IMF2分量进行故障选线[18-19]。当故障线路中故障零序电流IMF2分量的故障特征不明显时,则不能实现正确的故障选线。为了选择合适的模态分量,针对不同类型的单相接地故障,对VMD进行分析。

    在煤矿配电网模型中,针对以下2种故障类型进行分析:① 故障1发生在第4条线路S4的f1点,接地故障电阻Rg=5 Ω,故障初相角$\varphi $=90°。② 故障2发生在第4条线路S4的f1点,接地故障电阻Rg=2 000 Ω,故障初相角$\varphi $=90°。

    对于故障1,利用VMD对线路S2,S3,S4的故障零序电流if2if3if4进行分解,模态分量个数K=3和K=4的VMD结果如图3所示。图中,ifm为第m条线路Sm的故障零序电流,m=2,3,4;ikm为线路Sm故障零序电流的第k个模态分量。可看出:当K=3时,故障线路S4中故障零序电流IMF2分量(i24)的故障特征不明显,无法实现正确的故障选线;当K=4时,故障线路S4中故障零序电流IMF2分量的幅值比健全线路S2和S3大,且极性相反,故障特征明显,可以实现正确的故障选线。因此,VMD的模态分量个数K应取4。

    图  3  不同模态分量个数下的VMD结果
    Figure  3.  Variational pecomposition (VMD) decomposition results under different mode component number

    对于故障2,利用VMD对线路S2,S3,S4的故障零序电流if2if3if4进行分解,取K=4,发生高阻接地故障时VMD结果如图4所示。可看出:故障线路S4中故障零序电流IMF2分量(i24)的故障特征不明显,而IMF3分量(i34)的故障特征明显,应提取模态分量IMF3进行故障选线。因此,采用固定的模态分量,不能实现高准确度的故障选线。

    图  4  高阻接地故障时VMD结果
    Figure  4.  VMD decomposition results in case of high resistance grounding fault

    由上述分析可知,在利用VMD进行故障选线时,零序电流模态分量个数K=4,并根据模态分量的故障特征,在IMF2和IMF3中选择合适的模态分量进行故障选线。

    由于不同接地故障类型产生的故障零序电流模态分量IMFk的波形差别较大,幅值计算较为困难,因此,利用故障初始半个周波IMFk分量的暂态能量表征其幅值。

    $$ {E_{km}} = \sum\limits_{h = 1}^H {i_{km}^2\left( h \right)} $$ (10)

    式中:Ekm为母线第m条出线故障零序电流分量IMFk的暂态能量,k=2或k=3;h为采样点序号;H为半个周波采样点个数;ikm(h)为第m条出线故障零序电流IMFk分量的第h个采样值。

    不同线路故障零序电流IMFk分量的暂态能量在母线所有出线总能量中的占比不同。

    $$ {\varDelta _{km}} = \frac{{{E_{km}}}}{{{E_k}}} = \frac{{{E_{km}}}}{{ \displaystyle \sum\limits_{m = 1}^M {{E_{km}}} }} $$ (11)

    式中:Δkm为第m条出线故障零序电流IMFk分量的能量占比;Ek为母线所有出线故障零序电流IMFk分量的能量和;M为母线的出线数量。

    若第m条出线故障零序电流IMFk分量的暂态能量占比满足ΔkmΔthΔth为判据阈值),则第m条出线为故障线路,否则为非故障线路。

    受线路电气参数不同和线路分支的影响,利用暂态能量进行选线时,有时无法确定故障线路。此时,通过比较各出线半个周波故障零序电流IMFk分量的波形相似度进行故障选线。

    $$ {\rho _{kmq}} = \frac{{\displaystyle \sum\limits_{h = 1}^H {{i_{km}}(h){i_{kq}}(h)} }}{{\sqrt {\displaystyle \sum\limits_{h = 1}^H {{i}_{km}^2(h){i}_{kq}^2(h)} } }} $$ (12)

    式中:ρkmq为第m条和第q条出线故障零序电流IMFk分量的波形相似度;和ikq(h)为第q条出线故障零序电流IMFk分量的第h个采样值。

    为了实现故障选线,定义线路故障零序电流IMFk分量的波形相似度为

    $$ {\overline{\rho }}_{km}= \left\{ \begin{array}{l} \underset{q=[1,M],q\ne m}{\overset{}{\mathrm{max}}}{\rho }_{kmq} \qquad \exists {\rho }_{kmq} > 0 \\ \underset{q=[1,M],q\ne m}{\overset{}{\mathrm{min}}}{\rho }_{kmq} \qquad 其他 \end{array} \right. $$ (13)

    不同出线故障零序电流IMFk分量的波形相似度在母线所有出线总波形相似度中的占比不同。

    $$ {\eta _{km}} = \frac{{{{\bar \rho }_{km}}}}{{{\rho _k}}} $$ (14)
    $$ {\rho _k} = \sum\limits_{m = 1}^{M - 1} {{{\bar \rho }_{km}}} $$ (15)

    式中:ηkm为第m条出线故障零序电流IMFk分量的波形相似度占比;ρk为母线所有出线故障零序电流IMFk分量相似度总和。

    若第m条出线故障零序电流IMFk分量的波形相似度占比满足ηkmηthηth为判据阈值),则出线m为故障线路,否则为非故障线路。

    在不同配电网单相接地故障场景下,故障零序电流IMF2分量或IMF3分量可能存在故障特征不明显的情况;在同一个模态分量IMFkk=2,3)中,可能存在暂态能量或波形相似度故障特征不明显的情况。因此,基于单一模态分量和单一故障特征的故障选线方法的选线准确度较低。为了提高故障选线准确度,自动选取模态分量和故障特征进行故障选线,形成基于VMD的小电流接地故障融合选线方法。

    $$ \left\{ {\begin{array}{*{20}{l}} {{\varDelta _{km}} > \varDelta {}_{{\text{th}}{\text{.max}}}} \\ {{\eta _{km}} < {\eta _{{\text{th}}{\text{.max}}}}} \end{array}} \right. $$ (16)
    $$ \left\{ {\begin{array}{*{20}{l}} {{\varDelta _{km}} > \varDelta {}_{{\text{th}}{\text{.min}}}} \\ {{\eta _{km}} < {\eta _{{\text{th}}{\text{.min}}}}} \end{array}} \right. $$ (17)
    $$ \left\{ {\begin{array}{*{20}{l}} {\varDelta {}_{{\text{th}}{\text{.mid}}} > {\varDelta _{km}} > \varDelta {}_{{\text{th}}{\text{.max}}}} \\ {{\eta _{{\text{th}}{\text{.min}}}} < {\eta _{km}} < {\eta _{{\text{th}}{\text{.mid}}}}} \end{array}} \right. $$ (18)

    式中:Δth.maxΔth.minΔth.mid分别为暂态能量判据的最大阈值、最小阈值和中间阈值,通常取0.7<Δth.max<0.9,0.1<Δth.min<0.3,0.4<Δth.mid<0.6;ηth.maxηth.minηth.mid分别为波形相似度判据的最大阈值、最小阈值和中间阈值,通常取0.2<ηth.max<0.4,−0.6<ηth.min <−0.4,−0.1<ηth.mid<0.1。

    故障融合选线判据:当出线满足式(16)、式(17)、式(18)任意一个条件时,该出线为故障线路。当满足故障选线启动判据,但不满足式(16)、式(17)和式(18)时,则判定为母线故障。

    故障选线系统实时检测零序电压并判断是否发生单相接地故障,若发生单相接地故障,则执行故障融合选线算法,流程如图5所示。

    图  5  故障融合选线流程
    Figure  5.  Fault fusion line selection process

    1) 同步采集各出线故障初始半个周波的故障零序电流信号。

    2) 对故障零序电流进行VMD并提取IMF2和IMF3分量, 取VMD模态分量编号k=2,出线序号m=1。

    3) 计算第m条出线故障零序电流IMFk分量的暂态能量Ekm及能量占比Δkm

    4) 计算第m条出线故障零序电流IMFk分量的波形相似度ρkmq及波形相似度占比ηkm

    5) 根据式(16)、式(17)及式(18)判定第m条出线是否满足融合故障选线判据,若是,则该出线为故障线路,否则为健全线路。

    6) 循环判断其他出线是否为故障线路。若所有出线均为健全线路,则判断为母线故障。

    以煤矿配电网(图2)为例,利用电磁暂态仿真软件ATP/EMTP对单相接地故障进行仿真分析。故障发生在第4条线路S4的f1点,接地电阻Rg分别为0.5,10,1 000,2 000 Ω,故障初相角$\varphi $分别为0,45,90°。在上述各种故障类型下出线S2,S3,S4的暂态能量和波形相似度故障特征见表1。融合故障选线判据的阈值设置为Δth.max=0.8,Δth.min=0.2,Δth.mid=0.5,ηth.max=0.3,ηth.min=−0.5,ηth.mid=0,则利用公式(16)、(17)、(18)进行融合故障选线,结果见表1。可看出大部分故障选线结果正确,但在故障初相角为90°,接地电阻分别为1 000,2 000 Ω时,不能判断出故障线路。这是因为发生这2种故障时,故障线路S4故障零序电流IMF2分量的故障特征均不明显。对这2种故障选用IMF3分量进行融合故障选线,选线结果见表2。可看出故障选线结果均正确。因此,对于上述12种单相接地故障,所提出故障选线方法的选线正确率为100%。

    表  1  提取IMF2分量时的故障选线数据
    Table  1.  Fault line selection data when extracting IMF2 component
    Rgφ /(°)S2S3S4故障
    线路
    Δ22η22Δ23η23Δ24η24
    0.500.110 60.999 80.105 60.999 80.783 8−0.999 7S4
    450.113 21.064 10.114 11.064 10.772 6−1.128 2S4
    900.110 21.058 00.116 21.058 00.773 6−1.115 9S4
    1000.000 50.504 40.000 50.504 40.999 0−0.008 8S4
    450.150 21.571 60.132 31.571 60.717 5−2.143 2S4
    900.151 31.438 50.133 91.438 50.714 8−1.876 9S4
    1 00000.004 60.499 30.003 80.499 30.991 50.001 4S4
    450.000 5−1.00 710.000 00.003 50.999 50.003 5S4
    900.002 40.094 30.000 20.452 80.997 40.452 8
    2 00000.004 70.498 00.003 90.498 00.991 40.004 0S4
    450.000 4−1.002 30.000 00.001 10.999 50.001 1S4
    900.002 20.099 00.000 20.450 50.997 60.450 5
    下载: 导出CSV 
    | 显示表格
    表  2  提取IMF3分量时的故障选线数据
    Table  2.  Fault line selection data when extracting IMF3 component
    Rgφ/(°)S2S3S4故障
    线路
    Δ32η32Δ33η33Δ34η34
    1 000900.147 41.222 60.104 51.222 60.748 1−1.445 2S4
    2 000900.145 81.183 70.103 71.183 70.750 5−1.367 3S4
    下载: 导出CSV 
    | 显示表格

    为了验证故障融合选线方法对不同位置故障的适应性,在煤矿配电网第4条线路的f1点、第7条线路的f2点、第10条线路的f3点分别发生单相接地故障,接地电阻Rg=500 Ω,故障初相角φ=90°,选线结果见表3。可看出故障选线结果均正确,说明故障融合选线方法不受故障位置的影响。

    表  3  不同故障位置的故障选线数据
    Table  3.  Fault line selection data at different fault locations
    故障位置S2S3S4故障
    线路
    Δ32η32Δ33η33Δ34η34
    f10.313 90.020 70.060 60.020 70.625 4−1.041 4S4
    f20.132 61.065 90.067 01.065 90.800 4−1.131 8S7
    f30.143 40.998 40.067 60.998 40.789 0−0.996 9S10
    下载: 导出CSV 
    | 显示表格

    为了证明所提故障融合选线方法的优越性,将故障融合选线方法与能量法[18]和相关性聚类法[19]进行比较。利用能量法和相关性聚类法对上述12种故障进行选线,结果见表4。可看出能量法的故障选线正确率为83%,相关性聚类法的故障选线正确率为50%,故障融合选线方法的选线正确率较能量法和相关性聚类法分别提高了17%和50%。

    表  4  不同故障选线方法的选线结果
    Table  4.  Line selection results of different fault line selection methods
    Rgφ/(°)能量法相关性聚类法融合选线法
    故障线路选线结果故障线路选线结果故障线路选线结果
    0.50S4正确S4正确S4正确
    45S4正确S4正确S4正确
    90S4正确S4正确S4正确
    100S4正确S4正确S4正确
    45母线错误S4正确S4正确
    90母线错误S4正确S4正确
    1 0000S4正确母线错误S4正确
    45S4正确S2错误S4正确
    90S4正确母线错误S4正确
    2 0000S4正确母线错误S4正确
    45S4正确S2错误S4正确
    90S4正确母线错误S4正确
    下载: 导出CSV 
    | 显示表格

    为了进一步验证故障融合选线方法的有效性,在煤矿配电网模型中对多种单相接地故障类型进行仿真,接地故障电阻分别为0.5,5,10,20,50,100,200,500,1 000,2 000 Ω,故障初相角分别为0,30,45,60,90,135°,故障位置分别在10条线路S2−S11上,共进行600次仿真实验。利用所有仿真数据对基于VMD的融合暂态能量和波形相似度比较的故障选线方法进行验证,结果为故障选线正确率100%,表明该方法具有较高的故障选线准确度,可应用于故障选线装置。

    1) 在利用VMD提取煤矿配电网单相接地故障零序电流的暂态特征时,VMD的模态分量个数应为4,根据不同的故障类型选择模态分量IMF2或IMF3进行煤矿配电网故障选线。

    2) 用电磁暂态仿真软件ATP/EMTP搭建了煤矿配电网模型,在不同接地故障电阻、故障初相角和故障位置的单相接地故障场景下,对所提出的故障融合选线方法进行验证。结果表明:在配电网发生各种单相接地故障时,故障融合选线方法的选线正确率为100%,且不受故障位置的影响;故障融合选线方法具有较高的故障选线准确度,较能量法和相关性聚类法的故障选线正确率分别提高了17%和50%,可应用于故障选线装置。

  • 图  1   单相接地故障零序等效电路

    Figure  1.   Zero sequence equivalent circuit of single-phase ground fault

    图  2   煤矿配电网模型

    Figure  2.   Coal mine distribution network model

    图  3   不同模态分量个数下的VMD结果

    Figure  3.   Variational pecomposition (VMD) decomposition results under different mode component number

    图  4   高阻接地故障时VMD结果

    Figure  4.   VMD decomposition results in case of high resistance grounding fault

    图  5   故障融合选线流程

    Figure  5.   Fault fusion line selection process

    表  1   提取IMF2分量时的故障选线数据

    Table  1   Fault line selection data when extracting IMF2 component

    Rgφ /(°)S2S3S4故障
    线路
    Δ22η22Δ23η23Δ24η24
    0.500.110 60.999 80.105 60.999 80.783 8−0.999 7S4
    450.113 21.064 10.114 11.064 10.772 6−1.128 2S4
    900.110 21.058 00.116 21.058 00.773 6−1.115 9S4
    1000.000 50.504 40.000 50.504 40.999 0−0.008 8S4
    450.150 21.571 60.132 31.571 60.717 5−2.143 2S4
    900.151 31.438 50.133 91.438 50.714 8−1.876 9S4
    1 00000.004 60.499 30.003 80.499 30.991 50.001 4S4
    450.000 5−1.00 710.000 00.003 50.999 50.003 5S4
    900.002 40.094 30.000 20.452 80.997 40.452 8
    2 00000.004 70.498 00.003 90.498 00.991 40.004 0S4
    450.000 4−1.002 30.000 00.001 10.999 50.001 1S4
    900.002 20.099 00.000 20.450 50.997 60.450 5
    下载: 导出CSV

    表  2   提取IMF3分量时的故障选线数据

    Table  2   Fault line selection data when extracting IMF3 component

    Rgφ/(°)S2S3S4故障
    线路
    Δ32η32Δ33η33Δ34η34
    1 000900.147 41.222 60.104 51.222 60.748 1−1.445 2S4
    2 000900.145 81.183 70.103 71.183 70.750 5−1.367 3S4
    下载: 导出CSV

    表  3   不同故障位置的故障选线数据

    Table  3   Fault line selection data at different fault locations

    故障位置S2S3S4故障
    线路
    Δ32η32Δ33η33Δ34η34
    f10.313 90.020 70.060 60.020 70.625 4−1.041 4S4
    f20.132 61.065 90.067 01.065 90.800 4−1.131 8S7
    f30.143 40.998 40.067 60.998 40.789 0−0.996 9S10
    下载: 导出CSV

    表  4   不同故障选线方法的选线结果

    Table  4   Line selection results of different fault line selection methods

    Rgφ/(°)能量法相关性聚类法融合选线法
    故障线路选线结果故障线路选线结果故障线路选线结果
    0.50S4正确S4正确S4正确
    45S4正确S4正确S4正确
    90S4正确S4正确S4正确
    100S4正确S4正确S4正确
    45母线错误S4正确S4正确
    90母线错误S4正确S4正确
    1 0000S4正确母线错误S4正确
    45S4正确S2错误S4正确
    90S4正确母线错误S4正确
    2 0000S4正确母线错误S4正确
    45S4正确S2错误S4正确
    90S4正确母线错误S4正确
    下载: 导出CSV
  • [1] 罗超,耿蒲龙,曲兵妮,等. 基于小波包的矿井供电系统单相接地故障选线方法[J]. 工矿自动化,2018,44(2):68-74.

    LUO Chao,GENG Pulong,QU Bingni,et al. A line selection method for single phase ground fault in coal mine power supply system based on wavelet packet[J]. Industry and Mine Automation,2018,44(2):68-74.

    [2] 高宏杰,赵建文,郭秀才. 煤矿电网单相漏电故障区段自动定位探索[J]. 工矿自动化,2021,47(5):106-111.

    GAO Hongjie,ZHAO Jianwen,GUO Xiucai. Research on automatic location of single-phase leakage fault zone in coal mine power network[J]. Industry and Mine Automation,2021,47(5):106-111.

    [3] 董礼清,薛永端,李娟,等. 小电流接地故障零序电压分布特征与测距[J]. 电力系统自动化,2022,46(19):180-187.

    DONG Liqing,XUE Yongduan,LI Juan,et al. Zero-sequence voltage distribution characteristics and locatiaon of non-solidly grounding fault[J]. Automation of Electric Power Systems,2022,46(19):180-187.

    [4] 刘扬,刘建功,王毅颖,等. 煤矿坚强智能电网建设理论与技术探讨[J]. 煤炭学报,2020,45(6):2296-2307.

    LIU Yang,LIU Jiangong,WANG Yiying,et al. Discussion on theory and technology of building robust intelligent power grid in coal mine of China[J]. Journal of China Coal Society,2020,45(6):2296-2307.

    [5] 程路,陈乔夫. 小电流接地系统单相接地选线技术综述[J]. 电网技术,2009,33(18):219-224.

    CHENG Lu,CHEN Qiaofu. A survey on faulty line selection technology for single-phase grounded transmission line in small current neutral grounded system[J]. Power System Technology,2009,33(18):219-224.

    [6] 周宇,汤涛,曾祥君,等. 基于零序电流幅值比倍增系数的灵活接地系统故障选线方法[J]. 电力系统保护与控制,2022,50(23):112-120.

    ZHOU Yu,TANG Tao,ZENG Xiangjun,et al. Fault line detection method for a flexible grounding system based on a zero-sequence current amplitude ratio multiplication coefficient[J]. Power System Protection and Control,2022,50(23):112-120.

    [7] 方毅,薛永端,宋华茂,等. 谐振接地系统高阻接地故障暂态能量分析与选线[J]. 中国电机工程学报,2018,38(19):5636-5645,5921.

    FANG Yi,XUE Yongduan,SONG Huamao,et al. Transient energy analysis and faulty feeder identification method of high impedance fault in the resonant grounding system[J]. Proceedings of the CSEE,2018,38(19):5636-5645,5921.

    [8] 季鹏,陈芳芳,徐天奇,等. 基于Spearman相关系数法与有功分量法的高阻接地故障选线方法研究[J]. 山东电力技术,2022,49(12):8-13,31. DOI: 10.3969/j.issn.1007-9904.2022.12.002

    JI Peng,CHEN Fangfang,XU Tianqi,et al. Research on high-resistance grounding fault line selection method based on Spearman correlation coefficient and active component[J]. Shandong Electric Power,2022,49(12):8-13,31. DOI: 10.3969/j.issn.1007-9904.2022.12.002

    [9] 韦明杰,石访,张恒旭,等. 基于同步零序电流谐波群体比相的谐振接地系统高阻故障选线及区段定位方法[J]. 中国电机工程学报,2021,41(24):8358-8372.

    WEI Mingjie,SHI Fang,ZHANG Hengxu,et al. Feeder selection and section location of high impedance fault at resonant networks based on the phase differences between the synchronous harmonics of the zero sequence currents[J]. Proceedings of the CSEE,2021,41(24):8358-8372.

    [10] 樊淑娴,徐丙垠,张清周. 注入方波信号的经消弧线圈接地系统故障选线方法[J]. 电力系统自动化,2012,36(4):91-95.

    FAN Shuxian,XU Bingyin,ZHANG Qingzhou. A new method for fault line selection in distribution on system with arc suppression coil grounding with square-wave signal injection[J]. Automation of Electric Power Systems,2012,36(4):91-95.

    [11]

    SHI Shenxing,ZHU Beier,LEI Aoyu,et al. Fault location for radial distribution network via topology and reclosure-generating traveling waves[J]. IEEE Transactions on Smart Grid,2019,10(6):6404-6413. DOI: 10.1109/TSG.2019.2904210

    [12] 李聪聪,李玉敦,杨超,等. 基于序分量电压加权的双端故障定位算法研究[J]. 山东电力技术,2021,48(7):53-57,69.

    LI Congcong,LI Yudun,YANG Chao,et al. Research on two-terminal fault location algorithm based on weighted voltage sequence component[J]. Shandong Electric Power,2021,48(7):53-57,69.

    [13] 高锋阳,李昭君,袁成,等. 量子计算和免疫优化算法相结合的有源配电网故障定位[J]. 高电压技术,2021,47(2):396-406.

    GAO Fengyang,LI Zhaojun,YUAN Cheng,et al. Fault location for active distribution network based on quantum computing and immune optimization algorithm[J]. High Voltage Engineering,2021,47(2):396-406.

    [14]

    ZHAGN Chenyu,YUAN Xiaodong,SHI Mingming,et al. Fault location method based on SVM and similarity model matching[J]. Mathematical Problems in Engineering,2020(1):1-9.

    [15] 李科,随晓娜,张俊,等. 矿井电网故障选线方法研究[J]. 工矿自动化,2018,44(5):70-75.

    LI Ke,SUI Xiaona,ZHANG Jun,et al. Research on fault line selection method of mine power network[J]. Industry and Mine Automation,2018,44(5):70-75.

    [16] 原磊明,关瑞,石国栋,等. 基于EEMD与小波阈值去噪的单相接地故障选线方法[J]. 煤炭工程,2020,52(12):137-140.

    YUAN Leiming,GUAN Rui,SHI Guodong,et al. Line selection method for one-phase grounding fault based on EEMD and wavelet threshold denoising[J]. Coal Engineering,2020,52(12):137-140.

    [17] 毕胜,耿蒲龙,张建花,等. 基于CEEMD与自相关阈值去噪的单相接地故障选线方法研究[J]. 煤炭工程,2022,54(7):153-158.

    BI Sheng,GENG Pulong,ZHANG Jianhua,et al. Line selection method for single-phase ground fault based on CEEMD and autocorrelation threshold denoising[J]. Coal Engineering,2022,54(7):153-158.

    [18] 翟二杰,舒征宇,汪俊,等. 基于VMD−LSTM的小电流接地系统故障选线方法[J]. 电工电能新技术,2021,40(1):70-80.

    ZHAI Erjie,SHU Zhengyu,WANG Jun,et al. Fault line selection method of small current grounding system based on VMD-LSTM[J]. Advanced Technology of Electrical Engineering and Energy,2021,40(1):70-80.

    [19] 魏科文,张靖,何宇,等. 基于VMD和相关性聚类的谐振接地系统单相接地故障选线[J]. 电力系统保护与控制,2021,49(22):105-113.

    WEI Kewen,ZHANG Jing,HE Yu,et al. Single-phase grounding fault line selection in a resonant grounding system based on VMD and correlation clustering[J]. Power System Protection and Control,2021,49(22):105-113.

    [20] 薛永端,李娟,徐丙垠. 中性点经消弧线圈接地系统小电流接地故障暂态等值电路及暂态分析[J]. 中国电机工程学报,2015,35(22):5703-5714.

    XUE Yongduan,LI Juan,XU Bingyin. Transient equivalent circuit and transient analysis of single-phase earth fault in arc suppression coil grounded system[J]. Proceedings of the CSEE,2015,35(22):5703-5714.

    [21] 庄伟,牟龙华,童荣斌. 基于罗氏线圈二次信号的煤矿高压电网接地故障区段定位[J]. 煤炭学报,2014,39(6):1184-1190.

    ZHUANG Wei,MU Longhua,TONG Rongbin. Secondary signal of rogowski coil based ground fault section location in coalmine high-voltage distribution network[J]. Journal of China Coal Society,2014,39(6):1184-1190.

  • 期刊类型引用(3)

    1. 杨战社,张程,荣相. 矿井供电系统单相接地故障选线方法现状与发展趋势. 煤矿安全. 2025(01): 212-219 . 百度学术
    2. 程卫健,史小军,穆艳祥. 煤矿电网的单相接地故障选线研究. 煤矿机械. 2024(09): 160-163 . 百度学术
    3. 庄原,朱艺超,贺飞,汪若涵. 基于遗传算法的配电网小电流接地故障诊断研究. 电子元器件与信息技术. 2023(06): 69-71 . 百度学术

    其他类型引用(3)

图(5)  /  表(4)
计量
  • 文章访问数:  256
  • HTML全文浏览量:  60
  • PDF下载量:  22
  • 被引次数: 6
出版历程
  • 收稿日期:  2022-08-28
  • 修回日期:  2023-03-01
  • 网络出版日期:  2022-10-20
  • 刊出日期:  2023-03-24

目录

/

返回文章
返回