一种基于粒子滤波算法的锂电池SOC估算方法与系统与流程

一种基于粒子滤波算法的锂电池soc估算方法与系统
技术领域
1.本发明属于电池soc估算技术领域，更具体地说，是涉及一种基于粒子滤波算法的锂电池soc估算方法与系统。


背景技术：

2.随着能源紧缺和全球环境恶化的问题日益严峻，开发新能源等节能减排技术已经引起了全世界各国的关注。目前，具有自放电小、能量密度高和循环寿命长等优异性能的锂电池现已在光伏储能电站、电动汽车、通讯等领域得到广泛的应用。而为了满足应用对象的性能需求，通常将多块锂电池进行串并联组合，从而增大输出电压、输出电流及容量，但是这也给电池增加了安全隐患，在电池实际运行过程中，常常伴随着过充、过放、过热以及电池不一致性等问题的产生，导致电池的使用效率降低和电池的使用寿命缩短，且随着电池循环次数的增加，电池不一致性的问题会越来越明显，对于上述状况需要采取相应的措施进行应对。通常，借助电池管理系统实现对电池的均衡管理和荷电状态(state of charge，soc)、健康状态(state of health，soh)的估算。其中，锂电池soc是评估电池剩余可用能量的重要指标之一，快速准确的评估锂电池soc可以防止由于电池过充过放等导致的永久性内部恶化现象的发生，提升锂电池的安全性。
3.同时，由于锂电池内部存在复杂的物理化学反应，且在soc的估算过程中受到电流、电压、温度等多个因素的影响，很难通过直接检测得到电池soc值，只能通过这些影响因素进行推算间接得到，而电流、电压等参数在锂电池的实际运行过程中会不断发生变化，因此实时估算锂电池的soc值非常困难，如何实现锂电池soc值在线精确的估算成为锂电池技术的一大难点。
4.目前，国内外研究电池soc值的获取方法主要是从电池内部和电池外部两个方面入手，电池内部入手主要是通过对电池内部机理的研究，基于物质能量守恒定律结合电池的物理性质计算电池的剩余电量。从电池外部入手则是需要构建数学模型，在获取电池运行状态下外部输入参数的基础上，借助相应的算法对电池的soc值进行估算。有学者利用开路电压法准确估算出电池的soc，操作方法简单，但是该方法在估算过程中需要将电池进行长时间的静置，不能对电池soc进行在线估算。而安时积分法也可以在较短的时间内获取准确的电池soc值，但是由于难以获取准确的soc初始值，该方法在估算过程中会出现误差累计。另有学者使用改进的粒子滤波算法估算电池的soc值，但该方法在粒子重采样阶段，仅仅根据权值和阈值大小的对比对粒子进行重排，忽略了粒子集枯竭现象，从而对锂电池soc的估算精度造成影响。


技术实现要素：

