一种锂离子电池寿命早期预测方法

1.本发明涉及锂离子电池技术领域，更具体地说，涉及一种锂离子电池寿命早期预测方法。


背景技术：

2.随着能源危机和生态环境问题的日趋严重，新能源汽车已成为汽车产业的发展重点。然而，新能源汽车上的电池包的性能直接影响着新能源汽车的性能和运行情况。新能源汽车上用的电池包是由很多电池单体通过并联和串联而组成，这些电池单体的寿命如果相差较大，将会导致电池包使用一段时间后出现输出性能下降，甚至出现安全风险。因此，在电池单体成组前对单体寿命进行早期预测十分重要。
3.目前几乎所有锂离子电池负极材料都是采用的石墨，这些锂离子电池脱锂负极活性材料的损失，会导致放电电压发生变化，而容量并没有变化。这是因为脱锂负极材料的损失会改变锂离子的储存电位，而不会改变总容量。在脱锂负极活性物质损失的高比率下，负极容量最终将低于剩余的锂离子库存。此时，负极在充电过程中将没有足够的位点容纳锂离子，从而导致锂镀层。由于镀锂是不可逆的，所以容量损失会加速。
4.目前尚缺乏一种有效的锂离子电池寿命早期预测方法。


技术实现要素：

5.本发明针对锂离子电池寿命早期预测问题，提出一种锂离子电池寿命早期预测方法，实现对电池寿命进行早期预测。由于锂离子电池早期充放电循环中，脱锂负极活性物质损失会在不影响容量衰减曲线的情况下改变电压曲线，因而，本发明利用这一特征可实现锂离子电池寿命早期预测。
6.本发明提供的技术方案如下：
7.一种锂离子电池寿命早期预测方法，包括以下步骤：
8.(1)模型训练数据的获取：对n颗锂离子电池单体进行编号(编号分布为1,2,3,
···
,n)，采用恒流恒压方式进行充电，恒流方式进行放电，进行充放电循环直至电池有效容量衰减至标称容量80％；记录第n颗电池每一次充放电循环放电过程中等电压间隔采集的容量q
n,i
和电压v
n,i
的数据集，其中i为充放电循环次数，每个电压值对应一个容量值；q
n,10
为第n颗电池第10次放电过程中的容量数据，q
n,80
为第n颗电池第80次放电过程中的容量数据，以相同电压下容量数据差值q
n,10-80
＝q
n,10-q
n,80
为因变量，放电电压v范围2.2v～3.6v为自变量，绘制函数图像，并通过多项式拟合获得函数公式q
n,10-80
＝fn(v)；
9.(2)计算数值曲线的标准差mse(fn(v))作为输入一维向量x，并将电池单体寿命终止时的充放循环次数作为输出一维向量y；
10.(3)寿命早期预测模型训练：采用基于座头鲸算法优化的支持向量机作为预测算法训练预测模型，训练时采用输入向量x作为模型的输入，采用输出向量y作为模型的输出，对模型进行训练得到预测模型；
11.(4)电池寿命早期预测：在30℃环境温度下，对要预测的电池单体采用拟应用场景下的充电电流值进行恒流恒压充电，充满后，采用拟应用场景下的放电电流值进行放电，提取放电过程中的mse(fn(v))，作为输入值输入至步骤(3)所述的预测模型中，即可预测电池循环寿命。
12.进一步，所述步骤(1)中，n为200～500。优选的，n为300。
13.进一步，所述步骤(1)中，恒流充电阶段电流a设置为a＝4c-n/100c，放电采用3c恒流方式放电。
14.进一步，所述步骤(1)中，i为50～100。
15.进一步，所述步骤(2)中，mse(fn)的计算公式如下：
[0016][0017]
其中，fn(v)为第n颗电池的容量数据差值拟合曲线函数，为fn(v)的平均值，k为电压区间取值个数。
[0018]
更进一步，所述步骤(2)中，k为100～800。
[0019]
进一步，所述步骤(2)中，电池有效容量衰减至标称容量80％时，寿命终止。
[0020]
进一步，所述采用基于座头鲸算法优化的支持向量机作为预测算法训练预测模型的训练方法包括以下步骤：
[0021]
步骤
①
：读取数据：读取向量x，向量y；
[0022]
步骤
②
：初始化支持向量机svm：调用matlab7.10.0函数库中的svmtrain函数进行支持向量机模型训练；
[0023]
步骤
③
：初始化损失函数ε：损失函数ε决定了回归误差期望，其取值大小将会影响相应模型的支持向量数量以及泛化能力，ε取值越大，相应模型支持向量减小，预测精度越低，反之亦然；
[0024]
步骤
④
：初始化惩罚因子参数c、核函数参数σ：
[0025]
核函数参数σ通过座头鲸优化算法在1～300中获取：
[0026]
a＝2a
·
rand1-a
[0027]
s＝2
·
rand2
[0028]
d＝|s
·
σ
q，t-σ
i，t
|
[0029]d′
＝|σ
q，t-σ
i，t
|
[0030][0031]
式中：a和s为计算系数，d为当前鲸鱼个体与当前最优解之间的距离，d
′
为第i只鲸鱼到猎物之间的距离；σ
q，t
为当前最优参数，优选的取150；σ
i，t
为随机参数，取值范围内1～300；σ
i，t 1
为优化后的参数；rand和rand2为[0，1]之间的随机数；a＝2-2t/t为收敛因子，随着迭代的进行从2线性下降到0，t为当前迭代次数，t为最大迭代次数，优选的t为300；p为随机取值的数，其值为0或1；l是一个介于-1和1之间的随机数；b为常数，取值5。
[0032]
惩罚因子参数c主要用于平衡逼近误差以及模型复杂度，c的值越大，相应模型拟合误差则越小，对数据的拟合程度越高，但模型复杂程度亦越大，c一般取为(1～1000)；参
数σ通过座头鲸优化算法在1～1000中获取：
[0033]
a＝2a
·
rand1-a
[0034]
s＝2
·
rand2
[0035]
d＝|s
·cq，t-c
i，t
|
[0036]d′
＝|c
q，t-c
i，t
|
[0037][0038]
式中：c
q，t
为当前最优参数，优选的取500；c
i，t
为随机参数，取值范围内1～1000；c
i，t 1
为优化后的参数；rand和rand2为[0，1]之间的随机数；a＝2-2t/t为收敛因子，随着迭代的进行从2线性下降到0，t为当前迭代次数，t为最大迭代次数，优选的t为300；p为随机取值的数，其值为0或1；l是一个介于-1和1之间的随机数；b为常数，取值5；
[0039]
步骤
⑤
：优化的woa-svm预测模型：将计算的预测寿命与实测的寿命值进行对比，不断修正核函数参数σ、惩罚因子参数c以优化支持向量机，进而得到预测模型。
[0040]
进一步，所述训练方法步骤
③
中，ε为0.0001～0.01。优选的，ε取0.001。
[0041]
进一步，所述训练方法步骤
④
中，核函数包括多项式核函数、径向基函数(rbf)核函数和sigmoid核函数。优选的，核函数为径向基函数(rbf)核函数。
[0042]
本发明的有益效果
[0043]
(1)本发明通过衰减曲线的计算和变换，再结合预测算法进行训练，可实现对锂离子电池寿命进行早期预测，并根据其预测寿命可实现对电池进行有效的分类或者分级；
[0044]
(2)本发明可实现电池寿命早期预测，为电池的研发和制造工艺的改进提供有效依据；
[0045]
(3)本发明可实现对电池组寿命的有效管理，为电池组的维护更换提供依据；
[0046]
(4)本发明可促进电池管理优化，提升电池组的一致性，提高电池组的输出性能和安全性。
附图说明
[0047]
图1为电池寿命早期预测过程流程示意图；
[0048]
图2为放电过程的放电曲线；图2(a)是循环寿命为534周的电池的三条放电曲线，右边的实线为该电池第10次放电的放电曲线，中间实线是第80次放电的放电曲线，左边虚线为第534次放电的放电曲线；图2(b)是循环寿命为788周的电池的三条放电曲线，右边的实线为该电池第10次放电的放电曲线，中间实线是第80次放电的放电曲线，左边虚线为第788次放电的放电曲线；图2(c)是循环寿命为1054周的电池的三条放电曲线，右边的实线为该电池第10次放电的放电曲线，中间实线是第80次放电的放电曲线，左边虚线为第1054次放电的放电曲线；
[0049]
图3为放电电压v—容量差值q
10-80
数值曲线；
[0050]
图4为1号、2号、3号、4号电池循环寿命图。
具体实施方式
[0051]
下面结合具体实施例对本发明进一步说明，本发明的内容完全不限于此。
[0052]
实施例
[0053]
图1示出了寿命预测流程图，步骤如下：
[0054]
(1)模型训练数据获取：对300颗锂离子电池单体进行编号，编号n分别为1至300。在30℃环境温度下，对300颗电池单体采用恒流恒压方式进行充电，恒流充电阶段电流a设置为a＝4c-n/100c；放电采用3c恒流方式放电，直至达到电池放电截止电压停止放电。当电池有效容量衰减至标称容量80％时，代表寿命终止，该电池单体停止充放电试验。
[0055]
记录第n颗电池每一次充放电循环放电过程中等电压间隔采集的容量q
n,i
和电压v
n,i
的数据集，其中i为充放电循环次数，每个电压值对应一个容量值；q
n,10
为第n颗电池第10次放电过程中的容量数据，q
n,80
为第n颗电池第80次放电过程中的容量数据，以相同电压下容量数据差值q
n,10-80
＝q
n,10-q
n,80
为因变量，放电电压v范围2.2v～3.6v为自变量，绘制函数图像，并通过多项式拟合获得函数公式q
n,10-80
＝fn(v)。如图1所示，两条实线左边一条为第10次放电曲线，右边一条为第80次放电曲线，以放电2.5v电压做一条横线，以两条实线的交点的容量值的差值即为2.5v的容量数据差值。在放电电压v范围2.2v～3.6v的区间，每间隔0.01v获取实时的电压值和容量值。图2为放电容量差值-电压的曲线图。
[0056]
以其中随机抽取4颗电池作为测试数据集，4颗电池分别编号为1号、2号、3号、4号电池，其余296颗电池作为预测模型训练数据集为例。
[0057]
(2)计算数值曲线的标准差mse(fn(v))作为输入一维向量x，并将电池单体寿命终止时的充放循环次数作为输出一维向量y；
[0058]
mse(fn)的计算公式如下：
[0059][0060]
其中，fn(v)为第n颗电池的容量数据差值拟合曲线函数，为fn(v)的平均值，为fn(v)的总和除以k。k为电压区间取值个数，为100～800，优选的，k为500。
[0061]
(3)寿命早期预测模型训练：采用基于座头鲸算法优化的支持向量机作为预测算法训练预测模型，训练时采用随机选取的296颗电池，使用上述方法获取输入向量x作为模型的输入，输出向量y作为模型的输出，对模型进行训练得到预测模型。
[0062]
利用训练样本训练支持向量机模型，步骤如下：
[0063]
步骤
①
：读取数据：读取向量x,向量y；
[0064]
步骤
②
：初始化支持向量机svm：调用matlab7.10.0函数库中的svmtrain函数进行支持向量机模型训练；
[0065]
步骤
③
：初始化损失函数ε：损失函数ε决定了回归误差期望，其取值大小将会影响相应模型的支持向量数量以及泛化能力，ε取值越大，相应模型支持向量减小，预测精度越低，反之亦然，ε一般取为(0.0001～0.01)，优选的取0.001；
[0066]
步骤
④
：初始化惩罚因子参数c、核函数参数σ：
[0067]
支持向量机常用核函数包括多项式核函数、径向基函数(rbf)核函数、sigmoid核函数。本实施例选用rbf函数，核函数参数σ通过座头鲸优化算法在1～300中获取：
[0068]
a＝2a
·
rand1-a
[0069]
s＝2
·
rand2
[0070]
d＝|s
·
σ
q，t-σ
i，t
|
[0071]d′
＝|σ
q，t-σ
i，t
|
[0072][0073]
式中：a和s为计算系数，d为当前鲸鱼个体与当前最优解之间的距离，d
′
为第i只鲸鱼到猎物之间的距离；σ
q，t
为当前最优参数，优选的取150；σ
i，t
为随机参数，取值范围内1～300；σ
i，t 1
为优化后的参数；rand和rand2为[0，1]之间的随机数；a＝2-2t/t为收敛因子，随着迭代的进行从2线性下降到0，t为当前迭代次数，t为最大迭代次数，优选的t为300；p为随机取值的数，其值为0或1；l是一个介于-1和1之间的随机数；b为常数，取值5。
[0074]
惩罚因子参数c主要用于平衡逼近误差以及模型复杂度，c的值越大，相应模型拟合误差则越小，对数据的拟合程度越高，但模型复杂程度亦越大，c一般取为(1～1000)；参数σ通过座头鲸优化算法在1～1000中获取：
[0075]
a＝2a
·
rand1-a
[0076]
s＝2
·
rand2
[0077]
d＝|s
·cq，t-c
i，t
|
[0078]d′
＝|c
q，t-c
i，t
|
[0079][0080]
式中：c
q，t
为当前最优参数，优选的取500；c
i，t
为随机参数，取值范围内1～1000；c
i，t 1
为优化后的参数；rand和rand2为[0，1]之间的随机数；a＝2-2t/t为收敛因子，随着迭代的进行从2线性下降到0，t为当前迭代次数，t为最大迭代次数，优选的t为300；p为随机取值的数，其值为0或1；l是一个介于-1和1之间的随机数；b为常数，取值5；
[0081]
步骤
⑤
：优化的woa-svm预测模型：将计算的预测寿命与实测的寿命值进行对比，不断修正核函数参数σ、惩罚因子参数c以优化支持向量机，进而得到预测模型。
[0082]
(4)电池寿命早期预测：测试模型预测效果时采用剩余4颗电池，使用上述方法获取输入向量x作为模型的输入，模型的输出结果即为预测的4颗电池的循环寿命。提取的4颗电池的放电过程中容量变化值q
10-80
如图3所示，循环寿命预测结果如图4所示，本发明的电池寿命早期预测结果与实验结果非常接近，4颗电池预测结果误差均小于9％，表明本发明可以使用电池前80周的早期充放电循环数据有效预测电池循环寿命。
[0083]
对于相同材料配方的电池可实现有效预测：将待预测电池置于30℃环境温度下，对要预测的电池单体采用拟应用场景下的充电电流值进行恒流恒压充电，充满后，采用拟应用场景下的放电电流值进行放电，充放电循环80周后，停止实验。提取放电过程中的q
10
、q
80
，计算mse(fn)作为输入值输入至预测模型中，并可得到由模型预测的电池循环寿命。
[0084]
以上所述，仅为本发明较佳的具体实施方式，但本发明保护的范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内所做的任何修改，等同替换和改进等，均应包含在发明的保护范围之内。