采用NSGA-II算法的车载复合电源参数匹配
周放1, 宋传学1, 梁天唯2, 肖峰1
1.吉林大学 汽车仿真与控制国家重点实验室,长春 130022
2.吉林大学 机械科学与工程学院,长春 130022
通信作者:肖峰(1986-),男,助理研究员,博士.研究方向:汽车系统动力学.E-mail:xiaofeng4206@126.com

作者简介:周放(1988-),男,博士研究生.研究方向:车载复合电源关键技术.E-mail:zhoufang13@mails.jlu.edu.cn

摘要

提出了一种主动并联式车载复合电源的参数匹配/能量管理联合优化方法。首先,根据整车性能指标确定电池组和超级电容组的容量范围,并针对部件性能制定能量管理策略,然后,以复合电源替换成本和日均电耗最小为目标建立多目标优化函数,采用受控精英多目标遗传算法(Controlled elitist NSGA-II)在不同的日均里程下对复合电源参数和能量管理策略同时进行优化,得到对应的Pareto解集,作为备选。综合考虑备选集的初始成本、百公里总成本和电池组寿命里程,优选电池组和超级电容组参数。优化结果表明:日均里程在48 km以内时,匹配22 kW·h电池组,超过48 km后,里程越大,匹配电池组容量越大;在不同的日均里程下,不同容量的电池组匹配一个60 W·h左右的超级电容组可以达到综合性能最优。与电池单独优化相比,电源替换成本平均降低11.5 %,寿命里程平均提高19.8 %,日均总成本平均降低13.2 %。

关键词: 车辆工程; 复合电源; 受控精英多目标遗传算法; 替换成本; 多目标优化
中图分类号:U461.8 文献标志码:A 文章编号:1671-5497(2017)05-1336-08
Parameter matching of on-board hybrid energy storage system using NSGA-II algorithm
ZHOU Fang1, SONG Chuan-xue1, LIANG Tian-wei2, XIAO Feng1
1.State Key Laboratory of Automotive Simulation and Control, Jilin University, Changchun 130022, China
2.College of Mechanical Science and Engineering, Jilin University, Changchun 130022, China
Abstract

An integrated parameter matching/energy management optimization method is proposed for active on-board Hybrid Energy Storage System (HESS). The optimization model with two indicators, HESS replacement costs and energy consumption, was set up and solved using controlled elitist NSGA-Ⅱalgorithm to obtain Pareto sets as alternatives. After comparing the initial cost, daily total cost and battery pack mileage of lifetime of the alternatives, HESS parameters and corresponding energy management were selected. When the daily range is within 48 km, the optimal battery pack capacity is 22 kW·h, when the daily range exceeds 48 km, the optimal battery pack capacity increases with the daily range. An ultracapacitor pack around 60 W·h is adequate for various daily range. Compared with the optimization results of battery only, using the proposed optimization method, the replacement cost is reduced by 11.5% in average, the battery pack mileage of lifetime is raised by 19.8% in average, the daily total cost is reduced by 13.2%. The proposed optimization method ensures that the main indicators are optimal and the secondary indicators are satisfactory.

Key words: vehicle engineering; hybrid energy storage system(HESS); controlled elitist NSGA-II; replacement cost; multi-objective optimization
0 引 言

电动汽车能够降低燃油消耗, 并且具有零排放和环境友好的特点, 是应对能源和环境危机的有效解决方案。但电池作为电动汽车动力源往往存在功率密度低和循环寿命短的问题, 在车辆行驶过程中的峰值功率需求和瞬态功率会加速电池组的老化, 从而增加了电源在质保期内的替换成本。超级电容具有功率密度高和循环寿命长的优点, 由超级电容和电池组成的复合电源可以克服上述缺点, 充分发挥两者在功率密度和能量密度方面的优势。

