基于AO算法优化PNN的配电网线变关系精准辨识方法

技术领域

本发明属于配电网线变关系辨识技术领域，具体涉及一种基于AO算法优化PNN的配电网线变关系精准辨识方法。

背景技术

随着经济的快速发展，用户对能源的依赖程度越来越高，而电能作为不可替代的能源越来越凸显了其重要性，电网公司为了满足用户对电能的需求并且提高供电可靠性对配电网进行了大范围的改造，使得配电网的复杂度越来越高。

例如，在实际运用中，为了快速的改切负荷以及对停电区域的快速恢复送电，电网公司经常在配电网中安装大量的联络开关、断路器、环网柜等设备，当有需要时，只需对这些设备进行操作便可实现负荷改切以及恢复送电，但是由于在实际操作中，工作人员可能未及时将操作记录到GIS系统中，导致配电网实际运行的线变关系与GIS系统中的线变关系不一致。

目前对于配电网的线变关系识别主要采用有基于硬件的方式和基于数据驱动的方式。基于硬件的方式需要在配电网中安装大量的硬件设备，这种方式需要投入大量的财力去安装设备，而基于数据驱动的方式目前主要采用的是基于电压相关性的方式进行判断，这种方式的辨识阈值难以确定，导致辨识准确性难以得到保证。

发明内容

针对配电网目前的拓扑混乱，线变关系不清晰的难点，大量的专家学者提出了众多的方案进行辨识，但是所获得的方式并不能完美的解决所面临的难题，为了彻底解决这一难点，提升电网公司的运行服务质量。因此，本发明提出了一种基于AO（天鹰）算法优化PNN（概率神经网络）的配电网线变关系精准辨识方法，本方法通过提取配变的每日配变三相电压采样数据生成电压矩阵，首先通过电压值计算出配变之间的电压相关性矩阵以及配变和线路之间的相关性矩阵，其次将三个矩阵组合成为原始特征矩阵以后使用核主成分分析法（KPCA）进行降维获取新的特征矩阵，接着使用AO算法优化PNN网络，最后将新的特征矩阵作为输入优化后的PNN网络获取辨识结果。

本发明通过以下技术方案来实现。一种基于AO算法优化PNN的配电网线变关系精准辨识方法，包括以下步骤：

步骤1：提取线路下所有配变的三相电压数据并经过预处理后获取所有配变的电压矩阵；

步骤2：进行皮尔逊相关系数计算以获取配变间的电压相关性系数矩阵以及配变和线路之间的电压相关性系数矩阵；将配变间的电压相关性系数矩阵、配变和线路之间的电压相关性系数矩阵以及电压矩阵进行处理后合并组成配电网特征矩阵；

步骤3：使用核主成分分析法（KPCA）对归集后的配电网特征矩阵降维，得到降维后的配电网特征矩阵；

步骤4：使用AO算法优化PNN网络，使PNN网络的平滑因子的选取达到最优值；

步骤5：将实时监测的电压矩阵按步骤2和步骤3处理得到的降维后的配电网特征矩阵，送入到使用AO算法优化后的PNN网络进行分析以辨识发生线变关系错误的配变。

进一步优选，步骤4的具体过程如下：

步骤4.1：输入层，使用降维后的配电网特征矩阵作为PNN网络输入层的输入矩阵；

步骤4.2：模式层，计算输入矩阵和训练集之间的关系输出为：

式中，E表示输入矩阵，也是待识别样本，W

步骤4.3：求和层，将配变为正确或者错误的概率进行累加；

步骤4.4：输出层，输出每个神经元对应的配变的正确或者错误结论，表示为：

式中，

步骤4.5：使用AO算法优化PNN网络参数，选用均方根误差函数作为适应度函数进行计算：

式中，A

步骤4.6：初始化天鹰种群的位置，天鹰种群的初始位置为：

其中，

步骤4.7：更新天鹰位置，采用天鹰算法的垂直俯冲攻击、盘旋攻击、低空飞行和陆地行走捕食四种方法进行计算以实现位置的更新；

步骤4.8：再次使用适应度函数计算，判断是否满足适应度需求或者达到最大迭代次数，若不满足则返回步骤4.6，若满足则结束；

