一种T接输电线路故障支路诊断方法

一种t接输电线路故障支路诊断方法
技术领域
1.本发明属于t接输电线路故障支路诊断技术领域，具体涉及一种t接输电线路故障支路诊断方法。


背景技术：

2.电力系统的安全稳定运行是保障其它相关产业能够稳定发展的前提，对实现国民经济的稳健高速发展具有深远影响。输电线路是电力系统中重要的组成元件之一，在电力系统中承担电能的输送任务，其故障会严重影响电力系统的安全稳定运行，在国内外都曾发生过因输电线路故障而诱发电力系统崩溃的事故。所以，对输电线路故障诊断的研究具有重要意义。
3.t接输电线路被广泛应用于高压和超高压电力输电网中，常连接着大系统与大电厂，其线路传输的功率高、负荷重，因此当线路故障时，有可能导致大面积停电事故。为防止因输电线路故障引起事故的进一步扩大，以及减小因停电造成的经济损失，就要求能够快速、准确地诊断故障，进而排除故障。
4.目前对t接输电线路故障识别的研究主要有基于工频量和暂态量两种，基于工频量的t接线路故障识别方法主要利用区内三端近母线处保护单元测量的电压、电流工频信息以及输电线路分布参数信息建立保护判据，识别区内外故障。利用t接线路三端电压和电流各自故障分量矢量和的比值，识别区内外故障。利用t接线路三端电流故障分量之和以及三端电流故障分量中最大电流与另外两端电流之和的矢量差建立判据，识别区内外故障，但是判据中制动系数的选取会对故障识别的灵敏性和可靠性造成影响。在t接线路三端分别计算t节点处的正序电压，通过比较t节点正序电压叠加分量的最大幅值g1与三端正序电压叠加分量的最大幅度g2的关系识别区内外故障。当g1》g2时，发生区内故障；当g1《g2时，发生区外故障。利用t接线路三侧电压幅值差和测量阻抗特征等信息，综合电压幅值差主判据与自适应距离辅助判据识别区内外故障。基于工频量的故障诊断方法因计算数据窗较长而无法实现对故障的快速诊断。
5.基于暂态量的故障识别方法主要是对区内三端近母线处的行波保护单元测量的电压(电流)行波信息进行处理。将t接线路继电端测量得到的电压、电流信号提供给二阶泰勒-卡尔曼-傅里叶滤波器，以此估计电压、电流信号相量的瞬时值，然后通过瞬时电压、电流相量信息求取得到的正序阻抗识别区内外故障。利用t接线路三端暂态电流的余弦相似度建立判据，识别区内外故障。将小波变换应用于t接线路故障识别中，但高频噪声信号会对故障识别的效果产生影响。首先利用bior3.1小波对t接线路三端原始电流信号进行处理，然后通过对比三端运行电流与抑制电流的关系识别区内外故障。通过对比haar小波函数在t接线路每端检测到的故障电流极性来判别区内外故障。基于暂态量的故障诊断方法因其可以实现快速保护动作而受到人们的广泛关注，但基于暂态量的故障诊断方法精度不高。


技术实现要素：

