中性点不接地系统铁磁谐振故障检测方法
【专利摘要】本发明公开了一种中性点不接地系统铁磁谐振故障检测方法，本发明对母线三相电压和零序电压信号进行分析来检测系统的铁磁谐振故障，本发明以基于经验模态分解(EMD)的Hilbert变换(HHT)为核心，结合FFT变换，先对三相电压和零序电压信号进行FFT分析，再进行HHT分析，最后由三维Hilbert谱得出铁磁谐振发生的时刻及铁磁谐振过电压幅值大小，实现铁磁谐振类型、幅值、诱发时刻的准确定量分析。本发明采用HHT信号分析法彻底摆脱了线性和平稳性束缚，即可实现电压信号的频域分析，同时又可进行时域分析；进行HHT分析时采用RBF-点对称延拓法结合镜像延拓法改善HHT的端点效应，避免模态重叠。
【专利说明】中性点不接地系统铁磁谐振故障检测方法
【技术领域】
[0001] 本发明属于电力系统安全性分析领域，涉及一种中性点不接地系统铁磁谐振故障 检测方法。
【背景技术】
[0002] 中性点不接地系统中包含有许多电感、电容性元件（如变压器、互感器等的电感； 输电线路的对地电容及相间电容），可W构成一系列不同频率的振荡回路。在该些回路中， 电压互感器（简称PT)是铁芯电感元件，如果有某种大的扰动或操作，PT的非线性铁芯就 可能饱和，从而与线路和设备的对地电容形成特殊的单相或H相共振回路，激发起持续的 较高幅值的铁磁谐振过电压。谐振时出现的异常过电压和过电流引起绝缘闪络、避雷管爆 炸、设备损坏，严重时造成停电事故，严重威胁电网安全运行。铁磁谐振的理论分析和计算 主要有：
[0003] (1)图解法、相平面法，多在对铁磁谐振发生机理进行定性的分析；
[0004] (2)用非线性系统的分析法对谐振电路进行分析，如幅频法、平均法、谐波平衡法 等；
[0005] (3)国外学者把铁磁谐振与非线性动态系统和混沛分析结合起来，将分叉理论、奇 异和非奇异吸引子的概念引入铁磁谐振的研究领域，利用功率谱密度和鹿加莱映射的方法 和数字仿真技术对其进行动态分析；
[0006] (4)用数字仿真方法对铁磁谐振进行稳态和暂态计算。
[0007] 到目前为止，国内对于铁磁谐振的数值计算研究可分为两大类：
[0008] (1)在建立的数学模型基础上，用一些拟定的参数进行计算得出有关PT谐振的规 律；
[0009] (2)采用国外的电力系统电磁暂态计算程序（EMTP)，对实际系统进行仿真计算。 但该程序实际上没有专口针对铁磁谐振现象的计算，所W仿真效果并不理想，并且该程序 的使用和维护都很复杂。因此，有必要对配电网中的铁磁谐振故障进行研究，并提出中性点 不接地系统铁磁谐振故障检测方法，帮助运维检修人员进行分析。
【发明内容】
[0010] 为了解决现有技术中存在的问题，本发明提供了一种中性点不接地系统铁磁谐振 故障检测方法。
[0011] 本发明的上述技术问题主要是通过下述技术方案得W解决的：
[0012] 一种中性点不接地系统铁磁谐振故障检测方法，包括W下步骤：
[0013] 步骤一，故障录波仪记录故障时母线H相电压和零序电压信号U。；
[0014] 步骤二，对步骤一中记录下的母线H相电压和系统中性点U。电压进行快速傅氏变 换（FFT)，得到信号所含频率成分，构造出指定带宽的滤波器，对信号进行滤波，使各频率成 分分开，初步判断单相接地故障和铁磁谐振故障；
[0015] 步骤H，对经过步骤二处理的H相电压信号和零序电压信号通过经验模态分解 (EMD)分解为若干个固有模态函数IMF,再对每个IMF进行Hilbed变换，得到每个IMF随 时间变化的瞬时频率和瞬时幅值，构建信号的时间-频率-能量分布即H维化化ed谱，实 现铁磁谐振类型、幅值、诱发时刻的准确定量分析。
[0016] 更优的是，所述的步骤一中，至少要记录故障前tl砂到故障后t2砂期间的母线H 相电压和零序电压波形，其中tl G [0. 1,0. 2]，t2 G巧.0,3. 0]。
[0017] 更优的是，所述的步骤二中，初步判断单相接地故障和铁磁谐振故障的依据是中 性点不接地系统单相接地故障在故障发生的0. Is内只含有50化的频率分量，铁磁谐振的 0. Is内波形发生崎变，包含有很多相同频率的谐波，并能在W后的时间内逐渐进入稳定的 长久谐振或者谐振后自行消失。
[0018] 更优的是，所述的步骤H中，运用相似系数法从一组IMF中提取出能表征原信号 主要特征的主模态分量。
[0019] 更优的是，所述的步骤H中，对得到的每个IMF随时间变化的瞬时频率和瞬时幅 值进行最小二乘法拟合与均值运算，得到各分量的频率和幅值。
[0020] 更优的是，所述的步骤H中，H相电压信号和零序电压信号进行HHT分析的步骤 如下：
[0021] 1)记原始电压信号为X (t)，对于任意电压信号X (t)，经验模态分解EMD分解过程 如下：
[002引 （1)找到信号X (t)的所有极大值和极小值点，分别用曲线拟合，得到X (t)的上包 络线Ui (t)和下包络线U 2似；
[002引 似上、下包络线的平均值记为m(t)，则m的二非'1的+化的]，并令h(t)= X (t)-m (t)，则 h (t)为一个近似的 IMF ;
[0024] 做将h(t)作为新的x(t)，重复（1)、似操作，直到h(t)满足IMF条件停止，此 时得到第一阶IMF记作Cl (t);
[002引 （4)将r(t) =x(t)-ci(t)作为新的x(t)，重复W上过程，依次得到第二阶、第立阶 巧 IMF分量，……，最终可得电压信号二；其中，r(t)称为残余函数， k=l 代表信号的平均趋势；
[002引 2)对所有的IMF进行Hi化ert变换，记Xi (t)二Ci (t)，Xi (t)的Hi化ert变换 OO OO 为：?/佩二主/ 反变换则为正倘=i /葦學如，因化对原电压信号经 -OO -OO 过EMD分解得到的每个IMF分量做化化ed变换后，就能得到每个IMF分量的解析式： -啦）二馬的+J?/''的二a啦)巳抑制，Si (t)为每个IMF分量随时间变化的瞬时幅值， 白心）为相角，贝IJ每个IMF分量随时间变化的瞬时频率为乂'的=去^^。
[0027]扣根据每个IMF分量的瞬时频率和瞬时幅值，构建电压信号的时间-频率-能量 分布即H维Hilbed谱，实现铁磁谐振类型、幅值、诱发时刻的准确定量分析。
[0028] 更优的是，所述的步骤H中，HHT分析时通过RBF-点对称延拓法结合镜像延拓法 改善HHT的端点效应。
[0029] 综上所述，本发明与现有技术相比，具有W下优点：
[0030] 1、传统的铁磁谐振判别方法，如傅立叶变换只能处理线性、平稳的电压信号，小波 变换在实际算法实现中只能处理线性非平稳电压信号。本发明采用HHT信号分析法彻底摆 脱了线性和平稳性束缚，即可实现电压信号的频域分析，同时又可进行时域分析。
[0031] 2、先采用FFT分析法初步判别单相接地和铁磁谐振，再用基于经验模态分解 (EMD)的HHT分析法检测出基频谐振、分频谐振及高频谐振，并由H维化化ed谱得出铁磁 谐振发生的时刻及铁磁谐振过电压幅值大小。
[0032] 3、对电压信号进行HHT分析时采用RBF-点对称延拓法结合镜像延拓法改善HHT 的端点效应，避免模态重叠。
【专利附图】
【附图说明】
[0033] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W 根据该些附图获得其他的附图。
[0034] 图1为本发明的工作流程图；
[00巧]图2为本发明的算法流程图。
【具体实施方式】
[0036] 下面结合附图对本发明的优选实施例进行详细阐述，W使本发明的优点和特征能 更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。
[0037] 本发明的中性点不接地系统铁磁谐振故障检测方法的具有W下优点：
[0038] 1、传统的铁磁谐振判别方法，如傅立叶变换只能处理线性、平稳的电压信号，小波 变换在实际算法实现中只能处理线性非平稳电压信号。本发明采用HHT信号分析法彻底摆 脱了线性和平稳性束缚，即可实现电压信号的频域分析，同时又可进行时域分析。
[0039] 2、先采用FFT分析法初步判别单相接地和铁磁谐振，再用基于经验模态分解 (EMD)的HHT分析法检测出基频谐振、分频谐振及高频谐振，并由H维化化ed谱得出铁磁 谐振发生的时刻及铁磁谐振过电压幅值大小。
[0040] 3、对电压信号进行HHT分析时采用RBF-点对称延拓法结合镜像延拓法改善HHT 的端点效应，避免模态重叠。
[0041] 实施例：一种中性点不接地系统铁磁谐振故障检测方法，包括W下步骤：
[0042] 步骤一，故障录波仪记录故障时母线H相电压和零序电压信号U。；至少要记录故 障前tl砂到故障后t2砂期间的母线H相电压和零序电压波形，其中tl G [0. 1,0.2]， t2 G [2. 0,3. 0]；
[0043] 步骤二，对步骤一中记录下的母线H相电压和系统中性点U。电压进行快速傅氏变 换（FFT)，得到信号所含频率成分，构造出指定带宽的滤波器，对信号进行滤波，使各频率成 分分开，初步判断单相接地故障和铁磁谐振故障；初步判断单相接地故障和铁磁谐振故障 的依据是中性点不接地系统单相接地故障在故障发生的0. Is内只含有50化的频率分量， 铁磁谐振的0. Is内波形发生崎变，包含有很多相同频率的谐波，并能在W后的时间内逐渐 进入稳定的长久谐振或者谐振后自行消失。
[0044] 步骤H，对经过步骤二处理的H相电压信号和零序电压信号通过经验模态分解 (EMD)分解为若干个固有模态函数IMF,再对每个IMF进行Hilbed变换，得到每个IMF随 时间变化的瞬时频率和瞬时幅值，构建信号的时间-频率-能量分布即H维化化ed谱，实 现铁磁谐振类型、幅值、诱发时刻的准确定量分析；
[0045] 运用相似系数法从一组IMF中提取出能表征原信号主要特征的主模态分量；
[0046] 对得到的每个IMF随时间变化的瞬时频率和瞬时幅值进行最小二乘法拟合与均 值运算，得到各分量的频率和幅值；
[0047] H相电压信号和零序电压信号进行皿T分析的步骤如下：
[004引1)记原始电压信号为X (t)，对于任意电压信号X (t)，经验模态分解EMD分解过程 如下：
[004引 （1)找到信号X (t)的所有极大值和极小值点，分别用曲线拟合，得到X (t)的上包 络线Ui (t)和下包络线U 2似；
[0050] 似上、下包络线的平均值记为m(t)，则m的=|[巧的+巧削，并令h(t)= X (t)-m (t)，则 h (t)为一个近似的 IMF ;
[0051] 做将h(t)作为新的x(t)，重复（1)、似操作，直到h(t)满足IMF条件停止，此 时得到第一阶IMF记作Cl (t);
[005引 （4)将r(t) =x(t)-ci(t)作为新的x(t)，重复W上过程，依次得到第二阶、第；阶 巧 IMF分量，……，最终可得电压信号X的二玄邱的+ r闲，其中，r(t)称为残余函数， A-=I 代表信号的平均趋势；
[005引。对所有的IMF进行Hi化ert变换，记Xi (t) = Ci (t)，Xi (t)的Hi化ert变换 为：?/佩二主/ ￥^^如，反变换则为取'的=主7學耸dT，因此，对原电压信号经 -OO -OO 过EMD分解得到的每个IMF分量做HHbed变换后，就能得到每个IMF分量的解析式： 2；;的=正啦）+ = a'.(f)e州的，ai(t)为每个IMF分量随时间变化的瞬时幅值， 白心）为相角，贝IJ每个IMF分量随时间变化的瞬时频率为知亡）二去马驾!^。
[0054]扣根据每个IMF分量的瞬时频率和瞬时幅值，构建电压信号的时间-频率-能量 分布即H维Hilbed谱，实现铁磁谐振类型、幅值、诱发时刻的准确定量分析。
[00巧]HHT分析时通过RBF-点对称延拓法结合镜像延拓法改善HHT的端点效应。
[0056] W上仅仅W-个实施方式来说明本发明的设计思路，在系统允许的情况下，本发 明可W扩展为同时外接更多的功能模块，从而最大限度扩展其功能。
[0057] W上所述，仅为本发明的【具体实施方式】，但本发明的保护范围并不局限于此，任何 不经过创造性劳动想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的 保护范围应该W权利要求书所限定的保护范围为准。
【权利要求】
1. 一种中性点不接地系统铁磁谐振故障检测方法，其特征在于：包括以下步骤： 步骤一，故障录波仪记录故障时母线三相电压和零序电压信号Utl; 步骤二，对步骤一中记录下的母线三相电压和系统中性点Utl电压进行快速傅氏变换 (FFT)，得到信号所含频率成分，构造出指定带宽的滤波器，对信号进行滤波，使各频率成分 分开，初步判断单相接地故障和铁磁谐振故障； 步骤三，对经过步骤二处理的三相电压信号和零序电压信号通过经验模态分解（EMD) 分解为若干个固有模态函数MF，再对每个IMF进行Hilbert变换，得到每个IMF随时间变 化的瞬时频率和瞬时幅值，构建信号的时间-频率-能量分布即三维Hilbert谱，实现铁磁 谐振类型、幅值、诱发时刻的准确定量分析。
2. 根据权利要求1所述的中性点不接地系统铁磁谐振故障检测方法，其特征在于，所 述的步骤一中，至少要记录故障前tl秒到故障后t2秒期间的母线三相电压和零序电压波 形，其中tle[〇· 1，〇·2]，t2e[2·0,3·0]。
3. 根据权利要求1所述的中性点不接地系统铁磁谐振故障检测方法，其特征在于，所 述的步骤二中，初步判断单相接地故障和铁磁谐振故障的依据是中性点不接地系统单相接 地故障在故障发生的0.Is内只含有50Hz的频率分量，铁磁谐振的0.Is内波形发生畸变， 包含有很多相同频率的谐波，并能在以后的时间内逐渐进入稳定的长久谐振或者谐振后自 行消失。
4. 根据权利要求1所述的中性点不接地系统铁磁谐振故障检测方法，其特征在于，所 述的步骤三中，运用相似系数法从一组頂F中提取出能表征原信号主要特征的主模态分 量。
5. 根据权利要求1所述的中性点不接地系统铁磁谐振故障检测方法，其特征在于，所 述的步骤三中，对得到的每个IMF随时间变化的瞬时频率和瞬时幅值进行最小二乘法拟合 与均值运算，得到各分量的频率和幅值。
6. 根据权利要求1所述的中性点不接地系统铁磁谐振故障检测方法，其特征在于，所 述的步骤三中，三相电压信号和零序电压信号进行HHT分析的步骤如下： 1) 记原始电压信号为X(t)，对于任意电压信号X(t)，经验模态分解EMD分解过程如 下： (1) 找到信号x(t)的所有极大值和极小值点，分别用曲线拟合，得到x(t)的上包络线 Mt)和下包络线Mt); (2) 上、下包络线的平均值记为m(t)，则m(i) = +t'2(i)]，并令h(t)= x(t)-m(t)，则h(t)为一个近似的IMF; ⑶将h(t)作为新的x(t)，重复（1)、⑵操作，直到h(t)满足IMF条件停止，此时得 到第一阶IMF记作C1⑴； (4)将r(t)=X (t)-C1 (t)作为新的X(t)，重复以上过程，依次得到第二阶、第三阶MF η 分量，……，最终可得电压信号=ΣA⑷+Mi)，其中，r(t)称为残余函数，代表 k=l 信号的平均趋势； 2) 对所有的IMF进行Hilbert变换，记Xi (t) =Ci (t)，Xi⑴的Hilbert变换为： 识(?) = 士 7gdT·，反变换则为而·⑷7 @dr，因此，对原电压信号经 -OO-OO 过EMD分解得到的每个MF分量做Hilbert变换后，就能得到每个MF分量的解析式：々'(/:)=而⑷+加⑴=ajt)为每个IMF分量随时间变化的瞬时幅值， Si⑴为相角，贝IJ每个MF分量随时间变化的瞬时频率为上⑷=士^^〇 3)根据每个IMF分量的瞬时频率和瞬时幅值，构建电压信号的时间-频率-能量分布 即三维Hilbert谱，实现铁磁谐振类型、幅值、诱发时刻的准确定量分析。
7.根据权利要求1所述的中性点不接地系统铁磁谐振故障检测方法，其特征在于，所 述的步骤三中，HHT分析时通过RBF-点对称延拓法结合镜像延拓法改善HHT的端点效应。
【文档编号】G01R31/00GK104237683SQ201410441696
【公开日】2014年12月24日 申请日期:2014年9月1日 优先权日:2014年9月1日
【发明者】何菲, 孙志, 陈国珍, 王倩, 蒋浩 申请人:国家电网公司, 江苏省电力公司, 江苏省电力公司泰州供电公司, 东南大学