基于阶梯波的锂离子电池电化学阻抗谱低频段在线辨识方法

孙丙香1 王家驹2 苏晓佳1 张维戈1 赵鑫泽1

(1. 北京交通大学国家能源主动配电网技术研发中心 北京 100044 2. 国网宁夏电力有限公司银川供电公司 银川 750011)

摘要 电化学阻抗谱(EIS)是锂离子内部特性参数的外在表征,如何摆脱对专有设备的依赖,通过在线辨识得到EIS,尤其是低频段EIS对于准确地估计电池状态具有重要意义。该文以三元锂离子电池为研究对象,在电池管理系统和双向充放电机的协议框架下,设计了工况实验,利用阶梯信号和正弦信号的变换关系,提出了一种基于阶梯波的锂离子电池EIS低频段在线辨识方法。分析了不同荷电状态(SOC)和不同数据采集间隔对于阻抗谱辨识结果的影响,对比了该方法相对于方波傅里叶分解法的优势,给出了满足辨识精度的最低采样要求。该方法通过现有配备的电池外围系统即可实现,有效地解决了工程应用中无法实现锂离子电池EIS低频段在线辨识的难题。

关键词:锂离子电池 电化学阻抗谱 低频段 阶梯波 在线辨识

0 引言

锂离子电池的性能检测、状态评估和安全预警一直是使用过程中需要重点解决的问题[1-7],但在实际运行中,由于工况复杂、数据有限,人们难以准确地判断锂离子电池的状态,评估锂离子电池的性能。

电化学阻抗谱(Electrochemical Impedance Spectroscopy, EIS)是一种无损、有效的电化学检测工具,其测试频率很宽,被广泛应用于分析电极过程动力学、双电层和扩散现象,以及研究电极材料衰退、固体电解质形成、结构破损等老化机理[8-10]。锂离子电池在不同频率的激励下表现出不同的响应特性。如在高频段,电化学阻抗谱表现出电感特性,这主要是由电极的多孔性和表面不均匀等引起的,也是粘性系统的一个体现;在中低频段表现出电容特性,主要是由中频段的电荷传递过程和低频段的扩散过程引起的[11-13]。利用不同频段的阻抗谱信息对锂离子电池进行建模和参数辨识,能够充分挖掘电池的内部特征,更好地实现对电池的有效管理[14-19]。当前研究表明,EIS的中低频段数据能够有效实现锂离子电池健康状态估计[20-21],尤其是低频段能反映出的参数信息非常丰富,而且对采样频率的要求不苛刻。因此,在线获得EIS的低频段数据是实现在线主动健康状态评估的关键。

然而,传统的EIS辨识是通过专有设备输出正弦波激励测量得到,在实际应用中很难在线实现。目前工程上已有的在线辨识EIS的方法[22],主要是通过对方波进行傅里叶分解的方式得到相应的正弦波。然而,分解后得到的谐波幅值和基波幅值存在差异,其测量结果在EIS低频段的误差相对较大。因此,为进一步提高EIS低频段的辨识精度,在满足充放电机输出精度的前提下,结合我国现行的《电动汽车非车载传导式充电机与电池管理系统之间的通信协议》(GB/T 27930—2015)中规定的充电机与电池管理系统通信间隔50 ms的技术指标,本文提出一种基于阶梯波工况的在线辨识锂离子电池EIS低频段的方法,频率范围选定为2 Hz以下。

1 理论方法

本文提出的基于阶梯波的锂离子电池EIS低频段在线辨识方法示意图如图1所示,该方法涉及将正弦信号转换为阶梯信号的基础理论。

width=216.75,height=125.25

图1 在线辨识方法示意图

Fig.1 Schematic diagram of the identification method

首先,假设理想的正弦信号为

width=45.5,height=14.95 (1)

式中,width=45.5,height=12.9T为该正弦信号的周期;A为该正弦信号的幅值;width=10.2,height=12.25为原始正弦信号。将该正弦信号离散成每周期width=12.25,height=12.25个相等的间隔后,在离散点上,正弦信号的值可以表示为式(2)中的width=12.25,height=14.95[23]。同时,该正弦信号就可以用式(2)转换为阶梯信号。

width=107.3,height=29.2 (2)

式中,width=63.15,height=12.9,表示正弦信号离散后的第几个间隔,也表示阶梯信号中的第几个阶梯,其中width=12.25,height=12.25为阶梯数;width=12.25,height=14.95为第n个阶梯对应的幅值;width=12.9,height=12.25为参数,计算式为

width=35.3,height=26.5 (3)

运用上述公式,可以将EIS测试中的正弦波电流等效为阶梯波电流后,再进行EIS测试,具体操作方法如下。

首先,EIS测试时,先确定一个周期正弦波电流对应阶梯波电流的阶梯数width=12.25,height=12.25width=12.25,height=12.25一般可取为10,即将一个周期的正弦波分为10个相同的间隔。width=12.25,height=12.25取值越大,阶梯波电流等效为正弦波电流的效果越好。但阶梯数width=12.25,height=12.25越大,对应每一个阶梯的电流持续时间width=6.8,height=12.25就会越小,在实际工程中对设备精度的要求就会更高,不利于应用;若width=12.25,height=12.25取值太小,就无法有效地代替正弦波电流。

然后,根据所需等效的EIS测试的正弦波电流幅值,运用式(2)和式(3),确定阶梯波电流中每一个阶梯的电流幅值,即width=49.6,height=14.95。对锂离子电池施加阶梯波电流,根据实际需要确定施加阶梯波电流的频率范围,计算每一个频率下施加阶梯波电流的总时间width=10.2,height=12.25,进而确定不同频率下阶梯波电流中每一个阶梯的电流持续时间width=6.8,height=12.25,将每个频率下的阶梯波电流组合后得到变频阶梯波电流谱。若阶梯波电流的频率为f,则阶梯波电流的总时间width=10.2,height=12.25

width=29.2,height=27.85(4)

对应该频率下每一个阶梯的电流持续时间width=6.8,height=12.25

width=27.15,height=27.85(5)

对锂离子电池施加变频阶梯波电流谱,选取合适的采样频率,对阶梯波电流和电压响应进行采样,将采样得到的不同频率下的阶梯波电流和电压响应分别进行正弦拟合,得到不同频率下拟合的正弦电流和电压响应。各频率下的正弦拟合方式相同,频率为0.5 Hz下等效正弦波电流幅值为8 A的阶梯波电流和其正弦拟合电流如图2a所示,图2b为上述阶梯波电流激励下的阶梯波电压响应和正弦拟合电压,其拟合优度均达0.96。

width=206.25,height=342

图2 频率为0.5 Hz下等效正弦波电流幅值为8 A的阶梯波电流和电压响应

Fig.2 Current and voltage response with an equivalent sine wave current amplitude of 8 A at a frequency of 0.5 Hz

上述阶梯波电压响应为锂离子电池受到阶梯波电流激励后的电压变化值width=21.75,height=12.9,其表达式为

width=71.3,height=14.95 (6)

式中,UO为锂离子电池受到阶梯波电流激励后的端电压;UOCV为锂离子电池的开路电压。

最后,将对应频率下拟合后的锂离子电池正弦电压响应除以正弦电流,得到该频率下的复数阻抗,如式(7)~式(9)。进而组合不同频率下的阻抗值得到低频段的EIS。为了区别基于阶梯波和正弦波得到的两种EIS,本文将基于阶梯波得到的EIS简称为STEIS(Step Electrochemical Impedance Spectroscopy),将正弦波测试得到的EIS简称为SEIS(Sine Electrochemical Impedance Spectroscopy)。

width=65.9,height=27.85(7)

width=50.25,height=12.25 (8)

width=49.6,height=12.9 (9)

式中,width=19,height=12.25为拟合后的锂离子电池正弦电压响应的幅值;width=10.2,height=12.25为拟合后正弦响应电压的相角;width=8.15,height=12.25为拟合后正弦电流的幅值;width=12.25,height=10.2为拟合后正弦电流的相角;width=12.25,height=12.25为阻抗幅值;width=8.15,height=12.25为阻抗角;width=14.95,height=12.25为阻抗实部;Im为阻抗虚部。

2 实验平台及实验设计

2.1 测试对象及测试设备

本文测试对象为三元锰系的8 A·h软包电池,根据其产品手册列出样本电池的主要参数,见表1。

表1 电池参数

Tab.1 Battery parameters

参数数值 标称电压U/V3.60 额定容量Q/(A·h)8.00 充电截止电压Umax/V4.20 放电截止电压Umin/V2.75 标准充放电倍率1 C

本文采用的实验平台如图3所示,其中包括测试电化学阻抗谱的比奥罗杰VMP-300电化学工作站、控制环境温度的高低温试验箱、被测动力电池,以及用于控制和数据存储的计算机。

本次实验中进行的测试仅涉及基本的电池充放电实验和EIS实验,电化学工作站可以完全满足需求,因此无需使用其他电池充放电设备。

width=224.25,height=80.25

图3 实验平台

Fig.3 Laboratory platform

测试所选用的温箱为高低温试验箱,其温度范围为-20~80℃。为了防止外界温度变化对实验造成影响,在实验过程中均将电池置于温箱中进行测试。各设备的具体参数见表2和表3。

表2 VMP-300电化学工作站参数

Tab.2 The parameters of the VMP-300

参数数值 通道数/个4 控制电压/V±10 电压分辨率/μV1 最大电流/A±10 电流分辨率/fA760 阻抗测试10 μHz~3 MHz(1%,1°) 10 μHz~7 MHz(3%,3°)

表3 高低温试验箱参数

Tab.3 The parameters of thermostat

参数数值 内室尺寸/(mm×mm×mm)700´800´900 温度波动/℃<±0.5 温度均匀度/℃<±2 升温速率/(℃/min)1~3 降温速率/(℃/min)0.7~1 设定精度、指示精度/℃±0.1

2.2 实验设计

本文在研究过程中,主要涉及了电池的容量标定实验、SEIS实验及STEIS实验。

1)电池的容量标定实验具体实施方案如下:

(1)1C充电至电池单体上限电压,以该电压恒压充电至电流小于等于0.1 C,静置1 h。

(2)1C放电至电池单体下限电压,静置1 h。

(3)重复步骤(1)和步骤(2),共重复三次,取三次放电容量的均值作为电池的实际容量。

本次实验仅需在电池测试开始前进行一次,确定电池容量后,以标定实际容量作为后续测试依据。

2)为了能更好地分析本文提出的基于阶梯波的EIS低频段辨识方法,本次研究主要设计了三组实验,重点分析STEIS辨识方法的可行性、不同采样频率对STEIS低频段辨识结果的准确性,以及对比于方波傅里叶分解的EIS测试方法的优势性。测试频率均为10 mHz~2 Hz,并都在25℃下进行。

第一组实验是分别对20%荷电状态(State of Charge, SOC)、50%SOC、70%SOC、80%SOC、90%SOC和100%SOC下的电池,进行电流幅值均为8 A的SEIS测试和STEIS测试。在STEIS测试过程中,采样频率为1 kHz。通过该组实验,可以分析STEIS低频段在线辨识方法能否在不同SOC下实现低频段EIS的精确辨识。

第二组实验是对50%SOC下的电池,进行电流幅值为8 A的SEIS测试,并分别进行采样频率为 1 kHz、200 Hz、100 Hz、50 Hz的STEIS测试,电流幅值同样均为8 A。通过该组实验,可以对比分析在不同采集频率及误差下,STEIS低频段在线辨识方法的适用性。

第三组实验是在电池满电状态下,进行电流幅值为8 A的STEIS测试和基于方波傅里叶分解的EIS测试,并通过该实验对比这两种方法的辨识精度。方波傅里叶分解方法的具体步骤见文献[22]。

3 实验结果分析

3.1 STEIS测试方法在不同SOC下的可行性分析

电池在20%SOC、50%SOC、70%SOC、80%SOC、90%SOC和100%SOC下SEIS和STEIS的测试结果如图4所示。通过观察测试结果对比图可以发现,在不同SOC下,STEIS辨识方法在10 mHz~2 Hz频率范围内得到的STEIS均与SEIS吻合度极高。

