一种接地网腐蚀速率预测方法与流程

本发明涉及发电站技术领域，特别涉及一种接地网腐蚀速率预测方法。

背景技术：

影响接地网腐蚀速率的因素非常复杂，其中主要的因素包括ph值、含水量、电导率、有机质、全氮、的浓度、cl^-的浓度、的浓度、ca²⁺的浓度、mg²⁺的浓度、k⁺的浓度、na⁺的浓度。接地网腐蚀速率与以上12种主要影响因素之间存在着非常复杂的非线性关系，且各影响因素之间相互关联、相互制约，很难通过某一个精确的函数清楚地表达出它们之间模糊的本构关系。

目前虽然也有一些预测方法，例如基于k均值法的rbf网络的腐蚀速率预测方法，但是这类方法的预测准确率普遍不高，效果不够好。

技术实现要素：

本发明实施例提供了一种接地网腐蚀速率预测方法，用以解决现有技术中存在的问题。

一种接地网腐蚀速率预测方法，包括：

对采集得到的样本进行归一化处理，获得归一化后的样本；

对归一化后的样本进行聚类，输出最优的聚类解集，具体包括：

初始化蚁群参数；

根据归一化后的样本计算各个聚类中心，以每个样本到聚类中心的距离之和最小为目标函数，计算每个样本的目标函数值；

将计算得到的每个样本的目标函数值从小到大排序，在排序后的目标函数值中，取前l个目标函数值对应的样本作为需要交换的样本，对需要交换的样本进行局部搜索；

进行局部搜索后，对前l个目标函数值对应的样本进行信息素更新；

按照以上步骤进行循环迭代，直到达到最大迭代次数num，输出最优聚类解集；

基于蚁群聚类的rbf网络学习，包括：

样本训练，以所述最优聚类解集的蚁群聚类中心值为rbf网络的隐层中心值取，以最优聚类解集的聚类个数为隐层节点个数，计算隐节点数据中心间的距离，以隐节点间的最小距离作为扩展常数，先计算各归一化后的样本离各数据中心的距离，然后计算隐节点的输出阵，最后得到输出权值和偏移量；

样本测试，根据样本训练中得到的参数对归一化后的样本进行测试，对测试结果反归一化，得到腐蚀速率。

优选地，对采集得到的样本按照以下公式进行归一化处理：

其中，αi表示归一化后的样本，i＝(1,2,…,n)，x表示采集得到的样本集，其中样本集中有n个样本xi和m个聚类解集{sj,j＝1,2,…,m}。

优选地，按照以下公式计算每个样本的目标函数值：

其中，xip为第i个样本xi的第p个属性，cjp为第j个聚类解集sj的第p个属性，若第i个样本xi属于第j个聚类解集sj，则wij＝1，否则wij＝0。

优选地，对需要交换的样本进行局部搜索，具体包括：

对于每一个样本，产生一个随机数q'，若随机数q'小于q，则选择聚类中心与样本间距离最短的聚类解集作为样本需要归属的类；如果该随机数q'大于或等于q，根据设定的转换概率pij选择样本需要转换的聚类解集：

其中，τij为第i个样本xi和所属第j个聚类解集sj间的标准化信息素；

确定每个样本归属的聚类解集后重新聚类，重新计算每个样本的目标函数值，如果重新聚类后的样本的目标函数值小于重新聚类前的样本的目标函数值，则保留重新聚类后的聚类解集，否则保留重新聚类前的聚类解集。

优选地，按照以下公式对前l个目标函数值对应的样本进行信息素更新：

τij(t+1)＝(1-r)τij(t)+δτij(t)

其中，r为信息素衰减系数，τij(t)为第i个样本xi与第j个聚类解集sj在t时刻的信息素浓度；δτij表示本次循环信息素的增量；表示第k只蚂蚁在本次循环中所释放的信息素，设lk为蚂蚁的目标函数值，q为一参数常量值，若蚂蚁k中的第i个样本xi属于第j个聚类解集sj，则否则

本发明提供的一种接地网腐蚀速率预测方法，用蚁群聚类算法确定rbf网络的隐层中心值，对基本的蚁群算法进行了进一步改进，并且加入了局部搜索，然后把改进的蚁群聚类算法与rbf神经网络结合用于接地网腐蚀速率的预测。大量数据计算表明，该方法预测接地网腐蚀速率具有较高的精度和较好的泛化能力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种接地网腐蚀速率预测方法的预测结果与实测结果的对比图；

