基于最小二乘法的储能电池荷电和健康状态在线估计方法与流程

1.本发明属于储能电池技术领域，具体涉及一种基于最小二乘法的储能电池荷电和健康状态在线估计方法


背景技术：

2.储能是指通过介质或设备，利用化学或物理的方法把能量存储起来，根据应用需求以特定能量形式释放的过程，通常说的储能是指针对电能的储能。储能技术应用广泛，储能的作用主要是提高电力稳定性和可用性，储存的能量可以用做应急能源，也可以用于在电网负荷低的时候储能，在电网高负荷的时候输出能量，用于削峰填谷，减轻电网波动，同时还可以增强可再生能源利用。
3.随着电力系统、新能源发电(风能、太阳能等)、清洁能源动力汽车等行业的飞速发展，对储能技术尤其大规模储能技术提出了更高的要求。由于风能、太阳能等可再生能源发电具有不连续、不稳定、不可控的特性，可再生能源大规模并入电网会给电网的安全稳定运行带来严重的冲击，而大规模储能系统可有效实现可再生能源发电的调幅调频、平滑输出、跟踪计划发电，从而减小可再生能源发电并网对电网的冲击，提高电网对可再生能源发电的消纳能力，解决弃风、弃光问题。因此，储能技术在电力系统中的广泛应用将成为智能电网发展的一个必然趋势。
4.为了使储能电池与电力系统(配电网)达到更高的适配度，实现对电网电压、功率调整进行更有效的配合，对储能电池实时状态(特别是荷电状态和健康状态)的准确估计就显得尤为重要。只有掌握单体电池或电池组的实时状态，才能做出准确有效的调度决策，防止过充、过放、过温、短路等故障的发生。为此储能系统必须配备相应的电池管理系统，而荷电状态估算作为bms的核心部分，是操作人员控制电池充放电的依据，对于延长电池寿命和保证电池的安全有着重要意义。
5.随着可再生能源的大量并网，其不确定性给传统电网带来了诸多新挑战。储能电池由于其强大的缓冲能力，在电力系统中应用的比例不断提高。然而由于储能电池充放电过程中的出力特性不能准确估计，这给电力系统优化调度增加了不安全因素和不经济因素。调度过程中过充过放都会增加事故发生的可能性，也会减少电池的使用寿命。


技术实现要素：

