Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7616101B2 - SOC Estimation Method - Google Patents
[go: Go Back, main page]

JP7616101B2 - SOC Estimation Method - Google Patents

SOC Estimation Method Download PDF

Info

Publication number
JP7616101B2
JP7616101B2 JP2022006288A JP2022006288A JP7616101B2 JP 7616101 B2 JP7616101 B2 JP 7616101B2 JP 2022006288 A JP2022006288 A JP 2022006288A JP 2022006288 A JP2022006288 A JP 2022006288A JP 7616101 B2 JP7616101 B2 JP 7616101B2
Authority
JP
Japan
Prior art keywords
battery
soc
positive electrode
iron phosphate
ocv
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
Application number
JP2022006288A
Other languages
Japanese (ja)
Other versions
JP2023105454A (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.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
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 Toyota Motor Corp filed Critical Toyota Motor Corp
Priority to JP2022006288A priority Critical patent/JP7616101B2/en
Publication of JP2023105454A publication Critical patent/JP2023105454A/en
Application granted granted Critical
Publication of JP7616101B2 publication Critical patent/JP7616101B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries

Landscapes

  • Tests Of Electric Status Of Batteries (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)
  • Secondary Cells (AREA)
  • Battery Electrode And Active Subsutance (AREA)

Description

本開示は、リン酸鉄リチウムにより形成された正極を含むリン酸鉄リチウムイオン電池のSOC推定方法に関する。 This disclosure relates to a method for estimating the SOC of a lithium iron phosphate ion battery that includes a positive electrode formed from lithium iron phosphate.

近年、電気自動車(BEV)やハイブリッド車両(HEV)に搭載される二次電池や、車両に搭載される鉛蓄電池(12Vバッテリ)の代替品として、例えばLiFePO4等のリン酸鉄リチウムにより形成された正極(以下、適宜「LiFePO正極」という。)を含むリン酸鉄リチウムイオン電池が採用されつつある。リン酸鉄リチウムイオン電池のLiFePO正極は、メモリ効果や経路依存性、充電時と放電時とでOCV(開回路電圧)が変化する電圧ヒステリシスといった特異な挙動(分極挙動)を示すものである。また、LiFePO正極の挙動を説明するためのモデルとしては、LiFePO正極が互いに電子的およびイオン的に接続された粒子半径の異なる多数の球状粒子により構成されていると仮定した多粒子モデル(Many-Particle Model)が知られている(例えば、非特許文献1参照)。かかる多粒子モデルを用いることで、LiFePO正極の挙動を包括的に捉えることができる。 In recent years, lithium iron phosphate ion batteries including a positive electrode (hereinafter referred to as "LiFePO positive electrode") formed of lithium iron phosphate such as LiFePO4 are being adopted as a replacement for secondary batteries mounted on electric vehicles (BEVs) and hybrid vehicles (HEVs) and lead-acid batteries (12V batteries) mounted on vehicles. The LiFePO positive electrode of a lithium iron phosphate ion battery exhibits peculiar behavior (polarization behavior) such as memory effect, path dependency, and voltage hysteresis in which the OCV (open circuit voltage) changes between charging and discharging. In addition, as a model for explaining the behavior of the LiFePO positive electrode, a many-particle model is known in which the LiFePO positive electrode is assumed to be composed of a large number of spherical particles with different particle radii that are electronically and ionically connected to each other (see, for example, Non-Patent Document 1). By using such a multi-particle model, the behavior of the LiFePO positive electrode can be comprehensively grasped.

Comprehensive Study of the Polarization Behavior of LiFePO4 Electrodes Based on a Many-Particle Model Hiroki Kondo, Tsuyoshi Sasaki, Pallab Barai and Venkat Srinivasan Journal of The Electrochemical Society, Volume 165, Number 10 Published 6 July 2018Comprehensive Study of the Polarization Behavior of LiFePO4 Electrodes Based on a Many-Particle Model Hiroki Kondo, Tsuyoshi Sasaki, Pallab Barai and Venkat Srinivasan Journal of The Electrochemical Society, Volume 165, Number 10 Published 6 July 2018

ここで、三元系リチウムイオン電池等においては、一般に、電流の測定値を積算してSOCを推定する電流積算法と、SOCとOCVとの関係を用いた推定手法とを組み合わせてSOCが推定される。これに対して、リン酸鉄リチウムイオン電池では、高SOC域および低SOC域以外の広い範囲でSOCの変化に対するOCVの変化率が非常に小さくなり、しかも上記電圧ヒステリシス等が存在するので、SOCとOCVとの関係からSOCを精度よく推定することは困難である。このため、リン酸鉄リチウムイオン電池のSOCを精度よく推定する手法が求められているが、上記非特許文献1には、リン酸鉄リチウムイオン電池のSOCの推定手法が何ら記載されていない。 In ternary lithium-ion batteries and the like, the SOC is generally estimated by combining a current integration method, which estimates the SOC by integrating measured current values, with an estimation method using the relationship between SOC and OCV. In contrast, in lithium iron phosphate batteries, the rate of change in OCV relative to changes in SOC is very small over a wide range other than the high and low SOC ranges, and the above-mentioned voltage hysteresis and the like exist, making it difficult to accurately estimate the SOC from the relationship between SOC and OCV. For this reason, a method for accurately estimating the SOC of lithium iron phosphate batteries is required, but the above-mentioned non-patent document 1 does not disclose any method for estimating the SOC of lithium iron phosphate batteries.

そこで、本開示は、リン酸鉄リチウムにより形成された正極を含むリン酸鉄リチウムイオン電池のSOCを精度よく推定可能にすることを主目的とする。 Therefore, the primary objective of this disclosure is to enable accurate estimation of the SOC of a lithium iron phosphate ion battery that includes a positive electrode formed from lithium iron phosphate.

本開示のSOC推定方法は、リン酸鉄リチウムにより形成された正極を含むリン酸鉄リチウムイオン電池のSOCを推定するSOC推定方法であって、少なくとも、前記リン酸鉄リチウムイオン電池の電圧の測定値と、前記リン酸鉄リチウムイオン電池のSOCとOCVとの関係を規定するように予め定められた第1マップから取得されたOCVとに基づいて、前記正極の反応による過電圧を算出し、前記正極の反応による過電圧と、前記リン酸鉄リチウムのリチウム組成と前記正極のOCVとの関係に基づいて予め定められた第2マップとを用いて、前記正極がそれぞれリチウムを含む複数の粒子により構成されていると仮定したときの各粒子の過電圧を算出し、前記各粒子の過電圧に基づいて前記各粒子の電流密度を算出し、前記各粒子の電流密度に基づいて前記リン酸鉄リチウムイオン電池の推定電池電流を算出し、前記推定電池電流に基づいて前記リン酸鉄リチウムイオン電池のSOCを推定するものである。 The SOC estimation method disclosed herein is an SOC estimation method for estimating the SOC of a lithium iron phosphate ion battery including a positive electrode formed of lithium iron phosphate, and includes: calculating an overvoltage caused by a reaction of the positive electrode based on at least a measured voltage of the lithium iron phosphate ion battery and an OCV obtained from a first map that is predetermined to define the relationship between the SOC and OCV of the lithium iron phosphate ion battery; calculating an overvoltage of each particle when it is assumed that the positive electrode is composed of a plurality of particles each containing lithium, using the overvoltage caused by the reaction of the positive electrode and a second map that is predetermined based on the relationship between the lithium composition of the lithium iron phosphate and the OCV of the positive electrode; calculating a current density of each particle based on the overvoltage of each particle; calculating an estimated battery current of the lithium iron phosphate ion battery based on the current density of each particle; and estimating the SOC of the lithium iron phosphate ion battery based on the estimated battery current.

本開示のSOC推定方法は、リン酸鉄リチウムにより形成された正極がそれぞれリチウムを含む複数の粒子により構成されていると仮定する多粒子モデルを利用するものである。リン酸鉄リチウムイオン電池のSOCの推定に際しては、少なくとも、リン酸鉄リチウムイオン電池の電圧の測定値と、リン酸鉄リチウムイオン電池のSOCとOCVとの関係を規定するように予め定められた第1マップから取得されたOCVとに基づいて、正極の反応による過電圧が算出される。更に、正極の反応による過電圧と、リン酸鉄リチウムのリチウム組成と正極のOCVとの関係に基づいて予め定められた第2マップとを用いて、正極を構成する複数の粒子の各々の過電圧が算出される。そして、複数の粒子の各々の過電圧に基づいて複数の粒子の各々の電流密度が算出され、複数の粒子の各々の電流密度からリン酸鉄リチウムイオン電池の推定電池電流が算出される。これにより、第1マップから取得されたOCVと実際のOCVとのずれを相殺するように算出された推定電池電流に基づいてリン酸鉄リチウムイオン電池のSOCを推定することができる。この結果、メモリ効果や経路依存性、電圧ヒステリシスといった特異な挙動を示すリン酸鉄リチウムにより形成された正極を含むリン酸鉄リチウムイオン電池のSOCを精度よく推定することが可能となる。 The SOC estimation method of the present disclosure uses a multi-particle model that assumes that the positive electrode formed of lithium iron phosphate is composed of multiple particles each containing lithium. When estimating the SOC of the lithium iron phosphate ion battery, the overvoltage due to the reaction of the positive electrode is calculated based on at least the measured voltage of the lithium iron phosphate ion battery and the OCV obtained from a first map that is predetermined to define the relationship between the SOC and OCV of the lithium iron phosphate ion battery. Furthermore, the overvoltage of each of the multiple particles that constitute the positive electrode is calculated using the overvoltage due to the reaction of the positive electrode and a second map that is predetermined based on the relationship between the lithium composition of the lithium iron phosphate and the OCV of the positive electrode. Then, the current density of each of the multiple particles is calculated based on the overvoltage of each of the multiple particles, and the estimated battery current of the lithium iron phosphate ion battery is calculated from the current density of each of the multiple particles. As a result, the SOC of the lithium iron phosphate ion battery can be estimated based on the estimated battery current calculated to offset the deviation between the OCV obtained from the first map and the actual OCV. As a result, it is possible to accurately estimate the SOC of a lithium iron phosphate-ion battery that includes a positive electrode made of lithium iron phosphate, which exhibits peculiar behaviors such as memory effect, path dependency, and voltage hysteresis.

また、前記リン酸鉄リチウムイオン電池の電圧の測定値から、少なくとも前記第1マップから取得された前回のSOCに対応した前回のOCVを減じて前記正極の反応による過電圧を算出してもよく、前記正極の反応による過電圧と、前記第2マップから取得された前記各粒子の基準電圧とから、前記各粒子の過電圧を算出してもよく、前回のSOCと、前記推定電池電流と、電池容量とから前記リン酸鉄リチウムイオン電池の今回のSOCを推定すると共に、前記第1マップから今回のSOCに対応した今回のOCVを取得してもよい。これにより、リン酸鉄リチウムイオン電池のSOCの推定を繰り返し実行することで当該SOCの推定値を実際の値に収束させていくことが可能となる。 In addition, the overvoltage due to the reaction of the positive electrode may be calculated by subtracting from the measured voltage value of the lithium iron phosphate battery at least the previous OCV corresponding to the previous SOC obtained from the first map, and the overvoltage of each particle may be calculated from the overvoltage due to the reaction of the positive electrode and the reference voltage of each particle obtained from the second map. The current SOC of the lithium iron phosphate battery may be estimated from the previous SOC, the estimated battery current, and the battery capacity, and the current OCV corresponding to the current SOC may be obtained from the first map. This makes it possible to converge the estimated SOC value to an actual value by repeatedly estimating the SOC of the lithium iron phosphate battery.

更に、前記第1マップは、SOCごとに、前記正極の2つのスピノーダル点の平均電圧から前記リン酸鉄リチウムイオン電池の負極のOCVを減じることにより作成されてもよく、前記第2マップは、前記リン酸鉄リチウムのリチウム組成と前記正極のOCVとの関係を前記正極の2つのスピノーダル点の平均電圧が0ボルトになるようにオフセットすることにより作成されてもよい。これにより、正極の反応による過電圧と、第2マップとを用いて、複数の粒子の各々の過電圧をより適正に算出することが可能となる。 Furthermore, the first map may be created by subtracting the OCV of the negative electrode of the lithium iron phosphate battery from the average voltage of the two spinodal points of the positive electrode for each SOC, and the second map may be created by offsetting the relationship between the lithium composition of the lithium iron phosphate and the OCV of the positive electrode so that the average voltage of the two spinodal points of the positive electrode becomes 0 volts. This makes it possible to more accurately calculate the overvoltage of each of the multiple particles using the overvoltage due to the reaction of the positive electrode and the second map.

また、前記各粒子の電流密度に基づいて前記各粒子のリチウム濃度の変化量を算出すると共に、前回のリチウム濃度と前記変化量とから今回のリチウム濃度を算出してもよく、前記第2マップから前回のリチウム濃度に基づくリチウム組成に対応した前記各粒子の基準電圧が取得されてもよい。これにより、複数の粒子の各々の基準電圧ひいては過電圧をより適正に算出することが可能となる。 In addition, the amount of change in the lithium concentration of each particle may be calculated based on the current density of each particle, and the current lithium concentration may be calculated from the previous lithium concentration and the amount of change, and the reference voltage of each particle corresponding to the lithium composition based on the previous lithium concentration may be obtained from the second map. This makes it possible to more appropriately calculate the reference voltage and therefore the overvoltage of each of the multiple particles.

更に、前記リン酸鉄リチウムイオン電池の電圧の測定値から、前回のOCVと、塩濃度過電圧と、前記リン酸鉄リチウムイオン電池の電流の測定値と内部抵抗との積値とを減じて、前記正極の反応による過電圧を算出してもよい。これにより、演算負荷を低減しつつ、正極の反応による過電圧をより適正に算出することが可能となる。 Furthermore, the overvoltage due to the reaction of the positive electrode may be calculated by subtracting the previous OCV, the salt concentration overvoltage, and the product of the measured current and internal resistance of the lithium iron phosphate battery from the measured voltage value of the lithium iron phosphate battery. This makes it possible to more accurately calculate the overvoltage due to the reaction of the positive electrode while reducing the computational load.

本開示のSOC推定方法が適用されるリン酸鉄リチウムイオン電池を含むバッテリを搭載した車両の概略構成図である。1 is a schematic configuration diagram of a vehicle equipped with a battery including a lithium iron phosphate battery to which an SOC estimation method of the present disclosure is applied; リン酸鉄リチウムのリチウム組成と、LiFePO正極のOCVとの関係を示す図表である。1 is a diagram showing the relationship between the lithium composition of lithium iron phosphate and the OCV of a LiFePO positive electrode. (a)は、リン酸鉄リチウムイオン電池のSOCと正極のOCVの平均電圧との関係を示す図表であり、(b)は、リン酸鉄リチウムイオン電池のSOCと負極のOCVとの関係を示す図表であり、(c)は、リン酸鉄リチウムイオン電池のSOCとOCVとの関係を規定した第1マップの一例を示す図表である。1A is a graph showing the relationship between the SOC of a lithium iron phosphate ion battery and the average voltage of the OCV of the positive electrode, FIG. 1B is a graph showing the relationship between the SOC of a lithium iron phosphate ion battery and the OCV of the negative electrode, and FIG. 1C is a graph showing an example of a first map that specifies the relationship between the SOC and OCV of a lithium iron phosphate ion battery. 本開示のSOC推定方法を説明するためのフローチャートである。4 is a flowchart for explaining an SOC estimation method according to the present disclosure. 本開示のSOC推定方法において用いられる第2マップの一例を示す図表である。4 is a chart showing an example of a second map used in the SOC estimation method of the present disclosure.

次に、図面を参照しながら、本開示の発明を実施するための形態について説明する。 Next, the form for implementing the invention of this disclosure will be described with reference to the drawings.

図1は、本開示のSOC推定方法が適用されるリン酸鉄リチウムイオン電池を含むバッテリ1を搭載した車両Vを示す概略構成図である。同図に示す車両Vは、バッテリ1に加えて、インバータ等を含む電力制御装置(図示省略)を介してバッテリ1に接続されると共に当該バッテリ1と電力をやり取りして走行用の動力や回生制動力を出力可能なモータジェネレータ(三相交流電動機)MGを含む電気自動車(BEV)あるいはハイブリッド車両(HEV、PHEV)である。 Figure 1 is a schematic diagram showing a vehicle V equipped with a battery 1 including a lithium iron phosphate battery to which the SOC estimation method of the present disclosure is applied. The vehicle V shown in the figure is an electric vehicle (BEV) or hybrid vehicle (HEV, PHEV) that includes, in addition to the battery 1, a motor generator (three-phase AC motor) MG that is connected to the battery 1 via a power control device (not shown) including an inverter and the like and that can exchange power with the battery 1 to output driving power and regenerative braking force.

バッテリ1は、例えば直列に接続される複数の電池モジュールを含む電池スタックと、当該電池スタックを収容するケースとを含む、いわゆる高電圧バッテリであり、各電池モジュールは、例えば直列に接続されると共にモジュールケース内に収容される複数の電池セル2を含むものである。電池モジュールを構成する電池セル2は、本実施形態では、LiFePO4であるオリビン型の結晶構造を有するリン酸鉄リチウムにより形成された正極(LiFePO正極)と、黒鉛系炭素材料等により形成された負極とを含むリン酸鉄リチウムイオン電池であり、正極および負極は、セパレータや有機溶媒等である電解液と共に外装体の内部に収容される。 The battery 1 is a so-called high-voltage battery including a battery stack including a plurality of battery modules connected in series, for example, and a case that houses the battery stack, and each battery module includes a plurality of battery cells 2 that are connected in series and housed in a module case, for example. In this embodiment, the battery cells 2 constituting the battery module are lithium iron phosphate ion batteries that include a positive electrode (LiFePO positive electrode ) formed of lithium iron phosphate having an olivine crystal structure, which is LiFePO4, and a negative electrode formed of a graphite-based carbon material or the like, and the positive electrode and the negative electrode are housed inside an exterior body together with a separator and an electrolyte solution, which is an organic solvent or the like.

図2は、実験により得られた、リン酸鉄リチウムのリチウム組成x(LixFePO4における“x”)と、LiFePO正極のOCV(開回路電圧)との関係を示す図表である。同図に示すように、リン酸鉄リチウムイオン電池が充放電される際、LiFePO正極では、Li-rich相と、Li-poor相と、当該二相が共存する領域とが存在する状態で反応が進行する。また、LiFePO正極のOCVは、図示するように、Li-rich相側のスピノーダル点と、Li-poor相側のスピノーダル点との間でリチウム組成xの減少に概ね比例して減少し、2つのスピノーダル点間の電圧差は、およそ20mVとなる。更に、2つのスピノーダル点の平均電圧(Li-rich相およびLi-poor相の共存相における平衡電圧)は、およそ3.43Vとなる。 FIG. 2 is a diagram showing the relationship between the lithium composition x ("x" in LixFePO4 ) of lithium iron phosphate and the OCV (open circuit voltage) of the LiFePO positive electrode obtained by an experiment. As shown in the figure, when a lithium iron phosphate ion battery is charged and discharged, the reaction proceeds in the LiFePO positive electrode in a state in which the Li-rich phase, the Li-poor phase, and a region in which the two phases coexist are present. In addition, as shown in the figure, the OCV of the LiFePO positive electrode decreases approximately in proportion to the decrease in the lithium composition x between the spinodal point on the Li-rich phase side and the spinodal point on the Li-poor phase side, and the voltage difference between the two spinodal points is approximately 20 mV. Furthermore, the average voltage of the two spinodal points (the equilibrium voltage in the coexistence phase of the Li-rich phase and the Li-poor phase) is approximately 3.43 V.

ここで、メモリ効果や経路依存性、電圧ヒステリシスといったLiFePO正極の特異な挙動(分極挙動)を考慮して、LiFePO正極が互いに電子的およびイオン的に接続された粒子半径の異なるN個の球状粒子により構成(形成)されると仮定する多粒子モデル(Many-Particle Model)を導入した場合、図2より、複数の粒子の各々のリチウム組成xから各粒子のOCVを得ることができる。また、LiFePO正極の反応による過電圧が発生すると、当該LiFePO正極を構成する複数の粒子の各々は、リチウム組成に応じた過電圧ηjをもつ(ただし、“j”は、1からNまでの整数である。)。更に、過電圧ηjは、リン酸鉄リチウムイオン電池の充放電時におけるLiFePO正極の反応による過電圧を“ΔVMPM”とし、図2から得られるj番目の粒子のOCVを“OCVj”とし、上記2つのスピノーダル点の平均電圧を“VAVE”とすれば、次式(1)のように表され、N個の粒子の過電圧ηjの平均値は、過電圧ΔVMPMに一致する。 Here, in consideration of the peculiar behavior (polarization behavior) of the LiFePO positive electrode, such as the memory effect, path dependency, and voltage hysteresis, when a many-particle model is introduced that assumes that the LiFePO positive electrode is composed (formed) of N spherical particles with different particle radii that are electronically and ionically connected to each other, the OCV of each particle can be obtained from the lithium composition x of each of the particles, as shown in Fig. 2. Furthermore, when an overvoltage occurs due to the reaction of the LiFePO positive electrode, each of the particles that constitute the LiFePO positive electrode has an overvoltage η j according to the lithium composition (where "j" is an integer from 1 to N). Furthermore, if the overvoltage due to the reaction of the LiFePO positive electrode during charging and discharging of a lithium iron phosphate ion battery is "ΔV MPM ," the OCV of the jth particle obtained from Figure 2 is "OCV j ," and the average voltage of the above two spinodal points is "V AVE ," then the overvoltage η j can be expressed as in the following equation (1), and the average value of the overvoltage η j of N particles coincides with the overvoltage ΔV MPM .

Figure 0007616101000001
Figure 0007616101000001

また、過電圧ηjをもったj番目の粒子の電流密度を“ij”とすれば、電流密度ijは、次式(2)のように表され、LiFePO正極に印加される電流の推定値すなわちリン酸鉄リチウムイオン電池の推定電池電流を“Iapp”とすれば、推定電池電流Iappは、次式(3)のように表される。ただし、式(2)において、“F”は、ファラデー定数であり、“k”は、反応の速度定数であり、“cj”はj番目の粒子のリチウムイオン濃度であり、“ce”は、電解液のリチウムイオン濃度(本実施形態では、予め定められた一定値)であり、“cmax”は、活物質中の最大リチウム濃度であり、“R”は、気体定数であり、“T”は、電池温度(電池セルの温度)であり、陽極反応および陰極反応の電荷移動係数は、0.5である。また、式(3)における“Sj”は、j番目の粒子の表面積であり、j番目の粒子の粒子半径を“rj”とすれば、Sj=4πrj 2である。これにより、図2に示すリチウム組成xとLiFePO正極のOCVとの関係と、多粒子モデルとを利用することで、メモリ効果や経路依存性、電圧ヒステリシスといったLiFePO正極の特異な挙動を考慮したリン酸鉄リチウムイオン電池の推定電池電流Iappを算出することができる。 If the current density of the jth particle having an overvoltage η j is "i j ", the current density i j is expressed by the following formula (2), and if the estimated value of the current applied to the LiFePO positive electrode, i.e., the estimated battery current of the lithium iron phosphate ion battery, is "I app ", the estimated battery current I app is expressed by the following formula (3). In formula (2), "F" is the Faraday constant, "k" is the reaction rate constant, "c j " is the lithium ion concentration of the jth particle, "c e " is the lithium ion concentration of the electrolyte (a predetermined constant value in this embodiment), "c max " is the maximum lithium concentration in the active material, "R" is the gas constant, "T" is the battery temperature (temperature of the battery cell), and the charge transfer coefficient of the anode reaction and the cathode reaction is 0.5. Furthermore, "S j " in formula (3) is the surface area of the j-th particle, and if the particle radius of the j-th particle is "r j ", then S j = 4πr j 2. Thus, by utilizing the relationship between the lithium composition x and the OCV of the LiFePO positive electrode shown in Fig. 2 and the multi-particle model, it is possible to calculate the estimated battery current I app of a lithium iron phosphate ion battery taking into account the peculiar behavior of the LiFePO positive electrode, such as the memory effect, path dependency, and voltage hysteresis.

Figure 0007616101000002
Figure 0007616101000002

一方、リン酸鉄リチウムイオン電池の充放電時におけるLiFePO正極の反応による過電圧ΔVMPMは、リン酸鉄リチウムイオン電池(電池セル)の充放電時における電圧“V(t)”とし、リン酸鉄リチウムイオン電池のOCVの前回値を“OCV(t-1)”とし、電解液の塩濃度過電圧を“dVdce”とし、リン酸鉄リチウムイオン電池の充放電時における電流を“I(t)”とし、リン酸鉄リチウムイオン電池の内部抵抗を“Rd”とすれば、次式(4)のように表すことができる。また、式(4)における電圧V(t)および電流I(t)としては、それぞれ実測値(測定値)を用いることが可能であり、塩濃度過電圧dVdceは、次式(5)に基づいて算出することができる。 On the other hand, the overvoltage ΔV MPM due to the reaction of the LiFePO positive electrode when the lithium iron phosphate ion battery is charged and discharged can be expressed as in the following formula (4) if the voltage when the lithium iron phosphate ion battery (battery cell) is charged and discharged is "V(t)", the previous value of the OCV of the lithium iron phosphate ion battery is "OCV(t-1)", the salt concentration overvoltage of the electrolyte is "dV dce ", the current when the lithium iron phosphate ion battery is charged and discharged is "I(t)", and the internal resistance of the lithium iron phosphate ion battery is "Rd". In addition, the voltage V(t) and current I(t) in formula (4) can each be actual measured values (measured values), and the salt concentration overvoltage dV dce can be calculated based on the following formula (5).

Figure 0007616101000003
Figure 0007616101000003

従って、リン酸鉄リチウムイオン電池のOCVの前回値OCV(t-1)を得ることで、式(4)からLiFePO正極の反応による過電圧ΔVMPMを容易に算出することが可能となる。また、上述のように、図2に示すリチウム組成xとLiFePO正極のOCVとの関係と、多粒子モデルとを利用することで、LiFePO正極の特異な挙動を考慮した推定電池電流Iappを算出することができる。これらを踏まえて、本実施形態では、LiFePO正極のOCVに関して当該LiFePO正極の反応(分極)が無視され、図3(a)に示すように、LiFePO正極のOCVは、SOCに拘わらず上記平均電圧VAVE(一定値)になるとみなされる。そして、本実施形態では、当該平均電圧VAVEからリン酸鉄リチウムイオン電池の負極のOCV(図3(b)参照)を減じることにより、図3(c)に示すようなSOCとリン酸鉄リチウムイオン電池のOCVとの関係を規定する第1マップが用意され、当該第1マップからリン酸鉄リチウムイオン電池のOCVの前回値OCV(t-1)が取得される。 Therefore, by obtaining the previous value OCV(t-1) of the OCV of the lithium iron phosphate battery, it becomes possible to easily calculate the overvoltage ΔV MPM due to the reaction of the LiFePO positive electrode from formula (4). As described above, by utilizing the relationship between the lithium composition x and the OCV of the LiFePO positive electrode shown in FIG. 2 and the multi-particle model, it is possible to calculate the estimated battery current I app taking into account the peculiar behavior of the LiFePO positive electrode. In light of this, in this embodiment, the reaction (polarization) of the LiFePO positive electrode is ignored with respect to the OCV of the LiFePO positive electrode, and as shown in FIG. 3(a), the OCV of the LiFePO positive electrode is considered to be the above average voltage V AVE (constant value) regardless of the SOC. In this embodiment, a first map that defines the relationship between the SOC and the OCV of the lithium iron phosphate battery as shown in FIG. 3(c) is prepared by subtracting the OCV of the negative electrode of the lithium iron phosphate battery (see FIG. 3(b)) from the average voltage V AVE, and a previous value OCV(t-1) of the OCV of the lithium iron phosphate battery is obtained from the first map.

なお、本実施形態では、多粒子モデルにおける粒子の総数Nは、例えば200個から2000個であり、粒子半径の平均サイズは、例えば100nmであり、粒子半径は、標準偏差が例えば80nmの対数正規分布で分布するものとした。また、多粒子モデルでは、LiFePO正極を形成する複数の粒子の重心がリン酸鉄リチウムイオン電池のSOCに対応し、複数の粒子の重心に応じたSOCとLiFePO正極のOCVとの関係も、図2に示すものと同様の特徴を有するものとなり、2つのスピノーダル点間の電圧差がおよそ20mVとなり、かつ2つのスピノーダル点の平均電圧がおよそ3.43Vとなる。更に、多粒子モデルでは、SOCごとに粒子の分布状態が複数存在し得ることから、高SOC域および低SOC域では、正極のOCVが20mV程度の幅をもつことになる。従って、多粒子モデルにおいても、リン酸鉄リチウムイオン電池のSOCの変化に対するLiFePO正極のOCVの変化率は高SOC域および低SOC域以外の広い範囲で非常に小さくなり、SOCとOCVとの関係からリン酸鉄リチウムイオン電池のSOCを精度よく推定することは困難となる。 In this embodiment, the total number N of particles in the multi-particle model is, for example, 200 to 2000, the average size of the particle radius is, for example, 100 nm, and the particle radius is distributed in a log-normal distribution with a standard deviation of, for example, 80 nm. In the multi-particle model, the center of gravity of the particles forming the LiFePO positive electrode corresponds to the SOC of the lithium iron phosphate ion battery, and the relationship between the SOC according to the center of gravity of the particles and the OCV of the LiFePO positive electrode also has the same characteristics as those shown in FIG. 2, and the voltage difference between the two spinodal points is approximately 20 mV, and the average voltage of the two spinodal points is approximately 3.43 V. Furthermore, in the multi-particle model, since there can be multiple particle distribution states for each SOC, the OCV of the positive electrode has a width of about 20 mV in the high SOC region and the low SOC region. Therefore, even in the multi-particle model, the rate of change in the OCV of the LiFePO positive electrode relative to changes in the SOC of the lithium iron phosphate ion battery is very small over a wide range other than the high and low SOC regions, making it difficult to accurately estimate the SOC of the lithium iron phosphate ion battery from the relationship between SOC and OCV.

続いて、図4および図5を参照しながら、バッテリ1に含まれる各電池セル2のSOCを推定する手順について詳細に説明する。図4は、バッテリ1を管理するために車両Vに搭載された電子制御装置10(以下、「ECU10」という。)により実行されるSOC推定ルーチンの一例を示すフローチャートである。図4のSOC推定ルーチンは、車両Vがシステム起動されている間にECU10により例えば10-100msおきに繰り返し実行される。 Next, the procedure for estimating the SOC of each battery cell 2 included in the battery 1 will be described in detail with reference to Figures 4 and 5. Figure 4 is a flowchart showing an example of an SOC estimation routine executed by an electronic control unit 10 (hereinafter referred to as "ECU 10") mounted on the vehicle V to manage the battery 1. The SOC estimation routine in Figure 4 is repeatedly executed by the ECU 10, for example, every 10-100 ms while the system of the vehicle V is activated.

図4のSOC推定ルーチンの開始に際して、ECU10は、まず、変数mすなわちバッテリ1に含まれる複数の電池セル2の番号を値1に設定する(ステップS100)。また、ECU10は、図示しない電圧センサにより測定されたm番目の電池セル2(以下、「電池セル2m」という。)の電圧(端子間電圧)の今回(現在)の測定値V(t)、図示しない電流センサにより測定された電池セル2mの電流の今回(現在)の測定値I(t)、塩濃度過電圧dVdce、電池セル2mのOCVの前回値OCV(t-1)、電池セル2mのLiFePO正極が複数の粒子により構成されると仮定したときの各粒子のリチウム濃度の前回値Cj(t-1)、図示しない温度センサにより検出された電池セル2mの温度Tといった当該電池セル2mのSOCの推定に必要な情報を取得する(ステップS110)。塩濃度過電圧dVdceは、上記式(5)に基づいて別途算出されたものであり、電池セル2mのOCVの前回値OCV(t-1)と、各粒子のリチウム濃度の前回値Cj(t-1)は、SOC推定ルーチンの前回実行時に算出されたものである。 4, the ECU 10 first sets the variable m, i.e., the number of the battery cells 2 included in the battery 1, to a value of 1 (step S100). The ECU 10 also acquires information necessary for estimating the SOC of the battery cell 2 m, such as a current (present) measured value V(t) of the voltage (terminal voltage) of the m-th battery cell 2 (hereinafter referred to as "battery cell 2 m ") measured by a voltage sensor (not shown), a current (present) measured value I(t) of the current of the battery cell 2 m measured by a current sensor (not shown), the salt concentration overvoltage dV dce , the previous OCV value OCV(t-1) of the battery cell 2 m , the previous lithium concentration value C j (t-1) of each particle when it is assumed that the LiFePO positive electrode of the battery cell 2 m is composed of multiple particles, and the temperature T of the battery cell 2 m detected by a temperature sensor (not shown) (step S110). The salt concentration overvoltage dV dce is calculated separately based on the above formula (5), and the previous value OCV of the battery cell 2 m , OCV(t−1), and the previous value C j (t−1) of the lithium concentration of each particle are calculated the previous time the SOC estimation routine was executed.

次いで、ECU10は、上記式(4)に従って、ステップS110にて取得した電圧および電流の測定値V(t),I(t)および塩濃度過電圧dVdceに対応したLiFePO正極の反応による過電圧ΔVMPMを算出する(ステップS120)。また、ECU10は、ステップS120にて算出した過電圧ΔVMPMと、ステップS110にて取得したリチウム濃度の前回値Cj(t-1)と、図5に示す第2マップとを用いて、多粒子モデルにおける複数の粒子の各々の過電圧ηjを算出する(ステップS130)。ステップS130において、ECU10は、リチウム濃度の前回値Cj(t-1)に基づいて各粒子のリチウム組成xj(=Cj(t-1)/Cmax)を算出した上で、第2マップからリチウム組成xjに対応した基準電圧を取得する。図5に示す第2マップは、図2に示すリン酸鉄リチウム(LixFePO4)のリチウム組成xとLiFePO正極のOCVとの関係を平均電圧VAVEが0ボルトになるようにオフセットすることにより作成されたものである。すなわち、第2マップから取得される基準電圧は、各粒子のリチウム濃度の前回値Cj(t-1)に基づくリチウム組成xjに対応したOCVjから平均電圧VAVEを減じたものに相当する。更に、ステップS130において、ECU10は、取得した基準電圧を過電圧ΔVMPMから減じることにより多粒子モデルにおける複数の粒子の各々の過電圧ηjを算出する。 Next, the ECU 10 calculates the overvoltage ΔV MPM due to the reaction of the LiFePO positive electrode corresponding to the measured voltage and current values V(t), I(t) and the salt concentration overvoltage dV dce acquired in step S110 according to the above formula (4) (step S120). The ECU 10 also calculates the overvoltage η j of each of the multiple particles in the multi-particle model using the overvoltage ΔV MPM calculated in step S120, the previous value C j (t-1) of the lithium concentration acquired in step S110, and the second map shown in Fig. 5 (step S130). In step S130, the ECU 10 calculates the lithium composition x j (=C j (t-1 ) /C max ) of each particle based on the previous value C j (t-1) of the lithium concentration, and then acquires the reference voltage corresponding to the lithium composition x j from the second map. The second map shown in Fig. 5 was created by offsetting the relationship between the lithium composition x of lithium iron phosphate ( LixFePO4 ) shown in Fig. 2 and the OCV of the LiFePO positive electrode so that the average voltage V AVE becomes 0 volts. In other words, the reference voltage obtained from the second map corresponds to the OCV j corresponding to the lithium composition x j based on the previous value C j (t-1) of the lithium concentration of each particle minus the average voltage V AVE . Furthermore, in step S130, the ECU 10 calculates the overvoltage η j of each of the multiple particles in the multi-particle model by subtracting the obtained reference voltage from the overvoltage ΔV MPM .

ステップS130の処理の後、ECU10は、上記式(2)に従い、ステップS110にて取得した過電圧ηjとリチウム濃度の前回値Cj(t-1)と温度Tとに基づいて多粒子モデルにおける複数の粒子の各々の電流密度ijを算出する(ステップS140)。更に、ECU10は、上記式(3)に従い、ステップS140にて算出した各粒子の電流密度ijと各粒子の表面積Sjとに基づいて電池セル2mの推定電池電流Iappを算出する(ステップS150)。そして、ECU10は、算出した推定電池電流Iappを電池容量Cで除した値をSOC推定ルーチンの前回実行時に算出した電池セル2mのSOCである前回値SOC(t-1)に加算することにより、電池セル2mのSOCの今回値SOC(t)を算出する(ステップS160)。 After the process of step S130, the ECU 10 calculates the current density i j of each of the particles in the multi-particle model based on the overvoltage η j acquired in step S110, the previous lithium concentration value C j (t-1), and the temperature T in accordance with the above formula (2) (step S140). Furthermore, the ECU 10 calculates the estimated battery current I app of the battery cell 2 m in accordance with the current density i j of each particle calculated in step S140 and the surface area S j of each particle in accordance with the above formula (3) (step S150). Then, the ECU 10 calculates the current value SOC(t) of the SOC of the battery cell 2 m by adding a value obtained by dividing the calculated estimated battery current I app by the battery capacity C to the previous value SOC(t -1), which is the SOC of the battery cell 2 m calculated when the SOC estimation routine was executed the previous time (step S160).

また、ECU10は、図3(c)に示す第1マップから電池セル2mのSOCの今回値SOC(t)に対応した電池セル2mのOCVの今回値OCV(t)を取得する(ステップS170)。更に、ECU10は、次式(6)に従い、ステップS140にて算出した電流密度ijに基づいて各粒子のリチウム濃度の変化量∂Cj/∂tを算出すると共に、算出した変化量∂Cj/∂tをSOC推定ルーチンの前回実行時に算出されたリチウム濃度の前回値Cj(t-1)に加算して各粒子のリチウム濃度の今回値Cj(t)を算出する(ステップS180)。ただし、式(6)において、“Vj”は、 j番目の粒子の体積(=4/3・π・rj 3)である。ステップS170にて算出された電池セル2mのOCVの今回値OCV(t)は、SOC推定ルーチンの次回実行時にOCVの前回値OCV(t-1)として用いられ、ステップS180にて算出された各粒子のリチウム濃度の今回値Cj(t)は、SOC推定ルーチンの次回実行時にリチウム濃度の前回値Cj(t-1)として用いられる。 The ECU 10 also obtains the current value OCV(t) of the OCV of the battery cell 2 m corresponding to the current value SOC(t) of the SOC of the battery cell 2 m from the first map shown in Fig. 3(c) (step S170). Furthermore, the ECU 10 calculates the amount of change ∂C j /∂t in the lithium concentration of each particle based on the current density i j calculated in step S140 according to the following formula (6), and calculates the current value C j (t) of the lithium concentration of each particle by adding the calculated amount of change ∂C j /∂t to the previous value C j (t-1) of the lithium concentration calculated at the previous execution of the SOC estimation routine (step S180). Note that in formula (6), "V j " is the volume (= 4/3 x π x r j 3 ) of the j-th particle. The current OCV value OCV(t) of battery cell 2m calculated in step S170 is used as the previous OCV value OCV(t-1) the next time the SOC estimation routine is executed, and the current lithium concentration value C j (t) of each particle calculated in step S180 is used as the previous lithium concentration value C j (t-1) the next time the SOC estimation routine is executed.

Figure 0007616101000004
Figure 0007616101000004

ステップS180にて各粒子のリチウム濃度の今回値Cj(t)を算出した後、ECU10は、変数mをインクリメントし(ステップS190)、変数mがバッテリ1における電池セル2の総数Mを上回っているか否かを判定する(ステップS200)。ECU10は、変数mが総数M以下であると判定した場合(ステップS200:NO)、上記ステップS110-S190の処理を再度実行する。また、変数mが総数Mを上回っていると判定した場合(ステップS200:YES)、ECU10は、その時点で、図4のルーチンを一旦終了させる。 After calculating the current value C j (t) of the lithium concentration of each particle in step S180, the ECU 10 increments the variable m (step S190) and determines whether or not the variable m exceeds the total number M of the battery cells 2 in the battery 1 (step S200). If the ECU 10 determines that the variable m is equal to or smaller than the total number M (step S200: NO), it executes the processes of steps S110-S190 again. If the ECU 10 determines that the variable m exceeds the total number M (step S200: YES), it temporarily ends the routine of FIG. 4 at that point.

以上説明したように、リン酸鉄リチウムイオン電池であるバッテリ1の各電池セル2のSOCの推定に際しては、リン酸鉄リチウム(LiFePO4)により形成されたLiFePO正極がそれぞれリチウムを含む複数の粒子により構成されていると仮定する多粒子モデルが利用される。また、各電池セル2のSOCの推定に際しては、少なくとも、電池セル2mの電圧の測定値V(t)と、リン酸鉄リチウムイオン電池(電池セル2m)のSOCとOCVとの関係を規定するように予め定められた第1マップから取得された電池セル2mのOCVの前回値OCV(t-1)とに基づいて、LiFePO正極の反応による過電圧ΔVMPMが算出される(ステップS120)。更に、LiFePO正極の反応による過電圧ΔVMPMと、図2に示すリン酸鉄リチウム(LixFePO4)のリチウム組成xとLiFePO正極のOCVとの関係に基づいて予め定められた図5に示す第2マップとを用いて、LiFePO正極を構成する複数の粒子の各々の過電圧ηjが算出される(ステップS130)。そして、当該複数の粒子の各々の過電圧ηjから、複数の粒子の各々の電流密度ijが算出され(ステップS140)、当該複数の粒子の各々の電流密度ijから電池セル2mの推定電池電流Iappが算出される(ステップS150)。 As described above, when estimating the SOC of each battery cell 2 of the battery 1, which is a lithium iron phosphate ion battery, a multi-particle model is used, which assumes that the LiFePO positive electrode formed of lithium iron phosphate (LiFePO 4 ) is composed of multiple particles containing lithium. In addition, when estimating the SOC of each battery cell 2, an overvoltage ΔV MPM due to a reaction of the LiFePO positive electrode is calculated (step S120 ) based on at least the measured voltage value V(t) of the battery cell 2 m and the previous OCV value OCV(t-1) of the battery cell 2 m obtained from a first map that is predetermined to define the relationship between the SOC and OCV of the lithium iron phosphate ion battery (battery cell 2 m ). Furthermore, the overvoltage ηj of each of the particles constituting the LiFePO positive electrode is calculated ( step S130 ) using the overvoltage ΔV MPM due to the reaction of the LiFePO positive electrode and a second map shown in Fig. 5 that is determined in advance based on the relationship between the lithium composition x of lithium iron phosphate (LixFePO 4 ) and the OCV of the LiFePO positive electrode shown in Fig. 2. Then, the current density ij of each of the particles is calculated from the overvoltage ηj of each of the particles (step S140), and the estimated battery current I app of the battery cell 2 m is calculated from the current density ij of each of the particles (step S150).

より詳細には、電池セル2mの電圧の測定値V(t)から、少なくとも第1マップから取得されたSOCの前回値SOC(t-1)に対応したOCVの前回値OCV(t-1)を減じて正極の反応による過電圧ΔVMPMが算出される(ステップS120)。また、正極の反応による過電圧ΔVMPMと、第2マップから取得された複数の粒子の各々の基準電圧とから、複数の粒子の各々の過電圧ηjが算出される(ステップS130)。更に、SOCの前回値SOC(t-1)と、推定電池電流Iappと、電池容量Cとから電池セル2mのSOCの今回値SOC(t)が算出(推定)されると共に(ステップS140-S160)、第1マップから今回値SOC(t)に対応したOCVの今回値OCV(t)が取得される(ステップS170)。 More specifically, the overvoltage ΔV MPM due to the reaction of the positive electrode is calculated by subtracting the previous value OCV(t-1) of the OCV corresponding to at least the previous value SOC(t-1) of the SOC obtained from the first map from the measured value V(t) of the voltage of the battery cell 2 m (step S120). In addition, the overvoltage η j of each of the plurality of particles is calculated from the overvoltage ΔV MPM due to the reaction of the positive electrode and the reference voltage of each of the plurality of particles obtained from the second map (step S130). Furthermore, the current value SOC(t) of the SOC of the battery cell 2 m is calculated (estimated) from the previous value SOC(t-1), the estimated battery current I app , and the battery capacity C (steps S140-S160), and the current value OCV(t) of the OCV corresponding to the current value SOC(t) is obtained from the first map (step S170).

これにより、第1マップから取得された電池セル2mのOCVの前回値OCV(t-1)と実際の値とのずれを相殺するように算出された推定電池電流Iappに基づいて電池セル2mのSOCを推定することが可能となる。すなわち、第1マップから取得された前回値OCV(t-1)が実際のOCVよりも小さい場合、ステップS120にて算出される過電圧ΔVMPMが大きくなることで推定電池電流IappがOCVのずれを相殺するように大きく算出され、第1マップから取得された前回値OCV(t-1)が実際のOCVよりも大きい場合、ステップS120にて算出される過電圧ΔVMPMが小さくなることで推定電池電流IappがOCVのずれを相殺するように小さく算出される。従って、電池セル2mのSOCの推定が繰り返し実行されることで当該SOCの推定値すなわち今回値SOC(t)を実際の値に収束させていくことができる。この結果、メモリ効果や経路依存性、電圧ヒステリシスといった特異な挙動を示すリン酸鉄リチウム(LiFePO4)により形成されたLiFePO正極を含む各電池セル2のSOCを精度よく推定することが可能となる。 This makes it possible to estimate the SOC of the battery cell 2 m based on the estimated battery current I app calculated to offset the difference between the previous value OCV(t-1) of the OCV of the battery cell 2 m acquired from the first map and the actual value. That is, when the previous value OCV(t-1) acquired from the first map is smaller than the actual OCV, the overvoltage ΔV MPM calculated in step S120 is increased, so that the estimated battery current I app is calculated to be large so as to offset the difference in the OCV, and when the previous value OCV(t-1) acquired from the first map is larger than the actual OCV, the overvoltage ΔV MPM calculated in step S120 is decreased, so that the estimated battery current I app is calculated to be small so as to offset the difference in the OCV. Thus, by repeatedly estimating the SOC of the battery cell 2 m , the estimated value of the SOC, i.e., the current value SOC(t), can be converged to the actual value. As a result, it becomes possible to accurately estimate the SOC of each battery cell 2 including a LiFePO positive electrode formed from lithium iron phosphate (LiFePO 4 ), which exhibits peculiar behaviors such as memory effect, path dependency, and voltage hysteresis.

更に、上記実施形態において、第1マップは、SOCごとに、LiFePO正極の2つのスピノーダル点の平均電圧VAVEからリン酸鉄リチウムイオン電池の負極のOCVを減じることにより作成され、第2マップは、リン酸鉄リチウム(LixFePO4)のリチウム組成xと正極のOCVとの関係を平均電圧VAVEが0ボルトになるようにオフセットすることにより作成される。これにより、正極の反応による過電圧ΔVMPMと、第2マップとを用いて、複数の粒子の各々の過電圧ηjをより適正に算出することが可能となる。 Furthermore, in the above embodiment, the first map is created by subtracting the OCV of the negative electrode of the lithium iron phosphate battery from the average voltage V AVE of the two spinodal points of the LiFePO positive electrode for each SOC, and the second map is created by offsetting the relationship between the lithium composition x of lithium iron phosphate (LixFePO 4 ) and the OCV of the positive electrode so that the average voltage V AVE becomes 0 V. This makes it possible to more appropriately calculate the overvoltage η j of each of the multiple particles by using the overvoltage ΔV MPM due to the reaction of the positive electrode and the second map.

また、上記実施形態では、複数の粒子の各々の電流密度ijに基づいて各粒子のリチウム濃度の変化量∂Cj/∂tが算出されると共に、リチウム濃度の前回値Cj(t-1)と変化量∂Cj/∂tとからリチウム濃度の今回値Cj(t)が算出され(ステップS180)、図4のステップS130では、第2マップからリチウム濃度の前回値Cj(t-1)に基づくリチウム組成xjに対応した複数の粒子の各々の基準電圧が取得される。これにより、複数の粒子の各々の基準電圧ひいては過電圧ηjをより適正に算出することが可能となる。 In the above embodiment, the change amount ∂C j /∂t in the lithium concentration of each particle is calculated based on the current density i j of each of the plurality of particles, and the current value C j (t) of the lithium concentration is calculated from the previous value C j (t-1) of the lithium concentration and the change amount ∂C j /∂t (step S180), and in step S130 of Fig. 4, the reference voltage of each of the plurality of particles corresponding to the lithium composition x j based on the previous value C j (t-1) of the lithium concentration is obtained from the second map. This makes it possible to more appropriately calculate the reference voltage of each of the plurality of particles, and therefore the overvoltage η j .

更に、図2のステップS120では、電池セル2mの電圧の測定値V(t)から、OCVの前回値OCV(t-1)と、塩濃度過電圧dVdceと、電池セル2mの電流の測定値I(t)と内部抵抗Rdとの積値とを減じて、LiFePO正極の反応による過電圧ΔVMPMが算出される。これにより、演算負荷を低減しつつ、正極の反応による過電圧ΔVMPMをより適正に算出することが可能となる。ただし、図2のステップS120では、電池セル2mの電流の測定値I(t)を用いることなく、電池セル2mの電圧の測定値V(t)から、OCVの前回値OCV(t-1)と塩濃度過電圧dVdceとを減じて正極の反応による過電圧ΔVMPMと内部抵抗Rdによる過電圧との和ΔVMPM+IRが算出されてもよい。そして、この場合には、|ΔVMPM+IR-(ΔVMPM+Iapp・Rd)|を最小にする過電圧ΔVMPMと推定電池電流Iappとがイタレーションにより導出されてもよい。 2, the overvoltage ΔV MPM due to the reaction of the LiFePO positive electrode is calculated by subtracting the previous OCV value OCV(t-1), the salt concentration overvoltage dV dce , and the product of the measured current I(t) and the internal resistance Rd of the battery cell 2 m from the measured voltage V(t) of the battery cell 2 m . This makes it possible to more appropriately calculate the overvoltage ΔV MPM due to the reaction of the positive electrode while reducing the calculation load. However, in step S120 of FIG. 2, the sum ΔV MPM +IR of the overvoltage ΔV MPM due to the reaction of the positive electrode and the overvoltage due to the internal resistance Rd may be calculated by subtracting the previous OCV value OCV(t-1) and the salt concentration overvoltage dV dce from the measured voltage V( t ) of the battery cell 2 m without using the measured current I (t). In this case, the overvoltage ΔV MPM and the estimated battery current I app that minimize |ΔV MPM+IR - (ΔV MPM +I app · Rd)| may be derived by iteration.

なお、リン酸鉄リチウムイオン電池である電池セル2を複数含むバッテリ1は、モータジェネレータMGと電力をやり取りする高電圧バッテリに限られるものではない。すなわち、バッテリ1は、補機バッテリ(12Vバッテリ)といった低電圧バッテリとして構成されてもよい。 The battery 1, which includes multiple battery cells 2 that are lithium iron phosphate batteries, is not limited to being a high-voltage battery that exchanges power with the motor generator MG. In other words, the battery 1 may be configured as a low-voltage battery such as an auxiliary battery (12 V battery).

また、本開示の発明は上記実施形態に何ら限定されるものではなく、本開示の外延の範囲内において様々な変更をなし得ることはいうまでもない。更に、上記実施形態は、あくまで発明の概要の欄に記載された発明の具体的な一形態に過ぎず、発明の概要の欄に記載された発明の要素を限定するものではない。 Furthermore, the invention of this disclosure is in no way limited to the above embodiment, and it goes without saying that various modifications can be made within the scope of the present disclosure. Furthermore, the above embodiment is merely one specific form of the invention described in the Summary of the Invention section, and does not limit the elements of the invention described in the Summary of the Invention section.

本開示の発明は、リン酸鉄リチウムイオン電池の製造分野等において利用可能である。 The disclosed invention can be used in fields such as the manufacturing of lithium iron phosphate batteries.

1 バッテリ、2 電池セル、10 電子制御装置(ECU)、MG モータジェネレータ、V 車両。 1 battery, 2 battery cell, 10 electronic control unit (ECU), MG motor generator, V vehicle.

Claims (1)

リン酸鉄リチウムにより形成された正極を含むリン酸鉄リチウムイオン電池のSOCを推定するSOC推定方法であって、
少なくとも、前記リン酸鉄リチウムイオン電池の電圧の測定値と、前記リン酸鉄リチウムイオン電池のSOCとOCVとの関係を規定するように予め定められた第1マップから取得されたOCVとに基づいて、前記正極の反応による過電圧を算出し、
前記正極の反応による過電圧と、前記リン酸鉄リチウムのリチウム組成と前記正極のOCVとの関係に基づいて予め定められた第2マップとを用いて、前記正極がそれぞれリチウムを含む複数の粒子により構成されていると仮定したときの各粒子の過電圧を算出し、
前記各粒子の過電圧に基づいて前記各粒子の電流密度を算出し、
前記各粒子の電流密度に基づいて前記リン酸鉄リチウムイオン電池の推定電池電流を算出し、
前記推定電池電流に基づいて前記リン酸鉄リチウムイオン電池のSOCを推定するSOC推定方法。
1. A method for estimating an SOC of a lithium iron phosphate ion battery including a positive electrode formed of lithium iron phosphate, comprising:
calculating an overvoltage caused by a reaction of the positive electrode based on at least a measured value of the voltage of the lithium iron phosphate battery and an OCV obtained from a first map that is predetermined to define a relationship between an SOC and an OCV of the lithium iron phosphate battery;
calculating an overvoltage of each particle on the assumption that the positive electrode is composed of a plurality of particles each containing lithium, using an overvoltage due to a reaction of the positive electrode and a second map that is determined in advance based on a relationship between the lithium composition of the lithium iron phosphate and an OCV of the positive electrode;
Calculating a current density of each of the particles based on the overpotential of each of the particles;
Calculating an estimated battery current of the lithium iron phosphate battery based on the current density of each of the particles;
An SOC estimation method for estimating an SOC of the lithium iron phosphate battery based on the estimated battery current.
JP2022006288A 2022-01-19 2022-01-19 SOC Estimation Method Active JP7616101B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022006288A JP7616101B2 (en) 2022-01-19 2022-01-19 SOC Estimation Method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022006288A JP7616101B2 (en) 2022-01-19 2022-01-19 SOC Estimation Method

Publications (2)

Publication Number Publication Date
JP2023105454A JP2023105454A (en) 2023-07-31
JP7616101B2 true JP7616101B2 (en) 2025-01-17

Family

ID=87468893

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022006288A Active JP7616101B2 (en) 2022-01-19 2022-01-19 SOC Estimation Method

Country Status (1)

Country Link
JP (1) JP7616101B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020069137A (en) * 2018-10-31 2020-05-07 株式会社三洋物産 Game machine
JP2020069143A (en) * 2018-10-31 2020-05-07 株式会社三洋物産 Game machine
JP2020069140A (en) * 2018-10-31 2020-05-07 株式会社三洋物産 Game machine
JP2020069141A (en) * 2018-10-31 2020-05-07 株式会社三洋物産 Game machine
JP2020146312A (en) * 2019-03-14 2020-09-17 株式会社三洋物産 Game machine
JP2020146311A (en) * 2019-03-14 2020-09-17 株式会社三洋物産 Game machine

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015019482A1 (en) 2013-08-09 2015-02-12 株式会社日立製作所 Positive electrode material for lithium ion secondary batteries
JP2019144039A (en) 2018-02-19 2019-08-29 旭化成株式会社 Method, device and system for estimating soc of nonaqueous lithium type power storage element

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5737106B2 (en) * 2011-09-26 2015-06-17 トヨタ自動車株式会社 Secondary battery state estimation device
JP2018063191A (en) * 2016-10-13 2018-04-19 カルソニックカンセイ株式会社 Charge state estimation device
JP6969307B2 (en) * 2017-11-14 2021-11-24 株式会社Gsユアサ Management device, power storage system, method of equalizing the remaining capacity of the power storage element, method of estimating the internal state of the power storage element

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015019482A1 (en) 2013-08-09 2015-02-12 株式会社日立製作所 Positive electrode material for lithium ion secondary batteries
JP2019144039A (en) 2018-02-19 2019-08-29 旭化成株式会社 Method, device and system for estimating soc of nonaqueous lithium type power storage element

Also Published As

Publication number Publication date
JP2023105454A (en) 2023-07-31

Similar Documents

Publication Publication Date Title
JP7616101B2 (en) SOC Estimation Method
US12237477B2 (en) State of charge dependent plating estimation and prevention
US9728992B2 (en) Control apparatus and control method for secondary battery
CN102239597B (en) Method for determining charge state of secondary intercalation cell of rechargeable battery
US9399404B2 (en) Charging system for all-solid-state battery
US6700383B2 (en) Method of detecting and resolving memory effect
US9475480B2 (en) Battery charge/discharge control device and hybrid vehicle using the same
CN107179508B (en) Battery state of charge estimation based on reduced order electrochemical model
CN107179509B (en) Battery capacity estimation based on open-loop and closed-loop models
US20160023566A1 (en) Reduced order electrochemical battery model for vehicle control
US20160023567A1 (en) Temperature dependent electrochemical battery model for vehicle control
CN109643775B (en) Dual energy storage system and starter battery module
US10035426B2 (en) Battery power management in hybrid vehicles
US20080135315A1 (en) Fast search algorithm for finding initial diffusion voltage in electro-chemical systems
CN102771003B (en) Battery control system
JP4488426B2 (en) Storage device control device
US20160023568A1 (en) Interpolation of metal-ion concentrations in a battery model for vehicle control
US20080135316A1 (en) Method and Apparatus for Modeling Diffusion in an Electrochemical System
Sinha et al. Hybrid energy storage system for the life extension of lithium-ion batteries in electric vehicles
Mahmud et al. The impact of battery operating temperature and state of charge on the lithium-ion battery internal resistance
CN116472634A (en) Battery control device and battery system
Jarraya et al. Modeling and parameters estimation for lithium-ion cells in electric drive vehicle
Marcicki et al. A reduced-order electrochemical model of lithium-ion cells for system identification of battery aging
Samba et al. Statistical Sensitivity Analysis of Li-ion Pouch Battery Cell Dimension and Design

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240320

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20241125

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20241203

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241216

R150 Certificate of patent or registration of utility model

Ref document number: 7616101

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150