基于非线性观测器的超级电容器荷电状态在线估计

doi:10.16183/j.cnki.jsjtu.2021.210

基于非线性观测器的超级电容器荷电状态在线估计

杜宇石¹^,², 琚长江^,¹^,², 杨根科¹^,²

1.上海交通大学宁波人工智能研究院, 浙江宁波 315000

2.上海交通大学自动化系, 上海 200240

Online Estimation of Supercapacitor State of Charge Based on Nonlinear Observer

DU Yushi¹^,², JU Changjiang^,¹^,², YANG Genke¹^,²

1. Ningbo Artificial Intelligence Institute, Shanghai Jiao Tong University, Ningbo 315000, Zhejiang,China

2. Department of Automation, Shanghai Jiao Tong University, Shanghai 200240, China

通讯作者: 琚长江, 男, 副研究员; E-mail:juchangjiang@sjtu.edu.cn.

责任编辑: 王一凡

收稿日期: 2021-06-15

基金资助:

国家重点研发计划(2020YFB1711200)

Received: 2021-06-15

作者简介 About authors

杜宇石(1997-),男,甘肃省庆阳市人,硕士生,从事电池能量控制管理系统研究.

摘要

超级电容器具有快速充放、高功率密度和长寿命等优点, 被广泛用于新能源汽车的储能系统. 系统可靠运行需获取其剩余电量, 即对其荷电状态(SOC)进行估算. 依托超级电容单体的等效模拟电路模型, 建立了以模型中多电容端电压为状态, 电容器输入电流为控制输入, 电容器输出电压为观测输出的电容器二阶非线性系统的状态空间模型, 包含了自放电现象产生的泄漏电流的因素. 为提高模拟精度, 辨识不同的模型参数,分别刻画充电和放电工况. 采用非线性观测器算法来获取模型内部状态从而实现对SOC的估计. 充放电实验的结果表明, 考虑泄漏因素和建立不同参数下的充放电模型, 能够更好地模拟超级电容器的动态特性, 同时验证了非线性观测器算法具有稳定的跟踪能力.

关键词： 超级电容器; 荷电状态估计; 非线性观测器; 泄漏电流

Abstract

Supercapacitors have the advantages of fast charging and discharging, high power density, and long life, which are widely used in energy storage systems for new energy vehicles. Reliable operation of the system requires the acquisition of its remaining electricity, which is to estimate its state of charge (SOC). Relying on the equivalent analog circuit model of a single supercapacitor, this paper establishes a state-space model of the capacitor second-order nonlinear system with the multi-capacitor terminal voltage in the model as the state, the capacitor input current as the control input, and the capacitor output voltage as the observation output, and contains the leakage current caused by the self-discharge phenomenon. In order to improve the simulation accuracy, different model parameters were identified to characterize the charging and discharging conditions. In this paper, a nonlinear observer algorithm is used to obtain the internal state of the model to realize the estimation of SOC. The results of the charging and discharging experiment show that considering the leakage factor and establishing the charging and discharging model with different parameters, the dynamic characteristics of the supercapacitor can be better simulated, and the nonlinear observer algorithm has a stable tracking ability.

Keywords： supercapacitor; state of charge (SOC) estimation; nonlinear observer; leakage current

