一种梯级水库群优化调度新型群体搜索方法

一种梯级水库群优化调度新型群体搜索方法
【专利摘要】本发明涉及水电系统发电调度领域，公开了一种梯级水库群优化调度新型群体搜索方法，结合精英集合动态更新策略与邻域变异搜索机制，以均衡考虑方法的全局搜索与局部勘探，兼顾种群的多样性与方法的收敛速度。其技术方案为：将水电站水位作为状态变量，以调度期内梯级水库群综合发电量最大为优化目标，在初始化一定数量的蜘蛛个体后，逐代实行子蛛群的内部协作行为、异性个体的婚配行为、精英个体动态更新与邻域变异搜索策略，逐次逼近梯级水库群的最优调度策略。其中的精英个体动态更新策略可以保证精英蛛群引导种群进行有效进化，平衡方法的搜索能力与勘探能力；优秀个体邻域变异策略可以维持种群多样性，提升方法的计算效率与收敛速度，所提方法具有良好的推广价值与应用前景。
【专利说明】
-种梯级水库群优化调度新型群体搜索方法
技术领域
[0001] 本发明设及水电系统发电调度领域，特别设及一种梯级水库群优化调度新型群体 捜索方法。 技术背景
[0002] 伴随水电开发的平稳有序推进，中国逐步形成巨型梯级水库群新格局。梯级水库 群的协同统一调度有利于水能资源的合理高效利用和电力系统的安全稳定运行，可W产生 显著的经济、社会与生态等综合效益，使得梯级水库群优化调度是当前乃至W后水电能源 系统重要的理论与实践课题。但是，在梯级水库间复杂的水力与电力联系并存，且需综合考 虑水位、出库、出力等多种复杂时空约束，梯级水库群优化调度问题也成为典型的多阶段多 约束非线性最优控制问题。现有解决方案大致分为两类：W线性规划、动态规划、大系统分 解协调等为代表的传统优化方法；W布谷鸟算法、蛙跳算法、粒子群算法等为代表的人工智 能方法。尽管上述方法已取得了较为丰富的研究成果，但在工程实际应用中仍然受到不同 程度的限制，如动态规划的维数灾难题、粒子群算法的早熟收敛问题等。因此，梯级水库群 优化调度领域迫切需要研究新型方法W实现更快、更好的优化调度。
[0003] 受自然界蜘蛛群体的协同机制启发，已有学者提出了一种新型高效的智能优化方 法一一蛛群优化方法，该方法在多个标准测试函数均表现出更为优越的捜索能力，逐渐在 数值优化、网络设计等领域表现出强劲的生命力。但是在应用该方法求解梯级水库群优化 调度问题时，仍存在早熟收敛等问题。因此，本发明成果依托国家自然科学基金委重大国际 合作巧1210014)和国家十二五科技支撑计划项目（2013BAB06B04)，W欄沧江梯级水库群优 化调度问题为背景，W其中下游"二库五级"梯级水库群为主要对象，发明了具有高效性和 实用性的梯级水库群优化调度新型群体捜索方法。

【发明内容】

