特許 6363426 電池システム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6363426

(24)【登録日】2018年7月6日

(45)【発行日】2018年7月25日

(54)【発明の名称】電池システム

(51)【国際特許分類】

   H01M 10/48 20060101AFI20180712BHJP

   H01M 10/42 20060101ALI20180712BHJP

   H01M 2/10 20060101ALI20180712BHJP

   H01M 10/0566 20100101ALN20180712BHJP

【ＦＩ】

   H01M10/48 102

   H01M10/48 301

   H01M10/42 P

   H01M2/10 E

   !H01M10/0566

【請求項の数】1

【全頁数】19

(21)【出願番号】特願2014-162626(P2014-162626)

(22)【出願日】2014年8月8日

(65)【公開番号】特開2016-39074(P2016-39074A)

(43)【公開日】2016年3月22日

【審査請求日】2017年3月2日

(73)【特許権者】

【識別番号】000003207

【氏名又は名称】トヨタ自動車株式会社

(73)【特許権者】

【識別番号】000003609

【氏名又は名称】株式会社豊田中央研究所

(74)【代理人】

【識別番号】100087398

【弁理士】

【氏名又は名称】水野 勝文

(74)【代理人】

【識別番号】100128783

【弁理士】

【氏名又は名称】井出 真

(74)【代理人】

【識別番号】100128473

【弁理士】

【氏名又は名称】須澤 洋

(72)【発明者】

【氏名】田代 広規

(72)【発明者】

【氏名】戸村 修二

【審査官】 坂東 博司

(56)【参考文献】

【文献】 特開２０１５−２０７４７３（ＪＰ，Ａ）

【文献】 特開２０１３−８３５２５（ＪＰ，Ａ）

【文献】 特開２０１３−２１３６９１（ＪＰ，Ａ）

【文献】 特開２０１３−０４４５８０（ＪＰ，Ａ）

【文献】 特開２００９−１００５１３（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｈ０１Ｍ １０／４８

Ｈ０１Ｍ ２／１０

Ｈ０１Ｍ １０／４２

Ｈ０１Ｍ １０／０５６６

(57)【特許請求の範囲】

【請求項1】

活物質を含む電極板を有し、充放電を行う発電要素が、電解液とともに電池ケースに収容された二次電池と、
前記電池ケースの外面と接触する突起部を有し、前記突起部を介して前記二次電池に所定方向の拘束荷重を与える拘束部材と、
前記発電要素に含浸された前記電解液中の塩濃度の分布を算出するコントローラと、を有し、
前記コントローラは、
前記電解液の流れを規定する液流れ方程式を用いて、前記発電要素の内部から前記発電要素の外部に向かって前記電解液が流れるときの流速と、前記発電要素の外部から前記発電要素の内部に向かって前記電解液が流れるときの流速とを、前記電解液が流れる方向における前記発電要素内の位置毎に算出し、
前記位置毎の前記流速に基づいて、前記発電要素内の前記塩濃度の分布を算出し、
前記流速を算出するとき、前記電解液の圧力と、前記二次電池の温度に応じて変化する前記電解液の密度と、前記発電要素内における前記電解液の体積分率とをパラメータとして含む前記液流れ方程式を用い、
前記発電要素内で前記電解液が存在する空間の体積変化量から算出される前記体積分率の変化量に基づいて、前記体積分率を算出し、
前記電解液を除く前記発電要素の体積膨張率から算出され、前記電解液を除く前記発電要素の体積変化量と、前記電解液を含む前記発電要素の体積変化量とに基づいて、前記電解液が存在する空間の体積変化量を算出し、
前記所定方向で前記突起部と対向する前記発電要素内の第１領域と、前記所定方向で前記突起部と対向しない前記発電要素内の第２領域とにおける体積変化量を前記電解液の圧力に基づいてそれぞれ算出し、これらの体積変化量を加算して、前記電解液を含む前記発電要素の体積変化量を算出する、
ことを特徴とする電池システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、電解液中の塩濃度の偏りを把握できる電池システムに関する。

【背景技術】

【0002】

特許文献１に記載されているように、電解液中の塩濃度が偏ることにより、二次電池の内部抵抗値が上昇することが知られている。特許文献１では、電池モデルから推定される推定電流密度と、電流測定部（電流センサ）によって測定された測定電流密度との差が大きくなるときに、塩濃度の偏りが大きくなり、二次電池の内部抵抗値が高くなることを認識している。この認識の下において、推定電流密度および測定電流密度の比を電池抵抗上昇率として算出している。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１０−０６０４０６号公報（段落［０１２７］〜［０１３２］）

【発明の概要】

【発明が解決しようとする課題】

【0004】

特許文献１では、塩濃度の偏りが大きくなるときの条件（推定電流密度および測定電流密度の差）を認識しているだけであり、塩濃度の偏りが発生する要因については認識していない。また、特許文献１では、推定電流密度および測定電流密度の差に基づいて、塩濃度の偏りを間接的に把握しているだけである。本願発明者によれば、電解液の体積変化（膨張および収縮）によって、塩濃度の偏りが発生することが分かった。この知見によれば、電解液中の塩濃度の偏りを直接的に把握することができる。

【課題を解決するための手段】

【0005】

本発明の電池システムは、二次電池と、拘束部材と、コントローラとを有する。二次電池では、発電要素および電解液が電池ケースに収容されており、発電要素は、活物質を含む電極板を有し、充放電を行う。拘束部材は、電池ケースの外面と接触する突起部を有しており、突起部を介して二次電池に所定方向の拘束荷重を与える。コントローラは、発電要素に含浸された電解液中の塩濃度の分布を算出する。

【0006】

コントローラは、電解液の流れを規定する液流れ方程式を用いて、発電要素の内部から発電要素の外部に向かって電解液が流れるときの流速と、発電要素の外部から発電要素の内部に向かって電解液が流れるときの流速とを、電解液が流れる方向における発電要素内の位置毎に算出する。そして、コントローラは、発電要素内の位置毎の流速に基づいて、発電要素内の塩濃度の分布を算出する。流速を算出するときには、電解液の圧力と、二次電池の温度に応じて変化する電解液の密度と、発電要素内における電解液の体積分率とをパラメータとして含む液流れ方程式が用いられる。

【0007】

塩濃度の偏りは、発電要素の内部から発電要素の外部に向かう電解液の流れと、発電要素の外部から発電要素の内部に向かう電解液の流れとによって発生する。このような電解液の流れは、電解液の体積変化（膨張および収縮）と、活物質の体積変化（膨張および収縮）とによって発生することが分かった。電解液の体積変化は、二次電池の温度に依存し、二次電池の温度に応じて変化する電解液の密度によって規定される。また、活物質の体積変化が発生すると、発電要素内の電解液の体積分率が変化する。このため、電解液の密度や体積分率をパラメータとして含む液流れ方程式を用いて電解液の流速を算出すれば、電解液や活物質の体積変化によって発生する塩濃度の偏り（すなわち、塩濃度の分布）を把握できる。

【0008】

