一种矿井漏电保护选线方法

于群，尚雪丽

(山东科技大学 电气与自动化工程学院， 山东 青岛 266590)

0 引言

矿井供电系统中单相接地故障约占矿井电气故障的80%，若不及时检测并切除故障线路，可能会造成人身触电，甚至引起瓦斯爆炸事故[1]，对作业人员生命安全造成严重威胁。因此，要求漏电保护装置必须灵敏可靠，快速准确切除故障线路[2]。

目前，漏电保护方法分为基于稳态量的选线方法与基于暂态量的选线方法，主要有零序电压型保护法、零序功率方向保护法、五次谐波保护法、首半波法、附加直流检测式保护法等[3-4]。由于矿井条件限制及强电磁干扰等问题，上述方法在实际保护中的效果并不十分可靠，装置拒动、误动现象时有发生。文献[5]利用零序功率方向保护法对故障线路进行判别，但该方法受中性点接地方式的影响，当采用消弧线圈接地运行方式时，会导致误判；文献[6]利用小波算法获取特定频段的暂态分量，以故障线路零序电流小波系数对应的模极大值与波形相似趋势为判据进行故障选线，但当线路中出现强烈的干扰信号时，该方法不可靠，且不同的信号需要选择不同的小波函数，选择不当可能会导致选线错误。文献[7]提出了一种以零序电流波形相似性为判据的选线方法，采用互近似熵来反映零序电流波形之间的相似度，该方法具有较好的抗干扰能力，但当系统发生母线故障时，该方法不能有效识别。

本文提出了一种基于零序电流积分与改进Bhattacharyya距离算法的矿井漏电保护选线方法。该方法利用零序电流积分递推公式对零序电流波形进行预处理，减小了突变数据的影响；通过改进的Bhattacharyya距离算法求取各线路零序电流积分序列之间的Bhattacharyya系数，以Bhattacharyya系数累加和求得的综合累加系数作为选线判据，并以同一时刻各线路电流采样值的乘积和作为辅助判据对母线故障进行判别，可提高选线的准确性。

1 单相接地故障零序电流波形相似性分析

图1 单相接地故障零序网络Fig.1 Zero-sequence network of single-phase grounding fault

(1)

(2)

由于中性点经消弧线圈接地系统常处于过补偿运行状态，消弧线圈感性电流补偿容性电流，降低了故障线路零序电流幅值，且故障线路与非故障线路的零序电流方向相同，仅依靠零序电流稳态分量作为判据不再有效，因此，有必要研究故障后暂态零序电流的变化情况。为深入分析中性点经消弧线圈接地系统单相接地故障后的零序电流特征，将图1转换为单相接地故障暂态等效电路，如图2所示。u0，R0，L0分别为零序回路中等值电压、电阻、电感；C0∑为对地电容总和；i0L为暂态电感电流；i0C为暂态电容电流。

图2 单相接地故障暂态等效电路Fig.2 Transient equivalent circuit of single-phase grounding fault

u0=Umsin(ωt+φs)

(3)

式中：Um为相电压幅值；ω为角频率；t为时间；φs为故障时相电压相角。

(4)

故障线路与非故障线路流过的暂态零序电流主要由暂态电容电流i0C决定，由式(3)、式(4)可求出i0C。

由于消弧线圈的作用，故障线路直流分量衰减较快且振荡频率高，其高频分量为全部非故障线路高频分量之和。因此可通过线路之间零序电流波形的差异判断故障线路。

2 选线方法

2.1 零序电流积分递推

为表征零序电流在数据窗内的总体分布趋势，采用连续时间内零序电流积分递推的处理方法，并通过半周积分算法[8]来定义零序电流积分递推公式。

(5)

则电流有效值为

(6)

通过矩形面积求和近似代替积分，可得

(7)

式中：Ts为间隔时间；N为1个工频周期内采样点数；i(n)为第n个采样点零序电流。

(8)

