一种基于HMAQGA的水库群联合优化调度方法技术

本发明专利技术公开了一种基于HMAQGA的水库群联合优化调度方法，步骤如下：首先获取水库群的基本信息，对系统进行概化；其次分析水库群任务，构建水库群多目标联合优化调度数学模型；然后执行HMAQGA求解水库群多目标联合优化调度模型获取Pareto非劣解集；最后采用基于组合赋权的模糊优选法在Pareto解集中确定最佳水库群调度运行方案。本发明专利技术实现全局寻优，提高计算效率，满足自主选择水库群最佳调度方案要求。

【技术实现步骤摘要】
一种基于HMAQGA的水库群联合优化调度方法
本专利技术属于水利水电领域的水库
，特别是一种基于HMAQGA的水库群联合优化调度方法。
技术介绍
合理高效的水库调度模式，能在硬件设施不变的情况下，获得更可观的社会、经济效益。随着流域水库群规模的不断扩大，对于水库群联合运行方式的要求也不断提高。因此，如何通过模型的构建及求解，快速准确地找到大规模水库群的最优调度方案，为水库群联合调度提供理论基础，已成为亟待解决的问题。近年来，随着数学规划和计算机水平的日益提升，在传统智能算法基础上引入Pareto非劣解集理论的多目标进化算法，如多目标遗传算法、多目标粒子群算法、多目标蛙跳算法及多目标蚁群算法等被广泛应用于水库群联合优化调度中。但是部分智能优化算法存在寻优过程随机性较大，常出现多次求解结果不一致的现象，且易陷入局部搜索、易早熟。因此，针对复杂的水库群优化调度模型，如何选择合理高效的优化算法进行模型求解，值得进一步研究。上世纪量子力学的提出，使我们能从原子层面了解物质的基本结构，并于世纪末由Shor教授研究应用于信息科学，提出了量子算法。1996年，Narayanan结合了量子多宇宙的思想，将其引入遗传算法中，提出了一种全新的量子衍生算法，并在旅行商问题的应用中取得成功。2000年，K.H.Han等在遗传算法中引入量子的态矢量编码方式，开创了量子遗传算法，提出了量子旋转门更新的方式，实现了个体染色体的量子基因遗传，提高了量子遗传算法的可行性。多目标量子遗传算法因其计算高效、优化性能好、精度高等优点，很多学者加入到该算法的研究中来，使其成为解决多目标优化问题的有力工具。在水库调度领域中，量子遗传算法的运用较少，尤其在多目标问题中，量子遗传算法的应用方面的文献稀缺，现阶段提出的基于多目标量子遗传算法的水库群联合优化调度存在易陷于局部最优解、非劣解分布不均等缺陷。同时，采用多目标决策方法处理非劣解集时，需要基于主客观因素相结合确定评价指标权重。
技术实现思路
本专利技术的目的在于针对传统多目标量子遗传算法容易陷入局部收敛，外部档案集非劣解分布不均，个体难收敛到非劣解前沿等缺陷，提出一种基于调和距离多目标的自适应量子遗传算法(Harmonicdistancemulti-objectivequantum-inspiredAdaptiveGeneticAlgorithm，HMAQGA)的水库群联合优化调度方法，并在外部档案非劣解集的基础上采用基于组合权重的模糊优选法选择最佳水库群调度方案。技术方案：一种基于HMAQGA的水库群联合优化调度方法，包括以下步骤：Step1：获取水库群的基本信息数据，包括：各水库的基本特征参数、水位-库容关系曲线、下游水位-下泄流量关系曲线等，并根据各水库之间的水力、水文联系对水库群系统进行概化；Step2：分析水库群任务，构建以发电量最大、缺水量最小、弃水量最小、最小出力最大等的目标函数，以水量平衡、水库下泄流量、机组过流能力、下游防洪流量、水库水位为约束的水库群多目标联合优化调度模型；Step3：执行HMAQGA求解水库群多目标联合优化调度模型，获取Pareto(帕累托)非劣解集；Step4：采用基于组合赋权的模糊优选法在Pareto非劣解集中确定最佳水库群调度运行方案。进一步地，所述Step3：执行HMAQGA求解水库群多目标联合优化调度模型，获取Pareto非劣解集，具体为：Step3-1，选择以各时段各个水库的水位值作为决策变量，确定各个水库的水位值的上下限，划分水库群整体调度的周期时段T，确定种群规模N，最大迭代次数Ngen，决策变量个数Nd，量子比特位数Nq，交叉概率Pc，灾变因子gk，灾变个体数目Nk，外部归档集规模Np，适应度网格数目Ngrid；Step3-2，令迭代次数g＝0，对种群中的每一个个体进行初始化，其中，每个个体的概率幅对(αi,l,βi,l)初始值均为产生N个量子染色体，组成初始种群Q(g)；Step3-3，测量种群Q(g)中的所有个体得到相应的确定解，计算所有个体的非支配排序等级[Pig]rank(i＝1,2,...,N)、调和距离并对其进行基于调和距离的快速非支配排序，根据相应的原则，选取Np个个体作为初始的外部归档集Step3-4，基于自适应网格机制，从外部归档集中为每一个个体选取一个目标个体；Step3-5，基于量子旋转门U(θ)，根据量子调整策略，对种群Q(g)进行更新操作，并采用式(1)进行量子H门修正个体概率幅，得到新的种群Q(g+1)；Step3-5，采用量子旋转门U(θ)，根据表1量子调整策略，对种群Q(g)进行更新操作，并利用下式进行量子H门修正个体概率幅，得到新的种群Q(g+1)；式中，[α′iβ′i]为量子旋转门更新后基因，[α″iβ″i]为H门修正后的基因，Step3-6，将更新后种群Q(t+1)与外部归档集组合，将组合后的N+Np个个体实施测量，得到相应的确定解，计算所有个体的非支配排序等级调和距离并对其进行快速非支配排序，选取前Np个个体组成更新后的外部归档集Step3-7，产生随机概率值P(P∈[0,1])，当P≤Pi时进行量子交叉操作，利用二项式交叉使种群个体间的信息得到充分交流；当P＞Pi时则跳过量子交叉操作，进行下一步；Step3-8，每进行gk次迭代均实施一次量子灾变操作，选择调和距离最小的Nk个个体进行重生成操作，使Nk个个体初始化；Step3-9，判断迭代次数g是否已达到Ngen，若是，则终止运行并输出结果；若否，则g＝g+1，返回Step3-4继续进化。进一步地，所述Step3-3具体包括以下步骤：Step3-3-1对于每个个体i构造两个参数ni和Si，其中ni表示当前种群中支配个体i的个体数量，Si表示被个体i支配的全部个体集合；Step3-3-2检索所有ni＝0的个体，将他们存入子集F1，此即等级为1的非支配个体集合，对该集合内的个体赋予同一非支配等级P[i]rank，剩余个体集合为P1。接着检索P1中所有的nk＝1个体，将他们存入子集F2中，即第二级非支配个体集合，对该集合内的个体赋予一个相同的非支配等级P[i]rank+1。以此类推，直到种群中的所有个体都被分级为止，每一个个体均被赋予非支配等级；Step3-3-3按下式计算个体的调和距离di，以调和距离来替代原先快速非支配排序法中的拥挤距离，以更好地评价个体的拥挤度：式中，di,k表示个体i与其邻近的第k个个体的欧几里得距离，k表示总个体数减1；Step3-3-4非支配排序选取精英个体组成外部归档集。经上述操作，每个个体i均具有两个参数P[i]rank、di，即非支配等级与调和距离。此时，采用以下原则选取精英个体组成外部归档集：当个体间非支配等级不同时，选取等级较低的个体，以保证种群的收敛性；当个体间非支配等级相同时，则选取调和距离较大的个体，以保证种群的多样性。进本文档来自技高网...

