(57)【特許請求の範囲】
【請求項１】
電気化学蓄電池のエネルギ状態を決定するための方法において、
電力（Ｐ）と、温度（Ｔ）と、エネルギ状態ＳＯＥと、残留エネルギＥｎとを含む電気化学蓄電池の動作点、又は電力（Ｐ）と、温度（Ｔ）と、エネルギ状態ＳＯＥと、蓄電池のエネルギ状態の関数である蓄電池内の残留エネルギＥｎの傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎとを含む電気化学蓄電池の動作点に関連する４値の所定の集合を準備するステップと、
メモリから第１のエネルギ状態ＳＯＥ１を読み出すステップ（Ｅ１０１）と、
前記電気化学蓄電池の現在の動作を表す温度（Ｔ１）及び電力（Ｐ１）を測定するステップ（Ｅ１０２）と、
前記４値の所定の集合に基づいて、前記測定された温度及び前記測定された電力（Ｔ１，Ｐ１）の関数としての前記第１のエネルギ状態ＳＯＥ１を表す点において、前記電気化学蓄電池のエネルギ状態ＳＯＥの関数としての蓄電池の残留エネルギＥｎの傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎを決定するステップ（Ｅ１０３）と、
前記第１のエネルギ状態ＳＯＥ１の前記決定された傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎ及びエネルギ量（Ｐ．ｄｔ）の関数として、第２のエネルギ状態ＳＯＥ２を決定するステップ（Ｅ１０４）とを有する方法。
【請求項２】
前記傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎは、前記４値の所定の集合に含まれる関連値における前記測定された温度（Ｔ１）、前記測定された電力（Ｐ１）及び前記第１のエネルギ状態ＳＯＥ１の少なくとも１つの概算によって求められる請求項１記載の方法。
【請求項３】
前記概算は、前記所定の集合の最も近い値の関連値を選択すること又は前記所定の集合の直下の関連値を選択することに対応している請求項２記載の方法。
【請求項４】
前記第２のエネルギ状態ＳＯＥ２は、充電フェーズにおける所定の期間に前記電気化学蓄電池に供給される電力の正の値、又は放電フェーズにおける所定の期間に前記電気化学蓄電池から出力され又は前記電気化学蓄電池から取り出される電力の負の値に関連するエネルギ量をＰ．ｄｔとして、式ＳＯＥ２≒ＳＯＥ１＋Ｐ．ｄｔ＊ｄＳＯＥ／ｄＥｎを適用することによって求められる請求項１乃至３いずれか１項記載の方法。
【請求項５】
前記方法は、繰返し実行され、各繰返しの最後に、メモリに格納されている前記第１のエネルギ状態ＳＯＥ１の値を、前記第２のエネルギ状態ＳＯＥ２の値に置換し、この値を次の繰返しにおいて、第１のエネルギ状態として使用する請求項１乃至４いずれか１項記載の方法。
【請求項６】
前記所定の期間（ｄｔ）は、前記方法の繰返しの期間を表す請求項４又は５記載の方法。
【請求項７】
前記傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎは、前記４値の所定の集合からの単純な読出によって取得される請求項１乃至６いずれか１項記載の方法。
【請求項８】
前記傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎは、前記所定の集合に基づく演算によって取得される請求項１乃至６いずれか１項記載の方法。
【請求項９】
前記傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎは、
前記測定された温度（Ｔ１）、前記測定された電力（Ｐ１）及び前記第１のエネルギ状態ＳＯＥ１について、それぞれ概算された温度（Ｔ１ａｐｐ）、概算された電力（Ｐ１ａｐｐ）及び概算された第１のエネルギ状態ＳＯＥ１ａｐｐを選択することによって概算を行い、前記概算された温度及び電力の値は、前記所定の集合に存在する値であって、対応する測定値に等しい、最近の又は直下の値であり、前記概算されたエネルギ状態ＳＯＥ１ａｐｐの値は、前記所定の集合に存在する値であって、前記第１のエネルギ状態ＳＯＥ１に等しい、最近の又は直下の値であり、
前記概算された温度（Ｔ１ａｐｐ）、前記概算された電力（Ｐ１ａｐｐ）及び前記概算された第１のエネルギ状態ＳＯＥ１ａｐｐの関数として、前記所定の集合から第１の残留エネルギＥｎ１抽出し、
前記概算された温度（Ｔ１ａｐｐ）、前記概算された電力（Ｐ１ａｐｐ）及び前記所定の集合に含まれ、前記概算された第１のエネルギ状態ＳＯＥ１ａｐｐに関連する行の直上にインデクスされた第２の概算されたエネルギ状態ＳＯＥ２ａｐｐの関数として、前記所定の集合から第２の残留エネルギＥｎ２を抽出し、
式ｄＥｎ／ｄＳＯＥ＝（Ｅｎ１−Ｅｎ２）／（ＳＯＥ１ａｐｐ−ＳＯＥ２ａｐｐ）を用いて、傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎを算出することによって取得される請求項８記載の方法。
【請求項１０】
電気化学蓄電池と前記電気化学蓄電池の電力（Ｐ）及び温度（Ｔ）を測定するセンサとを備え、
前記センサから取得した値に基づいて請求項１乃至９いずれか１項記載の方法を実現するコンピュータ。
【請求項１１】
請求項１乃至９いずれか１項記載の各ステップをコンピュータに実行させる、コンピュータが読み取り可能なプログラム。
【請求項１２】
請求項１１に記載のプログラムが格納された記録媒体。
【請求項１３】
電力（Ｐ）と、温度（Ｔ）と、エネルギ状態（ＳＯＥ）と、残留エネルギ（Ｅｎ）とを含む電気化学蓄電池の動作点又は電力（Ｐ）と、温度（Ｔ）と、エネルギ状態（ＳＯＥ）と、前記電気化学蓄電池のエネルギ状態ＳＯＥに基づく前記電気化学蓄電池の残留エネルギＥｎの傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎとを含む前記電気化学蓄電池の動作点に関連する第１の４値の集合を処理する方法において、
前記第１の４値の集合に基づいて、補間によって、第２の４値の集合を生成するステップ（Ｅ１）であって、
ａ．温度（Ｔ）の補間、特に線形補間を行い（Ｅ１−１）、
ｂ．電力（Ｐ）及びエネルギ状態（ＳＯＥ）の補間、特に３次スプライン補間を行う（Ｅ１−２）ステップ（Ｅ１）と、
電力（Ｐ）と、温度（Ｔ）と、エネルギ状態（ＳＯＥ）と、少なくとも前記第２の４値の集合に基づく前記電気化学蓄電池のエネルギ状態ＳＯＥの関数としての前記電気化学蓄電池の残留エネルギＥｎの傾斜（ｄＥｎ／ｄＳＯＥ，ｄＳＯＥ／ｄＥｎ）とを含む前記電気化学蓄電池の動作点に関連する第３の４値の集合を構成するステップ（Ｅ２）とを有する方法。
【請求項１４】
前記第３の４値の集合を構成するステップ（Ｅ２）は、前記第１及び第２の４値の集合の併合として定義される請求項１３記載の方法。
【請求項１５】
前記第１の４値の集合は、標準蓄電池に対して測定された実験データから導出される請求項１３又は１４記載の方法。
【請求項１６】
前記第１の４値の集合に基づいて温度の補間を行い（Ｅ１−１）、続いて、前記温度の補間（Ｅ１−１）によって得られた中間集合に基づいて電力及びエネルギ状態の補間を行う（Ｅ１−２）請求項１３乃至１５いずれか１項記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電気化学蓄電池の分野に関する。
【０００２】
詳しくは、本発明は、蓄電池のエネルギに関連する動作点を補足するデータの使用に関する。
【０００３】
これに関連し、本発明は、特に、電力、温度、エネルギ状態及び残留エネルギ、又は電力、温度、エネルギ状態及び傾斜を含む電気化学蓄電池の動作点に関連する４値（４個一組の値：クワドラプレット：quadruplet）の第１の集合を処理する方法に関する。
【背景技術】
【０００４】
従来の蓄電池状態インジケータは、蓄電池に蓄積されている電荷の量の評価に基づいている。蓄電池から取り出される及び／又は蓄電池に供給される電流の強度の測定値を積分することによって、「充電状態（State Of Charge：ＳＯＣ）」インジケータを作成することができる。
【０００５】
換言すれば、以下の公式が適用される。
【０００６】
【数１】
0006305988.tif 000002
【０００７】
ここで、Ｑは、時刻ｔにおいてバッテリに蓄積されている電荷の量をクーロンで表し、Ｑ０は、バッテリに蓄積された初期の電荷の量をクーロンで表し、Ｑｍａｘは、バッテリの最大の（フル充電時の）電荷の量をクーロンで表し、ＳＯＣは、充電状態をパーセンテージで表している。
【０００８】
この従来の充電状態インジケータは、蓄電池における損失、特に蓄電池の内部抵抗に起因する損失を考慮できない点が不十分である。
【０００９】
実際、蓄電池の内部抵抗が大きい程、再生されるエネルギ量が小さくなる。この場合、蓄電池が非常に大きな量の電荷を蓄積していても、実際に使用可能な電荷が少なくなる。したがって、蓄電池の内部抵抗が高い程、充電状態の値が不正確になる。
【００１０】
したがって、蓄電池の実際の状態を表す他のインジケータが求められている。
【００１１】
特許文献ＦＲ２９４７６３７は、蓄電池のエネルギ状態を特徴付ける方法を開示している。
【００１２】
この方法では、図１に示すように、３次元空間にマッピングして表現することができるＳＯＥ（Ｗｈで表されるエネルギ状態）、Ｐ（Ｗで表される取り出された使用可能な電力）及びＥｎ（Ｗｈで表される残留エネルギ）の値の集合を定義して、蓄電池の挙動の幾つかの特徴点を判定する。
【００１３】
エネルギ状態は、基準電力において使用可能なエネルギに関連する。この基準電力は、使用可能なエネルギが最大となる電力であってもよい。エネルギ状態（State Of Energy：ＳＯＥ）は、０〜１、すなわち、０％〜１００％の範囲で変化する。例えば、基準電力における基準エネルギが１０Ｗｈである蓄電池について、ＳＯＥ＝５０％、Ｐ＝２０Ｗ、Ｅｎ＝３Ｗｈの実験点を選択すると、これは、蓄電池が実際に２０Ｗの電力で使用されている場合、（５Ｗｈではなく）３Ｗｈの残留エネルギを提供できることを意味する。
【００１４】
このような推論を可能にする図１のＳＯＥ値は、標準の蓄電池又はバッテリを構成する一組の蓄電池から判定でき、したがって、蓄電池の電力及び残留エネルギを制御可能な環境下で、実験室で標準化することができる。
【００１５】
本発明は、特に、演算リソースが限定的である実時間のオンボードアプリケーションのコンテキストにおいて、これらのデータを用いる際に生じる問題を解決するものである。
【発明の概要】
【発明が解決しようとする課題】
【００１６】
本発明の目的の１つは、上述した問題を解決し、蓄電池の速やかな較正／初期化を可能にする解決策を提供することである。
【課題を解決するための手段】
【００１７】
本発明に基づく判定方法は、請求項１によって定義される。
【００１８】
判定方法の異なる実施形態は、請求項２乃至９によって定義される。
【００１９】
本発明に基づく判定デバイスは、請求項１０によって定義される。
【００２０】
本発明に基づくデータ記録媒体は、請求項１１によって定義される。
【００２１】
本発明に基づくコンピュータプログラムは、請求項１２によって定義される。
【００２２】
本発明に基づく処理方法は、請求項１３によって定義される。
【００２３】
処理方法の異なる実施形態は、請求項１４乃至１８によって定義される。
【００２４】
本発明の他の利点及び特徴は、以下の添付の図面を用いて説明する非制限的な具体例である本発明の特定の実施形態から明らかとなる。
【図面の簡単な説明】
【００２５】
【図１】蓄電池の動作点の分布を残留エネルギ、電力及びエネルギ状態の関数として示す図である。
【図２】蓄電池の動作温度を考慮に入れることによる図１の改善例を示す図である。
【図３】実験的な特徴付けによって、例えば、図２に示すようなデータの第１の集合を処理する方法を図式的に示す図である。
【図４】図２のデータから開始し、図３の方法の適用に続く補間による精密化を示す図である。
【図５】第３の４値の集合を構成する特定の実施形態を示す図である。
【図６】固定の電力及び温度におけるエネルギ状態の関数としての残留エネルギの動作点のモデリングを表すグラフを示す図である。
【図７】蓄電池のエネルギ状態を推定する方法のステップを示す図である。
【図８】傾斜値の算出のための可能な異なるステップを示す図である。
【図９】エネルギ状態を判定するための方法の妥当性を検証するステップにおけるエネルギ状態を時間の関数として表す曲線及び蓄電池の電圧の変化を時間の関数として表す曲線を示す図である。
【発明を実施するための形態】
【００２６】
オンボードアプリケーションの一部としてエネルギ状態を判定する場合、リソースの管理は、無視できないパラメータである。このため、電気化学蓄電池の動作点に関連する第１の値の処理を集合から開始することが提案されている。各点は、温度Ｔである追加的な変数を考慮に入れることによって組み合わされた電力Ｐ、エネルギ状態ＳＯＥ及び残留エネルギＥｎ、又は電力Ｐ、エネルギ状態ＳＯＥ及び傾斜ｄＥｎ／ｄＳＯＥを含む。電気化学蓄電池の内部抵抗の挙動には、温度が影響するため、ここでは、温度を組み込んでいる。
【００２７】
換言すれば、第１の集合は、例えば、テーブルＥｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）又はテーブルｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）の形式で表される４値（４個一組の値：クワドラプレット：quadruplet）を表している。
【００２８】
これらのデータは、図２に示すように、マッピングの形式で表現することができる。図２は、それぞれ、Ｗｈで表される残留エネルギ、Ｗで表される電力及び％で表されるエネルギ状態ＳＯＥによって定義される３次元空間を示している。この空間の各層は、そのメッシュ化の関数として、温度（この具体例では、５つの温度）に関連する仮想表面を画定している。各層の各メッシュは、４値（ＳＯＥ、残留エネルギ、電力、温度）の関数として、実験的測定によって判定された動作点に対応している。実験は、自動化しても、時間が長くかかるため、動作点の数は極めて少なくしてある。したがって、これらの初期のマッピングを補足する必要がある。
【００２９】
本発明の実施形態である方法のステップを詳細に説明する前に、幾つかの定義を説明する。
【００３０】
「エネルギ状態（State Of Energy）」ＳＯＥは、蓄電池の公称条件下におけるエネルギの放電の仮定下で使用可能な残留エネルギＥｄ／ＰＮと、公称エネルギＥＮｏｍとの間の比として定義され、したがって、式ＳＯＥ＝Ｅｄ／ＰＮ／ＥＮｏｍとして表される。このＳＯＥの値は、０〜１の範囲であり、１の値は、フル充電されたエネルギ状態に対応し、０の値は、完全な放電状態に対応する。また、この値は、パーセンテージで表してもよい。
【００３１】
電力Ｐは、蓄電池の製造業者が推奨する電力使用範囲内の値であり、この範囲は、製造業者が直接的に指定したものであってもよく、又は例えば、この製造業者によって指定された電流の範囲を指定された公称電圧に乗算することによって推定したものであってもよい。この電力は、蓄電池の使用の状態、すなわち、充電又は放電の関数である。放電とは、蓄電池から電力Ｐが取り出されることを意味し、充電とは、蓄電池に電力Ｐが供給されることを意味する。時刻ｔにおいて使用可能な電力Ｐは、エネルギ状態及び温度に依存することがある。
【００３２】
充電状態及び放電状態は、蓄電池の技術に基づいて判定される。これらは、蓄電池製造業者が推奨する情報から判定してもよく、通常、閾値電圧に基づいて判定される。
【００３３】
残留エネルギＥｎは、蓄電池に実際に蓄積されている内部エネルギと、蓄電池の内部抵抗のジュール効果によって失われたエネルギとを考慮に入れたＷｈで表される蓄電池の使用可能なエネルギである。これは、以下の式で表される。
Ｅｎ＝Ｅｉ−Ｅｐ
【００３４】
【数２】
0006305988.tif 000003
【００３５】
傾斜の概念は、エネルギ状態の関数としての残留エネルギに関連する。これは、好ましくは、固定の電力及び温度について、残留エネルギ／エネルギ状態の対を起点とし、同じ電力及び同じ温度に関連する全ての点を通過する曲線の残留エネルギ／エネルギ状態の対によって形成された点において局所的に評価される傾斜の値に対応している。すなわち、傾斜は、ｄＥｎ／ｄＳＯＥ又はその逆数であるｄＳＯＥ／ｄＥｎに関連付けることができる。包括的に言えば、この説明において、傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎは、蓄電池のエネルギ状態の関数としての蓄電池の残留エネルギの傾斜と読み替えることができる。
【００３６】
第１の４値の集合は、仏国特許出願ＦＲ２９４７６３７に記載されている手法に、更に、温度を考慮に入れて生成してもよい（図２）。したがって、この第１の集合の生成については、ここでは詳細に記述しない。この仏国特許出願には、ｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）を得るためのテーブルについては説明されていないが、この仏国特許出願の開示に基づいて情報を処理して、このテーブルを取得する手法は、当業者にとって明らかである。通常、初期のテーブルの各点において、傾斜（Ｅｎ（ｉ＋１）−Ｅｎ（ｉ））／（ＳＯＥ（ｉ＋１）−ＳＯＥ（ｉ））を算出でき、ここで、ｉは、初期のテーブルの行を表している。
【００３７】
図３は、例えば、Ｅｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプの電力Ｐ、温度Ｔ、エネルギ状態ＳＯＥ及び残留エネルギＥｎを含む、若しくはｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）又はこの逆数のタイプの電力Ｐ、温度Ｔ、エネルギ状態ＳＯＥ及び傾斜を含む電気化学蓄電池の動作点に関連する第１の４値の集合を処理する方法の具体例を示している。
【００３８】
この処理方法は、第１の４値の集合に基づいて、補間によって、第２の４値の集合を生成するステップＥ１を含み、このステップＥ１は、以下のサブステップを含む。
・温度Ｔについて、補間、特に線形補間を実行するサブステップＥ１−１
・電力Ｐ及びエネルギ状態ＳＯＥについて、補間、特に３次スプライン補間を実行するサブステップＥ１−２
温度の線形補間の利点は、演算の簡潔性のために、実効が容易である点である。電力及びエネルギ状態における３次スプライン補間の場合、より定常的で単調（regular and monotonic）な結果を得ることができる。
【００３９】
更に、方法は、特に第１及び第２の４値の集合に基づいて、第３の４値の集合を構成するステップＥ２を含む。
【００４０】
好ましくは、第３の４値の集合は、電力（Ｐ）と、温度（Ｔ）と、エネルギ状態（ＳＯＥ）と、少なくとも第２の４値の集合に基づく蓄電池のエネルギ状態ＳＯＥの関数としての蓄電池の残留エネルギＥｎの傾斜（ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎ）とを含む電気化学蓄電池の動作点に関連する値から構成される。
【００４１】
実際、第３の４値の集合を構成するこのステップＥ２は、より包括的に言えば、少なくとも第２の４値の集合に基づいて実現することができる。換言すれば、第３の４値の集合は、実際には、第２の４値の集合のみから構成できる。この実施形態では、第２の４値の集合を生成するステップＥ１と、第３の４値の集合を構成するステップＥ２は、単一の共通のステップを構成している。
【００４２】
一実施形態では、第３の４値の集合を構成するステップは、第１及び第２の４値の集合の併合として定義できる。
【００４３】
実際、「補間」の概念は、簡潔な統計的数値、この場合、第１の４値の集合に基づく判定によって構成され、新しい値は、実験的測定が実行されていない中間値に対応する。換言すれば、第１及び第２の４値の集合は、好ましくは、別個である。この結果、第１及び第２の集合を併合することによって、可能な限り多数の値を含む第３の集合を得ることができる。
【００４４】
このように、蓄電池の既知の動作点は、通常、実験データに由来し、標準の蓄電池について測定された第１の４値の集合に基づいて補足される。これは、オンボードアプリケーションのアップストリーム側の高性能マシンによって実現でき、したがって、オンボードアプリケーションは、多くの複雑な演算を実行することなく、これらの結果を利用することができる。
【００４５】
図３の特定の具体例は、まず、Ｅｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）によって構成される第１の４値の集合に基づいて、温度の補間、好ましくは、線形補間を行うことができることを示している。次に、温度の補間によって得られた中間集合に基づいて、電力とエネルギ状態の補間、例えば、３次スプライン補間を行う。したがって、第２の集合は、中間集合と、この中間集合の電力及びエネルギ状態の補間によって得られた集合との併合に対応している。
【００４６】
したがって、中間集合は、第１の集合と、温度の補間の間に得られたデータとの併合に対応している。したがって、第３の集合は、中間集合と、この中間集合の電力及びエネルギ状態の補間によって得られた集合との併合に対応している。
【００４７】
図３の特定の具体例では、実験的発見に基づいて、Ｅｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプの元のテーブルが得られ、ＳＯＥの６つの異なる値と、電力Ｐの６つの異なる値と、温度Ｔの８つの異なる値との組合せについて、残留エネルギ値を求める。
【００４８】
温度の補間、好ましくは、線形補間では、１℃毎に一連のデータを取得し、すなわち、ステップＥ１−１において、６個の電力×６個のエネルギ状態について、９１個のデータを連続的に取得する。
【００４９】
そして、これらの一連のデータに対して、２０個の電力及び２０個のエネルギ状態の補間、好ましくは、３次スプライン補間を行う（ステップＥ１−２）。
【００５０】
図４は、電力及びエネルギ状態ＳＯＥにおける補間の表現の具体例を示している。図４では、説明を明瞭にするために、温度の補間の結果を示しておらず、補間後の９１個の層のうち、実験による温度に関連する５つの層のみを示している。
【００５１】
このように、ＳＯＥの６個の点、電力Ｐの６個の点及び温度Ｔの８個の点、すなわち、エネルギの２８８個の点の実験データに基づいて、ＳＯＥの２０個の点、電力Ｐの２０個の点及び温度Ｔの９１個の点の集合、すなわち、エネルギの３６４００個の点を生成する。
【００５２】
倍精度浮動小数点符号化すなわち、１つの値あたり８バイトによる符号化では、第３の集合は、３６４００×８＝２９１２００バイトで表現される。
【００５３】
９１個の温度に亘る補間は、蓄電池の通常の動作範囲に亘って、１℃毎に直接的に行うことができるため、実用的である。当業者は、蓄電池の用途に応じて、温度の補間を適応化できる。
【００５４】
図２及び図４の比較から、図４の層では、不規則性が滑らかになっていることがわかり、この結果、充電状態を正確に推定することができる。
【００５５】
この具体例のマッピングを格納するために必要なメモリ量は、合計で３００キロバイトの桁であり、特にオンボードアプリケーションにおいて、処理の集中の問題は生じない。
【００５６】
補間された電力の段階の数をＮｐをとし、補間されたエネルギ状態の段階の数をＮｓｏｅとし、補間された温度の段階の数をＮｔｅｍｐとすると、メモリ内に格納される４値の数は、Ｎｔｅｍｐ＊Ｎｐ＊Ｎｓｏｅとなる。
【００５７】
以下の表１は、Ｎｔｅｍｐ、Ｎｓｏｅ及びＮｐの値及び符号化ビット数の関数として、第３の集合を格納するメモリサイズ（キロビット又はキロバイト）を評価した結果を示している。
【００５８】
【表１】
0006305988.tif 000004
【００５９】
６４ビットで１００×１００×９１といった、先験的で高精度なデータ補間及び符号化を行ったとしても、メモリ量は合理的な大きさ、すなわち、８メガバイト以下であり、電子カードに容易に格納することができる。
【００６０】
補間の精度は、要求されるエネルギ状態の推定の精度に応じて最適化する必要がある。使用電力及び温度に関するエネルギ状態関数の不規則性が高い程（図２の領域Ｚ参照）、より多数のモデリング点／動作点が望まれる。
【００６１】
また、不規則性が存在する箇所のみにおいて、より多くの数の補間を行うことによって、動作点の数を増やすこともできる。これによって、アプリケーションが実行される際のメモリ内の探索の簡潔性を犠牲にして、マッピングを格納するメモリのサイズを小さくすることができる。
【００６２】
マッピング、すなわち、補間から派生する追加的な動作点は、数値解析ソフトウェアを用いて、コンピュータによって生成できる。このようなソフトウェアとしては、例えば、ｍａｔｌａｂ、ｍａｔｈｃａｄ、ｏｃｔａｖｅ、ｓｃｉｌａｂ等を用いることができる。
【００６３】
図５に示す１つの特定の実施形態では、第１の４値の集合及び第２の４値の集合は、電力Ｐ、温度Ｔ、エネルギ状態ＳＯＥ及び残留エネルギＥｎを含む電気化学蓄電池の動作点、すなわち、例えば、Ｅｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプに関連する。第１及び第２の４値の集合に基づいて第３の４値の集合を構成するステップは、電力Ｐ、温度Ｔ、エネルギ状態ＳＯＥ及び傾斜ｄＥｎ／ｄＳＯＥを含む電気化学蓄電池の動作点、すなわち、例えば、ｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプの第３の４値の集合を判定するステップを含む。これは、温度の補間、及び電力及びエネルギ状態の補間の後に、少なくとも第１の集合及び第２の集合の併合を表す４値の集合ついて行うことができる（ステップＥ２−１）。この場合、好ましくは、以下の手法によって、第３の４値の集合の温度Ｔ及び電力Ｐの組合せ毎に傾斜ｄＥｎ／ｄＳＯＥを得ることができる。
・それぞれが残留エネルギ及びエネルギ状態によって構成される座標の点を構成する残留エネルギ／エネルギ状態の対の集合を判定する。
・対の集合の各対について、処理された対の点及び異なる対に関連する異なる点の関数として関連する傾斜値を評価する（Ｅ２−３）。
【００６４】
換言すれば、ステップＥ２−２では、エネルギ状態の関数として、残留エネルギのグラフを求め、残留エネルギ／エネルギ状態の各対は、グラフ上の点を構成し、この後のステップＥ２−３では、グラフの各点において、グラフ上の全ての点を通過する曲線を表す傾斜を局所的に評価する。図５のステップＥ２−２及び図６のグラフは、この方法を説明するために示しているだけであり、この方法では、ステップＥ２−２を行わなくてもよい。
【００６５】
当業者は、曲線内の前の又は次の点を考慮に入れることによって、周知の手法で傾斜の局所的な評価を行うことができる。傾斜は、好ましくは、正の値であり、傾斜が負の値である場合、これは、補間のエッジ効果の影響である。実際、負の傾斜は、物理的意味を有さず、これは、蓄電池の放電中に蓄電池のエネルギ状態が高まることを意味する。実際には、残留エネルギ及びエネルギ状態は、同じ方向に変化する。しかしながら、電力及びエネルギ状態の補間の後に幾つかの傾斜が負になることがあり、この場合、これらの傾斜は、エネルギが取り出されている（又は供給されている）際の弱い放電（又は充電）を反映するように０より僅かに大きい正値となるように制限される。
【００６６】
１つの変形例においては、第１の集合は、電力Ｐ、温度Ｔ、エネルギ状態ＳＯＥ及び傾斜ｄＥｎ／ｄＳＯＥを含む電気化学蓄電池の動作点に関連する４値に関連し、方法は、電力Ｐ、温度Ｔ、エネルギ状態ＳＯＥ及び残留エネルギＥｎを含む電気化学蓄電池の動作点に関連する第４の４値の集合に基づいて第１の集合を判定する先行するステップを含む。そして、傾斜は、上述と同様に判定することができる。
【００６７】
図６は、固定された１０℃の温度、１５Ｗの電力について、エネルギ状態ＳＯＥの関数として残留エネルギをＷｈで示すグラフである。点線は、グラフの点の全てを通過する仮想曲線を示している。
【００６８】
各点（残留エネルギ／エネルギ状態）において局所的に傾斜を評価した後、ステップＥ２−４において、ｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）を好適に求めることができる温度、エネルギ状態及び傾斜を組み合わせたテーブルを容易に得ることができ、このテーブルは、第３の集合を表す。
【００６９】
全ての変形例を含む上述した方法によって構成される第３の集合は、電気化学蓄電池のエネルギ状態を判定するための方法において包括的に使用することができる。この方法は、オンボードアプリケーションにおいて実時間で好適に実行される。
【００７０】
このように、エネルギ状態を判定するための方法は、好ましくは、第３の４値の集合を使用する。より包括的に言えば、方法は、電力Ｐと、温度Ｔと、エネルギ状態ＳＯＥと、残留エネルギＥｎとを含む電気化学蓄電池の動作点又は電力Ｐと、温度Ｔと、エネルギ状態ＳＯＥと、蓄電池のエネルギ状態の関数としての蓄電池の残留エネルギの傾斜、特に傾斜ｄＥｎ／ｄＳＯＥ又はｄＳＯＥ／ｄＥｎとを含む電気化学蓄電池の動作点に関連する４値の所定の集合を使用するステップを含む。この所定の集合は、第３の４値の集合であってもよく、他の手法で取得してもよい。したがって、以下の説明では、第３の４値の集合は、４値の所定の集合によって置き換えてもよい。
【００７１】
上述したような第３の集合への分割によって、モデリング及び複雑な演算を実時間アプリケーションの外部で行うことができ、実時間アプリケーションの実行の間には、以下に説明するように、比較的簡単な繰り返し演算及び幾つかの測定を実行すればよく、これによって、エネルギ状態インジケータにアクセスすることができ、エネルギ状態インジケータは、電池管理システム（Battery Management System：ＢＭＳ）の電子回路に容易に統合できる。
【００７２】
図７は、蓄電池のエネルギ状態を判定する方法の特定の具体例を示している。この方法は、好ましくは繰り返し行われ、メモリから第１のエネルギ状態ＳＯＥ１を読み出すステップＥ１０１を有する。この方法が繰り返される場合、ＳＯＥ１の値は、前のステップで判定された蓄電池のエネルギ状態に対応する。換言すれば、この方法は、必要な場合、各繰返しの最後に、メモリに格納されている第１のエネルギ状態ＳＯＥ１の値を、繰返しの間に判定されたエネルギ状態（現在のエネルギ状態）を表す第２のエネルギ状態ＳＯＥ２の値に置換するステップを有することが好ましく、このＳＯＥ２の値は、次の繰返しにおいて、第１のエネルギ状態として使用される。
【００７３】
最初の初期状態では、蓄電池は、最大限に充電されていてもよく、充電が終了すると、メモリの値は、１００％を示す。逆に、蓄電池は、完全に放電されていてもよく、この場合、初期化時点でメモリに保存される値は、０％を示していてもよい。
【００７４】
そして、ステップＥ１０２において、蓄電池の現在の動作を表す温度Ｔ１及び電力Ｐ１を測定する。「現在」という用語は、繰返しの間の動作状態を意味する。「蓄電池を表す温度及び電力」とは、蓄電池から取り出されるエネルギ又は蓄電池に供給されるエネルギの電力及び蓄電池の動作温度を意味する。ここでは、電力は、符号付きであり、正の場合と、負の場合がある。正の電力は、蓄電池が充電フェーズであることを表し、負の電力は、蓄電池が放電フェーズであることを表す。したがって、４値の集合は、正の電力値及び負の電力値を含むことができる。蓄電池を表す温度は、内部温度に最も近い必要があり、したがって、センサが電解質に耐性を有していれば、蓄電池内に温度センサを配置することができる。温度センサは、第１の集合を生成する際に使用した標準の蓄電池の温度センサと同じ位置に配置することが望ましいことは明らかである。テーブルは、正の電力値及び負の電力値を含んでもよいが、単純化のために、テーブルには、電力の絶対値を格納してもよい。
【００７５】
ステップＥ１０３では、第３の４値の集合に基づいて、測定された温度Ｔ１及び測定された電力Ｐ１の関数である第１のエネルギ状態ＳＯＥ１を表す点における蓄電池のエネルギ状態の関数として、蓄電池内の残留エネルギの傾斜ｄＥｎ／ｄＳＯＥ判定する。この傾斜ｄＥｎ／ｄＳＯＥは、例えば、第３の４値の集合がｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプである場合、単なる読出によって判定でき、第３の４値の集合がＥｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプである場合、計算によって取得することができる。
【００７６】
最後に、ステップＥ１０４において、第１のエネルギ状態ＳＯＥ１の傾斜ｄＥｎ／ｄＳＯＥ及びエネルギ量の関数として、第２のエネルギ状態ＳＯＥ２を判定する。
【００７７】
この第２のエネルギ状態ＳＯＥ２は、判定中の蓄電池の現在のエネルギ状態に対応しており、その動作変数、すなわち、測定された電力Ｐ１及び測定された温度Ｔ１を考慮に入れている。
【００７８】
オンボードアプリケーションのコンテキストにおける演算を制限するために、エネルギ状態ＳＯＥを判定する方法の異なるステップの間に、概算を行うことが好ましい。
【００７９】
すなわち、傾斜ｄＥｎ／ｄＳＯＥは、測定された温度Ｔ１、測定された電力Ｐ１、及び第３の４値の集合に含まれる関連値における第１のエネルギ状態ＳＯＥ１の少なくとも１つの概算によって得られる。特に、これらの値が第３の集合の部分として含まれない場合、これらの値の全てを概算によって求めることが好ましい。
【００８０】
実際、測定センサの分解能によっては、温度Ｔ１及び電力Ｐ１の測定値、又は読み出されたエネルギ状態ＳＯＥ１は、第３の４値の集合に含まれない値であることがある。オンボードアプリケーションのリソース使用を制限するために、充電の状態の判定の場合、離散化された値、すなわち、第３の集合に含まれる値を考慮に入れることが有利である。これを実現するため、概算は、第３の集合内の最も近い関連値を選択すること、又は第３の集合における直下の関連値を選択することによって実行してもよい。もちろん、第３の集合に含まれている値が出現した場合は、概算を行う必要はない。直下の値を選択するより、最も近い値を選択した方がより正確な結果を得ることができる。
【００８１】
第２のエネルギ状態ＳＯＥ２は、好ましくは、Ｐ．ｄｔをエネルギ量として、式ＳＯＥ２≒ＳＯＥ１＋Ｐ．ｄｔ＊ｄＳＯＥ／ｄＥｎを適用することによって求められる。実際、エネルギ量は、充電フェーズにおける所定の期間に蓄電池に供給される電力の正の値、又は放電フェーズにおける所定の期間に蓄電池から出力され又は蓄電池から取り出される電力の負の値に関連する。
【００８２】
したがって、この方法は、エネルギ状態ＳＯＥ２を判定する前に、蓄電池が充電フェーズであるか放電フェーズであるかを判定するステップを有していてもよい。この判定は、当業者に周知の如何なる手法で行ってもよく、例えば、電流を測定することによって行ってもよい。
【００８３】
以上のように、エネルギ量は、取り出される又は供給されるエネルギに対応している。値に関しては、エネルギ量は、Ｑｅｎｅｒｇｙ＝（Ｅｎ２−Ｅｎ１）に対応し、ここで、Ｅｎ１は、ＳＯＥ１に関連する残留エネルギであり、Ｅｎ２は、ＳＯＥ２に関連する残留エネルギである（図６参照）。Ｑｅｎｅｒｇｙは、Ｑｅｎｅｒｇｙ＝Ｕ．Ｉ．ｄｔ＝Ｐ．ｄｔによって算出できる。時間導関数ｄｔは、判定された期間の長さに対応する。上述したように、方法は、繰り返してもよく、判定された期間とは、方法における繰返し期間を表し、又はこの期間に等しい。傾斜を算出する必要がある特定の場合の図７のステップＥ１０３の詳細を図８に示す。まず、ステップＥ１０３−１において、測定された温度Ｔ１、測定された電力Ｐ１及び第１のエネルギ状態ＳＯＥ１について、それぞれ概算された温度Ｔ１ａｐｐ、概算された電力Ｐ１ａｐｐ及び概算された第１はエネルギ状態ＳＯＥ１ａｐｐを選択することによって、概算を実行する。概算された温度及び概算された電力については、対応する値の概算は、例えば、Ｅｎ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプの第３の４値の集合内に存在する測定値に等しい、測定値の直下の、又は測定値に最も近い対応する値を選択することによって行うことができる。ＳＯＥ１ａｐｐについては、第３の集合に含まれるＳＯＥ１の値に等しい値又はＳＯＥ１の直下の値を選択することが好ましい。第３の４値の集合は、既知の保存された値の組合せを含んでいるため、ステップＥ１０３−２では、第３の４値の集合から、簡単な読出によって、概算された温度Ｔ１ａｐｐ、概算された電力Ｐ１ａｐｐ及び概算された第１のエネルギ状態ＳＯＥ１ａｐｐの関数である第１の残留エネルギＥｎ１を抽出することができる。更に、ステップＥ１０３−３では、上述と同じ理由により、第３の４値の集合から、簡単な読出によって、概算された温度Ｔ１ａｐｐと、概算された電力Ｐ１ａｐｐと、第３の集合に含まれる第１の概算されたエネルギ状態ＳＯＥ１ａｐｐに関連する行の直上の行にインデクスされた概算された第２のエネルギ状態ＳＯＥ２ａｐｐとの関数である第２の残留エネルギＥｎ２を抽出することができる。「上位のインデクス」は、所与のＴ１ａｐｐ及びＰ１ａｐｐの組合せについて、ＳＯＥ１ａｐｐの直上の値ＳＯＥ２ａｐｐを意味する。図８では、ステップＥ１０３−３の前にステップＥ１０３−２を実行しているが、傾斜値を算出するステップＥ１０３−４値の時点で、関連するＥｎ１、Ｅｎ２、ＳＯＥ１ａｐｐ及びＳＯＥ２ａｐｐが存在していれば、この順序は任意である。
【００８４】
最後に、ステップＥ１０３−４において、式ｄＥｎ／ｄＳＯＥ＝（Ｅｎ１−Ｅｎ２）／（ＳＯＥ１ａｐｐ−ＳＯＥ２ａｐｐ）を用いて、傾斜ｄＥｎ／ｄＳＯＥを算出する。そして、この傾斜値は、図７のステップＥ１０４で使用することができる。
【００８５】
以上のように、ｄＥｎ／ｄＳＯＥ＝ｆ（ＳＯＥ，Ｐ，Ｔ）のタイプの第３の集合を用いることにより、傾斜を直接的に取得できるため、オンボードアプリケーションのコンテキストにおいて有利である。これによって、エネルギ状態を判定する方法を実行するオンボードコンピュータのリソース使用を制限でき、蓄電池のエネルギ状態の現在の値をより速やかに取得することができる。
【００８６】
２つの連続する繰返し間で、電力又は温度が変化することもある。この特定の場合、繰返しの間、現在のエネルギ状態の判定は、前の繰返しからの温度及び／又は電力の変化を考慮に入れる。
【００８７】
実際、メモリを読み出して値ＳＯＥ１を判定することによって、及び第３の４値の集合に基づいて、測定された温度及び測定された電力の関数として傾斜を判定することによって、温度又は電力の変化が自動的に考慮に入れられる。前の繰返しの演算からは、エネルギ状態（新しいステップの間のＳＯＥ１）だけが引き継がれる。新たな測定値に基づいて、新たな電力及び新たな温度を判定し、これによって、テーブルにアクセスして、残留エネルギ又は傾斜を判定する判定し、新たなエネルギ状態ＳＯＥ２を算出することができる。このように、電力及び温度の変動は、原理的に考慮に入れられる。
【００８８】
式ＳＯＥ２≒ＳＯＥ１＋Ｐ．ｄｔ＊ｄＳＯＥ／ｄＥｎは、温度補正することもできる。この場合、この式は、ＳＯＥ２≒ＳＯＥ１＋ｅｆｆＣｈ＊Ｐ．ｄｔ＊ｄＳＯＥ／ｄＥｎに置き換えられ、ここで、ｅｆｆＣｈは、充電又は放電の関数である補正係数である。補正係数を使用することによって、充電に関連するテーブル及び放電に関連するテーブルの使用を回避でき、したがって、データサイズが小さくなるために有益である。
【００８９】
電力は、蓄電池電圧と、蓄電池を流れる電流の値に基づいて測定することができる。
【００９０】
コンピュータプログラムが記録されているコンピュータ可読データ記録媒体は、上述した処理及び／又は判定方法のステップ及び／又はサブステップを実行するコンピュータプログラムコード手段を含むことができる。
【００９１】
同様にコンピュータプログラムは、プログラムがコンピュータによって実行されると、この判定方法及び／又は処理方法のステップ及び／又はサブステップを実行するように適応化されたコンピュータプログラムコード手段を含むことができる。また、本発明は、蓄電池のエネルギ状態を判定するデバイスを提供し、このデバイスは、判定方法及び／又は処理方法のステップを実行する（すなわち、判定方法及び／又は処理方法を実行する）ハードウェア及び／又はソフトウェア手段を含む。デバイスは、通常、メモリを備え、メモリには、例えば、データベースの形式で第３の集合が格納され、このデータベースの主キー（primary key）は、単一の残留エネルギ値又は単一の傾斜値を定義する電力、温度及びエネルギ状態である。このデバイスは、上述した記録媒体及び／又はコンピュータプログラムを含んでいてもよい。
【００９２】
以上の説明から明らかなように、記録媒体のコンピュータプログラムは、判定方法及び／又は処理方法を実行するために、上述したデバイスのソフトウェア手段によって実行可能なコンピュータプログラムコード手段を含むことができる。
【００９３】
更に、コンピュータプログラムは、特に、プログラムがコンピュータによって実行されると、この判定方法及び／又は処理方法を実現する上述したデバイスのソフトウェア手段で実行可能なコンピュータプログラムコード手段を含むことができる。
【００９４】
したがって、デバイスは、蓄電池の正確な動作点を実時間で速やかに判定でき、正確で確実なエネルギ状態値を抽出することができる。
【００９５】
エネルギ状態を判定する方法を使用プロファイルについて検査し、有効性を確認した。検査条件は、以下の通りとした。
・充電フェーズ及び放電フェーズを有する電力プロファイルを蓄電池に印加した。
・蓄電池の端子における電圧及び蓄電池の温度を測定した。
【００９６】
更に、上述したエネルギ状態を推定するために、蓄電池の電力プロファイル及び動作温度を演算アルゴリズムのシミュレーションに入力した。
【００９７】
そして、エネルギ状態のシミュレート結果と、実際に測定された電圧とを比較した。
【００９８】
検査の結果を図９に示す。この図９は、時間の関数としてのエネルギ状態ＳＯＥ、及び時間の関数としての蓄電池の端子における電圧の変化を示している。この図９は、放電が終了し、電圧が２Ｖ近傍に達すると、エネルギ状態（ＳＯＥ）が実際に約０に近付く（１．４８％）ことを示している。誤差の幅は、非常に小さく、オンボードアプリケーションのコンテキストでは、十分である。
【００９９】
ここでは、電気化学蓄電池について説明した。蓄電池の定義は、広義に解釈され、基本素子としての蓄電池及びバッテリの形式に構成された複数の蓄電池素子の両方を含むものとする。参照用の蓄電池として、Ａ１２３システムズ社（A123systems）の製品番号ＡＮＲ２６６５０Ｍ１を検査のために使用した。
【図１】
0006305988.tif 000005
【図２】
0006305988.tif 000006
【図３】
0006305988.tif 000007
【図４】
0006305988.tif 000008
【図５】
0006305988.tif 000009
【図６】
0006305988.tif 000010
【図７】
0006305988.tif 000011
【図８】
0006305988.tif 000012
【図９】
0006305988.tif 000013