6.针对现有技术中的上述不足，本发明提供的一种t接输电线路故障支路诊断方法解决了基于工频量的故障诊断方法无法实现对故障的快速诊断，基于暂态量的故障诊断方法精度不高的问题。
7.为了达到上述发明目的，本发明采用的技术方案为：一种t接输电线路故障支路诊断方法，包括以下步骤：
8.s1、采集t接输电线路区内母线端各行波保护单元的初始电压行波数据和初始电流行波数据；
9.s2、对初始电压行波数据和初始电流行波数据进行解耦和模量组合处理，得到组合电压和组合电流；
10.s3、对组合电压和组合电流分别进行s变换，得到电压时频复矩阵和电流时频复矩阵；
11.s4、根据电压时频复矩阵和电流时频复矩阵，计算两两行波保护单元各频率对应时窗下无功功率均值和的比值；
12.s5、根据无功功率均值和的比值，构建t接输电线路故障特征向量样本集，结合随机森林构建故障识别模型；
13.s6、采用故障识别模型识别t接输电线路的故障支路。
14.进一步地，所述步骤s2包括以下分步骤：
15.s21、对初始电压行波数据和初始电流行波数据分别进行解耦处理，得到模电压和模电流；
16.s22、将模电压和模电流进行组合，得到组合电压和组合电流。
17.上述进一步地方案的有益效果为：在三相输电系统中，各相电压和各相电流之间的耦合会对电压电流造成影响，因此需要对相电压和相电流进行解耦处理，本发明采用克拉克相模变换对相电压和相电流进行解耦处理，再利用组合模量来反应t接线路的各种故障类型。
18.进一步地，所述步骤s22中组合公式为：
[0019][0020]
其中，δuz为组合电压，δu
α
为克拉克α模电压，δu
β
为克拉克β模电压，δiz为组合电流，δi
α
为克拉克α模电流，δi
β
为克拉克β模电流。
[0021]
进一步地，所述步骤s4包括以下分步骤：
[0022]
s41、根据电压时频复矩阵和电流时频复矩阵，得到各行波保护单元各频率下各采样点的无功功率q
mn
(l)，其中，m为行波保护单元编号，n为频率编号，l为采样点编号，采样点对应电压、电流故障后0.1ms时间段内的数据；
[0023]
s42、将各频率下各采样点的无功功率q
mn
(l)以11个采样点为固定数据窗长，1个采样点为滑移尺度因子进行截取，得到多个时窗的功率序列q
mnk
(x)，其中，k为时窗编号，x＝1,2,3,
…
,11为时窗的11个采样点编号；
[0024]
s43、根据功率序列q
mnk
(x)，分别计算各行波保护单元下在每个时窗内无功功率均
值
[0025][0026]
s44、根据各行波保护单元下在每个时窗内无功功率均值计算各行波保护单元在时窗k下无功功率均值和q
m_k
；
[0027][0028]
s45、根据各行波保护单元在时窗k下无功功率均值和q
m_k
，计算两两行波保护单元对应时窗下无功功率均值和的比值和其中，为行波保护单元tr1在时窗k下与行波保护单元tr2在时窗k下无功功率均值和的比值，为行波保护单元tr1在时窗k下与行波保护单元tr3在时窗k下无功功率均值和的比值，为行波保护单元tr2在时窗k下与行波保护单元tr3在时窗k下无功功率均值和的比值。
[0029]
上述进一步地方案的有益效果为：采用电压、电流行波数据进行故障诊断，有利于提升故障诊断的稳定性，采用多个频率多个时窗的数据，可进一步提升故障诊断的抗干扰能力。
[0030]
进一步地，所述步骤s41中无功功率q
mn
(l)公式为：
[0031][0032]
其中，q
mn
(l)为行波保护单元m在频率fn下第l采样点的无功功率，im[
·
]为求取虚部函数，为电压时频复矩阵中行波保护单元m在频率fn下第l采样点的电压相量，为的共轭相量，为电流时频复矩阵中行波保护单元m在频率fn下第l采样点的电流相量。
[0033]
进一步地，所述步骤s5包括以下分步骤：
[0034]
s51、将两两行波保护单元各频率对应时窗下无功功率均值和的比值组成t接输电线路故障特征向量；
[0035]
s52、将t接输电线路故障特征向量构成训练样本集；
[0036]
s53、采用训练样本集训练和测试随机森林故障识别模型，得到训练完成的故障识别模型。
[0037]
本发明的有益效果为：
[0038]
1、本发明利用t接输电线路区内三端测点的初始电压行波数据和初始电流行波数据，首先对电压、电流行波数据解耦并对解耦后的模量信息进行组合，再基于s变换计算t接输电线路多个频率下对应短时时窗序列的初始行波无功功率均值和，然后将两两行波保护单元对应短时序列的无功功率均值和作比，并将其组成t接输电线路故障特征向量样本集，
采用样本集对随机森林故障识别模型进行训练和测试，以此实现t接输电线路故障支路的诊断。
[0039]
2、利用区内外初始行波无功功率的差异反映不同支路故障特征，结合随机森林的学习及泛化能力对故障进行诊断，本发明相较于传统t接输电线路故障诊断方法，不仅能实现区内外故障的识别，还能准确识别出t接输电线路具体的故障支路，相比于暂态量的故障诊断方法精度更高。
[0040]
3、本发明故障诊断的数据窗长仅0.1ms，相较于传统故障诊断使用的全周(20ms)或半周(10ms)傅氏算法数据窗，大大缩短了算法计算时间，能实现t接输电线路故障的快速诊断。
附图说明
[0041]
图1为一种t接输电线路故障支路诊断方法的流程图；
[0042]
图2位本发明中500kv的t接输电线路示意图；
[0043]
图3为本发明中t接输电线路区内co支路故障时的彼得逊等效电路图；
[0044]
图4为本发明中t接线路区外cf支路故障时彼得逊等效电路图；
[0045]
图5为本发明中t接输电线路区内co支路故障时行波保护单元tr1相关波形图；
[0046]
图6为本发明中t接输电线路区内co支路故障时行波保护单元tr2相关波形图；
[0047]
图7为本发明中t接输电线路区内co支路故障时行波保护单元tr3相关波形图；
[0048]
图8为本发明中t接输电线路区外cf支路故障时行波保护单元tr1相关波形图；
[0049]
图9为本发明中t接输电线路区外cf支路故障时行波保护单元tr2相关波形图；
[0050]
图10为本发明中t接输电线路区外cf支路故障时行波保护单元tr3相关波形图；
[0051]
图11为本发明中训练模型的流程图。
具体实施方式
[0052]
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
[0053]
如图1所示，一种t接输电线路故障支路诊断方法，包括以下步骤：
[0054]
s1、采集t接输电线路区内近母线端各行波保护单元的初始电压行波数据和初始电流行波数据；
[0055]
图2为500kv的t接输电线路，定义图2中ao、bo、co三条支路为t接输电线路的区内支路，ad、be、cf三条支路为t接输电线路的区外支路。
[0056]
t接输电线路由区内支路ao、bo、co和区外支路ad、be、cf组成，电压等级为500kv，行波保护单元tr1～tr3分别安装在区内支路靠近a、b、c三端母线处。当故障发生在区内ao支路f1处时，行波沿线路从故障点向两侧传播，在线路波阻抗不连续处发生折反射。
[0057]
根据行波传播理论，设t
0m
(m＝1,2,3)分别为初始行波首次到达a、b、c3端的时刻，t
1m
(m＝1,2,3)为行波在线路波阻抗不连续处发生折反射后第2次到达a、b、c三端的时刻；在t
0m
～t
1m
时间段内，区内支路近a、b、c 3端处的行波保护单元trm(m＝1,2,3)获取的故障行波
称为初始电压行波数据δum(m＝1,2,3)和初始电流行波数据δim(m＝1,2,3)。
[0058]
s2、对初始电压行波数据和初始电流行波数据进行解耦和模量组合处理，得到组合电压和组合电流；
[0059]
所述步骤s2包括以下分步骤：
[0060]
s21、对初始电压行波数据和初始电流行波数据分别进行解耦处理，得到模电压和模电流；
[0061]
s22、将模电压和模电流进行组合，得到组合电压和组合电流。
[0062]
所述步骤s22中组合公式为：
[0063][0064]
其中，δuz为组合电压，δu
α
为克拉克α模电压，δu
β
为克拉克β模电压，δiz为组合电流，δi
α
为克拉克α模电流，δi
β
为克拉克β模电流。
[0065]
s3、对组合电压和组合电流分别进行s变换，得到电压时频复矩阵和电流时频复矩阵；
[0066]
s4、根据电压时频复矩阵和电流时频复矩阵，计算两两行波保护单元各频率对应时窗下无功功率均值和的比值；
[0067]
所述步骤s4包括以下分步骤：
[0068]
s41、根据电压时频复矩阵和电流时频复矩阵，得到各行波保护单元各频率下各采样点的无功功率q
mn
(l)，其中，m为行波保护单元编号，n为频率编号，l为采样点编号，采样点对应电压、电流故障后0.1ms时间段内的数据；
[0069]
所述步骤s41中无功功率q
mn
(l)公式为：
[0070][0071]
其中，q
mn
(l)为行波保护单元m在频率fn下第l采样点的无功功率，im[
·
]为求取虚部函数，为电压时频复矩阵中行波保护单元m在频率fn下第l采样点的电压相量，为的共轭相量，为电流时频复矩阵中行波保护单元m在频率fn下第l采样点的电流相量。
[0072]
和为：在电压时频复矩阵和电流时频复矩阵分别选择10个不同频率fn(n＝1,2,
…
,10)下故障后0.1ms时间段内的电压相量和电流相量。
[0073]
在本实施例中，频率fn具体为：f1＝6khz，f2＝7khz，f3＝8khz，f4＝9khz，f5＝10khz，f6＝11khz，f7＝12khz，f8＝13khz，f9＝14khz，f
10
＝15khz；采样点个数l＝1,2，
…
，20，采样点个数l为步骤s1中所指采样过程的采样点。
[0074]
s42、将各频率下各采样点的无功功率q
mn
(l)以11个采样点为固定数据窗长，1个采样点为滑移尺度因子进行截取，得到多个时窗的功率序列q
mnk
(x)，其中，k为时窗编号，x＝1,2,3,
…
,11为时窗的11个采样点编号；
[0075]
s43、根据功率序列q
mnk
(x)，分别计算各行波保护单元下在每个时窗内无功功率均值
[0076][0077]
s44、根据各行波保护单元下在每个时窗内无功功率均值计算各行波保护单元在时窗k下无功功率均值和q
m_k
；
[0078][0079]
s45、根据各行波保护单元在时窗k下无功功率均值和q
m_k
，计算两两行波保护单元对应时窗下无功功率均值和的比值和其中，为行波保护单元tr1在时窗k下与行波保护单元tr2在时窗k下无功功率均值和的比值，为行波保护单元tr1在时窗k下与行波保护单元tr3在时窗k下无功功率均值和的比值，为行波保护单元tr2在时窗k下与行波保护单元tr3在时窗k下无功功率均值和的比值。
[0080]
在本实施例中，步骤s4的具体原理过程为：
[0081]
设电流极性为流出母线为正，流进母线为负，功率的正负可根据各母线关联线路电流极性定义。
[0082]
当故障发生在t接线路区内co支路f1处时，t接线路的彼得逊等值电路如图3所示，其中为故障点附加网络电压，分别为测量到的母线a的初始电压、电流行波，线路oa、ob、oc、ad、be、cf的波阻抗分别为z
c1
～z
c6
，由于线路波阻抗近似为实常数，则令z
c1
＝z
c2
＝z
c3
＝z
c4
＝z
c5
＝z
c6
≈r；母线c对地等值电容阻抗为z
cc
。
[0083]
由初始行波复功率的定义可知该线路母线c端的初始行波复功率为：
[0084][0085]
区内oc故障时，由图3彼得逊等效电路可知：
[0086][0087]
c端行波保护单元tr3测量到的复功率为：
[0088][0089]
其中，pc为线路初始行波有功功率，qc为线路初始行波无功功率，则t接线路区内故障时：
[0090][0091]
当故障发生在t接线路区外cf支路f2处时，t接线路的彼得逊等值电路如图4所示，由图4可知：
[0092][0093]
则行波保护单元tr3测量到的复功率为：
[0094][0095]
式中，pc为线路初始行波有功功率，qc为线路初始行波无功功率，则t接线路区外故障时：
[0096][0097]
当t接输电线路区内co支路距o点110km处发生ac相接地故障，故障初始角为60度，过渡电阻为100ω时，以s变换后频率为15khz对应信号为例，各行波保护单元trm相应波形分别如图5～图7所示，其中
△im
、
△
um为trm(m＝1,2,3)的初始电流、电压行波，pm为行波保护单元trm(m＝1,2,3)初始行波无功功率分布波形。由图5～图7可知，当t接线路区内支路co发生故障时，行波保护单元trm(m＝1,2,3)初始电压、电流行波波头附近数据所求无功功率大小均为负值。
[0098]
当t接线路区外支路cf距o点220km处发生bc相间故障，故障初始角为5
°
，过渡电阻350ω，各行波保护单元的相关波形图分别如图8～图10所示(以s变换后15khz频率对应的信号为例)。由图8～图10可知，当t接线路区外支路cf发生故障时，行波保护单元tr3初始电压、电流行波波头附近数据所求无功功率大小均为正值，行波保护单元trm(m＝1,2)初始电压、电流行波波头附近数据所求无功功率大小均为负值。
[0099]
s5、根据无功功率均值和的比值，构建t接输电线路故障特征向量样本集，结合随机森林构建故障识别模型；
[0100]
所述步骤s5包括以下分步骤：
[0101]
s51、将两两行波保护单元各频率对应时窗下无功功率均值和的比值组成t接输电线路故障特征向量；
[0102]
特征向量q为
[0103]
s52、将t接输电线路故障特征向量构成训练样本集；
[0104]
s53、采用训练样本集训练和测试随机森林故障识别模型，得到训练完成的故障识别模型。
[0105]
在本实施例中，故障识别模型为多个决策树模型组成的分类器，如图11所示，步骤s53具体为：
[0106]
对训练样本集x进行bootstrap采样，采样j次后获得bootstrap训练子集xj(j＝1,2,
…
,j)。针对每个训练子集xj采用cart构建决策树模型hj，最终获得一组决策树模型组成的分类器{h1,
…
,hj,
…
,hj}。当测试样本输入训练好后的分类器时，通过建立好的j个决策
树模型进行类别投票，得票最高的类别作为测试样本的最后输出类别，分类决策如下：
[0107][0108]
其中，hj(x)为第j棵决策树模型；i{
·
}为示性函数，当满足表达式时为1，否则为0；y为类别标签，arg为取平均值，f为测试样本在分类器中的输出类别。
[0109]
为保障本发明在故障支路识别中的有效性与可靠性，选择在不同故障初始角、不同过渡电阻、不同故障类型以及不同故障距离4种情况下对t接输电线路6条支路分别仿真5组故障，共计得到120组故障特征向量，以此组成故障支路识别的训练样本集。
[0110]
s6、采用故障识别模型识别t接输电线路的故障支路。
[0111]
在不同故障初始角、不同过渡电阻、不同故障类型以及不同故障距离4类情况下对t接输电线路区内外6条支路分别仿真4组不同于训练样本的故障，得到4个测试样本集，每个样本集包含24组故障特征向量，分别将4组测试样本集输入训练好的故障识别模型，对故障支路进行识别。
[0112]
本发明利用t接输电线路区内三端测点的初始电压行波数据和初始电流行波数据，首先对电压、电流行波数据解耦并对解耦后的模量信息进行组合，再基于s变换计算t接输电线路多个频率下对应短时时窗序列的初始行波无功功率均值和，然后将两两行波保护单元对应短时序列的无功功率均值和作比，并将其组成t接输电线路故障特征向量样本集，采用样本集对随机森林故障识别模型进行训练和测试，以此实现t接输电线路故障支路的诊断。