一种基于变分模态分解多尺度排列熵的电能质量分析方法

摘要

本发明公开了一种基于变分模态分解多尺度排列熵的电能质量分析方法，包括：采集电能质量监测点在受到不同类型的扰动时的原始训练数据；采用变分模态分解对其分解得到K个IMF分量；计算每个IMF分量的多尺度排列熵，构建原始训练数据的特征向量；从特征向量中选择R个特征，组成原始训练数据的优化特征向量；以原始训练数据的优化特征向量和相应的扰动类型分别作为输入和输出数据，训练ELM神经网络模型，得到电能质量扰动分类器；对待检测电能质量监测点的电压信号，按上述相同方法获取其优化特征向量，并输入到电能质量扰动分类器中，得到待检测电能质量监测点的扰动类型。本发明方法对电能质量的扰动类型进行实时诊断的效率大大提升。

权利要求

1.一种基于变分模态分解多尺度排列熵的电能质量分析方法，其特征在于，包括以下步骤：步骤1，采集电能质量监测点在受到不同扰动类型的扰动时的原始训练数据；

采集电能质量监测点在历史时间段内受到扰动的电压信号U，并进行预处理，将得到的电压时序数据作为原始训练数据Ut；其中，针对每种扰动类型均获取若干个原始训练数据；

步骤2，采用变分模态分解对原始训练数据进行分解，每个原始训练数据均得到K个IMF分量uk(t)；

步骤3，计算每个IMF分量在NS个不同尺度因子下的排列熵，将每个IMF分量在每个尺度因子下的排列熵均作为原始训练数据的特征，由所有特征构建原始训练数据的特征向量V＝[v1,v2,...,vQ]，Q＝K×NS；

步骤4，从特征向量V中选择R个特征p1,p2,…,pR，组成原始训练数据的优化特征向量P＝[p1,p2,…,pR]；

步骤5，以原始训练数据的优化特征向量和相应的扰动类型分别作为输入数据和输出数据，训练ELM神经网络模型，得到电能质量扰动分类器；

步骤6，对待检测电能质量监测点的电压信号，按步骤1-4进行处理，得到待检测电能质量监测点的电压信号的优化特征向量，并输入到步骤5得到的电能质量扰动分类器中，电能质量扰动分类器输出待检测电能质量监测点的扰动类型。

2.根据权利要求1所述的方法，其特征在于，步骤4采用二进制遗传算法从原始训练数据的特征向量V中选择R个特征，具体过程为：步骤4.1，生成初始种群；

设置种群规模为MGA，最大迭代次数为d，并将种群中的每个个体均采用0/1编码方法初始化为1个字符串长度为Q的基因编码，每个字符与特征向量中的每个特征一一对应，字符为1表示选取相应位的特征，字符为0表示摈弃相应位的特征；

步骤4.2，确定适应度函数；

对于种群中的每个个体，将该个体所对应的特征作为电能质量扰动分类器的输入，根据相应的真实扰动类型计算识别准确率，将得到的识别准确率作为该个体的适应度函数；

步骤4.3，生成新种群；

计算当前种群的平均适应度值，并根据平均适应度值对当前通过对种群进行选择、交叉、变异，产生下一代种群；其中，所述种群的平均适应度值，通过对种群中所有个体的适应度值取平均计算得到；

步骤4.4，选择特征；

当满足最大迭代次数时，结束迭代，得到最终种群并选择最优个体，按最优个体从特征向量V中选择特征；其中最优个体是指种群中适应度值最高的个体。

3.根据权利要求1所述的方法，其特征在于，K个IMF分量为分别紧密围绕在各自中心频率ωk附近的调幅调频信号：Ak(t)和 分别为IMF分量uk(t)的瞬时幅值和瞬时相位，瞬时相位 的导数为瞬时频率ωk(t)；

步骤2的具体过程为：

步骤2.1，第一代参数初始化；

初始化第一代IMF分量为 第一代中心频率为

第一代拉格朗日算子λ1＝0；

步骤2.2，迭代更新IMF分量；

采用以下公式迭代更新IMF分量、中心频率和拉格朗日算子：其中， 为第n代IMF分量 的频域表示， 为第n代拉格朗日算子λn(t)的频域表示，x表示频率； 为原始训练数据Ut的频域表示；τ表示噪声容限参数；

步骤2.3，重复步骤2.2以执行迭代循环，当满足以下公式所述结束条件时，结束迭代循环：步骤2.4，将得到的K个IMF分量由频域 变换到时域uk(t)。

4.根据权利要求1所述的方法，其特征在于，步骤3中计算IMF分量在尺度因子为s下的排列熵的方法为：步骤3.1，时间序列粗粒化；

采用尺度因子s对IMF分量uk(t)按以下公式进行粗粒化处理，得到粗粒化后的时间序列式中， 表示N向下取整于s，N表示IMF分量的时间序列长度；

步骤3.2，时间序列相空间重构；

对步骤3.1得到的粗粒化后的时间序列 进行相空间重构，得到以下公式所示的多维度重构向量组：式中，m为嵌入维度，α为延迟时间；

步骤3.3，计算时间序列排列熵；

对于多维度重构向量组中的任一个向量 对其中m个元素根据数值大小进行升序排列得到重构向量，使用ID1,ID2,...,IDm表示重构向量中各元素所在列的索引，每一个向量 均得到对应一组符号序列： 多维度重构向量组中的所有向量共有g种序号排列方式，对应g组符号序列，其中g≤m！，m！表示m的阶乘；

计算多维度重构向量组的g组符号序列 出现的概率分别得到P1,P2,...,Pg，则多维度重构向量组的排列熵为：对多维度重构向量组的排列熵进行归一化处理，得到IMF分量在尺度因子为s下的排列熵：

5.根据权利要求1所述的方法，其特征在于，ELM神经网络模型的输入层的神经元数量为R+1，隐含层的神经元数量为14，输出层的神经元数量为7；对ELM神经网络模型的训练采用BP反向传播算法，最大迭代次数设置为1000，训练学习率为0.1，训练ELM神经网络模型的连接权值和阈值。

6.根据权利要求1所述的方法，其特征在于，扰动类型包括：应电压骤升、电压骤降、电压中断、电压振荡、电压缺口、谐波和无扰动。

说明书

一种基于变分模态分解多尺度排列熵的电能质量分析方法

技术领域

[0001] 本发明属于电能质量检测领域，具体涉及一种基于变分模态分解多尺度排列熵的电能质量分析方法。

背景技术

[0002] 随着科学技术的发展，电力负荷的变化和非线性负荷的增加，电网经常受到各种干扰的影响，从而导致波形失真和电能质量问题日益突出。电网中存在的电压暂升、暂降，谐波，衰减振荡以及电压闪变等电能质量问题得到电力部门和电网用户的广泛关注。电能质量是对电力系统中电能的质量的描述，用来衡量电能品质好坏。一般而言，电能的时序波形应该是幅值稳定，频率一致的稳定正弦波，但由于各类一次设备和二次设备对电力系统的影响，造成系统不对称性和非线性，以及各种自然或非自然因素导致的附加干扰及各类故障等原因，电能波形很难保持理想状态。所以电能质量的监测对于电力系统的稳定运行及故障预警有着极为重要的作用。

[0003] 要想实现电能质量的监测分析，最为重要的就是实现对电能质量暂态扰动进行检测与分类，而其中的首要步骤就是通过对原始电压进行处理分析提取电能质量扰动特征量。目前主流的特征量提取与分析的方法主要有：数学形态学算法、傅里叶变换、小波变换、S变换、EMD分解等，一系列的研究表明这些提取方法提取的特征量对于电能质量扰动具有一定的表征性，但同时也具有对噪声敏感，缺乏自适应性，计算过程复杂、效率低，端点效应严重，分解结果模态混叠严重，提取特征用于分类时分类精度低等问题。

[0004] 而另一方面，现有研究采用的分类手段主要为模板匹配或各类传统神经网络，缺乏优化，输入特征过多时识别模型的复杂度过高，分类器的效率低下。同时，现有研究难以实现电能质量的实时在线测量，以及分类器根据目标环境的不同做出相应调整。

发明内容

