(57)【要約】【課題】２次電池の容量維持率を高精度で推定することができる推定方法を提供する。【解決手段】２次電池１の容量維持率の時間的変化率である劣化速度ｙ’を時間間隔をおいて次式（Ａ）により逐次推定する第１ステップ（ＳＴＥＰ１〜４）と、劣化速度ｙ’の推定値を用いて２次電池１の容量維持率ｙを推定する第２ステップ（ＳＴＥＰ５）とを備える。  ｙ’＝(ｃ／(１−ｙ)(1/c)-1)・(ａ・ｂ＋ｚ)    ……（Ａ）【選択図】図５
2014052186.tif 000002
【特許請求の範囲】
【請求項１】
２次電池の容量維持率を推定する方法であって、
前記２次電池の容量維持率の時間的変化率である劣化速度ｙ’を時間間隔をおいて次式（Ａ）により逐次推定する第１ステップと、該劣化速度ｙ’の推定値を用いて前記２次電池の容量維持率を推定する第２ステップとを備えることを特徴とする２次電池の容量維持率推定方法。
ｙ’＝(ｃ／(１−ｙ)(1/c)-1)・(ａ・ｂ＋ｚ)    ……（Ａ）
ただし、
ｃ：あらかじめ定めた所定値
ｙ：劣化速度ｙ’の今回の推定時点での容量維持率の推定値
ａ，ｚ：劣化速度ｙ’の今回の推定時点又はその直前における２次電池の使用条件に応じて値が決定される変数
ｂ：２次電池の初期状態からの経過時間ｔの増加に伴い、値が増加していくように決定される変数
【請求項２】
請求項１記載の２次電池の容量維持率推定方法において、
前記劣化速度ｙ’の各推定時点において、前記２次電池に単一又は複数の周波数の交流電流を通電し、該交流電流の各周波数における前記２次電池のインピーダンスを計測するステップと、該インピーダンスの計測データに基づいて前記変数ｂの値を決定するステップとをさらに備えることを特徴とする２次電池の容量維持率推定方法。
【請求項３】
請求項２記載の２次電池の容量維持率推定方法において、
前記変数ｂは、前記２次電池に複数の周波数の交流電流を通電し、該交流電流の各周波数における前記２次電池のインピーダンスを計測する実験を事前に実施しておくことによりあらかじめ得られた前記２次電池のインピーダンスの周波数特性の計測データに基づいて、前記劣化速度ｙ’との相関性を有するものとして選出された所定種類の特徴量であり、且つ、該所定種類の特徴量は、その値が、１つ以上の特定の周波数における前記２次電池のインピーダンスにより規定される特徴量であり、
前記劣化速度ｙ’の各推定時点において、前記２次電池に通電する交流電流は、前記特定の周波数の成分を少なくとも含む交流電流であることを特徴とする２次電池の容量維持率推定方法。。
【請求項４】
請求項３記載の２次電池の容量維持率推定方法において、
前記２次電池のインピーダンスの実数成分、虚数成分、絶対値、及び位相のそれぞれを基本パラメータと定義したとき、前記変数ｂは、単一の周波数におけるいずれかの前記基本パラメータにより値が規定される特徴量と、２つの周波数の一方の周波数におけるいずれかの前記基本パラメータ及び他方の周波数におけるいずれかの前記基本パラメータから成る２つの基本パラメータの値の和、差、積及び比のいずれかによって値が規定される特徴量とのうちから選出されたいずれかの特徴量であることを特徴とする２次電池の容量維持率推定方法。
【請求項５】
請求項１記載の２次電池の容量維持率推定方法において、
前記劣化速度ｙ’の各推定時点において、前記変数ｂの値の経時変化を近似する次式（Ｂ１）又は（Ｂ２）に基づいて、当該推定時点における前記変数ｂの値を推定するステップをさらに備えることを特徴とする２次電池の容量維持率推定方法。
ｂ(t)＝exp(ｕ(t))−ｐ    ……（Ｂ１）
ｂ(t)＝(ｕ(t))1/q−ｐ    ……（Ｂ２）
ただし、
ｔ：２次電池の初期状態からの経過時間
ｂ(t)：経過時間がｔである時刻での変数ｂの値
ｕ(t)：ｕ(t)＝α・ｔ＋βにより表される線形関数
exp(  )：自然対数の底の指数関数
α：経過時間がｔである時刻又はその直前での２次電池の使用条件に応じて値が決定される変数
β：２次電池の初期状態での使用条件に応じて値が決定される変数
ｐ，ｑ：あらかじめ定められた所定値
【請求項６】
２次電池の容量維持率を推定する方法であって、
前記２次電池に複数の周波数の交流電流を通電し、該交流電流の各周波数における前記２次電池のインピーダンスを計測する実験を事前に実施することにより得られた該インピーダンスの周波数特性の計測データに基づいて、前記容量維持率の時間的変化率である劣化速度と相関性を有する特徴量をあらかじめ選出しておき、前記２次電池の容量維持率を推定するとき、前記特徴量の値に基づいて前記劣化速度を推定し、該劣化速度の推定値を用いて前記容量維持率を推定することを特徴とする２次電池の容量維持率推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、２次電池の容量維持率を推定する方法に関する
【背景技術】
【０００２】
２次電池の容量劣化（２次電池が満充電状態で貯蔵し得る電気量の低下）の度合いを示す指標としての容量維持率を推定する手法としては、所謂、ルート則（１／２乗則）を用いる手法が一般に知られている。また、例えば特許文献１、２に見られる手法も提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−２２０９００号公報
【特許文献２】特開２０１１−３８８５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
リチウムイオン電池等の２次電池の容量劣化の要因の１つは、内部の電極表面に界面被膜が生成されることである。この界面被膜の生成速度は経過時間の平方根にほぼ比例し、これを利用して、２次電池の容量維持率を推定する手法がルート則である。
【０００５】
この手法では、任意の時点における２次電池の容量維持率の時間的変化率（単位時間当たりの容量維持率の低下量。以降、「劣化速度」又は「容量維持率の低下速度」ということがある）は、２次電池の使用条件（温度、充電率（ＳＯＣ：State of Charge）等）と、当該時点での容量維持率とに応じて規定されるものとなる。
【０００６】
一方、２次電池の劣化速度は、界面被膜の生成の他の要因の影響も受ける。特に、本願発明者の各種実験、検討によって、２次電池の内部状態の劣化の程度が劣化速度に大きな影響を及ぼすことが判明した。
【０００７】
従って、２次電池の容量維持率を高精度に推定するためには、２次電池の内部状態の劣化の程度を適切に反映させて劣化速度を推定し、この劣化速度の推定値を用いて容量維持率を推定することが望ましいと考えられる。
【０００８】
しかしながら、従来のルート則を用いる手法、あるいは、特許文献１に見られ手法では、２次電池の内部状態の劣化の程度に起因する劣化速度の変化の影響が配慮されていないため、容量維持率の推定精度をより一層高めることが困難なものとなっていた。
【０００９】
また、特許文献２に見られる手法では、容量維持率の推定を行なう時点で、２次電池のインピーダンスの周波数特性を計測し、その計測データを利用して容量維持率の推定を行なうことは行なわれている。
【００１０】
しかるに、この特許文献２の手法では、劣化速度を推定することなく、上記計測データに基づき決定した単一の特徴周波数に応じて直接的に容量維持率を推定するようにしている。このため、２次電池の仕様や使用条件等によって、２次電池の内部状態の劣化の程度に起因する劣化速度の変化の影響が、容量維持率の推定値に適切に反映されずに、該推定値の精度が低下する虞がある。
【００１１】
本発明は、かかる背景に鑑みてなされたものであり、２次電池の容量維持率を高精度で推定することができる推定方法を提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明の２次電池の容量維持率推定方法は、２次電池の容量維持率を推定する方法であって、前記２次電池の容量維持率の時間的変化率である劣化速度ｙ’を時間間隔をおいて次式（Ａ）により逐次推定する第１ステップと、該劣化速度ｙ’の推定値を用いて前記２次電池の容量維持率を推定する第２ステップとを備えることを特徴とする（第１発明）。
【００１３】
ｙ’＝(ｃ／(１−ｙ)(1/c)-1)・(ａ・ｂ＋ｚ)    ……（Ａ）
ただし、
ｃ：あらかじめ定めた所定値
ｙ：劣化速度ｙ’の今回の推定時点での容量維持率の推定値
ａ，ｚ：劣化速度ｙ’の今回の推定時点又はその直前における２次電池の使用条件に応じて値が決定される変数
ｂ：２次電池の初期状態からの経過時間ｔの増加に伴い、値が増加していくように決定される変数
なお、本発明において、「今回の推定時点」というのは、前記劣化速度ｙ’を新たに推定しようとしている時点を意味する。
【００１４】
また、２次電池の「初期状態」は、２次電池の容量劣化が未だ生じていない（２次電池が新品である）とみなせる状態を意味し、例えば２次電池の製造直後、工場出荷時等の状態である。
【００１５】
また、２次電池の「使用条件」というのは、２次電池の容量劣化に影響を及ぼす１つ又は複数種類の状態量（例えば、２次電池の温度、ＳＯＣ、通電電流等）が、どのような状態であるか（例えば各種類の状態量の実際の値がどのような範囲内の値であるか等）を示す条件を意味する。この場合、「使用条件」に関する「状態量」の種別は、２次電池の容量劣化への影響度合いが高いものを、実験等に基づいてあらかじめ選定しておけばよい。
【００１６】
本願発明者の各種実験、検討によって、前記式（Ａ）の変数ｃ，ａ，ｚ，ｂを、前記式（Ａ）のただし書きで規定する如く決定するようにすることで、前記劣化速度ｙ’（容量維持率の低下速度）を、式（Ａ）の演算によって精度よく推定できることが見出された。
【００１７】
ここで、式（Ａ）の導出について説明しておく。
【００１８】
従来のルート則によれば、２次電池の初期状態からの経過時間がｔであるときの容量維持率ｙは、次式（１）により与えられる。
【００１９】
ｙ＝１−ｋ・ｔ0.5    ……（１）
本願発明者は、まず、２次電池の容量劣化に対する界面被膜以外の各種要因を容量維持率ｙに反映させるために、式（１）をより一般化して、容量維持率ｙが、式（１）のｔの指数「０．５」を一般変数「ｃ」に置き換えて得られる次式（２）の形で与えられるものとみなした。なお、「ｃ」の値は、２次電池に固有の所定の定数値である。
【００２０】
ｙ＝１−ｋ・ｔc    ……（２）
この式（２）を時間微分することによって、次式（３）が得られる。また、式（２）を変形することで、次式（４）が得られる、。
【００２１】
ｙ’＝−ｋ・ｃ・ｔc-1    ……（３）
ｔ＝((１−ｙ)／ｋ)1/c    ……（４）
上記式（４）を式（３）に代入して整理することで、次式（５）が得られる。
【００２２】
ｙ’＝−(ｃ／(１−ｙ)(1/c)-1)・ｋ1/c    ……（５）
本願発明者は、この式（５）に、２次電池の内部状態の劣化の程度の影響と、２次電池の使用条件の影響とを反映させるために、式（５）のｋ1/cが、当該内部状態の劣化の程度と使用条件とに依存するものとみなして、ｋ1/cを次式（６）により置き換えた。
【００２３】
ｋ1/c＝ａ・ｂ＋ｚ    ……（６）
この場合、変数ａ，ｚは、２次電池の使用条件に依存する値を有する変数であり、変数ｂは２次電池の内部状態の劣化の程度に依存する値を有する変数である。
【００２４】
この式（６）を前記式（５）に代入することで、前記式（Ａ）が得られる。この式（Ａ）を下記に改めて示す。
【００２５】
ｙ’＝(ｃ／(１−ｙ)(1/c)-1)・(ａ・ｂ＋ｚ)    ……（Ａ）
この式（Ａ）により、任意の時点での劣化速度ｙ’は、その時点での容量維持率と、変数ａ，ｂ，ｚの値とから推定できることとなる。
【００２６】
この場合、ａ，ｂ，ｚのうち、ａ，ｚの値を、２次電池の使用条件に適切に依存させるように設定し、また、ｂの値を、２次電池の内部状態の劣化の程度に適切に依存させるように設定することで、式（Ａ）に基づいて、高精度に劣化速度ｙ’を推定することができることが本願発明者の各種実験、検討により確認された。
【００２７】
さらに、ｂの値を、２次電池の内部状態の劣化の程度に適切に依存させるように設定した場合、結果的には、ｂの値は、２次電池の初期状態からの経過時間ｔの増加に伴い、増加していく傾向があることも本願発明者の各種実験、検討により確認された。
【００２８】
そこで、第１発明では、前記式（Ａ）により劣化速度ｙ’を推定することとした。
【００２９】
かかる第１発明によれば、前記第１ステップにおいて、式（Ａ）により劣化速度ｙ’を推定することで、その推定値を高精度に実際の劣化速度に合致させることができる。そして、前記第２ステップでは、この劣化速度ｙ’の推定値を用いて容量維持率を推定することで、すなわち、容量維持率の推定値が、劣化速度ｙ’の推定値に一致する低下速度で低下していくように該容量維持率を推定することで、該容量維持率を高精度に推定することができる。
【００３０】
前記第１発明では、前記変数ｂの値は、より具体的には、例えば次のような手法によって適切に決定できる。
【００３１】
その１つの手法では、前記劣化速度ｙ’の各推定時点において、前記２次電池に単一又は複数の周波数の交流電流を通電し、該交流電流の各周波数における前記２次電池のインピーダンスを計測するステップと、該インピーダンスの計測データに基づいて前記変数ｂの値を決定するステップとをさらに備える（第２発明）。
【００３２】
すなわち、変数ｂの値は、一般に、ある単一又は複数の周波数（交流周波数）における２次電池のインピーダンスと密接に関連する。そこで、第２発明では、前記インピーダンスを計測して、このインピーダンスの計測データに基づいて前記変数ｂの値を決定する。
【００３３】
これにより、劣化速度を精度よく推定する上で適切な変数ｂの値を高い信頼性で決定できる。
【００３４】
なお、第２発明において、２次電池に通電する「交流電流」は、正弦波の交流電流でよいことはもちろん、複数の周波数成分を含む交流電流（例えば矩形波電流、三角波電流等）であってもよい。
【００３５】
上記第２発明では、前記変数ｂは、前記２次電池に複数の周波数の交流電流を通電し、該交流電流の各周波数における前記２次電池のインピーダンスを計測する実験を事前に実施しておくことによりあらかじめ得られた前記２次電池のインピーダンスの周波数特性の計測データに基づいて、前記劣化速度ｙ’との相関性を有するものとして選出された所定種類の特徴量であり、且つ、該所定種類の特徴量は、その値が、１つ以上の特定の周波数における前記２次電池のインピーダンスにより規定される特徴量であり、前記劣化速度ｙ’の各推定時点において、前記２次電池に通電する交流電流は、前記特定の周波数の成分を少なくとも含む交流電流であることが好ましい（第３発明）。
【００３６】
かかる第３発明によれば、前記変数ｂは、事前の実験により得られた２次電池のインピーンダンスの周波数特性の計測データに基づいて、上記の如く選出された所定種類の特徴量であるので、劣化速度ｙ’と相関性が極力高いものとなるような当該変数ｂの種別を適切に選出できる。
【００３７】
この場合、該所定種類の特徴量は、その値が、１つ以上の特定の周波数における前記２次電池のインピーダンスにより規定される特徴量であるので、前記劣化速度ｙ’の各推定時点での変数ｂの値は、その推定時点での上記特定の周波数におけるインピーダンスの計測データに応じて規定されるものとなる。
【００３８】
従って、前記劣化速度ｙ’の各推定時点において、前記２次電池に通電する交流電流の周波数を、前記特定の周波数とすることで、該推定時点でのインピーダンスの計測データに基づいて、劣化速度ｙ’を推定するために好適な変数ｂの値を適切に決定することができる。
【００３９】
上記第３発明では、前記２次電池のインピーダンスの実数成分、虚数成分、絶対値、及び位相のそれぞれを基本パラメータと定義したとき、前記変数ｂは、単一の周波数におけるいずれかの前記基本パラメータにより値が規定される特徴量と、２つの周波数の一方の周波数におけるいずれかの前記基本パラメータ及び他方の周波数におけるいずれかの前記基本パラメータから成る２つの基本パラメータの値の和、差、積及び比のいずれかによって値が規定される特徴量とのうちから選出されたいずれかの特徴量であることが好ましい（第４発明）。
【００４０】
なお、この第４発明において、前記一方の周波数におけるいずれかの前記基本パラメータと、前記他方の周波数におけるいずれかの前記基本パラメータとは、互いに異なる種類の基本パラメータであってもよい。
【００４１】
この第４発明によれば、種々様々な仕様の２次電池に関して、劣化速度ｙ’との相関性が高い好適な種類の特徴量を、前記変数ｂとして選出することができる。
【００４２】
前記第１発明では、前記変数ｂの値を決定する手法として、前記第２〜第４発明の手法以外に、例えば次のような手法を採用することもできる。
【００４３】
すなわち、前記劣化速度ｙ’の各推定時点において、前記変数ｂの値の経時変化を近似する次式（Ｂ１）又は（Ｂ２）に基づいて、当該推定時点における前記変数ｂの値を推定するステップをさらに備える（第５発明）。
ｂ(t)＝exp(ｕ(t))−ｐ    ……（Ｂ１）
ｂ(t)＝(ｕ(t))1/q−ｐ    ……（Ｂ２）
ただし、
ｔ：２次電池の初期状態からの経過時間
ｂ(t)：経過時間がｔである時刻での変数ｂの値
ｕ(t)：ｕ(t)＝α・ｔ＋βにより表される線形関数
exp(  )：自然対数の底の指数関数
α：経過時間がｔである時刻又はその直前での２次電池の使用条件に応じて値が決定される変数
β：２次電池の初期状態での使用条件に応じて値が決定される変数
ｐ，ｑ：あらかじめ定められた所定値
ここで、本願発明者による実験、検討によれば、前記式（Ａ）における変数ｂの値は、２次電池の使用条件が一定に維持される場合、経過時間ｔの増加に伴い、指数関数又はべき乗関数に近似する形態で増加していく。
【００４４】
従って、変数ｂの値の経時変化を、前記式（Ｂ１）又は（Ｂ２）により精度よく近似することができる。この場合、前記式（Ｂ１）における変数ｐ、あるいは、前記式（Ｂ２）における変数ｐ，ｑを所定値とすると共に、線形関数ｕ(t)に関する変数α，βの値を２次電池の使用条件に適切に依存させるように設定することで、式（Ｂ１）又は（Ｂ２）に基づいて、劣化速度ｙ’を精度よく推定することを可能とする好適な変数ｂの値を決定できる。
【００４５】
このため、第５発明によれば、劣化速度ｙ’の各推定時点での２次電池のインピーダンスを計測せずとも、変数ｂの値を適切に決定して、精度よく劣化速度ｙ’を推定できる。
【００４６】
補足すると、前記変数αは、前記線形関数ｕ(t)の値の単位時間当たりの変化量を意味し、この変数αの値に応じて変数ｂの値の単位時間当たりの増加量が規定されるものとなる。また、前記変数βは、２次電池の初期状態、すなわち、ｔ＝０の時刻での状態における線形関数ｕ(t)の値を意味し、このβの値に応じて２次電池の初期状態での変数ｂの値が規定されるものとなる。
【００４７】
従って、２次電池の初期状態の時刻（ｔ＝０の時刻）においては、該初期状態での２次電池の使用条件に応じて変数βの値を決定することで、そのβの値が、該初期状態での線形関数ｕ(t)の値として得られることとなる。そして、その線形関数ｕ(t)の値から、前記式（Ｂ１）又は（Ｂ２）により該初期状態での変数ｂの値を決定できる。
【００４８】
また、２次電池の初期状態の時刻以降の各推定時点においては、該推定時点又はその直前における２次電池の使用条件に応じて変数αの値を決定することで、そのαの値が、該推定時点での線形関数ｕ(t)の単に時間当たりの変化量（変化速度）として得られることとなる。
【００４９】
このため、決定したαの値に、今回の推定時点とその１つ前の推定時点（前回の推定時点）との間の時間間隔を乗じることで、その時間間隔での線形関数ｕ(t)の値の変化量を決定でき、さらにその変化量を、前回の推定時点でのｕ(t)の値に加算することで、今回の推定時点での線形関数ｕ(t)の値を決定できる。そして、その線形関数ｕ(t)の値から、前記式（Ｂ１）又は（Ｂ２）により該初期状態での変数ｂの値を決定できる。
【００５０】
また、本発明の２次電池の容量維持率推定方法は、２次電池の容量維持率を推定する方法であって、前記２次電池に複数の周波数の交流電流を通電し、該交流電流の各周波数における前記２次電池のインピーダンスを計測する実験を事前に実施することにより得られた該インピーダンスの周波数特性の計測データに基づいて、前記容量維持率の時間的変化率である劣化速度と相関性を有する特徴量をあらかじめ選出しておき、前記２次電池の容量維持率を推定するとき、前記特徴量の値に基づいて前記劣化速度を推定し、該劣化速度の推定値を用いて前記容量維持率を推定することを特徴とする（第６発明）。
【００５１】
この第６発明によれば、事前の実験により得られた２次電池のインピーンダンスの周波数特性の計測データに基づいて、前記特徴量を選出するので、劣化速度と相関性が極力高いものとなるような種類の特徴量を選出できる。そして、２次電池の容量維持率を推定するにあたっては、この特徴量の値に基づいて前記劣化速度を推定するので、該劣化速度を精度よく推定することができる。
【００５２】
さらに、この劣化速度の推定値を用いて容量維持率を推定することで、すなわち、容量維持率の推定値が、劣化速度ｙ’の推定値に一致する低下速度で低下していくように該容量維持率を推定することで、該容量維持率を高精度に推定することができる。
【図面の簡単な説明】
【００５３】
【図１】本発明の第１実施形態において、２次電池の容量維持率を推定するためのシステムの構成を示す図。
【図２】２次電池のインピーダンスの周波数特性の例を示す図。
【図３】変数ｂの値を規定する特徴量の候補の例を示す図。
【図４】２次電池の容量維持率の経時変化とインピーダンスの周波数特性との関係を例示する図。
【図５】図１に示す測定器の劣化速度推定部及び容量維持率推定部による処理を示すフローチャート。
【図６】第１実施形態の実施例における２次電池のインピーダンスの周波数特性の例を示す図。
【図７】第１実施形態の実施例における容量維持率の推定値の経時変化を示す図。
【図８】本発明の第２実施形態において、２次電池の容量維持率を推定するためのシステムの構成を示す図。
【図９】変数ｂの値の経時変化の例を示す図。
【図１０】変数ｂの値から得られらる線形関数ｕの値の経時変化を示す図。
【図１１】図８に示す測定器の劣化速度推定部及び容量維持率推定部による処理を示すフローチャート。
【図１２】図１１のＳＴＥＰ１３の処理を説明するための図。
【発明を実施するための形態】
【００５４】
［第１実施形態］
本発明の第１実施形態を図１〜図７を参照して以下に説明する。
【００５５】
図１を参照して、本実施形態において、２次電池１（以下、バッテリ１という）の容量維持率を推定するためのシステム（以下、容量維持率推定システム）は、バッテリ１の一対の出力端子（正極及び負極）間の電圧（以降、バッテリ電圧という）を検出する電圧センサ２と、バッテリ１に流れる電流（充電電流又は放電電流。以降、バッテリ電流という）を検出する電流センサ３と、バッテリ１の温度（以降、バッテリ温度という）を検出する温度センサ４と、これらのセンサ２〜４による検出データを用いてバッテリ１の容量維持率等の推定処理を行う測定器５とを備える。
【００５６】
また、バッテリ１の出力端子間には、スイッチ６を介して交流電源７が接続され、該交流電源７からバッテリ１に、適宜、交流電流を通電することが可能となっている。
【００５７】
さらに、バッテリ１の出力端子間には、スイッチ８，１０をそれぞれ介して充電器９、負荷１１が接続され、充電器９によりバッテリ１を適宜、充電することと、バッテリ１から適宜、負荷１１に放電することとが可能となっている。
【００５８】
測定器５は、ＣＰＵ、ＲＡＭ、ＲＯＭ等を含む電子回路ユニット、あるいは、パーソナルコンピュータ等のコンピュータにより構成される。この測定器５は、実装されるプログラムを実行することにより実現される主な機能として、バッテリ１のＳＯＣ（State of Charge）を推定するＳＯＣ推定部２１と、バッテリ１の劣化速度を推定する劣化速度推定部２２と、バッテリ１の容量維持率を推定する容量維持率推定部２３とを備える。
【００５９】
ＳＯＣ推定部２１は、例えば電流センサ３によるバッテリ電流の検出データ（検出値）を逐次積算した値（詳しくは、バッテリ１の満充電状態等、ＳＯＣが既知の値となっている状態からのバッテリ電流の積算値）を基に、ＳＯＣを逐次推定する。なお、バッテリ１のＳＯＣの推定手法は、公知の任意の手法を採用することができる。
【００６０】
劣化速度推定部２２は、本発明を適用した手法によって、後述の如く劣化速度を推定する。
【００６１】
容量維持率推定部２３は、劣化速度推定部２２により決定された劣化速度の推定値を用いて、バッテリ１の容量維持率を推定する。
【００６２】
<<事前の準備処理>>
本実施形態の容量維持率推定システムによりバッテリ１の容量維持率の推定を行なうにあたって、以下のように、事前に準備処理が実行される。
【００６３】
その準備処理では、バッテリ１の初期状態から任意に設定した時間間隔で、バッテリ１の容量維持率を実測することと、バッテリ１に交流電流を通電して該バッテリ１のインピーダンスを実測することを含む測定処理が行なわれる。また、この測定処理は、バッテリ温度、バッテリ電流、ＳＯＣ等の種々様々の使用条件の基で行なわれる。
【００６４】
容量維持率及びインピーダンスの実測を行なう各時点（以降、当該時点をバッテリ状態測定時という）におけるバッテリ１の容量維持率の実測は次のように行なわれる。すなわち、バッテリ１の放電及び充電を既定の仕方で行うことで、バッテリ１の容量（満充電容量）を計測し、その計測値の、初期状態での容量値に対する比率を求めることで、当該比率が容量維持率として実測される。なお、バッテリ１の初期状態は、バッテリ１の新品状態（バッテリ１の製造直後や工場出荷時の状態等、バッテリ１の劣化が生じていないとみなし得る状態）である。
【００６５】
バッテリ１の容量劣化は、少なくとも初期状態からの経過時間の増加に伴い進行するので、実測される容量維持率は、経過時間の増加に伴い低下していく。
【００６６】
また、各バッテリ状態測定時におけるバッテリ１のインピーダンスの実測は次のように行なわれる。すなわち、あらかじめ定めた複数の周波数のそれぞれの周波数で、交流電源からバッテリ１に交流電流を通電した状態で、バッテリ電圧及びバッテリ電流が測定される。
【００６７】
そして、該バッテリ電圧及びバッテリ電流の振幅値の測定値と、該バッテリ電圧及びバッテリ電流の間の位相差の測定値とから、実数成分（抵抗成分）と虚数成分（リアクタンス成分）とから成るインピーダンス（複素インピーダンス）が測定される。以降、各バッテリ状態測定時における複数の周波数のそれぞれに対応するインピーダンスの測定値の全体をインピーダンスの周波数特性の測定データという。
【００６８】
なお、バッテリ１に通電する交流電流は、正弦波に限らず、矩形波の交流電流等、複数の周波数成分を含む交流電流であってもよい。
【００６９】
図２は、上記測定により得られるインピーダンスの周波数特性の測定データを、Cole-Coleプロットの形式で例示した図である。図２の実軸（横軸）は、インピーダンスの実数成分の値を示す軸であり、虚軸（縦軸）は、インピーダンスの虚数成分の値を示す軸である。図２中の白丸もしくは黒丸の各丸点が、各周波数でのインピーダンスの測定点を示している。
【００７０】
ここで、バッテリ１の容量維持率の推移と、インピーダンスの周波数特性との関係等に関して図４を参照して説明しておく。
【００７１】
図４中の４Ａの図は、同一仕様の２つのバッテリ１，１についての容量維持率の実測値の経時的な推移を例示したものである。黒点が一方のバッテリ１の容量維持率の実測値、白点が他方のバッテリ１の容量維持率の実測値である。この図示例では、バッテリ１の初期状態からの経過時間ｔがｔ＝ｔxとなる時刻以前では、２つのバッテリ１，１の使用条件を互いに異らせ、ｔ＝ｔxの時刻以後では、２つのバッテリ１，１の使用条件を同じにしている。
【００７２】
また、図４中の４Ｂの図は、バッテリ１の初期状態の時刻である初期時刻（ｔ＝０の時刻）でのインピーダンスの周波数特性の測定データを示し、４Ｃの図はｔ＝ｔxの時刻（例えば３０００時間程度の経過時間の時刻）でのインピーダンスの周波数特性の測定データを示している。この場合、４Ｂ、４Ｃのそれぞれの図中の実線が、４Ａの図の黒点の実測値に対応するバッテリ１のインピーダンスの周波数特性の測定データを示し、４Ｂ、４Ｃのそれぞれの図中の破線が、４Ａの図の白点の実測値に対応するバッテリ１のインピーダンスの周波数特性の測定データを示している。
【００７３】
なお、４Ｂ、４Ｃの図では、インピーダンスの周波数特性の測定データを、インピーダンスの測定点を連結した状態で示している。
【００７４】
４Ａの図を参照して判るように、２つのバッテリ１，１のいずれについても、容量維持率は、経過時間ｔの増加に伴い、低下していく。すなわち、経過時間ｔの増加に伴い、バッテリ１の容量劣化が進行していく。
【００７５】
この場合、図４に示す例では、ｔ＝ｔxの時刻以前での２つのバッテリ１，１の使用条件は互いに異なるものの、ｔ＝ｔxの時刻での容量維持率が２つのバッテリ１，１で互いにほぼ同一となっている。また、ｔ＝ｔxの時刻以後は、バッテリ１，１の使用条件は互いに同じであるものの、容量維持率の時間的変化率（低下速度）が２つのバッテリ１，１で互いに異なるものとなっている。
【００７６】
このように、ある時刻でのバッテリ１，１の容量維持率が同じであっても、バッテリ１，１のその後の容量維持率の低下速度が異なるものとなる場合がある。
【００７７】
一方、インピーダンスの周波数特性は、初期時刻（ｔ＝０の時刻）では、２つのバッテリ１，１で互いにほぼ一致する。また、いずれのバッテリ１，１についても、ｔ＝ｔxの時刻でのインピーダンスの周波数特性は、初期状態から変化している。
【００７８】
このように、バッテリ１の内部状態を示すインピーダンスは、経過時間ｔの増加に伴い（バッテリ１の容量劣化の進行に伴い）変化する。
【００７９】
この場合、図４に示す例では、ｔ＝ｔxの時刻でのバッテリ１，１のインピーダンスの周波数特性は互いに顕著に異なるものとなっている。これは、時刻ｔx以前でのバッテリ１，１の使用条件が互いに異なることに起因するものと考えられる。
【００８０】
そして、ｔ＝ｔxの時刻以後の２つのバッテリ１，１の使用条件は互いに同じであるものの、ｔ＝ｔxの時刻以後の容量維持率の低下速度は、２つのバッテリ１，１で互いに異なっている。
【００８１】
このことから、任意の経過時間の時点でのバッテリ１のインピーダンスの周波数特性は、その後のバッテリ１の容量維持率の低下速度、すなわち、劣化速度ｙ’と高い相関性を有すると考えられる。
【００８２】
事前の準備処理では、以上の事項を踏まえて、各バッテリ状態測定時におけるバッテリ１の容量維持率の推移の測定データと、インピーダンスの周波数特性の測定データとに基づいて、バッテリ１の容量維持率の低下速度（劣化速度）と高い相関性を有する所定種類の特徴量（インピーダンスの周波数特性に係わる特徴量）を、前記式（Ａ）の変数ｂとして選出することが行なわれる。
【００８３】
ここで、インピーダンスの実数成分、虚数成分、絶対値（＝実数成分の２乗値と虚数成分の２乗値との和の平方根）、及び位相をそれぞれ、ｘ、ｙ、Ｉ、θと表記し、これらを基本パラメータと呼ぶ。そして、図２に示す如く、任意の経過時間ｔの時刻でのインピーダンスの周波数特性の測定データにおいて、任意の１つの周波数に対応するインピーダンスの測定値（例えば図２の測定点Ｐのインピーダンス）の各基本パラメータｘ、ｙ、Ｉ、θの値をそれぞれｘ(t)、ｙ(t)、Ｉ(t)、θ(t)と表記する。
【００８４】
また、任意の２つの周波数にそれぞれ対応するインピーダンスの測定値のうちの一方の測定値（例えば図２の測定点Ｐ１のインピーダンス）の各基本パラメータｘ、ｙ、Ｉ、θの値をそれぞれｘ1(t)、ｙ1(t)、Ｉ1(t)、θ1(t)と表記し、他方の測定値（例えば図２の測定点Ｐ２のインピーダンス）の各基本パラメータｘ、ｙ、Ｉ、θの値をそれぞれｘ2(t)、ｙ2(t)、Ｉ2(t)、θ2(t)と表記する。
【００８５】
本実施形態では、任意の経過時間ｔにおいて、１つの周波数でのインピーダンスについてのｘ(t)、ｙ(t)、Ｉ(t)、θ(t)のいずれかにより値が規定されると特徴量と、２つの周波数でのインピーダンスについてのｘ1(t)，ｙ1(t)、Ｉ1(t)、θ1(t)のいずれかと、ｘ2(t)，ｙ2(t)、Ｉ2(t)、θ2(t)のいずれかとの組から、和、差、積及び比のいずれかの演算（四則演算）により値が規定される特徴量とのそれぞれが、変数ｂの値を規定する所定種類の特徴量の候補とされる。すなわち、変数ｂの値は、上記候補のうちのいずれかの特徴量の値の関数値として決定されるものとされる。
【００８６】
そして、上記候補のうち、ある１つの特定周波数に対応する当該特徴量、又は、ある２つの特定周波数に対応する当該特徴量が、変数ｂの値を規定する所定種類の特徴量として選出される。
【００８７】
具体的には、本実施形態では、例えば図３に示すように、インピーダンスの各基本パラメータｘ、ｙ、Ｉ、θのうちの任意の１つを代表的にＳと表記したとき、単一の周波数でのインピーダンスのＳの値Ｓ(t)の初期状態からの変化量（＝Ｓ(t)−Ｓ(0)。以降、特徴量候補１という）と、２つの周波数でのＳの値Ｓ1(t)，Ｓ2(t)の和の初期状態からの変化量（＝(Ｓ1(t)＋Ｓ2(t))−(Ｓ1(0)＋Ｓ2(0))。以降、特徴量候補２という）と、２つの周波数でのＳの値Ｓ1(t)，Ｓ2(t)の差の初期状態からの変化量（＝(Ｓ1(t)−Ｓ2(t))−(Ｓ1(0)−Ｓ2(0))。以降、特徴量候補３という）と、２つの周波数でのＳの値Ｓ1(t)，Ｓ2(t)の積の初期状態からの変化量（＝Ｓ1(t)・Ｓ2(t)−Ｓ1(0)・Ｓ2(0)。以降、特徴量候補４という）と、２つの周波数でのＳの値Ｓ1(t)，Ｓ2(t)のそれぞれの初期状態の変化量の積（＝（Ｓ1(t)−Ｓ1(0)）・（Ｓ2(t)−Ｓ2(0)）。以降、特徴量候補５という）と、２つの周波数でのＳの値Ｓ1(t)，Ｓ2(t)のそれぞれの初期状態からの変化量の比（＝（Ｓ1(t)−Ｓ1(0)）／（Ｓ2(t)−Ｓ2(0)）。以降、特徴量候補６という）と、２つの周波数でのＳの値Ｓ1(t)，Ｓ2(t)のそれぞれの初期状態に対する比率の差（＝（Ｓ1(t)／Ｓ1(0)）−（Ｓ2(t)／Ｓ2(0)。以降、特徴量候補７という）とが、変数ｂの値を規定する所定種類の特徴量の候補とされる。
【００８８】
ただし、この場合において、２つの周波数に関しては、Ｓ1とＳ2とが、互いに異なる種類の基本パラメータであってもよい。例えば、Ｓ1がｘ1、Ｓ2がｙ2であってもよい。
【００８９】
なお、単一の周波数に対応する特徴量候補１の種類数は、基本パラメータの種類数（４個）と、バッテリ１に通電する交流電流の周波数（インピーダンスの測定に用いる周波数）の種類数との積なる。また、２つの周波数に対応する特徴量候補２〜７のそれぞれの種類数は、S1、Ｓ2のそれぞれの種類の組み合せの総数と、バッテリ１に通電する交流電流の周波数（インピーダンスの測定に用いる周波数）のうちの２つの周波数の組合せの総数との積となる。
【００９０】
そして、種々様々な使用条件の基で得られたインピーダンスの周波数特性の全ての測定データに基づいて、上記の各候補（特徴量候補１〜７）の値と、容量維持率の低下速度（劣化速度）との相関性が解析される。
【００９１】
さらにその解析結果に基づいて、容量維持率の低下速度（劣化速度）に対して最も相関性の高い候補（特徴量候補１〜７のいずれか）が、変数ｂの値を規定する所定種類の特徴量として選出される。そして、該所定種類の特徴量の値から変数ｂの値を決定するための変換関数（変換の規則）が決定される。
【００９２】
該変換関数は、例えば、所定種類の特徴量の値を、そのまま変数ｂの値とするように決定される。あるいは、該変換関数は、変数ｂの値を既定の範囲内の値に制限する（例えば０から１の範囲内の値に正規化する）ために、当該所定種類の特徴量の値に、所定値のゲインを乗じたもの、あるいは、所定のオフセット値を加算したもの、あるいは、該ゲインの乗算及びオフセット値の加算の両方を施したものを変数ｂの値とするように決定される。
【００９３】
事前の準備処理では、以上のようにして前記式（Ａ）の変数ｂの値を規定する所定種類の特徴量が選出される。そして、このように選出した特徴量の値を算出するために必要な１つ又は２つの特定周波数の値と、該算出の演算式と、該特徴量の値を変数ｂの値に対応付ける変換関数とが前記測定器５のメモリに記憶保持される。
【００９４】
準備処理ではさらに、前記式（Ａ）の変数ｃの値（固定値）が決定されると共に、任意の経過時間ｔの時点での前記式（Ａ）の変数ａ、ｚのそれぞれの値を、当該時点（又はその直前）のバッテリ１の使用条件に応じて決定するためのマップが作成される。該マップは、変数ａ，ｚのそれぞれ毎に各別のマップであり、ａ，ｚのそれぞれの値が、バッテリ温度、ＳＯＣ、バッテリ電流等の使用条件の関数となるように作成される。
【００９５】
この場合、変数ｃの値と、上記マップとは、任意の経過時間ｔの時点において、当該ｃの値の設定値と、前記マップにより決定されるａ，ｚの値と、当該時点でのインピーダンスの周波数特性の測定データから特定されるｂの値とを用いて前記式（Ａ）により劣化速度ｙ’を算出したときに、その算出値が実際の劣化速度（容量維持率の推移の測定データにより示される劣化速度）にできるだけ精度よく合致するように決定される。そして、これらのｃの値の設定値と、マップとが前記測定器５のメモリに記憶保持される。
【００９６】
以上が、事前の準備処理である。
【００９７】
<<容量維持率の推定>>
次に、本実施形態の容量維持率推定システムによるバッテリ１の容量維持率の推定処理を説明する。
【００９８】
測定器５は、バッテリ１の初期状態での時刻（ｔ＝０の初期時刻）から電圧センサ２、電流センサ３及び温度センサ４の各センサによるバッテリ電圧、バッテリ電流、バッテリ温度の検出データを逐次取得すると共に、取得したバッテリ電流の検出データ（バッテリ電流の検出値）を逐次積算した値に基づいて、ＳＯＣ推定部２１によりＳＯＣの推定値を逐次算出する。
【００９９】
また、充電器９によるバッテリ１の充電、あるいは、負荷１１によるバッテリ１の放電が適宜行なわれる。
【０１００】
このような状況において、測定器５は、初期時刻から所定の時間間隔で、劣化速度推定部２２及び容量維持率推定部２３の処理を図５のフローチャートで示す如く実行することで、バッテリ１の劣化速度及び容量維持率を、上記時間間隔で逐次推定する。なお、バッテリ１の劣化速度及び容量維持率の推定を行なう時間間隔は一定でなくてもよい。
【０１０１】
バッテリ１の劣化速度及び容量維持率の各推定時点（図５のフローチャートの処理が実行される各時点）において、まず、ＳＴＥＰ１〜５の処理が劣化速度推定部２２により実行される。
【０１０２】
ＳＴＥＰ１では、劣化速度推定部２２は、前記式（Ａ）の変数ａ，ｚの値を、今回の推定時点（又はその直前）でのバッテリ電流及びバッテリ温度の検出データと、ＳＯＣの推定値とから、前記準備処理であらかじめ作成されたマップに基づいて決定する。
【０１０３】
次いで、ＳＴＥＰ２において、劣化速度推定部２２は、交流電源７及びスイッチ６を制御して、交流電源７からバッテリ１に交流電流を一時的に通電させる。その交流電流は、前記準備処理によって決定された特定周波数（１つ又は２つの特定周波数）の正弦波交流、又は、該特定周波数の成分を少なくとも含む交流電流である。
【０１０４】
そして、劣化速度推定部２２は、バッテリ１に交流電流を通電した状態でのバッテリ電圧及びバッテリ電流の検出データに基づいて、上記特定周波数でのバッテリ１のインピーダンスを計測する。
【０１０５】
次いで、ＳＴＥＰ３において、劣化速度推定部２２は、特定周波数でのインピーダンスの計測値から、前記準備処理により選出された特徴量の値を算出し、さらにその算出値を変数ｂの値に変換することで、今回の推定時点での変数ｂの値を決定する。なお、前記準備処理で、特徴量の値をそのまま変数ｂの値とするようにした場合には、特徴量の算出値を変数ｂの値に変換する必要なない。
【０１０６】
次いで、ＳＴＥＰ４において、劣化速度推定部２２は、ＳＴＥＰ１で決定した変数ａ，ｚの値と、ＳＴＥＰ３で決定した変数ｂの値と、前記準備処理であらかじめ決定されたｃの値と、今回の推定時点でのバッテリ１の容量維持率ｙの値と用いて、前記式（Ａ）の右辺の演算を行なうことで、今回の推定時点での劣化速度ｙ’の推定値を算出する。
【０１０７】
この場合、式（Ａ）の右辺の演算に用いる容量維持率ｙの値は、前回の推定時点（今回の推定時点の１つ前の推定時点）で決定された容量維持率ｙの推定値である。
【０１０８】
なお、劣化速度ｙ’の推定値は、負の値である。また、ｔ＝０の時点（初期時刻）での劣化速度ｙ’の推定値は、例えば、初期時刻でのＳＯＣの推定値、バッテリ温度の検出データ、バッテリ電流の検出データ等から、あらかじめ作成されたマップ（前記準備処理であらかじめ作成されたマップ）により決定される。
【０１０９】
次いで、ＳＴＥＰ５において、容量維持率推定部２３により、今回の推定時点での容量維持率ｙの推定値が算出される。
【０１１０】
この場合、次式（７）で示す如く、前回の推定時点で決定された劣化速度ｙ’の推定値ｙ'_preに、前回の推定時点と今回の推定時点との間の時間間隔Δｔを乗算してなる値（これは前回の推定時点と今回の推定時点との間の期間での容量維持率の変化量を意味する）を、前回の推定時点での容量維持率ｙ_preに加算することで、今回の推定時点での容量維持率ｙの推定値が算出される。
【０１１１】
ｙ＝ｙ_pre＋ｙ'_pre・Δｔ    ……（７）
これにより、各推定時点での容量維持率ｙの推定値は、初期時刻（ｔ＝０の時刻）での容量維持率ｙ(0)に、該推定時点までの劣化速度ｙ’の推定値を積算（積分）した値を加えた値として算出される。なお、初期時刻（ｔ＝０の時刻）での容量維持率ｙ(0)は、“１”（１００％）である。
【０１１２】
以上の処理により、各推定時点において、バッテリ１の劣化速度ｙ’と容量維持率ｙとが推定される。
【０１１３】
なお、本実施形態では、ＳＴＥＰ１〜４の処理が、本発明における第１ステップに相当し、ＳＴＥＰ５の処理が、本発明における第２ステップに相当する。
【０１１４】
以上説明した本実施形態によれば、各推定時点でのインピーダンスの計測データに基づいて決定した変数ｂの値を用いて前記式（Ａ）により劣化速度ｙ’が推定されるので、高精度に劣化速度ｙ’を推定できる。そして、この劣化速度ｙ’の推定値を用いて上記式（７）の演算により容量維持率ｙを推定するので、精度よく容量維持率ｙを推定できる。
【０１１５】
ここで、本実施形態による容量維持率の推定精度の検証結果について図６及び図７を参照して説明する。本願発明者は、本発明の一実施例として、例えば三元正極と、ハードカーボンの負極を有するバッテリセルについて、前記実施形態の手法によって、劣化速度及び容量維持率の推定を行なった。
【０１１６】
この場合、前記実施形態で説明した如く、事前の準備処理において、式（Ａ）の変数ｂの値を規定する特徴量の選出と、変数ｃの値の設定と、変数ａ，ｚの値を決定するためのマップの作成とを行なった。
【０１１７】
この実施例における変数ｂの値を規定する特徴量として、図３に示した特徴量候補のうち、インピーダンスの実数成分を基本パラメータとして用いる特徴量候補３を選出した。すなわち、図６に示す如く、２つの特定周波数ｆ1，ｆ2におけるインピーダンスの実数成分ｘ1，ｘ2の差Δｘの初期状態からの変化量を、実施例のバッテリセルに関して、変数ｂの値を規定する所定種類の特徴量として選出した。なお、図６は、実施例のバッテリセルのインピーダンスの周波数特性の測定データの一例を示している。
【０１１８】
また、実施例における式（Ａ）の変数ｃの値は、０．７に設定した。
【０１１９】
そして、バッテリセルの初期状態から所定の時間間隔で図５のフローチャートに示す処理を実行することで、劣化速度ｙ’の推定を行い、さらに、この劣化速度ｙ’の推定値を用いて容量維持率ｙの推定を行なった。
【０１２０】
図７は、当該実施例における容量維持率の推定値のデータと、該容量維持率の実測値のデータとをそれぞれ、実線、破線で示している。なお、図７には、比較のために、従来のルート則によって推定した容量維持率のデータを一点鎖線で併記した。
【０１２１】
図７に見られるように、実施例における容量維持率の推定値は、ルート則による推定値に比して、高精度に実測値に一致することが判る。
【０１２２】
［第２実施形態］
次に、本発明の第２実施形態を図８〜図１２を参照して説明する。
【０１２３】
まず、本実施形態の概略を説明する。前記第１実施形態では、劣化速度ｙ’及び容量維持率ｙの各推定時点において、特定周波数でのバッテリ１のインピーダンスを計測し、その計測データを用いて式（Ａ）の変数ｂの値を決定するようにした。
【０１２４】
これに対して本実施形態では、バッテリ１のインピーダンスを計測せずに、式（Ａ）の変数ｂの値を推定し、その推定値を用いて劣化速度ｙ’（容量維持率ｙの低下速度）を推定する。
【０１２５】
ここで、変数ｂの値は、バッテリ１の使用条件を一定に維持した場合、一般に、バッテリ１の初期状態からの経過時間ｔの増加に伴い、図９のグラフで示すように、単調に増加していく。その増加の形態は、指数関数又はべき乗関数の形態に近似する。
【０１２６】
ただし、変数ｂの値の増加率（単位時間当たりのｂの値の増加量）は、一般にバッテリ１の使用条件に応じて変化するので、それを加味して変数ｂの値を推定する必要がある。
【０１２７】
そこで、本願発明者は、変数ｂの推移を線形化することを試みた。具体的には、任意の経過時間ｔの時刻での変数ｂの値ｂ(t)に対して、次式（８ａ）又は（８ｂ）の関係を有する関数ｕ(t)を想定する。
【０１２８】
ｕ(t)＝log((ｂ(t)＋p)    ……（８ａ）
ｕ(t)＝(ｂ(t)＋ｐ)q    ……（８ｂ）
なお、log(  )は、自然対数である。また、ｐ，ｑは、ある適当な所定値（固定値）である。
【０１２９】
このとき、式（８ａ）のｐ、あるいは、式（８ｂ）のｐ，ｑの値を適切な所定値（固定値）に設定することで、変数ｂの値ｂ(t)が図９に示すような形態で経過時間ｔの増加に伴い増加すると、関数ｕ(t)は、近似的に、図１０に示すような線形関数（一次関数）、すなわち、次式（９）の形の線形関数として表される。従って、変数ｂは、式（８ａ）又は式（８ｂ）によって、線形関数ｕ(t)に変換される。
【０１３０】
ｕ(t)＝α・ｔ＋β    ……（９）
そして、式（９）を式（８ａ），（８ｂ）に代入して変形することで、前記式（Ｂ１）、（Ｂ２）が得られる。これらの式（Ｂ１），（Ｂ２）を下記に改めて示す。
【０１３１】
ｂ(t)＝exp(ｕ(t))−ｐ＝exp(α・ｔ＋β)−ｐ    ……（Ｂ１）
ｂ(t)＝(ｕ(t))1/q−ｐ＝(α・ｔ＋β)1/q−ｐ    ……（Ｂ２）
この場合、式（９）のα、βを、その値がバッテリ１の使用条件に依存する変数として、該変数α、βの値を用いて関数ｕ(t)＝α・ｔ＋βの値を決定し、そのｕ(t)の値から、式（Ｂ１）又は（Ｂ２）によりｂ(t)を決定することで、バッテリ１の使用条件に応じた変数ｂの増加率の変化を適切に加味して、該変数ｂの値を決定できる。
【０１３２】
そこで、本実施形態では、変数ｂの値を、式（Ｂ１）及び（Ｂ２）の一方、例えば（Ｂ１）に基づき推定することとした。
【０１３３】
以上を踏まえて、本実施形態を以下に具体的に説明する。なお、本実施形態の説明では、第１実施形態と同一の構成要素については、第１実施形態と同一の参照符号を用いる。
【０１３４】
図８を参照して、本実施形態における容量維持率推定システムは、劣化速度ｙ’（容量維持率の低下速度）や容量維持率の推定を行なうための検出データを取得するセンサとして、バッテリ電流を検出する電流センサ３とバッテリ温度を検出する温度センサ４とを備えると共に、第１実施形態と同様にバッテリ１の出力端子間にスイッチ８，１０をそれぞれ介して接続された充電器９、負荷１１とを備える。
【０１３５】
さらに、本実施形態の容量維持率推定システムは、電流センサ３及び温度センサ４による検出データを用いて劣化速度ｙ’及び容量維持率の推定処理を行う測定器１５とを備える。
【０１３６】
この測定器１５は、第１実施形態の測定器５と同様に、ＣＰＵ、ＲＡＭ、ＲＯＭ等を含む電子回路ユニット、あるいは、パーソナルコンピュータ等のコンピュータにより構成される。そして、測定器１５は、実装されるプログラムを実行することにより実現される主な機能として、第１実施形態と同様に、バッテリ１のＳＯＣを推定するＳＯＣ推定部２１と、バッテリ１の劣化速度を推定する劣化速度推定部２４と、バッテリ１の容量維持率を推定する容量維持率推定部２３とを備える。
【０１３７】
この場合、ＳＯＣ推定部２１及び容量維持率推定部２３の処理は、第１実施形態と同じである。一方、本実施形態における劣化速度推定部２４は、その一部の処理（式（Ａ）の変数ｂの値を決定する処理）が第１実施形態と相違する。
【０１３８】
<<事前の準備処理>>
本実施形態においても、バッテリ１の容量維持率の推定を行なうにあたって、事前に準備処理が実行される。
【０１３９】
この準備処理では、まず、第１実施形態と同様に、バッテリ温度、バッテリ電流、ＳＯＣ等の種々様々の使用条件の基で、バッテリ１の初期状態から任意に設定した時間間隔をおいて、バッテリ１の容量維持率を実測することと、バッテリ１に交流電流を通電して該バッテリ１のインピーダンスの周波数特性を実測することとが行なわれる。
【０１４０】
そして、バッテリ１の容量維持率の推移の測定データと、インピーダンスの周波数特性の測定データとに基づいて、バッテリ１の容量維持率の低下速度（劣化速度）と高い相関性を有する所定種類の特徴量を前記式（Ａ）の変数ｂの値を規定する特徴量として選出することと、式（Ａ）の変数ｃの値を決定することと、式（Ａ）の変数ａ，ｚの値をバッテリ１の使用条件に応じて決定するためのマップを作成することとが第１実施形態と同様に行なわれる。
【０１４１】
なお、変数ｃの値と、変数ａ，ｚにそれぞれ対応するマップとは、第１実施形態と同じでよい。
【０１４２】
加えて、本実施形態における準備処理では、前記式（Ｂ１）の変数ｐの値（固定値）が決定されると共に、前記式（Ｂ１）の変数α、βのそれぞれの値を、バッテリ１の使用条件に応じて決定するためのマップが作成される。該マップは、変数α、βのそれぞれ毎に各別のマップであり、α、βのそれぞれの値が、バッテリ温度、ＳＯＣ、バッテリ電流等の使用条件の関数となるように作成される。
【０１４３】
この場合、変数ｐの値と、変数α，βにそれぞれ対応するマップとは、バッテリ１の使用条件が一定に維持された状態で、当該ｐの値の設定値と、当該マップにより決定されるα，βの値とを用いて、前記式（Ｂ１）により任意の経過時間ｔの時刻での変数ｂの値ｂ(t)を算出したときに、その算出値が、バッテリ１のインピーダンス（特定周波数でのインピーダンス）の測定データから特定される変数ｂの値にできるだけ精度よく合致するように決定される。
【０１４４】
なお、変数βの値は、ｔ＝０の初期時刻での線形関数ｕ(t)の値ｕ(0)に一致し、該初期時刻での変数ｂの値ｂ(0)を規定するものである。従って、βに対応するマップは、ｔ＝０の初期時刻でのβの値（＝ｕ(0)）を決定するために用いるマップである。
【０１４５】
また、変数αの値は、任意の経過時間ｔの時刻での線形関数ｕ(t)の値の時間的変化率（増加速度）、ひいては、変数ｂの増加率を規定するものである。従って、αに対応するマップは、任意の経過時間ｔの時刻でのαの値（＝ｕ(t)の増加速度）決定するために用いるマップである。
【０１４６】
本実施形態における事前の準備処理では、以上の如く決定した変数ｃ、ｐの値の設定値と、変数ａ，ｚにそれぞれ対応するマップと、変数α，βにそれぞれ対応するマップとが、前記測定器１５のメモリに記憶保持される。
【０１４７】
以上が、本実施形態での事前の準備処理である。
【０１４８】
<<容量維持率の推定>>
次に、本実施形態の容量維持率推定システムによるバッテリ１の容量維持率の推定処理を説明する。
【０１４９】
測定器１５は、バッテリ１の初期状態での時刻（ｔ＝０の初期時刻）から電流センサ３及び温度センサ４の各センサによるバッテリ電流、バッテリ温度の検出データを逐次取得すると共に、ＳＯＣ推定部２１により第１実施形態と同様に、ＳＯＣの推定値を逐次算出する。
【０１５０】
また、充電器９によるバッテリ１の充電、あるいは、負荷１１によるバッテリ１の放電が適宜行なわれる。
【０１５１】
このような状況において、測定器１５は、初期時刻から所定の時間間隔で、劣化速度推定部２４及び容量維持率推定部２３の処理を図１１のフローチャートで示す如く実行することで、バッテリ１の劣化速度及び容量維持率を、上記時間間隔で逐次推定する。なお、第１実施形態の場合と同様に、バッテリ１の劣化速度及び容量維持率の推定を行なう時間間隔は一定でなくてもよい。
【０１５２】
バッテリ１の劣化速度及び容量維持率の各推定時点（図１１のフローチャートの処理が実行される各時点）において、まず、ＳＴＥＰ１１〜１５の処理が劣化速度推定部２４により実行される。
【０１５３】
ＳＴＥＰ１１では、劣化速度推定部２４は、前記式（Ａ）の変数ａ，ｚの値を、今回の推定時点（又はその直前）でのバッテリ電流及びバッテリ温度の検出データと、ＳＯＣの推定値とから、前記準備処理であらかじめ作成されたマップに基づいて決定する。このＳＴＥＰ１１の処理は、第１実施形態のＳＴＥＰ１の処理と同じである。
【０１５４】
次いで、ＳＴＥＰ１２において、劣化速度推定部２４は、前記式（Ｂ１）の変数α，βの値を、バッテリ電流及びバッテリ温度の検出データと、ＳＯＣの推定値とから、前記準備処理であらかじめ作成されたマップに基づいて決定する。
【０１５５】
この場合、変数αの値は、今回の推定時点（又はその直前）でのバッテリ電流及びバッテリ温度の検出データと、ＳＯＣの推定値とから、前記準備処理であらかじめ作成されたマップに基づいて決定される。
【０１５６】
一方、変数βの値（＝ｕ(0)）は、ｔ＝０の初期時刻でのバッテリ電流及びバッテリ温度の検出データと、ＳＯＣの推定値とから、前記準備処理であらかじめ作成されたマップに基づいて決定される。
【０１５７】
なお、初期時刻以外の推定時点での変数βの値（＝ｕ(0)）は、マップを使用せずに、初期時刻での変数βの値と同じ値に決定してもよい。
【０１５８】
また、ｔ＝０の初期時刻では、変数αの値によらずに、α・ｔ＝０であるから、初期時刻でのαの値は、任意の値（例えばあらかじめ定めた固定値等）でもよい。
【０１５９】
次いで、ＳＴＥＰ１３において、劣化速度推定部２４は、ＳＴＥＰ１２で決定した変数α、βの値を用いて線形関数ｕ(t)の値を決定する。
【０１６０】
具体的には、劣化速度推定部２４は、今回の推定時点が初期時刻である場合には、βの値を今回の推定時点でのu(t)の値（＝ｕ(0)）として決定する。
【０１６１】
また、今回の推定時点が初期時刻でない場合には、劣化速度推定部２４は、今回の推定時点でのＳＴＥＰ１２で決定した変数αの値に、前回の推定時点と今回の推定時点との間の時間間隔Δｔを乗算してなる値（＝α・Δｔ）を、前回の推定時点と今回の推定時点との間での線形関数ｕの値の変化量として算出する。そして、劣化速度推定部２４は、この変化量の算出値を、前回の推定時点で決定した線形関数ｕの値ｕ(t−Δｔ)に加算することにより、今回の推定時点での線形関数ｕ(t)の値を決定する。
【０１６２】
このようにｕ(t)の値を決定することで、バッテリ１の使用条件が一定に維持されている場合（各推定時点（又はその直前）での使用条件が同じである場合）には、ｕ(t)の値は、図１２のいずれかの破線ラインで示す如く、一定の増加率（＝α）で増加していくこととなる。
【０１６３】
一方、バッテリ１の使用条件がある推定時点（又はその直前）で切替わった場合には、図１２に示すように、その切替わり時点でのu(t)の値を保持したまま、変数αの値（ｕ(t)の増加率）が切替わり後の使用条件に応じた値に変更され、その変更後のαの値により示される増加率で、ｕ(t)の値が増加していくこととなる。
【０１６４】
次いで、ＳＴＥＰ１４において、劣化速度推定部２４は、ＳＴＥＰ１３で決定した線形関数ｕ(t)の値から、前記式（Ｂ１）により変数ｂの値を算出する。なお、式（Ｂ１）の変数ｐの値は、事前の準備処理で決定された設定値である。
【０１６５】
次いで、ＳＴＥＰ１５において、劣化速度推定部２４は、ＳＴＥＰ１１で決定した変数ａ，ｚの値と、ＳＴＥＰ１４で決定した変数ｂの値と、前記準備処理であらかじめ決定されたｃの値と、今回の推定時点でのバッテリ１の容量維持率ｙの値と用いて、前記式（Ａ）の右辺の演算を行なうことで、今回の推定時点での劣化速度ｙ’の推定値（＜０）を算出する。
【０１６６】
この場合、第１実施形態におけるＳＴＥＰ４の処理と同様に、式（Ａ）の右辺の演算に用いる容量維持率ｙの値は、前回の推定時点で求められた容量維持率ｙの推定値である。
【０１６７】
また、ｔ＝０の時点（初期時刻）での劣化速度ｙ’の推定値は、該初期時刻でのＳＯＣの推定値、バッテリ温度の検出データ、バッテリ電流の検出データ等からあらかじめ作成されたマップにより決定された値とされる。
【０１６８】
次いで、ＳＴＥＰ１６において、容量維持率推定部２３により、今回の推定時点での容量維持率ｙの推定値が算出される。この処理は、第１実施形態におけるＳＴＥＰ５の処理と同様に行なわれる。
【０１６９】
以上の処理により、各推定時点において、バッテリ１の劣化速度ｙ’と容量維持率ｙとが推定される。
【０１７０】
なお、本実施形態では、ＳＴＥＰ１１〜１５の処理が、本発明における第１ステップに相当し、ＳＴＥＰ１６の処理が、本発明における第２ステップに相当する。
【０１７１】
以上説明した本実施形態によれば、前記式（Ｂ１）に基づき、各推定時点での変数ｂの値を推定することで、各推定時点で、インピーダンスの計測を行なうことなく、前記式（Ａ）により劣化速度ｙ’を精度よく推定することができる。そして、この劣化速度ｙ’の推定値を用いて上記式（７）の演算により容量維持率ｙを推定するので、精度よく容量維持率ｙを推定できる。
【０１７２】
なお、第２実施形態では、前記式（Ｂ１）に基づいて、変数ｂの値を推定するようにしたが、図１１のＳＴＥＰ１４に括弧書きで示すように、前記式（Ｂ２）により変数ｂの値を算出するようにしてもよい。なお、この場合においては、式（Ｂ２）の変数ｐ，ｑの値は、事前の準備処理において、式（Ｂ１）の変数ｐと同様に所定値（固定値）に決定される。
【０１７３】
また、前記第２実施形態では、事前の準備処理において、バッテリ１のインピーダンスの計測を行い、その計測データに基づいて、変数ｂの値を規定する特徴量を選出した上で、前記式（Ｂ１）の変数ｐの値（又は式（Ｂ２）の変数ｐ，ｑの値）と、式（９）の変数α，βの値を決定するためのマップとを決定するようにしたが、インピーダンスの計測データを用いずに、式（Ｂ１）の変数ｐの値（又は式（Ｂ２）の変数ｐ，ｑの値）と、変数α，βの値を決定するためのマップとを決定するようにしてもよい。
【０１７４】
この場合には、式（Ｂ１）（又は式（Ｂ２））に基づき推定される変数ｂの値を用いて前記式（Ａ）により劣化速度ｙ’を推定したときに、その推定値が、容量維持率の推移の計測データにより示される実際の劣化速度にできるだけ合致するように、式（Ｂ１）の変数ｐの値（又は式（Ｂ２）の変数ｐ，ｑの値）と、式（９）の変数α，βの値を決定するためのマップと、式（Ａ）の変数ｃの値と、式（Ａ）の変数ａ，ｚを決定するためのマップとの全体を決定するようにすればよい。
【０１７５】
また、前記各実施形態では、式（Ａ）の変数ａ，ｚ、あるいは、式（９）の変数α，βの値をバッテリ温度、バッテリ電流及びＳＯＣの３種類の状態量の値に関する使用条件に応じて決定するようにしたが、これらの状態量のうちの１つ又は２つの状態量に関する使用条件に応じて式（Ａ）の変数ａ，ｚ、あるいは、式（９）の変数α，βの値を決定するようにしてもよい。
【符号の説明】
【０１７６】
１…バッテリ（２次電池）、２２，２４…劣化速度推定部、２３…容量維持率推定部、ＳＴＥＰ１〜４，ＳＴＥＰ１１〜１５…第１ステップ、ＳＴＥＰ５，ＳＴＥＰ１６…第２ステップ。
【図１】
2014052186.tif 000003
【図２】
2014052186.tif 000004
【図３】
2014052186.tif 000005
【図４】
2014052186.tif 000006
【図５】
2014052186.tif 000007
【図６】
2014052186.tif 000008
【図７】
2014052186.tif 000009
【図８】
2014052186.tif 000010
【図９】
2014052186.tif 000011
【図１０】
2014052186.tif 000012
【図１１】
2014052186.tif 000013
【図１２】
2014052186.tif 000014
【手続補正書】
【提出日】平成24年10月16日(2012.10.16)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】０００８
【補正方法】変更
【補正の内容】
【０００８】
しかしながら、従来のルート則を用いる手法、あるいは、特許文献１に見られる手法では、２次電池の内部状態の劣化の程度に起因する劣化速度の変化の影響が配慮されていないため、容量維持率の推定精度をより一層高めることが困難なものとなっていた。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０１０５
【補正方法】変更
【補正の内容】
【０１０５】
次いで、ＳＴＥＰ３において、劣化速度推定部２２は、特定周波数でのインピーダンスの計測値から、前記準備処理により選出された特徴量の値を算出し、さらにその算出値を変数ｂの値に変換することで、今回の推定時点での変数ｂの値を決定する。なお、前記準備処理で、特徴量の値をそのまま変数ｂの値とするようにした場合には、特徴量の算出値を変数ｂの値に変換する必要はない。