PDF (1181KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

杜宇石, 琚长江, 杨根科. 基于非线性观测器的超级电容器荷电状态在线估计[J]. 上海交通大学学报, 2022, 56(12): 1630-1637 doi:10.16183/j.cnki.jsjtu.2021.210

DU Yushi, JU Changjiang, YANG Genke. Online Estimation of Supercapacitor State of Charge Based on Nonlinear Observer[J]. Journal of Shanghai Jiaotong University, 2022, 56(12): 1630-1637 doi:10.16183/j.cnki.jsjtu.2021.210

超级电容器作为一种新型的储能元件, 具有充放电迅速、功率密度高、循环寿命长以及能量转换效率高等优点, 在混合动力汽车等领域有越来越广泛的应用. 为了汽车安全可靠运行, 判断可行驶路程等信息, 需要对其剩余电量, 即荷电状态(State of Charge, SOC)进行估算. SOC估计方法分安时法和电压软测量两类. 安时法^[1]是一种基于直接电流测试的SOC估计方法, 具有较高的精度, 是实验室中最常用的SOC估计方法. 但从长时运维来看, 由于自放电现象和电流传感器精度的影响, 如果不定期进行校正, 误差会越来越大; 基于模型的电压软测量方法是机理或半机理或无机理的, 具有精度高和鲁棒性强等优点, 被广泛用于SOC 的实时估计^[2]. 本文主要研究基于模型的SOC估算方法.

超级电容器的软测量模型基本可以分为4种: 电化学模型、等效电路模型、智能模型和分数阶模型. 电化学模型能够描述超级电容器内部真实的反应过程, 具有精度高的优点, 但是耦合偏微分方程带来了很大的计算量, 某些参数也是不可测量的^[3]. 等效电路模型是根据经验和实验数据建立的, 具有良好的准确性和实用性, 但也存在物理意义缺失和易受老化影响等不足^[4]. 智能模型能够描述自变量和因变量之间的复杂非线性关系, 而不需要详细了解其潜在机制和原理, 具有很好的建模能力, 但是对训练数据的质量和数量要求高, 非机理模型鲁棒性也不强^[5]. 分数阶模型引入了分数阶微积分, 有较好的拟合实验数据的能力, 但计算成本大, 嵌入式实施难度高^[6]. 等效电路模型结构简单且建模精度高, 在实时能量管理系统中得到了广泛的应用. 因此, 本文深入研究基于等效电路模型的SOC估计算法.

Pavkovic等^[7]基于一阶非线性模型提出了基于扩展卡尔曼滤波的SOC估计算法. 此外, 还有如无迹卡尔曼滤波^[8]等一系列改进方法. 基于卡尔曼滤波的SOC估计器虽然具有较高的精度, 但是此类方法需要精确的模型以及噪声的统计信息, 并且计算成本较高. Ceraolo等^[9]基于线性RC串联电路模型采用Luenberger型观测器来计算SOC, 但单一的线性RC 串联电路不能完全模拟超级电容器的充放电特性. 以上算法虽然原理不同, 但在估计SOC时, 皆用等效电路模型中的非线性电容存储的电荷近似代替总电荷, 忽略了模型中其余的定值电容存储的电荷, 这在一定程度上会降低算法的精度.

为了得到更加准确的SOC估计值, Fan等^[10]将SOC的定义扩充到所有电容存储的电荷, 利用滑模观测器来估计SOC, 在降低计算量的同时取得了较高的精度. 但是文中讨论的模型未考虑由于自放电现象产生的泄漏电流的影响, 这在短期内忽略是合理的, 而在长期使用或静置时可能会带来一定的累积误差. 另外, 超级电容器在充电和放电工况下的特性并不一致, 在不同的工况下建立不同的模型参数来表征其充放电特性可提高模拟精度^[11].

本文在二分支等效模拟电路模型^[12]的基础上, 建立内部多电容端电压的状态空间模型, 将超级电容器的电流和电压分别作为可控制输入和可观测输出, 来计算内部多电容的端电压值; 再根据非线性观测器^[13]的思想, 用内部多电容的端电压来重构SOC. 针对充电和放电两种工况分别辨识模型参数, 每个模型包括有泄漏和无泄漏两种情况. 最后, 通过充放电实验讨论区分充放电不同工况以及泄漏电流对SOC估计精度的影响, 并验证非线性观测器的性能.

1 超级电容器两分支等效电路建模

超级电容器有一对金属箔电极, 每个电极表面都有活性炭, 多孔的活性炭表面在电解液中与离子形成电双层, 可以决定超级电容器的容性. 超级电容器的物理化学特性比较复杂, 单一的电路器件不足以描述其充放电行为. 因此, 研究人员通常利用一些简化的等效电路模型来代替超级电容器^[14].

本文综合考虑了超级电容单体等效电路模型的准确性和实用性, 选用Faranda^[12]提出的二分支等效模拟电路模型, 如图1所示. 图中: $R_{0}$ 为等效串联电阻; $R_{2}$ 为电荷再分配电阻; $C_{1}$ 为微分电容; $C_{2}$ 为小孔隙电容(电荷再分配电容); $R_{l}$ 为泄漏电阻; $i_{1}$ 为短期支路电流; $i_{2}$ 为中长期支路电流; $i_{3}$ 为泄漏电流;i为输入的充放电电流; $v_{1}$ 为微分电容 $C_{1}$ 两端的电压; $v_{2}$ 为小孔隙电容 $C_{2}$ 两端的电压;v为单体的端电压.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 二分支等效模拟电路模型

Fig.1 Two-branch equivalent analog circuit model

利用Kirchhoff定律, 建立了基于图1所示的二阶非线性系统的状态空间模型:

(1)

\begin{matrix} \begin{array}{l} {\overset{\cdot}{x}}_{1} = ϕ (x, u) = \frac{a_{11} x_{1}}{1 + a x_{1}} + \\ \frac{a_{12} x_{2}}{1 + a x_{1}} + \frac{b_{1}}{1 + a x_{1}} u \\ {\overset{\cdot}{x}}_{2} = a_{21} x_{1} + a_{22} x_{2} + b_{2} u \\ y = d_{1} x_{1} + d_{2} x_{2} + d u \end{array}\} \end{matrix}

式中:状态变量 $x_{1} = v_{1} 和 x_{2} = v_{2} 为电容器内部多电容的端电压; 电容器的可控输入 u = i 为充放电电流, 可测输出 y = v$ 为端电压.

(2)

\begin{matrix} \begin{array}{l} a_{11} = \frac{- R_{2} - R_{l}}{R_{e} C_{0}}, a_{12} = \frac{R_{l}}{R_{e} C_{0}} \\ a_{21} = \frac{R_{l}}{C_{2} R_{e}}, a_{22} = \frac{- R_{0} - R_{l}}{C_{2} R_{e}} \\ a = \frac{2 k}{C_{0}}, b_{1} = \frac{R_{2} R_{l}}{R_{e} C_{0}} \\ b_{2} = \frac{R_{0} R_{l}}{C_{2} R_{e}}, d_{1} = \frac{R_{2} R_{l}}{R_{e}} \\ d_{2} = \frac{R_{0} R_{l}}{R_{e}}, d = \frac{R_{0} R_{2} R_{l}}{R_{e}} \\ R_{e} \overset{Δ}{=} R_{0} R_{2} + R_{2} R_{l} + R_{0} R_{l} \end{array}\} \end{matrix}

式中: $k 为微分电容 C_{1} 随两端电压 v_{1} 线性变化的斜率; C_{0} 为微分电容 C_{1}$ 两端电压为0时的初始电容. 这些参数因超级电容器的型号及规格而不同.

