CN102854437B - 应用时频原子分解理论的小电流接地系统故障选线方法 - Google Patents

Abstract

本发明提出一种基于时频原子分解理论的小电流接地系统故障选线方法。该方法主要基于时频原子分解理论，将零序电流数据在Gabor过完备原子库中进行稀疏分解，再通过相关参数的优化、求解获得匹配的衰减正弦量原子。时频原子分解法能够准确得到基波和各次谐波的起止时刻、幅值、频率和变化规律等扰动特征，并能够有效的滤除干扰信号。根据时频原子分解后原子的能量熵是按从大到小排列的，除去零序暂态电流基波原子，比较每条线路零序电流频率相近原子的相角（极性）：如果线路零序暂态频率相近的原子相角（极性）与其他线路相反，则为故障线路；如果每条线路对应原子相角（极性）相同，则为母线故障，综合各频率原子相角的比较结果来确定故障线路。

Description

应用时频原子分解理论的小电流接地系统故障选线方法

技术领域

本发明涉及一种配电网的小电流接地系统故障选线方法，特别是涉及一种应用时频原子分解理论进行小电流接地系统故障选线方法。

背景技术

我国大多数配电网均采用中性点不直接接地系统(NUGS)，即小接地电流系统，它包括中性点不接地系统(NUS)，中性点经消弧线圈接地系统(NES，也称谐振接地系统)，中性点经电阻接地系统(NRS)。

我国在小电流接地故障选线方面做了大量的研究，提出了多种选线方法，取得了一定的成效，但是仍然不能完全做到准确可靠的选线，这会阻碍配电网自动化顺利发展，威胁电网的安全稳定运行。

中性点不直接接地系统(NUGS)发生单相接地故障的几率最高，发生单相故障时，系统会产生零序电流，这时供电系统仍能保证线电压的对称性，且故障电流较小，不影响对负荷的连续供电，故不必立即跳闸，规程规定可以继续运行1~2h。但是接地点的出现使得故障相对地电压大幅度降低，非故障相对地电压升高为接近线电压，很容易在电网的绝缘薄弱处引起另一点的接地，从而导致两点或多点接地短路。弧光接地还会引起全系统过电压，进而损坏设备，破坏系统安全运行，所以必须及时找到故障线路予以切除。

国内外学者提出了很多种故障选线的方法。在中性点经消弧线圈接地的系统，发生金属性单相接地时，由于消弧线圈通常处于过补偿状态，故障线路与非故障线路的基波零序电流在数值和方向上都很难区分，已有的稳态量的选线方法^[1-3]很难满足现场运行要求、而已有的暂态量的选线方法^[4-10]中仍存在许多问题，文献[4]在研究S变换提取信号幅频特性和相频特性的基础上，提出了一种基于S变换的融合多个采样点投票结果的配电网故障选线方法，这种方法适用的前提是需要采集到正确的馈线相角和频率信息，文献[5][6]利用小波变换提取故障后的行波信息，构造判据以实现故障选线，小波变换具有良好的时域–频域局部化特性，能提供信号在不同尺度的特征，但易受噪声影响，应用效果不好。文献[7]引入故障测度概念，用Dempster-Shafer证据理论实现了融合的选线方法。文献[8]通过比较暂态零序电流的幅值捕捉特征频带，进而滤波得到特征频带内的信号。文献[9]利用S变换处理各馈线的零序电流，通过比较不同频率点的暂态能量确定容性电流的主导频率，并根据能量的大小选出故障线路。S变换是对连续小波变换和短时傅里叶变换的发展，具有良好的时频特性，但分解后信息量太多。同时还有把稳态量和暂态量两者结合的选线方法^[10]，如神经网络算法，但此算法存在局部最优问题，收敛性较差，训练时间较长，可靠性相对较低。

发明内容

本发明所要解决的技术问题，就是提供一种应用时频原子分解理论的小电流接地系统故障选线方法，其可满足现场运行要求，不易受噪声影响，应用效果较好且分解后信息量不多，收敛性较好，可靠性相对较高。

解决上述技术问题，本发明采用的技术方案如下：

一种应用时频原子分解理论的小电流接地系统故障选线方法，包括以下步骤：

S1建立配电系统发生小电流接地故障时的零序电流数据库：

以母线零序电压瞬时值u(t)大于K_uU_n作为故障启动条件，其中K_u取值为0.15，U_n为母线额定电压，通过选线装置记录故障启动前后2个周波的各馈线零序电流，建立零序电流数据库；

S2对零序电流数据库数据进行时频原子分解，挑选特征量原子：

应用匹配追踪（Matching Pursuits,MP）算法采用离散的Gabor原子库将零序电流数据库数据在Gabor过完备原子库中进行稀疏分解，在满足下式（1）条件下,得到最匹配Gabor原子也即特征量原子（式8），并得到离散的原子参量[s,ξ,τ]；

