一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法

摘要

本发明提供一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法，属于电力系统继电保护技术领域。当同塔双回输电线路某回线路发生故障时，检测并记录各条线路电流行波数据；利用相模变换矩阵分别求取各出线α模电流行波并进行小波变换，选取高频第一尺度小波系数作为研究对象；计算出短时窗内各条线路第一尺度小波系数的相对极性，并形成故障选线判别矩阵R，根据R中的元素即可判断出故障线路。仿真分析表明，对于同塔双回线中单回线故障，该方法可以准确判断出故障线。

权利要求

1.一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法，其特征在于：当同塔双回输电线路某回线路发生故障时，检测并记录各条线路电流行波数据；利用相模变换矩阵分别求取各条线路α模电流行波并进行小波变换，选取高频第一尺度小波系数作为研究对象；计算出短时窗内各条线路第一尺度小波系数的相对极性，并形成故障选线判别矩阵R，根据R中的元素判断出故障线路；

具体实施步骤为：

(1)当同塔双回输电线路发生故障时，利用各出线侧的电流表记录各条线路的电流行波数据ia、ib、ic，利用相模变换矩阵Q分别求取各条线路α模电流行波iα，

(2)将各条线路α模电流行波iα进行小波分解，得到各条出线的第一尺度小波系数d1_n，其中n＝1,2,3…；

(3)根据式(2)计算各条出线第一尺度小波系数间的相对极性；

式中，d1_i表示线路i高频第一尺度小波系数；d1_j表示线路j高频第一尺度小波系数；rij表示d1_i和d1_j极性比较的结果，i＝1,2,3…；j＝1,2,3…；rij＝-1表示d1_i和d1_j极性相反，rij＝1表示d1_i和d1_j极性相同，k表示小波系数长度；

(4)根据式(3)计算，形成同塔双回输电线路故障选线矩阵R：

根据每行元素为“-1”的个数判断故障线路：若R中所有的元素都为“1”，则判断为母线故障；

若第i行中为“-1”的元素为n-1，则第i条线路为故障电路。

说明书

一种基于电流行波短时窗小波系数极性比较的同塔双回线选

线方法

技术领域

[0001] 本发明涉及一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法，属于电力系统继电保护技术领域。

背景技术

[0002] 同塔双回输电线路因其工程造价低、占用走廊宽度小、建设周期短、运行维护简单、经济效益显著等优点，因此得到了广泛的应用。同塔双回线输电方式有120种故障类型，其中，跨线故障有98种，虽然，从种类上来讲，跨线故障占多数；但是，从故障发生的概率上来讲，同塔双回输电线路发生单回线故障的概率占到了70％-80％。

[0003] 所以，在同塔双回线输电的过程中，如果能在线路发生故障的情况下，快速、准确的选出故障线路，不仅能通过快速切除故障线路，以免故障的扩大化。而且也可以减少停电时间，提高系统的可靠性。本发明利用电流行波短时窗小波系数极性的比较矩阵判断故障线路，不仅原理简单，而且可以准确判断故障位于母线还是某一条出线。

发明内容

[0004] 本发明要解决的技术问题是针对同塔双回线发生单回线故障的情况，提出一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法。

[0005] 本发明的技术方案是：一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法，当同塔双回输电线路某回线路发生故障时，检测并记录各条线路电流行波数据；利用相模变换矩阵分别求取各条线路α模电流行波并进行小波变换，选取高频第一尺度小波系数作为研究对象；计算出短时窗内各条线路第一尺度小波系数的相对极性，并形成故障选线判别矩阵R，根据R中的元素判断出故障线路。

[0006] 具体实施步骤为：

[0007] (1)当同塔双回输电线路发生故障时，利用各出线侧的电流表记录各条线路的电流行波数据ia、ib、ic，利用相模变换矩阵Q分别求取各条线路α模电流行波iα，[0008]

[0009] (2)将各条线路α模电流行波iα进行小波分解，得到各条出线的第一尺度小波系数d1_n，其中n＝1,2,3…；

[0010] (3)根据式(2)计算各条出线第一尺度小波系数间的相对极性；

[0011]