流过泄漏电阻 $R_{l}$ 的电流是自放电现象产生的泄漏电流,一般都比较小.由于泄漏电阻 $R_{l}$ 远大于 $R_{0} 和 R_{2},$ 所以泄漏电流比较小, 在某些情况下也可以忽略. 电容器的状态空间模型式(1)在无泄漏电阻情况下^[15]的各参数为

(3)

\begin{matrix} \begin{array}{l} a_{11} = \frac{- 1}{R_{e} C_{0}}, a_{12} = \frac{1}{R_{e} C_{0}} \\ a_{21} = \frac{1}{C_{2} R_{e}}, a_{22} = \frac{- 1}{C_{2} R_{e}}, a = \frac{k}{C_{0}} \\ b_{1} = \frac{R_{2}}{R_{e} C_{0}}, b_{2} = \frac{R_{0}}{C_{2} R_{e}} \\ d_{1} = \frac{R_{2}}{R_{e}}, d_{2} = \frac{R_{0}}{R_{e}} \\ d = \frac{R_{0} R_{2}}{R_{e}}, R_{e} \overset{Δ}{=} R_{0} + R_{2} \end{array}\} \end{matrix}

注意到无泄漏效应对应模型式(3)并非式(2)的特例.

文献[11]针对模型参数选取不准给扩展卡尔曼滤波SOC估计带来误差的问题, 设计了变电流放电工况和变电流充电工况实验, SOC估计结果显示: 对于变电流充电工况, 使用充电模型参数的扩展卡尔曼滤波SOC 估计值更接近真实值, 而使用放电模型参数的SOC估计值则偏离真实值; 反之亦然. 因此, 电池充电模型参数和放电模型参数是不同的, 如果单独使用一种模型的参数进行SOC估计就会增大误差.

本文针对充电和放电两种不同的工况, 分别辨识不同的模型参数来表征超级电容器不同的物理特性.

2 充放电实验与模型参数辨识

以某型号为UCK42V28000的超级电容器单体为测试研究对象. 对超级电容单体进行100 A恒定电流的充放电实验, 规定充电电流为正, 放电电流为负, 直到超级电容端电压达到额定电压3.8 V为止, 分别记录电流i和端电压v的测试数据, 如图2和图3所示. 图中: t为时间.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 充电工况的输入电流和端电压

Fig.2 Input current and terminal voltage of charging conditions

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 放电工况的输入电流和端电压

Fig.3 Input current and terminal voltage of discharging conditions

参考Liu等^[16]使用遗传算法辨识等效电路模型相关参数的方法, 分别对充电和放电两种不同场景下的模型参数进行辨识. 对比不区分充放电工况^[8], 即仅一组参数来表征超级电容器的充放电特性, 辨识了单一模型的参数. 各模型参数如表1所示.

表1 辨识的模型参数

Tab.1 Model parameters identified

模型类别	R₀/Ω	R₂/Ω	R_l/Ω	C₀/F	k/ (F·V^-1)	C₂/F
单一模型	0.000 63	0.5	3 000	28 000	350	500
充电模型	0.001 5	0.000 05	3 000	28 000	440	2 000
放电模型	0.000 53	0.8	3 000	28 000	400	0.1

新窗口打开| 下载CSV

在MATLAB R2020a/Simulink环境中, 保持和实验一致的条件, 对式(1)进行仿真. 设置带噪声的100 A电流作为输入, 仿真时间与试验数据一致, 即272.6 s. 分别记录有泄漏和无泄漏情况下输出的端电压值, 并与实验得到的端电压值作对比, 如图4所示.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 充电时端电压的仿真值和实测值比较

Fig.4 Comparison of simulated value and measured value of terminal voltage during charging

不引入泄漏电阻情况下的端电压仿真值与实测值间的误差会越来越大, 整体最大误差为 0.075 V. 而引入泄漏电阻情况下的端电压仿真值与实测值在中间部分的误差比较大, 但在两端的误差都比较小, 整体最大误差为0.048 V. 在充电起始阶段和接近结束时误差变小, 而在充电中间过程误差较大, 这是由于超级电容器端电压的非线性变化以及模型中RC电路的滞后特性引起的. 在起始阶段和结束阶段, 模型能够比较准确地反映端电压的变化情况;而在中间过程, 模型准确性相对较低. 在实际工程应用场景中, 超级电容器模型在接近充电结束时的准确度更加重要, 这对精确估计SOC值从而判断超级电容器是否充满电具有非常重要的意义. 在充电场景下, 考虑泄漏电流因素, 即引入泄漏电阻情况下的模型能够提高约0.03 V的精度.

类似地, 设置带噪声的 -100 A 电流作为输入(规定充电电流为正,放电电流为负), 测试时间为272.4 s. 分别记录有泄漏和无泄漏情况下输出的端电压值, 并与实验得到的端电压值作对比, 如图5所示.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 放电时端电压的仿真值和实测值比较

Fig.5 Comparison of simulated value and measured value of terminal voltage during discharging