图2为现有技术中基于k均值法的rbf网络的腐蚀速率预测结果与实测结果的对比图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提供了一种接地网腐蚀速率预测方法，该方法包括：

对采集得到的样本按照公式(1)进行归一化处理：

其中，αi表示归一化的样本，i＝(1,2,…,n)，x表示采集得到的样本集，其中样本集中有n个样本xi和m个聚类解集{sj,j＝1,2,…,m}；

对归一化的样本进行聚类，输出最优的聚类解集。具体步骤如下：

初始化蚁群参数，具体需要初始化的蚁群参数包括：蚂蚁数量ant、类的数量class、转换规则参数q(0<q<1)、信息素衰减系数r、局部搜索阈值h、迭代次数num、初始信息素tau和初始聚类解集solution；

根据归一化后的样本计算各个聚类中心，以每个样本到聚类中心的距离之和最小为目标函数，计算每个样本的目标函数值：

其中，xip为第i个样本xi的第p个属性，cjp为第j个聚类解集sj的第p个属性，若第i个样本xi属于第j个聚类解集sj，则wij＝1，否则wij＝0；

将计算得到的每个样本的目标函数值按照从小到大的顺序排序，在排序后的目标函数值中，取前l个目标函数值对应的样本作为需要交换的样本，对需要交换的样本进行局部搜索，具体为：

对于每一个样本，产生一个随机数q'，若随机数q'小于q，则选择聚类中心与样本间距离最短的聚类解集作为样本需要归属的类；如果该随机数q'大于或等于q，根据设定的转换概率pij选择样本需要转换的聚类解集：

其中，τij为第i个样本xi和所属第j个聚类解集sj间的标准化信息素；

确定每个样本归属的聚类解集后重新聚类，根据式(2)重新计算每个样本的目标函数值，如果重新聚类后的样本的目标函数值小于重新聚类前的样本的目标函数值，则保留重新聚类后的聚类解集，否则保留重新聚类前的聚类解集。

进行局部搜索后，对前l个目标函数值对应的样本按照下式进行信息素更新：

τij(t+1)＝(1-r)τij(t)+δτij(t)(5)

其中，τij(t)为第i个样本xi与第j个聚类解集sj在t时刻的信息素浓度；δτij表示本次循环信息素的增量；表示第k只蚂蚁在本次循环中所释放的信息素，设lk为蚂蚁的目标函数值，q为一参数常量值，若蚂蚁k中的第i个样本xi属于第j个聚类解集sj，则否则

按照以上步骤进行循环迭代，直到达到最大迭代次数num，输出最优聚类解集。

基于蚁群聚类的rbf网络学习，主要包括：

样本训练。隐层中心值和隐层节点个数的确定，这里rbf网络的隐层中心值取上述最优聚类解集的蚁群聚类中心值，隐层节点个数取最优聚类解集的聚类个数。计算隐节点数据中心间的距离，以隐节点间的最小距离作为扩展常数。先计算各归一化后的样本离各数据中心的距离，然后计算隐节点的输出阵，最后得到输出权值和偏移量。

样本测试。根据样本训练中得到的隐层中心值、隐层节点个数、扩展常数、输出权值和偏移量等参数对归一化后的样本进行测试，对测试结果反归一化，获得腐蚀速率。

为了检验该模型的正确性，我们用全国不同地方的变电站腐蚀速率及影响因素数据来验证。共有101组数据，如表1所示。

表1接地网腐蚀速率及影响因素

这些数据中，前93组用来训练，后8组用来测试。在蚁群聚类中参数设置如下：蚂蚁数目ant＝120，类的个数class＝20，迭代次数num＝2000，转换规则参数q＝0.9，信息素衰减系数r＝0.05，局部搜索阈值h＝0.03，l＝2。采用本发明的方法预测结果如图1所示。图2是用k均值与rbf网络结合的预测结果。

后8组数据的实测值为y，用本发明方法预测的结果用y1表示，基于k均值法的rbf网络的腐蚀速率预测结果用y2表示，则y、y1、y2如表2所示。

表2腐蚀速率的实测值及两种不同方法的预测值

把y、y1、y2分别看做一个向量，则y1到y的欧氏距离为||y1-y||2＝0.3438，y2到y的欧氏距离为||y2-y||2＝0.3973。很明显，||y1-y2＜||y2-y||2。因此，本发明的预测效果比基于k均值法的rbf网络的预测结果好，准确率高。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

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

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