一种考虑单体不一致的锂离子动力电池组在线参数估计方法
【专利摘要】本发明涉及一种考虑单体不一致的锂离子动力电池组在线参数估计方法，用以获取锂离子动力电池组中各个单体的参数，包括欧姆内阻、极化电阻和时间常数，包括以下步骤：1)获取锂离子动力电池组的总体平均参数；2)分别构造欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C并根据总体平均参数分别获取锂离子动力电池组中各个单体的参数。与现有技术相比，本发明具有考虑单体不一致性、计算简单等优点。
【专利说明】
一种考虑单体不一致的锂离子动力电池组在线参数估计方法
技术领域
[0001]本发明涉及动力电池组领域，尤其是涉及一种考虑单体不一致的锂离子动力电池 组在线参数估计方法。
【背景技术】
[0002] 目前的电池管理系统(Battery Management System，BMS)对电池内部参数进行估 计时，往往将整个电池组看成均匀的。实际上，电池在制作和装配过程中，由于工艺和材质 的不均匀等，使得同一类型、规格、型号的电池在出厂时的参数不可能完全一致。而在电池 的使用过程中，每个电池单体的充放电过程、自放电程度、通风条件和温度环境等的差异在 一定程度上也会增加电池参数的不一致性电池组内部各单体之间存在不一致。对于成组的 电池来说，由于其内阻不一致性，在充电过程中，相同的充电电流下，内阻大的电池升压高， 容易导致过充;而在放电过程中，相同的放电电流下，内阻大的电池压降高，容易导致过放。 因此，以平均参数代替各个单体参数进行管理是不合理的。在电池的使用过程中，随着循环 次数的增加，电池的不一致性会加剧，电池组很有可能会因为某一节单体的损坏而提前报 废，甚至引起安全隐患，导致系统使用和维护成本增加。这一问题已经成为制约电动汽车产 业化发展和应用的瓶颈问题之一。
[0003] 由于单体不一致的问题，我们必须要对每一节单体电池进行参数估计，对于BMS来 说，电池信息显示必须实时。因此，如果我们直接将传统参数辨识方法在每节单体电池上实 施，将导致计算量过大，从而就对处理器提出较高的要求，势必会增加 BMS的成本。因此，本 发明所提的方法是对现有估计算法进行改进和优化，使之能在不增加计算量或者只增加少 量计算量的前提下估计出考虑一致性的参数值。
【发明内容】
[0004] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑单体不一 致性、计算简单的考虑单体不一致的锂离子动力电池组在线参数估计方法。
[0005] 本发明的目的可以通过以下技术方案来实现：
[0006] -种考虑单体不一致的锂离子动力电池组在线参数估计方法，用以获取锂离子动 力电池组中各个单体的参数，包括欧姆内阻、极化电阻和时间常数，包括以下步骤：
[0007] 1)获取锂离子动力电池组的总体平均参数；
[0008] 2)分别构造欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C并根据总体 平均参数分别获取锂离子动力电池组中各个单体的参数。
[0009] 所述的步骤1)中通过最小二乘法或自适应滤波法获取锂离子动力电池组的总体 平均参数。
[0010]当通过最小二乘法获取锂离子动力电池组的总体平均参数时包括以下步骤：
[0011] 11)根据多个电池单体串联的电池组结构建立电池平均模型，该模型为：
[0013] 其中，UL为端电压，U。。为开路电压，UThm为平均极化电压，II为电流，R〇 m为平均欧姆 内阻，RThm为平均极化电阻，CThm为平均极化电容，为，
[0014] 12)构建电池平均模型的差分方程，并采用递推公式，获取动力电池组的总体平均 参数，包括平均欧姆内阻、平均极化电阻和平均时间常数，
[0015] 所述的差分方程为：
[0023] 动力电池组的总体平均参数为：
[0025] 其中，Ud,k、Ud,k+i分别为电池组k时刻和k+Ι时刻的平均动态电压向量，IL,k、lL,k+i分 别为电池组k时刻和k+Ι时刻的电流，T m为平均时间常数且Tm=RThmCTt?，Δ t为采样周期，Pk、 Pk+i分别为k和k+1时刻的协方差矩阵，Lk+i为k+1时刻的增益矩阵。
[0026] 所述的步骤2)中，欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C的表达 式为：
[0029] 其中，Udn为电池组中第η个单体的平均动态电压，υΒη为电池组中第η个单体基于权 向量Β的动态电压，U Cn为电池组中第η个单体基于权向量C的动态电压。
[0030] 所述的步骤2)中，
[0031]欧姆内阻权向量Α的递推公式为：
[0032] UP,k = Ud mean, kAk-i
[0033 ] Ak = Ak-ι+gAUd mean, k(Ud,k-UP,k)
[0034] 极化电阻权向量B的递推公式为：
[0035] UBp,k = UBmean,kBk-l
[0036] Bk = Bk-l+gBUBmean, k ( Ub, k~UBp, k )
[0037] 时间常数权向量C的递推公式为：
[0038] UcP>k = Ucmean,kCk-l
[0039] Ck = Ck-l+gcUcmean, k ( UC> k'Ucp, k )
[0040] 其中，UP,k为电池组k时刻的预测动态电压，Udm^k为k时刻所有电池组单体动态电 压的算数平均值，Ud, k为电池组k时刻的平均动态电压向量，gA为欧姆增益，UBp, k为电池组k时 亥IJ的UB的预测值，UBmean, k为k时刻电池组各单体对应UB的算数平均值，gB为极化增益，UCp, k为 电池组k时刻的UC的预测值，UCmean,Ak时刻电池组各单体对应U B的算数平均值，gc为时间增 益。
[0041 ]与现有技术相比，本发明具有以下优点：
[0042] -、考虑单体不一致性:本发明分别根据每节单体的端电压，荷电状态等参数估计 出每节电池的参数，而不是把整包电池看成均匀状态来求出平均参数。
[0043]二、计算简单:本发明从原理出发，推导出考虑单体不一致性的权向量，计算时只 需将平均参数与权向量相乘即可得到每节单体的参数，与目前的方法相比，大大减小了计 舁里。
【附图说明】
[0044] 图1为本发明算法流程图。
[0045] 图2为本发明实施例涉及的电池模型图。
[0046] 图3为本发明实施例的两种电池模型示意图，其中，图（3a)为本发明实施例的电池 串联模型示意图，图（3b)为本发明实施例的电池平均模型示意图。
[0047] 图4为本发明实施例的电池参数曲线，其中，图（4a)为本发明实施例中两节电池动 态电压曲线，图（4b)为本发明实施例中两节电池欧姆内阻曲线。
【具体实施方式】
[0048] 下面结合附图和具体实施例对本发明进行详细说明。
[0049] 实施例：
[0050] 为了使本发明的目的、技术方案及创新点更加清晰，以下结合附图及实施例对本 发明作进一步的阐释。应当理解，此处所描述的具体实施例仅用于解释本发明，并不用于限 定本发明。
[0051 ]如图1所示，图为本发明的算法流程图。
[0052]图2所示为本发明实施例涉及的电池模型图，其中U。。为开路电压，仇为端电压，II 为电流Rq为欧姆内阻，RTh和Orh分别为极化电阻和极化电容，I Th为极化电容的电流，UTh为极 化电压，模型的数学关系为：
[0054]为简化起见，假设一个简单的电池模组，其中有两节电池串联，其阻抗模型如图3a 所示。首先，可以构建一个上述两节电池的平均阻抗模型，如图3b所示，RQm为平均欧姆内阻， RThm为平均极化电阻，CThm平均极化电容。通过一般的参数估计方法，我们可以求出电池组的 平均参数，此处以最小二乘法为例。电池组的差分方程如下式：
[0055] Ud,k+i = aUd,k+blL,k+i+clL,k
[0056] 上式中，Ud,k为k时刻平均端电压和平均开路电压之差，即电池组平均动态电压， 1^为袖寸刻的电流， &、13、(3为：
[0058] 上式中，RQm为平均欧姆内阻，RThm为平均极化电阻，τ η为平均时间常数，e为自然对 数，At为采样周期。最小二乘法的递推公式如下：
[0060]上式中，堯为k时刻的参数估计值，zk+i = Ud,k+i，hk+i= [Ud,k+i，lL,k+i， IL,k]。设定好初值之后可以递推出a、b、c的值，然后根据下式可以求出平均参数：
[0062] 求出平均参数之后，下面求单体参数，如图3a所示，有：
[0063] Roi: R〇2 = Uoi : U〇2
[0064] 上式中Uoi和U〇2分别为两节电池欧姆内阻Roi和R〇2上的电压。但是在实际计算中Uoi 和UQ2的求解比较困难，而实际情况中发现欧姆内阻之比与动态电压之比比较接近。本发明 的实施例，如图4a和4b所示，在动态电压曲线上随机选取了两个时刻，两节电池动态电压的 比值为1.06和1.07,两节电池欧姆内阻的比值为1.05。
[0065] 因此，我们可以推断：
[0066] R〇i：R〇2^Udi：Ud2
[0067] 根据上述比例关系，我们可以构造一个权向量A:
[0071] 根据平均欧姆内阻和向量A，我们可以求出两节电池的欧姆内阻：
[0072] [ Roi, R02 ] = R〇mA
[0073] 这样就可以求出两节电池的欧姆内阻，然后根据动态电压Ud和电流II可以求出极 化电压UTh:
[0075]根据模型的数学关系可以推出：
[0077]由于采样周期At-般远远小于时间常数τ，因此有：
[0081 ] 上式整理得：
[0083] 令：
[0084] UB = UTh,k+l/UTh,k
[0085] 这样我们就可以构造一个权向量B:
[0089]根据平均时间常数~和权向量B，我们可以求出：
[0091]这样就求出了两节电池的时间常数τ#Ρτ2，由电池模型的数学关系得：
[0093]因此有两节电池的极化电阻之比：
[0097]这样我们就可以构造一个权向量C:
[0099] 根据RTh4P权向量C，我们可以求出：
[0100] [RThl，RTh2]=RThmC
[0101] 然后，根据两节电池的时间常数和极化电阻的值可以求出两节电池的极化电容： [0103]以上是两节电池的的推导，同理可以推出多节电池的权向量：
[0106] 如果直接利用上述权向量进行求解，由于电池电压的波动比较大，在参数估计过 程中会造成参数比较大的波动。首先，我们必须对权向量进行归一化处理，权向量实质上反 应了所有电池参数的比率，当参数全部相同时，A = B = C=[1，1，1，......，1]，权向量中的元 素大于1说明对应的参数大于平均参数，小于1说明对应的参数小于平均参数。
[0107] 其次，为了减缓波动，我们可以采用递推的方式得到权向量。下面以基于动态电压 Ud求权向量Ak为例，说明权向量的求解方法。k为当前时刻，首先利用上一时刻的权向量预测 当前的时刻的动态电压值：
[0108] UP,k = Ud mean, kAk-i
[0109] 上式中udmean，ks当前时刻的平均动态电压值，然后利用预测的动态电压值和真实 动态电压值之差来修正当前时刻的权向量：
[0110] Ak = Ak-l+gAUd mean, k(Ud,k-UP,k)
[0111] 上式中，gA为增益，增益越大，修正量越大，算法越快，也可能出现过修调的情况， 因此精度越低。
[0112] 同理可以求出权向量Bk和Ck的递推公式：
[0115] 上两式中，UBmean,k和UCmean,k分别为k时刻U B和Uc的平均值，gB和gc分别为权向量B和C 的增益。
【主权项】
1. 一种考虑单体不一致的裡离子动力电池组在线参数估计方法，用W获取裡离子动力 电池组中各个单体的参数，包括欧姆内阻、极化电阻和时间常数，其特征在于，包括W下步 骤： 1) 获取裡离子动力电池组的总体平均参数； 2) 分别构造欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C并根据总体平均 参数分别获取裡离子动力电池组中各个单体的参数。2. 根据权利要求1所述的一种考虑单体不一致的裡离子动力电池组在线参数估计方 法，其特征在于，所述的步骤1)中通过最小二乘法或自适应滤波法获取裡离子动力电池组 的总体平均参数。3. 根据权利要求2所述的一种考虑单体不一致的裡离子动力电池组在线参数估计方 法，其特征在于，当通过最小二乘法获取裡离子动力电池组的总体平均参数时包括W下步 骤： 11) 根据多个电池单体串联的电池组结构建立电池平均模型，该模型为：其中，化为端电压，Uoc为开路电压，UThm为平均极化电压，II为电流，Rom为平均欧姆内阻， RThm为平均极化电阻，CThm为平均极化电容，馬/,,,，为， 12) 构建电池平均模型的差分方程，并采用递推公式，获取动力电池组的总体平均参 数，包括平均欧姆内阻、平均极化电阻和平均时间常数， 所述的差分方程为： Ud,k+1 -过 Ud, k+b II, k+1+C II, k动力电池组的总体平均参数为：其中，Ud,k、Ud,k+i分别为电池组k时刻和k+1时刻的平均动态电压向量，lL,k、lL,k+i分别为 电池组k时刻和k+1时刻的电流，Tm为平均时间常数且Tm = RT〇：Tto，At为采样周期，Pk、Pk+l分 别为k和k+1时刻的协方差矩阵，Lk+1为k+1时刻的增益矩阵。4. 根据权利要求1所述的一种考虑单体不一致的裡离子动力电池组在线参数估计方 法，其特征在于，所述的步骤2)中，欧姆内阻权向量A、极化电阻权向量B和时间常数权向量C 的表达式为：其中，Udn为电池组中第η个单体的平均动态电压，化η为电池组中第η个单体基于权向量B 的动态电压，化η为电池组中第η个单体基于权向量C的动态电压。5. 根据权利要求4所述的一种考虑单体不一致的裡离子动力电池组在线参数估计方 法，其特征在于，所述的步骤2)中， 欧姆内阻权向量A的递推公式为： Up,k = Udmean,kAk-l Ak 二 Ak-l+gAUdmean，k ( Ud，k_Up，k ) 极化电阻权向量B的递推公式为： U化，k = UBmean，kBk-l Bk 二 Bk-1+卽 UBmean，k ( Ub，k_U化，k ) 时间常数权向量c的递推公式为： Ucp, k - Ucmean, kCk-1 Ck 二 Ck-1+容cUcmean，k ( Uc，k-Ucp，k ) 其中，Up, k为电池组k时刻的预测动态电压，Udmean，k为k时刻所有电池组单体动态电压的 算数平均值，Ud，k为电池组k时刻的平均动态电压向量，gA为欧姆增益，UBp，k为电池组k时刻的 化的预测值，UBmean,k为k时刻电池组各单体对应Ub的算数平均值，即为极化增益，U也k为电池 组k时刻的化的预测值，Ucmean,k为k时刻电池组各单体对应化的算数平均值，gC为时间增益。
【文档编号】G06F19/00GK106096281SQ201610417242
【公开日】2016年11月9日
【申请日】2016年6月15日 公开号201610417242.X, CN 106096281 A, CN 106096281A, CN 201610417242, CN-A-106096281, CN106096281 A, CN106096281A, CN201610417242, CN201610417242.X
【发明人】戴海峰, 孙泽昌, 魏学哲, 曾雷
【申请人】同济大学