CN104951662B - 一种磷酸铁锂电池能量状态soe的估算方法 - Google Patents
一种磷酸铁锂电池能量状态soe的估算方法 Download PDFInfo
- Publication number
- CN104951662B CN104951662B CN201510418388.1A CN201510418388A CN104951662B CN 104951662 B CN104951662 B CN 104951662B CN 201510418388 A CN201510418388 A CN 201510418388A CN 104951662 B CN104951662 B CN 104951662B
- Authority
- CN
- China
- Prior art keywords
- msub
- energy
- mrow
- battery
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Battery Electrode And Active Subsutance (AREA)
- Secondary Cells (AREA)
Abstract
本发明提供了一种磷酸铁锂电池能量状态(State of Energy,SOE)的估算方法,即根据磷酸铁锂电池内部能量组成形式,建立包含电能和热能的SOE数学估算模型。通过在不同放电倍率下对电池放电过程的外部消耗电能以及电池内部欧姆热能、极化热能和熵产热能进行测算,获取各种放电倍率下最大可用能量,并模拟得到最大理论总能和各种放电倍率下的效率函数关系式,对估算模型中的总能进行实时修正,从而有效提高电池能量状态SOE的估算精度。本发明引入热能形式完善了电池放电过程中能量消耗的种类,同时以最大理论总能及能量释放效率对估算过程的释放总能进行实时调整,使模型更贴近实际工况,具有物理概念清晰,估算准确度高等优点。
Description
技术领域
本发明属于磷酸铁锂电池技术领域,特别涉及一种磷酸铁锂电池能量状态SOE(State Of Energy)的估算方法。
背景技术
雾霾现象的发生,特别是PM2.5粒子引发的健康问题与新能源应用中所固有的“间歇性、随机性”等特性问题进一步推动了电池储能系统取代化石能源在汽车上的应用,随着储能设备的大规模应用,单一电池荷电状态参数不能满足现实要求,大量有关于电池能量状态估算的研究兴起。
电池剩余能量的估计一般是通过两种方式,即电池容量状态SOC(State ofCharge,即荷电状态,也叫剩余电量)间接表示和通过电池的SOC模型与相对应的电压模型相互耦合得到剩余能量,由于锂离子电池的大规模应用,储电设备必须配备更为复杂且精准的电池管理系统,单纯依靠SOC并不能线性表征电池剩余能量,而现实运行工况的多样性、两模型耦合的复杂度以及SOC所表征的电池能量非线性均限制了耦合能量表征方式的进一步发展。
本发明首先通过研究电池充放电过程中的能量转换关系,引入热能完善能量消耗形式,同时提出电池最大允许能量和能量释放效率概念,建立该估算方法的数学模型,更有利于直接表征电池不同运行工况下的剩余能量,有效预测电池剩余运行时间。
发明内容
本发明的目的在于提供一种磷酸铁锂电池能量状态SOE的估算方法,以解决当前电池能量状态SOE估算不精确的问题。
为了实现上述目的,本发明通过下列技术方案来实现:
一种磷酸铁锂电池能量状态SOE的估算方法,其包括以下步骤:
步骤1、建立磷酸铁锂电池包含电能和热能两种能量损耗形式的SOE数学估算模型,所述电能为电池对外提供的能量总和,热能为电池内部产生的欧姆热能、极化热能和熵产热能的总和;所述SOE数学估算模型为
其中,SOE0为电池初始能量状态,U为电池端电压,I为电流,R为电阻,T为温度,ΔS为电池反应熵变,Emax为电池最大理论总能,η(T,I)为电池能量释放效率;
步骤2、以电池的截止电压a1V,安全温度b1℃作为其放电截止的条件,在不同放电倍率下对电池进行放电,由测算得到的电能和热能计算获取最大可用能量Ex,Ex为在放电倍率为x时测量得到的电池外部电能和内部热能之和;
步骤3、根据不同放电倍率下的最大可用能量Ex拟合得到最大理论总能Emax;
步骤4、通过不同放电倍率下最大可用能量Ex与最大理论总能Emax的比值获得电池能量释放效率的函数关系式:
其中,c1-c6为参量,通过联立多个放电倍率x与其对应的最大可用能量Ex以及最大理论总能Emax之间的关系进行求解;
步骤5、将求解后的参量c1-c6代入式(2)中,并将式(2)获得的关于电池能量释放效率的函数关系式以及步骤3中得到的最大理论总能Emax代入式(1)中,作为修正因子对SOE数学估算模型在积分段t0-t1内电池能量状态进行实时修正:
在式(1)中,电阻R=R0+Rp,其中,R0为欧姆电阻,Rp为极化电阻;在电池不同荷电状态下通过采用若干放电倍率对电池在一定脉冲电流下进行脉冲放电,脉冲放电时刻,电池电压会有一个急速下降和缓慢下降的阶段,其中,急速下降阶段的电压变化值与该脉冲电流的比值为欧姆电阻R0的值,缓慢下降阶段的电压变化值与该脉冲电流的比值为极化电阻Rp的值。
所述公式(1)中,电池反应熵变其中,n为电池交换电子数,为常数1,F为法拉第常数,数据通过将不同荷电状态的电池静置放置在依次从0℃—60℃变化温度通过电压的变化梯度下求得。
在上述的一种磷酸铁锂电池能量状态SOE的估算方法中,所述的数学模型中对最大理论总能Emax的确定,以磷酸铁锂电池的截止电压a1V,安全温度b1℃作为电池放电截止的条件,选取若干放电倍率电流对磷酸铁锂电池进行放电直至截止条件,计算出该过程消耗的电能与热能,电能与热能之和即为最大可用能量,拟合若干放电倍率下总能得到一条横坐标为放电倍率,纵坐标为最大可用能量的曲线,该曲线的最大值可认为是电池的最大理论总能Emax。
在上述的一种磷酸铁锂电池能量状态SOE的估算方法中,所述的数学模型中对能量释放效率的确定,电池最大理论总能是一定的,但在不同放电倍率下电池所能放出的最大可用能量是不相同的,这其中就存在能量释放效率问题,此时计算若干放电倍率下电池释放的最大可用能量与最大理论总能Emax的比值,将得到的比值拟合,得到一条横坐标为放电倍率,纵坐标为能量释放效率曲线,该曲线即为能量释放效率的函数关系式。
本发明与现有模型相比,具有如下优点:添加了电池放电过程中的热量形式,完善了能量损耗类型,同时提出了最大理论总能Emax及能量释放效率两个概念,并加入到SOE数学模型中去,提高了现有数学模型的估算精度。
附图说明
图1是本磷酸铁锂电池能量状态SOE的估算方法结构图;
图2是本磷酸铁锂电池热能随放电倍率变化的关系图;
图3是本磷酸铁锂电池电能随放电倍率变化的关系图;
图4是本磷酸铁锂电池总能随放电倍率变化的关系图。
具体实施方式
下面结合附图和具体实施方式对本发明的内容做进一步详细说明。
如图1所示,本磷酸铁锂电池能量状态SOE的估算方法,包括以下步骤:
(1)建立包含电能和热能两种能量损耗形式的SOE数学估算模型,电能为电池对外提供的能量总和,热能为电池内部产生的欧姆热能、极化热能和熵产热能的总和;
(2)以磷酸铁锂电池的截止电压a1V,安全温度b1℃作为电池放电截止的条件,在不同倍率下对电池进行放电,由测算得到的电能和热能得到各个放电倍率下的最大可用能量;
(3)模拟得到电池最大理论总能和各放电倍率下的效率函数关系式,作为修正因子对估算过程每个积分段内的总能参数进行实时修正,从而提高SOE估算的准确度。
首先,选取了一款额定电压、容量为3.65V、20Ah的磷酸铁锂电池作为实验对象,其放电截止电压为2V,工作温度为-20℃—60℃,允许最大放电电流为3C。通过以上参数可知,其放电截止条件为截止电压到达2V或达到温度上限60℃。
然后,选取0.1C—3C中10个不同的放电倍率针对满电状态下的电池进行放电,直至达到截止条件。分别计算10个放电倍率下电池的热能与电能消耗。热能和电能随放电倍率变化的数据如图2、3所示。
其次,将上述不同放电倍率下电能与热能相加,即得到不同倍率放电到截止条件的最大可用能量,将最大可用能量拟合,得到当倍率为2.6C时其能量最高点为321562J,如图4所示,则321562J为该电池的最大理论总能Emax。
最后,计算上述10个不同放电倍率下释放的最大可用能量与最大理论总能321562J相比,得到10个对应的数值,再将其进行拟合,得到一条随放电倍率变化的能量释放效率曲线,且通过拟合,得到曲线关系式为:
η=-0.0465x5+0.3145x4-0.7591x3+0.8012x2-0.2693x+0.7815 (4)
其中,η为电池能量释放效率,x为电池放电倍率。
需要说明的是:能量释放效率曲线关系式中参量c1-c6可以采用上述拟合曲线方式获得,也可以通过联立六个以上的方程(即取6个不同的放电倍率,其能量释放效率则可以确定)。最终将拟合或求解得到的电池能量释放效率曲线关系式(4)代入SOE数学估算模型中进行修正。
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所做出的等效的变化或修饰,都应涵盖在本发明的保护范围内。
Claims (3)
1.一种磷酸铁锂电池能量状态SOE的估算方法,其特征在于,其包括以下步骤:
步骤1、建立磷酸铁锂电池包含电能和热能两种能量损耗形式的SOE数学估算模型,所述电能为电池对外提供的能量总和,热能为电池内部产生的欧姆热能、极化热能和熵产热能的总和;所述SOE数学估算模型为
<mrow>
<mi>S</mi>
<mi>O</mi>
<mi>E</mi>
<mo>=</mo>
<msub>
<mi>SOE</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mi>U</mi>
<mi>I</mi>
<mo>+</mo>
<msup>
<mi>I</mi>
<mn>2</mn>
</msup>
<mi>R</mi>
<mo>+</mo>
<mi>T</mi>
<mi>&Delta;</mi>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>,</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,SOE0为电池初始能量状态,U为电池端电压,I为电流,R为电阻,T为温度,ΔS为电池反应熵变,Emax为电池最大理论总能,η(T,I)为电池能量释放效率;
步骤2、以电池的截止电压a1V,安全温度b1℃作为其放电截止的条件,在不同放电倍率下对电池进行放电,由测算得到的电能和热能计算获取最大可用能量Ex,Ex为在放电倍率为x时测量得到的电池外部电能和内部热能之和;
步骤3、根据不同放电倍率下的最大可用能量Ex拟合得到最大理论总能Emax;
步骤4、通过不同放电倍率下最大可用能量Ex与最大理论总能Emax的比值获得电池能量释放效率的函数关系式:
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>E</mi>
<mi>x</mi>
</msub>
<msub>
<mi>E</mi>
<mi>max</mi>
</msub>
</mfrac>
<mo>=</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msup>
<mi>x</mi>
<mn>5</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msup>
<mi>x</mi>
<mn>4</mn>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msup>
<mi>x</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>5</mn>
</msub>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,c1-c6为参量,通过联立多个放电倍率x与其对应的最大可用能量Ex以及最大理论总能Emax之间的关系进行求解;
步骤5、将求解后的参量c1-c6代入式(2)中,并将式(2)获得的关于电池能量释放效率的函数关系式以及步骤3中得到的最大理论总能Emax代入式(1)中,作为修正因子对SOE数学估算模型在积分段t0-t1内电池能量状态进行实时修正:
<mrow>
<mi>S</mi>
<mi>O</mi>
<mi>E</mi>
<mo>=</mo>
<msub>
<mi>SOE</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<mi>U</mi>
<mi>I</mi>
<mo>+</mo>
<msup>
<mi>I</mi>
<mn>2</mn>
</msup>
<mi>R</mi>
<mo>+</mo>
<mi>T</mi>
<mi>&Delta;</mi>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msup>
<mi>x</mi>
<mn>5</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msup>
<mi>x</mi>
<mn>4</mn>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msup>
<mi>x</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>5</mn>
</msub>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>6</mn>
</msub>
<mo>)</mo>
<mo>*</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2.根据权利要求1所述的磷酸铁锂电池能量状态SOE的估算方法,其特征在于,在式(1)中,电阻R=R0+Rp,其中,R0为欧姆电阻,Rp为极化电阻;在电池不同荷电状态下通过采用若干放电倍率对电池在一定脉冲电流下进行脉冲放电,脉冲放电时刻,电池电压会有一个急速下降和缓慢下降的阶段,其中,急速下降阶段的电压变化值与该脉冲电流的比值为欧姆电阻R0的值,缓慢下降阶段的电压变化值与该脉冲电流的比值为极化电阻Rp的值。
3.根据权利要求1或2所述的磷酸铁锂电池能量状态SOE的估算方法,其特征在于,所述公式(1)中,电池反应熵变其中,n为电池交换电子数,为常数1,F为法拉第常数,数据通过将不同荷电状态的电池静置放置在依次从0℃—60℃变化温度通过电压的变化梯度下求得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510418388.1A CN104951662B (zh) | 2015-07-16 | 2015-07-16 | 一种磷酸铁锂电池能量状态soe的估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510418388.1A CN104951662B (zh) | 2015-07-16 | 2015-07-16 | 一种磷酸铁锂电池能量状态soe的估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104951662A CN104951662A (zh) | 2015-09-30 |
CN104951662B true CN104951662B (zh) | 2017-11-07 |
Family
ID=54166315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510418388.1A Active CN104951662B (zh) | 2015-07-16 | 2015-07-16 | 一种磷酸铁锂电池能量状态soe的估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104951662B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597292A (zh) * | 2016-11-01 | 2017-04-26 | 深圳市沃特玛电池有限公司 | 一种电动汽车剩余能量估算方法 |
CN106707180B (zh) * | 2016-12-01 | 2020-09-11 | 深圳市麦澜创新科技有限公司 | 一种并联电池组故障检测方法 |
CN107091993B (zh) * | 2017-06-22 | 2019-08-20 | 安徽锐能科技有限公司 | 用于估计电池能量状态的装置 |
CN107015166B (zh) * | 2017-06-22 | 2019-05-07 | 安徽锐能科技有限公司 | 用于估计电池能量状态的方法及计算机可读存储介质 |
CN107238803B (zh) * | 2017-07-10 | 2019-11-01 | 安徽锐能科技有限公司 | 用于估计电池soe的方法和计算机可读存储介质 |
CN107346002B (zh) * | 2017-07-10 | 2019-11-01 | 安徽锐能科技有限公司 | 用于估计电池soe的方法和计算机可读存储介质 |
CN109507599A (zh) * | 2017-09-12 | 2019-03-22 | 北京奔驰汽车有限公司 | 一种动力电池soe的优化算法 |
CN107677965B (zh) * | 2017-09-21 | 2019-10-18 | 中国科学院广州能源研究所 | 一种锂电池能量状态估算方法 |
CN109342967A (zh) * | 2018-09-26 | 2019-02-15 | 合肥国轩高科动力能源有限公司 | 一种动力锂离子电池组低温可用能量的估算方法 |
CN109299552A (zh) * | 2018-09-29 | 2019-02-01 | 清华大学 | 一种电池能量状态的评估方法及其评估系统 |
CN110596594A (zh) * | 2019-09-23 | 2019-12-20 | 广东毓秀科技有限公司 | 一种通过大数据预测轨交锂电池soe的方法 |
CN111273177B (zh) * | 2019-10-23 | 2022-01-04 | 浙江零跑科技股份有限公司 | 一种估计电池剩余可用能量的方法 |
CN115784529B (zh) * | 2022-12-23 | 2023-08-18 | 浙江海拓环境技术有限公司 | 一种磷酸铁生产废水的资源化处理方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459551A (zh) * | 2014-11-28 | 2015-03-25 | 山东理工大学 | 一种电动汽车动力电池能量状态估算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2991105B1 (fr) * | 2012-05-24 | 2016-12-09 | Commissariat Energie Atomique | Procede de determination d'un etat d'energie d'un accumulateur electrochimique, dispositif, support et programme informatique |
-
2015
- 2015-07-16 CN CN201510418388.1A patent/CN104951662B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459551A (zh) * | 2014-11-28 | 2015-03-25 | 山东理工大学 | 一种电动汽车动力电池能量状态估算方法 |
Non-Patent Citations (3)
Title |
---|
A method for joint estimation of state-of-charge and avaliable energy of LiFePO4 batteries;Wang Yujie等;《Applied energy》;20141215;第135卷;81-87 * |
A method for state of energy estimation of lithium-ion batteries at dynamic currents and temperatures;Liu Xingtao等;《Journal of power sources》;20141215;第270卷;151-157 * |
纯电动汽车锂动力电池能量状态估算算法研究;王海峰;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20121015;C042-1321 * |
Also Published As
Publication number | Publication date |
---|---|
CN104951662A (zh) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104951662B (zh) | 一种磷酸铁锂电池能量状态soe的估算方法 | |
CN106909716B (zh) | 计及容量损耗的磷酸铁锂电池建模及soc估计方法 | |
CN106585422B (zh) | 一种动力电池soh估算方法 | |
CN105319515B (zh) | 锂离子电池荷电状态和健康状态联合估算方法 | |
CN106650134B (zh) | 一种锂离子电池系统热仿真方法 | |
Lu et al. | A method of cell-to-cell variation evaluation for battery packs in electric vehicles with charging cloud data | |
CN105866700B (zh) | 一种锂离子电池快速筛选的方法 | |
CN108732503A (zh) | 一种电池健康状态与电池容量检测方法及装置 | |
CN104051810B (zh) | 一种锂离子储能电池系统soc估算快速修正方法 | |
CN102608535A (zh) | 一种锂离子电池容量的预测方法 | |
CN105353316B (zh) | 动力电池充电时soc变化量及充电电量折算系数测量方法 | |
CN109900937B (zh) | 一种具有温度补偿功能的锂电池电荷状态估算方法 | |
CN105467328A (zh) | 一种锂离子电池荷电状态估计方法 | |
CN105116350B (zh) | 动力电池放电时soc变化量及放电电量折算系数测量方法 | |
CN109839599B (zh) | 基于二阶ekf算法的锂离子电池soc估计方法 | |
Feng et al. | A graphical model for evaluating the status of series‐connected lithium‐ion battery pack | |
CN107677965A (zh) | 一种锂电池能量状态估算方法 | |
CN103728570B (zh) | 一种基于电池热特性的健康状态检测方法 | |
CN103424712A (zh) | 基于粒子群优化的电池剩余容量在线测量方法 | |
CN112098849A (zh) | 基于求积分卡尔曼滤波的锂电池剩余电量估计方法 | |
CN112363077A (zh) | 一种电池健康度容量预估方法 | |
CN109557478A (zh) | 一种锂离子电池安全性能评估方法 | |
CN108829911A (zh) | 一种开路电压与soc函数关系优化方法 | |
CN105974320A (zh) | 一种液态或半液态金属电池荷电状态估计方法 | |
Zhang et al. | An improved Peukert battery model of nonlinear capacity considering temperature effect |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |