一种基于免疫RBF神经网络的电网故障诊断方法与流程

本发明属于电网故障诊断技术领域，具体涉及一种基于免疫rbf神经网络的电网故障诊断方法。

背景技术：

随着电网规模的不断扩大，在电网发生故障时，若不能及时准确地定位故障区域，处理故障设备，恢复供电，很可能造成更加严重的大面积停电事故并造成巨大的经济损失。当电网发生故障时，调度监测人员会在短时间内收到大量的故障报警信息，若仅靠人工处理这些故障信息并快速准确地判断故障区域是十分困难的。因此，及时准确诊断出故障元件对于电网恢复供电和稳定运行具有十分重要的意义。

电网故障诊断一般是指通过scada采集的保护和断路器动作的开关量变化信息来识别出故障的元件位置、类型等。目前，许多学者已对电网故障诊断领域进行了广泛的研究，并提出了许多智能诊断算法，如专家系统、优化技术、贝叶斯理论、petri网、神经网络等。rbf神经网络是一种以函数逼近为基础的前馈型神经网络，由于rbf神经网络具有结构简单、学习速度快、泛化能力强等特点，并且具有任意函数逼近能力，使之在电网故障领域获得了广泛的运用。然而rbf中心点的选择一直是rbf神经网络算法的关键之一，目前传统的rbf神经网络通常采用k-means聚类确定rbf中心点，但k-means算法对聚类的初始中心敏感，会造成局部极小值的出现，网络的收敛速度也会变慢，从而直接降低rbf神经网络的性能。因此有必要研究一种方法更好的确定rbf中心点，以此增强rbf神经网络的性能，提高在故障诊断时的容错性，从而更好的辅助调度人员进行事故的分析和处理。

技术实现要素：

发明目的：针对传统的rbf神经网络采用k-means聚类确定rbf中心点时，会造成局部极小值的出现，网络的收敛速度也会变慢，从而直接降低rbf神经网络的性能，降低在故障诊断时的容错性问题，本发明提供一种基于免疫rbf神经网络的电网故障诊断方法。

技术方案：一种基于免疫rbf神经网络的电网故障诊断方法，包括以下步骤：

(1)采集样本数据，并将样本数据划分为训练集和测试集；

(2)构建基于免疫rbf神经网络的电网故障网络结构，包括确定rbf神经网络输入向量、输出向量的维数，确定隐含层节点数；

(3)构建三层前向网络拓扑结构，第一层为输入层，由采集的保护、断路器报警数据构成；第二层为隐含层，采用高斯函数进行非线性变换；第三层为输出层，对隐含层进行加权线性组合完成线性映射，以其输出结果确定电网故障区域；

(4)网络训练，输入训练集，先根据免疫算法原理，对电网故障诊断问题进行可行性分析，提取故障先验知识，构造适应度函数，并制定各种约束条件，把训练样本中的电网报警数据作为抗原，rbf中心点作为抗体，在进行克隆、变异等免疫操作后，将得到的样本数据中心作为rbf神经网络的中心点，再确定扩展常数，并利用最小二乘法确定隐含层到输出层的权值，得到电网故障诊断模型；

(5)电网故障诊断模型，将测试集输入到训练好的电网故障诊断模型中进行故障诊断的效果验证；

(6)通过电网故障诊断模型进行电网故障区域定位。

进一步的，步骤(2)所述的rbf神经网络中输入向量的特征维数由判断故障区域的继电保护和断路器的数量决定；隐含层节点数量设置成与输入向量的特征维数相同，可根据电网规模和诊断模型测试效果进行适当调整。更进一步的，确定rbf神经网络输入向量、输出向量的维数，确定隐含层节点数。其中：输入层含有m个神经元(x1,x2,…,xm)，对应着输入向量的特征维数，由判断故障区域的继电保护和断路器的数量决定；隐含层含有n个神经元(h1,h2,…,hn)；输出层含有p个神经元(y1,y2,…yp)，对应于电网故障的p个区域。

步骤(3)所述的隐含层的激活函数为高斯函数，其表达式如下：

hi(x)＝exp(-||x-ci||²/qi²)

上式中：x＝[x1,x2,…,xm]^t为m维输入向量；ci＝[ci1,ci2,…,cim]^t为第i个隐含层神经元径向基函数的中心；||x-ci||表示向量x与向量ci之间的欧几里得距离；qi为为扩展常数，又称神经元的“宽度”，表示高斯函数的影响范围；i＝1,2,…,n。

