一种综合无功优化与网络重构的配网运行优化方法

文档序号:10572097阅读:408来源:国知局
一种综合无功优化与网络重构的配网运行优化方法
【专利摘要】本发明公布了一种综合无功优化与网络重构的配网运行优化方法,该方法充分考虑了配网优化运行中无功优化与网络重构两种技术的协调与配合。本发明首先定义了以年综合费用作为目标函数的优化模型,接着给出新的编码策略以及解的辨识方法,最后,采用混合蛙跳算法求得最佳综合优化方案。本发明能充分发掘无功优化与网络重构两项技术相互配合存在的潜在经济价值,具有一定的实用参考价值。
【专利说明】
一种综合无功优化与网络重构的配网运行优化方法
技术领域
[0001] 本发明属于电力系统运行分析和控制技术领域,涉及一种综合无功优化与网络重 构的配网运行优化方法。 技术背景
[0002] 配网重构和配网无功优化作为配网优化运行的两项重要技术手段,是配电自动化 的关键技术,其在保障电能质量、降低网络损耗等方面有着重要作用。配网重构通过改变网 络开关的闭合来获得最佳优化目标值下的网络拓扑结构;配网无功优化一般通过定点投切 无功补偿设备以实现有功损耗最小化且保证较高的电压水平。实质上,配网重构是非线性 组合优化问题,配网无功优化是非线性整数规划问题,二者的综合优化使得问题的求解更 加复杂,针对此问题,有学者分别提出采用先重构后补偿交替迭代的方法、分别以重构和无 功优化为主进行优化,后者采用了智能算法求解,相比前一种交替迭代法,提高了计算精 度,它们并不是真正意义上同步进行重构与无功优化,此外,现有研究均以网损作为优化目 标,并没有考虑无功补偿的经济性,综合优化的潜在价值没有得到充分研究与挖掘。因而本 发明公布了一种综合无功优化与网络重构的配网运行优化方法。

【发明内容】

