一种动力锂电池寿命数据后台挖掘分析方法与流程

本发明涉及动力锂电池的健康管理技术领域，具体涉及动力锂电池的衰退程度及寿命状态的分析与判断。

背景技术：

电动汽车产业的快速发展，使得锂离子电池得到了更大规模的应用和深入的研究。目前，锂离子电池已从早期能充放500次循环的钻酸锂电池发展到了满区间充放电1500～5000次循环的三元电池。虽然锂离子电池具有重量轻、电压高以及寿命长等性能，但是其性能又会受运行环境、工作状态和老化程度等因素的影响而改变。

因此，亟需寻找到一种方法，能够通过挖掘、分析动力锂电池后台运行数据，掌握电池真实的衰退程度和寿命状态，从而评估电池的充放电性能，优化电池充放电策略，以避免滥用给电池带来永久性的损害。

锂离子电池在充放循环的过程中，存在荷电状态(soc)相同时，开路电压(ocv)却收敛于不同值的现象，研究人员称之为迟滞现象。由于迟滞效应的存在，电池的ocv-soc函数关系并不是一一对应的，根据电流大小、环境温度、充放电历史等因素的不同，同一soc点，可能会得到不同的与之对应的开路电压值。

技术实现要素：

本发明的目的在于针对现有技术的不足或缺陷，提供了一种动力锂电池寿命数据后台挖掘分析方法。该方法通过锂电池在投入使用过程中，对锂电池的后台数据挖掘、分析，得到迟滞电压曲线，并将得到的迟滞电压曲线与初始迟滞电压曲线进行比对分析，从而判断并掌握动力锂电池的真实衰退程度和所处寿命状态。

本发明的目的通过如下技术方案实现。

一种动力锂电池寿命数据后台挖掘分析方法，包括：

s1、在动力锂电池投入使用前，通过间歇性充放电，得到动力锂电池的ocv-soc曲线；

s2、标定初始迟滞电压，得到初始迟滞电压曲线；

s3、动力锂电池投入使用过程中，挖掘、分析后台数据，得到实时迟滞电压曲线，并与所述初始迟滞电压曲线比对分析，得出动力锂电池衰退程度及寿命状态。

作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s1中，所述间歇性充放电，包括如下步骤：

s11、使用cccv法(恒压恒流法)将动力锂电池充满，静置；静置完毕后放电，每次放电2～4％soc后静置，如此循环持续放电，直至电压达到最低截止电压vmin，静置；

s12、静置完毕后，再次开始充电，每次充2～4％soc后静置，如此循环持续充电，直至电压达到最高截止电压vmax；

s13、恒压充电，直至充电电流≤0.05c，静置；静置完毕后，结束充放电。

更优选的，s11、s12、s13中，所述静置均为1～4小时，更进一步优选为4小时。

更优选的，s11、s12、s13中，所述充电和放电均以0.05～1c倍率充电和放电，更进一步优选为0.5～1c；而在更优选的实施例中，所述充电和放电均以1c倍率充电和放电。

在具体的间歇性充放电实施例中，充放电倍率越小，需要静置的时间也就越小。如此，可根据实际实验条件适当减小充放电倍率，而同时，相应的静置时间也可以减小。

作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s2、s3中，采用二阶rc驰豫模型标定迟滞电压，迟滞电压的最大值表示如下：

式中，z表示soc，ocvc(z)表示充电曲线的开路电压测量值，ocvd(z)表示放电曲线的开路电压测量值，g(z)表示由迟滞引起的最大偏差，亦为soc的函数。

作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s3中，所述后台数据包括电池的内阻r值、电池端电压ut、充电电流i1、充电极化电压up1、放电电流i2、放电极化电压up2和实时的soc值。

更优选的，s3中，动力锂电池投入使用过程中的充电公式为：

ut(z)＝ocvc(z)+up1+up2+i1r，(2)

动力锂电池投入使用过程中的放电公式为：

ut(z)＝ocvd(z)-up1-up2-i2r，(3)

式中，z表示soc，ocvc(z)表示充电曲线的开路电压测量值，ocvd(z)表示放电曲线的开路电压测量值。

