JP2016004726A - 電池システム - Google Patents

Abstract

【課題】 リチウム析出量を推定するとき、リチウムイオン二次電池の負極における反応抵抗を考慮する必要がある。
【解決手段】 析出反応および溶解反応における過電圧および交換電流密度から算出される反応電流に基づいて、活性リチウムの量を算出する。不活性化反応における過電圧および交換電流密度から算出される反応電流に基づいて、不活性リチウムの量を算出する。リチウム析出量は、活性リチウムの量および不活性リチウムの量との総量となる。析出反応および溶解反応における交換電流密度は、反応抵抗変化率に基づいて補正される。反応抵抗変化率は、リチウムイオン二次電池に与えられる拘束荷重およびリチウムイオン二次電池の温度のうちの少なくとも一方から算出する。反応抵抗変化率を算出すれば、リチウムイオン二次電池における現在の反応抵抗を把握できる。これにより、反応抵抗に依存する活性リチウムの量を把握した上で、リチウム析出量を算出できる。
【選択図】 図１

Description

本発明は、リチウムイオン二次電池におけるリチウム析出量を算出する技術に関する。

リチウムイオン二次電池では、負極板の表面に金属リチウムが析出することが知られている。そこで、従来、リチウム析出量を推定する技術が提案されている。例えば、特許文献１では、二次電池の端子間電圧を段階的に低下させたときの電流値に基づいて、リチウム析出量を算出している。

特開２０１０−０８６９０１号公報

金属リチウムの析出は、リチウムイオン二次電池の負極電位だけでなく、金属リチウムの析出反応における反応抵抗にも依存する。この反応抵抗が増加するほど、金属リチウムが析出しやすくなる。そこで、金属リチウムの析出量を算出するときには、析出反応における反応抵抗を考慮する必要がある。

本発明の電池システムは、充放電を行うリチウムイオン二次電池と、リチウム析出量を算出するコントローラと、を有する。リチウム析出量とは、活性リチウムを析出させる析出反応および、活性リチウムを溶解させる溶解反応が行われたときの活性リチウムの量と、活性リチウムの不活性化反応によって生成された不活性リチウムの量との総量を示す。

コントローラは、析出反応および溶解反応における過電圧および交換電流密度から算出される反応電流に基づいて、活性リチウムの量を算出する。また、コントローラは、不活性化反応における過電圧および交換電流密度から算出される反応電流に基づいて、不活性リチウムの量を算出する。

ここで、析出反応および溶解反応における交換電流密度は、反応抵抗変化率に基づいて補正される。反応抵抗変化率とは、リチウムイオン二次電池の負極における反応抵抗に関して、基準値および現在値の比率を示す。反応抵抗変化率と、リチウムイオン二次電池に与えられる拘束荷重およびリチウムイオン二次電池の温度のうちの少なくとも一方であるパラメータとの対応関係を用いれば、パラメータを取得することにより、取得したパラメータに応じた反応抵抗変化率を算出できる。

本発明によれば、リチウムイオン二次電池における現在の反応抵抗変化率を考慮して、活性リチウムの量を算出している。反応抵抗変化率を算出すれば、リチウムイオン二次電池における現在の反応抵抗を把握できる。これにより、反応抵抗に依存する活性リチウムの量を把握した上で、リチウム析出量を算出できる。

電池システムの構成を示す図である。 析出反応、溶解反応および不活性化反応を示す図である。 反応抵抗変化率、拘束荷重および電池温度の対応関係を示す図である。 リチウム析出量を算出する処理を示すフローチャートである。

以下、本発明の実施例について説明する。

実施例１である電池システムについて、図１を用いて説明する。二次電池１０は、正極ラインＰＬおよび負極ラインＮＬを介して負荷２０に接続されている。二次電池１０としては、リチウムイオン二次電池が用いられる。二次電池１０には、二次電池１０を挟む方向において、所定の拘束荷重が与えられている。具体的には、二次電池１０を挟む位置に一対のエンドプレーを配置し、これらのエンドプレートを互いに連結することにより、二次電池１０に拘束荷重を与えることができる。

正極ラインＰＬには、システムメインリレーＳＭＲ−Ｂが設けられ、負極ラインＮＬには、システムメインリレーＳＭＲ−Ｇが設けられている。システムメインリレーＳＭＲ−Ｂ，ＳＭＲ−Ｇは、コントローラ４０からの制御信号を受けて、オンおよびオフの間で切り替わる。

負荷２０は、二次電池１０の放電電力を受けて動作したり、二次電池１０に電力（充電電力）を供給したりする。本実施例の電池システムを車両に搭載したとき、負荷２０としては、モータ・ジェネレータを用いることができる。モータ・ジェネレータは、二次電池１０の放電電力を受けて、車両を走行させるための運動エネルギを生成する。また、モータ・ジェネレータは、車両の制動時に発生する運動エネルギを電力に変換し、この電力（回生電力）を二次電池１０に出力することができる。なお、二次電池１０を車両に搭載するときには、複数の二次電池１０を直列に接続することによって構成された組電池を車両に搭載することができる。

電圧センサ３１は、二次電池１０の電圧値Ｖｂを検出し、検出結果をコントローラ４０に出力する。電流センサ３２は、二次電池１０の電流値Ｉｂを検出し、検出結果をコントローラ４０に出力する。本実施例において、二次電池１０を放電しているときの電流値Ｉｂを正の値とし、二次電池１０を充電しているときの電流値Ｉｂを負の値とする。温度センサ３３は、二次電池１０の温度（電池温度）Ｔｂを検出し、検出結果をコントローラ４０に出力する。コントローラ４０は、メモリ４１を有する。メモリ４１は、所定の情報を記憶している。

