一种基于经验模态分解的电压异常特征识别方法
【技术领域】
[0001] 本发明涉及电力系统电能质量监测分析领域，具体是涉及一种基于经验模态分解 的电压异常特征识别方法。
【背景技术】
[0002] 在现代电力系统中，电力电子设备的应用日益广泛，随之带来各种随机性、非线 性、冲击性、波动性以及暂态性的电力信号，电力系统电能质量日益遭受越来越为严重的污 染。鉴于此，负载端（用户端)对供电可靠性和电能质量分析的要求也日趋增加。另一方面， 电压异常信号的存在，给电力系统带来严重的不良影响。例如:伺服电机的不正常运行，减 少电力设备的使用寿命，增加电能损耗，严重时可能导致供电故障等。
[0003] 针对电力系统电压异常扰动信号问题，国内外对此进行了深入细致的研究工作， 提出了众多提取特征信息的方法，如:傅里叶变换法、小波变换法、加窗傅里叶变换、S变换 等，取得了一些效果的同时也存在一定不足之处。例如：
[0004] 基于傅里叶变换及其一系列改进的分析方法，存在一定的频谱泄漏现象和栅栏效 应，只能分析平稳信号，对暂态和突变扰动信号往往采用加窗算法对其进行修正，然而窗函 数的选取需要依据信号特点，窗函数选取的不适当可能会造成信号失真；
[0005] 基于小波变换的分析方法，在电能质量扰动处理方面得到较为广泛的应用，然而 其极易受信号噪声影响，不适应对时域变化为主的扰动信号；
[0006] S变换是一种新型的时频域分析方法，近年来在电力信号特征信息提取方面取得 长足发展，然而其算法结构复杂，分析长信号时对系统硬件要求高。
【发明内容】
[0007] 为解决上述技术问题，本发明提供了一种基于经验模态分解的电压异常特征识别 方法，为实现从海量电压扰动信号中自动提取并正确分类电能质量检测分析。
[0008] 为实现上述目的，本发明采用了以下技术方案：
[0009] -种基于经验模态分解的电压异常特征识别方法，包括以下步骤：
[0010]步骤A:信号处理器对电力系统电压信号设定采样周期为T，并进行实时采样和量 化得到原始信号v(t)，其中? = 1、2、3···πι，所述m根据实际设备采样率设定，所述采样周期单 位为s，所述电压信号单位为V;
[0011] 步骤B:对原始信号v(t)采用经验模态分解得到各阶頂F模态分量；
[0012] 步骤C:对原始信号v(t)的第一个MF模态分量采用希尔伯特-黄变换得到瞬时频 率图和瞬时幅值图。
[0013] 进一步的，在步骤C后面设有步骤D，其中：
[0014] 步骤D:依据希尔伯特-黄变换图中的突变点、幅值变化、瞬时频率走势特征信息对 电压异常信号进行检测和识别。
[0015] 进一步的技术方案:所述步骤B:对原始信号v(t)采用经验模态分解得到各阶MF 模态分量的具体步骤为：
[0016] S101:第1次筛选:对原始信号v(t)所有局部极大值点和所有局部极小值点用三次 样条函数进行插值，并拟合上、下包络线；
[0017] S102:求取上、下包络线的平均值曲线MKt)，则原始信号奴〇与％(〇之差即为卩丄 (t)J/ft^Pi(t) = v(t)-Mi(t);
[0018] S103:如果Pi(t)同时满足下述MF模态分量的两个条件，则其为第一个頂F模态分 量，否则将其作为新的原始信号重复步骤S101到S102,得到Pn(t)，所述PiKOiPKO-Mn (t)，其中:MWOSPKt)的上、下包络线的平均曲线；
[0019] 所述IMF模态分量满足的两个条件为：（1)整个时间历程内，穿越零点次数与极值 点数相等或至多相差1; (2)且信号上任意一点，由局部极大值定义的上包络线和局部极小 值点定义的下包络线的均值为〇，即信号关于时间轴局部对称；
[0020] S104:重复上述步骤筛选，直到第k次筛选时由式（1)得到的plk(t)满足頂F模态分 量的两个条件：
[0021] Pik(t) =Pi(i-k)(t)-Mik(t) (1);
[0022] S105:在实际计算时可以通过式(2)求取门限值SD来判断每次筛选结果是否为IMF 模态分量：
[0023]
[0024]其中:m为电力系统信号的采样点数，门限值SD取0.2到0.3;
[0025] 3106:令(：1(〇=?11{(〇，则(：1(〇即为第一个頂?模态分量，其包含了原始信号奴〇 中周期最短的MF模态分量;将(^&)从v(t)中分离出来，得到RKt)，所述 ⑴；
[0026] SlOTd^RKt)作为新的原始信号重复以上步骤S101至S105n次，可获得原始信号v (t)的 η 个頂 F 分量 Ci(t)，i = l，2,3…n，且有公式
[0027] S108:iRn(t)为单调函数从原始信号v(t)不能再分解出其他分量时，整个分解过 程结束，此时有如下公式：
[0028]
[0029] 进一步的技术方案:所述步骤C:对原始信号v(t)的第一个IMF模态分量采用希尔 伯特-黄变换得到瞬时频率图和瞬时幅值图的具体步骤为：
[0030] S109:将通过经验模态分解后获得的所有MF模态分量Ci(t)进行希尔伯特-黄变 换，其中i = l，2,3…n，给定C( t)的希尔伯特-黄变换形式为：
[0031]
L其中，λ为积分变量，C( t)为所有IMF模态分 量Ci(t)的统称；
[0032] S110:构造一个解析信号Z(t):
[0033] Z(t) = C(t)+iH(t) = A(t)e10(t)，其中 i为单位虚数，i2 = _l。
[0034] 上式中：幅值函数：
，所述幅值单位为V;
[0035] 幅角函数：
述幅角单位为rad;
[0036] 即可得到，瞬时频率：
〒述频率单位为Hz。
[0037] 进一步的技术方案:所述步骤D:依据希尔伯特-黄变换图中的突变点、幅值变化、 瞬时频率走势特征信息对电压异常信号进行检测和识别的具体步骤为：
[0038] S111:对电力系统中电压异常信号进行分类，所述电压异常信号包括电压骤降、电 压骤升、电压骤停、电压谐波、电压脉冲、电压暂态；
[0039] S112:依据希尔伯特-黄变换图中针对不同的电压异常信号的突变点、幅值变化以 及瞬时频率走势的特性信息不同，分类识别出电力系统中电压发生的异常情况，具体识别 方法为：
[0040] 正常时刻电压瞬时频率维持在50Hz，若电压瞬时频率经历先上升再下降后维持一 段平稳时间，然后再上升再下降；相对应的，电压幅值先上升后经历同样一段平稳时间过程 再下降，由此表明此种异常信号为电压骤降；
[0041 ] 正常时刻电压瞬时频率维持在50Hz，若电压瞬时频率经历先下降再上升后维持一 段平稳时间，然后再下降再上升;相对应的，电压幅值先下降后经历同样一段平稳时间再上 升，由此表明此种异常信号为电压骤升；
[0042]电压骤停异常信号与电压骤降异常信号特征信息相似，区别在于无论是电压瞬时 频率还是电压幅值异常曲线变化是否收敛到零，电压骤停异常瞬时频率和电压幅值收敛到 零，电压骤降不收敛到零；
[0043] 若电压瞬时频率和电压幅值在某一时刻出现尖峰信息，由此表明此种异常信号为 电压脉冲；
[0044] 若电压瞬时频率在某一段时间内经历暂态变化过程后稳定到信号基频附近，并且 电压幅值在同一时间段内经历从零逐渐增加到信号幅值后趋于稳定，由此表明此种异常信 号为电压暂态；
[0045] 若电压瞬时频率在某一时间段内经历剧烈抖动，加之电压幅值在同一时间段内先 经历暂升，后电压幅值剧烈变化后回归稳态，由此表明此种异常信号为电压谐波。
[0046] 进一步的技术方案:所述m根据实际设备采样率设定为1000。
[0047] 本发明的有益效果主要体现在:该方法可以针对电力系统非线性、非平稳信号进 行处理，从希尔伯特-黄变换图谱中可以清晰辨析出信号的不同特征，实现了从海量电压扰 动信号中自动提取并正确分类电能质量检测分析问题。
【附图说明】
[0048] 图1是本发明电压异常信号检测方法的流程图。
[0049] 图2是本发明模拟电压骤降信号IMF1模态分量的希尔伯特-黄变换瞬时频率和幅 值图。
[0050] 图3是本发明模拟电压骤升信号IMF1模态分量的希尔伯特-黄变换瞬时频率和幅 值图。
[0051]图4是本发明模拟电压骤停信号IMF1模态分量的希尔伯特-黄变换瞬时频率和幅 值图。
[0052]图5是本发明模拟电压脉冲信号MF1模态分量的希尔伯特-黄变换瞬时频率和幅 值图。
[0053]图6是本发明模拟电压暂态信号MF1模态分量的希尔伯特-黄变换瞬时频率和幅 值图。
[0054] 图7是本发明模拟电压谐波信号IMF1模态分量的希尔伯特-黄变换瞬时频率和幅 值图。
【具体实施方式】
[0055] 下面结合实例和附图，对本发明的实施方式进行详细的解释和说明。
[0056] 实施案例：
[0057]如图1所示，一种基于经验模态分解的电压异常信号检测方法，包括以下步骤： [0058]步骤A:信号处理器对电力系统信号设定采样周期为T = 0.001s，并进行实时采样 和量化得到原始信号v(t)，t = 1000;
[0059] 步骤B:对原始信号v(t)采用经验模态分解得到各阶頂F模态分量；
[0060] 对原始信号V(t)采用经验模态分解得到各阶頂F模态分量的具体步骤为：
[0061] S101:第1次筛选:对原始信号v(t)所有局部极大值点和所有局部极小值点用三次 样条函数进行插值，并拟合上、下包络线；
[0062] S102:求取上、下包络线的平均值曲线MKt)，则原始信号奴〇与％(〇之差即为卩丄 (t)J/ft^Pi(t) = v(t)-Mi(t);
[0063] S103:如果Pi(t)同时满足下述MF模态分量的两个条件，则其为第一个頂F模态分 量，否则将其作为新的原始信号重复步骤S101到S102,得到Pn(t)，所述PiKOiPKO-Mn (t)，其中:MWOSPKt)的上、下包络线的平均曲线；
[0064] 所述IMF模态分量满足的两个条件为：（1)整个时间历程内，穿越零点次数与极值 点数相等或至多相差1; (2)且信号上任意一点，由局部极大值定义的上包络线和局部极小 值点定义的下包络线的均值为〇，即信号关于时间轴局部对称；
[0065] S104:重复上述步骤筛选，直到第k次筛选时由式（1)得到的Plk(t)满足頂F模态分 量的两个条件，所述式⑴为:Plk(t) =Pl(l-k)⑴-Mlk(t);
[0066] S105:在实际计算时可以通过式(2)求取门限值SD来判断每次筛选结果是否为IMF 模态分量：
[0067] 所述式(2)为
[0068]其中:m为电力系统信号的采样点数，门限值SD取0.2到0.3;
[0069] 3106:令&(〇=?11{(〇，则(：1(〇即为第一个頂?模态分量，其包含了原始