JP7616101B2 - SOC Estimation Method - Google Patents
SOC Estimation Method Download PDFInfo
- 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
Links
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy 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.
ここで、三元系リチウムイオン電池等においては、一般に、電流の測定値を積算して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
そこで、本開示は、リン酸鉄リチウムにより形成された正極を含むリン酸鉄リチウムイオン電池の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.
次に、図面を参照しながら、本開示の発明を実施するための形態について説明する。 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
バッテリ1は、例えば直列に接続される複数の電池モジュールを含む電池スタックと、当該電池スタックを収容するケースとを含む、いわゆる高電圧バッテリであり、各電池モジュールは、例えば直列に接続されると共にモジュールケース内に収容される複数の電池セル2を含むものである。電池モジュールを構成する電池セル2は、本実施形態では、LiFePO4であるオリビン型の結晶構造を有するリン酸鉄リチウムにより形成された正極(LiFePO正極)と、黒鉛系炭素材料等により形成された負極とを含むリン酸鉄リチウムイオン電池であり、正極および負極は、セパレータや有機溶媒等である電解液と共に外装体の内部に収容される。
The
図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 .
また、過電圧η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.
一方、リン酸鉄リチウムイオン電池の充放電時における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).
従って、リン酸鉄リチウムイオン電池の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
図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
次いで、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
ステップ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
また、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
ステップ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
以上説明したように、リン酸鉄リチウムイオン電池であるバッテリ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
より詳細には、電池セル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
これにより、第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
更に、上記実施形態において、第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
なお、リン酸鉄リチウムイオン電池である電池セル2を複数含むバッテリ1は、モータジェネレータMGと電力をやり取りする高電圧バッテリに限られるものではない。すなわち、バッテリ1は、補機バッテリ(12Vバッテリ)といった低電圧バッテリとして構成されてもよい。
The
また、本開示の発明は上記実施形態に何ら限定されるものではなく、本開示の外延の範囲内において様々な変更をなし得ることはいうまでもない。更に、上記実施形態は、あくまで発明の概要の欄に記載された発明の具体的な一形態に過ぎず、発明の概要の欄に記載された発明の要素を限定するものではない。 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と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.
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)
| 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)
| 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)
| 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 |
-
2022
- 2022-01-19 JP JP2022006288A patent/JP7616101B2/en active Active
Patent Citations (2)
| 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 |