一种基于变分贝叶斯滤波的电池荷电状态估计方法技术

本发明专利技术提供了一种基于变分贝叶斯滤波的电池荷电状态估计方法，通过采集静置的电池端电压标定OCV与SOC关系曲线，建立系统状态方程和量测方程，建立电池模型参数的状态空间方程和量测方程，利用双变分贝叶斯滤波算法在线联合估计电池的SOC和电池模型参数。本发明专利技术通过交替使用两个变分贝叶斯滤波器在线估计电池状态和参数，同时实时估计量测噪声协方差来补偿电池模型的不确定因素，使电池模型能较好的描述电池的动态特性，具有了较好的适应性，因而具有精度高、实时性好、鲁棒性强等优点，实验证明，本发明专利技术的SOC估计精度与鲁棒性都优于传统的双扩展卡尔曼滤波算法。

A Battery Charge State Estimation Method Based on Variational Bayesian Filtering

The invention provides a battery charging state estimation method based on variational Bayesian filtering. The relationship curve between OCV and SOC is calibrated by collecting static battery terminal voltage, the system state equation and measurement equation are established, the state space equation and measurement equation of battery model parameters are established, and the battery SOC and battery model parameters are estimated online jointly by using double variational Bayesian filtering algorithm. The invention uses two variational Bayesian filters alternately to estimate the battery state and parameters online, and estimates the covariance of measurement noise in real time to compensate for the uncertainties of the battery model, so that the battery model can better describe the dynamic characteristics of the battery, has better adaptability, and has the advantages of high precision, good real-time performance and strong robustness. The experiment proves that the invention has the advantages of high accuracy, good real-time performance and strong robustness. The accuracy and robustness of SOC estimation are better than the traditional double extended Kalman filter algorithm.

