(57)【要約】【課題】二次電池の劣化状態およびそれに基づく寿命予測を正確に行う二次電池の劣化率算出方法、二次電池の寿命予測方法、二次電池の劣化率算出システムおよび二次電池の寿命予測システムを提供する。【解決手段】電池パック１０の放電量、放電レートおよび放電温度の履歴に基づき、現時点での積算放電量Ｐ、積算平均放電レートＣおよび積算平均放電温度Ｔｐを算出し、電池パック１０の放置時の放置時間、充電状態および放置温度の履歴に基づき、現時点での積算放置時間Ｕ、積算平均充電状態Ｓおよび積算平均放置温度Ｔｓを算出し、算出したこれらのデータを用い、電池パック１０の現時点劣化量Ωを算出し、算出した現時点劣化量Ωと寿命時劣化量Ωｍａｘとの比から現時点劣化率Ωｒａｔｅを算出する。【選択図】図１
2013181875.tif 000002
【特許請求の範囲】
【請求項１】
二次電池の放電量、放電レートおよび放電温度を計測するステップと、
計測した前記放電量、放電レートおよび放電温度の履歴に基づき、所定時点での積算放電量、積算平均放電レートおよび積算平均放電温度を算出するステップと、
前記二次電池の放置時の放置時間、充電状態および放置温度を計測するステップと、
計測した前記放置時間、前記充電状態および前記放置温度の履歴に基づき、前記所定時点での積算放置時間、積算平均充電状態および積算平均放置温度を算出するステップと、
算出した前記積算放電量、前記積算平均放電レート、前記積算平均放電温度、前記積算放置時間、前記積算平均充電状態および前記積算平均放置温度を用い、前記所定時点での前記二次電池の劣化量を算出するステップと、
算出した前記所定時点での前記二次電池の劣化量と予め設定した寿命時の前記二次電池の劣化量との比から前記所定時点での前記二次電池の劣化率を算出するステップと、を含むことを特徴とする二次電池の劣化率算出方法。
【請求項２】
前記二次電池の劣化量を算出するステップでは、劣化量をΩとし、前記積算放電量をＰとし、前記積算平均放電レートをＣとし、前記積算平均放電温度をＴｐとし、前記積算放置時間をＵとし、前記積算平均充電状態をＳとし、前記積算平均放置温度をＴｓとすると、下記式（１）
【数１】
2013181875.tif 000010
により、前記二次電池の劣化量を算出することを特徴とする請求項１記載の二次電池の劣化率算出方法。
【請求項３】
請求項１又は２記載の二次電池の劣化率算出方法によって現時点での前記二次電池の劣化量を算出するステップと、
計測した前記放電レート、放電温度、充電状態および放置温度の履歴に基づき、過去の任意のサンプリング期間の期間平均放電レート、期間平均放電温度、期間平均充電状態および期間平均放置温度を算出するステップと、
算出した前記期間平均放電レート、前記期間平均放電温度、前記期間平均充電状態および前記期間平均放置温度を用い、前記二次電池の劣化量の予測関数を導出するステップと、
算出した現時点の前記二次電池の劣化量および前記予測関数を用い、現時点の前記二次電池の劣化量から寿命時の前記二次電池の劣化量に至るまでの余命期間を算出するステップと、を含み、前記二次電池の寿命予測を行うことを特徴とする二次電池の寿命予測方法。
【請求項４】
前記二次電池の前記余命期間を表示するステップを更に含み、
前記サンプリング期間として複数の異なる前記サンプリング期間を設定し、
前記複数の異なるサンプリング期間に応じたそれぞれの期間平均放電レート、期間平均放電温度、期間平均充電状態および期間平均放置温度を算出し、算出した前記それぞれの期間平均放電レート、期間平均放電温度、期間平均充電状態および期間平均放置温度を用い、前記複数の異なるサンプリング期間に応じたそれぞれの予測関数を導出し、算出した現時点の前記二次電池の劣化量および導出した前記それぞれの予測関数を用い、現時点の前記二次電池の劣化量から寿命時の前記二次電池の劣化量に至るまでのそれぞれの余命期間を算出し、前記それぞれの余命期間を同時に表示することを特徴とする請求項３記載の二次電池の寿命予測方法。
【請求項５】
二次電池の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する放電履歴記録部と、
前記二次電池の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する放置履歴記録部と、
前記放電履歴記録部および前記放置履歴記録部に記録されたデータに基づき、所定時点の前記二次電池の劣化量を算出すると共に、算出した前記所定時点での前記二次電池の劣化量と予め設定した寿命時の前記二次電池の劣化量との比から前記所定時点での前記二次電池の劣化率を算出する劣化算出部と、を備えたことを特徴とする二次電池の劣化率算出システム。
【請求項６】
二次電池の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する放電履歴記録部と、
前記二次電池の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する放置履歴記録部と、
前記放電履歴記録部および前記放置履歴記録部に記録されたデータに基づき、現時点の前記二次電池の劣化量を算出する劣化算出部と、
前記放電履歴記録部および前記放置履歴記録部に記録されたデータに基づいて算出された、過去の任意のサンプリング期間の期間平均放電レート、期間平均放電温度、期間平均充電状態および期間平均放置温度を用い、前記二次電池の劣化量の予測関数を導出し、前記二次電池の劣化変遷を予測する劣化変遷予測部と、
前記劣化変遷予測部の予測に基づき、現時点の前記二次電池の劣化量から寿命時の前記二次電池の劣化量に至るまでの余命期間を算出する寿命算出部と、を備えたことを特徴とする二次電池の寿命予測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、二次電池の劣化率算出方法、二次電池の寿命予測方法、二次電池の劣化率算出システムおよび二次電池の寿命予測システムに関する。詳しくは、二次電池の劣化状態や寿命を正確に把握するための二次電池の劣化率算出方法、二次電池の寿命予測方法、二次電池の劣化率算出システムおよび二次電池の寿命予測システムに関する。
【背景技術】
【０００２】
近年、モータとエンジンとを組み合わせたハイブリッド車両（ＨＥＶ）やプラグインハイブリッド車両（ＰＨＥＶ）、モータを搭載した電気自動車（ＢＥＶ）などの電動車両には、Ｎｉ−ＭＨ電池やリチウムイオン二次電池などを直並列に並べた大容量・高出力の組電池が搭載されている。
また、電力平準化・緊急対応の観点から、太陽光発電や風力発電などの電力プラント、ビルや住宅などに二次電池が使用される。大規模サーバにおいては、無停電電源装置として二次電池が利用される。
【０００３】
上記のように利用されている二次電池は、使用に伴い劣化するものである。二次電池は、劣化すると容量・出力が低下していき、最終的には製品仕様を満足する性能が発揮できなくなり、交換する必要が生じる。
このため、二次電池の劣化を正確に計測して寿命を判定し、二次電池の交換時期を正確に把握することが求められる。
特許文献１には、二次電池の寿命を判定するために、二次電池の内部抵抗値を利用する方法が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−１６４３２９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、二次電池の内部抵抗値は、充電状態（ＳＯＣ）や電池温度によって変動するので、二次電池が一定の劣化状態でも異なる値となる場合がある。このため、二次電池の内部抵抗値を用いた場合には、二次電池の劣化を正確に算出できなかった。
【０００６】
本発明は、上記課題に鑑みたもので、その目的は、二次電池の劣化状態を正確に把握し、それに基づく寿命を正確に予測する二次電池の劣化率算出方法、二次電池の寿命予測方法、二次電池の劣化率算出システムおよび二次電池の寿命予測システムを提供することにある。
【課題を解決するための手段】
【０００７】
（１）  二次電池（例えば、後述の電池パック１０）の放電量、放電レートおよび放電温度を計測するステップ（例えば、後述のステップＳ４）と、計測した前記放電量、放電レートおよび放電温度の履歴に基づき、所定時点での積算放電量（例えば、後述の積算放電量Ｐ）、積算平均放電レート（例えば、後述の積算平均放電レートＣ）および積算平均放電温度（例えば、後述の積算平均放電温度Ｔｐ）を算出するステップ（例えば、後述のステップＳ５）と、前記二次電池の放置時の放置時間、充電状態および放置温度を計測するステップ（例えば、後述のステップＳ７）と、計測した前記放置時間、充電状態および放置温度の履歴に基づき、前記所定時点での積算放置時間（例えば、後述の積算放置時間Ｕ）、積算平均充電状態（例えば、後述の積算平均充電状態Ｓ）および積算平均放置温度（例えば、後述の積算平均放置温度Ｔｓ）を算出するステップ（例えば、後述のステップＳ８）と、算出した前記積算放電量、前記積算平均放電レート、前記積算平均放電温度、前記積算放置時間、前記積算平均充電状態および前記積算平均放置温度を用い、前記所定時点での前記二次電池の劣化量（例えば、後述の現時点劣化量Ω）を算出するステップ（例えば、後述のステップＳ９）と、算出した前記所定時点での前記二次電池の劣化量と予め設定した寿命時の前記二次電池の劣化量（例えば、後述の寿命時劣化量Ωｍａｘ）との比から前記所定時点での前記二次電池の劣化率（例えば、後述の現時点劣化率Ωｒａｔｅ）を算出するステップ（例えば、後述のステップＳ１０）と、を含むことを特徴とする二次電池の劣化率算出方法。
【０００８】
本発明者らは、二次電池の電極表面に付着する不純物量が、二次電池の劣化量に対して相関関係を示すことを明らかにした。
そこで（１）の発明は、二次電池の劣化量を、電池放電時の二次電池の電極表面に付着する不純物量に相関する劣化量分と、電池放置時の二次電池の電極表面に付着する不純物量に相関する劣化量分との和として求める。
ここで、電池放電時の二次電池の電極表面に付着する不純物量に相関する劣化量分は、積算放電量、積算平均放電レートおよび積算平均放電温度をパラメータとする関係式により求まる。
また、電池放置時の二次電池の電極表面に付着する不純物量に相関する劣化量分は、積算放置時間、積算平均充電状態および積算平均放置温度をパラメータとする関係式により求まる。
したがって、二次電池の使用状態の履歴、すなわち放電量、放電レート、放電温度、放置時間、充電状態および放置温度の履歴によって、二次電池の劣化率を算出できる。
【０００９】
（１）の発明によると、上記各種パラメータを用いて二次電池の内部抵抗値に基づかずに二次電池の劣化率を正確に算出できる。このため、二次電池の劣化状態を、劣化状態測定時の充電状態や電池温度によって変動しない真の状態として正確に把握できる。
【００１０】
（２）  前記二次電池の劣化量を算出するステップでは、劣化量をΩとし、前記積算放電量をＰとし、前記積算平均放電レートをＣとし、前記積算平均放電温度をＴｐとし、前記積算放置時間をＵとし、前記積算平均充電状態をＳとし、前記積算平均放置温度をＴｓとすると、下記式（１）
【数１】
2013181875.tif 000003
により、前記二次電池の劣化量を算出することを特徴とする（１）記載の二次電池の劣化率算出方法。
【００１１】
（２）の発明によると、二次電池の劣化量を、電池放電時の二次電池の電極表面に付着する不純物量に相関する劣化量分と、電池放置時の二次電池の電極表面に付着する不純物量に相関する劣化量分との和として、数式（１）に基づき定量的に算出できる。
【００１２】
（３）  （１）又は（２）記載の二次電池の劣化率算出方法によって現時点での前記二次電池の劣化量を算出するステップ（例えば、後述のステップＳ９）と、計測した前記放電レート、放電温度、充電状態および放置温度の履歴に基づき、過去の任意のサンプリング期間の期間平均放電レート（例えば、後述の期間平均放電レートＣ１）、期間平均放電温度（例えば、後述の期間平均放電温度Ｔｐ１）、期間平均充電状態（例えば、後述の期間平均充電状態Ｓ１）および期間平均放置温度（例えば、後述の期間平均放置温度Ｔｓ１）を算出するステップ（例えば、後述のステップＳ１１）と、算出した前記期間平均放電レート、前記期間平均放電温度、前記期間平均充電状態および前記期間平均放置温度を用い、前記二次電池の劣化量の予測関数（例えば、後述の予測関数Ωｆ（ｔ））を導出するステップ（例えば、後述のステップＳ１２）と、算出した現時点の前記二次電池の劣化量および導出した前記予測関数を用い、現時点の前記二次電池の劣化量から寿命時の前記二次電池の劣化量に至るまでの余命期間（例えば、後述の余命期間Ｔ）を算出するステップ（例えば、後述のステップＳ１３）と、を含み、前記二次電池の寿命予測を行うことを特徴とする二次電池の寿命予測方法。
【００１３】
（３）の発明では、まず（１）又は（２）の発明を用いて現時点における二次電池の劣化量を正確に算出し、次にサンプリング期間についての二次電池劣化量の予測関数を導出する。そして導出した二次電池の劣化量の予測関数を用い、現時点の二次電池の劣化量から寿命時の二次電池の劣化量に至るまでの正確な余命期間を算出する。
したがって、サンプリング期間における二次電池の使用状態に応じた二次電池の寿命を正確に予測できる。
【００１４】
（４）  前記二次電池の前記余命期間を表示するステップ（例えば、後述のステップＳ１４）を更に含み、前記サンプリング期間として複数の異なるサンプリング期間（例えば、後述のＳＰ１、ＳＰ２）を設定し、前記複数の異なるサンプリング期間に応じたそれぞれの期間平均放電レート、期間平均放電温度、期間平均充電状態および期間平均放置温度を算出し、算出した前記それぞれの期間平均放電レート、期間平均放電温度、期間平均充電状態および期間平均放置温度を用い、前記複数の異なるサンプリング期間に応じたそれぞれの予測関数を導出し、算出した現時点の前記二次電池の劣化量および導出した前記それぞれの予測関数を用い、現時点の前記二次電池の劣化量から寿命時の前記二次電池の劣化量に至るまでのそれぞれの余命期間（例えば、後述のＴ１、Ｔ２）を算出し、前記それぞれの余命期間を同時に表示することを特徴とする（３）記載の二次電池の寿命予測方法。
【００１５】
（４）の発明によると、複数の異なるサンプリング期間における二次電池の使用状態に基づき、それぞれの余命期間として二次電池の寿命を正確に予測できる。このため、指定した種々のサンプリング期間の二次電池の使用状態に応じて、二次電池の寿命を様々に予測できる。
また、それぞれの余命期間を同時に表示することにより、どのサンプリング期間における二次電池の使用状態が好ましいか一目で視認できる。
【００１６】
（５）  二次電池の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する放電履歴記録部（例えば、後述のサイクル履歴部２１）と、前記二次電池の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する放置履歴記録部（例えば、後述の放置履歴部２２）と、前記放電履歴記録部および前記放置履歴記録部に記録されたデータに基づき、所定時点の前記二次電池の劣化量（例えば、後述の現時点劣化量Ω）を算出すると共に、算出した前記所定時点での前記二次電池の劣化量と予め設定した寿命時の前記二次電池の劣化量（例えば、後述の寿命時劣化量Ωｍａｘ）との比から前記所定時点での前記二次電池の劣化率（例えば、後述の現時点劣化率Ωｒａｔｅ）を算出する劣化算出部（例えば、後述の劣化算出部３）と、を備えたことを特徴とする二次電池の劣化率算出システム。
【００１７】
（５）の発明によると、（１）の発明と同様の作用、効果を奏する。
【００１８】
（６）  二次電池の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する放電履歴記録部（例えば、後述のサイクル履歴部２１）と、前記二次電池の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する放置履歴記録部（例えば、後述の放置履歴部２２）と、前記放電履歴記録部および前記放置履歴記録部に記録されたデータに基づき、現時点の前記二次電池の劣化量（例えば、後述の現時点劣化量Ω）を算出する劣化算出部（例えば、後述の劣化算出部３）と、前記放電履歴記録部および前記放置履歴記録部に記録されたデータに基づいて算出された、過去の任意のサンプリング期間の期間平均放電レート（例えば、後述の期間平均放電レートＣ１）、期間平均放電温度（例えば、後述の期間平均放電温度Ｔｐ１）、期間平均充電状態（例えば、後述の期間平均充電状態Ｓ１）および期間平均放置温度（例えば、後述の期間平均放置温度Ｔｓ１）を用い、前記二次電池の劣化量の予測関数（例えば、後述の予測関数Ωｆ（ｔ））を導出し、前記二次電池の劣化変遷を予測する劣化変遷予測部（例えば、後述の劣化変遷予測部４）と、前記劣化変遷予測部の予測に基づき、現時点の前記二次電池の劣化量から寿命時の前記二次電池の劣化量（例えば、後述の寿命時劣化量Ωｍａｘ）に至るまでの余命期間（例えば、後述の余命期間Ｔ）を算出する寿命算出部（例えば、後述の寿命算出部５）と、を備えたことを特徴とする二次電池の寿命予測システム。
【００１９】
（６）の発明によると、（３）の発明と同様の作用、効果を奏する。
【発明の効果】
【００２０】
本発明によれば、二次電池の劣化状態を正確に把握し、それに基づく寿命を正確に予測する二次電池の劣化率算出方法、二次電池の寿命予測方法、二次電池の劣化率算出システムおよび二次電池の寿命予測システムを提供できる。
【図面の簡単な説明】
【００２１】
【図１】本発明の実施形態に係る電池寿命予測システムの概略構成を示すブロック図である。
【図２】上記実施形態に係る積算放電量と電池放電時劣化量との関係を示す図である。
【図３】上記実施形態に係る積算放置時間と電池放置時劣化量との関係を示す図である。
【図４】上記実施形態に係る電池パックの劣化量の予測関数を示す図である。
【図５】上記実施形態に係る表示部に表示される任意のサンプリング期間に応じた余命期間を示す図である。
【図６】上記実施形態に係る表示部に表示される２種のサンプリング期間に応じた余命期間を示す図である。
【図７】上記実施形態に係る電池パックの寿命予測方法を示すフローチャートである。
【発明を実施するための形態】
【００２２】
以下に図面を参照して本発明の実施形態を説明する。
【００２３】
＜実施形態＞
図１は、本実施形態に係る二次電池の劣化率算出方法、二次電池の寿命予測方法、二次電池の劣化率算出システムおよび二次電池の寿命予測システムを適用する電池寿命予測システム１の概略構成を示すブロック図である。
【００２４】
（電池寿命予測システム）
本実施形態に係る電池寿命予測システム１は、二次電池の電池パック１０に接続され、電池パック１０の劣化状態を把握すると共に電池パック１０の寿命を予測する。
【００２５】
二次電池の電池パック１０は、ハイブリッド車両（ＨＥＶ）やプラグインハイブリッド車両（ＰＨＥＶ）、モータを搭載した電気自動車（ＢＥＶ）などの電動車両あるいは太陽光発電や風力発電などの電力プラント、ビルや住宅などに使用される。
電池パック１０の蓄電池は、２極式セルであり、正極材料には、ＬｉＣｏＯ２を用い、負極材料には、グラファイトを用いる。電池パックのセパレータには、ポリプロピレン製樹脂フィルムを用い、電解液は、ＬｉＰＦ６／ＥＣ・ＤＥＣであってよい。
【００２６】
電池パック１０は、電流センサ１１と、温度センサ１２と、電圧センサ１３とを有する。
電流センサ１１および温度センサ１２は、電池パック１０の劣化状態を把握するために電池パック１０の蓄電池の放電時および放置時の電流および電池温度をそれぞれ検出する。
電圧センサ１３は、電池パック１０が使用時に性能を発揮できるか否かを判別するために電池パック１０の蓄電池の両極間の電圧を検出する。
電流センサ１１、温度センサ１２および電圧センサ１３は、それぞれ電池寿命予測システム１に接続されており、これらの出力は、電池寿命予測システム１に送信される。
【００２７】
電池寿命予測システム１は、電池使用履歴算出部２と、劣化算出部３と、劣化変遷予測部４と、寿命算出部５と、表示部６とを備える。
【００２８】
電池使用履歴算出部２は、サイクル履歴部２１と放置履歴部２２とから構成され、電池パック１０の使用履歴を算出および記録する。
電池使用履歴算出部２には、電流センサ１１および温度センサ１２の検出値が送信されると共に、電池パック１０の電池セルの情報などが予め記憶される。
サイクル履歴部２１は、電流センサ１１および温度センサ１２が検出した電池パック１０の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する。
放置履歴部２２は、電流センサ１１および温度センサ１２が検出した電池パック１０の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する。
【００２９】
劣化算出部３は、電池使用履歴算出部２に記録されたデータに基づき、電池パック１０の現時点の劣化量（以下、現時点劣化量Ωという）を算出する。また劣化算出部３は、算出した現時点劣化量Ωと、予め設定した電池パック１０の寿命時の劣化量（以下、寿命時劣化量Ωｍａｘという）との比から現時点の劣化率（以下、現時点劣化率Ωｒａｔｅという）を算出する。
【００３０】
劣化変遷予測部４は、電池使用履歴算出部２に記録されたデータに基づいて算出された、過去の任意のサンプリング期間の期間平均放電レートＣ１、期間平均放電温度Ｔｐ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１を用い、電池パック１０の劣化量の予測関数Ωｆ（ｔ）を算出し、電池パック１０の劣化変遷を予測する。
【００３１】
寿命算出部５は、劣化変遷予測部４の予測に基づき、電池パック１０の現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔを算出する。
【００３２】
表示部６は、劣化算出部３で算出された電池パック１０の現時点劣化率Ωｒａｔｅや、寿命算出部５で算出された電池パックの現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔ（寿命）などを表示する。
【００３３】
上記構成の本実施形態に係る電池寿命予測システム１は、電池パック１０の劣化量を、電池放電時の電池パック１０の電極表面に付着する不純物量に相関する劣化量分と、電池放置時の電池パック１０の電極表面に付着する不純物量に相関する劣化量分との和として求める。
ここで、電池パック１０の電極表面に付着する不純物量と、電池パック１０の劣化量とに相関関係が生じることは、本発明者らが見出したものである。
そして、上記方法により電池パック１０の現時点劣化量Ωを正確に算出する。また過去の任意のサンプリング期間についての電池パック１０の劣化量の予測関数Ωｆ（ｔ）を算出する。そして算出した電池パック１０の劣化量の予測関数Ωｆ（ｔ）を用い、電池パック１０の現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの正確な余命期間Ｔを算出する。
ここで、電池パック１０の劣化量の予測関数Ωｆ（ｔ）を用いて算出された現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔは、サンプリング期間に基づくので、サンプリング期間における電池パック１０の使用状態に応じて当該余命期間Ｔが変化する。つまり、設定するサンプリング期間における電池パック１０の使用状態が、予測される電池パック１０の寿命に反映される。
【００３４】
以下、電池寿命予測システム１の各構成を用いて、電池パック１０の寿命予測方法について詳述する。
【００３５】
（サイクル履歴部）
サイクル履歴部２１は、電流センサ１１および温度センサ１２が検出した電池パック１０の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する。
放電量、放電レートおよび放電温度を計測する際の基準となる放電期間は、電流センサ１１によって放電開始から放電終了（すなわち、充電又は放置への遷移）までの期間としてよい。放電量は、放電時電流と放電期間の積により求めればよい。放電レートは、放電期間における放電電流に対応するＣ値と当該蓄電池の基準Ｃ値との比で求めればよい。放電温度は、放電期間における電池温度を求めればよい。
サイクル履歴部２１は、電池放電時の放電量、放電レートおよび放電温度の履歴に基づき、現時点での積算放電量Ｐ、積算平均放電レートＣおよび積算平均放電温度Ｔｐを算出する。
積算放電量Ｐは、電池パック１０の使用開始から現時点までの電池放電時の放電量を積算して求めればよい。積算平均放電レートＣは、電池パック１０の使用開始から現時点までの電池放電時の放電レートの平均値とすればよい。積算平均放電温度Ｔｐは、電池パック１０の使用開始から現時点までの放電期間の温度の平均値とすればよい。
【００３６】
（放置履歴部）
放置履歴部２２は、電流センサ１１および温度センサ１２が検出した電池パック１０の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する。
放置時間は、電流センサ１１によって電流がゼロになってから放置終了（すなわち、充電又は放電への遷移）までの期間（放置期間）としてよい。放置期間の充電状態は、直近の満充電状態を基準にそれ以降であって当該放置期間に入るまでの積算放電量および積算充電量から求めればよい。放置温度は、放置期間における電池温度を求めればよい。
放置履歴部２２は、電池放置時の放置時間、充電状態および放置温度の履歴に基づき、現時点での積算放置時間Ｕ、積算平均充電状態Ｓおよび積算平均放置温度Ｔｓを算出する。
積算放置時間Ｕは、電池パック１０の使用開始から現時点までの電池放置時の放置時間を積算して求めればよい。積算平均充電状態Ｓは、電池パック１０の使用開始から現時点までの電池放置時の充電状態の平均値とすればよい。積算平均放置温度Ｔｓは、電池パック１０の使用開始から現時点までの放置期間の温度の平均値とすればよい。
【００３７】
（劣化算出部）
劣化算出部３は、電池使用履歴算出部２に記録された各種データに基づき、電池パック１０の現時点劣化量Ωを算出する。
【００３８】
ここで、電池放電時の電池パック１０の電極表面に付着する不純物量に相関する劣化量（以下、電池放電時劣化量Ωａという）は、積算放電量Ｐ、積算平均放電レートＣおよび積算平均放電温度Ｔｐをパラメータとする関係式により求めることができる。
【００３９】
図２は、本実施形態に係る積算放電量Ｐと電池放電時劣化量Ωａとの関係を示す図である。
図２に示すように、電池放電時劣化量Ωａは、積算放電量Ｐの１／２累乗に対し、以下の式（２）の関係を示す。
【数２】
2013181875.tif 000004
なお、Ωａは電池放電時劣化量、Ｐは積算放電量、Ｃは積算平均レート、Ｔｐは積算平均放電温度である。
また、ａ、ｂ、Ｌは電池材料に固有の定数であり、予め電池セル情報として劣化算出部３に記憶されている。
【００４０】
また、電池放置時の電池パック１０の電極表面に付着する不純物量に相関する劣化量（以下、電池放置時劣化量Ωｂという）は、積算放置時間Ｕ、積算平均充電状態Ｓおよび積算平均放置温度Ｔｓをパラメータとする関係式により求めることができる。
【００４１】
図３は、本実施形態に係る積算放置時間Ｕと電池放置時劣化量Ωｂとの関係を示す図である。
図３に示すように、電池放置時劣化量Ωｂは、積算放置時間Ｕの１／２累乗に対し、以下の式（３）の関係を示す。
【数３】
2013181875.tif 000005
なお、Ωｂは電池放置時劣化量、Ｕは積算放置時間、Ｓは積算平均充電状態、Ｔｓは積算平均放置温度である。
また、ｄ、ｅ、Ｍは電池材料に固有の定数であり、予め電池セル情報として劣化算出部３に記憶されている。
【００４２】
以上の原理を用い、劣化算出部３は、現時点劣化量Ωを、電池放電時劣化量Ωａと電池放置時劣化量Ωｂとの和として算出する。
すなわち、下記式（１）により、電池パック１０の現時点劣化量Ωを算出する。
【数４】
2013181875.tif 000006
なお、Ωは現時点劣化量、Ｐは積算放電量、Ｃは積算平均レート、Ｔｐは積算平均放電温度、Ｕは積算放置時間、Ｓは積算平均充電状態、Ｔｓは積算平均放置温度である。
【００４３】
劣化算出部３は、算出した現時点劣化量Ωと、予め設定した電池パック１０の寿命時劣化量Ωｍａｘとの比から現時点劣化率Ωｒａｔｅを算出する。
すなわち、下記式（４）により、電池パック１０の現時点劣化率Ωｒａｔｅを算出する。
【数５】
2013181875.tif 000007
ここで、Ωｒａｔｅは現時点劣化率、Ωは現時点劣化量、Ωｍａｘは寿命時劣化量である。
【００４４】
（劣化変遷予測部）
劣化変遷予測部４は、電池使用履歴算出部２に記録されたデータに基づき、電池パック１０の劣化量の予測関数Ωｆ（ｔ）を算出し、電池パック１０の劣化変遷を予測する。
【００４５】
まず、ユーザにより過去の任意のサンプリング期間を予め設定する。サンプリング期間は、例えば、現時点を含む１日や１週間や１月などの期間、前日のみ、一昨日のみ、先週、先月などの期間として指定できる。またサンプリング期間は、ユーザを変更した時点から現時点までの期間や、ユーザごとの期間などとして指定してもよい。
サンプリング期間は、サンプリング開始時点とサンプリング終了時点とを指定することで設定するものであってよい。
サンプリング期間は、劣化変遷予測部４に記憶しておく。
【００４６】
次に、サンプリング期間における、期間積算放電量Ｐ１、期間平均放電レートＣ１、期間平均放電温度Ｔｐ１、期間放置時間Ｕ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１を算出する。
期間積算放電量Ｐ１、期間平均放電レートＣ１および期間平均放電温度Ｔｐ１は、サイクル履歴部２１の電池使用時の放電量、放電レートおよび放電温度の履歴に基づく。
期間積算放電量Ｐ１は、サンプリング期間における放電量を積算して求めればよい。期間平均放電レートＣ１は、サンプリング期間における放電レートの平均値とすればよい。期間平均放電温度ＴＰ１は、サンプリング期間における放電温度の平均値とすればよい。
期間放置時間Ｕ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１は、放置履歴部２２の電池放置時の放置時間、充電状態および放置温度の履歴に基づく。
期間放置時間Ｕ１は、サンプリング期間における放置時間を積算して求めればよい。期間平均充電状態Ｓ１は、サンプリング期間における充電状態の平均値とすればよい。期間平均放置温度Ｔｓ１は、サンプリング期間における放置温度の平均値とすればよい。
【００４７】
次に、サンプリング期間における、期間平均放電レートＣ１、期間平均放電温度Ｔｐ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１を用い、式（１）におけるサンプリング期間における劣化量の係数Ｋ１および係数Ｒ１を算出する。
すなわち、下記式（５）、（６）のように算出される。
【数６】
2013181875.tif 000008
なお、Ｃ１は期間平均放電レート、Ｔｐ１は期間平均放電温度、Ｓ１は期間平均充電状態、Ｔｓ１は期間平均放置温度である。
【００４８】
次に、係数Ｋ１および係数Ｒ１を用いた予測関数Ωｆ（ｔ）を導出する。
すなわち、予測関数Ωｆ（ｔ）は、式（１）を用い、下記式（７）のように導出される。
【数７】
2013181875.tif 000009
ここで、Ωｆ（ｔ）は劣化量の予測関数、Ｐ（ｔ）は放電量の変数、Ｕ（ｔ）は放置時間の変数である。
【００４９】
図４は、本実施形態に係る電池パック１０の劣化量の予測関数Ωｆ（ｔ）を示す図である。図４では、横軸に使用時間をとり、縦軸に劣化量をとって予測関数Ωｆ（ｔ）を示す。図４に示す黒丸は、現時点に至るまでの記録された劣化量のプロットであり、白丸は、予測関数Ωｆ（ｔ）に従った将来予測される劣化量のプロットである。
図４に示すように、電池パック１０の劣化量の予測関数Ωｆ（ｔ）により、係数Ｋ１および係数Ｒ１を用いた場合、すなわちサンプリング期間の使用状態を維持した場合の劣化量の変遷が予測できる。
【００５０】
（寿命算出部）
寿命算出部５は、劣化算出部３が算出した現時点劣化量Ωおよび劣化変遷予測部４が導出した予測関数Ωｆ（ｔ）により、現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの期間を、現時点と、図４に示す予測関数Ωｆ（ｔ）が寿命時劣化量Ωｍａｘに到達する寿命時点との間の余命期間Ｔとして算出する。
電池パック１０の寿命を定める電池セル寿命などの情報は、寿命算出部５に予め記憶されている。
【００５１】
（表示部）
表示部６は、劣化算出部３で算出された電池パック１０の現時点劣化率Ωｒａｔｅや、寿命算出部５で算出された電池パック１０の現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔなどを表示する。
【００５２】
図５〜図６は、表示部６に表示されるデータの例である。
図５は、本実施形態に係る表示部６に表示される任意のサンプリング期間ＳＰ１に応じた余命期間Ｔ１を示す図である。図５に示す黒三角は、サンプリング期間ＳＰ１で記録された劣化量のプロットである。また、白三角は、サンプリング期間ＳＰ１の予測関数Ωｆ（ｔ）に従った将来予測される劣化量のプロットである。これらプロットは、電池パック１０の劣化量の変化を表す。
また、Ｔ１は、サンプリング期間ＳＰ１の予測関数Ωｆ（ｔ）に従った現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔである。
【００５３】
図６は、本実施形態に係る表示部６に表示される２種のサンプリング期間ＳＰ１，ＳＰ２に応じた余命期間Ｔ１，Ｔ２を示す図である。図６に示す黒三角は、サンプリング期間ＳＰ１で記録された劣化量のプロットであり、黒四角は、サンプリング期間ＳＰ２で記録された劣化量のプロットである。また、白三角は、サンプリング期間ＳＰ１の予測関数Ωｆ（ｔ）１に従った将来予測される劣化量のプロットであり、白四角は、サンプリング期間ＳＰ２の予測関数Ωｆ（ｔ）２に従った将来予測される劣化量のプロットである。これらプロットは、電池パック１０の劣化量の変化を表す。
また、Ｔ１は、サンプリング期間ＳＰ１の予測関数Ωｆ（ｔ）１に従った現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔであり、Ｔ２は、サンプリング期間ＳＰ２の予測関数Ωｆ（ｔ）２に従った現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔである。
図６に示すように、サンプリング期間（ＳＰ１，ＳＰ２）を異ならせて複数設定することで、複数の異なるサンプリング期間（ＳＰ１，ＳＰ２）における電池パック１０の使用状態のそれぞれに基づき、複数の余命期間Ｔ（Ｔ１，Ｔ２）として電池パック１０の寿命をそれぞれ正確に予測して同時に表示できる。
図６に示す表示により、ユーザは、どのサンプリング期間における電池パック１０の使用状態が好ましいか一目で視認できる。
【００５４】
以下、フローチャートを用いて電池パック１０の寿命予測方法について説明する。
図７は、本実施形態に係る電池パック１０の寿命予測方法を示すフローチャートである。
【００５５】
ステップＳ１では、サンプリング期間の設定を行う。ユーザが電池寿命予測システム１に対してサンプリング期間を入力し、劣化変遷予測部４にサンプリング期間が記憶される。サンプリング期間は、サンプリング開始時点および終了時点を指定することで入力されるものであってよい。
なお、サンプリング期間は、異なる期間（ＳＰ１，ＳＰ２など）で複数設定されるものであってよい。
【００５６】
ステップＳ２では、電池パック１０を稼動（放電）させた状態か否かを判別する。Ｓ２において、電池パック１０を稼動させた状態（放電）の場合には、Ｓ３へ移行する。Ｓ２において、電池パック１０を稼動させない状態（放置）の場合には、Ｓ６へ移行する。
【００５７】
ステップＳ３では、電流センサ１１と、温度センサ１２と、電圧センサ１３とによって、電池放電時の電流、電池温度および電圧を検出する。
【００５８】
ステップＳ４では、サイクル履歴部２１により、電流センサ１１および温度センサ１２が検出した電池放電時の電池パック１０の放電量、放電レートおよび放電温度を計測してこれらの履歴を記録する。
【００５９】
ステップＳ５では、サイクル履歴部２１により、ステップＳ４で求めた電池放電時の放電量、放電レートおよび放電温度の履歴に基づき、現時点での積算放電量Ｐ、積算平均放電レートＣおよび積算平均放電温度Ｔｐを算出する。
【００６０】
一方ステップＳ６では、電流センサ１１と、温度センサ１２と、電圧センサ１３とによって、電池放置時の電流、電池温度および電圧を検出する。
【００６１】
ステップＳ７では、放置履歴部２２により、電流センサ１１および温度センサ１２が検出した電池放置時の電池パック１０の放置時間、充電状態および放置温度を計測してこれらの履歴を記録する。
【００６２】
ステップＳ８では、放置履歴部２２により、電池放置時の放置時間、充電状態および放置温度の履歴に基づき、現時点での積算放置時間Ｕ、積算平均充電状態Ｓおよび積算平均放置温度Ｔｓを算出する。
【００６３】
ステップＳ９では、劣化算出部３により、サイクル履歴部２１と放置履歴部２２とから構成される電池使用履歴算出部２に記録された各種データに基づき、電池パック１０の現時点劣化量Ωを算出する。すなわち、式（１）により電池パック１０の現時点劣化量Ωを算出する。
【００６４】
ステップＳ１０では、劣化算出部３により、式（４）により電池パック１０の現時点劣化率Ωｒａｔｅを算出する。
【００６５】
ステップＳ１１では、劣化変遷予測部４により、電池使用履歴算出部２に記録されたデータに基づき、サンプリング期間の期間平均放電レートＣ１、期間平均放電温度Ｔｐ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１を算出する。
サンプリング期間が複数存在する場合には、それぞれの期間平均放電レートＣ１、期間平均放電温度Ｔｐ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１を算出する。
【００６６】
ステップＳ１２では、劣化変遷予測部４により、サンプリング期間の期間平均放電レートＣ１、期間平均放電温度Ｔｐ１、期間平均充電状態Ｓ１および期間平均放置温度Ｔｓ１を用い、電池パック１０の劣化量の予測関数Ωｆ（ｔ）を導出し、電池パック１０の劣化変遷を予測する。すなわち、式（７）により予測関数Ωｆ（ｔ）を導出し、図４に示すように電池パック１０の劣化変遷を予測する。
サンプリング期間が複数存在する場合には、それぞれの予測関数（Ωｆ（ｔ）１，Ωｆ（ｔ）２など）を導出する。
【００６７】
ステップＳ１３では、寿命算出部５により、図４に示すように、劣化変遷予測部４の予測に基づき、現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔを算出する。
サンプリング期間（ＳＰ１，ＳＰ２）が複数存在する場合には、それぞれの余命期間Ｔ（Ｔ１，Ｔ２）を算出する。
【００６８】
ステップＳ１４では、表示部６により、劣化算出部３で算出された電池パック１０の現時点劣化率Ωｒａｔｅや、図５〜図６に示す寿命算出部５で算出された電池パック１０の現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの余命期間Ｔなどを表示する。サンプリング期間（ＳＰ１，ＳＰ２）が複数存在する場合には、それぞれの余命期間Ｔ（Ｔ１，Ｔ２）を表示する。
【００６９】
以上の本実施形態に係る電池寿命予測システム１によれば、以下の効果を奏する。
（１）上記した積算放電量Ｐ、積算平均放電レートＣ、積算平均放電温度Ｔｐ、積算放置時間Ｕ、積算平均充電状態Ｓおよび積算平均放置温度Ｔｓという各種パラメータを用いて電池パック１０の内部抵抗値に基づかずに電池パック１０の現時点劣化率Ωｒａｔｅを正確に算出できる。このため、電池パック１０の劣化状態を、現時点の充電状態や電池温度によって変動しない真の状態として正確に把握できる。
【００７０】
（２）電池パック１０の現時点劣化量Ωを、電池使用時の電池パック１０の電極表面に付着する不純物量に相関する劣化量分（Ωａ）と、電池放置時の電池パック１０の電極表面に付着する不純物量に相関する劣化量分（Ωｂ）との和として、数式（１）に基づき定量的に算出できる。
【００７１】
（３）まず式（１）を用いて電池パック１０の現時点劣化量Ωを正確に算出し、次にサンプリング期間についての電池パック１０の劣化量の予測関数Ωｆ（ｔ）を導出する。そして導出した電池パック１０の劣化量の予測関数Ωｆ（ｔ）を用い、現時点劣化量Ωから寿命時劣化量Ωｍａｘに至るまでの正確な余命期間Ｔを算出する。
したがって、サンプリング期間における電池パック１０の使用状態に応じた電池パック１０の寿命を正確に予測できる。
【００７２】
（４）複数の異なるサンプリング期間（ＳＰ１，ＳＰ２）における電池パック１０の使用状態に基づき、それぞれの余命期間Ｔ（Ｔ１，Ｔ２）として電池パック１０の寿命を正確に予測できる。このため、採用する種々のサンプリング期間の電池パック１０の使用状態に応じて、電池パック１０の寿命を様々に予測できる。
また、それぞれの余命期間Ｔ（Ｔ１，Ｔ２）を同時に表示することにより、どのサンプリング期間における電池パック１０の使用状態が好ましいか一目で視認できる。
【００７３】
なお、本発明は上記実施形態に限定されず、本発明の目的を達成できる範囲での変形、改良などは本発明に包含される。
上記実施形態における電池寿命予測システム１では、電池パック１０の余命期間Ｔを求めることが前提のため、劣化算出部３で電池パック１０の現時点劣化量Ωを算出するものであったが、これに限られない。例えば、劣化算出部では電池パックの所定時点での劣化量を算出するものであってもよい。
また、本発明は、二次電池の寿命予測システムとしてだけでなく、二次電池の劣化率算出システムとしても構成できる。
【符号の説明】
【００７４】
１…電池寿命予測システム
２１…サイクル履歴部（放電履歴記録部）
２２…放置履歴部（放置履歴記録部）
３…劣化算出部
４…劣化変遷予測部
５…寿命算出部
１０…電池パック
【図１】
2013181875.tif 000011
【図２】
2013181875.tif 000012
【図３】
2013181875.tif 000013
【図４】
2013181875.tif 000014
【図５】
2013181875.tif 000015
【図６】
2013181875.tif 000016
【図７】
2013181875.tif 000017