作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s3中，对所述后台数据挖掘、分析过程中，采用统计学方法识别错误的异常值，并通过包括估算或删除的方法处理错误或重复的数据。

与现有技术相比，本发明具有如下优点和有益效果：

本发明方法通过锂电池在投入使用过程中，对锂电池的后台数据挖掘、分析，得到迟滞电压曲线，并将得到的迟滞电压曲线与初始迟滞电压曲线进行比对分析，从而判断并掌握动力锂电池的真实衰退程度和所处寿命状态，进而优化动力锂电池的充放电策略，以避免滥用给动力锂电池带来永久性的损害。

附图说明

图1为实施例1中动力锂电池进行间歇性充放电的流程图；

图2为实施例1中动力锂电池的充放电ocv-soc曲线图；

图3为实施例1中动力锂电池进行迟滞电压标定采用的二阶rc驰豫模型的示意图；

图4为实施例1中动力锂电池的主回路迟滞曲线的函数拟合图。

具体实施方式

以下结合具体实施例及附图对本发明的技术方案作进一步详细的描述，但本发明的保护范围及实施方式不限于此。

本发明通过对动力锂电池的极化效应原理和迟滞效应原理进行分析，并结合改进的间歇性充放电方法及计算，得到不同衰退程度电池的迟滞电压曲线，从而证明了迟滞电压能够反映动力锂电池的衰退程度。并由此提出本发明的动力锂电池寿命数据后台挖掘分析方法，通过对锂电池的后台数据挖掘、分析，得到迟滞电压曲线，并将得到的迟滞电压曲线与初始迟滞电压曲线进行比对分析，从而判断并掌握动力锂电池的真实衰退程度和所处寿命状态。

本发明的一种动力锂电池寿命数据后台挖掘分析方法，包括：

s1、在动力锂电池投入使用前，通过间歇性充放电，得到动力锂电池的ocv-soc曲线；

s2、标定初始迟滞电压，得到初始迟滞电压曲线；

s3、动力锂电池投入使用过程中，挖掘、分析后台数据，得到实时迟滞电压曲线，并与所述初始迟滞电压曲线比对分析，得出动力锂电池衰退程度及寿命状态。

而作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s1中，所述间歇性充放电为改进的间歇性充电方法，其包括如下步骤：

s11、使用cccv法(恒压恒流法)将动力锂电池充满，静置；静置完毕后放电，每次放电2～4％soc后静置，如此循环持续放电，直至电压达到最低截止电压vmin，静置；

s12、静置完毕后，再次开始充电，每次充2～4％soc后静置，如此循环持续充电，直至电压达到最高截止电压vmax；

s13、恒压充电，直至充电电流≤0.05c，静置；静置完毕后，结束充放电。

采用间歇性充放电的方法可通过充分的静置消除极化，其中，在电流变为零的瞬间(毫秒级或微秒级)，电池的欧姆极化和电化学极化被消除；而在静置的过程中，电池内部的锂离子会继续进行扩散运动，从而消除浓差极化。而且，静置结束后的端电压即可认为是电池在该荷电状态(soc)下的开路电压。如此，通过该改进的间歇性充放电方法进行充放电，不仅可以得到更为准确的开路电压(ocv)，并且可以尽可能的将充放电曲线的soc点对齐，减小了传动测试方法中因对充放电曲线soc＝0％标定方式的不同而产生的误差。

而作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s2、s3中，采用二阶rc驰豫模型标定迟滞电压。二阶rc驰豫模型的示意图参见图3所示，其中，图3的左半部分电容cn代表电池容量，该模型假定当电容满放满充时，其两端电压在0v～1v范围内变化，所以soc与电容两端电压在数值上相等；而图3的右半部分将ocv等效为一个随z(即soc)变化的受控电压源，r为欧姆内阻，vrelaxation代表驰豫效应，且vrelaxation的内部为两个串联的rc并联电路模块，分别代表电化学计划和浓差极化，而vt代表电池端电压。