5.本发明的目的在于提供一种基于粒子滤波算法的锂电池soc估算方法与系统，旨在解决现有对锂电池soc的估算精度低的问题。
6.为实现上述目的，本发明采用的技术方案是：
7.一种基于粒子滤波算法的锂电池soc估算方法，包括以下步骤：
8.步骤1：根据锂电池内部电路构建等效电路模型；
9.步骤2：对所述锂电池的ocv曲线进行拟合得到等效电路模型中锂电池的开路电压；
10.步骤3：使用人工蜂群算法对等效电路模型中作为输出参数的锂电池的端电压进行辨识，得到最优阻抗参数；
11.步骤4：根据所述最优阻抗参数，使用粒子滤波算法对锂电池的soc值进行估算。
12.优选的，所述步骤1中等效电路模型的微分方程为：
[0013][0014]
式中，r1表示电化学极化的电阻，r2表示浓差极化的电阻，c1表示电化学极化的电容，c2表示浓差极化的电容，r0表示欧姆内阻，u
t
表示端电压，u1表示电化学极化电压，u2表示浓差极化电压，u
oc
表示开路电压，i表示锂电池工作电流。
[0015]
优选的，所述步骤3：使用人工蜂群算法对等效电路模型的参数进行辨识得到最优阻抗参数，包括：
[0016]
步骤3.1：确定人工蜂群算法的输入/输出参数；其中，锂电池的工作电流作为输入参数并视为侦察蜂，锂电池的端电压作为输出参数并视为最优蜜源；
[0017]
步骤3.2：使用侦察蜂在近邻搜索蜜源；
[0018]
步骤3.3：根据蜜源花蜜量所对应的适应度函数判断蜜源的好坏，并确定出最优蜜源位置；
[0019]
步骤3.4：根据最优蜜源位置得到最优阻抗参数。
[0020]
优选的，所述步骤3.2：使用侦察蜂在近邻搜索蜜源，包括：
[0021]
采用公式：
[0022][0023]
在近邻搜索蜜源；其中，x
id
＝[r
0 r
1 r
2 c
1 c2]
t
，为[-1，1]的随机数，i＝k+1(k＝1，2，
…
，99)，k表示一个不等于i的蜜源，d＝1，2，
…
，5，id表示当前蜜源的的位置，x
kd
表示已选择的蜜源附近的蜜源，x
′
id
表示选择的蜜源。
[0024]
优选的，所述适应度函数为：
[0025][0026]
其中，ek表示适应度函数，表示模型端电压，u
t
表示实测的端电压。
[0027]
优选的，所述步骤4：根据所述最优阻抗参数使用粒子滤波算法对锂电池的soc值进行估算，包括：
[0028]
步骤4.1：根据锂电池的soc值，初始化粒子数和粒子权重；
[0029]
步骤4.2：对每一个粒子的权重进行更新并归一化，得到归一化后的粒子权重；
[0030]
步骤4.3：对每一个粒子进行重采样，得到新的粒子集；
[0031]
步骤4.4：使用状态方程对锂电池的soc进行估算，得到估算完成的soc值。
[0032]
优选的，所述状态方程表达式为：
[0033][0034]
其中，vk表示系统的观测噪声，wk表示系统的过程噪声，δt表示系统的采样周期，x
k+1
表示粒子的下一个状态，ik表示放电电流，ηi表示放电效率，η
t
表示充电效率，ηn表示转换效率，qn表示电池的额定容量，y
k+1
表示状态观测值，k0、k1、k3、k4表示常数项系数，r表示电池电阻值，n表示电池循环次数。
[0035]
本发明还提供了一种基于粒子滤波算法的锂电池soc估算系统，包括：
[0036]
等效电路模型构建模块，用于根据锂电池内部电路构建等效电路模型；
[0037]
曲线拟合模块，用于对所述锂电池的ocv曲线进行拟合得到等效电路模型中锂电池的开路电压；
[0038]
参数辨识模块，用于使用人工蜂群算法对等效电路模型的参数进行辨识得到最优阻抗参数；
[0039]
soc值估算模块，用于根据所述最优阻抗参数使用粒子滤波算法对锂电池的soc值进行估算。
[0040]
本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现上述的一种基于粒子滤波算法的锂电池soc估算中的步骤。
[0041]
本发明提供的一种基于粒子滤波算法的锂电池soc估算方法与系统的有益效果在于：与现有技术相比，本发明通过构建等效电路模型，并使用人工蜂群算法实现等效电路模型参数的辨识，可以在获得精准的模型参数基础上，采用粒子滤波算法对锂电池的soc值进行估算，不仅简化了计算过程，还可以得到准确的估算结果。
附图说明
[0042]
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0043]
图1为本发明实施例提供的一种基于粒子滤波算法的锂电池soc估算方法流程图；
[0044]
图2为本发明实施例提供的锂电池等效电路模型示意图。
具体实施方式
[0045]
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
[0046]
本发明的目的在于提供一种基于粒子滤波算法的锂电池soc估算方法与系统，旨在解决现有对锂电池soc的估算精度低的问题。
[0047]
请参阅图1，一种基于粒子滤波算法的锂电池soc估算方法，包括以下步骤：
[0048]
步骤1：根据锂电池内部电路构建等效电路模型；
[0049]
图2为本发明根据锂电池内部电路构建的等效电路模型(即二阶等效电路模型)。二阶等效电路模型，即双极化(dual polarization，dp)模型在计算过程较为便利的同时也保障了估算结果的准确性。因此本发明选用dp模型，来描述锂电池的动态特性。其dp模型的微分方程为：
[0050][0051]
式中，r1表示电化学极化的电阻，r2表示浓差极化的电阻，c1表示电化学极化的电容，c2表示浓差极化的电容，r0表示欧姆内阻，u
t
表示端电压，u1表示电化学极化电压，u2表示浓差极化电压，u
oc
表示开路电压，i表示锂电池工作电流。
[0052]
步骤2：对所述锂电池的ocv曲线进行拟合得到等效电路模型中锂电池的开路电压；
[0053]
在实际应用中，本发明可根据供应商给的电池运行数据得到ocv曲线，并采用lsqcurvefit函数对该曲线进行拟合，其拟合公式如下：
[0054]uoc
＝a1soc6+a2soc5+a3soc4+a4soc3+a5soc2+a6soc+b
[0055]
其中，a1、a2、a3、a4、a5、a6、b表示常数项系数，u
oc
表示开路电压。
[0056]
步骤3：使用人工蜂群算法对等效电路模型的参数进行辨识得到最优阻抗参数；
[0057]
进一步的，步骤3包括：
[0058]
步骤3.1：确定人工蜂群算法的输入/输出参数；其中，锂电池的工作电流作为输入参数并视为侦察蜂，锂电池的端电压作为输出参数并视为最优蜜源。
[0059]
进一步的，在本发明实施例中，本发明还要设定种群的参数初始值，由于人工蜂群算法(abc算法)是在5维空间进行求解进而对问题进行优化的方法，因此本发明假设引领蜂、观察蜂和蜜源的数量相同，其中，且以模型端电压与实测电压u
t
的差值的绝对值作为适应度函数进而对蜜源的优劣进行评判。
[0060]
步骤3.2：使用侦察蜂在近邻搜索蜜源；
[0061]
需要说明的是，侦察蜂在近邻搜索蜜源，可以找到蜜源的称为引领蜂，否则称为观察蜂，引领蜂在近邻搜索蜜源，其表达式为：
[0062][0063]
其中，x
id
＝[r
0 r
1 r
2 c
1 c2]t，为[-1，1]的随机数，i＝k+1(k＝1，2，
…
，99)，d＝1，2，
…
，5，x
kd
表示已选择的蜜源附近的蜜源，x
′
id
表示选择的蜜源。
[0064]
步骤3.3：根据蜜源花蜜量所对应的适应度函数判断蜜源的好坏，并确定出最优蜜源位置；其中，适应度函数为：
[0065][0066]
其中，ek表示适应度函数，表示模型端电压，u
t
表示实测的端电压。
[0067]
具体的，当新蜜源优于原蜜源时，更新蜜源的位置，并计算出观察蜂的跟随概率(观察蜂通过跟随概率选择蜜源)，计算公式如下：(观察蜂通过跟随概率选择蜜源)，计算公式如下：其中，pi表示跟随概率，fi表示第i个蜜源的适应度值，fitj表示第j个种群的适应度值，abs(fi)表示第i个蜜源的适应值的绝对值。
[0068]
步骤3.4：根据最优蜜源位置得到最优阻抗参数。
[0069]
经过步骤3.1-3.3的迭代过程，即可确定出最优蜜源位置，即可精确辨识到dp模型在最小电压误差下的最优阻抗参数值。
[0070]
步骤4：根据所述最优阻抗参数使用粒子滤波算法对锂电池的soc值进行估算。
[0071]
进一步的，所述步骤4包括：
[0072]
步骤4.1：根据锂电池的soc值初始化粒子数和粒子权重；
[0073]
在本发明实施例中，可采用先验概率p(x0)产生n个soc初始粒子其相应的粒子权重{q
0i
}
ni＝1
＝1/n。
[0074]
步骤4.2：对每一个粒子的权重进行更新并归一化得到归一化后的粒子权重；
[0075]
需要说明的是，在步骤4.2中，本发明更新粒子权重采用公式为：需要说明的是，在步骤4.2中，本发明更新粒子权重采用公式为：归一化权值的公式为：
[0076]
表示先验条件概率，表示粒子权重值，表示粒子权重值，表示目标状态的后验概率分布。
[0077]
步骤4.3：对每一个粒子进行重采样得到新的粒子集；
[0078]
进一步的，在重新采样之前，本发明需要计算最小均方估计，计算公式为：然后计算有效粒子数n
eff
，计算公式为：n
eff
＝1/当n
eff
≤ns时，ns表示设定的阈值，可以得到若干粒子集
[0079]
步骤4.4：使用状态方程对锂电池的soc进行估算得到估算完成的soc值。
[0080]
进一步的，本发明的步骤4.4包括：
[0081]
(1)采用状态方程对未知参数进行预测，表示电池soc值。
[0082]
(2)根据程序结束条件判断程序是否结束，若未结束，则回到更新粒子权重步骤继续执行。
[0083]
需要说明的是，在粒子滤波算法对电池soc进行估算过程中，其状态方程表达式如下：
[0084][0085]
其中，vk表示系统的观测噪声，wk表示系统的过程噪声，δt表示系统的采样周期。x
k+1
表示粒子的下一个状态，ik表示放电电流，ηi表示放电效率，η
t
表示充电效率，ηn表示转换效率，qn表示电池的额定容量，y
k+1
表示状态观测值，k0、k1、k3、k4表示常数项系数，r表示电池电阻值，n表示电池循环次数。
[0086]
本发明通过构建等效电路模型，并使用人工蜂群算法实现等效电路模型参数的辨识，可以在获得精准的模型参数基础上，采用粒子滤波算法对锂电池的soc值进行估算，不仅简化了计算过程，还可以得到准确的估算结果，有利于锂电池得到更广泛的发展与应用。
[0087]
本发明还提供了一种基于粒子滤波算法的锂电池soc估算系统，包括：
[0088]
等效电路模型构建模块，用于根据锂电池内部电路构建等效电路模型；
[0089]
曲线拟合模块，用于对所述锂电池的ocv曲线进行拟合得到等效电路模型中锂电池的开路电压；
[0090]
参数辨识模块，用于使用人工蜂群算法对等效电路模型的参数进行辨识得到最优阻抗参数；
[0091]
soc值估算模块，用于根据所述最优阻抗参数使用粒子滤波算法对锂电池的soc值进行估算。
[0092]
本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现上述的一种基于粒子滤波算法的锂电池soc估算中的步骤。
[0093]
本发明提供的一种基于粒子滤波算法的锂电池soc估算方法与系统的有益效果在于：与现有技术相比，本发明通过构建等效电路模型，并使用人工蜂群算法实现等效电路模型参数的辨识，可以在获得精准的模型参数基础上，采用粒子滤波算法对锂电池的soc值进行估算，不仅简化了计算过程，还可以得到准确的估算结果。
[0094]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。