Model predictive control of battery-supercapacitor hybrid energy storage system
-
摘要: 为匹配中国环流三号装置磁体线圈功率和能量需求,基于模型预测控制理论,设计电池-超级电容混合储能系统放电策略。以环向场线圈为储能系统输出端负载,建立系统预测模型,根据电池、超级电容特性和负载能量需求设计目标函数,实时求解最优开关序列。对电池储能子系统、超级电容储能子系统分别采用长周期慢控、短周期快控,实现电池稳定放电和超级电容瞬态响应。基于MATLAB/Simulink平台进行仿真验证,混合储能系统稳定输出满足负载需求的平顶电流,其电流纹波为0.22%,验证控制策略有效性。Abstract: A discharge strategy for a battery-supercapacitor hybrid energy storage system is designed based on model predictive control theory to match the power and energy requirement of the magnet coil of HL-3. Using toroidal field coil as the load of the energy storage system, the mathematical model of the system and the objective function based on battery/supercapacitor characteristics and energy demands of the load is established. The optimal switching sequence is solved in real time. Long cycle control is applied on battery energy storage system to achieve stable discharge of the battery, while short cycle control is applied on supercapacitor energy storage system to achieve transient response of the supercapacitor. Simulation experiments are conducted using MATLAB/Simulink. The hybrid energy storage system stably outputs a flat top current that meets the load demand, with a current ripple of 0.22%. The simulation results verify the effectiveness of the proposed control method.
-
中国环流三号(HL-3)装置是我国自主建造的新一代磁约束聚变实验装置,是目前中国设计参数最高、规模最大的常规托卡马克装置[1]。环向场(Toroidal Field, TF)线圈是HL-3装置的重要组成部分,用来形成环形的主磁场。为满足装置运行需求,TF线圈电源应提供平顶达7 s的140 kA电流,电流纹波小于1%,需要200 MV·A容量和
1300 MJ能量[2-3]。目前,HL-3装置TF线圈通过两套90 MV·A飞轮储能式脉冲发电机组供电,单套机组最大释能约为500 MJ(转速范围1650 ~1200 r/min),无论是瞬时最大输出功率还是能量均无法满足环向场7 s平顶140 kA放电的需求。为了匹配HL-3环向场功率和能量的物理需求,提供更高的约束能量,必须展开新的能量供给系统研究。储能介质实现能量的转化、储存和释放。以蓄电池、锂离子电池等为代表的能量型储能介质具有能量密度大的优点,但功率密度小;以超级电容、飞轮储能等为代表的功率型储能介质,具有功率密度大、响应速度快的优点,但能量密度小[4]。TF线圈是一个强感性负载,放电时长为秒级,需求高功率、大容量能量输入,要求储能介质具有快响应速度、高能量密度和高功率密度,单一的储能介质缺乏经济性。因此,采用电池-超级电容混合储能方式,电池组稳定提供大容量能量输出,超级电容组提供瞬时高功率,两种储能介质配合工作满足负载需求。
目前,混合储能控制策略研究集中于电动汽车和新能源发电领域[5-12],缺乏对大电流、高功率、大容量能量需求的强感性负载的放电控制策略研究。模型预测控制(Model Predictive Control, MPC)适用于非线性系统的实时优化,其核心思想是根据当前状态,结合系统模型预测系统未来状态,根据优化目标实时更新控制策略[13]。文献[14]采用MPC方法结合基于规则的控制策略,对应用于电动汽车的电池-超级电容混合储能系统进行功率分配控制,降低电池峰值放电、延长电池寿命;文献[11,15]将基于MPC控制的混合储能系统应用于微电网,有效抑制电压波动,调节电力系统频率,提高系统暂态性能。MPC方法具有实时性且物理意义明确,适用于非线性的、具有变优化目标和多优化目标的系统控制。因此,针对HL-3装置的TF线圈放电需求,基于MPC方法开展电池-超级电容混合储能系统放电控制策略研究,以实现储能子系统间合理的功率分配。
1. 混合储能系统结构组成及模型建立
电池-超级电容混合储能系统结构可分为被动式、半控式、全控式,其中,只有全控式结构可实现对多组储能装置输出功率的主动控制,且各个子系统的控制相对独立,因此选择采用全控式结构[16-17]。如图1所示为混合储能系统结构示意图,电池组和超级电容组经DC-DC变换器连接至直流母线为负载供电。一个电池储能系统与一个超级电容储能系统构成一个子模块,如图1(a)所示;n个子模块连接至直流母线,共同为负载供电,如图1(b)所示。DC-DC变换器采用Cuk电路,其具有输入输出电流连续、输出电压范围灵活的优点,满足装置实验运行需求[18-19]。
1.1 电池组模型
使用规格统一的电池单体、超级电容单体构成储能单元。由于储能装置具有电压均衡电路,因此为便于分析,假定每个单体的电压、电流和荷电状态(State of Charge, SOC)相等,储能装置组工作于电流、电压均衡的理想状态。
针对电池单体,采用简化电池模型;为了获得所需的电压、电流,将电池单体经串并联后得到电池组,其等效模型由一个恒定电压源VB串联一个等效内阻RB构成。设vBT为电池组输出电压,iB为电池组输出电流,得到电池组数学模型
vBT(t)=VB−iB(t)RB (1) 电池组荷电状态SSOC_B采用安时积分法计算
SSOC\_B(t)=SSOC\_B(0)−1CB∫t0iB(τ)dτ (2) 式中:CB为电池组容量(A·s)。
1.2 超级电容组模型
针对超级电容,采用经典简化模型。同样经过超级电容单体的串并联得到超级电容组,其等效模型由等效电容C串联等效内阻RC构成。设vC为等效电容电压,vSC为超级电容组输出电压,iC为超级电容组输出电流,得到超级电容组数学模型
iC(t)=−CdvC(t)dtvSC(t)=vC(t)−RCiC(t) (3) 超级电容组荷电状态SSOC_SC按照其存储能量与最大标称储能上限之比计算,即
SSOC\_SC(t)=12Cv2C(t)12Cv2C_max=v2C(t)v2C_max (4) 1.3 TF线圈模型
HL-3装置的磁体线圈是一类强感性负载。表1所示为TF线圈参数,其中,LTF表示TF线圈电感,RTF表示TF线圈电阻,tflat表示平顶放电时间,iTF_flat表示TF线圈平顶段电流。
表 1 TF线圈参数Table 1. Parameters of TF coilLTF/mH RTF/mΩ tflat/s iTF_flat/kA 32 5.6 7 140 TF线圈模型由LTF、RTF串联构成,vTF、iTF分别表示TF线圈电压、电流,其数学模型为
vTF(t)=LTFdiTF(t)dt+RTFiTF(t) (5) 1.4 Cuk变换器模型
图2所示为Cuk变换器电路结构,vin为输入电压,L1、i1为输入电感及其电流,L2、i2为输出电感及其电流,C1、v1为中间电容及其端电压,C2、v2为输出电容及其电压。假设C1足够大,使Cuk电路工作于中间电容电压连续模式(Continuous Capacitor Voltage Mode, CCVM)。
Cuk电路工作于连续导通模式(Continuous Conduction Mode, CCM)时,存在两种工作状态。
工作状态1:开关S导通,二极管D截止,输入端经开关管S向电感L1充能,中间电容C1经开关管S向电感L2及负载放电,电感电流i1与i2均线性上升。
工作状态2:开关S关断,二极管D导通,电感L1向中间电容C1充能,电感L2继续经由D续流,电感电流i1、i2均下降,此时,i1或i2可能为负值,但二者之和始终为正。
Cuk电路工作于断续导通模式(Discontinuous Conduction Mode, DCM)时,在一个开关周期内存在电感电流i1、i2之和始终为0的阶段,在CCVM状态下,该种工况仅出现在开关S关断时,输入端、输出端直接贯通,若i1=−i2≠0,则输入端电流反向。为防止电池进入瞬时反向充电状态,应避免电池储能系统的Cuk电路输入端电流反向,通过对电池组输出电流的控制可实现该要求;超级电容具有吸收或释放瞬时功率的能力,且耐过充、放电能力优良,因此,对于超级电容储能系统,允许该种工况出现,同时应监测超级电容组荷电状态SSOC_SC,保证超级电容处于正常工作状态。
依据电路原理,得到CCM Cuk变换器的状态空间方程如式(6)所示,其中g表示控制信号,当g=1时,开关S开通,当g=0时,开关S关断。
{L1di1(t)dt=vin(t)−(1−g)v1(t)L2di2(t)dt=gv1(t)−v2(t)C1dv1(t)dt=−gi2(t)+(1−g)i1(t)C2dv2(t)dt=i2(t)−iout(t) (6) 2. 基于模型预测控制的放电策略
2.1 预测模型
以图1(a)所示的单个子模块为研究对象,以TF线圈作为系统负载。混合储能系统由10个子模块构成,每个子模块需要为负载提供7 s平顶的14 kA电流。图3所示为系统主电路原理图,电池储能系统与超级电容储能系统经过输出滤波电容C3为TF线圈供电。
设采样时间间隔为T。根据式(6)得到储能系统状态预测方程如式(7)、式(8)。其中,G1(k)、G2(k)分别表示电池储能系统开关S1状态、超级电容储能系统开关S2状态,且
G′1(k)=1−G1(k) ,G′2(k)=1−G2(k) 。当G1(k)=1 时,开关S1闭合;当G1(k)=0 时,开关S1断开;G2(k)、S2同理。[i1(k+1)i2(k+1)v1(k+1)]=[10−G′1(k)TL101G1(k)TL2G′1(k)TC1−G1(k)TC11][i1(k)i2(k)v1(k)]+[TL100−TL200][vin1(k)vTF(k)] (7) [i3(k+1)i4(k+1)v2(k+1)]=[10−G′2(k)TL301G2(k)TL4G′2(k)TC2−G2(k)TC21][i3(k)i4(k)v2(k)]+[TL300−TL400][vin2(k)vTF(k)] (8) 由式(1)得到电池组端电压预测方程
vin1(k+1)=VB−RBi1(k+1) (9) 由式(3)得到超级电容组端电压预测方程
vin2(k+1)=(RC−TC)i3(k)−RCi3(k+1)+vin2(k) (10) 根据图3所示电路拓扑,得到输出电容C3支路微分方程
C3dvTF(t)dt=i2(t)+i4(t)−iTF(t) (11) 对式(5)、式(11)使用前向差分法,得到输出电压电流预测方程
{iTF(k+1)=TLTFvTF(k)+(1−TRTFLTF)iTF(k)vTF(k+1)=TC3[i2(k)+i4(k)−iTF(k)]+vTF(k) (12) 2.2 反馈校正
由于存在电路元件杂散参数、外界干扰等不确定因素,引起预测值与实际值的偏离,因此通过对输出误差加权的方式对式(7)~(12)的预测输出进行修正,得到
xm(k+1)=x(k+1)+fe(k)=x(k+1)+f[x(k)−xm(k)] (13) 式中:e(k)为第k周期采样值x(k)与第(k−1)采样周期时对第k周期的预测值xm(k)之差,f为误差校正系数。
2.3 滚动优化
设电池储能系统和超级电容储能系统的预测步数分别为N1、N2,控制周期分别为T1、T2,且T1≥N2T2。
电池储能系统优化目标是使电池组稳定输出电流并为TF线圈提供主要功率,其代价函数为
J1=α1N1∑j=1[i1(k+j)−i1ref]2+α2N1∑j=1[i2(k+j)−i2ref]2 (14) 式中:i1ref、i2ref分别为电池组端电流和输出电感电流参考值,α1、α2为权重系数,且α1+α2=1。采用混合储能方式的主要目的是保持电池组稳定输出功率,负载端的瞬时功率由超级电容储能系统吸收或释放,对电池储能系统的控制着重于电池组端电流的控制,因此权重系数应有α1>α2。对于i1ref,可取电池组额定输出电流IBN;对于i2ref,TF线圈从电流上升阶段到平顶放电阶段存在电路工况的变化,为提高电路对工况变化的动态响应能力,增强输出稳定性,在i2ref的基础上增加电压反馈补偿项εv,得到输出电感电流参考值i2ref*为
{i2ref*=i2ref+εv=i2ref+[KPvΔvTF+KIv∫ΔvTFdt]ΔvTF=vTF*−vTF (15) 式中:KPv为输出电压比例补偿系数,KIv为输出电压积分补偿系数,vTF*为输入指令电压。
设负载电流参考值为iTFref。超级电容储能系统优化目标:(1)根据电池储能系统输出,补充负载期望功率,使负载电流快速上升,并在平顶段保持在iTFref;(2)吸收或释放瞬时功率,降低输出电压波动。由电池储能系统MPC得到其实时最优控制信号G1(k),由于T1≥N2T2,可在N2步预测中认为G1(k)不变,由此得到i2在每一步的预测值。
为实现对负载电流的控制,建立目标函数J21如下
J21=N2∑j=1[i2(k+j)+i4(k+j)−iTFref]2 (16) 由于负载电感较大,当输出电感电流变化较大时,容易引起输出电压的剧烈变化,不利于储能系统的正常工作,因此需要限制输出电感电流的控制效果,即
J22=N2∑j=1[i2(k+j)+i4(k+j)−iTF(k+j−1)]2 (17) 得到超级电容储能系统MPC控制目标函数为
J2=β1N2∑j=1[i2(k+j)+i4(k+j)−iTFref]2+β2N2∑j=1[i2(k+j)+i4(k+j)−iTF(k+j−1)]2 (18) 式中:β1、β2为权重系数,且β1+β2=1。式(16)和式(17)的控制效果存在冲突,需要在不同放电阶段设定权重系数β1、β2的值。在电流上升阶段,需要超级电容组快速建立输出电压电流,因此在此阶段β1占主导;输出电流上升至参考值iTFref附近,需要稳定输出电压电流,限制控制效果,在此阶段适当增大β2的值。为防止超级电容电流过大,设定最大放电电流ISCmax,当i3(k) > ISCmax时,强制令G2(k)=0。
同样,电路工况受放电阶段影响,为提高超级电容储能系统对电流的补偿作用,在iTFref的基础上增加电流反馈补偿项εi,得到电流指令值iTFref为
{iTFref*=iTFref+εi=iTFref+(KPiΔiTF+KIi∫ΔiTFdt)ΔiTF=iTFref−iTF (19) 式中:KPi为输出电流比例补偿系数,KIi为输出电流积分补偿系数。
图4所示为混合储能系统控制原理图,其中,图4(a)为电池储能系统控制原理图,图4(b)为超级电容储能系统控制原理图。
3. 仿真分析
按照图3所示的系统原理图搭建仿真模型,仿真模型使用参数如表2所示。其中,rL、rC分别为Cuk电路电感i1~i4、电容C1~C3的寄生电阻;β1_up、β1_flat分别表示电流上升段、平顶段的权重系数β1。
表 2 仿真模型参数Table 2. Parameters of simulation modelvB/V RB/mΩ RC/mΩ C/F L2, L4/mH vC(0)/V SSOC_B(0) 100 0.102 0.955 26400 0.6 114 0.8 C1, C2/F C3/μF L1, L3/mH rC/mΩ rL/mΩ LTF/mH RTF/mΩ 10 470×20 0.5 1 0.2 32 5.6 N1, N2 T1/ms T2/ms i1ref/kA i2ref/kA ISCmax/kA iTFref/kA 2 0.25 0.1 12 10 10 14 α1 β1_up β1_flat KPi KIi KPv KIv 0.9 0.9 0.5 2 0 5 10 在Simulink中建立混合储能系统仿真模型。图5所示为TF线圈仿真波形,图5(a)为负载电压vTF波形,图5(b)为负载电流iTF波形,图5(c)为系统瞬时输出功率pTF波形,图5(d)为系统输出能量WTF波形。TF电流超调量为0.5%,在7 s的平顶段内电流纹波为0.22%,满足线圈电流纹波小于1%的需求。系统输出最大瞬时功率2.6 MV·A,输出能量12.2 MJ。
图6所示为电池储能系统电压电流波形。图6(a)为电池组端电流即输入电感电流i1波形,图6(b)为电池组端电压vin1波形,图6(c)为输出电感电流i2波形,图6(d)为中间电容电压v1波形。在0.15 s时,完成电池组输出电压、电流建立;平顶放电阶段,i1波动为1.67%,所设计的控制策略有效实现电池组稳定放电的目标。输出电感电流i2提供主要TF电流输出,具有对应于负载电压波动的脉动,同时起到稳定电压的作用;Cuk变换器输入输出端通过C1进行能量交换,C1对i2的脉动进行能量的吸收与释放,使电池组输出电压电流保持相对稳定。通过中间电容C1的能量缓冲作用降低负载侧功率变化对电池组的影响,这就是Cuk电路应用于电池储能的优势。
图7所示为超级电容储能系统电压电流波形。图7(a)为超级电容组端电流即输入电感电流i3波形,图7(b)为超级电容组端电压vin2波形,图7(c)为输出电感电流i4波形,图7(d)为中间电容电压v2波形。仿真开始后,输入端电流快速上升,在0.12 s到达峰值;在2.5 s之前,超级电容组以上限电流ISCmax放电,实现负载电流快速上升;之后,系统进入平顶放电阶段,超级电容储能系统开始工作于补充负载电流、吸收或释放强感性负载导致的功率波动的工作状态。
图8为储能装置荷电状态波形。图8(a)为电池组SSOC_B曲线,图8(b)为超级电容组SSOC_SC曲线。图中表明,电池组实现平滑放电,超级电容组吸收、释放电路中的功率波动,混合储能系统功率分配合理。
为体现本文提出控制算法的优势,将使用MPC与使用闭环PI控制得到的输出电流iTF进行对比,如图9所示,从图中可以看出MPC方法使输出电流具有更短的上升时间,平顶段输出电流iTF符合要求,且能有效补偿因器件杂散参数带来的输出偏差。
4. 结 论
针对HL-3装置TF线圈这类强感性、大电流、高功率、大容量能量输出的负载,基于Cuk电路拓扑和模型预测控制方法,设计电池-超级电容混合储能系统的放电控制策略。仿真实验结果表明,所提出的控制策略实现电池组稳定放电,避免电池组在放电过程中出现电压电流放电尖峰,并且充分利用超级电容高功率密度、强耐充放电能力特性,吸收或释放负载端因强感性造成的瞬时功率,使储能系统输出满足负载需求的7 s平顶段14 kA稳定电流,电流纹波为0.22%。所设计的控制方法有效达到所有控制目标,具有一定的参考价值。
-
表 1 TF线圈参数
Table 1. Parameters of TF coil
LTF/mH RTF/mΩ tflat/s iTF_flat/kA 32 5.6 7 140 表 2 仿真模型参数
Table 2. Parameters of simulation model
vB/V RB/mΩ RC/mΩ C/F L2, L4/mH vC(0)/V SSOC_B(0) 100 0.102 0.955 26400 0.6 114 0.8 C1, C2/F C3/μF L1, L3/mH rC/mΩ rL/mΩ LTF/mH RTF/mΩ 10 470×20 0.5 1 0.2 32 5.6 N1, N2 T1/ms T2/ms i1ref/kA i2ref/kA ISCmax/kA iTFref/kA 2 0.25 0.1 12 10 10 14 α1 β1_up β1_flat KPi KIi KPv KIv 0.9 0.9 0.5 2 0 5 10 -
[1] Zhong Wulyu, HL-3 Team. China's HL-3 tokamak achieves H-mode operation with 1 MA plasma current[J]. The Innovation, 2024, 5: 100555. [2] 宣伟民, 王英翘, 李华俊, 等. HL-2M装置供电系统研制[J]. 核聚变与等离子体物理, 2021, 41(s2):437-442Xuan Weimin, Wang Yingqiao, Li Huajun, et al. Development of HL-2M power supply system[J]. Nuclear Fusion and Plasma Physics, 2021, 41(s2): 437-442 [3] 袁保山, 龙永兴, 邱银, 等. HL-2M环向场线圈电磁场和受力的计算分析[J]. 核聚变与等离子体物理, 2017, 37(3):249-256Yuan Baoshan, Long Yongxing, Qiu Yin, et al. Electromagnetic calculation and load analysis of the HL-2M TF coils[J]. Nuclear Fusion and Plasma Physics, 2017, 37(3): 249-256 [4] 李逢兵. 含锂电池和超级电容混合储能系统的控制与优化研究[D]. 重庆: 重庆大学, 2015Li Fengbing. Control and optimization of hybrid energy storage systems containing lithium-ion batteries and super-capacitors[D]. Chongqing: Chongqing University, 2015 [5] Kollimalla S K, Mishra M K, Narasamma N L. Design and analysis of novel control strategy for battery and supercapacitor storage system[J]. IEEE Transactions on Sustainable Energy, 2014, 5(4): 1137-1144. doi: 10.1109/TSTE.2014.2336896 [6] 王红艳, 张文倩. 改进型逻辑门限混合储能系统控制策略研究[J]. 智慧电力, 2020, 48(5):41-46 doi: 10.3969/j.issn.1673-7598.2020.05.007Wang Hongyan, Zhang Wenqian. Control strategy of improved logic threshold hybrid energy storage system[J]. Smart Power, 2020, 48(5): 41-46 doi: 10.3969/j.issn.1673-7598.2020.05.007 [7] Essoufi M, Hajji B, Rabhi A. Fuzzy logic based energy management strategy for fuel cell hybrid electric vehicle[C]//Proceedings of 2020 International Conference on Electrical and Information Technologies. 2020: 1-7. [8] 吴铁洲, 薛竹山, 向富超, 等. 基于遗传算法的电动汽车HESS功率优化分配[J]. 电气传动, 2019, 49(2):49-55Wu Tiezhou, Xue Zhushan, Xiang Fuchao, et al. Optimization of HESS power allocation for electric vehicle based on genetic algorithm[J]. Electric Drive, 2019, 49(2): 49-55 [9] Javed M S, Ma Tao, Jurasz J, et al. Economic analysis and optimization of a renewable energy based power supply system with different energy storages for a remote island[J]. Renewable Energy, 2021, 164: 1376-1394. doi: 10.1016/j.renene.2020.10.063 [10] Khoobi Arani S, Niasar A H, Zadeh A H. Energy management of dual-source propelled electric vehicle using fuzzy controller optimized via genetic algorithm[C]//Proceedings of the 7th Power Electronics and Drive Systems Technologies Conference. 2016: 338-343. [11] Raveendran Nair U, Costa-Castelló R. A model predictive control-based energy management scheme for hybrid storage system in islanded microgrids[J]. IEEE Access, 2020, 8: 97809-97822. doi: 10.1109/ACCESS.2020.2996434 [12] Xiong Rui, Cao Jiayi, Yu Quanqing. Reinforcement learning-based real-time power management for hybrid energy storage system in the plug-in hybrid electric vehicle[J]. Applied Energy, 2018, 211: 538-548. doi: 10.1016/j.apenergy.2017.11.072 [13] Xiong Rui, Chen Huan, Wang Chun, et al. Towards a smarter hybrid energy storage system based on battery and ultracapacitor - A critical review on topology and energy management[J]. Journal of Cleaner Production, 2018, 202: 1228-1240. doi: 10.1016/j.jclepro.2018.08.134 [14] Santucci A, Sorniotti A, Lekakou C. Power split strategies for hybrid energy storage systems for vehicular applications[J]. Journal of Power Sources, 2014, 258: 395-407. doi: 10.1016/j.jpowsour.2014.01.118 [15] Kumar K, Bae S. Coordinated dynamic power management for renewable energy-based grid-connected microgrids using model predictive control[J]. IEEE Transactions on Industrial Informatics, 2023, 19(9): 9596-9608. doi: 10.1109/TII.2022.3231409 [16] 宋元明, 刘亚杰, 金光, 等. 锂离子电池/超级电容器混合储能系统能量管理方法综述[J]. 储能科学与技术, 2024, 13(2):652-668Song Yuanming, Liu Yajie, Jin Guang, et al. Review of energy management methods for lithium-ion battery/supercapacitor hybrid energy storage systems[J]. Energy Storage Science and Technology, 2024, 13(2): 652-668 [17] Song Ziyou, Hofmann H, Li Jiaqiu, et al. A comparison study of different semi-active hybrid energy storage system topologies for electric vehicles[J]. Journal of Power Sources, 2015, 274: 400-411. doi: 10.1016/j.jpowsour.2014.10.061 [18] Gorji S A, Sahebi H G, Ektesabi M, et al. Topologies and control schemes of bidirectional DC-DC power converters: an overview[J]. IEEE Access, 2019, 7: 117997-118019. doi: 10.1109/ACCESS.2019.2937239 [19] Jaga O P, Gupta R, Jena B, et al. Bi-directional DC/DC converters used in interfacing ESSs for RESs and EVs: a review[J]. IETE Technical Review, 2023, 40(3): 334-370. doi: 10.1080/02564602.2022.2116362 -