一种考虑多点故障的高压输电线路行波测距方法
【专利摘要】本发明公开一种考虑多点故障的高压输电线路行波测距方法，通过在高压输电线路上装设若干组测距装置，利用故障电流检测装置采集各检测点的故障电压电流信号，运用小波变换求取模极大值的方法获取检测点的突变量信号，采用电流偏离度与故障行波到达时刻相结合方法判断故障区间，在上述获得的故障区间及获取的突变量基础上，依靠故障发生后故障点两侧行波波头所到达的位置拟合行波波头传播特性曲线实现对故障测距的准确定位。本发明能够解决常规行波测距方法不能够完全对输电线路所有故障准确定位的问题，有效解决针对输电线路同时发生多处故障的精确定位问题，保障电力线路故障准确高效修复以及保证电力系统安全稳定运行。
【专利说明】
一种考虑多点故障的高压输电线路行波测距方法
技术领域
[0001] 本发明属于电力系统继电保护技术的领域，具体涉及一种运用获取有效区间内的 故障行波到达故障点时刻并结合检测点故障行波的时序特征来构造行波故障特性曲线获 取多点故障位置距离的方法。
【背景技术】
[0002] 目前，对于输电线路的在故障定位主要有故障分析法和行波法两大类。其中行波 法因具有相对简单的模型，抗干扰性强，有较低的测距误差等特点得到广泛运用。行波法主 要分为单端测距法以及双端测距法。
[0003] 单端测距方法只需要获取线路一侧的电压或电流的行波信号，设备简单，但要区 分行波来自故障点反射还是对端母线反射。然而对于行波而言，极性和幅值是其最主要的 识别特征，并且线路结构和故障情况不同，很多情况下采用单端测距普适性不强。双端测距 法利用故障行波到达两侧的时间差来确定故障距离，记录输电线路两端的故障行波到达时 亥IJ，模型相对简单，可靠性和测距精度均很高，但需要实现输电线路两端信息的同步传输， 成本较高。
[0004] 当输电线路某一处发生故障时，尽管在靠近线路端点、或过渡电阻较大的情况下， 常规行波测距方法的测量精度会受到一定的影响，但是误差范围都不是很大，基本能够满 足故障测距的要求，但当线路发生两处故障时，由于线路上叠加了两个行波波头从线路不 同地点朝两端传输，同时在两个故障点都将发生行波的折反射，常规的行波测距在特定的 故障下无法正确获取靠近线路首段故障的距离，不能对故障进行正确定位。因此，迫切需要 提出一种考虑多点故障的高压输电线路行波测距方法。
【发明内容】
[0005] 本发明的目的在于解决目前输电线路发生多处故障时，常规的行波测距在特定的 故障下无法正确获取距离故障点的初始行波和反射行波到达线路首端的时间，从而不能对 故障进行正确定位问题。
[0006] 为解决上述技术问题，发明人采用了如下的技术方案：
[0007] -种考虑多点故障的高压输电线路行波测距方法，包括以下步骤：
[0008] 采用分布在输电线路上的若干故障检测装置得到故障电压、电流信号，对其进行 相模解耦变换，接着运用小波变换获取模极大值的方法得到各个检测点的故障电压、电流 突变量信息，利用获得的故障电压、电流信号以及对输电线路各个检测点电流之间的偏离 度得到故障点所处的区间，利用获取的电压、电流突变量和故障所处的有效区间，得到故障 行波波头传输特性曲线，从而得到故障位置。
[0009] 具体步骤为：
[0010] 第一步，基于分布式的行波故障测距系统读取高压输电线路上的各检测装置是否 发出故障信号，判断线路是否发生故障。若发生故障则得到故障暂态电压、电流信号，若未 发生故障则系统继续等待故障信号。
[0011] 第二步，得到检测点信号突变量信息。
[0012] 步骤2-1)对于高压输电线路，直接测量得到的各相电流、电压之间存在耦合关系， 因此，首先采用克拉克变换对输电线路电压电流量进行相模变换，使其解耦，得到独立的电 压、电流量，即将保护安装处的电流、电压值乘以克拉克变换矩阵得到解耦电压、电流分量。
[0014]步骤2-2)对于高压输电线路故障时产生的行波信号，其频率会随着传递时间的变 化而发生变化，因此，采用离散的小波变换对故障行波进行分析。表示如下：
[0016] 式中，b为小波变换平移因子，a为其尺度因子，供为小波变换基小波，t，a，b取离散 值，f (t)为解耦之后故障电流信号。
[0017] 采用小波变换模极大值方法分析故障信号突变点，将时间to代入式(2)得到小波 变换后的值，并对te(t『S，t〇+S)所有变换后的值进行检验，当满足式：|w sf(t)|彡|Wsf(t〇) (3)
[0018] 得到时间to为故障信号突变点的时间，Wsf (to)为模极大值。
[0019]第三步，得到故障点所处的区间。
[0020]步骤3-1)首先，检测点Ah和4"故障行波波首到达时间，若满足如下条件：
[0022] 初步判断得到故障所在的区间为之间。其中，分别为故障行波波首 到达An-^An的时间，1^、1^分别为输电线路故障测距检测点Aw和AIM端的距离，v为行波 波速。若不满足式(4)则可判断有一故障在检测点n-1和η之间。
[0023] 步骤3-2)利用欧式空间计算位于故障同一侧两个检测点iAk和IAg的电流偏离度dkg 和位于故障点不同侧的两个检测点1?和I Aq的电流偏离度dM (Kk，g<n)。其中lAk，k=l，2, 3···η为各检测点的工频故障电流。设每组检测装置采集的数据长度为L，偏离度计算公式如 下：
[0025]将得到的dkg、dPq进行处理比较，由（5)得到若满足如下关系
[0026] dPq>Kidkg，（Kk，g或i<k，g<n)且（p>i>q或p<i<q) (6)
[0027] 其中L为偏离系数，且KOI，得到区段故障区间位于检测点i和i+1之间。
[0028]第四步，得到故障行波波头传输特性曲线。
[0029]步骤4-1)获取故障行波时序。根据附图1，利用输电线路检测装置得到各个检测点 检测到的故障行波到达时间，公式表示如下：
[0032] tzl，tz2，…，tzn为从故障点F1发出的正向行波经折反射后依次到达检测点的时间 与故障发生时刻的差值，tfl，tf2，…，tfn为反向行波经折反射后依次到达的时间与故障发生 时刻的差值。. . .，νζη为从故障点F2发出的正向行波经折反射后依次到达检测点 的时间与故障发生时刻的差值，...为反向行波经折反射后依次到达的时间 与故障发生时刻的差值。
[0033] 又因为Lx2-Lxi^：0且j，可进一步得到检测点Αι，A2，…，Ai和Aj+i，Aj+2，…，A n测得 的初始故障行波时序：
[0036] 其中，tzk为各检测点检测到的第一个行波到达时间。接着确定两故障点之间 A1+1，…，、行波到达时间，分以下三种情况：当线路发生两处故障区间相同，该测量点行波 到达时间满足式(9);当线路发生两处故障区间相邻，若L 1+1-LX1>LX2-L1+1即tWfG+υ时， 1/ z(i+i)>tf(i+i)。若Li+i_Lxi>Lx2-Li+i即tzi〉!/ f(j+i)时，1/ z(i+i)>tf(i+i)，根据测量点Ai，Aj+i检测 到的第一个行波时间判断测量点A 1+1的行波时序;当线路发生两处故障区间不相邻，对于测 量点Ak(k= i+1，i+2，…，j)，若Lk-Lxi>Lx2-Lk，1/ zk〈i/ fk。若Lk-Lxi>Lx2-Lk，1/ zk>i/ fk。因此，最 终得到整条输电线路每个故障点行波波头到达时间。
[0037] 步骤4-2)根据得到的故障所在的区间及结合步骤4-1)故障行波的时序拟合各时 亥IJ行波波头达到位置距离线路首端Μ和线路尾端N的距离曲线，即得到曲线XM=nu(t)，XM = m2 (t)，XN = m (ti)，XN = n2 (t)。接着统一两曲线的纵坐标，即都以距线路首段Μ的距离为纵坐 标，根据输电线路满足Xm+Xn = L，分别将XN=m (t)，XN = n2 (t)转化为曲线L-Xn = f! (t)，L-XN =f2(t);得到的故障行波波头传输特性曲线如附图。
[0038] 第五步，得到故障位置。
[0039]联立XMiimUhL-XNzfKt)，两条曲线的交点所在的位置纵坐标就是故障点距输 电线路首端Μ的距离，所在位置的横坐标是故障时刻。XM = m2(t)，L-XN = f2(t)同理，从而完 成了多点故障定位。
【附图说明】
[0040]图1是本发明分布式故障测距检验系统示意图。
[0041 ]图2是本发明故障测距系统结构框图。
[0042]图3是本发明分布式的改进故障测距算法流程图。
[0043]图4是本发明故障行波波头传输特性曲线拟合流程图。
[0044] 图5是本发明系统发生两处故障的仿真模型图。
【具体实施方式】
[0045] 结合图1、图2、图3以及图4所示，本发明公开了一种考虑多点故障的高压输电线路 行波测距方法，通过在高压输电线路上装设若干组测距装置，利用故障电流检测装置采集 各检测点的故障电压电流信号，运用小波变换求取模极大值的方法获取检测点的突变量信 号，采用电流偏离度与故障行波到达时刻相结合方法判断故障区间，在上述获得的故障区 间及获取的突变量基础上，依靠故障发生后某一时刻故障点两侧行波波头所到达的位置拟 合行波波头传播特性曲线实现对故障测距的准确定位。
[0046]故障测距系统结构图如附图2,得到正确的故障测距结果方法如下：
[0047]第一步，通过装设在输电线路上的若干三相故障检测装置，为了提高故障暂态信 号的准确性采用罗氏线圈进行采集，得到故障暂态电压、电流信号Ua，Ub，Uc，ia，ib，ic。
[0048]第二步，得到检测点信号突变量信息：
[0049] 步骤2-1)高压输电线路，直接测量得到的各相电流ia，ib,ic、电压U a，Ub，Uc之间存 在耦合关系，因此，首先采用克拉克变换对输电线路电压电流量进行相模变换，使其解耦， 得到独立的电压Uo，山，U2、电流量i〇，i!，i 2，即将保护安装处的电流、电压值乘以克拉克变换 矩阵得到解耦电压、电流分量。
[0051]步骤2-2)高压输电线路故障时产生的电流信号i2,其频率会随着传递时间 的变化而发生变化，因此，采用离散的小波变换对故障行波进行分析。表示如下：
[0053] 式中，b为小波变换平移因子，a为其尺度因子，供为小波变换基小波，t，a，b取离散 值，f(t)为解耦之后故障电流信号
[0054] 采用小波变换模极大值方法分析故障信号突变点，将时间to代入式(2)得到小波 变换后的值，并对te(t『S，t〇+S)所有变换后的值进行检验，当满足式：|w sf(t)|彡|Wsf(t〇) (3)
[0055] 得到时间to为故障信号突变点的时间，Wsf (to)为模极大值。
[0056]第三步，得到故障点所处的区间：
[0057]步骤3-1)首先，检测点和六"故障行波波首到达时间，若满足如下条件：
[0059] 初步判断得到故障所在的区间为之间。其中，分别为故障行波波首 到达An-^An的时间，1^、1^分别为输电线路故障测距检测点Aw和AIM端的距离，v为行波 波速。若不满足式(4)则可判断有一故障在检测点n-1和η之间。
[0060] 步骤3-2)利用欧式空间计算位于故障同一侧两个检测点iAk和IAg的电流偏离度dkg 和位于故障点不同侧的两个检测点1?和I Aq的电流偏离度dM (Kk，g<n)。其中lAk，k=l，2, 3···η为各检测点的工频故障电流。设每组检测装置采集的数据长度为L，偏离度计算公式如 下：
[0062]将得到的dkg、dPq进行处理比较，由（5)得到若满足如下关系 dM>Kidkg, (Kk，g或i彡k，g彡η)且(p>i>q或p<i<q) (6)
[0063] 其中L为偏离系数，且KOI，得到区段故障区间位于检测点i和i+1之间。
[0064] 第四步，得到故障行波波头传输特性曲线：
[0065] 步骤4-1)获取故障行波时序。根据附图1，利用输电线路检测装置得到各个检测点 检测到的故障行波到达时间，公式表示如下：
[0067] tzl，tz2，…，tzn为从故障点F1发出的正向行波经折反射后依次到达检测点的时间 与故障发生时刻的差值，tfl，tf2，…，tfn为反向行波经折反射后依次到达的时间与故障发生 时刻的差值。匕，ξ 2，.. .，t为从故障点F2发出的正向行波经折反射后依次到达检测点的 时间与故障发生时刻的差值，^:，<c，_··，'为反向行波经折反射后依次到达的时间与故 障发生时刻的差值。
[0068] 又因为Lx2-Lxi^：0且j，可进一步得到检测点Αι，A2，…，Ai和Aj+i，Aj+2，…，A n测得 的初始故障行波时序：
[0071] 其中，tzk为各检测点检测到的第一个行波到达时间。接着确定两故障点之间 A1+1，…，、行波到达时间，分以下三种情况：当线路发生两处故障区间相同，该测量点行波 到达时间满足式(9);当线路发生两处故障区间相邻，若L 1+1-LX1>LX2-L1+1即tWfG+υ时， 1/ z(i+i)>tf(i+i)。若Li+l_Lxi>Lx2-Li+i即 tzi〉!/ f(j+i)时，1/ z(i+i)>tf(i+i)，根据测量点Ai，Aj+i检测 到的第一个行波时间判断测量点Ai+1的行波时序；当线路发生两处故障区间不相邻，对于 测量点Ak(k= i+1，i+2，…，j)，若Lk-Lxi>Lx2_Lk，1/ zk〈i/ fk。若Lk_Lxi>LX2-Lk，1/ zk>i/ fk。因此， 最终得到整条输电线路每个故障点行波波头到达时间。
[0072] 步骤4-2)根据得到的故障所在的区间及结合步骤4-1)故障行波的时序拟合各时 亥IJ行波波头达到位置距离线路首端Μ和线路尾端N的距离曲线，即得到曲线X M=nu(t)，XM = m2(t)，XN=m(t)，XN=n2(t)。接着统一两曲线的纵坐标，即都以距线路首段Μ的距离为纵坐 标，根据输电线路满足Xm+Xn = L，分别将XN = m(t)，XN = n2(t)转化为曲线L-Xn = f i(t)，L-Xn = f2(t);得到故障行波波头传输特性曲线流程图如附图4。
[0073]第六步，得到故障位置：
[0074]联立XMiimUhL-XNzfKt)，两条曲线的交点所在的位置纵坐标就是故障点距输 电线路首端Μ的距离，所在位置的横坐标是故障时刻。XM = m2(t)，L-XN = f2(t)同理，从而完 成了多点故障定位。
[0075]通过仿真来分析本发明提出的基于分布式的改进故障行波测距算法对于多点故 障的有效性。假设在输电线路上距离Μ端40km和170km处发生A相接地故障，仿真系统图见附 图 5。接着分别设置故障点在 65km、240km;140km、260km;50km、80km;150km、170km;240km、 275km处。并采用传统的故障行波测距方法和本文所提的基于分布式的改进故障行波测距 方法进行仿真实验，测距结果如表1和表2。
[0076] 表1不同故障类型及不同故障距离下的传统故障行波测距结果
[0079] 表2不同故障类型及不同故障距离下的分布式行波测距结果
[0082]根据表1可知，由于传统的故障测距只能在特定的范围内即（LX1<1/3LX2,且L X1彡 1/3LX-1/3LX2)可以正确区分故障点发出的故障行波波首到达时间，在特定范围（L X1<1/ 3LX2，且LX1 < 1/3LX-1/3LX2)之外无论是单端测距还是双端测距均不能准确得到两个故障点 的距离。
[0083]由表2可知，运用基于分布式的改进故障行波测距方法对故障位置、故障类型都有 一定的适应性，并且检测误差基本小于0.3%，满足工程需要。由上可知，本发明的方法能够 有效地计算输电线路多点故障距离位置。
【主权项】
1. 一种考虑多点故障的高压输电线路行波测距方法，其特征在于，包括W下步骤： 1) 利用分布在输电线路上的若干故障检测装置，对电力线路各个检测点的故障电压、 电流信号进行采集处理； 2) 获取检测点信号突变量信息，对保护安装处的电压、电流进行相模解禪变换，然后运 用小波变换获取模极大值的方法得到各个检测点的故障电压、电流突变量信息； 3) 判别故障区间，利用输电线路上故障电流行波的传输特征，即获得的故障电压、电流 突变量信号W及对输电线路各个检测点电流之间的偏离度得到故障点所处的区间； 4) 对故障行波波头进行曲线拟合，利用获取的电压、电流突变量和故障所处的有效区 间，得到故障行波波头传输特性曲线； 5) 确定故障位置，根据行波波头传输特性对故障进行测距，确定故障位置。2. 根据权利要求1所述的考虑多点故障的高压输电线路行波测距方法，其特征在于，在 所述步骤1)中，通过装设在输电线路上的若干Ξ相故障检测装置进行各个检测点的故障电 压、电流信号的采集，Ξ相故障检测装置采用罗氏线圈进行采集，得到故障暂态电压、电流 信号。3. 根据权利要求1所述的考虑多点故障的高压输电线路行波测距方法，其特征在于，在 所述步骤2)中，得到各个检测点的故障电压、电流突变量信息，具体实现步骤如下： 步骤2-1)对于高压输电线路，直接测量得到的各相电流、电压之间存在禪合关系，首先 采用克拉克变换对输电线路电压电流量进行相模变换，使其解禪，得到独立的电压、电流 量，即将保护安装处的电流、电压值乘W克拉克变换逆矩阵得到解禪电压、电流分量，所采 用的变换逆矩阵如下：步骤2-2)对于高压输电线路故障时产生的行波信号，其频率会随着传递时间的变化而 发生变化的特性，采用离散的小波变换对故障行波进行分析，表示如下：式中，b为小波变换平移因子，a为其尺度因子，<55为小波变换基小波，t，a，b取离散值，f (t)为解禪之后故障电流信号； 采用小波变换模极大值方法分析故障信号突变点，将时间to代入式(2)得到小波变换后 的值，并对te(t〇-S，t〇+S)所有变换后的值进行检验，当满足式(3)时： Wsf(t)|^|Wsf(to) (3) 得到时间to为故障信号突变点的时间，Wsf (to)为模极大值。4. 根据权利要求1所述的考虑多点故障的高压输电线路行波测距方法，其特征在于:在 所述步骤3)中，判别故障区间具体步骤如下： 步骤3-1)设线路中设有Ai，A2，A3, 一,Α。个测量点，假设线路发生两处故障时，故障点分 别为F1和F2,其中F1位于测量点Ai和Ai+i间，F2位于测量点Aj和Aj+i间，且j ; 首先检测点An-l和An故障行波波首到达时间，若满足如下条件：则初步判断得到故障所在的区间为Al与An-込间，其中，tn-l、tn分别为故障行波波首到 达An-1和An的时间，Ln-l、Ln分别为输电线路故障测距检测点An-1和An距输电线路首段Μ端的距 离，V为行波波速； 步骤3-2)利用欧式空间偏离度公式计算位于故障同一侧两个检测点lAk和lAg的电流偏 离度dkg和位于故障点不同侧的两个检测点Iap和lAq的电流偏离度(1。9，1《1^，邑《11，其中141<，4 =1，2，3 ···!!为各检测点的工频故障电流； 设每组检测装置采集的数据长度为L，偏离度计算公式如下：将得到的dkg、dpq进行处理比较，由（5)得到:若满足如下关系： dpq>Kldkg，（l《k,g或i《k,g《n)且(p>i>q或p<i<q)(6) 得到区段故障区间位于巧日i+1区间，其中Κι为偏离系数，且Κι〉1。5.根据权利要求4所述的考虑多点故障的高压输电线路行波测距方法，其特征在于:在 所述步骤4)中，拟合故障行波波头传输特性曲线具体步骤如下： 步骤4-1)获取故障行波时序，利用输电线路检测装置得到各个检测点检测到的故障行 波到达时间，公式表示如下：tzl，tz2，-，，tzn为从故障点F1发出的正向行波经折反射后依次到达检测点的时间与故 障发生时刻的差值，tfl，tf2，···，tfn为反向行波经折反射后依次到达的时间与故障发生时刻 的差值，玲* ?,'2，...，4为从故障点F2发出的正向行波经折反射后依次到达检测点的时间 与故障发生时刻的差值，《'/节…，为反向行波经折反射后依次到达的时间与故障 发生时刻的差值； 由于Lx2-Lxi > 0且i《j，可进一步得到检测点Ai，A2，…，Ai和Aj+i，Aj+2，…，An测得的初始 故障行波时序：其中，tzk为各检测点检测到的第一个行波到达时间； 步骤4-2)根据前述得到的故障所在的区间W及步骤4-1)故障行波的时序拟合各时刻 行波波头达到位置距离线路首端Μ和线路尾端N的距离曲线，即得到曲线XM=mi(t)，XM = m2 (t) ,XN=ni(t) ,Χν = Π 2(?)； 接着统一两曲线的纵坐标，即都W距线路首段Μ的距离为纵坐标，根据输电线路满足Χμ+ Xn=L，分别将 XN = ni(t)，XN=n2(t)转化为曲线kXN = fi(t)，kXN = f2(t)。6.根据权利要求4所述的考虑多点故障的高压输电线路行波测距方法，其特征在于:在 所述步骤5)中，确定故障位置具体步骤如下： 联立XM=mi(t)，kXN = fi(t)，两条曲线的交点所在的位置纵坐标就是故障点距输电线 路首端Μ的距离，所在位置的横坐标是故障时刻； 联立XM = m2(t) ,L-XN = f2(t)，同理可完成多点故障定位。
【文档编号】G01R31/11GK106093702SQ201610394721
【公开日】2016年11月9日
【申请日】2016年6月6日 公开号201610394721.4, CN 106093702 A, CN 106093702A, CN 201610394721, CN-A-106093702, CN106093702 A, CN106093702A, CN201610394721, CN201610394721.4
【发明人】卜京, 孙震宇, 马迪, 刘馨月, 王俊, 吉晓筱, 杨阳, 张梦月, 殷明慧, 谢云云
【申请人】南京理工大学