发明名称 一种梯级水库群优化调度新型群体搜索方法
摘要 本发明涉及水电系统发电调度领域,公开了一种梯级水库群优化调度新型群体搜索方法,结合精英集合动态更新策略与邻域变异搜索机制,以均衡考虑方法的全局搜索与局部勘探,兼顾种群的多样性与方法的收敛速度。其技术方案为:将水电站水位作为状态变量,以调度期内梯级水库群综合发电量最大为优化目标,在初始化一定数量的蜘蛛个体后,逐代实行子蛛群的内部协作行为、异性个体的婚配行为、精英个体动态更新与邻域变异搜索策略,逐次逼近梯级水库群的最优调度策略。其中的精英个体动态更新策略可以保证精英蛛群引导种群进行有效进化,平衡方法的搜索能力与勘探能力;优秀个体邻域变异策略可以维持种群多样性,提升方法的计算效率与收敛速度,所提方法具有良好的推广价值与应用前景。
申请公布号 CN105956714A 申请公布日期 2016.09.21
申请号 CN201610340489.6 申请日期 2016.05.21
申请人 华能澜沧江水电股份有限公司;大连理工大学 发明人 吴英;程春田;李江;周毅;冯仲恺;申建建;牛文静
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉
主权项 一种梯级水库群优化调度新型群体搜索方法,其特征包括如下步骤:(1)设置蛛群个体数目为Np,最大迭代次数K、精英集合个体数目ω等计算参数;(2)记迭代次数k=1,利用式①计算雄、雌两类性征子蛛群的个体数量,然后利用式②随机产生雌性初始蛛群F和雄性初始蛛群M,实现蛛群所有个体的初始化;<img file="FDA0000996198750000011.GIF" wi="518" he="151" /><img file="FDA0000996198750000012.GIF" wi="518" he="166" />式中,N<sub>f</sub>、N<sub>m</sub>分别为雌性、雄性子蛛群个体数目;r<sub>1</sub>为[0,1]区间均匀分布的随机数;[]表示取整函数;j表示维度标号,j=1,2,…,D;<img file="FDA0000996198750000013.GIF" wi="68" he="70" />表示第k轮迭代时第i个雌性蜘蛛的第j维取值,i=1,2,…,N<sub>f</sub>;<img file="FDA0000996198750000014.GIF" wi="83" he="70" />表示第k轮迭代时第l个雄性蜘蛛的第j维取值,l=1,2,…,N<sub>m</sub>;r<sub>2</sub>,r<sub>3</sub>为[0,1]区间均匀分布的随机数;<img file="FDA0000996198750000015.GIF" wi="86" he="67" />X<sub>j</sub>分别为第j维变量取值上、下限;此时蛛群S<sup>k</sup>由N<sub>p</sub>个蜘蛛构成,包含了N<sub>f</sub>个雌性蜘蛛构成的群F<sup>k</sup>与N<sub>m</sub>个雄性蜘蛛构成的群M<sup>k</sup>,记作:<img file="FDA0000996198750000016.GIF" wi="1614" he="79" />(3)将各蜘蛛个体修正至可行域并评估计算其适应度,利用式③获得个体相应权重;<img file="FDA0000996198750000017.GIF" wi="558" he="151" />式中,J(S<sub>i</sub>)、w<sub>i</sub>分别表示第i个蜘蛛S<sub>i</sub>的适应度及其权重;(4)采用精英集合动态更新策略动态更新精英集合<img file="FDA0000996198750000018.GIF" wi="74" he="62" />具体操作如下:设定<img file="FDA0000996198750000019.GIF" wi="56" he="60" />的精英个体数目取为ω,精英集合<img file="FDA00009961987500000110.GIF" wi="52" he="61" />动态更新具体方案为:将当前种群中适应度排在前ω名的个体复制至中间集合ν中,若<img file="FDA00009961987500000111.GIF" wi="163" he="62" />则令<img file="FDA00009961987500000112.GIF" wi="147" he="61" />否则取<img file="FDA00009961987500000113.GIF" wi="51" he="61" />与ν中较为优秀的一半个体构成新的精英集合<img file="FDA00009961987500000114.GIF" wi="73" he="63" />同时,雌、雄蜘蛛进化公式分别如下:④雌性蜘蛛进化公式为:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>F</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>+</mo><msub><mi>r</mi><mn>4</mn></msub><msub><mi>V</mi><msub><mi>c</mi><mi>i</mi></msub></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>c</mi></msub><mo>-</mo><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mn>5</mn></msub><msub><mi>V</mi><msub><mi>b</mi><mi>i</mi></msub></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>b</mi></msub><mo>-</mo><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mn>6</mn></msub><mrow><mo>(</mo><msub><mi>r</mi><mn>7</mn></msub><mo>-</mo><mn>0.5</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mn>11</mn></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>&alpha;</mi></msub><mo>-</mo><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>r</mi><mn>7</mn></msub><mo>&le;</mo><msub><mi>P</mi><mi>f</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>-</mo><msub><mi>r</mi><mn>4</mn></msub><msub><mi>V</mi><msub><mi>c</mi><mi>i</mi></msub></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>c</mi></msub><mo>-</mo><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mn>5</mn></msub><msub><mi>V</mi><msub><mi>b</mi><mi>i</mi></msub></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>b</mi></msub><mo>-</mo><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mn>6</mn></msub><mrow><mo>(</mo><msub><mi>r</mi><mn>7</mn></msub><mo>-</mo><mn>0.5</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mn>11</mn></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>&alpha;</mi></msub><mo>-</mo><msubsup><mi>F</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>r</mi><mn>7</mn></msub><mo>&gt;</mo><msub><mi>P</mi><mi>f</mi></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00009961987500000115.GIF" wi="1418" he="167" /></maths>式中,S<sub>α</sub>表示从<img file="FDA00009961987500000116.GIF" wi="52" he="61" />所选第α个精英个体,α=[r<sub>12</sub>ω];r<sub>11</sub>,r<sub>12</sub>为[0,1]区间均匀分布的随机数;⑤雄性蜘蛛进化公式为:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>M</mi><mi>l</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>M</mi><mi>l</mi><mi>k</mi></msubsup><mo>+</mo><msub><mi>r</mi><mn>8</mn></msub><msub><mi>V</mi><msub><mi>g</mi><mi>l</mi></msub></msub><mrow><mo>(</mo><msub><mi>F</mi><mi>g</mi></msub><mo>-</mo><msubsup><mi>M</mi><mi>l</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mn>9</mn></msub><mrow><mo>(</mo><msub><mi>r</mi><mn>10</mn></msub><mo>-</mo><mn>0.5</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mn>13</mn></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>&beta;</mi></msub><mo>-</mo><msubsup><mi>M</mi><mi>l</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>w</mi><mrow><msub><mi>N</mi><mi>f</mi></msub><mo>+</mo><mi>l</mi></mrow></msub><mo>&gt;</mo><msub><mi>w</mi><mrow><msub><mi>N</mi><mi>f</mi></msub><mo>+</mo><mi>m</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>M</mi><mi>l</mi><mi>k</mi></msubsup><mo>+</mo><msub><mi>r</mi><mn>8</mn></msub><mrow><mo>(</mo><msup><mi>C</mi><mi>k</mi></msup><mo>-</mo><msubsup><mi>M</mi><mi>l</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mn>13</mn></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>&beta;</mi></msub><mo>-</mo><msubsup><mi>M</mi><mi>l</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>w</mi><mrow><msub><mi>N</mi><mi>f</mi></msub><mo>+</mo><mi>l</mi></mrow></msub><mo>&le;</mo><msub><mi>w</mi><mrow><msub><mi>N</mi><mi>f</mi></msub><mo>+</mo><mi>m</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000996198750000021.GIF" wi="1318" he="174" /></maths>式中;S<sub>β</sub>表示从<img file="FDA0000996198750000022.GIF" wi="49" he="55" />所选第β个精英个体,β=[r<sub>14</sub>ω];r<sub>13</sub>,r<sub>14</sub>为[0,1]区间均匀分布的随机数;(5)利用式⑥和⑦实行种群邻域变异搜索策略以增强方法的勘探能力;<img file="FDA0000996198750000023.GIF" wi="1102" he="94" />⑦X<sub>i</sub>=argmax{F(X<sub>i</sub>'),F(X<sub>i</sub>)}式中,r<sub>15</sub>为[0,1]区间均匀分布的随机数;γ为服从N(0,1)正态分布的随机数;(6)分别采用式④、⑤产生新的雌性与雄性子蛛群,以实现子蛛群的内部协作行为;(7)采用式⑧~⑩对雄性子蛛群中支配个体实行婚配行为,以提升种群的多样性;<img file="FDA0000996198750000024.GIF" wi="438" he="127" /><img file="FDA0000996198750000025.GIF" wi="1006" he="98" /><img file="FDA0000996198750000026.GIF" wi="462" he="111" />式中,R为婚配半径;<img file="FDA0000996198750000027.GIF" wi="69" he="61" />为第k轮迭代时第l个雄性蜘蛛个体;T<sub>M</sub>为所有在<img file="FDA0000996198750000028.GIF" wi="66" he="63" />婚配半径R内的雌性个体组成子蛛群;P<sub>s</sub>为T<sub>M</sub>重所有个体的分配概率;(8)令k=k+1,若k≤K,则返回步骤(3);反之,停止计算,并输出最优个体。
地址 650000 云南省昆明市昆明世纪城中路1号