在放电初始阶段, 无泄漏电阻和有泄漏电阻两种情况的端电压基本一致. 随着放电过程的进行, 无泄漏电阻情况下的端电压误差先增大再减小, 最大误差为0.06 V, 在接近结束的时候误差会变小, 但与实测值仍然有0.038 V的误差. 有泄漏电阻情况下的端电压仿真值与实测值的误差较小, 仅在放电中间阶段误差比较大, 最大误差为0.036 V, 在接近结束时误差会变得很小, 最终为0.004 V. 类似地, 端电压在放电中间过程误差较大, 在起始和结束时的误差较小, 这也是由于超级电容器端电压的非线性变化和模型中RC电路的滞后特性导致的. 与无泄漏电阻的模型相比, 引入泄漏电阻能够把端电压最大误差从1.60%降低到0.96%, 将接近放电结束时的端电压误差从1.02%降低到0.11%. 另外, 超级电容器内部存在自放电现象, 故需要在建模时引入泄漏电阻来描述这一特性. 因此, 在放电场景下, 考虑泄漏电流, 即引入泄漏电阻情况下的模型有更好的准确性.

综上所述, 在充电和放电不同场景下, 考虑泄漏效应的模型有更高的模拟精度.

3 非线性观测器及性能验证

引入文献[10]的SOC, 定义为

(4)

\begin{matrix} S O C = \frac{1}{Q_{n}} [(C_{0} + k v_{1}) v_{1} + C_{2} v_{2}] \end{matrix}

式中: $Q_{n} 为超级电容器的额定容量, 即存储的最大电荷量 . 本文的 Q_{n}$ 为 25 200 C. 系统模型式(1)是非线性的, 因此利用非线性观测器来观测内部电压, 进而估计SOC.

对式(1)构造非线性状态观测器:

(5)

\begin{matrix} \begin{array}{l} {\overset{\cdot}{\hat{x}}}_{1} = ϕ (\hat{x}, u) + λ_{1} (y - \hat{y}) \\ {\overset{\cdot}{\hat{x}}}_{2} = a_{21} {\hat{x}}_{1} + a_{22} {\hat{x}}_{2} + b_{2} u + λ_{2} (y - \hat{y}) \\ \hat{y} = d_{1} {\hat{x}}_{1} + d_{2} {\hat{x}}_{2} + d u \end{array}\} \end{matrix}

式中: $λ_{1}, λ_{2}$ 为观测器的增益 $; {\overset{^{}}{x}}_{1}$ 为模型内部电容C₁ 端电压v₁ 的估计值; ${\overset{^{}}{x}}_{2}$ 为模型内部电容C₂ 端电压v₂ 的估计值; $\overset{^{}}{y}$ 为电容器终端电压v的估计值; ${\overset{\cdot}{\overset{^{}}{x}}}_{1} 和 {\overset{\cdot}{\overset{^{}}{x}}}_{2}$ 分别为 ${\hat{x}}_{1} 和 {\hat{x}}_{2}$ 的一阶导数.

非线性观测器式(5)的全局渐近稳定性的证明过程如下.

设控制输入信号是有界的, 即存在 $u_{M} > 0,$ 对于任意的t,有 $|u (t)| \leq u_{M} .$

定义饱和函数:

(6)

\begin{matrix} σ (x) = \{\begin{array}{l} x_{M}, & x \geq x_{M} \\ x, & 0 < x < x_{M} \\ 0, & x \leq 0 \end{array} \end{matrix}

式中:x_M是充分大的正常数.

由于状态变量电压x_k非负, 将其代入式(6)可得:

\begin{matrix} σ (x_{k}) = x_{k}, k = 1, 2 \end{matrix}

且所得函数σ(x_k)具有Lipschiz性质, 即对于任意的 $n 维向量 x, y \in R^{n},$ 有

|σ (x) - σ (y)| \leq |x - y|

定义函数:

\begin{array}{l} ϕ_{σ} (x, u) = ϕ_{1} (x_{1}) ϕ_{2} (x, u) \\ ϕ_{1} (x_{1}) = \frac{1}{1 + a σ (x_{1})} \\ ϕ_{2} (x, u) = a_{11} σ (x_{1}) + a_{12} σ (x_{2}) + b_{1} u \end{array}

注意到函数 $ϕ_{σ} (x, u) 对 x$ 满足Lipschiz性质, 即存在实数γ>0, 有

|ϕ_{σ} (x, u) - ϕ_{σ} (\hat{x}, u)| \leq γ |x - \hat{x}|

其中: $\hat{x} 为向量 x$ 的估计值.

对于式(5), 估计误差 $\tilde{x} = x - \hat{x} .$ 根据式(1)和(5)可得:

(7)

\begin{array}{l} {\overset{\cdot}{\tilde{x}}}_{1} = ϕ (x, u) - ϕ_{σ} (\tilde{x}, u) - λ_{1} (d_{1} {\tilde{x}}_{1} + d_{2} {\tilde{x}}_{2}) \\ {\overset{\cdot}{\tilde{x}}}_{2} = (a_{21} - d_{1} λ_{2}) {\tilde{x}}_{1} + (a_{22} - d_{2} λ_{2}) {\hat{x}}_{2} \end{array}\}

对任意的t≥0, 定义Lyapunov函数:

\begin{matrix} V (t) = \frac{1}{2} ({\tilde{x}}_{1}^{2} (t) + {\tilde{x}}_{2}^{2} (t)) \end{matrix}

此函数关于t求导, 得:

(8)

\begin{matrix} \overset{\cdot}{V} (t) = \frac{d V}{d t} = {\tilde{x}}_{1} {\overset{\cdot}{\tilde{x}}}_{1} + {\tilde{x}}_{2} {\overset{\cdot}{\tilde{x}}}_{2} \end{matrix}

将(7)式代入(8)式得:

(9)

