基于在线宽频带阻抗的电池模型参数和SoC联合估计方法

基于在线宽频带阻抗的电池模型参数和soc联合估计方法
技术领域
1.本发明涉及储能电池技术领域，具体涉及一种基于在线宽频带阻抗的电池模型参数和soc联合估计方法。


背景技术：

2.规模化的储能电池系统所使用的单体电池容量大、电池簇单体数量多、充放电深度深且电堆电流大，在使用的过程中存在着巨大的安全隐患，因此对电池能量管系统(bms)提出了更高的要求。荷电状态(soc)是bms中最关键的参数之一，精确的soc估算可有效避免单体电芯的过充、过放；精确的电池模型及模型参数能够通过数学的形式模拟电池，便于通过算法实现电池的最优控制及管理。
3.常见的soc估计方法有：安时积分法、开路电压法和人工智能算法等。其中，安时积分法高度依赖于传感器的进度和初始soc的准确性；开路电压法是一种可以通过ocv-soc关系曲线确定soc的简单算法，但该方法需要电池经过长时间的静置来获取开路电压，不适合应用于实际现场；此外，随着技术的发展，一些人工智能算法应用与电池soc也应用于电池soc估算中，其性能取决于大量的训练数据。
4.基于模型的方法通常以一种电池模型为基础，由于电化学模型的复杂性，使得基于等效电路模型的拓展卡尔曼滤波器广泛用于soc的估计中，此类方法收敛速度快，能够对soc进行实时估计，但其估计精度受模型影响较大。双拓展卡尔曼滤波(dekf)可同时对系统模型、电池soc进行估计，提高了soc估计的精度，但该算法在初始值不确定、电池电压电流变化缓慢时，容易出现发散，造成结果不准确。
5.综上所述，基于卡尔曼滤波的soc在线估算方法存在的问题是：
6.soc估算精度不高，收敛速度不快，滤波器稳定性差；
7.电池模型参数初始值不确定、运行过程中模型参数时变，导致卡尔曼滤波算法发散或不稳定。


技术实现要素：

8.鉴于以上技术问题，本发明的目的是针对现有技术的缺陷，提供了一种基于在线宽频带阻抗的电池模型参数和soc联合估计方法，其能够同时对电池模型参数和soc进行估计，同时准确度较高，收敛速度较快。
9.本发明采用以下技术方案为：
10.一种基于在线宽频带阻抗的电池模型参数和soc联合估计方法，包括以下步骤：
11.s1、建立双极化分数阶等效电路模型，所述双极化分数阶等效电路模型包括串联的电源、电阻r0、r
sei
和r
ct
，同时还包括与r
sei
并联的恒相角元件cpe1，与r
ct
并联的恒相角元件cpe2，建立所述双极化分数阶等效电路模型的频域阻抗方程，同时建立双极化分数阶等效电路模型的双拓展卡尔曼滤波状态方程；
12.s2、获取电池的soc-ocv的关系曲线以及电化学阻抗谱；
13.s3、根据频域阻抗方程和电化学阻抗谱获取双极化分数阶等效电路模型的初始参数，并将所述参数作为双拓展卡尔曼滤波状态方程的初始值，同时设定初始值t＝1，且t为正整数；
14.s4、获取模型参数先验估计值以及模型参数先验协方差矩阵，同时获取电池状态先验估计值以及电池状态先验协方差矩阵；
15.s5、计算电池状态卡尔曼增益，同时获取电池状态后验估计值以及电池状态后验协方差矩阵；
16.s6、判断和s
×
t的大小关系，其中，s为大于0小于1的常数，s越小计算量越大，但其计算结果相对精确；s越大状态估计结果越不精确，但其计算量相对较小，计算速度相对较快：本领域技术人员可根据实际需求选定合适的值。
17.若则在向电池在线注入幅值为i
amp
、频率范围为[f
min
,f
max
]的宽频电流，其中i
amp
、f
min
、f
max
均为大于0的常量，同时获取电池的宽频电压和电流，得到电池的阻抗谱。随后采用s3的方法获得此时双极化分数阶等效电路模型的参数，并令新获得参数等于此时模型参数的后验估计值，初始化模型参数协方差矩阵；
[0018]
若|χ
est,k
(1)-χ0(1)|＜s
×
t，则计算此时模型参数的卡尔曼增益，随后计算模型参数的后验估计值，最后更新模型参数的协方差矩阵；
[0019]
s7、将s6获得的参数估计值带入步骤1进行递归运算即得。
[0020]
本发明的有益效果是：本发明针对双拓展卡尔曼滤波器在初始值不确定时，易出现不收敛问题，引入电池在线宽频带阻抗测量，将测量结果作为双拓展卡尔曼滤波器的初始值进行递归运算。在此基础上，选取一定soc间隔测量电池的宽频带阻抗，获取准确的模型参数，并将其作为估计器下一次的参数估计值。所提方法将电池在线宽频带阻抗测量方法与双拓展卡尔曼滤波器进行有机融合，一方面提高了双拓展卡尔曼滤波器的收敛速度及稳定性，另一方面提高了电池模型参数及soc估计的准确度。
附图说明
[0021]
图1为本发明的流程示意图；
[0022]
图2为本发明建立的双极化分数阶等效电路模型的结构示意图；
[0023]
图3为本发明一种基于在线宽频带阻抗的电池模型参数和soc联合估计方法的soc估算结果对比图；
[0024]
图4为本发明一种基于在线宽频带阻抗的电池模型参数和soc联合估计方法的soc估算误差对比图。
具体实施方式
[0025]
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现结合附图和实施例对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。
[0026]
如图1所示，一种基于在线宽频带阻抗的电池模型参数和soc联合估计方法，包括以下步骤：
[0027]
s1、建立双极化分数阶等效电路模型，所述双极化分数阶等效电路模型包括串联的电源、电阻r0、r
sei
和r
ct
，同时还包括与r
sei
并联的恒相角元件cpe1，与r
ct
并联的恒相角元件cpe2，建立所述双极化分数阶等效电路模型的频域阻抗方程，同时建立双极化分数阶等效电路模型的双拓展卡尔曼滤波状态方程；
[0028]
对于双极化分数阶等效电路模型而言，其包括串联的电源、电阻r0、r
sei
和r
ct
，同时还包括与r
sei
并联的恒相角元件cpe1，与r
ct
并联的恒相角元件cpe2，具体结构参见图2，图2中，u
oc
为开路电压，u
l
表示电池端电压；i表示流过整个电路的电流。其中，cpe1和cpe2的频域阻抗如下所示：
[0029][0030]
式中，α、β分别为恒相位角元件cpe1、cpe2的分数阶阶数；c1、c2分别为cpe1、cpe2的电容值；ω为角频率；
[0031]
同时根据双极化分数阶等效电路模型中各元件，建立该模型的频域阻抗方程：
[0032][0033]
式中，z(jω)为双极化等效电路模型的频域阻抗，z1(jω)为恒相角元件cpe1的频域阻抗，z2(jω)为恒相角元件cpe2的频域阻抗。
[0034]
同时建立该模型的双拓展卡尔曼滤波状态方程：
[0035][0036]ul,k
＝u
oc
(sock)-r
0,kik-u
1,k-u
2,k
+vk[0037]
式中，χk＝[sock,u
1,k
,u
2,k
]
t
；；θk＝[r
0,k
,r
sei,k
,c
1,k
,αk,r
ct,k
,c
2,k
,βk]
t
；ωk是均值为0，方差为qw的过程噪声；rk是均值为0，方差为qr的过程噪声；vk是均值为0，方差为qv的测量噪声；sock为k时刻电池的荷电状态；cq为电池的容量；u
1,k
、u
2,k
分别是k时刻电阻r
sei
、r
ct
两端电压，αk、βk分别为k时刻恒相位角元件cpe1、cpe2的分数阶阶数；ts为采样周期；k为大于等于2的正整数。
[0038]
s2、获取电池的soc-ocv的关系曲线以及电化学阻抗谱；
[0039]
soc-ocv的关系曲线采用以下方法获得：每间隔10％soc测量一次电池的开路电压，并利用6次多项式进行拟合获得，事实上，采用更多次的多项式比如7次多项式均可，但是其次数越高，计算难度越大，次数越低，则准确性越低，因此综合考虑计算准确性和计算难度，采用6次多项式进行拟合。所述6次多项式如下所示：
[0040]uoc
＝p0+p1·
soc+p2·
soc2+p3·
soc3+p4·
soc4+p5·
soc5+p6·
soc6[0041]
式中，pi(i＝0,1,
…
,6)为多项式系数，soc为电池剩余电量状态，u
oc
为电池开路电压。
[0042]
电化学阻抗谱采用以下方法获得：电池工作前，向电池注入幅值为i
amp
、频率范围
为[f
min
,f
max
]、直流偏置为0的宽频电流，其中i
amp
＞0.2a、f
min
≤10hz、f
max
≥800hz，获取电池的宽频电压和电流，得到电池的阻抗谱。此处所指的电池工作前，是指电池在投入使用前；其中，宽频电流可通过dc/dc变换器注入电池。在本实施例中，注入幅值i
amp
＝0.2a，频率范围为1hz～1khz。
[0043]
s3、根据频域阻抗方程对s2中的阻抗谱进行最小二乘拟合获取模型中的各个参数：r0的初始电阻值r
0,0
，r
sei
的初始电阻值r
sei,0
，r
ct
的初始电阻值r
ct,0
，恒相角元件cpe1的初始电容值c
1,0
，恒相角元件cpe2的初始电容值c
2,0
，恒相角元件cpe1的初始分数阶阶数α0，恒相角元件cpe2的初始分数阶阶数β0，并将这些参数作为双拓展卡尔曼滤波器的初始值，同时设定初始值t＝1，且t为正整数；
[0044]
s4、获取模型参数先验估计值以及模型参数先验协方差矩阵，同时获取电池状态先验估计值以及电池状态先验协方差矩阵；
[0045]
其中，模型参数先验估计值和模型参数先验协方差矩阵如下所示：
[0046][0047]
式中，为k时刻的模型参数先验估计值；θ
est,k-1
为k-1时刻的模型参数后验估计值；为k时刻的模型参数先验协方差矩阵；p
θ,k-1
为k-1时刻的模型参数后验协方差矩阵。
[0048]
电池状态先验估计值以及电池状态协方差矩阵如下所示：
[0049][0050]
式中，为k时刻的电池状态先验估计值；χ
est,k-1
为k-1时刻的电池状态后验估计值；为k时刻的电池状态先验协方差矩阵；p
χ,k-1
为k-1时刻的电池状态后验协方差矩阵。
[0051]
s5、计算电池状态卡尔曼增益，同时获取电池状态后验估计值以及电池状态后验协方差矩阵；
[0052]
其中，电池状态卡尔曼增益如下所示：
[0053][0054]
式中，为k时刻的状态卡尔曼增益；为k时刻的电池状态先验协方差矩阵；
[0055]
电池状态后验估计值和电池状态后验协方差矩阵如下所示：
[0056][0057]
式中，χ
est,k
为k时刻的电池状态后验估计值，为k时刻的状态先验估计值；χ
est,k-1
为k-1时刻的状态后验估计值；p
χ,k-1
为k时刻的状态后验协方差矩阵。
[0058]
s6、判断和s
×
t的大小关系，s为大于0小于1的常数，s越小计算量越大，但其计算结果相对精确；s越大状态估计结果越不精确，但其计算量相对较小，计算速度相对较快：本领域技术人员可根据实际需求选定合适的值，本实施例中s＝0.05。
[0059]
若则在电池原有工作电流基础上叠加幅值为0.2a、频率范围为1hz～1khz的宽频电流，同时获取电池新的电压、电流，随后采用s3的方法获得此时双极化分数阶等效电路模型的参数：并令此时模型参数的后验估计值等于新获得参数，即θ
est,k
＝θ
deis
，同时初始化模型参数协方差矩阵；
[0060]
若|χ
est,k
(1)-χ0(1)|＜s
×
t，则计算此时模型参数的卡尔曼增益，随后计算模型参数的后验估计值，最后更新模型参数的协方差矩阵：
[0061]
首先计算模型参数的卡尔曼增益：
[0062][0063]
式中，为k时刻的参数卡尔曼增益；为k时刻的参数卡尔曼增益；为k时刻的模型参数的先验协方差矩阵；
[0064]
随后计算模型参数的后验估计值：式中，θ
est,k
为k时刻的模型参数后验估计值。
[0065]
最后更新模型参数协方差矩阵：式中，p
θ,k
为k时刻的模型参数后验协方差矩阵
[0066]
s7、将s6获得的模型参数后验估计值θ
est,k
带入s4进行递归运算。在运算过程中，若想知道k时刻的电池状态和模型参数，调取该时刻电池状态估计结果χ
est,k
、模型参数θ
deis,k
的结果即得。
[0067]
为验证本发明中soc估计算法的准确性，将本发明提出方法与传统dekf方法进行比较，结果如图3所示，误差结果如图4所示。从图中可以看出所提方法对soc估算的精度较高，并且收敛速度较快。
[0068]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明实施例揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。