一种锂离子电池荷电状态与开路电压模型建立方法与流程

本发明涉及新能源汽车电池应用领域，具体涉及用于锂离子电池荷电状态(soc)估算之前，建立合理的锂离子电池荷电状态soc与开路电压(ocv)之间关系模型。

背景技术：

电池的荷电状态soc反映了电池的剩余电量，是电池的管理系统中是一个非常基本又重要的监测指标。由于电池soc同电池其他参数关系表现出高度的非线性，不得不通过中间的物理量，迂回获得soc；加之外部环境和内部环境参数变化的随机性，使系统基于数学模型的预估方法不够准确，因此必须对电池荷电状态估计的抗干扰能力和自适应能力进行抑制和提高。

目前用于soc估算的方法有：安时计量法、开路电压法、阻抗分析法、卡尔曼滤波法和神经网络法等。无论采用什么方法，都离不开soc-ocv关系曲线模型；因而准确对soc-ocv建模，是提高soc估算精度的重要前提。当前所有soc-ocv之间的模型建立都是在某一定条件下建立的，而实际应用中电池环境条件和电池本身性能都随时间和场景而变化，单一条件下的soc-ocv模型已无法满足后续soc的估算要求。

本发明的方法是给出具有较强的自适应能力的soc-ocv关系模型，其方法计算量小于目前大多方法，且容易实现、有利于估算提高电池soc估算精度。

技术实现要素：

本发明的目的在于提供一种锂离子电池荷电状态与开路电压模型建立方法，解决传统锂离子电池soc估测算法中，soc-ocv模型复杂、计算量大，估算精度低，实现困难的问题。

实现本发明目的的技术方案主要包括以下步骤：

步骤1：在常温下，实测电池soc、ocv数据；

步骤2：将soc重叠分成三段；

步骤3：对三段分别用低价多项式进行拟合

步骤4：对重合部分进行线性插值处理；

步骤5：考虑soh变化，对拟合公式进行修正；

步骤6：考虑温度变化，对拟合公式进行修正；

步骤7：综合考虑soh、温度变化，对拟合公式进行最终修正；

步骤8：建立soc-ocv关系模型。

将soc重叠分成三段，第一段在(0,0.2+δ1)范围、第二段在(0.2-δ1，0.75+δ2)范围、第三段在(0.75-δ2,1)。

分别对三段拟合，拟合多项式用低价的(4阶)。

将二处重合部分，在低价多项式拟合基础上进行线性插值。

考虑soh变化影响时，将拟合公式进行soc用soc/soh修正。

考虑温度变化影响时，将拟合公式用soce^0.002(t-25)代替soc修正。

当同时考虑soh和温度变化影响时，将拟合公式中自变量用(soce^0.002(t-25)/soh)修正。

仅需在常温下，实测电池soc、ocv数据，即可建立同时考虑soh和温度影响的模型。

本发明与现有技术相比，其显著优点为：(1)将soc重叠分段再进行低价多项式拟合，其阶次低、计算量小，而整体精度高；(2)拟合模型考虑到电池健康度的变化对模型的影响。同时也考虑到温度变化影响，拟合结果没有增加模型复杂程度；(3)本发明具有较好的通用性，不仅仅局限于某一种电池；(4)本发明具有较好的适应性，适合于电池长期应用后环境变化和电池特性下降情况。

下面结合附图对本发明作进一步详细描述。

附图说明

图1是本发明的数据处理流程示意图

图2是不同温度不同soh下实测soc-ocv数据示意图

图3是本方法拟合后示意图

具体实施方式

下面结合附图对本发明算法做进一步的说明。

如图1所示，本发明的一种锂离子电池荷电状态与开路电压模型建立方法，用于提高锂离子电池荷电状态soc估测精度的能力。方法主要包括以下过程和处理手段，首先将soc-ocv分段处理，降低拟合阶次、提高拟合精度；在分段重叠区，用线性插值方法，保证拟合的平滑性；其次考虑电池健康度soh和温度t对soc-ocv模型的影响，使soc-ocv模型更符合实际运行规律。本发明给出的soc-ocv建模方法，解决soc估测过程中基础模型不匹配问题，有利于电池后续管理工作。该方法可用于各种电池soc估测方法。

结合图1，说明soc-ocv关系建立的详细步骤。

步骤一，在常温下，实测电池soc、ocv数据。

一般电池在常温(25^o)下恒流恒压充电下，可以获取电池开路电压(ocv)与电池soc的变化变化关系见图2中soh＝1的曲线。

实验过程分为以下几步：实现以将电池放电至其下限截止电压，放电电流为1/2c并静置3h以上；其次以充电电流为1/2c充入电池标称容量5％的电量，并静置2h以上，记录静置后的电池端电压作为当前soc处的ocv；重复以上步骤2，直至充电过程中电池的电压达到其上限截止电压，实验结束。