二次電池１０を充放電したとき、負極板および正極板の間において、リチウムイオンが移動する。具体的には、二次電池１０を放電したとき、正極板では、リチウムイオンおよび電子を吸収する反応が行われ、負極板では、リチウムイオンおよび電子を放出する反応が行われる。一方、二次電池１０を充電したとき、正極板では、リチウムイオンおよび電子を放出する反応が行われ、負極板では、リチウムイオンおよび電子を吸収する反応が行われる。

一方、負極板では、図２に示す副反応が行われる。析出反応（副反応）では、リチウムイオンおよび電子が反応することにより、活性リチウムが生成される。この析出反応は、上述した二次電池１０の充放電に応じて、リチウムイオンが放出されたり、リチウムイオンが吸収されたりする反応とは異なる。溶解反応（副反応）は、活性リチウムがリチウムイオンおよび電子に分離する反応であり、溶解反応の発生は、二次電池１０の負極電位に依存する。

不活性化反応（副反応）は、活性リチウムが電解液と反応して不活性リチウム（リチウム化合物）を生成する反応（不可逆反応）である。不活性リチウムは、負極活物質層において、被膜として形成され、二次電池１０の充放電に関与しなくなる。このため、不活性リチウムの量が増えるほど、二次電池１０の満充電容量が低下してしまう。なお、活性リチウムも二次電池１０の充放電に関与しなくなるが、溶解反応が行われれば、活性リチウムをリチウムイオンに戻して充放電を行うことができる。

析出反応および溶解反応が行われたときの活性リチウムの量と、不活性化反応によって生成された不活性リチウムの量との総量を、リチウム析出量Ｑ_Ｌｉとして規定する。析出反応および溶解反応が行われたときの活性リチウムの量とは、析出反応によって生成された活性リチウムであって、溶解反応が行われた後に残存している活性リチウムの量である。

リチウム析出量Ｑ_Ｌｉを算出する方法について説明する。本実施例において、析出反応および溶解反応が行われたときの活性リチウムの量を、活性リチウムの電荷量Ｑ_３（単位［Ａｈ］）として規定する。また、不活性化反応によって生成された不活性リチウムの量を、不活性リチウムの電荷量Ｑ_４（単位［Ａｈ］）として規定する。なお、以下の各式（１）〜（３１）で説明した記号と重複する記号については、繰り返して説明しない。

リチウム析出量Ｑ_Ｌｉは、下記式（１）に基づいて算出される。また、活性リチウムの量Ｑ_３は、下記式（２）に基づいて算出され、不活性リチウムの量Ｑ_４は、下記式（３）に基づいて算出される。

活性リチウムの量Ｑ_３は、析出反応および溶解反応に依存する。また、活性リチウムの量Ｑ_３を電荷量として規定することにより、活性リチウムの量Ｑ_３は、析出反応および溶解反応の両方向を考慮した電流値を、初期状態（０秒）から現在の時間ｔ（単位［秒］）までの間で積算した値となる。初期状態とは、二次電池１０を初めて使用し始めた状態である。この電流値は、析出反応および溶解反応における反応電流密度Ｉ_３’（単位［Ａ／ｃｍ^２］）に対して、電極の反応面積Ｓ（単位［ｃｍ^２］）を乗算することによって算出される。したがって、活性リチウムの量Ｑ_３は、上記式（２）によって表すことができる。反応面積Ｓとは、二次電池１０において、正極板（正極活物質層）および負極板（負極活物質層）が対向する領域における電極板（正極板や負極板）上の面積であり、予め求めておくことができる。

不活性リチウムの量Ｑ_４は、不活性化反応に依存する。また、不活性リチウムの量Ｑ_４を電荷量として規定することにより、不活性リチウムの量Ｑ_４は、不活性化反応における電流値を、上述した初期状態（０秒）から現在の時間ｔ（単位［秒］）までの間で積算した値となる。この電流値は、不活性化反応における反応電流密度Ｉ_４’ （単位［Ａ／ｃｍ^２］）に対して、電極の反応面積Ｓを乗算することによって算出される。したがって、不活性リチウムの量Ｑ_４は、上記式（３）によって表すことができる。

反応電流密度Ｉ_３’，Ｉ_４’は、下記式（４），（５）に示すように規定することができる。反応電流密度Ｉ_３’，Ｉ_４’は、後述するように、所定時間（演算周期）Δｔの間で算出される。

上記式（４）において、ｉ_３’は、析出反応および溶解反応において、二次電池１０の負極の界面における反応電流密度（単位［Ａ／ｃｍ^２］）である。Ａ_３は、析出反応および溶解反応において、単位電極反応面積当たりの反応表面積を示す面積比（単位［ｃｍ^２／ｃｍ^２］）である。すなわち、面積比Ａ_３は、反応表面積を反応面積Ｓで除算した値である。反応表面積とは、負極板の表面のうち、活性リチウムが析出する領域の総面積である。

上記式（５）において、ｉ_４’は、不活性化反応において、二次電池１０の負極の界面における反応電流密度（単位［Ａ／ｃｍ^２］）である。Ａ_４は、単位電極反応面積当たりの反応表面積を示す面積比（単位［ｃｍ^２／ｃｍ^２］）である。すなわち、面積比Ａ_４は、反応表面積を反応面積Ｓで除算した値である。ここで、面積比Ａ_４は、面積比Ａ_３に等しいとみなすことができる。

反応電流密度ｉ_３’および面積比Ａ_３を算出すれば、活性リチウムの量Ｑ_３を算出でき、反応電流密度ｉ_４’および面積比Ａ_４（面積比Ａ_３）を算出すれば、不活性リチウムの量Ｑ_４を算出できる。そして、算出した量Ｑ_３，Ｑ_４に基づいて、リチウム析出量Ｑ_Ｌｉを算出できる。

析出反応および溶解反応における反応電流密度ｉ_３’を算出する方法について、以下に説明する。