[0003] 发明目的:本发明的目的在于针对现有技术的不足,提出了一种综合无功优化与 网络重构的配网运行优化方法,建立以年综合费用为目标函数的综合优化模型,并在求解 过程中,简化网络,提高解的辨识度,最后通过蛙跳算法获得综合优化方案,为运行人员提 供决策支持。
[0004] 技术方案:本发明提供一种综合无功优化与网络重构的配网运行优化方法,包括 以下步骤:
[0005] 步骤1:建立以年综合费用作为目标函数的综合优化模型;
[0006] 步骤2:给出新的编码策略以及解的识别方法;
[0007] 步骤3:采用混合蛙跳算法求得最佳综合优化方案。
[0008] 进一步,所述步骤1中综合优化模型为:
[0009]
[0010] 式中:CT为年综合费用,λ为电价;Tmax为年最大负荷损耗小时数;ki为补偿设备的 年维护费用率;k 2为投资回收系数;Wi为第i个节点的无功补偿量,Cl为无功补偿的价格;C2 为单个补偿点的安装费用,k3为无功补偿点个数,Pi〇ss为网络的有功损耗,其大小受网络拓 扑结构与无功补偿量的影响,其值为每条线路有功损耗的总和,表达式如下:
[0011;
[0012] ΑΨ,Nb衣不文?数;Kk衣示支路k的电阻;Hk表示支路闭合状态,1表示闭合,0表示 打开;Pk表不支路k的有功功率;Qk表不支路k的无功功率;Vk表不支路k的末端电压;
[0013] 潮流计算中,除保持网络功率平衡外,还需要考虑以下约束条件:
[0014] Vmin^Vj^Vmax
[0015] \Sk\<Sr
[0016] 〇^ffi^ffi;max
[0017] 式中,Vmin,vmax分别表示配网正常运行时节点电压的上下限;S k表示支路k的载流 量,Smax表示线路k的最大载流量;Wi,max表示第i个补偿点补偿容量上限。
[0018] 进一步,所述步骤2包括以下步骤:
[0019] 步骤201:网络简化及编码:
[0020] 配网重构通过改变网络开关的闭合来获得最佳优化目标值下的网络拓扑结构,为 了提高求解的搜索效率,以33节点配网为例,对网络做如下简化:根节点并入1节点;由于节 点2和节点20之间只允许断开一个开关,因此可以将这两个节点间的所有支路视为一个支 路组,以此类推,简化网络,简化后的网络由8个节点,12个支路组组成。基于此,本文采用基 于独立环路的编码策略可以有效地减少变量维数,能够提高求解的效率;
[0021] 在混合蛙跳算法中,每只蛙相当于一个优化方案,第i只蛙编码为Ziinxi X2,......,Χη?} ; {yi,y2,......,yn2}],其中,{xi,X2,··· · · ,Χη?}表示重构断开的支路在支路组中 的编号。{yi,y2,......,yn2}表示选定节点的无功补偿量,nl、n2分别表示支路组和无功补偿 点的个数;
[0022] 步骤202:解的辨识:
[0023] 利用图论中连通度理论对不可行解进行辨识,改进使网络呈孤岛或环网的解,使 得解满足网络辐射状约束条件,等效后的配电网络可以用一个图表示,其各节点间的连接 关系可以用邻接矩阵A表示:
[0024]
[0025] 其中,Μ为等效网络的节点个数。若节点i与j相连,则ai,j为1,否则为0,构造图的 laplacian 矩阵 B:
[0026] B = diag(sum(A) )-A
[0027] 当rank(B)=M_l等式成立,并且当断开的支路等于网络的独立环路时,网络的拓 扑结构则满足辐射状要求。
[0028] 进一步,所述步骤3中蛙跳算法求解过程包括以下两个步骤:
[0029]步骤301:全局搜索
[0030] 步骤a:初始化参数,包括:蛙群的数量F;族群的数量m;族群中青蛙的数量η;最大 允许跳动步长Smax;全局最优解Ρ ζ;局部最优解Pb;局部最差解Pw;全局迭代进化次数Ng,局部 迭代进化次数Λ,各补偿点无功补偿上限W 1>max;
[0031] 步骤b:随机生成初始蛙群,由目标函数计算每个蛙的评价值;
[0032] 步骤c :按照评价值大小进行升序排序,记录下最优解Pz,并且将蛙群按以下方式 分成族群:第1只蛙放入第1个族群,第2只蛙放入第2个族群,第m只蛙放入第m个族群,第m+1 只蛙放入第1个族群,以此类推,直至所有蛙被放入入指定位置;
[0033] 步骤d:按照下式对每个族群进行进化操作
[0034] SL = ceil(Rand() X (Pw-Pb))
[0035] NewPw= Pw+SL , -SLmin^ SL^ SLmax
[0036] 式中,ceil表示取整,rand()表示产生0~1的随机数,SL表示蛙跳的步长,SLmax, Lmin为娃跳的步长限制,NewPw表不更新后的Pw;
[0037] 步骤e:所有族群更新完毕后,计算蛙群中所有蛙的评价值;
[0038] 步骤f:判断是否满足停止条件。如果满足则停止搜索,否则转到步骤c;
[0039]步骤302:局部搜索,是指对上述步骤d展开的具体描述,过程如下:
[0040] 步骤di:设Ιμ= Ιν = 0,Im表示族群进化的计数器,In表示局部进化计数器;
[0041 ] 步骤d2:选出当前族群的Pb和Pw,Im加1;
[0042]步骤d3:IN加 1;
[0043]步骤d4:根据步骤d中两个式子改进族群中最差蛙;
[0044] 步骤d5:如果上步改进了最差蛙,则用该新蛙取代最差蛙,否则用Pz替代式步骤d中 的Pb,重新进化;
[0045] 步骤d6:如果上步仍没有改进最差蛙,则随机产生一个可行解来代替最差蛙;
[0046] 步骤d7:如果In小于局部进化次数Ln,则转入步骤d3;
[0047] 步骤d8:如果Im小于族群数m,则转入步骤d2,否则进入全局搜索的步骤d 5。
[0048] 工作原理:本发明首先定义了以年综合费用作为目标函数的优化模型,接着给出 新的编码策略以及解的辨识方法,最后,采用混合蛙跳算法求得最佳综合优化方案。本发明 能充分发掘无功优化与网络重构两项技术相互配合存在的潜在经济价值,具有一定的实用 参考价值。
[0049]有益效果:与现有综合优化技术相比,本发明具有如下优点和技术效果:
[0050] 1)该方法公布的模型不仅能够进一步降低网损、提升节点电压,还充分考虑了无 功补偿的经济性,使得配网综合优化更加符合经济性原则,即充分挖掘了综合优化的经济 价值,也为综合优化调度提供了新的评价指标。
[0051] 2)除了提出新的优化目标,本发明采用了易于理解、易于编程实现及寻优能力强 的混合蛙跳算法,并结合基于独立支路法的网络简化策略,能够快速搜索到最优解,保证了 模型求解效率。
【附图说明】
[0052]图1为本发明配电网络简化图;
[0053]图2为本发明总流程图;
[0054]图3为本发明中蛙跳算法流程图;
[0055]图4为本发明节点电压分布图;
[0056]图5为本发明无功优化过程目标值收敛曲线;
[0057]图6为本发明重构过程目标值收敛曲线;
[0058]图7为本发明综合优化过程目标值收敛曲线。
【具体实施方式】:
[0059] 以下结合附图和实例对本发明的实施作进一步说明,但本发明的实施和包含不限 于此。
[0060] 一种综合无功优化与网络重构的配网运行优化方法,包括以下步骤:
[0061 ]步骤1:建立以年综合费用作为目标函数的综合优化模型;
[0062] 步骤2:给出新的编码策略以及解的识别方法;
[0063] 步骤3:采用混合蛙跳算法求得最佳综合优化方案。
[0064] 进一步,所述步骤1中综合优化模型为:
[0065]
[0066] 式中:CT为年综合费用,λ为电价;!'max为年最大负荷损耗小时数;lu为补偿设备的 年维护费用率;k 2为投资回收系数;Wi为第i个节点的无功补偿量,Ci为无功补偿的价格;C2 为单个补偿点的安装费用,k3为无功补偿点个数,Pi〇ss为网络的有功损耗,其大小受网络拓 扑结构与无功补偿量的影响,其值为每条线路有功损耗的总和,表达式如下:
[0067]
[0068] 式中,Nb表示支路数;Rk表示支路k的电阻;Hk表示支路闭合状态,1表示闭合,0表示 打开;Pk表不支路k的有功功率;Qk表不支路k的无功功率;Vk表不支路k的末端电压;
[0069] 潮流计算中,除保持网络功率平衡外,还需要考虑以下约束条件:
[0070] Vmin^Vj^Vmax
[0071] \sk\<sr
[0072] 〇^ffi^ffi;max
[0073] 式中,Vmin,Vmax分别表示配网正常运行时节点电压的上下限;S k表示支路k的载流 量,SmaX表示线路k的最大载流量;Wi,max表示第i个补偿点补偿容量上限。
[0074] 进一步,所述步骤2包括以下步骤:
[0075] 步骤201:网络简化及编码:
[0076] 配网重构通过改变网络开关的闭合来获得最佳优化目标值下的网络拓扑结构,为 了提高求解的搜索效率,以33节点配网为例,对网络做如下简化:根节点并入1节点;由于节 点2和节点20之间只允许断开一个开关,因此可以将这两个节点间的所有支路视为一个支 路组,以此类推,简化网络,简化后的网络由8个节点,12个支路组组成。基于此,本文采用基 于独立环路的编码策略可以有效地减少变量维数,能够提高求解的效率;
[0077] 在混合蛙跳算法中,每只蛙相当于一个优化方案,第i只蛙编码为Z i =[{ X i, X2,......,Χη?} ; {yi,y2,......,yn2}],其中,{xi,X2,··· · · ,Χη?}表示重构断开的支路在支路组中 的编号。{yi,y2,......,yn2}表示选定节点的无功补偿量,nl、n2分别表示支路组和无功补偿 点的个数;
[0078] 步骤202:解的辨识:
[0079] 利用图论中连通度理论对不可行解进行辨识,改进使网络呈孤岛或环网的解,使 得解满足网络辐射状约束条件,等效后的配电网络可以用一个图表示,其各节点间的连接 关系可以用邻接矩阵A表示:
[0080]
[0081] 其中,Μ为等效网络的节点个数。若节点i与j相连,则ai,j为1,否则为0,构造图的 laplacian 矩阵 B:
[0082] B = diag(sum(A) )-A
[0083] 当rank(B)=M_l等式成立,并且当断开的支路等于网络的独立环路时,网络的拓 扑结构则满足辐射状要求。
[0084] 进一步,所述步骤3中蛙跳算法求解过程包括以下两个步骤:
[0085]步骤301:全局搜索
[0086] 步骤a:初始化参数,包括:蛙群的数量F;族群的数量m;族群中青蛙的数量η;最大 允许跳动步长Smax;全局最优解Ρ ζ;局部最优解Pb;局部最差解Pw;全局迭代进化次数Ng,局部 迭代进化次数Λ,各补偿点无功补偿上限W 1>max;
[0087] 步骤b:随机生成初始蛙群,由目标函数计算每个蛙的评价值;
[0088] 步骤c :按照评价值大小进行升序排序,记录下最优解Pz,并且将蛙群按以下方式 分成族群:第1只蛙放入第1个族群,第2只蛙放入第2个族群,第m只蛙放入第m个族群,第m+1 只蛙放入第1个族群,以此类推,直至所有蛙被放入入指定位置;
[0089] 步骤d:按照下式对每个族群进行进化操作
[0090] SL = ceil(Rand() X (Pw-Pb))
[0091 ] NewPw= Pw+SL , -SLmin^ SL^ SLmax
[0092] 式中,ceil表示取整,rand()表示产生0~1的随机数,SL表示蛙跳的步长,SLmax, Lmin为娃跳的步长限制,NewPw表不更新后的Pw;
[0093] 步骤e:所有族群更新完毕后,计算蛙群中所有蛙的评价值;
[0094] 步骤f:判断是否满足停止条件。如果满足则停止搜索,否则转到步骤c;
[0095]步骤302:局部搜索,是指对上述步骤d展开的具体描述,过程如下:
[0096] 步骤di:设Ιμ= Ιν = 0,Im表示族群进化的计数器,In表示局部进化计数器;
[0097] 步骤d2:选出当前族群的Pb和Pw,Im加 1;
[0098] 步骤d3:1^1;
[0099] 步骤d4:根据步骤d中两个式子改进族群中最差蛙;
[0100] 步骤d5:如果上步改进了最差蛙,则用该新蛙取代最差蛙,否则用Pz替代式步骤d中 的Pb,重新进化;
[0101]步骤d6:如果上步仍没有改进最差蛙,则随机产生一个可行解来代替最差蛙;
[0102] 步骤d7:如果In小于局部进化次数Ln,则转入步骤d3;
[0103] 步骤d8:如果Im小于族群数m,则转入步骤d2,否则进入全局搜索的步骤d5。
[0104] 工作原理:本发明首先定义了以年综合费用作为目标函数的优化模型,接着给出 新的编码策略以及解的辨识方法,最后,采用混合蛙跳算法求得最佳综合优化方案。本发明 能充分发掘无功优化与网络重构两项技术相互配合存在的潜在经济价值,具有一定的实用 参考价值。
[0105] 实施例
[0106] 本发明以IEEE33节点系统作为说明案例,总流程图见附图2,该系统电压等级为 12.66kV,有功总负荷3715kW,无功总负荷2300kvar。针对网络重构、无功优化、综合优化三 种优化方案进行对比分析。
[0107] 蛙跳算法参数设置如下:蛙群大小为80,族群数为20,全局进化次数为50,局部进 化次数为3,流程图见附图3。无功补偿点选择无功负荷最重的23、24、29节点,一方面符合无 功就地补偿原则。无功补偿量以系统总无功负荷的1.2倍作为无功补偿的上限,以lOkvar作 为最小搜索步长,设71£[0,60], 72£[0,60],73£[0,120]。其它参数4 = 0.5元/妍·!!,!^ = 5000h,ki = 0 · 13,k2 = 0 · 1,Ci = 60元/kvar,C2 = 5000元/节点。
[0108] 三种方案优化结果对比如下表1,各优化方案下系统各节点电压分布如附图4,求 解迭代过程目标函数值收敛情况如附图5至附图7。
[0109]表1本发明三种方案优化结果对比表
[0110]
[0111]
[0112]表1中各方案优化结果对比表明,综合优化方案相比单纯的重构与无功优化,能够 进一步提高电压质量:最低节点电压0.9575P.U,相比单独重构与无功优化的结果分别提高 了2.1 %、3.5% ;能够进一步减少网络有功损耗:有功损耗102.09kW,相比单独重构与无功 优化的结果降低了 26.7%;进一步降低年综合费用:年综合费用为27.68万,相比单独重构 与无功优化的结果分别降低了20.5%、25.8%。此外,综合优化方案中的无功补偿相比单独 的无功补偿,补偿量降低了 15.0%。验证了本发明所提综合优化模型的合理性与有效性。
[0113] 结合表1与图4可知,单独定点无功优化与重构都能够有效降低网络有功损耗,但 在提升系统各节点电压水平方面,重构效果较无功优化更好。综合优化方案相比重构,系统 各节点电压整体有进一步提升,从而进一步验证了表1的结论。
[0114] 由图5至图7可知,SFLA在求解三个优化模型过程中,具有良好的收敛效果,求解时 间均在25秒以内,能够_效求解本文所提t旲型。
[0115] 综上,案例说明了本发明的正确性与实用性。
【主权项】
1. 一种综合无功优化与网络重构的配网运行优化方法,其特征在于:包括以下步骤: 步骤1:建立以年综合费用作为目标函数的综合优化模型; 步骤2:给出新的编码策略以及解的辨识方法; 步骤3:采用混合蛙跳算法求得最佳综合优化方案。2. 根据权利要求1所述的综合无功优化与网络重构的配网运行优化方法,其特征在于: 所述步骤1中综合优化模型为:式中:CT为年综合费用,λ为电价;!'max为年最大负荷损耗小时数;Iu为补偿设备的年维护 费用率;k2为投资回收系数;Wi为第i个节点的无功补偿量,Ci为无功补偿的价格;C2为单个 补偿点的安装费用,k 3为无功补偿点个数;Plciss为网络的有功损耗,其大小受网络拓扑结构 与无功补偿量的影响,其值为每条线路有功损耗的总和,表达式如下:式中,Nb表示支路数;Rk表示支路k的电阻;Hk表示支路闭合状态,1表示闭合,O表示打 开;Pk表不支路k的有功功率;Qk表不支路k的无功功率;Vk表不支路k的末端电压; 潮流计算中,除保持网络功率平衡外,还需要考虑以下约束条件: Vmin^Vj ^VmaxO^ffi^ffi.max 式中,V?in,V-分别表示配网正常运行时节点电压的上下限;&表示支路k的载流量,S-表示线路k的最大载流量;Wi,max表示第i个补偿点补偿容量上限。3. 根据权利要求1所述的种综合无功优化与网络重构的配网运行优化方法,其特征在 于:所述步骤2包括以下步骤: 步骤201:网络简化及编码: 配网重构通过改变网络开关的闭合来获得最佳优化目标值下的网络拓扑结构,为了提 高求解的搜索效率,以33节点配网为例,对网络做如下简化:根节点并入1节点;由于节点2 和节点20之间只允许断开一个开关,因此可以将这两个节点间的所有支路视为一个支路 组,以此类推,简化网络,简化后的网络由8个节点,12个支路组组成;基于此,本发明采用基 于独立环路的编码策略可以有效地减少变量维数,能够提高求解的效率; 在混合娃跳算法中,每只娃相当于一个优化方案,第i只娃编码为Zi = [ {χι,Χ2,......, Xnl} ; {yi,Y2,......,yn2}],其中,{xi,X2,··· · · ,Xnl}表示重构断开的支路在支路组中的编号。 {yi,y2,......,yn2}表示选定节点的无功补偿量,nl、n2分别表示支路组和无功补偿点的个 数; 步骤202:解的辨识: 利用图论中连通度理论对不可行解进行辨识,改进使网络呈孤岛或环网的解,使得解 满足网络辐射状约束条件,等效后的配电网络可以用一个图表示,其各节点间的连接关系 可以用邻接矩阵A表示:其中,M为等效网络的节点个数;若节点i与j相连,则ai, j为I,否则为O ;构造图的 Iaplacian 矩阵 B: B = diag(sum(A) )_A 当rank(B) =M-I等式成立,并且当断开的支路等于网络的独立环路时,网络的拓扑结 构则满足辐射状要求。4.根据权利要求1所述的一种综合无功优化与网络重构的配网运行优化方法,其特征 在于:所述步骤3蛙跳算法求解过程包括以下两个步骤: 步骤301:全局搜索,包括以下分步骤: 步骤a:初始化参数,包括蛙群的数量F;族群的数量m;族群中青蛙的数量η;最大允许跳 动步长SLmax;全局最优解Pz;局部最优解Pb;局部最差解P w;全局迭代进化次数Ng,局部迭代 进化次数阶,各补偿点无功补偿上限W1>max; 步骤b:随机生成初始蛙群,由目标函数计算每个蛙的评价值; 步骤c:按照评价值大小进行升序排序,记录下最优解Pz,并且将蛙群按以下方式分成族 群:第1只蛙放入第1个族群,第2只蛙放入第2个族群,第m只蛙放入第m个族群,第m+1只蛙放 入第1个族群,以此类推,直至所有蛙被放入指定位置; 步骤d:按照下式对每个族群进行进化操作 SL = ceil(Rand() X (Pw-Pb)) NewPw = Pw+SL , -SLmin^ SL^ SLmax 其中,ce i 1表示取整,rand ()表示产生0~1的随机数,SL表示蛙跳的步长,SLmax,SLmin为 娃跳的步长限制,NewPw表不更新后的Pw; 步骤e:所有族群更新完毕后,计算蛙群中所有蛙的评价值; 步骤f:判断是否满足停止条件,如果满足则停止搜索,否则转到步骤c; 步骤302:局部搜索,是指对上述步骤d展开的具体描述,分步骤如下: 步骤di:设Im= In=O,Im表示族群进化的计数器,In表示局部进化计数器; 步骤d2:选出当前族群的Pb和Pw,Im加1; 步骤d3: In加1; 步骤d4:根据步骤d中两个式子改进族群中最差蛙; 步骤d5:如果上步改进了最差蛙,则用该新蛙取代最差蛙,否则用Pz替代式步骤d中的 Pb,重新进化; 步骤d6:如果上步仍没有改进最差蛙,则随机产生一个可行解来代替最差蛙; 步骤d7:如果In小于局部进化次数Ln,则转入步骤d3; 步骤d8:如果Im小于族群数m,则转入步骤d2,否则进入全局搜索的步骤d5。
【文档编号】G06F17/50GK105932690SQ201610335137
【公开日】2016年9月7日
【申请日】2016年5月19日
【发明人】卫志农, 王薪苹, 孙国强, 李逸驰, 臧海祥
【申请人】河海大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1