(57)【要約】  電池などの電気エネルギーの電気化学的貯蔵用の装置の内部状態を予測する方法。  同じ種類の複数の電池が分析されるときの電池の様々な内部状態について、複数の電池を通過する電気信号を電流に加えることによってインピーダンス計測が実施され、それからこれらのインピーダンスをモデル化するようにＲＣ回路が定義される。次に、多変量統計解析を使用して、ＳｏＣ（とＳｏＨの少なくともいずれか）とＲＣ回路の複数のパラメータとの間の関係が較正される。分析される電池のインピーダンスの計測が実施され、これはＲＣ回路を使用してモデル化される。最後に、この関係を分析される電池について定義された等価電気回路の複数のパラメータに適用することによって、電池の内部状態が予測される。
2013519893.tif 000002
【特許請求の範囲】
【請求項１】
電池のような電気エネルギーの貯蔵のための第１の電気化学装置の内部状態を予測する方法であって、前記第１の電気化学装置の前記内部状態に関する少なくとも１つの特性が、インピーダンス分光法によって得られた電気計測から予測される方法において、
前記第１の電気化学装置と同じ種類の少なくとも１つの第２の電気化学装置のさまざまな内部状態について、前記第２の電気化学装置の前記内部状態に関連している前記特性を計測し、インピーダンス分光法を使用してさまざまな周波数で前記第２の電気化学装置の電気計測を実施するステップと、
前記第２の電気化学装置の複数の電気応答をモデル化する少なくとも１つのパラメータを有する等価電気回路を定めるステップと、
前記複数の内部状態について求められた前記特性および前記パラメータの複数の値の統計解析を使用して前記特性と前記等価電気回路の前記パラメータとの間の関係を較正するステップと、
前記等価電気回路の電気応答が前記第１の電気化学装置の前記電気応答と等しくなるように前記パラメータを求めることで前記等価電気回路を使ってモデル化されている様々な周波数について前記第１の電気化学装置の電気応答を求めるステップと、
前記関係を使用して前記電気化学装置の前記内部状態に関連している前記特性を計算することによって前記第１の電気化学装置の前記内部状態を予測するステップと、
を有することを特徴とする方法。
【請求項２】
前記様々な内部状態は前記第１の電気化学装置と同じ種類の前記電気エネルギーの貯蔵のための前記第２の電気化学装置の加速されたエージングを実施することによって得られる、請求項１に記載の方法。
【請求項３】
前記様々な内部状態は、前記第１の電気化学装置と同じ種類の一式の第２の電気化学装置を選択することによって得られ、前記一式の第２の電気化学装置は異なる内部状態を有している、請求項１に記載の方法。
【請求項４】
前記電気化学装置の前記内部状態に関連している前記特性である前記電気化学装置の充電状態（ＳｏＣ）と前記電気化学装置の劣化状態（ＳｏＨ）のうち少なくとも１つが計算される、請求項１ないし３のいずれか１項に記載の方法。
【請求項５】
前記等価電気回路は、以下のパラメータ、すなわち、抵抗、容量、温度、または前記パラメータの任意の組み合わせから選択された複数のパラメータによって定義されている、請求項１ないし４のいずれか１項に記載の方法。
【請求項６】
さまざまな周波数に対する電気応答は、前記電気化学装置を通過する電流に電気信号を加えることによって得られる複数の電気インピーダンス線図を計測することによって求められる、請求項１ないし５のいずれか１項に記載の方法。
【請求項７】
前記複数の電気インピーダンス線図は、正弦波電流の外乱を前記電気化学装置に加えて、前記電気化学装置の複数の端子の位置に誘起される正弦波電圧を計測することによって計測される、請求項６に記載の方法。
【請求項８】
前記複数の電気インピーダンス線図は、複数の正弦曲線の重ね合わせの形態またはホワイトノイズの形態の外乱を前記電気化学装置に加えて、前記電気化学装置の複数の端子の位置に誘起される正弦波電圧を計測することによって計測される、請求項６に記載の方法。
【請求項９】
前記電気化学装置は休止しているか動作している、請求項１ないし８のいずれか１項に記載の方法。
【請求項１０】
電気エネルギーの貯蔵のための電気化学装置の内部状態を推定する装置であって、
インピーダンス分光法によって前記電気化学装置の電気インピーダンスを計測する手段を有するセンサ（Ｇ）と、
等価電気回路と、前記電気化学装置の前記内部状態に関する特性と前記等価電気回路の複数のパラメータとの間の関係であって前記電気化学装置と同じ種類の少なくとも１つの第２の電気化学装置のさまざまな内部状態の複数の計測によって較正されている関係と、を保存するのに使用可能なメモリと、
前記電気化学装置の電気応答をモデル化する前記等価電気回路の前記複数のパラメータを定義する手段と、
前記関係を使用して前記電気化学装置の前記内部状態に関する特性を計算する手段と、
を有する装置。
【請求項１１】
電気インピーダンスを計測する前記手段は、
正弦波電流の外乱、複数の正弦曲線の重ね合わせの形態の外乱、またはホワイトノイズの形態の外乱を前記電気化学装置に加えるガルバノスタットと、
前記電気化学装置の複数の端子の位置に誘起される正弦波電圧を計測する手段と、
を有する、請求項１０に記載の装置。
【請求項１２】
請求項１０または１１に記載の電池の内部状態を推定する装置を有する、電池管理装置。
【請求項１３】
電池と、請求項１２に記載の電池管理装置を有する乗り物。
【請求項１４】
請求項１０または１１に記載の内部状態を推定する装置を有している電気エネルギーを蓄える光起電力装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電池（鉛、Ｎｉ−ＭＨ、リチウムイオン等）などの電気エネルギーの貯蔵のための電気化学装置の内部状態を予測する方法に関する。本方法は、静止または搭載の用途において、特に動作中に電池を管理するために使用することができる。
【背景技術】
【０００２】
電池は、ハイブリッド乗り物または電気乗り物の用途の場合に、または光起電力太陽エネルギーの貯蔵にとって、最も重要な構成要素の１つである。これらの用途の適切な動作は、様々な動的な負荷レベルの間の最良の妥協点で電池を動作させることに関する電池管理装置（ＢＭＳ）に依存している。そのＢＭＳには、充電状態（ＳｏＣ）と劣化状態（ＳｏＨ）の正確で信頼性のある知識が必要である。
【０００３】
電池の充電状態（ＳｏＣ）は、その利用可能な容量に該当し、製造者によって示されているその呼び容量のパーセントとして、または計測が可能なときに所与の条件下で計測されたその総容量のパーセントとして表される。ＳｏＣを知ることは、電池の次の再充電前に電池が所与の電流で継続してエネルギーを供給できる時間や、電池の次の放電前に電池がエネルギーを吸収できるまでの時間を予測することができることを意味する。この情報は、電池を使用する装置の動作の条件を決定する。
【０００４】
電池の寿命の間、電池の性能は、電池が使用できなくなるまで、使用中に発生する物理的および化学的変動のせいで徐々に低下する傾向がある。劣化状態（ＳｏＨ）は電池の消耗の状態を表している。このパラメータは、電池の寿命の間の時間ｔでの電池の総容量に該当しており、製造者によって示されている呼び容量または所与の複数の条件下で寿命の開始時に計測された容量に等しい、寿命の開始時に定められる全容量のパーセントで表されている。
【０００５】
乗り物についてのＳｏＣとＳｏＨの正確で信頼性のある予測値は、たとえば、乗り物の運転者が電池のエネルギーの潜在力を使用することに過度に慎重になったりその反対になったりしなくてもよいことを意味している。充電状態の不十分な診断によって、駆動できるキロメートル数が過度に評価されて、運転者にとって問題を引き起こすかもしれない。これらの指標の良好な予測は、電池を安全上の理由で過度に大きくする必要がなく、搭載される重量を削減し、その結果、燃料消費を削減することも意味する。ＳｏＣとＳｏＨの予測は、乗り物の総コストを削減することもできる。したがって、正確な予測器は、乗り物の動作の全範囲にわたる効率が高く安全な電池の容量の使用を保証する。
【０００６】
電池の充電状態（ＳｏＣ）と劣化状態（ＳｏＨ）とを予測するために多数の方法が公知である。
【０００７】
公知の方法の例には「クーロン計数」または「記帳」の方法がある。しかし、そのような方法は、自己放電などの現象を無視するため、予測に誤差が生じる。開放電圧がＳｏＣの指標として計測される方法が公知である。他の指標、例えば内部抵抗（米国特許第６１９１５９０号明細書、欧州特許出願公開第１８３５２９７号明細書）の使用も知られている。
【０００８】
これら２つの方法は、ＳｏＣは、スタティックチャートまたは解析的な機能依存性によって１つまたは２つ以上の計測可能なまたは容易に予測される量（電位、内部抵抗）に関連付けられていることを特徴としている。しかし、実際には、そのような依存性は、ＢＭＳにおいて通常考慮されるものよりもはるかに複雑であって、ＳｏＣの推定誤差につながることが多い。
【０００９】
潜在的により有望な方法は、インピーダンス分光法（ＥＩＳ）を使用して、ＳｏＣによって左右される量を計測することに基づいている。例として、米国特許出願公開第２００７／００９０８４３号明細書は、ＥＩＳによって容量性／誘導性転移に関連付けられている周波数ｆ±を求めることを提案している。周波数ｆ±とＳｏＣとの間の相関は、鉛電池だけでなくニッカド電池、およびニッケル水素電池についても示されている。同様なアプローチは、米国特許第６７７８９１３号明細書に記載されているように、構成要素がＳｏＣによって左右される等価電気回路によるＥＩＳ分光のモデル化に基づいており、それは、酸−鉛の対についての複数の周波数の電気化学インピーダンス分光法に基づいている自動車用バッテリーテスター、Ｓｐｅｃｔｒｏ  ＣＡ−１２（Ｃａｄｅｘ  Ｅｌｅｃｔｒｏｎｉｃｓ  Ｉｎｃ、カナダ）の開発につながった。ＥＩＳスペクトルは、等価電気回路によって検討され、構成要素の変化は、ＳｏＣによって左右される。同様に、米国特許第６０３７７７７号明細書において、電池の充電状態と他の特性は、鉛電池または他の装置について、複素インピーダンス／アドミタンスの実部と虚部とを計測することによって求められる。ＲＣモデルの使用は、欧州特許出願公開第０８８０７１０号明細書にも記載されており、電極の位置と電解質内での電気化学的および物理的な現象の記述がＲＣモデルの開発を支援する役割を果たし、外部の計測値に対しての精度を上げるために、電池の温度がモデルによってシミュレーションされている。
【００１０】
代替のアプローチは、他の分野で公知の推定技法を使用するために、電池の数学的モデルに基づいている。米国特許出願公開第２００７／００３５３０７号明細書は、数学的な電池モデルを使用して動作データ（電圧Ｕ、電流Ｉ、Ｔ）から電池の複数のパラメータと複数の状態変数とを推定する方法を特に記述している。数学的モデルは、複数の数学的サブモデルを有し、より高速に応答することを特徴としている。サブモデルは、限定された周波数範囲に対応しているＲＣモデルとも呼ばれる等価電気回路型のモデルである。
【００１１】
文献（Ｇｕ、Ｗｈｉｌｔｅ他）において知られている他のＳｏＣ推定方法は、電気化学装置の反応の数学的記述に基づいている。ＳｏＣは、装置の複数の状態変数から計算される。その記述は、材料、変化、エネルギー等のバランスだけでなく半経験的な相関に基づいている。
【００１２】
文献において知られているＳｏＨを予測する方法に関して、国際公開第２００９／０３６４４４号の著者等は、電極の劣化反応を観察するために基準電極を業務用の要素に導入している。しかし、この方法は、大量の計装、特に基準電極を要素に導入するための大量の計装と、電池のより複雑な電子管理を要求する。
【００１３】
仏国特許出願公開第２８７４７０１号明細書は、得られた応答を基準応答と比較するために時間的な電気的外乱を使用する方法を記述している。しかし、その方法は、そのような種類の外乱に応答して非常に弱い変動を示し従ってＳｏＨを正確に計測できないリチウムイオン型の複数の要素について実装がより困難である。
【００１４】
インピーダンス分析も、文献に記載されている。Ｕ．Ｔｒｏｌｔｚｓｃｈら（Ｅｌｅｃｔｒｏｃｈｉｍｉｃａ  Ａｃｔａ  ５１、２００６年、１６６４−１６７２頁）は、要素の劣化状態を求めるために、電気モデルに従ったインピーダンスの調整と組み合わされているインピーダンス分光法を使用する方法を記載している。しかし、その技法には、計測を行うために電気要素の使用を停止することが必要である。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】米国特許第６１９１５９０号明細書
【特許文献２】欧州特許出願公開第１８３５２９７号明細書
【特許文献３】米国特許出願公開第２００７／００９０８４３号明細書
【特許文献４】米国特許第６７７８９１３号明細書
【特許文献５】米国特許第６０３７７７７号明細書
【特許文献６】欧州特許出願公開第０８８０７１０号明細書
【特許文献７】米国特許出願公開第２００７／００３５３０７号明細書
【特許文献８】国際公開第２００９／０３６４４４号
【特許文献９】仏国特許出願公開第２８７４７０１号明細書
【非特許文献】
【００１６】
【非特許文献１】Ｇｕ、Ｗｈｉｌｔｅ他
【非特許文献２】Ｕ  Ｔｒｏｌｔｚｓｃｈ他、Ｅｌｅｃｔｒｏｃｈｉｍｉｃａ  Ａｃｔａ（電気化学会報）  ５１、２００６年、１６６４−１６７２頁
【発明の概要】
【発明が解決しようとする課題】
【００１７】
したがって、本発明は、電池などの電気エネルギーの貯蔵のための電気化学装置の内部状態を予測する代替の方法の提供に関する。本方法は、電池とその用途のモデルの機能として、所与の統計的なモデルによって装置の内部状態を再構築するための装置のインピーダンスの計測に基づいている。特に、本方法は、電気化学電池の充電状態（ＳｏＣ）と劣化状態（ＳｏＨ）の予測に使用することが可能で、これらは、静止しているか搭載されているかに関わらず電池を使用する大部分の用途において、複数の最も重要な内部特性である。
【課題を解決するための手段】
【００１８】
本発明は、第１の電気化学装置の内部状態に関する少なくとも１つの特性がインピーダンス分光学によって得られた電気計測から予測される、電池のような電気エネルギーの貯蔵のための第１の電気化学装置の内部状態を予測する方法に関する。本方法は、以下のステップ、すなわち、
第１の電気化学装置と同じ種類の少なくとも１つの第２の電気化学装置のさまざまな内部状態について、第２の装置の内部状態に関連している特性を計測し、インピーダンス分光法を使用してさまざまな周波数で第２の電気化学装置の電気計測を実施するステップと、
第２の装置の複数の電気応答をモデル化する少なくとも１つのパラメータを有する等価電気回路を定めるステップと、
複数の内部状態について求められた特性およびパラメータの複数の値の統計解析を使用して特性と等価回路のパラメータとの間の関係を較正するステップと、
等価電気回路の電気応答が第１の電気化学装置の電気応答と等しくなるように、パラメータを求めることで等価電気回路を使ってモデル化されている様々な周波数について第１の電気化学装置の電気応答を求めるステップと、
関係を使用して電気化学装置の内部状態に関連している特性を計算することによって第１の電気化学装置の内部状態を予測するステップと、
を有している。
【００１９】
本発明によれば、様々な内部状態は第１の電気化学装置と同じ種類の電気エネルギーの貯蔵のための第２の電気化学装置の加速されたエージングを実施することによって得られる。様々な内部状態は、第１の電気化学装置と同じ種類であって、異なる内部状態を有している一式の第２の電気化学装置を選択することによって得られてもよい。
【００２０】
電気化学装置の内部状態に関連している以下の特性、すなわち装置の充電状態（ＳｏＣ）および装置の劣化状態（ＳｏＨ）の少なくとも１つを計算することができる。
【００２１】
等価電気回路は、以下のパラメータ、すなわち、抵抗、容量、温度、または前記パラメータの任意の組み合わせから選択された複数のパラメータによって定義されてもよい。
【００２２】
本発明によれば、電気化学装置を通過する電流に電気信号を加えることによって得られる複数の電気インピーダンス線図を計測することによってさまざまな周波数に対する電気応答を求めることができる。複数の電気インピーダンス線図は、正弦波電流の外乱を電気化学装置に加えて、電気化学装置の複数の端子の位置に誘起される正弦波電圧を計測することによって計測してもよい。これらの電気インピーダンス線図は、複数の正弦曲線の重ね合わせの形態またはホワイトノイズの形態の外乱を電気化学装置に加えて、電気化学装置の複数の端子の位置に誘起される正弦波電圧を計測することによって計測してもよい。
【００２３】
本発明によれば、電気化学装置は休止（乗り物が停止または静止）していてもよいし、動作していてもよい。
【００２４】
本発明は電気エネルギーの貯蔵のための電気化学装置の内部状態を推定する装置であって、
インピーダンス分光法によって電気化学装置の電気インピーダンスを計測する手段を有するセンサ（Ｇ）と、
等価電気回路と、電気化学装置の内部状態に関する特性と等価電気回路の複数のパラメータとの間の関係であって電気化学装置と同じ種類の少なくとも１つの第２の電気化学装置のさまざまな内部状態の複数の計測によって較正されている関係と、を保存するのに使用可能なメモリと、
電気化学装置の電気応答をモデル化する等価電気回路の複数のパラメータを定義する手段と、
前記関係を使用して電気化学装置の内部状態に関する特性を計算する手段と、
を有する装置にも関する。
【００２５】
本発明によれば、電気インピーダンスを計測する手段は、
正弦波電流の外乱、複数の正弦曲線の重ね合わせの形態の外乱、またはホワイトノイズの形態の外乱を電気化学装置に加えるガルバノスタットと、
電気化学装置の複数の端子の位置に誘起された正弦波電圧を計測する手段と、
を有している。
【００２６】
本発明は、本発明の電池の内部状態を推定する装置を有する電池管理装置にも関する。
【００２７】
本発明は、電池と、本発明の電池管理装置とを有する乗り物にも関する。
【００２８】
本発明は、本発明の内部状態を予測する装置を有する、電気エネルギーを貯蔵する光起電力装置にも関する。
【００２９】
本発明の方法のさらなる特徴と利点とは、以下に説明する添付の図面を参照して行う実施形態の非限定的な例の以下の説明から明らかになるであろう。
【図面の簡単な説明】
【００３０】
【図１】本発明の方法の論理図である。
【図２】インピーダンス計測値を使用した検査手続きの図である。
【図３】Ｌｉ４Ｔｉ５Ｏ１２／ＬｉＦｅＰＯ４型の電池についての２０％の充電状態でのＶＥＨ用途を表している５〜６個のエージングの状態について得られたインピーダンスの比較の図である。
【図４】電気化学蓄電池を表している等価電気回路の例の図である。
【図５】ナイキスト表現ａ）とボード表現ｂ）における、そして図４の等価回路モデルを使用した２０％ＳｏＣでのＬｉ４Ｔｉ５Ｏ１２／ＬｉＦｅＰＯ４電池についての６５ｋＨｚと０．１Ｈｚとの間のインピーダンスのモデルついての調整の例の図である。
【図６】正弦波信号（ＳＳ）を印加することによって得られたインピーダンスとホワイトノイズ（ＢＢ）を使用して得られたインピーダンスとの比較の図である。
【図７】電池の計測された容量に対してＳｏＨを予測するための関係に基づいて電池の容量について計算されている直線ａ）と、複数のインピーダンス線図から計算された容量と電池の計測された容量との差を表している残差ｂ）の図である。
【図８】２０℃での検査（ＣＫ）中の完全なサイクルとエージング（ＶＩ）中の５０℃でのインピーダンスとによる電池の容量の測定値を示している図である。
【図９】ＳｏＣの計測値に対してＳｏＣを予測するための関係に基づいて電池のＳｏＣについて計算されている直線ａ）と、複数のインピーダンスの図から計算されたＳｏＣと計測されたＳｏＣとの間の差を表している残差ｂ）との図である。
【発明を実施するための形態】
【００３１】
本発明の方法は、予め特定されているモデルを使用して電池の充電状態または劣化状態の計測器を作るために使用可能で、また運輸用途（牽引電池）においてまたは再生可能なエネルギーを蓄えるために、その計測器の使用のための技術とすることができる。提案されている原理は、ＢＭＳによるＳｏＣとＳｏＨとの予測値を統合するが、これらのデータは直接計測可能ではない。
【００３２】
本発明は、潜在的に乗り物において搭載状態であるか、送電網に接続されている光起電力太陽装置の場合のエネルギーの蓄え時に使用され、装置の複数の電極の複数の端子の位置での電気インピーダンスの計測に基づいて電池、特にリチウムイオン電池の充電状態（ＳｏＣ）と劣化状態（ＳｏＨ）を定量的に求めるために使用することが可能で、計測は非侵入的で制御された温度で行われる。
【００３３】
本発明の論理線図を図１に示す。本発明の方法は以下のステップを有している。
【００３４】
ステップＥ１：ＳｏＣ、ＳｏＨ、およびＴの関数として複数のインピーダンス（Ｚ）の線図を計測するために、一群の電池（Ｂａｔ．）について複数の実験室用試験のプログラムを実行する。
【００３５】
ステップＥ２：ＳｏＣ、ＳｏＨおよびＴの関数である一式のパラメータ（ｐａｒａ．）を求めるために、計測されたインピーダンス線図（Ｚ）を使用して、選択されたモデル（ＲＣ回路）（ｍｏｄ．）を調整する。
【００３６】
ステップＥ３：複数のパラメータの多変量の組み合わせから量ＳｏＣとＳｏＨとを計算する。ＳｏＣの計算のための関係とＳｏＨの計算のための関係の少なくとも一方が得られる（Ｒｅｌ．１およびＲｅｌ．２）。
【００３７】
ステップＥ４：選択されたモデルと計算された複数の関係は、電気信号を調査中の電池（ＢａｔＥ．）に印加してインピーダンスＺを計測する計器（ＩＭＩ）と、計測されたインピーダンスＺへの選択されたモデル（ｍｏｄ．）の調整のためと、それから、得られたパラメータ（ｐａｒａ．）と以前に計算された関係からＳｏＣとＳｏＨの少なくとも一方の計算（ＣＡＬＣ）のためのソフトウェア部分（ＬＯＧ）とから構成されている計測器（Ｇ）において使用される。
【００３８】
（１  ＳｏＣ、ＳｏＨの関数としての電気化学インピーダンス線図の計測）
ＳｏＣ、ＳｏＨ、および場合によっては温度の関数として複数の電気化学インピーダンス線図を記録するために実験室用試験のプログラムが実施される。一般に、第１の電気化学装置と同じ種類の少なくとも１つの第２の電気化学装置のさまざまな内部状態について、第２の装置の内部状態に関連している特性が計測され（ＳｏＣ、ＳｏＨ）、第２の電気化学装置の電気応答が様々な周波数で計測される。
【００３９】
一実施形態において、所与の種類の電池（ＢａｔＥ）とその電池の所与の用途について、同じ種類の電池が使用される（Ｂａｔ．）。次に、複数の電気応答の計測がその電池の様々な充電状態と劣化状態について行われる。その電池の様々な劣化状態を得るために、考えられる用途を表している加速されたエージングを実施してもよい。例として、実験室において、電池はハイブリッド乗り物の種類の搭載用途をシミュレーションする加速されたエージングのプロトコル、または送電網に接続されている光起電力由来のエネルギーを蓄える用途をシミュレーションする加速されたエージングのプロトコルの対象とされる。
【００４０】
ガルバノスタットを用いて正弦の外乱、好ましくは電流の外乱を電池に加え、複数の端子の位置に誘起された正弦波電圧を計測することよって複数のインピーダンス線図を計測してもよい。他の実施形態において、外乱は、単純な正弦波外乱の形態ではなく、複数の正弦曲線の重ね合わせの形態、またはホワイトノイズ（全ての周波数が同じ信号内で重ね合わされている）の形態で加えられてさえもよく、これは５〜６個のまたは全ての周波数応答をそれから同時に分析できることを意味している。
【００４１】
ＳｏＣの関数としての複数のインピーダンス線図の計測は、ＳｏＣの全範囲または用途で使用される範囲に対応しているＳｏＣの範囲にわたって行われてもよい。
【００４２】
用途の動作の温度範囲にわたる温度を備えた複数のインピーダンス線図の変形例も計測する。
【００４３】
変化の各状態とエージングの状態との少なくとも一方で、ガルバノスタットを使用して電流外乱を加えることによって電気化学装置の電気インピーダンスＺが計測される。
【００４４】
複素量Ｚ（実部ＲｅＺと虚部ＩｍＺを備えている）は、ナイキスト線図の形態で表してもよく、Ｉｍ（Ｚ）はＲｅＺの関数であって、各点は１つの周波数に対応している。そのような線図を図３に示している。それから、高速の現象（高周波数での内部抵抗）に対する応答を、複数の電極での反応などの中速の現象や、低速の現象（ワールブルクインピーダンスによって示される低周波数での媒体中のイオンの拡散）から区別することができる。これらのさまざまな現象は、ＳｏＣとＳｏＨに対して様々な程度で影響を受けやすい。したがって、インピーダンス応答は、充電とエージングの状態の関数として変化し、これらの効果を分離するのが困難である。
【００４５】
対象の電池と同じ種類の第２の電池の使用を説明しており、同じ種類の一式の電池の使用も可能であって、これらの電池の各々は充電状態と劣化状態の少なくとも一方が異なっている。
【００４６】
（２  等価電気回路によるインピーダンス線図のモデリング）
全ての状態（ＳｏＣ、ＳｏＨ、および温度）について得られたナイキスト線図が、好ましくは、複数の抵抗器と複数のコンデンサがＳｏＣとＳｏＨに依存しているが、単純な比例関係で依存しているということを知って、等価電気回路（直列と並列の少なくとも一方の複数の抵抗器と複数のコンデンサの構成）に基づいてモデル化される。
【００４７】
図４は、電気化学蓄電池を表している等価電気回路の例を示している。Ｒ０は、高周波抵抗または素子の直列抵抗を表しており、Ｒ１は電荷移動抵抗を表しており、Ｑ１は電気化学２重層現象を表している一定位相素子を表しており、Ｗは分散現象を表しているワールブルクインピーダンスを表している。
【００４８】
等価回路は、構成要素の数を限定し、物理的な意味をできるだけ維持しながら、電池の全ての状態について装置のインピーダンスの最良のモデルを実現するように選択される。
【００４９】
選択されたモデル（ｍｏｄ．）は、モデルの複数のパラメータを変化させることによって、電池のＳｏＣ、ＳｏＨ、および温度（Ｔ）の各々に対応している試験プログラムの各インピーダンス線図に対して調整される。幾何学的なアプローチによるモデリングは粗いが、より高速である（たとえば、半円の直径と、線形な低周波拡散部分の勾配を求める場合）。
【００５０】
両方共、複数のモデルを記述する複数の量は、ＳｏＣ、ＳｏＨ、および温度に左右される。
【００５１】
（３  ＳｏＣまたはＳｏＨとモデルパラメータとの間の関係の決定）
このステップの間、特性（ＳｏＣ、ＳｏＨ）と等価電気回路（モデル）のパラメータとの間の関係が、各内部状態について得られた特性およびパラメータの複数の値の統計解析を使用して較正される。
【００５２】
多変量の組み合わせの形式の方程式が、ＳｏＣまたはＳｏＨと、考慮しているモデルの記述との間で求められる。
【００５３】
その結果、一方のＳｏＣまたはＳｏＨと他方のモデルの複数のパラメータ（および場合によっては電池の温度と電圧）との間で多変量解析が実行される。したがって、ＳｏＣおよびＳｏＨの少なくとも一方は、互いに独立している複数の等価電気回路の様々なパラメータの値の変化のみから予測されない。それに対して、本発明によれば、これらのパラメータ全ての組み合わせに依存している法則が、この多変量解析を使用して定められる。これは、ＳｏＣまたはＳｏＨの最適な予測を保証する最適な多変量の法則を求められることを意味している。
【００５４】
例として、複数の電気パラメータの主成分分析の形式の処理を使用してもよい。
【００５５】
電気モデルの例として、ＳｏＣと抵抗、ワールブルク係数との間に以下の関係が得られる。
【００５６】
ＳｏＨ＝ａ×Ｃ１＋ｂ×ａｌ１２＋ｃ×Ｗ＋ｄ×Ｗ２＋ｅ×Ｒ１＋ｆ×Ｌ０．Ｒ０＋ｇ×Ｃ１Ｒ１
有用であることが証明された場合、温度をパラメータとしてモデルパラメータに追加してもよい。同様に、電池の電圧をパラメータとして複数のモデルパラメータに追加してもよい。したがって、この関係は、パラメータＴ、ＳｏＣ、およびＳｏＨを制御することによって、またモデルの複数のパラメータの主成分分析などの数学的な処理によって、選択された電池の種類および考えられる用途についての実験室用試験のプログラムの結果として確立される。
【００５７】
（４  前記関係を使用した電気化学装置の内部状態の予測）
検討中の電気化学装置の電気応答が様々な周波数について求められる。この応答は、等価電気回路の電気応答が所定の電気応答と等しくなるパラメータを求めることで等価電気回路を使ってモデル化される。次に、電気化学装置の内部状態が、前記関係を使用して電気化学装置の内部状態に関する特性を計算することによって予測される。
【００５８】
実際には、前のステップで得られた関係が、ステップ１で説明した任意の方法を使用するインピーダンス計測装置（ＩＭＩ）によって構成されているセンサ（Ｇ）において一方で使用され、ソフトウェア部分として他方で使用される。ソフトウェア部分は、
計測装置ａ）がインピーダンスを有していない場合、インピーダンスを計算し、
計測されたインピーダンスにステップ２で選択されたモデルを自動的に当てはめ（ＬＯＧ）、
ステップ３で求められ、かつ以前に調整されている複数のモデルの複数のパラメータに基づいた、場合によっては電池の温度と電圧とに関連付けられた複数の関係を用いてＳｏＣまたはＳｏＨを計算する（ＣＡＬＣ）。
【００５９】
［例］
例として、本発明の方法のステップは、材料の対が異なる２つの電池（リチウムイオン蓄電池）、例えば、
リン酸鉄リチウム（ＬｉＦｅＰＯ４）の正電極への使用とチタン酸リチウム（Ｌｉ４Ｔｉ５Ｏ１２）の負電極への使用に基づいており、使用が開始されつつある対を使用したプロトタイプの蓄電池、
リン酸鉄リチウム（ＬｉＦｅＰＯ４）の正電極への使用と黒鉛、Ｃ６の負電極への使用に基づいている、より従来型の対を使用している商用蓄電池、
に適用される。
【００６０】
（実験室で実施された加速されたエージング）
場合によって、電池はハイブリッド乗り物の種類の搭載用途をシミュレーションする加速されたエージングのプロトコル、または送電網に接続されている光起電力由来のエネルギーを蓄える用途をシミュレーションする加速されたエージングのプロトコルの対象とされた。
【００６１】
（インピーダンス計測手続き）
ＳｏＣとＳｏＨの両方の電池診断のための方法の有効性を確認するために、「検査」試験手続きが定義された。この手続きは、エージングの前後で、通常４週間ごとに、雰囲気温度で電池を特徴付けるために使用された。
【００６２】
この試験は、図２に示しているように、４つの連続しているサイクルで構成された。この図において、サイクル番号は、ＮＣｙが前に付いている数字によって示されており、複数の曲線が充電の状態を表している。
【００６３】
第１のサイクル（ＮＣｙ１）は、電池が完全に充電されていることを保証するため、残りの放電とのそれに続く完全な充電とから構成された。第２のサイクル（ＮＣｙ２）は、容量の喪失、したがって電池の劣化状態を評価する試験であった。この試験は、次の２つのサイクル中の充電−放電電流を適合させるためにも使用することができる。インピーダンス計測に対する充電状態の影響は、最後の２つのサイクル中に一連の計測によって検討された。第３のサイクル（ＮＣｙ３）の目的は、休止期間の後で定電位型のインピーダンス分光法（図でＳＩＰと表されている）を使用することであって、第４のサイクル（ＮＣｙ４）の目的は、電流を中断することなくインピーダンスを計測することであったが、これは、充電及び放電段階の間の定電流モード（図でＳＩＧと表されている）でのインピーダンスの計測を意味する。しかし、定電位モードは、電圧調整充電の最後に必然的に使用される。
【００６４】
（インピーダンスと調整）
エージングの様々な程度について得られた複数のインピーダンスは、電池の総インピーダンスに対するエージングの様々な影響を観察するために、同じナイキスト線図（図３の例）上で表される。図３は、Ｌｉ４Ｔｉ５Ｏ１２／ＬｉＦｅＰＯ４型の電池についての２０％の充電状態でのＶＥＨ用途を表している５〜６個のエージングの状態について得られたインピーダンスの比較を表している。Ｖ１：初期状態、Ｖ２：２週間のエージング後、Ｖ３：４週間後、Ｖ４：６週間後、Ｖ５：８週間後。
【００６５】
注目すべきことは、インピーダンスは重ならず、エージングされた電池のインピーダンスは、より大きな半径の半円となることがわかるであろう。これらの差は、本発明によって、Ｒ０＋Ｒ１／Ｑ１＋Ｗ（図４）の電気モデルに対してインピーダンスを調整することによって定量化され、調整の品質をナイキスト線図またはボード線図（図５）において試験することができる。図５は、ナイキスト表現ａ）とボード表現ｂ）における、そして図４の等価回路モデルを使用した２０％ＳｏＣでのＬｉ４Ｔｉ５Ｏ１２／ＬｉＦｅＰＯ４型の電池についての６５ｋＨｚと０．１Ｈｚとの間でのインピーダンスのモデルついての調整の例の図である。（ＥＸ）：実験計測と（ＭＡ）：調整済みのモデル。
【００６６】
調整によって得られた複数の抵抗を比較すると、それらはエージングの関数として増加することが観察されることが注目される。この例では、この特性は、インピーダンスが計測された温度に関わらず同一である。したがって、電池のエージングの電気回路の複数の構成要素の値に対する影響は示され、電池の劣化状態（ＳｏＨ）を計測するために使用される。
【００６７】
前述のインピーダンスは、周波数の異なる連続している複数の正弦信号を使用することによって得られた。インピーダンスは、ホワイトノイズを電池の充電／放電信号に重ね合わせるなど、異なる態様で得ることもできる。図６は、従来の経路（正弦波信号（ＳＳ））によって計測されたインピーダンスだけでなく、ホワイトノイズ（ＢＢ）によって計測されたインピーダンスを示している。ホワイトノイズを使用してはるかに多くの点が得られていることに注目すべきであり、これは調整がより正確にできることを意味している。
【００６８】
（電池の劣化状態（ＳｏＨ）の決定）
これらの試験は、Ｌｉ４Ｔｉ５Ｏ１２／ＬｉＦｅＰＯ４型のプロトタイプにおいて実施された。実験のプロトコルは、容量試験と５〜６個のインピーダンス試験（図２）とが実施された、２５℃での複数の検査期間を含む、５０℃でのハイブリッド乗り物型のエージングに基づいていた。
【００６９】
プロトタイプの総容量は、検査手続きの各ステップについて既知であった。
【００７０】
本発明によれば、インピーダンスは、抵抗、コンデンサ（つまり一定位相素子、ＣＰＥ）およびワールブルク素子（例えば図４）などの単純な電気要素によって構成されている非線形モデルを使用して調整された。これらの実験中に、複数のエージング期間の間、複数のインピーダンス計測も実施された。
【００７１】
我々は、電池の異なる充電状態と５つの異なる劣化状態とで実施された２９個のインピーダンス計測を得ている。また、これらの計測値の各々について、単純な等価電気回路モデルに対して、Ｒ０（直列抵抗）、Ｒ１（電荷移動抵抗）、Ｃ１（Ｃ１、ＣＰＥの量Ｑ１）、ａｌ１（ａｌ１、ＣＰＥの指数部）、およびＷ（ワールブルクインピーダンス）の５つの因子を調整することができる。
【００７２】
これらの因子の統計的な処理（重回帰分析）後、以下の関係が保持された。
【００７３】
ＳｏＨ＝ａ＋ｂ×Ｒ１＋ｃ×Ｒ１２＋ｄ×Ｒ１３＋ｅ×ａｌ１＋ｆ×Ｗ＋ｇ×Ｒ１＋ｈ×Ｗ２
ここで、Ｒ１、ａｌ１およびＷは、調整された電気インピーダンスパラメータを表している。
【００７４】
この関係を使用したエージング中の複数の予測された容量の変化が図７のプロトタイプの既知の複数の容量と比較された。
【００７５】
図７は、上のグラフａ）において、ＳｏＨの予測関係から予測された容量（Ｑｃａｌｃ）が縦軸上にあって、電池の計測された容量（Ｑ）が横軸上にあることを表している。直線は線形回帰に該当しており、下のグラフｂ）上では、複数のインピーダンス線図から計算された容量の間の差を表している残差（ΔＱｃａｌｃ）は縦軸上にあって、計測された電池の容量（Ｑ）は横軸に沿っている。
【００７６】
統計的には、グラフａ）上の複数の点は、モデルの信頼性を確認するためには、ｙ＝ｘの種類の直線に近くなければならない。残差については、（本明細書の場合のように）無作為の分散とならなければならない。非無作為の分散は、関係に妥当性がないことを示している。
【００７７】
計算された値と実際の値との間の差を表している残差の変化は、予測が正確に動作していることを示している。また、モデルによって説明される分散を示している相関係数Ｒ２は、０．９９９９に等しかった（Ｒ２＝０の場合、相関がなく、Ｒ２＝１の場合、完全な相関が存在する）。モデルによる標準誤差は０．２５％であって、これは非常に低い値であって、モデルの正確さを表している。分散分析も、調整因子がすべてモデルを表していることを示した。コルモロゴフ−スミルノフ検定も実施された。この検定は、モデルによって計算された複数の値と計測された複数の値とが同じ法則に従うことを確認することができる。この検定は、０．９５のＰ値（Ｐ＝０ならば、法則が異なる、Ｐ＝１ならば、法則は同じである）となり、これはモデルに対しては非常によい。したがって、実施された統計検定の観点から、得られたモデルは有効である。
【００７８】
（エージング中に計測されたインピーダンスの検証）
検査段階（図２）の間と同じ複数のインピーダンス計測がエージング期間中に実施され、以前に求められたモデルをこれらのインピーダンスに対して調整されている複数のパラメータの複数の値に適用することが可能であった。ただし、エージングインピーダンスは、検査段階と同じように２５℃ではなく５０℃で計測され、したがって、電池のインピーダンスはその温度に依存するので計測には歪みが存在した。
【００７９】
図８は、２つの方法、すなわち、２０℃での検査（ＣＫ）中のサイクルと、エージング（ＶＩ）中に計測された５０℃でのインピーダンスの計測と、によって求められた容量（Ｑ）を時（ｈｏｕｒ）で表した時間の関数として示している。時間の関数として２つの種類の結果の間の非常に強い類似性を示している。温度の差からくる歪みは規則的であって、計測された容量の値よりも予測される容量の値は常に高かった。この結果は意味があり、それは電池の容量は常に温度と共に上昇するからである。エージング中のＳｏＨ（容量によって表されている）は、温度の違いにも関わらず、エージング中に計測された複数のインピーダンス線図に適用された方法によって予測することができる。実際に、使用する例において、温度パラメータは検討されていない。このパラメータの積分は、予測の精度を改善することができることを意味することになる。
【００８０】
（電池の充電状態（ＳｏＣ）の決定）
実験のプロトコルは、２．３Ａｈの容量を有し完全に充電された黒鉛／リン酸鉄リチウムの商用リチウムイオン電池の使用と、５％の充電状態のステップでのその電池の放電からなる。各充電状態において、電池は安定化させるために休止させ、それからインピーダンス計測が定電流モードで実施された。データ処理は、劣化状態を判断するために適用されたデータ処理に類似している。
【００８１】
統計的な処理の後、保持された関係は次のようになる。
【００８２】
ＳｏＣ＝ａ×Ｃ１＋ｂ×ａｌ１２＋ｃ×Ｗ＋ｄ×Ｗ２＋ｅ×Ｒ１＋ｆ×Ｌ０．Ｒ０十ｇ×Ｃ１Ｒ１
ここで、Ｒ０、Ｒ１、Ｃ１、ａｌ１、Ｌ０およびＷは前述のように調整された電気インピーダンスパラメータを表している（Ｒ０は電解質の抵抗、Ｒ１は伝達抵抗、Ｃ１はＣＰＥの量Ｑ１、ａｌ１はＣＰＥの指数部、Ｌ０は高周波数インダクタンス、Ｗはワールブルクインピーダンス）。
【００８３】
この関係を使用して予測されたＳｏＣの変化は、電池の既知のＳｏＣと比較された。
【００８４】
図９は、上のグラフａ）において、ＳｏＣ予測関係から予測されたＳｏＣ（ＳｏＣ  ｃａｌｃ）が縦軸上にあって、電池の計測されたＳｏＣ（ＳｏＣ）が横軸上にあることを表している。直線は線形回帰に該当しており、下のグラフｂ）上では、複数のインピーダンス線図から計算された値の間の差を表している残差（ΔＳｏＣ  ｃａｌｃ）は縦軸上にあって、電池の計測された充電状態（ＳｏＣ）は横軸に沿っている。
【００８５】
複数の残差の変化は、統計的にモデルが機能していることを示している。また、モデルによって説明される分散を示している相関係数Ｒ２は、０．９９７に等しかった（Ｒ２＝０の場合、相関がなく、Ｒ２＝１の場合、完全な相関が存在する）。モデルによる標準誤差は４％であり、この値は非常に小さく（電池のＳｏＣについて４％の不確かさ）、モデルの精度を表している。分散分析も、調整因子がすべてモデルを表していることを示した。コルモロゴフ−スミルノフ検定も実施された。この検定は、モデルによって計算された複数の値と計測された複数の値とが同じ法則に従うことを確認することができる。この検定は、１のＰ値（Ｐ＝０ならば、法則が異なる、Ｐ＝１ならば、法則は同じである）となり、これはモデルに対しては非常によい。したがって、実施された統計検定の観点から、得られたモデルは有効である。
【図１】
2013519893.tif 000003
【図２】
2013519893.tif 000004
【図３】
2013519893.tif 000005
【図４】
2013519893.tif 000006
【図５ａ】
2013519893.tif 000007
【図５ｂ】
2013519893.tif 000008
【図６】
2013519893.tif 000009
【図７】
2013519893.tif 000010
【図８】
2013519893.tif 000011
【図９】
2013519893.tif 000012
【国際調査報告】
2013519893.tif 000013
2013519893.tif 000014
2013519893.tif 000015
2013519893.tif 000016
2013519893.tif 000017
2013519893.tif 000018