$\left\{\begin{array}{c}{f_x}^0=f\\ {f_x}^m={f_x}^{m-1}-<{f_x}^{m-1},{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>{g_{\mathrm{\&gamma;}}}^{\left(m\right)}\\ {g_{\mathrm{\&gamma;}}}^{\left(m\right)}=\arg \underset{{g_{\mathrm{\&gamma;}}}^{\left(i\right)}\&Element;D}{\max }|<{f_x}^{m-1},{g_{\mathrm{\&gamma;}}}^{\left(i\right)}>|\end{array}\right.---\left(1\right)$

S3用伪牛顿算法(Pseudo-Newton)将离散的原子参量[s,ξ,τ]连续化，并根据得到的连续化的原子参量[s,ξ,τ]，求出此时的最佳相角φ；

S4根据原子四个参量[s,ξ,τ,φ]，推导对应的衰减正弦量原子，包括以下子步骤：

S4-1检查正弦量原子是衰减还是发散：寻求Gabor原子与当前残余信号具有较大内积的半平面，如果是右半平面(高斯窗函数中心τ的右半部分)，则正弦量原子是衰减的；如果是左半平面(高斯窗函数中心τ的左半部分)，则正弦量原子是发散的；

S4-2由四个参量[s,ξ,τ,φ]中的尺度因子s计算初始衰减因子ρ的估计值：由上一步S4-1的结果，当正弦量原子为衰减时，当正弦量原子为发散时

S4-3确定起始和终止时间t_sq与t_eq：定义t_sq＝m_s和t_eq＝m_e分别为衰减正弦量原子的起始时间和结束时间；当正弦原子衰减时：m_s=τ、m_e＝N-1，此时起始时间已经确定，只需要确定终止时间m_e；在第n次匹配追踪过程中，如果原子与当前信号的内积满足 $<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)>\&GreaterEqual;<f_{x\left(t\right)}^{\left(n\right)},{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>;$ 其中，

$<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)>=\frac{<f_{x\left(t\right)}^{\left(n\right)}{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>-f_{X\left(m_e\right)}^{\left(n\right)}{g_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(m_e\right)}{\sqrt{1-{g_{\mathrm{\&gamma;}}}^{\left(n\right)2}\left(m_e\right)}}---\left(2\right)$

则m_e＝m_e-1，重复上述过程，当 $<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)><<f_{x\left(t\right)}^{\left(n\right)},{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>$ 时，迭代结束；求得最终的终止时间m_e，同理可求得当正弦原子发散时的m_s、m_e；

S4-4根据所得各参数得到衰减正弦量原子

$g\left(t\right)=A_q\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{eq}}))---\left(3\right);$

S4-5，利用伪牛顿算法对衰减正弦量原子的初始衰减因子ρ和频率因子ξ进行优化，并利用优化后衰减正弦量原子再次计算最佳相角，得到最终参变量和最优衰减正弦量原子，存储最终原子参量和原始信号中除去最优衰减正弦量原子的残余信号，并求取残余能量，进行下一次迭代得到第二个原子；当残余信号能量小于原始信号能量的0.01时，迭代结束；

式(1)表明了最佳匹配原子需满足的条件，即提取的原子与当前残余信号具备内积值最大，据此，定义故障选线原子分解能量熵Ei：

E_i=max(|<f_x ^(i-1),g_γ ⁽ⁱ⁾>|),i＝1，2,...,D                 (4)

由式(4)可知，原子分解能量熵是根据原子能量(即故障零序电流的能量)特征自适应求取，内积值最大即表明该原子所含能量最大；

S5设计选线方案

获得每条线路零序暂态电流经原子稀疏分解后按能量熵从大到小排列的原子，除去零序暂态电流基波原子，随后比较每条线路零序电流频率相近原子的相角（极性），如果线路零序暂态频率相近的原子相角（极性）与其他线路相反，则为故障线路，如果每条线路对应原子相角（极性）相同，则为母线故障。

所述的步骤S3包括以下子步骤：

S3-1依次增加参量[s,ξ,τ]的值，增加量为其自身的一半；

S3-2以新参量为标准构成新原子，将新原子与当前残余信号作内积；如果内积值增加，原子参变量以当前值代替初始值，继续相同步骤；反之如果内积值不增加，则将增加值乘以负0.5与自身相加，再进行内积并进行判断；

S3-3在每次迭代过程中，衰减正弦量原子的三参量[s,ξ,τ]按顺序进行优化；

S3-4在判断过程中，如果下面两个条件其中之一得到满足，那么当前参量优化过程终止：(1)内积的增加值不足当前值的1%；(2)当前参变量的增加值不足自身值的10%；

S3-5由求得的原子参量[s,ξ,τ]构造实原子和虚原子，定义G_γ(t)＝g_γ(t)e^jξt+φ，||G_γ(t)||＝1；G_γ(t)为复原子，φ初始值为随机角度，取值范围为[0,2π]；实原子P_γ(t)为复原子G_γ(t)的实部，虚原子Q_γ(t)为复原子G_γ(t)的虚部；用MP算法对零序电流数据进行第m次原子分解迭代后，当前残余值为m次迭代的实原子为P_γ ^(m)(t)，虚原子为Q_γ ^(m)(t)；

1）如果ξ≠0且a≠0，

当 $<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|>0$ 时，φ₀＝arctan(-b/a)；

当 $<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|<0$ 时，φ₀＝arctan(-b/a)+π。

2）如果ξ＝0,

当 $-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|>0$ 时，φ₀＝0；

当 $-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|<0$ 时，φ₀＝φ。

3）如果a＝0，

当 $(f_{x\left(t\right)}^{\left(m\right)}<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>a+<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>b)/\left(a\right|\left|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right||+b|\left|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right|\left|\right)>0$ 时，φ₀＝π/2；

当 $(<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>a+<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>b)/\left(a\right|\left|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|+b\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|\right)<0$ 时，φ₀＝3π/2。

其中 $a=<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>{\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|}^2-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)><{p_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right),{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>,$

$b=<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>{\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|}^2-<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)><{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right),{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>.$

由上述过程可以求得每个原子的最佳相角。

所述的步骤S4-4的根据所得各参数得到衰减正弦量原子包括以下子步骤：

S4-4-1首先构造归一化的衰减正弦量原子：

$g_{\mathrm{\&gamma;}}\left(t\right)=K_{\mathrm{\&gamma;}}\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{eq}}))---\left(5\right)$

