用于估计电池SOE的装置的制作方法

本发明涉及电池管理技术领域，具体地，涉及用于估计电池soe的装置。

背景技术：

受能源危机与环境危机的影响，电动汽车得到了快速发展。针对电动汽车而言，续航里程估计是一项难题，目前常用锂电池荷电状态(state-of-charge，soc)来对其续航里程进行评判，但是soc与电动汽车续航流程并不呈线性变化，在电池组放电末期，soc的变化要快过可续航里程，这就容易导致用户对电动汽车续航里程的误判。

由于锂电池能量状态(state-of-energy，soe)与电动汽车续航里程密切相关，用soe代替soc对续航里程进行判断是一种更好的方式。当电动汽车行驶状态不变时，其soe与续航里程基本为线性关系，可使用户对续航里程的判断更为准确。

技术实现要素：

本发明的目的是提供一种用于估计电池soe的装置，以此更好地表征电动汽车的续航里程。

为了实现上述目的，本发明的实施方式提供一种用于估计电池soe的装置，该装置包括处理器，该处理器被配置成：建立电池的等效电路模型；根据等效电路模型得到电池系统离散空间模型的状态方程和观测方程；根据接收到的电池两端的检测电压和电池的检测电流、状态方程和观测方程，使用扩展卡尔曼滤波ekf算法来估计电池soe值。

优选地，等效电路模型为戴维南thevenin等效电路模型。

优选地，处理器被配置成：

根据thevenin等效电路模型得到等式(1)：

v(t)＝voc(soe)-v1-rohm*i(t)等式(1)

其中：

v(t)为t时刻的电池的端电压，voc(soe)为t时刻电池的开路电压，v1为电池的极化电压，rohm为电池的欧姆内阻，i(t)为t时刻电池的电流，rpol为电池的极化内阻，cpol为电池的极化电容。

优选地，处理器被配置成：

根据soe的定义得到等式(2)：

其中，soe(t)为t时刻电池的soe值，soe(t0)为电池的初始soe值，η为电池充放电倍率的影响参数，voc(t)为t时刻电池的开路电压，voc(0)为电池满电时的开路电压，c0为电池额定容量；

令voc＝f(soe)，xk＝soek，yk＝vk，uk＝ik，分别将等式(1)和等式(2)离散，得到状态方程和观测方程：

状态方程：

观测方程：

其中，xk为当前状态电池的soe值，xk-1为上一状态电池的soe值，f(xk-1)为上一状态电池的开路电压，uk-1为上一状态电池的观测电流，δt为采样周期，wk为过程噪声，yk为当前状态电池的端电压的观测值，f(xk)为当前状态电池的开路电压，uk为当前状态电池的观测电流，v1(uk)为当前状态电池的极化电压，vk为观测噪声，k＝1,2,3…n，n是正整数。

优选地，处理器被配置成：确定估计次数n和采样周期δt；设定初始状态的误差协方差初始值和初始状态的soe初始值；进行误差协方差时间更新；进行soe状态时间更新；计算卡尔曼增益；进行误差协方差观测更新；根据卡尔曼增益对soe进行状态观测估计更新，以得到估计的soe值。

优选地，该装置还包括：与处理器连接的电压传感器，用于检测电池的电压；与处理器连接的电流传感器，用于检测电池的电流。

优选地，该装置还包括与处理器连接的显示器，显示器被配置成显示电池的soe值。

优选地，处理器包括卡尔曼滤波器。

优选地，上述装置是电池管理系统bms。

上述技术方案基于电池的等效电路模型建立电池系统离散空间模型，并采用扩展卡尔曼滤波(ekf)算法实现电池soe的估计，为用户判断续航里程提供更好的依据。

本发明的其它特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本发明的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：

图1是根据本发明的一实施方式的用于估计电池soe的方法的流程图；

图2是根据本发明的一实施方式的用于估计电池soe的方法的流程图；

图3示出了使用扩展卡尔曼算法估计电池soe值的流程图；

图4是根据本发明的一实施方式的用于估计电池soe的装置的示意框图；

图5示出了电池的等效电路模型thevenin等效电路模型的示意图。具体实施方式

以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。

在本申请中当前状态与第k次估计相对应，上一状态与第k-1次估计相对应。

图1是根据本发明的一实施方式的用于估计电池soe的方法的流程图。如图1所示，在本发明的一实施方式中，提供了一种用于估计电池soe的方法，该方法可以包括：

在步骤s11中，检测电池两端的电压和电池的电流；

在步骤s12中，建立电池的等效电路模型；

在步骤s13中，根据等效电路模型得到电池系统离散空间模型的状态方程和观测方程；

在步骤s14中，根据检测到的电压和电流、状态方程和观测方程，使用ekf算法来估计电池soe值。

电池可以例如是锂电池。