超级电容的加入旨在延长电池组寿命、降低成本, 同时保证经济性, 参数匹配和能量管理策略制定是以上目标能否达到预期的核心问题, 也是近年来的研究热点。Rui等[1]以复合电源成本和使用成本加权和最小为目标对电源参数进行优化匹配。王庆年等[2]基于目标工况对功率需求进行解耦, 并以经济性、效率、寿命和成本为参考, 对HEV的发动机、电池和超级电容进行了匹配。Hu等[3]以复合电源初始成本和使用成本加权和为优化目标, 采用凸优化对电源参数和功率分配进行优化。Shen等[4]以复合电源质量和电池工况衰减加权最低为优化目标, 采用DIRECT算法对复合电源参数和控制策略参数进行优化。Song等[5]以复合电源成本和工况下电池容量衰减最小为目标, 采用多次动态规划对超级电容参数和功率分配策略进行了优化。

本文以车载主动并联式复合电源构型为研究对象, 以整车性能为约束, 以复合电源的替换成本和车辆日均电耗为优化目标建立优化函数, 针对不同的日均里程工况, 采用受控精英快速非支配遗传算法(Controlled elitist NSGA-II)进行优化求解, 并根据电池组寿命里程等性能指标对求解结果进一步优选, 进行复合电源的多目标优化匹配设计。

1 复合电源构型选择

复合电源构型包括直接并联和主动并联。直接并联结构中电流不可控; 电池组主动并联结构中, 电机控制器直接和超级电容组连接, 电压范围波动大, 控制难度增加[6]。超级电容组主动并联结构中, 母线电压波动范围小, DC/DC工作时间短, 负荷较小, 不易出现发热问题, 选为本文目标构型, 如图1所示。

图1 复合电源构型Fig.1 HESS topology

2 参数匹配目标

复合电源的构型确定后, 需要确定电池组串、并联数nbsnbp, 超级电容组串、并联数nucsnucp, 以及DC/DC的功率等级PDC。复合电源的参数优化匹配任务就是在满足整车性能需求的前提下, 寻找一组参数集X=[ nbs, nbp, nucs, nucp, PDC], 使目标性能指标达到最优[7]

本文选择复合电源替换成本和整车经济性作为待优化指标, 建立多目标优化框架, 寻找最优参数集X

3 参数优化匹配流程

用于匹配优化的电池类型选取松下18650B, 超级电容选取Maxwell的BCAP 3000, 单体参数如表1所示。

电机峰值功率55 kW, 峰值转矩190 N· m, 最大转速为9000 r/min, 电压等级312 V, DC/DC功率不小于电机峰值功率, PDC定为55 kW。超级电容单体容量较大, 并联数选为1, 即nucp为1。

表1 单体参数 Table 1 Cell parameter
3.1 容量范围

车载电源需要满足车辆能量需求, 受续驶里程限制。电池组容量应满足车辆在许用荷电状态(SOC)范围内达到续驶里程, 能量需求为:

EbE100L100(SOCH-SOCL)(1)

式中:E100=12.5 kW· h, 为车辆百公里电耗; L=120 km, 为续驶里程; SOCH=90%, 为电池组SOC上限; SOCL=20%, 为电池SOC下限。

根据式(1)计算得到:电池组能量不小于21.4 kW· h; nbp≥ 21.4× 1000/nbs/Q=21.19, 其中Q为电池单体容量, 11.88 W· h; nbp不小于22。

整车能量需求分为循环能量和耗散能量[8], 其中循环能量包括车辆的动能和势能, 这部分能量可以通过再生制动回收, 并在超级电容组中循环, 超级电容应该满足循环流动的能量需求。图2为UDDS工况下车辆的循环功率和耗散功率需求, 其中循环功率在子工况(相邻停车之间的工况)内积分为0。

图2 功率需求分离Fig.2 Power demand seperation

定义车辆制动回收能量与循环能量的比值为制动效率, 电池组提供平均充放电功率, 剩余功率由超级电容承担, 据此求出超级电容组的充电和放电总能量, 然后求出每个子工况下的超级电容组平均能量需求, 最后根据许用电压范围以及超级电容组和DC/DC效率计算超级电容组的容量范围:

EUCtcyclcPUC(t)dtnsubcycle(UH2-UL2Umax2)ηUCηDC(2)

式中:PUC为超级电容功率, kW; nsubcycle为子工况数; UHUL分别为超级电容电压上限和下限, V; Umax为超级电容最大电压, V; η UC为超级电容组效率; η DC为DC/DC平均效率。

根据式(2)计算得到:超级电容组容量不小于44.3 W· h, nucs不小于15。

3.2 能量管理策略

本文制定了逻辑门限策略并对其参数进行优化, 以保证结果能够近似最优且能实际应用。

控制逻辑如下:再生制动功率不超过PUC, limit时, 电池充电, 超过后, 超级电容充电, 同时电池为超级电容以恒功率Pch供电; 功率需求不超过Plow时, 电池单独供电; 功率需求在PminPhigh之间时, 电池以Pmin供电, 超级电容补足剩余需求; 功率需求超过Phigh之后, 电池以Pmax供电, 超级电容补足剩余需求; 功率需求为正, 如果USOC小于0.5, 则电池单独供电, 并且向超级电容充电; 在电池功率输出上加入一阶低通滤波防止振荡, 滤波常数为tau。控制策略共有6个待定参数:ParaEMS=[PUC, limit, Pmin, Pmax, Phigh, Pch, tau]。

再生制动策略表述如下:在满足ECE法规和I线限制的要求下, 前、后制动力分配尽量使前轮制动力最大; 机电制动分配中, 电机转速低于160 r/min或电池组SOC大于0.9时, 电机不回收能量[9]

3.3 行驶-充电工况

以UDDS工况为目标工况, 根据NHTS统计数据[10], 选择12、24、36、48、72和96 km为不同的日均里程, 并与充电工况结合, 组成行驶消耗-充电工况。工况初始SOC为0.9, 充电方式为每天晚上充电, 交流电220 V, 30 A, SOC达到0.9时, 充电停止。

3.4 多目标优化

3.4.1 优化模型

优化目标是质保期内的替换成本Ctotal和日均电耗Etotal, 优化模型表示为:

F(x)=min{Ctotal(x), Etotal(x)}x=[nbp, nucs, ParaEMS]x[xL, xH](3)

式中:x为优化参数集; xLxH分别为优化参数集下限和上限。

3.4.2 复合电源替换成本

复合电源的替换成本定义为复合电源各部件的初始成本和质保期内替换次数的乘积。其中DC/DC和超级电容组循环寿命很长, 在质保期内不用替换, 只计算电池组替换成本:

Ctotal=CUCnucs+CDCPDC+nreplaceCBnbsnbp(4)

式中:CUC为超级电容单体价格, 万元; CDC为DC/DC单位功率价格, 万元; CB为电池单体价格, 万元; nreplace为电池组替换次数。

根据工信部《关于免征新能源汽车车辆购置税的公告》规定[11], 对新能源汽车动力电池、电机、电控等关键零部件提供不低于5年或10万公里(以先到者为准)质保, 因此, 对车辆行驶10万公里的复合电源替换成本进行计算。

当电池组中单体可用容量衰减至初值的80%时, 电池组进行替换, 用加权电荷流通量表示电池容量衰减, 当单体的加权电荷流通量达到标准工况下的可用安时流通量时, 电池组进行替换, nreplace根据下式计算:

nreplace=100000Ah, eff, dayLdaynbpAh, life(5)

式中: Ah, eff, day为电池组加权安时流通量, A· h; Ah, life=1492 A· h, 为电池安时寿命; Lday为日均里程, km。

根据式(6)计算加权安时[12, 13], 包括循环和充电工况:

Ah, eff, day=tcyclew(t)i(t)dt+tchargew(t)i(t)dt(6)

