基于遗忘因子自适应反馈修正的电池能量状态评估方法与流程

1.本发明涉及新能源电池测控领域，尤其涉及基于遗忘因子自适应反馈修正的电池能量状态评估方法。


背景技术：

2.在锂离子电池组的整个生命周期中，电池管理系统（battery management system，bms）对核心参数soe（state of energy，能量状态）的监控和调节将影响应急动力输出的效果和安全性；因此，实时监测该参数的变化，并基于此保障锂离子电池组的工作性能是非常必要的。由于bms中的成组soe估算技术尚不成熟，使用过程中存在的安全隐患严重制约了锂离子电池组的发展。
3.对于锂离子电池组而言，可靠的bms管理依靠准确的soe值；在该值已知的情况下，不仅能够对其进行可靠能量管理和安全控制，而且还能避免锂离子电池组的提前损坏，延长其使用寿命；因此，精确估算soe值，对保障锂离子电池组的工作性能及其能量和安全管理至关重要。锂离子电池组的soc（state of charge，荷电状态）估算模型构建和精确估算值得获取，已成为其能量和安全管理的核心问题；锂离子电池组由具有高能量密度和闭路电压的钴酸锂电池单体组合构成，其安全性受到所处工作状态的影响；soe表征了锂电池在当前状态下可持续放出的能量，其单位为wh，为电池管理系统最基本也是最重要的的一个关键参数。此外，锂离子电池组的充放电过程包含复杂电能、化学能和热能转换等环节，过充电和过放电现象易引发安全事故，精确soe估算在防止过充电和过放电中起着重要作用；在锂离子电池组的应用中，其安全性依然是最为关注的问题，soe估算是其安全使用的基础和前提；锂离子电池组采用电池单体级联结构，满足了辅助动力供能过程中的容量和电压需求；然而，由于无法避免的材料和工艺差异，单体间不一致现象客观存在且无法避免；并且，该现象会随着循环次数的增加越来越明显，这就使得单体间不一致性的表达与修正成为成组soe估算的重要组成部分，同时也给成组soe精确估算带来了巨大的挑战。
4.针对soe估算的必要性和紧迫性需求，相关研究机构和高校，如麻省理工学院、宾州州立大学、美国南卡大学、英国利兹大学、英国罗伯特高登大学、美国国家可再生能源室、美国莱登能源公司、德国英飞凌科技公司、清华大学、北京航空航天大学、北京理工大学、北京交通大学、重庆大学、中国科学技术大学和哈尔滨工业大学等，针对soe估算展开了大量研究并进行了深入的探索。
5.目前对于soe的估算方法包括积分法、开路电压法、基于数据驱动的方法以及基于模型的方法等。其中，基于数据驱动的方法有神经网络算法、主元分析和支持向量机等方法；基于模型的方法包括卡尔曼滤波算法与粒子滤波等算法。基于数据驱动的方法往往需要大量精确的实验数据作为训练样本，计算量大，且不同工况下的训练结果只能用于相同的情况，泛化能力差。积分法计算方便，但对soe的初始值要求较高，抗干扰能力差，在存在干扰时造成误差后误差不能修复。


技术实现要素：

6.本发明的目的就在于提供一种基于遗忘因子自适应反馈修正的电池能量状态评估方法，以解决上述问题。为了实现上述目的，本发明采用的技术方案是这样的：一种基于遗忘因子自适应反馈修正的电池能量状态评估方法，包括下述步骤：s1：在线辨识锂离子电池模型参数，将锂离子电池模型离散化处理；s2：采用自适应遗忘因子递推最小二乘算法估算参数并对参数进行分离；s3：根据模型误差运用模糊逻辑控制器实现遗忘因子的自适应整定；s4: 根据分离得到的参数，计算状态空间变量的一步预测及其方差矩阵，根据状态方程可以直接得到如下所示预测方程，即进行先验估算；s5: 计算卡尔曼增益k(k)的表达式；s6:根据卡尔曼增益k(k)来修正初始预测值并更新过程噪音；得到该时刻的soe、极化电压以及系统残差，即进行后验估算。
7.作为优选的技术方案，步骤s1中，其方法为:
ꢀꢀ
（1）上式中，为开路电压；为工作电压输出观测变量；为系统输入控制变量；c1、c2、c3、c4和c5为相应的常数系数；为k 1时刻方程差值；下标k表示为第k时刻；下标k 1为第k 1时刻；下标k-1为第k-1时刻；下标k-2为第k-2时刻。
8.作为优选的技术方案，步骤s2中，其方法为：(2) (3)上式中，为待辨识系统的变量；为待辨识系统变量的估算值；为k时刻系统的观测噪声即噪声矩阵；为k时刻预测误差协方差矩阵；t为动力锂离子电池电压与电流的采样时间；为遗忘因子；为参数矩阵；下标k表示为第k时刻；下标k 1为第k 1时刻；下标k-1为第k-1时刻；下标k-2为第k-2时刻。
9.作为优选的技术方案，步骤s4中，其方法为：根据基尔霍夫电压定律和基尔霍夫电
流定律列写出锂离子电池能量状态观测方程和观测方程。
10.作为优选的技术方案，步骤s5中，其方法为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）上式中，为k 1时刻预测误差协方差矩阵；k
k 1
为k 1时刻卡尔曼增益；c为系统观测矩阵；为k时刻观测协方差；t为动力锂离子电池电压与电流的采样时间；下标k表示为第k时刻；下标k 1表示为第k 1时刻。
11.作为优选的技术方案，步骤s6中，其方法为： （5）上式中， 为k 1时刻观测变量预测误差；为k 1时刻状态变量预测量；为基于k时刻状态变量预测量预测的k 1时刻状态变量；u
k 1
为k 1时刻系统输入控制变量；c、d为系统观测矩阵; k
k 1
为k 1时刻卡尔曼增益；为k 1时刻误差协方差矩阵，即；为k时刻预测的k 1时刻误差协方差矩阵；e为单位矩阵；下标k表示为第k时刻；下标k 1表示为第k 1时刻；为k 1时刻观测变量。式中的与，一个是基于k即可的预测值，一个是k 1时刻的真值，一个是先验估计值，一个是后验估计值。
12.现有锂离子电池组bms应用中，基于瓦时积分和开路电压的soe估算方法，未能准确表征soe估算中存在的累积误差，并且不能结合当前状态进行参数修正；通过现有soe估算方法分析，基于affrls（遗忘因子递推最小二乘法）和ekf（扩展卡尔曼滤波算法）的锂电池soe估计研究，把闭路电压、电流为实时输入参量，在soe估算过程中考虑锂离子电池组的工况信息，克服了传统soe估算方法实时修正不足所造成的误差较大和逐渐累积等缺点；针对锂电池的极化特性表征难题，构建锂离子电池二阶rc等效电路模型；针对ffrls在线参数辨识方法不能很好适应工况变化难题，提出使用模糊逻辑控制器实现遗忘因子的在线自适应整定；在电池等效电路模型基础上运用affrls和ekf的算法实现锂离子电池组soe估算模型的建立和soe值的数学迭代运算；提出affrls和ekf的锂电池soe估计研究，实现了soe估算模型的构建与实验验证。
13.本文运用放电端电压残差和残差变化趋势作为模糊逻辑控制器的输入，设计模糊逻辑控制策略自适应在线调节ffrls的遗忘因子。封居强等人所使用的方法是基于ffrls和aekf的锂离子电池soc在线估计研究，其在线参数辨识算法ffrls不涉及遗忘因子的自适应调节，该文献使用的方法是在扩展卡尔曼（ekf）算法的基础上运用了传统的sage-husa自适应滤波方法对噪声估计特性进行自适应在线调节。
14.本发明针对锂离子电池成组应用时的soe估算问题，提出了一种锂离子电池组遗忘因子自适应反馈修正的电池能量状态评估方法，通过间歇式老化程度测定和实时校准计算处理过程，实现了锂离子电池成组soe估算的有效表征；基于遗忘因子自适应反馈修正的电池能量状态评估方法在容量归一化表征的基础上，通过老化状态对能量的影响系数计算获得老化因素影响的数学表达；基于遗忘因子自适应反馈修正的电池能量状态评估方法在定期测定校准基础上，通过额定能量状态与循环次数相关值的同步获取和修正，获得叠加循环次数修正的函数关系；基于遗忘因子自适应反馈修正的电池能量状态评估方法在老化影响系数和循环次数修正计算的基础上，结合两因素影响的叠加计算处理，获得老化过程对额定能量影响修正的计算方法；该方法在充分考虑锂离子电池成组工作基础上，结合soe估算的建立，实现对锂离子电池组老化过程特性的数学表达，构造基于遗忘因子自适应反馈修正的电池能量状态评估方法；与现有技术相比，本发明的优点在于：本发明结合锂离子电池组的状态空间模型，实现soe值的迭代计算，在用于跟踪锂离子电池组输出电压时，平均估算误差为0.01v，最大估算误差为0.05v，显著提高了电池soe估算精度。
附图说明
15.图1为本发明锂离子电池组soe估算模型结构示意图；图2为本发明根据模型误差运用模糊逻辑控制器实现遗忘因子的自适应整定；图3为采用不同方法的锂离子电池soe估算结果图；图4为采用不同方法的锂离子电池soe估算误差图。
具体实施方式
16.为了更好地解释本发明，在本实施例中仅以锂离子电池组为例进行说明，但本领域技术人员应该熟知，根据本发明的技术思想可以实现多种锂离子电池组的基于遗忘因子自适应反馈修正的电池能量状态评估方法；以下对锂离子电池组基于平方根扩展卡尔曼的锂离子电池组soc估算方法的实现步骤进行详细说明。
[0017] 针对提高soe估算精度目标，基于泰勒级数展开对锂离子电池组的非线性特征进行描述，以便于运用自适应遗忘因子递推最小二乘算法联合扩展卡尔曼滤波算法进行soe估计；结合锂离子电池组的状态空间模型，实现soe值的迭代计算，在用于跟踪锂离子电池组输出电压时，平均估算误差为0.01v，最大估算误差为0.05v；通过把soe作为其状态方程中的变量，输出闭路电压作为观测方程的变量，构建状态方程和观测方程表达式；soe(k)为状态变量，是k时刻的soe值；u
l
(k)为工作电压输出观测变量；状态方程系数 a为系统矩阵，b为控制输入矩阵；h为观测矩阵，初始值为[0 0 1]；系统噪声参数w(k)和观测噪声参数v(k)均为高斯白噪声，协方差分别为q和r；u
l
(k)为考虑测量误差v(k)影响的工作电压输出观测变量；通过迭代计算，从上一个状态值soe(k-1)、输入信号i(k)和工作电压输出观测变量u
l
(k)计算出卡尔曼模型的估算值soe(k)；利用无迹变换代替状态变量统计特性线性化变换，对于不同时刻的k值，具有高斯白噪声w(k)的随机向量soe和具有高斯白噪声v(k)的观测变量u
l
(k)构成离散时间非线性系统；通过把该估算框架应用于估算过程中，构建锂离
子电池组soe估算模型如图1所示。
[0018]
针对不同时刻k，该soe估算过程包括融合高斯白噪声的状态矩阵w(k)的随机状态变量soe，以及融入高斯白噪声的观测误差矩阵v(k)的工作电压输出观测变量u
l
(k)；f(*)是一个非线性状态方程，用于描述锂离子电池组的soe状态；g(*)是一个非线性观测方程，用于描述输出闭路电压的特征；噪声矩阵w(k)的方差使用q进行描述，噪声矩阵v(k)的方差使用r进行描述；在随机噪声的影响下，针对锂离子电池组soe精确估算目标，不同时刻k的估算通过以下步骤实现。
[0019]
图1为soe估算的全过程，包含锂离子电池参数获取过程，soe与极化电压计算过程；将扩展卡尔曼滤波算法用于极化电压估算过程，能降低初始极化电压未知时对锂离子电池soe估算带来的误差；具体步骤如下：s1：在线辨识锂离子电池模型参数，将锂离子电池模型离散化处理：
ꢀꢀ
（1）上式中，为开路电压；为工作电压输出观测变量；为系统输入控制变量；c1、c2、c3、c4和c5为相应的常数系数；为k 1时刻方程差值；下标k表示为第k时刻；下标k 1表示为第k 1时刻；下标k-1表示为第k-1时刻；s2：采用模糊逻辑自适应遗忘因子递推最小二乘算法估算参数并对参数进行分离：(2) (3)上式中，为待辨识系统的变量；为待辨识系统变量的估算值；为k时刻系统的观测噪声即噪声矩阵；为k时刻预测误差协方差矩阵；t为动力锂离子电池电压与电流的采样时间；为遗忘因子；为参数矩阵；下标k表示为第k时刻；下标k 1表示为第k 1时刻；下标k-1表示为第k-1时刻；s3：根据模型误差运用模糊逻辑控制器实现遗忘因子的自适应整定，如图2所示。
[0020]
s4：根据分离得到的参数，计算状态空间变量的一步预测及其方差矩阵，根据基尔
霍夫电压定律和基尔霍夫电流定律列写出锂离子电池能量状态观测方程和观测方程。
[0021]
s5：计算卡尔曼增益k(k)的表达式： （4）上式中，为k 1时刻预测误差协方差矩阵；k
k 1
为k 1时刻卡尔曼增益；c为系统观测矩阵；为k时刻观测协方差；下标k表示为第k时刻；下标k 1表示为第k 1时刻；s6：根据卡尔曼增益k(k)来修正初始预测值并更新过程噪音；得到该时刻的soe、极化电压以及系统残差，即进行后验估算：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（5）上式中， 为k 1时刻观测变量预测误差；为k 1时刻状态变量预测量；为基于k时刻状态变量预测量预测的k 1时刻状态变量；u
k 1
为k 1时刻系统输入控制变量；c、d为系统观测矩阵; k
k 1
为k 1时刻卡尔曼增益；为k 1时刻误差协方差矩阵；为k时刻预测的k 1时刻误差协方差矩阵；e为单位矩阵；下标k表示为第k时刻；下标k 1表示为第k 1时刻；为k 1时刻观测变量。
[0022]
在锂离子电池soe估算过程中，其迭代过程如图1所示，通过以上一系列公式进行迭代，输出的xk为估算的soe与极化电压值，通过使用泰勒公式展开来解决非线性转换问题，通过扩展卡尔曼滤波算法实现极化电压估算；通过以上迭代计算过程，基于自适应遗忘因子递推最小二乘（affrls）和扩展卡尔曼滤波算法（ekf）的锂电池能量状态（soe）估算方法，实现了锂离子电池组的soe的估算模型构建。
[0023]
采用上述实施例的自适应遗忘因子递推最小二乘结合扩展卡尔曼滤波算法（affrls-ekf）法与现有的遗忘因子递推最小二乘结合扩展卡尔曼滤波算法（ffrls-ekf）法的锂离子电池soe估算结果如图3所示，在图3中real为实验室条件下测得的能量状态真实值，从图3中可以看出；在估算过程中，affrls-ekf法与ffrls-ekf法均可以是实现锂离子电池soe的有效在线估算，affrls-ekf法的估算结果曲线相比于ffrls-ekf法更接近真实值，表明affrls-ekf法相比ffrls-ekf法有更高的估算精度。
[0024]
采用上述实施例的affrls-ekf法与现有的ffrls-ekf法的锂离子电池soe估算误差如图4所示，从图4中可以看出：ffrls-ekf法的估算最大误差在6%左右，再放电末期误差增大到7.5%左右；affrls-ekf法估算误差收敛在3%以内，表明affrls-ekf法有更好的鲁棒性和估算精度。
[0025]
本发明的以上实施例仅以锂离子电池组为例进行了基于遗忘因子自适应反馈修
正的电池能量状态评估方法的说明，但可以理解的是，在不脱离本发明精神和范围下本领域技术人员可以对其进行任意的改变和变化。