一种基于分解文化进化算法的多目标生态调度方法和系统与流程
本发明属于水电能源优化领域，更具体地，涉及一种基于分解文化进化算法的多目标生态调度方法和系统。
背景技术：
传统的多目标调度侧重于以发电量、生态调度或者最小出力为单一主要目标，其他次要目标作为约束条件，缺乏将三者均列为主要目标的考虑和研究，实际应用较少。传统调度忽略三者间的内在联系，导致弃水现象频发，降低经济效益，增加运营成本。
现有多目标粒子群算法有基于分解的多目标进化算法(moea/d)和文化算法(ca)，这两种算法均存在许多问题。例如在一种基于分解的多目标进化算法(moea/d)中，解决方案的比较完全基于惩罚函数的边界交叉点(pbi)这样的分解方法，这种单一的评估方式在进化过程中可能会损失部分精英个体。文化算法(ca)是一种基本框架，在进化过程中提取信息，经处理的信息反馈促进进化过程，但是单独使用时有较大局限性。
由此可见，现有技术存在评估方式单一，局限性较大，经济效益低的技术问题。
技术实现要素：
针对现有技术的以上缺陷或改进需求，本发明提供了一种基于分解文化进化算法的多目标生态调度方法和系统，由此解决现有技术存在评估方式单一，局限性较大，经济效益低的技术问题。
为实现上述目的，按照本发明的一个方面，提供了一种基于分解文化进化算法的多目标生态调度方法，依次包括如下步骤：
(1)以水位为自变量，以发电量最大、生态流量改变度最小和最小出力最大化为目标建立多目标函数模型；
(2)为多目标函数模型生成多组权重向量，基于多组权重向量初始化邻域索引集和归档集，随机生成包含多个个体的初始种群，所述个体为水库调度方案；
(3)对于每个个体，从该个体的邻域索引集中随机选取个体作为父代个体进行交叉变异得到子代个体，将父代个体和子代个体进行比较并保留较优个体；
(4)当初始种群中每个个体都完成步骤(3)后则初始种群进化完成，从初始种群中选取较优个体作为较优水库调度方案保存至归档集；
(5)当进化次数满足预设次数，得到最终的归档集，最终的归档集是以发电量最大、生态流量改变度最小和最小出力最大化为目标的最优水库调度方案，否则进入步骤(2)。
进一步地，初始化邻域索引集的具体实现方式为：
对于初始种群中每个个体，选取与其距离最近的部分个体，构成邻域索引集。
进一步地，初始化的归档集中较优水库调度方案的个数为0。
进一步地，步骤(3)包括如下子步骤：
(31)对于每个个体，从0-1中产生一个随机数与进化参数进行比较，大于进化参数则进行步骤(32)，否则舍弃该个体；
(32)对于大于进化参数的个体，从邻域索引集中随机选取个体作为父代个体进行交叉变异，得到子代个体；
(33)将父代个体和子代个体进行比较，保留较优个体。
进一步地，步骤(33)包括如下子步骤：
(331)判断父代个体和子代个体的支配关系，若父代个体支配子代个体，则将父代个体保留，若子代个体支配父代个体，则将子代个体保留；若父代个体和子代个体之间不具有支配关系，则进入步骤(332)；
(332)将父代个体及其对应的权重向量带入多目标函数模型，得到父代个体的多目标函数值，将子代个体及其对应的权重向量带入多目标函数模型，得到子代个体的多目标函数值，保留多目标函数值较大者。
进一步地，支配关系的判断方式为：
当父代个体对应的发电量、生态流量改变度和最小出力均大于子代个体对应的发电量、生态流量改变度和最小出力，则父代个体支配子代个体；
当父代个体对应的发电量、生态流量改变度和最小出力均小于子代个体对应的发电量、生态流量改变度和最小出力，则子代个体支配父代个体；
否则父代个体和子代个体之间不具有支配关系。
进一步地，步骤(4)还包括：当归档集中较优水库调度方案的个数大于m个，对归档集进行裁剪，使得归档集中的较优水库调度方案的个数为m个。
进一步地，对归档集进行裁剪的具体实现方式为：
判断第m+1个较优水库调度方案与归档集中前m个较优水库调度方案之间是否有支配关系，若有则根据支配关系对第m+1个较优水库调度方案进行取舍；
否则，计算归档集中m+1个较优水库调度方案两两之间的距离，计算距离最小值对应的两个较优水库调度方案的多目标函数值，保留多目标函数值较大者。
按照本发明的另一方面，提供了一种基于分解文化进化算法的多目标生态调度系统，包括如下模块：
模型建立模块，用于以水位为自变量，以发电量最大、生态流量改变度最小和最小出力最大化为目标建立多目标函数模型；
初始化模块，用于为多目标函数模型生成多组权重向量，基于多组权重向量初始化邻域索引集和归档集，随机生成包含多个个体的初始种群，所述个体为水库调度方案；
交叉变异模块，用于对每个个体，从该个体的邻域索引集中随机选取个体作为父代个体进行交叉变异得到子代个体，将父代个体和子代个体进行比较并保留较优个体；
选取模块，用于当初始种群中每个个体都完成交叉变异模块后则初始种群进化完成，从初始种群中选取较优个体作为较优水库调度方案保存至归档集；
最优解模块，用于当进化次数满足预设次数，得到最终的归档集，最终的归档集是以发电量最大、生态流量改变度最小和最小出力最大化为目标的最优水库调度方案，否则执行初始化模块。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
(1)本发明以水位为自变量，以发电量最大、生态流量改变度最小和最小出力最大化为目标建立多目标函数模型，对多个目标进行评估，然后通过生成权重向量、初始化邻域索引集、种群进化得到最终的归档集，进而最终得到以发电量最大、生态流量改变度最小和最小出力最大化为目标的经济可靠的水电站优化调度方案，对解决分布不均引起的重大生态问题和选择实际水库运行控制方案具有重要指导意义。
(2)与传统的moea/d算法和ca算法不同，本发明提出的分解文化算法由于首先进行支配关系判断，然后进行多目标函数值比较，评估方式多样，不易陷入局部最优，具有较好的收敛性，可得到非支配解集，适用于多目标模型的求解，且有利于开展水电站调度决策和发电可靠性评估等工作。
附图说明
图1是本发明实施例提供的一种基于分解文化进化算法的多目标生态调度方法的流程图；
图2是本发明实施例1提供的三峡水库适宜生态流量图；
图3是本发明实施例1提供的在方案9、方案15、方案42状态下，三峡水库水位变化过程；
图4是本发明实施例1提供的在方案9、方案15、方案42状态下，三峡水库下泄流量变化过程。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示，一种基于分解文化进化算法的多目标生态调度方法，依次包括如下步骤：
(1)以水位为自变量，以发电量最大、生态流量改变度最小和最小出力最大化为目标建立多目标函数模型，并设置多目标函数模型的模型约束。
根据水电站发电系数、水电站每个时段的净水头以及水电站每个时段的发电流量得到水电站每个时段的出力，利用水电站每个时段的出力建立以水电站发电量最大为目标的第一函数：
其中，f1是水电站发电量，t代表时段总数，t是时段序列号，t∈[1，t]，δt是时间间隔，nt是水电站第t个时段的出力，a为水电站发电系数，ht为水电站第t个时段的净水头，qt为水电站第t个时段的发电流量。
净水头是通过本时段水位加上下一时段水位的平均值减去下游水位和水头损失得到的，下游水位的计算为：依据本时段水位和下一时段水位，可以算出本时段下泄流量，通过多个时段的下泄流量绘制尾水流量曲线，通过尾水流量曲线可以得到下游水位。
为了达到最佳的生态效益，生态供水量应至少满足生态基流量的需求，即两者之间的差距应最小化。因此，以水电站的生态供水量与水电站的生态基础流量之间的差距最小为目标建立第二函数：
其中，f2为水电站的生态供水量与水电站的生态基础流量之间的差距即为生态流量改变度，qt为水电站第t个时段的生态供水量(即为下泄流量)，q′t为水电站第t个时段的生态基础流量。
以水电站每个时段的下泄流量最小且水电站整个时段的泄流量最大为目标建立第三函数：
其中，rt为水电站在第t个时段的泄流量(即为出力)，f3为水电站在t时段的泄流量。
第一函数、第二函数和第三函数组合形成多目标函数模型：
f＝maxf1-minf2+maxf3
设置多目标函数模型的模型约束。模型约束包括：水量平衡约束、水位约束、库容约束、下泄流量约束和出力约束。
水量平衡约束为：
vt＝vt-1+it-rt，t∈[1，t]
其中，vt是水电站第t个时段的库容，vt-1是水电站第t-1个时段的库容，it是水电站第t个时段的入库流量，rt为水电站在第t个时段的泄流量。
水位约束为：
其中，zt、
库容约束为：
其中，
下泄流量约束为：
其中，
出力约束为：
其中，
(2)为多目标函数模型生成多组权重向量，基于多组权重向量初始化邻域索引集和归档集，随机生成包含多个个体的初始种群，所述个体为水库调度方案；对于初始种群中每个个体，选取与其距离最近的部分个体，构成邻域索引集。初始化的归档集中较优水库调度方案的个数为0。
(3)对于每个个体，从该个体的邻域索引集中随机选取个体作为父代个体进行交叉变异得到子代个体，将父代个体和子代个体进行比较并保留较优个体；
步骤(3)包括如下子步骤：
(31)对于每个个体，从0-1中产生一个随机数与进化参数进行比较，大于进化参数则进行步骤(32)，否则舍弃该个体；
(32)对于大于进化参数的个体，从邻域索引集中随机选取个体作为父代个体进行交叉变异，得到子代个体；
(33)将父代个体和子代个体进行比较，保留较优个体。
具体地，进化参数为de算法中的进化参数。
进一步地，步骤(33)包括如下子步骤：
(331)判断父代个体和子代个体的支配关系，若父代个体支配子代个体，则将父代个体保留，若子代个体支配父代个体，则将子代个体保留；若父代个体和子代个体之间不具有支配关系，则进入步骤(332)；
(332)将父代个体及其对应的权重向量带入多目标函数模型，得到父代个体的多目标函数值，将子代个体及其对应的权重向量带入多目标函数模型，得到子代个体的多目标函数值，保留多目标函数值较大者。
进一步地，支配关系的判断方式为：
当父代个体对应的发电量、生态流量改变度和最小出力均大于子代个体对应的发电量、生态流量改变度和最小出力，则父代个体支配子代个体；
当父代个体对应的发电量、生态流量改变度和最小出力均小于子代个体对应的发电量、生态流量改变度和最小出力，则子代个体支配父代个体；
否则父代个体和子代个体之间不具有支配关系。
(4)当初始种群中每个个体都完成步骤(3)后则初始种群进化完成，从初始种群中选取较优个体作为较优水库调度方案保存至归档集；
步骤(4)还包括：当归档集中较优水库调度方案的个数大于m个，对归档集进行裁剪，使得归档集中的较优水库调度方案的个数为m个；
所述对归档集进行裁剪的具体实现方式为：
判断第m+1个较优水库调度方案与归档集中前m个较优水库调度方案之间是否有支配关系，若有则根据支配关系对第m+1个较优水库调度方案进行取舍；
否则，计算归档集中m+1个较优水库调度方案两两之间的距离，计算距离最小值对应的两个较优水库调度方案的多目标函数值，保留多目标函数值较大者。
(5)当进化次数满足预设次数，得到最终的归档集，最终的归档集是以发电量最大、生态流量改变度最小和最小出力最大化为目标的最优水库调度方案，否则进入步骤(2)。
实施例1
本发明实施例1以2003年三峡大坝的实际入流量作为调度实例进行计算。三峡工程是当今世界最大的水利发电工程，蓄水高度175米，坝高185米。装机容量2240万千瓦，最大排水量10万立方米/秒，总库容393亿立方米。
调度周期和计算周期划分如下：
根据初步设计，三峡水库调度采用“蓄水排水”原则。在6月至9月的汛期，水库将以145.0m的洪水位运行，蓄水则从10月1日开始。本发明时间长度为一年。每个计算周期以10天为步长，每个月分为三个计算周期，一年计算周期总数为36。本发明采用tennant方法计算三峡水库的生态基流。处理结果如图2所示。
分解文化(mocea/d)算法参数设置如下：进化参数包括：交叉概率和变异概率，根据mocea/d算法的优化思想和过程，程序由java编写，目的是优化设计。种群大小选择为42，最大进化次数为5000，交叉概率为0.8，变异概率为0.03。通过mocea/d算法迭代计算获得42组pareto最优解。42组pareto前沿与最优解集分布如表1所示。
本发明利用2003年三峡大坝的实际入流量作为最优调度实例进行计算。在满足发电效益，保证最小出力的前提下，兼顾考虑生态需水量。表1显示了三峡水库多目标优化调度的42种优化方案。从表1可以看出，生态流量的改变度与发电量成反比例。单方面提高生态调度效益将不可避免地导致发电效益降低。同时可看出，mocea/d算法适用于求解多目标优化模型问题。
表1三峡水库多目标优化调度方案
为进一步分析各种调度方案之间的差异，选取了三种方案进行比较：方案9(最大发电量)，方案42(生态流量改变度最小)和方案15(最小出力最大化)。(1)当多目标优化模型以发电量为主时，可选取方案9；(2)当多目标优化模型以生态流量改变度为主时，可选取方案42；(3)当多目标优化模型以最小出力为主时，可选取方案15；因此，可根据电站的实际情况选取对应方案。
图3和图4显示了三种典型优化方案下三峡水库的水位和泄流量变化过程。从图中可看出诸多差异，部分差异不甚明显。多目标优化模型需深入研究三个目标内在联系。为进一步分析这种联系，故选取三种最具代表性方案进行对比：方案9(最大发电量)，方案42(生态流量改变度最小)和方案15(最小出力最大化)。
三峡水库全年计算可通过数据分析比较三个目标之间的内在联系。(1)仅考虑发电效益。相比方案15，方案9的发电量增加0.31％；相比方案42，方案9的发电量增加0.31％。(2)仅考虑生态调度效益。相比方案9，方案42的生态流量改变度减小14.69％；相比方案15，方案42的生态流量改变度减小10.35％。(3)仅考虑最小出力最大化时，相比方案9，方案15的最小出力增加1.02％；相比方案42，方案15的最小出力增加0.93％。以上计算分析结果如表2所示。
表2三种典型方案的各类指标比较
上述数据充分说明略微减小电站电量或电站的功率输出的措施是合理的。例如，发电量减少约0.31％，生态流量改变度减少14.69％。其他请参考表格。此外，本发明实施例1选取三峡水库正常运行年度进行计算，结果对汛期前后对防洪效果无显著影响。既不会增加防洪压力，又能提高水资源的利用程度。本发明为实际多目标优化调度的决策提供了参考。
与传统的moea/d算法和ca算法不同，本发明提出的mocea/d算法不易陷入局部最优，具有较好的收敛性，可得到非支配解集，适用于多目标模型的求解，且有利于开展水电站调度决策和发电可靠性评估等工作。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
- 还没有人留言评论。精彩留言会获得点赞！