一种基于改进型人工蜂群算法的配电网无功优化方法
【专利摘要】本发明提出一种基于改进型人工蜂群算法的配电网无功优化方法，具体包括：建立配电网无功优化模型，并设置约束条件；设置控制变量和状态变量；初始化种群；随机产生种群的初始位置，并计算其适应度值；局部寻优；更新蜜源位置并计算其适应度值，保留最优解；判断终止条件，输出最优解。本发明将差分算法中的交叉和变异机制引入到人工蜂群算法中，适用配电系统的无功优化，加快了收敛速度，提高人工蜂群算法的局部搜索能力，降低网络损耗的效果，并降低电网的运行成本。
【专利说明】
-种基于改进型人工蜂群算法的配电网无功优化方法
技术领域
[0001] 本发明属于电力信息技术领域，尤其是一种基于改进型人工蜂群算法的配电网无 功优化方法。
【背景技术】
[0002] 电力系统无功优化指给定已知量包括有功负荷、有功电源、潮流分布，通过优化计 算确定系统中能够满足所有约束条件的系统最小有功网损值W及质量电压最优值等。系统 的无序流动，将会增加其有功损耗，当系统一旦发生无功不足W及扰动现象时，将使电压水 平低下，当低于临界电压时，会产生电压崩溃，从而导致系统因失去同步而崩溃的灾难性事 故。
[0003] 无功优化问题是一个多变量、多约束的非线性数学规划问题。它的操作系统包括 连续变量和离散变量，其中连续变量包括节点电压;离散变量按有载调压分接头档位、补偿 电容器的投切组数。进行无功优化潮流即电压控制一般要对发电机端电压、可调变压器变 比和节点步长无功做综合调整，W网损最小作为目标函数，对无功功率进行优化控制，可W 改善系统电压质量，降低有功功率损耗，节约运行成本。在满足系统各种约束的条件下，来 实现无功功率在电网中的优化配置，使电网的有功损耗最小。
[0004] 人工蜂群算法是一种新型的群智能算法，它具有原理简单、控制参数少、鲁棒性强 而被广泛应用。但是该算法在捜索阶段存在解的利用不足，收敛速度慢，容易陷入局部最优 等问题，为了解决运一问题。
【发明内容】
[0005] 本发明所解决的技术问题在于提供一种基于改进型人工蜂群算法的配电网无功 优化方法，将差分算法中的交叉和变异机制引入到人工蜂群算法中，加快了收敛速度，提高 人工蜂群算法的局部捜索能力，降低网络损耗的效果，并降低电网的运行成本。
[0006] 实现本发明目的的技术解决方案为：
[0007] -种基于改进型人工蜂群算法的配电网无功优化方法，包括W下步骤：
[000引步骤1:建立配电网无功优化模型，并设置约束条件，约束条件包括功率约束和变 量约束；
[0009] 步骤2:设置无功优化模型的控制变量和状态变量；
[0010] 步骤3:初始化种群，设置种群参数；
[0011] 步骤4:采用第一位置更新公式随机产生种群的初始位置，并计算其适应度值，用 贪婪选择机制更新最优蜜源，所述第一位置更新公式为：
[0012] xi,d=mind+ran山（0,1) X (maxd-mind)，
[0013]其中，每一个解Xi,d都用D维向量表示，mind、maxd分别表示为维数的下界、上界， randi(0,l)表示0至Ijl之间均匀分布的随机值，i = l,2,. . . ,NP，NP表示种群大小；
[0014]步骤5:跟随蜂根据概率公式并采用轮盘赌原则选择某个雇佣蜂跟随，并在该雇佣 蜂所在的蜜源附近进行寻优；
[0015] 步骤6:生成一个均匀分布的随机值rand2(0,1)，
[0016] 若rand2(0，l)<prob(i)，则使用第二位置更新公式来更新vi,d，所述第二位置更 新公式为：
[0017] vi,d = Xa,d+L(A) X (xb,d-Xc,d),
[001引其中，VM表示更新后的蜜源位置，a，b，cG[l，N門且互不相等，L(A)表示为服从参 数入的莱维分布产生一个随机捜索向量，即
'如为第t次迭代保留的的最
优蜜源，U和V均服从正态分布 曰V=I，r表示为Gamma函数；
[0019] 若^11(12(0，1)>9的6。），则使用第^位置更新公式来更新巾1,<1，所述第^位置更 新公式为：
[0020] vi,d = Xi,d+ri,dX (xk,d-Xi,d),
[0021] 其中，XM是当前个体蜜源位置，xk,d是随机选择的蜜源位置，k= I，2…NP, i声k， Ti,d是随机数，Ti,dG [-1，1]，d=l ,2.Jmax表示最大 迭代次数，t表示目前的迭代次数；
[0022] 步骤7:计算每一个种群的适应度值，并采用贪婪选择机制更新最优蜜源；
[0023] 步骤8:根据限定值limit判断种群是否陷入局部最优，若种群适应度值在经过 limit次迭代后仍然没有得到改进，则判断该种群陷入局部最优，采用第一位置更新公式进 行随机捜索并计算适应度值，采用贪婪选择机制更新最优蜜源;若否，转到步骤9;
[0024] 步骤9:判断是否达到最大迭代次数，若是，则输出最优解;若否，则转到步骤4。
[0025] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤1中无 功优化模型的目标函数为：
[0026]
[0027] 其中，Ui为节点i处的电压，Uj为节点j处的电压，0ij为节点i、j之间的相角差，Gij为 节点i、j之间的电导，N为节点总数，Pioss表示无功损耗功率。
[0028] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤1中的 功率约束为：
[0029；
[0030；
[0031] 其中，Pi、Qi、Ui分别为节点i处注入的有功功率、无功功率和电压，Uj为节点j处的 电压，Gu、Bu、0i汾别为节点i、j处的电导、电纳和相角差，N为节点总数。
[0032] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤1中的 变量约束为：
[0033]
[0034] 其中，UGk为发电机的电压，UGk,min、UGk,max为可调发电机的上下限，Gk= I ,2 ,…Ng，Ng 为所有可调发电机节点总数;Cl为不成电容器的投切组数，Cl,min、Ci,max为补偿电容器投切组 数的上下限，il=l，2,…Ne, Nc为所有无功补偿节点总数;T功变压器的分接头档位，L,min、 Tj,max为可调变压器分接头档位的上下限，j = 1，2，…化，化为所有变压器支路总数;Ui为PQ节 点的电压，化如。、化,。3功?9节点电压的上下限，12=1，2，一化，化为?9节点总数;舶功发电机 的无功出力，0时,。1。、9时,。3、为发电机无功出力的上下限，13=1，2^''抓。
[0035] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤2中的 控制变量包括发电机端电压、有载调压变压器分接头位置、并联电容器和电抗器投切组 数。
[0036] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤2中的 状态变量包括除平衡节点外的所有节点的电压相角、除发电机或具有无功补偿的设备外的 节点的电压模值。
[0037] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤3中的 种群参数包括种群大小NP和控制变量数D，其中种群数的一半为雇佣蜂数量，另一半为跟随 蜂数量，控制变量数即维数，包括发电机端电压、有载调压变压器变比、无功补偿容量、最大 迭代次数MAX 口邸和限定值limit。
[0038] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤4和步 骤7中计算适应度值的公式为：
[0039]
[0040] 其中，fi表示目标函数，即配电系统的网损函数公式，fiti表示种群的适应度值。
[0041] 进一步的，本发明的基于改进型人工蜂群算法的配电网无功优化方法，步骤5中的 概率公式为：
[0042]
[0043] 其中，fiti表示种群的适应度值，NP表示种群大小。
[0044] 本发明采用W上技术方案与现有技术相比，具有W下技术效果：
[0045] 1、本发明将差分算法中的交叉和变异机制引入到人工蜂群算法中，提高了种群的 多样性，实现了多路径捜索，防止算法陷入局部最优.
[0046] 2、本发明增强了算法的局部捜索能力，加快了收敛速度，提高全局最优值概率；
[0047] 3、本发明适用于配电系统的无功优化，能够更大程度地减少电网的运行成本。
【附图说明】
[0048] 图1是本发明的基于改进型人工蜂群算法的配电网无功优化方法流程图；
[0049] 图2是本发明的IEEE30节点电力系统结构图。
【具体实施方式】
[0050] 下面详细描述本发明的实施方式，所述实施方式的示例在附图中示出，其中自始 至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过 参考附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
[0051] 本发明提出了一种基于改进型人工蜂群算法的配电网无功优化方法，其流程图如 图1所示，具体包括W下步骤：
[0052] 步骤1:建立配电网无功优化模型，并设置约束条件，约束条件包括功率约束和变 量约束。
[0053] 配电网无功优化模型是研究无功优化问题的基础，该配电网由发电机、变压器、输 电线和负荷等组成。
[0054] 无功优化模型的目标函数为：
[0化5]
[0056] 其中，Ui为节点i处的电压，Uj为节点j处的电压，0ij为节点i、j之间的相角差，Gij为 节点i、j之间的电导，N为节点总数，Pioss表示无功损耗功率。
[0057] 功率约束为：
[0化引 [0化9]
[0060] 其中，Pi、化、Ui分别为节点i处注入的有功功率、无功功率和电压，Gij、Bij、0u分别 为节点i、j处的电导、电纳和相角差，N为节点总数。
[0061] 变量约束为：
[0062]
[006引片:下，化止，UGk,min、UGk,max为可调发电机的上下限，Gk = 1,2 ,…Ng ,Ng 为所有可调发电机节点总数;Cl为补偿电容器的投切组数，Cl,min、Ci,max为补偿电容器投切组 数的上下限，il=l，2,…Ne, Nc为所有无功补偿节点总数;T功变压器的分接头档位，L,min、 Tj,max为可调变压器分接头档位的上下限，j = 1，2，…化，Nt为所有变压器支路总数;Ui为PQ 节点的电压，化如。、化,。3功?9节点电压的上下限，12=1，2，一化，化为?9节点总数;舶功发电 机的无功出力，Qci, min、QGi, max为发电机无功出力的上下限，13=1,2,…Ng。
[0064] 步骤2:设置无功优化模型的控制变量和状态变量。
[0065] 所述控制变量包括发电机端电压、有载调压变压器分接头位置、并联电容器和电 抗器投切组数。
[0066] 所述状态变量包括除平衡节点外的所有节点的电压相角、除发电机或具有无功补 偿的设备外的节点的电压模值。
[0067] 本实施例中采用的是IE邸30系统，如图2所示，该系统包括PV节点、PQ节点和平衡 节点、可调变压器的抽头、发电厂的无功出力、无功补偿装置。其中包括41条支路、21个负荷 节点、6台发电机、4台可调变压器及2个并联电容无功补偿点。
[0068] 步骤3:初始化种群，设置种群参数。
[0069] 人工蜂群算法有=种蜂，分别为雇佣蜂、跟随蜂、侦查蜂;其中它们的任务是在蜂 巢周围随机捜索蜜源，并比较捜索前后食物源的丰富程度;选择丰富的蜜源作为捜索目标 后会返回蜂巢与其他蜜蜂分享食物源信息，雇佣蜂会通过跳舞的形式传达蜜源的信息;跟 随蜂会根据雇佣蜂提供的信息按照概率选择优质的蜜源进行开采;越丰富的蜜源被选择的 概率越大，然后跟随蜂和雇佣蜂一样进行领域捜索，并选择较好的解。当一个蜜源经过多次 开采后仍没有改进，雇佣蜂将会放弃该蜜源并转成侦查蜂进行跳跃式随机捜索。
[0070] 所述种群参数包括种群大小NP和控制变量数D，其中，种群数的一半为雇佣蜂数 量，另一半为跟随蜂数量;控制变量数即维数，表示为待优化的参数，包括发电机端电压、有 载调压变压器变比、无功补偿容量和最大迭代次数。
[0071] 步骤4:采用第一位置更新公式随机产生种群的初始位置，并计算其适应度值，用 贪婪选择机制更新最优蜜源，所述第一位置更新公式为：
[0072] Xi,d=mind+randi(0,1) X (maxd-mind)，
[0073] 其中，每一个解xi,d都用D维向量表示，mind、maxd分别表示为维数的下界、上界， randi(0,l)表示0至Ijl之间均匀分布的随机值，i = l,2,. . . ,NP，NP表示种群大小。
[0074] 其中，计算适应度值的公式为：
[0075]
[0076] 其中，fi表示目标函数，即配电网的网损函数公式，fiti表示种群的适应度值。
[0077] 步骤5:跟随蜂根据概率公式并采用轮盘赌原则选择某个雇佣蜂跟随并在该蜜源 附近进行局部寻优；
[0078] 加:杰/.V 气为；
[0079]
[0080] 其中，fiti表示种群的适应度值，NP表示种群大小。
[0081 ] 步骤6:生成一个均匀分布的随机值rand2(0，1)，
[0082] 若rand2(0，l)<prob(i)，则使用第二位置更新公式来更新vi,d，所述第二位置更 新公式为：
[008；3] vi,d = Xa,d+L(入）X (xb,d-Xc,d)，
[0084] 其中，VM表示更新后的蜜源位置，a，b，cG [I，N門且互不相等，L(A)表示为服从参 数入的莱维分布产生一个随机捜索向量，目[
'xbt为第t次迭代保留的的 最优蜜源，U和V均服从正态分布：
， Ov=I，r表示为Gamma函数；
[0085] 若rand2(0，1)>prob(i)，则使用第S位置更新公式来更新Vi,d，所述第S位置更 新公式为：
[0086] Vi,d = Xi,d+ri,dX (xk,d-xi,d)，
[0087] 其中，Xi,d是当前个体蜜源位置，Xk,d是随机选择的蜜源位置，k= I，2…NP，i声k， ri,d是随机数，ri,dG [-1，1]，d=l ,2...D;p;rob(i)=a+(l-a)*(Tmax-t)/Tmax，Tmax表示最大 迭代次数，t表示目前的迭代次数。
[0088] 步骤7:计算每一个种群的适应度值，并采用贪婪选择机制，比较Xi,d和Vi,d的适应 度值，如果Xi, d优于Vi, d，则用Xi, d代替Vi, d，更新最优蜜源。
[0089] 步骤8:根据限定值limit判断种群是否陷入局部最优，若种群适应度值在经过 limit次迭代后仍然没有得到改进，则判断该种群陷入局部最优，采用第一位置更新公式进 行随机捜索并计算适应度值，采用贪婪选择机制更新最优蜜源;若否，转到步骤9。
[0090] 步骤9:判断是否达到最大迭代次数，若是，则输出最优解;若否，则转到步骤4。
[0091] 在本实施例中，采用如图2所示的IE邸30节点系统进行仿真。I邸E30节点为国际的 一个标准网络系统，应用官方给定的具体参数进行计算。I邸E30节点标准系统有41条支路、 21个负荷节点;6台发电机、4台可调变压器及2个并联电容无功补偿点。节点1，2,5,8,11，13 为发电机节点，在发电机节点中，节点1设为平衡节点；节点2，5，8，11，13为PV节点；其余的 为PQ节点。发电机端电压的取值范围在0.95~1.10之间连续取值;有载调压变压器调节范 围为0.9~1.1，分16个档，调节步长为0.0125;节点10、24为并联电容器，初始状态各补偿点 的并联电容器全部投入，无功补偿装置范围为[0,0.5]，调节补偿为0.5。
[0092] 本实施例取人工蜂群种群大小NP为20,雇佣蜂和跟随蜂各取10,维度D为12,最大 迭代次数MAX 口邸=150,限定值Iimit = IOO,采用粒子群优化算法所得的系统网络损耗值 为5.1734MW，采用传统人工蜂群优化算法所得的系统网络损耗值为4.2510MW，采用本发明 的改进型人工蜂群优化算法所得的系统网络损耗值为4.0326MW。通过比较发现，本发明所 采用的改进型人工蜂群算法具有更优的捜寻功能，能够有效的降低有功功率网络损耗值。
[0093] W上所述仅是本发明的部分实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可W做出若干改进，运些改进应视为本发明的保 护范围。
【主权项】
1. 一种基于改进型人工蜂群算法的配电网无功优化方法，其特征在于，包括以下步骤： 步骤1:建立配电网无功优化模型，并设置约束条件，约束条件包括功率约束和变量约 束； 步骤2:设置无功优化模型的控制变量和状态变量； 步骤3:初始化种群，设置种群参数； 步骤4:采用第一位置更新公式随机产生种群的初始位置，并计算其适应度值，用贪婪 选择机制更新最优蜜源，所述第一位置更新公式为： xi,d=mind+randi(0,1) X (maxa-mina), 其中，每一个解xi,d都用D维向量表示，mind、maxd分别表示为维数的下界、上界，randi (〇，1)表示0到1之间均匀分布的随机值，i = l，2,...，NP，NP表示种群大小； 步骤5:跟随蜂根据概率公式并采用轮盘赌原则选择某个雇佣蜂跟随，并在该雇佣蜂所 在的蜜源附近进行寻优； 步骤6:生成一个均勾分布的随机值rand2(0,1)， 若瓜11(12(0，1)<?仰13(1)，则使用第二位置更新公式来更新^,(1，所述第二位置更新公式 为： Vi, d 一 Xa, d+L (λ ) X ( Xb, d_Xc, d )， 其中，Vl,d表示更新后的蜜源位置，a，b，ce[l，NP]且互不相等，L(A)表示为服从参数λ 的莱维分布产生一个随机搜索向量，SP为第t次迭代保留的最优蜜 源，u和v均服从正态分布:u~N(0,〇u2)，V~Ν(0，σν 2))ν=1，Γ 表示为Ga_a函数； 若rand2(0，l)多prob(i)，则使用第三位置更新公式来更新Vl,d，所述第三位置更新公式 为： Vi, d - Xi，d+Γ?，d X ( Xk, d_Xi，d )， 其中，Xi,d是当前个体蜜源位置，Xk,d是随机选择的蜜源位置，k = l，2···ΝΡ，;?辛k，ri,d是 随机数，ri,de [-1,1]，d = 1,2'"D;prob(i) = α+(l-a)*(Tmax_t)/Tmax，Tmax表示最大迭代 次数，t表示目前的迭代次数； 步骤7:计算每一个种群的适应度值，并采用贪婪选择机制更新最优蜜源； 步骤8:根据限定值limit判断种群是否陷入局部最优，若种群适应度值在经过limit次 迭代后仍然没有得到改进，则判断该种群陷入局部最优，采用第一位置更新公式进行随机 搜索并计算适应度值，采用贪婪选择机制更新最优蜜源;若否，转到步骤9; 步骤9:判断是否达到最大迭代次数，若是，则输出最优解;若否，则转到步骤4。2. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤1中无功优化模型的目标函数为：其中，Ui为节点i处的电压，Uj为节点j处的电压，Θ i」为节点i、j之间的相角差，Gi j为节点 i、j之间的电导，N为节点总数，Plciss表示无功损耗功率。3. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤1中的功率约束为：其中，Pi、Qi、Ui分别为节点i处注入的有功功率、无功功率和电压，Uj为节点j处的电压， 6"』^、01汾别为节点1、」处的电导、电纳和相角差4为节点总数。4. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤1中的变量约束为：其中，UGk为发电机的电压，UGk,min、UGk,max为可调发电机的上下限，Gk= 1，2，…Ng，Ng为所 有可调发电机节点总数为不成电容器的投切组数，为补偿电容器投切组数的 上下限，11 = 1，2，-教，价为所有无功补偿节点总数;1'」为变压器的分接头档位，1'」,_、1'」,_ 为可调变压器分接头档位的上下限，j = 1，2，"·Ντ，Ντ为所有变压器支路总数;U^PQ节点的 电压，Ui,min、Ui,max为PQ节点电压的上下限，i 2=l，2,…Nl，Nl为PQ节点总数;QGi为发电机的无 功出力，QGi,min、QGi,max为发电机无功出力的上下限，i3 = 1，2,…Ng。5. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤2中的控制变量包括发电机端电压、有载调压变压器分接头位置、并联电容器和电 抗器投切组数。6. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤2中的状态变量包括除平衡节点外的所有节点的电压相角、除发电机或具有无功补 偿的设备外的节点的电压模值。7. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤3中的种群参数包括种群大小NP和控制变量数D，其中种群数的一半为雇佣蜂数量， 另一半为跟随蜂数量，控制变量数即维数，包括发电机端电压、有载调压变压器变比、无功 补偿容量、最大迭代次数MAXITER和限定值1 imi t。8. 根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤4和步骤7中计算适应度值的公式为：其中，fi表示目标函数，即配电系统的网损函数公式，行^表示种群的适应度值。9.根据权利要求1所述的基于改进型人工蜂群算法的配电网无功优化方法，其特征在 于，步骤5中的概率公式为：其中，fiti表示种群的适应度值，NP表示种群大小。
【文档编号】G06Q10/04GK105956706SQ201610300670
【公开日】2016年9月21日
【申请日】2016年5月9日
【发明人】岳东, 高 浩, 秦媛
【申请人】苏州泛能电力科技有限公司