基于脉冲放电压差的储能电池模型参数辨识方法与流程

本发明涉及电气工程领域，具体涉及一种基于脉冲放电压差的储能电池模型参数辨识方法。

背景技术：

大规模铅炭电池储能系统是由单体电池串并联构成，各单体电池伏安特性及运行工况的差异性，导致铅炭电池储能系统因个别电池首先达到运行极限而制约其他单体电池调节能力的利用，故储能系统是由“短板电池”制约的受限调节电源。因此，设计一款更加贴近真实电池充放电效果的电池充放电模型以及相关的精确计算方法，对于蓄电池的合理利用以及相关数学模型的建立具有重大意义。

传统计算方法是基于充放电曲线的纯数学计算，计算过程涉及高阶方程等超越方程，需要软件仿真与计算，手算参数难以实现。且传统计算方法没有考虑电路模型的基本电路知识，不利于对电池模型本身的理解。

技术实现要素：

为解决上述问题，本发明提供了一种基于脉冲放电压差的储能电池模型参数辨识方法，利用了基本电路特性，避免了参数方程的拟合，简化了计算过程，使铅碳电池充放电模型的外特性得到更好地理解。使用该方法计算铅碳电池充放电模型阻容参数的计算量小、过程简单，计算结果精确。

为实现上述目的，本发明采取的技术方案为：

基于脉冲放电压差的储能电池模型参数辨识方法，包括如下步骤：

s1、ro辨识：

基于铅炭电池等效电路模型，放电开始瞬间，电流i在内阻ro上产生电压降，形成电压曲线ab段，此时电容电压up等于零，因此，内阻ro可由ab段电压差除以流过ro的电流得到，用ua表示电池放电前的开路电压，ub表示放电瞬间电流流过ro时的电池端电压，可以得到ro的求解表达式如下：

s2、rp辨识：

恒流放电达到稳态时，流过电容cp的电流等于零，即电流i全部通过电阻ro，电容cp的电压等于电阻rp两端的电压，即

up＝i*rp(6)

放电结束瞬间，电流i等于零，电阻ro上的电压等于零，根据ul＝uoc-ir0-up可知，电池端电压u1的表达式为：

ul＝uoc-up(7)

用ud表示放电结束瞬间，电流i等于零时的电池端电压，根据公式(7)，结合前述假设uoc不变，可得电容电压求解表达式如下：

up＝ua-ud(8)

根据公式(6)和(8)可求得rp的表达式如下：

s3、cp辨识：

放电过程：

a点之前，电池处于静止状态，cp两端电压等于零，b点电池开始放电，可等效为零状态响应，设b点为零状态响应起点，则up等于：

其中，τdis＝rp*cpdis；

将式(10)代入ul＝uoc-ir0-up，可得电池端电压u1的表达式为：

用uc表示放电终止时刻的电池端电压，用tbc表示b、c两点的时间间隔，则式(11)可表示为：

由式(12)可求得电容cpdis的表达式为：

静止过程：

放电结束后，电容cp存储的能量通过rp释放，可等效为零输入响应，设d点为零输入响应的起点，则up等于：

其中，τsta＝rp*cpsta；

将式(10)代入ul＝uoc-ir0-up，可得电池端电压u1的表达式为：

ul＝uoc-ude^-t/τ(15)

用ue表示电池在e点的端电压，用tde表示d、e两点的时间间隔，根据式(15)可求得电容cpsta的表达式为：

根据公式(5)、(9)、(13)、(16)和放电曲线，求得的铅炭电池模型参数：ro；rp；cpdis；cpsta。

本发明具有以下有益效果：

该计算方法对于电池模型辨识以及电池快速soc计算具有重要作用，相较传统计算方法，该模型辨识方法计算量小，计算过程结合电池物理模型，使计算具有更加灵活的环境适应性以及计算快速性。

附图说明

图1为thevenin模型。

图2为铅炭电池脉冲放电曲线。

图3为铅炭电池simulink仿真模型。

图4为thevenin模型输出电压与实际电压对比图。

图5为仿真结果与实际电压差值。

具体实施方式

为了使本发明的目的及优点更加清楚明白，以下结合实施例对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

本发明的基于脉冲放电压差的储能电池模型参数辨识方法是基于已被提出并广泛应用的电池充放电模型thevenin模型提出的。thevenin模型由理想电压源e、并联rc网络和串联电阻ro组成，如图1所示。理想电压源为电池开路电压，与soc有关。并联rc网络用于表示电池的极化反应过程。串联电阻ro用于表示电池极板、电解液和间隔板的电阻。thevenin模型可以很好地模拟电池静态和暂态特性，且模型精度和计算复杂度均比较适中，是模拟电池特性较好的选择。

基于thevinin等效电路的铅炭电池模型参数与电池内部结构、电解液浓度以及电池运行工况有关，无法通过测量直接获取，需要根据电池外部特性，通过特定的实验和数据处理方法得到。电池模型需要辨识的参数有开路电压uoc、欧姆内阻ro、极化电阻rp和极化电容cp。

