发明名称 基于多目标离散粒子群的智能配电网故障恢复方法
摘要 基于多目标离散粒子群智能优化算法的智能故障恢复方法,包括以下步骤1)初始化参数;2)初始化离散粒子群优化算法中的位置、速度,并根据智能配电网孤岛划分以及负荷恢复的算法,对每个粒子计算相对于多目标函数的各个适应度值fitness<sub>k</sub>;3)基于适应度支配的概念,进行支配种群和非支配种群的归类;4)根据相应的更新规则,对支配种群中的粒子位置和速度进行更新;5)进行支配种群和非支配种群间的动态交换;6)检测是否达到最大迭代次数,若达到最大迭代次数,则跳转步骤(4),否则进入步骤7);7)输出最终的优化结果。
申请公布号 CN104112165A 申请公布日期 2014.10.22
申请号 CN201410209704.X 申请日期 2014.05.19
申请人 浙江工业大学 发明人 陈明军;冯杰;陈俊宇
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 基于多目标离散粒子群智能优化算法的智能故障恢复方法,包括以下步骤:步骤1,初始化参数。1.1确定离散粒子群优化算法中的最大迭代次数、粒子数、维度数、学习因子、惯性权重的取值以及各类目标函数的计算公式和相应参数,包含了三个方面的目标函数:1.1.1尽可能多的恢复失电负荷<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>max</mi><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><mi>D</mi></mrow></munder><msub><mi>&lambda;</mi><mi>i</mi></msub><msub><mi>P</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506928630000011.GIF" wi="1034" he="116" /></maths>其中,P<sub>i</sub>表示失电负荷的大小;λ<sub>i</sub>表示失电负荷i的权重系数,即表示负荷的优先等级;D代表整个系统的负荷集合;1.1.2尽可能较少开关的操作次数:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>min</mi><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Element;</mo><msub><mi>T</mi><mi>s</mi></msub></mrow></munder><msub><mi>K</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506928630000012.GIF" wi="1048" he="118" /></maths>其中,K<sub>k</sub>表示开关的状态,1表示闭合,0表示打开;T<sub>s</sub>表示整个开关的集合;1.1.3恢复后,整个系统的网损尽可能的小:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>min</mi><msub><mi>f</mi><mn>3</mn></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>N</mi><mi>l</mi></msub></mrow></munder><msubsup><mi>I</mi><mi>l</mi><mn>2</mn></msubsup><msub><mi>R</mi><mi>l</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506928630000013.GIF" wi="1047" he="123" /></maths>其中,f<sub>3</sub>表示形成供电孤岛的有功功率损耗,I<sub>l</sub>表示支路的电流,N<sub>l</sub>表示整个系统中的供电支路;R<sub>l</sub>表示支路的阻抗参数。步骤2,离散粒子的位置和速度初始化;2.1初始化离散粒子的位置值各个粒子的位置值即为各类开关的状态值,1表示闭合,0表示打开;因此,在约束条件的范围内,随机定义各类粒子的位置值,即随机的开断各个开关以形成相应的孤岛;约束条件主要包括:2.1.1节点的电压约束,即各个节点的电压应该保持在允许的范围内,介于电压容许的最大值和最小值之内。2.1.2支路的功率约束,即各条支路上的功率不应该超过其最大的允许功率。2.1.3一类负荷不失电约束,即故障后所形成的各个孤岛内,负荷等级为一类的负荷必须保证不失电。2.1.4孤岛的容量约束,即各个孤岛内的总负荷和总损耗之和不能超过该孤岛内所有的DG发电总量。2.2初始化离散粒子群的速度;随机对离散粒子群进行速度的初始化,并将速度的初始值限制在0到1之间;2.3计算各个粒子的适应值;根据各类粒子随机后得到的位置值,即开关的闭合或打开状态,输入到整个智能配电网系统中,形成相应的孤岛,并对满足约束条件的孤岛进行潮流计算,得出相应的网络损耗。并根据这些得到的结果计算是目标函数(1)‑(3)的值,得出各个粒子的适应度值fitness<sub>k</sub>;步骤3,种群的归类;基于适应度支配的概念,将整个种群分为两个子群,分别为非支配子群P_set和支配子群NP_set;其中P_set子群中粒子的个数为n<sub>1</sub>,NP_set中粒子的个数为n<sub>2</sub>;步骤4,粒子位置和速度的更新;4.1根据下面公式对粒子的速度值更新:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><msub><mi>v</mi><mi>k</mi></msub><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>=</mo><msup><msub><mi>v</mi><mi>k</mi></msub><mi>t</mi></msup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>p</mi><mi>best</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>k</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>g</mi><mi>best</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>k</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506928630000031.GIF" wi="1196" he="81" /></maths>其中,x<sub>k</sub>表示粒子的位置值,v<sub>k</sub><sup>t</sup>表示粒子k在第t次迭代计算时的速度值,p<sub>best</sub>表示每个粒子的历史最优值;g<sub>best</sub>为随机在非支配子群中选出的最优值;4.2根据下面的公式对粒子的位置值更新:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><msub><mi>x</mi><mi>k</mi></msub><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo><mi>ifrand</mi><mo>(</mo><mo>)</mo><mo>&lt;</mo><mi>sig</mi><mrow><mo>(</mo><msup><msub><mi>v</mi><mi>k</mi></msub><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo><mi>else</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506928630000032.GIF" wi="1134" he="166" /></maths>其中,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>sig</mi><mrow><mo>(</mo><msup><msub><mi>v</mi><mi>k</mi></msub><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msubsup><mi>v</mi><mi>k</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506928630000033.GIF" wi="1105" he="141" /></maths>步骤5,支配种群和非支配种群间的动态交换;对NP_set子群中的每个粒子与P_set子群中的每个粒子逐一进行比较,并进行相应的交换操作,以使得最终NP_set子群中的所有粒子为支配粒子,即适应值均占优于P_set子群中的粒子的适应值;并且,最后删除P_set子群和NP_set子群中重复的粒子;步骤6,方法的收敛性检验;检验整个离散粒子群优化算法是否达到最大的迭代次数,若达到迭代次数已经达到所设定的最大次数,则跳出迭代步骤,进入下一步骤(7);若还未达到最大迭代次数,则回到步骤(4);步骤7,输出最终的优化结果;根据最终所得到的结果,即粒子的状态值(开关状态值),得出最终的孤岛方案和负荷恢复方案,并输出最终的三维帕累托最优曲面,以供操作者选择相应的开关操作量以及负荷恢复方案。
地址 310014 浙江省杭州市下城区潮王路18号
您可能感兴趣的专利