【技术保护点】
1.一种基于HMAQGA的水库群联合优化调度方法，其特征在于，包括以下步骤：/nStep 1：获取水库群的基本信息数据，包括：各水库的基本特征参数、水位-库容关系曲线、下游水位-下泄流量关系曲线，并根据各水库之间的水力、水文联系对水库群系统进行概化；/nStep 2：分析水库群任务，构建以发电量最大、缺水量最小、弃水量最小、最小出力最大的目标函数，以水量平衡、水库下泄流量、机组过流能力、下游防洪流量、水库水位为约束的水库群多目标联合优化调度数学模型；/nStep 3：执行HMAQGA求解水库群多目标联合优化调度模型，获取Pareto非劣解集；/nStep 4：采用基于组合赋权的模糊优选法在Pareto非劣解集中确定最佳水库群调度运行方案。/n

【技术特征摘要】
1.一种基于HMAQGA的水库群联合优化调度方法，其特征在于，包括以下步骤：
Step1：获取水库群的基本信息数据，包括：各水库的基本特征参数、水位-库容关系曲线、下游水位-下泄流量关系曲线，并根据各水库之间的水力、水文联系对水库群系统进行概化；
Step2：分析水库群任务，构建以发电量最大、缺水量最小、弃水量最小、最小出力最大的目标函数，以水量平衡、水库下泄流量、机组过流能力、下游防洪流量、水库水位为约束的水库群多目标联合优化调度数学模型；
Step3：执行HMAQGA求解水库群多目标联合优化调度模型，获取Pareto非劣解集；
Step4：采用基于组合赋权的模糊优选法在Pareto非劣解集中确定最佳水库群调度运行方案。


