1 单一电能质量事件的监测

电压暂降、暂升、中断的监测

信号的采样频率为 1 kHz(本文信号中的噪声、采样频率均如此设置)。分别采用全周 FT 算法、基于 Morlet复小波的 CWT(中心频率设置为 50 Hz)和 TFT 方法(TFA 的中心频率为 50 Hz，频域窗口半径为 13 Hz)测量波形中基波分量的瞬时幅值并比较其动态特性 。

基于 TFT 的电压暂降监测结果可见：对于单一电压暂降事件，上述 3 种方法都可准确检测出暂降事件开始和结束的时刻；暂降的幅值和暂降过程中的频率与实际情况相吻合；TFT 的动态响应速度也要优于连续小波变换方法。类似地， TFT 可对单一电压暂升和电压中断 2 种电能质量事件进行有效监测 。

TFT 的动态特性为 2 个周波，因此 TFT 无法精确检测 0.5~2 个周波的瞬时电压暂降(暂升)、电压中断的幅值。实验证明，当电压暂降(暂升)的持续时间小于 2 个周波时，基于 TFT 的测量结果会比实际值略大(小)，误差约为 10%，因此本文采用 TFT 和半波 FT 算法相结合的方法，当TFT 检测到电压暂降(升)、电压中断的持续时间小于 2 个周波时，依据 TFT 测量的频率来设置半波FT 算法的参数实现其同步采样，该软件自动启动半波 FT 算法检测电压的幅值，并将其作为 TFT 的修正值，用于上述瞬时电压暂降(升)、电压中断的特征提取 。

频率偏移的监测

用仿真软件中生成频率为 50.5 Hz，幅值为 1 pu的频率偏移扰动波形，即 s(t) = sin(2 × 50.5πt) n(t)，t = 0~0.6 s。利用 TFT 得到的频率可看出，测量结果与实际情况完全相符， 说明 TFT 可以在非同步采样条件下准确监测电能扰动波形中的基波频率偏移，基波幅值测量精度不受频率偏移的影响 。

电压波动的监测

用仿真软件生成电压波动的扰动波形，基波频率为 50 Hz，幅值为 1 pu 且波动频率为 4 Hz，电压波动值为 0.32 pu，即 s = sin(2 × 50πt)(1 0.16sin(2 ×4πt)) n(t) 。

2 多重电能质量事件的监测

电能扰动波形中往往同时存在多种电能质量事件，需要对每一个电能质量事件进行精确监测和准确识别。0~0.6 s 期间电能扰动波形中存在频率偏移、 3次谐波和间谐波 3 种电能质量事件； 0.2~0.4 s 期间另有电压暂升。

采用 TFT 逐采样间隔监测依次测量基波和谐波、间谐波分量的参数，如图 5、 6 所示。 TFT 不受基波频率波动的影响和各信号分量相互之间的干扰，能够准确检测出基波分量的 2 种电能质量事件(电压暂升和频率偏移)并精确测量出基波分量的幅值和频率，且动态响应速度也优于复连续小波变换方法。而全周 FT 算法受频率偏移和谐波、间谐波分量的影响，幅值和频率测量误差较大。同时TFT 能够在多种电能质量事件同时存在的情况下准确地检测出扰动波形中的谐波和间谐波，且幅值和频率的监测精度远优于复连续小波变换方法 。

利用 TFT 提取基波幅值标么值 U1、基波频率标么值 f、谐波畸变率 THD1、间谐波畸变率THD2、基波幅值波动值的百分数 d(以下简称为波动值)，将其组成一个 5 维向量，用于上述 7 种电能扰动现象的类型识别。考虑到对每个检测数据进行识别计算量太大，故对各扰动分量的幅值和频率值从初始时刻对应的数据开始每隔 0.1 s 求一次平均。其 中 波 动 值 d = (Umax − Umin) /UN × 100% ， UN 与Umax、 Umin 分别为电压额定值和 0.1 s 时间窗口内基波幅值的最大值、最小值， pu。从电能扰动波形中提取的特征参数构成的 5 维模式特征向量可以为单一或同时存在的多重电能质量事件的模式识别提供准确、可靠的模式信息 。