步骤4.9：输出满足条件的天鹰位置，通过解码以后得到最优平滑因子，天鹰的位置会随着AO算法的运行而不断进行迭代更新，直到满足所需条件；

步骤4.10：将使用AO算法得到的最优平滑因子送入到PNN网络进行计算。

进一步优选，步骤2的具体过程为：

步骤2.1：计算地理信息系统中保存的同一条线路下所有配变之间的皮尔逊相关系数组成配变间的电压相关性系数矩阵；

步骤2.2：计算所有配变和所属线路的皮尔逊相关系数，采用皮尔逊相关系数公式计算，组建配变和线路之间的电压相关性系数矩阵；

步骤2.3：分别将电压矩阵、配变间的电压相关性系数矩阵以及配变和线路之间的电压相关性系数矩阵进行归一化；

步骤2.4：将归一化后的电压矩阵、配变间的电压相关性系数矩阵以及配变间的电压相关性系数矩阵组合成原始配电网特征矩阵；将原始配电网特征矩阵的数据以配变为单位重新归集形成各台配变的特征矩阵，得到归集后的配电网特征矩阵。

进一步优选，步骤3的具体过程为：

步骤3.1：针对归集后的配电网特征矩阵，计算核矩阵，并将其中心化得到中心化核矩阵；

步骤3.2：计算中心化核矩阵的特征值和特征向量；

步骤3.3：将中心化核矩阵的特征值，并进行降序排序，取中心化核矩阵的前r个特征值和对应的特征向量；

步骤3.4：利用核矩阵和中心化核矩阵的前r个的特征值和对应的特征向量计算降维结果；

步骤3.5：对降维结果计算主成分累计贡献率；当前P个降维结果的主成分累计贡献率大于设定百分比时，选取前P个降维结果的特征值及特征向量组建成为降维后的配电网特征矩阵。

具体地，所述电压矩阵为：

式中，n表示每日的电压采样次数，m表示同一条线路下所辖配变数量，V

具体地，配变之间的皮尔逊相关系数计算方式为：

式中，ρ

具体地，所述配变间的电压相关性系数矩阵为：

式中，ρ

配变和线路之间的电压相关性系数矩阵为

本发明使用了核主成分分析法（KPCA）对采集到的非线性数据组合成的原始特征矩阵进行降维，在保留最大特征的基础下运用了PNN（概率神经网络）对配变的正确与否进行分析，在运用过程中，针对PNN的平滑因子

附图说明

图1是本发明的方法流程图。

图2是一条线路的线路拓扑图。

具体实施方式

下面结合附图和实施例对本发明进一步详细阐明。

参照图1，一种基于AO算法优化PNN的配电网线变关系精准辨识方法，包括以下步骤：

步骤1：提取线路下所有配变的三相电压数据并经过预处理后获取所有配变的电压矩阵；

步骤1.1：从配电网运行调度中心的生产管理系统（PMS系统）、数据采集与监视控制系统（SCADA系统）以及地理信息系统（GIS）中导出所需的配变的每日96次采样的电压数据如表1以及图2所示的线路拓扑图，图2中ZT1、ZT2、ZT3、YS1、YS2、YS3、TZ1、TZ2、TZ3、MDC1、MDC2、MDC3、MD1、MD2、MD3、MD4、ML5、ML6、ML7、ML8、GD1为各配变的配变名称；表1中，u1表示第1次采样的电压，u2表示第2次采样的电压，u3表示第3次采样的电压，u95表示第95次采样的电压，u96表示第96次采样的电压。

表1 一条线路下变压器日电压采样数据，单位：伏

步骤1.2：采用数学方法对数据进行处理，包括：缺失值填补、异常值替换、电压干扰装置检测、零采样值删除以及三相电压归一化等，将线路下的所有配变的预处理完的电压数据处理成电压矩阵：

式中，n表示每日的电压采样次数，m表示同一条线路下所辖配变数量，V

步骤2：进行皮尔逊相关系数计算以获取配变间的电压相关性系数矩阵以及配变和线路之间的电压相关性系数矩阵。将配变间的电压相关性系数矩阵、配变和线路之间的电压相关性系数矩阵以及电压矩阵进行处理后合并组成配电网特征矩阵；

步骤2.1：计算地理信息系统中保存的同一条线路下所有配变之间的皮尔逊相关系数组成配变间的电压相关性系数矩阵；

