基于TWP SVR的锂离子电池健康状态估计.doc

《基于TWP SVR的锂离子电池健康状态估计.doc》由会员分享，可在线阅读，更多相关《基于TWP SVR的锂离子电池健康状态估计.doc（22页珍藏版）》请在第一文库网上搜索。

1、基于TWP SVR的锂离子电池健康状态估计目前我国真正达到低碳节能标准意义的建筑不到10%，其他的90%都是耗能建筑，而且新建房屋建筑中95%以上仍为高能耗建筑，但我国资源占有量不到世界平均水平的20%，单位建筑面积能耗却是气候相近发达国家的35倍。因此针对超低能耗、近零能耗建筑，开发安全、高效、绿色、经济、耐久的保温装饰一体化板保温系统具有极其重要的价值，而储能电池是智能传感应用的关键。现今，锂离子电池(lithium-ion batteries，LIBs)已成为储能系统电源的主要选择。为提高电池系统实际运行的安全性、稳定性和可靠性，电池SOH估计成为至关重要的研究热点。SOH通常根据容量衰

2、减和内阻变化来对电池老化程度进行定量评估。然而，目前无论是容量还是内阻都不能用市售的传感器直接测量。因此，间接分析方法是开发SOH估计方法的一个重要方向。间接分析的关键在于如何能从低成本传感器采集的数据中提取出间接健康特征参数来表征SOH。从电压、温度曲线中提取间接健康特征的分析方法已经获得了广泛的研究。ICA通过测量电池在不同循环周期中电压与电池增量容量之间关系的变化趋势，提取出电池的退化特征参数。差分电压分析(differential voltage analysis，DVA)通过分析电压对容量的导数与电池循环退化过程中的峰谷变化，能够检测和量化退化机制。考虑到电池使用中电池温度的变化，差

3、热伏安法(differential thermal voltammetry，DTV)通过电压和温度的差分比率计算参数，参数曲线的峰值变化反映了电池容量衰减信息。曲线相似性分析是一种常用的聚类分析方法，在大规模电池组故障分析研究中，通过对不同电池单体的电压曲线聚类，能够分析出电池组的健康状况。对同一电池单体不同充放循环的电压曲线进行相似性分析可以估计出电池单体SOH。参考文献11引入平均弗雷歇距离(average Frchet distance，AFD)方法来确定充放电曲线的相似性，实验表明AFD是良好的SOH估计参数。参考文献12为解决不同长度的循环曲线序列，使用动态时间规整(dynamic

4、time warping，DTW)算法和长短期记忆网络(long short term memory networks，LSTM)的组合模型估计电池的SOH。基于曲线相似性分析思想，本研究提出了一种使用时间规整图TWP提取间接健康特征参数，使用支持SVR模型估计SOH的方法。流程如图1所示。图 1TWP-SVR方法流程图首先，通过等间隔采样获取LIBs初始循环的放电电压曲线作为参考序列，依次获取不同循环的放电电压曲线作为采样序列。然后，使用DTW算法找出最优规划路径，通过坐标转换获得TWP并提取出4个间接健康特征。其次，对健康特征数据和容量数据进行标准化处理，获得标准化的训练数据集。接着，用处

5、理后的数据集训练基于线性核的SVR模型。最后，通过对新循环的电压曲线进行采样获得间接健康特征，使用训练好的SVR模型估计SOH。1 方法1.1 特征提取算法1.1.1 动态时间规整动态时间规整是一种准确率高、鲁棒性强的时间序列相似性度量方法，能够通过弯曲时间序列的时域对时间序列的数据点进行匹配，不仅能够得到更好的形态度量效果，而且能够度量两条非等长的时间序列相似性。考虑两条不同长度的时间序列和；其中S1是采样序列，S2是参考序列。根据两条序列任意两点的欧氏距离构造距离矩阵A为(1)其中。如图2所示，动态时间规整距离D(S1,S2)可以定义为：在距离矩阵A中找到一条最优的规整路径W使得S1和S2