近年来，电能质量扰动检测是电力工程界的热点和难点问题。由于扰动信号大多属于非线性信号，一些新的非线性信号分析方法被引入电能质量扰动检测领域。如小波变换、S 变换、改进 S变换、数学形态学、分形理论、时频原子算法 、 原 子 分 解 算 法和 希 尔 伯 特 黄 变 换(Hilbert-Huang transform，HHT)等方法。这些非线性分析方法虽在电能质量扰动检测领域取得了较好的检测结果，但也各自存在一些问题。如小波变换的检测效果不仅受 Heisenberg 测不准原理制约，而且检测的效果取决于基函数的选择和分解尺度，无法保证最优的分解效果；S 变换根据 S 矩阵的幅值矩阵很难考察扰动信号频率随时间的分布，检测暂态扰动信号时效果也不太理想；改进 S 变换中如何选取高斯窗调节因子缺乏理论依据；用多刻度形态学中的形态谱可表征电能质量中的各种扰动，但各刻度下形态谱的大小受结构函数幅值和形状、采样频率的影响很大；分形理论应用于电能质量扰动检测处于起步阶段， 文献利用小波变换和分形指数提取动态电能质量的扰动特征量，但小波与分形之间的定量关系在进一步地研究中；时频原子算法无法精确地检测出小于 2 个周波的瞬时电压暂降(暂升)的特征参数；原子分解算法对如何构造波动和闪变的电能质量扰动相关原子库还需要进一步研究；HHT 首先用经验模态分解(empirical mode decomposition， EMD)将复杂信号分解为若干固有模态函数(intrinsic mode function，IMF) 分 量 之 和 ， 然 后 用 希 尔 伯 特 变 换 (Hilberttransform，HT)求取每个 IMF 的频率和幅值，根据频率突变点定位扰动时刻，该方法有很好的自适应，克服了小波变换、S 变换等传统时频分析方法的局限性。但 HHT 存在一些理论问题，如端点效应、IMF 判据、没有快速算法和过包络等问题。另外， 因受 Bedrosian 和 Nuttall 理论的限制， 通过 HT获取的幅值在端部明显失真。这些理论问题的存在导致应用 HHT 分析电能质量扰动时不仅在端点处的分析效果较差，而且瞬时幅值函数波动较严重 。

2005 年，Smith 等人提出了局部均值分解算法(local mean decomposition，LMD)，LMD 可将复杂信号分解为乘积函数(product function，PF)之和。每个 PF 由包络函数和纯调频函数之积组成，包络函数是 PF 的瞬时幅值，纯调频函数的频率即为 PF的瞬时频率。LMD 和 HHT 类似，也是根据信号固有特征尺度分解复杂信号， 但 LMD 获取 PF 分量的迭代过程采用除法运算，而 EMD 获取 IMF 的迭代过程采用减法，较之 EMD 获取一个 IMF 分量的迭代次数，LMD 获取一个 PF 分量的迭代次数明显较少，而迭代次数越少，端点效应污染数据序列的程度就越轻，幅值与频率检测结果较为准确 。

HHT 中采用 HT 获取 IMF 分量的幅值， 由于 HT 的边缘效应，在端点处的幅值和频率信息会出现部分失真。而 LMD 中将包络估计函数相乘得到幅值信息，端部失真较小。最初 LMD 用于脑电信号分析，程军圣等人将其应用于机械故障诊断，杨世锡将其应用于信号瞬时频率的提取，最近唐巍[16]将 LMD 应用于电力系统低频振荡分析，但将LMD 用于分析含有高频暂态、脉冲等电能质量扰动信号的研究工作还未见报道 。