[0005] 为了在现有研究的基础上对电能质量暂态扰动识别方法进行改进，本发明提供一种基于变分模态分解多尺度排列熵的电能质量分析方法，采用基于VMD多尺度排列熵作为提取的特征训练基于ELM神经网络模型的电能质量扰动分类器，对不同的电能质量暂态扰动进行分类诊断，效率大大提升。

[0006] 为实现上述技术目的，本发明采用如下技术方案：

[0007] 一种基于变分模态分解多尺度排列熵的电能质量分析方法，包括以下步骤：

[0008] 步骤1，采集电能质量监测点在受到不同扰动类型的扰动时的原始训练数据；

[0009] 采集电能质量监测点在历史时间段内受到扰动的电压信号U，并进行预处理，将得到的电压时序数据作为原始训练数据Ut；其中，针对每种扰动类型均获取若干个原始训练数据；

[0010] 步骤2，采用变分模态分解对原始训练数据进行分解，每个原始训练数据均得到K个IMF分量uk(t)；

[0011] 步骤3，计算每个IMF分量在NS个不同尺度因子下的排列熵，将每个IMF分量在每个尺度因子下的排列熵均作为原始训练数据的特征，由所有特征构建原始训练数据的特征向量V＝[v1,v2,...,vQ]，Q＝K×NS；

[0012] 步骤4，从特征向量V中选择R个特征p1,p2,…,pR，组成原始训练数据的优化特征向量P＝[p1,p2,…,pR]；

[0013] 步骤5，以原始训练数据的优化特征向量和相应的扰动类型分别作为输入数据和输出数据，训练ELM神经网络模型，得到电能质量扰动分类器；

[0014] 步骤6，对待检测电能质量监测点的电压信号，按步骤1-4进行处理，得到待检测电能质量监测点的电压信号的优化特征向量，并输入到步骤5得到的电能质量扰动分类器中，电能质量扰动分类器输出待检测电能质量监测点的扰动类型。

[0015] 本方案通过采集电能质量监测点的电压时序数据，并采用变分模态分解及多尺度排列熵计算，提取并选择有效特征量，从而建立一种电能质量扰动分类器，实现对电能质量的实时监控。

[0016] 本方案所采用的变分模态分解作为一种新颖的信号分解处理方法，与传统的经验模态分解等递归式分解方法相比，能够根据预设数量K对电能质量监测点的电压时序数据进行分解得到K个IMF分量，具有受噪声影响较小，模态混叠不明显，分解误差小等特征，且能够有效分离频率接近的纯谐波信号，即IMF分量uk(t)，使得基于各纯谐波信号提取的特征所训练的电能质量扰动分类器的分类诊断效率大大提升。

[0017] 同时，该方案采用了计算多尺度排列熵的方式提取特征，多尺度排列熵是一种度量时间序列复杂性方法，具有通用性好，计算简单快捷、抗噪声干扰能力强等优点，适用于本发明机械故障诊断领域对非线性、非平稳时间序列的电压信号提取特征以进行机械故障诊断。

[0018] 进一步地，步骤4采用二进制遗传算法从原始训练数据的特征向量V中选择R个特征，具体过程为：

[0019] 步骤4.1，生成初始种群；

[0020] 设置种群规模为MGA，最大迭代次数为d，并将种群中的每个个体均采用0/1编码方法初始化为1个字符串长度为Q的基因编码，每个字符与特征向量中的每个特征一一对应，字符为1表示选取相应位的特征，字符为0表示摈弃相应位的特征；

[0021] 步骤4.2，确定适应度函数；

[0022] 对于种群中的每个个体，将该个体所对应的特征作为电能质量扰动分类器的输入，根据相应的真实扰动类型计算识别准确率，将得到的识别准确率作为该个体的适应度函数；

[0023] 步骤4.3，生成新种群；

[0024] 计算当前种群的平均适应度值，并根据平均适应度值对当前通过对种群进行选择、交叉、变异，产生下一代种群；其中，所述种群的平均适应度值，通过对种群中所有个体的适应度值取平均计算得到；

[0025] 步骤4.4，选择特征；

[0026] 当满足最大迭代次数时，结束迭代，得到最终种群并选择最优个体，按最优个体从特征向量V中选择特征；其中最优个体是指种群中适应度值最高的个体。