6.本发明解决的技术问题：提供一种通过电池内部模型参数的构建与辨识来实现荷电状态估计，对多种储能电池均有很好的普适性的基于最小二乘法的储能电池荷电和健康状态在线估计方法。
7.技术方案：为了解决上述技术问题，本发明采用的技术方案如下：
8.一种基于最小二乘法的储能电池荷电和健康状态在线估计方法，包括：首先，根据储能电池内特性将储能电池的内部结构建模为二阶电路模型；采用最小二乘方法，利用历史数据对储能电池内部二阶电路模型中的参数进行辨识，实现在一定误差前提下，获取最
满足端口电路特性的充放电仿真曲线；然后通过建立soc、soh估计公式，基于充放电仿真曲线实现soc和soh状态估计。
9.具体包括以下步骤：
10.s1：收集额定工作条件下储能电池的电流和电压数据；
11.s2：构建储能电池二阶电路内部模型；
12.s3：建立二阶电路的网络拓扑频域方程，经z变换得到差分方程，然后利用最小二乘法的辨识模型，基于充放电数据识别模型参数；
13.s4：根据获得的电池二阶模型进行充放电曲线仿真。
14.s5：soc和soh状态估计
15.作为优选，步骤s1中，历史数据包括充电和放电两个过程的中的端电压电压、端电流电流数据，数据量的大小应满足下列步骤中求解待定系数的要求。
16.作为优选，步骤s3中，二阶电路的网络拓扑频域方程为：
[0017][0018]
把电流i当成系统输入，电压u当成系统输出，利用z变换，把上式变换成差分方程的模式，具体公式为：
[0019]
u(k)＝(1-k
1-k2)u
oc
(k)+k1u(k-1)+k2u(k-2)+k3i(k-1)+k5i(k-2)
[0020]
写成最小二乘的形式，可以得到：
[0021][0022][0023]
θ＝[(1-k
1-k2)u
oc
，k1，k2，k3，k4，k5]
t
[0024]
式中：表示数据向量，表示数据向量的转置，k1，k2，k3，k4，k5为待定系数、θ为待定系数向量，可由与实际数据比对获取；
[0025]
k1，k2，k3，k4，k5和电池模型的参数关系，如下：
[0026][0027][0028][0029][0030][0031]
通过求解得到二阶电路模型的元件参数，式中：t为采样时间。
[0032]
作为优选，步骤s4中，根据已获取的电池内部二阶模型，对电池充放电曲线进行仿真，包括充放电过程中的电压、电流数据。
[0033]
作为优选，步骤s5中，历史研究表明soc与v
oc
有很好的线性相关性，但由于v
oc
的测
量需要一段时间才能达到热力学平衡下的稳定状态，凭借v
oc
估计soc只适合离线状态下的状态估计，不能满足在线估计的要求。因此在本步骤中引入电压变化率的概念v
′
oc
，经数据比对发现soc与1/v
′
oc
也有很好的线性相关性，可建立soc状态估计公式如下：
[0034]
soc(％)＝α
·voc
+β
·
(1/v
′
oc
)+γ
[0035]
其中v
oc
为开路电压；为充放电电压变化仿真曲线的斜率，α，β，γ为待定系数，可由与实际数据比对获取。
[0036]
作为优选，步骤s5中，为满足soh在线状态估计的要求，建立soh状态估计公式如下：
[0037]
soh(％)＝λ(soc)
·
[m
·
(1/v
′
oc
)+n]
[0038]
λ(soc)＝a+b
·
(soc)+c
·
(soc)2+d
·
(soc)3[0039]
其中m，n，a，b，c，d为待定系数，可由与实际数据比对获取，λ为与soc呈三次相关的修正系数，可以看出soh的估计公式并不依赖循环次数的测量。
[0040]
有益效果：与现有技术相比，本发明具有以下优点：
[0041]
本发明的基于最小二乘法的在线估计储能电池荷电状态和健康状态的建模方法，基于最小二乘法对电池内部建立假想二阶内部模型，使用最小二乘法进行参数辨识。根据二阶内部模型进行充放电曲线仿真，利用新颖的估计公式能够更为准确地在线估计电池的荷电状态和健康状态。这将有利于储能电池在电力系统中的并网，更加准确有效地发挥应对不确定性地缓冲作用。与以往的评估方法相比，进一步缩小了仿真模型和实际模型之间的差距，提高了评估精度，并满足在线估计的要求，增强了应用于实际的可能性。经软件模拟预测仿真，和电池实际充放电数据相比较，本模型能够准确的估计储能电池的荷电状态和健康状态。
[0042]
本发明可以较为准确的估计电池充放电过程中电池的荷电状态变化以及健康状态，调度员可以清楚的根据电池当前状态和指令后状态变化下发调度指令，避免过充过放，这将有利于充分发挥储能电池在电力系统中的应用潜力，提升电力系统的经济性和安全性。
附图说明
[0043]
图1是本发明方法的流程框图；
[0044]
图2是储能电池二阶电路内部模型；
[0045]
图3是锂电池恒流充电过程电压变化仿真曲线；
[0046]
图4是锂电池恒流充电过程soc-v
oc
对应关系曲线。
具体实施方式
[0047]
下面结合具体实施例，进一步阐明本发明，实施例在以本发明技术方案为前提下进行实施，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围。
[0048]
本发明提出的基于最小二乘法的储能电池的荷电状态估计方法，首先根据储能电池内特性将储能电池的内部结构建模为含有电容、电阻的二阶电路模型。综合考虑电池充放电时端口电压电流及功率曲线和耦合关系，模拟与实际电路特性相比较为精准的电路结
构。采用最小二乘方法，利用历史数据和对储能电池内部二阶电路模型中的参数进行辨识，实现在一定误差前提下，实现在一定误差前提下，获取最满足端口电路特性的充放电仿真曲线；然后通过建立soc、soh估计公式，基于充放电仿真曲线实现soc和soh状态估计。
[0049]
具体包括以下步骤：
[0050]
步骤1：收集额定工作条件(额定温度、额定负荷等)下目标型号储能电池的历史数据，包括充电和放电两个过程的中的端电压电压、端电流电流数据，数据量的大小应满足下列步骤中求解待定系数的最低要求。
[0051]
步骤2：基于储能电池充放电特性构建储能电池二阶电路内部模型；
[0052]
基于储能电池充放电特性构建含有电阻、电容的电池内部模型，经仿真分析选择建立储能电池二阶电路内部模型，具体如图2所示。r0为电池内部的电阻内阻；r1、r2为极化的电阻；u
oc
为整体开路电压；c1以及c2为极化电容；而i为流过内阻的电流；u为电池端的电压。
[0053]
步骤3：建立二阶电路的网络拓扑频域方程，经z变换得到一差分方程。根据历史数据，利用最小二乘法求取二阶电路模型中的元件参数，使二阶电路模型在给定输入条件下输出与历史数据相同或近似相同，这样最终得到的仿真模型端口特性与实际储能电池相拟合。
[0054]
二阶电路的网络拓扑频域方程为：
[0055][0056]
把电流i当成系统输入，电压u当成系统输出，利用z变换，把上式变换成差分方程的模式，具体公式为：
[0057]
u(k)＝(1-k
1-k2)u
oc
(k)+k1u(k-1)+k2u(k-2)+k3i(k-1)+k5i(k-2)
[0058]
写成最小二乘的形式，可以得到：
[0059][0060][0061]
θ＝[(1-k
1-k2)u
oc
，k1，k2，k5，k4，k5]
t
[0062]
式中：表示数据向量，表示数据向量的转置，k1，k2，k3，k4，k5为待定系数、θ为待定系数向量，可由与实际数据比对获取；
[0063]
k1，k2，k3，k4，k5和电池模型的参数关系，如下：
[0064][0065][0066][0067]
[0068][0069]
通过求解得到二阶电路模型的元件参数。式中：t为采样时间。
[0070]
以容量为2.55ah，充放电截止电压分别为4.2v和2.75v的三元锂电池为例进行实验，经过最小二乘拟合获取二阶电路模型的元件参数。
[0071]
步骤4：根据已获取的电池内部二阶模型，对电池充放电曲线进行仿真，可以获取充放电过程中的电压、电流数据，这里只根据需求获取恒流充电过程电压变化仿真曲线如图3所示。
[0072]
步骤5：进行荷电状态估计和健康状态估计
[0073]
历史研究表明soc与v
oc
有很好的线性相关性，但由于v
oc
的测量需要一段时间才能达到热力学平衡下的稳定状态，凭借v
oc
估计soc只适合离线状态下的状态估计，不能满足在线估计的要求。因此在本步骤中引入电压变化率的概念v
′
oc
，经数据比对发现soc与1/v
′
oc
也有很好的线性相关性。可建立soc状态估计公式如下：
[0074]
soc(％)＝α
·voc
+β
·
(1/v
′
oc
)+γ
[0075]
其中v
oc
为开路电压；为充放电电压变化仿真曲线的斜率，α，β，γ为待定系数，可由与实际数据比对获取：α＝0.9，β＝-0.04，γ＝-2.5。将其代入soc状态估计公式有：
[0076]
soc(％)＝0.9
·voc-0.04
·
(1/v
′
oc
)-2.5。
[0077]
为满足soh在线状态估计的要求，建立soh状态估计公式如下：
[0078]
soh(％)＝λ(soc)
·
[m
·
(1/v
′
oc
)+n]
[0079]
λ(soc)＝a+b
·
(soc)+c
·
(soc)2+d
·
(soc)3[0080]
其中m，n，a，b，c，d为待定系数，可由与实际数据比对获取。λ为与soc呈三次相关的修正系数，可以看出soh的估计公式并不依赖循环次数的测量，满足在线测量得要求。
[0081]
基于图3中的曲线可以获取每一充电时刻v
oc
、v
′
oc
等数据，将多组数据代入soc状态估计公式与soh状态估计公式，与离线实际值比对可以求取m，n，a，b，c，d等参数，本实施例中获取的参数为：m＝0.16，n＝0.04，a＝0.7，b＝-1.5，c＝3.3，d＝-0.9。
[0082]
将其代入soh状态估计公式有：
[0083]
soh(％)＝λ(soc)
·
[0.16
·
(1/v
′
oc
)+0.04]
[0084]
λ(soc)＝0.7-1.5
·
(soc)+3.3
·
(soc)
2-0.09
·
(soc)3[0085]
如以soc为50％的典型状态为例得到soh状态估计公式为：soh(％)＝0.76375
·
[0.16
·
(1/v
′
oc
)+0.041]，结合图3中v
′
oc
的数据可以求取soh实时状态。
[0086]
由soc状态估计公式绘制soc-v
oc
对应关系曲线如图4所示。
[0087]
本发明的基于最小二乘法的储能电池的荷电状态和健康状态估计方法，在构建储能电池二阶电路模型的基础上，采用最小二乘方法，利用历史数据对储能电池内部二阶电路模型中的参数进行辨识，使电路模型与实际模型的充放电特性相一致，获取最满足端口电路特性的充放电仿真曲线。
[0088]
首先需要量测目标储能电池充放电过程中端口的电压、电流数据。然后根据储能电池内部特性建立含电阻、电容的二阶电路模型。基于确定好的二阶电路模型，列写电路频
域拓扑方程，并经z变换转化为差分方程形式。将差分方程化为最小二乘形式，并代入充放电过程中量测好的历史数据，求取最符合充放电特性时的二阶电路元件参数。根据获得的电池二阶模型进行充放电曲线仿真，通过建立依靠于电压变化率的soc、soh
[0089]
状态估计表达式，实现储能电池的在线状态估计。
[0090]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。