\begin{array}{l} \overset{\cdot}{V} = {\tilde{x}}_{1} \tilde{ϕ} - λ_{1} d_{1} {\tilde{x}}_{1}^{2} + (a_{22} - λ_{2} d_{2}) {\tilde{x}}_{2}^{2} + \\ (a_{21} - λ_{1} d_{2} - λ_{2} d_{1}) {\tilde{x}}_{1} {\tilde{x}}_{2} \end{array}

式中: $\tilde{ϕ} = ϕ_{σ} (x, u) - ϕ_{σ} (\hat{x}, u) \leq γ |x - \hat{x}| .$

由于ϕ_σ具有Lipschitz性, 所以存在γ≥0, 有

(10)

\begin{array}{l} {\tilde{x}}_{1} \tilde{ϕ} \leq |{\tilde{x}}_{1}| |\tilde{ϕ}| = |{\tilde{x}}_{1}| |ϕ_{σ} (x, u) - ϕ_{σ} (\hat{x}, u)| \leq \\ γ |{\tilde{x}}_{1}| |x - \hat{x}| = γ |{\tilde{x}}_{1}| |\tilde{x}| = \\ γ ({\tilde{x}}_{1}^{2} + |{\tilde{x}}_{1}| |{\tilde{x}}_{2}|) \end{array}

根据杨氏不等式, 存在设计常数ε>0, 使得

(11)

\begin{matrix} |{\tilde{x}}_{1}| |{\tilde{x}}_{2}| \leq \frac{{\tilde{x}}_{1}^{2}}{2 ε} + \frac{ε {\tilde{x}}_{2}^{2}}{2} \end{matrix}

将式(11)代入式(10)得:

(12)

\begin{array}{l} {\tilde{x}}_{1} \tilde{ϕ} \leq γ ({\tilde{x}}_{1}^{2} + |{\tilde{x}}_{1}| |{\tilde{x}}_{2}|) \leq \\ (γ + \frac{γ}{2 ε}) {\tilde{x}}_{1}^{2} + \frac{γ ε}{2} {\tilde{x}}_{2}^{2} \end{array}

式(12)代入式(9), 并整理可得:

\begin{array}{l} \overset{\cdot}{V} \leq - β_{1} {\tilde{x}}_{1}^{2} - β_{2} {\tilde{x}}_{2}^{2} \\ β_{1} = - γ + λ_{1} d_{1} - \frac{γ + a_{21} - λ_{1} d_{2} - λ_{2} d_{1}}{2 ε} \\ β_{2} = - a_{22} + λ_{2} d_{2} - \frac{ε (γ + a_{21} - λ_{1} d_{2} - λ_{2} d_{1})}{2 ε} \end{array}

选择合适的观测器增益 $λ_{1}, λ_{2} 和设计常数 ε, 使得 β_{1} > 0, β_{2} > 0, 则 \overset{\cdot}{V} < 0 .$

根据Lyapunov定理, 平衡点 $({\tilde{x}}_{1}, {\tilde{x}}_{2}) = (0, 0)$ 是全局渐进稳定的.

因此, 可观测到系统内部状态 $x_{1} = v_{1}, x_{2} = v_{2}$ 的信息, 代入式(4)从而求得SOC的估计值.

通过实验数据来验证式(5)的有效性. 参考文献[13], 考虑到非线性观测器的鲁棒性和收敛速度, 选择观测器增益为 $λ_{1} = 7, λ_{2} = 9,$ 充电电流和放电电流均与第2节中的保持一致. 在恒流的模式下进行充放电, 采用安时法得到的SOC作为实际参照值.

在充电工况下, 为了研究区分充放电不同工况对SOC估计的影响, 把基于充电模型的SOC估计值与参照值之间的误差和基于充放电单一模型^[8]的SOC估计值与参照值之间的误差作对比, 如图6所示. 图中:SOCE为SOC估计误差. 计算充电模型和单一模型的最大误差(ME)、平均绝对误差(MAE)和均方根误差(RMSE), 如表2所示.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 充电模型与单一模型的SOCE比较

Fig.6 Comparison of SOCE between charging model and single model

表2 单一模型和充电模型的SOC估计误差

Tab.2 SOC estimation error of single model and charging model%

模型类别	ME	MAE	RMSE
单一模型	2.493	0.748	1.065
充电模型	0.184	0.107	0.122

新窗口打开| 下载CSV

由图6和表2可知, 不区分充放电的模型得到的SOC估计误差是逐渐增大的, 并不收敛.本文建立的充电模型得到的估计误差波动范围很小. 因此在充电时, 基于充电模型估计的SOC值比基于不区分充放电模型估计的SOC值更加准确.

把仿真得到的观测器估计的SOC值(包括基于有泄漏的模型和基于无泄漏的模型^[10]两种情况) 与参照值进行比较, 如图7所示. 计算有泄漏模型和无泄漏模型的ME、MAE和RMSE, 结果如表3所示.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 充电时有无泄漏估计的SOC值与参照值比较

Fig.7 Comparison of estimated SOC value with reference value when charging with or without leakage

表3 充电工况下不同模型的SOC估计误差

Tab.3 SOC estimation errors of different models under charging conditions%

模型类别	ME	MAE	RMSE
不考虑泄漏	7.776	3.653	4.290
考虑泄漏	0.184	0.108	0.123

新窗口打开| 下载CSV

无泄漏电流情况下非线性观测器估计的SOC值与实际参照值的误差会越来越大, 到充电后期, 甚至会超过1, 这是很不合理的. 而在考虑有泄漏电阻情况下估计的SOC值与实际参照值吻合程度非常好, 均方根误差降低了4%. 因此, 考虑泄漏会提高SOC的估计精度.