[0027] 本方案通过采用二进制遗传算法排除特征向量中比较无用的特征，提升后续电能质量扰动分类器的分类效率，选择有效特征则可提高电能质量扰动分类器的分类准确性，提升整体系统的运行效率。通过选用电能质量扰动分类器的识别准确率作为二进制遗传算法的适应度评价指标，相当于以电能质量扰动分类器的识别准确率为优化目标，来选择合适的特征，从而可进一步提高电能质量扰动分类器的分类准确率。

[0028] 进一步地，K个IMF分量为分别紧密围绕在各自中心频率ωk附近的调幅调频信号：

[0029]

[0030] Ak(t)和 分别为IMF分量uk(t)的瞬时幅值和瞬时相位，瞬时相位 的导数为瞬时频率ωk(t)；

[0031] 步骤2的具体过程为：

[0032] 步骤2.1，第一代参数初始化；

[0033] 初 始化 第一 代I MF 分量 为 第一 代中 心频 率 为第一代拉格朗日算子λ1＝0；

[0034] 步骤2.2，迭代更新IMF分量；

[0035] 采用以下公式迭代更新IMF分量、中心频率和拉格朗日算子：

[0036]

[0037]

[0038]

[0039] 其中， 为第n代IMF分量 的频域表示， 为第n代拉格朗日算子λn(t)的频域表示，x表示频率； 为原始训练数据Ut的频域表示；τ表示噪声容限参数；

[0040] 步骤2.3，重复步骤2.2以执行迭代循环，当满足以下公式所述结束条件时，结束迭代循环：

[0041]

[0042] 步骤2.4，将得到的K个IMF分量由频域 变换到时域uk(t)。

[0043] 进一步地，步骤3中计算IMF分量在尺度因子为s下的排列熵的方法为：

[0044] 步骤3.1，时间序列粗粒化；

[0045] 采用尺度因子s对IMF分量uk(t)按以下公式进行粗粒化处理，得到粗粒化后的时间序列

[0046]

[0047] 式中， 表示N向下取整于s，N表示IMF分量的时间序列长度；

[0048] 步骤3.2，时间序列相空间重构；

[0049] 对步骤3.1得到的粗粒化后的时间序列 进行相空间重构，得到以下公式所示的多维度重构向量组：

[0050]

[0051] 式中，m为嵌入维度，α为延迟时间；

[0052] 步骤3.3，计算时间序列排列熵；

[0053] 对于多维度重构向量组中的任一个向量 对其中m个元素根据数值大小进行升序排列得到重构向量，使用ID1,ID2,...,IDm表示重构向量中各元素所在列的索引，每一个向量 均得到对应一组符号序列：
多维度重构向量组中的所有向量共有g种序号排列方式，对应g组符号序列，其中g≤m！，m！表示m的阶乘；

[0054] 计算多维度重构向量组的g组符号序列 出现的概率分别得到P1,P2,...,Pg，则多维度重构向量组的排列熵为：

[0055]

[0056] 对多维度重构向量组的排列熵进行归一化处理，得到IMF分量在尺度因子为s下的排列熵：

[0057]

[0058] 进一步地，ELM神经网络模型的输入层的神经元数量为R+1，隐含层的神经元数量为14，输出层的神经元数量为7；对ELM神经网络模型的训练采用BP反向传播算法，最大迭代次数设置为1000，训练学习率为0.1，训练ELM神经网络模型的连接权值和阈值。

[0059] 进一步地，扰动类型包括：应电压骤升、电压骤降、电压中断、电压振荡、电压缺口、谐波和无扰动。

[0060] 有益效果

[0061] 本方案通过采集电能质量监测点的电压时序数据，并采用变分模态分解及多尺度排列熵计算，提取有效特征量，并导入基于遗传算法优化的ELM神经网络模型中分类，从而建立了一种电能质量扰动分类器，实现对电能质量的实时监控。

[0062] 本方案所采用的变分模态分解作为一种新颖的信号分解处理方法，与传统的经验模态分解等递归式分解方法相比，能够根据预设模态分量个数对电能质量监测点的电压时序数据进行分解，具有受噪声影响较小，模态混叠不明显，分解误差小等特征，且能够有效分离频率接近的纯谐波信号，使得基于各纯谐波信号提取的特征所训练的电能质量扰动分类器的分类诊断效率大大提升。