ここで、電解液の体積分率は、体積分率の変化量に基づいて算出される。また、体積分率の変化量は、発電要素内で電解液が存在する空間の体積変化量から算出される。電解液が存在する空間の体積変化量は、電解液を除く発電要素の体積変化量と、電解液を含む発電要素の体積変化量とに基づいて算出される。電解液を除く発電要素の体積変化量は、電解液を除く発電要素の体積膨張率から算出される。電解液を含む発電要素の体積変化量は、所定方向で突起部と対向する発電要素内の第１領域と、所定方向で突起部と対向しない発電要素内の第２領域とにおける体積変化量を加算した量となる。第１領域および第２領域における体積変化量は、電解液の圧力に基づいて算出される。

【0009】

発電要素は拘束部材からの拘束荷重を受けるが、第１領域および第２領域が受ける拘束荷重は互いに異なる。これにより、第１領域の体積変化量と、第２領域の体積変化量とは互いに異なる。そこで、本発明では、第１領域および第２領域における体積変化量をそれぞれ算出した上で、電解液を含む発電要素の体積変化量を算出している。これにより、第１領域および第２領域における体積変化の相違を考慮して、電解液を含む発電要素の体積変化量を算出できる。

【図面の簡単な説明】

【0010】

【図1】電池システムの構成を示す図である。

【図2】二次電池の構成を示す図である。

【図3】発電要素の展開図である。

【図4】発電要素の外観図である。

【図5】二次電池および仕切り部材を示す上面図である。

【図6】負極板および正極板が対向する方向における塩濃度の偏りを説明する図である。

【図7】負極板および正極板のそれぞれの表面上における塩濃度の偏りを説明する図である。

【図8】電解液の膨張および収縮によって電解液が流れる方向を説明する図である。

【図9】発電要素の体積の変化量と、電解液の圧力との関係を示す図である。

【図10】抵抗上昇量および塩濃度差の関係を示す図である。

【図11】塩濃度分布を算出する処理を説明するフローチャートである。

【図12】二次電池の充放電を制御する処理を説明するフローチャートである。

【発明を実施するための形態】

【0011】

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

【実施例1】

【0012】

本発明の電池システムについて、図１を用いて説明する。組電池１は、正極ラインＰＬおよび負極ラインＮＬを介して、負荷２と接続されている。組電池１は、直列に接続された複数の二次電池（単電池）１０を有する。なお、組電池１には、並列に接続された複数の二次電池１０が含まれていてもよい。負荷２は、組電池１から出力された電力を受けて動作する。また、負荷２は、発電を行うこともでき、負荷２によって生成された電力は、組電池１に供給される。これにより、組電池１が充電される。

【0013】

図１に示す電池システムは、例えば、車両に搭載できる。負荷２としては、モータ・ジェネレータを用いることができる。モータ・ジェネレータは、組電池１から出力された電力を受けて、車両を走行させるための動力を生成する。モータ・ジェネレータが生成した動力は、車輪に伝達される。モータ・ジェネレータは、車両の制動時に発生する運動エネルギを電力に変換し、この電力を組電池１に供給できる。

【0014】

監視ユニット３１は、各二次電池１０の電圧値Ｖｂを検出し、検出結果をコントローラ４０に出力する。電流センサ３２は、二次電池１０の電流値Ｉｂを検出し、検出結果をコントローラ４０に出力する。本実施例において、二次電池１０を放電しているときの電流値Ｉｂを正の値とし、二次電池１０を充電しているときの電流値Ｉｂを負の値としている。

【0015】

コントローラ４０は、電圧値Ｖｂおよび電流値Ｉｂに基づいて、二次電池１０の充放電を制御できる。また、コントローラ４０は、電圧値Ｖｂや電流値Ｉｂに基づいて、二次電池１０のＳＯＣ（State of Charge）を算出できる。ＳＯＣとは、満充電容量に対する現在の充電容量の割合である。ＳＯＣの算出方法としては、公知の方法を適宜採用できるため、ＳＯＣの算出方法に関する詳細な説明は省略する。

【0016】

温度センサ３３は、組電池１（二次電池１０）の温度（電池温度）Ｔｂを検出し、検出結果をコントローラ４０に出力する。ここで、複数の温度センサ３３を配置することもできる。複数の温度センサ３３を用いることにより、互いに異なる位置に配置された複数の二次電池１０の温度Ｔｂをそれぞれ検出することができる。コントローラ４０は、メモリ４１を有する。メモリ４１は、コントローラ４０が所定の処理（特に、本実施例で説明する処理）を行うときに用いられる情報を記憶する。なお、メモリ４１は、コントローラ４０の外部に設けることもできる。

【0017】

次に、二次電池１０の構造について、図２を用いて説明する。図２において、Ｘ軸およびＺ軸は、互いに直交する軸である。本実施例では、鉛直方向に相当する軸をＺ軸としている。なお、Ｘ軸およびＺ軸と直交する軸をＹ軸とする。

【0018】

二次電池１０は、電池ケース１１０および発電要素１２０を有する。電池ケース１１０は、発電要素１２０を収容している。電池ケース１１０は密閉状態となっており、電池ケース１１０の内部には電解液が注入されている。電池ケース１１０には負極端子１１１および正極端子１１２が固定されている。負極端子１１１および正極端子１１２は、発電要素１２０と電気的に接続されている。

【0019】

発電要素１２０は、充放電を行う要素であり、図３に示すように、負極板（本発明の電極板に相当する）１２１と、正極板（本発明の電極板に相当する）１２２と、セパレータ１２３とを有する。図３は、発電要素１２０の一部を展開した図である。負極板１２１は、集電箔１２１ａと、集電箔１２１ａの表面に形成された負極活物質層１２１ｂとを有する。負極活物質層１２１ｂは、負極活物質、導電剤、バインダーなどを含んでいる。負極活物質層１２１ｂは、集電箔１２１ａの一部の領域に形成されており、集電箔１２１ａの残りの領域には、負極活物質層１２１ｂが形成されていない。

【0020】

正極板１２２は、集電箔１２２ａと、集電箔１２２ａの表面に形成された正極活物質層１２２ｂとを有する。正極活物質層１２２ｂは、正極活物質、導電剤、バインダーなどを含んでいる。正極活物質層１２２ｂは、集電箔１２２ａの一部の領域に形成されており、集電箔１２２ａの残りの領域には、正極活物質層１２２ｂが形成されていない。

【0021】

負極活物質層１２１ｂ、正極活物質層１２２ｂおよびセパレータ１２３には、電解液が含浸している。この電解液は、発電要素１２０の内部に存在する。一方、発電要素１２０の外部、言い換えれば、発電要素１２０および電池ケース１１０の間に形成されたスペースにも、余剰液としての電解液が存在している。

【0022】

図３に示す順番で、負極板１２１、正極板１２２およびセパレータ１２３を積層し、この積層体をＸ軸の周りで図４に示す矢印Ｄの方向に巻くことにより、発電要素１２０が構成される。ここで、負極板１２１および正極板１２２の間には、セパレータ１２３が配置される。

【0023】

Ｘ軸が延びる方向（Ｘ方向という）における発電要素１２０の一端では、負極板１２１の集電箔１２１ａだけが巻かれている。集電箔１２１ａだけが巻かれた部分は、図２に示す負極端子１１１と電気的に接続される。また、Ｘ方向における発電要素１２０の他端では、正極板１２２の集電箔１２２ａだけが巻かれている。集電箔１２２ａだけが巻かれた部分は、図２に示す正極端子１１２と電気的に接続される。