図２に示すように、析出反応および溶解反応は、互いに逆方向の反応となる。析出反応および溶解反応の両者を考慮したときの反応電流密度ｉ_３’は、バトラー・ボルマー式に基づいて下記式（６）によって表すことができる。下記式（６）に示す右辺第１項は、析出反応における反応電流密度を示し、下記式（６）に示す右辺第２項は、溶解反応における反応電流密度を示す。

上記式（６）において、ｉ_０３’は、単位反応表面積当たりの交換電流密度（単位［Ａ／ｃｍ^２］）であり、α_ａ３は酸化反応（すなわち、析出反応）の移動係数であり、α_ｃ３は還元反応（すなわち、溶解反応）の移動係数である。移動係数α_ａ３，α_ｃ３は、予め求めておくことができる。Ｆはファラデー定数、Ｒはガス定数、Ｔは温度（すなわち、電池温度Ｔｂ）、η_３’は反応過電圧である。

反応過電圧η_３’は、負極板の界面における負極電位Ｖ_２’と、析出反応および溶解反応における平衡電位（定数）Ｕ_３との差になるため、下記式（７）によって表すことができる。下記式（７）では、負極板の表面に形成された被膜の抵抗値（被膜抵抗という）Ｒ_ｆ３’に伴う電圧降下量も考慮している。この被膜は、二次電池１０を初めて充電したときに負極板（負極活物質層）の表面に形成されるＳＥＩ（Solid Electrolyte Interface）被膜である。

上記式（７）に示す反応過電圧η_３’を上記式（６）に代入することにより、反応電流密度ｉ_３’を算出することができる。

上記式（６）に示す交換電流密度ｉ_０３’は、下記式（８）に基づいて算出される。下記式（８）に示すｋ_ｒａｔｅは、析出反応および溶解反応における反応抵抗変化率である。また、下記式（８）に示す交換電流密度ｉ_０３は、下記式（９）に基づいて算出される。交換電流密度ｉ_０３は、負極板（負極活物質層）に含まれる電解液中の塩濃度ｃ_ｅ２および温度Ｔ（すなわち、電池温度Ｔｂ）に依存するため、下記式（９）によって表すことができる。後述するように塩濃度ｃ_ｅ２を算出し、温度Ｔ（電池温度Ｔｂ）を検出することにより、下記式（９）に基づいて、交換電流密度ｉ_０３を算出できる。

上記式（９）において、塩濃度ｃ_{ｅ，ｉｎｉ}は、電解液中の平均塩濃度の初期値である。Ｅ_ｉ０３は、温度依存性を表す活性化エネルギである。参照温度Ｔ_ｒｅｆにおける交換電流密度ｉ_０３を予め求めておき、温度Ｔの依存性も考慮して、アレニウス式を適用すると、交換電流密度ｉ_０３は、上記式（９）によって表すことができる。上記式（９）に示す交換電流密度ｉ_０３（Ｔ_ｒｅｆ）は、参照温度（所定値）Ｔ_ｒｅｆにおいて、交流インピーダンス測定から同定される。

本実施例では、上記式（８）から分かるように、上記式（９）から算出された交換電流密度ｉ_０３を反応抵抗変化率ｋ_ｒａｔｅによって補正している。反応抵抗変化率ｋ_ｒａｔｅとは、予め設定した基準値（固定値）としての反応抵抗Ｒｒ＿ｒｅｆと、現在値としての反応抵抗Ｒｒ＿ｃｕｒとの比率（Ｒｒ＿ｃｕｒ／Ｒｒ＿ｒｅｆ）を示す。反応抵抗変化率ｋ_ｒａｔｅを算出すれば、二次電池１０における現在の反応抵抗Ｒｒ＿ｃｕｒを把握できる。

反応抵抗変化率ｋ_ｒａｔｅは、二次電池１０が受ける拘束荷重や電池温度Ｔｂに依存する。このため、反応抵抗変化率ｋ_ｒａｔｅと、拘束荷重および電池温度Ｔｂの少なくとも一方との対応関係（マップ又は演算式）を実験に基づいて予め求めておくことができる。この対応関係を用いれば、拘束荷重や電池温度Ｔｂを特定（取得）することにより、反応抵抗変化率ｋ_ｒａｔｅを算出することができる。すなわち、反応抵抗変化率ｋ_ｒａｔｅを算出するときには、拘束荷重および電池温度Ｔｂの少なくとも一方を考慮することができる。

図３には、反応抵抗変化率ｋ_ｒａｔｅ、拘束荷重および電池温度Ｔｂの関係（一例）を示す。図３において、縦軸は反応抵抗変化率ｋ_ｒａｔｅを示し、横軸は、二次電池１０に作用する拘束荷重を示す。図３に示す実線および一点鎖線は、互いに異なる２つの電池温度Ｔｂにおいて、反応抵抗変化率ｋ_ｒａｔｅおよび拘束荷重の対応関係を示している。図３に示す高温および低温は、２つの電池温度Ｔｂを比較したときに、高い側の電池温度Ｔｂと、低い側の電池温度Ｔｂとをそれぞれ示している。

図３に示すように、所定の拘束荷重Ｌｄを境界として、反応抵抗変化率ｋ_ｒａｔｅおよび拘束荷重の対応関係が異なっている。具体的には、図３の実線で示す対応関係では、所定の拘束荷重Ｌｄよりも高い範囲では、拘束荷重が変化しても、反応抵抗変化率ｋ_ｒａｔｅは変化しにくくなる。一方、所定の拘束荷重Ｌｄ以下の範囲では、拘束荷重の変化に応じて、反応抵抗変化率ｋ_ｒａｔｅが変化しやすくなる。具体的には、所定の拘束荷重Ｌｄ以下の範囲内において、拘束荷重が低いほど、反応抵抗変化率ｋ_ｒａｔｅが高くなる。言い換えれば、拘束荷重が高いほど、反応抵抗変化率ｋ_ｒａｔｅが低くなる。

図３の一点鎖線で示す対応関係でも、所定の拘束荷重Ｌｄ以下の範囲では、所定の拘束荷重Ｌｄよりも高い範囲に比べて、拘束荷重の変化に応じて、反応抵抗変化率ｋ_ｒａｔｅが変化しやすくなる。一方、図３に示すように、電池温度Ｔｂが異なっていると、拘束荷重が等しくても、反応抵抗変化率ｋ_ｒａｔｅが異なる。

二次電池１０の拘束荷重を測定するためには、圧力センサを用いることができる。具体的には、二次電池１０の外装を構成するケース（内壁面又は外壁面）に対して圧力センサを配置すれば、二次電池１０に作用する拘束荷重を測定できる。一方、二次電池１０のＳＯＣや電池温度Ｔｂに応じて拘束荷重が変化するため、ＳＯＣや電池温度Ｔｂに基づいて、拘束荷重を推定することもできる。

具体的には、ＳＯＣおよび電池温度Ｔｂの少なくとも一方と、拘束荷重との対応関係（マップ又は演算式）を実験に基づいて予め求めておけば、ＳＯＣや電池温度Ｔｂを特定することにより、拘束荷重を算出（推定）することができる。すなわち、拘束荷重を推定するときには、ＳＯＣおよび電池温度Ｔｂの少なくとも一方を考慮すればよい。なお、二次電池１０のＳＯＣを推定する方法としては、公知の方法を適宜採用することができる。

上記式（７）に示す被膜抵抗Ｒ_ｆ３’は、析出反応における被膜抵抗Ｒ_{ｆ３ｄｅｐｏ}と、溶解反応における被膜抵抗Ｒ_{ｆ３ｄｉｓｓ}とに分けられる。析出反応では、反応過電圧η_３’が０よりも小さくなるため、反応過電圧η_３’が０よりも小さいときには、被膜抵抗Ｒ_ｆ３’として被膜抵抗Ｒ_{ｆ３ｄｅｐｏ}が用いられる。溶解反応では、反応過電圧η_３’が０以上となるため、反応過電圧η_３’が０以上であるときには、被膜抵抗Ｒ_ｆ３’として被膜抵抗Ｒ_{ｆ３ｄｉｓｓ}が用いられる。このため、析出反応においては、下記式（１０）に基づいて被膜抵抗Ｒ_ｆ３’を算出できる。また、溶解反応においては、下記式（１１）に基づいて被膜抵抗Ｒ_ｆ３’を算出できる。

析出反応における被膜抵抗Ｒ_{ｆ３ｄｅｐｏ}は、負極板の被膜におけるイオン伝導の抵抗値（被膜抵抗という）Ｒ_ｆｉ２に相当する。イオン伝導とは、被膜中におけるリチウムイオンの伝導である。被膜抵抗Ｒ_ｆｉ２は、負極活物質の単位表面積当たりの被膜抵抗（単位［Ωｃｍ^２］）である。ここで、参照温度Ｔ_ｒｅｆにおける被膜抵抗Ｒ_ｆｉ２を予め求めておき、温度Ｔの依存性も考慮して、アレニウス式を適用すると、析出反応における被膜抵抗Ｒ_{ｆ３ｄｅｐｏ}は、上記式（１０）によって表すことができる。上記式（１０）に示すＥ_Ｒｆｉ２は、温度依存性を有する活性化エネルギである。温度Ｔ（電池温度Ｔｂ）を検出すれば、上記式（１０）に基づいて、析出反応における被膜抵抗Ｒ_{ｆ３ｄｅｐｏ}を算出できる。

一方、溶解反応における被膜抵抗Ｒ_{ｆ３ｄｉｓｓ}は、単位反応表面積当たりの活性リチウムの量（Ｑ_３／Ａ_３Ｓ）に比例して増加し、活性リチウムの量Ｑ_３がゼロであるときには、負極板の被膜における電子伝導の抵抗値（被膜抵抗という）Ｒ_ｆｅ２に相当する。電子伝導とは、被膜中における電子の伝導である。したがって、参照温度Ｔ_ｒｅｆにおける被膜抵抗Ｒ_ｆｅ２を予め求めておき、温度Ｔの依存性を考慮してアレニウス式を適用すると、溶解反応における被膜抵抗Ｒ_ｆ３’は、上記式（１１）によって表すことができる。上記式（１１）に示すＥ_Ｒｆｅ２は、温度依存性を有する活性化エネルギである。上記式（１１）に示すｋ_ｒｑは、被膜抵抗の増加係数（単位［ｃｍ^２／Ａｈ］）であり、この増加係数ｋ_ｒｑは、単位反応表面積当たりの活性リチウムの量に対する値である。増加係数ｋ_ｒｐは予め求めておくことができる。量Ｑ_３および面積比Ａ_３を更新し、温度Ｔ（電池温度Ｔｂ）を検出すれば、上記式（１１）に基づいて、溶解反応における被膜抵抗Ｒ_{ｆ３ｄｉｓｓ}を算出できる。

面積比Ａ_３は、下記式（１２）に基づいて算出することができる。

上記式（１２）において、Ａ_３０は面積比Ａ_３の初期値、ｋ_３は活性リチウムによる反応表面積の増加係数、ｋ_４は不活性リチウムによる反応表面積の増加係数である。初期値Ａ_３０および増加係数ｋ_３，ｋ_４は、予め定めておくことができる。また、量Ｑ_３，Ｑ_４の初期値としては、例えば、０に設定することができる。量Ｑ_３，Ｑ_４を算出するたびに、面積比Ａ_３を更新することができる。

次に、不活性化反応における反応電流密度ｉ_４’を算出する方法について説明する。

