一种串联型电池系统的荷电状态估计方法

本发明属于智能电网中大容量电池储能系统设计与控制，涉及一种串联型电池系统荷电状态估计方法。

背景技术：

1、随着人们对电动汽车的动力学和续航里程要求的提高，单个电池的电极电势和材料技术水平的限制已经不能满足需求，因此需要将多个电池单体串联组成大容量电池系统，即串联电池系统。但随着电池数量的增多，电池不一致性的问题是不可避免的，将影响着串联型电池系统荷电状态(state of charge,soc)的准确估计，整个串联电池系统的荷电状态估计影响着电池管理系统的稳定运行，关系着电动汽车的行驶安全和使用寿命。因此准确的电池系统荷电状态估计有着重要的意义。

2、目前，国内外关于电池荷电状态估计的方法主要分为安时积分法、开路电压法、基于模型的滤波法和数据驱动法，但这些方法都有不足之处。其中安时积分法难以克服累计误差和测量误差造成的影响；开路电压法因需要长时间静置不能应用于实时估计；基于模型的滤波法需要建立精确的模型；数据驱动法需要大量数据进行训练。目前主要的估计方法为基于滤波法中的卡尔曼滤波族法，但是由于系统噪声未知，会导致滤波算法的估计精度不高、稳定性差的缺点，所以常用噪声估计器获取系统噪声信息去代替滤波中的噪声统计信息。如文献(cn107831448a)公开的一种并联型电池系统的荷电状态估计方法，该方法基于极大似然准则和最大期望算法的思想，构建出带次优递推的噪声估计器，提高了估计精度、减小了计算量。但是其中自适应系数的取值决定着噪声估计器计算的权重值，本发明通过使用滑动窗口的方法对权重取值不确定造成误差的问题提出了解决方法，从而提高估算精度；此外，目前关于串联型电池系统方面的文献不是很多，文献(cn105183934b)公开的一种基于参数校正器串联型电池参数建模方法，该方法通过结合soc估计模块与参数校正器模块对串联型电池系统的荷电状态进行修正，然后通过修正值去对串联型电池系统模型进行更新，建立准确的串联型电池系统模型，但是该方法没有考虑到随着时间推移串联型电池系统内部电池不一致性是无法避免的，本发明针对此问题提出了对通过不一致性检测器来对电池不一致性进行检测，从而获得精确的电池系统内部参数。

技术实现思路

1、本发明的目的在于针对上述问题，提出一种串联型电池系统荷电状态估计方法，一是通过带滑动窗口的噪声估计器对噪声统计信息进行平稳化处理，解决噪声统计信息估计中权重值对模型精度影响较大的问题；二是通过计算适应度因子来控制无迹卡尔曼滤波法ukf迭代次数，解决无迹卡尔曼滤波法ukf迭代速度不可控的问题；三是通过对欧姆内阻导致压降幅度不一致性进行检测，从而对串联型电池系统中不一致性情况进行检测，获取精确的电池系统参数，解决串联型电池系统中不一致性导致电池系统参数难以被准确测量的问题，最终得到精确的串联型电池系统荷电状态估计值。

2、本发明目的是通过以下技术方案来实现：

3、步骤一：结合电池单体模型参数与串联电路工作特性确定串联型电池系统模型(1)并建立一个二阶rc等效电路模型；步骤二：根据等效电路模型获取空进状态方程结合在线检测电压值与预测电压值作为输入量，经带滑动窗口的自适应无迹卡尔曼滤波法aukf模块(2)，即首先通过带滑动窗口的噪声估计器(3)获取系统噪声信息(4)去代替无迹卡尔曼滤波法ukf中噪声系统信息值，然后经无迹卡尔曼滤波法ukf(5)计算出采样点的测量值与预测值并结合实际观测值构建适应度函数计算出适应度因子(6)，根据比值大小判断是否需要重新迭代，然后经状态估计更新，最终得到状态估计值socb；步骤三：将电池单体离线数据经不一致性检测器模块(7)获取修正电压值，即首先通过数据均值筛选模块(8)进行均值计算初步对电池参数不一致性进行筛选，然后经欧姆内阻不一致性判别模块(9)根据欧姆内阻不一致会影响放电压降幅度不一致性进行判断是否出现不一致性现象，剔除出现不一致性的电池数据，最后得到筛选后的电压u；步骤四：结合修正后电压值与在线检测电压值经参数校准器(10)得到荷电状态补偿值δsocb；步骤五：将估计值socb与补偿值δsocb迭加得到校准值socr，再进一步更新电池系统模型，并得到下一时刻的电池系统socb，如此循环迭代。

4、所述串联型电池系统是由n个电池单体串联而成，n为大于的自然数，如图2所示。所述串联型电池系统等效电路模型由两个并联rc电路、受控电压源u0(soc)及电池内阻r组成，电路图如图3所示。根据基尔霍夫定律kvc得电池模型表达式为：ub(t)＝ub0[soc(t)]-ib(t)zb(t),式中ub0[soc(t)]为电池系统的开路电压，zb(t)为电池系统的阻抗，其计算如下：

5、

