(57)【要約】  本発明の実施の形態に係るエネルギー貯蔵システムの充電容量算出する方法は、前記エネルギー貯蔵システムに一定の値を有する所定の電流を印加する電流印加ステップと、前記エネルギー貯蔵システムに一定の値を有する所定の電流が印加される間の前記エネルギー貯蔵システムの電圧を測定する第１の電圧測定ステップと、前記エネルギー貯蔵システムに印加されていた、一定の値を有する所定の電流を切った後、エネルギー貯蔵システムの電圧を測定する第２の電圧測定ステップと、前記第１の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、第２の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、及び前記一定の値を有する所定の電流に基づいて、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する抵抗（Ｒ）及びキャパシタンス（Ｃ）算出ステップと、を含んでなる。
2020530905.tif 000002
【特許請求の範囲】
【請求項１】
エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する方法であって、
前記エネルギー貯蔵システムに一定の値を有する所定の電流を印加する電流印加ステップと、
前記エネルギー貯蔵システムに一定の値を有する所定の電流が印加される間の前記エネルギー貯蔵システムの電圧を測定する第１の電圧測定ステップと、
前記エネルギー貯蔵システムに印加されていた、一定の値を有する所定の電流を切った後、エネルギー貯蔵システムの電圧を測定する第２の電圧測定ステップと、
前記第１の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、第２の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、及び前記一定の値を有する所定の電流に基づいて、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する抵抗（Ｒ）及びキャパシタンス（Ｃ）算出ステップと、
を含んでなることを特徴とするエネルギー貯蔵システムの充電容量算出方法。
【請求項２】
前記電流印加ステップにおいて印加される一定の値を有する所定の電流は、前記エネルギー貯蔵システムの抵抗及びキャパシタンス値が測定可能な最小の電流値から２ｃ−ｒａｔｅまでの電流であることを特徴とする請求項１に記載のエネルギー貯蔵システムの充電容量算出方法。
【請求項３】
前記電流印加ステップにおいて印加される一定の値を有する所定の電流は、前記エネルギー貯蔵システムの抵抗及びキャパシタンス値が測定可能な最小の電流値から０．１ｃ−ｒａｔｅまでの電流であることを特徴とする請求項１又は２に記載のエネルギー貯蔵システムの充電容量算出方法。
【請求項４】
前記エネルギー貯蔵システムの抵抗及びキャパシタンスを用いて、拡張カルマンフィルターの状態変数を算出する状態変数算出ステップと、
前記状態変数を前記拡張カルマンフィルターに代入して、前記エネルギー貯蔵システムの充電容量を算出するエネルギー貯蔵システム充電容量算出ステップと、
をさらに含んでなることを特徴とする請求項１〜３のいずれか一項に記載のエネルギー貯蔵システムの充電容量算出方法。
【請求項５】
前記拡張カルマンフィルターは、
前記状態変数を時間更新するステップと、
前記状態変数の誤差共分散を時間更新するステップと、
前記拡張カルマンフィルターのカルマンゲインを算出するステップと、
前記カルマンゲインを用いて、状態変数を推定するステップと、
前記カルマンゲインを用いて、前記状態変数の誤差共分散を補正するステップと、
を繰り返し行うことを特徴とする請求項４に記載のエネルギー貯蔵システムの充電容量算出方法。
【請求項６】
前記状態変数算出ステップは、
前記エネルギー貯蔵システムの充電容量に対する状態変数を生成する充電容量状態変数算出ステップと、
前記エネルギー貯蔵システムの電圧に対する状態変数を生成する電圧状態変数算出ステップと、
を含んでなることを特徴とする請求項４又は５に記載のエネルギー貯蔵システムの充電容量算出方法。
【請求項７】
前記充電容量状態変数算出ステップは、
前記エネルギー貯蔵システムにおいて測定された電流を積算する電流積算法を用いて、充電容量状態変数を算出し、
前記算出された充電容量状態変数は、下記の数式１に従って時間更新されることを特徴とする請求項６に記載のエネルギー貯蔵システムの充電容量算出方法。
【数１】
2020530905.tif 000014
（但し、Ｑｃａｐａｃｉｔｙ：二次電池の容量、ｋ：時間インデックス、i［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【請求項８】
前記電圧状態変数算出ステップは、
電圧モデル回路を用いて、前記エネルギー貯蔵システムの電圧状態変数を算出し、
前記算出された電圧状態変数は、下記の数式２に従って時間更新されることを特徴とする請求項６又は７に記載のエネルギー貯蔵システムの充電容量算出方法。
【数２】
2020530905.tif 000015
（但し、Ｒ１及びＣ１は、図５の回路モデルに含まれている抵抗及びキャパシタンスの値、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【請求項９】
前記電圧モデル回路は、
前記エネルギー貯蔵システムに流れる電流値に応じて所定の電圧を出力する開放電圧源と、
一つ以上の抵抗と、
一つ以上のキャパシターと、
を備えてなることを特徴とする請求項８に記載のエネルギー貯蔵システムの充電容量算出方法。
【請求項１０】
複数のバッテリーラックと、
前記複数のバッテリーラックを制御するバッテリーセクションコントローラー（ＢＳＣ）と、
を備えてなるエネルギー貯蔵システムにおいて、
前記ＢＳＣは、
エネルギー貯蔵システムの出力電圧を測定する電圧測定部と、
エネルギー貯蔵システムの出力電流を測定する電流測定部と、
一定の値を有する所定の電流により算出された温度別のエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）のルックアップテーブルが保存されている保存部と、
前記ルックアップテーブルに基づいて、前記エネルギー貯蔵システムの充電容量を算出する充電容量算出部と、
を備えてなることを特徴とするエネルギー貯蔵システム。
【請求項１１】
前記一定の値を有する所定の電流は、
前記エネルギー貯蔵システムの抵抗及びキャパシタンス値が測定可能な最小の電流値から２ｃ−ｒａｔｅまでの電流であることを特徴とする請求項１０に記載のエネルギー貯蔵システム。
【請求項１２】
前記一定の値を有する所定の電流は、
前記エネルギー貯蔵システムの抵抗及びキャパシタンス値が測定可能な最小の電流値から０．１ｃ−ｒａｔｅまでの電流であることを特徴とする請求項１０又は１１に記載のエネルギー貯蔵システム。
【請求項１３】
前記電圧測定部は、
前記エネルギー貯蔵システムに所定の電流が印加される間の前記エネルギー貯蔵システムの電圧である第１の電圧と、
前記エネルギー貯蔵システムに印加される電流が切れた後の前記エネルギー貯蔵システムの電圧である第２の電圧と、
を測定することを特徴とする請求項１０〜１２のいずれか一項に記載のエネルギー貯蔵システム。
【請求項１４】
前記充電容量算出部は、
前記ルックアップテーブルの抵抗及びキャパシタンスに基づいて、エネルギー貯蔵システムの電圧に対する状態変数を算出する電圧状態変数算出部と、
前記測定された出力電流を積算する電流積算法を用いて、エネルギー貯蔵システムの充電容量に対する状態変数を算出する充電容量状態変数算出部と、
を備えてなり、
前記算出された電圧状態変数及び充電容量状態変数をカルマンフィルターに代入して、エネルギー貯蔵システムの充電容量を算出することを特徴とする請求項１０〜１３のいずれか一項に記載のエネルギー貯蔵システム。
【請求項１５】
前記充電容量状態変数算出部において算出された充電容量状態変数は、
下記の数式１に従って時間更新されることを特徴とする請求項１４に記載のエネルギー貯蔵システム。
【数３】
2020530905.tif 000016
（但し、Ｑｃａｐａｃｉｔｙ：二次電池の容量、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【請求項１６】
前記電圧状態変数算出部において算出される前記算出された電圧状態変数は、
下記の数式２に従って時間更新されることを特徴とする請求項１４又は１５に記載のエネルギー貯蔵システム。
【数４】
2020530905.tif 000017
（但し、Ｒ１及びＣ１は、図５の回路モデルに含まれている抵抗及びキャパシタンスの値、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エネルギー貯蔵システムの充電容量を算出する装置及び方法に関する。
【０００２】
具体的には、拡張カルマンフィルターに用いられる電圧モデルの正確度を高めて、エネルギー貯蔵システムの充電容量を正確に算出する装置及び方法に関する。
【背景技術】
【０００３】
産業の発達と伴い、電力の需要が高まりつつあり、昼夜間、季節間、日別間の電力使用量の格差が次第に広がってきている。近年、この理由から、系統の余剰の電力を活用して、ピーク負荷を削減するための数多くの技術が急速に開発されている。
【０００４】
これらの技術に代表されるものが、系統の余剰の電力をバッテリーに貯蔵したり、系統の不足した電力をバッテリーから供給したりするエネルギー貯蔵システムである。
【０００５】
かようなエネルギー貯蔵システムの寿命を長時間にわたって維持し、かつ、安全に使用するためには、エネルギー貯蔵システムを適正な充電容量の範囲内で駆動せねばならず、エネルギー貯蔵システムの寿命は、充電／放電回数に応じて大きく異なる。
【０００６】
したがって、エネルギー貯蔵システムの電流、電圧及び温度を測定して、エネルギー貯蔵システムの充電容量（ＳＯＣ）を正確に算出することにより、エネルギー貯蔵システムの状態を把握することが重要である。
【０００７】
従来のエネルギー貯蔵システムの充電容量（ＳＯＣ）を算出する方法としては、電流積算法があった。しかしながら、従来の電流積算法は、電流を測定する電流センサーに誤差が存在し、電流を積算する過程において電流センサーの誤差もまた積算されてしまい、正確な充電容量を算出することが困難であった。
【０００８】
このような従来の問題を解消するための別の従来の技術としては、エネルギー貯蔵システムの充電容量を算出するとき、前記電流センサーにおいて測定される電流に基づいて、充電容量を算出する電流積算方式と、エネルギー貯蔵システムの電圧モデルを拡張カルマンフィルターに代入して充電容量を算出する方法とがあった。
【０００９】
しかしながら、従来の拡張カルマンフィルターに電流積算方式とエネルギー貯蔵システムの電圧モデルを使用する技術においても、エネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）、キャパシタンス（Ｃ）の正確度が低下して、最終的に算出される充電容量に誤差が生じてしまうという問題が依然として存在していた。
【００１０】
したがって、本発明では、エネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）、キャパシタンス（Ｃ）を正確に算出し、電圧モデルの比重が高いカルマンフィルターを用いて、最終的に算出されるエネルギー貯蔵システムの充電容量の誤差を減らす装置及び方法を提案する。
【先行技術文献】
【特許文献】
【００１１】
【特許文献１】大韓民国特許公報第１０−２０１３−０１０５１２３号
【発明の概要】
【発明が解決しようとする課題】
【００１２】
本発明は、従来の問題を解消するためのものであって、エネルギー貯蔵システムの電圧モデルを正確に算出して、エネルギー貯蔵システムの正確な充電容量を算出する装置及び方法を提供する。
【００１３】
より具体的には、エネルギー貯蔵システムの電圧モデルに用いられる抵抗（Ｒ）、キャパシタンス（Ｃ）を正確に算出して、正確な電圧モデルを拡張カルマンフィルターに代入することにより、正確な充電容量を算出する装置及び方法を提供する。
【課題を解決するための手段】
【００１４】
本発明の実施の形態に係るエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する方法は、前記エネルギー貯蔵システムに一定の値を有する所定の電流を印加する電流印加ステップと、前記エネルギー貯蔵システムに一定の値を有する所定の電流が印加される間の前記エネルギー貯蔵システムの電圧を測定する第１の電圧測定ステップと、前記エネルギー貯蔵システムに印加されていた、一定の値を有する所定の電流を切った後、エネルギー貯蔵システムの電圧を測定する第２の電圧測定ステップと、前記第１の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、第２の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、及び前記一定の値を有する所定の電流に基づいて、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する抵抗（Ｒ）及びキャパシタンス（Ｃ）算出ステップと、を含んでなってもよい。
【００１５】
前記電流印加ステップにおいて印加される一定の値を有する所定の電流は、ｃ−ｒａｔｅ（Ｃレート）の低い電流である。
【００１６】
本発明の他の実施の形態に係る低いｃ−ｒａｔｅにより算出された抵抗（Ｒ）及びキャパシタンス（Ｃ）に基づいて算出されるエネルギー貯蔵システムの電圧モデル及び拡張カルマンフィルターを用いてエネルギー貯蔵システムの充電容量を算出する方法は、前記拡張カルマンフィルターの状態変数を算出する状態変数算出ステップ、前記状態変数を前記拡張カルマンフィルターに代入してエネルギー貯蔵システムの充電容量を算出するエネルギー貯蔵システム充電容量算出ステップを含んでなってもよい。
【００１７】
前記拡張カルマンフィルターは、前記状態変数を時間更新するステップと、前記状態変数の誤差共分散を時間更新するステップと、前記拡張カルマンフィルターのカルマンゲイン（ｇａｉｎ）を算出するステップと、前記カルマンゲイン（ｇａｉｎ）を用いて、状態変数を推定するステップと、前記カルマンゲイン（ｇａｉｎ）を用いて、前記状態変数の誤差共分散を補正するステップと、を繰り返し行ってもよい。
【００１８】
前記状態変数算出ステップは、前記エネルギー貯蔵システムの充電容量に対する状態変数を生成する充電容量状態変数算出ステップ及び前記エネルギー貯蔵システムの電圧に対する状態変数を生成する電圧状態変数算出ステップを含んでなってもよい。
【００１９】
前記充電容量状態変数算出ステップにおいては、前記エネルギー貯蔵システムの電流を積算して充電容量状態変数を算出し、算出された充電容量状態変数は、拡張カルマンフィルターにおいて数式１に基づいて前記エネルギー貯蔵システムの電流を積算して算出された充電容量状態変数を時間更新してもよい。
【００２０】
【数１】
2020530905.tif 000003
（但し、Ｑｃａｐａｃｉｔｙ：二次電池の容量、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【００２１】
前記電圧状態変数算出ステップにおいては、電圧モデル回路を用いて、前記エネルギー貯蔵システムの電圧状態変数を算出し、前記算出された電圧状態変数は、拡張カルマンフィルターにおいて、数式２に基づいて、前記エネルギー貯蔵システムの電圧状態変数を時間更新してもよい。
【００２２】
【数２】
2020530905.tif 000004
（但し、Ｒ１及びＣ１は、図５の回路モデルに含まれている抵抗及びキャパシタンスの値、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【００２３】
一方、本発明の実施の形態に係るエネルギー貯蔵システムは、複数のバッテリーラックと、前記複数のバッテリーラックを制御するバッテリーセクションコントローラー（ＢＳＣ；Ｂａｔｔｅｒｙ  Ｓｅｃｔｉｏｎ  Ｃｏｎｔｒｏｌｌｅｒ）と、を備えてなり、
前記ＢＳＣは、エネルギー貯蔵システムの出力電圧を測定する電圧測定部と、エネルギー貯蔵システムの出力電流を測定する電流測定部と、０．１ｃ−ｒａｔｅの電流により算出された温度別のエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）のルックアップテーブルが保存されている保存部と、前記測定された電圧と電流及び算出された抵抗及びキャパシタンスに基づいて、エネルギー貯蔵システムの充電容量を算出する充電容量算出部と、を備えてなってもよい。
【００２４】
前記電圧測定部は、前記エネルギー貯蔵システムに所定の電流が印加される間のエネルギー貯蔵システムの電圧である第１の電圧及び前記エネルギー貯蔵システムに印加される電流が切れた後に、エネルギー貯蔵システムの電圧である第２の電圧を測定してもよい。
【００２５】
前記充電容量算出部は、前記ルックアップテーブルの抵抗及びキャパシタンスに基づいて、エネルギー貯蔵システムの電圧に対する状態変数を算出する電圧状態変数算出部、及び前記測定されたエネルギー貯蔵システムの出力電流を積算する電流積算法を用いて、エネルギー貯蔵システムの充電容量に対する状態変数を算出する充電容量状態変数算出部を備えてなり、前記算出された電圧状態変数及び充電容量状態変数をカルマンフィルターに代入して、エネルギー貯蔵システムの充電容量を算出してもよい。
【発明の効果】
【００２６】
本発明は、エネルギー貯蔵システムの電圧モデルに用いられる抵抗（Ｒ）、キャパシタンス（Ｃ）を正確に算出してエネルギー貯蔵システムの電圧モデルを正確に算出することができる。
【００２７】
また、本発明は、正確なエネルギー貯蔵システムの電圧モデルを算出することにより、エネルギー貯蔵システムの充電容量を正確に算出することができる。
【図面の簡単な説明】
【００２８】
【図１】本発明のエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）、キャパシタンス（Ｃ）を算出する方法を示すフローチャートである。
【図２】本発明の実施の形態に係る拡張カルマンフィルターの動作を示すフローチャートである。
【図３】本発明の実施の形態に係るエネルギー貯蔵システムの電圧モデル算出するのに用いられる抵抗（Ｒ）、キャパシタンス（Ｃ）を算出する実際の実験方法（ＤＣＰＲ）を示す図である。
【図４】従来の電圧モデルに基づいて算出された放電後の充電容量の変化と、本発明の電圧モデルの算出方法により算出された充電容量の変化とを比較した図である。
【図５】本発明の実施の形態に係るエネルギー貯蔵システムの電圧を算出する電圧モデル回路である。
【図６】本発明のエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗及びキャパシタンス算出方法を用いてエネルギー貯蔵システムの充電容量を算出する方法を用いて算出された充電容量と、従来の技術（電流積算方法のみを用いる）により算出された充電容量とを比較する図である。
【図７】誤差条件を反映して、本発明のエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗及びキャパシタンスの算出方法を用いてエネルギー貯蔵システムの充電容量を算出する方法を用いて算出された充電容量と、既存のアルゴリズム（高いｃ−ｒａｔｅにより算出された電圧モデル及びカルマンフィルターにおいて電圧モデルの比重が低いアルゴリズムを用いる従来の技術）により算出された充電容量とを比較する図である。
【図８】既存のアルゴリズムと、本発明の充電容量推定方法を用いて容量テストを行った結果とを比較した図である。
【図９】本発明の実施の形態に係るエネルギー貯蔵システムの構成を示す図である。
【発明を実施するための形態】
【００２９】
以下では、添付図面に基づいて、本発明が属する技術分野において通常の知識を有する者が容易に実施できるように本発明の実施の形態について詳しく説明する。しかしながら、本発明は、種々の異なる形態に具体化可能であり、ここで説明する実施の形態に何ら限定されるものではない。なお、図中、本発明を明確に説明するために、説明とは無関係な部分は省略し、明細書の全般に亘って、類似の部分には類似の図面符号を付している。
【００３０】
「第１の」、「第２の」などのように序数を含む言い回しは、様々な構成要素を説明するうえで使用可能であるが、前記構成要素は、前記言い回しによって何等限定されない。前記言い回しは、ある構成要素を他の構成要素から区別する目的でしか使えない。例えば、本発明の権利範囲を逸脱しない範囲内において第１の構成要素は第２の構成要素と命名されてもよく、同様に、第２の構成要素もまた第１の構成要素と命名されてもよい。本出願において用いた用語は、単に特定の実施の形態を説明するために用いられたものであり、本発明を限定しようとする意図はない。単数の表現は、文脈からみて明らかに他の意味を有さない限り、複数の言い回しを含む。
【００３１】
明細書の全般に亘って、ある部分が他の部分と「連結」されているとか、「接続」されているとか、と言及された場合、これは、前記ある部分が前記他の部分に「直接的に連結されたり接続」されたりする場合だけではなく、これらの間に他の素子を間に挟んで「電気的に連結されたり接続」されたりする場合をも含む。なお、ある部分がある構成要素を「備える」としたとき、これは、特に断りのない限り、他の構成要素を除外するわけではなく、他の構成要素をさらに備えていてもよいことを意味する。本願の明細書の全般に亘って用いられる度合いの言い回しである「〜するステップ」又は「〜のステップ」は、「〜のためのステップ」を意味するものではない。
【００３２】
本発明において用いられる用語としては、本発明における機能を考慮しつつ、できる限り現在汎広く用いられている一般的な用語を選択したが、これは、当分野に携わっている技術者の意図又は判例、新たな技術の出現などによって異なる。なお、特定の場合は、出願人が任意に選定した用語もあり、この場合に、当該する発明の説明の部分の欄において詳しくその意味を記載する。よって、本発明において用いられる用語は、単なる用語の名称ではなく、その用語が有する意味と本発明の全般に亘っての内容を踏まえて定義されるべきである。
【００３３】
１．本発明の実施の形態に係るエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する方法
従来のエネルギー貯蔵システムの充電容量を算出する方法であって、電流積算方法を用いてエネルギー貯蔵システムの充電容量を算出していた。電流積算方法は、充電電流及び放電電流を時間により積分して充電容量を決定する。しかしながら、電流積算方法は、充電電流及び放電電流を測定する電流センサーに誤差が存在し、このような誤差は、電流を積算する時間が累積されることにつれて一緒に累積されて、最終的に電流積算方法により算出される充電容量の正確度が低下してしまう。
【００３４】
従来には、このような問題を解消するために、エネルギー貯蔵システムの電圧モデルと電流積算方法の両方を拡張カルマンフィルターに代入して、より正確な充電容量を算出していた。
【００３５】
しかし、従来には、拡張カルマンフィルターに代入されるエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）が正確に算出されないため、エネルギー貯蔵システムの充電容量に誤差が依然として存在していた。
【００３６】
したがって、本発明では、エネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）を正確に算出する方法を提案する。
【００３７】
図１は、本発明のエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する方法である。
【００３８】
本発明は、エネルギー貯蔵システムの電圧モデルを算出するのにエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）が用いられる。
【００３９】
本発明のエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）の算出方法は、低いｃ−ｒａｔｅの電流を使用するところに特徴がある。
【００４０】
本発明において、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）の算出方法は、図３に示すＤＣＰＲ方法を用いて算出可能である。
【００４１】
以下では、図１及び図３に基づいて、本発明のエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する方法について説明する。
【００４２】
１−１．エネルギー貯蔵システムに一定の値を有する所定の電流を印加する電流印加ステップ（Ｓ１１０）
本発明は、エネルギー貯蔵システムに一定の値を有する所定の電流を印加して、エネルギー貯蔵システムの電圧の変化を測定し、これに基づいて、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）値を算出する。
【００４３】
前記電流印加ステップは、所定の充電容量で充電されている状態のエネルギー貯蔵システムに一定の値を有する所定の電流を印加するステップである。このように、一定の値を有する所定の電流をエネルギー貯蔵システムに印加した後、次の第１の電圧測定ステップを行う。このとき、印加される一定の値を有する所定の電流は、低いｃ−ｒａｔｅを有する電流である。ｃ−ｒａｔｅは、電流の強さを意味してもよい。
【００４４】
印加される一定の値を有する電流としては、通算して２〜３ｃ−ｒａｔｅ以上を有する電流を用いる。ところが、本発明の発明者は、２〜３ｃ−ｒａｔｅよりも低いｃ−ｒａｔｅの電流を印加して、エネルギー貯蔵システムの抵抗及びキャパシタンス値をより正確に算出することを想到した。
【００４５】
したがって、本発明における前記低いｃ−ｒａｔｅの電流は、エネルギー貯蔵システムの抵抗及びキャパシタンスを測定可能な最小の電流値から２ｃ−ｒａｔｅまでの値を有する電流であってもよい。
【００４６】
好ましくは、前記低いｃ−ｒａｔｅを有する電流は、０．１ｃ−ｒａｔｅの電流であってもよく、これは、１ｃ−ｒａｔｅの電流の１／１０に相当する値である。
【００４７】
高いｃ−ｒａｔｅを有する電流を印加すれば、電流の印加が終わった後、エネルギー貯蔵システムの電圧に変動が生じる。これに対し、低いｃ−ｒａｔｅを有する電流を印加すれば、電流の印加が終わった後、エネルギー貯蔵システムの電圧が素早く安定してより正確なエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出することができる。
【００４８】
１−２．エネルギー貯蔵システムに一定の値を有する所定の電流が印加される間のエネルギー貯蔵システムの電圧を測定する第１の電圧測定ステップ（Ｓ１２０）
本発明の第１の電圧測定ステップは、エネルギー貯蔵システムに一定の値を有する所定の電流が流れる間にエネルギー貯蔵システムの電圧を測定するステップである。
【００４９】
具体的に、前記電圧印加ステップにおいて所定の充電容量で充電されている状態のエネルギー貯蔵システムに一定の値を有する所定の電流を印加すれば、前記エネルギー貯蔵システムの電圧は、所定の値だけ変動する。このとき、測定されるエネルギー貯蔵システムの電圧は、後述する抵抗（Ｒ）及びキャパシタンス（Ｃ）算出ステップにおいて用いられる。
【００５０】
１−３．エネルギー貯蔵システムに印加されていた一定の値を有する所定の電流を切った後、エネルギー貯蔵システムの電圧を測定する第２の電圧測定ステップ（Ｓ１３０）
本発明の第２の電圧測定ステップは、前記第１の電圧測定ステップが終わった後、一定の値を有する所定の電流の印加を終了し、エネルギー貯蔵システムの電圧を測定するステップである。
【００５１】
具体的に、エネルギー貯蔵システムに一定の値を有する所定の電流が印加されてから切れると、エネルギー貯蔵システムの電圧は変動する。例えば、充電方向に電流が印加されてから電流が切れる場合には、エネルギー貯蔵システムの電圧が下がり、放電方向に電流が印加されてから電流が切れる場合には、エネルギー貯蔵システムの電圧が上がってもよい。このときに測定されるエネルギー貯蔵システムの電圧の変動は、後述する抵抗（Ｒ）及びキャパシタンス（Ｃ）算出ステップにおいて用いられる。
【００５２】
１−４．前記第１の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧、第２の電圧測定ステップにおいて測定されたエネルギー貯蔵システムの電圧及び前記一定の値を有する所定の電流に基づいて、エネルギー貯蔵システムの抵抗（Ｒ）、キャパシタンス（Ｃ）を算出する抵抗（Ｒ）、キャパシタンス（Ｃ）算出ステップ（Ｓ１４０）
本発明の抵抗（Ｒ）、キャパシタンス（Ｃ）算出ステップは、前記第１の電圧測定ステップ（Ｓ１２０）において測定されたエネルギー貯蔵システムの電圧、第２の電圧測定ステップ（Ｓ１３０）において測定されたエネルギー貯蔵システムの電圧及び前記一定の値を有する所定の電流に基づいて、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出するステップである。具体的に、エネルギー貯蔵システムに印加される一定の値を有する所定の電流は既知であり、前記一定の値を有する所定の電流の流れの有無による電圧をも既知であるため、エネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出することができる。
【００５３】
例えば、図５の回路において、前記Ｖｃｅｌｌは、ＯＣＶ電圧とＲ０に印加される電圧及びＣ１に印加される電圧を合計した値である。一方、ＯＣＶ値は、バッテリーセルの充電容量に応じて既に設定されている値であり、Ｒ０の抵抗値は、時間ｔに測定されるＶｃｅｌｌの電圧（Ｖｔ）と時間ｔ＋１において測定されるＶｃｅｌｌの電圧（Ｖｔ＋１）との差分を回路に流れる低いｃ−ｒａｔｅの電流で割った値であり、Ｒ０に印加される電圧もまた既知である。したがって、Ｃ１に印加される電圧値及び電流値が既知であるため 、Ｒ１値とＣ１値を算出することができる。
【００５４】
このような方法でエネルギー貯蔵システムの温度を変更しながら繰り返し算出されたエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を温度別のエネルギー貯蔵システムの抵抗及びキャパシタンスルックアップテーブルの形式で作成して保存されてもよい。
【００５５】
一方、本発明では、低いｃ−ｒａｔｅ（０．１ｃ−ｒａｔｅ）を有する電流を用いているため、前記第２の電圧測定ステップにおいて測定されるエネルギー貯蔵システムの電圧は電流が流れない後に素早く安定することから、電流が流れるときのエネルギー貯蔵システムの電圧と、電流が流れないときのエネルギー貯蔵システムの電圧とが正確に測定される。
【００５６】
したがって、前述したエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）の算出方法を用いて、エネルギー貯蔵システムの電圧モデルを算出し、これに基づいて、エネルギー貯蔵システムの充電容量を算出すれば、エネルギー貯蔵システムの充電容量を正確に算出することができる。
【００５７】
一方、本発明の電流印加ステップと第１の電圧測定ステップとの間、及び第１の電圧測定ステップと第２の電圧測定ステップとの間には所定の休止（ｒｅｓｔ）時間があってもよい。
【００５８】
以下では、前述した本発明のエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）の算出方法を用いて、エネルギー貯蔵システムの電圧モデルを正確に算出し、これに基づいて、エネルギー貯蔵システムの充電容量を算出する方法について説明する。
【００５９】
２．本発明のエネルギー貯蔵システムの充電容量を算出する方法。
従来のエネルギー貯蔵システムの充電容量は、電圧モデルと電流積算モデルを拡張カルマンフィルターの状態変数として代入して算出していた。しかしながら、従来には、電圧モデルを算出するときに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）値が正確に算出されないため、算出される充電容量に誤差が生じていた。
【００６０】
具体的に、従来には、電圧モデルを算出するときに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出するとき、約２〜３ｃ−ｒａｔｅの高いｃ−ｒａｔｅの電流を用いて算出していた。
【００６１】
しかし、本発明では、電圧モデルを算出するときに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出するとき、低いｃ−ｒａｔｅの電流を用いることにより、抵抗（Ｒ）及びキャパシタンス（Ｃ）の算出正確度を高めている。
【００６２】
一方、前記低いｃ−ｒａｔｅの電流は、エネルギー貯蔵システムの抵抗及びキャパシタンスを測定可能な最小の電流値から２ｃ−ｒａｔｅまでの値であってもよい。
【００６３】
好ましくは、前記低いｃ−ｒａｔｅを有する電流は、０．１ｃ−ｒａｔｅの電流であってもよく、これは、１ｃ−ｒａｔｅの電流の１／１０に相当する値である。図４は、従来の電圧モデルに基づいて算出された放電後の充電容量の変化と、本発明の電圧モデルの算出方法により算出された充電容量の変化とを比較した図である。
【００６４】
図４を参照すると、従来の方式により算出された抵抗（Ｒ）及びキャパシタンス（Ｃ）に基づいて電圧モデルを算出すれば、放電後に拡張カルマンフィルターにより算出される充電容量の変化率（３％）が大きいことを確認することができる。
【００６５】
放電後には充電容量の変化がないことが正常であるが、従来の方式は３％の変化率があるため、従来の方式は不正確である。
【００６６】
しかし、前述した本発明のように、低いｃ−ｒａｔｅの電流を用いて算出された抵抗（Ｒ）及びキャパシタンス（Ｃ）に基づいて電圧モデルを算出すれば、放電後に拡張カルマンフィルターにより算出される充電容量の変化率（０．５％）が非常に小さいことを確認することができる。
【００６７】
すなわち、本発明は、放電後における充電容量の変化率が小さいことから、より正確に充電容量を算出する効果を奏するということが分かる。
【００６８】
したがって、本発明のエネルギー貯蔵システムの充電容量算出方法を用いれば、従来よりも正確な充電容量を算出することができる。
【００６９】
図２は、本発明の実施の形態に係る拡張カルマンフィルターの動作を示すフローチャートである。
【００７０】
一方、拡張カルマンフィルターは、動的なシステムに対して外部において測定可能な変数とシステムの外乱（外部誤差）を考慮してシステムの状態を確率統計的に推定可能な適応的なソフトウェアアルゴリズムである。
【００７１】
具体的に、拡張カルマンフィルターは、状態変数を時間更新するステップと、前記状態変数の誤差共分散を時間更新するステップと、前記拡張カルマンフィルターのカルマンゲイン（ｇａｉｎ）を算出するステップと、前記算出されたカルマンゲイン（ｇａｉｎ）を用いて状態変数を推定するステップと、前記カルマンゲイン（ｇａｉｎ）を用いて前記状態変数の誤差共分散を補正するステップと、を繰り返し行ってもよい。
【００７２】
すなわち、拡張カルマンフィルターは、推定→補正→推定→補正を繰り返し行って誤差を減らす方法である。
【００７３】
本発明において、拡張カルマンフィルターの状態方程式は、状態変数を含んでなり、状態変数を時間に応じて更新する。状態変数としては、エネルギー貯蔵システムの充電容量状態変数とエネルギー貯蔵システムの電圧状態変数が用いられる。
【００７４】
２．１  状態変数の一つであるエネルギー貯蔵システムの充電容量状態変数を算出する方法
本発明の状態変数のうち、エネルギー貯蔵システムの充電容量算出方法は、エネルギー貯蔵システムの電流を測定し、測定された電流を積算することにより、エネルギー貯蔵システムの充電容量を算出することができる。
【００７５】
エネルギー貯蔵システムの電流を測定する方法としては、エネルギー貯蔵システムを制御及び管理するエネルギー貯蔵システム制御部において測定されるエネルギー貯蔵システムの電流情報を用いる方法が挙げられる。
【００７６】
すなわち、充電容量状態変数の初期値は、エネルギー貯蔵システムの電流を測定し、これを電流積算法により積算することにより、ＳＯＣ１値が算出可能である。
【００７７】
前述した方法により算出されるエネルギー貯蔵システムの充電容量状態変数は、拡張カルマンフィルターにおいて下記の数式１で時間更新されてもよい。
【００７８】
【数３】
2020530905.tif 000005
（但し、Ｑｃａｐａｃｉｔｙ：二次電池の容量、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【００７９】
２．２  状態変数のうち、エネルギー貯蔵システムの電圧状態変数を算出する方法
本発明の状態変数のうち、エネルギー貯蔵システムの電圧状態変数を算出する方法は、図５の電圧モデル回路を用いて算出可能である。
【００８０】
図５の電圧モデル回路は、エネルギー貯蔵システムに流れる所定の電流値に応じて可変となる開放電圧源２１０と、少なくとも一つの抵抗（Ｒ）０  ２２０と、Ｒ１  ２３０と、少なくとも一つのキャパシターＣ１と、を備えてなる。前記電圧モデル回路を構成する抵抗（Ｒ）０  ２２０と、Ｒ１  ２３０とキャパシターＣ１との値は、上述した本発明の実施の形態に従って低いｃ−ｒａｔｅの一定の値を有する電流により算出されて保存される温度別のエネルギー貯蔵システムの抵抗及びキャパシタンスルックアップテーブルから検出してもよい。
【００８１】
このように、温度別のエネルギー貯蔵システムの抵抗及びキャパシタンスルックアップテーブルから前記電圧モデル回路を構成する抵抗及びキャパシタンス値を検出すれば、電圧状態変数の初期値であるＶ１値が算出可能である。
【００８２】
具体的に、Ｒ１、Ｃ１、及びＲ０のそれぞれの値とＲ１、Ｃ１、及びＲ０のそれぞれに流れる電流が既知であるため、Ｃ１に印加される電圧であるＶ１値を算出することができる。
【００８３】
一方、前述した方法により算出されるエネルギー貯蔵システムの電圧状態変数は、下記の数式２を用いて、拡張カルマンフィルターにおいて時間更新されてもよい。
【００８４】
【数４】
2020530905.tif 000006
（但し、Ｒ１及びＣ１は、図５の回路モデルに含まれている抵抗及びキャパシタンスの値、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【００８５】
２．３     エネルギー貯蔵システムの充電容量を算出するエネルギー貯蔵システムの充電容量算出ステップ
一方、数式３は、本発明の拡張カルマンフィルターに状態変数として用いられる前記数式１（充電容量状態変数を時間更新する数式）及び数式２（電圧状態変数を時間更新する数式）をベクトル状態方程式で示す数式である。
【００８６】
【数５】
2020530905.tif 000007
（但し、Ｃａｐａｃｉｔｙ：二次電池の容量、Ｒ１及びＣ１は、図５の回路モデルに含まれている抵抗及びキャパシタンスの値、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【００８７】
前記数式３を拡張カルマンフィルターの初期値（ｋ＝０のとき）に取り込んでエネルギー貯蔵システムの充電容量を正確に算出することができる。
【００８８】
本発明において用いられる拡張カルマンフィルターは、次の手順（時間更新（推定）→測定更新（補正）→時間更新（推定）→測定更新（補正））を繰り返し行って算出される充電容量の誤差を減らすことができる。
【００８９】
【数６】
2020530905.tif 000008
（Ｐ：誤差共分散、Ｈ：変換係数、Ｋ：カルマンゲイン、Ｑ：真値に対する標準偏差Ａ：状態変数（数式３）、ｚ：観測値、Ｒ：観測真値との誤差、ｋ：ステップの次数、ｕ：追加入力値、Ｉ：単位行列）
【００９０】
そして、拡張カルマンフィルターの誤差は、入力される状態変数の誤差が小さいほど、さらに正確な値を算出することができる。
【００９１】
換言すれば、本発明は、前記低いｃ−ｒａｔｅの電流により算出された抵抗（Ｒ）及びキャパシタンス（Ｃ）の値で正確な電圧モデルを算出して拡張カルマンフィルターに代入することから、従来の技術により算出された電圧モデルに比べて誤差が減る。
【００９２】
＜具体的な実験データ＞
以下では、本発明の特徴である低いｃ−ｒａｔｅを用いて測定されたエネルギー貯蔵システムの抵抗（Ｒ）及びキャパシタンス（Ｃ）を用いてエネルギー貯蔵システムの電圧モデルを生成し、これに基づいてエネルギー貯蔵システムの充電容量を算出した結果データと、従来の技術により算出された充電容量の結果データとを比較して説明する。
【００９３】
表１及び図６は、本発明のエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗及びキャパシタンス算出方法を用いてエネルギー貯蔵システムの充電容量を算出する方法を用いて算出された充電容量と、従来の技術（電流積算方法のみを用いる）により算出された充電容量とを比較する表及び図である。
【００９４】
【表１】
2020530905.tif 000009
【００９５】
表１において、前記外部誤差は、電流センサーによる誤差である。一方、前記外部誤差の有無に応じた温度別の誤差値は、実際に測定した充電容量値と、本発明または従来の方式による算出方法により算出された充電容量値との差分である。
【００９６】
表１及び図６を参照すると、外部誤差が存在しない場合（電流センサーに誤差がない場合）には、従来の技術により算出されたＳＯＣ（充電容量）の誤差が本発明の方法により算出されたＳＯＣ誤差よりも小さいことを確認することができる。これに対し、外部誤差が存在する場合には、本発明の方法により算出されたＳＯＣ誤差が従来の技術により算出されたＳＯＣ誤差よりも小さいことを確認することができる。なお、外部誤差が存在するとき、ＳＯＣ誤差の増加率を見ても、本発明の方法により算出されたＳＯＣの誤差は最大で約３％であるのに対し、従来の技術により算出されたＳＣＯ誤差は最大で７％を超えるため、２倍以上の差分が出る。
【００９７】
したがって、外部誤差がない場合には従来の技術の方がより正確であるが、実際にエネルギー貯蔵システムに用いられる環境下では電流センサーには誤差が生じざるを得ないため、実際の使用環境下でエネルギー貯蔵システムの充電容量を算出する方法は、本発明の方法の方が従来の技術の方法よりも正確であるといえる。
【００９８】
一方、表２及び図７は、誤差条件を反映して、本発明のエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗及びキャパシタンス算出方法を用いてエネルギー貯蔵システムの充電容量を算出する方法を用いて算出された充電容量と、既存のアルゴリズム（高いｃ−ｒａｔｅにより算出された電圧モデル及びカルマンフィルターにおいて電圧モデルの比重が低いアルゴリズムを用いる従来の技術）により算出された充電容量とを比較する表及び図である。
【００９９】
【表２】
2020530905.tif 000010
【０１００】
表２及び図７を参照すると、電流センサーに＋−３％の誤差（外部誤差）がある場合（図７の（ａ））に、既存のアルゴリズムにより算出された充電容量の最大の誤差は２．５％であるのに対し、本発明の充電容量（ＳＯＣ）の算出方法を用いる場合（図７の（ｂ））、充電容量の最大の誤差は１％であって、既存のアルゴリズムよりも本発明の充電容量算出方法の方が、誤差がさらに小さいことを確認することができる。なお、エネルギー貯蔵システムの退化度（ＳＯＨ）誤差（外部誤差）が＋−５％であり、固定で発生する工程容量の許容誤差（外部誤差）の範囲が＋−２．５％である場合には、既存のアルゴリズムにより算出された充電容量の最大の誤差は９％であるのに対し、本発明の充電容量算出方法を用いる場合、充電容量の誤差は最大で１．５％であることを確認することができる。
【０１０１】
結果的に、外部の誤差が大きいほど、既存のアルゴリズムを用いて充電容量を測定する場合よりも本発明の充電容量算出方法を用いる場合の方が、誤差が少なく生じることを確認することができる。
【０１０２】
一方、表３及び図８は、既存のアルゴリズムと本発明の充電容量推定方法を用いて容量テストを行った結果とを比較した表及び図である。
【０１０３】
【表３】
2020530905.tif 000011
【０１０４】
表３及び図８を参照すると、既存のアルゴリズムを用いて容量テストを行ったところ、充電方向では最大で３．９％の誤差が、かつ、放電方向では最大で５％の誤差が生じたことを確認することができる。これに対し、本発明の充電容量算出方法を用いて容量テストを行ったところ、充電方向では最大で２％、放電方向では最大で２、３％の誤差が生じることを確認することができる。
【０１０５】
結果的に、本発明の充電容量算出方法を用いる方が、充放電の両方に対して容量テストを行ったところ、誤差が小さいことを確認することができる。
【０１０６】
一方、既存のアルゴリズムは、前術したように、最大の誤差が生じると、現在の充電容量を正確に判断することができないため、最大で約１分ほどエネルギー貯蔵システムが充電または放電を続けることがある。このように、充電または放電を続けると、エネルギー貯蔵システムの最大の充電容量または最大の放電容量を外れた範囲内においても充電または放電が続くため、エネルギー貯蔵システムにダメージを与える虞がある。
【０１０７】
３．本発明の実施の形態に係るエネルギー貯蔵システム
図９は、本発明の実施の形態に係るエネルギー貯蔵システムの構成を示す図である。
【０１０８】
以下では、図９に基づき、本発明の実施の形態に係るエネルギー貯蔵システムについて説明する。
【０１０９】
本発明の実施の形態に係るエネルギー貯蔵システム１０は、複数のバッテリーラック１１及び前記複数のバッテリーラックを制御するＢＳＣ（Ｂａｔｔｅｒｙ  Ｓｅｃｔｉｏｎ  Ｃｏｎｔｒｏｌｌｅｒ）１１０を備えてなってもよい。
【０１１０】
具体的に、前記ＢＳＣ  １１０は、エネルギー貯蔵システムの出力電圧を測定する電圧測定部１２０と、エネルギー貯蔵システムの出力電流を測定する電流測定部１３０と、一定の値を有する所定の電流により算出された温度別のエネルギー貯蔵システムの抵抗及びキャパシタンスルックアップテーブルが保存されている保存部１４０及び前記測定された電圧と電流及び算出された抵抗及びキャパシタンスに基づいて、エネルギー貯蔵システムの充電容量を算出する充電容量算出部１５０を備えてなってもよい。
【０１１１】
一定の値を有する電流としては、通算して、２〜３ｃ−ｒａｔｅ以上を有する電流を用いる。ところが、本発明の発明者は、２〜３ｃ−ｒａｔｅよりも低いｃ−ｒａｔｅの電流を印加して、エネルギー貯蔵システムの抵抗及びキャパシタンス値をより正確に算出することを想到した。
【０１１２】
したがって、本発明における前記低いｃ−ｒａｔｅの電流は、エネルギー貯蔵システムの抵抗及びキャパシタンスを測定可能な最小の電流値から２ｃ−ｒａｔｅまでの値を有する電流であってもよい。
【０１１３】
好ましくは、前記低いｃ−ｒａｔｅを有する電流は、０．１ｃ−ｒａｔｅの電流であってもよく、これは、１ｃ−ｒａｔｅの電流の１／１０に相当する値である。
【０１１４】
一方、前記電圧測定部１２０は、エネルギー貯蔵システムに所定の電流が印加される間のエネルギー貯蔵システムの電圧である第１の電圧及び前記エネルギー貯蔵システムに印加される電流が切れた後に、エネルギー貯蔵システムの電圧である第２の電圧を測定してもよい。
【０１１５】
一方、前記電流測定部１３０は、エネルギー貯蔵システムの出力電流を測定し、測定されたエネルギー貯蔵システムの出力電流は、後述するエネルギー貯蔵システムの充電容量状態変数値を算出するのに用いられる。
【０１１６】
前記保存部１４０に保存されているルックアップテーブルは、上述した本発明の実施の形態に係るエネルギー貯蔵システムの電圧モデルを算出するのに用いられる抵抗（Ｒ）及びキャパシタンス（Ｃ）を算出する方法に従って作成されて保存されてもよい。
【０１１７】
一方、前記充電容量算出部１５０は、前記ルックアップテーブルの抵抗及びキャパシタンスに基づいて、エネルギー貯蔵システムの電圧に対する状態変数を算出する電圧状態変数算出部１５１及び前記測定されたエネルギー貯蔵システムの出力電流を積算する電流積算法を用いて、エネルギー貯蔵システムの充電容量に対する状態変数を算出する充電容量状態変数算出部１５２を備えてなり、前記算出された電圧状態変数及び充電容量状態変数をカルマンフィルターに代入して、エネルギー貯蔵システムの充電容量を算出してもよい。
【０１１８】
具体的なエネルギー貯蔵システムの充電容量は、上述した本発明のエネルギー貯蔵システムの充電容量を算出する方法の手順を行って算出されてもよい。
【０１１９】
例えば、前記充電容量状態変数算出部において算出された充電容量状態変数は、下記の数式１に従って時間更新されてもよい。
【０１２０】
【数７】
2020530905.tif 000012
（但し、Ｑｃａｐａｃｉｔｙ：二次電池の容量、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【０１２１】
一方、前記電圧状態変数算出部において算出される前記算出された電圧状態変数は、下記の数式２に従って時間更新されてもよい。
【０１２２】
【数８】
2020530905.tif 000013
（但し、Ｒ１及びＣ１は、図５の回路モデルに含まれている抵抗及びキャパシタンスの値、ｋ：時間インデックス、Ｉ［ｋ］：時間インデックスｋから測定された電流、ｔ：時間）
【０１２３】
一方、本発明の技術的思想は、前記実施形態に基づいて具体的に記述されたが、前記実施形態はその説明のためのものであり、その制限のためのものではないということに留意すべきである。なお、本発明の技術分野における当業者であれば、本発明の技術思想の範囲内において種々の実施形態が実施可能であるということが理解できる筈である。
【符号の説明】
【０１２４】
１０  エネルギー貯蔵システム
１１  バッテリーラック
１１０  ＢＳＣ（Ｂａｔｔｅｒｙ  Ｓｅｃｔｉｏｎ  Ｃｏｎｔｒｏｌｌｅｒ）
１２０  電圧測定部
１３０  電流測定部
１４０  保存部
１５０  充電容量算出部
１５１  電圧状態変数算出部
１５２  充電容量状態変数算出部
２１０  開放電圧源
【図１】
2020530905.tif 000018
【図２】
2020530905.tif 000019
【図３】
2020530905.tif 000020
【図４】
2020530905.tif 000021
【図５】
2020530905.tif 000022
【図６】
2020530905.tif 000023
【図７】
2020530905.tif 000024
【図８】
2020530905.tif 000025
【図９】
2020530905.tif 000026
【国際調査報告】
2020530905.tif 000027
2020530905.tif 000028
2020530905.tif 000029
2020530905.tif 000030