本发明属于电力工程领域，涉及复杂环境下多因素电力系统的电力负荷分配方法；具体来说是一种基于信息交流策略连续域蚁群算法的电力系统经济调度方法，可用于火电经济调度计划、水火电协调、梯级调度和燃料约束调度等领域。
背景技术：
电力系统经济调度是能源生产和分配的核心。随着能源需求和燃料成本的增加，找到降低能源生产成本的方法变得越来越重要。电力系统经济调度需要找出电力系统中发电机的优化负荷，以满足整个电力需求，同时产生最小可能的发电成本。因此，需要制定和应用能够考虑发电机组特性成本函数曲线的数学优化模型，以找到最优的负荷分配。
电力系统的成本函数可以简单地用单调递增的二阶非线性函数来近似。然而，这样一个简单的成本函数没有适当地模拟出发电机的实际情况，首先电力系统涉及到阀点效应和多种燃料使用等约束。如果将阀点效应引入成本函数，会导致多个局部最小值产生。通常，在进行发电机负荷分配时，有一些区域的值是禁止运行的，禁止运行的区域使成本函数不连续。除此之外，还应考虑诸如爬坡率限制和输电线路损耗等约束条件，以便近似估计发电机的实际行为。如今，由于人们对环境的关注日益增加，燃烧燃料排放的有毒气体也被纳入成本函数。这些限制会使电力系统经济调度的任务更加复杂。因此，电力系统经济调度中的优化问题由多个目标组成。其任务是找到发电机组的负荷分配，使发电量满足整个系统电力需求，同时产生最小可能的发电成本。许多经典的优化技术，如线性规划、齐次线性规划、梯度搜索算法、lambda迭代法、模式搜索法、二次规划方法、拉格朗日乘子法和非线性规划法被用来解决电力系统经济调度问题。这些经典优化技术的主要问题在于它们返回的最终解决方案高度依赖于初始值，通常容易陷入局部最优值。成本函数的不连续性也严重影响了这些优化技术的性能。虽然动态规划技术可以处理不连续的非线性成本函数，但优化问题的维数对它来说是一个非常重要的影响因素。
技术实现要素：
针对上述问题，本发明提出一种基于信息交流策略连续域蚁群算法的电力系统经济调度方法，该方法以系统总运行费用为目标函数建立模型，同时考虑阀点负载效应、输电线路损失、爬坡率限制、禁止区域等限制。本发明方法不仅能够降低电能的费用，而且还能保障电力系统的电压质量，满足实际工程的设计需求。
本发明采用的技术方案：
一种基于信息交流策略连续域蚁群算法的电力系统经济调度方法，包括以下步骤：
(1)建立电力系统经济调度模型：
其中，n为系统发电机组的总个数，f为发电总成本，fi为机组i的发电成本，ai，bi，ci是机组i的成本系数，pi为机组i的有功输出功率，vi为阀点效应带来的变化值，di，ei为常数，
电力系统经济调度模型的限制条件为：
(1.1)传输损耗：
其中，bi，j为耗量系数矩阵b中第i行第j列的系数，ti为耗量系数矩阵t中第i行的系数，l为燃料耗量系数；
(1.2)运行区域：
其中，
(1.3)发电机爬升率：
其中，dri为机组i的向下爬坡速率值，uri为机组i的向上爬坡速率值，
(1.4)禁用区域：
其中，
(2)初始化各项参数：解档案个数k，关联度q＝0.0001，权值参数ξ＝0.85，蚂蚁个数s，优解率r，最大迭代次数iter_max，随机生成k个解；
(3)利用电力系统经济调度模型计算解档案中解的适应度值，并根据适应度值将解进行排序；
(4)采用标准蚁群算法更新解档案中的解，具体步骤为：
(4.1)计算解档案中每个解的权重，wi为第i个解在解档案中的权重，定义为：
其中，q是算法中的一个参数，被用来平衡局部最优和全局最优；
(4.2)计算选取每个解作为指导解的概率，选取第i个解的概率为：
(4.3)每只蚂蚁根据解档案中解的权重选取一个解xi，i＝1，2，…，k，并分别通过高斯函数gi，j(x)对这个解的d个解组件xi，j，j＝1，2，…d，进行取样；每一个高斯函数的均值ui，j＝xi，j，标准差σi，j为
(4.4)s只蚂蚁都产生新解之后，把这k+s个解重新排序，然后选取k个最优的解存储到解档案中。
(5)从更新后的解档案中选取m个优质解，m个优质解根据优解率r计算得出，则有m＝ceil(r×k)，ceil表示向上取整；
(6)在m个优质解中，每一个解都选择另外一个解来进行信息交换，得到候选解，采用贪婪的方式接受能够改善解的质量的候选解，具体步骤为：
(6.1)在m个优质解中，信息交流策略通过如下公式为每个解xi，i＝1，2，...m确定另一个解xh，h∈(1，m)进行维变量互换。
其中，d表示解的维数；
(6.2)每个xi解确定好进行维变量互换的解之后，则产生的新解
其中，rand表示产生一个(0,1)之间的随机数，pj为：
(6.3)通过贪婪的方式确定xi和
(7)当前迭代次数加1，判断是否到达最大迭代次数，当前迭代次数超过最大迭代次数iter_max时，则停止迭代，输出解档案中第一个解作为全局最优值；当前迭代次数未超过最大迭代次数iter_max时，则跳转到步骤(3)；
(8)根据得到的解进行电力负荷分配，得到最优电力系统经济调度。
本发明方法是一种基于信息交流策略连续域蚁群算法的电力系统经济调度方法，该发明方法在acor算法中引入信息交流策略，避免出现指导解的选择过于单一的情况。信息交流策略在解档案中选取一部分较优的解，在这部分较优的解里每一个解都选择另外一个解来进行信息交换，这样可以增加解的多样性，避免陷入局部最优。同时，信息交流策略通过对解的信息交换，融合不同解中较优的维变量，得到较优的子解，而且只处理了部分较优解,这样能加快算法的收敛速度，提高算法的收敛精度，使得该发明方法能得到较优的电力系统经济调度方法。
附图说明
图1为本发明方法的流程图；
具体实施方式
以下结合实例对本发明作进一步详细说明，但不作为对本发明的限定。
实施例1
步骤1.建立电力系统经济调度模型。以电力系统相关的总燃料成本为目标函数，发电机的燃料成本函数可以被认为是由发电机产生功率的二阶二次函数，但是，考虑实际情况，热发电机通常使用蒸汽涡轮机，其由于蒸汽阀的连续打开而出现阀点效应。当考虑这些影响时，发电机的成本函数不再是线性的。考虑阀点效应后，其耗量特性函数为：
式中，n为系统发电机组的总个数，f为发电总成本，fi为机组i的发电成本，ai，bi，ci是机组i的成本系数，pi为机组i的有功输出功率，vi为阀点效应带来的变化值，di，ei为常数，
在解决多发电机系统的经济调度问题时，还需要考虑一些系统限制因素，这些因素将作为模型的限制条件：
1)传输损耗：总传输损耗是电力系统经济调度问题的一个重要制约因素。不仅希望在系统中发生的损失与总燃料成本一起降至最低，而且系统还必须产生足够的电力以满足负载需求并补偿传输损耗。系统中的总损失可以使用以下关系计算：
其中bi，j为耗量系数矩阵b中第i行第j列的系数，ti为耗量系数矩阵t中第i行的系数，l为燃料耗量系数；
系统运行时要满足平衡约束，即整个发电机系统产生的功率必须等于总功率需求和系统内部以及由于传输造成的损失之和，即：
其中pd为总功率需求，pl为系统内部以及由于传输造成的损失。
2)运行区域：每台热电发电机都有一个指定的范围，其运行稳定。因此，希望发电机在该范围内运行以保持系统稳定性。
其中
3)发电机爬升率：对于负载的任何突然变化，发电机会增加或减少其供电量，以便将功率失配降至零。但是，发电机组只能按照上升和下降速率确定的一定速率改变电源。如果发电机在特定点运行，则其运行点可以改变，直至由斜坡率确定的某个极限。因此，为了更接近实际情况，在处理电力系统经济调度问题时，也需要考虑加速率限制。
式中，dri为机组i的向下爬坡速率值，uri为机组i的向上爬坡速率值
4)禁用区域：可能存在指定范围内的操作区域，系统可能会失去稳定性。这些被称为禁止区和实际操作；这些地区通常在发电过程中被避免。
其中，
步骤2.初始化各项参数：解档案个数k，关联度q＝0.0001，权值参数ξ＝0.85，蚂蚁个数s，优解率r，最大迭代次数iter_max，随机生成k个解；
步骤3.利用电力系统经济调度模型计算解档案中解的适应度值，并根据适应度的值将解进行排序；
步骤4.采用标准蚁群算法更新解档案中的解，具体步骤为：
(4.1)计算解档案中每个解的权重，wi为第i个解在解档案中的权重，定义为：
其中q是算法中的一个参数，被用来平衡局部最优和全局最优。
(4.2)计算选取每个解作为指导解的概率，选取第i个解的概率为：
(4.3)每只蚂蚁根据解档案中解的权重选取一个解xi＝(i＝1，2，…，k)，分别通过一个高斯函数gi，j(x)对这个解的d个解组件xi，j(j＝1，2，…d)进行取样。每一个高斯函数的均值ui，j＝xi，j，标准差σi，j为
(4.4)s只蚂蚁都产生新解之后，把这k+s个解重新排序然后选取k个最优的解存储到解档案中。
步骤5.根据参数r计算出m个较优解，具体步骤为：
(5.1)根据适应度值的质量来评价解的质量，只选取一部分质量较优的解进行处理，设这部分解的个数为m，则有m＝ceil(r×k)。
步骤6.在m个较优的解里每一个解都选择另外一个解来进行信息交换得到候选解，采用贪婪的方式接受能够改善解的质量的候选解，具体步骤为：
(6.1)在m个较优解中，信息交流策略通过如下公式为每个解xi，i＝1，2，...m确定另一个解xh，h∈(1，m)进行维变量互换。
(6.2)每个xi解确定好进行维变量互换的解之后，则产生的新解
其中，rand表示产生一个(0,1)之间的随机数，pj为：
(6.3)通过贪婪的方式确定xi和
步骤7.当前迭代次数加1，判断是否到达最大迭代次数，如果当前迭代次数超过最大迭代次数iter_max，则停止迭代，输出解档案中第一个解作为全局最优值，如果当前迭代次数未超过最大迭代次数iter_max，则跳转到步骤(3)；
步骤8.根据得到的解进行电力负荷分配，得到最优电力系统经济调度。
本发明的效果可以通过以下仿真实验进一步说明。
1.仿真条件及参数
本发明方法和基于acor的方法、基于daco的方法、基于haco的方法、基于abc的方法、基于pso的方法做实验对比。acor表是标准的连续域蚁群算法，daco算法和haco算法是目前较为常用的连续域蚁群算法的改进算法，abc表示蜂群算法，pso表示粒子群算法。采用ieee-30节点标准测试系统。该系统为10机组电力系统，系统总负荷为300mw。具体参数如表1所示：
表1系统燃料耗量的参数和有功功率上下限
燃料耗量系数矩阵如下：
t＝0.001*[0.2870.0120.08960.14710.00870.31210.2330.11230.09120.1121]
l＝0.038
实验中基于acor的方法参数设置为：蚂蚁个数m＝2，解档案个数k＝50，q＝0.0001，ξ＝0.85，基于haco的方法中学习速率α＝0.9。基于daco的方法参数设置为：m＝2，k＝50，q＝0.9，ξ＝0.85。本发明方法参数设置为：m＝50，k＝50，q＝0.0001，ξ＝0.85，r＝0.2。为了消除因初始值不同而导致结果产生差异，所有算法对每个测试函数运行之前都使用同一组随机生成的初始值。同时，所有算法都独立运行15次，每次运行的最大迭代次数为30000次。
2.仿真内容及结果分析
实验结果如表2所示：
表2电力系统发电机功率输出和总发电成本(300mw)
表2的优化结果表明，在计算电力系统经济调度问题上，本发明方法得到的最优调度方案为：系统总功率300.94mw，网损0.94mw，总费用5988.58($/h)，基于acor的方法最优调度方案为：系统总功率303.99mw，网损3.99mw，总费用6326.44($/h)，基于daco的方法最优调度方案为：系统总功率306.46mw，网损6.46mw，总费用6645.12($/h)，基于haco的方法最优调度方案为：系统总功率303.17mw，网损3.17mw，总费用6254.75($/h)，，基于abc的方法最优调度方案为：系统总功率301.26mw，网损1.26mw，总费用6114.75($/h)，基于pso的方法最优调度方案为：系统总功率308.19mw，网损8.19mw，总费用6763.32($/h)，本发明方法跟其他方法相比，既降低了总功率又降低了总费用，说明改进的蚁群算法在计算电力系统经济调度问题上更有效。
为了测试本发明方法的鲁棒性，将本发明方法和其他方法应用到系统负载为500mw的电力系统中，实验参数同上，实验结果如表3所示。
表3电力系统发电机功率输出和总发电成本(500mw)
从表3可以看出，在500mw负载，本发明方法跟其他方法相比，既降低了总功率又降低了总费用，说明改进的蚁群算法在计算电力系统经济调度问题上更有效。