反応電流密度ｉ_４’は、バトラー・ボルマー式に基づいて、下記式（１３）によって表すことができる。図２に示すように、不活性化反応は不可逆的な還元反応であるため、下記式（１３）では、還元反応の電流だけを考慮している。

上記式（１３）において、ｉ_０４は、単位反応表面積当たりの交換電流密度（単位［Ａ／ｃｍ^２］）である。α_ｃ４は還元反応（すなわち、不活性化反応）の移動係数であり、予め求めておくことができる。η_４’は反応過電圧である。反応過電圧η_４’は、負極電位Ｖ_２’と、不活性化反応における平衡電位（定数）Ｕ_４との差になるため、下記式（１４）によって表すことができる。下記式（１４）では、被膜抵抗Ｒ_ｆ４’に伴う電圧降下量も考慮している。

上記式（１４）に示す反応過電圧η_４’を上記式（１３）に代入することにより、反応電流密度ｉ_４’を算出することができる。

上記式（１３）に示す交換電流密度ｉ_０４は、下記式（１５）に基づいて算出される。交換電流密度ｉ_０４は、温度Ｔ（すなわち、電池温度Ｔｂ）に依存する。このため、参照温度Ｔ_ｒｅｆにおける交換電流密度ｉ_０４を予め求めておき、温度Ｔの依存性も考慮して、アレニウス式を適用すると、交換電流密度ｉ_０４は、下記式（１５）によって表すことができる。ここで、交換電流密度ｉ_０４は、電解液中の塩濃度ｃ_ｅ２に依存しないものとしている。

上記式（１５）において、交換電流密度ｉ_０４（Ｔ_ｒｅｆ）は、参照温度Ｔ_ｒｅｆにおいて、交流インピーダンス測定から同定される。Ｅ_ｉ０４は、温度依存性を表す活性化エネルギである。温度Ｔ（電池温度Ｔｂ）を検出すれば、上記式（１５）に基づいて、交換電流密度ｉ_０４を算出できる。

ここで、下記式（１６）に示すように、上述した反応抵抗変化率ｋ_ｒａｔｅに基づいて、上記式（１５）から算出される交換電流密度ｉ_０４（Ｔ）を補正することもできる。この場合には、上記式（１３）に示す交換電流密度ｉ_０４として、下記式（１６）に基づいて算出される交換電流密度ｉ_０４’が用いられる。

上記式（１６）に示す反応抵抗変化率ｋ_ｒａｔｅは、不活性化反応における反応抵抗変化率ｋ_ｒａｔｅである。本実施例では、不活性化反応における反応抵抗変化率ｋ_ｒａｔｅを、析出反応および溶解反応における反応抵抗変化率ｋ_ｒａｔｅとみなしている。上述したように、反応抵抗変化率ｋ_ｒａｔｅは、拘束荷重および電池温度Ｔｂの少なくとも一方に基づいて算出できる。

上記式（１４）に示す被膜抵抗Ｒ_ｆ４’は、析出反応時に発生する不活性化反応における被膜抵抗Ｒ_{ｆ４ｄｅｐｏ}と、溶解反応時に発生する不活性化反応における被膜抵抗Ｒ_{ｆ４ｄｉｓｓ}とに分けられる。被膜抵抗Ｒ_{ｆ４ｄｅｐｏ}については、下記式（１７）に示すように、被膜抵抗Ｒ_{ｆ３ｄｅｐｏ}とみなすことができる。被膜抵抗Ｒ_{ｆ４ｄｉｓｓ}については、下記式（１８）に示すように、被膜抵抗Ｒ_{ｆ３ｄｉｓｓ}とみなすことができる。

二次電池１０の負極板では、下記式（１９）に示す電荷保存則が成り立つ。

上記式（１９）において、Ｉは電流センサ３２によって検出される電流値（電流値Ｉｂ）であり、Ｌ_１は正極活物質層の厚さであり、Ｌ_２は負極活物質層の厚さである。ｉ_１は、単位正極活物質表面積当たりの電流密度であり、ｉ_２は、単位活物質表面積当たりの電流密度である。電流密度ｉ_１，ｉ_２については後述する。ａ_ｓ１は、単位正極体積当たりの正極活物質の表面積を示す比表面積（単位［ｃｍ^２／ｃｍ^３］）である。すなわち、比表面積ａ_ｓ１は、正極活物質の表面積を、正極活物質層の体積で除算した値である。ａ_ｓ２は、単位負極体積当たりの負極活物質の表面積を示す比表面積（単位［ｃｍ^２／ｃｍ^３］）である。すなわち、比表面積ａ_ｓ２は、負極活物質の表面積を、負極活物質層の体積で除算した値である。厚さＬ_１，Ｌ_２および比表面積ａ_ｓ１，ａ_ｓ２は、予め求めておくことができる。

上記式（７），（１４），（１９）を満足するように、上記式（７），（１４），（１９）を連立して解くことにより、反応過電圧η_３’，η_４’を算出できる。反応過電圧η_３’ ，η_４’を算出すれば、上記式（６），（１３）に基づいて反応電流密度ｉ_３’，ｉ_４’を算出できる。上記式（４）に基づいて反応電流密度Ｉ_３’を算出すれば、上記式（２）に基づいて活性リチウムの量Ｑ_３を算出できる。上記式（５）に基づいて反応電流密度Ｉ_４’を算出すれば、上記式（３）に基づいて不活性リチウムの量Ｑ_４を算出できる。量Ｑ_３，Ｑ_４を算出すれば、上記式（１）に基づいて、リチウム析出量Ｑ_Ｌｉを算出できる。

ここで、反応過電圧η_３’，η_４’を算出するためには、上記式（７），（１４）に示す負極電位Ｖ_２’を推定する必要がある。負極電位Ｖ_２’は、下記式（２０）に基づいて算出（推定）できる。下記式（２０）では、負極活物質層の厚さ方向における電流分布が一定の電流値を示すものと仮定している。この電流分布とは、負極活物質層の厚さ方向における位置と、この位置に応じた電流値との関係を示す。

上記式（２０）において、Ｖ_２は負極電位（平均値）である。κ_２ ^ｅｆｆは負極板における実効イオン伝導率であり、塩濃度ｃ_ｅ２および温度Ｔに依存する。このため、実効イオン伝導率κ_２ ^ｅｆｆ、塩濃度ｃ_ｅ２および温度Ｔの対応関係（マップ又は演算式）を予め求めておけば、塩濃度ｃ_ｅ２を算出し、温度Ｔ（電池温度Ｔｂ）を検出することにより、実効イオン伝導率κ_２ ^ｅｆｆを算出できる。σ_２ ^ｅｆｆは負極活物質内の実効電子伝導率であり、予め決めておくことができる。

上記式（２０）に示す負極電位（平均値）Ｖ_２は、下記式（２１）に示す関係を有する。

上記式（２１）において、Ｖは電圧センサ３１によって検出された電圧値（電圧値Ｖｂ）、Ｖ_１は正極電位（平均値）である。Ｒ_ｄは二次電池１０の直流抵抗、Δφ_ｅは、正極板および負極板の間における濃度過電圧である。

上記式（２１）に示す直流抵抗Ｒ_ｄは、下記式（２２）に基づいて算出される。

上記式（２２）において、Ｒ_{ｄ，ｍａｐ}は、塩濃度（初期値）ｃ_{ｅ，ｉｎｉ}において温度Ｔに依存する直流抵抗成分である。直流抵抗成分Ｒ_{ｄ，ｍａｐ}は、温度Ｔ毎に交流インピーダンス測定から同定される。このため、温度Ｔ（電池温度Ｔｂ）を検出することにより、直流抵抗成分Ｒ_{ｄ，ｍａｐ}を算出できる。

上記式（２２）の右辺第２項および第３項では、塩濃度の変化に伴う直流抵抗成分の変化を算出している。具体的には、塩濃度ｃ_ｅ１，ｃ_ｅ２，ｃ_ｅｓおよび温度Ｔに依存する直流抵抗成分Ｒ_ｄ，ｅｓから、塩濃度（初期値）ｃ_{ｅ，ｉｎｉ}および温度Ｔに依存する直流抵抗成分Ｒ_ｄ，ｅｓを減算している。ここで、塩濃度ｃ_ｅ１は、正極板（正極活物質層）に含まれる電解液中の塩濃度であり、塩濃度ｃ_ｅｓは、セパレータに含まれる電解液中の塩濃度である。

上記式（２２）の右辺第２項に示す直流抵抗成分Ｒ_ｄ，ｅｓは、上記式（２３）に基づいて算出できる。正極板の実効イオン伝導率κ_１ ^ｅｆｆは、塩濃度ｃ_ｅ１および温度Ｔに依存するため、実効イオン伝導率κ_１ ^ｅｆｆ、塩濃度ｃ_ｅ１および温度Ｔの対応関係（マップ又は演算式）を予め求めておけば、塩濃度ｃ_ｅ１を算出し、温度Ｔ（電池温度Ｔｂ）を検出することにより、実効イオン伝導率κ_１ ^ｅｆｆを算出できる。

また、セパレータの実効イオン伝導率κ_ｓ ^ｅｆｆは、塩濃度ｃ_ｅｓおよび温度Ｔに依存するため、実効イオン伝導率κ_ｓ ^ｅｆｆ、塩濃度ｃ_ｅｓおよび温度Ｔの対応関係（マップ又は演算式）を予め求めておけば、塩濃度ｃ_ｅｓを算出し、温度Ｔ（電池温度Ｔｂ）を検出することにより、実効イオン伝導率κ_ｓ ^ｅｆｆを算出できる。正極の実効電子伝導率σ_１ ^ｅｆｆおよびセパレータの厚さＬ_ｓは、予め求めておくことができる。

上記式（２２）の右辺第３項に示す直流抵抗成分Ｒ_ｄ，ｅｓおよび温度Ｔの対応関係（マップ又は演算式）を予め求めておけば、温度Ｔ（電池温度Ｔｂ）を検出することにより、上記式（２２）の右辺第３項に示す直流抵抗成分Ｒ_ｄ，ｅｓを算出できる。

上記式（２１）に示す濃度過電圧Δφ_ｅは、下記式（２４）に基づいて算出できる。下記式（２４）において、ｔ_＋ ^０はカチオンの輸率、「１＋ｄｌｎｆ_±／ｄｌｎｃ_ｅ」は、塩濃度ｃ_ｅの変化量に対する濃度過電圧の平均の活量係数に関する相関係数である。塩濃度ｃ_ｅ１，ｃ_ｅ２を算出し、温度Ｔ（電池温度Ｔｂ）を検出すれば、下記式（２４）に基づいて、濃度過電圧Δφ_ｅを算出できる。

塩濃度ｃ_ｅｓを塩濃度（初期値）ｃ_{ｅ，ｉｎｉ}とすると、塩濃度ｃ_ｅ１，ｃ_ｅ２は、下記式（２５）に基づいて算出できる。

上記式（２５）において、ｔは時間、Δｔは、塩濃度の変化量Δｃ_ｅを算出する周期（所定時間）である。Ｄ_ｅは電解液中のリチウムの化学拡散係数である。化学拡散係数Ｄ_ｅは、温度Ｔに依存する。このため、化学拡散係数Ｄ_ｅおよび温度Ｔの対応関係（マップ又は演算式）を予め求めておけば、温度Ｔ（電池温度Ｔｂ）を検出することにより、化学拡散係数Ｄ_ｅを算出できる。係数α，βは、予め定めておくことができる。塩濃度の変化量（初期値）Δｃ_ｅ（０）としては、０が設定される。そして、所定時間Δｔが経過するたびに、今回の電流値Ｉ（電流値Ｉｂ）と、前回の塩濃度の変化量Δｃ_ｅ（ｔ）に基づいて、今回の塩濃度の変化量Δｃ_ｅ（ｔ＋Δｔ）を算出できる。また、塩濃度の変化量Δｃ_ｅに基づいて、塩濃度ｃ_ｅ１，ｃ_ｅ２を算出できる。