6、式中rb(t)为电池系统内阻，rbs(t)、rbl(t)和cbs(t)、cbl(t)分别描述电池系统暂态响应特性的电阻、电容，其相关计算如下：ub0[soc(t)]＝nu0(soc)、rb(t)＝nr(t)、rbs(t)＝nrs(t)、rbl(t)＝nrl(t)、其中r(t)为电池单体内阻，rs(t)、rl(t)和cs(t)、cl(t)为电池单体暂态响应特性。以上参数均与荷电状态soc有关，soc的定义为：

7、

8、式中soc0为电池单体soc初始值，一般为0～1的常数，c0为电池单体额定容量，所述电池单体参数与soc的计算关系为：

9、

10、

11、

12、

13、

14、

15、其中a0～a5、b0～b5、c0～c2、d0～d2、e0～e2、f0～f2均为模型系数，可由电池测量数据经拟合而得。

16、所述空间状态方程以电池socb及两个rc端电压作为系统状态变量xk,以ub,ib分别作为系统量测变量yk以及系统输入变量，设计如下：

17、

18、式中ubs和ubl为两个并联rc电路的端电压，τ1和τ2为时间常数，ωk为系统噪声，δt为采样周期，根据霍尔夫电压定律，结合电池等效电路模型yk＝[ub,k]＝ub0,k-rb,kib,k-ubs,k-ubl,k+vk,vk为系统观测噪声，k为大于1的自然数。

19、所述带滑动窗口的自适应无迹卡尔曼滤波法aukf的步骤如下：

20、(1)系统初始化状态估计:

21、(2)计算采样点：

22、式中λ＝α2(n+h)-n,n为状态变量的维度；ωm、ωc分别表示方差及均值的权重，算子为对称阵的cholesky分解，α、β、h均为常数；

23、(3)时间更新：通过公式将所有sigma点向后更新，可以得到系统状态的先验估计：均方误差时间更新为：为状态方程噪声协方差。系统输出时间更新：采样点对应的测量值预测测量值式中g(·)为测量方程。

24、(4)计算适应度因子：根据预测测量值与实际观测的适应度函数以及采样点对应的测量值与实际观测值的适应度函数比值求得适应度因子ρ：式中f1和f2分别为预测测量值与实际观测的适应度函数以及采样点对应的测量值与实际观测值的适应度函数，ωi为采样点对应权值，r为测量噪声协方差。若ρ<1，表示采样点有效逼近真实估计值，则不需要迭代计算通过无迹卡尔曼滤波法ukf进行状态估计更新得到电池荷电状态估计值socb；若ρ≥1，表示采样点与真实估计偏差较大，则需要重新计算采样点进行迭代更新。

25、(5)状态估计更新：式中lk为卡尔曼滤波增益，py,k为自协方差，pxy,k为自互协方差，rk为状态方程噪声方差。

26、状态估计与状态协方差修正：

27、所述带滑动窗口的噪声估计器为式中为状态噪声协方差qk的观测值；δxk是量测更新前后的状态变化；pk为状态估计协方差；为k时刻的状态噪声估计值,n是滑动窗口宽度。

28、所述数据均值筛选模块设计如下：以电池单体离线数据中电压、电容、欧姆内阻为输入，对其相关数据均值计算，相关数据均值包括:uc、qc、rc分别为单体电池的电压、可用容量和欧姆内阻的平均值，ui,qi,ri分别为第i个电池的电压、可用容量和欧姆内阻，从而对池组不一致性检测进行初步筛选。

29、将筛选后的数据进入欧姆内阻不一致性判别模块，所述模块是根据放电起始阶段，电池电压存在一个瞬间的下降，该电压下降值δu是因为欧姆内阻分压导致的，其数学关系为：δu＝i*r0,式中i为放电电流，r0为欧姆内阻，根据压降幅值的大小可间接反应欧姆内阻大小情况；放电阶段由电池单体电压计算得到电池组电压其公式为：式中δur为电池组欧姆内阻压降幅值；δur,i为第i个电池欧姆内阻压降幅值；δur为单体电池的压降幅度；n为电池组中单体电池数量；将电池组与单体电池压降幅度关系进行线性拟合,并对其求导得到线性关系变化速率k′，其表达式如下：根据开始电池不一致性影响小的情况下所得的变化速率k′和随时间变化所得到电池数据得到的变化速率k′进行比较，若出现变化速率幅度差距较大的情况，则说明电池组之间出现不一致性现象；

30、所述参数校准器设计如下：计算得到补偿值δsocb，ksoc为校正系数。

31、最后将状态更新值socb与补偿值δsocb迭加计算得校准值socr，进而更新电池系统等效模型，如此循环。

32、与公开文献相比，本发明具有以下有益的技术效果：第一，考虑到在噪声统计信息估计中权重值对模型精度影响很大，通过带滑动窗口的噪声估计器对噪声统计信息进行平稳化处理，提高了估计精度；第二，针对常规无迹卡尔曼滤波法ukf迭代速度不可控的问题，提出通过计算自适应因子实现迭代次数的自适应控制，从而提高了ukf迭代速度；第三，考虑到电池系统中不一致性问题无法避免的问题，首先通过对电容、电压、欧姆内阻数据进行均值处理对检测不一致性先进行初步筛选，然后针对欧姆内阻导致压降幅度变化的问题对其变化幅度进行检测，获取精度的电池参数，最终提高荷电状态估计精度。