【技术实现步骤摘要】
一种基于变分贝叶斯滤波的电池荷电状态估计方法
本专利技术涉及电池荷电领域，尤其是一种电池荷电状态估计方法。
技术介绍
电池荷电状态(stateofcharge，SOC)表征了电池的剩余可用电量，对其准确估计是电池管理系统的核心功能之一，是实现电池均衡、故障诊断、电池充/放电控制等其他功能的前提，其估计精度直接影响整个电池管理系统的效率。然而，SOC通常难以直接测量获取，只能基于其它可测参数(电压、电流、温度)间的关系进行估算。由于这种关系存在强烈的非线性特性，而且往往会受到工况、温度、老化程度等诸多因素影响而实时地发生变化，所以SOC的精准在线估计十分困难。目前常用的SOC估计方法分为四类：开路电压法、安时积分法、数据驱动的方法和基于模型的方法。现有方法的不足之处在于：1)开路电压法需要长时间静置电池，不适于实时估计。2)安时积分法对准确的SOC初值依赖性强，对电流采集设备的精度要求很高。若SOC初值不准确或者电流表发生偏移，则会产生较大的累积误差。3)神经网络法、模糊逻辑法和支持向量机法等数据驱动类方法，需要大量的实验数据作为先验知识，若样本数据不能较为全面地反映电池特性则估计精度也难以保证。4)基于模型的方法，以扩展卡尔曼滤波(EKF)及其衍生算法为主。基于模型的方法克服了安时积分法对SOC初值要求严格的问题，同时并不需要大量的训练数据，具有研究与发展前景。但其精度受电池模型、传感器测量误差等因素的影响，在实际应用过程中，电池模型参数随SOC、温度以及电池老化程度的不同而变化，而且传感器误差的统计信息可能未知或时变，这将导致传统的EKF类方法的估计精度不高、鲁棒性差。因而需要开发鲁棒且高精度的SOC估计算法。
技术实现思路
为了克服现有技术的不足，本专利技术提供基于变分贝叶斯滤波方法的电池荷电状态估计方法，通过利用两个变分贝叶斯滤波器来实时估计电池的SOC、电池模型参数以及量测噪声协方差，具有精度高、实时性好、鲁棒性强等优点。本专利技术解决其技术问题所采用的技术方案包括以下步骤：步骤1：选择二阶RC网络作为电池等效电路模型，通过电池混合脉冲功率特性(HPPC)实验对电池模型的各个参数进行辨识，包括电池的欧姆内阻、两个极化电阻和两个极化电容；步骤2：通过采集静置的电池端电压标定开路电压(OCV)与SOC关系曲线，采集步骤为：(1)：标准充电：以0.5C电流对电池持续充电至截止电压；(2)：静置：将电池静置2分钟；(3)：持续充电：以0.2C电流对电池持续充电至截止电压；(4)：静置：将电池静置3小时，测量电池此时端电压，此时电池SOC为100％；(5)：间断放电：以0.2C电流对电池放电，直至电池SOC下降10％；(6)：静置：将电池静置3小时，测量电池此时端电压；(7)：重复步骤(5)-(6)九次，直至SOC下降到0％；步骤3：建立系统状态方程和量测方程，表示如下：式中，k表示时刻，Xk表示k时刻系统状态变量，Xk＝[SOCk,U1,k,U2,k]T，SOCk为k时刻的电池荷电状态，U1,k和U2,k为两个RC并联电路在k时刻的端电压；Ik为电池充放电电流；θk为时变的电池模型参数向量，θk＝[R0,R1,τ1,R2,τ2]T，其中，R0为电池的欧姆内阻，R1、R2分别表示电池模型中两个RC并联电路的内阻，τ1＝R1C1、τ2＝R2C2分别表示电池模型中两个RC并联电路的时间常数；yk表示k时刻的模型输出量，这里为电池端电压Uk；wk为系统过程噪声，服从零均值，协方差为的高斯分布；vk为量测噪声，服从零均值，方差为Σx,k的高斯分布，Σx,k是未知的，其后验分布用逆伽马分布来表示，即p(Σx,k|y1:k)＝Inv-Gamma(Σx,k|αx,k,βx,k)，其中αx,k和βx,k为逆伽马分布的参数；f(·)和h(·)分别表示关于状态变量X、输入量I以及模型参数变量θ的非线性函数，其数学表达式为：h(·)＝Uk＝UOC(SOCk)-U1,k-U2,k-IkR0其中，η为库伦效率，Δt为采样周期，Qmax为电池额定容量，UOC为电池开路电压；利用步骤2的OCV与SOC关系曲线，通过最小二乘法曲线拟合得到开路电压UOC与SOC的关系函数UOC(SOCk)；步骤4：根据电池模型参数缓慢变化的特点，建立电池模型参数的状态空间方程和量测方程如下：式中，rk表示协方差为的零均值高斯白噪声，反映了电池参数随时间的缓慢变化；量测方程dk是关于θk的量测方程，白噪声ek表征了传感器噪声和模型的不确定因素，其方差Σθ,k未知，Σθ,k的后验分布用逆伽马分布来表示，即p(Σθ,k|y1:k)＝Inv-Gamma(Σθ,k|αθ,k,βθ,k)，其中αθ,k和βθ,k为逆伽马分布的参数；步骤5：利用双变分贝叶斯(DualVariationalBayes，DVB)滤波算法在线联合估计电池的SOC和电池模型参数，具体步骤如下：1)初始化如下参数：X和θ的均值分别为和协方差Px,0和Pθ,0，过程噪声协方差以及量测噪声协方差概率分布参数其中，电池模型参数的初始值在步骤1中通过电池充放电实验得到，其余参数根据实际应用中电池的初始荷电状态和测量设备精度进行设置；2)计算系统状态估计、电池参数估计及其协方差，以及量测噪声协方差概率分布参数的预测值：其中，表示k-1时刻的系统状态估计值；表示k-1时刻的电池参数估计值；Px,k-1表示k-1时刻系统状态估计协方差；Pθ,k-1表示k-1时刻电池参数估计协方差；和表示k-1时刻量测噪声方差概率分布参数的估计值。(·)-表示相应变量的一步预测值，ρx和ρθ为衰减因子，ρx∈(0,1]，ρθ∈(0,1]；3)计算系统状态估计及其协方差Px,k、电池参数估计及其协方差Pθ,k，以及量测噪声协方差概率分布参数的测量更新；通过N次循环迭代实现，主要包括以下步骤：第一步，初始化，令第二步，计算量测协方差的估计值：第三步，计算系统状态估计及其协方差的测量更新：其中，为量测方程关于系统状态Xk的雅克比矩阵；第四步，计算电池参数估计及其协方差的测量更新：其中，Hθ,k为量测方程关于电池参数θk的雅克比矩阵，表示为：第五步，计算量测噪声协方差概率分布参数的测量更新：第六步，将第二步～第五步循环迭代N次，N取2～10；第七步，将第N次循环迭代得到的估计值作为k时刻的最终估计值，即：令得到的估计值即为k时刻的SOC估计和电池参数估计，SOC的估计值反映了电池剩余电量的估计，从而为电池管理系统提供参数。本专利技术有益效果在于通过交替使用两个变分贝叶斯滤波器在线估计电池状态和参数，同时实时估计量测噪声协方差来补偿电池模型的不确定因素，使电池模型能较好的描述电池的动态特性，具有了较好的适应性，因而具有精度高、实时性好、鲁棒性强等优点。实验证明，本专利技术的SOC估计精度与鲁棒性都优于传统的双扩展卡尔曼滤波算法。附图说明图1为本专利技术的电池二阶等效电路模型。图2为本专利技术的方法流程图。图3为本专利技术的SOC0＝0.8时，电池恒流放电工况下的SOC估计情况。图4为本专利技术的SOC0＝0.8且设定的量测噪声协方差偏离实际值时，电池恒流放电工况下的SOC估计情况。具体实施方式下面结合附图和实施例对本专利技术进一步说明。步骤1：为了更准确地模拟电池本文档来自技高网...

