CN113702836A - 一种基于emd-gru锂离子电池荷电状态估计方法 - Google Patents

Abstract

本发明公开了一种基于EMD‑GRU锂离子电池荷电状态估计方法，具体为：步骤1，对锂离子电池进行脉冲电流放电，直至电池电压降到放电截止电压时停止放电，采集放电电流、电池电压、温度及SOC(t)时间序列；步骤2，采用经验模态分解EMD算法对采集放电电流进行分解，将电流时间序列分解成不同频率的子电流集时间序列和残差电流时间序列；步骤3，将子电流集时间序列和残差电流时间序列以及电池电压、温度进行归一化处理；步骤4，基于门控循环单元GRU建立子电流集、电压和温度时间序列的SOC估计模型并进行锂离子电池SOC估计。本发明解决了现有技术中存在的循环神经网络锂离子电池SOC估计丢失中长期电流时间序列的问题。

Description

一种基于EMD-GRU锂离子电池荷电状态估计方法

技术领域

本发明属于锂电池状态估计方法技术领域，涉及一种基于EMD-GRU锂离子电池荷电状态估计方法。

背景技术

锂离子电池因其能量密度高，自放电率低，无记忆效应等优点，被广泛电动汽车及各个领域，逐渐成为未来诸多重要领域的关键和支撑。然而锂离子电池作为电动汽车核心技术仍面临者诸多瓶颈，续航里程及安全性能成为用户关注的重点。为了提高电动汽车续航里程及安全性能，国内外的高校、公司及科研机构都将电池管理***(BatteryManagement System，BMS)技术作为重要的发展方向。

锂离子电池荷电状态(State of Charge,SOC)精准评估及预测是BMS的关键功能之一。精确的SOC值能够表征电池电量使用情况和充放电程度，制定理想的电池充放电策略；单体电池由于制造工艺的差异，在使用过程中的电量变化不一致，进而影响整个电池组的工作效率，需要通过控制每个电池单体的SOC值平衡单体电池间的差异；此外，SOC影响电池管理决策的制定，是电池管理***分配电量的依据。然而电池在使用过程中，其负载变化大，使用环境多变，并且SOC是电池的内部状态参数，无法直接测量得到，使得SOC估计问题一直是研究热点。

工程应用中通常采用开路电压法、安时积分法(库伦积分法)或二者相结合的方法估计SOC，但是估计精度不高。近几年，学者提出循环神经网络的方法(RNN、LSTM和GRU)应用于锂离子电池SOC估计当中，神经网络法高度非线性特点和电池相契合，故可以很好的应用。输入层的输入量一般为电池的电压，电流，电池温度，电池内阻，循环次数等，输出为电池的SOC，但循环神经网络的锂离子电池SOC估计在面对较长时间序列时，会出现丢失部分中长期序列信息的情况，影响模型预测精度。

发明内容

本发明的目的是提供一种基于EMD-GRU锂离子电池荷电状态估计方法，解决了现有技术中存在的循环神经网络锂离子电池SOC估计丢失中长期电流时间序列的问题。

本发明所采用的技术方案是，一种基于EMD-GRU锂离子电池荷电状态估计方法，具体按照如下步骤实施：

步骤1，对锂离子电池进行脉冲电流放电，直至电池电压降到放电截止电压时停止放电，采集放电电流i(t)、电池电压u(t)、温度T(t)及SOC(t)时间序列；

步骤2，采用经验模态分解EMD算法对步骤1采集的采集放电电流i(t)进行分解，将电流时间序列分解成不同频率的子电流集时间序列和残差电流时间序列；

步骤3，将步骤2得到的子电流集时间序列和残差电流时间序列以及电池电压u(t)、温度T(t)进行归一化处理；

步骤4，基于门控循环单元GRU建立子电流集、电压和温度时间序列的SOC估计模型并进行锂离子电池SOC估计。

本发明的特征还在于，

步骤1对锂离子电池进行脉冲电流放电前，将锂离子电池电量充满，即SOC＝100％，锂离子电池静置在恒温恒湿箱中，静置时间为1-3小时，然后进行脉冲电流放电。

步骤1中进行采样的采样频率为脉冲电流最大频率的十倍及以上。

步骤2具体为：

步骤2.1，找到放电电流i(t)的所有极大值和极小值点，采用三次样条插值得到上下包络线e_max(t)和e_min(t)，计算上下包络线的均值m(t)：

步骤2.2，将放电电流i(t)与上下包络平均值m(t)进行差值，获得剩余分量h₁(t)：

h₁(t)＝i(t)-m(t) (2)

步骤2.3，判断h₁(t)是否满足固有模态分量的两个条件，不满足的话，用h₁(t)作为新的时间序列返回步骤2.1，直到剩余分量h₁(t)满足固有模态分量的两个条件为止，固有模态分量c₁(t)记为：