width=219.75,height=176.25

图4 不同SOC下SEIS和STEIS测试结果

Fig.4 Test results of SEIS and STEIS under different SOCs

为了能更直观地分析测试结果,本文将不同SOC下SEIS和STEIS的实部、负虚部分别进行了对比,如图5和图6所示。显然,图5和图6中STEIS实部、负虚部与SEIS实部、负虚部都存在细微的偏差。通过对数据进一步计算,得到了不同SOC下STEIS相比SEIS实部与负虚部的相对误差,如图7、图8所示。

width=221.25,height=180.75

图5 不同SOC下的SEIS和STEIS实部

Fig.5 Real part of SEIS and STEIS under different SOCs

width=204.75,height=167.25

图6 不同SOC下的SEIS and STEIS负虚部

Fig.6 Negative imaginary part of SEIS and STEIS under different SOCs

从图7和图8中可以得到,STEIS测试方法在10 mHz~2 Hz频率范围内,无论电池SOC如何变化,STEIS实部的相对误差可以保持在1.75%以下,STEIS负虚部的相对误差基本可以保持在3%以下。由该数据也能得出,STEIS实部的相对误差较STEIS负虚部的相对误差更小一些。

同时,STEIS实部的相对误差随着测试频率width=12.25,height=14.95的改变并没有出现明显的规律性变化,而STEIS负虚部的相对误差随着测试频率width=12.25,height=14.95变大呈增大趋势。引起这个结果的主要原因是STEIS负虚部的值比STEIS实部的值小一个数量级。当STEIS负虚部出现误差时,误差相对其真实值就会占比更多。但从总体来看,无论是STEIS实部的相对误差,还是STEIS负虚部的相对误差,其大小都在可接受的范围之内。因此,本文提出的基于阶梯波的EIS低频段的在线辨识方法可以应用于电池的不同SOC状态。

width=200.25,height=164.25

图7 不同SOC下STEIS实部的相对误差

Fig.7 Relative error of the real part of STEIS under different SOCs

width=198,height=162

图8 不同SOC下STEIS负虚部的相对误差

Fig.8 Relative error of the negative imaginary part of STEIS under different SOCs

3.2 采样频率对STEIS低频段辨识结果的准确性影响

为了验证采样精度对该方法的影响,本文在电池为50%SOC的状态下,分别完成了采样频率width=12.25,height=14.95为1 000、200、100、50 Hz的STEIS测试,测试频率为10 mHz~2 Hz。同时,将上述采样频率下得到的STEIS与SEIS测试结果进行对比,得到结果如图9所示。

width=207,height=171.75

图9 不同采样频率下SEIS and STEIS测试结果

Fig.9 Test results of SEIS and STEIS at different sampling frequencies

从图9的对比结果可以发现,采样精度会影响STEIS的测试结果,但在采样频率为50~1 000 Hz范围内的测试结果均能接受。采样频率越高,STEIS测试方法得到的结果就越接近SEIS的结果。当采样频率width=12.25,height=14.95=1 000 Hz时,STEIS和SEIS基本重合。但当采样频率width=12.25,height=14.95=50 Hz时,STEIS中在不同测试频率点下的阻抗均出现了一些偏移。为了更好地分析STEIS测试方法受采样频率的影响,将不同采样频率下的STEIS实部、负虚部和SEIS的实部、负虚部分别进行了对比,如图10和图11所示。

width=207.75,height=171

图10 不同采样频率下SEIS and STEIS实部

Fig.10 Real part of SEIS and STEIS at different sampling frequencies