2.如权利要求1所述的基于HMAQGA的水库群联合优化调度方法，其特征在于，Step3包括以下步骤：
Step3-1，选择以各时段各个水库的水位值作为决策变量，确定各个水库的水位值的上下限，划分水库群整体调度的周期时段T，确定种群规模N，最大迭代次数Ngen，决策变量个数Nd，量子比特位数Nq，交叉概率Pc，灾变因子gk，灾变个体数目Nk，外部归档集规模Np，适应度网格Ngrid；
Step3-2，令迭代次数g＝0，对种群中的每一个个体进行初始化，其中，每个个体的概率幅对(αi,l,βi,l)初始值均为产生N个量子染色体，组成初始种群Q(g)；
Step3-3，测量种群Q(g)中的所有个体计算所有个体的非支配排序等级调和距离并对其进行基于调和距离的快速非支配排序，根据相应的原则，选取Np个个体作为初始的外部归档集
Step3-4，基于自适应网格机制，从外部归档集中为每一个个体选取一个目标个体；
Step3-5，基于量子旋转门U(θ)，根据量子调整策略，对种群Q(g)进行更新操作，并采用式(1)进行量子H门修正个体概率幅，得到新的种群Q(g+1)；



式中，[α′iβ′i]为量子旋转门更新后基因，[α″iβ″i]为H门修正后的基因，
Step3-6，将更新后种群Q(g+1)与外部归档集组合，将组合后的N+Np个个体实施测量，计算所有个体的非支配排序等级调和距离并对其进行快速非支配排序，选取前Np个个体组成更新后的外部归档集
Step3-7，产生随机概率值P(P∈[0,1])，当P≤Pc时进行量子交叉操作，利用二项式交叉，使种群个体间的信息得到充分交流；当P＞Pc时则跳过量子交叉操作，进行下一步；
Step3-8，每进行gk次迭代均实施一次量子灾变操作，选择调和距离最小的Nk个个体进行重生成操作，使Nk个个体初始化；
Step3-9，判断迭代次数g是否已达到Ngen，若是，则终止运行并输出结果；若否，则g＝g+1，返回Step3-4继续进化。


3.如权利要求2所述的基于HMAQGA的水库群联合优化调度方法，其特征在于，所述Step3-3包括以下步骤：
Step3-3-1对于每个个体i构造两个参数ni和Si，其中ni表示当前种群中支配个体i的个体数量，Si表示被个体i支配的全部个体集合；
Step3-3-2检索所有ni＝0的个体，将他们存入子集F1，此即等级为1的非支配个体集合，对该集合...

【专利技术属性】
技术研发人员：郭玉雪，许月萍，尤延锋，谢京凯，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33