式中，ρ

式中，ρ

步骤2.2：计算所有配变和所属线路的皮尔逊相关系数，采用皮尔逊相关系数公式计算，组建配变和线路之间的电压相关性系数矩阵

步骤2.3：分别将电压矩阵、配变间的电压相关性系数矩阵以及配变和线路之间的电压相关性系数矩阵进行归一化；

式中，g

步骤2.4：将归一化后的电压矩阵、配变间的电压相关性系数矩阵以及配变和线路之间的电压相关性系数矩阵组合成原始配电网特征矩阵；

其中，X为归一化后的原始特征矩阵，

式中，T表示转置，x

步骤3：使用核主成分分析法（KPCA）对归集后的配电网特征矩阵降维，得到降维后的配电网特征矩阵，降维后的配电网特征矩阵最大程度保留了原始的特征。

步骤3.1：针对归集后的配电网特征矩阵

式中，a为核参数，x

式中，

步骤3.2：计算中心化核矩阵的特征值和特征向量；

步骤3.3：将中心化核矩阵的特征值，并进行降序排序，取中心化核矩阵的前r个特征值

步骤3.4：利用核矩阵和中心化核矩阵的前r个的特征值和对应的特征向量计算降维结果Z

步骤3.5：对降维结果计算主成分累计贡献率C；

当前P个降维结果的主成分累计贡献率C大于85%时，选取前P个降维结果的特征值及特征向量组建成为降维后的配电网特征矩阵。

步骤4：使用AO算法优化PNN网络，使PNN网络的平滑因子的选取达到最优值，更适应本发明中的线变关系精准辨识的实际需要。

步骤4.1：输入层，使用降维后的配电网特征矩阵作为PNN网络输入层的输入矩阵；

步骤4.2：模式层，计算输入矩阵和训练集之间的关系输出为：

式中，E表示输入矩阵，也是待识别样本，W

步骤4.3：求和层，将配变为正确或者错误的概率进行累加；

步骤4.4：输出层，输出每个神经元对应的配变的正确或者错误结论，表示为：

式中，

步骤4.5：针对PNN网络的平滑因子

式中，A

步骤4.6：初始化天鹰种群的位置，天鹰种群的初始位置为：

其中，

式中，L

步骤4.7：更新天鹰位置，采用天鹰算法的四种方法进行计算以实现位置的更新，分别是：①垂直俯冲攻击；②盘旋攻击；③低空飞行；④陆地行走捕食；

步骤4.8：再次使用适应度函数计算，判断是否满足适应度需求或者达到最大迭代次数，若不满足则返回步骤4.6，若满足则结束；

步骤4.9：输出满足条件的天鹰位置，通过解码以后得到最优平滑因子。值得注意的是，在这步中，天鹰的位置会随着AO算法的运行而不断进行迭代更新，直到满足所需条件。

步骤4.10：将使用AO算法得到的最优平滑因子送入到PNN网络进行计算。

步骤5：将实时监测的电压矩阵按步骤2和步骤3处理得到的降维后的配电网特征矩阵，送入到使用AO算法优化后的PNN网络进行分析以辨识发生线变关系错误的配变，得到表2所示的PNN输出结果。

表2 PNN输出结果

在本实施例中，在地理信息系统中所有的配变均是使用的MLX线进行供电的，经过使用本发明方法进行分析识别以后，发现PNN网络的输出结果中GD、ML系列、MD系列的配变和MLX线下的配变并不属于同一集群，判断其已经出现了线变关系错误，后经过实际的现场确认发现两条线路连接处的断路器已经断开，GD、ML系列、MD系列的配变实际上是有35kVGC变电站进行供电，与地理信息系统的记录不符造成了线变关系的错误，最终通过本发明所提方法成功辨识出了所有的错误配变。

对于线变关系的辨识事实上是一种分类问题，分类的结果只有正确和错误，若使用原始特征矩阵需要巨大的计算量，并且原始特征矩阵中的数值并不在同一量纲，经过数据的归一化及降维以后才能使其可用，最后，由于人工选取的PNN网络平滑因子

尽管已描述了本申请的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。

显然，本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样，倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内，则本申请也意图包含这些改动和变型在内。