上海交通大学学报, 2022, 56(10): 1294-1307 doi: 10.16183/j.cnki.jsjtu.2021.371

新型电力系统与综合能源

基于孔雀优化算法的配电网储能系统双层多目标优化配置

杨博1, 王俊婷1, 俞磊1, 曹璞璘,1, 束洪春1, 余涛2,3

1.昆明理工大学 电力工程学院,昆明 650500

2.华南理工大学 电力学院,广州 510640

3.广东省电网智能量测与先进计量企业重点实验室,广州 510640

Peafowl Optimization Algorithm Based Bi-Level Multi-Objective Optimal Allocation of Energy Storage Systems in Distribution Network

YANG Bo1, WANG Junting1, YU Lei1, CAO Pulin,1, SHU Hongchun1, YU Tao2,3

1. Faculty of Electric Power Engineering, Kunming University of Science and Technology, Kunming 650500, China

2. College of Electric Power, South China University of Technology, Guangzhou 510640, China

3. Guangdong Provincial Key Laboratory of Intelligent Measurement and Advanced Metering of Power Grid, Guangzhou 510640, China

通讯作者: 曹璞璘,男,副教授;E-mail:pulincao_kust@sina.com.

责任编辑: 孙伟

收稿日期: 2021-09-24  

基金资助: 国家自然科学基金(61963020)
国家自然科学基金-国家电网公司智能电网联合基金(U2066212)
云南省自然科学基金(202001AT070096)
云南省重大科技专项计划(202002AF080001)

Received: 2021-09-24  

作者简介 About authors

杨博(1988-),男,云南省昆明市人,教授,从事新能源发电/储能系统优化与控制,以及人工智能在智能电网中的应用研究.

摘要

考虑电池储能系统(BESSs)规划与运行之间的联系,建立兼顾经济性和技术性要求的BESSs多目标优化配置模型并进行双层架构,保证BESSs规划的有效性和运行的高效性.内层以BESSs运营收益最大为目标,提出孔雀优化算法求解BESSs充放电运行策略的最优解;外层以BESSs投资运维成本、配电网电压波动和负荷波动最小为目标,设计多目标孔雀优化算法求解选址定容规划方案的Pareto非支配解集.考虑配电网运行条件的不确定性,采用聚类算法获得典型场景集,并基于IEEE-33节点系统进行仿真.结果表明:所提算法实现了局部探索和全局搜索的平衡,有效获得高质量解;与传统多目标优化算法相比,其能够获得分布更广泛且均匀的Pareto前沿,实现BESSs投资效益最优,显著提升配电网电压质量和功率稳定性.

关键词: 电池储能系统; 双层优化配置模型; 孔雀优化算法; Pareto多目标优化

Abstract

Based on the relation between battery energy storage systems (BESSs) planning and operation, a multi-objective optimal allocation model that takes into account both economic and technical requirements is established, and a bi-level optimization structure is constructed to ensure effective planning and high-efficient operation of BESSs. In the inner layer, a peafowl optimization algorithm (POA) is employed to solve the BESSs charge-discharge operation strategy with the purpose of BESSs operation benefit maximization. In the outer layer, a multi-objective peafowl optimization algorithm (MOPOA) is devised to solve the Pareto solution set of BESSs siting and sizing scheme, which aims at minimizing BESSs cost, as well as voltage fluctuation and load fluctuation in distribution network. Furthermore, a typical scenario set is obtained via the clustering algorithm considering uncertain operating conditions. The simulation is performed based on the extended IEEE-33 bus system. The results show that the proposed algorithm achieves a trade-off between local search and global search, thus obtains a high-quality solution. It can obtain a more widely distributed and uniform Pareto front, which not only achieves the best investment benefit, but also improves voltage quality and power stability.

Keywords: battery energy storage systems (BESSs); bi-level optimal allocation model; peafowl optimization algorithm (POA); Pareto based multi-objective optimization

PDF (2519KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

杨博, 王俊婷, 俞磊, 曹璞璘, 束洪春, 余涛. 基于孔雀优化算法的配电网储能系统双层多目标优化配置[J]. 上海交通大学学报, 2022, 56(10): 1294-1307 doi:10.16183/j.cnki.jsjtu.2021.371

YANG Bo, WANG Junting, YU Lei, CAO Pulin, SHU Hongchun, YU Tao. Peafowl Optimization Algorithm Based Bi-Level Multi-Objective Optimal Allocation of Energy Storage Systems in Distribution Network[J]. Journal of Shanghai Jiaotong University, 2022, 56(10): 1294-1307 doi:10.16183/j.cnki.jsjtu.2021.371

电池储能系统(Battery Energy Storage Systems, BESSs)凭借其快速的功率调节和灵活的能量管理能力,可有效解决配电网弃风弃光、电压越限、潮流反向和调峰能力不足等问题[1-6].合理配置BESSs对配电网高效经济运行具有重要意义[7].

目前,国内外诸多学者已对BESSs优化配置模型展开研究.例如,文献[8-9]分别以配电网年净收益最大、年综合成本最小为优化目标,建立了BESSs单目标规划模型.然而,在进行BESSs规划时,需要同时考虑安全性、可靠性、经济性等相互冲突和影响的多个目标,单一目标下的优化难以满足实际工程需要.进而,不少学者利用线性加权等方式将BESSs规划的多个优化目标加权为单个目标进行优化求解[10-11],但未能客观地分配各目标权重.文献[12]提出了BESSs容量配置和优化布点的双层优化模型,上层以投资经济性确定容量,下层以负荷方差最小进行选址,即在上层优化的结果中再进行下层优化,但实际上未考虑两层优化模型之间的耦合关系.文献[13]提出了一种含峰值功率能力的源-网-荷-储双层协同规划模型,有效实现了各层目标函数和决策变量在上下层的优化传递.

此外,鉴于BESSs的优化配置是一个多变量、且含离散变量的复杂非线性多目标优化问题,多目标进化算法能够比传统的数学优化方法[14]更有效地获得全局最优解[15].文献[16]采用非支配解排序遗传算法(Non-Dominated Sorting Genetic Algorithm-II, NSGAII)求解BESSs多目标规划模型,但该算法交叉变异操作的随机性较强.文献[17]提出改进的多目标粒子群优化(Multi-Objective Particle Swarm Optimization, MOPSO)算法求解BESSs安装位置和容量,然而该算法惯性权重的取值缺乏指导,寻优性能稍差.

针对上述问题,本文建立BESSs双层多目标优化配置模型,并设计性能良好的优化算法求解BESSs优化配置方案:

(1) 建立基于Pareto的多目标优化模型,以实现BESSs投资成本、配电网电压波动和负荷波动3个目标的最佳权衡.相较于单目标优化[7-8]和基于线性加权的多目标优化[9-10],所提模型不仅综合考虑了BESSs投资经济性和配电网运行可靠性,而且能客观地分配各目标权重,避免因采取主观权重而导致对某一目标产生较大偏好.

(2) 对BESSs优化配置模型进行双层架构,将建立的多目标优化模型作为外层模型,以求解BESSs选址定容方案;内层模型则以最大化BESSs运营效益为目标求解其充放电运行策略.相较于单层优化和上下层独立循环的双层优化[11],所构建的双层模型考虑了BESSs规划和运行之间的联系,通过内外层的嵌套循环迭代,保证了优化配置方案的有效性.

(3) 提出一种孔雀优化算法(Peafowl Optimization Algorithm, POA),并基于Pareto理论设计多目标孔雀优化算法(Multi-Objective Peafowl Optimization Algorithm, MOPOA),分别求解内层和外层模型.在自适应搜索和逼近机制中设计可平衡局部探索和全局搜索的算子,使内层优化有效逼近高质量的最优解,外层优化获得分布广泛且均匀的Pareto非支配解集,从而避免了数学优化方法[13]和传统多目标进化算法[14-15]存在的易陷入局部最优等固有缺陷.

(4) 采用聚类算法对全年的负荷曲线和风、光出力曲线进行典型日聚类,基于其时序特性组合得到4种典型场景.在不同典型场景下分别进行仿真计算,得到各目标函数的全年总和,更好地反映了负荷和新能源在全年的不确定性.以扩展的IEEE-33节点配电系统为例进行仿真测试,验证了所提方法的有效性.

1 BESSs双层多目标优化配置模型

鉴于BESSs的后期运行与初期规划相互影响,并考虑BESSs选址定容方案与其充放电功率之间的耦合关系,建立BESSs双层多目标优化配置模型.一方面,将长时间尺度的BESSs规划问题放在外层模型,短时间尺度的BESSs运行问题放在内层模型,通过内外层优化的有效衔接,使得BESSs规划方案更合理,且更符合实际运行,有助于提高BESSs的利用率.另一方面,将BESSs充放电功率放于内层进行优化,使得外层模型的决策变量维度呈指数式缩减,可有效提高Pareto解集的分布性,并加快收敛速度.

1.1 BESSs最优选址定容规划

外层以BESSs年投资运维成本、配电网年电压波动和年负荷波动最小为目标,建立考虑BESSs投资成本效益和配电网运行可靠性的Pareto多目标优化模型:

minXFh(X,x), h=1,2,3s.t.g(X,x)0

式中:$F_{h}(\boldsymbol{X},\boldsymbol{x})$为外层模型的第h个目标函数;X和x分别为外层和内层模型的决策变量向量;$g(\boldsymbol{X}, \boldsymbol{x})$为外层模型需满足的约束条件.

1.1.1 目标函数

BESSs的年投资运维成本.投资方往往以成本效益为原则配置BESSs,投资效益主要来自于BESSs的低储高发套利.BESSs年投资运维成本为

1.1.1 目标函数

BESSs的年投资运维成本.投资方往往以成本效益为原则配置BESSs,投资效益主要来自于BESSs的低储高发套利.BESSs年投资运维成本为

F1=CTCC+COM-IS-P-Isub

式中:CTCCCOM分别为BESSs的等年值投资成本和年运行维护成本;IS-PIsub分别为每年的低储高发套利和发售电量的政府补贴.

对于CTCC,有

CTCC=[CinvNBESS+   n=1NBESS(aPBESS,n+bEBESS,n)]μCRFμCRF=r(1+r)y(1+r)y-1

式中:Cinv为单台BESS的固定投资成本;NBESS为BESSs的安装数量;EBESS,nPBESS,n分别为第n台BESSs的配置容量和功率;ab分别为BESSs的单位功率成本和不同充电倍率(Crate)下的单位容量成本;μCRF为年资金回收率;y为BESSs的使用年限,本文取y=10;r为折现率,采用加权平均资金成本[18]计算.需说明的是,在目标函数计算中,把一年分为多个典型场景计算,同时考虑BESSs的投资运行年限,将所有资金折算成现值.

对于COM,有

COM=ρomn=1NBESS(aPBESS,n+bEBESS,n)

式中:ρom为运维系数,即BESSs年运维费用占设备投资成本的比率,约为1%~10%[19].

对于IS-P,有

IS-P=m-1MdDmn=1NBESSt=1T(ρsell(t)Pdis,n(t)-ρpur(t)Pcha,n(t))

式中:Md为场景数量;Dm为第m种场景对应的天数;ρpur(t)ρsell(t)分别为BESSs在

t时段的购电和售电电价[20];Pcha,n(t)Pdis,n(t)分别为第n台BESSs在t时段内的充电和放电功率;T=24 h.

对于Isub,有

Isub=m=1MdDmn=1NBESSt=1T(λPdis,n(t))

式中:λ为政府给予BESSs发售电量的运营补贴.

配电网的年电压波动.BESSs可通过功率控制减小节点电压波动.在BESSs的选址定容规划中考虑配电网的电压质量,用全年的电压波动总和表示:

F2=m=1MdDmj=1Nnodest=1TVj(t)-V-j

式中:Nnodes为系统节点总数;Vj(t)j节点在t时刻的电压;V-jj节点在T时段内的平均电压.

配电网的年负荷波动.BESSs具有快速能量吞吐能力,可对负荷波动进行平抑.在BESSs的选址定容规划中同时考虑了配电网的功率稳定性,用全年的负荷标准差之和表示:

 F3=m=1MdDmt=1T[(Pload(t)-PPV(t)-Pwind(t)-Pcha/dis(t)-P-L)2/T]

式中: Pload(t)PPV(t)(t)和Pwind(t)分别为t时段内的系统负荷功率、光伏出力和风力发电(以下简称风电)出力;Pcha/dis(t)为BESSs在t时段的充放电功率;P-L为负荷功率的平均值,且

P-L=t=1T(Pload(t)-PPV(t)-Pwind(t)-Pcha/dis(t))T

1.1.2 约束条件

外层模型的约束条件主要包括两类:一类是系统运行约束,即节点功率平衡约束、节点电压约束、并网点功率约束、弃风弃光约束,该类约束以网络潮流的计算结果为依据,以确保整个配电网运行状态的安全性和可靠性,并通过满足弃风弃光率来尽可能促进新能源的就地消纳;另一类是BESSs配置约束,即安装位置约束、功率和容量约束,以满足外层决策变量范围的限制约束.

节点功率平衡约束为

Pi(t)-(Gijcosθij(t)+ Bijsinθij(t))Vi(t)j=1NnodesVj(t)=0Qi(t)-(Gijsinθij(t)- Bijcosθij(t))Vi(t)j=1NnodesVj(t)=0

式中: Pi(t)Qi(t)分别为t时刻节点i的注入有功功率和无功功率;Vi(t)Vj(t)分别为t时段节点i和节点j的电压;θij(t)t时段节点ij之间的电压相位差;GijBij分别节点ij之间线路的电导和电纳.

节点电压约束为

ViminVi(t)Vimax

式中:

VimaxVimin和分别为节点i的电压上下限.根据功率型前推回代法计算,有

Vi(t)={[Vj(t)-(rijPij(t)+xijQij(t))]2+(rijPij(t)+xijQij(t))2}12

式中:Pij(t)Qij(t)分别为流经节点ij之间的有功功率和无功功率;rijxij分别为连接节点ij之间输电线的电阻和电抗.

并网点功率约束为

PgridminPgrid(t)PgridmaxQgridminQgrid(t)Qgridmax

式中:Pgrid(t)Qgrid(t)分别为并网点有功功率和无功功率;PgridmaxPgridminQgridmaxQgridmin分别为其上下限.

弃风弃光约束为

RcurRcur,max

式中:Rcur,max为最大年弃风弃光率,设为5%[21];Rcur为年弃风弃光率,且

Rcur=m=1MdDmt=1T(PPV(t)+Pwind(t)+ Pcha/dis(t)-Pload(t)-Ploss(t))/[m=1MdDmt=1T(PPV(t)+Pwind(t))]

BESSs安装位置约束为

LBESS,nNnodesLBESS,nLgridLBESS,nLBESS,n+1

式中:LBESS,n为第n台BESS的安装节点;Lgrid为配网与主网的联络点.需说明的是,BESSs可以安装在除联络点以外的任一节点,但不能安装在同一个节点.

BESSs功率和容量约束为

EBESSminEBESS,nEBESSmaxPBESSminPBESS,nPBESSmax

式中:EBESSmaxEBESSminPBESSmaxPBESSmin为在安装场地、并网功率和总负荷量等条件下允许的BESSs配置容量的上下限和配置功率的上下限,分别设置为0.3715、3.375 MW·h和 0.3715、 3.375 MW.需说明的是,配置功率的上下限分别按照系统总负荷功率的90%和10%设置,配置容量和配置功率的上下限按照1∶1设置.

1.2 BESSs最优充放电运行

内层模型以BESSs年运营收益最大为目标,建立BESSs充放电运行优化模型:

maxxf(X,x)s.t.g'(X,x)0

式中:f(X,x);g'(X,x)为内层模型需满足的约束条件.

1.2.1 目标函数

在外层规划保证了配电网可靠运行的前提下,BESSs运营方旨在最大化年运营效益,主要包括低储高发套利、降损效益和环境效益.BESSs年运营效益定义为

f=IS-P+Isub+ΔIPloss+ΔIENV

式中:ΔIPloss和ΔIENV分别为配置BESSs后每年的降损效益和环境效益.

对于ΔIPloss,有

  ΔIPloss=m=1MdDmt=1T(ρsell(t)|ΔPloss(t)|)

式中:ΔPloss(t)为配置BESSs后相比配置前配电网在t时段内减少的功率损耗.

对于ΔIENV,有

ΔIENV=m=1MdDmt=1TΔPgrid(t)p=1P(Upup)s.t.Pgrid(t)>0

式中:ΔPgrid(t)为BESSs配置后配电网在t时段内减少的向上级电网购电电量;P为污染物的种类数量;Up为第p种污染物的排放系数;up为第p种污染物的污染罚款[22].

1.2.2 约束条件

内层模型的约束条件除系统运行约束和BESSs配置约束以外,还设置了BESSs充放电约束、荷电状态(State of Charge, SOC)约束和BESSs充电倍率约束,旨在降低BESSs的寿命损耗和运行成本.

BESSs充放电约束为

0Pcha,n(t)PBESS,nηcha-PBESS,n/ηdisPdis,n(t)0

式中:ηchaηdis分别为BESSs的充电效率和放电效率.

SOC约束为

SOCminSOCn(t)SOCmax

式中:SOCmin=20%、SOCmax=90%;SOCn(t)为第n台BESSs在t时段的SOC,有

SOCn(t)=(1-δΔt)SOCn(t-1)+     Pcha,n(t)Δtηcha,SOCn(t)=(1-δΔt)SOCn(t-1)-     Pdis,n(t)Δt/ηdis,

式中:δ为BESSs的自放电率;Δt为BESSs的充放电时长.

BESSs充电倍率约束为

Crate4CCrate=k'C, k'=0.5,1,2,4

式中:C为BESSs在1 h内完全放电时的电流强度.

2 求解方案设计

2.1 基于孔雀优化算法的内层模型求解

在模型内层,采用寻优性能良好的POA求解最优的BESSs充放电运行策略,并在算法内部调用Matpower工具箱采用快速解耦算法进行配电网潮流计算.POA的算法原理和验证参见文献[23],本文仅对算法的寻优机制进行简述.

2.1.1 孔雀优化算法

POA通过雄孔雀求偶、雌孔雀自适应接近雄孔雀、幼孔雀自适应搜索食物源和雄孔雀交互4种寻优机制完成迭代更新,以逐渐逼近问题最优解.过程中,所有孔雀根据其适应度值的排序进行角色分配,具有最高适合度值的前5只为雄孔雀,剩下的前30%为雌孔雀,其余为幼孔雀.

雄孔雀求偶.雄孔雀的位置更新机制可描述为

xc,1=xc,1(t)+Rsxr,1xr,1xc,n=xc,N(t)+σRsxr,Nxr,N,rN<εxc,N(t),

式中:xc,1xc,N分别为第1只和第N只雄孔雀的位置,N=2,3,4,5;σε为决定雄孔雀位置更新的算子,当N=2,3,4,5时,σ=1.5,2,3,5ε=0.9,0.8,0.6,0.3;rN为随机生成的实数;xr,1xr,N为一组随机向量;为雄孔雀围绕食物源旋转的半径,有

  xr=2rand(1,Dim)-1
Rs(t)=Rs0-(Rs0-0)(k/kmax)0.01Rs0=Cv(xub-xlb)

式中:Dim为决策变量的数量;0表示孔雀绕圈旋转的圆心;kkmax分别为当前迭代次数和最大迭代次数;Rs0为初始旋转半径;Cv为雄孔雀旋转因子,设为0.2;xubxlb为决策变量的上下限.

在该机制下,适应度值越高的雄孔雀围绕食物源旋转的概率越大,且绕圈半径越小,因此更趋近于局部最优解.可见,雄孔雀位置代表的决策变量解趋近最优解的能力与其适应度值成正相关,与绕圈半径成负相关.

雌孔雀自适应接近雄孔雀.其位置更新机制为

$\boldsymbol{x}_{\mathrm{h}}=\boldsymbol{x}_{\mathrm{h}}(t)+3 \theta\left(\boldsymbol{x}_{\mathrm{c}, N}-\boldsymbol{x}_{\mathrm{h}}(t)\right), \quad \text { 若 } r_{5} \in A$

式中:xh为雌孔雀的位置;r5为[0,1]范围内的随机数;A为决定雌孔雀位置更新的算子,当N=1,2,3,4,5时,A分别为[0.6,1)[0.4,0.6)[0.2,0.4)[0.1,0.2)[0,0.1);θ探索和全局搜索的算子,有

θ=θ0+(θ1-θ0)k/kmax

式中:θ0=0.1θ1=1.

在该机制下,当θ<1/3时(迭代初期),雌孔雀趋向于所选择的雄孔雀,进行局部勘测;当θ>1/3时(迭代中后期),雌孔雀倾于向所选雄孔雀相对的位置移动,进行全局搜索.因此,较小的θ值有利于雌孔雀在局部勘测过程中寻找高质量的决策变量解;较大的θ值有利于增强算法的随机性和全局搜索能力,避免陷入局部最优.

幼孔雀自适应搜索食物源.幼孔雀向雄孔雀移动的同时借助Levy飞行机制[24]在搜索空间进行随机搜索:

 xcu=xcu(t)+αLevy(xc,1(t)-xcu(t))+  β(xpu(t)-xcu(t))xpu=xc,N(t), r6B

式中:xcuxpu分别为幼孔雀位置和幼孔雀跟随的雄孔雀位置;r6为[0,1]范围内的随机数;B为决定幼孔雀位置更新的算子,当N=1,2,3,4,5时,B分别为[0.8,1)、[0.6,0.8)、[0.4,0.6)、[0.2,0.4)、[0,0.2);αα为随迭代次数动态变化的算子,有

α=α0(α0-α1)(k/kmax)2β=β0+(β1-β0)(k/kmax)0.5

式中:α0=0.9α1=0.4β0=0.1β1=1.

在该机制下,当α>β时(迭代初期),幼孔雀主要进行随机搜索;当β>α时(迭代中后期),幼孔雀逐渐向5只雄孔雀收敛.可见,αβ共同指导幼孔雀的位置更新,引导劣势解向最优解移动,并加快收敛速度.

雄孔雀交互行为.拥有最佳食物源的第1只雄孔雀被视为领导者,第2~4只雄孔雀逐渐向第1只雄孔雀移动:

xc,N=xc,N(t)+θdN+r'NDNDN

式中:

dN=xc,1-xc,NDN=x'r,N-x'r,NdNdNdNdN

r'n为[0,1]范围内的随机数;x'r,N为随机向量,可根据式(26)计算.

2.1.2 BESSs充放电校验

内层优化实际上是一个时序耦合优化问题,即BESSs的运行约束与多个时段关联.在常规网络潮流的基础上进行多时段的扩展,以保证BESSs在各时段的充放电运行均最优,即采用POA求解BESSs在24个时段的充放电功率,基于Matpower工具箱对每个时段单独进行潮流计算.另外,鉴于BESSs在充放电过程中可能会出现超出SOC安全裕度的情况,导致BESSs容量空间未被充分利用或是超出安全运行范围.因此,设置SOC安全约束,对BESSs的充放电功率进行校验和修正:

Pcha,n(t)=min{PBESSmax,[(1-δ)En(t)-     SOCminEBESS,n]/ηcha},s.t.En(t+1)<SOCminEBESS,nPcha,n(t)=min{PBESSmax,[(1-δ)En(t)-     SOCmaxEBESS,n]/ηcha},s.t.En(t+1)>SOCmaxEBESS,nEn(t+1)=(1-δ)En(t)-ηchaPcha,n(t)
Pdis,n(t)=min{-PBESSmax,[(1-δ)En(t)-     SOCminEBESS,n]ηdis},s.t.En(t+1)<SOCminEBESS,nPdis,n(t)=min{-PBESSmax,[(1-δ)En(t)-     SOCmaxEBESS,n]ηdis},s.t.En(t+1)>SOCmaxEBESS,nEn(t+1)=(1-δ)En(t)-Pcha,n(t)/ηdis

式中En(t+1)En(t)分别为第n台BESSst+1时刻和t时刻的剩余容量.另外,为保证BESSs的持续工作能力,设定BESSs每天初始时刻和结束时刻的SOC为50%.

2.2 基于多目标孔雀优化算法的外层模型求解

在模型外层,设计基于Pareto的MOPOA求解BESSs的最佳选址定容方案.首先,根据式(25)~(32)逐代更新外层孔雀种群的位置,以不断更新BESSs选址定容方案的Pareto解集.进而,采用改进理想点决策(Improved Ideal Point Based Decision, IIPBD)[25]方法决策BESSs最佳选址定容方案.

2.2.1 Pareto解集存储与筛选

在算法迭代过程中,采用外部归档集来存储Pareto非支配解集.MOPOA-POA在每次迭代中获得的新非支配解集,须逐一与原非支配解集进行比较并更新外部归档集[26]:① 若新解支配归档集里的某一个或多个解,则将归档集中被支配的解替换为新解;② 若归档集里至少有一个解支配新解,则放弃新解;③ 若新解与归档集里的所有解均互不支配,则将新解放入外部归档集.

为增加解的多样性,当归档集中非支配解的数量超过上限时,需剔除多余的解.采用拥挤距离排序法选出分布较密集的一组非支配解集,并通过轮盘赌法剔除多余解集.非支配解的相邻解数量越多,被剔除的概率就越大.

2.2.2 改进理想点决策 首先,对所有Pareto非支配解的适应度值进行归一化处理:

  yh(χi)=Fh(xi)-FhminFhmax-Fhmin, h=1,2,3

式中: yh(χi)为第i个非支配解下的第h个目标函数的归一化值;FhPareto非支配解的适应度函数;FhmaxFhmin分别为第h个目标函数的最大值和最小值.归一化后的目标理想点为 (0, 0, 0).

其次,计算各个非支配解到目标理想点的趋近度,即欧氏距离平方:

EU(χi)=h=13(yh(χi)-0)2ωh2

式中:ωhh个目标的权重系数.

最后,以所有Pareto非支配解在各个目标上的欧氏距离平方之和最小为原则来客观地确定BESSs选址定容的最佳折中决策方案:

$\chi_{\text {best }}=\arg \min _{i=1,2, \cdots, n_{\mathrm{r} h=1}}^{3}\left(y_{h}\left(\chi_{i}\right)-0\right)^{2} \omega h^{2}$

式中:nr为外部归档集中Pareto非支配解的数量.

2.3 算法应用设计

构建的BESSs双层多目标优化配置模型是一个含多维变量、多个目标、多项约束的高度非线性复杂优化模型,分别应用POA和MOPOA对内层和外层模型进行求解,在算法内部调用Matpower工具箱进行基于快速解耦算法的配电网潮流计算.

(1) 决策变量设置.内层优化的决策变量为由BESSs分别在24个时段的充放电运行功率构成的优化向量,即为POA算法中每只孔雀的位置向量,所有孔雀随迭代过程不断进行位置更新.BESSs外层优化的决策变量包括BESSs的安装位置、配置功率和容量.特别地,BESSs安装位置为离散变量,通过在连续空间取整即可.其中,外层所得BESSs配置功率即为内层BESSs充放电功率的上下限,即式(21);BESSs配置功率和配置容量在设定变量范围内还应满足充电倍率约束,即式(24).

(2) 适应度函数处理.算法的适应度函数须有效结合BESSs优化配置模型的目标和约束条件.其中,不等式约束采用罚函数的形式将其加到原目标函数上进行处理.第h个目标函数对应的算法适应度函数设计为

Fh(xi)=Fh(xi)+ηq

式中:η为惩罚系数,取正无穷;q为违背约束条件的个数.

(3) 算法求解流程:① 根据配电网参数和时序运行数据,初始化外层孔雀群的位置和各适应度值,获得初始化的BESSs安装位置、容量和功率;② 将外层得到的BESSs选址定容方案输入内层模型,并初始化内层孔雀群的位置和适应度值;③ 在内层循环迭代过程中,采用POA对内层孔雀群的位置和适应度值进行更新,直到内层最大迭代次数,则输出内层最优解,即BESSs在24个时段的充放电功率;④ 将内层迭代得到的BESSs最优充放电运行策略输入外层模型,计算外层孔雀群的各适应度值;⑤ 在外层循环迭代过程中,采用MOPOA更新外层孔雀群的位置和外部归档集,并筛选Pareto非支配解集;⑥重复步骤④~⑤,直到外层最大迭代次数,则输出一组完整且均匀分布的外层Pareto非支配解集;⑦利用IIPBD从外层Pareto非支配解集中挑选BESSs选址定容的最佳折中决策方案,即最优的BESSs安装位置、容量和功率;⑧将BESSs最优选址定容方案输入内层模型,重复步骤③~⑦.

综上,所提BESSs双层多目标优化配置的具体流程如图1所示.其中,K为外层优化的迭代次数.

图1

图1   BESSs双层多目标优化配置的流程图

Fig.1   Flow chart of bi-level multi-objective optimal allocation of BESSs


3 运行场景聚类

考虑负荷、风电和光伏出力的不确定性,采用多个典型日场景下的负荷、风电和光伏出力曲线代表全年,更好地反映负荷和新能源在全年的不确定性.首先,采用一种模糊核C-均值聚类算法[27]处理云南省某地的年历史数据,分别得到负荷、风电和光伏的典型日曲线,如图2所示.进而,综合考虑负3种曲线的相关性,通过场景组合得到4种典型场景,每种场景的组合方式及对应天数如表1所示.

图2

图2   负荷、风电和光伏的典型日曲线

Fig.2   Typical diurnal curves of load, wind, and PV power


表1   组合场景聚类结果

Tab.1  Clustering results of composition scenes

场景序号组合方式t/d
1负荷2/风电4/光伏1120
2负荷2/风电4/光伏384
3负荷4/风电4/光伏283
4负荷2/风电3/光伏278

新窗口打开| 下载CSV


4 算例分析

4.1 仿真模型

采用扩展的IEEE-33节点配电网测试系统[28]进行仿真分析,以验证所提方法进行BESSs优化配置的有效性,其拓扑结构如图3所示.其中,G为平衡节点,可向配电网提供缺额功率,也可吸收配电网的多余功率,以保证功率平衡.考虑高比率新能源的渗透,在节点9、20、25分别接入风电1号、2号、3号,在节点28接入光伏1号.仅考虑新能源场站的有功输出,不计无功.电网主要参数为系统基准容量10 MV·A,总负荷功率(3.715+j2.3) MV·A,风电1号、2号、3号有功功率1.2 MW,光伏1号有功功率2.5 MW.选择锂电池作为BESSs的储能元件,相关参数如表2所示[29].为验证所提MOPOA-POA的优越性,仿真算例引入基于MOPSO[17]、NSGAII[16]和MOGWO[30]设计的双层MOPSO-PSO、NSGAII-GA和MOGWO-GWO进行对比.4种算法的基本参数设为一致,即内、外层最大迭代次数均设为100,内、外层种群规模均设为100,外部归档集规模设为20.此外,所有算法均在主频为2.80 GHz 的Intel® CoreTM i5-8400 CPU、内存为8 GB 的计算机和MATLAB 2018b环境下进行优化计算,并调用Matpower工具箱进行潮流计算.

图3

图3   扩展的IEEE-33节点系统拓扑图

Fig.3   Topology of extended IEEE-33 bus system


表2   BESSs的相关参数[29]

Tab.2  Parameters of BESSs[29]

成本参数数值
Cinv/(美元·台-1)1×106
a/(美元·MW-1)137×106
b/(美元·MW-1·h-1)98×105, s.t. Crate=0.5C
190×106, s.t. Crate =C
373×106, s.t. Crate =2C
734×106, s.t. Crate =4C
λ/(美元·kW-1·h-1)0.1
ρom/%5
r/%6.33
ηcha/%95
ηdis/%95
δ/%1

新窗口打开| 下载CSV


4.2 算法性能比较

图4图5(a)分别为4种算法在IEEE-33节点系统获得的两目标和三目标Pareto前沿.需说明的是,不同算法的Pareto前沿由其Pareto非支配解集中的解构成.从图中可知,在同样的迭代次数、种群数量及外部归档集规模条件下,MOPOA-POA获得的两目标和三目标Pareto前沿明显比其他3种算法延伸范围更广、分布更均匀.另外,将4种算法分别得到的Pareto非支配解集对应的目标值全部组合起来,确定该组合解集的Pareto前沿,以此作为所有算法求解该多目标优化问题的近似Pareto最优前沿[31].比较不同算法的Pareto前沿与近似Pareto最优前沿的差距,可知由MOPOA-POA求得的Pareto前沿与该多目标优化问题的近似Pareto最优前沿重合度更高,因此本文所提算法的多目标优化性能较优.图5(b)为采用IIPBD方法从

图4

图4   不同算法的两目标Pareto前沿

Fig.4   Bi-objective Pareto front obtained by different algorithms


图5

图5   三目标Pareto寻优结果

Fig.5   Results of three-objective Pareto front


MOPOA-POA获得的Pareto非支配解集中选取BESSs选址定容方案的最佳折中解.图中数值为各量的归一化值,ω1ω2ω3分别为目标函数F1F2F3的权重,满足ω1+ω2+ω3=1.

表3给为4种算法在外层模型获得的Pareto优化统计结果.可知,MOPOA-POA的Pareto解集在各个目标下的最好(小)值和最差(大)值的差均是最大的,进一步说明其Pareto解集分布广泛,算法的全局搜索能力最强;同时,在BESSs年投资运维成本和配电网年负荷波动两个目标下,MOPOA-POA的最好值是最小的,而配电网年电压波动仅次于MOGWO-GWO,可见算法的局部搜索能力也较强.因此,MOPOA-POA可以获得一组更合理的Pareto非支配解集,为决策者提供多个不同偏好的候选方案.

表3   不同算法下的外层Pareto优化结果

Tab.3  Pareto optimization results of the outer model obtained by different algorithms

算法标准F1/元F2(p.u.)F3/MW
MOPSO-PSO最差值9.6331×10515.7440335.49
最好值4.3513×1050.3392325.22
平均值4.7894×1050.6860335.08
折中解4.5364×1050.6549334.83
NSGAII-GA最差值6.9207×1053.3892338.61
最好值2.5263×1050.3431321.21
平均值4.1647×1050.7902330.98
折中解3.4864×1050.5931328.84
MOGWO-GWO最差值3.6801×10533.6600508.55
最好值-3.7789×1050.2488298.90
平均值-8.3822×1048.2495379.74
折中解-2.1783×1030.4436298.90
MOPOA-POA最差值5.8130×10535.5390538.81
最好值-1.7104×1060.3138282.56
平均值-2.9754×1059.3221387.86
折中解-3.5192×1050.3736328.64

新窗口打开| 下载CSV


另外,采用世代距离(Generational Distance, GD)、反转世代距离(Inverted Generational distance, IGD)、空间分布(Spacing, SP)、广泛性、纯粹多样性(Pure Diversity, PD)、分布度指标(Diversity Metric, DM)、超体积(Hyper Volume, HV)定量评估算法的Pareto解集性能[32],结果如表4所示.可知, MOPOA-POA的GD和IGD指标最低,

表4   不同算法下的Pareto解集性能指标

Tab.4  Performance indexes of Pareto solution sets obtained by different algorithms

算法GDIGDSP广泛性PDDMHV
MOPSO-PSO63142.798×1052.798×1051.117.918×1070.2670.1878
NSGAII-GA29942.266×1052.266×1041.0127.691×1070.3440.1015
MOGWO-GWO94691.038×1051.038×1051.3585.468×1080.6980.2961
MOPOA-POA872.91.568×1041.988×1040.9697.47×1080.8890.2855

新窗口打开| 下载CSV


故收敛性能最佳;具有最小的SP和广泛性指标,故Pareto解集的分布最广泛且均匀;PD和DM指标最高,故分布度最高;而HV指标仅低于MOGWO-GWO,故多样性较好.综上, MOPOA-POA拥有优良的最优解搜索能力和Pareto解集寻优能力.

4.3 结果分析

表5为基于IIPBD的MOPSO-PSO、NSGAII-GA、MOGWO-GWO、MOPOA-POA(算法1~4),以及基于加权求和(Weighted Sum, WS)的POA-POA(算法5)共5种算法获得的BESSs优化配置方案和双层优化结果.其中,为验证IIPBD的有效性,基于WS的POA-POA算法通过简化处理BESSs外层多目标优化模型,即主观地设定3个目标函数的权重均为1/3,从而将多目标问题转化为单目标问题.可知,相比其他算法,基于IIPBD的MOPOA-POA算法得到的BESSs配置方案具有最小的BESSs年投资运维成本和配电网年电压波动,而配电网年负荷波动值仅高于MOGWO-GWO,并且BESSs年运营收益最高.因此,MOPOA-POA算法在保证经济效益的同时可以给配电网带来最佳的综合优化效果,而且IIPBD可以获得更为客观的目标权重系数和折中解,避免采取主观权重而导致对某一目标产生较大偏好,从而获得了较合理的BESSs优化配置方案.

表5   不同算法获得的BESSs配置方案和优化结果

Tab.5  BESSs allocation schemes and optimization results obtained by different algorithms

算法BESSs优化配置方案内层优化结果
f/元
ω1ω2ω3外层优化结果
位置节点容量/(MW·h)功率/MWF1/元F2(p.u.)F3/MW
1(27, 33)(0.375, 0.375)(0.093, 0.09)3.814×1040.5540.4290.0174.536×1050.655334.83
2(25, 21)(0.646, 0.612)(0.162, 0.15)3.319×1050.2960.5530.1513.486×1050.593328.84
3(23, 2)(1.426, 3.375)(0.476, 1.222)1.343×1060.1670.4990.334-2.178×1030.444298.80
4(23, 7)(1.827, 2.445)(0.457, 0.611)1.582×1060.2440.4760.28-3.519×1050.374328.64
5(8, 3)(0.701, 0.484)(0.175, 0.242)4.067×1051/31/31/3-2.825×1055.87328.88

新窗口打开| 下载CSV


图6为基于4种双层多目标优化算法配置BESSs后的IEEE-33节点配电网电压分布(U-)和负荷水平(P-).可知,与未配置BESSs相比,采用4种算法配置BESSs均能在一定程度上改善电压分布和负荷水平.特别地,算法4可显著提高配电网的节点电压分布,并改善配电网负荷波动情况.另外,算法1~4的BESSs优化配置方案均满足配电网弃风弃光约束,分别为3.27%、3.92%、3.64%、2.71%,相较于其他算法,算法4拥有最小弃风弃光率.综上,本文所提方法在解决配电网BESSs规划问题具有一定的有效性和优越性.

图6

图6   IEEE-33节点系统下4种算法的配电网优化结果

Fig.6   Optimization results of distribution network obtained by four algorithms in IEEE 33-bus system


由IEEE-33节点系统参数[26]可知,节点6、7、23、24、29为系统重负荷节点,所提算法将2台BESSs分别配置在节点7和节点23,以满足重负荷地区的电能需要,并尽可能消纳新能源出力.特别地,安装在线路上游的BESS容量最大,可以更好地提升节点电压、降低网损.另外,所配置BESSs的充电倍率均为4C,虽然增加了容量成本,但高充电倍率下电池的恒流充电时间更短,更能满足电网对BESSs的响应要求,也能较大发挥BESSs的充放电能力,从而获得较高的运营收益.然而,较高的充电倍率也会加速电池的容量衰减快慢,影响电池寿命,而本文未考虑BESSs的循环寿命损耗,乃不妥之处.

5 结语

本文设计了一种BESSs优化配置方法.首先,兼顾BESSs投资效益和配电网运行质量,建立BESSs双层多目标优化配置模型;其次,设计寻优性能良好的POA和MOPOA,分别对内层模型和外层模型进行可靠求解.通过内外层的嵌套循环迭代,最终形成综合考虑长期规划经济性和短期运行高效性的BESSs优化配置方案.IEEE-33节点系统的仿真结果显示,所提方法不仅可以实现BESSs投资经济性、配电网电压质量和功率稳定性的最佳权衡,而且可以实现BESSs年运营收益的最大化.

未来研究将提高算法收敛性能并增加算法的种群规模,同时考虑BESSs的循环寿命损耗,以提高配置的合理性.

参考文献

YANG B, WANG J T, CHEN Y X, et al.

Optimal sizing and placement of energy storage system in power grids: A state-of-the-art one-stop handbook

[J]. Journal of Energy Storage, 2020, 32: 101814.

[本文引用: 1]

YANG B, YU L, CHEN Y X, et al.

Modelling, applications, and evaluations of optimal sizing and placement of distributed generations: A critical state-of-the-art survey

[J]. International Journal of Energy Research, 2021, 45(3): 3615-3642.

DOI:10.1002/er.6104      URL     [本文引用: 1]

孙立明, 杨博.

蓄电池/超导混合储能系统非线性鲁棒分数阶控制

[J]. 电力系统保护与控制, 2020, 48(22): 76-83.

[本文引用: 1]

SUN Liming, YANG Bo.

Nonlinear robust fractional-order control of battery/SMES hybrid energy storage systems

[J]. Power System Protection and Control, 2020, 48(22): 76-83.

[本文引用: 1]

YANG B, WANG J B, WANG J T, et al.

Robust fractional-order PID control of supercapacitor energy storage systems for distribution network applications: A perturbation compensation based approach

[J]. Journal of Cleaner Production, 2021, 279: 123362.

DOI:10.1016/j.jclepro.2020.123362      URL     [本文引用: 1]

YANG B, WANG J T, ZHANG X S, et al.

Control of SMES systems in distribution networks with renewable energy integration: A perturbation estimation approach

[J]. Energy, 2020, 202: 117753.

DOI:10.1016/j.energy.2020.117753      URL     [本文引用: 1]

YANG B, WANG J B, ZHANG X S, et al.

Applications of battery/supercapacitor hybrid energy storage systems for electric vehicles using perturbation observer based robust control

[J]. Journal of Power Sources, 2020, 448: 227444.

DOI:10.1016/j.jpowsour.2019.227444      URL     [本文引用: 1]

杨博, 俞磊, 王俊婷, .

基于自适应蝠鲼觅食优化算法的分布式电源选址定容

[J]. 上海交通大学学报, 2021, 55(12): 1673-1688.

[本文引用: 2]

YANG Bo, YU Lei, WANG Junting, et al.

Optimal sizing and placement of distributed generation based on adaptive manta ray foraging optimization

[J]. Journal of Shanghai Jiao Tong University, 2021, 55(12): 1673-1688.

[本文引用: 2]

颜志敏, 王承民, 郑健, .

配电网中蓄电池储能系统的价值评估模型

[J]. 电力自动化设备, 2013, 33(2): 57-61.

[本文引用: 2]

YAN Zhimin, WANG Chengmin, ZHENG Jian, et al.

Value assessment model of battery energy storage system in distribution network

[J]. Electric Power Automation Equipment, 2013, 33(2): 57-61.

[本文引用: 2]

苏海锋, 胡梦锦, 梁志瑞.

基于时序特性含储能装置的分布式电源规划

[J]. 电力自动化设备, 2016, 36(6): 56-63.

[本文引用: 2]

SU Haifeng, HU Mengjin, LIANG Zhirui.

Distributed generation & energy storage planning based on timing characteristics

[J]. Electric Power Automation Equipment, 2016, 36(6): 56-63.

[本文引用: 2]

章美丹, 宋晓喆, 辛焕海, .

计及网损的配电网电池储能站优化运行策略

[J]. 电网技术, 2013, 37(8): 2123-2128.

[本文引用: 2]

ZHANG Meidan, SONG Xiaozhe, XIN Huanhai, et al.

Optimal operation strategy of battery energy storagesystem in distribution networks with consideration of power losses

[J]. Power System Technology, 2013, 37(8): 2123-2128.

[本文引用: 2]

LI R, WANG W, CHEN Z, et al.

Optimal planning of energy storage system in active distribution system based on fuzzy multi-objective bi-level optimization

[J]. Journal of Modern Power Systems and Clean Energy, 2018, 6(2): 342-355.

DOI:10.1007/s40565-017-0332-x      URL     [本文引用: 2]

贾雨龙, 米增强, 刘力卿, .

分布式储能系统接入配电网的容量配置和有序布点综合优化方法

[J]. 电力自动化设备, 2019, 39(4): 1-7.

[本文引用: 1]

JIA Yulong, MI Zengqiang, LIU Liqing, et al.

Comprehensive optimization method of capacity configuration and ordered installation for distributed energy storage system accessing distribution network

[J]. Electric Power Automation Equipment, 2019, 39(4): 1-7.

[本文引用: 1]

张忠会, 雷大勇, 李俊, .

基于自适应ε-支配多目标粒子群算法的含SOP的主动配电网源-网-荷-储双层协同规划模型

[J]. 电网技术, 2022, 6, 2199-2209.

[本文引用: 2]

ZHANG Zhonghui, LEI Dayong, LI Jun, et al.

Source-network-load-storage bi-level collaborative planning model of active distribution network with sop based on adaptive ε-dominating multi-objective particle swarm optimization algorithm

[J]. Power System Technology, 2022, 6, 2199-2209.

[本文引用: 2]

ABU-MOUTI F S, EL-HAWARY M E.

Optimal distributed generation allocation and sizing in distribution systems via artificial bee colony algorithm

[J]. IEEE Transactions on Power Delivery, 2011, 26(4): 2090-2101.

DOI:10.1109/TPWRD.2011.2158246      URL     [本文引用: 2]

叶鹏, 刘思奇, 关多娇, .

基于自适应均衡技术的分布式储能聚合模型及评估方法

[J]. 上海交通大学学报, 2021, 55(12): 1689-1699.

[本文引用: 2]

YE Peng, LIU Siqi, GUAN Duojiao, et al.

An aggregation model and evaluation method of distributed energy storage based on adaptive equalization technology

[J]. Journal of Shanghai Jiao Tong University, 2021, 55(12): 1689-1699.

[本文引用: 2]

SABER H, MOEINI-AGHTAIE M, EHSAN M.

Developing a multi-objective framework for expansion planning studies of distributed energy storage systems (DESSs)

[J]. Energy, 2018, 157: 1079-1089.

DOI:10.1016/j.energy.2018.06.081      URL     [本文引用: 2]

陆立民, 褚国伟, 张涛, .

基于改进多目标粒子群算法的微电网储能优化配置

[J]. 电力系统保护与控制, 2020, 48(15): 116-124.

[本文引用: 2]

LU Limin, CHU Guowei, ZHANG Tao, et al.

Optimal configuration of energy storage in a microgrid based on improved multi-objective particle swarm optimization

[J]. Power System Protection and Control, 2020, 48(15): 116-124.

[本文引用: 2]

马溪原, 吴耀文, 方华亮, .

采用改进细菌觅食算法的风/光/储混合微电网电源优化配置

[J]. 中国电机工程学报, 2011, 31(25): 17-25.

[本文引用: 1]

MA Xiyuan, WU Yaowen, FANG Hualiang, et al.

Optimal sizing of hybrid solar-wind distributed generation in an islanded microgrid using improved bacterial foraging algorithm

[J]. Proceedings of the CSEE, 2011, 31(25): 17-25.

[本文引用: 1]

何颖源, 陈永翀, 刘勇, .

储能的度电成本和里程成本分析

[J]. 电工电能新技术, 2019, 38(9): 1-10.

[本文引用: 1]

HE Yingyuan, CHEN Yongchong, LIU Yong, et al.

Analysis of cost per kilowatt-hour and cost per mileage for energy storage technologies

[J]. Advanced Technology of Electrical Engineering and Energy, 2019, 38(9): 1-10.

[本文引用: 1]

张志义, 余涛, 王德志, .

基于集成学习的含电气热商业楼宇群的分时电价求解

[J]. 中国电机工程学报, 2019, 39(1): 112-125.

[本文引用: 1]

ZHANG Zhiyi, YU Tao, WANG Dezhi, et al.

Optimal solution of time-of-use price based on ensemble learning for electricity-gas-heat commercial building

[J]. Proceedings of the CSEE, 2019, 39(1): 112-125.

[本文引用: 1]

中华人民共和国国家发展和改革委员会, 国家能源局. 解决弃水弃风弃光问题实施方案[R]. 北京: 中华人民共和国国家发展和改革委员会, 国家能源局, 2017.

[本文引用: 1]

National Development and Reform Commission, National Energy Administration. Implementation plan to solve the problem of water, air and light abandonment[R]. Beijing: National Development and Reform Commission, National Energy Administration, 2017.

[本文引用: 1]

林君豪, 张焰, 陈思, .

考虑可控负荷影响的主动配电系统分布式电源优化配置

[J]. 电力自动化设备, 2016, 36(9): 46-53.

[本文引用: 1]

LIN Junhao, ZHANG Yan, CHEN Si, et al.

Optimal DG allocation considering effect of controllable load for active distribution system

[J]. Electric Power Automation Equipment, 2016, 36(9): 46-53.

[本文引用: 1]

GANDOMI A H, YANG X S, ALAVI A H.

Cuckoo search algorithm: A metaheuristic approach to solve structural optimization problems

[J]. Engineering with Computers, 2013, 29(1): 17-35.

DOI:10.1007/s00366-011-0241-y      URL     [本文引用: 1]

YANG B, WANG J B, ZHANG X S, et al.

Applications of battery/supercapacitor hybrid energy storage systems for electric vehicles using perturbation observer based robust control

[J]. Journal of Power Sources, 2020, 448: 227444.

DOI:10.1016/j.jpowsour.2019.227444      URL     [本文引用: 1]

杨蕾, 吴琛, 黄伟, .

含高比例风光新能源电网的多目标无功优化算法

[J]. 电力建设, 2020, 41(7): 100-109.

DOI:10.12204/j.issn.1000-7229.2020.07.013      [本文引用: 1]

为适应新能源大量接入电网的趋势,基于不同时刻的风速、光照强度、温度等气象条件信息,评估出风光新能源的无功调节容量,搭建了含高比例风光新能源参与调控的电网多目标无功优化模型。为快速获得电网中变压器分接头档位调节、无功补偿设备投切、传统发电机组电压调节以及风光的无功输出等控制措施的帕累托最优解集,采用寻优性能高效的多目标樽海鞘群算法(multi-objective salp swarm algorithm, MSSA)进行无功优化求解。为更客观找出电网线损、电压偏差、静态电压稳定裕度等不同目标之间的折中解,采用改进的理想点法进行多目标最优解集决策。最后,利用扩展的IEEE标准9节点和39节点算例进行仿真分析,并引入传统多目标智能优化算法来进行比较验证。仿真结果表明:与其他2种传统多目标智能优化算法相比,所提算法获得的帕累托前沿分布更广、更均匀;利用改进理想点法进行决策之后,可有效降低电网的线损和电压偏差,同时提高了电网的静态电压稳定裕度。

YANG Lei, WU Chen, HUANG Wei, et al.

Pareto-based multi-objective reactive power optimization for power grid with high-penetration wind and solar renewable energies

[J]. Electric Power Construction, 2020, 41(7): 100-109.

DOI:10.12204/j.issn.1000-7229.2020.07.013      [本文引用: 1]

To adapt the trend of high-penetration renewable energies paralleled in power grid, this paper constructs a multi-objective reactive power optimization for power grid with the controlled participation of high-penetration wind and solar renewable energies. Particularly, the reactive power regulation capacities of renewable energies are evaluated according to the wind speed, solar irradiation, and temperature in different time. To obtain the optimal dispatch scheme of transformer taps, shunt capacitor states, voltage outputs of generators, and reactive power outputs of renewable energies, a multi-objective salp swarm algorithm (MSSA) is employed for the multi-objective reactive power optimization. Then an improved ideal-point based decision method is designed to select a compromise solution among multiple non-dominated points, thus three objectives of power loss, voltage deviation, and static voltage stability margin can be properly balanced. Finally, an extended IEEE 9-bus system and an extended IEEE 39-bus system are used to evaluate the performance of the proposed algorithm compared with conventional multi-objective intelligent optimization algorithms. Simulation results demonstrate that the proposed algorithm can obtain a widely spread and well-distributed Pareto front compared with conventional multi-objective optimization algorithms. Moreover, the improved ideal-point based decision method not only can effectively reduce the power loss and voltage deviation, but also can improve the static voltage stability margin.

FARAMARZI A, HEIDARINEJAD M, STEPHENS B, et al.

Equilibrium optimizer: A novel optimization algorithm

[J]. Knowledge-Based Systems, 2020, 191: 105190.

DOI:10.1016/j.knosys.2019.105190      URL     [本文引用: 2]

伍忠东, 高新波, 谢维信.

基于核方法的模糊聚类算法

[J]. 西安电子科技大学学报, 2004, 31(4): 533-537.

[本文引用: 1]

WU Zhongdong, GAO Xinbo, XIE Weixin.

A study of a new fuzzy clustering algorithm based on the kernel method

[J]. Journal of Xidian University, 2004, 31(4): 533-537.

[本文引用: 1]

GOSWAMI S K, BASU S K.

A new algorithm for the reconfiguration of distribution feeders for loss minimization

[J]. IEEE Transactions on Power Delivery, 1992, 7(3): 1484-1491.

DOI:10.1109/61.141868      URL     [本文引用: 1]

RAN F, TIMOTHY R, ROBERT M.

Utility-scale photovoltaics-plus-energy storage system costs benchmark

[R]. USA: National Renewable Energy Laboratory, 2018.

[本文引用: 3]

MIRJALILI S, SAREMI S, MIRJALILI S M, et al.

Multi-objective grey wolf optimizer: A novel algorithm for multi-criterion optimization

[J]. Expert Systems with Applications, 2016, 47: 106-119.

DOI:10.1016/j.eswa.2015.10.039      URL     [本文引用: 1]

DURILLO J J, NEBRO A J, COELLO COELLO C A, et al.

A comparative study of the effect of parameter scalability in multi-objective metaheuristics

[C]//2008 IEEE Congress on Evolutionary Computation (IEEE World Congress on Computational Intelligence). Hong Kong, China: IEEE, 2008: 1893-1900.

[本文引用: 1]

王丽萍, 邱飞岳. 复杂多目标问题的优化方法及应用[M]. 北京: 科学出版社, 2018.

[本文引用: 1]

WANG Liping, QIU Feiyue. Complex multi-objective optimization and it’s applications[M]. Beijing: Science Press, 2018.

[本文引用: 1]

/