[0004] 本发明要解决的技术问题是提供一种梯级水库群优化调度新型群体捜索方法，结 合精英集合动态更新策略与邻域变异捜索机制，W均衡考虑方法的全局捜索与局部勘探， 兼顾种群的多样性与方法的收敛速度。
[0005] 本发明的技术方案为：本发明掲示了一种梯级水库群优化调度新型群体捜索方 法，按照下述步骤(1 )-(7)完成梯级库群协调优化过程：
[0006] (1)设置蛛群个体数目为化，最大迭代次数K、精英集合个体数目CO等计算参数；
[0007] (2)记迭代次数k=l，利用式①计算雄、雌两类性征子蛛群的个体数量，然后利用 式②随机产生雌性初始蛛群F和雄性初始蛛群M，实现蛛群所有个体的初始化；
[000引
[0009]
[0010]式中，Nf、Nm分别为雌性、雄性子蛛群个体数目术为[(U]区间均匀分布的随机数； []表示取整函数;j表示维度标号，j = 1，2，…，D;巧;^表示第k轮迭代时第i个雌性蜘蛛的第j 维取值，i = l，2,…，Nf; AC表示第k轮迭代时第1个雄性蜘蛛的第j维取值，1 = 1，2,…，Nm; 化〇为[(U]区间均匀分布的随机数；式、&分别为第j维变量取值上、下限；
[001U 此时蛛群Sk由Np个蜘蛛构成，包含了Nf个雌性蜘蛛构成的群Fk与Nm个雄性蜘蛛构成 的群Mk,记作
[0012] (3)将各蜘蛛个体修正至可行域并评估计算其适应度，利用式③获得个体相应权 重；
[0013]
[0014] 式中，J(Si)、wi分别表示第i个蜘蛛Si的适应度及其权重；
[0015] (4)采用精英集合动态更新策略动态更新精英集合骑，具体操作如下：
[0016] 设定起的精英个体数目取为《，精英集合的动态更新具体方案为:将当前种群中 适应度排在前《名的个体复制至中间集合V中，若的=0，则令怒=V ;否则取媒与V中较为 优秀的一半个体构成新的精英集合巧。同时，雌、雄蜘蛛进化公式分别如下：
[0017] ④雌性蜘蛛进化公式为：
[001 引
[0019]式中，Sa表示从呜所选第a个精英个体，日=[郎"]术1，印为[0山区间均匀分布 的随机数。
[0020] 间雄忡咖晚讲化公立责.
[0021]
[002^ 式中；站表示从礎所选第0个精英个体，0= [ru ? ];ri3，ri4为[0，U区间均匀分布 的随机数。
[0023] 化）肃I巧式减巧仍连行种群邻城巧择榜泰策胳Pi增輔力?法的勘探能力；
[0024]
[00巧]⑦Xi = argmax{F(Xi，），F(Xi)I
[0026] 式中，ri5为[0，1]区间均匀分布的随机数；丫为服从N(0，1)正态分布的随机数；
[0027] (6)分别采用式④、⑤产生新的雌性与雄性子蛛群，W实现子蛛群的内部协作行 为；
[0028] (7)采用式⑧~⑩对雄性子蛛群中支配个体实行婚配行为，W提升种群的多样性；
[0029]
[0030]
[0031]
[0032] 式中，R为婚配半径;Mf为第k轮迭代时第I个雄性蜘蛛个体;Tm为所有在婚配半 径R内的雌性个体组成子蛛群;Ps为Tm重所有个体的分配概率；
[0033] (8)令k = k+l，若k《K，则返回步骤(3);反之，停止计算，并输出最优个体。
[0034] 本发明对比现有技术有如下有益效果:本发明一种梯级水库群优化调度新型群体 捜索方法，将水电站水位作为状态变量，W调度期内梯级水库群综合发电量最大为优化目 标，在初始化一定数量的蜘蛛个体后，逐代实行子蛛群的内部协作行为、异性个体的婚配行 为、精英个体动态更新与邻域变异捜索策略，逐次逼近梯级水库群的最优调度策略。对比现 有技术，本发明中精英个体动态更新策略可W保证精英蛛群引导种群进行有效进化，平衡 方法的捜索能力与勘探能力；优秀个体邻域变异策略可W维持种群多样性，提升方法的计 算效率与收敛速度;本发明方法能够有效提升蛛群优化方法的捜索性能在实际工程中具有 良好的推广价值与应用前景。
【附图说明】
[0035] 图1是本发明方法基本原理示意图；
[0036] 图2(a)是小湾电站的优化调度计算结果；
[0037] 图2(b)是漫湾电站的优化调度计算结果；
[0038] 图2(c)是大朝山电站的优化调度计算结果；
[0039] 图2(d)是懦扎渡电站的优化调度计算结果；
[0040] 图2(e)是景洪电站的优化调度计算结果。
【具体实施方式】
[0041] 下面结合附图和实施例对本发明作进一步的描述。
[0042] 梯级水库群的协同统一调度有利于水能资源的合理高效利用和电力系统的安全 稳定运行，可W产生显著的经济、社会与生态等综合效益，但是，在梯级水库间复杂的水力 与电力联系并存，且需综合考虑水位、出库、出力等多种复杂时空约束，使得梯级水库群优 化调度面临着"维数灾"、"早熟收敛"等严重问题，尤其是求解大规模梯级水库群优化调度 问题非常困难。本发明掲示一种梯级水库群优化调度新型群体捜索方法，将水电站水位作 为状态变量，W调度期内梯级水库群综合发电量最大为优化目标，在初始化一定数量的蜘 蛛个体后，逐代实行子蛛群的内部协作行为、异性个体的婚配行为、精英个体动态更新与邻 域变异捜索策略，逐次逼近梯级水库群的最优调度策略。其中的精英个体动态更新策略可 W保证精英蛛群引导种群进行有效进化，平衡方法的捜索能力与勘探能力;优秀个体邻域 变异策略可W维持种群多样性，提升方法的计算效率与收敛速度，所提方法具有良好的推 广价值与应用前景。
[0043] 本发明的一种梯级水库群优化调度新型群体捜索方法，W梯级水库群优化调度中 常见的发电量最大模型开展研究，目标函数见下式：
[0044]
[0045] 式中，E为调度期内总发电量，kW ? h; N为电站数目；i为电站序号，i = l，2,…，N;T 为调度周期;t为时段序号，t = l，2,…，T;Pi,t为水库i在时段t的出力，kW;Ai为水库i的出力 系数;Qi,t为水库i在时段t的发电流量，m3/s;出,t为水库i在时段t的水头，m;At为时段t的小 时数，h。
[0046] 为保证优化结果的可行和可用性，上述目标函数需要考虑梯级水库群优化调度问 题求解过程中要满足的各类运行约束，下面将一一介绍。
[0047] (1)始末水位约束：
[004引
[0049] 式中，4为水库i的初始水位，m;远为水库i的末水位，m。
[0050] (2)水量平衡约束：
[0化1 ]
[0052] 式中，Vi,t为水库i在时段t的库容，m3; Qi表示水库i的上游水库集合，对龙头水库 有^=0;qi,t、Oi,t、di,t分别为水库i在时段t的区间流量、出库流量和弃水流量，m3/s。
[0053] (3)水位约束：
[0化4] ,,
[00对 C库i在时段t的水位，m毎,t、挺为水库i在时段t的水位上、下限。
[0化6] 约束；
[0057] t
[0化引 t为水库i在时段t的发电流量上、下限。
[0化9] 约束；
[0060] ,
[0061] t为水库i在时段t的出库流量上、下限。
[0062] ；
[0063]
[0064] 为水库i在时段t的出力上、下限。
[00化] 约束：
[0066]
[0067] 巧甲，々,、心分别为时段t的水电系统带宽上、下限，kW。
[0068] 根据上述思想，一次完整的优化调度过程，按照下述步骤(1)-(8)予W实现：
[0069] (1)设置蛛群个体数目为化，最大迭代次数K、精英集合个体数目CO等计算参数；
[0070] (2)记迭代次数k=l，利用式①计算雄、雌两类性征子蛛群的个体数量，然后利用 式②随机产生雌性初始蛛群F和雄性初始蛛群M，实现蛛群所有个体的初始化；
[0071] (3)将各蜘蛛个体修正至可行域并评估计算其适应度，利用式③获得个体相应权 重；
[0072] (4)采用精英集合动态更新策略动态更新精英集合畴;
[0073] (5)利用式⑥和⑦实行种群邻域变异捜索策略W增强方法的勘探能力；
[0074] (6)分别采用式④、⑤产生新的雌性与雄性子蛛群，W实现子蛛群的内部协作行 为；
[0075] (7)采用式⑧~⑩对雄性子蛛群中支配个体实行婚配行为，W提升种群的多样性；
[0076] (8)令k = k+l，若k《K，则返回步骤(3);反之，停止计算，并输出最优个体。
[0077] 现W欄沧江中下游梯级水库群为研究对象，采用本发明方法制作梯级水库群优化 调度方案。欄沧江流域是中国13大水电基地之一，干流规划设计14级水电站，总装机规模约 25750MW，现阶段仅有部分电站投产运行。本文研究对象包括已投运的5座电站，分别为小 湾、漫湾、大朝山、懦扎渡和景洪，流域拓扑图和电站特征参数分别如表1所示。各电站优化 调度结果由图2(a)-(e)所示，不同方法优化调度结果对比如表2所示，不同方案分别代表不 同运行次数。由图2(a)-(e)分析可知，本发明方法计算结果中，小湾作为龙头水库，在蓄水 期迅速抬升水位，枯水期逐渐消落至设定水位，充分发挥了梯级补偿效应;其余水库在讯期 快速降低水位迎接大流量来水，其余时段尽可能维持高水头运行W降低水耗增发电量;各 水库优化结果均满足水位、出力等各项约束条件。运充分表明本发明方法所得调度过程合 理可行，能够服务于梯级水库群优化调度实际作业。由表2可知，从发电量看，本发明方法能 够获得多于粒子群和蛛群优化方法的发电量，在10中方案下，分别多于粒子群方法14.63、 23.82、19.06、15.61、27.67、12.08、17.43、16.71、9.67、17.41 亿 kWh,分别多于蛛群优化方 法 6.19、2.69、4.57、7.91、15.35、4.13、4.81、1.94、3.33、2.43 亿 kWh;从发电量标准差看，本 发明方法不同方案的标准差远小于其他两类对比方法，说明本发明方法捜索机制优于粒子 群方法捜索机制，同时所提策略能够有效完善蛛群优化方法捜索机理，提升了寻优结果的 稳定性与鲁棒性;从计算耗时看，本发明方法优于粒子群，与蛛群优化方法相近，说明本发 明方法可W有效提高梯级水电站群优化调度求解效率。
[007引 表1
[0079]
[0080]
[0081]
【主权项】
1. 一种梯级水库群优化调度新型群体搜索方法，其特征包括如下步骤： (1) 设置蛛群个体数目为Np，最大迭代次数K、精英集合个体数目ω等计算参数； (2) 记迭代次数k=l，利用式①计算雄、雌两类性征子蛛群的个体数量，然后利用式② 随机产生雌性初始蛛群F和雄性初始蛛群Μ，实现蛛群所有个体的初始化；式中，Nf、Nm分别为雌性、雄性子蛛群个体数目；^为[0，1]区间均匀分布的随机数；□表 示取整函数;j表示维度标号，j = l，2，-_，D;i^表示第k轮迭代时第i个雌性蜘蛛的第j维取 值，i = 1，2，…，Nf; 表示第k轮迭代时第1个雄性蜘蛛的第j维取值，1 = 1，2，…，Nm; r2，r3 为[〇，1 ]区间均匀分布的随机数；尤、心分别为第j维变量取值上、下限； 此时蛛群Sk由NP个蜘蛛构成，包含了Nf个雌性蜘蛛构成的群Fk与仏个雄性蜘蛛构成的群 Mk，记作= P =碎，…名,+2 p (3)将 各蜘蛛个体修正至可行域并评估计算其适应度，利用式③获得个体相应权重；1<7<,V v ' ^ .]</<^· 、·， 式中，·KSO、·分别表示第i个蜘蛛Si的适应度及其权重； (4) 采用精英集合动态更新策略动态更新精英集合4，具体操作如下： 设定允的精英个体数目取为ω，精英集合4动态更新具体方案为:将当前种群中适应 度排在前ω名的个体复制至中间集合ν中，若4 =必，则令圮=ν;否则取祀与ν中较为优秀的 一半个体构成新的精英集合4;伺时，雌、雄蜘蛛进化公式分别如下： ④ 雌性蜘蛛进化公式为：式中，Sa表示从巧所选第α个精英个体，α=[；τ?2ω ] ;rii，ri2为[0,1]区间均勾分布的随机 数； ⑤ 雄性蜘蛛进化公式为：式中；Sfs表示从#.所选第β个精英个体，β=[η4ω ] ;ri3，ri4为[0,1]区间均勾分布的随机 数； (5) 利用式⑥和⑦实行种群邻域变异搜索策略以增强方法的勘探能力；⑦Xi = argmax{F(Xi，)，F(Xi)} 式中，r15为[Ο，1 ]区间均匀分布的随机数；γ为服从N( Ο，1)正态分布的随机数； (6) 分别采用式④、⑤产生新的雌性与雄性子蛛群，以实现子蛛群的内部协作行为； (7) 采用式⑧~⑩对雄性子蛛群中支配个体实行婚配行为，以提升种群的多样性；式中，R为婚配半径; < 为第k轮迭代时第1个雄性蜘蛛个体;TM为所有在Mf婚配半径R 内的雌性个体组成子蛛群;Ps为TM重所有个体的分配概率； (8) 令k = k+Ι，若k彡K，则返回步骤(3);反之，停止计算，并输出最优个体。
【文档编号】G06Q50/06GK105956714SQ201610340489
【公开日】2016年9月21日
【申请日】2016年5月21日
【发明人】吴英, 程春田, 李江, 周毅, 冯仲恺, 申建建, 牛文静
【申请人】华能澜沧江水电股份有限公司, 大连理工大学