CN113848493A - 基于机器学习的三元锂离子电池早期加速老化诊断方法 - Google Patents
基于机器学习的三元锂离子电池早期加速老化诊断方法 Download PDFInfo
- Publication number
- CN113848493A CN113848493A CN202111042948.XA CN202111042948A CN113848493A CN 113848493 A CN113848493 A CN 113848493A CN 202111042948 A CN202111042948 A CN 202111042948A CN 113848493 A CN113848493 A CN 113848493A
- Authority
- CN
- China
- Prior art keywords
- curve
- discharge
- early
- cycle
- lithium ion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000032683 aging Effects 0.000 title claims abstract description 90
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 title claims abstract description 65
- 229910001416 lithium ion Inorganic materials 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000003745 diagnosis Methods 0.000 title claims abstract description 20
- 238000010801 machine learning Methods 0.000 title claims abstract description 19
- 238000007477 logistic regression Methods 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000007637 random forest analysis Methods 0.000 claims abstract description 9
- 238000010219 correlation analysis Methods 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims description 43
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 7
- 238000005562 fading Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 2
- 238000002405 diagnostic procedure Methods 0.000 claims 3
- 230000036541 health Effects 0.000 abstract description 7
- 238000011156 evaluation Methods 0.000 abstract 1
- 238000007599 discharging Methods 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 5
- 238000013399 early diagnosis Methods 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 3
- 239000010405 anode material Substances 0.000 description 2
- 239000010406 cathode material Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000032677 cell aging Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000007773 negative electrode material Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 239000007774 positive electrode material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/392—Determining battery ageing or deterioration, e.g. state of health
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种基于机器学习的三元锂离子电池早期加速老化诊断方法，从三元锂离子电池放电容量‑电压曲线、放电IC曲线和放电DV曲线上提取表征三元锂离子电池的健康状态的17个老化特征参数，然后利用新的组合算法对三元锂离子电池的加速老化进行早期地准确诊断，首先通过随机森林选择重要特征，然后再通过线性相关分析降低重要特征线性相关性，最后通过逻辑回归模型判断加速老化，实现三元锂离子电池的加速老化早期准确诊断，从而在早期判断三元锂离子电池是否会发生加速老化，为锂离子电池的健康状态管理与健康状态评估提供重要的信息。
Description
技术领域
本发明涉及三元电池健康状态管理技术领域，特别涉及一种机器学习的三元锂离子电池早期加速老化诊断方法。
背景技术
锂离子电池在运行过程中可能会经历加速老化。锂离子电池的加速老化会降低电池性能并导致安全问题。因此锂离子电池早期的加速老化诊断对于锂离子电池健康状态管理十分重要。锂离子电池容量加速衰退会伴随着容量快速衰减，以及内部锂离子和正负极材料的大量损失。锂离子电池目前加速老化诊断相关的方法有基于逻辑回归的方法。各国研究人员现在越来越关注锂离子电池的加速老化现象，并从锂离子电池加速老化机理和拐点识别等方面开展研究，但是加速老化提前诊断方法依然欠缺。传统的机器学习模型可以进行各种故障的诊断，但是三元锂离子电池的衰退机理复杂多变，因此机器学习模型应用到三元锂离子电池的加速老化诊断需要结合电池的相关知识。
发明内容
针对现有技术中存在的缺陷，本发明的目的在于提供一种基于机器学习的三元锂离子电池早期加速老化诊断方法，从放电容量-电压曲线，放电IC曲线和放电DV曲线上提取表征三元锂离子电池的健康状态参数，然后利用机器学习算法对三元锂离子电池进行早期加速老化诊断。
为达到以上目的，本发明采取的技术方案是：
一种基于机器学习的三元锂离子电池早期加速老化诊断方法，包括以下步骤：
步骤S1、选取三元锂离子电池样本，对不同三元锂离子电池样本在不同温度和不同放电倍率条件下进行电池循环衰退试验，得到三元锂离子电池样本在不同循环时期的放电容量-电压曲线、放电IC曲线和放电DV曲线；
步骤S2、利用第n次循环的放电容量-电压曲线、放电IC曲线和放电DV曲线分别与第15次的放电容量-电压曲线、放电IC曲线和放电DV曲线做差，得到放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线；
步骤S3、分别提取放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线在电压区间2.9V-4V上的最大值，最小值，平均值，标准差，偏度以及第n次循环相对于第15次循环的放电容量变化和放电时间变化，一共提取出17个老化特征参数；
步骤S4、对提取的17个老化特征参数进行标准化，得到三元锂离子电池样本的特征矩阵；
步骤S5、通过随机森林算法评估17个老化特征参数的重要性，筛选对加速老化诊断重要的老化特征参数，作为重要特征，然后通过线性相关性分析，降低重要特征之间的相关性，选择出的重要特征的重要性系数总和＞70％，重要特征之间的标准差膨胀因子(VIF)＜10；
步骤S6、利用步骤S5得到的重要特征训练逻辑回归模型，并用训练好的逻辑回归模型进行三元锂离子电池加速老化判断。
在上述方案的基础上，步骤S3中分别提取放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线在电压区间2.9V-4V上的最大值、最小值、平均值、标准差、偏度以及第n次循环相对于第15次循环的放电容量变化和放电时间变化的具体过程为：
ΔQn(Vi)＝Qn(Vi)-Q15(Vi) (1)
ΔdV/dtn(Vi)＝dV/dtn(Vi)-dV/dt15(Vi) (2)
ΔdQ/dVn(Vi)＝dQ/dVn(Vi)-dQ/dV15(Vi) (3)
feature2＝ln(|max(ΔQn(Vi)|) (8)
feature3＝ln(|min(ΔQn(Vi)|) (9)
feature7＝ln(|max(ΔdQ/dVn(Vi)|) (13)
feature8＝ln(|min(ΔdQ/dVn(Vi)|) (14)
feature12＝ln(|max(ΔdV/dtn(Vi)|) (18)
feature13＝ln(|min(ΔdV/dtn(Vi)|) (19)
feature16＝tn-t15 (22)
feature17＝Qdisn-Qdis15 (23)
其中，Qn(Vi)为第n圈放电容量-电压曲线，Q15(Vi)为第15圈放电容量-电压曲线，ΔQn(Vi)为第n圈放电容量-电压曲线与第15圈放电容量-电压曲线之差，作为放电容量-电压曲线的早期变化曲线，dV/dtn(Vi)为第n圈放电DV曲线，dV/dt15(Vi)为第15圈放电DV曲线，ΔdV/dtn(Vi)为第n圈放电DV曲线与第15圈放电DV曲线之差，作为放电DV曲线的早期变化曲线，dQ/dVn(Vi)为第n圈放电IC曲线，dQ/dV15(Vi)为第15圈放电IC曲线，ΔdQ/dVn(Vi)为第n圈放电IC曲线与第15圈放电IC曲线之差，作为放电IC曲线的早期变化曲线；tn为第n次循环的总放电时间，t15为第15次循环的总放电时间，Qdisn为n次循环的总放电容量，Qdis15为第15次循环的总放电容量，
为放电容量-电压曲线的早期变化曲线的平均值，
为放电DV曲线的早期变化曲线的平均值，
为放电IC曲线的早期变化曲线的平均值，k为放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线的上的电压个数，i＝1，2，3，…，k，feature1为放电容量-电压曲线的早期变化曲线的平均值的绝对值的对数，feature2为放电容量-电压曲线的早期变化曲线的最大值的绝对值的对数，feature3为放电容量-电压曲线的早期变化曲线的最小值的绝对值的对数，feature4为放电容量-电压曲线的早期变化曲线的标准差的绝对值的对数，feature5为放电容量-电压曲线的早期变化曲线的偏度，feature6为放电IC曲线的早期变化曲线的平均值的绝对值的对数，feature7为放电IC曲线的早期变化曲线的最大值的绝对值的对数，feature8为放电IC曲线的早期变化曲线的最小值的绝对值的对数，feature9为放电IC曲线的早期变化曲线的标准差的绝对值的对数，feature10为放电IC曲线的早期变化曲线的偏度，feature11为放电DV曲线的早期变化曲线的平均值的绝对值的对数，feature12为放电DV曲线的早期变化曲线的最大值的绝对值的对数，feature13为放电DV曲线的早期变化曲线的最小值的绝对值的对数，feature14为放电DV曲线的早期变化曲线的标准差的绝对值的对数，feature15为放电DV曲线的早期变化曲线的偏度，feature16为第n次循环相对于第15次循环的放电时间变化，feature17为第n次循环相对于第15次循环的放电容量变化，Vi代表第n次循环的放电DV曲线、放电IC曲线和放电容量-电压曲线的电压，Vi∈(2.9V-4V)。
在上述方案的基础上，所述三元锂离子电池样本的特征矩阵具体如下式所示：
其中，
表示第w个电池第n次循环试验提取的第k1个特征，n为循环次数，k1＝1，2，...，17，labelwn表示第n次循环的第w个电池的特征样本所对应的标签，对于加速老化电池的特征样本标签为1，正常老化电池的特征样本标签为0，特征样本为17个老化特征参数组成的样本。
在上述方案的基础上，所述重要特征包括：feature1、feature3、feature6、feature7、feature10、feature12和feature17。
在上述方案的基础上，逻辑回归模型可以输出加速老化的发生概率，能够很好的契合加速老化的不确定性。
三元锂离子电池加速老化的不确定性来源于生产工艺的误差。
在上述方案的基础上，本发明所述方法先通过随机森林选择重要特征，然后再通过线性相关分析降低重要特征之间的线性相关性，最后通过逻辑回归模型判断加速老化。
在上述方案的基础上，采用网格搜索的方法寻找训练好的逻辑回归模型的参数值。
本发明所达到的有益效果为：
本发明从三元锂离子电池放电容量-电压曲线，放电IC曲线，放电DV曲线的早期变化曲线上提取了表征三元锂离子电池的健康状态的17个老化特征。这17个老化特征能反映三元电池在不同温度区间衰退的锂离子损失和内部的正负极材料损失。基于随机森林算法以及线性相关性分析，提取了相关性较低，重要的7个老化特征。本发明提出了一种新的组合方法，组合方法流程是先通过RF(随机森林)选择重要特征，然后再通过线性相关分析降低重要特征线性相关性，最后通过LR(逻辑回归模型)判断加速老化。本发明能够实现三元锂离子电池的加速老化早期准确诊断，从而在早期判断三元锂离子电池是否会发生加速老化。
附图说明
本发明有如下附图：
图1为一种基于机器学习的三元锂离子电池早期加速老化诊断方法流程图；
图2为100次循环的17个老化特征重要性示意图；
图3为200次循环的17个老化特征重要性示意图；
图4为其他老化时期的17个老化特征重要性示意图；
图5为200次循环的17个老化特征的VIF示意图；
图6为去除高VIF特征后，200次循环的重要特征示意图。
具体实施方式
以下结合附图1-6对本发明作进一步详细说明。
机器学习用于各种分类以及模式识别问题，因此机器学习算法适合进行三元电池加速老化早期诊断，在已有大量关于锂离子电池加速老化和正常老化的历史数据的前提下，利用历史数据训练机器学习模型实现三元锂离子电池的早期加速老化诊断。
基于此，本发明提供一种基于机器学习的三元锂离子电池早期加速老化诊断方法，参见图1，包括如下步骤：
步骤1：选取三元锂离子电池样本，对不同三元锂离子电池样本在不同温度和不同放电倍率条件下进行电池循环衰退试验，得到三元锂离子电池样本在不同循环时期的放电容量-电压曲线、放电IC曲线和放电DV曲线；
IC曲线可以表示为dQ/dV关于电压V的函数，IC曲线各个峰以及对应面积的变化对应了三元电池内部不同的衰退模式。
DV曲线可以表示为dV/dQ关于容量Q的函数，DV曲线各个谷以及对应面积的变化对应了三元电池内部不同的衰退模式。
放电曲线反映了电池的电化学过程。
常见的三元锂离子电池内部衰退模式包括锂离子损失，正极材料损失以及负极材料损失。
步骤2、利用第n次循环的放电容量-电压曲线、放电IC曲线和放电DV曲线分别与第15次的放电容量-电压曲线、放电IC曲线和放电DV曲线做差，得到放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线；
步骤3、分别提取三元锂离子电池的放电容量-电压曲线，放电IC曲线和放电DV曲线的早期变化曲线在电压区间2.9V-4V上的老化特征参数，一共提取出17个老化特征参数，老化特征参数能够反映三元电池内部老化机理；
步骤4、对提取完的老化特征参数进行标准化，得到三元锂离子电池样本的特征矩阵。
步骤5：通过随机森林算法评估17个老化特征参数的重要性，筛选对加速老化诊断重要的老化特征参数，作为重要特征，然后通过线性相关性分析，降低重要特征之间的相关性。选择的重要特征的重要性系数总和>70％，重要特征之间的标准差膨胀因子(VIF)<10；
步骤6：利用步骤5得到的重要特征训练逻辑回归模型，并用训练好的逻辑回归模型进行三元锂离子电池加速老化判断。
进一步的，由于放电曲线，DV曲线和IC曲线的主要反应区间可以反映三元电池的内部衰退情况，因此提取放电容量-电压曲线，放电DV曲线和放电IC曲线的早期变化曲线在电压区间2.9V-4V上的老化特征参数，具体公式如下所示：
ΔQn(Vi)＝Qn(Vi)-Q15(Vi) (1)
ΔdV/dtn(Vi)＝dV/dtn(Vi)-dV/dt15(Vi) (2)
ΔdQ/dVn(Vi)＝dQ/dVn(Vi)-dQ/dV15(Vi) (3)
feature2＝ln(|max(ΔQn(Vi)|) (8)
feature3＝ln(|min(ΔQn(Vi)|) (9)
feature7＝ln(|max(ΔdQ/dVn(Vi)|) (13)
feature8＝ln(|min(ΔdQ/dVn(Vi)|) (14)
feature12＝ln(|max(ΔdV/dtn(Vi)|) (18)
feature13＝ln(|min(ΔdV/dtn(Vi)|) (19)
feature16＝tn-t15 (22)
feature17＝Qdisn-Qdis15 (23)
其中，Qn(Vi)为第n圈放电容量-电压曲线，Q15(Vi)为第15圈放电容量-电压曲线，ΔQn(Vi)为第n圈放电容量-电压曲线与第15圈放电容量-电压曲线之差，作为放电容量-电压曲线的早期变化曲线，dV/dtn(Vi)为第n圈放电DV曲线，dV/dt15(Vi)为第15圈放电DV曲线，ΔdV/dtn(Vi)为第n圈放电DV曲线与第15圈放电DV曲线之差，作为放电DV曲线的早期变化曲线，dQ/dVn(Vi)为第n圈放电IC曲线，dQ/dV15(Vi)为第15圈放电IC曲线，ΔdQ/dVn(Vi)为第n圈放电IC曲线与第15圈放电IC曲线之差，作为放电IC曲线的早期变化曲线；tn为第n次循环的总放电时间，t15为第15次循环的总放电时间，Qdisn为n次循环的总放电容量，Qdis15为第15次循环的总放电容量，
为放电容量-电压曲线的早期变化曲线的平均值，
为放电DV曲线的早期变化曲线的平均值，
为放电IC曲线的早期变化曲线的平均值，k为放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线的上的电压个数，i＝1，2，3，…，k，feature1为放电容量-电压曲线的早期变化曲线的平均值的绝对值的对数，feature2为放电容量-电压曲线的早期变化曲线的最大值的绝对值的对数，feature3为放电容量-电压曲线的早期变化曲线的最小值的绝对值的对数，feature4为放电容量-电压曲线的早期变化曲线的标准差的绝对值的对数，feature5为放电容量-电压曲线的早期变化曲线的偏度，feature6为放电IC曲线的早期变化曲线的平均值的绝对值的对数，feature7为放电IC曲线的早期变化曲线的最大值的绝对值的对数，feature8为放电IC曲线的早期变化曲线的最小值的绝对值的对数，feature9为放电IC曲线的早期变化曲线的标准差的绝对值的对数，feature10为放电IC曲线的早期变化曲线的偏度，feature11为放电DV曲线的早期变化曲线的平均值的绝对值的对数，feature12为放电DV曲线的早期变化曲线的最大值的绝对值的对数，feature13为放电DV曲线的早期变化曲线的最小值的绝对值的对数，feature14为放电DV曲线的早期变化曲线的标准差的绝对值的对数，feature15为放电DV曲线的早期变化曲线的偏度，feature16为第n次循环相对于第15次循环的放电时间变化，feature17为第n次循环相对于第15次循环的放电容量变化，Vi代表第n次循环的放电DV曲线、放电IC曲线和放电容量-电压曲线的电压，Vi∈(2.9V-4V)。
进一步的，三元锂离子电池样本的特征矩阵具体如下式所示：
其中，
表示第w个电池第n次循环试验提取的第k1个特征，n为循环次数，k1＝1，2，...，17，labelwn表示第n次循环的第w个电池的特征样本所对应的标签，对于加速老化电池的特征样本标签为1，正常老化电池的特征样本标签为0，特征样本为17个老化特征参数组成的样本。
进一步的，重要特征包括：feature1、feature3、feature6、feature7、feature10、feature12和feature17。
实施例
现对国内某厂家额定116Ah和38Ah的三元锂离子电池进行试验并进行早期加速老化判断，具体实施步骤如下：
步骤1：开展三元锂离子电池在不同温度、不同放电倍率下的电池循环衰退试验，一共36块电池。温度范围是10℃-55℃，放电倍率是0.5C-2C。
步骤2：从三元锂离子电池的放电容量-电压曲线，放电IC曲线和放电DV曲线的早期变化曲线上提取老化特征参数，老化特征参数能够反映三元电池内部老化机理，放电容量-电压曲线，放电IC曲线和放电DV曲线的早期变化曲线具体为：第n次循环的放电容量-电压曲线，放电IC曲线和放电DV曲线与第15次的放电容量-电压曲线，放电IC曲线和放电DV曲线的做差后得到放电容量-电压曲线，放电IC曲线和放电DV曲线的早期变化曲线。
老化特征参数包括放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线的最大值、最小值、平均值、方差、偏度以及第n次循环相对于第15次循环的放电容量变化和放电时间变化。一共有17个老化特征参数反映电池老化机理。
其中，选取放电容量-电压曲线，放电IC曲线和放电DV曲线的早期变化曲线在电压区间2.9V-4V上的老化特征参数。
进一步标准化后得到36块三元电池在不同循环次数n下的特征矩阵：
式中，w为电池编号。
步骤3：利用随机森林算法评估17个老化特征参数的重要性，利用线性相关性分析去除高度共线性的重要特征，feature1(在横坐标轴用1表示)、feature2(在横坐标轴用2表示)、feature3(在横坐标轴用3表示)、feature4(在横坐标轴用4表示)、feature5(在横坐标轴用5表示)、feature6(在横坐标轴用6表示)、feature7(在横坐标轴用7表示)、feature8(在横坐标轴用8表示)、feature9(在横坐标轴用9表示)、feature10(在横坐标轴用10表示)、feature11(在横坐标轴用11表示)、feature12(在横坐标轴用12表示)、feature13(在横坐标轴用13表示)、feature14(在横坐标轴用14表示)、feature15(在横坐标轴用15表示)、feature16(在横坐标轴用16表示)和feature17(在横坐标轴用17表示)，如图2-图4所示，特征1,2,3,4,6,7,8,9,10,12,17具有一定的重要性。如图5-图6所示，特征3，6，10，12和17的VIF都较小，特征2，4，8和9的VIF都大于100。因此综合考虑特征重要性以及特征的多重共线性，我们选择特征1,3,6,7,10,12,17作为最终的重要特征。
步骤4：基于逻辑回归模型，将步骤3中选取的对加速老化诊断重要且线性相关性较低的重要特征作为逻辑回归模型的输入。
将36块电池200次循环的特征矩阵分为两个部分，80％特征矩阵作为训练集，用于训练逻辑回归模型，采用网格搜索的方法寻找合适的参数值。20％特征矩阵作为验证集，用于验证逻辑回归模型，进行加速老化早期诊断的准确度。
本发明所述的老化诊断方法中逻辑回归模型在80％的训练样本上的准确度为100％。在20％的验证样本上的加速老化判断准确度为100％。本发明所述方法加速老化诊断结果展示在表1中。同时作为对比逻辑回归模型结合主成分分析的加速老化诊断结果也展示在表1中作为对比，主成分分析提取解释度为80％的4个主成分。结合主成分分析的逻辑回归模型在训练集上的准确度为100％，但是在验证集上的准确度只有86％。
表1基于组合算法和逻辑回归的加速老化早期诊断结果
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (5)
1.一种基于机器学习的三元锂离子电池早期加速老化诊断方法，其特征在于，包括以下步骤：
步骤S1、选取三元锂离子电池样本，对不同三元锂离子电池样本在不同温度和不同放电倍率条件下进行电池循环衰退试验，得到三元锂离子电池样本在不同循环时期的放电容量-电压曲线、放电IC曲线和放电DV曲线；
步骤S2、利用第n次循环的放电容量-电压曲线、放电IC曲线和放电DV曲线分别与第15次的放电容量-电压曲线、放电IC曲线和放电DV曲线做差，得到放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线；
步骤S3、分别提取放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线在电压区间2.9V-4V上的最大值，最小值，平均值，标准差，偏度以及第n次循环相对于第15次循环的放电容量变化和放电时间变化，一共提取出17个老化特征参数；
步骤S4、对提取的17个老化特征参数进行标准化，得到三元锂离子电池样本的特征矩阵；
步骤S5、通过随机森林算法评估17个老化特征参数的重要性，筛选对加速老化诊断重要的老化特征参数，作为重要特征，然后通过线性相关性分析，降低重要特征之间的相关性，选择出的重要特征的重要性系数总和>70％，重要特征之间的标准差膨胀因子<10；
步骤S6、利用步骤S5得到的重要特征训练逻辑回归模型，并用训练好的逻辑回归模型进行三元锂离子电池加速老化判断。
2.如权利要求1所述的基于机器学习的三元锂离子电池早期加速老化诊断方法，其特征在于，步骤S3中分别提取放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线在电压区间2.9V-4V上的最大值、最小值、平均值、标准差、偏度以及第n次循环相对于第15次循环的放电容量变化和放电时间变化的具体过程为：
ΔQn(Vi)＝Qn(Vi)-Q15(Vi) (1)
ΔdV/dtn(Vi)＝dV/dtn(Vi)-dV/dt15(Vi) (2)
ΔdQ/dVn(Vi)＝dQ/dVn(Vi)-dQ/dV15(Vi) (3)
feature2＝ln(|max(ΔQn(Vi)|) (8)
feature3＝ln(|min(ΔQn(Vi)|) (9)
feature7＝ln(|m2x(ΔdQ/dVn(Vi)|) (13)
feature8＝ln(|min(ΔdQ/dVn(Vi)|) (14)
feature12＝ln(|max(ΔdV/dtn(Vi)|) (18)
feature13＝ln(|min(ΔdV/dtn(Vi)|) (19)
feature16＝tn-t15 (22)
feature17＝Qdisn-Qdis15 (23)
其中，Qn(Vi)为第n圈放电容量-电压曲线，Q15(Vi)为第15圈放电容量-电压曲线，ΔQn(Vi)为第n圈放电容量-电压曲线与第15圈放电容量-电压曲线之差，作为放电容量-电压曲线的早期变化曲线，dV/dtn(Vi)为第n圈放电DV曲线，dV/dt15(Vi)为第15圈放电DV曲线，ΔdV/dtn(Vi)为第n圈放电DV曲线与第15圈放电DV曲线之差，作为放电DV曲线的早期变化曲线，dQ/dVn(Vi)为第n圈放电IC曲线，dQ/dV15(Vi)为第15圈放电IC曲线，ΔdQ/dVn(Vi)为第n圈放电IC曲线与第15圈放电IC曲线之差，作为放电IC曲线的早期变化曲线；tn为第n次循环的总放电时间，t15为第15次循环的总放电时间，Qdisn为n次循环的总放电容量，Qdis15为第15次循环的总放电容量，
为放电容量-电压曲线的早期变化曲线的平均值，
为放电DV曲线的早期变化曲线的平均值，
为放电IC曲线的早期变化曲线的平均值，k为放电容量-电压曲线、放电IC曲线和放电DV曲线的早期变化曲线的上的电压个数，i＝1，2，3，…，k，feature1为放电容量-电压曲线的早期变化曲线的平均值的绝对值的对数，feature2为放电容量-电压曲线的早期变化曲线的最大值的绝对值的对数，feature3为放电容量-电压曲线的早期变化曲线的最小值的绝对值的对数，feature4为放电容量-电压曲线的早期变化曲线的标准差的绝对值的对数，feature5为放电容量-电压曲线的早期变化曲线的偏度，feature6为放电IC曲线的早期变化曲线的平均值的绝对值的对数，feature7为放电IC曲线的早期变化曲线的最大值的绝对值的对数，feature8为放电IC曲线的早期变化曲线的最小值的绝对值的对数，feature9为放电IC曲线的早期变化曲线的标准差的绝对值的对数，feature10为放电IC曲线的早期变化曲线的偏度，feature11为放电DV曲线的早期变化曲线的平均值的绝对值的对数，feature12为放电DV曲线的早期变化曲线的最大值的绝对值的对数，feature13为放电DV曲线的早期变化曲线的最小值的绝对值的对数，feature14为放电DV曲线的早期变化曲线的标准差的绝对值的对数，feature15为放电DV曲线的早期变化曲线的偏度，feature16为第n次循环相对于第15次循环的放电时间变化，feature17为第n次循环相对于第15次循环的放电容量变化，Vi代表第n次循环的放电DV曲线、放电IC曲线和放电容量-电压曲线的电压，Vi∈(2.9V-4V)。
3.如权利要求2所述的基于机器学习的三元锂离子电池早期加速老化诊断方法，其特征在于，所述三元锂离子电池样本的特征矩阵具体如下式所示：
其中，
表示第w个电池第n次循环试验提取的第k1个特征，n为循环次数，k1＝1,2,…,17，labelwn表示第n次循环的第w个电池的特征样本所对应的标签，对于加速老化电池的特征样本标签为1，正常老化电池的特征样本标签为0，特征样本为17个老化特征参数组成的样本。
4.如权利要求3所述的基于机器学习的三元锂离子电池早期加速老化诊断方法，其特征在于，所述重要特征包括：feature1、feature3、feature6、feature7、feature10、feature12和feature17。
5.如权利要求1所述的基于机器学习的三元锂离子电池早期加速老化诊断方法，其特征在于，采用网格搜索的方法寻找训练好的逻辑回归模型的参数值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111042948.XA CN113848493A (zh) | 2021-09-07 | 2021-09-07 | 基于机器学习的三元锂离子电池早期加速老化诊断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111042948.XA CN113848493A (zh) | 2021-09-07 | 2021-09-07 | 基于机器学习的三元锂离子电池早期加速老化诊断方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113848493A true CN113848493A (zh) | 2021-12-28 |
Family
ID=78973127
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111042948.XA Pending CN113848493A (zh) | 2021-09-07 | 2021-09-07 | 基于机器学习的三元锂离子电池早期加速老化诊断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113848493A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118033461A (zh) * | 2024-02-29 | 2024-05-14 | 广东电网有限责任公司 | 一种电池健康状态的评估方法、装置及电子设备 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108051751A (zh) * | 2017-11-25 | 2018-05-18 | 合肥国轩高科动力能源有限公司 | 一种锂离子动力电池分选方法和系统 |
| CN108732508A (zh) * | 2018-05-23 | 2018-11-02 | 北京航空航天大学 | 一种锂离子电池容量的实时估计方法 |
| CN108846227A (zh) * | 2017-12-05 | 2018-11-20 | 北京航空航天大学 | 一种基于随机森林及容量自恢复效应分析的锂离子电池容量退化预测评估方法 |
| CN109946612A (zh) * | 2019-04-01 | 2019-06-28 | 北京交通大学 | 一种三元锂离子电池容量加速衰退转折点识别方法 |
| CN110133527A (zh) * | 2019-05-08 | 2019-08-16 | 深圳市比克动力电池有限公司 | 一种基于三电极锂离子电池分析容量衰减的方法 |
| CN110927607A (zh) * | 2019-11-22 | 2020-03-27 | 武汉理工大学 | 辨识和量化分析锂离子电池衰退机理的方法及系统 |
| CN111077465A (zh) * | 2019-12-25 | 2020-04-28 | 欣旺达电动汽车电池有限公司 | 电池特征参数提取方法、装置、计算机设备及存储介质 |
| CN112382798A (zh) * | 2020-11-12 | 2021-02-19 | 湖南立方新能源科技有限责任公司 | 一种电池循环失效的判断方法及判断系统 |
-
2021
- 2021-09-07 CN CN202111042948.XA patent/CN113848493A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108051751A (zh) * | 2017-11-25 | 2018-05-18 | 合肥国轩高科动力能源有限公司 | 一种锂离子动力电池分选方法和系统 |
| CN108846227A (zh) * | 2017-12-05 | 2018-11-20 | 北京航空航天大学 | 一种基于随机森林及容量自恢复效应分析的锂离子电池容量退化预测评估方法 |
| CN108732508A (zh) * | 2018-05-23 | 2018-11-02 | 北京航空航天大学 | 一种锂离子电池容量的实时估计方法 |
| CN109946612A (zh) * | 2019-04-01 | 2019-06-28 | 北京交通大学 | 一种三元锂离子电池容量加速衰退转折点识别方法 |
| CN110133527A (zh) * | 2019-05-08 | 2019-08-16 | 深圳市比克动力电池有限公司 | 一种基于三电极锂离子电池分析容量衰减的方法 |
| CN110927607A (zh) * | 2019-11-22 | 2020-03-27 | 武汉理工大学 | 辨识和量化分析锂离子电池衰退机理的方法及系统 |
| CN111077465A (zh) * | 2019-12-25 | 2020-04-28 | 欣旺达电动汽车电池有限公司 | 电池特征参数提取方法、装置、计算机设备及存储介质 |
| CN112382798A (zh) * | 2020-11-12 | 2021-02-19 | 湖南立方新能源科技有限责任公司 | 一种电池循环失效的判断方法及判断系统 |
Non-Patent Citations (1)
| Title |
|---|
| 高洋: "三元材料锂离子电池老化诊断评估与建模方法", 《北京交通大学博士论文》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118033461A (zh) * | 2024-02-29 | 2024-05-14 | 广东电网有限责任公司 | 一种电池健康状态的评估方法、装置及电子设备 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114114049B (zh) | 一种基于样本迁移的锂离子电池寿命预测方法 | |
| CN109782190B (zh) | 用于估计单颗电池或单批电池的剩余使用寿命的方法 | |
| CN110058165A (zh) | 锂离子电池容量衰减和寿命预测的数据驱动模型 | |
| CN104569844B (zh) | 阀控密封式铅酸蓄电池健康状态监测方法 | |
| CN101504443B (zh) | 一种锂离子电池放电容量的预测方法 | |
| CN111036575B (zh) | 一种基于温度变化分析的锂离子电池分选方法 | |
| CN110988695A (zh) | 动力电池健康状态评估方法、装置、存储介质及电子设备 | |
| CN113533988B (zh) | 一种锂离子电池长期循环容量衰减分析方法 | |
| CN113075554A (zh) | 一种基于运行数据的锂离子电池组不一致性辨识方法 | |
| CN114280479A (zh) | 一种基于电化学阻抗谱的退役电池快速分选方法 | |
| CN113866655B (zh) | 融合车联网运行数据与测试数据的动力电池容量评估方法 | |
| CN113884929B (zh) | 一种磷酸铁锂电池循环寿命预测方法 | |
| CN112098866B (zh) | 一种判断电池在循环过程中是否发生析锂的无损分析方法 | |
| CN113406502B (zh) | 基于机器学习的缺陷电池筛选方法、设备及介质 | |
| CN111580003A (zh) | 一种基于阻抗谱的二次电池不一致性鉴别方法及装置 | |
| CN111832221A (zh) | 基于特征筛选的锂电池寿命预测方法 | |
| CN110850297A (zh) | 一种通过大数据预测轨交锂电池soh的方法 | |
| CN114720899A (zh) | 退役电池梯次利用分选方法及系统、电子设备、存储介质 | |
| CN111077457A (zh) | 一种梯次利用磷酸铁锂电池加速衰减的评价方法及装置 | |
| CN115494394A (zh) | 鉴定电池质量的电化学方法 | |
| CN115586452A (zh) | 基于新型健康特征的锂离子电池健康状态估计方法 | |
| CN113848493A (zh) | 基于机器学习的三元锂离子电池早期加速老化诊断方法 | |
| CN115236528A (zh) | 锂离子电池循环寿命预测方法 | |
| CN115201692A (zh) | 一种锂电池剩余寿命分析预测方法 | |
| CN114094043A (zh) | 评估锂电池正极材料循环性能的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |