梁建英 ,陈维荣
(1.中车青岛四方机车车辆股份有限公司,山东 青岛 266000;2.西南交通大学电气工程学院,四川 成都 610031)
随着社会对城市内景观要求的不断提高,且传统的化石燃料也在开始逐渐枯竭,开发各种可再生能源开始受到各国的支持[1], 其中,氢能的应用前景最为广阔[2].在此背景下,燃料电池随之出现在众多研究者的眼中,作为可将氢中的化学能转化为电能的潜力发电装置备受到国内外研究者喜爱[3].在不同的燃料电池类型中,质子交换膜燃料电池 (proton exchange membrane fuel cell,PEMFC)具有低温运行和启动快等优势[4].因此,PEMFC系统在分布式发电、新能源汽车移动设备等领域得到广泛运用.
如今,需要用大功率驱动的交通车辆,如现代城市的有轨电车和高速列车,因其舒适、快捷、方便等优点,作为一种省与省间的常用交通运输工具,将PEMFC发电系统应用于这类运输领域是一种具有巨大发展潜力的新型运营方式[5].但是由于单堆技术、材料等无法满足大功率列车长时间运行的缺点,导致燃料电池 (fuel cell, FC)发电系统在大功率机车行业的大规模应用成为巨大的挑战.不过,使用多堆燃料电池系统 (multi-stack fuel cell system,MFCS)可以一定程度上地代替大功率PEMFC的开发.MFCS是由几套低功率FC系统组成,而不是只由一个高功率的堆单独构成,这使得对FC的替换更加方便[6],从而大大提高了系统的稳定性和耐久性.但是,如今所用PEMFC是一个强耦合的多物理系统,其电气性能取决于其性能退化程度和运行参数,如气体湿度、压力和温度.虽然有大量的研究人员正在研究电磁系统,但这些现有策略的主要缺点是将燃料电池视为静态模型,没有考虑燃料电池的性能会随着运行状态而变化.因此,有必要在MFCS控制中考虑这些参数变化对FC运行效率区间变化的影响.在文献[7]中,提出了自适应递归最小二乘 (adaptive recursive least squares,ARLS)算法来识别燃料电池的性能,以找到其最大效率和功率点.文献[8]提出了一种遗忘因子递推最小二乘 (forgetting factor recursive least squares,FFRLS)在线辨识方法,以提高在变化负载下光纤通道参数估计精度.结果表明,该在线辨识方法可以在不考虑多物理量波动的情况下就估计出光纤通道参数.文献[9]提出了一种效率优化策略,在双堆FC系统中进行功率分配,提高FC运行时长,仿真结果验证了该策略的有效性.此外,为了实现串联式MFCS的最优功率分配,文献[10]提出了一种最大功率跟踪方法.
综上所述,需考虑PEMFC的一个“最大效率范围”(maximum efficiency range,MER),在这个范围内控制PEMFC功率输出,有助于优化系统效率和PEMFC性能[11].此外,根据前面的描述,该MER会随着外部参数(运行条件、温度、湿度等)的变化而移动[12].然而,很少有研究将MFCS中的功率分配与单个PEMFC的MER识别方法联系起来共同优化FC系统效率.因此,本文提出了一种基于在线辨识的MFCS改进型链式功率分配方法,用以提高FC系统效率和延长堆的使用寿命,促进FC在高功率场合的大规模应用.此外,为说明此方法的优势所在,本文将所提方法与平均功率分配方法和传统链式功率分配方法做了对比分析.最后,通过搭建的RT-LAB半实物平台,测试了所搭建的MFCS的性能以及所提出的基于在线辨识的改进型链式功率分配方法的有效性和实用性.
由于并联拓扑结构在系统复杂度和系统性能方面都优于串联拓扑结构,因此,本文构建的MFCS采用前一种拓扑,如图1所示.本研究开发的大功率MFCS由3套独立的PEMFC系统、控制器以及单向DC/DC变换器组成.所用3个PEMFC虽然是同一型号,但由于使用时间不同性能上有一定差异,其主要参数如表1中所示,3个PEMFC的极化曲线如图2所示.
表1 PEMFC电堆参数Tab.1 PEMFC stack parameters
图1 多堆燃料电池系统结构Fig.1 Multi-stack fuel cell system structure
建立PEMFC仿真模型,如式(1)所示[13].
式中:UFC为电堆输出电压;Eoc为开路电压;Uact为极化电压降;Uohmic为欧姆电压降.
这些参数可进一步表示[13]如下:
式(2)、(3)中:Kc为额定电压参数;T为 PEMFC的运行温度;E0为标准压力下的电动势;z为电子转移数;F为法拉第常数;R为气体常数;PH2和PO2分别为氢气和氧气压力;N为单体片数;RFC为PEMFC内阻;iFC为PEMFC模型测试电流;A和i0分别为埃菲尔斜率和交流电流,如式(4)所示.
式中:k为玻尔兹曼常数;h普朗克常数;-ΔG为吉布斯能量;α为转移系数.
根据建立的模型和实验测试数据校准相应参数[14],并将实验测试数据和模型仿真结果对比分析,如图2所示,本文所建立模型能很好地反映PEMFC的极化曲线变化.
图2 燃料电池极化曲线Fig.2 Fuel cell polarization curve
本研究的主要目的是优化燃料电池系统的效率,使其各电堆尽可能在MER中运行.因此,PEMFC系统的效率是需要分析的重要指标.
PEMFC效率定义为电池组产生的电能与供应的氢气化学能之间的比值.PEMFC的效率计算如式(5)所示[15].
式中:PSFC为单片电池产生的功率;Vcell为PEMFC单片电压;IFC为PEMFC输出电流; ΔHHV为氢气高热值.
为确保PEMFC系统正常运行,一般需要一些辅助设备例如单向DC/DC变换器、控制器、燃料电池辅机等,这些均会影响燃料电池系统效率[16].因此,每个PEMFC系统的效率 ηsys可由式(6)计算得到.
式中:ηaux为所有辅助设备的效率;ηFC为氢气利用率,通常该值为 99% ~ 100%[17],本文中ηFC的值为100%;Paux为PEMFC控制器消耗的功率、DC/DC变换器消耗功率以及DSP控制器消耗功率的总和.
根据式(4)可得PEMFC系统的效率曲线如图3所示,图中:PFCj_L和PFCj_H分别为PEMFCj的低功率值和高功率值,j= 1,2,3.
如图3所示,依据所用3个不同使用时长的堆所测效率/功率曲线数据,将效率高于EFF1、EFF2以及EFF3的区域定义为各自的高效率区域,而相对应的功率区间即为本文研究中的MER.此外,根据前面的描述,MER将随着外部条件的改变而移动,图3中各个 PEMFC的边界功率值(PFCj_L和PFCj_H)会随着堆长期的运行而发生变化.因此,为更好地动态调节PEMFC高效运行区域,使用FFRLS在线参数识别算法来估计MER的边界功率值,使各PEMFC系统尽可能在MER中运行,进而实现延长堆的使用寿命,提高MFCS平均效率,减少氢气消耗.
图3 燃料电池系统效率-功率曲线Fig.3 Fuel cell system efficiency-power curves
FFRLS算法已经在许多参数识别应用中得到采用,该算法在一些研究中有详细描述[18].根据图3的曲线形式,本文是通过六阶多项式来描述PEMFC系统输出功率与效率的数学模型, 表达式如(7)所示.
式中:PFC为燃料电池输出功率矩阵;a0~a6为系统的未知参数.
对于单输入单输出(single input single output,SISO)系统,可以采用自回归移动平均(auto regressive moving average, ARMA)模型[19]来描述被控对象.模型可用式(8)表示.
式中:y和θ分别为系统输出和输入向量;ϕ为系统未知参数向量;e为辨识误差;k为离散时间.
FFRLS算法由误差e定义,性能标准用J表示,如式(9)所示.
式中:μ(0<μ< 1)为遗忘因子,用于增强对最近识别误差的敏感性,并降低前期估计误差的权重[20];L为数据迭代时间长度.
可以利用式(10)来实现FFRLS算法.
式中:P为预测的协方差矩阵;L为卡尔曼增益矩阵.
链式功率分配方法是一种顺序方法[9].提供功率的燃料电池数量取决于功率基准.只有当前PEMFC达到最大输出功率后,下一个PEMFC才能继续提供剩余的功率,直到达到最大功率.该过程一直持续到所有燃料电池都被使用或达到所需功率为止,表达式如式(11)所示.
式中:Pload为系统负载功率;Pmax为燃料电池最大功率.
本文为提高MFCS效率,提出了一种基于在线辨识的改进型链式功率分配方法.该方法利用FFRLS在线辨识算法实时获取MER的值(PFCj_L和PFCj_H),并将其作为链式功率输出约束条件,以优化PEMFC运行过程中效率不佳和性能差异造成的功率分配问题,增加性能好的堆出力,减少性能差的堆出力以及输出功率波动以便提升PEMFC的耐久性.该算法具体执行过程如图4所示.
图4 改进链式功率分配算法结构Fig.4 Structure of improved chain power distribution algorithm
加拿大Opal-RT公司开发的RT-LAB 实时仿真系统是专门针对电力电子、电力系统、电力拖动等系统实时仿真的运行平台,可以对基于MATLAB的Simulink模型进行实时仿真.本文利用RT-LAB的半实物平台对Simulink模型中的模型和控制器中的算法进行测试,首先将在线辨识和功率分配算法写入DSP中,在RT-LAB中搭建仿真模型,最后在DSP控制下完成实时在线验证.整个系统如图5所示.
图5 半实物平台Fig.5 Semi-physical platform
功率分配方法是在组成系统的不同PEMFC之间分配参考功率,则
式中:PFCj为PEMFCj的功率.
分配系数K1~K3由式(13)引入和定义.
通过在3个PEMFC之间平均分配参考功率来达到供应负载需求的目的,则
按照式(14)的分配策略,其试验结果如图6、图7所示.从图6中可看出:所有PEMFC都具有相同的功率需求,说明本文所用的平均功率分配方法验证有效.但由于3个PEMFC性能上有一定差异,全部以相同功率输出最终会造成PEMFC3先停运,不利于整个MFCS稳定运行,并且3个PEMFC效率特性曲线不同,以相同功率长期运行也不利于提升MFCS整体效率.从图7可知:此时MFCS平均效率为48.98%.
图6 平均功率分配方法Fig.6 Average power distribution method
图7 平均功率分配的多堆燃料电池系统效率Fig.7 Efficiency of multi-stack fuel cell system for average power distribution
传统链式功率分配是依据任意MFCS中PEMFC的排列顺序,依次启动.当第1个PEMFC输出功率达到最大功率时,第2个PEMFC补充输出,直至满足负载总需求功率为止.各个PEMFC输出功率和MFCS效率的实验测试结果如图8、9所示.
图8 传统链式功率分配方法Fig.8 Traditional chain power distribution method
从图8中不难看出:MFCS系统出力的方式是PEMFC1先输出功率直到输出达到限幅值,然后PEMFC2和PEMFC3依次接力输出负载剩余需求功率,说明所用传统链式功率分配方法验证有效.但是,结合图3对各PEMFC的效率/功率曲线的分析,这种传统链式功率分配方法没有对PEMFC的运行区间约束,以至于PEMFC1、PEMFC2和PEMFC3长期在最大或最低限幅功率处运行,导致所有PEMFC效率都较低,从图9中可知:MFCS平均效率只有47.96%,造成MFCS运行成本上升.
图9 传统链式功率的多堆燃料电池系统效率Fig.9 Efficiency of multi-stack fuel cell system for traditional chain power distribution
4.4.1 FFRLS在线辨识
通过FFRLS在线辨识算法确定MER是最关键的步骤,本文利用FFRLS强大的在线拟合能力,通过获取每个PEMFC系统运行时的效率/功率曲线,然后利用高阶方程求解设定目标效率的对应功率解.具体实验结果如图10所示.
图10 PEMFC系统效率和功率Fig.10 PEMFC system efficiency and power
4.4.2 基于FFRLS在线辨识的链式功率分配
改进链式功率分配方法是利用FFRLS在线辨识算法实时数据更新约束值,具体实验结果如图11、图12所示.
图11 基于在线辨识的改进链式功率分配方法Fig.11 Improved chain power distribution method based on online identification
从图11中不难看出:该方法能有效控制PEMFC的效率区间,保证单体PEMFC系统尽量在最优效率区间运行,不仅可以提升MFCS效率,减少系统氢气消耗,还有利于减缓性能较差的PEMFC老化速率,增强了系统稳定性.从图12可知:此时MFCS平均效率为 49.91%.
图12 多堆燃料电池系统效率Fig.12 Efficiency of multi-stack fuel cell system for Improved chain power distribution
从图7、9、12的对比中不难看出:本文所提出的基于在线辨识的改进链式功率分配方法相比平均功率分配方法和传统链式功率分配方法而言,对MFCS效率分别提升了0.93%和1.95%,不仅提高了MFCS运行效率,减少了系统氢耗,也一定程度上优化了MFCS中每个PEMFC运行区间的功率波动,减缓了PEMFC老化速率.具体MFCS效率变化数据比较结果如表2所示.
表2 3种功率分配方法下MFCS效率变化Tab.2 Changes in MFCS efficiency under three power distribution methods %
为优化各个PEMFC发电系统运行效率区间,提高MFCS平均效率延长PEMFC使用寿命,本文在分析MFCS功率分配方法的基础上,提出了一种基于FFRLS的在线辨识算法对各个PEMFC系统的MER寻找,并将其作为功率分配的实时约束条件,以优化单个PEMFC系统运行效率区间.为验证所提方法的有效性和实用性,本文利用RT-LAB搭建了快速原型系统,并在测试平台上完成了功率分配方法和效率区间优化相结合的验证和分析,结论如下:
1) 利用FFRLS在线辨识算法实现了对系统效率/功率曲线的在线辨识.
2) 通过对高阶函数的实时求解,估算出了运行中PEMFC的MER,确定出了链式功率分配策略的边界区间.
3) 对比分析了平均功率分配和传统链式功率分配方法,本文所提方法弥补了PEMFC性能不同时功率分配的不合理和MFCS平均效率较低的缺点.