从图10和图11的对比结果可以发现,采样频率的降低对STEIS的实部和负虚部的影响都十分有限。不同采样频率下STEIS实部和负虚部相比于SEIS实部和负虚部的相对误差分别如图12和图13所示。通过计算STEIS相比SEIS实部的相对误差,可知,采样频率在50~1 000 Hz范围内时,随着采样频率从1 000 Hz降低到50 Hz,STEIS中每个测试频率点下的阻抗实部的相对误差会小幅增加,但STEIS中阻抗实部的相对误差均能保持在1.4%以内。同时,通过计算STEIS相比SEIS负虚部的相对误差,可知采样频率width=12.25,height=14.95=1 000 Hz时,STEIS的负虚部相对误差可以保持在3%以内。当采样频率width=12.25,height=14.95降低为200、100、50 Hz时,STEIS负虚部的最大相对误差可以达到6%~7%之间。但10 mHz~0.8 Hz频率范围内的阻抗虚部误差在采样频率为50 Hz时,仍能保持在3%以内。由于STEIS中阻抗实部占总阻抗比重在95%以上,因此STEIS中阻抗负虚部的误差变化可以忽略。对于在线辨识EIS,上述误差仍处于可接受范围内。

width=207.75,height=170.25

图11 不同采样频率下SEIS and STEIS负虚部对比

Fig.11 Negative imaginary part of SEIS and STEIS at different sampling frequencies

width=213,height=173.25

图12 不同采样频率下STEIS实部相对误差

Fig.12 Relative error of the real part of STEIS at different sampling frequencies

width=200.25,height=164.25

图13 不同采样频率下STEIS负虚部相对误差

Fig.13 Relative error of the negative imaginary part of STEIS at different sampling frequencies

综上所述,STEIS辨识方法受限于采样频率,即采样频率越高,STEIS结果越好。这种情况的出现是因为在进行STEIS测试时,所测阻抗对应的频率越高,相同采样频率下得到的采样点就越少。这个结果就导致阶梯波电流和电压响应在进行正弦拟合时,采样数据点的减少使得拟合后的正弦电流和电压误差增大,其中拟合后的正弦电压误差变化更为明显,致使最终计算得到的阻抗误差也随之变大。所以在同一采样频率下,阻抗对应的频率越低,STEIS辨识方法效果就越好。

因此,在使用阶梯波进行EIS低频段辨识时,其采样频率需要根据所测频率点的变化,选取不同的采样频率。所测频率点越高,采样频率点就要相应提高。当前测试结果表明,在进行2 Hz频率点的阻抗测试时,采样频率至少要达到50 Hz。

3.3 对比STEIS方法和方波傅里叶分解EIS的辨识精度

图14给出了STEIS辨识方法和方波傅里叶分解法EIS低频段的测量结果,阶梯波信号和方波信号的采样频率分别为1 000 Hz和100 Hz。通过计算得到STEIS实部和负虚部的平均相对误差分别为1.25%和2.8%,方波傅里叶分解法EIS实部和负虚部的平均相对误差分别为5.29%和11.62%。经计算,STEIS的平均相对误差比方波的测量结果小6.4%。所以即使方波具有较高的采样频率,但阶梯波得到的结果明显优于方波。STEIS方法测量阻抗谱时各个频率的幅值一致,而方波傅里叶分解的方法在不同频率处的幅值存在差异,因此STEIS方法优于方波方法。

width=213,height=171.75

图14 方波EIS和STEIS测试结果

Fig.14 Test results of Square EIS and STEIS

4 结论

采用十阶梯的阶梯波能够实现正弦波的精准拟合,拟合优度为0.96,进而能够实现锂离子电池EIS的在线辨识。从辨识结果分析可以发现,在不同的SOC下,只要采样频率足够高,就可以实现较宽频域的EIS测试。实验结果表明,采样频率为50 Hz时,在10 mHz~2 Hz测试频率范围内可得到高精度的STEIS结果,平均相对误差在1.12%以内。STEIS方法与方波傅里叶分解的方法进行对比,测量频域的平均相对误差比方波小6.4%。由此可见,STEIS方法明显优于方波。本文提出的方法能够摆脱对正弦波设备的依赖,在现有的通信协议框架内,通过电池管理系统控制双向充放电设备就可以实现,为电池的运行状态在线主动评估提供了重要支撑。

参考文献