[0063] 同时，该方案采用了计算多尺度排列熵的方式提取特征，多尺度排列熵是一种度量时间序列复杂性方法，具有通用性好，计算简单快捷、抗噪声干扰能力强等优点，适用于任何类型的非线性、非平稳时间序列，广泛运用于生物医学信号研究，机械故障诊断等领域。最后，本方案采用遗传算法对特征进行选择，优化了电能质量扰动分类器的分类准确性，提升了整体系统的运行效率。

附图说明

[0064] 图1为本发明所述方法的流程示意图。

具体实施方式

[0065] 下面对本发明的实施例作详细说明，本实施例以本发明的技术方案为依据开展，给出了详细的实施方式和具体的操作过程，对本发明的技术方案作进一步解释说明。

[0066] 如图1所示，本发明提供一种基于变分模态分解多尺度排列熵的电能质量分析方法，通过提取电信号多频段多尺度排列熵信息，采用神经网络分类模型，实现对电能质量的实时监测与故障类型判断。主要内容包括以下步骤：

[0067] 步骤1，采集电能质量监测点在受到不同类型的扰动时的原始训练数据；

[0068] 采集电能质量监测点在历史时间段内受到暂态扰动的电压信号U，并进行预处理，将得到的电压时序数据作为原始训练数据Ut；其中，针对每种扰动类型均获取若干个原始训练数据；

[0069] 具体地，采用电压互感器采集电网中电能质量监测点的电压信号U，作为电能质量扰动分析的输入数据，由于电压互感器采样频率较高，采用原始的采样电压数据进行后续处理的计算成本较大，效率低下，故本发明对原始的10kHz采样频率进行预处理，通过求取固定步长内的电压有效值将10kHz的原始数据处理为基频信号频率50Hz的电压时序数据，并按照滑动窗的方式，以步长为80个周期基频信号长度向前移动，每次移动后窗口截取100个周期基频信号作为系统所需的训练数据，获得原始训练数据Ut。

[0070] 步骤2，采用变分模态分解对原始训练数据进行分解，得到K个IMF分量uk(t)，为分别紧密围绕在各自中心频率ωk附近的调幅调频信号：

[0071] 在语音处理、风速预测等领域中，经验模态分解(Empirical  Mode Decomposition,EMD)广泛的应用于时间序列的分解和分析。但是，由于它缺乏理论定义和递归结构，经验模态分解(EMD)不具有反向纠错能力，高度依赖于极值点的搜寻方法、载波包络线的插值和终止条件，使其具有模态混叠、端点效应、分解停止准则等问题。而变分模态分解(Variational Mode Decomposition,VMD)采用变分模型自适应地确定相关频段，对含有噪声的采样数据进行处理时具有更好的鲁棒性。故本发明采用变分模态分解方法对步骤1获取到的原始训练数据进行分解，使得到的K个IMF分量的带宽之和最小，降低噪声对电能质量扰动类型检测分类的干扰，提高分类精确。

[0072] 具体地，本发明采用变分模态分解对原始训练数据进行分解的过程为：

[0073] 步骤2.1，第一代参数初始化；

[0074] VMD是通过自适应的将信号Ut的分解成K个紧密围绕在各自中心频率ωk附近的IMF分量uk(t)，分解后的每一IMF分量uk(t)被重新定义为一个调幅-调频信号，表示如下：

[0075]

[0076] 其中Ak(t)和 分别为IMF分量uk(t)的瞬时幅值和瞬时相位，瞬时相位 的导数为瞬时频率ωk(t)。每个IMF分量对应一个中心频率，而每个IMF分量的瞬时频率ωk(t)与时间有关，本发明只用到了各IMF分量的中心频率ωk。

[0077] 由于VMD是采用迭代的方式更新每一层IMF的uk(t)，故初始化第一代调幅-调频信号uk(t)为 及第一代中心频率 同时初始化第一代拉格朗日算子λ1＝0及迭代数n＝1。

