1、DOI:10.11991/yykj.202111013网络出版地址：https:/ IEEE33 节点系统为例进行仿真验证，仿真结果证明本文改进算法具有寻优速度快、寻优成功率高的优点。关键词：重构；多目标；改进烟花算法；最小网损；最少开关动作；配电网；优化；柯西变异中图分类号：TM732文献标志码：A文章编号：1009671X(2023)04008406Optimal reconfiguration of distribution network based on improvedfireworks algorithmPENGBo1,ZOULe1,LIXin1,ZHANGYuchen1,CHE
2、NMule1,TANDashuai21.StateGridBeijingHaidianElectricPowerSupplyCompany,Beijing100089,China2.DongfangElectronicsCorporation,Yantai264000,ChinaAbstract:Inordertoimprovetheefficiencyofdistributionnetworkreconfiguration,amulti-objectiveoptimizationfunctionisconstructedtoreducenetworklossandminimizeswitch
3、ingtimes.Inordertosolvetheproblemsoflowaccuracyandslowspeedoftheoriginalfireworksalgorithmandothertraditionalintelligentoptimizationalgorithms,animprovedfireworksalgorithmisproposed.Thealgorithmadoptsrealnumbercodingtoreducethevariabledimension,andemploystheloopswitcharraytodeterminethefireworksindi
4、vidualexplosionspace,soastoreducetheinvalidsearchintheoptimization.Atthesametime,byintroducingCauchymutationtoreplacetraditionalGaussianmutation,anadaptiveexplosionradiusandaneliteopposition-basedlearningselectionstrategyareproposed,whichhasgreatlyimprovedthesolutionspeedandaccuracyoftheoptimization
5、algorithm.Finally,theIEEE33bussystemistakenasanexampletoverifythattheimprovedalgorithmhastheadvantagesoffastoptimizationspeedandhighoptimizationsuccessrate.Keywords:reconfiguration;multiple targets;improved fireworks algorithm;minimum network loss;minimumswitchingaction;distributionnetwork;optimizat
6、ion;Cauchyvariation网络拓扑重构是配电网优化的重要手段之一，它是在满足一定的约束条件下，通过改变联络开关与分段开关的开闭状态从而改变网络拓扑。配电网重构是 Merlin 等1在 1975 年首次提出的概念。目前配电网重构方法主要有数学规化法、启发式搜索法和智能算法。数学规化法23面对复杂配电网时会产生维数灾问题，导致求解量和计算用时大量增加。启发式算法主要包括最优流模式算法45和支路交换法67：最优流模式算法具有速度快的优点，但其在求解过程容易产生孤立节点；支路交换法则容易受到配电网初始状态的影响且全局寻优能力不强。随着人工智能技术的不断发展，智能算法被用于配电网重构，并收稿
7、日期：20211104.网络出版日期：20230613.基金项目：山东省自然科学基金项目（ZR202103030510）；国网北 京 市 电 力 公 司 科 技 项 目（SGBJHD00DQJS2100459）.作者简介：彭博，男，工程师.谭大帅，男，高级工程师.通信作者：谭大帅，E-mail：.第50卷第4期应用科技Vol.50No.42023年7月AppliedScienceandTechnologyJul.2023取得一系列的成果，如遗传算法、粒子群算法、模拟退火法、和声搜索算法和禁忌搜索等。文献 8构造了网络损耗和电压偏差最小的多目标配电网优化模型，并利用烟花算法（fireworksa
8、lgorithm，FWA）进行求解。文献 9 提出了一种改进型二进制粒子群优化（binaryparticleswarmoptimization，BPSO）算法来快速调整拓扑结构，以显著提高系统在风险状态下的稳定性，但没考虑算法的寻优成功率，使其难以满足实际中故障重构对限定时间内的寻优成功率要求。为了避免或减少寻优过程中产生不可行解，文献 10 提出了基本环矩阵的概念，研究结果表明它能够简化编码方式、提高搜索效率。但是其将已确定元素所在基本环路与未确定维所在基本环路的相同元素删除的操作会剔除部分可行解。针对上述问题，本文以降低网损和开关次数最少为优化目标，并对原始烟花算法进行改进，采用实数编码来
9、减少变量维数，以基本环路矩阵划分而成的环路数组作为烟花算子的爆炸空间，来实现在大幅减少不可行解数量的同时保证其最优解不被丢失。通过算例验证表明，本文提出的改进后的 FWA 配电网重构算法，在降低配电网运行网损的同时，可快速得到多目标最优解。1配电网优化重构模型1.1目标函数配电网重构是根据配电网现有参数，通过改变分段开关和联络开关的通断状态，使原有配电网络既满足运行约束的要求，又达到某些性能的最优。因此，配电网重构是一个多目标多约束的优化求解问题。为便于算法性能的比较，本文仍构建降低网损和减少开关次数的多目标优化模型。1）网络损耗子目标函数的表达式11为min f1=ni=1P2i+Q2iU2
10、iRiKiniRiiKiiPiiQiiUii式中：为电网分支总数；为分支编号；为分支的电阻；为分支 的状态量（0 为断开，1 为闭合）；为分支 末端的有功功率；为分支 末端的无功功率；为分支 末端的电压。2）开关动作次数子目标函数的表达式为min f2=ND+NKNDNK式中：为网络中分段开关的动作次数，为网络中联络开关动作次数。将上述 2 个目标函数通过赋予不同的权重得到加权混合优化模型：min f=1f1+2f2121+2=1式中、分别为 2 个目标函数的权重系数，且。1.2约束条件配电网重构还需要满足配电网各种运行约束条件，如潮流方程约束、电压越限约束、辐射状约束3。Pi=PL,i+Ui
11、Ni=1Uj(Gijcosij+Bijsinij)Qi=QL,i+UiNi=1Uj(Gijcosij+Bijsinij)PiQiLiPL,iQL,iiUiUjLGijBijijLLN式中：、分别为在电网分支 的 节点注入的有功和无功功率，、分别为 节点的负荷有功和无功功率，和为分支 的 2 个端部节点的电压幅值，、分别为电网分支 的电导、电纳和分支 的 2 个端部节点的相角差，是电网的总节点数。Ui,min Ui Ui,maxUi,miniUi,maxi式中：为节点 的电压下限，为节点 的电压上限。gkGgkG式中：为重构后的网络结构，为所有满足辐射状要求的网络集合。2环路开关数组在配电网重构
12、研究中，多采用二进制编码来为分段开关和联络开关的状态。随着配电网规模的增大，这会产生维数泛滥，且其中的大量不满足辐射状约束的不可行解会降低算法的寻优效率。因此，本文采用实数编码并利用环路开关数组来限制寻优的解空间，限制了搜索空间，提高了搜索效率。生成环路开关数组需要注意的是与电源点相连的开关以及不处在闭合回路中的开关不能断开。具体实施步骤如下：1）对环路进行编码，每个联络开关对应一个环路，即有多少联络开关就相应的有多少环路开关集合，每一个开关集合将作为对应变量的可行空间；2）从编号为 1 的环路任意节点顺时针搜索，直至回到原点，中间所搜索到的开关作为环路1 集合的元素；3）对所有环路重复步骤
13、2），得到 n 个环路开关组集合，n 为配电网络对应的环路个数。以图 1 所示的 IEEE33 节点配电网为例，生成第4期彭博，等：基于改进烟花算法的配电网优化重构85的环路开关数组为L1=2,3,4,5,6,7,33,20,19,18L2=2,3,4,5,6,7,8,9,10,11,35,21,20,19,18L3=3,4,5,25,26,27,28,37,24,23,22L4=6,7,8,9,10,11,12,13,14,15,16,17,36,31,30,32,29,28,27,26,25L5=9,10,11,12,13,14,34101511217=242 550235在基于环路开关数
14、组确定的解空间中，最大的待选解个数为个，而用常规方法编码的待选解个数为个。可见，用环路开关数组限制解空间，能极大地提高算法搜索效率。3改进烟花算法3.1烟花算法原理烟花算法是一种模拟烟花爆炸过程的寻优算法，它通过烟花在爆炸过程中产生的火花实现对特定点的局部空间搜索12。当需要找到一个满足特定条件的点时，可以在潜在的空间不断地燃放“烟花”。这个特定的潜在空间相当于其他优化方法的可行解空间。通过不停地爆炸，直到所产生的“火花”相当接近于目标点。FWA 的具体流程如图 2 所示。烟花种群初始化通过爆炸和高斯变异产生火花对超过边界的火花采用映射规则使用选择策略选择下一代种群满足精度或者最大迭代次数开始
15、结束YN图2烟花算法框图1）初始化烟花种群进行参数的设定和烟花种群的初始化。参数主要包括初始种群数 P、变量维数 D、变异火花数 M、爆炸数目限制因子 a 和 b、爆炸半径调节常数 A、爆炸火花数调节常数 S 和最大迭代次数 N。2）产生爆炸火花xiSiAi对每一个烟花个体进行爆炸操作，设其爆炸数目为、爆炸半径为，对应的爆炸计算公式为Si=Symax f(xi)+Ni=1(ymax f(xi)+Ai=Af(xi)ymin+Ni=1(f(xi)ymin)+(1)yminymaxf(xi)i式中：为当前烟花种群中最小的适应度，为当前烟花种群中最大的适应度，为第 个烟花的适应度，为一个机器最小量值，
16、用来防止分母为零的情况。为了避免个体产生的爆炸火花数过多或过少，因此还需对其爆炸数目进行限制：Si=round(aS),Si aS,a b 1round(Si),其他round式中为取整函数。3）生成变异火花为保证种群的多样性，需要额外产生变异火花。在种群中随机选择 m 个烟花和 z 个维度进行高斯变异操作，具体为xki=xkiGxkixikG式中：为烟花维度为 的变量，为方差均值均为 1 的高斯分布。4）映射规则经爆炸操作和高斯变异操作，产生的新火花可能会在某个维度超出可行域限制范围。用映射规则来对越界的火花位置修正，将其映射到可行123456789101112131415161718192
17、0212223242526272829303132 33S1S2S3S4S5S6S7S8S9S10 S11S12S13 S14 S15 S16 S17S18S19 S20S21S35S33S22S23S24S37S25S26S27S28 S29 S30S31S32S36S3413245图1IEEE33 节点配电网86应用科技第50卷域中的新位置，这个新位置为xki=XLB+?xki?mod(XUBXLB)xkixikXUBXLB式中：为烟花维度 的值，和分别为可行域的上、下边界，mod 为模运算。5）选择策略xiP(xi)对当代烟花个体、爆炸火花和变异火花都进行适应度计算，选择最佳适应度个体作
18、为最优个体，其将被选择进入下一代种群中。通过距离选择剩余 N-1 个个体。对于，被选中的概率为P(xi)=R(xi)Pj=1R(xi)式中R(xi)=Pj=1d(xi,xj)=Pj=1xixkR(xi)xid(xi,xj)xixj其中：为与其他个体的欧氏距离，为和的欧式距离。3.2用于配电网重构的改进烟花算法本文针对配电网重构优化算法提出了改进：用搜索能力更强的柯西变异代替传统的高斯变异；利用迭代次数和烟花适应度信息自适应地改变爆炸半径，以更好地平衡全局搜索和局部搜索能力；在映射规则上采用随机映射规则；在选择策略上用反向学习增加选择种群的多样性进一步加强搜索能力。具体实现如下。1）爆炸半径的改
19、进由式（1）可知，在适应度值均较小的初期，若产生了一个局部最优解，该烟花会被认为优质烟花，其爆炸半径等于或接近于零，但是爆炸产生的火花数较多，会占据着大量的资源，并最终容易陷入局部收敛。为此，对烟花算法爆炸半径作如下改变：在爆炸初期时，要使其爆炸半径尽可能的大，在爆炸后期因逐渐趋于收敛要注意缩减爆炸半径，改进后的爆炸半径为At,i=(TmaxtTmax)Af(xi)fmax(x)At,iif(xi)ifmax(x)ATmaxt式中：为第 t 代第 个烟花的爆炸半径，为第个烟花的适应度值，为当前种群的最大适应度值，是常数，为最大迭代次数，为当前迭代次数。2）变异规则的改进柯西分布图形的两翼宽于高
20、斯分布，且其在零点附近的波峰低于标准高斯分布、两边趋向于零的速度比高斯分布要慢，由此可知柯西变异的扰动能力比高斯变异强13，变异范围更广，较容易跳出局部最优，因此采用柯西变异来替代传统高斯变异，其计算公式为xki=xkiCayuchy(0,1),rand(0,1)pxki,其他Cayuchy(0,1)p式中：是标准柯西分布函数，为随机变异概率。3）映射规则的改进对于爆炸后和变异后产生的子代火花不一定都分布在可行域中，传统映射规则容易将解映射到搜索空间原点附近位置，因此提出与重构相适应的映射规则。当第 k 维上的位置越限时，采用映射到可行域的一个随机位置：xki=xLB+round(rand(x
21、UBxLB)xkixixUBxLBrand式中：为在 k 维上的变量，和为优化空间的下边界和上边界，为取（0，1）的随机数。4）选择策略的改进为了扩大种群的多样性，运用精英反向学习选择策略，其内容如下：从候选集 1（烟花、爆炸火花、变异火花）中选取适应度最好的个体作为下一代烟花的“精英”；然后对其生成反向种群候选集 2，增加其种群搜索范围以及选择下一代种群的多样性；然后从 2 个候选集中采用轮盘赌的形式选取下一代个体，直到达到原来的种群规模。通过对没选择的个体生成反向候选集，一般动态反向14学习定义如下：xi,j=k(daj+dbj)xi,jdajdbjk 0,1式中：、分别为当前种群搜索空间
22、中第 j 维度上的最小值和最大值，d 反映距离；。改进后的烟花算法重构流程如图 3 所示。开始采用实数编码生成初始种群和初始化参数计算自适应爆炸半径和爆炸火花数目，并通过位移操作生成爆炸火花通过柯西变异产生变异火花应用精英反向学习选择策略选择下一代烟花种群达到最大迭代次数迭代次数+1采用随机映射规则获取网络参数，建立目标函数和约束条件，构成多目标函数输出最优个体结束火花越界火花越界NY图3改进烟花算法重构流程第4期彭博，等：基于改进烟花算法的配电网优化重构874算例分析采用图 1 所示的 IEEE33 节点配电网15来验证本文提出的算法。该配电网系统共有 37 条线路和 1 个电源节点，其中
23、S33、S34、S35、S36、S37为线路配置联络开关。经多轮测试，算法参数设置如下：种群数 P=10，变量维数 D=5，变异火花数M=5，爆炸数目下限限制因子 a=0.04，爆炸数目上限限制因子 b=0.8，爆炸半径调节常数 A=10，爆炸火花数调节常数为 S=50，最大迭代次数 N=30。重构前后的统计结果见表 1，图 4 给出了 33 个节点在重构前后的电压标幺值。由表 1 可知，重构后的网络有功损耗由重构前的 202.667kW 减小到144.537kW，减小了 28.7%，降低了网络损失，提高了网络运行的经济性；配电网的最低节点电压（p.u.）从重构前的 0.913 增加到 0.9
24、34，按照配电网 10%的电压最大允许偏差来看，明显提高了配电系统的供电质量和安全边际。表1IEEE33 节点配电网重构前后对比结果配电网结构断开开关编号 网络有功损耗/kW最低节点电压(p.u.)重构前33、34、35、36、37202.6670.913重构后7、11、34、36、37144.5370.9341.000.990.980.970.960.950.940.930.920.91电压幅值05101520253035节点重构前重构后图433 个节点的电压分布情况由图 4 可见重构后各节点电压比重构前得到整体改善，这表明本文算法得到的重构方案能够提高系统的稳定性。为体现本文所提方法的性能
25、，引入文献 16提出的双重混合粒子群优化算法（doublehybridparticleswarmoptimizationalgorithm，DHPSO）做对比，算例和其参数均按其建议设置，进行 30 轮独立测试，统计结果如表 2 所示。表23 种算法的收敛性能对比方法最优值最劣值平均值获最优值次数FWA40.13444.91643.5352DHPSO40.13442.96241.02318本文算法40.13440.50940.17226表 2 中，最优值和最劣值分别是 30 轮独立测试中各算法获得最佳优化解和最劣优化解，平均值则是 30 轮独立测试结果的均值。通过表 2 中3 种算法的结果对比
26、可知，3 种算法均能搜寻到最优方案；但在最劣目标值与平均目标值方面，相比于另外 2 种算法，本文算法能够保证在未搜寻到最佳重构方案时尽可能地接近最优。在 30 次独立计算中获得的最优值次数方面，本文算法获得最佳解的概率为 86.7%，高于原始 FWA 的6.7%和 DHPSO 算法的 63.3%，这表明本文算法有助于提高全局搜索能力，增加最优获取比率。图 5 为目标值随迭代次数的变化曲线。从图 5 中可见，FWA 算法在第 4 次迭代中陷入局部最优，在第 24 次迭代左右跳出陷阱达到全局最优；本文算法在第 46 次迭代中到达局部最优值，仅在 2 次左右迭代就跳出局部最优，在第 9 次目标函数值
27、获得最优值 40.134，显著优于 FWA 算法。这表明本文算法能够快速跳出局部最优，有效防止早熟，易获得最优值。44.043.543.042.542.041.541.040.540.0目标值051015202530迭代次数FWA本文算法图5算法迭代对比曲线5结论本文介绍了配电网重构的一种新方法。重构目标兼顾最小网损和最少开关次数，提出环路开关数组概念，并利用环路数组确定算法的解空间，基于烟花算法原理并根据其不足提出了一种改进型的 FWA。通过 IEEE33 节点配电网标准算例测试并与其他算法比对。结果表明本文提出的重构优化方法在寻优速度和寻优质量方面均获得了提高，得出结论如下。1）本文提出的
28、配电网重构模型和解空间压缩方法，能够有效解决配电网优化重构问题。2）本 文 提 出 的 改 进 算 法 与 传 统 FWA、DHPSO 算法相比，具有更加优秀的寻优质量和更快的寻优速度，有助于提高在限定时间内获取的“最优解”质量，助力重构决策的质量。88应用科技第50卷参考文献：MERLINA,BACKH.Searchforaminimal-lossoperatingspanningtreeconfigurationinanurbanpowerdistributionsystemC/TheFifthPowerSystemsComputerConference.Cambridge:s.n.197
29、5:118.1许向雷.含分布式电源的配电网重构研究 D.沈阳:东北大学,2017.2杨明,翟鹤峰,马嘉翼,等.计及分布式电源发电不平衡度约束的三相不对称配电网动态重构 J.中国电机工程学报,2019,39(12):34863499.3SHIRMOHAMMADID,HONGHW.Reconfigurationofelectric distribution networks for resistive line lossesreductionJ.IEEEtransactionsonpowerdelivery,1989,4(2):14921498.4刘蔚,韩祯祥.基于最优流法和遗传算法的配电网重构J
30、.电网技术,2004,28(19):2933.5CIVANLARS,GRAINGERJJ,YINH,etal.DistributionfeederreconfigurationforlossreductionJ.IEEEtransactionsonpowerdelivery,1988,3(3):12171223.6张栋,张刘春,傅正财.配电网络重构的快速支路交换算法 J.电网技术,2005,29(9):8285.7MOHAMEDIMRANA,KOWSALYAM.Anewpowersystemreconfigurationschemeforpowerlossminimizationandvolt
31、ageprofileenhancementusingFireworks8AlgorithmJ.International journal of electrical power&energysystems,2014,62:312322.孙旭日,李延真.配电网在风险状态下的重构优化方法 J.应用科技,2020,47(4):100105.9陈春,汪沨,刘蓓,等.基于基本环矩阵与改进和声搜索算法的配电网重构 J.电力系统自动化,2014,38(6):5560.10徐泽,杨伟,张文强,等.基于连锁环网与改进离散粒子群算法的多目标配电网重构 J.电力系统保护与控制,2021,49(6):114123.1
32、1TAN Ying,Zhu Yuanchun.Fireworks algorithm foroptimizationC/Proceedings of First InternationalConferenceonAdvancesinSwarmIntelligence.Beijing:Springer-Verlag,2010:355364.12曾敏,赵志刚,李智梅.带柯西变异的自学习改进烟花算法 J.小型微型计算机系统,2020,41(2):264270.13毛清华,张强.融合柯西变异和反向学习的改进麻雀算法 J.计算机科学与探索,2021,15(6):11551164.14刘路宁,彭春华,温泽
33、之,等.基于配电网动态重构的分布式光伏消纳策略 J.电力自动化设备,2019,39(12):5662.15马草原,孙展展,尹志超,等.基于双重混合粒子群算法的 配 电 网 重 构 J.电 工 技 术 学 报,2016,31(11):120128.16本文引用格式：彭博,邹乐,李昕,等.基于改进烟花算法的配电网优化重构 J.应用科技,2023,50(4):8489.PENGBo,ZOULe,LIXin,etal.OptimalreconfigurationofdistributionnetworkbasedonimprovedfireworksalgorithmJ.Appliedscienceandtechnology,2023,50(4):8489.第4期彭博，等：基于改进烟花算法的配电网优化重构89