通过二阶rc驰豫模型标定迟滞电压，得到迟滞电压的最大值表示如下：

式中，z表示soc，ocvc(z)表示充电曲线的开路电压测量值，ocvd(z)表示放电曲线的开路电压测量值；g(z)为soc的函数，是充电所得ocv值与放电ocv值的差值的一半，表示由迟滞引起的最大偏差，亦表示为迟滞电压的最大值。

并且，在具体实施例中，可采用迟滞电压的主回路迟滞曲线的拟合函数表示迟滞电压曲线。

进一步的，作为本发明动力锂电池寿命数据后台挖掘分析方法的一个优选实施例，s3中，所述后台数据包括电池的内阻r值、电池端电压ut、充电电流i1、充电极化电压up1、放电电流i2、放电极化电压up2和实时的soc值。

在动力锂电池投入使用过程中，电池的运行数据会由bms或测试设备收集并上传到监控后台。对于内阻r值及实时的soc值的获取，可通过对电池实时监测得到电池连续的电压电流及运行环境温度等数据，而后，再通过电池管理系统(bms)的内置算法，直接得到电池的内阻r值和实时的soc值；或者，通过电池测试设备，以r＝δu/δi方式计算得到内阻，而实时的soc值由设备的内置算法读取。

更优选的，s3中，动力锂电池投入使用过程中的充电公式为：

ut(z)＝ocvc(z)+up1+up2+i1r，(2)；

动力锂电池投入使用过程中的放电公式为：

ut(z)＝ocvd(z)-up1-up2-i2r，(3)；

式中，z表示soc，ocvc(z)表示充电曲线的开路电压测量值，ocvd(z)表示放电曲线的开路电压测量值。

而在对所述后台数据挖掘、分析过程中，由于采样误差、延时、干扰或设备精度等问题造成上传到后台的数据会存在错误、丢失及重复等问题，分析前首先需要进行数据清洗。可采用包括标准差、四分位距、最大值或最小值的统计学方法识别错误的异常值，并通过包括估算或删除的方法处理错误或重复的数据；如根据电流的方向判断此时电池为连续的充电状态，端电压应随着soc的增大而升高，此时若出现电压降低的数据则证明采集错误，要对该数据进行删除操作，对于重复数据的处理也是如此；若出现部分数据缺失，可以根据其前后数据值，使用平均值、最大值、最小值或更为复杂的概率估计代替缺失的值，从而对数据完成清洗。

对数据完成清洗后，再将数据分为充电方向数据和放电方向数据两组，并分别将两个方向的端电压数据和其对应的soc值对应绘制成ut-soc曲线，即得到式(2)和式(3)所示曲线。在动力锂电池投入使用过程中，内阻r值、电池端电压ut、充电电流i1、充电极化电压up1、放电电流i2和放电极化电压up2均可在后台数据中直接找到，结合式(1)、式(2)和式(3)，即可得到实时的迟滞电压曲线g(z)；再将实时的迟滞电压曲线与标定的初始迟滞电压曲线对比，即可分析得出电池的衰退程度和所处寿命状态。

实施例1

本实施例中采用48ah盟固利三元软包电池作为实验对象。

s1、首先，对新电池进行容量标定。实用额定容量的1c电流按照国标方式对新电池进行充放电循环，相邻循环3次的容量误差小于0.01，取三次平均值视为实际容量。其中，新电池的放电容量略小于充电容量，实际容量标定和soc计算均以放电容量为标准。通过测量得到放电容量为49.58ah。

在动力锂电池投入使用前，通过间歇性充放电，得到动力锂电池的ocv-soc曲线，间歇性充放电流程如图1所示，具体包括如下步骤：

s11、使用cc-cv法以1c倍率(本实施例中为48a电流)将动力锂电池充满，静置4小时；每次以1c放电，每次(本实施例中，每次用时2分24秒)放电4％soc，静置4小时，如此循环持续放电，直至电压达到最低截止电压vmin(本实施例中，最低截止电压3v)；

s12、静置4小时后，开始1c倍率(本实施例中为48a电流)充电，每充4％soc静置4小时，如此循环持续充电，直至电压达到最高截止电压vmax(本实施例中，最高截止电压4.2v)；