步骤(4)基于免疫算法确定rbf神经网络的中心点具体步骤如下：

(41)首先进行抗原识别，将训练样本中的电网报警数据作为抗原，然后产生初始抗体群，把中心点的可行取值表示成样本空间中的抗体，从样本空间内随机产生一个初始种群；同时初始化免疫算法的其他参数，如种群刷新比例、免疫选择比例、最大进化迭代次数等。

(42)对抗体进行解码确定rbf神经网络中心点，采用最小二乘法确定rbf神经网络的输出权值，计算网络适应度，具体网络适应度计算公式如下：

式中，yk为网络输出层第k个节点的实际输出；为网络输出层第k个节点的期望输出；p为网络输出层的总节点个数。

(43)根据网络的适应度计算每个抗体的亲和力，具体计算公式如下：

式中，e(i)为种群中的第i个抗体的网络适应度。

(44)判断是否满足终止条件，若满足终止条件则终止算法寻优过程，输出结果；若不满足终止条件，则继续执行(5)。

(45)计算抗体的浓度和激励度，具体计算公式如下：

sim(i)＝a·aff(i)-b·den(i)

式中，den(i)为抗体的浓度；n为种群规模；δs为相似度阈值；bff(i,j)为抗体间的亲和力；h(i,j)为抗体i和抗体j间的距离；sim(i)为抗体的激励度；a和b为激励度系数。

(46)进行免疫选择、克隆、变异、克隆抑制。

免疫选择：根据亲和力和抗体浓度选择优质抗体，并使其活化。

克隆：对活化后的抗体进行复制，得到多个副本。

变异：对克隆后的副本进行变异，使其亲和度发生突变。

克隆抑制：对变异后的抗体进行选择，抑制亲和度低的抗体，保留亲和度高的抗体。

(47)种群刷新，按激励度大小对抗体进行排序，用随机新生成的抗体替代激励度较小的抗体，形成新的抗体种群，并跳转(42)。

更进一步的，步骤(4)所述的扩展常数由数据中心的散布确定，为避免径向基函数图像太尖或太平，按如下公式计算：

式中，dmax为所选聚类中心的最大距离；m为聚类中心的个数。

更进一步的，步骤(4)所述的隐含层到输出层的权值确定，采用最小二乘法计算隐含层到输出层的权值矩阵w，具体为：

隐含层的输出为矩阵b，对所有的训练集样本，神经网络的输出为：

则神经网络的实际输出与训练集的误差为：

利用最小二乘法使得ε最小化，由此可以求得隐含层到输出层的权值矩阵：

w＝b⁺y

b⁺＝(b^tb)^-1b^t

式中，b⁺为b的伪逆。

有益效果：与现有技术相比，本发明所述方法可以更加精确的确定rbf神经网络中基函数的中心点，有效提高了网络的收敛速度，改善了rbf神经网络的性能，使得所构建的电网故障诊断模型具有较强的容错性，能够有效地实现已知故障的诊断。

附图说明

图1为实施例中配电系统结构示意图；

图2为本发明所述的电网故障诊断方法的流程示意图；

图3为本发明利用免疫算法寻找rbf中心点的流程示意图。

具体实施方式

为详细的说明本发明所公开的技术方案，下面结合说明书附图和具体实施例做进一步的阐述。

图1所示的是某个配电系统，以此配电系统为例，该系统共分为5个区域(sec1～sec5),每个区域均配有断路器保护开关(qf1～qf5)和过电流保护(co1～co5),其中，sec1和sec3配有距离保护(rr1、rr3)。

本发明所述方法的流程如图2所示。具体实施步骤如下：

步骤1：样本选取

选取断路器保护、过电流保护和距离保护的状态作为判断电网故障区域的决策属性，将其作为神经网络输入，并将故障区域作为神经网络的输出，以此进行神经网络的训练。

将已知的电网故障样本数据分成训练集和测试集，训练集样本见表1，一共有13组样本。表2为含有错误信息的测试集样本，一共有6组样本，和表1的训练集样本对比可知，样本1和样本6是qf2断路器误动，样本2是co1过电流保护误动，样本3和样本5是qf3断路器拒动，样本4是qf4断路器拒动。

表1.故障训练集样本

表2.故障测试集样本

步骤2：网络结构设计

从故障样本可以看出每一组故障数据有12个判断条件，对应于神经网络的12个输入层节点，即可确定神经网络的输入向量的维数为12；故障样本中有5个故障区域和1个无故障区域，对应于神经网络的6个输出层节点，即可确定神经网络的输入向量的维数为6。神经网络的隐含层节点数取11。

步骤3：网络拓扑设计

选用高斯函数作为隐含层的激活函数，函数如下：

hi(x)＝exp(-||x-ci||²/qi²)

式中：x＝[x1,x2,…,xm]^t为m维输入向量；ci＝[ci1,ci2,…,cim]^t为第i个隐含层神经元径向基函数的中心；||x-ci||表示向量x与向量ci之间的欧几里得距离；qi为对应神经元的“宽度”，又称为扩展常数，表示高斯函数的影响范围；i＝1,2,…,n。

步骤4：网络训练

rbf神经网络的训练主要由训练样本来调整中心点、扩展常数、输出权值这三个参数，通过对这三个参数的调整，完成rbf神经网络的训练。

首先，采用免疫算法确定rbf神经网络的中心点，流程如图3所示，具体过程如下：

(41)首先进行抗原识别，将输入数据作为抗原，然后产生初始抗体群，把中心点的可行取值表示成样本空间中的抗体，从样本空间内随机产生一个初始种群。同时初始化免疫算法的其他参数，如种群刷新比例、免疫选择比例、最大进化迭代次数等；

(42)对抗体进行解码确定rbf神经网络中心点，采用最小二乘法确定rbf神经网络的输出权值，计算网络适应度，具体网络适应度计算公式如下：

式中，yk为网络输出层第k个节点的实际输出；为网络输出层第k个节点的期望输出；p为网络输出层的总节点个数。

(43)根据网络的适应度计算每个抗体的亲和力，具体计算公式如下：

式中，e(i)为种群中的第i个抗体的网络适应度。

(44)判断是否满足终止条件，若满足终止条件则终止算法寻优过程，输出结果；若不满足终止条件，则继续执行(5)。

(45)计算抗体的浓度和激励度，具体计算公式如下：

sim(i)＝a·aff(i)-b·den(i)

式中，den(i)为抗体的浓度；n为种群规模；δs为相似度阈值；bff(i,j)为抗体间的亲和力；h(i,j)为抗体i和抗体j间的距离；sim(i)为抗体的激励度；a和b为激励度系数。

(46)进行免疫选择、克隆、变异、克隆抑制。

免疫选择：根据亲和力和抗体浓度选择优质抗体，并使其活化。

克隆：对活化后的抗体进行复制，得到多个副本。

变异：对克隆后的副本进行变异，使其亲和度发生突变。

克隆抑制：对变异后的抗体进行选择，抑制亲和度低的抗体，保留亲和度高的抗体。

(47)种群刷新，按激励度大小对抗体进行排序，用随机新生成的抗体替代激励度较小的抗体，形成新的抗体种群，并跳转步骤(42)。

然后，由数据中心的散布确定扩展常数，为避免径向基函数太尖或太平，按如下公式计算：

式中，dmax为所选聚类中心的最大距离；m为聚类中心的个数。

最后，采用最小二乘法确定隐含层到输出层的权值矩阵w，具体为：

隐含层的输出为矩阵b，对所有的训练集样本，神经网络的输出为：

则神经网络的实际输出与训练集的误差为：

利用最小二乘法使得ε最小化，由此可以求得隐含层到输出层的权值矩阵：

w＝b⁺y

b⁺＝(b^tb)^-1b^t

式中，b⁺为b的伪逆。

步骤5：网络测试

将测试集输入到训练好的电网故障诊断模型中进行故障诊断的效果验证，为了验证本发明的有效性，在将测试集输入到本发明采用免疫算法确定中心点的故障诊断模型的同时，也将测试集输入到采用k-means算法确定中心点的故障诊断模型进行对比分析，诊断结果如表3所示。

表3.故障诊断结果

从表3故障诊断的结果来看，本发明采用免疫算法确定中心点的rbf神经网络要比传统的采用k-means算法确定中心点的rbf神经网络有更好的分类效果，传统的rbf神经网络在诊断样本2和样本4时出现诊断错误，而采用免疫算法确定中心点的rbf神经网络对6个样本的诊断均正确。

综上，本发明所公开的一种基于免疫rbf神经网络的电网故障诊断方法，采用免疫算法进行rbf中心点的确定，弥补了传统采用k-means聚类确定rbf中心点会出现局部极小值的问题，提高了rbf神经网络的性能。使用免疫rbf神经网络所构建的电网故障诊断模型具有较强的容错性，能够有效地实现已知故障的诊断，进而可以很好的辅助调度人员进行事故的分析和处理。