[1] 王榘, 熊瑞, 穆浩. 温度和老化意识融合驱动的电动车辆锂离子动力电池电量和容量协同估计[J]. 电工技术学报, 2020, 35(23): 4980-4987. Wang Ju, Xiong Rui, Mu Hao. Co-estimation of lithium-ion battery state-of-charge and capacity through the temperature and aging awareness model for electric vehicles[J]. Transactions of China Electrotechnical Society, 2020, 35(23): 4980-4987.

[2] 赵伟, 闵婕, 李章溢, 等. 基于一致性模型的梯次利用锂离子电池组能量利用率估计方法[J]. 电工技术学报, 2021, 36(10): 2190-2198. Zhao Wei, Min Jie, Li Zhangyi, et al. Energy utilization efficiency estimation method for second-use lithium-ion battery packs based on a battery consistency model[J]. Transactions of China Electrotechnical Society, 2021, 36(10): 2190-2198.

[3] 冯飞, 逯仁贵, 朱春波. 一种锂离子电池低温SOC估计算法[J]. 电工技术学报, 2014, 29(7): 53-58. Feng Fei, Lu Rengui, Zhu Chunbo. State of charge estimation of Li-ion battery at low temperature[J]. Transactions of China Electrotechnical Society, 2014, 29(7): 53-58.

[4] 刘伟龙, 王丽芳, 王立业. 基于电动汽车工况识别预测的锂离子电池SOE估计[J]. 电工技术学报, 2018, 33(1): 17-25. Liu Weilong, Wang Lifang, Wang Liye. Estimation of state-of-energy for electric vehicles based on the identification and prediction of driving condition[J]. Transactions of China Electrotechnical Society, 2018, 33(1): 17-25.

[5] Dai Haifeng, Yu Chenchen, Wei Xuezhe, et al. State of charge estimation for lithium-ion pouch batteries based on stress measurement[J]. Energy, 2017, 129: 16-27.

[6] Hu Xiaosong, Xu Le, Lin Xianke, et al. Battery lifetime prognostics[J]. Joule, 2020, 4(2): 310-346.

[7] Xiong Rui, Sun Wanzhou, Yu Quanqing, et al. Research progress, challenges and prospects of fault diagnosis on battery system of electric vehicles[J]. Applied Energy, 2020, 279: 115855.

[8] 吴磊, 吕桃林, 陈启忠, 等. 电化学阻抗谱测量与应用研究综述[J]. 电源技术, 2021, 45(9): 1227-1230. Wu Lei, Lü Taolin, Chen Qizhong, et al. Review of measurement and application of electrochemical impedance spectroscopy[J]. Chinese Journal of Power Sources, 2021, 45(9): 1227-1230.

[9] 庄全超, 杨梓, 张蕾, 等. 锂离子电池的电化学阻抗谱分析研究进展[J]. 化学进展, 2020, 32(6): 761-791. Zhuang Quanchao, Yang Zi, Zhang Lei, et al. Research progress on diagnosis of electrochemical impedance spectroscopy in lithium ion batteries[J]. Progress in Chemistry, 2020, 32(6): 761-791.

[10] 黄秋安, 李伟恒, 汤哲鹏, 等. 电化学阻抗谱基础[J].自然杂志, 2020, 42(1): 12-26. Huang Qiuan, Li Weiheng, Tang Zhepeng, et al. Fundamentals of electrochemical impedance spectroscopy[J]. Chinese Journal of Nature, 2020, 42(1): 12-26.

[11] 韩慧颖, 冯素蕊, 陶钰禧. 电化学阻抗谱的原理及其应用[J]. 缔客世界, 2020, 6(6): 121. Han Huiying, Feng Surui, Tao Yuxi. Principle and application of electrochemical impedance spectroscopy[J]. Maker world, 2020, 6(6): 121.

[12] 王盼. 电化学阻抗谱在锂离子电池中的应用[J]. 电源技术, 2020, 44(12): 1847-1850, 1854. Wang Pan. Application of electrochemical impedance spectroscopy in lithium ion batteries[J]. Chinese Journal of Power Sources, 2020, 44(12): 1847-1850, 1854.