c₁(t)＝h₁(t) (3)；

其中，固有模态分量的两个基本条件：1)h₁(t)函数局部极值点和过零点的个数相等或者相差不超过一个；2)在任意时刻范围内，由局部最大值和最小值所形成的包络曲线的平均值必须为零；

步骤2.4，将放电电流i(t)减去固有模态分量c₁(t)，求出残差r₁(t)：

r₁(t)＝i(t)-c₁(t) (4)；

步骤2.5，将残差r₁(t)作为新的信号重复步骤2.1～2.4，分别得到c₂(t)和r₂(t)，如此不断循环，当残差满足终止条件后结束，得到了固有模态分量c₁(t)、c₂(t)、…、c_n(t)即为不同频率的子电流集，最后得到的r_n(t)即为残差电流；

其中，终止条件为：当最后的残差r_n(t)是单调递增或者单调递减信号或者残差幅值不大于信号幅度的百分之一。

步骤3具体为：

子电流集c₁(t)、c₂(t)、…、c_n(t)、残差电流r_n(t)、电压u(t)、温度T(t)构成输入集x(t)，SOC(t)构成输出集y(t)，分别对输入集各个元素采用公式(5)进行归一化，归一化公式如下：

式中，x′(t)为各元素归一化后输入集，a为各个元素的平均值，b为各个元素的标准差。

步骤4具体为：

步骤4.1，采用输入集x′(t)作为GRU网络的输入，输出集y(t)作为GRU网络的输出，GRU模型公式如下：

式中，x′(t)为第t时刻隐含层的输入，h(t)为第t时刻当前层的输出，h(t-1)为第t-1时刻当前层的输出，z(t)和r(t)分别为更新门和重置门，

是输入x′(t)和前一时刻输出h(t-1)的汇总，σ为sigmoid激活函数，tanh为双曲正切函数，W_r为重置门训练参数矩阵，W_z为更新门训练参数矩阵，其中[]表示两个向量相连，*表示矩阵的乘积；

步骤4.2，采用反向传播算法BPTT算法对GRU模型进行训练，得到误差最小的参数矩阵即最优GRU模型，误差函数如下：

式中，y(t)为第t时刻实际输出，W_O为输出层参数矩阵；

步骤4.3，重新按照步骤1、步骤2和步骤3方法得到后一段时间锂离子电池子电流集c′₁(t)、c′₂(t)、…、c′_n(t)、残差电流r′_n(t)、电压u′(t)、温度T′(t)作为锂离子电池SOC估计输入集，结合步骤4.2中最优GRU模型，得到模型输出即为SOC估计值。

本发明的有益效果是：

本发明首先在脉冲放电下采集锂离子电池电压、电流、温度和SOC数据；其次采用EMD对电流数据进行预处理，把电流数据分解为子电流集和残差电流，残差电流呈现平稳趋势可以提高中长期时间序列预测的精度；最后基于GRU建立子电流、电压和温度时间序列的SOC估计模型，对锂离子电池SOC估计精度高，实现动态估计，解决循环神经网络锂离子电池SOC估计丢失中长期电流时间序列的问题。

附图说明

图1是本发明一种基于EMD-GRU锂离子电池SOC估计方法的总流程图；

图2是本发明一种基于EMD-GRU锂离子电池SOC估计方法脉冲电流放电图；

图3是本发明一种基于EMD-GRU锂离子电池SOC估计方法EMD电流时间序列分解流程图；

图4是本发明一种基于EMD-GRU锂离子电池SOC估计方法GRU结构框图；

图5是本发明一种基于EMD-GRU锂离子电池SOC估计方法电流分布和测量电压实例测量图；

图6是本发明一种基于EMD-GRU锂离子电池SOC估计方法估计结果图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

本发明一种基于EMD-GRU锂离子电池荷电状态估计方法，其流程如图1所示，具体按照如下步骤实施：

步骤1，对锂离子电池进行脉冲电流放电前，将锂离子电池电量充满，即SOC＝100％，锂离子电池静置在恒温恒湿箱中，静置时间为1-3小时，然后进行脉冲电流放电，直至电池电压降到放电截止电压时停止放电，采集放电电流i(t)、电池电压u(t)、温度T(t)及SOC(t)时间序列，采样频率为脉冲电流最大频率的十倍及以上；脉冲激励信号的选取原则是对被测对象充分激励，同时保证不会对电池造成损害。因此采用脉冲电流对锂离子电池进行循环放电，脉冲电流如图2所示，实验结束后，采用恒流恒压将电池充满电，以备后续使用。

步骤2，采用经验模态分解EMD算法对步骤1采集的采集放电电流i(t)进行分解，将电流时间序列分解成不同频率的子电流集时间序列和残差电流时间序列，电流时间序列是影响SOC最重要的成分，故本发明分解电流时间序列；如图3所示，具体为：