在放电工况下, 将基于原充放电单一模型和放电模型得到的SOCE进行对比, 如图8所示. 计算其ME、MAE和RMSE, 如表4所示.

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 放电模型与单一模型的SOCE比较

Fig.8 SOCE comparison of discharging model and single model

表4 单一模型和放电模型的SOC估计误差

Tab.4 SOC estimation error of single model and discharging model%

模型类别	ME	MAE	RMSE
单一模型	2.087	0.731	0.962
放电模型	0.303	0.190	0.214

新窗口打开| 下载CSV

由图8和表4可知, 基于单一模型估计的SOC范围并没有包含0%~2%之间的部分, 这会损失一部分的SOC信息. SOC估计误差是逐渐增大的, 并不收敛. 而基于放电模型估计的SOC在合理的范围内, 均方根误差比基于单一模型的小4倍左右. 因此, 基于放电模型的估计效果更好.

将基于无泄漏和有泄漏的模型仿真得到的SOC估计值进行对比, 如图9所示. 计算其ME、MAE和RMSE, 如表5所示.

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 放电时有无泄漏估计的SOC值与参照值比较

Fig.9 When discharging with or without leakage, the estimated SOC value is compared with the reference value

表5 放电工况下不同模型的SOC估计误差

Tab.5 SOC estimation errors of different models under discharging conditions%

模型类别	ME	MAE	RMSE
不考虑泄漏	3.744	2.097	2.365
考虑泄漏	0.303	0.195	0.217

新窗口打开| 下载CSV

无泄漏情况下非线性观测器估计的SOC值甚至会低于0, 这在实际应用中不可能发生, 且SOC估计误差逐渐增大. 而有泄漏情况下估计的SOC最大误差和均方根误差都比无泄漏情况下的小一个数量级. 所以, 考虑泄漏会提高SOC估计结果的可靠性.

综上所述, 非线性观测器算法具有良好的估计能力, 考虑泄漏因素会提高模型的准确度, 而区分充放电不同工况进行建模可以更加准确地表征超级电容器的充放电特性, 从而显著降低SOC的估计误差.

4 结语

以超级电容器单体的二分支等效模拟电路模型为基础, 推导了以电容器模型中多电容端电压为状态的二阶非线性系统的状态空间方程. 针对充电和放电工况不同, 辨识了充电和放电各自的模型参数集. 讨论了自放电现象产生的泄漏电流对模型精度的影响. 实际的充放电实验数据表明, 与不考虑泄漏相比, 考虑泄漏能够降低约1%的电容器端电压模型预测误差. 依托电容器模型中多电容端电压, 设计了非线性观测器来估计SOC的值, 并通过仿真和实验对观测器进行了验证. 结果表明, 非线性观测器算法具有稳定的跟踪能力. 与不区分充电放电相比, 区分充电放电能够降低SOC的估计误差, 最大误差降低了一个数量级, 故基于区分充电放电的模型能够更好模拟实际充放电特性, 显著提高了SOC的估计精度. 不考虑泄漏得到的SOC估计值在充放电后期会偏离实际参照值, 甚至不在合理范围内, 所以考虑泄漏会降低SOC估计值的误差, 提高结果的可靠性. 本文提出的考虑泄漏因素且区分充电放电不同工况估计SOC的方法具有一定的理论价值和实用性, 后续将进一步在车辆加速、匀速、减速等不同工况下开展深入研究和应用.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

WANG

, CAO

, CHEN

, et al.

Combined state of charge estimator for electric vehicle battery pack

[J]. Control Engineering Practice, 2007, 15(12): 1569-1576.

DOI:10.1016/j.conengprac.2007.03.004 URL [本文引用: 1]

[2]

ZHANG

, HU

, WANG

, et al.

A review of supercapacitor modeling, estimation, and applications: A control/management perspective

[J]. Renewable and Sustainable Energy Reviews, 2018, 81: 1868-1878.

DOI:10.1016/j.rser.2017.05.283 URL [本文引用: 1]

[3]

ALLU

, ASOKAN

B V

, SHELTON

W A

, et al.

A generalized multi-dimensional mathematical model for charging and discharging processes in a supercapacitor

[J]. Journal of Power Sources, 2014, 256: 369-382.

DOI:10.1016/j.jpowsour.2014.01.054 URL [本文引用: 1]

[4]

PARVINI

, SIEGEL

J B

, STEFANOPOULOU

A G

, et al.

Supercapacitor electrical and thermal modeling, identification, and validation for a wide range of temperature and power applications

[J]. IEEE Transactions on Industrial Electronics, 2015, 63(3): 1574-1585.

DOI:10.1109/TIE.2015.2494868 URL [本文引用: 1]

[5]

, LI

S E

, YANG

Advanced machine learning approach for lithium-ion battery state estimation in electric vehicles

[J]. IEEE Transactions on Transportation electrification, 2015, 2(2): 140-149.

DOI:10.1109/TTE.2015.2512237 URL [本文引用: 1]

[6]

WANG

, LI

S E

, PENG

, et al.

Fractional-order modeling and parameter identification for lithium-ion batteries

[J]. Journal of Power Sources, 2015, 293: 151-161.

DOI:10.1016/j.jpowsour.2015.05.059 URL [本文引用: 1]

[7]

PAVKOVIC

, KOMLJENOVI

, HRGETI

, et al.

Design of EKF-based SOC estimator for an ultracapacitor module

[C]//IEEE EUROCON 2015—International Conference on Computer as a Tool. Salamanca, Spain:IEEE, 2015: 1-6.

[本文引用: 1]

[8]

SAHA

, DEY

, KHANRA

Accurate estimation of state-of-charge of supercapacitor under uncertain leakage and open circuit voltage map

[J]. Journal of Power Sources, 2019, 434(9): 226696.

DOI:10.1016/j.jpowsour.2019.226696 URL [本文引用: 3]

[9]

CERAOLO

, LUTZEMBERGER

, POLI

State-of-charge evaluation of supercapacitors

[J]. Journal of Energy Storage, 2017, 11: 211-218.

DOI:10.1016/j.est.2017.03.001 URL [本文引用: 1]

[10]

FAN

, DUAN

, SUN

, et al.

State of charge estimate for super-capacitor based on sliding mode observer

[C]//2017 IEEE Transportation Electrification Conference and Expo, Asia-Pacific. Harbin, China: IEEE, 2017: 1-5.

[本文引用: 3]

[11]

武国良.

电动汽车用镍氢电池剩余电量估计方法研究[D]. 哈尔滨: 哈尔滨工业大学, 2010.

[本文引用: 2]

Guoliang

Research on state of charge estimation of Ni-MH battery used in electric vehicles[D]. Harbin:Harbin Institute of Technology, 2010.

[本文引用: 2]

[12]

FARANDA

A new parameters identification procedure for simplified double layer capacitor two-branch model

[J]. Electric Power Systems Research, 2010, 80(4): 363-371.

DOI:10.1016/j.epsr.2009.10.024 URL [本文引用: 2]

[13]

XIA

, CHEN

, TIAN

, et al.

A novel method for state of charge estimation of lithium-ion batteries using a nonlinear observer

[J]. Journal of Power Sources, 2014, 270: 359-366.

DOI:10.1016/j.jpowsour.2014.07.103 URL [本文引用: 1]

[14]

ZUBIETA

, BONERT

Characterization of double-layer capacitors for power electronics applications

[J]. IEEE Transactions on Industry Applications, 2000, 36(1): 199-205.

DOI:10.1109/28.821816 URL [本文引用: 1]

[15]

WANG

, HE

, ZHANG

, et al.

A comparative study on the applicability of ultracapacitor models for electric vehicles under different temperatures

[J]. Applied Energy, 2017, 196: 268-278.

DOI:10.1016/j.apenergy.2017.03.060 URL [本文引用: 1]

[16]

LIU

, WANG

, CHEN

, et al.

A variable capacitance based modeling and power capability predicting method for ultracapacitor

[J]. Journal of Power Sources, 2018, 374: 121-133.

DOI:10.1016/j.jpowsour.2017.11.033 URL [本文引用: 1]

Combined state of charge estimator for electric vehicle battery pack

2007

... 超级电容器作为一种新型的储能元件, 具有充放电迅速、功率密度高、循环寿命长以及能量转换效率高等优点, 在混合动力汽车等领域有越来越广泛的应用. 为了汽车安全可靠运行, 判断可行驶路程等信息, 需要对其剩余电量, 即荷电状态(State of Charge, SOC)进行估算. SOC估计方法分安时法和电压软测量两类. 安时法^[1]是一种基于直接电流测试的SOC估计方法, 具有较高的精度, 是实验室中最常用的SOC估计方法. 但从长时运维来看, 由于自放电现象和电流传感器精度的影响, 如果不定期进行校正, 误差会越来越大; 基于模型的电压软测量方法是机理或半机理或无机理的, 具有精度高和鲁棒性强等优点, 被广泛用于SOC 的实时估计^[2]. 本文主要研究基于模型的SOC估算方法. ...

A review of supercapacitor modeling, estimation, and applications: A control/management perspective

2018

A generalized multi-dimensional mathematical model for charging and discharging processes in a supercapacitor

2014

... 超级电容器的软测量模型基本可以分为4种: 电化学模型、等效电路模型、智能模型和分数阶模型. 电化学模型能够描述超级电容器内部真实的反应过程, 具有精度高的优点, 但是耦合偏微分方程带来了很大的计算量, 某些参数也是不可测量的^[3]. 等效电路模型是根据经验和实验数据建立的, 具有良好的准确性和实用性, 但也存在物理意义缺失和易受老化影响等不足^[4]. 智能模型能够描述自变量和因变量之间的复杂非线性关系, 而不需要详细了解其潜在机制和原理, 具有很好的建模能力, 但是对训练数据的质量和数量要求高, 非机理模型鲁棒性也不强^[5]. 分数阶模型引入了分数阶微积分, 有较好的拟合实验数据的能力, 但计算成本大, 嵌入式实施难度高^[6]. 等效电路模型结构简单且建模精度高, 在实时能量管理系统中得到了广泛的应用. 因此, 本文深入研究基于等效电路模型的SOC估计算法. ...

Supercapacitor electrical and thermal modeling, identification, and validation for a wide range of temperature and power applications

2015

Advanced machine learning approach for lithium-ion battery state estimation in electric vehicles

2015

Fractional-order modeling and parameter identification for lithium-ion batteries

2015

Design of EKF-based SOC estimator for an ultracapacitor module

2015