s13、转为恒压充电，直至充电电流≤0.05c(本实施例中，充电至电流小于247ma)，静置4小时，结束充电。

其中，静置末端电压视为开路电压，通过上述充放电过程绘制的ocv-soc曲线如图2所示，soc计算所用分母为实际放电容量49.58ah。

s2、对s1过程中的充放电过程得到的ocv-soc曲线采用式(1)的二阶rc驰豫模型(模型图如图3所示)标定初始迟滞电压，得到初始迟滞电压曲线；

具体的，本实施例中，采用初始迟滞电压的主回路迟滞曲线的拟合函数表示初始迟滞电压曲线，如下：

该拟合函数为三段函数，其中前两段为二次函数，第三段为五次函数；式中，a取值0.01653±0.00114，b取值-9.42011e-5±1.07215e-4，c取值-5.69664e-6±2.10883e-6，m取值-0.02579±0.00217，n取值0.00996±4.57471e-4，e取值-4.86034e-4±2.22448e-5，而d5～d0依次分别取值为0.17787±0.02066、-0.0285±0.00356、0.00187±2.37808e-4、-5.92904e-5±7.68054e-6、9.1014e-7±1.20382e-7、-5.44263e-9±7.34321e-10。

s3、动力锂电池投入使用，100个充放循环后进行挖掘分析，本实施例中，电池的运行数据会由bms收集并上传到监控后台。对直接监测的实时的soc值、电池端电压ut、充电电流i1和放电电流i2的数据从后台直接获取；而对于内阻r值，则通过电池管理系统(bms)的内置算法得到。

充电极化电压up1和放电极化电压up2则通过如下公式计算得到：

式(4)为n阶阻容模型的极化电压公式，本实施例中n为2；式中，vp代表极化电压，t为时间常数，rp1cp1、rp2cp2分别为模型中代表电化学极化和浓差极化的电容电阻。v01、v02分别为零输入时的电压。

对后台数据整理如下表1和表2所示。

表1挖掘的后台放电数据

表2挖掘的后台充电数据

将表1、表2根据soc进行数据采用插值法对齐，再结合式(2)以及式(3)，并采用式(1)的二阶rc驰豫模型(模型图如图3所示)标定实时迟滞电压，如表3所示。

表3挖掘的后台迟滞电压数据

得到实时迟滞电压曲线，并对实时迟滞电压曲线进行函数拟合，拟合结果如图4所示。并且，具体的，本实施例中，采用实时迟滞电压的主回路迟滞曲线的拟合函数表示实时迟滞电压曲线，如下：

该拟合函数为三段函数，其中前两段为二次函数(图4中拟合曲线1表示第一段函数，拟合曲线2表示第二段函数)，第三段为五次函数(图4中的拟合曲线3)；式中，a取值-1.11079±0.97716，b取值-0.01488±0.05427，c取值0.00822±6.29907e-4，m取值-0.63823±0.02949，n取值0.21272±0.00988，e取值-0.00485±7.44496e-4，而d5～d0依次分别取值为-0.98801±0.09911、3.10223±0.30234、-3.78007±0.35797、2.21642±0.20501、-0.62671±0.05663、0.07301±0.00602。

将得到的实时迟滞电压曲线与初始迟滞电压曲线对比，发现电池循环次数增多，电池容量开始衰退，其对应的迟滞电压曲线也随之发生变化。衰退后的迟滞曲线的第一个峰值所对应的soc向左偏移，且数值减小，有消失的趋势。由此，通过电池迟滞电压曲线变化趋势可以反映电容衰退程度，从而进行寿命分析和预测，进而优化动力锂电池的充放电策略，以避免滥用给动力锂电池带来永久性的损害。

以上实施例仅为本发明的较优实施例，仅在于对本发明的技术方案作进一步详细的描述，但本发明的保护范围及实施方式不限于此，任何未脱离本发明精神实质及原理下所做的变更、组合、删除、替换或修改等均将包含在本发明的保护范围内。