为准确监测电能扰动波形中基波分量的特征参数同时保证良好的动态特性， 设置 TFA 的中心频率为 50 Hz，调整尺度参数使频率窗口半径为 13 Hz。TFA 在时域上逐采样间隔滑动等效于(50± 13) Hz 的带通滤波过程， 根据式(4)—(7)可以得到基波分量在每个采样时刻的特征参数。类似地，将 TFA 的中心频率分别设置为 1、 3 Hz、直至采样频率的 1/2(在实际应用中可以根据需要灵活调整范围)， 调整尺度参数使频率窗口半径为 1 Hz，应用时频原子变换可以得到谐波和间谐波分量在每个采样时刻的特征参数 。

TFT 不仅适用于单相电压、电流波形的监测，也能应用于对称分量(序分量)情况，具有很好的适应能力和很高的测量精度，适合于电能扰动在线监测 。

自 组 织 映 射 网 络 (self-organizing map ，SOM)模拟大脑神经系统自组织映射的功能，是一种能无监督地进行自组织学习的人工神经网络。 SOM 能够将高维模式特征包含的复杂非线性关系转换成 2 维平面上的简单几何关系，并将输入的相似模式特征汇聚在输出层特定区域进行直观表达。不同于常规 SOM， ISOM 是有监督学习的人工神经网络，即学习过程中训练样本所属类别是已知的；输出层平面被划分为互不重叠的多个区域，分别对应于训练样本中的不同模式特征。训练完成后，输入测试样本，可根据其在输出层中最佳匹配点(best match unit， BMU)所属区域确定测试样本的模式类别。与常规 SOM 相比， ISOM 能够提高模式分类的准确性 。

将 ISOM 作为分类器实现电能质量事件的分类识别。首先，采用时频原子变换方法提取的典型电能质量事件的模式特征作为训练样本来训练ISOM，使其输出层的不同区域对应于不同的模式特征；进一步，采用测试样本对训练好的 ISOM 进行测试；最后，将其作为分类器识别电能扰动波形单一或同时存在的多个电能质量事件，并在输出层直观表达这些电能质量事件各自的发生、变化和结束过程，实现电能质量事件的直观表达和准确分类/识别 。

随着大量电力电子和精密仪器的使用，电能质量问题越来越引起人们的关注。准确、完善的电能质量监测是研究和治理电能质量的前提条件。小波变换具有良好的时频局部化特性，能够同时提取信号的时频特性，是一种良好的时频分析工具，已被广泛地应用于电能质量监测领域 。

电力系统中不可避免地存在着噪声，电气信号往往被淹没在大量的噪声信号中。这些噪声会降低小波的检测性能，甚至会造成检测失效。在噪声环境中及时准确地提取各种扰动信息已经成为电能质量监测的一个重要问题。

电能作为一种高效、清洁、可控的能源，已经得到最广泛的应用。但随着科技的发展，电力系统中用电负荷类型的改变和非线性负载的使用使得波形畸变、电能质量下降；而带有基于微处理机的控制器和功率电子器件的现代用电设备的使用，又提高了对电能质量的要求。因此，电能质量的研究成为一个热点问题 。

电能质量问题主要为电能质量扰动和谐波问题。电能质量扰动主要包括短时电压波动(电压暂升、电压暂降和中断)、电压下陷、电磁暂态脉冲和振荡。在电能质量检测中，常用的方法包括傅里叶变换法、小波变换法、数学形态学等。傅里叶变换具有计算量小、适用范围广、计算稳定可靠等优点，因此得到了广泛的应用，但傅里叶变换在应用中也经常遇到一些问题，如对工频稳态分量进行处理时，大大削弱了对畸变信号的检测能力，无法满足现代电力系统快速保护的要求。由于实际应用中不可能用无限窗，故存在频谱泄露等问题。小波分析法是近些年兴起的一种方法， 用小波分析方法可准确检测出信号的奇异点，从而检测出扰动的位置，但其计算十分复杂、实时性不强，难以用价格低廉的硬件实现。数学形态学是近年来提出的一种从图像处理演变而来的新方法，它利用图像处理理论提取信号的主要特征，而不改变其大体形状。这种基于时域的非线性数学方法具有计算简单、实时性强的特点，有较好的应用前景 。

小波包变换能够实现信号频带的均匀划分，在任意频率聚焦，是分析暂态电能质量扰动时频特性的良好工具。但是电气信号中的电磁噪声严重影响了小波包的检测特性 。

有文献提出了一种基于小波变换的除噪算法，能够在噪声环境中检测各种扰动信号。大量的仿真实验证明，该算法除噪效果良好。有文献指出，当信噪比 SNR 降低到 35dB 时，监测的准确度保持在 74%。但是作者在应用上述算法时，发现该算法存在 2 个问题：① 当信噪比 SNR 大于 35dB时，检测精度高，但当信噪比降低至 30dB 时，检测的准确度大幅下降。② 当干扰噪声是白噪声时， 检测效果良好， 但是当干扰噪声是有色噪声时，检测精度大幅下降。小波包变换建立在小波变换的基础上，可以实现信号频带的均匀划分，具有更好的时频特性 。

改进的 Brownian 桥式经验公式

设置参数 c 的目的是把小波包分解系数中的扰动分量与噪声分量分离。 c 的合理设置是选取阈值λ 的关键。对于噪声信号，如果令 c≥dmax，则可以实现信号与噪声分量的理想分离。因此，可以把噪声信号的 dmax看作是理想的参数值 c。在此基础上，对上述算法进行分析和改进 。

改进的概率公式

由表 1 可以看出，参数 c 虽然与噪声的方差无关，但是与采样点数 N 有关。有文献中利用概率的观点求出 c N a = −1/ 2 ln / 2 。这样设定 c 值在小波除噪算法中得到较为满意的检测结果。但是由于小波包分解系数与小波分解系数的分布不同，这样的 c 值不适用于小波包除噪。对具有相同分布特性，而采样点数不同的噪声信号进行分析，得到采样点数与 dmax 的关系曲线。其中， d1max=dmax×N。可以看出， d1max曲线的幅值保持在[160， 210]之间，而且不随采样点数的变化而变。因此，可以判断，dmax与采样点数成反比。因为可以把噪声信号的 dmax看作是理想的参数值 c，因此可以设c=dmax=k/N。

当信号发生扰动时，其 WPT(xn)的绝对值会远大于白噪声信号，因此可以通过设置阈值λ来实现噪声与扰动信号的分离。如果| (W x PT n ′ > λ) |，则认为是扰动信号分量，保留；如果| (W x PT n ′ < λ) |，则认为是噪声信号，不保留。 λ 的合理设置是信噪分离的关键 。

小波包除噪算法中阈值λ 的设置可以借鉴小波除噪的阈值设置算法。有文献给出了一种小波除噪的阈值选取算法，有较好的除噪效果， 介绍如下：令叠加了噪声的信号在尺度 n 的小波变换系数为 Cn，并假设

H1: C1=C2=…=CN

H2: C1=C2=…=Cn≠Cn 1≠…=CN-1=CN

假设 H1 正确，则小波变换就表示这是不含有扰动的检测信号，即只包含有平稳信号或噪声。另一方面，假设 H2 正确，则扰动有可能发生在 n 1的位置 。

小波变换阈值λ 的设置步骤如下：（ 1）验证假设的 H1；（ 2） H1 不成立，而 H2 成立，除去绝对值最大的小波变换系数，令 N=N-1 重新回到第（ 1）步。否则，进行第（ 3）步。（ 3） H1 成立，绝对值最大的小波变换系数就设置为阈值λ ；（ 4）用阈值λ 处理原来的小波变换系数。为验证假设的 H1，利用 Brownian 桥式经验公式得到一个序列 B(n/N) 。

如果把上述算法中的小波变换系数 Cn 用小波包变换系数 C[m n]（ m 为小波包分解层数， n 为小波包分解的结点位置）来替代，就得到了小波包除噪的阈值选取算法。由于小波变换系数和小波包变换系数的分布不完全一样，这样直接得到的小波包阈值选取算法除噪效果不理想，因此，需要对算法进行改进 。