... Pavkovic等^[7]基于一阶非线性模型提出了基于扩展卡尔曼滤波的SOC估计算法. 此外, 还有如无迹卡尔曼滤波^[8]等一系列改进方法. 基于卡尔曼滤波的SOC估计器虽然具有较高的精度, 但是此类方法需要精确的模型以及噪声的统计信息, 并且计算成本较高. Ceraolo等^[9]基于线性RC串联电路模型采用Luenberger型观测器来计算SOC, 但单一的线性RC 串联电路不能完全模拟超级电容器的充放电特性. 以上算法虽然原理不同, 但在估计SOC时, 皆用等效电路模型中的非线性电容存储的电荷近似代替总电荷, 忽略了模型中其余的定值电容存储的电荷, 这在一定程度上会降低算法的精度. ...

Accurate estimation of state-of-charge of supercapacitor under uncertain leakage and open circuit voltage map

2019

... 参考Liu等^[16]使用遗传算法辨识等效电路模型相关参数的方法, 分别对充电和放电两种不同场景下的模型参数进行辨识. 对比不区分充放电工况^[8], 即仅一组参数来表征超级电容器的充放电特性, 辨识了单一模型的参数. 各模型参数如表1所示. ...

... 在充电工况下, 为了研究区分充放电不同工况对SOC估计的影响, 把基于充电模型的SOC估计值与参照值之间的误差和基于充放电单一模型^[8]的SOC估计值与参照值之间的误差作对比, 如图6所示. 图中:SOCE为SOC估计误差. 计算充电模型和单一模型的最大误差(ME)、平均绝对误差(MAE)和均方根误差(RMSE), 如表2所示. ...

State-of-charge evaluation of supercapacitors

2017

State of charge estimate for super-capacitor based on sliding mode observer

2017

... 为了得到更加准确的SOC估计值, Fan等^[10]将SOC的定义扩充到所有电容存储的电荷, 利用滑模观测器来估计SOC, 在降低计算量的同时取得了较高的精度. 但是文中讨论的模型未考虑由于自放电现象产生的泄漏电流的影响, 这在短期内忽略是合理的, 而在长期使用或静置时可能会带来一定的累积误差. 另外, 超级电容器在充电和放电工况下的特性并不一致, 在不同的工况下建立不同的模型参数来表征其充放电特性可提高模拟精度^[11]. ...

... 引入文献[10]的SOC, 定义为 ...

... 把仿真得到的观测器估计的SOC值(包括基于有泄漏的模型和基于无泄漏的模型^[10]两种情况) 与参照值进行比较, 如图7所示. 计算有泄漏模型和无泄漏模型的ME、MAE和RMSE, 结果如表3所示. ...

2010

... 文献[11]针对模型参数选取不准给扩展卡尔曼滤波SOC估计带来误差的问题, 设计了变电流放电工况和变电流充电工况实验, SOC估计结果显示: 对于变电流充电工况, 使用充电模型参数的扩展卡尔曼滤波SOC 估计值更接近真实值, 而使用放电模型参数的SOC估计值则偏离真实值; 反之亦然. 因此, 电池充电模型参数和放电模型参数是不同的, 如果单独使用一种模型的参数进行SOC估计就会增大误差. ...

2010

A new parameters identification procedure for simplified double layer capacitor two-branch model

2010

... 本文在二分支等效模拟电路模型^[12]的基础上, 建立内部多电容端电压的状态空间模型, 将超级电容器的电流和电压分别作为可控制输入和可观测输出, 来计算内部多电容的端电压值; 再根据非线性观测器^[13]的思想, 用内部多电容的端电压来重构SOC. 针对充电和放电两种工况分别辨识模型参数, 每个模型包括有泄漏和无泄漏两种情况. 最后, 通过充放电实验讨论区分充放电不同工况以及泄漏电流对SOC估计精度的影响, 并验证非线性观测器的性能. ...

... 本文综合考虑了超级电容单体等效电路模型的准确性和实用性, 选用Faranda^[12]提出的二分支等效模拟电路模型, 如图1所示. 图中:

R_{0}

为等效串联电阻;

R_{2}

为电荷再分配电阻;

C_{1}

为微分电容;

C_{2}

为小孔隙电容(电荷再分配电容);

R_{l}

为泄漏电阻;

i_{1}

为短期支路电流;

i_{2}

为中长期支路电流;

i_{3}

为泄漏电流;i为输入的充放电电流;

v_{1}

为微分电容

C_{1}

两端的电压;

v_{2}

为小孔隙电容

C_{2}

两端的电压;v为单体的端电压. ...

A novel method for state of charge estimation of lithium-ion batteries using a nonlinear observer

2014

Characterization of double-layer capacitors for power electronics applications

2000

... 超级电容器有一对金属箔电极, 每个电极表面都有活性炭, 多孔的活性炭表面在电解液中与离子形成电双层, 可以决定超级电容器的容性. 超级电容器的物理化学特性比较复杂, 单一的电路器件不足以描述其充放电行为. 因此, 研究人员通常利用一些简化的等效电路模型来代替超级电容器^[14]. ...

A comparative study on the applicability of ultracapacitor models for electric vehicles under different temperatures

2017

... 流过泄漏电阻

R_{l}

的电流是自放电现象产生的泄漏电流,一般都比较小.由于泄漏电阻

R_{l}

远大于

R_{0} 和 R_{2},

所以泄漏电流比较小, 在某些情况下也可以忽略. 电容器的状态空间模型式(1)在无泄漏电阻情况下^[15]的各参数为 ...

A variable capacitance based modeling and power capability predicting method for ultracapacitor

2018

〈

〉