[13] 冷晓伟, 戴作强, 郑莉莉, 等. 锂离子电池电化学阻抗谱研究综述[J]. 电源技术, 2018, 42(11): 1749-1752. Leng Xiaowei, Dai Zuoqiang, Zheng Lili, et al. Review on electrochemical impedance spectroscopy of lithium-ion batteries[J]. Chinese Journal of Power Sources, 2018, 42(11): 1749-1752.

[14] 龚敏明, 卞景季, 孙丙香, 等. 锂离子电池分数阶等效电路模型低频参数演变规律研究[J]. 重庆理工大学学报(自然科学), 2020, 34(2): 6-14. Gong Minming, Bian Jingji, Sun Bingxiang, et al. Study on parameter evolution of fractional order equivalent circuit model for Li-ion batteries in low frequency area[J]. Journal of Chongqing University of Technology (Natural Science), 2020, 34(2): 6-14.

[15] 孙丙香, 刘佳, 韩智强, 等. 不同区间衰退路径下锂离子电池的性能相关性及温度适用性分析[J]. 电工技术学报, 2020, 35(9): 2063-2073. Sun Bingxiang, Liu Jia, Han Zhiqiang, et al. Performance correlation and temperature applicability of Li-ion batteries under different range degradation paths[J]. Transactions of China Electrotechnical Society, 2020, 35(9): 2063-2073.

[16] 戴海峰, 王冬晨, 姜波. 基于电化学阻抗谱的电池荷电状态估计[J]. 同济大学学报(自然科学版), 2019, 47(S1): 95-98. Dai Haifeng, Wang Dongchen, Jiang Bo. Estimation of battery state of charge based on electrochemical impedance spectroscopy[J]. Journal of Tongji University (Natural Science), 2019, 47(S1): 95-98.

[17] 吴健, 尹泽, 李豪, 等. 基于分数阶理论的锂离子电池高频等效电路模型[J]. 电工技术学报, 2021, 36(18): 3902-3910. Wu Jian, Yin Ze, Li Hao, et al. High-frequency equivalent circuit model of lithium-ion battery based on fractional order theory[J]. Transactions of China Electrotechnical Society, 2021, 36(18): 3902-3910.

[18] 李晓宇, 朱春波, 魏国, 等. 基于分数阶联合卡尔曼滤波的磷酸铁锂电池简化阻抗谱模型参数在线估计[J]. 电工技术学报, 2016, 31(24): 141-149. Li Xiaoyu, Zhu Chunbo, Wei Guo, et al. Online parameter estimation of a simplified impedance spectroscopy model based on the fractional joint Kalman filter for LiFePO4 battery[J]. Transactions of China Electrotechnical Society, 2016, 31(24): 141-149.

[19] 张彩萍, 姜久春, 张维戈, 等. 梯次利用锂离子电池电化学阻抗模型及特性参数分析[J]. 电力系统自动化, 2013, 37(1): 54-58. Zhang Caiping, Jiang Jiuchun, Zhang Weige, et al. Characterization of electrochemical impedance equivalent model and parameters for Li-ion batteries echelon use[J]. Automation of Electric Power Systems, 2013, 37(1): 54-58.

[20] 张利中, 穆苗苗, 赵书奇, 等. 再利用退役锂动力电池的性能评估[J]. 电源技术, 2018, 42(7): 964-967. Zhang Lizhong, Mu Miaomiao, Zhao Shuqi, et al. Performance assessments of retired lithium-ion power batteries for reuse[J]. Chinese Journal of Power Sources, 2018, 42(7): 964-967.

[21] 胡国荣, 彭清远, 彭忠东, 等. 两种方法制备的磷酸铁锂/石墨烯复合材料的性能对比[J]. 无机化学学报, 2015, 31(6): 1153-1158. Hu Guorong, Peng Qingyuan, Peng Zhongdong, et al. Comparison on properties of lithium iron phosphate/graphene composite prepared by two methods[J]. Chinese Journal of Inorganic Chemistry, 2015, 31(6): 1153-1158.

[22] 沈迪, 阮海军, 姜久春, 等. 一种锂离子电池的EIS快速测量方法: CN106970266A[P]. 2017-07-21.