步骤2.1，找到放电电流i(t)的所有极大值和极小值点，采用三次样条插值得到上下包络线e_max(t)和e_min(t)，计算上下包络线的均值m(t)：

步骤2.2，将放电电流i(t)与上下包络平均值m(t)进行差值，获得剩余分量h₁(t)：

h₁(t)＝i(t)-m(t) (2)

步骤2.3，判断h₁(t)是否满足固有模态分量的两个条件，不满足的话，用h₁(t)作为新的时间序列返回步骤2.1，直到剩余分量h₁(t)满足固有模态分量的两个条件为止，固有模态分量c₁(t)记为：

c₁(t)＝h₁(t) (3)；

其中，固有模态分量的两个基本条件：1)h₁(t)函数局部极值点和过零点的个数相等或者相差不超过一个；2)在任意时刻范围内，由局部最大值和最小值所形成的包络曲线的平均值必须为零；

步骤2.4，将放电电流i(t)减去固有模态分量c₁(t)，求出残差r₁(t)：

r₁(t)＝i(t)-c₁(t) (4)；

步骤2.5，将残差r₁(t)作为新的信号重复步骤2.1～2.4，分别得到c₂(t)和r₂(t)，如此不断循环，当残差满足终止条件后结束，得到了固有模态分量c₁(t)、c₂(t)、…、c_n(t)即为不同频率的子电流集，最后得到的r_n(t)即为残差电流；

其中，终止条件为：当最后的残差r_n(t)是单调递增或者单调递减信号或者残差幅值不大于信号幅度的百分之一；

步骤3，将步骤2得到的子电流集时间序列和残差电流时间序列以及电池电压u(t)、温度T(t)进行归一化处理；具体为：

子电流集c₁(t)、c₂(t)、…、c_n(t)、残差电流r_n(t)、电压u(t)、温度T(t)构成输入集x(t)，SOC(t)构成输出集y(t)，分别对输入集各个元素采用公式(5)进行归一化，归一化公式如下：

式中，x′(t)为各元素归一化后输入集，a为各个元素的平均值，b为各个元素的标准差；

步骤4，基于门控循环单元GRU建立子电流集、电压和温度时间序列的SOC估计模型并进行锂离子电池SOC估计，如图4所示，具体为：

步骤4.1，采用输入集x′(t)作为GRU网络的输入，输出集y(t)作为GRU网络的输出，GRU模型公式如下：

式中，x′(t)为第t时刻隐含层的输入，h(t)为第t时刻当前层的输出，h(t-1)为第t-1时刻当前层的输出，z(t)和r(t)分别为更新门和重置门，

是输入x′(t)和前一时刻输出h(t-1)的汇总，σ为sigmoid激活函数，tanh为双曲正切函数，W_r为重置门训练参数矩阵，W_z为更新门训练参数矩阵，其中[]表示两个向量相连，*表示矩阵的乘积；

步骤4.2，采用反向传播算法BPTT算法对GRU模型进行训练，得到误差最小的参数矩阵即最优GRU模型，误差函数如下：

式中，y(t)为第t时刻实际输出，W_O为输出层参数矩阵；

步骤4.3，重新按照步骤1、步骤2和步骤3方法得到后一段时间锂离子电池子电流集c′₁(t)、c′₂(t)、…、c′_n(t)、残差电流r′_n(t)、电压u′(t)、温度T′(t)作为锂离子电池SOC估计输入集，结合步骤4.2中最优GRU模型，得到模型输出即为SOC估计值。

实施例

本发明采用美国马里兰大学高级生命周期工程中心的数据仓库中的公共数据集。试验中，将电池A123放置在温度室中，测量电池的温度。电流分布和测量电压如图5所示。很明显，动态应力测试DST数据集(图5顶部)与联邦城市驾驶时间表FUDS数据集(图5中部)和US06数据集(图5底部)都不同。US06和FUDS与DST在放电电流和电压上有明显的差异，因此以DST和FUDS数据集作为训练数据集，US06和作为测试数据集，可以很好地验证GRU网络的泛化能力。DST配置文件的温度为0℃、10℃、25℃、30℃、40℃、50℃，与FUDS和US06配置文件相同。

采用GRU模型，输入维度为8，输出维度为1，隐含层有两层节点个数为32，batch_size为100。如图6所示，为EMD-GRU锂离子电池SOC估计图，左边是训练，右边是测试，测试最终MSE误差<0.5％，满足估计需求。

Claims (6)

1.一种基于EMD-GRU锂离子电池荷电状态估计方法，其特征在于，具体按照如下步骤实施：