式中:w(t)为老化因子; i(t)为电池实时电流, A。

老化因子代表温度和工况变化引起的老化加速速率, 根据电池组实时温度和电流查表得出, 温度根据均质对流换热模型计算得出[14]:

Q·bat=kAbat(Tbat-Tamb)+cbatmbat dTbatdt(7)

式中:Qbat为电池产热率(内阻消耗功率), W; k为对流传热系数; Abat为电池表面积, m2; Tbat为电池温度, ℃; Tamb为环境温度, ℃; cbat为电池比热容; mbat为电池质量, kg。

3.4.3 复合电源日均电耗

对行驶消耗工况下的电池组和超级电容组功率积分, 计算日均电耗, 如下式所示:

Etotal=PB(t)dt+PUC(t)dt(8)

式中:PB(t)为电池组功率, kW。

3.4.4 优化算法

上述优化问题具有多目标、多变量以及高度非线性的特点, 多目标遗传算法非常适合于求解此类问题。Controlled elitist NSGA-II算法在NSGA-II算法的基础上对精英保留策略进行了改进。受控精英保留策略保证每次迭代中精英个体和非精英个体按一定的比例保留下来, 改善了种群在非支配集垂向的探索能力, 提高了全局收敛性[15, 16]

受控精英保留策略表示如下:在通过精英策略从合并种群Rt选择新种群Pt+1个体时, 第i个前沿Fi允许的个体数量ni和第i-1个前沿Fi-1允许的个体数量ni-1有下面的关系:

ni=rni-19

式中:r为几何分配参数。

假设Rt(个体数量为2N)共有K个前沿, ni的计算式为:

ni=N1-r1-rKri-1(10)

通过几何分配, 新种群Pt+1在每个Fi(i=1, 2, …, K)上都以一定比例分布了个体, 保证了解集多样性, 提高了优化向最优解集收敛的能力。

Controlled elitist NSGA-II算法参数如表2所示, 流程及Simulink模型调用步骤如图3所示。

表2 NSGA-II参数设置 Table 2 NSGA-II parameters

图3 NSGA-II流程Fig.3 NSGA-II process

3.5 优化结果分析

3.5.1 Pareto解集

优化结果用非支配集(Pareto集)表示, 48 km日均里程对应的Pareto集如图4所示, 图中每个点都表示一组备选方案。替换成本和日均电耗呈现了明显的冲突, 两者无法同时达到最优。由于降低替换成本倾向于延长电池寿命, 电池的使用被最大限度的限制, 超级电容组应对更多的电流负荷。由于超级电容组的低压大电流和DC/DC的高损耗, 总体增加了电耗量, 从而引起了两个目标之间的冲突。

图4 48 km日均里程的Pareto集Fig.4 Pareto sets with 48 km daily range

Pareto集中各个方案之间无法量化地比较哪个更好, 只能根据经验选择。当倾向于提高经济性时, 选择备选集中右侧的参数; 当倾向于降低电源替换成本时, 选择备选集中左侧参数。

3.5.2 备选方案优选

进一步考察备选方案的电源初始成本Cinit、日均总成本Ctotal, day(日均替换成本和电耗成本的和)以及池组寿命里程(电池组寿命终结时的车辆总里程)LB, life三者的大小:

Cinit=CUCnucs+CDCPDC+CBnbsnbpCtotal, day=Ctotal100000/Lday+PelecEtotalLB, life=10nreplace(11)

式中:Pelec为电价, 元/(kW· h)。

图5(a)(b)(c)分别为初始成本、日均总成本和电池组里程寿命与替换成本的关系。初始成本直接影响车辆售价, 应该选择尽可能低的方案, 图5(a)明显分为两部分, 在左侧垂直部分的方案, 初始成本基本不变, 替换成本大幅度下降, 这是能量管理策略调控的结果, 右侧部分的方案, 随着初始成本增加, 电源容量增加, 替换成本缓慢下降; 日均总成本应选择尽可能小的方案; 电池组寿命里程显著影响消费者对整车的印象, 随初始成本增加而增加, 应该尽可能选择大的方案。由于初始成本低和寿命里程高两个目标相互冲突, 首先选择初始成本分布中拐点右侧部分作为备选方案, 然后根据尽量少替换电池组的原则, 在寿命里程大于5万公里的备选集中选择初始成本最低的方案。

