CN103728567A - Charge state estimation method based on initial value optimization - Google Patents

Charge state estimation method based on initial value optimization Download PDF

Info

Publication number
CN103728567A
CN103728567A CN201310753696.0A CN201310753696A CN103728567A CN 103728567 A CN103728567 A CN 103728567A CN 201310753696 A CN201310753696 A CN 201310753696A CN 103728567 A CN103728567 A CN 103728567A
Authority
CN
China
Prior art keywords
soc
battery
initial value
charge
draw
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.)
Granted
Application number
CN201310753696.0A
Other languages
Chinese (zh)
Other versions
CN103728567B (en
Inventor
向勇
范正航
毕闯
向俊杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Angao Special Electric Technology Co ltd
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201310753696.0A priority Critical patent/CN103728567B/en
Publication of CN103728567A publication Critical patent/CN103728567A/en
Application granted granted Critical
Publication of CN103728567B publication Critical patent/CN103728567B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention belongs to the technical field of automobile power batteries, and particularly relates to an SOC estimation method of an automobile power battery. Charging and discharging are conducted on the battery over and over again, then standing is performed, the voltage value of the battery in the process from charging and discharging to standing is detected, stable voltage time t(off) is recorded, and an initial value SOC(t0+delta t) is obtained. The battery is used within recovery time, a precise SOC initial value can be provided for the battery, and thus the SOC value estimated by expanding a Kalman filtering method can be converged more quickly.

Description

A kind of state-of-charge evaluation method based on optimizing initial value
Technical field
The invention belongs to automobile power cell technical field, relate in particular to the method for automobile power cell SOC estimation.
Background technology
SOC is one of important parameter of battery, refers to the dump energy of battery, and the running status that SOC is battery provides important evidence.But SOC is affected by several factors, its value can only be estimated by measuring these factor values.
The method of estimation of electrokinetic cell SOC mainly contains discharge test method, open-circuit voltage method, Ah counting method, internal resistance measurement method, Kalman filtering method, load method and neural network etc. both at home and abroad at present.The most frequently used: open-circuit voltage method, Ah counting method, Kalman filtering method.
Open-circuit voltage method: determine SOC size by the funtcional relationship of open circuit electromotive force and SOC.In electric discharge latter stage, the effect of test is better, but uses and to have a recovery stage later at battery, at the corresponding relation of this stage voltage and SOC, is not very obvious.
Ah counting method: Ah counting method is to estimate the SOC of battery by calculating the accumulation electric weight of battery when the charge or discharge, but inaccurate current measurement will increase SOC evaluated error, and through accumulating for a long time, it is increasing that this error can become.
Kalman filtering method: Kalman filtering method is a kind of minimum variance estimate method of recursion, utilized the value in a upper moment and the parameter value of measurement in real time to estimate.When Kalman filtering method is estimated for electrokinetic cell SOC, be a discrete system, need an accurate model, and calculated amount be very large.In the research of electrokinetic cell SOC, there is the correction algorithm of a lot of Kalman filtering methods, these algorithms have been optimized estimation precision to a great extent.
Because Kalman filtering method has error correcting action, so the SOC algorithm of great majority based on Kalman filtering all tends to its correction factor improve or use in conjunction with other algorithms and use, and ignore the precise decreasing that initial value error is brought.The SOC algorithm (application number 201110142292.9) of speeding scientific and technological such as Zhenjiang perseverance is exactly to carry out Integrated using by Ah counting method and Kalman filtering method, to improve its robustness, convergence, and adaptivity.In the patent of LG, (application number 200680016312.5) revised efficiency factor.But in the arranging of SOC initial value, or use experience method is weighted estimation.Its precision is not very high.
Summary of the invention
The object of the present invention is to provide a kind of state-of-charge evaluation method based on optimizing initial value, comprise the steps:
S1, the parameters initialization to battery;
S2, estimate battery in a upper operational phase because the quantity of electric charge of polarization effect accumulation;
S3, battery is is repeatedly discharged and recharged, then carry out standing, and detect battery described discharge and recharge standing process in magnitude of voltage, the time t (off) that recording voltage is stable;
S4, t described in S3 (off) is carried out to linearization process;
S5, draw initial value SOC ( t 0 + Δt ) = SOC ( t 0 ) + ∫ t 0 t 0 + Δt [ Q ( up ) / t ( off ) ] dt .
The quantity of electric charge computing method of accumulating described in S2 further, are as follows:
Under steady temperature:
Charging: power consumption=SOC increment+stored charge amount
?
Figure BDA0000451730190000022
by Δ SOC corresponding to different current values, draw
Figure BDA0000451730190000023
by matlab to different current values and α corresponding to different current value 1(I) carry out matching, draw α 1(I)=p 1i 4+ p 2i 3+ p 3i 2+ p 4i 1+ p 5, α 1(I) be that charge efficiency is about the coefficient of electric current;
Electric discharge: SOC reduction=power consumption+stored charge amount
?
Figure BDA0000451730190000024
by Δ SOC corresponding to different current values, draw
Figure BDA0000451730190000025
by matlab to different current values and α corresponding to different current value 2(I) carry out matching, draw α 2(I)=p 1i 4+ p 2i 3+ p 3i 2+ p 4i 1+ p 5, α 2(I) be discharging efficiency about the coefficient of electric current,
Wherein
Figure BDA0000451730190000026
for emitting the quantity of electric charge, Δ SOC is SOC reduction, and Q (up) is stored charge amount, p 1i 4for being matching parameter value out, current i is independent variable, and p is exactly the coefficient of independent variable,
Because η=α (I) β (T), so, in the charging stage at discharge regime
Figure BDA0000451730190000028
wherein, η is the function about electric current I and temperature T.
Further, Δ t≤t (off).
The invention has the beneficial effects as follows: within release time, use battery, can provide a comparatively accurate SOC initial value to battery, make more rapid convergence of SOC value that EKF method estimates.
Accompanying drawing explanation
Fig. 1 is evaluation method process flow diagram of the present invention.
Fig. 2 is model of the present invention.
When Fig. 3 is charging, efficiency is about the figure of electric current
Initial value estimation figure when Fig. 4 is charging.
Fig. 5 be under discharge scenario α about the function of I.
Wherein, E is standard voltage source, and R is ohmic internal resistance, and U (up) is polarizing voltage.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is described:
As shown in Fig. 2 model, set up Kalman filtering algorithm.Wherein, E is standard voltage source, and R is ohmic internal resistance, and U (up) is polarizing voltage.
S1, the parameters initialization to battery;
S2, estimate battery in a upper operational phase because the quantity of electric charge of polarization effect accumulation comprises:
Under steady temperature:
Charging: power consumption=SOC increment+stored charge amount
?
Figure BDA0000451730190000031
by Δ SOC corresponding to different current values, draw
Figure BDA0000451730190000032
by matlab to different current values and α corresponding to different current value 1(I) carry out matching, draw α 1(I)=p 1i 4+ p 2i 3+ p 3i 2+ p 4i 1+ p 5, α 1(I) be that charge efficiency is about the coefficient of electric current;
Electric discharge: SOC reduction=power consumption+stored charge amount
?
Figure BDA0000451730190000033
by Δ SOC corresponding to different current values, draw by matlab to different current values and α corresponding to different current value 2(I) carry out matching, draw α 2(I)=p 1i 4+ p 2i 3+ p 3i 2+ p 4i 1+ p 5, α 2(I) be discharging efficiency about the coefficient of electric current,
Wherein
Figure BDA0000451730190000035
for emitting the quantity of electric charge, Δ SOC is SOC reduction, and Q (up) is stored charge amount, p 1i 4for being matching parameter value out, current i is independent variable, and p is exactly the coefficient of independent variable,
Because η=α (I) β (T), so, in the charging stage
Figure BDA0000451730190000036
at discharge regime
Figure BDA0000451730190000037
wherein, η is the function about electric current I and temperature T.
S3, battery is is repeatedly discharged and recharged, then carry out standing, and detect battery described discharge and recharge standing process in magnitude of voltage, the time t (off) that recording voltage is stable;
S4, t described in S3 (off) is carried out to linearization process;
S5, draw initial value SOC ( t 0 + Δt ) = SOC ( t 0 ) + ∫ t 0 t 0 + Δt [ Q ( up ) / t ( off ) ] dt , Δt≤t(off);
S6, by initial value SOC (t described in S5 0+ Δ t) for expanded Kalman filtration algorithm model and carry out SOC estimation.
According to matlab models fitting, show that under discharge scenario α is about the function of electric current I:
F (x)=p1*i^4+p2*i^3+p3*i^2+p4*i+p5, wherein:
p1=-5.23e+011(-5.469e+011,-4.992e+011)
p2=7.233e+011(6.904e+011,7.563e+011)
p3=-3.751e+011(-3.922e+011,-3.58e+011)
p4=8.647e+010(8.253e+010,9.041e+010)
p5=-7.474e+009(-7.815e+009,-7.133e+009)
According to matlab models fitting, show that under charge condition α is about the function of electric current I:
F (x)=p1*x^4+p2*x^3+p3*x^2+p4*x+p5, wherein:
p1=-0.09137(-0.4133,0.2306)
p2=0.979(-2.161,4.119)
p3=-3.527(-13.52,6.468)
p4=4.833(-7.766,17.43)
P5=-1.144 (6.546,4.258), wherein, because fitting function selection is 4 rank, so most significant digit is 4 rank, also can use 5 rank, 6 rank, but conventionally use 4 rank.

Claims (3)

1. the state-of-charge evaluation method based on optimizing initial value, is characterized in that: comprise the steps:
S1, the parameters initialization to battery;
S2, estimate battery in a upper operational phase because the quantity of electric charge of polarization effect accumulation;
S3, battery is is repeatedly discharged and recharged, then carry out standing, and detect battery described discharge and recharge standing process in magnitude of voltage, the time t (off) that recording voltage is stable;
S4, t described in S3 (off) is carried out to linearization process;
S5, draw initial value SOC ( t 0 + Δt ) = SOC ( t 0 ) + ∫ t 0 t 0 + Δt [ Q ( up ) / t ( off ) ] dt .
2. a kind of state-of-charge evaluation method based on optimizing initial value according to claim 1, is characterized in that: the quantity of electric charge computing method of accumulating described in S2 are as follows:
Under steady temperature:
Charging: power consumption=SOC increment+stored charge amount
? by Δ SOC corresponding to different current values, draw
Figure FDA0000451730180000013
by matlab to different current values and α corresponding to different current value 1(I) carry out matching, draw α 1(I)=p 1i 4+ p 2i 3+ p 3i 2+ p 4i 1+ p 5, α 1(I) be that charge efficiency is about the coefficient of electric current;
Electric discharge: SOC reduction=power consumption+stored charge amount
?
Figure FDA0000451730180000014
by Δ SOC corresponding to different current values, draw
Figure FDA0000451730180000015
by matlab to different current values and α corresponding to different current value 2(I) carry out matching, draw α 2(I)=p 1i 4+ p 2i 3+ p 3i 2+ p 4i 1+ p 5, α 2(I) be discharging efficiency about the coefficient of electric current,
Wherein
Figure FDA0000451730180000016
for emitting the quantity of electric charge, Δ SOC is SOC reduction, and Q (up) is stored charge amount, p 1i 4for being matching parameter value out, current i is independent variable, and p is exactly the coefficient of independent variable,
Because η=α (I) β (T), so, in the charging stage
Figure FDA0000451730180000017
at discharge regime
Figure FDA0000451730180000018
wherein, η is the function about electric current I and temperature T.
3. a kind of state-of-charge evaluation method based on optimizing initial value according to claim 1, is characterized in that: Δ t≤t (off).
CN201310753696.0A 2013-12-31 2013-12-31 A kind of charge state estimation method based on optimized initial value Active CN103728567B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310753696.0A CN103728567B (en) 2013-12-31 2013-12-31 A kind of charge state estimation method based on optimized initial value

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310753696.0A CN103728567B (en) 2013-12-31 2013-12-31 A kind of charge state estimation method based on optimized initial value

Publications (2)

Publication Number Publication Date
CN103728567A true CN103728567A (en) 2014-04-16
CN103728567B CN103728567B (en) 2016-06-08

Family

ID=50452725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310753696.0A Active CN103728567B (en) 2013-12-31 2013-12-31 A kind of charge state estimation method based on optimized initial value

Country Status (1)

Country Link
CN (1) CN103728567B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105891717A (en) * 2015-07-01 2016-08-24 乐视移动智能信息技术(北京)有限公司 Battery electric quantity obtaining method and device

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5119011A (en) * 1990-08-08 1992-06-02 General Electric Company Battery state of charge indicator
JPH10319100A (en) * 1997-05-15 1998-12-04 Nippon Soken Inc Battery charging state detecting device
CN101138142A (en) * 2005-03-09 2008-03-05 Lg化学株式会社 Method of setting initial value of soc of battery using OCV temperature hysteresis
WO2008026477A1 (en) * 2006-08-29 2008-03-06 Nec Corporation Method and device for estimating soc value of secondary battery and degradation judging method and device
CN101176234A (en) * 2005-05-11 2008-05-07 Lg化学株式会社 Method of estimating SOC of battery for hybrid electric vehicle
CN102788957A (en) * 2011-05-20 2012-11-21 镇江恒驰科技有限公司 Estimating method of charge state of power battery
CN103439668A (en) * 2013-09-05 2013-12-11 桂林电子科技大学 Charge state evaluation method and system of power lithium ion battery

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5119011A (en) * 1990-08-08 1992-06-02 General Electric Company Battery state of charge indicator
JPH10319100A (en) * 1997-05-15 1998-12-04 Nippon Soken Inc Battery charging state detecting device
CN101138142A (en) * 2005-03-09 2008-03-05 Lg化学株式会社 Method of setting initial value of soc of battery using OCV temperature hysteresis
CN101176234A (en) * 2005-05-11 2008-05-07 Lg化学株式会社 Method of estimating SOC of battery for hybrid electric vehicle
WO2008026477A1 (en) * 2006-08-29 2008-03-06 Nec Corporation Method and device for estimating soc value of secondary battery and degradation judging method and device
CN102788957A (en) * 2011-05-20 2012-11-21 镇江恒驰科技有限公司 Estimating method of charge state of power battery
CN103439668A (en) * 2013-09-05 2013-12-11 桂林电子科技大学 Charge state evaluation method and system of power lithium ion battery

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘浩: "纯电动汽车用锂离子电池SOC估算方案的研究", 《电子技术·冶金电气》 *
华贵山等: "电动汽车锂离子电池SOC估算策略研究", 《电源技术》 *
李秉宇等: "基于卡尔曼滤波器的蓄电池剩余容量估算法", 《电源技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105891717A (en) * 2015-07-01 2016-08-24 乐视移动智能信息技术(北京)有限公司 Battery electric quantity obtaining method and device

Also Published As

Publication number Publication date
CN103728567B (en) 2016-06-08

Similar Documents

Publication Publication Date Title
CN107991623B (en) Battery ampere-hour integral SOC estimation method considering temperature and aging degree
CN103472403B (en) A kind of electrokinetic cell SOC compound method of estimation based on PNGV equivalent-circuit model
CN103424710B (en) For monitoring the method and system that the performance of the aged monomer in set of cells changes
CN103675706B (en) A kind of power battery electric charge quantity estimation method
CN105891729B (en) The condition detection method and device of battery and battery pack
JP5466564B2 (en) Battery degradation estimation method, battery capacity estimation method, battery capacity equalization method, and battery degradation estimation apparatus
Ng et al. Enhanced coulomb counting method for estimating state-of-charge and state-of-health of lithium-ion batteries
CN102645637B (en) Method for estimating SOC (state of charge) of equalized batteries
WO2015106691A1 (en) Soc estimation method for power battery for hybrid electric vehicle
CN104502855B (en) Lead-acid accumulator SOH detection methods and system
CN103823189B (en) Method for calculating residual capacity of power battery pack
CN102749588B (en) Method for fault diagnosis on basis of storage battery state of charge (SOC) and state of health (SOH)
CN104678316A (en) Method and device for estimating charge state of lithium ion battery
JP5419831B2 (en) Battery degradation degree estimation device
CN104360286B (en) Lithium ion battery charge state estimation modification method
CN104459551A (en) Electric vehicle power battery state-of-energy estimation method
CN105068008A (en) Battery SOC (state of charge) estimation method by utilizing vehicle-mounted charging machine identification battery parameter
CN107247235A (en) A kind of battery capacity evaluation method for considering batteries in parallel connection difference
CN103328997A (en) Device for estimating state of charge of battery
CN105353313A (en) Estimation method of battery state of charge and estimation device of battery state of charge
CN108717164A (en) The state-of-charge SOC scaling methods and system of battery
CN104535932A (en) Lithium ion battery charge state estimating method
CN103364736A (en) Method for calculating RAC (residual available capacity) of lithium ion battery pack
CN103675683A (en) Lithium battery state of charge (SOC) estimation method
CN104407298A (en) Lithium ion battery pack available surplus capacity calculation method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240204

Address after: 610000 Chengdu Future Science and Technology City, Chengdu, Sichuan Province (No. 1-2, Group 8, Dingshui Village, Futian Street, Eastern New Area)

Patentee after: Sichuan Angao Special Electric Technology Co.,Ltd.

Country or region after: China

Address before: 611731, No. 2006, West Avenue, Chengdu hi tech Zone (West District, Sichuan)

Patentee before: University of Electronic Science and Technology of China

Country or region before: China

TR01 Transfer of patent right