随着全球能源危机日益严重及社会对环境保护需求的增长,人们迫切需要更高效、更环保、更经济的电力供应[1],因此,提出了智能电网的概念,通过智能化管理电力的生产、运输和零售过程来保证电力系统安全、可靠、经济的运行[2].
目前,智能电网通过大力发展可再生资源发电来减少对火力发电的依赖[3],但可再生资源发电具有随机、波动和间歇性强的缺点,将其并入电网将会对电网产生冲击并影响电力系统的调峰能力[4].随着电池储能和超级电容储能等新型储能技术的快速发展和产业化[5-6],结合可再生资源发电功率的预测技术与海量储能单元功率双向流动的特性来平抑大规模可再生资源发电并网所产生的功率波动,已逐渐成为研究及应用的热点[7-8].
由于太阳能和风能等可再生能源具有不确定性,故需要实时调整智能电网储能系统的输出功率来补偿实际功率与计划功率间的偏差[9].然而,海量储能单元存在地理位置分散难以获取同步和全局信息的问题,传统的集中式调度方法灵活性低,无法满足智能电网对储能单元“即插即用”的需求[9].而分布式调度方法仅需局部信息的交互即可实现储能单元的实时功率分配,具有更好的灵活性和鲁棒性[10].分布式调度策略通过应用一致性算法来实现储能单元的实时功率分配,如文献[11]通过仿真分析揭示了不平衡量的反馈系数、通信拓扑结构和网络规模对分布式调度算法收敛速度的影响;文献[12]利用矩阵论和图论分析了通信时延与一致性算法收敛速度间的关系.
本文基于一致性算法[13]实现储能单元的调度,从而抑制可再生能源发电所引起的功率波动,该调度考虑储能单元的内阻对充放电效率的影响,通过引入一致性算法实现充放电效率的最大化.
假设有n个储能单元分散配置在配电网中,构成一个分布式储能网络,每个储能单元均有一个计划输出功率,该功率由调度策略决定.当可再生能源的输出功率发生波动时,需要根据一定的原则调整储能单元的实时功率.为了实现抑制可再生资源发电功率波动的效果,本文使用一致性算法实现功率在储能单元中的公平分配,即在保证储能单元充放电状态不变的同时,确保每个储能单元的输出值与计划输出功率成比例.
智能电网的有功平衡式可表示为
(1)
式中:PB,i为储能单元的充放电功率;PG,k、PD,j分别为发电机与负荷的需求功率;SG、SD、SB分别为智能电网中的发电机、负荷和储能单元.电网中储能单元需要提供的总有功功率为
(2)
当PL>0时,表示储能单元充电;当PL<0时,表示储能单元放电.为了保证智能电网功率的平衡,本文使用调度策略来控制储能单元的充放电功率PB,i(i∈SB)来调整PL,即
(3)
由于储能单元也存在一定的内阻,故在充放电时会存在功率损耗,即
(4)
(5)
式中:分别为储能单元的充放电功率与充放电效率;ηC,i、ηD,i为损耗系数,0<ηC,i<1<ηD,i.
由于放电状态可以由充电状态表示出来,本文为了简化分析,仅考虑储能单元运行在充电状态时的情况.储能单元的实际充电功率为
(6)
为了使智能电网更经济运行,需要最小化网络损失,即通过调整储能单元的充电参数最大化式(6).本文通过协调控制目标函数来保证智能电网在快速变化时仍可进行快速、稳定地调控,即
(7)
式中,为储能单元的最小、最大充电功率.
由文献[13]可知,储能单元的充电速率αi与充电功率PB,i线性表达式为
αi=ai-biPB,i
(8)
式中,ai、bi为常系数.
将式(8)代入式(7)可得目标函数为
(9)
文中定义储能单元的边际成本(λi)为式(9)相对于充电功率PB,i的偏导,即
λi=ai-2biPB,i
(10)
一致性算法即连通图中所有顶点收敛到同一个与各顶点初始值大小无关的值,所有顶点值的和不变,即对于离散系统则有
其中,δj(k)为顶点i第k次迭代时的值,Wi,j为顶点i与j所连成边的权重.
智能电网各储能单元B之间的通信拓扑构成一个连通图g=(B,E,A),其中,E为相邻两个储能单元的连接线,A为其对应的权重.第i个储能单元的入邻居为表示储能单元i能从其他的入邻居接收到信息.假设网络中存在一个虚拟顶点向B中其他点发送功率偏差信息,该虚拟节点称为一致性算法的“Leader”,而其他顶点用于分配总的功率偏差称为一致性算法的“Follower”.用顶点0表示虚拟顶点“Leader”,用
表示“Follower”,“Leader”向“Follower”平均分配的总功率为ΔPΣ.使用一致性算法求解式(9)的最优解,则有
(11)
式中:Pi(0)为储能单元i的初始功率,本文实验设置为40 kW;ΔPi(0)为需要调整的功率偏差值.
实时功率分配策略为
(12)
式中:ε为扰动系数,其值为足够小的常数;为一致性部分,表示所有储能单元均达到同一个值;ci=uiPi(0),ui为第i个储能单元的充放电状态,ui=1表示放电,ui=-1表示充电.由于智能电网的通信拓扑网络为强连通图,所以系统中各储能单元均能得到其他储能单元的信息(ΔPi(k),xi(k)).整合这些信息可以得到全局信息xi(k),根据xi(k)可以调整各储能单元的输出功率ΔPi(k),从而平抑可再生能源引入的发电功率.
将式(12)转化为矩阵的形式来分析分布式调度策略的性质和收敛性,即
(13)
式中:λ、ΔP、x、C、Q分别为λi、ΔPi、xi、Ci、q的向量形式;ΔP(k+1)-ΔP(k)为控制所有储能单元达到最优值的反馈机制.由式(13)则有
x(k+1)=Qx(k)-(ΔP(k+1)-ΔP(k))
=x(k)-(ΔP(k+1)-ΔP(k))
⟹x(k+1)+ΔP(k+1)
=x(k)+ΔP(k)
(14)
式中,x(k)+ΔP(k)取任意的k均为一常数.由初始值xi(0)和ΔPi(0)可得到∑xi(0)+ΔPi(0)=ΔPΣ,因此,x(k)=ΔPΣ-ΔP(k)为需要调节的功率偏差量,其值可由储能单元系统的分布式调度获得.
将式(13)代入式(14)则有
x(k+1)=(Q-εC)x(k)+C(I-P)λ(k)
(15)
将式(14)和式(15)矩阵化则有
(16)
定义
(17)
则式(17)所示的系统矩阵可以看做被εΔ扰动的下三角M矩阵,构造向量为
(18)
式中:μ为变换系数;且有
(19)
式中,VT、φ满足VTφ=I.
当扰动ε>0时,随着ε的增大,矩阵M的特征值k2逐渐减小.假设ε≤ξ1,则|k2|<1.因为M的特征值与ε有关,故存在上界ξ2,使得当ε≤ξ2时,|kj|<1,j=3,4,…,N.综上所述,当ε≤min(ξ1,ξ2)时,M具有特征值k1=1,且其他特征值位于单位圆上.此时,特征向量为[1,0]T,当k→∞时,[λ(k),x(k)]T收敛到[1,0]T,所以当ε足够小时,系统能稳定运行且所有储能单元能收敛到最优特征值λ*.
为了验证本文所提基于一致性算法的储能单元分布式调度算法的有效性,使用如图1所示的IEEE27节点系统进行仿真测试,图中虚线表示通信链路,实线表示配电线路.该系统的712、718和725母线为风力发电系统,且将虚拟命令节点连接704和713母线,即有本文设置储能单元的总输入功率为2 000 kW,各储能单元初始功率为40 kW,扰动系数ε为0.000 1,由各储能单元的铭牌可得其常系数,如表1所示,由表1和式(8)可以得到各储能单元的充电速率.
图1 IEEE27节点系统拓扑图
Fig.1 Topological diagram of IEEE27 node system
27节点仿真结果如图2所示,图2a、b、c展示了储能单元1、13和27的λ、PB和ΔPB的变化情况,图2d为27条母线总PL变化情况.根据图2可知,27个储能单元能收敛到相同值且总储能功率为2 000 kW.
表1 储能单元常系数
Tab.1 Constant coefficients of energy storage units
节点编号aibi10.8730.008320.8550.007630.7600.006840.9260.006550.8550.007160.7800.007170.8270.007680.6310.006390.7320.0053100.8280.0085110.9190.0096120.9300.0045130.8440.0102140.8080.0058节点编号aibi150.9870.0016160.8930.0081170.8780.0108180.8720.0019190.8300.0087200.8170.0048210.9060.0067220.9890.0057230.9980.0064240.9720.0065250.8520.0063260.9990.0025270.9270.0124
图2 分配策略有效性仿真结果
Fig.2 Simulation results for effectiveness of distributed strategy
本文在图1系统保证总储能值不变的情况下接入第28个储能单元,接入后的系统拓扑图如图3所示,第28个储能单元的ai为0.874,bi为0.008 5.仿真实验结果如图4所示,图4a、b、c展示了储能单元1、13、27和28的λ、PB和ΔPB的变化情况,由图4a、b、c可以看出,在接入节点28后,三项指标值虽然有一定的变化,但最终均收敛于一致;图4d为28条母线总PL变化情况,从图4d中可以看出,在接入新的储能单元后,系统边界成本先出现短暂波动,然后收敛并趋于一致,且未引起PL较大变化,表明系统能满足对储能单元即插即用特性的需求,可以针对太阳能和风能的不确定接入作出良好的反应,更适合用于实际工况.
图3 IEEE28节点系统拓扑图
Fig.3 Topological diagram of IEEE28 node system
为了实时调整智能电网储能系统的输出功率,以补偿实际功率与计划功率间的偏差,本文建立了储能单元分布式调度模型,并提出了一种基于一致性算法的智能电网储能单元分布式调度策略.该调度考虑了储能单元内阻对充放电效率的影响,通过引入一致性算法实现充放电效率的最大化.IEEE27节点系统和即插即用特性仿真结果表明,所提出的分布式调度策略能保证所有储能单元可以收敛到同一个最优值,克服了可再生资源发电具有随机、波动和间歇性强的缺点,满足太阳能及风能等不确定工况需要.
图4 即插即用特性仿真结果
Fig.4 Simulation results of plug and play characteristics
[1]Kang J,Park J,Jung L,et al.Development of a systematic,self-consistent algorithm for the K-DEMO steady-state operation scenario [J].Nuclear Fusion,2017(6):513-525.
[2]赵红毅,刘利坚.一种分布式系统进程调度方法研究 [J].电子科技,2010,23(6):55-58.
(ZHAO Hong-yi,LIU Li-jian.A process scheduling method for distributed systems [J].Electronic Science and Technology,2010,23(6):55-58.)
[3]Guan Y,Meng L,Li C,et al.A dynamic consensus algorithm to adjust virtual impedance loops for discharge rate balancing of AC microgrid energy storage units [J].IEEE Transactions on Smart Grid,2017(9):11-19.
[4]Kiaei I,Lotfifard S.Tube-based model predictive control of energy storage systems for enhancing transient stability of power systems [J].IEEE Transactions on Smart Grid,2017(9):42-50.
[5]杜家兵,陈衍鹏,梁满发.基于分布式实时调度策略的智能电网控制系统的设计与实现 [J].电子设计工程,2018,26(7):119-122.
(DU Jia-bing,CHEN Yan-peng,LIANG Man-fa.Design and implementation of smart grid control system based on distributed real-time scheduling strategy [J].Electronic Design Engineering,2018,26(7):119-122.)
[6]米阳,吴彦伟,纪宏澎,等.基于多组储能动态调节的独立直流微电网协调控制 [J].电力自动化设备,2017,37(5):170-176.
(MI Yang,WU Yan-wei,JI Hong-peng,et al.Coordinated control of independent DC microgrid based on multi group energy storage dynamic regulation [J].Power Automation Equipment,2017,37(5):170-176.)
[7]吴晓刚,侯维祥,帅志斌,等.电动汽车复合储能系统的功率分配优化研究 [J].电机与控制学报,2017,21(11):114-124.
(WU Xiao-gang,HOU Wei-xiang,SHUAI Zhi-bin,et al.Power distribution optimization of electric vehicle hybrid energy storage system [J].Electric Machines and Control,2017,21(11):114-124.)
[8]雷为民,苗友忠,罗玮,等.电网电压不平衡情况下双馈风力发电变流器控制 [J].沈阳工业大学学报,2017,39(6):612-616.
(LEI Wei-min,MIAO You-zhong,LUO Wei,et al.Control of doubly-fed wind power generation conver-ter under unbalanced grid voltage [J].Journal of Shen-yang University of Technology,2017,39(6):612-616.)
[9]从旭光,任永峰,胡志帅,等.基于两级全钒液流电池的风电储能系统的功率优化分配控制 [J].水电能源科学,2017(7):210-214.
(CONG Xu-guang,REN Yong-feng,HU Zhi-shuai,et al.Storage of wind power system optimal power allocation control based on two of all vanadium flow battery [J].Water Resources and Power,2017(7):210-214.)
[10]王帅,赵克,安群涛,等.钻井直流微电网冲击功率的混合储能平衡技术 [J].电机与控制学报,2017,21(4):53-61.
(WANG Shuai,ZHAO Ke,AN Qun-tao,et al.Hybrid energy storage technology for drilling DC microgrid impact power [J].Electric Machines and Control,2017,21(4):53-61.)
[11]杨毅,刘洁,于建成,等.分布式储能优化配置影响因素的评价方法研究 [J].现代电子技术,2015,38(2):152-157.
(YANG Yi,LIU Jie,YU Jian-cheng,et al.Evaluation method of distributed energy storage system in consi-deration of optimization configuration influence factors [J].Modern Electronics Technique,2015,38(2):152-157.)
[12]吴杰,丁明.采用自适应小波包分解的混合储能平抑风电波动控制策略 [J].电力系统自动化,2017,41(3):7-12.
(WU Jie,DING Ming.Control strategy of hybrid energy storage to suppress wind power fluctuation using adaptive wavelet packet decomposition [J].Automation of Electric Power Systems,2017,41(3):7-12.)
[13]杨家豪.基于一致性算法的孤岛型微电网群实时协同功率分配 [J].电力系统自动化,2017,41(5):8-15.
(YANG Jia-hao.Real-time cooperative power allocation of islanding micro grid based on consistency algorithm [J].Automation of Electric Power Systems,2017,41(5):8-15.)