[发明专利]一种高压直流输电线路雷击干扰识别方法有效
申请号: | 202010852572.8 | 申请日: | 2020-08-21 |
公开(公告)号: | CN111985412B | 公开(公告)日: | 2022-08-05 |
发明(设计)人: | 郝治国;叶冬萌;谢凡 | 申请(专利权)人: | 西安交通大学 |
主分类号: | G06K9/00 | 分类号: | G06K9/00;G06F17/14;G06F17/16;G01R31/58;G01R31/52 |
代理公司: | 西安智大知识产权代理事务所 61215 | 代理人: | 何会侠 |
地址: | 710049 陕*** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 高压 直流 输电 线路 雷击 干扰 识别 方法 | ||
1.一种高压直流输电线路雷击干扰识别方法,其特征在于:基于短数据窗,对不同类型的故障、干扰适应性较好,不受过渡电阻和雷电流参数影响,包括以下步骤:
步骤1:以一定采样频率采集整流侧正、负极电压信号和电流信号,分别对电压信号和电流信号进行相模变换,获取电压1模分量和电流1模分量;
步骤2:设保护启动的时刻为t,雷击干扰识别系统自动记录启动后[t ms,(t+3)ms]内的电压1模分量和电流1模分量;
步骤3:基于CEEMD-HHT计算电压1模分量时频矩阵EU1,具体实现方法如下:
a.通过CEEMD算法将3ms时窗内的电压1模分量u1(n)分解为多个固有模态函数IMF:
其中,ck(n)表示分解得到的第k个IMF,n为时间序列,其中k=1,2,3…K,K为IMF总数,R(n)表示残余函数;
b.定义IMF分量ck(n)的希尔伯特变换HT为:
其中,为第k个IMF ck(n)的希尔伯特变换,m为新引入的积分变量,m∈Z;
得到各个IMF的瞬时幅值ak(n)、瞬时相位φk(n)和瞬时频率fk(n),计算方法如下:
c.根据各IMF的瞬时幅值和瞬时频率,求得电压1模分量u1(n)的时频矩阵EU1,矩阵中的元素计算方法如下:
其中,f为频率变量,f∈[0,fs],fs为信号的采样频率,单位kHz;
步骤4:对各个频率的电压平方值在3ms时窗内进行求和,计算得到频率边际谱EU1(f):
步骤5:以1kHz作为高频、低频的阈值,计算电压1模分量的低频能量ELF_U1和高频能量EHF_U1,并计算其比值ρU1:
步骤6:判断比值ρU1是否超过阈值ρset1,若超过,则判定为故障,进入步骤7;若未超过,则判定为雷击干扰,保护复归;
步骤7:基于CEEMD-HHT计算电流1模分量时频矩阵EI1,具体实现方法如下:
a.通过CEEMD算法将1ms时窗内的电流1模分量i1(n)分解为多个固有模态函数(IMF):
其中,ck(n)表示分解得到的第k个IMF,其中k=1,2,3…K,K为IMF总数,R(n)表示残余函数;
b.根据式(2)得到IMF分量ck(n)的希尔伯特变换,根据式(3)得到各个IMF的瞬时幅值ak(n)、瞬时相位φk(n)和瞬时频率fk(n);
c.根据各IMF的瞬时幅值和瞬时频率,求得电流1模分量i1(n)的时频矩阵EI1,矩阵中的元素计算方法如下:
其中,f为频率变量,f∈[0,fs],fs为信号的采样频率,单位kHz;
步骤8:对各个频率的电流平方值在1ms时窗内进行求和,计算得到频率边际谱EI1(f):
步骤9:以1kHz作为高频、低频的阈值,计算电流1模分量的低频能量ELF_I1和高频能量EHF_I1,并计算其比值ρI1:
步骤10:判断比值ρI1是否超过阈值ρset2,若超过,则判定为普通短路故障,若未超过,则判定为雷击故障,算法结束,保护正常出口。
2.根据权利要求1所述的一种高压直流输电线路雷击干扰识别方法,步骤1所述的采样频率设定为1MHz。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安交通大学,未经西安交通大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010852572.8/1.html,转载请声明来源钻瓜专利网。