哈尔滨理工大学学报
JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY
第26卷第2期
2021年4月
Vol. 26 No. 2
Apr. 2021
嵇文路- 赵晓龙S 张 明- 杨红磊S 翁嘉明'
(1.国网江苏省电力有限公司南京供电分公司，南京210019；
2.中国电力科学研究院有限公司，北京100192；
3.上海交通大学电力传输与功率变换控制教育部重点实验室，上海200240)
摘要：为提高配电网单相接地故障特征提取的有效性以及实现对故障选线的准确判别，提出 了一种基于小波包全频带分析和在线序贯极限学习机(OS-ELM)的新的故障选线分析方法。本文 灵活运用小波包的多分辨率分频特性，提取了小电流单相接地故障暂态信号的特征量；利用OS- ELM 将批处理和逐次迭代相结合，不 断更新训练数据和网 络结构，实现了对单相接地故障线路的 快速选取。仿真实验结果表明：利用小波包全频带分析能够获得更有效的故障特征,OS-ELM 在测
试时间和诊断精度上都优于该反向传播(BP)网络、径向基网络(RBF)和支持向量机(SVM)，说明
了所提方法在单相接地故障选线的有效性，更能满足在线应用需求。
关键词：小波包全频带分析；在线序贯极限学习机；故障选线
DOI ：10. 15938/j. jhust. 2021. 02. 014
中图分类号：TM711 文献标志码：A 文章编号：1007-2683(2021 )02-0110-08
Fault Line Selection of Single Phase Grounding Based on Wavelet Packet Full Frequency A
nalysis and OS-ELM
JI Wen-lu , ZHAO Xiao-long , ZHANG Ming 1, YANG Hong-lei 2, WENG ]ia-ming
(1. Slate Grid Jiangsu Electric Power Co. , Ltd. Nanjing Power Supply Branch , Nanjing Jiangsu , 210019, China ；
2. China Electric Power Research Institute, Beijing 100192, China ；
3. Key Laboratory of Control of Power Transmission and Conversion , Ministry of Education , Shanghai Jiao Tong University , Shanghai 200240, China)
Abstract : In order to enhance effectiveness of single-phase ground fault feature extraction and to achieve exactly identification of fault line selection , a new fault line selection method of single phase grounding fault based on wavelet packet and online sequential extreme learning machine (OS-ELM ) is proposed ・ In this paper , the fault characteristics of fault transient signals are extracted by using the multi-resolution frequency division characteristics of wavelet packets ・ Combining batch and successive iterative , and constantly updating the training data and network structure , to achieve quickly fault line selection of single-phase gro
unding ・ The simulation results show that the more effective fault characteristics can be obtained using wavelet packet full frequency analysis , and 0S ・ELM is superior to the back propagation ( BP) network , radial basis function network ( RBF) and support vector machine (SVM ) in test time and diagnostic accuracy. The proposed method is effective in the single phase grounding fault line, can better meet the needs of online application ・
Keywords : wavelet packet full frequency analysis ； OS-ELM ； single-phase ground fault ； fault line selection
收稿日期：2019 -09 -24
基金项目：国家自然科学基金(51677116)；江苏省电力有限公司科技项目(J2018110-2).作者简介：嵇文路(1974—)，男，博士，高级工稈师；
翁嘉明(1986—)，男，博士，讲师.
通信作者：赵晓龙(1988—)，男，硕士，E ・mail ：314344808@qq. com.
第2期嵇文路等：基于小波包全频带分析和OS-ELM的小电流单相接地故障选线111
0引言
我国配电网具有点多面广的特点，配网故障也以单相接地故障为主E，然而在小电流接地系统发生单相接地故障时，由于中性点接地方式与大电网的不同，产生的短路电流信号微弱，不利于准确可靠地检测出故障线路，如何提高小电流单相接地故障选线的成功率，对配网可靠运行具有重要意义。
目前，国内外学者已开展的传统单相接地故障研究方法主要包括基于暂态零序电流突变信号、能量信号、幅值信号的方法等S3。文［5］用小波变换的方法提取故障后的行波信息，构造判据以实现故障线路的判定。文［6］利用小波变换的方法提取出各馈线零序电流的暂态量，通过比较电流极性实现故障选线。在配电网单相接地故障中，由于暂态信号受到过渡电阻等多种噪声因素的干扰，易造成误选的情况，严重影响了算法的准确性。由于小波变换具有对信号局部提取和滤除噪声信号的良好性能，因此更适应于暂态信号杂乱、故障状态复杂的配网单相接地故障分析⑺。
在配电网单相接地故障分析中，具有分布并行处理、非线性映射、自适应学习和鲁棒容错等特性的神经网络方法尤为适用。国内外学者利用不同的神经网络分析方法，开展了小电流单相接地故障选线分析的研究，并取得了显著进展⑻。传统的神经网络分析方法存在学习速度慢、泛化能力不强、易陷入局部最优等缺陷，文［9］提出一种新型单隐层前馈神经网络算法，称为极限学习机(extreme learning machine,ELM)，其具有学习速度快、泛化性能优良等特点，已被应用于电力系统优化分析求解当中。在ELM基础上，文［10］提出了一种在线序贯ELM (online sequential-ELM,OS-ELM)算法，其特点是历史数据分批进行训练并且支持可变化的样本数，每一轮训练过程中训练算法仅输入当前批次数据并
更新网络权值，无需重复扫描历史数据，使得该方法具有更强的泛化能力。考虑到实际配电网单相接地故障选线系统中所有的数据可能不是一次性添加到网络中，且实际配网数据存在错误或缺失，因此os-elm算法更能满足配电网单相故障接地选线的需要。基于极限学习机的众多优点，本文提出了一种小波包全频带分析和OS-ELM相结合的配电网单相接地故障选线新方法。2单相接地故障特征提取
2.1配电网单相接地故障暂态分析
谐振接地配电网发生单相接地故障时，非故障相电压升高，故障相电压降低，故障相电容通过故障线路向故障点放电，由于线路具有分布电容和电感参数，使得故障电流具有衰减震荡特性〔⑴。暂态电容电流是由故障相电压突然降低而引起的放电电流和非故障相电压突然升高而引起的充电电容电流之和。在故障点有两个暂态分量电流流过，其中暂态电容电流衰减较快，暂态电感电流衰减较慢。流过。配电网谐振接地系统模型如图1所示，由文［12］知 故障接地电流表达式为：
i d=(/Cm cos(a/+宰)+/(；,”［(ftysin^sintiX)/o)-cos(pcosa)jt］e_l/TC+/厶,”cos°e-‘仇 (1)式中:/c”，为电容电流稳态幅值；兮为暂态自由分量振荡分量的角频率；仇为电容回路的时间常数；匚为电感电流的稳态值；九为电感回路的时间常数。
图110kV谐振接地系统模型
Fig.1Model of a10kV resonant grounding system
故障暂态能量主要集中在故障零序电流的首半波，因此故障信息主要从故障后首半波获得。故障发生瞬间，零序电流的暂态特征明显，但持续时间较短。暂态电流大小在很大程度上取决于暂态电容电流的大小，频率大小亦是如此。其幅值跟初始相位角相关，故障出现在相电压处在峰值时刻时，暂态电容电流幅值达到最大值，而故障岀现在相电压处在过零位置时，暂态电容幅值很小，衰减直流分量幅值、能量较大。
手机包装
2.2基于小波包全频带分析的故障特征提取
在谐振接地配电网的单相接地故障中，其故障信号的暂态特征受故障发生时刻、接地点过度阻抗、线路参数等随机条件的影响很大。为反映故障条件的随机性及其对故障暂态信号特征的影响，利用小
112哈尔滨理工大学学报第26卷
自动抄表水表
波包良好的时频局部化能力与能量无损性特点，对故障暂态电流信号做小波包全频带分析。
小波包函数碣心)定义为
昧心)(2)式中：71=0,1,2,…为振荡参数J w Z和keZ分别是尺度参数和平移参数。
当n=0,1；j=k=0时，初始的2个小波包函数定义为
如0)=W(t),"1Q)=/(/)(3)式中呼Q)和&(/)分别是正交尺度函数和正交小波函数。初始小波函数满足以下双尺度方程：
=72X A(A;)/Z,o(2z-k)■
蛀z(4) Mi(0⑵一力).
KE.Z
式中山(巧和gQ)分别为对应的多尺度分析中低通滤波系数和高通滤波系数。则当“=2,3…时，其他的小波包函数满足：
如Q)=工力(％)“”⑵-％)
詔(5)如+i(')=血Yg仏)⑵一町.
KG.Z
那么，由式(5)所决定的函数集合1^(01,71=0,1, 2,…，就称为关于正交尺度函数卩(/)的小波包。对于一组离散信号%(0，小波包分解重构算法如式(6)和式(7)：
(6)
丐：；5)=*g(k-2t)d；(k)1
d；(k)=2[》h(k-2T)d；：J(k)+
弓g(—2T)硝心)] (7)式中:d'；仏)表示经小波包分解后节点(j,n)所对应的第力个系数，节点0，小表示第j层的第n个频带。
利用上述小波包过程进行单相接地故障特征提取，对各线路零序电流进行小波包分解。第i个频带对应的能量值为
车流量
E严愆烤仏)|2；21,2,...,2)(8)式中:^；)(n)为小波包分解第(_//)子频带下的系数，每个子频带下共有n个系数。
通过对暂态信号的故障特征的分析可知：
在电压相角＜p=0附近发生任意过渡电阻大小的接地故障和在任意电压相角大小发生过渡电阻较大的接地故障两种情况下，各线路的零序暂态信号能量集中频带都是在低频带上。因此可以统一在低频带上按照各线路暂态信号的能量大小来确定故障选线。在低频带上，故障线路所具有的能量要大于其他线路所具有的能量值。
在电压相角＜p=&2附近发生过渡电阻较小的单相接地故障情况下，各线路能量集中的频带可能不同，能量集中的频带也可能有多个，依特征频带选线较困难，但此时过渡电阻较小，暂态零序电流的幅值较大，谐波的五次分量相对也较大。因此，可以按照各线路暂态信号在五次谐波所在频带的能量大小来判别出故障线路。
因为故障接地线路的暂态零序电流信号的幅值为系统网络内的其他所有健全线路的暂态零序电流信号的幅值总和，零序暂态信号在各频带的信号能量表示为小波包系数的平方和，所以，综合考虑各线路在全频带上的能量统计和分布特性，在绝大多数子频带上满足故障线路所具有的能量大于其他线路所具有的能量值。
据此，定义线路厶在(几仍频带上的故障权重因子为：
F
G二严，心1,2,…,N(9)其中：N为线路的总条数;j为小波包的分解层数; Eg为线路在(j問子频带上所具有的能量；£(川)为所有线路在(_//)子频带上所具有的能量总和。
在式(9)的基础上，定义L，线路的综合全频带故障权重因子为
C：、叫=£,1=\,2,N(10)
其中:G为线路厶在所有的子频带当中故障权重因子为最大的子频带的个数;C为子频带的总个数。
故障权重因子的大小反映了各线路的故障的可能性大小，故障线路的故障权重因子大于健全线路的故障权重因子，据此可以作为选线的判据。
3OS-ELM理论
与传统神经网络参数调整机制不同，在极限学习机理论中，由线性方程求解出参数权值，进而随机生成权值、偏差等隐层学习参数。这使得极端学习机的训练速度、泛化性能均较传统神经网络法有大幅提升O
在线序贯极限学习机算法对于单隐层神经网络输出权重的学习过程主要分为两个部分，第一部分为初始步骤，通过少量样本得到单隐层前馈神经网络的输出权重0，第二部分为在线学习步骤，即利用单个样本或样本数据块更新在初始阶段学习到的单隐层前馈神经网络的输出权重0。
第2期嵇文路等：基于小波包全频带分析和OS-ELM的小电流单相接地故障选线113 3.1初始化阶段
给定N个训练样本0=｛（叫山）｝驚，其中输入样本X：=[%n,x a,•••,%]'e R"，输出样本t,=，如，…，俎]丁丘疋，具有厶个隐藏节点和激活函数为G（g,bjE的SLFNs的输出为：
九（乞）=£BjG（Wj,bj,x）=和=\,…,N（11）
随机生成输入层与隐藏层的连接权值吗和隐含层阈值巧，为是连接第丿•个隐层神经元和输出神经元的权值。
从n中选取部分数据集仏=｛（亿山）1芻，其中他代表初始训练样本的个数。计算初始隐藏层输出矩阵：
'G（W\,b\,%!）…G（w L,b L,x｝）_
H a=::
_G3，仇，％）…G（叫，伉,％）_|叫从
（12）已知目标输出：
丁0=[上1，…,％]hxm（13）计算初始权值仇也就是计算||H°0-T o||最小值问题。
常规条件下，隐层节点数目比起训练样本数目可以忽略，即皿＞〉厶，此时H。是非方阵，根据广义逆引理，上述线性系统的最小范数最小二乘解为仇=叫（14）式中：H'是矩阵H（）的Moore-Penrose广义逆。如果非奇异，则式（14）可写为
仇=日咒=仇民几（15）式中P0=（HjH0）_1o
工艺品加工设备
3.2在线学习阶段
对于新加入的第％+1批训练数据Gw=｛（乞，4）｝加，其中M+i为第k+1批训练样本的个数，计算隐层输出矩阵H k+｝，更新P中和0中:
P k+i=Pk-卩同+心+PEI.）-'H,+l P,
0k+i=Bk*Pk+\h：+\（Tk+\-Hk+、Bk）（16）式中验+严[耳，…，％+,]T。
当逐个获得训练数据时,M+I=1，式（16）由如下简化形式：
p_p_+
尸"“一1+仏+必心
Pk+i=伤+P*+i屁+i（4+1-hk+、Bk）（17）式中h k=[G（W],6]冋+i）,…，G（叫，饥冋+i）]。4小波包全频带分析和OS-ELM故障诊断
在某些故障情况（故障时电压初相角在峰值附近且故障接地电阻较小的情况）下，各线路的特征频带不同，能量较集中的频带也不止一个，此时如果仅仅依靠各线路在其特征频带上的故障特征信息进行选线，就不能有效的利用该线路的暂态故障特征，甚至会出现误判的情况（健全线路在其特征频带上的能量大于故障线路，健全线路的小波包系数模极大值大于故障线路），降低了选线的可靠性；实际配电网单相接地故障选线问题中，仍存在数据非一次性添加、实际配网数据存在错误或缺失等问题。针对小波选线存在的上述问题，论文给出一种小波包全频带分析和OS-ELM网络相结合的故障选线新方法。
先用小波包对暂态零序电流进行适当频带宽度的分解，然后分别计算各线路在低频带（4,0）、五次谐波所在频带（4,1）、综合全频带上的故障权重因子作为各线路的故障表征值构造故障特征向量Z，再
将Z输入到OS-ELM神经网络进行智能融合，输出一个由各线路的故障权重因子构成的故障特征向量T，综合故障权重因子最大的线路即为故障线路；同时在线选线得到的在线样本更新到样本集中，更新OS-ELM神经网络的输出权值。该方法综合利用了全频带的故障信息，可有效地解决小波选线存在的能量集中频带不同引起误判的问题；利用了选线系统的历史信息和在线信息，提高了故障选线算法的鲁棒性，流程图如图2所示。
图2小波在线序贯极限学习机诊断流程图Fig.2Diagnosis flow chart of wavelet
OS-ELM
114哈尔滨理工大学学报第26卷5诊断实例
5.1仿真模型
利用MATLAB软件搭建图1消弧线圈接地配
电网系统的故障选线的仿真模型，如图3所示。
图3仿真模型
Fig.3Simulation model
该配电仿真系统具有4条主馈线，架空线路长
度分别为厶二8km,£2=10km，电缆线路长度为
厶=5km,L4=3km架空线路、电缆线路参数如表1
所示。消弧线圈接地系统的补偿度设为p=10%,
由线路参数求出系统对地分布电容5，进而计算出
消弧线圈等效电感L=(1/1.1)x1/(3w2C x)
=0.7421H o
表1线路参数
Tab.1Line Parameters
线路类型相序
电阻/
(Q•km'1)
电感/电容/
(mH•km_1)(jxF•km'1)
架空线路
正序0.1700 1.20000.0097
角钉
零序0.2300 5.480 00.0060电缆线路
全息3d智能炫屏正序0.26500.25500.1700
零序 2.5400 1.01900.1530利用上述建立的10kV接地配电网系统的故障选线的仿真模型，分别对4条线路在电压初相角为0°、15。、30。、45。、60。、75。、90。，故障位置为线路距离母线的15%.30%.45%.60%、75%.90%长度处，过渡电阻为00.20Q、100(1.500(1.10000. 3000Q时的不同故障情况做单相接地仿真试验，总共得到4x7x6x6=l008组样本，随机选取600组样本进行训练样，其中300组用于初始化阶段训练样本,300个进行增量学习，剩余408组样本作为测试样本。
本文选取dB10小波包来分析，信号的采样频率取为4kHz，即有效频率为2kHz，选线的频带宽度定为125Hz，即分解的层数为4层，暂态信号分解为24=16个子频带，其中(4,0)代表0〜150Hz频带，(4,1)代表125〜250Hz频带，即五次谐波所在频带，以此类推,(4,15)代表1875-2000Hz频带。5.2实例分析
5.2.1故障合闸角为0。时的仿真验证
算例1：线路人故障，故障合闸角0。，过渡电阻20Q，故障点距离母线80%长度处。
系统发生上述条件下的单相接地故障后，仿真可得线路厶〜厶4的零序电流的波形，如图4所示。
time/s
图4各线路零序电流波形
Fig.4The zero sequence current waveform of each line 将暂态零序电流信号经dB10小波包4层分解后，可按照式(8)计算得到线路厶〜人的暂态零序电流信号在各子频带上的能量比，如图5所示。图中横轴坐标表示子频带的编号，纵轴坐标表示小波包分解在各子频带的能量所占所有频带总能量的比值。
0.80.8
频段号
图5暂态零序电流频谱能量分布图
Fig.5Transient zero-sequence current
spectrum energy
distribution
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议