CN105445618B - 一种小电流接地系统故障选线方法及装置 - Google Patents
一种小电流接地系统故障选线方法及装置 Download PDFInfo
- Publication number
- CN105445618B CN105445618B CN201510810748.2A CN201510810748A CN105445618B CN 105445618 B CN105445618 B CN 105445618B CN 201510810748 A CN201510810748 A CN 201510810748A CN 105445618 B CN105445618 B CN 105445618B
- Authority
- CN
- China
- Prior art keywords
- line
- fault
- calculating
- atoms
- singular
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000007935 neutral effect Effects 0.000 title abstract description 4
- 230000001052 transient effect Effects 0.000 claims abstract description 102
- 238000005070 sampling Methods 0.000 claims description 65
- 238000004364 calculation method Methods 0.000 claims description 49
- 238000010606 normalization Methods 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 23
- 238000010187 selection method Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 12
- 230000001629 suppression Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 2
- 238000012986 modification Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 238000000354 decomposition reaction Methods 0.000 description 18
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000010219 correlation analysis Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Emergency Protection Circuit Devices (AREA)
- Locating Faults (AREA)
Abstract
本发明提供一种小电流接地系统故障选线方法及装置，所述方法包括：根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；根据各条线路的暂态零序电流，计算各条线路的最佳原子；根据各条线路的最佳原子，计算各条线路的原子奇异熵；根据各条线路的最佳原子，计算各条线路的原子实际能量；根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；将各条线路的故障综合测度值进行比较，进行故障选线。该方法提高了故障选线的可靠性和准确性，能满足各种复杂电网结构的故障工况要求。
Description
技术领域
本发明涉及电力系统故障选线技术领域，具体设计一种小电流接地系统故障选线方法及装置。
背景技术
我国配电网主要采用中性点非有效接地系统，即小电流接地系统，小电流接地故障又称为单向接地故障，是配电网最主要的故障形式，随着用户对供电可靠性的要求越来越高，学者对小电流接地故障选线进行了大量的研究，因为在发生单向接地故障时，暂态信号特征量非常丰富，使得基于暂态信号的分析方法得到了广泛的关注，常用的暂态信号故障分析方法有：暂态能量法、普罗尼算法、相关分析法、小波分解法、S变换法以及经验模态分解法等。
但是，现有的暂态信号故障分析方法对故障信号进行处理时仍存在许多缺陷，比如，利用暂态能量法进行选线时，当发生大电阻接地故障时，由于消弧线圈感性电流对故障瞬间零序电流的补偿作用，使得线路间的暂态零序电流相差较小，从而导致误判；普罗尼算法虽然对低频暂态信号的拟合能力较好，但是对高频暂态信号的拟合能力欠佳；相关分析法在流过电缆健全线路与架空故障线路放入零序电流的幅度相近时失效；小波变换法虽然对故障信号处理具有良好的时频特性，但在发生高阻接地故障时，故障分量较小，难以判断故障线路；S变换法虽然具有良好的时频特性，但是变换分解后的信息量太多，无法合理的利用相角信息；经验模态分解法虽然对非线性、非平稳信号的处理具有优势，但是会造成模态混叠现象。
因此，现有的小电流接地系统故障选线方法具有以下缺陷：无法对暂态信号蕴含丰富的故障信息进行量化处理，无法合理有效的利用故障特征信息，导致选线结果不准确；对故障特征量的提取不具有自 适用性，导致出现了无法解释的分解项；仅是将单一的故障暂态信息作为进行选线的判断依据，没有将多种故障暂态信息进行融合，进行综合选线，因此，不能覆盖所有的接地工况，故障选线准确性不高。所以，现有的小电流接地系统故障选线装置的选线结果准确度较低，而且难以适用各种复杂电网结构的故障工况要求。
发明内容
针对现有技术存在的缺陷，本发明提出一种小电流接地系统故障选线方法及装置，以解决现有技术存在的选线结果准确性较低，而且难以适用各种复杂电网结构的故障工况要求的问题。
为此目的，第一方面，本发明提供一种小电流接地系统故障选线方法，所述方法包括：
根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；
根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数；
根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵；
根据各条线路的Q个最佳原子，计算各条线路的特征原子的实际能量；
根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；
根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；
将各条线路的故障综合测度值进行比较，选取故障综合测度值最大的线路为故障线路。
其中，所述根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵，包括：
计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇 异值；
对各条线路的Q个奇异值进行归一化处理，得到各条线路归一化后的Q个奇异值；
根据各条线路归一化后的Q个奇异值，计算各条线路的原子奇异熵。
其中，在所述根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子之后，所述方法还包括：
根据预设的采样工频周期对各条线路的Q个最佳原子进行采样，得到各条线路的Q个最佳采样原子；
相应地，所述计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值，包括：
计算各条线路的Q个最佳采样原子的奇异值，得到各条线路的Q个奇异值。
其中，在所述根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵之后，所述方法还包括：
对各条线路的原子奇异熵进行归一化处理，得到各条线路归一化后的原子奇异熵；
相应地，根据各条线路归一化后的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值。
其中，所述根据各条线路的Q个最佳原子，计算各条线路的特征原子的实际能量，包括：
对各条线路的Q个最佳原子求和，得到各条线路的特征原子；
根据各条线路的特征原子，基于预设的特征原子能量计算公式，计算各条线路的特征原子能量；
根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量；
根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量。
其中，所述各条线路的故障综合测度函数为：
Xj＝EjCj，
其中，Xj为各条线路的故障综合测度函数，Ej为各线路的特征原子的实际能量，Cj为各条线路的故障信任度，j为小电流接地系统中线路的编号。
其中，在所述根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流之前，所述方法还包括：
获取小电流接地系统当前的零序电压和母线额定电压；
将所述获取零序电压和所述母线额定电压进行比较，若零序电压小于或等于预设的母线额定电压值，执行所述获取小电流接地系统当前的零序电压和母线额定电压的步骤；
若零序电压大于预设的母线额定电压值，检测小电流接地系统中电压互感器的接通状态，若所述接通状态为断线，发出电压互感器断线警告信息；
若所述接通状态正常，检测小电流接地系统中消弧线圈的当前状态，若所述当前状态为串联谐振，调节消弧线圈远离谐振点，若所述当前状态不是串联谐振，采集故障时刻起小电流接地系统中各条线路的暂态零序电流。
第二方面，本发明提供一种小电流接地系统故障选线装置，所述装置包括：
采集模块，用于根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；
第一计算模块，用于根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数；
采样模块，用于根据预设的采样工频周期对各条线路的Q个最佳原子进行采样，得到各条线路的Q个最佳采样原子；
第二计算模块，用于根据各条线路的Q个最佳采样原子，计算各 条线路的原子奇异熵；
第三计算模块，用于根据各条线路的Q个最佳采样原子，计算各条线路的特征原子的实际能量；
第四计算模块，用于根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；
修正模块，用于根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；
选取模块，用于将各条线路的故障综合测度值进行比较，选取故障综合测度值最大的线路为故障线路。
其中，所述第二计算模块，包括：
奇异值计算单元，用于计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值；
归一化单元，用于对各条线路的Q个奇异值进行归一化处理，得到各条线路归一化后的Q个奇异值；
原子奇异熵计算单元，用于根据各条线路归一化后的Q个奇异值，计算各条线路的原子奇异熵。
其中，所述修正模块，包括：
特征原子计算单元，用于对各条线路的Q个最佳原子求和，得到各条线路的特征原子；
特征原子能量计算单元，用于根据各条线路的特征原子，基于预设的特征原子能量计算公式，计算各条线路的特征原子能量；
电流能量计算单元，用于根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量；
实际能量计算单元，用于根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量。
本发明提供的一种小电流接地系统故障选线方法及装置，具有如下优点：通过用小电流接地系统中各条线路的原子奇异熵来表征各条线路的暂态故障信息，量化出各条线路的原子包含的故障暂态信息的 复杂程度，提高选线准确性；采用匹配追踪算法对暂态零序电流分解得到的最佳原子，不仅表征了故障暂态信息，还降低了故障信息维度，而多个最佳原子组成的序列能从时频域表征暂态零序电流的故障特征和能量分布；利用多个最佳原子重构特征原子，将不同的最佳原子的故障暂态信息进行融合，涵盖了多种故障暂态信息，利用特征原子的实际能量来修正故障信任度值，得到故障综合测度函数不仅保留了各条线路的原始暂态故障信息，还涵盖了特征原子的暂态故障信息，提高了小电流接地系统故障选线的可靠性和准确性，能满足各种复杂电网结构的故障工况要求。
附图说明
为了更清楚地说明本公开实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些图获得其他的附图。
图1为本发明一实施例提供的一种小电流接地系统故障选线的方法流程图；
图2为本发明一实施例提供的辐射状缆线混合配电网络结构示意图；
图3为本发明一实施例提供的一种小电流接地系统故障选线方法总体流程图；
图4为本发明一实施例提供的小电流接地系统故障瞬间，各条线路2个工频周期内的暂态零序电流示意图；
图5为本发明一实施例提供的线路1的4个最佳原子波形图；
图6为本发明一实施例提供的线路2的4个最佳原子波形图；
图7为本发明一实施例提供的线路3的4个最佳原子波形图；
图8为本发明一实施例提供的线路4的4个最佳原子波形图；
图9为本发明一实施例提供的小电流接地系统中各条线路的特征原子波形图；
图10为本发明一实施例提供的一种小电流接地系统故障选线装置的结构示意图。
具体实施方式
下面将结合本公开实施例中的附图，对本公开实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。
如图1所示，本公开一实施例提供了一种小电流接地系统故障选线方法，该方法包括如下步骤S1至S7：
S1、根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；
以采样频率为100kHz为例：根据100kHz的采样频率确定采集小电流接地系统故障时刻起2个工频周期内的各条线路的暂态零序电流。
应需说明的是，本实施例中的采样频率仅为举例说明，本实施例不限定采样频率的具体值，本领域技术人员可根据实际情况设定采样频率。
具体地，根据预设采样点对每个工频周期内的各条线路的暂态零序电流进行采集，得到各条线路的暂态零序电流为Ij(n)，其中，n为采样点，n＝1,2,…,N，j为各条线路的编号，所述预设采样点由预设的采样频率确定。
S2、根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数；
S3、根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵；
S4、根据各条线路的Q个最佳原子，计算各条线路的特征原子的 实际能量；
S5、根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；
S6、根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；
S7、将各条线路的故障综合测度值进行比较，选取故障综合测度值最大的线路为故障线路。
本实施例公开的一种小电流接地系统故障选线方法，通过预先构建的故障信任度函数和小电流接地系统中各条线路的原子奇异熵计算各条线路的故障信任度值，量化暂态故障信息的复杂程度，并通过特征原子的实际能量对故障信任度值进行修正，将多种暂态故障信息进行融合，得到故障综合测度值，通过根据故障信任度值进行故障选线，提高了选线结果的可靠性准确性。
在本实施例中，步骤S2：“根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数”，包括图中未示出的如下步骤S2`：
S2`、采用匹配追踪算法对各条线路的暂态零序电流进行Q次迭代分解产生Q个最佳原子波形，其中，Q为预设常数。
具体地，本实施例中的步骤S2`中，迭代分解过程如下：
采用匹配追踪算法对最佳原子波形信号I(n)进行分析时，若迭代分解次数为m，最佳原子波形信号I(n)分解为：
其中，为第k次迭代时求得的最佳匹配原子，Rkf(n)为当前信号，表示内积运算，RmI(n)为残余信号。
同时，最佳原子波形信号I(n)也可以分解为||I(n)||2，公式为：
其中，满足如下公式：
其中，0≤α≤1，如果分解m次后达到设定的精确度，则停止分解。
通过采用匹配追踪算法对小电流接地系统中的各条线路的暂态零序电流进行迭代分解，具有自适用性，使得选取的原子能最佳的匹配各条线路的原始暂态故障信号的结构特征。
应需说明的是，本实施例中可以通过设定精确度来确定迭代分解的次数，也可以直接设定迭代分解次数。
在本实施例中，步骤S3：“根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵”，包括图1中未示出的如下细分步骤S31至S33：
S31、计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值；
S32、对各条线路的Q个奇异值进行归一化处理，得到各条线路归一化后的Q个奇异值；
S33、根据各条线路归一化后的Q个奇异值，计算各条线路的原子奇异熵。
具体地，本实施例的步骤S31中，计算奇异值的过程如下：
根据小电流接地系统中各条线路的最佳原子波形的系数构成c×h阶矩阵A；
根据分别为c×c和h×h阶的正交矩阵U和V对所述c×h阶矩阵A进行奇异值分解，公式为：
A＝UΛVT，
其中，Λ＝diag(λ1,λ2,…,λp)，p＝min(c,h)，且存在λ1>>λ2>…>λp，λ1,λ2,…,λp即为A的奇异值。
具体地，在故障选线中，由于线路类型、故障电阻、故障位置、故障初相角以及现场噪声强度的影响，使得根据暂态零序电流的时频状态不能准确的识别故障线路，通过采用奇异值分解理论将各条线路的暂态零序电流分解成最佳原子序列矩阵，得到一组能反映原子序列 矩阵基本特征的奇异值，并对各条线路的奇异值组进行归一化处理，简单、准确的定量结果，使得奇异值受环境的影响减小。
在本实施例的步骤S32中，奇异值归一化公式为：
其中，ujm为线路j的第m个归一化后的奇异值，λjm为线路j的第m个奇异值，J为小电流接地系统中所有的线路数，m为迭代次数，也即最佳原子个数，j＝1,2,…,J，m＝1,2,3,4。
在本实施例的步骤S33中，原子奇异熵计算公式为：
其中，Sj小电流接地系统中线路j的奇异熵。
在本实施例中，在步骤S3：“根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵”之后，所述方法还包括图1中未示出的如下步骤S3`：
S3`、对各条线路的原子奇异熵进行归一化处理，得到各条线路归一化后的原子奇异熵；
相应地，步骤S5包括：根据各条线路归一化后的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值。
具体地，以对各条线路的4个原子奇异熵进行归一化处理为例，原子奇异熵归一化公式为：
其中，j为小电流接地系统中各条线路的编号，J为小电流接地系统中所有的线路数，ej为线路j归一化后的奇异熵，j＝1,2,…,J。
通过对小电流接地系统中各条线路的原子奇异熵进行归一化处 理，对各条线路的原子奇异熵进行简单、准确的定量，防止原子奇异熵受小电流接地系统环境的影响，使得原子奇异熵值发生较大波动。
在本实施例中，所述预先构建的故障信任度函数为：
其中，C(ej)为故障信任度函数，ej为各条线路的奇异熵，j为小电流接地系统的线路的编号，a、b均为预设常数，0≤k1≤k2≤1。
具体地，在本实施例中，故障信任度函数可为：
应需说明的是，本实施例的故障信任度函数仅为举例说明，本实施例不限定a、b、k1、k2的具体值，本领域技术人员可根据实际情况通过多次实验验证得到a、b、k1、k2的具体值。
在本实施例中，在步骤S2：“根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子”之后，还包括图1中未示出的如下步骤S21：
S21、根据预设的采样工频周期对各条线路的Q个最佳原子进行采样，得到各条线路的Q个最佳采样原子；
可选地，设定采样工频周期为1/4工频周期。应需说明的是，本实施例中的采样工频周期仅为举例说明，本实施例不限定采样工频周期的具体值，本领域技术人员可根据实际情况设定采样工频周期。
相应地，步骤S31：“计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值”，包括：计算各条线路的Q个最佳采样原子的奇异值，得到各条线路的Q个奇异值。
在本实施例的步骤S4：“根据各条线路的Q个最佳原子，计算各条线路的特征原子的实际能量”，包括图1中未示出的如下步骤S41至 S44：
S41、对各条线路的Q个最佳原子求和，得到各条线路的特征原子；
具体地，特征原子计算公式为：
Yj(n)＝Rj1(n)+Rj2(n)+···+Rjm(n)+···+RjQ(n)，
其中，j为小电流接地系统中各条线路的编号，m为迭代分解次数，Rjm(n)为线路j的第m次迭代分解的最佳原子。
S42、根据各条线路的特征原子，基于预设的特征原子能量计算公式，计算各条线路的特征原子能量；
具体地，所述预设的特征原子能量计算公式为：
其中，Hj为线路j的特征原子能量，n＝1,2,…,N，N为根据预设的采样频率确定的采样点数。
S43、根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量；
具体地，所述预设的暂态零序电流能量计算公式为：
其中，Zj为线路j的暂态零序电流能量，n＝1,2,…,N，N为根据预设的采样频率确定的采样点数。
S44、根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量。
具体地，特征原子的实际能量计算公式为：
Ej＝Zj·Hj，
其中，Hj和Zj分别为线路j的特征原子能量和暂态零序电流能量， Ej为线路j的特征原子的实际能量，j＝1,2,…,J，J为小电流接地系统中所有的线路数。
在本实施例中，在步骤S6中，所述各条线路的故障综合测度值计算公式为：
Xj＝EjCj，
其中，Xj为各条线路的故障综合测度值，Ej为各线路的特征原子的实际能量，Cj为各条线路的故障信任度值。
本实施例公开的一种小电流接地系统故障选线方法，通过根据小电流系统中各条线路的特征原子的实际能量来修正预先构建故障信任度值，得到故障综合测度值，既保留了故障信任度值中的暂态故障信息，同时，还涵盖了线路特征原子的实际能量，将基于故障信任度函数中的原始暂态故障信息选线的优势和基于线路特征原子的实际能量选线的优势进行互补，提高了小电流接地系统故障选线的可靠性和准确性。
在本实施例中，在步骤S1：“根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流”之前，所述方法还包括图1中未示出的如下步骤S01至S06：
S01、获取小电流接地系统当前的零序电压和母线额定电压；
S02、将所述获取零序电压和所述母线额定电压进行比较，若零序电压小于或等于预设的母线额定电压值，执行步骤S01；若零序电压大于预设的母线额定电压值，则执行步骤S03；
S03、检测小电流接地系统中电压互感器的接通状态，若所述接通状态为断线，则执行步骤S04；若所述接通状态正常，则执行步骤S05；
S04、发出电压互感器断线警告信息；
S05、检测小电流接地系统中消弧线圈的当前状态，若所述当前状态为串联谐振，则执行步骤S06；若所述当前状态不是串联谐振，则执行步骤S1；
S06、调节消弧线圈远离谐振点。
具体地，如图2所示，以辐射状缆线混合配电网络为例，其中，架空线路1为线路1，长度为13.5km，架空线路2为线路2，长度为24km，缆线混合线路3为线路3，其中，电缆线长度为5km，架空线路长度为12km；电缆线路4为线路4，线路正序参数为R1＝0.17Ω/km，L1＝1.2mH/km，C1＝9.697nF/km，零序参数为R0＝0.23Ω/km，L0＝5.48mH/km，C0＝6nF/km；线路4为电缆线路，长度10km，线路正序参数为R11＝0.193Ω/km，L11＝0.442mH/km，C11＝143nF/km，零序参数为R00＝1.93Ω/km，L00＝5.48mH/km，C00＝143nF/km，消弧线圈Lx的过补偿度为10％，消弧线圈的电感经计算为1.574H，采样频率100kHz，变压器接线方式为Δ/Y，变比为220kV/35kV。上述小电流接地系统故障选线流程描述如下：
图3示出了上述小电流接地系统故障选线方法的总体流程示意图：
故障选线启动判断阶段，包括步骤C1至C6：
C1、获取小电流接地系统当前的零序电压U0(t)和母线额定电压U；
C2、将所述获取零序电压U0(t)和所述母线额定电压U进行比较，当U0(t)≤0.15U时，执行步骤C1，当U0(t)>0.15U时，则执行步骤C3；
C3、检测小电流接地系统中电压互感器的接通状态，若所述接通状态为断线，则执行步骤C4；若所述接通状态正常，则执行步骤C5；
C4、发出电压互感器断线警告信息；
C5、检测小电流接地系统中消弧线圈的当前状态，若所述当前状态为串联谐振，则执行步骤C6；若所述当前状态不是串联谐振，则执行步骤S1；
C6、调节消弧线圈远离谐振点。
小电流接地系统故障选线阶段，包括步骤D1至D14：
D1、根据100kHz的采样频率分别采集故障时刻起2个工频周期内小电流接地系统中各条线路的暂态零序电流：I1(n),I2(n),I3(n),I4(n)，其中，n为采样点，n＝1,2,…,N，N为2个工频周期内的总采样点数；
图4示出了小电流接地系统发生故障瞬间，各条线路在2个工频周期内的暂态零序电流波形示意图，其中，图4(a)为分支线路1的暂态零序电流波形，图4(b)为分支线路2的暂态零序电流波形，图4(c)为分支线路3的暂态零序电流波形，图4(d)为分支线路4的暂态零序电流波形。
D2、采用匹配追踪算法对各条线路的暂态零序电流进行4次迭代分解产生4个最佳原子波形：R1m(n)，R2m(n)，R3m(n)，R4m(n)，其中，R1m(n)为第一条线路经m次迭代分解产生的m个最佳原子波形，其中，R1m(n)＝[R11(n)，R12(n)，R13(n)，R14(n)]，其中，m＝1,2,3,4；
图5示出了小电流接地系统发生故障瞬间，采用匹配追踪算法4次迭代分解线路1的暂态零序电流得到的4个最佳原子波形，其中，图5(a)最佳原子1的波形，图5(b)为最佳原子2的波形，图5(c)为最佳原子3的波形，图5(d)为最佳原子4的波形。
图6示出了小电流接地系统发生故障瞬间，采用匹配追踪算法4次迭代分解线路2的暂态零序电流得到的4个最佳原子波形，其中，图6(a)最佳原子1的波形，图6(b)为最佳原子2的波形，图6(c)为最佳原子3的波形，图6(d)为最佳原子4的波形。
图7示出了小电流接地系统发生故障瞬间，采用匹配追踪算法4次迭代分解线路3的暂态零序电流得到的4个最佳原子波形，其中，图7(a)最佳原子1的波形，图7(b)为最佳原子2的波形，图7(c)为最佳原子3的波形，图7(d)为最佳原子4的波形。
图8示出了小电流接地系统发生故障瞬间，采用匹配追踪算法4次迭代分解线路4的暂态零序电流得到的4个最佳原子波形，其中，图8(a)最佳原子1的波形，图8(b)为最佳原子2的波形，图8(c)为最佳原子3的波形，图8(d)为最佳原子4的波形。
D3、根据1/4工频周期的采样工频周期对各条线路的4个最佳原子进行采样，得到各条线路的4个最佳采样原子序列：R`1m(n)，R`2m(n)， R`3m(n)，R`4m(n)，其中，R`1m(n)＝[R`11(n)，R`12(n)，R`13(n)，R`14(n)]， 表示不超过的最大整数；
D4、计算各条线路的4个最佳采样原子序列的奇异值，得到各条线路的奇异值组：λ1m,λ2m,λ3m,λ4m，其中，λ1m为线路1的奇异值组，其中，λ1m＝[λ11,λ12,λ13,λ14]，λ11，为线路1的最佳采样原子1的奇异值；
得到4条线路的最佳采样原子的奇异值如表1所示：
表1
D5、对各条线路的奇异值组进行归一化处理，得到各条线路归一化后的奇异值组：u1m、u2m、u3m、u4m，其中，u1m为线路1的归一化后的奇异值组，u1m＝[u11,u12,u13,u14]，u11为线路1的最佳采样原子1的归一化的奇异值，其中，奇异值归一化公式为：
其中，ujm为线路j的第m个归一化后的奇异值，j＝1,2,3,4，m＝1,2,3,4。
D6、根据各条线路归一化后的奇异值组：u1m、u2m、u3m、u4m，计算各条线路的原子奇异熵：S1,S2,S3,S4，其中，原子奇异熵的计算公式为：
其中，Sj为线路j的原子奇异熵，j＝1,2,3,4，m＝1,2,3,4。
计算得到4条线路的原子奇异熵为：Sj＝[0.3884 1.023 0.8152 0.7439]。
D7、对各条线路的原子奇异熵：S1,S2,S3,S4，进行归一化处理，得到各条线路归一化后的奇异熵：e1,e2,e3,e4，其中，原子奇异熵归一化公式为：
其中，ej为线路j的归一化后的原子奇异熵。
计算得到4条线路的归一化后的原子奇异熵为：ej＝[0.4241 0.161 0.19350.2214]。
D8、根据各条线路归一化后的原子奇异熵：e1,e2,e3,e4，基于预先构建的故障信任度函数计算各条线路的故障信任度值：C(e1)，C(e2)，C(e3)，C(e4)，其中，故障信任度函数公式为：
计算得到4条线路的故障信任度值为：Cj＝[1 0.4442 0.5741 0.6857]。
D9、对各条线路的4个最佳采样原子求和，得到各条线路的特征原子，其中，特征原子的计算公式为：
Yj(n)＝R`j1(n)+R`j2(n)+R`j3(n)+R`j4(n)，
其中，Yj(n)为线路j的特征原子，R`j1(n)为线路j的最佳采样原子1。
在本实施例中，各条线路的特征原子波形如图9所示。
D10、根据各条线路的特征原子：Y1(n),Y2(n),Y3(n),,Y4(n)，基于预 设的特征原子能量计算公式，计算各条线路的特征原子能量：H1,H2,H3,H4，其中，特征原子能量计算公式为：
其中，Hj为线路j的特征原子能量，
D11、根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量：Z1,Z2,Z3,Z4，其中，暂态零序电流能量计算公式为：
其中，Zj为线路j的1/4采样工频周期内的暂态零序电流的能量，
D12、根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量：E1,E2,E3,E4，公式为：
Ej＝Zj·Hj，
其中，Hj和Zj分别为线路j的特征原子能量和暂态零序电流能量，Ej为线路j的特征原子的实际能量，j＝1,2,3,4。
计算得到4条线路的特征原子的实际能量为：Ej＝[1648.7 84.2 460.7 909.8]。
D13、根据各条线路的特征原子的实际能量：E1,E2,E3,E4，修正所述预先构建的故障信任度值，得到各条线路的故障综合测度值为：
Xj＝EjCj。
计算得到4条线路的故障综合测度值为：Xj＝[1648.7 37.383 264.45 623.84]。
D14、将各条线路的故障综合测度值进行比较，选取故障综合测度 值最大的线路为故障线路。
由于最大故障综合测度值Xmax＝1648.7，因此，线路1被选为故障线路，选线过程结束。
如图10所示，本实施例提供一种小电流接地系统故障选线装置，可包括以下单元：采集模块11、第一计算模块12、采样模块13、第二计算模块14、第三计算模块15、第四计算模块16、修正模块17以及选取模块18。
采集模块11，用于根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；
第一计算模块12，用于根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数；
采样模块13，用于根据预设的采样工频周期对各条线路的Q个最佳原子进行采样，得到各条线路的Q个最佳采样原子；
第二计算模块14，用于根据各条线路的Q个最佳采样原子，计算各条线路的原子奇异熵；
第三计算模块15，用于根据各条线路的Q个最佳采样原子，计算各条线路的特征原子的实际能量；
第四计算模块16，用于根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；
修正模块17，用于根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；
选取模块18，用于将各条线路的故障综合测度值进行比较，选取故障综合测度值最大的线路为故障线路。
在本实施例中，第二计算模块14，包括图中未示出的：奇异值计算单元141、归一化单元142以及原子奇异熵计算单元143；
奇异值计算单元141，用于计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值；
归一化单元142，用于对各条线路的Q个奇异值进行归一化处理，得到各条线路归一化后的Q个奇异值；
原子奇异熵计算单元143，用于根据各条线路归一化后的Q个奇异值，计算各条线路的原子奇异熵。
在本实施例中，修正模块17，包括图中未示出的：特征原子计算单元171、特征原子能量计算单元172、电流能量计算单元173以及实际能量计算单元174；
特征原子计算单元171，用于对各条线路的Q个最佳原子求和，得到各条线路的特征原子；
特征原子能量计算单元172，用于根据各条线路的特征原子，基于预设的特征原子能量计算公式，计算各条线路的特征原子能量；
电流能量计算单173元，用于根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量；
实际能量计算单元174，用于根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量。
本实施例提供的小电流接地系统故障选线装置可执行上述小电流接地系统故障选线方法实施例的步骤，本实施例不再赘述。
本实施例公开的一种小电流接地系统故障选线装置，通过用小电流接地系统中各条线路的原子奇异熵来表征各条线路的暂态故障信息，量化出各条线路的原子包含的故障暂态信息的复杂程度，提高选线准确性；而且多个最佳原子组成的序列能从时频域表征暂态零序电流的故障特征和能量分布；利用多个最佳原子重构特征原子，将不同的最佳原子的故障暂态信息进行融合，涵盖了多种故障暂态信息，利用特征原子的实际能量来修正故障信任度值，得到故障综合测度函数不仅保留了各条线路的原始暂态故障信息，还涵盖了特征原子的暂态故障信息，提高了小电流接地系统故障选线的可靠性和准确性，能满足各种复杂电网结构的故障工况要求。
需要说明的是，本文中“第一”、“第二”、“第三”和“第四”仅 仅用来区分名称相同的实体或操作，并不暗示这些实体或操作之间顺序或关系。
本领域普通技术人员可以理解：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (10)
1.一种小电流接地系统故障选线方法，其特征在于，所述方法包括：
根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；
根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数；
根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵；
根据各条线路的Q个最佳原子，计算各条线路的特征原子的实际能量；
根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；
根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；
将各条线路的故障综合测度值进行比较，选取故障综合测度值最大的线路为故障线路。
2.根据权利要求1所述的方法，其特征在于，所述根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵，包括：
计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值；
对各条线路的Q个奇异值进行归一化处理，得到各条线路归一化后的Q个奇异值；
根据各条线路归一化后的Q个奇异值，计算各条线路的原子奇异熵。
3.根据权利要求2所述方法，其特征在于，在所述根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子之后，所述方法还包括：
根据预设的采样工频周期对各条线路的Q个最佳原子进行采样，得到各条线路的Q个最佳采样原子；
相应地，所述计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值，包括：
计算各条线路的Q个最佳采样原子的奇异值，得到各条线路的Q个奇异值。
4.根据权利要求1所述的方法，其特征在于，在所述根据各条线路的Q个最佳原子，计算各条线路的原子奇异熵之后，所述方法还包括：
对各条线路的原子奇异熵进行归一化处理，得到各条线路归一化后的原子奇异熵；
相应地，根据各条线路归一化后的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值。
5.根据权利要求2所述的方法，其特征在于，所述根据各条线路的Q个最佳原子，计算各条线路的特征原子的实际能量，包括：
对各条线路的Q个最佳原子求和，得到各条线路的特征原子；
根据各条线路的特征原子，基于预设的特征原子能量计算公式，计算各条线路的特征原子能量；
根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量；
根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量。
6.根据权利要求1所述的方法，其特征在于，各条线路的故障综合测度函数为：
Xj＝EjCj，
其中，Xj为各条线路的故障综合测度函数，Ej为各线路的特征原子的实际能量，Cj为各条线路的故障信任度，j为小电流接地系统中线路的编号。
7.根据权利要求1所述的方法，其特征在于，在所述根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流之前，所述方法还包括：
获取小电流接地系统当前的零序电压和母线额定电压；
将所述获取的零序电压和所述母线额定电压进行比较，若零序电压小于或等于预设的母线额定电压值，执行所述获取小电流接地系统当前的零序电压和母线额定电压的步骤；
若零序电压大于预设的母线额定电压值，检测小电流接地系统中电压互感器的接通状态，若所述接通状态为断线，发出电压互感器断线警告信息；
若所述接通状态正常，检测小电流接地系统中消弧线圈的当前状态，若所述当前状态为串联谐振，调节消弧线圈远离谐振点，若所述当前状态不是串联谐振，采集故障时刻起小电流接地系统中各条线路的暂态零序电流。
8.一种小电流接地系统故障选线装置，其特征在于，所述装置包括：
采集模块，用于根据预设的采样频率分别采集故障时刻起小电流接地系统中各条线路的暂态零序电流；
第一计算模块，用于根据各条线路的暂态零序电流，计算各条线路的Q个最佳原子，其中，Q为预设常数；
采样模块，用于根据预设的采样工频周期对各条线路的Q个最佳原子进行采样，得到各条线路的Q个最佳采样原子；
第二计算模块，用于根据各条线路的Q个最佳采样原子，计算各条线路的原子奇异熵；
第三计算模块，用于根据各条线路的Q个最佳原子，计算各条线路的特征原子的实际能量；
第四计算模块，用于根据各条线路的原子奇异熵，基于预先构建的故障信任度函数计算各条线路的故障信任度值；
修正模块，用于根据各条线路的原子实际能量，修正所述各条线路的故障信任度值，得到各条线路的故障综合测度值；
选取模块，用于将各条线路的故障综合测度值进行比较，选取故障综合测度值最大的线路为故障线路。
9.根据权利要求8所述的装置，其特征在于，所述第二计算模块，包括：
奇异值计算单元，用于计算各条线路的Q个最佳原子的奇异值，得到各条线路的Q个奇异值；
归一化单元，用于对各条线路的Q个奇异值进行归一化处理，得到各条线路归一化后的Q个奇异值；
原子奇异熵计算单元，用于根据各条线路归一化后的Q个奇异值，计算各条线路的原子奇异熵。
10.根据权利要求8所述的装置，其特征在于，所述修正模块，包括：
特征原子计算单元，用于对各条线路的Q个最佳原子求和，得到各条线路的特征原子；
特征原子能量计算单元，用于根据各条线路的特征原子，基于预设的特征原子能量计算公式，计算各条线路的特征原子能量；
电流能量计算单元，用于根据预设的暂态零序电流能量计算公式，计算各条线路的暂态零序电流能量；
实际能量计算单元，用于根据各条线路的特征原子能量和各条线路的暂态零序电流能量，计算各条线路的特征原子的实际能量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510810748.2A CN105445618B (zh) | 2015-11-20 | 2015-11-20 | 一种小电流接地系统故障选线方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510810748.2A CN105445618B (zh) | 2015-11-20 | 2015-11-20 | 一种小电流接地系统故障选线方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105445618A CN105445618A (zh) | 2016-03-30 |
| CN105445618B true CN105445618B (zh) | 2018-06-05 |
Family
ID=55556045
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510810748.2A Expired - Fee Related CN105445618B (zh) | 2015-11-20 | 2015-11-20 | 一种小电流接地系统故障选线方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105445618B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106501678B (zh) * | 2016-11-02 | 2019-06-14 | 李晓明 | 一种小电流接地故障选线方法及系统 |
| CN106526429B (zh) * | 2016-12-06 | 2019-08-13 | 李晓明 | 一种具有纠错功能的接地故障选线方法 |
| CN107045093B (zh) * | 2017-03-16 | 2019-09-27 | 西安理工大学 | 基于快速s变换的小电流单相接地故障选线方法 |
| CN112415324A (zh) * | 2020-10-15 | 2021-02-26 | 西安理工大学 | 一种小电流接地系统故障选线方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854437A (zh) * | 2012-08-22 | 2013-01-02 | 广东电网公司电力科学研究院 | 应用时频原子分解理论的小电流接地系统故障选线方法 |
| CN103499769A (zh) * | 2013-09-23 | 2014-01-08 | 武汉大学 | 一种谐振接地系统单相接地故障自适应选线方法 |
| US8744637B1 (en) * | 2013-09-12 | 2014-06-03 | Deepak Maragal | Methods and apparatus for detecting the clearance of fault in shunt reactor compensated transmission lines |
| CN104035002A (zh) * | 2014-04-16 | 2014-09-10 | 武汉大学 | 一种利用行波固有频率和原子分解能量熵的故障选相法 |
| CN104297635A (zh) * | 2014-10-14 | 2015-01-21 | 河南理工大学 | 基于原子稀疏分解与极限学习机的配电网故障选线方法 |
-
2015
- 2015-11-20 CN CN201510810748.2A patent/CN105445618B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854437A (zh) * | 2012-08-22 | 2013-01-02 | 广东电网公司电力科学研究院 | 应用时频原子分解理论的小电流接地系统故障选线方法 |
| US8744637B1 (en) * | 2013-09-12 | 2014-06-03 | Deepak Maragal | Methods and apparatus for detecting the clearance of fault in shunt reactor compensated transmission lines |
| CN103499769A (zh) * | 2013-09-23 | 2014-01-08 | 武汉大学 | 一种谐振接地系统单相接地故障自适应选线方法 |
| CN104035002A (zh) * | 2014-04-16 | 2014-09-10 | 武汉大学 | 一种利用行波固有频率和原子分解能量熵的故障选相法 |
| CN104297635A (zh) * | 2014-10-14 | 2015-01-21 | 河南理工大学 | 基于原子稀疏分解与极限学习机的配电网故障选线方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于原子分解的谐振接地系统自适应故障选线方法;吴靓 等;《高电压技术》;20141130;第40卷(第11期);第3554-3562页 * |
| 小波奇异熵及其在高压输电线路故障选相中的应用;何正友 等;《中国电机工程学报》;20070131;第27卷(第1期);第31-36页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105445618A (zh) | 2016-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Suonan et al. | Distance protection for HVDC transmission lines considering frequency-dependent parameters | |
| CN103344875B (zh) | 谐振接地系统单相接地故障分类选线方法 | |
| CN105445618B (zh) | 一种小电流接地系统故障选线方法及装置 | |
| CN106990324B (zh) | 一种配电网接地故障检测定位方法 | |
| CN108957225B (zh) | 一种计及电缆分布电容的直流配电线路单端故障测距方法 | |
| CN111007364B (zh) | 一种电缆早期自恢复故障的识别方法 | |
| CN106405285A (zh) | 一种电力系统故障录波数据突变时刻检测方法及系统 | |
| CN110954743B (zh) | 一种分布式录波装置及小电流接地选线方法 | |
| CN104316768B (zh) | 一种三相不平衡扰动源定位的负序阻抗参数估算方法 | |
| CN112098886B (zh) | 一种小电流接地系统单相接地选线方法及系统 | |
| CN112098889B (zh) | 一种基于神经网络和特征矩阵的单相接地故障定位方法 | |
| CN109061414A (zh) | 光伏系统直流故障电弧检测方法 | |
| CN103427405B (zh) | 基于高阶累积量的输电线路差动保护方法 | |
| CN110579684A (zh) | 一种基于融合算法的小电流接地系统选线方法 | |
| CN107091958A (zh) | 一种输电线路参数在线辨识系统及其辨识方法 | |
| CN115327300A (zh) | 一种基于录波文件分析的故障起始点定位方法及系统 | |
| CN112345968A (zh) | 一种直流系统接地故障综合检测方法 | |
| CN110646703B (zh) | 一种配电网单相接地故障选线方法及系统 | |
| CN113933647A (zh) | 一种基于首半波功率方向的小电流接地系统故障选线方法 | |
| CN110703134B (zh) | 一种基于故障序分量的小电流接地选线选相方法 | |
| CN117937368A (zh) | 基于频域能量矩阵相似度的柔直电网纵联保护方法 | |
| CN115219787A (zh) | 基于改进矩阵束的电网相量移动测量方法、系统及介质 | |
| Zhu et al. | High impedance grounding fault detection in resonance grounding system based on nonlinear distortion of zero-sequence current | |
| CN112649694B (zh) | 一种小电流接地系统单相接地故障的判定方法 | |
| CN114966324A (zh) | 一种基于改进变分模态分解的单相接地故障定位方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180605 Termination date: 20201120 |