[23] 陆祖良, 杨雁, 黄璐, 等. 阶梯波性质的进一步探讨: 阶梯波研究之一[J]. 计量学报, 2018, 39(6): 759-767. Lu Zuliang, Yang Yan, Huang Lu, et al. Further discussion on characteristics of staircase waveform[J]. Acta Metrologica Sinica, 2018, 39(6): 759-767.

Study on Online Identification Method of Low Frequency Electrochemical Impedance Spectroscopy for Lithium-Ion Battery Based on Step Wave

Sun Bingxiang1 Wang Jiaju2 Su Xiaojia1 Zhang Weige1 Zhao Xinze1

(1. National Active Distribution Network Technology Research Center Beijing Jiaotong University Beijing 100044 China 2. Yinchuan Power Supply Company State Grid Ningxia Power Co. Ltd Yinchuan 750011 China)

Abstract Electrochemical impedance spectroscopy (EIS) is an external characterization of the internal performance parameters of lithium-ion battery. How to get rid of the dependence on proprietary equipment for online identification of EIS, especially low-frequency EIS identification is of great significance for accurately estimating battery status.

At present, the existing methods of online identification of EIS in engineering mainly obtain the corresponding sine wave by Fourier decomposition of the square wave. However, there are differences between the harmonic amplitudes obtained after decomposition and the fundamental amplitudes, and the identification results have relatively large errors in the low frequency band of EIS. Therefore, in order to further improve the identification accuracy of the EIS low frequency band, this paper takes NCM lithium-ion battery as the research object in the framework of the communication protocols between battery management system (BMS) and bi-directional charging-discharging device. Further, an online identification method for lithium-ion batteries low-frequency EIS based on step waves is proposed via the transformation relationship of step signals and sine signals, the frequency range is selected to be below 2 Hz.

In order to better analyze the EIS low-frequency identification method based on step wave proposed in this paper, three groups of experiments are mainly designed. The feasibility of the step wave EIS (STEIS) identification method, the accuracy of the STEIS low-frequency identification results with different sampling frequencies, and the advantages of the EIS test method compared with the square wave Fourier decomposition are mainly analyzed. The test frequency is 10 mHz~2 Hz, and the temperature is all carried out at 25℃. The equipment and components used in the experiment include Bio-Logic VMP-300 electrochemical workstation for electrochemical impedance spectroscopy, high and low temperature test chamber for controlling ambient temperature, power battery under test, and computer for control and data storage.

The experimental results show that the ten-step step wave can realize the accurate fitting of the sine wave, and the goodness of fit is 0.96, which can realize the online identification of the EIS of the lithium-ion battery. From the analysis of the identification results, it can be found that at different SOC points, as long as the sampling frequency is high enough, the EIS test in a wider frequency domain can be realized. When the sampling frequency is 50 Hz, high-precision STEIS results can be obtained in the test frequency range of 10 mHz~2 Hz, and the average relative error is within 1.12%. Comparing the STEIS method with the square wave Fourier decomposition method, the average relative error in the measurement frequency domain is 6.4% smaller than that of the square wave. It can be seen that the STEIS method is significantly better than the square wave. The method proposed in this paper can get rid of the dependence on the sine wave equipment. Within the framework of the existing communication protocol, it can be realized by controlling the bidirectional charging and discharging equipment through the battery management system, which provides an important method support for the online active evaluation of the operating status of the battery.

keywords:Lithium-ion battery, electrochemical impedance spectroscopy (EIS), low frequency, step wave, online identification

DOI:10.19595/j.cnki.1000-6753.tces.220280

中图分类号:TM911

国家自然科学基金资助项目(52177206,51907005)。

收稿日期 2022-02-18

改稿日期 2022-07-16

作者简介

孙丙香 女,1979年生,博士,教授,研究方向为锂离子电池成组应用技术,包括建模与仿真、热特性与低温加热策略、优化充电、储能电池性能评估与经济性分析等。E-mail:bxsun@bjtu.edu.cn(通信作者)

王家驹 男,1995年生,硕士,研究方向为锂离子电池状态估计。E-mail:wangjiaju@bjtu.edu.cn

(编辑 李冰)