[0078] 步骤2.2，迭代更新IMF分量；

[0079] 采用以下公式迭代更新IMF分量、中心频率和拉格朗日算子：

[0080]

[0081]

[0082]

[0083] 其中，， 为第n代IMF分量 的频域表示， 为第n代拉格朗日算子λn(t)的频域表示，x表示频率； 为原始训练数据Ut的频域表示；τ表示噪声容限参数；

[0084] 步骤2.3，重复步骤2.2以执行迭代循环，当满足以下公式所述结束条件时，结束迭代循环：：

[0085]

[0086] 步骤2.4，将得到的K个IMF分量由频域 变换到时域uk(t)。

[0087] 通过傅里叶反变换得到时域内K个带宽之和最小的IMF分量，即VMD分解结果，并将分解后每个分量下的时间序列表示为uk(t)。

[0088] 步骤3，提取多尺度排列熵信息；

[0089] 多尺度排列熵(multi-scale permutation entropy，MPE)指的是在不同尺度下的排列熵，对原时间序列在多个尺度上进行粗粒化处理，然后计算时间序列在每个尺度粗粒化后的排列熵值。

[0090] 通过计算每个IMF分量在NS个不同尺度因子下的排列熵，将每个IMF分量的每个尺度因子下的排列熵均作为原始训练数据的特征，由所有特征构建原始训练数据的特征向量V＝[v1,v2,...,vQ]，Q＝K×NS。

[0091] 其中，计算IMF分量在尺度因子为s下的排列熵的具体过程为：

[0092] 步骤3.1，时间序列粗粒化；

[0093] 采用尺度因子s对IMF分量uk(t)按以下公式进行粗粒化处理，即使用长度为s的窗口对原时间序列按序划分，然后计算每个时间窗内的均值，得到粗粒化后的时间序列[0094]

[0095] 式中， 表示N向下取整于s，N表示IMF分量的时间序列长度；由于本实施例在步骤1中截取的是100个周期的基频信号作为原始训练数据，故此处N＝100。

[0096] 步骤3.2，时间序列相空间重构；

[0097] 对步骤3.1得到的粗粒化后的时间序列 进行相空间重构，得到以下公式所示的多维度重构向量组：

[0098]

[0099] 式中，m为嵌入维度，α为延迟时间；

[0100] 步骤3.3，计算时间序列排列熵；

[0101] 对于多维度重构向量组中的任一个向量 对其中m个元素根据数值大小进行升序排列得到重构向量，使用ID1,ID2,...,IDm表示重构向量中各元素所在列的索引，每一个向量 均得到对应一组符号序列：
多维度重构向量组中的所有向量共有g种序号排列方式，对应g组符号序列，其中g≤m！；

[0102] 计算多维度重构向量组的g组符号序列 出现的概率P1,P2,...,Pg，则多维度重构向量组的排列熵为：

[0103]

[0104] 多维度重构向量组的排列熵进行归一化处理，得到第k个IMF分量在尺度因子为s下的排列熵：

[0105]

[0106] MPE值得大小表示时间序列的随机程度，即在本实施例中，MPE表示IMF分量uk(t)在尺度因子为s时的随机程度，其值越小说明时间序列越规则，反之时间序列越接近随机。

[0107] 最后，计算每一IMF分量下的多尺度排列熵值作为特征，构建特征向量V＝[v1,v2,...,vQ]，Q＝K×NS。其中K和NS分别表示VMD分解后的IMF分量总数和时间序列粗粒化处理的尺度数。

[0108] 步骤4，从特征向量V中选择R个特征p1,p2,…,pR，组成原始训练数据的优化特征向量P＝[p1,p2,…,pR]；

[0109] 在将特征向量导入分类模型之前需要对特征向量进行特征选择，提取出其中的有效特征，从而提升分类器的分类准确性。本发明采用二进制遗传算法对其进行特征选择，具体步骤如下：

[0110] 步骤4.1，生成初始种群；

[0111] 设置种群规模为MGA，最大迭代次数为d，并将种群中的每个个体均采用0/1编码方法随机初始化为1个字符串长度为Q的基因编码，基因编码中每个字符与特征向量中的每个特征一一对应，字符为1表示选取相应位的特征，字符为0表示摈弃相应位的特征；