其中，K_γ为使||f_γ（t）||=1的系数；频率因子f_q，相位因子φ_q分别与优化后的连续化的原子参量ξ和φ对应，步骤S4-2和S4-3中已经求解出衰减因子ρ，起始和终止时间t_sq与t_eq；

S4-4-2利用MP算法求得归一化后的衰减正弦量原子与当前残余信号的最大内积值 $|<f,{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>|=\underset{\mathrm{\&gamma;}\&Element;\mathrm{\&Gamma;}}{\max }|<f,g_{\mathrm{\&gamma;}}>|,$ 幅值 $A_q=\frac{|<f,g_{\mathrm{\&gamma;}}^{\left(m\right)}>|}{K_{\mathrm{\&gamma;}}};$

S4-4-3由A_q，f_q，ρ_q，φ_q，t_sq，t_eq共6个参数可求得衰减正弦量原子g(t)。（式3所示）。

本发明的理论分析

1、单相接地故障分析

1.1、故障稳态信号特征分析

1.1.1中性点不接地系统（参见图1）

NUS发生单相接地故障，当A相接地时，完好的B、C两相的对地电压升高倍，即升高为线电压；B、C两相对地电容电流也增大为原来电容电流的倍，发生单相接地故障后，虽然相电压不再对称，但是线电压仍然三相对称，三相负荷电流对称。因此，小电流接地故障的稳态电气量具有以下特点：

A.流过故障点的电流数值是正常运行状态下电网三相对地电容电流之和。

B.流过故障元件的零序电流在数值上等于所有非故障元件对地电容电流之和，即故障线路上的零序电流最大。

C.故障线路的零序电流方向为由线路流向母线，非故障线路的零序电流方向为母线流向线路。

1.1.2中性点经消弧线圈接地系统（参见图2）

NES谐振接地电网的单相接地故障，相当于在系统的中性点处，通过消弧线圈注入了一个感性电流来抵消接地点的容性电流，根据补偿度的不同可分为全补偿、欠补偿和过补偿。电力系统中普遍采用过补偿。中性点经消弧线圈接地的电力系统，在单相接地故障时，与中性点不接地的系统一样，其他两相对地电压也要升高到线电压，即升高为原来对地电压的倍。与中性点不接地系统不同的是，经消弧线圈接地的系统允许在系统发生单相接地故障后短时间（一般规定2h）继续运行。在这种情况下，基波中的零序电流最大法和比相法都不能使用，实际中一般采用五次谐波分量法、基于小波分析的故障暂态分析法等，零序电流五次谐波分量在NES 中有着与NUS 中零序电流基波相同的特点，再利用针对NUS的方法，即可解决NES 的选线问题。但负荷中的五次谐波源、CT 不平衡电流和过渡电阻大小都会影响选线精度。

1.2故障暂态信号分析

单相接地时，故障电压和电流的暂态过程持续时间短，并且含有丰富的特征量，因此如果选用一种适合分析暂态分量的方法，将有利于故障选线。

以NES为例，分析故障的暂态过程。

故障瞬间，流过故障点的暂态接地电流存在工频和高频振荡分量，流过故障点的接地电流是由故障相对地电容的放电电流、非故障相对地电容的充电电流和消弧线圈的暂态电感电流叠加而成的。

暂态接地电流表达式为：

I_Cm为电容稳态幅值；ω_f为暂态自由分量振荡分量的角频率；τ_C为电容分量的时间常数；I_Lm为电感电流的稳态幅值；τ_L为电感回路的时间常数。

接地电流的暂态量等于电容电流的暂态自由分量与电感电流的暂态直流分量之和，两者幅值不仅能相互抵消，还会叠加，暂态电流幅值将更大。当故障发生在相电压接近于最大值瞬间，暂态电容电流比暂态电感电流大得多，而且暂态的频率也高，所以在故障初期，电感电流与电容电流是不能相互补偿的，其暂态接地电流的特性主要是由暂态电容电流特性决定的。

中性点不接地系统没有消弧线圈，暂态接地电流就是暂态电容电流，因此可以用暂态电流来进行谐振接地或者是不接地系统的故障选线。

小电流接地故障的暂态电气量具有以下特点：

A.故障线路与非故障线路出现零序电流，非故障线路零序电流等于该线路本身的电容电流，故障线路零序电流为全系统非故障线路对地电容电流总和。

B.非故障线路零序电流超前零序电压90度，故障线路零序电流滞后零序电压90度，因此故障与非故障零序电流相差180度。

因此，可以利用暂态零序电流的幅值和极性来实现选线。

2、时频原子分解法（参见图3）

时频原子分解算法由Mallat和Zhang于1993年提出，该算法将信号在一组过完备的非正交基上分解，分解结果高度稀疏。为了得到信号的稀疏表示，时频原子分解法需要构建一个过完备的展开函数集合，对于某个特定的信号，可根据信号的特征自适应的从过完备集合中选择最佳的展开函数，这样信号分解联系更加紧密，能用更少的函数更准确的表示信号，这个过完备的集合中的展开函数就被称为原子，由原子组成的过完备展开函数集合称为原子库。

2.1、Gabor过完备原子库

为了达到信号稀疏分解的目的，原子库需要是高度冗余的，这样才能保证任意信号都可以从中自适应的选择一组最佳的原子来表示该信号。目前原子稀疏分解中应用最多的是Gabor原子库。表达式为：

$g_{\mathrm{\&gamma;}}\left(t\right)=\frac{1}{\sqrt{s}}g\left(\frac{t-\mathrm{\&tau;}}{s}\right)e^{\mathrm{j\&xi;t}}---\left(7\right)$

其所对应的实Gabor原子为：

$g_{\mathrm{\&gamma;}}\left(t\right)=\frac{K_{\mathrm{\&gamma;}}}{\sqrt{s}}g\left(\frac{t-\mathrm{\&tau;}}{s}\right)\cos (\mathrm{\&xi;t}+\mathrm{\&phi;})---\left(8\right)$

$g\left(t\right)=2^{\frac{1}{4}}{e}^{-\mathrm{\&pi;}{t}^2}---\left(9\right)$

其中(9)式为高斯窗函数，γ＝(s,τ,ξ,φ)，γ为g_γ(t)的索引，s为尺度参量、τ为位移参量、ξ为频率参量、φ为相位参量，参数是为了使||g_γ(t)||＝1。这样的原子空间是无穷的，在实际中不可能搜索一个无穷的空间，所以对原子库进行离散处理。γ＝(s,τ,ξ)离散处理后为γ＝(a^j,pa^jΔτ,ka^-jΔξ)，其中a＝2，Δτ＝1/2，Δξ＝π，

 0＜j＜log₂N                                         (10)

 0＜p＜N2^-j+1                                        (11)

 0＜k＜2^j+1                                          (12)

g_rd(n)＝g_j(n-p2^j)cos(nkπ2^1-j+φ)n＝0,1,…N-1        (13)

$g_j\left(n\right)=\left\{\begin{array}{cc}\mathrm{\&delta;}\left(j\right)& j=0\\ K_{\mathrm{\&gamma;d}}g\left({n2}^{-j}\right)& j\&Element;[1,L)\\ \frac{1}{\sqrt{N}}& j=L\end{array}\right.---\left(14\right)$

其中，L＝log₂N。根据上述Gabor原子库的构建方法，任选参数范围内的参数，可以得到一个Gabor原子，它以横轴u为中心，能量集中在u附近，能量的大小与尺度参数s成比例。由以上离散化过程可知，分频基数为2，随着倍频程j值的递增（即尺度参数s的增加），位移参数的扫描间隔在不断稀疏，而频率参数的扫描间隔在不断密集。

2.2、匹配追踪算法

信号的原子分解过程采用信号的原子分解过程采用匹配追踪(matching pursuits，MP)算法。MP 算法是一种贪婪迭代算法，在每次迭代计算中,依照索引方式扫描原子库从而得到此次迭代中与分析信号最为相关的原子，然后从分析信号中抽取出该最佳原子成分，形成新的残余信号。一般通过控制残余信号的能量形成匹配追踪的结束条件。设D为过完备原子库，f为待分析信号，g_γ为原子库D中的原子。在第m次分解迭代过程中，MP算法寻求原子g_γ ^(m)∈D，使它与当前残余信号f_x ^(m-1)具有最大的内积。

$|<f,{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>|=\underset{\mathrm{\&gamma;}\&Element;\mathrm{\&Gamma;}}{\max }|<f,g_{\mathrm{\&gamma;}}>|---\left(15\right)$

在每次迭代完成后，将当前所得的最佳原子g_γ ^(m)从残余信号中抽取，形成新的残余信号。令最初的残余信号为f_x ⁽⁰⁾＝f，迭代关系式为：

f_x ^(m)＝f_x ^(m-1)-a_mg_γ ^(m)                  (16)

其中，a_m＝<f_x ^(m-1),g_γ ^(m)>，进行m次迭代后，当前残余值为f_x ^(m)，则原信号x可表示为：

$f=\underset{n=1}{\overset{m}{\mathrm{\&Sigma;}}}{a}_n{g_{\mathrm{\&gamma;}}}^{\left(n\right)}+{f_x}^{\left(m\right)}---\left(17\right)$

其中，a_n＝<f_x ^(n-1),g_γ ⁽ⁿ⁾>。因此，经过m次迭代后，信号x可用m个原子的线性组合来表示，其误差为第m次迭代计算后的残差。

2.3、伪牛顿(Pseudo-Newton)法

这种方法的目的是使原始离散的参量实际上更趋于连续变量，弥补了离散的过完备原子库无法实现高密度分解的不足。

具体的优化过程如下：

1）依次增加参量[s,ξ,τ]的值，增加量为其自身的一半；

2）以新参量为标准构成新原子，将新原子与当前残余信号作内积。如果内积值增加，原子参变量以当前值代替初始值，继续相同步骤；反之如果内积值不增加，则将增加值乘以负0.5与自身相加，再进行内积并进行判断；

3）在每次迭代过程中，衰减正弦量原子的三参量[s,ξ,τ]按顺序进行优化；

4）在判断过程中，如果下面两个条件其中之一得到满足，那么当前参量优化过程终止：(1)内积的增加值不足当前值的1%；(2)当前参变量的增加值不足自身值的10%。

2.4、正弦衰减原子

对于线性电力系统模型，信号通常用衰减的正弦量模型来表示，但这时还需要考虑如开关动作等导致的非连续分量，因此电力系统中常采用如下的信号模型来表示电流波形中的扰动信号：

$G\left(t\right)=\underset{q=0}{\overset{Q-1}{\mathrm{\&Sigma;}}}{A}_q\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{sq}}))---\left(15\right)$