图5 Pareto集成本、里程对比Fig.5 Pareto sets costs and life range comparison

同理, 从12、24、36、72和96 km日均里程的备选方案中优选出参数集, 列于表3中。可以看出, 日均里程在48 km以内时, 电池组最优匹配为22 kW· h左右, 日均里程超出48 km后, 随着里程增加, 电池组容量增加; 超级电容组容量基本维持在60 W· h左右。

表3 不同日均里程优选结果 Table 3 Optimal results with different daily range
4 对比分析

与电池组单独供电优化结果进行对比, 图6为替换成本对比。在不同的日均里程下, 复合电源相较于电池组替换成本减小2~3万元。随着日均里程的增加, 替换成本减少的比例增加, 这是因为日均里程越大, 电池单独工作时温度上升越高, 加速电池老化, 而复合电源中超级电容减缓了电池组的温升趋势, 且里程越长, 效果越明显, 进而缓解电池组老化趋势, 发挥更大的优势。

图6 替换成本对比Fig.6 Replacement costs comparison

图7为日均电耗对比, 复合电源日均电耗相较于电池组单独供电基本没有变化, 这说明超级电容组的低压和DC/DC的效率引起的电耗增加被最大限度的削弱, 优化在两个目标上都达到了令人满意的结果。

图7 日均电耗对比Fig.7 Daily power consumption

计算优化对应的电池组初始成本、寿命里程和日均总成本, 列于表4中。可以看出, 匹配方案在多个目标之间取得了良好的平衡, 主要性能指标替换成本平均降低11.5%, 寿命里程平均增加19.8%, 日均总成本平均降低13.2%。此外, 电池单独供能时电池组除了在日均里程48 km时替换一次, 其他里程下都替换两次, 而复合电源在所有里程下都只替换一次电池组。

根据对比结果, 所匹配的复合电源使车辆在多个性能指标上得到了优化, 多目标遗传算法在降低替换成本的同时保证了整车经济性, 基于Pareto集的优选方法既保证了主要性能指标最优, 也使其他性能指标达到理想结果, 证明了匹配方法的有效性。

表4 性能指标对比 Table 4 Performance comparison
5 结束语

提出了一种电动汽车复合电源的参数优化匹配方法。在性能约束下, 以复合电源替换成本和整车经济性为目标建立多目标优化函数, 并在不同的日均里程工况下对多目标问题进行求解, 得到Pareto集, 作为备选方案。然后, 以电源初始成本、日均总成本和电池组寿命里程为参考依据, 从备选方案中优选得到匹配结果。结果表明, 日均里程在50 km以内时, 匹配22 kW· h的电池组, 里程超过50 km时, 电池组容量随里程增加; 不同里程下不同容量电池组匹配60 W· h左右的超级电容组, 综合性能达到最优:替换成本平均降低11.5%, 寿命里程平均增加19.8%, 日均总成本平均降低13.2%。从成本角度来看, 复合电源与电池组单独供电相比具有明显的优势, 非常适合应用于电动汽车, 且日均里程越大, 复合电源的优势越大。后期控制策略的开发可以进一步改进复合电源的性能。

The authors have declared that no competing interests exist.