【技术保护点】
1.一种基于变分贝叶斯滤波的电池荷电状态估计方法，其特征在于包括下述步骤：步骤1：选择二阶RC网络作为电池等效电路模型，通过电池混合脉冲功率特性实验对电池模型的各个参数进行辨识，包括电池的欧姆内阻、两个极化电阻和两个极化电容；步骤2：通过采集静置的电池端电压标定开路电压OCV与SOC关系曲线，采集步骤为：(1)：标准充电：以0.5C电流对电池持续充电至截止电压；(2)：静置：将电池静置2分钟；(3)：持续充电：以0.2C电流对电池持续充电至截止电压；(4)：静置：将电池静置3小时，测量电池此时端电压，此时电池SOC为100％；(5)：间断放电：以0.2C电流对电池放电，直至电池SOC下降10％；(6)：静置：将电池静置3小时，测量电池此时端电压；(7)：重复步骤(5)‑(6)九次，直至SOC下降到0％；步骤3：建立系统状态方程和量测方程，表示如下：

【技术特征摘要】
1.一种基于变分贝叶斯滤波的电池荷电状态估计方法，其特征在于包括下述步骤：步骤1：选择二阶RC网络作为电池等效电路模型，通过电池混合脉冲功率特性实验对电池模型的各个参数进行辨识，包括电池的欧姆内阻、两个极化电阻和两个极化电容；步骤2：通过采集静置的电池端电压标定开路电压OCV与SOC关系曲线，采集步骤为：(1)：标准充电：以0.5C电流对电池持续充电至截止电压；(2)：静置：将电池静置2分钟；(3)：持续充电：以0.2C电流对电池持续充电至截止电压；(4)：静置：将电池静置3小时，测量电池此时端电压，此时电池SOC为100％；(5)：间断放电：以0.2C电流对电池放电，直至电池SOC下降10％；(6)：静置：将电池静置3小时，测量电池此时端电压；(7)：重复步骤(5)-(6)九次，直至SOC下降到0％；步骤3：建立系统状态方程和量测方程，表示如下：式中，k表示时刻，Xk表示k时刻系统状态变量，Xk＝[SOCk,U1,k,U2,k]T，SOCk为k时刻的电池荷电状态，U1,k和U2,k为两个RC并联电路在k时刻的端电压；Ik为电池充放电电流；θk为时变的电池模型参数向量，θk＝[R0,R1,τ1,R2,τ2]T，其中，R0为电池的欧姆内阻，R1、R2分别表示电池模型中两个RC并联电路的内阻，τ1＝R1C1、τ2＝R2C2分别表示电池模型中两个RC并联电路的时间常数；yk表示k时刻的模型输出量，这里为电池端电压Uk；wk为系统过程噪声，服从零均值，协方差为的高斯分布；vk为量测噪声，服从零均值，方差为Σx,k的高斯分布，Σx,k是未知的，其后验分布用逆伽马分布来表示，即p(Σx,k|y1:k)＝Inv-Gamma(Σx,k|αx,k,βx,k)，其中αx,k和βx,k为逆伽马分布的参数；f(·)和h(·)分别表示关于状态变量X、输入量I以及模型参数变量θ的非线性函数，其数学表达式为：h(·)＝Uk＝UOC(SOCk)-U1,k-U2,k-IkR0其中，η为库伦效率，Δt为采样周期，Qmax为电池额定容量，UOC为电池开路电压；利用步骤2的OCV与SOC关系曲线，通过最小二乘法曲线拟合得到开路电压UOC与SOC的关系函数UOC(SOCk)；步...

【专利技术属性】
技术研发人员：侯静，羊彦，高田，张一凡，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61