式中：u(t)为单位阶跃函数；A_q，f_q，φ_q分别为衰减正弦量的幅值、频率和相位；ρ_q为衰减因子；t_sq和t_eq分别为衰减正弦量的起始时间和结束时间。衰减正弦量原子可由Gabor原子求得，每个衰减正弦量原子包含A_q，f_q，ρ_q，φ_q，t_sq，t_eq共6个参数。

这个模型与Prony分析中的模型相似，

$G\left(t\right)=\underset{q=0}{\overset{Q-1}{\mathrm{\&Sigma;}}}{A}_q\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}---\left(16\right)$

但Prony分析中并未考虑到各分量可能具有不同起始和终止时间。因此这里的信号模型相对于Prony分析增加了时间量定位功能，这也是这种信号模型的优势之处。

本发明与现有技术相比，具有以下优点和效果：

传统的信号展开方法有：傅立叶变换和小波变换等。

傅立叶变换的基函数是频率不同的正弦波集合，它在频率域是全部局部化的，在时间域是随时间周期波动的。傅立叶变换的基函数主要用来分解信号中正弦波，当用它来展开一个非平稳信号时，不能够得到信号的局部成分。另一个重要的基元函数是小波函数，它在时域和频域都具有良好的局部化性质，但是它不能把同一频带中的多个正弦信号分量分解。总之，由于受到展开函数固定有限的制约，基函数表示信号的能力和范围是有限的，如果试图用一个有限的函数或向量集合来表示任意的信号，自然无法较好地适应信号自身的特点；想要得到能够适应信号自身特点的基函数，就需要基函数高度冗余。

现有的实用的小电流接地系统故障选线算法有：

基于傅立叶变换零序电流基波的比幅法、比相法、以及群体比幅比相法。但以上方法不能排除CT不平衡电流及过渡电阻大小的影响，“时针效应”可能存在。

基于小波分析的零序电流暂态过程选线法是利用小波变换把故障发生后的暂态信号分解成不同尺度和位置的小波之和，从而很容易提取出故障选线暂态特征量。例如小波包单频带选线算法，根据小电流系统发生单相接地故障后，非故障线路零序电流中暂态高频分量在特征频带SFB上，图形变换特征基本相同，而故障线路呈现几乎相反的变化特征，构造出选线判据。但实际上由于暂态过程的复杂性，很可能在SFB频段内找不到一个各条线路暂态电容电流都比较集中的频带，从而无法利用SFB频段下信号的极性来选线。

本文采用的基于时频原子分解的选线算法则弥补了这一不足，由于时频原子库是高度冗余的，分解的策略是贪婪自适应的，能从复杂信号和具有特殊结构的信号（例如暂态零序电流）中按能量熵Ei的大小提取特征原子。选取的原子是具有良好的时频局部特性，可准确定量地得到故障成分的起止时刻、幅值、频率和变化规律等扰动特征，获得小电流接地系统中各种成分参量化的原子解析表示，适用于暂态扰动、稳态扰动和多重扰动，而且由该方法得到的能量密度没有Wigner和Cohen类等时频方法的交叉干扰项。按原子分解的特点，可以获得每条线路零序暂态电流经时频原子稀疏分解后按能量熵从大到小排列的原子，除去零序暂态电流基波原子，随后比较每条线路零序电流谐波频率相近原子的相角（极性），综合各频率下的比较结果即可完成选线。仿真结果验证了该方法的可行性和有效性。

附图说明

图1为中性点不接地小电流系统单相接地时电容电流分布图；

图2为中性点接消弧线圈接地小电流系统单相接地时电容电流分布图；

图3为时频原子分解法程序图；

图4为故障线路零序暂态电流原子分解和重构波形的三次迭代分解的原子图；

图5为故障线路零序暂态电流原子分解和重构波形的原始信号、重构信号与残余信号图；

图6为非故障线路零序暂态电流原子分解和重构波形的三次迭代分解的原子图；

图7为非故障线路零序暂态电流原子分解和重构波形的原始信号、重构信号与残余信号图；

图8为时频原子故障选线算法流程图。

具体实施方式

本发明提出了一种应用时频原子分解法的配电网故障选线方法，其实施流程图如图8所示。

本发明的具体实现包括以下步骤：

S1建立配电系统发生小电流接地故障时的零序电流数据库：

以母线零序电压瞬时值u(t)大于K_uU_n作为故障启动条件，其中K_u取值为0.15，U_n为母线额定电压，通过选线装置记录故障启动前后2个周波的各馈线零序电流，建立零序电流数据库；

S2对零序电流数据库数据进行时频原子分解，挑选特征量原子：

应用匹配追踪（Matching Pursuits,MP）算法采用离散的Gabor原子库将零序电流数据库数据在Gabor过完备原子库中进行稀疏分解，在满足下式（1）条件下,得到最匹配Gabor原子也即特征量原子（式8），并得到离散的原子参量[s,ξ,τ]；