步骤二，根据步骤一的数据，将soc重叠分成三段，第一段在(0,0.2+δ1)范围、第二段在(0.2-δ1，0.75+δ2)范围、第三段在(0.75-δ2,1)。

由于电池开路电压与荷电状态之间的关系高度非线性，目前多采用六阶或以上多项式函数来拟合ocv与soc关系模型，或者用更为复杂的函数来拟合建模。在不降低拟合精度下，为了降低拟合的计算复杂度和运算量；将soc从0到1分成重叠的三段：(0,0.2+δ1)、(0.2-δ1，0.75+δ2)、(0.75-δ2,1)。

步骤三，分别对以上三段进行拟合，拟合多项式用低价的(4阶)。

每一段可以用四阶或更低的多项式函数拟合如式(1)、(2)和(3)，且达到更高的拟合精度，每一段拟合仍采用常规的最小二乘法；其中δ1取0.03附近，其中δ2取0.05附近。

结合图1，soc在第一段(0,0.2+δ1)的拟合关系如下：

ocv1(soc)＝a11soc⁴+a12soc³+a13soc²+a14soc+a15(1)

结合图1，soc在第二段(0.2-δ1，0.75+δ2)的拟合关系如下：

ocv2(soc)＝a21soc⁴+a22soc³+a23soc²+a24soc+a25(2)

结合图1，soc在第三段(0.75-δ2,1)的拟合关系如下：

ocv3(soc)＝a31soc⁴+a32soc³+a33soc²+a34soc+a35(3)

步骤四，将第一段与第二段、第二段与第三段二处重合部分，在低价多项式拟合基础上进行线性插值。

结合图1，由于分段拟合，三段之间不易保证平滑连续；下面给出一种简单实用的处理方法。第一段与第二段的soc在重叠区(0.2-δ1，0.2+δ1)模型为ocv12(soc)，由式(1)和式(2)进行线性差值：

第二段与第三段的soc在重叠区(0.75-δ2，0.75+δ2)模型为ocv23(soc)，由式(2)和式(3)进行线性插值处理：

结合图1，通过上面处理，实质上将soc从0到1分成的五段：(0,0.2-δ1)、(0.2-δ1，0.2+δ1)、(0.2+δ1，0.75-δ2)、(0.75-δ2，0.75+δ2)和(0.75+δ2,1),它们使用的模型分别是式(1)、(2)、(3)、(4)和(5)。

步骤五，考虑soh变化影响时，将拟合公式进行soc用soc/soh修正。

结合图1，以及图2中soh＝0.8曲线的可知，当电池经过多次循环工作后其健康度soh值要下降，其ocv与soc关系已经发生改变。若仍然按上面的拟合模型进行soc相关估算，必然导致很大的误差。从图2中soh＝0.8曲线，可见随着电池健康状态(soh)值下降，ocv与soc关系在横坐标方向往原点收缩；因而可以将上面拟合的公式ocvi(soc)，进一步修改为下式：

其中i＝1,2,3

步骤六，考虑温度变化影响时，将拟合公式用soce^0.002(t-25)代替soc修正。

结合图1，以及其次图2中soh＝1和soh＝0.8曲线(虚线和点划线)的可知，在运行或环境温度发生变化时，电池的温度(t)也相应发生改变，若仍然按上面的拟合式进行soc相关估算，也必然导致很大的误差。从图2中soh＝1和soh＝0.8曲线(虚线和点划线)，可见随着温度下降，ocv与soc关系在纵向坐标方向往下偏移，反之亦然；因而可以将拟合的公式ocvi(soc)，最终修改为下式：

ocvi⁺(soc)＝ai1(soce^0.002(t-25))⁴+ai2(soce^0.002(t-25))³+ai3(soce^0.002(t-25))²+ai4(soce^0.002(t-25))+ai5

其中i＝1,2,3

步骤七，仅需在常温下，实测电池soc、ocv数据，即可建立同时考虑soh和温度影响的模型；最终分成的五段：(0,0.2-δ1)、(0.2-δ1，0.2+δ1)、(0.2+δ1，0.75-δ2)、(0.75-δ2，0.75+δ2)和(0.75+δ2,1)。

重叠段(0.2-δ1，0.2+δ1)和(0.75-δ2，0.75+δ2)模型分别如下：

其他段(0,0.2-δ1)、(0.2+δ1，0.75-δ2)和(0.75+δ2,1)模型分别如下：

其中i＝1,2,3

为了验证该方法的有效性，图3给出按上述设计过程的拟合结果。实验结果表明这种模型能适应温度和soh的变化。虽然是分段处理，由于采用重叠部分线性插值，仍然保持良好的平滑性，有利于后续工作。

显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无法对所有的实施方式予以穷举。凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。