正弦波对任意Δt(Δt

(9)

式中：I(x)为电流积分结果；Nt为Δt时间内采样点数。

考虑到故障选线对速度的要求，本文将数据窗缩短，取Nt为3，即对连续3个采样数对应时间间隔的零序电流求积分平均值得出零序电流积分序列。

2.2 改进Bhattacharyya距离算法

在统计学中，Bhattacharyya距离算法用于测量2个离散概率分布的相似性[9]。通过对2个离散概率分布的重叠量进行近似计算，可分析两者之间的相似程度[10-11]。

在定义域X中，定义2个离散概率分布函数p(x)和q(x)的Bhattacharyya距离为

DB(p,q)=-lnBC(p,q)

(10)

对2条线路零序电流波形的采样值进行处理，得到2个离散概率分布函数p(x)和q(x)，然后利用Bhattacharyya距离算法计算零序电流相似度。相同区间内，直方图重叠程度越高，表明2条线路零序电流波形的相似度越高；直方图重叠程度越低，表明2条线路零序电流波形的相似度越低；直方图完全没有重叠时，表明2条线路零序电流波形在该区间具有很大差异，呈不相关性。

原始的Bhattacharyya距离算法存在以下问题：① 没有考虑2个离散概率分布函数在不同区间所占的权重比；② 当2个离散概率分布函数在相同区间的频数同时为0时，直方图高度重叠，而原始的Bhattacharyya距离算法忽略了该种特殊情况，认为p(x)q(x)=0，导致零序电流波形的相似度提高，影响选线准确性。

针对上述问题，引入权重wh对Bhattacharyya距离系数进行修正，以放大离散概率分布函数在不同区间的差异。若2个离散概率分布函数在同一区间的频数均为0，令p(x)q(x)=1。

(11)

(12)

式中Cph和Cqh分别为p(x)和q(x)在区间的频数。

2.3 选线判据

为更明显地判断出故障线路，将第i条线路的综合累加系数Pi作为选线判据。

(13)

对比各线路的综合累加系数并选出最小值P，P对应的线路即为故障线路。

P=min(P1,P2,…,Pm)

(14)

当母线发生故障时，系统中各线路零序电流波形相似程度变大，求得的综合累加系数差值变小，仅采用综合累加系数作为选线判据无法有效判断母线故障，因此利用同一时刻各线路电流采样值的乘积和Ssgn作为辅助判据来区分线路故障和母线故障。

(15)

式中：Ns为1/4工频周期内采样点数；ii(n)为第i条线路第n个采样点零序电流。

若Ssgn=-1，表示系统发生线路故障；若Ssgn=1，表示系统发生母线故障。

矿井电网发生单相接地故障时，选线流程如图3所示。

图3 单相接地故障选线流程Fig.3 Line selection process of single-phase grounding fault

3 仿真验证

3.1 算例分析

利用Matlab/Simulink搭建单相接地故障仿真模型[12-13]，如图4所示，其中母线电压为10.5 kV，具体线路参数见表1。

图4 单相接地故障仿真模型Fig.4 Simulation model of single-phase grounding fault

表1 仿真线路参数Table 1 Simulation line parameters

以中性点不接地系统为例，设置采样频率为4 kHz，仿真时间为0.3 s，线路l1发生单相接地故障，故障点距离母线2 km，故障初相角为0°，接地电阻为2 kΩ。通过对故障发生后各线路首个1/4工频周期的零序电流数据进行积分递推处理，得到各线路零序电流积分序列，如图5所示。

图5 零序电流积分序列Fig.5 Zero-sequence current integral sequence

提取各线路零序电流积分序列在数据窗内的最小值a与最大值c，将区间[a，c]分为5个子区间并依次标号为1—5，建立零序电流积分序列在各子区间的分布直方图，如图6所示。

图6 零序电流积分序列分布直方图Fig.6 Distribution histogram of zero-sequence current integral sequence

通过零序电流积分序列分布直方图得出概率分布函数，根据式(12)计算各线路零序电流积分序列之间的Bhattacharyya系数：

BC(li,lj)=

(16)

各线路综合累加系数：P1=0.78，P2=14.37，P3=13.95，P4=13.71，P5=10.37，P6=14.37。辅助判据：Ssgn=-1。由此可得线路l1为故障线路，仿真结果正确。

3.2 方法适应性分析

3.2.1 不同中性点接地方式

分别在中性点不接地方式与中性点经消弧线圈接地方式下设置线路l2发生单相接地故障，故障点距离母线1 km，故障初相角为90°，接地电阻为500 Ω。故障选线结果见表2。

表2 不同中性点接地方式下故障选线结果Table 2 Fault line selection results under different neutral grounding modes

从表2可看出，不同中性点接地方式下，故障线路l2的综合累加系数始终最小，表明本文方法能够在不同中性点接地方式下正确选出故障线路，选线结果不受中性点接地方式的影响。

3.2.2 噪声影响

考虑矿井环境噪声干扰大，会对选线结果造成一定的影响，仿真测试在高斯白噪声干扰下本文方法的性能。以中性点经消弧线圈接地系统为例，设置线路l3为故障线路，故障点距离母线4 km，故障初相角为30°，在仿真过程中分别加入10，50 dB的高斯白噪声，在不同接地电阻条件下进行仿真。故障选线结果见表3。

表3 噪声干扰下故障选线结果Table 3 Fault line selection results under noise interference

从表3可看出，在噪声干扰下，故障线路l3的综合累加系数明显小于其他线路，且故障线路与非故障线路的综合累加系数差异明显，不同噪声干扰与不同接地电阻对综合累加系数的影响小，表明本文方法具有较强的抗噪声干扰能力。

3.2.3 不同故障工况

当矿井供电系统发生单相接地故障时，零序电流波形随着故障初相角、过渡电阻、故障距离的改变而改变[14-15]。以中性点不接地系统为例，设置不同故障工况进行仿真，故障选线结果见表4。

表4 不同故障工况下故障选线结果Table 4 Fault line selection results under different fault conditions

从表4可看出：当母线发生故障时，可通过辅助判据进行有效区分；当线路发生单相接地故障时，不同条件下各线路综合累加系数发生变动，但故障线路的综合累加系数始终最小；过渡电阻与故障距离对综合累加系数的影响不明显，故障初相角对综合累加系数的影响较大，当故障初相角大于90°时，各线路综合累加系数的差异减小，但仍能正确选出故障线路，表明本文方法在不同故障工况下选线的可靠性较高。

4 结语

基于矿井电网发生单相接地故障时故障线路与非故障线路零序电流波形存在差异的原理，提出了一种基于零序电流积分与改进Bhattacharyya距离算法的矿井漏电保护选线方法。通过对故障后首个1/4工频周期的零序电流数据进行积分递推处理，得出零序电流积分序列的分布趋势，有效避免了突变数据对选线准确率的影响；采用改进的Bhattacharyya距离算法使故障线路与非故障线路之间的故障特征指标差别明显增大，并以同一时刻各线路电流采样值的乘积和作为辅助判据对母线故障进行判别，从而准确区分线路故障和母线故障。仿真结果表明，该方法适用于不同中性点接地方式，具有较强的抗噪声干扰能力，受不同故障工况的影响小，选线准确率高。