$\left\{\begin{array}{c}{f_x}^0=f\\ {f_x}^m={f_x}^{m-1}-<{f_x}^{m-1},{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>{g_{\mathrm{\&gamma;}}}^{\left(m\right)}\\ {g_{\mathrm{\&gamma;}}}^{\left(m\right)}=\arg \underset{{g_{\mathrm{\&gamma;}}}^{\left(i\right)}\&Element;D}{\max }|<{f_x}^{m-1},{g_{\mathrm{\&gamma;}}}^{\left(i\right)}>|\end{array}\right.---\left(1\right)$

S3用伪牛顿算法(Pseudo-Newton)将离散的原子参量[s,ξ,τ]连续化，并根据得到的连续化的原子参量[s,ξ,τ]，求出此时的最佳相角φ；

具体包括以下子步骤：

S3-1依次增加参量[s,ξ,τ]的值，增加量为其自身的一半；

S3-2以新参量为标准构成新原子，将新原子与当前残余信号作内积；如果内积值增加，原子参变量以当前值代替初始值，继续相同步骤；反之如果内积值不增加，则将增加值乘以负0.5与自身相加，再进行内积并进行判断；

S3-3在每次迭代过程中，衰减正弦量原子的三参量[s,ξ,τ]按顺序进行优化；

S3-4在判断过程中，如果下面两个条件其中之一得到满足，那么当前参量优化过程终止：(1)内积的增加值不足当前值的1%；(2)当前参变量的增加值不足自身值的10%；

S3-5由求得的原子参量[s,ξ,τ]构造实原子和虚原子，定义G_γ(t)＝g_γ(t)e^jξt+φ，||G_γ(t)||＝1；G_γ(t)为复原子，φ初始值为随机角度，取值范围为[0,2π]；实原子P_γ(t)为复原子G_γ(t)的实部，虚原子Q_γ(t)为复原子G_γ(t)的虚部；用MP算法对零序电流数据进行第m次原子分解迭代后，当前残余值为m次迭代的实原子为P_γ ^(m)(t)，虚原子为Q_γ ^(m)(t)；

1）如果ξ≠0且a≠0，

当 $<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|>0$ 时，φ₀＝arctan(-b/a)；

当 $<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|<0$ 时，φ₀＝arctan(-b/a)+π。

2）如果ξ＝0,

当 $-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|>0$ 时，φ₀＝0；

当 $-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|<0$ 时，φ₀＝φ。

3）如果a＝0，

当 $(f_{x\left(t\right)}^{\left(m\right)}<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>a+<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>b)/\left(a\right|\left|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right||+b|\left|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right|\left|\right)>0$ 时，φ₀＝π/2；

当 $(<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>a+<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>b)/\left(a\right|\left|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|+b\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|\right)<0$ 时，φ₀＝3π/2。

其中 $a=<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>{\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|}^2-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)><{p_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right),{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>,$

$b=<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>{\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|}^2-<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)><{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right),{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>.$

由上述过程可以求得每个原子的最佳相角。

S4根据原子四个参量[s,ξ,τ,φ]，推导对应的衰减正弦量原子，包括以下子步骤：

S4-1检查正弦量原子是衰减还是发散：寻求Gabor原子与当前残余信号具有较大内积的半平面，如果是右半平面(高斯窗函数中心τ的右半部分)，则正弦量原子是衰减的；如果是左半平面(高斯窗函数中心τ的左半部分)，则正弦量原子是发散的；

S4-2由四个参量[s,ξ,τ,φ]中的尺度因子s计算初始衰减因子ρ的估计值：由上一步S4-1的结果，当正弦量原子为衰减时，当正弦量原子为发散时

S4-3确定起始和终止时间t_sq与t_eq：定义t_sq＝m_s和t_eq＝m_e分别为衰减正弦量原子的起始时间和结束时间；当正弦原子衰减时：m_s=τ、m_e＝N-1，此时起始时间已经确定，只需要确定终止时间m_e；在第n次匹配追踪过程中，如果原子与当前信号的内积满足 $<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)>\&GreaterEqual;<f_{x\left(t\right)}^{\left(n\right)},{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>;$ 其中，

$<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)>=\frac{<f_{x\left(t\right)}^{\left(n\right)}{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>-f_{X\left(m_e\right)}^{\left(n\right)}{g_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(m_e\right)}{\sqrt{1-{g_{\mathrm{\&gamma;}}}^{\left(n\right)2}\left(m_e\right)}}---\left(2\right)$

则m_e＝m_e-1，重复上述过程，当 $<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)><<f_{x\left(t\right)}^{\left(n\right)},{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>$ 时，迭代结束；求得最终的终止时间m_e，同理可求得当正弦原子发散时的m_s、m_e；

S4-4根据所得各参数得到衰减正弦量原子

$g\left(t\right)=A_q\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{eq}}))---\left(3\right);$

具体包括以下子步骤：

S4-4-1首先构造归一化的衰减正弦量原子：

$g_{\mathrm{\&gamma;}}\left(t\right)=K_{\mathrm{\&gamma;}}\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{eq}}))---\left(5\right)$

其中，K_γ为使||f_γ（t）||=1的系数；频率因子f_q，相位因子φ_q分别与优化后的连续化的原子参量ξ和φ对应，步骤S4-2和S4-3中已经求解出衰减因子ρ，起始和终止时间t_sq与t_eq；

