CN112487597A - 一种基于改进ekf算法的动力电池soc估算方法 - Google Patents
一种基于改进ekf算法的动力电池soc估算方法 Download PDFInfo
- Publication number
- CN112487597A CN112487597A CN201910806080.2A CN201910806080A CN112487597A CN 112487597 A CN112487597 A CN 112487597A CN 201910806080 A CN201910806080 A CN 201910806080A CN 112487597 A CN112487597 A CN 112487597A
- Authority
- CN
- China
- Prior art keywords
- battery
- model
- capacity
- soc
- soc estimation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/70—Energy storage systems for electromobility, e.g. batteries
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tests Of Electric Status Of Batteries (AREA)
- Secondary Cells (AREA)
Abstract
本发明提出了一种基于改进EKF算法的动力电池SOC估算方法。采用了PNGV电池等效模型,改进了模型精度,并以扩展卡尔曼滤波算法为基础,引入了与电池表面温度有关的容量校正模型以及和库伦效率有关的电流效率模型,更准确地表征电池的动态特性,实现对SOC估算的校正。本发明具有估算精度较高、误差较小、且具有来回震荡的调节过程、抗干扰能力较强的优点,能够使SOC估算在整个充放电区间里和较为复杂的使用条件下保持较高的精度。
Description
技术领域
本发明属于电池管理系统技术领域,具体涉及一种基于改进EKF算法的动力电池SOC估算方法。
背景技术
动力电池是电动汽车的动力源,是电动汽车的重要组成部分。动力电池及管理系统是制约电动汽车发展的核心技术之一,其中SOC估算是电池管理系统的关键技术,准确地估算动力电池SOC,不仅可以提高动力电池的使用性能,还可以延长动力电池的使用寿命。
SOC是电池管理系统最重要的参数之一,电池管理系统的管理和控制都以准确估算SOC为前提和基础。目前,卡尔曼滤波算法是应用最广泛的SOC估算方法,它是一种递推的线性最小估计算法,在系统噪声统计特性已知的情况下,其估算精度较高。传统的卡尔曼滤波是在标准条件下获得的,是一种无偏的线性最小方差估计算法,在动态系统的数学模型和噪声统计特性已知的情况下,卡尔曼滤波通过测量值对预测估计进行修正,可以得到状态的精确估计。但是由于动力电池是一个高度非线性系统,在实际应用中很难得到精确的数学模型和噪声的统计特性,所以实际中通常利用扩展卡尔曼滤波算法(EKF)估算SOC。
目前现有技术的SOC估算方法主要存在模型过于简化和噪声统计特性不准确等问题。为满足运算速度的要求,大多数情况下都采用简化模型,其中的模型参数往往是固定不变的,而动力电池系统是一个高度非线性系统,其容量、内阻、充放电效率等参数都会随着温度、充放电电流的改变而发生变化,这就可能造成应用模型与实际系统之间的不匹配,模型精度降低,从而导致SOC估算误差增大。因此,设计一种改进的SOC估算方法,在估算过程中充分考虑电池温度和充放电倍率对电池特性的影响,更准确地表征动力电池的动态特性,对提高SOC估算精度具有重大的现实意义。
发明内容
本发明提出了一种基于改进EKF算法的动力电池SOC估算方法,在EKF算法基础上加入与电池表面温度有关的容量校正模型以及和库伦效率有关的电流效率模型,更准确地表征动力电池的动态特性,提高SOC估算精度。
本发明的具体实施技术方案如下:
本发明提供一种基于改进EKF算法的动力电池SOC估算方法,包括如下步骤:
步骤S1:建立动力电池等效模型;
步骤S2:建立与电池表面温度有关的容量校正模型;
步骤S3:建立和库伦效率有关的电流效率模型;
步骤S4:以EKF算法为基础,并结合S2和S3,采用改进的EKF算法对动力电池SOC进行估算。
更进一步的,所述步骤S1的具体过程包括:
采用PNGV电池等效模型,并在模型上添加一个RC环节进行改进,使模型更加符合电池的内部结构和动态变化情况,并且更精确。改进模型的数学关系式为:
UL(t)=Uocv(SOC,t)-IL(t)R0-Up1(t)-Up2(t) (1)
Up(t)=Up1(t)+Up2(t)=U1exp(-t/τ1)+U2exp(-t/τ2) (2)
其中UL(t)为电池端电压,Uocv(SOC,t)为开路电压,R0为电池的欧姆内阻,IL为电池环路的电流,τ1、τ2分别为两个RC环节的时间常数。
更进一步的,所述步骤S2的具体过程包括:
将电池分别置于-20、-10、0、25、40、60℃的温度下,以固定放电倍率对电池进行放电到截止电压,通过实验得到各温度下电池的实际容量,通过与标称容量比较,得到不同温度下的容量补偿系数ηT:
ηT=Qi/QT (3)
上式中:Qi为实际容量,QT为标称容量;
将所测实际容量与对应温度进行曲线拟合,得到电池内部温度对电池容量校正拟合曲线;
由于电池内部温度和电池极柱温度有差异,对所测电池温度进行自校,并对所测数据进行曲线拟合,得到温度自校系数拟合曲线。
更进一步的,所述步骤S3的具体过程包括:
设置充、放电条件:在25℃条件下,充电倍率为0.5C,分别以0.1C、0.5C、1C、1.5C、2C、3C放电倍率进行放电至截止电压,得到不同放电倍率下的容量补偿系数,即库仑效率ηI:
ηI=Qd/Qc (4)
式中:Qd为实际放电容量,Qc为实际充电容量;
将测量点数据进行曲线拟合,得到库伦效率拟合曲线。
更进一步的,所述步骤S4的具体过程包括:
S41:建立电池系统离散状态空间模型
根据改进的PNGV电池等效模型的数学方程式,联合电池安时积分法原理,将电池的SOC与PNGV电池模型的极化电压作为电池的状态变量,选取测量的电池端电压UL作为观测量,并结合S2、S3的容量校正模型和电流效率模型,得到状态方程与观测方程如式(5)和(6)所示:
UL,k=Uocv(SOCk)-IL,kR0-Up1,k-Up2,k (6)
令:
S42:基于改进EKF算法的SOC估算过程:
1)k=0,选择初始值X0,0,P0,0;
2)计算预测值Xk,k-1
Xk,k-1=Ak-1Xk-1,k-1+Bk-1Ik-1
3)预测协方差矩阵Pk,k-1
4)更新卡尔曼增益矩阵Kk
5)更新状态参数Xk,k
Xk,k=Xk,k-1+Kk[Uk-(Ak-1Xk-1,k-1+Bk-1Ik-1)]
6)更新协方差矩阵Pk,k
Pk,k=(I-KkCk)Pk,k-1
7)判断滤波过程是否完成,如果否,令k=k+1,返回步骤2);如果是,结束算法。
由上述本发明的技术方案可以看出,本发明具有如下技术效果:
在EKF算法的基础上,本发明充分考虑电池温度和充放电倍率对SOC估算的影响,引入与电池表面温度有关的容量校正模型以及和库伦效率有关的电流效率模型,更准确地表征动力电池的动态特性,相比传统的EKF算法具有更高的估算精度。
附图说明
图1为本发明的电池等效模型示意图;
图2为本发明的SOC算法流程图;
图3为25℃,1C充放电倍率下SOC估算结果对比;
图4为25℃,3C充放电倍率下SOC估算结果对比;
图5为-5℃,3C充放电倍率下SOC估算结果对比;
具体实施方式
下面结合附图,通过实施例对本发明作进一步的说明。
本发明提供一种基于改进EKF算法的动力电池SOC估算方法,包括如下步骤:
步骤S1:建立动力电池等效模型。
考虑到电池模型的动态特性及复杂性,本发明采用PNGV电池等效模型,并在模型上添加一个RC环节进行改进,使模型更加符合电池的内部结构和动态变化情况,并且更精确,其模型结构图如图1所示。
其中UL为电池端电压,Uocv为开路电压,R0为电池的欧姆内阻,Rp1、Rp2和Cp1、Cp2为电池两个RC环节的内阻和电容,两个RC环节反映了电池内部的极化现象,IL为电池环路的电流,模型的数学关系式如下:
UL(t)=Uocv(SOC,t)-IL(t)R0-Up1(t)-Up2(t) (1)
Up(t)=Up1(t)+Up2(t)=U1exp(-t/τ1)+U2exp(-t/τ2) (2)
其中τ1、τ2分别为两个RC环节的时间常数,τ1=Rp1Cp1,τ2=Rp2Cp2。
步骤S2:建立与电池表面温度有关的容量校正模型。
将电池分别置于-20、-10、0、25、40、60℃的温度下,以固定放电倍率对电池进行放电到截止电压,通过实验得到各温度下电池的实际容量,通过与标称容量比较,得到不同温度下的容量补偿系数ηT:
ηT=Qi/QT (3)
上式中:Qi为实际容量,QT为标称容量;
将所测实际容量与对应温度进行曲线拟合,得到电池内部温度对电池容量校正拟合曲线;
由于电池内部温度和电池极柱温度有差异,对所测电池温度进行自校,并对所测数据进行曲线拟合,得到温度自校系数拟合曲线。
步骤S3:建立和库伦效率有关的电流效率模型。
设置充、放电条件:在25℃条件下,充电倍率为0.5C,分别以0.1C、0.5C、1C、1.5C、2C、3C放电倍率进行放电至截止电压,得到不同放电倍率下的容量补偿系数,即库仑效率ηI:
ηI=Qd/Qc (4)
式中:Qd为实际放电容量,Qc为实际充电容量;
将测量点数据进行曲线拟合,得到库伦效率拟合曲线。
步骤S4:以EKF算法为基础,并结合S2和S3,采用改进的EKF算法对动力电池SOC进行估算。
S41:建立电池系统离散状态空间模型
根据改进的PNGV电池等效模型的数学方程式,联合电池安时积分法原理,将电池的SOC与PNGV电池模型的极化电压作为电池的状态变量,选取测量的电池端电压UL作为观测量,并结合S2、S3的容量校正模型和电流效率模型,得到状态方程与观测方程如式(5)和(6)所示:
UL,k=Uocv(SOCk)-IL,kR0-Up1,k-Up2,k (6)
令:
S42:基于改进EKF算法的SOC估算过程,其流程图如图2所示。
1)k=0,选择初始值X0,0,P0,0;
2)计算预测值Xk,k-1
Xk,k-1=Ak-1Xk-1,k-1+Bk-1Ik-1
3)预测协方差矩阵Pk,k-1
4)更新卡尔曼增益矩阵Kk
5)更新状态参数Xk,k
Xk,k=Xk,k-1+Kk[Uk-(Ak-1Xk-1,k-1+Bk-1Ik-1)]
6)更新协方差矩阵Pk,k
Pk,k=(I-KkCk)Pk,k-1
7)判断滤波过程是否完成,如果否,令k=k+1,返回步骤2);如果是,结束算法。
为了验证本发明提出的SOC估算算法的准确性,在Matlab/Simulink中建立算法仿真模型,对改进前后的SOC估算算法模型进行仿真对比分析。为考虑充、放电倍率对SOC估算的影响,选择在室温(25℃)下分别对电池以1C、3C进行恒流放电,静置一段时间后,再分别以1C、3C进行恒流充电、放电,SOC估算结果如图3、图4所示。其中参考曲线为通过安时积分法得到的曲线,因充、放电过程中电流均保持不变,故以这条曲线为参考曲线,且赋予SOC初值SOC0为0.8。曲线2是经PNGV电池模型+未改进EKF估算的SOC曲线;曲线3是经改进EKF算法估算的SOC曲线。
从不同充、放电倍率下SOC估算结果对比可以看出,随着充放电倍率的增大,未改进EKF估算的SOC距离参考曲线的误差逐渐增大,最高达到3.5%左右,随着时间的延长,误差一直存在。而无论是1C充放电倍率,还是3C充放电倍率,改进EKF算法估算的SOC曲线始终在参考曲线上下波动,误差明显较小。在两种充放电倍率下,改进EKF估算的SOC值比未改进EKF估算的SOC值精确度要高,体现了本发明算法的有效性。
为验证温度对SOC估算的影响,选择在3C充放电倍率下,分别以-5℃、25℃进行充放电,其结果如图4、图5所示。通过图4和图5可以明显看出,温度对SOC估算影响很大。-5℃下未改进EKF估算的SOC曲线与参考曲线最大误差达到近4%,而改进EKF算法得到的SOC曲线误差较小,最大约为0.9%,还存在围绕参考曲线来回波动的调节过程,说明改进EKF在极端环境下SOC估算精度更高,算法鲁棒性较好。
Claims (5)
1.一种基于改进EKF算法的动力电池SOC估算方法,其特征在于:包括如下步骤:
步骤S1:建立动力电池等效模型;
步骤S2:建立与电池表面温度有关的容量校正模型;
步骤S3:建立和库伦效率有关的电流效率模型;
步骤S4:以EKF算法为基础,并结合S2和S3,采用改进的EKF算法对动力电池SOC进行估算。
2.根据权利要求1所述的一种基于改进EKF算法的动力电池SOC估算方法,其特征在于,所述步骤S1的具体过程包括:
采用PNGV电池等效模型,并在模型上添加一个RC环节进行改进,使模型更加符合电池的内部结构和动态变化情况,并且更精确。改进模型的数学关系式为:
UL(t)=Uocv(SOC,t)-IL(t)R0-Up1(t)-Up2(t) (1)
Up(t)=Up1(t)+Up2(t)=U1exp(-t/τ1)+U2exp(-t/τ2) (2)
其中UL(t)为电池端电压,Uocv(SOC,t)为开路电压,R0为电池的欧姆内阻,IL为电池环路的电流,τ1、τ2分别为两个RC环节的时间常数。
3.根据权利要求1所述的一种基于改进EKF算法的动力电池SOC估算方法,其特征在于,所述步骤S2的具体过程如下:
将电池分别置于-20、-10、0、25、40、60℃的温度下,以固定放电倍率对电池进行放电到截止电压,通过实验得到各温度下电池的实际容量,通过与标称容量比较,得到不同温度下的容量补偿系数ηT:
ηT=Qi/QT (3)
上式中:Qi为实际容量,QT为标称容量;
将所测实际容量与对应温度进行曲线拟合,得到电池内部温度对电池容量校正拟合曲线;
由于电池内部温度和电池极柱温度有差异,对所测电池温度进行自校,并对所测数据进行曲线拟合,得到温度自校系数拟合曲线。
4.根据权利要求1所述的一种基于改进EKF算法的动力电池SOC估算方法,其特征在于,所述步骤S3的具体过程如下:
设置充、放电条件:在25℃条件下,充电倍率为0.5C,分别以0.1C、0.5C、1C、1.5C、2C、3C放电倍率进行放电至截止电压,得到不同放电倍率下的容量补偿系数,即库仑效率ηI:
ηI=Qd/Qc (4)
式中:Qd为实际放电容量,Qc为实际充电容量;
将测量点数据进行曲线拟合,得到库伦效率拟合曲线。
5.根据权利要求1所述的一种基于改进EKF算法的动力电池SOC估算方法,其特征在于,所述步骤S4的具体过程如下:
S41:建立电池系统离散状态空间模型
根据改进的PNGV电池等效模型的数学方程式,联合电池安时积分法原理,将电池的SOC与PNGV电池模型的极化电压作为电池的状态变量,选取测量的电池端电压UL作为观测量,并结合S2、S3的容量校正模型和电流效率模型,得到状态方程与观测方程如式(5)和(6)所示:
UL,k=Uocv(SOCk)-IL,kR0-Up1,k-Up2,k (6)
令:
S42:基于改进EKF算法的SOC估算过程:
1)k=0,选择初始值X0,0,P0,0;
2)计算预测值Xk,k-1
Xk,k-1=Ak-1Xk-1,k-1+Bk-1Ik-1
3)预测协方差矩阵Pk,k-1
4)更新卡尔曼增益矩阵Kk
5)更新状态参数Xk,k
Xk,k=Xk,k-1+Kk[Uk-(Ak-1Xk-1,k-1+Bk-1Ik-1)]
6)更新协方差矩阵Pk,k
Pk,k=(I-KkCk)Pk,k-1
7)判断滤波过程是否完成,如果否,令k=k+1,返回步骤2);如果是,结束算法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910806080.2A CN112487597B (zh) | 2019-08-20 | 2019-08-20 | 一种基于改进ekf算法的动力电池soc估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910806080.2A CN112487597B (zh) | 2019-08-20 | 2019-08-20 | 一种基于改进ekf算法的动力电池soc估算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112487597A true CN112487597A (zh) | 2021-03-12 |
| CN112487597B CN112487597B (zh) | 2023-05-05 |
Family
ID=74920158
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910806080.2A Expired - Fee Related CN112487597B (zh) | 2019-08-20 | 2019-08-20 | 一种基于改进ekf算法的动力电池soc估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112487597B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113919145A (zh) * | 2021-09-27 | 2022-01-11 | 珠海格力电器股份有限公司 | 电池soc的估算方法、装置、电池、存储介质及处理器 |
| CN115754742A (zh) * | 2022-12-07 | 2023-03-07 | 昆明理工大学 | 储能系统的soc估算方法及装置、设备及存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018096953A (ja) * | 2016-12-16 | 2018-06-21 | 三菱自動車工業株式会社 | 電池状態推定装置 |
| CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
-
2019
- 2019-08-20 CN CN201910806080.2A patent/CN112487597B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018096953A (ja) * | 2016-12-16 | 2018-06-21 | 三菱自動車工業株式会社 | 電池状態推定装置 |
| CN108490356A (zh) * | 2018-03-06 | 2018-09-04 | 福建工程学院 | 一种改进ekf算法的锂电池soc估算方法 |
Non-Patent Citations (3)
| Title |
|---|
| 刘晓悦,等: "锂离子电池组SOC估算方法", 《计算机仿真》 * |
| 樊波等: "基于改进PNGV模型的电池SOC估计算法研究", 《电测与仪表》 * |
| 雷津等: "基于EKF的船舶动力锂电池剩余容量估算", 《电子测量技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113919145A (zh) * | 2021-09-27 | 2022-01-11 | 珠海格力电器股份有限公司 | 电池soc的估算方法、装置、电池、存储介质及处理器 |
| CN115754742A (zh) * | 2022-12-07 | 2023-03-07 | 昆明理工大学 | 储能系统的soc估算方法及装置、设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112487597B (zh) | 2023-05-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107402353B (zh) | 一种对锂离子电池的荷电状态进行滤波估计的方法及系统 | |
| CN111060834A (zh) | 一种动力电池健康状态估算方法 | |
| CN111985154B (zh) | 一种自适应模糊卡尔曼估计soc算法 | |
| CN111323705A (zh) | 基于鲁棒递归最小二乘的电池参数辨识方法及系统 | |
| CN112327182B (zh) | 基于量测值残差序列的自适应h无穷滤波soc估计方法 | |
| WO2025016154A1 (zh) | 一种电池soc估算方法、装置及储能设备 | |
| WO2024152592A1 (zh) | 电池soc估计方法、装置、设备、电池模组及存储介质 | |
| CN115656848A (zh) | 一种基于容量修正的锂电池soc估算方法 | |
| CN113125969B (zh) | 基于aukf的电池数据处理方法、设备和介质 | |
| CN110907835B (zh) | 一种具有噪声免疫特性的电池模型参数辨识与soc估计方法 | |
| CN112964992A (zh) | 基于aukf的电池内部温度信息处理方法、设备和介质 | |
| CN109839596B (zh) | 基于ud分解的自适应扩展卡尔曼滤波的soc估算方法 | |
| CN107153163A (zh) | 一种基于自适应ukf的锂电池soc估算方法 | |
| CN112433154A (zh) | 基于ffrls和ekf的锂离子电池soc估测算法 | |
| CN115248381B (zh) | 基于安时校正的扩展卡尔曼滤波锂电池soc估算方法 | |
| CN109828215A (zh) | 一种提升电池单体soc估算精度的方法和系统 | |
| CN112415412A (zh) | 估算电池soc值的方法和装置及车辆、存储介质 | |
| CN112601968A (zh) | 用于48v轻混汽车锂离子电池的荷电状态在线估算系统 | |
| CN112816873A (zh) | 一种基于改进电池模型的iekf锂电池soc估算方法 | |
| CN116990691A (zh) | 一种电池剩余满充时间的评估方法、装置、设备和介质 | |
| CN114779093A (zh) | 基于最小误差熵扩展卡尔曼滤波的电池荷电状态估计方法 | |
| CN113671377A (zh) | 一种基于改进ekf算法的锂电池soc估算方法 | |
| CN114089193A (zh) | 电池的温度和负极电位在线估计方法、装置和计算机设备 | |
| CN110082683A (zh) | 抑制安时积分soc估计误差的闭环补偿方法 | |
| CN116794517A (zh) | 基于分数阶卡尔曼滤波的锂离子电池soc估计方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230505 |




