6、的累积距离值达到最小。即：(2)其中，表示规整路径W在距离矩阵A中的位置。为了保证公式(2)的解唯一性，给出3条约束性要求：边界线：单调性：对于给定的和，必有连续性：对于给定的和，必有可以使用动态规划方法求解公式(2)，即：(3)其中初始条件D(1,1)=d11。图 2DTW算法示意图1.1.2 时间规整图相对于欧氏距离，动态时间规整距离通过时间规整，能够较好地描述时间序列的相似性。但这主要是通过时间序列幅值的差异性表征的，忽视了时间序列中时域上的相位差异。为解决这一问题，参考文献15基于DTW算法获得的最优规整路径，提出了一种基于坐标变换的时间规整图(time warp profile)方法

7、估计不同时间序列之间相位变化的差异性。首先，通过DTW算法获得两条时间序列的距离矩阵，将距离矩阵的对角线和反对角线作为相位参考系(t,)，并将原点固定在(1，1)处，如图3(a)中红线坐标所示。如图3(a)所示将规整路径上的每一点Wk都进行坐标变换。具体方法如以下公式(4)(6)：(4)(5)(6)其中。利用公式(4)(6)将规整路径W转换成时间规整图TWP中的相位差异序列(t)，如图3(b)所示。在TWP中(t)序列值表示时间序列之间的局部滞后，也即是相位差异性，索引指数t对应于两个序列同采样时刻的时间。图 3TWP坐标变换示意图1.1.3 间接健康特征TWP将DTW中的二维的最佳规整路径W

8、压缩成一维的相位差异序列，序列的统计性质能反映出参考序列S0和采样序列S1的相似性，为了更方便地比较S0和S1的相似性，选取4个序列的统计指标。(7)其中，反映了序列与0值的偏差，越大序列越远离0值，表明S0和S1幅值相似性越低。反映了序列变化率的大小，越大随着时间变化越陡峭，表明S0和S1相位差异性越大。是序列平均值，越大序列离0值越大，表明S0和S1幅值相似性越低。反映了序列波动性，越大幅值变化越大，说明S0和S1相位差异性越大。因此，利用以上4个统计特征，能够反映出S0和S1的相似性。在LIBs循环老化实验中，考虑电池充放循环中放电时间较长，且曲线长度不一致。为了缩小模型规模和统一参数数

9、目。选择初始循环的放电电压曲线作为参考序列S0，选择其余循环的放电电压曲线作为采样序列Si，依次使用上述算法获得采样序列与参考序列的相位差异序列。利用以上4个特征参数计算初始循环充放电压曲线和其余循环充放电压曲线的相似性，利用不同循环电压曲线的相似性间接表征LIBs的老化特征。1.2 SOH拟合模型目前，有多种使用间接健康特征估计SOH的模型。神经网络模型有：LSTM，门控循环单元(gated recurrent unit，GRU)，注意力机制模型和混合模型。统计回归模型有GPR模型。机器学习模型有支持向量机(support vector machine，SVM)模型。本研究选择SVR模型作为

10、SOH拟合模型。SVR原理如下：对于数据集(X1,Y1),(X2,Y2)(XN,YN),考虑线性关系的时候，有回归函数：(8)其中W和b是模型的权重参数和偏置参数。支持向量机回归SVR将拟合公式(8)转换为求解凸二次优化问题，SVR引入不敏感损失参数和松弛变量，求解问题转变为：(9)其中，C为惩罚参数，用于控制对位于边界之外的观测值施加的惩罚，这有助于防止过度拟合。松弛变量和用于描述样本不满足约束的程度。引入拉格朗日乘子算子求解公式(9)的对偶问题：(10)其中和是拉格朗日乘子。根据KKT条件，得到回归函数：(11)引入核函数之后可得回归函数：(12)核函数有多种选择，包括线性核、多项式核、高

11、斯核以及混合核。1.3 评价指标为了评估选取的4个间接健康指标能否反映电池容量衰减特性。选择三个相关性指标进行评价：皮尔森(Pearson)系数、斯皮尔曼(Spearman)系数、肯德尔(Kendall)系数。(13)(14)(15)其中，表示第i次循环某一个间接健康特征；表示多次循环的某一个间接健康特征值的平均值；表示第i次循环的电池容量；表示多次循环的电池容量平均值；表示排列在第i位的F和Q的秩差；sgn是符号函数。三个相关系数可以说明电池容量衰减曲线分别和4个间接健康特征的相关性，相关性系数绝对值越高，说明两个参数的相关性越大，可以评估4个间接健康特征是否反映出电池容量的衰减特性。为了更