根据基尔霍夫电压、电流定律和电容电压与电流的关系，可以得到thevenin模型的数学表达式如下：

ul＝uoc-ir0-up(1)

根据soc定义，采用thevenin模型时，电池soc的表达式如下：

式中：ce为电池额定容量，单位为ah，η为充放电效率，i为电池电流，以放电方向为正。

开路电压uoc可表示为soc的函数，即uoc＝f(ssoc)。

对铅炭电池进行脉冲放电，得到放电曲线如图2所示，红色曲线为电流曲线，蓝色曲线为电压曲线，因为放电时间很短，可以认为电池soc不变，从而可以推出与soc对应的电池开路电压不变，即放电前和放电后电池开路电压uoc不变。

基于脉冲放电压差的储能电池模型参数辨识方法，包括如下步骤：

(1)ro辨识：

结合图1所示的铅炭电池等效电路模型，放电开始瞬间，电流i在内阻ro上产生电压降，形成电压曲线ab段，此时电容电压up等于零，因此，内阻ro可由ab段电压差除以流过ro的电流得到，用ua表示电池放电前的开路电压，ub表示放电瞬间电流流过ro时的电池端电压，可以得到ro的求解表达式如下：

(2)rp辨识：

恒流放电达到稳态时，流过电容cp的电流等于零，即电流i全部通过电阻ro，电容cp的电压等于电阻rp两端的电压，即

up＝i*rp(6)

放电结束瞬间，电流i等于零，电阻ro上的电压等于零，根据公式(1)可知，电池端电压u1的表达式为：

ul＝uoc-up(7)

用ud表示放电结束瞬间，电流i等于零时的电池端电压，根据公式(7)，结合前述假设uoc不变，可得电容电压求解表达式如下：

up＝ua-ud(8)

根据公式(6)和(8)可求得rp的表达式如下：

(3)cp辨识：

观察电池端电压，bc和de两段电压曲线的下降和上升斜率差别较大，用一个cp不能准确描述电池特性，因此需要根据电池运行条件分别进行cp辨识。

放电过程：

a点之前，电池处于静止状态，cp两端电压等于零，b点电池开始放电，可等效为零状态响应，设b点为零状态响应起点，则up等于：

其中，τdis＝rp*cpdis。

将式(10)代入公式(1)，可得电池端电压u1的表达式为：

用uc表示放电终止时刻的电池端电压，用tbc表示b、c两点的时间间隔，则式(11)可表示为：

由式(12)可求得电容cpdis的表达式为：

静止过程：

放电结束后，电容cp存储的能量通过rp释放，可等效为零输入响应，设d点为零输入响应的起点，则up等于：

其中，τsta＝rp*cpsta。

将式(10)代入公式(1)，可得电池端电压u1的表达式为：

ul＝uoc-ude^-t/τ(15)

用ue表示电池在e点的端电压，用tde表示d、e两点的时间间隔，根据式(15)可求得电容cpsta的表达式为：

根据公式(5)、(9)、(13)、(16)和图2的放电曲线，求得的铅炭电池模型参数：ro；rp；cpdis；cpsta。

在辨识得到的电池模型参数基础上，将电流和soc作为电池模型的输入量，实时输出等效电路模型的端电压，对电池模型输出的电压和实际电池电压比较，进行误差分析，验证等效电路模型的合理性和参数辨识的准确性。

模型应用之前，需要对其离散化，对thevenin等效电路模型离散化得到以下表达式：

其中，u1(k)、uoc(k)、up(k)和i(k)分别为k时刻电池端电压、开路电压、极化电压和电池工作电流，up(k+1)为k+1时刻电池极化电压，t为采样周期，τ＝rpcp，为时间常数。

在simulink仿真平台搭建的铅炭电池模型如图3所示，其中，模型输入为电池实际端电压、电池实际工作电流、电池初始soc值、电流效率、电池额定容量、采样周期、极化电阻、电池内阻、放电时的极化电容、静止时的极化电容和k时刻的电容电压，模型输出为电池端电压、通过模型得到的电池端电压与实际端电压的差值和k+1时刻的电容电压。

室温下，以0.1c的电流幅值对铅炭电池进行脉冲放电，记录电池电流和端电压曲线。将电池电流曲线作为thevenin电池模型输入，记录模型输出电压曲线，并与实际电池端电压曲线进行对比，如图4所示，图中，红色曲线为实际电池端电压曲线，蓝色曲线为电池模型输出曲线。

图5为仿真结果与实际电压的差值，可以看出，在放电开始和放电结束瞬间，仿真结果与实际电压差值较大，其他时间电压差值均很小，计算得到的平均误差为0.6mv。

模型仿真结果与实际电池电压曲线一致，除了放电开始和结束瞬间，仿真误差小于1.5mv，平均误差为0.6mv，能准确地描述电池原型，验证了模型的合理性和参数的准确性。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。