参考文献
[1] Rui E A, Castro R D, Pinto C, et al. Combined sizing and energy management in EVs with batteries and supercapacitors[J]. IEEE Transactions on Vehicular Technology, 2014, 63(7): 3062-3076. [本文引用:1]
[2] 王庆年, 曲晓冬, 于远彬. 基于目标工况的复合电源混合动力客车优化匹配[J]. 吉林大学学报: 工学版, 2013, 43(5): 1153-1159.
Wang Qing-nian, Qu Xiao-dong, Yu Yuan-bin. Optimal matching of HEV with composite electric power supply[J]. Journal of Jilin University(Engineering and Technology Edition), 2013, 43(5): 1153-1159. [本文引用:1]
[3] Hu X S, Murgovki N, Johannesson L M, et al. Comparison of three electrochemical energy buffers applied to a hybrid bus powertrain with simultaneous optimal sizing and energy management[J]. IEEE Transactions on Intelligent Transportation Systems, 2014, 15(3): 1193-1205. [本文引用:1]
[4] Shen J Y, Dusmez S, Khaligh A. Optimization of sizing and battery cycle life in battery/ultracapacitor hybrid energy storage systems for electric vehicle applications[J]. IEEE Transactions on Industrial Informatics, 2014, 10(4): 2112-2121. [本文引用:1]
[5] Song Z Y, Hofmann H, Li J Q, et al. Optimization for a hybrid energy storage system in electric vehicles using dynamic programming approach[J]. Applied Energy, 2015, 139: 151-162. [本文引用:1]
[6] 米勒. 超级电容的应用[M]. 韩晓娟, 李建林, 田春光, 译. 北京: 机械工业出版社, 2014: 97-98. [本文引用:1]
[7] 张抗抗, 徐梁飞, 华剑锋, . 基于多目标优化的纯电动车动力系统参数匹配方法[J]. 汽车工程, 2015(7): 757-765.
Zhang Kang-kang, Xu Liang-fei, Hua Jian-feng, et al. A parameter matching method for the powertrain of battery electric vehicle based on multi-objective optimization[J]. Automotive Engineering, 2015(7): 757-765. [本文引用:1]
[8] Ott T, Zurbriggen F, Onder C, et al. Cycle-averaged efficiency of hybrid electric vehicles[J]. Proceedings of the Institution of Mechanical Engineers Part D Journal of Automobile Engineering, 2013, 227(1): 78-86. [本文引用:1]
[9] 叶敏, 郭金刚. 电动汽车再生制动及其控制技术[M]. 北京: 人民交通出版社, 2013: 128-135. [本文引用:1]
[10] Department of Transportation. National Household Travel Survey[EB/OL]. [2016-04-02]. http://nhts.ornl.gov/. [本文引用:1]
[11] 中华人民共和国工业和信息化部. 关于免征新能源汽车车辆购置税的公告[EB/OL]. [2016-03-05]. http://www.miit.gov.cn/n1146295/n1652858/n1652930/n4509607/c4512240/content.html. [本文引用:1]
[12] Onori S. A new life estimation method for lithium-ion batteries in plug-in hybrid electric vehicles applications[J]. International Journal of Power Electronics, 2012, 4(3): 302-319. [本文引用:1]
[13] Marano V, Onori S, Guezennec Y, et al. Lithium-ion batteries life estimation for plug-in hybrid electric vehicles[C]//Vehicle Power and Propulsion Conference, IEEE, 2009: 536-543. [本文引用:1]
[14] Huria T, Ceraolo M, Gazzarri J, et al. High fidelity electrical model with thermal dependence for characterization and simulation of high power lithium battery cells[C]//Electric Vehicle Conference(IEVC), IEEE, 2012: 1-8. [本文引用:1]
[15] Deb K, Goel T. Controlled Elitist Non-dominated Sorting Genetic Algorithm for Better Convergence[M]. Berlin : Springer, 2001: 67-81. [本文引用:1]
[16] Deb K, Agrawal S, Pratap A, et al. A fast elitist non-dominated sorting genetic algorithm for multi-objective optimization: NSGA-II[C]//International Conference in Parallel Problem Solving from Nature. Berlin : Springer-Verlag, 2000: 849-858. [本文引用:]