12、全面地评估所提方法的准确性，给出三个常用的评价指标：RMSE、MAE、MAPE。(16)(17)(18)其中，表示真实值；表示模型的预测值。RMSE、MAE和MAPE值越小说明模型准确性越高。为了评估所提方法应用于不同电池的稳定性，使用不同电池实验的同类型准确性指标的样本标准差和四分位差两个指标说明方法的稳定性。(19)(20)其中，表示不同LIBs实验获得的同类型准确性指标，表示所有电池实验的准确性指标平均值。SSD表示样本标准差(sample standard deviation)，反映了对于不同的电池进行实验，所提方法准确性指标的波动性，SSD越小说明所提方法越稳定。表示准确性指标序列的

13、第75个百位数，也叫第3个四分位数。表示准确性指标序列的第25个百位数，也叫第1个四分位数。IQR是准确指标序列的四分位距(inter-quartile range)，IQR越小说明准确性指标分布越集中，波动性越小，可以表明所提方法越稳定。2 实验验证与分析2.1 数据描述与实验过程为了验证所提出方法估计SOH的准确性及其适用性，本工作使用了2个开源数据和实际储能电站运行数据进行实验，数据集的相关信息见表1。本工作所有的实验在CPU型号为i7-6500U2.50 GHz，RAM内存为12 G的设备上使用Matlab R2021a软件进行分析。所提算法由特征提取和拟合模型两部分组成。实验使用NA

14、SA数据集，首先，根据相关系数验证TWP算法提取出的4个间接健康特征能否反映电池容量衰减特性。然后，评估TWP类型的特征提取算法分别接入SVR模型和GPR模型后的准确性。使用UL-PUR数据集验证TWP-SVR和TWP-GPR两种方法的准确性和稳定性，最后使用储能电站的数据验证所提TWP-SVR方法的准确性和稳定性。表1实验数据集信息2.2 特征提取算法的验证为了验证TWP提取4个间接健康特征与电池容量衰减的相关性，本工作使用NASA数据集中5、6、7、18号电池进行分析。在每一次充放循环实验中，选取2 A(1 C)恒定电流(CC)放电阶段的电压曲线作为实验数据。注意到电压数据是非等间隔采样获

15、得，所以在提取特征之前先进行等间隔插值处理。将第1次循环的放电电压曲线作为参考序列，不同循环的放电电压曲线作为采样序列，利用公式(7)从TWP中提取4个间接健康特征。使用公式(3)、(4)、(5)的相关系数评估本研究所提出的4个间接健康特征和电池衰减容量的相关性。从图4可以看出，间接健康特征和衰减容量的相关系数均大于0.9，具有强相关性，因此所提的间接健康特征可以作为表征电池容量衰减和估计SOH的参数。图4中相关性系数值是4个电池的平均值，详细的实验数据见表2。图 4间接健康特征和衰减容量的相关性表 2间接健康特征和衰减容量的相关系数为了验证本工作所提特征提取算法的准确性，本工作使用上文提到的

16、NASA开源的4个2000 mAh的电池数据与2种算法进行比较分析。这些算法是由特征提取和拟合模型组合的方法，算法详细说明见表3。表 3实验所用算法信息使用表3中的算法对4个电池使用前15%的数据作为训练数据进行SOH估计，实验结果表明TWP作为特征提取的算法具有较高的准确性。图5是使用两种算法对4个电池数据进行SOH估计的结果，图6是准确性指标平均值统计图。由图6可以看出，TWP类型算法在小规模的训练数据的条件下具有较高的准确性，在取前15%的数据作为训练数据的条件下，TWP-SVR算法和TWP-GPR算法的平均RMSE分别为0.031和0.018，平均MAE为0.027和0.015，平均MAPE为3.905和