电池的电压可以例如采用电压传感器来检测。电池的电流可以例如采用电流传感器来检测。可替换地，电池的电压和电流可以例如通过电池管理系统(batterymanagementsystem，bms)来检测。

上述的等效电路模型可以例如是thevenin等效电路模型(如图5所示)。该等效电路模型具有物理意义明确、模型参数辨识实验容易执行等优点。上述的等效电路模型还可以是本领域技术人员所知的任意一种电池等效电路模型。

上述的状态方程例如可以根据电池能量状态soe的定义获得。上述的测量方程例如可以根据电池的等效电路模型的电池端电压与电流的关系式获得。

图2是根据本发明的一实施方式的用于估计电池soe的方法的流程图。如图2所示，在本发明的一实施方式中，提供了一种用于估计电池soe的方法，该方法可以包括：

在步骤s21中，检测电池两端的电压和电池的电流；

在步骤s22中，建立电池的等效电路模型；

在步骤s23中，根据等效电路模型得到电池的端电压与soe的关系式；

在步骤s24中，将电池的端电压与soe的关系式离散，得到电池系统离散空间模型的观测方程；

在步骤s25中，将soe的定义式离散，得到电池系统离散空间模型的状态方程；

在步骤s26中，根据检测到的电压和电流、状态方程和观测方程，使用扩展卡尔曼滤波ekf算法来估计电池soe值。

在参考图2所示的实施方式中，步骤s21、s22和s26可以分别与参考图1所示的实施方式中的步骤s11、s12和s14相同或相似。

在本发明的一实施方式中，等效电路模型采用thevenin等效电路模型(如图5所示)，thevenin等效电路模型可以由等式(5)来表达：

v(t)＝voc(soc)-v1-rohm*i(t)等式(5)

其中

v(t)为t时刻的电池的端电压，voc(soc)为t时刻电池的开路电压，表示成soc的函数，v1为电池的极化电压，rohm为电池的欧姆内阻，i(t)为t时刻电池的电流，rpol为电池的极化内阻，cpol为电池的极化电容；

电池的soe和电池的soc的关系可以由等式(6)表达：

其中，soc(t)为t时刻的电池的soc值，ct为t时刻的电池的剩余容量，c0为电池的额定容量，soe(t)为t时刻电池的soe值，e(t)为t时刻的电池的剩余能量，e0为电池的额定能量，voc(t)为t时刻电池的开路电压，voc(0)为电池满电时的开路电压。

由于电池的soc与电池的开路电压voc(t)之间成单调函数关系，同时根据等式(6)，可以得出结论：电池的soe与电池的开路电压voc(t)之间也存在单调函数关系。因此等式(5)可以转化为上述的等式(1)：

v(t)＝voc(soe)-v1-rohm*i(t)等式(1)

其中voc(soe)是t时刻电池的开路电压，表示为soe的函数。

电池的soe与电池的开路电压voc(t)可以由等式(7)来表示：

voc＝f(soe)等式(7)

等式(7)所表达的电池的soe与电池的开路电压voc(t)之间的对应关系(函数关系)可以通过查表获得，也可以通过对电池进行多次充放电实验，记录不同时刻电池的soe值和电池的开路电压值，进而拟合出电池的soe与电池的开路电压voc(t)之间的关系式。

从等式(6)还可以看出，随着电池放电的进行，t时刻电池的开路电压voc(t)与电池满电时的开路电压voc(0)之间的差距越来越大，电池的soe与电池的soc之间并不呈线性关系，因此使用电池的soc来反应电动汽车的续航里程并不准确。

上述的电池的soe的定义式可以由等式(2)来表达：

其中，soe(t)为t时刻电池的soe值，soe(t0)为初始时刻电池的soe值，η为电池充放电倍率的影响参数，voc(t)为t时刻电池的开路电压，voc(0)为电池满电时的开路电压，c0为电池的额定容量。

在本发明的一实施方式中，上述的将电池的端电压与soe的关系式离散得到电池系统离散空间模型的观测方程、将soe的定义式离散得到电池系统离散空间模型的状态方程可以包括：

令xk＝soek，yk＝vk，uk＝ik，分别将等式(1)和等式(2)离散，得到上述的状态方程和观测方程：

状态方程：

观测方程：

其中，xk为当前状态电池的soe值，xk-1为上一状态电池的soe值，f(xk-1)为上一状态电池的开路电压，uk-1为上一状态电池的观测电流，δt为采样周期，wk为过程噪声，yk为当前状态电池的端电压的观测值，f(xk)为当前状态电池的开路电压，uk为当前状态电池的观测电流，v1(uk)为当前状态电池的极化电压，vk为观测噪声，k＝1,2,3…n,n是正整数。

在本发明的一实施方式中，过程噪声wk和测量噪声vk可以例如是高斯白噪声。一般情况下，过程噪声wk和测量噪声vk可以取零值。