S4-4-2利用MP算法求得归一化后的衰减正弦量原子与当前残余信号的最大内积值 $|<f,{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>|=\underset{\mathrm{\&gamma;}\&Element;\mathrm{\&Gamma;}}{\max }|<f,g_{\mathrm{\&gamma;}}>|,$ 幅值 $A_q=\frac{|<f,g_{\mathrm{\&gamma;}}^{\left(m\right)}>|}{K_{\mathrm{\&gamma;}}};$

S4-4-3由A_q，f_q，ρ_q，φ_q，t_sq，t_eq共6个参数可求得衰减正弦量原子g(t)。（式3所示）。

S4-5，利用伪牛顿算法对衰减正弦量原子的初始衰减因子ρ和频率因子ξ进行优化，并利用优化后衰减正弦量原子再次计算最佳相角，得到最终参变量和最优衰减正弦量原子，存储最终原子参量和原始信号中除去最优衰减正弦量原子的残余信号，并求取残余能量，进行下一次迭代得到第二个原子；当残余信号能量小于原始信号能量的0.01时，迭代结束；

式(1)表明了最佳匹配原子需满足的条件，即提取的原子与当前残余信号具备内积值最大，据此，定义故障选线原子分解能量熵Ei：

E_i=max(|<f_x ^(i-1),g_γ ⁽ⁱ⁾>|),i＝1，2,...,D                            (4)

由式(4)可知，原子分解能量熵是根据原子能量(即故障零序电流的能量)特征自适应求取，内积值最大即表明该原子所含能量最大；

S5设计选线方案

获得每条线路零序暂态电流经原子稀疏分解后按能量熵从大到小排列的原子，除去零序暂态电流基波原子，随后比较每条线路零序电流频率相近原子的相角（极性），如果线路零序暂态频率相近的原子相角（极性）与其他线路相反，则为故障线路，如果每条线路对应原子相角（极性）相同，则为母线故障。

参考文献

[1]曾祥君，尹项根，张哲，等.零序导纳法馈线接地保护的研究[J].中国电机工程学报，2001，21(4)：5-10

[2]齐郑，杨以涵.中性点非有效接地系统单相接地选线技术分析[J].电力系统自动化，2004，28（14）：1-5.

[3]束洪春，刘娟，司大军，等.自适应消弧线圈接地系统故障选线实用新方法[J].电力系统自动化，2005，29（13）：64-68.

[4]张仲孝，苗世洪，林湘宁，刘沛.基于多孔算法的小电流接地系统故障选线算法[J].电力系统自动化，2011，35(1)：67-70.

[5]董新洲，毕见广.配电线路暂态行波的分析和接地选线研究[J].中国电机工程学报，2005，25(4)：1-6．

[6]张帆，潘贞存，张慧芬，等.基于方向行波的小电流接地系统故障选线[J].中国电机工程学报，2007，27(34)：70-75．

[7]贾清泉，杨奇逊，杨以涵.基于故障测度概念与证据理论的配电网单相接地故障多判据融合[J].中国电机工程学报,2003,23(12):6-11.

[8]张保会，赵慧梅，张文豪，等.基于特征频带内暂态零序电流特点的配电网单相接地故障选线方法[J].电力系统保护与控制，2008，36(13)：5-10．

[9]束洪春，彭仕欣.基于短窗数据S变换能量的缆-线混合配电网络故障选线方法[J].电工技术学报，2009，24(10)：152-159．

[10]张海平，何正友，张钧，等.基于量子神经网络和证据融合的小电流接地选线方法[J].电工技术学报，2009，24(12)：171-178。

Claims (3)

1.一种应用时频原子分解理论的小电流接地系统故障选线方法，包括以下步骤：

S1建立配电系统发生小电流接地故障时的零序电流数据库：

以母线零序电压瞬时值u(t)大于K_uU_n作为故障启动条件，其中K_u取值为0.15，U_n为母线额定电压，通过选线装置记录故障启动前后2个周波的各馈线零序电流，建立零序电流数据库；

S2对零序电流数据库数据进行时频原子分解，挑选特征量原子：

应用匹配追踪算法采用离散的Gabor原子库将零序电流数据库数据在Gabor过完备原子库中进行稀疏分解，在满足下式(1)条件下,得到最匹配Gabor原子也即特征量原子，并得到离散的原子参量[s,ξ,τ]；

$\left\{\begin{array}{c}{f_x}^0=f\\ {f_x}^m={f_x}^{m-1}-<{f_x}^{m-1},{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>{g_{\mathrm{\&gamma;}}}^{\left(m\right)}\\ {g_{\mathrm{\&gamma;}}}^{\left(m\right)}=\arg \underset{{g_{\mathrm{\&gamma;}}}^{\left(i\right)}\&Element;D}{\max }|<{f_x}^{m-1},{g_{\mathrm{\&gamma;}}}^{\left(i\right)}>|\end{array}\right.---\left(1\right)$

S3用伪牛顿算法将离散的原子参量[s,ξ,τ]连续化，并根据得到的连续化的原子参量[s,ξ,τ]，求出此时的最佳相角φ；

S4根据原子四个参量[s,ξ,τ,φ]，推导对应的衰减正弦量原子，包括以下子步骤：

S4-1检查正弦量原子是衰减还是发散：寻求Gabor原子与当前残余信号具有较大内积的半平面，如果是右半平面，则正弦量原子是衰减的；如果是左半平面，则正弦量原子是发散的；

S4-2由四个参量[s,ξ,τ,φ]中的尺度因子s计算初始衰减因子ρ的估计值：由上一步S4-1的结果，当正弦量原子为衰减时，当正弦量原子为发散时

S4-3确定起始和终止时间t_sq与t_eq：定义t_sq＝m_s和t_eq＝m_e分别为衰减正弦量原子的起始时间和结束时间；当正弦量原子衰减时：m_s＝τ、m_e＝N-1，此时起始时间已经确定，只需要确定终止时间m_e；在第n次匹配追踪过程中，如果原子与当前信号的内积满足 $<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)>\&GreaterEqual;<f_{x\left(t\right)}^{\left(n\right)},{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>;$ 其中，

$<f_{x\left(t\right)}^{\left(n\right)},{P_{{\mathrm{\&gamma;}}^{\&prime;}}}^{\left(n\right)}\left(t\right)>=\frac{<f_{x\left(t\right)}^{\left(n\right)}{P_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(t\right)>-f_{x\left(m_e\right)}^{\left(n\right)}{g_{\mathrm{\&gamma;}}}^{\left(n\right)}\left(m_e\right)}{\sqrt{1-{g_{\mathrm{\&gamma;}}}^{\left(n\right)2}\left(m_e\right)}}---\left(2\right)$

则m_e＝m_e-1，重复上述过程，当时，迭代结束；求得最终的终止时间m_e，同理可求得当正弦量原子发散时的m_s、m_e；

S4-4根据所得各参数得到衰减正弦量原子

$g\left(t\right)=A_q\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{eq}}))---\left(3\right);$

S4-5，利用伪牛顿算法对衰减正弦量原子的初始衰减因子ρ和频率因子ξ进行优化，并利用优化后衰减正弦量原子再次计算最佳相角，得到最终原子参量和最优衰减正弦量原子，存储最终原子参量和原始信号中除去最优衰减正弦量原子的残余信号，并求取残余信号能量，进行下一次迭代得到第二个原子；当残余信号能量小于原始信号能量的0.01时，迭代结束；

式(1)表明了最佳匹配原子需满足的条件，即提取的原子与当前残余信号具备内积值最大，据此，定义故障选线原子分解能量熵Ei：

E_i＝max(|<f_x ^(i-1),g_γ ⁽ⁱ⁾>|),i＝1，2,...,D                 (4)

由式(4)可知，原子分解能量熵是根据原子能量特征自适应求取，内积值最大即表明该原子所含能量最大；

S5设计选线方案

获得每条线路零序暂态电流经原子稀疏分解后按能量熵从大到小排列的原子，除去零序暂态电流基波原子，随后比较每条线路零序电流频率相近原子的相角，如果线路零序暂态频率相近的原子相角与其他线路相反，则为故障线路，如果每条线路对应原子相角相同，则为母线故障。

2.根据权利要求1所述的应用时频原子分解理论的小电流接地系统故障选线方法，其特征是：所述的步骤S3包括以下子步骤：

S3-1依次增加参量[s,ξ,τ]的值，增加量为其自身的一半；

S3-2以新参量为标准构成新原子，将新原子与当前残余信号作内积；如果内积值增加，原子参量以当前值代替初始值，继续相同步骤；反之如果内积值不增加，则将增加值乘以负0.5与自身相加，再进行内积并进行判断；

S3-3在每次迭代过程中，原子参量[s,ξ,τ]按顺序进行优化；

S3-4在判断过程中，如果下面两个条件其中之一得到满足，那么当前原子参量优化过程终止：(1)内积的增加值不足当前值的1％；(2)当前原子参量的增加值不足自身值的10％；

S3-5由求得的原子参量[s,ξ,τ]构造实原子和虚原子，定义G_γ(t)＝g_γ(t)e^jξt+φ，||G_γ(t)||＝1；G_γ(t)为复原子，φ初始值为随机角度，取值范围为[0,2π]；实原子P_γ(t)为复原子G_γ(t)的实部，虚原子Q_γ(t)为复原子G_γ(t)的虚部；用MP算法对零序电流数据进行第m次原子分解迭代后，当前残余值为m次迭代的实原子为P_γ ^(m)(t)，虚原子为Q_γ ^(m)(t)；

1)如果ξ≠0且a≠0，