[0012] 式中，d1_i表示线路i高频第一尺度小波系数；d1_j表示线路j高频第一尺度小波系数；rij表示d1_i和d1_j极性比较的结果，i＝1,2,3…；j＝1,2,3…；rij＝-1表示d1_i和d1_j极性相反，rij＝1表示d1_i和d1_j极性相同，k表示小波系数长度；

[0013] (4)根据式(3)计算，形成同塔双回输电线路故障选线矩阵R：

[0014]

[0015] 根据每行元素为“-1”的个数判断故障线路：

[0016] 若R中所有的元素都为“1”，则判断为母线故障；

[0017] 若第i行中为“-1”的元素为n-1，则第i条线路为故障电路。

[0018] 本发明的原理是：……

[0019] 本发明的有益效果是：

[0020] 1、本发明比较的是小波系数在短时窗内的相对极性，与小波变换模极大值相比较，可靠性高。

[0021] 2、本发明仅需利用电流量即可实现故障选线，不需引入电压量。

附图说明

[0022] 图1为本发明同塔双回输电线路结构示意图：图中，P、M、N、Q为系统的母线，CE为母线对地的杂散电容，F1、F2、F3、F4分别为线路PM段、双回线MN段Ⅰ、Ⅱ回和母线M处仿真故障点。PM段长200km，MN段长140km，NQ段长260km。

[0023] 图2为本发明的选线流程图；

[0024] 图3为本发明实施例1中，F1故障时1、2和3保护安装处α模电流行波的高频第一尺度小波系数曲线图；

[0025] 图4为本发明实施例2中，F2故障时1、2和3保护安装处α模电流行波的高频第一尺度小波系数曲线图；

[0026] 图5为本发明实施例3中，F3故障时1、2和3保护安装处α模电流行波的高频第一尺度小波系数曲线图；

[0027] 图6为本发明实施例4中，F4故障时1、2和3保护安装处α模电流行波的高频第一尺度小波系数曲线图。

具体实施方式

[0028] 下面结合附图和具体实施方式，对本发明作进一步说明。

[0029] 一种基于电流行波短时窗小波系数极性比较的同塔双回线选线方法，当同塔双回输电线路某回线路发生故障时，检测并记录各条线路电流行波数据；利用相模变换矩阵分别求取各条线路α模电流行波并进行小波变换，选取高频第一尺度小波系数作为研究对象；计算出短时窗内各条线路第一尺度小波系数的相对极性，并形成故障选线判别矩阵R，根据R中的元素判断出故障线路。

[0030] 具体实施步骤为：

[0031] (1)当同塔双回输电线路发生故障时，利用各出线侧的电流表记录各条线路的电流行波数据ia、ib、ic，利用相模变换矩阵Q分别求取各条线路α模电流行波iα，[0032]

[0033] (2)将各条线路α模电流行波iα进行小波分解，得到各条出线的第一尺度小波系数d1_n，其中n＝1,2,3…；

[0034] (3)根据式(2)计算各条出线第一尺度小波系数间的相对极性；

[0035]

[0036] 式中，d1_i表示线路i高频第一尺度小波系数；d1_j表示线路j高频第一尺度小波系数；rij表示d1_i和d1_j极性比较的结果，i＝1,2,3…；j＝1,2,3…；rij＝-1表示d1_i和d1_j极性相反，rij＝1表示d1_i和d1_j极性相同，k表示小波系数长度；

[0037] (4)根据式(3)计算，形成同塔双回输电线路故障选线矩阵R：

[0038]

[0039] 根据每行元素为“-1”的个数判断故障线路：

[0040] 若R中所有的元素都为“1”，则判断为母线故障；

[0041] 若第i行中为“-1”的元素为n-1，则第i条线路为故障电路。

[0042] 实施例1：500kV双回输电线路系统图如图1所示。其线路参数如下：线路全长PM段200km，MN段140km，NQ端260km。故障位置：PM段距M端10km处发生A相接地故障。接地阻抗0Ω，初始故障角90°,采样率为1MHz。

[0043] (1)根据说明书中的步骤(1)得到3条线路的电流行波数据ia、ib、ic，并利用相模变换矩阵Q分别求取各条线路α模电流行波iα。

[0044] 其中：

[0045] (2)将(1)中3条线路α模电流行波iα进行小波分解，得到各条出线的频率范围在250KHz～500KHz的第一尺度小波系数d1_n，其曲线图如3所示。