一方、電極電位Ｖ_ｊは、下記式（２６），（２７）に基づいて算出される。下記式（２６），（２７）において、添字ｊが１であるときには正極に関する値を示し、添字ｊが２であるときには負極に関する値を示す。例えば、Ｖ_１は正極電位を示し、Ｖ_２は負極電位を示す。なお、下記式（２８）〜（３０），（３２）〜（３４）に付された添字ｊについても同様である。

上記式（２６）において、Ｕ_ｊは電極（正極又は負極）の開放電位であり、後述する局所的ＳＯＣθ_ｊに依存する。開放電位Ｕ_ｊおよび局所的ＳＯＣθ_ｊの対応関係（マップ又は演算式）を予め求めておけば、局所的ＳＯＣθ_ｊを算出することにより、開放電位Ｕ_ｊを算出できる。η_ｊは電極の反応過電圧である。ｉ_ｃｔｊは、活物質に対してリチウムイオンが挿入・脱離する反応における反応電流密度であり、上記式（２７）に基づいて算出される。Ｒ_ｆｉｊは、電極板の被膜におけるイオン伝導の抵抗値（被膜抵抗という）である。

上記式（２７）において、α_ａｊは、電極板（正極板又は負極板）の酸化反応における移動係数であり、α_ｃｊは、電極の還元反応における移動係数である。ここで、負極の酸化反応は、二次電池１０の放電時の反応に相当し、負極の還元反応は、二次電池１０の充電時の反応に相当する。正極の酸化反応は、二次電池１０の充電時の反応に相当し、正極の還元反応は、二次電池１０の放電時の反応に相当する。ｉ_０ｊは、電極の交換電流密度であり、下記式（２８）に基づいて算出できる。

交換電流密度ｉ_０ｊは、後述する局所的ＳＯＣθ_ｊ、塩濃度ｃ_ｅｊおよび温度Ｔに依存する。したがって、上記式（２８）に基づいて、交換電流密度ｉ_０ｊを算出できる。Ｅ_ｉ０ｊは温度依存性を有する活性化エネルギである。局所的ＳＯＣθ_ｊおよび塩濃度ｃ_ｅｊを算出し、温度Ｔ（電池温度Ｔｂ）を検出すれば、上記式（２８）に基づいて、交換電流密度ｉ_０ｊを算出できる。

ここで、下記式（２９）に示すように、上述した反応抵抗変化率ｋ_ｒａｔｅに基づいて、上記式（２８）から算出される交換電流密度ｉ_０ｊを補正することもできる。この場合には、上記式（２７）に示す交換電流密度ｉ_０ｊとして、下記式（２９）に基づいて算出された交換電流密度ｉ_０ｊ’が用いられる。

上記式（２９）に示す反応抵抗変化率ｋ_ｒａｔｅは、リチウムイオンが挿入・脱離する反応における反応抵抗変化率ｋ_ｒａｔｅである。本実施例では、リチウムイオンが挿入・脱離する反応における反応抵抗変化率ｋ_ｒａｔｅを、析出反応および溶解反応における反応抵抗変化率ｋ_ｒａｔｅとみなしている。上述したように、反応抵抗変化率ｋ_ｒａｔｅは、拘束荷重および電池温度Ｔｂの少なくとも一方に基づいて算出できる。

上記式（２６）に示す被膜抵抗Ｒ_ｆｉｊは、上記式（１０）と同様に、下記式（３０）に基づいて算出できる。下記式（３０）に示すＥ_Ｒｆｉｊは、温度依存性を有する活性化エネルギである。

二次電池１０を充放電するときには、下記式（３１）に示す電荷保存則が成り立つ。なお、下記式（３１）では、反応電流密度ｉ_ｃｔｊだけを考慮しているが、これに限るものではない。二次電池１０の内部では、電気二重層キャパシタンス成分が存在するため、電気二重層キャパシタンス成分に関して、単位活物質表面積当たりの電流密度ｉ_ｄｌｊを考慮することもできる。ここで、電荷保存則に関して、電流密度ｉ_ｄｌｊは支配的ではないため、下記式（３１）では、電流密度ｉ_ｄｌｊを省略している。

ここで、上記式（１９）に示す電流密度ｉ_１としては、反応電流密度ｉ_ｃｔ１とすることもできるし、電流密度ｉ_ｃｔ１，ｉ_ｄｌ１の合計値とすることもできる。また、上記式（１９）に示す電流密度ｉ_２としては、反応電流密度ｉ_ｃｔ２とすることもできるし、電流密度ｉ_ｃｔ２，ｉ_ｄｌ２の合計値とすることもできる。

電圧値Ｖ（電圧値Ｖｂ）および電流値Ｉ（電流値Ｉｂ）を検出すれば、上記式（３１）を満足する条件の下で、上記式（２１），（２６）に基づいて、反応過電圧η_１，η_２を算出できる。これにより、負極電位Ｖ_２，Ｖ_２’を算出できる。

ここで、電極の局所的ＳＯＣθ_ｊ（θ_１又はθ_２）は、下記式（３２）に基づいて算出できる。局所的ＳＯＣθ_ｊとは、活物質の界面における塩濃度によって特定されるＳＯＣである。