步骤1，对锂离子电池进行脉冲电流放电，直至电池电压降到放电截止电压时停止放电，采集放电电流i(t)、电池电压u(t)、温度T(t)及SOC(t)时间序列；

步骤2，采用经验模态分解EMD算法对步骤1采集的采集放电电流i(t)进行分解，将电流时间序列分解成不同频率的子电流集时间序列和残差电流时间序列；

步骤3，将步骤2得到的子电流集时间序列和残差电流时间序列以及电池电压u(t)、温度T(t)进行归一化处理；

步骤4，基于门控循环单元GRU建立子电流集、电压和温度时间序列的SOC估计模型并进行锂离子电池SOC估计。

2.根据权利要求1所述的一种基于EMD-GRU锂离子电池荷电状态估计方法，其特征在于，所述步骤1对锂离子电池进行脉冲电流放电前，将锂离子电池电量充满，即SOC＝100％，锂离子电池静置在恒温恒湿箱中，静置时间为1-3小时，然后进行脉冲电流放电。

3.根据权利要求2所述的一种基于EMD-GRU锂离子电池荷电状态估计方法，其特征在于，所述步骤1中进行采样的采样频率为脉冲电流最大频率的十倍及以上。

4.根据权利要求1所述的一种基于EMD-GRU锂离子电池荷电状态估计方法，其特征在于，所述步骤2具体为：

步骤2.1，找到放电电流i(t)的所有极大值和极小值点，采用三次样条插值得到上下包络线e_max(t)和e_min(t)，计算上下包络线的均值m(t)：

步骤2.2，将放电电流i(t)与上下包络平均值m(t)进行差值，获得剩余分量h₁(t)：

h₁(t)＝i(t)-m(t) (2)

步骤2.3，判断h₁(t)是否满足固有模态分量的两个条件，不满足的话，用h₁(t)作为新的时间序列返回步骤2.1，直到剩余分量h₁(t)满足固有模态分量的两个条件为止，固有模态分量c₁(t)记为：

c₁(t)＝h₁(t) (3)；

其中，固有模态分量的两个基本条件：1)h₁(t)函数局部极值点和过零点的个数相等或者相差不超过一个；2)在任意时刻范围内，由局部最大值和最小值所形成的包络曲线的平均值必须为零；

步骤2.4，将放电电流i(t)减去固有模态分量c₁(t)，求出残差r₁(t)：

r₁(t)＝i(t)-c₁(t) (4)；

步骤2.5，将残差r₁(t)作为新的信号重复步骤2.1～2.4，分别得到c₂(t)和r₂(t)，如此不断循环，当残差满足终止条件后结束，得到了固有模态分量c₁(t)、c₂(t)、…、c_n(t)即为不同频率的子电流集，最后得到的r_n(t)即为残差电流；

其中，终止条件为：当最后的残差r_n(t)是单调递增或者单调递减信号或者残差幅值不大于信号幅度的百分之一。

5.根据权利要求4所述的一种基于EMD-GRU锂离子电池荷电状态估计方法，其特征在于，所述步骤3具体为：

子电流集c₁(t)、c₂(t)、…、c_n(t)、残差电流r_n(t)、电压u(t)、温度T(t)构成输入集x(t)，SOC(t)构成输出集y(t)，分别对输入集各个元素采用公式(5)进行归一化，归一化公式如下：

式中，x′(t)为各元素归一化后输入集，a为各个元素的平均值，b为各个元素的标准差。

6.根据权利要求5所述的一种基于EMD-GRU锂离子电池荷电状态估计方法，其特征在于，所述步骤4具体为：

步骤4.1，采用输入集x′(t)作为GRU网络的输入，输出集y(t)作为GRU网络的输出，GRU模型公式如下：

式中，x′(t)为第t时刻隐含层的输入，h(t)为第t时刻当前层的输出，h(t-1)为第t-1时刻当前层的输出，z(t)和r(t)分别为更新门和重置门，

是输入x′(t)和前一时刻输出h(t-1)的汇总，σ为sigmoid激活函数，tanh为双曲正切函数，W_r为重置门训练参数矩阵，W_z为更新门训练参数矩阵，其中[]表示两个向量相连，*表示矩阵的乘积；

步骤4.2，采用反向传播算法BPTT算法对GRU模型进行训练，得到误差最小的参数矩阵即最优GRU模型，误差函数如下：

式中，y(t)为第t时刻实际输出，W_O为输出层参数矩阵；

步骤4.3，重新按照步骤1、步骤2和步骤3方法得到后一段时间锂离子电池子电流集c′₁(t)、c′₂(t)、…、c′_n(t)、残差电流r′_n(t)、电压u′(t)、温度T′(t)作为锂离子电池SOC估计输入集，结合步骤4.2中最优GRU模型，得到模型输出即为SOC估计值。

Images