[0046] (3)根据说明书中的步骤(3)计算各条出线第一尺度小波系数间的相对极性；

[0047]

[0048] 得r11＝r22＝r33＝1，r12＝r21＝-1，r13＝r31＝-1，r23＝r32＝1。

[0049] (4)根据说明书中的步骤(4)得到同塔双回输电线路故障选线矩阵R如下式：

[0050]

[0051] 从上式可以看出只有第1行元素的“-1”个数为2，按照图2所示的流程图，结合图1中保护装置的位置可知，线路PM段为故障线路。

[0052] 实施例2：500kV双回输电线路系统图如图1所示。其线路参数如下：线路全长PM段200km，MN段140km，NQ端260km。故障位置：MN段Ⅰ回线距M端60km处发生A相接地故障。接地阻抗10Ω，初始故障角90°,采样率为1MHz。

[0053] (1)根据说明书中的步骤(1)得到3条线路的电流行波数据ia、ib、ic，并利用相模变换矩阵Q分别求取各条线路α模电流行波iα。

[0054] 其中：

[0055] (2)将(1)中3条线路α模电流行波iα进行小波分解，得到各条出线的频率范围在250KHz～500KHz的第一尺度小波系数d1_n，其曲线图如3所示。

[0056] (3)根据说明书中的步骤(3)计算各条出线第一尺度小波系数间的相对极性；

[0057]

[0058] 得r11＝r22＝r33＝1，r12＝r21＝-1，r13＝r31＝1，r23＝r32＝-1。

[0059] (4)根据说明书中的步骤(4)得到同塔双回输电线路故障选线矩阵R如下式：

[0060]

[0061] 从上式可以看出只有第2行元素的“-1”个数为2，按照图2所示的流程图，结合图1中保护装置的位置可知，线路MN段Ⅰ回线为故障线路。

[0062] 实施例3：500kV双回输电线路系统图如图1所示。其线路参数如下：线路全长PM段200km，MN段140km，NQ端260km。故障位置：MN段Ⅱ回线距M端90km处发生A相接地故障。接地阻抗0Ω，初始故障角90°,采样率为1MHz。

[0063] (1)根据说明书中的步骤(1)得到3条线路的电流行波数据ia、ib、ic，并利用相模变换矩阵Q分别求取各条线路α模电流行波iα。

[0064] 其中：

[0065] (2)将(1)中3条线路α模电流行波iα进行小波分解，得到各条出线的频率范围在250KHz～500KHz的第一尺度小波系数d1_n，其曲线图如3所示。

[0066] (3)根据说明书中的步骤(3)计算各条出线第一尺度小波系数间的相对极性；

[0067]

[0068] 得r11＝r22＝r33＝1，r12＝r21＝1，r13＝r31＝-1，r23＝r32＝-1。

[0069] (4)根据说明书中的步骤(4)得到同塔双回输电线路故障选线矩阵R如下式：

[0070]

[0071] 从上式可以看出只有第3行元素的“-1”个数为2，按照图2所示的流程图，结合图1中保护装置的位置可知，线路MN段Ⅱ回线为故障线路。

[0072] 实施例4：500kV双回输电线路系统图如图1所示。其线路参数如下：线路全长PM段200km，MN段140km，NQ端260km。故障位置：母线M发生故障，采样率为1MHz。

[0073] (1)根据说明书中的步骤(1)得到3条线路的电流行波数据ia、ib、ic，并利用相模变换矩阵Q分别求取各条线路α模电流行波iα。

[0074] 其中：

[0075] (2)将(1)中3条线路α模电流行波iα进行小波分解，得到各条出线的频率范围在250KHz～500KHz的第一尺度小波系数d1_n，其曲线图如3所示。

[0076] (3)根据说明书中的步骤(3)计算各条出线第一尺度小波系数间的相对极性；

[0077]

[0078] 得r11＝r12＝r21＝r22＝r13＝r31＝r23＝r32＝r33＝1。

[0079] (4)根据说明书中的步骤(4)得到同塔双回输电线路故障选线矩阵R如下式：

[0080]

[0081] 从上式可以看出R中所有的元素都为“1”，按照图2所示的流程图，结合图1中保护装置的位置可知，故障发生在母线M。

[0082] 上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。