发明名称 基于解析图求解与SVM的小电流接地系统故障选线方法
摘要 本发明涉及一种基于解析图求解与SVM的小电流接地系统故障选线方法。首先,判断故障发生瞬间的故障初相角是否属于[90°,270°),若故障初相角属于[90°,270°),则应用EEMD对1/6工频周期的暂态零序电流进行分解,通过相关系数获得零序特征电流,并对其进行Hilbert变换;然后,选取特征解析点,并分别求取虚相位特征角并构造虚相位特征角向量;最后,利用虚相位特征角向量训练SVM<sub>1</sub>分类器,输入测试集,输出选线结果。若初相角不属于[90°,270°),则利用相同的处理方法求得虚相位特征角向量,进而训练SVM<sub>2</sub>分类器,输入测试集,输出选线结果。理论分析和仿真结果表明,该选线方法原理简单,选线准确率高。
申请公布号 CN103760464B 申请公布日期 2016.05.18
申请号 CN201410028877.1 申请日期 2014.01.07
申请人 河南理工大学 发明人 王晓卫;魏向向;张涛;高杰;侯雅晓
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 代理人
主权项 一种基于解析图求解与SVM的小电流接地系统故障选线方法,其特征在于,故障发生瞬间,判断故障初相角是否属于[90°,270°),若故障初相角属于[90°,270°),应用EEMD(Ensemble Empirical Mode Decomposition,集合经验模态分解)对各分支线路故障时刻起1/6工频周期的暂态零序电流进行分解,通过相关系数获得零序特征电流,并对其进行Hilbert变换;然后,选取特征解析点,并分别求取虚相位特征角并构造虚相位特征角向量;最后,利用虚相位特征角向量训练SVM<sub>1</sub>(Support Vector Machine,支持向量机)分类器,输入测试集,输出选线结果;若初相角不属于[90°,270°),则利用同样的方法求得在该故障初相角范围内的虚相位特征角向量并训练SVM<sub>2</sub>分类器,输入测试集,输出选线结果;该方法具体步骤如下:步骤1当小电流接地系统发生单相接地故障时,故障选线装置立即启动,判断故障初相角δ是否属于[90°,270°),若此时的故障初相角δ属于[90°,270°),则进行以下步骤:步骤①采用100KHz的采样频率记录故障时刻起1/6工频周期内各条分支线路暂态零序电流<img file="FSB0000149833650000011.GIF" wi="217" he="80" />为分支线路的编号,k=1,2,…,q;n为采样点,n=1,2,…,N;步骤②将采集到的1/6工频周期内第k条分支线路暂态零序电流<img file="FSB0000149833650000012.GIF" wi="153" he="81" />进行EEMD分解,得到IMF分量<img file="FSB0000149833650000013.GIF" wi="178" he="80" />h为IMF分量的编号,h=1,2,…U,U为分解得到的IMF分量组数;步骤③遴选零序特征电流<img file="FSB0000149833650000014.GIF" wi="164" he="81" />求取<img file="FSB0000149833650000015.GIF" wi="148" he="81" />与原始暂态零序电流<img file="FSB0000149833650000016.GIF" wi="156" he="81" />的相关系数<img file="FSB0000149833650000017.GIF" wi="118" he="80" />利用相关系数最大的2组IMF分量<img file="FSB0000149833650000018.GIF" wi="223" he="81" />和<img file="FSB0000149833650000019.GIF" wi="228" he="80" />进行重构得到<img file="FSB00001498336500000110.GIF" wi="162" he="81" />求取相关系数<img file="FSB00001498336500000111.GIF" wi="87" he="81" />的计算式如下:<img file="FSB0000149833650000021.GIF" wi="1163" he="323" />其中,<img file="FSB0000149833650000022.GIF" wi="77" he="72" />为第k条暂态零序电流<img file="FSB0000149833650000023.GIF" wi="154" he="71" />的平均值,<img file="FSB0000149833650000024.GIF" wi="85" he="80" />为<img file="FSB0000149833650000025.GIF" wi="154" he="79" />经EEMD分解后得到的第h组IMF分量<img file="FSB0000149833650000026.GIF" wi="151" he="80" />的平均值;零序特征电流<img file="FSB0000149833650000027.GIF" wi="139" he="79" />的计算式如下:<img file="FSB0000149833650000028.GIF" wi="671" he="102" />步骤④对零序特征电流<img file="FSB0000149833650000029.GIF" wi="133" he="71" />进行Hilbert变换,得到变换后的电流<img file="FSB00001498336500000210.GIF" wi="162" he="71" />进而计算零序解析电流<img file="FSB00001498336500000211.GIF" wi="171" he="80" />获取解析点<img file="FSB00001498336500000212.GIF" wi="409" he="72" />零序解析电流<img file="FSB00001498336500000213.GIF" wi="156" he="81" />的计算过程如下:<img file="FSB00001498336500000214.GIF" wi="530" he="150" /><img file="FSB00001498336500000215.GIF" wi="732" he="81" />步骤⑤对Hilbert变换后的电流<img file="FSB00001498336500000216.GIF" wi="144" he="70" />作一阶差分,选出差分绝对值最大且零序特征电流<img file="FSB00001498336500000217.GIF" wi="143" he="81" />最小的解析点,也即特征解析点<img file="FSB00001498336500000218.GIF" wi="444" he="81" />对<img file="FSB00001498336500000219.GIF" wi="144" he="80" />作一阶差分计算式如下:<img file="FSB00001498336500000220.GIF" wi="610" he="81" />步骤⑥计算虚相位特征角<img file="FSB00001498336500000221.GIF" wi="109" he="80" />构建虚相位特征角向量θ<sup>(1</sup><sup>)</sup>:求取各特征解析点<img file="FSB00001498336500000222.GIF" wi="423" he="81" />所对应的虚相位特征角<img file="FSB00001498336500000223.GIF" wi="109" he="79" />也即<img file="FSB00001498336500000224.GIF" wi="421" he="82" />与解析原点G(p,0)构成的虚相位角φ<sub>k</sub>,所求得的<img file="FSB00001498336500000225.GIF" wi="83" he="81" />组成<img file="FSB00001498336500000226.GIF" wi="540" he="88" />虚相位特征角<img file="FSB00001498336500000227.GIF" wi="83" he="80" />的计算式如下:<img file="FSB00001498336500000228.GIF" wi="542" he="179" />步骤⑦训练SVM<sub>1</sub>分类器:在故障初相角δ∈[90°,270°)的范围内,改变δ和 接地电阻R的数值,求得C组虚相位特征角向量θ<sup>(1)</sup>,取L(L<C)组θ<sup>(1)</sup>训练SVM<sub>1</sub>分类器;训练过程中,将输出结果k设定为分支线路k故障,将输出结果0设定为母线故障;步骤⑧利用步骤⑦中剩余的C‑L组虚相位特征角向量θ<sup>(1)</sup>作为测试集,输入SVM<sub>1</sub>分类器中,当输出结果为k时,判定为分支线路k故障,当输出结果为0时,则判定母线故障;步骤2若故障初相角δ不属于[90°,270°),则进行以下步骤:步骤I采用100KHz的采样频率记录故障时刻起1/6工频周期内各条分支线路暂态零序电流<img file="FSB0000149833650000031.GIF" wi="185" he="81" />步骤II将采集到的1/6工频周期内第k条分支线路暂态零序电流<img file="FSB0000149833650000032.GIF" wi="159" he="79" />进行EEMD分解,得到IMF分量<img file="FSB0000149833650000033.GIF" wi="173" he="80" />步骤III遴选零序特征电流<img file="FSB0000149833650000034.GIF" wi="175" he="80" />求取<img file="FSB0000149833650000035.GIF" wi="156" he="81" />与原始暂态零序电流<img file="FSB0000149833650000036.GIF" wi="157" he="80" />的相关系数<img file="FSB0000149833650000037.GIF" wi="122" he="82" />利用相关系数最大与次大的2组IMF分量<img file="FSB0000149833650000038.GIF" wi="222" he="72" />和<img file="FSB0000149833650000039.GIF" wi="227" he="72" />进行重构得到<img file="FSB00001498336500000310.GIF" wi="170" he="80" />求取相关系数<img file="FSB00001498336500000311.GIF" wi="96" he="81" />的计算式如下:<img file="FSB00001498336500000312.GIF" wi="1047" he="323" />其中,<img file="FSB00001498336500000313.GIF" wi="91" he="73" />为第k条暂态零序电流<img file="FSB00001498336500000314.GIF" wi="163" he="80" />的平均值,<img file="FSB00001498336500000315.GIF" wi="91" he="82" />为<img file="FSB00001498336500000316.GIF" wi="162" he="81" />经EEMD分解后得到的第h组IMF分量<img file="FSB00001498336500000317.GIF" wi="155" he="72" />的平均值;零序特征电流<img file="FSB00001498336500000318.GIF" wi="149" he="79" />的计算式如下:<img file="FSB00001498336500000319.GIF" wi="686" he="84" />步骤IV对零序特征电流<img file="FSB00001498336500000320.GIF" wi="148" he="81" />进行Hilbert变换,得到变换后的电流<img file="FSB00001498336500000321.GIF" wi="171" he="81" />进而计算零序解析电流<img file="FSB0000149833650000041.GIF" wi="185" he="81" />获取解析点<img file="FSB0000149833650000042.GIF" wi="425" he="81" />零序解析电流<img file="FSB0000149833650000043.GIF" wi="162" he="80" />的计算过程如下:<img file="FSB0000149833650000044.GIF" wi="541" he="146" /><img file="FSB0000149833650000045.GIF" wi="576" he="84" />步骤V对Hilbert变换后的电流<img file="FSB0000149833650000046.GIF" wi="149" he="80" />作一阶差分,选出差分绝对值最大且零序特征电流<img file="FSB0000149833650000047.GIF" wi="150" he="80" />最小的解析点,也即特征解析点<img file="FSB0000149833650000048.GIF" wi="460" he="81" />对<img file="FSB0000149833650000049.GIF" wi="147" he="80" />作一阶差分计算式如下:<img file="FSB00001498336500000410.GIF" wi="616" he="72" />步骤VI计算虚相位特征角<img file="FSB00001498336500000411.GIF" wi="117" he="79" />构建虚相位特征角向量θ<sup>(2)</sup>:求取各特征解析点<img file="FSB00001498336500000412.GIF" wi="435" he="80" />所对应的虚相位特征角<img file="FSB00001498336500000413.GIF" wi="113" he="80" />也即<img file="FSB00001498336500000414.GIF" wi="435" he="81" />与解析原点G(p,0)构成的虚相位角,所求得的<img file="FSB00001498336500000415.GIF" wi="89" he="80" />组成<img file="FSB00001498336500000416.GIF" wi="565" he="80" />虚相位特征角<img file="FSB00001498336500000417.GIF" wi="90" he="80" />计算式如下:<img file="FSB00001498336500000418.GIF" wi="555" he="178" />步骤VII训练SVM<sub>2</sub>分类器:在故障初相角δ∈[0°,90°)∪[270°,360°]的范围内,改变δ和接地电阻R的数值,求得C组虚相位特征角向量θ<sup>(2)</sup>,取L(L<C)组θ<sup>(2)</sup>训练SVM<sub>2</sub>分类器;训练过程中,将输出结果k设定为分支线路k故障,将输出结果0设定为母线故障;步骤VIII利用步骤VII中剩余的C‑L组虚相位特征角向量θ<sup>(2)</sup>作为测试集,输入SVM<sub>2</sub>分类器中,当输出结果为k时,判定为分支线路k故障,当输出结果为0时,则判定母线故障。
地址 454003 河南省焦作市高新区世纪大道2001号河南理工大学