[0112] 步骤4.2，确定适应度函数；

[0113] 对于种群中的每个个体，将该个体所对应的特征作为电能质量扰动分类器的输入，根据相应的真实扰动类型计算识别准确率，将得到的识别准确率作为该个体的适应度函数，适应度函数值越大则代表该个体越优秀。

[0114] 步骤4.3，生成新种群；

[0115] 计算当前种群的平均适应度值，并根据平均适应度值通过对当前种群进行选择、交叉、变异，产生下一代种群；其中，所述种群的平均适应度值，通过对种群中所有个体的适应度值取平均计算得到；

[0116] 步骤4.4，选择特征；

[0117] 当满足最大迭代次数时结束迭代，得到最终种群并选择最优个体，按最优个体从特征向量V中选择特征，得到优化特征向量P＝[p1,p2,...,pR]。其中最优个体是指种群中适应度值最高的个体，R为最优个体中为1的字符个数。

[0118] 其中，选择是指以一定概率选择较优的父代复制并遗传至下一代；交叉是指以一定概率把两个父代的基因进行交叉重组，获得新的下一代；变异是指以一定概率选择父代中的某个基因进行更改。

[0119] 步骤5，以原始训练数据的优化特征向量和相应的扰动类型分别作为输入数据和输出数据，训练ELM神经网络模型，得到电能质量扰动分类器。具体过程为：

[0120] 步骤5.1，初始化ELM参数

[0121] ELM神经网络模型的输入层的神经元数量，根据步骤4选择所得的特征的数量R而得到，具体根据不同的场景而有所变动，输入向量Pi为：Pi＝[1,p(1),p(2),…,p(R)]。隐含层的神经元数量为14，隐含层输出的向量Ph为：Ph＝[1,ph(1),ph(2),ph(3)...ph(14)]。输出层的神经元数量，根据扰动类型的数量而设为7，输出向量Po为：Po＝[po(1),po(2),po(3)...po(7)]，分别对应电压骤升(swell)，电压骤降(sag)，电压中断(interruption)，电压振荡(oscillatory transients)，电压缺口(notch)，谐波(harmonic)六种事件，及无扰动事件。

[0122] 其中，定义输入层与隐含层之间的权值矩阵A如下：

[0123]

[0124] 定义隐藏层与输出层之间的权值矩阵B如下：

[0125]

[0126] 权值矩阵A和B中的第一列分别代表着输入层和隐含层的阈值，将其与权值相融合有助于后续计算机运算。

[0127] 选择广泛使用的Sigmoid函数作为各神经元的激活函数，其表达式如下：

[0128]

[0129] 步骤5.2，ELM神经网络模型训练及验证

[0130] 采用常见的BP反向传播算法对ELM神经网络分类模型进行训练，将最大迭代次数设置为1000，训练学习率设为0.1。

[0131] 在电能质量监测点，针对每种扰动类型均采集若干个电压信号，并对所有电压信号均按照上述步骤1-4将电压信号处理转化为带标签(已知扰动类型，将扰动类型作为标签)的优化特征向量，所有优化特征向量组成特征向量矩阵；然后将矩阵按3:1的比例分为训练集和测试集。训练集用于训练ELM分类器，并且用于上述二进制遗传算法选择特征的适应度值计算。测试集用于计算分类器模型的误差，评估分类器模型的精确性和鲁棒性。综合以上过程便得到了最终的电能质量在线故障诊断模型，便可用来执行步骤6对待检测电能质量监测点进行实时扰动检测，判断具体的暂态扰动的类型。

[0132] 步骤6，实时扰动监测；

[0133] 采集待检测电能质量监测点的电压信号，按步骤1-4进行处理，得到待检测电能质量监测点的电压信号的优化特征向量，并输入到步骤5得到的电能质量扰动分类器中，电能质量扰动分类器输出待检测电能质量监测点的扰动类型。

[0134] 以上实施例为本申请的优选实施例，本领域的普通技术人员还可以在此基础上进行各种变换或改进，在不脱离本申请总的构思的前提下，这些变换或改进都应当属于本申请要求保护的范围之内。