一种并联型电池系统荷电状态估计方法
【技术领域】
[0001] 本发明属于智能电网中丽级电池储能系统设计与控制技术领域，涉及一种并联 型电池系统荷电状态估计方法。
【背景技术】
[0002] 随着风电、光伏发电等可再生能源及电网智能化的大力发展，电池系统作为电 池储能系统能量存储的主要载体，已越来越多地受到世界各国的关注和应用。同时可 再生能源规模的不断扩大及用电负荷的快速增长，也将促使电池系统向大容量化（MW 级）方向发展，通过多个电池单体的并联可实现电池容量的扩大，即并联型电池系统 (Parallel-connected Battery System, PBS)。然而，由于应用环境的复杂性（如秒级波动 功率平滑、一次高频等高动态场合）及电池电量不能直接测量等因素，准确估计电池系统 荷电状态（State of Charge, S0C)不仅直接决定电池系统能否安全、可靠、高效运行，且对 电池系统优化配置、设计与控制等至关重要。
[0003] 传统的S0C估计算法主要有：安时法、阻抗法、开路电压法等，近年来相继出现了 神经网络、模糊逻辑法、支持向量机及标准卡尔曼滤波法、EKF等高级算法。安时法因其算法 简单、易行等优点，已得到广泛应用，但存在自身开环、误差时间积累等缺点，其精度受限； 开路电压法适合稳态下S0C估计，不宜于在线估计；神经网络、模糊逻辑等高级算法适宜于 恒负载、恒流充放电状态下不同类型电池S0C估计，但存在训练数据量大、训练方法对估计 误差影响大的局限；标准卡尔曼滤波法具有鲁棒性高、抗扰动能力强等优点，适宜于线性系 统的S0C估计，然而电池系统是一种非线性时变系统，其精度仍受限；为此，针对非线性时 变的电池系统，目前常采用EKF，并取得良好的效果，然而由于EKF存在自身计算复杂、忽略 高阶项等问题，必会产生一定误差，使电池的S0C估计精度仍待进一步研究。
【发明内容】
[0004] 本发明解决的问题是在于提供一种基于无迹卡尔曼滤波（Unscented Kalman Filter, UKF)的并联型电池系统荷电状态估计方法，解决并联型电池系统性能参数受SOC 影响、扩展卡尔曼滤波法计算复杂、精度不高而导致电池系统S0C难以被准确测量、估算的 问题，达到准确估计并联型电池系统S0C的目的。
[0005] 本发明目的是通过以下技术方案来实现：
[0006] 本发明提供一种并联型电池系统，该系统由N个电池单体并联而成，其中N均为大 于1的自然数。
[0007] -种并联型电池系统荷电状态估计方法如下：根据已知锂离子电池单体性能参 数，利用并联电路工作特性及其充放电工作特性确定并联型电池系统性能参数与电池单体 性能参数的关系，再结合基尔霍夫定律KVC确定电池系统输出端电压方程，建立并联型电 池系统等效模型（1);将并联型电池系统的荷电状态S0C及等效模型中2个RC并联电路的 端电压作为状态变量，以电池系统的电流及输出电压分别作为系统输入量与输出量，结合 并联型电池系统等效电路模型，得并联型电池系统空间状态方程（2);将并联型电池系统 空间状态方程（2)中的电池系统S0C、2个RC并联电路的端电压作为无迹卡尔曼滤波算法 UKF的状态变量；以并联型电池系统空间状态方程（2)的输入状态空间方程、输出电压状态 空间方程分别作为UKF算法的非线性状态方程及测量方程；通过电压传感器测量并联型电 池系统端电压（4)的实际值与UKF算法获得的电池端电压估计值来更新增益矩阵（5)，最后 由UKF算法经循环迭代，从而实时得到电池系统SOC的估计值。
[0008] 所述并联型电池系统等效电路模型（1)为二阶等效电路模型，模型主电路由2个 RC并联电路、受控电压源UbQ(S0C)及电池内阻Rb等组成。建立准确的电池系统等效电路模 型关键在于如何根据电池工作特性来确定电池系统性能参数与电池单体性能参数的关系。 本发明中并联型电池系统性能参数与电池单体性能参数关系式为：
[0010] 式中，1^、1^1、(；；3、(；1分别表示电池系统模型中2个此并联电路的电阻和电容；下 标i表示第i个电池单体；分别表示电池单体的开路电压、内阻；^^^、(^、(^分别 表示电池单体模型中2个此并联电路的电阻和电容；1]1。、1?1、1^、1?11、(： 1；3、(：11均与30(：有关， S0C的定义为：
，其中，S0C。为电池单体S0C初始值，一般 为0~1的常数；Qu(t)为电池单体不可用容量，Q。为电池单体额定容量。U 1(](S0C)、Rls、Ru 和(：1；3、(：11、1?1的计算分别如下：
d。~d 2、e。~e 2、f。~f 2、bQ~b 5均为模型系数，可由电池测量数据经拟合而得。
[0011] 所述并联型电池系统空间状态方程（2)的建立如下：a、以电池系统的荷电状态 S0Cp及等效模型中2个RC并联电路的端电压作为状态变量，以电池系统的电流I p为系统 输入量，根据等效电路模型建立电池系统空间状态方程为
[0013] 式中，Ups、Upl为2个RC并联电路端电压，R ps、Rpl为2个RC并联电路的电阻，Q 电池系统额定电量，τρ τ2为时间常数，Wk为系统观过程噪声，At为采样周期，k为大于 1的自然数；b、根据基尔霍夫电压定律，结合电池系统等效电路模型，可得电池系统输出电 压方程空间状态方程为：
，式中，Up为电 池系统端电压，Rp为电池系统内阻。
[0015] 所述无迹卡尔曼滤波算法UKF的主要步骤为：1)初始化状态变量X均值 E0和均方误差P。：
2)获取采样 点Xi及对应权重ω :
式中，λ = a2(n+h)-n，〇Γ、〇^分别表示方差及均值 的权重，a、β分别表示采样点中粒子分布距离及高阶项误差大小；3)状态估计及均方 误差的时间更新：状态估计时间更新为
，均 方误差时间更新为
，系统输出时间 更新为 式中，gkl(·)为测量方程；4)计算增益矩阵（5):
;5)状态估计 及均方误差的测量更新：状态估计测量更新为
，均方误差测 量更新为，
[0016] 与采用扩展卡尔曼滤波算法EKF进行并联型电池系统S0C估计相比，本发明具有 以下有益的技术效果：一是整个放电过程，本发明所采用的UKF算法比EKF算法进行并联型 电池系统S0C估计时UKF估计精度更高，尤其是放电初期和末期效果更明显；二是所采用的 UKF算法比EKF算法能更快收敛于实验数据，鲁棒性更好。
【附图说明】
[0017] 图1为一种并联型电池系统荷电状态估计方法流程图；
[0018] 图2为并联型电池系统结构示意图；
[0019] 图3为含2个电池单体的并联型电池系统结构示意图；
[0020] 图4为含2个RC并联电路的并联型电池系统等效电路模型图；
[0021] 图5为无迹卡尔曼滤波算法流程图；
[0022] 图6-1~图6-4为S0C。不同时电池恒流放电特性，其中图6-1为S0CQ= 1时S0C 变化情况，图6-2为S0CQ= 1时电池系统端电压变化情况，图6-3为S0C。= 0. 8时S0C变 化情况，图6-4为SOQi 0. 8时电池系统端电压变化情况；
[0023] 图7-1~图7-4为S0C。不同时电池脉冲放电特性，其中图7-1为S0CQ= 1时S0C 变化情况，图7-2为S0CQ= 1时电池系统端电压变化情况，图7-3为S0C。= 0. 8时S0C变 化情况，图7-4为SOQi 0. 8时电池系统端电压变化情况。
【具体实施方式】
[0024] 下面结合具体的实例对本发明作进一步的详细说明，所述为对本发明的解释而不 是限定。
[0025] 根据本发明实施例，如图1、图2、图3、图4和图5所示，提供了一种并联型电池系 统荷电状态估计方法，实施例的流程图如图1所示，主要包括以下几个步骤：
[0026] 1、建立并联型电池系统等效电路模型
[0027] 1)并联型电池系统
[0028] 并联型电池系统是由由N个电池单体并联而成，其结构图如图2所示。为便于分 析，本实例中假设并联型电池系统由2个电池单体经并联而成，即1 X 2并联型电池系统，如 图3所示。并联型电池系统中每个电池单体的额定电压为3. 2V，额定容量为25Ah，放电截 止电压为2. 5V。
[0029] 2)建立1