本发明涉及到锂电池技术领域，具体涉及到一种锂电池能量状态估算方法。
背景技术：
电池状态一般以电池剩余电量进行表征，可以从容量或者能量的角度进行描述，前者表示剩余容量占额定容量的比值，即电池荷电状态(soc)；后者表示剩余可用能量占总能量的比值，即电池能量状态(soe)。由于soc描述的是电量变化的情况，同时与端电压之间存在较强的耦合关系，采用soc预测与功率直接相关的参数时易于导致较大的预测误差。而soe则与功率参数线性相关，因此在实际应用中，采用soe对电池的续驶里程、续航时间等进行衡量和估算，将是一种较为精确的方法。
但是目前电池能量状态估算精度仍然不高，一方面是计算模型不精确，电池变工况运行时，其电学、热学参数均发生变化，上述变化对电池所带来的损失能量暂时没有得到较为完善的考虑；另一方面，电池初始能量状态的估算也未形成较好的方案，导致目前整体估算精度仍然较低。
技术实现要素：
本发明的目的是提供一种锂电池能量状态估算方法，以解决锂电池能量状态估算不精确的问题。
为实现上述目的，本发明采取的技术方案是：
一种锂电池能量状态估算方法，其包括以下步骤：
步骤1、将锂电池放电过程的能量分为电能、热能以及损失能量三部分，在考虑电能、热能和损失能量的基础上，建立锂电池能量状态的数学估算模型为：
其中，soe(k)为当前时刻电池能量状态，soe(k-1)为上一时刻电池能量状态，tk、tk-1分别为当前时刻和上一时刻时间值，u为电池端电压，i为放电电流且取正值，r为电池内阻，emax为电池最大可用能量，eloss(i,t)为电池的损失能量，esum(i,t)为电池释放的电能和热能的总和；
步骤2、采用电池全过程放电法获得电池的电能和热能参数：以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，放电过程分为两个阶段，第一阶段：电池在非标准放电工况下放电直至达到截止条件，记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至额定温度，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数；
步骤3、利用上述电能和热能的参数，分别得到至少两组电池在放电过程中第一阶段消耗的电能we1(i,t)和热能wh1(i,t)，第二阶段消耗的电能we2(i,t)和热能wh2(i,t)，则
we1(i,t)＝u1×i1×δt(2)
wh1(i,t)＝i12×r1×δt(3)
并得到电池第一阶段释放的能量：
we2(i,t)＝u2×i2×δt(5)
wh2(i,t)＝i22×r2×δt(6)
并得到电池第二阶段释放的能量：
电池的总放电电能为：e(i,t)＝esum(i,t)+e'sum(i,t)
其中，δt＝tk-tk-1，u1和i1分别为第一阶段时电池端电压和放电电流，u2和i2分别为第二阶段时电池端电压和放电电流，电池内阻r为电池欧姆内阻和极化内阻的总和；
步骤4、利用至少两组e(i,t)数据进行函数拟合，并获得拟合函数的最大值e(i,t)max，
从而电池最大可用能量即为：emax＝e(i,t)max(8)
则可知电池在t时刻损失能量为：
eloss(i,t1)＝emax-esum(i,t1)(9)；
步骤5、将式(4)、(8)和(9)代入式(1)中，即可得到t时刻电池的能量状态soe(k)。
实践表明，锂电池在不同放电工况下，其释放的电能和热能并不相同，而电能和热能的总和也不一致，这一现象说明在特定放电工况下电池存在能量损失，导致其最大可用能量无法全部转化为电能和热能。本发明通过在估算模型中引入损失能量的参数，可反映不同放电电流、不同环境温度对能量状态的影响，综合考虑了电能、热能和损失能量对电池能量状态估算的影响，从而获得准确的锂电池能量状态。同时，电池放电采用全过程放电法，并以两阶段的电能和热能拟合得到电池最大可用能量，对传统额定能量参数进行修正，进一步提高了估算的准确度。
soe为100％的状态定义为电池满充状态，也即在标准充电工况下，采用先恒流后恒压的方式对电池进行充电直至充电电流小于0.02c，此时认为soe为100％；soe为0％的状态定义为电池最大可用能量emax全部消耗的状态。
作为本发明的一种改进，电池初始能量状态soe(0)通过电池当前的荷电状态soc(0)得到，其关系式为：
soe(0)＝m×soc(0)+n(10)
其中，m、n为soe和soc曲线的耦合系数。
进一步地，所述电池全过程放电的第一阶段的具体过程为：在额定环境温度下，以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，在不同放电倍率下对电池进行放电，放电直至达到截止条件。
进一步地，所述电池全过程放电的第一阶段的具体过程为：在额定放电倍率下，以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，在不同环境温度下对电池进行放电，放电直至达到截止条件。
作为本发明的一种改进，所述步骤2中采用电池全过程放电法获得电池的电能和热能参数：第一组电池的电能和热能参数获得方式为：以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，放电过程分为两个阶段，第一阶段：在额定环境温度下，在不同放电倍率下对电池进行放电，放电直至达到截止条件，记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至额定温度，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数；
第二组电池的电能和热能参数获得方式为：以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，放电过程分为两个阶段，第一阶段：在额定放电倍率下，在不同环境温度下对电池进行放电，放电直至达到截止条件，记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至额定温度，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数。
与现有技术相比，本发明具有以下优点：
本发明通过在估算模型中引入损失能量的参数，可反映不同放电电流、不同环境温度对能量状态的影响，在综合考虑了电能、热能和损失能量对电池能量状态估算影响的情况下，获得准确的锂电池能量状态；同时，电池放电采用全过程放电法，并以两阶段的电能和热能拟合得到电池最大可用能量，对传统额定能量参数进行修正，进一步提高了估算的准确度。
附图说明
图1为本发明锂电池能量状态估算方法的流程图；
图2为soe和soc的耦合曲线图。
具体实施方式
下面结合附图和具体实施方式对本发明的内容做进一步详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部内容。
本发明通过在估算模型中引入损失能量的参数，可反映不同放电电流、不同环境温度对能量状态的影响，在综合考虑了电能、热能和损失能量对电池能量状态估算影响的情况下，获得准确的锂电池能量状态。
soe为100％的状态定义为电池满充状态，也即在标准充电工况下，采用先恒流后恒压的方式对电池进行充电直至充电电流小于0.02c，此时认为soe为100％；soe为0％的状态定义为电池最大可用能量emax全部消耗的状态。
实施例1
请参考图1，一种锂电池能量状态估算方法，其包括以下步骤：
步骤1、将锂电池放电过程的能量分为电能、热能以及损失能量三部分，在考虑电能、热能和损失能量的基础上，建立锂电池能量状态的数学估算模型为：
其中，soe(k)为当前时刻电池能量状态，soe(k-1)为上一时刻电池能量状态，tk、tk-1分别为当前时刻和上一时刻时间值，u为电池端电压，i为放电电流且取正值，r为电池内阻，emax为电池最大可用能量，eloss(i,t)为电池的损失能量，esum(i,t)为电池释放的电能和热能的总和；
其中，在本实施例中，电池初始能量状态soe(0)通过电池当前的荷电状态soc(0)得到，其关系式为：
soe(0)＝m×soc(0)+n(10)
m、n为soe和soc曲线的耦合系数。
步骤2、采用电池全过程放电法获得电池的电能和热能参数：以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，放电过程分为两个阶段，第一阶段：电池在非标准放电工况下放电直至达到截止条件，记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至额定温度，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数；
在本实施例中，所述电池全过程放电的第一阶段的具体过程为：在额定环境温度下，以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，在不同放电倍率下对电池进行放电，放电直至达到截止条件。
步骤3、利用上述电能和热能的参数，分别得到至少两组电池在放电过程中第一阶段消耗的电能we1(i,t)和热能wh1(i,t)，第而阶段消耗的电能we2(i,t)和热能wh2(i,t)，则
we1(i,t)＝u1×i1×δt(2)
wh1(i,t)＝i12×r1×δt(3)
并得到电池第一阶段释放的能量：
we2(i,t)＝u2×i2×δt(5)
wh2(i,t)＝i22×r2×δt(6)
并得到电池第二阶段释放的能量：
电池的总放电电能为：e(i,t)＝esum(i,t)+e'sum(i,t)
其中，δt＝tk-tk-1，u1和i1分别为第一阶段时电池端电压和放电电流，u2和i2分别为第二阶段时电池端电压和放电电流，电池内阻r为电池欧姆内阻和极化内阻的总和。
步骤4、利用至少两组e(i,t)数据进行函数拟合，并获得拟合函数的最大值e(i,t)max，
从而电池最大可用能量即为：emax＝e(i,t)max(8)
则可知电池在t时刻损失能量为：
eloss(i,t1)＝emax-esum(i,t1)(9)；
步骤5、将式(4)、(8)和(9)代入式(1)中，即可得到t时刻电池的能量状态soe(k)。
以下通过具体例子对上述过程进行讲解：选取一款额定电压/容量为3.2v/20ah的磷酸铁锂电池为实验对象，通过bms对其进行监控，设置其放电截止电压为2v，上限截止温度为60℃，允许其最大放电电流为3c。
为得到电池初始能量状态soe(0)，需要首先确定该锂电池soe与soc的对应关系，其方法为：采用标准充电方法对电池进行充电直至达到满充状态，规定此时电池soe为100％、soc为100％。接着在标准放电工况下(放电倍率为1/3c、环境温度为25℃)，对电池进行放电直至达到放电截止条件。通过对放电过程的容量和能量数据进行百分比处理，以soc作为自变量、soe作为因变量，得到soe和soc的对应曲线如图2所示，拟合得到耦合系数m为1.0374，n为-3.7508。因此，根据当前静置状态下锂电池soc数值，代入初始能量状态计算公式，得到初始能量soe的当前值，该计算结果代表百分比数值。电池初始能量状态由静置状态下的电池荷电状态soc得到，有效提高了电池能量状态soe的估算精度。
接着，以放电倍率1/3c、环境温度25℃作为标准放电工况，采用电池全过程放电法，放电过程分为两个阶段，第一阶段：设置恒温箱环境温度为25℃，以放电截止电压为2v，上限截止温度为60℃作为放电截止条件，分别采用2.5c、2c、1.5c、1c、0.5c、1/3c的放电倍率进行放电，并记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至25℃，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数。
基于上述放电过程得到的实验数据，结合采样时间间隔以及电池内阻r(r取电池欧姆内阻和极化内阻的总和)，根据公式：
we1(i,t)＝u1×i1×δt(2)
wh1(i,t)＝i12×r1×δt(3)
得到电池第一阶段释放的能量：
根据公式：
we2(i,t)＝u2×i2×δt(5)
wh2(i,t)＝i22×r2×δt(6)
得到电池第二阶段释放的能量：
将上述电池第一阶段和第二阶段释放的能量相加即得到电池的总放电电能为：e(i,t)＝esum(i,t)+e'sum(i,t)。利用至少两组e(i,t)数据进行函数拟合，并获得拟合函数的最大值e(i,t)max，从而电池最大可用能量即为：
emax＝e(i,t)max(8)
结合得到的emax以及任一工况下esum(i,t1)的数值，计算得到该工况对应的电池损失能量eloss(i,t1)，则可知电池在t时刻损失能量为：
eloss(i,t1)＝emax-esum(i,t1)(9)
将初始能量状态soe(0)、(4)、(8)和(9)代入式(1)中，在积分时间段段[tk,tk-1]内，通过bms采集得到的电压u、电流i和温度t，对电池能量状态soe进行实时估算，即可得到t时刻电池的能量状态soe(k)。
实施例2
本实施例与实施例1不同之处在于，在步骤2中，所述电池全过程放电的第一阶段的具体过程为：在额定放电倍率下，以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，在不同环境温度下对电池进行放电，放电直至达到截止条件。
以下通过具体例子对上述过程进行讲解：选取一款额定电压/容量为3.2v/20ah的磷酸铁锂电池为实验对象，通过bms对其进行监控，设置其放电截止电压为2v，上限截止温度为60℃，允许其最大放电电流为3c。
为得到电池初始能量状态soe(0)，需要首先确定该锂电池soe与soc的对应关系，其方法为：采用标准充电方法对电池进行充电直至达到满充状态，规定此时电池soe为100％、soc为100％。接着在标准放电工况下(放电倍率为1/3c、环境温度为25℃)，对电池进行放电直至达到放电截止条件。通过对放电过程的容量和能量数据进行百分比处理，以soc作为自变量、soe作为因变量，得到soe和soc的对应曲线如图2所示，拟合得到耦合系数m为1.0374，n为-3.7508。因此，根据当前静置状态下锂电池soc数值，代入初始能量状态计算公式，得到初始能量soe的当前值，该计算结果代表百分比数值。
接着，以放电倍率1/3c、环境温度25℃作为标准放电工况，采用电池全过程放电法，放电过程分为两个阶段，第一阶段：设置电池的放电倍率为3c，以放电截止电压为2v，上限截止温度为60℃作为放电截止条件，分别采用35℃、25℃、15℃、5℃、0℃的环境温度进行放电，并记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至25℃，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数。
基于上述放电过程得到的实验数据，结合采样时间间隔以及电池内阻r(r取电池欧姆内阻和极化内阻的总和)，根据公式：
we1(i,t)＝u1×i1×δt(2)
wh1(i,t)＝i12×r1×δt(3)
得到电池第一阶段释放的能量：
根据公式：
we2(i,t)＝u2×i2×δt(5)
wh2(i,t)＝i22×r2×δt(6)
得到电池第二阶段释放的能量：
将上述电池第一阶段和第二阶段释放的能量相加即得到电池的总放电电能为：e(i,t)＝esum(i,t)+e'sum(i,t)。利用至少两组e(i,t)数据进行函数拟合，并获得拟合函数的最大值e(i,t)max，从而电池最大可用能量即为：
emax＝e(i,t)max(8)
结合得到的emax以及任一工况下esum(i,t1)的数值，计算得到该工况对应的电池损失能量eloss(i,t1)，则可知电池在t时刻损失能量为：
eloss(i,t1)＝emax-esum(i,t1)(9)
将初始能量状态soe(0)、(4)、(8)和(9)代入式(1)中，在积分时间段段[tk,tk-1]内，通过bms采集得到的电压u、电流i和温度t，对电池能量状态soe进行实时估算，即可得到t时刻电池的能量状态soe(k)。
实施例3
其与实施例1不同的是：所述步骤2中采用电池全过程放电法获得电池的电能和热能参数：第一组电池的电能和热能参数获得方式为：以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，放电过程分为两个阶段，第一阶段：在额定环境温度下，在不同放电倍率下对电池进行放电，放电直至达到截止条件，记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至额定温度，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数；
第二组电池的电能和热能参数获得方式为：以放电截止电压a1v，上限截止温度b1℃作为放电截止条件，放电过程分为两个阶段，第一阶段：在额定放电倍率下，在不同环境温度下对电池进行放电，放电直至达到截止条件，记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至额定温度，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数。
两组电池的电能和热能参数分别通过在不同非标准放电工况下获得，然后利用上述两组e(i,t)数据进行函数拟合，并获得拟合函数的最大值e(i,t)max，此时得到的电池最大可用能量比只通过在单一非标准放电工况下拟合得到的电池最大可用能量更加准确。
以下通过具体例子对上述过程进行讲解：选取一款额定电压/容量为3.2v/20ah的磷酸铁锂电池为实验对象，通过bms对其进行监控，设置其放电截止电压为2v，上限截止温度为60℃，允许其最大放电电流为3c。
为得到电池初始能量状态soe(0)，需要首先确定该锂电池soe与soc的对应关系，其方法为：采用标准充电方法对电池进行充电直至达到满充状态，规定此时电池soe为100％、soc为100％。接着在标准放电工况下(放电倍率为1/3c、环境温度为25℃)，对电池进行放电直至达到放电截止条件。通过对放电过程的容量和能量数据进行百分比处理，以soc作为自变量、soe作为因变量，得到soe和soc的对应曲线如图2所示，拟合得到耦合系数m为1.0374，n为-3.7508。因此，根据当前静置状态下锂电池soc数值，代入初始能量状态计算公式，得到初始能量soe的当前值，该计算结果代表百分比数值。
接着，第一组电池的电能和热能参数获得方式为：以放电倍率1/3c、环境温度25℃作为标准放电工况，采用电池全过程放电法，放电过程分为两个阶段，第一阶段：设置恒温箱环境温度为25℃，以放电截止电压为2v，上限截止温度为60℃作为放电截止条件，分别采用2.5c、2c、1.5c、1c、0.5c、1/3c的放电倍率进行放电，并记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至25℃，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数。
第二组电池的电能和热能参数获得方式为：以放电倍率1/3c、环境温度25℃作为标准放电工况，采用电池全过程放电法，放电过程分为两个阶段，第一阶段：设置电池的放电倍率为3c，以放电截止电压为2v，上限截止温度为60℃作为放电截止条件，分别采用35℃、25℃、15℃、5℃、0℃的环境温度进行放电，并记录第一阶段电池的电能和热能的参数；第二阶段：第一阶段结束后，电池静置，待电池温度恢复至25℃，采用标准放电工况对电池进行放电直至达到截止条件，记录第二阶段电池的电能和热能参数。
基于上述放电过程得到的实验数据，结合采样时间间隔以及电池内阻r(r取电池欧姆内阻和极化内阻的总和)，根据公式：
we1(i,t)＝u1×i1×δt(2)
wh1(i,t)＝i12×r1×δt(3)
得到电池第一阶段释放的能量：
根据公式：
we2(i,t)＝u2×i2×δt(5)
wh2(i,t)＝i22×r2×δt(6)
得到电池第二阶段释放的能量：
将上述电池第一阶段和第二阶段释放的能量相加即得到电池的总放电电能为：e(i,t)＝esum(i,t)+e'sum(i,t)。利用两组e(i,t)数据进行函数拟合，并获得拟合函数的最大值e(i,t)max，从而电池最大可用能量即为：
emax＝e(i,t)max(8)
结合得到的emax以及任一工况下esum(i,t1)的数值，计算得到该工况对应的电池损失能量eloss(i,t1)，则可知电池在t时刻损失能量为：
eloss(i,t1)＝emax-esum(i,t1)(9)
将初始能量状态soe(0)、(4)、(8)和(9)代入式(1)中，在积分时间段段[tk,tk-1]内，通过bms采集得到的电压u、电流i和温度t，对电池能量状态soe进行实时估算，即可得到t时刻电池的能量状态soe(k)。
上述实施例仅用以说明本专利而并非限制本专利所描述的技术方案；因此，尽管本说明书参照上述的各个实施例对本专利已进行了详细的说明，但是，本领域的普通技术人员应当理解，仍然可以对本专利进行修改或者等同替换；而一切不脱离本专利的精神和范围的技术方案及其改进，其均应涵盖在本专利的权利要求范围当中。