上記式（３２）に示すｒは、活物質を球体と見なしたときの半径を示し、ｒ_ｓｊは、球体の活物質の表面（界面）を示す。局所的ＳＯＣθ_ｊを算出するときには、球体の活物質の表面における塩濃度ｃ_ｓｊが用いられる。ｃ_{ｓｊ，ｍａｘ}は、活物質の内部における限界塩濃度であり、予め求めておくことができる。活物質内の塩濃度分布は、下記式（３３）に示す拡散方程式に基づいて算出できる。活物質内の塩濃度分布とは、球体の活物質の半径方向における塩濃度ｃ_ｓｊの分布である。

上記式（３３）に示すＤ_ｓｊは、活物質内におけるリチウムの拡散係数であり、温度Ｔに依存する。拡散係数Ｄ_ｓｊおよび温度Ｔの対応関係（マップ又は演算式）を予め求めておけば、温度Ｔ（電池温度Ｔｂ）を検出することにより、拡散係数Ｄ_ｓｊを算出できる。上記式（３３）を適用するとき、下記式（３４）に示す境界条件が設定される。

上記式（３３），（３４）に基づいて、活物質の界面における塩濃度ｃ_ｓｊを算出すれば、上記式（３２）に基づいて、局所的ＳＯＣθ_ｊを算出できる。

次に、リチウム析出量Ｑ_Ｌｉを算出する処理について、図４に示すフローチャートを用いて説明する。図４に示す処理は、コントローラ４０によって実行される。図４に示す処理は、所定時間Δｔを周期として繰り返し行われる。なお、上述したように、所定時間Δｔを周期として、塩濃度ｃ_ｓｊ、局所的ＳＯＣθ_ｊや塩濃度の変化量Δｃ_ｅが更新されるとともに、これに依存するパラメータも更新される。

ステップＳ１０１において、コントローラ４０は、電圧センサ３１、電流センサ３２および温度センサ３３を用いて、電圧値Ｖｂ、電流値Ｉｂおよび電池温度Ｔｂを検出する。ステップＳ１０２において、コントローラ４０は、上記式（８）と、上記式（１５）又は（１６）とに基づいて、交換電流密度ｉ_０３，ｉ_０４を算出する。ステップＳ１０３において、コントローラ４０は、上記式（１０）又は（１１）と、上記式（１７）又は（１８）とに基づいて、被膜抵抗Ｒ_ｆ３’，Ｒ_ｆ４’を算出する。ステップＳ１０４において、コントローラ４０は、上記式（１２）に基づいて、面積比Ａ_３を算出する。ここで、面積比Ａ_４は、面積比Ａ_３と等しいとみなしているため、面積比Ａ_３を算出することにより、面積比Ａ_４も算出できる。

ステップＳ１０５において、コントローラ４０は、上記式（２０）に基づいて、負極電位Ｖ_２’を算出する。ステップＳ１０６において、コントローラ４０は、反応過電圧η_３’，η_４’を更新する。ステップＳ１０７において、コントローラ４０は、ステップＳ１０６の処理で更新された反応過電圧η_３’，η_４’に基づいて、反応電流密度Ｉ_３’，Ｉ_４’を算出する。具体的には、上記式（４），（６）に基づいて、反応電流密度Ｉ_３’が算出されるとともに、上記式（５），（１３）に基づいて、反応電流密度Ｉ_４’が算出される。

ステップＳ１０８において、コントローラ４０は、ステップＳ１０７の処理で算出された反応電流密度Ｉ_３’，Ｉ_４’に基づいて、上記式（２０）に示す電荷保存則が成立するか否かを判別する。電荷保存則が成立しないとき、コントローラ４０は、ステップＳ１０６の処理に戻り、反応過電圧η_３’，η_４’を更新する。このように、電荷保存則が成立するまで、反応過電圧η_３’，η_４’が更新される。

電荷保存則が成立するとき、コントローラ４０は、ステップＳ１０９において、上記式（２），（３）に基づいて、活性リチウムの量Ｑ_３および不活性リチウムの量Ｑ_４を算出する。また、ステップＳ１０８において、コントローラ４０は、上記式（１）に基づいて、リチウム析出量Ｑ_Ｌｉを算出する。

本実施例によれば、反応抵抗変化率ｋ_ｒａｔｅに基づいて交換電流密度ｉ_０３を補正し、補正後の交換電流密度ｉ_０３’に基づいて、反応電流密度ｉ_３’，Ｉ_３’を算出している。これにより、反応抵抗変化率ｋ_ｒａｔｅを考慮した量Ｑ_３やリチウム析出量Ｑ_Ｌｉを算出することができる。

１０：二次電池（リチウムイオン二次電池）、２０：負荷、３１：電圧センサ、
３２：電流センサ、３３：温度センサ、４０：コントローラ、４１：メモリ

Claims (1)

1. 充放電を行うリチウムイオン二次電池と、
   活性リチウムを析出させる析出反応および、活性リチウムを溶解させる溶解反応が行われたときの活性リチウムの量と、活性リチウムの不活性化反応によって生成された不活性リチウムの量との総量を示すリチウム析出量を算出するコントローラと、を有し、
   前記コントローラは、
   前記析出反応および前記溶解反応における過電圧および交換電流密度から算出される反応電流に基づいて、前記活性リチウムの量を算出し、
   前記不活性化反応における過電圧および交換電流密度から算出される反応電流に基づいて、前記不活性リチウムの量を算出し、
   前記リチウムイオン二次電池の負極における反応抵抗に関して、基準値および現在値の比率を示す反応抵抗変化率と、前記リチウムイオン二次電池に与えられる拘束荷重および前記リチウムイオン二次電池の温度のうちの少なくとも一方であるパラメータとの対応関係を用いて、取得した前記パラメータに応じた前記反応抵抗変化率を算出し、
   算出した前記反応抵抗変化率に基づいて、前記析出反応および前記溶解反応における交換電流密度を補正する、
   ことを特徴とする電池システム。