当 $<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|>0$ 时，φ₀＝arctan(-b/a)；

当 $<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|<0$ 时，φ₀＝arctan(-b/a)+π；

2)如果ξ＝0,

当 $-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|>0$ 时，φ₀＝0；

当 $-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>/\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|<0$ 时，φ₀＝φ；

3)如果a＝0，

当 $(f_{x\left(t\right)}^{\left(m\right)}<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>a+<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>b)/\left(a\right|\left|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right||+b|\left|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right|\left|\right)>0$ 时，φ₀＝π/2；

当 $(<f_{x\left(t\right)}^{\left(m\right)}{,P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>a+<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>b)/\left(a\right|\left|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\right||+b\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|)<0$ 时，φ₀＝3π/2；

其中 $a=<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>{\left|\right|{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|}^2-<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)><{p_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right),{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>,$

$b=<f_{x\left(t\right)}^{\left(m\right)},{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>{\left|\right|{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)\left|\right|}^2-<f_{x\left(t\right)}^{\left(m\right)},{P_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)><{P_r}^{\left(m\right)}\left(t\right),{Q_{\mathrm{\&gamma;}}}^{\left(m\right)}\left(t\right)>;$

由上述过程可以求得每个原子的最佳相角。

3.根据权利要求1所述的应用时频原子分解理论的小电流接地系统故障选线方法，其特征是：所述的步骤S4-4的根据所得各参数得到衰减正弦量原子包括以下子步骤：

S4-4-1首先构造归一化的衰减正弦量原子：

$g_{\mathrm{\&gamma;}}\left(t\right)=K_{\mathrm{\&gamma;}}\cos (2\mathrm{\&pi;}{f}_{q}t+{\mathrm{\&phi;}}_q)e^{-{\mathrm{\&rho;}}_q(t-t_{\mathrm{sq}})}\&times;(u(t-t_{\mathrm{sq}})-u(t-t_{\mathrm{eq}}))---\left(5\right)$

其中，K_γ为使||f_γ(t)||＝1的系数；频率因子f_q，相位因子φ_q分别与优化后的连续化的原子参量ξ和φ对应，步骤S4-2和S4-3中已经求解出衰减因子ρ，起始和终止时间t_sq与t_eq；

S4-4-2利用MP算法求得归一化后的衰减正弦量原子与当前残余信号的最大内积值 $|<f,{g_{\mathrm{\&gamma;}}}^{\left(m\right)}>|\underset{\mathrm{\&gamma;}\&Element;\mathrm{\&Gamma;}}{\max }|<f,g_{\mathrm{\&gamma;}}>|,$ 幅值 $A_q=\frac{|<f,g_{\mathrm{\&gamma;}}^{\left(m\right)}>|}{K_{\mathrm{\&gamma;}}};$

S4-4-3由A_q，f_q，ρ_q，φ_q，t_sq，t_eq共6个参数按式(3)可求得衰减正弦量原子g(t)。