【0024】

本実施例では、上述したように、積層体を巻くことにより、発電要素１２０を構成しているが、これに限るものではない。具体的には、積層体を巻かずに、負極板１２１、正極板１２２およびセパレータ１２３を積層しただけで、発電要素１２０を構成することもできる。

【0025】

図４に示す領域Ａは、負極活物質層１２１ｂおよび正極活物質層１２２ｂがセパレータ１２３を挟んで互いに向かい合う領域である。領域Ａにおいて、二次電池１０（発電要素１２０）の充放電に応じた化学反応が行われる。

【0026】

組電池１を構成する複数の二次電池１０は、Ｙ軸が延びる方向（Ｙ方向という）に並べて配置されており、Ｙ方向で隣り合う２つの二次電池１０の間には、図５に示す仕切り部材（本発明の拘束部材に相当する）２０が配置されている。Ｙ方向における仕切り部材２０の一端面には、複数の突起部２１が設けられている。各突起部２１は、Ｙ方向に突出しており、突起部２１の先端は、電池ケース１１０の外面に接触している。Ｙ方向における仕切り部材２０の他端面は、平坦な面で構成されており、図５に示すように電池ケース１１０に接触している。

【0027】

突起部２１を電池ケース１１０に接触させることにより、仕切り部材２０および電池ケース１１０の間には、スペースが形成され、このスペースは、二次電池１０の温度を調節するための熱交換媒体（空気など）が移動する通路となる。仕切り部材２０および電池ケース１１０の間にスペースを形成できればよく、突起部２１の形状は適宜設計できる。

【0028】

Ｙ方向における組電池１の両端には、一対のエンドプレート（図示せず）が配置され、一対のエンドプレートには、Ｙ方向に延びる連結部材の両端が固定される。一対のエンドプレートを互いに近づく方向（Ｙ方向）に変位させることにより、二次電池１０は、仕切り部材２０からの拘束荷重を受ける。この拘束荷重は、Ｙ方向（本発明の所定方向に相当する）において、二次電池１０を拘束する荷重である。

【0029】

図５に示すように、発電要素１２０の領域Ａには、Ｙ方向において突起部２１と対向する第１領域Ａ１１と、Ｙ方向において突起部２１と対向しない第２領域Ａ１２とが含まれる。第１領域Ａ１１は、突起部２１からの荷重を受けやすく、第２領域Ａ１２は、突起部２１からの荷重を受けにくい。このため、第１領域Ａ１１が受ける荷重は、第２領域Ａ１２が受ける荷重よりも高くなる。

【0030】

二次電池１０では、電解液中の塩濃度に偏りが発生することにより、二次電池１０の内部抵抗値が上昇してしまう。このような内部抵抗値の上昇量を抵抗上昇量Ｒｈとする。抵抗上昇量Ｒｈは、二次電池１０の劣化に伴う内部抵抗値の上昇量とは異なる。劣化に伴う内部抵抗値の上昇量は増加するだけであり、減少することはない。一方、抵抗上昇量Ｒｈは、塩濃度の偏りに依存するため、塩濃度が偏るほど、抵抗上昇量Ｒｈが増加し、塩濃度の偏りが緩和されるほど、抵抗上昇量Ｒｈが減少する。

【0031】

塩濃度の偏りの状態としては、図６および図７に示す状態がある。図６は、負極板１２１および正極板１２２が対向する方向（Ｙ方向）において、塩濃度の偏りが発生する状態を示す。図６では、負極板１２１、正極板１２２およびセパレータ１２３の位置関係を示す概略図（図６の上側の図）と、塩濃度の分布（一例）を示す図（図６の下側の図）とを表している。

【0032】

塩濃度分布を示す図において、縦軸は塩濃度であり、横軸はＹ方向における位置である。図６（上側の図）では、負極板１２１および正極板１２２がセパレータ１２３から離れているが、実際には、負極板１２１および正極板１２２がセパレータ１２３に接触している。図６（下側の図）に示すように、二次電池１０の充電時では、塩濃度の偏りとして、実線で示す塩濃度分布が発生することがある。また、二次電池１０の放電時では、塩濃度の偏りとして、一点鎖線で示す塩濃度分布が発生することがある。

【0033】

なお、図６では、Ｙ方向における塩濃度の偏りを示しているが、これに限るものではない。上述したように、本実施例の発電要素１２０では、負極板１２１および正極板１２２がＸ軸の周りで巻かれているため、負極板１２１（負極活物質層１２１ｂ）および正極板１２２（正極活物質層１２２ｂ）が対向する方向において、図６と同様の塩濃度の偏りが発生する。

【0034】

二次電池１０の充放電を行うときには、負極板１２１および正極板１２２の間において、負極板１２１および正極板１２２が対向する方向（例えば、図６に示すＹ方向）に塩が移動する。二次電池１０がリチウムイオン二次電池であるとき、この塩はリチウム塩となる。負極板１２１および正極板１２２が対向する方向に塩が移動することによって、負極板１２１および正極板１２２が対向する方向において、塩濃度の偏りが発生する。

【0035】

図７は、負極板１２１および正極板１２２のそれぞれの表面上（領域Ａ内）において、塩濃度の偏りが発生する状態を示す。図６に示す塩濃度の偏りが発生することに応じて、図７に示す塩濃度の偏りが発生する。図７では、領域Ａを含む負極板１２１の一部と、領域Ａを含む正極板１２２の一部とを上下に分けて示している。ここで、負極板１２１の一部と、正極板１２２の一部とは、セパレータ１２３を挟んで対向している。

【0036】

図７の矢印で示すように、塩濃度の偏りは、領域Ａ内において、Ｘ方向に発生しやすい。図７では、負極板１２１および正極板１２２のそれぞれにおいて、領域Ａ内における塩濃度の分布（一例）も示している。塩濃度分布を示す図において、縦軸は塩濃度であり、横軸はＸ方向における位置である。ここで、二次電池１０の充電時では、塩濃度の偏りとして、実線で示す塩濃度分布が発生することがある。また、二次電池１０の放電時では、塩濃度の偏りとして、一点鎖線で示す塩濃度分布が発生することがある。

【0037】

上述したように、Ｘ方向における発電要素１２０の両端部では、負極板１２１（集電箔１２１ａ）や正極板１２２（集電箔１２２ａ）がＸ軸の周りで巻かれているだけである。このため、Ｘ方向における発電要素１２０の両端部では、電解液が通過しやすい。言い換えれば、発電要素１２０の内部から発電要素１２０の外部に向かって電解液が移動したり、発電要素１２０の外部から発電要素１２０の内部に向かって電解液が移動したりしやすい。

【0038】

これにより、図７に示すように、領域Ａ内のＸ方向において、塩濃度の偏りが発生しやすくなる。上述したように、負極板１２１、正極板１２２およびセパレータ１２３を積層しただけの構成であっても、発電要素１２０の内部から発電要素１２０の外部に向かって電解液が移動したり、発電要素１２０の外部から発電要素１２０の内部に向かって電解液が移動したりしやすい。

【0039】

図７に示す塩濃度の偏りは、電解液の流れによって発生することが分かった。また、この電解液の流れは、電解液の体積変化（膨張および収縮）によって発生することが分かった。具体的には、電解液が膨張すると、Ｘ方向において、発電要素１２０の内部から発電要素１２０の外部に向かう電解液の流れが発生する。すなわち、図８において、Ｘ方向における領域Ａの中心Ｃを基準として、矢印Ｘ１で示す方向に電解液が流れる。一方、電解液が収縮すると、Ｘ方向において、発電要素１２０の外部から発電要素１２０の内部に向かう電解液の流れが発生する。すなわち、図８において、領域Ａの中心Ｃに向かう方向（矢印Ｘ２で示す方向）に電解液が流れる。矢印Ｘ１，Ｘ２で示す方向に電解液が流れることにより、図７に示す塩濃度の偏りが発生する。

【0040】

そこで、本実施例では、電解液の膨張および収縮によって発生する電解液の流れ（流速）を算出することにより、図７に示す塩濃度の偏り（塩濃度分布）を把握するようにしている。図７に示す塩濃度の偏りを把握すれば、例えば、この塩濃度の偏りによって発生する抵抗上昇量Ｒｈを把握できる。本発明における塩濃度の分布とは、図７に示す塩濃度分布である。

【0041】

以下、電解液の流れ（流速）を算出する方法について説明する。

【0042】

下記式（１）は電解液の流れを規定する方程式であり、Brinkman-Navier-Stokes方程式として知られている。下記式（２）は電解液の流れに関する連続式であり、質量保存則から導かれる式である。下記式（１），（２）は、本発明の液流れ方程式に相当する。

【0043】

【数1】

【0044】

上記式（１），（２）において、ｕ_ｊは電解液の流速、ρは電解液の密度、ε_ｅ，ｊは電解液の体積分率、ｔは時刻である。また、上記式（１）において、μは電解液の粘度、Ｋ_ｊは透過係数、ｐは電解液の圧力である。

【0045】

ここで、添字ｊは、負極板１２１、正極板１２２およびセパレータ１２３を区別するために用いられ、添字ｊには「ｎ」、「ｐ」および「ｓ」が含まれる。添字ｊが「ｎ」であるときには、負極板１２１に関する値を示し、添字ｊが「ｐ」であるときには、正極板１２２に関する値を示し、添字ｊが「ｓ」であるときには、セパレータ１２３に関する値を示す。

【0046】

上述したように、電解液は、負極板１２１（負極活物質層１２１ｂ）、正極板１２２（正極活物質層１２２ｂ）およびセパレータ１２３のそれぞれに含浸している。このため、電解液に関するパラメータ（流速ｕ_ｊ、体積分率ε_ｅ，ｊ、透過係数Ｋ_ｊ）としては、負極板１２１、正極板１２２およびセパレータ１２３のそれぞれにおいて規定される。なお、本明細書では、上記式（１），（２）に示すパラメータ以外にも添字ｊを用いることがある。

【0047】

粘度μとしては予め定めた固定値を用いたり、電解液の温度に応じて粘度μを変更したりすることができる。電解液の温度としては、温度センサ３３によって検出される電池温度Ｔｂが用いられる。粘度μおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂを検出することにより粘度μを特定できる。粘度μおよび電池温度Ｔｂの対応関係を示す情報はメモリ４１に記憶される。

【0048】

密度ρは、電解液の膨張および収縮を規定するパラメータであり、電解液の膨張および収縮に応じた値を示す。電解液の膨張および収縮は電解液の温度（すなわち、電池温度Ｔｂ）に依存し、密度ρも電解液の温度（電池温度Ｔｂ）に依存する。したがって、密度ρおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂを検出することにより密度ρを特定できる。密度ρを特定することにより、電解液の膨張および収縮を把握できる。密度ρおよび電池温度Ｔｂの対応関係を示す情報はメモリ４１に記憶される。透過係数Ｋ_ｊは予め定めた固定値を用いることができる。体積分率ε_ｅ，ｊの算出方法については、後述する。

【0049】

上記式（１）では、流速ｕ_ｊおよび圧力ｐが未知数となるため、上記式（２）に示す連続式を規定して上記式（１），（２）の連立方程式を解くことにより、流速ｕ_ｊおよび圧力ｐを算出できる。上記式（１），（２）には、電解液の膨張および収縮を規定する密度ρが含まれているため、上記式（１），（２）の連立方程式を解くことにより、電解液の膨張および収縮に応じた流速ｕ_ｊを算出できる。この流速ｕ_ｊには、図８の矢印Ｘ１で示す方向の流速と、図８の矢印Ｘ２で示す方向の流速とが含まれる。流速ｕ_ｊを算出するときには、例えば、上記式（１），（２）を用いて収束計算を行うことができる。また、上記式（１），（２）を用いた演算は所定の周期で行われるが、前回の演算周期で算出された値を今回の演算周期で用いることにより、流速ｕ_ｊを算出できる。

【0050】

負極板１２１では、負極活物質層１２１ｂの内部において電解液が移動する。このため、流速ｕ_ｊ（すなわち、流速ｕ_ｎ）は、負極活物質層１２１ｂの内部における位置毎に算出される。図７に示す塩濃度の偏りを把握するときには、Ｘ方向で互いに異なる位置において、流速ｕ_ｎが算出される。

【0051】

正極板１２２では、正極活物質層１２２ｂの内部において電解液が移動する。このため、流速ｕ_ｊ（すなわち、流速ｕ_ｐ）は、正極活物質層１２２ｂの内部における位置毎に算出される。図７に示す塩濃度の偏りを把握するときには、Ｘ方向で互いに異なる位置において、流速ｕ_ｐが算出される。セパレータ１２３では、この内部において電解液が移動する。このため、流速ｕ_ｊ（すなわち、流速ｕ_ｓ）は、セパレータ１２３の内部における位置毎に算出される。図７に示す塩濃度の偏りを把握するときには、Ｘ方向で互いに異なる位置において、流速ｕ_ｓが算出される。

【0052】

一方、上記式（１），（２）において、様々な条件を仮定した上で、上記式（１）を簡易化することもできる。以下、上記式（１）を簡易化するときの手法（一例）について説明する。

【0053】

電解液の密度ρが流速ｕ_ｊを算出する位置（Ｘ方向の位置を含む）に関わらず一定であると仮定すると、上記式（２）は下記式（３）で表される。

【0054】

【数2】

【0055】

上記式（３）から下記式（４）が導き出せる。下記式（４）に示すｖ_ｐ，ｖ_ｎ，ｖ_ｓは、正極板１２２、負極板１２１およびセパレータ１２３のそれぞれにおける電解液の動粘性係数（ｖ＝μ／ρ）である。下記式（４）に示すｘ，ｙは、Ｘ方向およびＹ方向における位置をそれぞれ示す。ここで、Ｙ方向とは、負極板１２１および正極板１２２がセパレータ１２３を挟んで対向する方向（図８の上下方向）である。

【0056】

【数3】

【0057】

電解液の連続性を考慮すると、上記式（４）は下記式（５）で表される。

【0058】

【数4】

【0059】

一方、上記式（１）において、下記式（６）に示す仮定を行うと、下記式（７）が得られる。

【0060】

【数5】

【0061】

Ｘ方向における電解液の圧力分布に関して、負極板１２１、正極板１２２およびセパレータ１２３における圧力分布が互いに等しいと仮定すると、下記式（８）が得られる。

【0062】

【数6】

【0063】

上記式（５），（８）によれば、下記式（９）が導き出せる。下記式（９）は、負極板１２１において、Ｘ方向の位置に応じた流速ｕ_ｎを示している。下記式（９）に示すｘはＸ方向の位置を示している。

【0064】

【数7】

【0065】

上記式（９）をｘで積分し、ｘが０であるときの流速ｕ_ｎを０と仮定すると、下記式（１０）が得られる。ここで、ｘが０であるときの位置は、Ｘ方向における領域Ａの一端を示している。また、Ｘ方向における領域Ａの長さをＬとすると、Ｘ方向における領域Ａの他端については、ｘがＬになる。

【0066】

【数8】

【0067】

流速ｕ_ｎと同様に、流速ｕ_ｐ，ｕ_ｓは下記式（１１），（１２）で表される。

【0068】

【数9】

【0069】

上記式（１０）〜（１２）は上記式（１）を簡易化した式となり、負極板１２１、正極板１２２およびセパレータ１２３のそれぞれにおいて、Ｘ方向の位置に応じた流速ｕ_ｎ，ｕ_ｐ，ｕ_ｓを算出できる。上記式（１０）〜（１２）は、本発明における液流れ方程式に相当する。なお、電解液の流速ｕ_ｊを算出するための方程式は、上記（１），（２），（１０）〜（１２）に限るものではない。電解液の流れは電解液の体積変化によって発生するため、電解液の体積変化を規定する電解液の密度をパラメータ（変数）として含み、電解液の流れを規定できる方程式であれば、本発明を適用できる。

【0070】

流速ｕ_ｊを算出すれば、下記式（１３）に基づいて、電解液中の塩濃度ｃ_ｅ，ｊを算出できる。ここで、塩濃度ｃ_ｅ，ｊを算出するときには、流速ｕ_ｎ，ｕ_ｐ，ｕ_ｓのすべてを考慮してもよいし、流速ｕ_ｎ，ｕ_ｐ，ｕ_ｓの一部（例えば、流速ｕ_ｎ）だけを考慮してもよい。

【0071】

【数10】

【0072】

上記式（１３）において、Ｄ_ｅ，ｊ^ｅｆｆは、電解液の実効拡散係数であり、ｔ_＋^０は電解液中の塩の輸率である。Ｆはファラデー定数であり、ｊ_ｊは、単位体積および単位時間において、電解液中の塩の生成量である。

【0073】

上記式（１３）の左辺第１項は、所定時間Δｔにおける塩濃度の変化を規定している。上記式（１３）の左辺第２項は、電解液の流れ（流速ｕ_ｊ）に依存する塩濃度の変化を規定している。上記式（１３）の右辺第１項は、電解液中の塩の拡散状態を規定している。上記式（１３）の右辺第２項は、塩の生成量を規定している。ここで、二次電池１０の放電時では、負極板１２１の表面（負極活物質層１２１ｂ）において塩が生成され、二次電池１０の充電時では、正極板１２２の表面（正極活物質層１２２ｂ）において塩が生成される。

【0074】

次に、体積分率ε_ｅ，ｊの算出方法について説明する。電解液の流れは、電解液の膨張および収縮に加えて、活物質（正極活物質や負極活物質）の膨張および収縮によって発生する。活物質の膨張や収縮によって各活物質層１２１ｂ，１２２ｂの体積が変化する。各活物質層１２１ｂ，１２２ｂの体積が変化すると、各活物質層１２１ｂ，１２２ｂにおいて、電解液が存在する空間の体積が変化することに応じて体積分率ε_ｅ，ｊが変化する。

【0075】

正極活物質層１２２ｂ、負極活物質層１２１ｂおよびセパレータ１２３のそれぞれにおいて、電解液が存在する空間の体積Ｖ_ｅ，ｊの変化量ΔＶ_ｅ，ｊは、下記式（１４）に基づいて算出される。

【0076】

【数11】

【0077】

上記式（１４）において、Ｖ_ｓ，ｊは、各活物質層１２１ｂ，１２２ｂでは、各活物質層１２１ｂ，１２２ｂに含まれる活物質の体積（Ｖ_ｓ，ｎ、Ｖ_ｓ，ｐ）を示し、セパレータ１２３では、セパレータ１２３自体の体積（Ｖ_ｓ，ｓ）を示す。ここで、体積Ｖ_ｓ，ｊは、電解液を除く発電要素１２０の体積となる。そして、体積Ｖ_ｓ，ｊの変化量がΔＶ_ｓ，ｊとなる。Ｖ_{ａｌｌ，ｊ}は体積Ｖ_ｓ，ｊおよび体積Ｖ_ｅ，ｊの総和である。ここで、体積Ｖ_{ａｌｌ，ｊ}は、電解液を含む発電要素１２０の体積となる。そして、体積Ｖ_{ａｌｌ，ｊ}の変化量がΔＶ_{ａｌｌ，ｊ}となる。

【0078】

変化量ΔＶ_{ａｌｌ，ｊ}および変化量ΔＶ_ｓ，ｊは、下記式（１５），（１６）によって表される。

【0079】

【数12】

【0080】

体積Ｖ_{ａｌｌ，ｊ}は、第１領域Ａ１１および第２領域Ａ１２における体積Ｖ_ｓ，ｊ，Ｖ_ｅ，ｊの総和となる。このため、上記式（１５）に示すように、変化量ΔＶ_{ａｌｌ，ｊ}は、図５に示す第１領域Ａ１１における体積（体積Ｖ_ｓ，ｊ，Ｖ_ｅ，ｊの合計）の変化量ΔＶ_１１，ｊと、図５に示す第２領域Ａ１２における体積（体積Ｖ_ｓ，ｊ，Ｖ_ｅ，ｊの合計）の変化量ΔＶ_１２，ｊとの和になる。ここで、図５に示すように、第１領域Ａ１１が複数設けられているため、変化量ΔＶ_１１，ｊは、すべての第１領域Ａ１１における体積の変化量の総和となる。同様に、変化量ΔＶ_１２，ｊは、すべての第２領域Ａ１２における体積の変化量の総和となる。

【0081】

各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊは、電解液の圧力ｐに依存する。このため、各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊおよび圧力ｐの対応関係（マップ又は演算式）を実験などに基づいて予め求めておくことができる。図９には、各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊおよび圧力ｐの対応関係（一例）を示す。図９に示すように、圧力ｐの変化に応じて、各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊが変化する。ただし、圧力ｐが高いほど、各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊが変化しにくくなる。

【0082】

ここで、圧力ｐが等しいとき、変化量ΔＶ_１２，ｊは、変化量ΔＶ_１１，ｊよりも大きくなる。第２領域Ａ１２が受ける荷重は、第１領域Ａ１１が受ける荷重よりも低いため、変化量ΔＶ_１２，ｊは、変化量ΔＶ_１１，ｊよりも大きくなる。各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊおよび圧力ｐの対応関係に基づいて、各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊを算出するときには、圧力ｐとして、前回の演算周期で算出された圧力ｐの値を用いることができる。

【0083】

上述したように、第１領域Ａ１１および第２領域Ａ１２が受ける荷重は互いに異なる。このため、第１領域Ａ１１における体積の変化量と、第２領域Ａ１２における体積の変化量とは互いに異なる。そこで、上述したように、第１領域Ａ１１および第２領域Ａ１２において、各変化量ΔＶ_１１，ｊ，ΔＶ_１２，ｊを算出することにより、第１領域Ａ１１および第２領域Ａ１２における体積変化の相違を考慮して、変化量ΔＶ_{ａｌｌ，ｊ}を算出できる。

【0084】

上記式（１６）において、β_ｊ（具体的にはβ_ｎ又はβ_ｐ）は活物質（負極活物質又は正極活物質）の体積膨張率である。上記式（１６）に示す体積分率ε_ｅ，ｊおよび体積Ｖ_{ａｌｌ，ｊ}としては、前回の演算周期で算出された値が用いられる。

【0085】

ここで、変化量Δε_ｅ，ｊを最初に算出するとき、前回の体積分率ε_ｅ，ｊとしては、電池温度Ｔｂや二次電池１０のＳＯＣに応じた値を用いることができる。電池温度ＴｂおよびＳＯＣの少なくとも一方と、体積分率ε_ｅ，ｊとの対応関係を示す情報を実験などによって予め求めておけば、電池温度ＴｂやＳＯＣを特定することにより、体積分率ε_ｅ，ｊを特定できる。

【0086】

体積膨張率β_ｊは線膨張率α_ｊから算出できる。すなわち、線膨張率α_ｊの３倍の値が体積膨張率β_ｊとなる。ここで、線膨張率α_ｊは、電池温度Ｔｂや二次電池１０のＳＯＣに依存するため、電池温度ＴｂおよびＳＯＣの少なくとも一方から線膨張率α_ｊを特定できる。例えば、線膨張率α_ｊおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂの検出又は推定によって線膨張率α_ｊを特定できる。また、線膨張率α_ｊおよびＳＯＣの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、ＳＯＣの推定によって線膨張率α_ｊを特定できる。上述した対応関係を示す情報はメモリ４１に記憶しておくことができる。なお、線膨張率α_ｊを算出せずに、電池温度ＴｂおよびＳＯＣの少なくとも一方から体積膨張率β_ｊを直接算出することもできる。

【0087】

なお、上述した説明では、透過係数Ｋ_ｊを固定値としているが、体積膨張率β_ｊに基づいて透過係数Ｋ_ｊを変更することもできる。具体的には、透過係数Ｋ_ｊおよび体積膨張率β_ｊの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、体積膨張率β_ｊに応じた透過係数Ｋ_ｊを特定できる。このように特定した透過係数Ｋ_ｊは、上記式（１）又は上記式（１０）〜（１２）において用いることができる。

【0088】

上記式（１５），（１６）に基づいて変化量ΔＶ_{ａｌｌ，ｊ}，ΔＶ_ｓ，ｊを算出すれば、上記式（１４）に基づいて変化量ΔＶ_ｅ，ｊを算出できる。ここで、上記式（１４）には体積膨張率β_ｊが含まれているため、体積膨張率β_ｊに基づいて変化量ΔＶ_ｅ，ｊを算出できる。

【0089】

電解液の流速ｕ_ｊを算出する式（上記式（１），（１０）〜（１２））では、変化量ΔＶ_ｅ，ｊがパラメータとして含まれていない。ただし、電解液が存在する空間の体積が変化するときには、電解液の体積分率ε_ｅ，ｊが変化するため、変化量ΔＶ_ｅ，ｊを変化量Δε_ｅ，ｊに変換することができる。具体的には、体積Ｖ_{ａｌｌ，ｊ}が予め定めた一定の体積Ｖ_{ａｌｌ＿ｍｏｄｅｌ，ｊ}であると仮定したとき、下記式（１７）に示すように、体積Ｖ_{ａｌｌ＿ｍｏｄｅｌ，ｊ}および変化量ΔＶ_ｅ，ｊから変化量Δε_ｅ，ｊを算出できる。

【0090】

【数13】

【0091】

変化量Δε_ｅ，ｊを算出すれば、前回の演算周期で算出された体積分率ε_ｅ，ｊに変化量Δε_ｅ，ｊを加算することにより、今回の体積分率ε_ｅ，ｊを算出できる。今回算出した体積分率ε_ｅ，ｊを用いれば、上記式（１），（２）に基づいて、又は、上記式（１０）〜（１２）に基づいて、流速ｕ_ｊを算出できる。また、上記式（１３）に基づいて塩濃度ｃ_ｅ，ｊを算出するときにも、今回算出した体積分率ε_ｅ，ｊを用いることができる。

【0092】

なお、各活物質層１２１ｂ，１２２ｂに関する変化量Δε_ｅ，ｎ，Δε_ｅ，ｐに加えて、セパレータ１２３に関する変化量Δε_ｅ，ｓを算出してもよい。変化量Δε_ｅ，ｓを算出するとき、上記式（１６）に示す体積膨張率β_ｊ（すなわち、β_ｓ）は、セパレータ１２３を構成する材料の体積膨張率となる。体積膨張率β_ｓは電池温度Ｔｂに依存する。このため、体積膨張率β_ｓおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂの検出又は推定によって体積膨張率β_ｓを特定できる。

【0093】

上記式（１３）によれば、図６および図７に示す塩濃度の偏りを把握することができる。上記式（１３）を解くことにより塩濃度ｃ_ｅ，ｊを算出できる。ここで、流速ｕ_ｊとしては、Ｘ方向の位置に応じた流速ｕ_ｊが用いられるため、上記式（１３）を解くことにより、Ｘ方向の位置に応じた塩濃度ｃ_ｅ，ｊを算出できる。これにより、Ｘ方向における塩濃度ｃ_ｅ，ｊの分布（図７参照）を算出できる。

【0094】

なお、塩濃度ｃ_ｅ，ｊの算出方法は、上記式（１３）を用いた方法に限るものではない。Ｘ方向において塩濃度ｃ_ｅ，ｊのバラツキが発生している状態であれば、流速ｕ_ｊに基づいて、Ｘ方向における塩濃度ｃ_ｅ，ｊの分布を算出できる。上記式（１３）では、塩濃度ｃ_ｅ，ｊのバラツキを把握するために、塩の拡散状態および塩の生成量を規定している。例えば、塩濃度ｃ_ｅ，ｊのバラツキを予め設定すれば、塩の拡散状態や生成量を考慮せずに、流速ｕ_ｊに基づいて塩濃度ｃ_ｅ，ｊの分布を算出できる。

【0095】

また、上記式（１３）によれば、電解液中の塩の拡散状態を規定しているが、塩の拡散状態を考慮しなくてもよい。塩の拡散に関する時定数は、塩の生成に関する時定数よりも大きいため、塩の拡散が発生していないこともある。そこで、塩濃度ｃ_ｅ，ｊの分布を算出するときに、塩の拡散状態を考慮しないようにしてもよい。

【0096】

上述したように塩濃度ｃ_ｅ，ｊの分布を算出すれば、塩濃度ｃ_ｅ，ｊの最大差（塩濃度差）Δｃ_ｅ，ｊ＿ｍａｘを算出できる。塩濃度差Δｃ_ｅ，ｊ＿ｍａｘは、塩濃度（最大値）ｃ_ｅおよび塩濃度（最小値）ｃ_ｅの差である。

【0097】

図１０に示すように、抵抗上昇量Ｒｈおよび塩濃度差Δｃ_ｅ＿ｍａｘの対応関係を実験などによって予め求めておけば、塩濃度差Δｃ_ｅ＿ｍａｘを算出することにより、この塩濃度差Δｃ_ｅ＿ｍａｘに対応した抵抗上昇量Ｒｈを算出できる。図１０に示すように、塩濃度差Δｃ_ｅ＿ｍａｘが大きくなるほど、抵抗上昇量Ｒｈが大きくなる。言い換えれば、塩濃度差Δｃ_ｅ＿ｍａｘが小さくなるほど、抵抗上昇量Ｒｈが小さくなる。上記式（１３）によれば、塩濃度ｃ_ｅは、負極板１２１および正極板１２２のそれぞれで算出される。ここで、図１０に示す塩濃度差Δｃ_ｅ＿ｍａｘを算出するときには、まず、負極板１２１および正極板１２２における塩濃度ｃ_ｅの分布を合算する。具体的には、負極板１２１および正極板１２２において、互いに対向する位置における塩濃度ｃ_ｅを合算する。そして、合算した塩濃度ｃ_ｅの分布において、塩濃度（最大値）ｃ_ｅおよび塩濃度（最小値）ｃ_ｅの差を塩濃度差Δｃ_ｅ＿ｍａｘとして算出する。

【0098】

抵抗上昇量Ｒｈおよび塩濃度差Δｃ_ｅ＿ｍａｘの対応関係は、マップ又は演算式として表すことができる。そして、この対応関係を特定する情報はメモリ４１に記憶しておくことができる。上記式（１３）によれば、所定時間（演算周期）Δｔが経過するたびに塩濃度ｃ_ｅ，ｊが算出されて塩濃度ｃ_ｅ，ｊの分布を把握できるため、所定時間Δｔが経過するたびに抵抗上昇量Ｒｈが算出される。これに伴い、抵抗上昇量Ｒｈの変化を把握できる。

【0099】

一方、抵抗上昇量Ｒｈは、平均塩濃度ｃ_ｅ＿ａｖｅに依存することがある。そこで、平均塩濃度ｃ_ｅ＿ａｖｅに基づいて、抵抗上昇量Ｒｈを算出することもできる。平均塩濃度ｃ_ｅ＿ａｖｅとは、負極板１２１および正極板１２２における塩濃度（上述した合算値）ｃ_ｅの分布を平均化した塩濃度（平均値）である。塩濃度差Δｃ_ｅ＿ｍａｘと同様に、抵抗上昇量Ｒｈおよび平均塩濃度ｃ_ｅ＿ａｖｅの対応関係（マップ又は演算式）を実験などによって予め求めておけば、平均塩濃度ｃ_ｅ＿ａｖｅを算出することにより、この平均塩濃度ｃ_ｅ＿ａｖｅに対応した抵抗上昇量Ｒｈを算出できる。

【0100】

また、塩濃度差Δｃ_ｅ＿ｍａｘおよび平均塩濃度ｃ_ｅ＿ａｖｅに基づいて、抵抗上昇量Ｒｈを算出することもできる。この場合には、塩濃度差Δｃ_ｅ＿ｍａｘ、平均塩濃度ｃ_ｅ＿ａｖｅおよび抵抗上昇量Ｒｈの対応関係（マップ又は演算式）を予め求めておけばよい。塩濃度差Δｃ_ｅ＿ｍａｘおよび平均塩濃度ｃ_ｅ＿ａｖｅを算出すれば、算出した塩濃度差Δｃ_ｅ＿ｍａｘおよび平均塩濃度ｃ_ｅ＿ａｖｅに対応する抵抗上昇量Ｒｈを算出できる。

【0101】

図１１は、塩濃度分布を算出する処理を示すフローチャートである。図１１に示す処理は、コントローラ４０によって実行される。

【0102】

ステップＳ１０１において、コントローラ４０は、温度センサ３３を用いて電池温度Ｔｂを検出する。ステップＳ１０２において、コントローラ４０は、電解液の流速ｕ_ｊを算出する。具体的には、上述したように、コントローラ４０は、Ｘ方向の位置に応じた流速ｕ_ｊを算出する。流速ｕ_ｊを算出するときには、ステップＳ１０１の処理で検出された電池温度Ｔｂが用いられる。

【0103】

ステップＳ１０３において、コントローラ４０は、ステップＳ１０２の処理で算出された流速ｕ_ｊに基づいて塩濃度ｃ_ｅ，ｊの分布を算出する。塩濃度ｃ_ｅ，ｊの分布を算出する方法は上述した通りである。

【0104】

本実施例では、ステップＳ１０１の処理において、温度センサ３３を用いて電池温度Ｔｂを検出しているが、これに限るものではない。具体的には、ステップＳ１０１の処理において、二次電池１０を充放電したときの発熱量（電池温度Ｔｂの上昇量）と、二次電池１０の放熱量（電池温度Ｔｂの低下量）とを考慮することにより、電池温度Ｔｂを推定できる。

【0105】

二次電池１０の発熱量は、二次電池１０の電流値Ｉｂおよび内部抵抗値から算出できる。二次電池１０の内部抵抗値は、電流値Ｉｂおよび電圧値Ｖｂから算出できる。二次電池１０の放熱量は、電池温度Ｔｂと、二次電池１０の周囲に存在する大気中の温度（環境温度）とに基づいて算出できる。ここで、電池温度Ｔｂが環境温度よりも高いときにおいて、電池温度Ｔｂおよび環境温度の差が広がるほど、放熱量が大きくなりやすい。

【0106】

二次電池１０を充放電していなければ、電池温度Ｔｂは環境温度と等しくなる。そこで、温度センサ３３とは異なる温度センサを用いて環境温度を検出しておき、二次電池１０を充放電するたびに算出される発熱量と、二次電池１０の放熱量とに基づいて、現在の電池温度Ｔｂを算出（推定）できる。発熱量および放熱量に基づいて電池温度Ｔｂを推定する方法は、公知の方法を適宜採用できる。二次電池１０を車両に搭載したとき、車両には、車両の外部における温度を検出する温度センサが設けられている。この温度センサを用いて、環境温度を検出できる。

【0107】

図１１に示す処理によって塩濃度分布を算出したときには、この塩濃度分布に基づいて抵抗上昇量Ｒｈを算出し、この抵抗上昇量Ｒｈに基づいて二次電池１０の充放電を制御できる。ここで、二次電池１０の充放電を制御するときの処理（一例）について、図１２に示すフローチャートを用いて説明する。図１２に示す処理は、コントローラ４０によって実行される。

【0108】

ステップＳ２０１において、コントローラ４０は、図１１に示す処理で算出した塩濃度ｃ_ｅ，ｊの分布に基づいて抵抗上昇量Ｒｈを算出する。抵抗上昇量Ｒｈを算出する方法は上述した通りである。ステップＳ２０２において、コントローラ４０は、抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるか否かを判別する。閾値Ｒｈ＿ｔｈは、抵抗上昇量Ｒｈの上限値であり、二次電池１０の劣化を抑制する観点に基づいて適宜設定できる。閾値Ｒｈ＿ｔｈを特定する情報は、メモリ４１に記憶しておくことができる。

【0109】

抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるとき、コントローラ４０は、ステップＳ２０３において、充電電力許容値Ｗｉｎや放電電力許容値Ｗｏｕｔを低下させる。充電電力許容値Ｗｉｎは、二次電池１０の充電を許容する上限の電力値であり、放電電力許容値Ｗｏｕｔは、二次電池１０の放電を許容する上限の電力値である。

【0110】

上述したように、二次電池１０を充電したときの電流値Ｉｂは負の値となるため、充電電力値も負の値となる。一方、二次電池１０を放電したときの電流値Ｉｂは正の値となるため、放電電力値は正の値となる。二次電池１０を充電するときには、充電電力値（絶対値）が充電電力許容値（絶対値）Ｗｉｎよりも高くならないように充電が制御される。また、二次電池１０を放電するときには、放電電力値が放電電力許容値Ｗｏｕｔよりも高くならないように放電が制御される。

【0111】

電池温度Ｔｂや二次電池１０のＳＯＣに基づいて、基準値としての充電電力許容値Ｗｉｎ＿ｒｅｆや放電電力許容値Ｗｏｕｔ＿ｒｅｆが設定される。ステップＳ２０３の処理では、充電電力許容値（絶対値）Ｗｉｎを充電電力許容値（絶対値）Ｗｉｎ＿ｒｅｆよりも低下させたり、放電電力許容値Ｗｏｕｔを放電電力許容値Ｗｏｕｔ＿ｒｅｆよりも低下させたりする。充電電力許容値Ｗｉｎや放電電力許容値Ｗｏｕｔを低下させることにより、抵抗上昇量Ｒｈの増加を抑制できる。

【0112】

抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈよりも小さいとき、コントローラ４０は、図１２に示す処理を終了する。このとき、充電電力許容値Ｗｉｎとしては、上述した充電電力許容値Ｗｉｎ＿ｒｅｆが設定され、放電電力許容値Ｗｏｕｔとしては、上述した放電電力許容値Ｗｏｕｔ＿ｒｅｆが設定される。

【0113】

図１２に示す処理では、抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるとき、充電電力許容値（絶対値）Ｗｉｎおよび放電電力許容値Ｗｏｕｔを低下させているが、これに限るものではない。具体的には、充電電力許容値（絶対値）Ｗｉｎおよび放電電力許容値Ｗｏｕｔの一方だけを低下させることができる。

【0114】

抵抗上昇量Ｒｈとしては、二次電池１０の充電に起因した抵抗上昇量Ｒｈと、二次電池１０の放電に起因した抵抗上昇量Ｒｈとがある。図７に示すように、充電時および放電時において、塩濃度ｃ_ｅ，ｊの分布が異なる。したがって、塩濃度ｃ_ｅ，ｊの分布を把握すれば、二次電池１０の充電に起因した抵抗上昇量Ｒｈと、二次電池１０の放電に起因した抵抗上昇量Ｒｈとを区別できる。

【0115】

ここで、二次電池１０の充電によって塩濃度差Δｃ_ｅ＿ｍａｘが発生し、この塩濃度差Δｃ_ｅ＿ｍａｘから算出された抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるときには、充電電力許容値（絶対値）Ｗｉｎだけを低下させることができる。これにより、充電に起因した抵抗上昇量Ｒｈに関して、抵抗上昇量Ｒｈの増加を抑制したり、抵抗上昇量Ｒｈを減少させたりすることができる。

【0116】

一方、二次電池１０の放電によって塩濃度差Δｃ_ｅ＿ｍａｘが発生し、この塩濃度差Δｃ_ｅ＿ｍａｘから算出された抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるときには、放電電力許容値Ｗｏｕｔだけを低下させることができる。これにより、放電に起因した抵抗上昇量Ｒｈに関して、抵抗上昇量Ｒｈの増加を抑制したり、抵抗上昇量Ｒｈを減少させたりすることができる。

【0117】

本実施例では、塩濃度ｃ_ｅの分布に基づいて抵抗上昇量Ｒｈを算出しているが、これに限るものではない。すなわち、本実施例のように塩濃度ｃ_ｅの分布を算出すれば、この塩濃度ｃ_ｅの分布を、二次電池１０の内部状態を把握するための様々な技術に応用することができる。二次電池１０の内部状態としては、例えば、以下に説明するように、ＳＯＣや電流値の分布がある。

【0118】

例えば、特開２０１３−０８３５２５号公報に記載の技術によれば、二次電池１０のＯＣＶ（Open Circuit Voltage）からＳＯＣ（推定値）を算出し、ＳＯＣ（補正量）を用いてＳＯＣ（推定値）を補正している。ＳＯＣ（補正量）は、ＳＯＣ（推定値）および電池温度Ｔｂから特定される基本値にゲインを乗算することによって算出される。ゲインは、塩濃度の偏りに応じて設定される。

【0119】

特開２０１３−０８３５２５号公報では、二次電池１０の電流値を積算した電流積算値と電池温度Ｔｂとに基づいて、ゲインを算出している。すなわち、電流積算値および電池温度Ｔｂを用いて塩濃度の分布を間接的に把握した上でゲインを算出している。ここで、本実施例のように、塩濃度ｃ_ｅの分布を算出すれば、この塩濃度ｃ_ｅの分布に基づいてゲインを算出することができる。

【0120】

一方、塩濃度ｃ_ｅの分布を把握できれば、図７に示すＸ方向の位置に応じた電流値の分布を把握できる。塩濃度ｃ_ｅが変化すれば、電流の流れやすさ（すなわち、電流値）が変化する。このため、塩濃度ｃ_ｅが互いに異なる位置（図７に示すＸ方向の位置）では、電流値が互いに異なることになる。二次電池１０がリチウムイオン二次電池であるとき、電流が流れやすくなるほど、リチウムが析出されやすくなることが知られている。このため、塩濃度ｃ_ｅの分布から把握される電流値の分布を用いて、リチウムの析出状態を把握することも可能と考えられる。

【符号の説明】

【0121】

１：組電池、１０：二次電池、２：負荷、３１：監視ユニット、３２：電流センサ、３３：温度センサ、４０：コントローラ、４１：メモリ、１１０：電池ケース、
１１１：負極端子、１１２：正極端子、１２０：発電要素、１２１：負極板、
１２１ａ：集電箔、１２１ｂ：負極活物質層、１２２：正極板、１２２ａ：集電箔、
１２２ｂ：正極活物質層、１２３：セパレータ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】