如图3所示，在本发明的一实施方式中，上述步骤s14或步骤s26中的根据检测到的电压和电流、状态方程和观测方程，使用ekf算法来估计电池soe值可以包括：

确定优化次数n和采样周期δt；

设定初始状态的误差协方差初始值p0和初始状态的soe初始值x0；

根据上一状态进行误差协方差观测更新后的误差协方差值pk′-1，对当前状态的误差协方差pk进行误差协方差时间更新，k＝1,2,3…n；

根据上一状态进行状态观测估计更新后的电池的soe值x′k-1，对当前状态的电池的soe值xk进行状态时间更新，k＝1,2,3…n；

计算当前状态下电池系统的卡尔曼增益kk；

根据卡尔曼增益kk对误差协方差pk进行误差协方差观测更新，得到观测更新后的误差协方差p′k；

根据卡尔曼增益kk对soe值xk进行状态观测估计更新，以得到估计的soe值x′k。

应该说明的是，第一次对误差协方差进行时间更新时(即k＝1时)，采用初始状态的误差协方差初始值p0作为上一状态进行误差协方差观测更新后的误差协方差值(即p′0＝p0)；第一次对电池的soe值进行状态时间更新时，采用初始状态的soe初始值x0作为上一状态进行状态观测估计更新后的电池的soe值(即x′0＝x0)。

对于上述的优化次数n和采样周期δt，本领域的技术人员应该理解优化次数n的取值越大、采样周期δt的取值越小，最终获得的soe估计值越精确。

上述的初始状态的误差协方差初始值p0和初始状态的soe初始值x0可以是本领域技术人员根据经验和实际需要给定，也可以采用如下方式获取：

对被测电池进行多次充放电实验，记录每次对应的电池soe初值和电池开路电压初值，并建立被测电池的soe初始值的多项式拟合函数；

求解多项式拟合函数，得到soe初始值；

通过梯度仿真实验，确定误差协方差初始值。

通过上述方法可以实现电池soe的估计，为用户判断续航里程提供更加准确的依据。

图4是根据本发明的一实施方式的用于估计电池soe的装置的示意框图。该装置可以用于执行上述实施方式中的用于估计电池soe的方法。具体来说，如图4所示，用于估计电池soe的装置包括处理器10。

在本发明的一个实施方式中，处理器10可以被配置成：

建立电池的等效电路模型；

根据等效电路模型得到电池系统离散空间模型的状态方程和观测方程；

根据接收到的电池两端的检测电压和电池的检测电流、状态方程和观测方程，使用ekf算法来估计电池soe值。

在本发明的一个实施方式中，该装置还可以包括电压传感器20，用于检测电池的端电压；电流传感器30，用于检测电池的电流。

该装置还可以包括显示器40，被配置成显示电池的soe值。

在本发明的另一实施方式中，处理器10可以被配置成：

建立电池的等效电路模型；

根据等效电路模型得到电池的端电压与soe的关系式；

将电池的端电压与soe的关系式离散，得到电池系统离散空间模型的状态方程；

将soe的定义式离散，得到电池系统离散空间模型的观测方程；

根据接收到的电池两端的检测电压和电池的检测电流、状态方程和观测方程，使用ekf算法来估计电池soe值。

在本发明的一个实施方式中，使用ekf算法来估计电池soe值时，处理器10可以被配置成：

确定优化次数n和采样周期δt；

设定初始状态的误差协方差初始值p0和初始状态的soe初始值x0；

根据上一状态进行误差协方差观测更新后的误差协方差值p′k-1，对当前状态的误差协方差pk进行误差协方差时间更新，k＝1,2,3…n；

根据上一状态进行状态观测估计更新后的电池的soe值x′k-1，对当前状态的电池的soe值xk进行状态时间更新，k＝1,2,3…n；

计算当前状态下电池系统的卡尔曼增益kk；

根据卡尔曼增益kk对误差协方差pk进行误差协方差观测更新，得到观测更新后的误差协方差p′k；

根据卡尔曼增益kk对soe值xk进行状态观测估计更新，以得到估计的soe值x′k。

在本发明的一个实施方式中，处理器10可以包括卡尔曼滤波器。

在本发明的一个实施方式中，上述的装置可以例如是电池管理系统bms。

在本发明的一个实施方式中，提供了一种计算机可读存储介质，在该存储介质上存储有指令，该指令被处理器执行时使得该处理器执行如图1、图2和图3中任意所示的用于估计电池soe的方法。

上述方法或装置实现了对电池的soe值的估计，使得电池可以直接采用soe值来表征电池的续航里程，为用户提供更为直接、更为准确的判断依据，防止电池提前终止放电，造成汽车半路“趴窝”等严重后果。

以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。

此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。