(57)【要約】【課題】比較的簡易かつ短時間で過充電メモリー効果に起因する電圧低下を推定可能な電圧低下推定方法を提供する。【解決手段】メモリー効果を有する電池の電圧低下推定方法であって、トリクル充電された電池を放電させる第１ステップと、電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する第２ステップと、電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を過渡電圧波形にフィッティングさせ、パラメータＲ１’を導出する第３ステップと、パラメータＲ１’を用いて電池の高ＳｏＣ領域における電圧低下を推定する第４ステップと、を含むことを特徴とする。【選択図】図１１
2020128910.tif 000002
【特許請求の範囲】
【請求項１】
メモリー効果を有する電池の電圧低下推定方法であって、
トリクル充電された前記電池を放電させる第１ステップと、
前記電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する第２ステップと、
前記電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を前記過渡電圧波形にフィッティングさせ、前記パラメータＲ１’を導出する第３ステップと、
前記パラメータＲ１’を用いて前記電池の高ＳｏＣ領域における電圧低下を推定する第４ステップと、
を含むことを特徴とする電圧低下推定方法。
【請求項２】
前記第４ステップでは、前記電池の高ＳｏＣ領域における最低電圧ＶＬＯＷと前記パラメータＲ１’との関係式を用いて、前記第３ステップで導出した前記パラメータＲ１’から前記最低電圧ＶＬＯＷを推定することを特徴とする請求項１に記載の電圧低下推定方法。
【請求項３】
メモリー効果を有する電池の電圧低下を推定する電圧低下推定装置であって、
トリクル充電された前記電池を放電させる放電部と、
前記電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する電圧測定部と、
前記電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を前記過渡電圧波形にフィッティングさせ、前記パラメータＲ１’を導出する演算部と、
前記パラメータＲ１’を用いて前記電池の高ＳｏＣ領域における電圧低下を推定する電圧低下推定部と、
を備えることを特徴とする電圧低下推定装置。
【請求項４】
前記電池の高ＳｏＣ領域における最低電圧ＶＬＯＷと前記パラメータＲ１’との関係式を記憶した記億部をさらに備え、
前記電圧低下推定部は、前記関係式を参照して、前記演算部で導出した前記パラメータＲ１’から前記最低電圧ＶＬＯＷを推定することを特徴とする請求項３に記載の電圧低下推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、メモリー効果を有する電池の電圧低下を推定する電圧低下推定方法および電圧低下推定装置に関する。
【背景技術】
【０００２】
メモリー効果を有する電池としては、例えば、Ｎｉ−ＭＨ電池（ニッケル水素二次電池）が知られている。Ｎｉ−ＭＨ電池は、可燃性の低い水溶性電解液が用いられていることや反応物質の溶解・析出反応が無いことから、防災機器、医療用機器、自動車などの高い安全性が要求される機器に多く採用されている。中でも誘導灯、非常用照明、無停電電源（ＵＰＳ）などの防災機器では、Ｎｉ−ＭＨ電池はトリクル充電により常に完全充電で待機状態にある（例えば、特許文献１参照）。
【０００３】
Ｎｉ−ＭＨ電池は、過充電によるメモリー効果の影響により、放電時に期待する電圧が得られず、防災機器が動作しない場合もある。このため、従来から、防災機器搭載のＮｉ−ＭＨ電池に対して定電流放電試験を行い、Ｎｉ−ＭＨ電池の端子間電圧が下限電圧（例えば、１．０Ｖ）に達するか否かで、非常時において防災機器が動作するか否かを診断している。
【０００４】
しかしながら、放電試験中に予測不能な災害が発生することも考えられ、放電試験はリスクを伴う。また放電試験には長時間（例えば、３０分程度）を要するため、特に大規模施設では多大な労力が必要となるばかりでなく、保守費用がかさむ。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１１−１７１１６１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は上記事情に鑑みてなされたものであって、その課題とするところは、過充電メモリー効果に起因する電圧低下を瞬時かつ簡易に推定可能な電圧低下推定方法および電圧低下推定装置を提供することにある。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明に係る電圧低下推定方法は、
メモリー効果を有する電池の電圧低下推定方法であって、
トリクル充電された前記電池を放電させる第１ステップと、
前記電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する第２ステップと、
前記電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を前記過渡電圧波形にフィッティングさせ、前記パラメータＲ１’を導出する第３ステップと、
前記パラメータＲ１’を用いて前記電池の高ＳｏＣ領域における電圧低下を推定する第４ステップと、
を含むことを特徴とする。
【０００８】
上記電圧低下推定方法は、
前記第４ステップでは、前記電池の高ＳｏＣ領域における最低電圧ＶＬＯＷと前記パラメータＲ１’との関係式を用いて、前記第３ステップで導出した前記パラメータＲ１’から前記最低電圧ＶＬＯＷを推定することが好ましい。
【０００９】
また、上記課題を解決するために、本発明に係る電圧低下推定装置は、
メモリー効果を有する電池の電圧低下を推定する電圧低下推定装置であって、
トリクル充電された前記電池を放電させる放電部と、
前記電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する電圧測定部と、
前記電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を前記過渡電圧波形にフィッティングさせ、前記パラメータＲ１’を導出する演算部と、
前記パラメータＲ１’を用いて前記電池の高ＳｏＣ領域における電圧低下を推定する電圧低下推定部と、
を備えることを特徴とする。
【００１０】
上記電圧低下推定装置は、
前記電池の高ＳｏＣ領域における最低電圧ＶＬＯＷと前記パラメータＲ１’との関係式を記憶した記億部をさらに備え、
前記電圧低下推定部は、前記関係式を参照して、前記演算部で導出した前記パラメータＲ１’から前記最低電圧ＶＬＯＷを推定することが好ましい。
【発明の効果】
【００１１】
本発明によれば、過充電メモリー効果に起因する電圧低下を瞬時かつ簡易に推定可能な電圧低下推定方法および電圧低下推定装置を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】未使用のＮｉ−ＭＨ電池の内部インピーダンスベクトル軌跡を示す図である。
【図２】Ｎｉ−ＭＨ電池の放電電圧波形のサイクル劣化特性を示す図である。
【図３】Ｎｉ−ＭＨ電池の内部インピーダンスのサイクル劣化特性を示す図である。
【図４】Ｎｉ−ＭＨ電池の放電電圧波形の過充電特性を示す図である。
【図５】過充電２００時間毎に測定したＮｉ−ＭＨ電池の内部インピーダンスベクトル軌跡を示す図である。
【図６】図４の各Ｎｉ−ＭＨ電池の内部インピーダンスベクトル軌跡を示す図である。
【図７】リコンディショニングによるＮｉ−ＭＨ電池の内部インピーダンスの変化を示す図である。
【図８】図７の各Ｎｉ−ＭＨ電池の放電電圧波形を示す図である。
【図９】Ｎｉ−ＭＨ電池の等価回路モデルである。
【図１０】（Ａ）は未使用のＮｉ−ＭＨ電池の過渡電圧波形に対する従来モデルのフィッティング結果を示す図である。（Ｂ）は過充電１０００時間のＮｉ−ＭＨ電池の過渡電圧波形に対する従来モデルのフィッティング結果を示す図である。
【図１１】（Ａ）は未使用のＮｉ−ＭＨ電池の過渡電圧波形に対する提案モデルのフィッティング結果を示す図である。（Ｂ）は過充電１０００時間のＮｉ−ＭＨ電池の過渡電圧波形に対する提案モデルのフィッティング結果を示す図である。
【図１２】ＫＷＷ指数αを変化させた際の内部インピーダンスベクトル軌跡の円弧成分の変化を示す図である。
【図１３】パラメータＲ１’およびＫＷＷ指数αのサイクル劣化依存性を示す図である。
【図１４】パラメータＲ１’およびＫＷＷ指数αの過充電時間依存性を示す図である。
【図１５】高ＳｏＣ領域における最低電圧ＶＬＯＷとパラメータＲ１’との関係を示す図である。
【図１６】本発明の一実施形態に係る電圧低下推定装置のブロック図である。
【発明を実施するための形態】
【００１３】
以下、添付図面を参照して、本発明に係る電圧低下推定方法および電圧低下推定装置の実施形態について説明する。また、メモリー効果を有する電池として、Ｎｉ−ＭＨ電池（ニッケル水素二次電池）を例に挙げて説明する。
【００１４】
［電圧低下推定方法］
本実施形態では、メモリー効果を「リコンディショニングにより回復可能である一時的な内部特性変化であり、放電電圧が低下する特性」と定義し、過充電により生じるメモリー効果を過充電メモリー効果と呼ぶ。
【００１５】
本実施形態では、公称電圧１．２Ｖ、公称容量３７００ｍＡｈのＮｉ−ＭＨ電池を用いる。また本実施形態において、リアクタンスの正負を逆転させて複素平面に描画した内部インピーダンスベクトル軌跡（Cole-Cole plot）は、特に記載がなければＳｏＣ（State of Charge）１００％にて交流重畳法を用いて測定したものとする。
【００１６】
図１に、未使用（Unused）のＮｉ−ＭＨ電池の内部インピーダンスベクトル軌跡を示す。図１の内部インピーダンスベクトル軌跡は、５００Ｈｚにおける実軸との交点が示す内部抵抗成分（Ａ）と、１００ｍＨｚから５００Ｈｚ間の円弧成分（Ｂ、Ｃ）と、１００ｍＨｚ以下の右上がりに増加していく拡散インピーダンス成分（Ｄ）からなる。円弧成分は、正円ではなく、ＢとＣのそれぞれの領域の楕円が重なり合った形となる。内部抵抗成分は、活物質間の接触抵抗および電池セパレータ内の電解液抵抗を示す。
【００１７】
図２に、Ｎｉ−ＭＨ電池の放電電圧波形のサイクル劣化特性を示す。Ｎｉ−ＭＨ電池をＳｏＣ１００％からＳｏＣ０％まで放電し、その後ＳｏＣ１００％まで充電することを劣化試験の１サイクルとする。また、定電流充電時の最大許容電流が２Ａ（０．５４Ｃ）であることから、劣化試験時の充電・放電は定電流０．５Ｃで実施する。
【００１８】
未使用のＮｉ−ＭＨ電池と、上記劣化試験を１００サイクル、２００サイクル、３００サイクル、４００サイクル行ったＮｉ−ＭＨ電池とを準備し、これらのＮｉ−ＭＨ電池に対して、ＳｏＣ１００％から端子間電圧が０．９Ｖに到達するまで、電流０．１Ｃで定電流放電を行う。その結果得られた電圧波形が、図２に示す放電電圧波形である。
【００１９】
図２より、Ｎｉ−ＭＨ電池は、サイクル劣化によっては高ＳｏＣ領域で放電電圧が電池搭載機器（例えば、防災機器）を使用不能とする下限電圧（例えば、１．０Ｖ）まで低下しないことが分かる。なお、本発明において高ＳｏＣ領域とは、ＳｏＣ８０％以上の領域、好ましくはＳｏＣ９０％以上の領域である。
【００２０】
図３に、図２で用いた各Ｎｉ−ＭＨ電池の内部インピーダンスのサイクル劣化特性を示す。図３に示すとおり、内部インピーダンスの内部抵抗はサイクル劣化の進行に伴い一様に大となる一方、内部インピーダンスの円弧高さ（虚部最大値）はわずかに小となった後に再び大となる。
【００２１】
すなわち、サイクル劣化の進行に伴い、実軸切片が示す内部抵抗は一様に大となることから、Ｎｉ−ＭＨ電池の劣化は内部抵抗と強い相関性を持つことが分かる。一方、サイクル劣化の進行に伴う円弧の変化は一様で無いことから、内部インピーダンスの円弧は劣化との顕著な相関性を有しないことが分かる。
【００２２】
Ｎｉ−ＭＨ電池は、使用環境・使用条件によってはメモリー効果が生じ、電池搭載機器の動作不良を引き起こす。例えば、Ｎｉ−ＭＨ電池は、トリクル充電による過充電状態におくと、過充電メモリー効果を発生させ、電池搭載機器の動作時間を大幅に減少させる。
【００２３】
Ｎｉ−ＭＨ電池のトリクル充電による過充電状態を模擬するために、Ｎｉ−ＭＨ電池を電流０．２Ｃにて連続充電する過充電加速試験を実施する。過充電加速試験における過充電時間は、ＳｏＣ１００％到達後からの充電時間とする。
【００２４】
図４に、Ｎｉ−ＭＨ電池の放電電圧波形の過充電特性を示す。未使用のＮｉ−ＭＨ電池と、過充電時間が１００時間、２５０時間、４５０時間、６５０時間、８００時間、１０００時間のＮｉ−ＭＨ電池とを準備し、これら７つのＮｉ−ＭＨ電池に対して、ＳｏＣ１００％から端子間電圧が０．９Ｖに到達するまで、電流０．１Ｃで定電流放電を行う。その結果得られたものが、図４に示す放電電圧波形である。
【００２５】
図４より、過充電時間の増加に伴い高ＳｏＣ領域（ＳｏＣ９０％以上の領域）で放電電圧が大きく低下していることが分かる。特に、過充電時間が４５０時間、６５０時間、１０００時間のＮｉ−ＭＨ電池は、高ＳｏＣ領域での放電電圧が一般的なＮｉ−ＭＨ電池の下限電圧である１．０Ｖまで低下することから、電池搭載機器の動作時間は想定値より大幅に短くなる。
【００２６】
過充電時間が４５０時間、６５０時間、１０００時間のＮｉ−ＭＨ電池は、高ＳｏＣ領域で放電電圧が下限電圧に達しているが、過充電時間が８００時間のＮｉ−ＭＨ電池は、高ＳｏＣ領域で放電電圧が下限電圧に達していない。つまり、Ｎｉ−ＭＨ電池の個体差によって、高ＳｏＣ領域での放電電圧、すなわち過充電メモリー効果の蓄積量に差が生じる可能性が高い（従来の動作時間による電池交換は不適当である）。
【００２７】
図５に、あるＮｉ−ＭＨ電池において、過充電２００時間毎に内部インピーダンスを測定した際の内部インピーダンスベクトル軌跡を示す。図５に示すとおり、過充電時間の増加に伴い、実軸切片が示す内部インピーダンスの内部抵抗および内部インピーダンスの円弧高さ（虚部最大値）は大となる。これより、過充電時間は内部インピーダンスに影響を与えることが明らかである。
【００２８】
図６に、図４の各Ｎｉ−ＭＨ電池の内部インピーダンスベクトル軌跡を示す。放電電圧の低下が大であった過充電４５０時間、６５０時間、１０００時間のＮｉ−ＭＨ電池は、内部抵抗が約３００ｍΩ以上かつ、内部インピーダンスの円弧の虚部最大値が約４．５ｍΩ以上まで大となっている。一方、放電電圧の低下が中程度であった過充電８００時間のＮｉ−ＭＨ電池は、内部抵抗が約１２０ｍΩ程度かつ、内部インピーダンスの円弧の虚部最大値が約４．０ｍΩ程度となった。
【００２９】
すなわち、放電電圧の低下が大であると、内部インピーダンスの円弧の虚部最大値も大となる。換言すれば、図６を見ると、高ＳｏＣ領域で放電電圧の低下が大となるＮｉ−ＭＨ電池を推定することができる。
【００３０】
図７に、過充電１０００時間のＮｉ−ＭＨ電池、１度のリコンディショニングを行ったＮｉ−ＭＨ電池、５度のリコンディショニングを行ったＮｉ−ＭＨ電池、未使用のＮｉ−ＭＨ電池の内部インピーダンスベクトル軌跡を示す。リコンディショニングとは、Ｎｉ−ＭＨ電池に対してＳｏＣ０％まで放電後、再度ＳｏＣ１００％まで充電することをいう。図７に示すとおり、過充電時間の増加に伴い大となった内部インピーダンスの内部抵抗および円弧の虚部最大値は、リコンディショニングに伴い小となった。
【００３１】
図８に、図７の各Ｎｉ−ＭＨ電池に対して、ＳｏＣ１００％から端子間電圧が０．９Ｖに到達するまで電流０．１Ｃで定電流放電を行い、その結果得られた放電電圧波形を示す。図８に示すとおり、高ＳｏＣ領域での放電電圧の低下はリコンディショニングによって大きく改善されている。このことから、内部インピーダンスが大となり、高ＳｏＣ領域で放電電圧が低下した現象は、過充電メモリー効果によるものであることが分かる。
【００３２】
内部インピーダンスの内部抵抗の増加は、過充電メモリー効果と劣化のそれぞれの影響を受ける。一方、内部インピーダンスの円弧高さ（虚部最大値）の変化は、過充電メモリー効果のみの影響を受ける。このため、内部インピーダンスの円弧成分の測定により，過充電メモリー効果による放電電圧の低下が推定可能となる。
【００３３】
なお、一般的に内部抵抗の測定によりＮｉ−ＭＨ電池の劣化診断が可能とされていたが、過充電メモリー効果によっても内部抵抗が増加するため、Ｎｉ−ＭＨ電池の劣化診断はリコンディショニング後に実施することが好ましい。
【００３４】
ところで、Ｎｉ−ＭＨ電池の端子間電圧は、充電・放電終了後に急峻な変化を示し、その後、時間経過に伴い開放電圧に緩やかに収束するという電圧過渡特性を有する。また、Ｎｉ−ＭＨ電池の内部インピーダンスは、図１に示すとおり容量性を示すことから、図９に示すようにＲＣ並列回路をｎ段直列接続した等価回路モデルで一般的に模擬される。
【００３５】
ここで、直列内部抵抗Ｒ０は内部インピーダンスの内部抵抗成分に対応し、ｎ段直列接続したＲＣ並列回路は内部インピーダンスの円弧成分に対応する。また、ＩはＮｉ−ＭＨ電池の電流、Ｖ０はＮｉ−ＭＨ電池の内部電圧、ＶＺは内部インピーダンスによる電圧降下、ＶＢはＮｉ−ＭＨ電池の端子間電圧に対応する。
【００３６】
上記のとおり、放電電圧の低下が大であると内部インピーダンスの増加も大となり、内部インピーダンスの円弧の虚部最大値の増加は過充電メモリー効果のみの影響を受けるので、対応する電圧過渡特性のみ抽出すれば、放電電圧の低下が推定可能である。この抽出は、過渡電圧波形に電池モデルをフィッティングし、その電池モデルのパラメータを導出することにより実現する。なお、フィッティングに際しては、内部インピーダンスの内部抵抗に対応する直列内部抵抗Ｒ０は対象とせず急峻な過渡応答成分を無視し、内部インピーダンスの円弧に対応する成分のみを検討する。
【００３７】
従来から、ＲＣ並列回路１段からなる等価回路モデル（従来モデル）を用いた模擬法も広く利用されており、階段状電流に対する電圧過渡応答は下記（１）式により示される。下記（１）式において、τ１＝Ｒ１Ｃ１である。
【数１】
2020128910.tif 000003
【００３８】
図１０に、Ｎｉ−ＭＨ電池（（Ａ）は未使用の電池、（Ｂ）は過充電１０００時間の電池）を定電流０．２Ｃで５秒間のパルス放電した直後の過渡電圧波形に対して、従来モデルの（１）式を非線形最小二乗法によりフィッティングした結果を示す。なお、図１より内部インピーダンスの円弧成分の最大周波数が５００Ｈｚであることから、サンプリング定理よりサンプリング周波数を１ｋＨｚ、電圧分解能ΔＶは０．２ｍＶとした。
【００３９】
図１０（特に、図１０（Ｂ））に示すとおり、（１）式の従来モデルは、内部インピーダンスの過渡特性、すなわち周波数特性を観測するには模擬精度が不十分である。ＲＣ並列回路を２段以上とする多段回路モデルでは、より高精度に過渡応答を再現可能であるが、計算負荷が大となるので短時間かつ簡易な診断実現の視点から適切でない。
【００４０】
本実施形態では、過渡電圧波形立ち上がり時の時定数補正を可能とする提案モデル、すなわちＫＷＷ（Kohlausch-Williams-Watts）関数を含んだ下記（２）式（本発明の「電圧過渡応答表現式」に相当）を提案する。下記（２）式のαをＫＷＷ指数という。
【数２】
2020128910.tif 000004
【００４１】
（２）式は、第１項のパラメータＲ０が内部インピーダンスの内部抵抗成分に対応し、第２項のパラメータＲ１’およびＫＷＷ関数が内部インピーダンスの円弧成分に対応する。さらに、ＫＷＷ指数αは円弧高さに対応し、パラメータＲ１’は円弧直径および円弧高さに対応する。すなわち、内部インピーダンスの円弧高さ（虚部最大値）は、ＫＷＷ指数αとパラメータＲ１’の２つで定まる。
【００４２】
図１１に、Ｎｉ−ＭＨ電池（（Ａ）は未使用の電池、（Ｂ）は過充電１０００時間の電池）を定電流０．２Ｃで５秒間のパルス放電した直後の過渡電圧波形に対して、提案モデルの（２）式を非線形最小二乗法によりフィッティングした結果を示す。なお、図１０と同様に、サンプリング周波数を１ｋＨｚ、電圧分解能ΔＶを０．２ｍＶとした。
【００４３】
図１１に示すとおり、提案モデルの（２）式は、過充電時間によらず過渡特性を良好に表現する。提案モデルの（２）式によれば、従来モデルの（１）式に一つのパラメータαを追加するのみで、Ｎｉ−ＭＨ電池の電圧過渡特性を高精度に表すことができる。
【００４４】
（２）式を微分することにより得られるインパルス応答関数を級数展開し、ラプラス変換することにより解析的に周波数領域へ変換する。これにより得られる内部インピーダンスの円弧成分を表す（２）式第２項の周波数関数を、下記（３）式に示す。
【数３】
2020128910.tif 000005
【００４５】
図１２に、（２）式においてＫＷＷ指数αのみを０＜α≦１の範囲で変化させた際の内部インピーダンスベクトル軌跡の変化を示す。ただし、ＫＷＷ指数α以外のパラメータは、例えば、パラメータＲ１’＝５．３ｍΩ、τ１＝６７．８ｍｓとする。
【００４６】
パラメータＲ１’を固定した場合、ＫＷＷ指数αを小にすると、内部インピーダンスの円弧直径は変化せず、円弧高さのみ小となる。すなわち、（２）式に示す提案モデルのＫＷＷ指数αは、内部インピーダンスの円弧高さの歪みを表す。なお、ＫＷＷ指数αを固定した場合、パラメータＲ１’が小となれば、円弧高さも小となる。
【００４７】
したがって、パルス放電終了後の過渡電圧波形に対して提案モデル（２）式をフィッティングすることにより、内部インピーダンスの円弧成分を抽出することが可能となり、放電電圧低下の推定が可能となる。
【００４８】
本実施形態に係る電圧低下推定方法では、まずトリクル充電されたＮｉ−ＭＨ電池を放電させる（本発明の「第１ステップ」に相当）。例えば、ＳｏＣ１００％のＮｉ−ＭＨ電池を定電流０．２Ｃで５秒間のパルス放電を行う。
【００４９】
次いで、Ｎｉ−ＭＨ電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する（本発明の「第２ステップ」に相当）。測定時間は、過渡電圧波形を取得できる時間が確保されていれば、電池に応じて適宜変更できる。本実施形態では、例えば、測定時間を１秒とする。
【００５０】
次いで、内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式（例えば、（２）式）を、第２ステップで取得した過渡電圧波形にフィッティングさせて、少なくともパラメータＲ１’を含む電圧過渡応答表現式のパラメータを導出する（本発明の「第３ステップ」に相当）。フィッティングの手法としては、例えば、非線形最小二乗法を用いる。
【００５１】
図１３に、図３の各Ｎｉ−ＭＨ電池に関する、第３ステップで導出したパラメータＲ１’およびＫＷＷ指数αのサイクル劣化依存性を示す。図１３に示すとおり、サイクル劣化の進行に伴うパラメータＲ１’およびＫＷＷ指数αの変化はわずかである。
【００５２】
図１４に、図４の各Ｎｉ−ＭＨ電池に関する、第３ステップで導出したパラメータＲ１’およびＫＷＷ指数αの過充電時間依存性を示す。図１４に示すとおり、パラメータＲ１’は過充電時間の増加に伴い大となる傾向がある。過充電８００時間のＮｉ−ＭＨ電池におけるパラメータＲ１’のみ小となったが、これは図４に示す高ＳｏＣ領域の放電電圧の低下がわずかであったためである。なお、ＫＷＷ指数αの変化は、パラメータＲ１’の変化と比べて小である。
【００５３】
パラメータＲ１’は、サイクル劣化が進行した場合における変化がわずかである一方、過充電時間の増加に伴い大となる。このため、パラメータＲ１’は、過充電メモリー効果のみに依存することが明らかである。すなわち、内部インピーダンスの円弧高さは、パラメータＲ１’とＫＷＷ指数αの組み合わせで定まるが、過渡応答領域における円弧高さ（過充電メモリー効果に伴う円弧高さ）は、主にパラメータＲ１’で定まる。
【００５４】
図１５に、図４の各Ｎｉ−ＭＨ電池の高ＳｏＣ領域（ＳｏＣ９０％以上の領域）における最低電圧ＶＬＯＷとパラメータＲ１’との関係を示す。図１５に示すとおり、最低電圧ＶＬＯＷとパラメータＲ１’は強い相関性を有しており、その特性は下記（４）式で表すことができる。下記（４）式は、例えば、線形最小二乗法で求めることができる。
【数４】
2020128910.tif 000006
【００５５】
すなわち、本実施形態に係る電圧低下推定方法では、Ｎｉ−ＭＨ電池の高ＳｏＣ領域における最低電圧ＶＬＯＷとパラメータＲ１’との関係式（例えば、（４）式）を用いて、第３ステップで導出したパラメータＲ１’から最低電圧ＶＬＯＷを推定する（本発明の「第４ステップ」に相当）。
【００５６】
結局、本実施形態に係る電圧低下推定方法によれば、診断のための電流波形測定は不要かつＳｏＣ低下はわずかであり、数秒の過渡電圧波形の測定のみにより内部インピーダンスを抽出可能であるため、過充電メモリー効果に起因する電圧低下を瞬時かつ簡易に推定できる。その結果、予測不能な災害時のリスクを減らし、保守費の削減にも貢献することができる。
【００５７】
また、（２）式で表される提案モデルは、ＲＣ並列回路１段からなる等価回路モデル（従来モデル）の欠点である低模擬精度、多段回路モデルの計算負荷が大となる問題を解決し、高精度かつ簡易に電池特性を模擬可能にする。しかも、提案モデルは、防災設備に搭載されているＮｉ−ＭＨ電池をはじめ、トリクル充電により長期間待機状態にあるＮｉ−ＭＨ電池に適用可能である。
【００５８】
［電圧低下推定装置］
図１６に、本実施形態に係る電圧低下推定装置１を示す。電圧低下推定装置１は、放電部２と、電圧測定部３と、演算部４と、記億部５と、電圧低下推定部６とを備え、Ｎｉ−ＭＨ電池の高ＳｏＣ領域における電圧低下を推定する。
【００５９】
放電部２は、任意の間隔で電圧低下推定方法の第１ステップを実行し、トリクル充電されたＮｉ−ＭＨ電池を放電させる。
【００６０】
電圧測定部３は、電圧低下推定方法の第２ステップを実行し、Ｎｉ−ＭＨ電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する。
【００６１】
演算部４は、内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式（例えば、（２）式）を、電圧測定部３で取得した過渡電圧波形にフィッティングさせて、少なくともパラメータＲ１’を含む電圧過渡応答表現式のパラメータを導出する。すなわち、演算部４は、電圧低下推定方法の第３ステップを実行する。
【００６２】
記億部５には、Ｎｉ−ＭＨ電池の高ＳｏＣ領域における最低電圧ＶＬＯＷとパラメータＲ１’との関係式（例えば、（４）式）が予め記憶されている。
【００６３】
電圧低下推定部６は、記億部５の関係式を参照して、演算部４で導出したパラメータＲ１’から最低電圧ＶＬＯＷを推定する。すなわち、電圧低下推定部６は、電圧低下推定方法の第４ステップを実行する。
【００６４】
結局、本実施形態に係る電圧低下推定装置１は、任意の間隔で本実施形態に係る電圧低下推定を実行し、過充電メモリー効果に起因する電圧低下を瞬時かつ簡易に推定することができる。その結果、予測不能な災害時のリスクを減らし、保守費の削減にも貢献することができる。
【００６５】
以上、本発明に係る電圧低下推定方法および電圧低下推定装置の実施形態について説明したが、本発明は上記実施形態に限定されるものではない。
【００６６】
本発明に係る電圧低下推定方法は、トリクル充電された電池を放電させる第１ステップと、電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する第２ステップと、電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を過渡電圧波形にフィッティングさせ、パラメータＲ１’を導出する第３ステップと、パラメータＲ１’を用いて電池の高ＳｏＣ領域における電圧低下を推定する第４ステップとを含むのであれば、適宜構成を変更できる。
【００６７】
本発明における電圧低下推定装置は、トリクル充電された電池を放電させる放電部と、電池の放電直後の端子間電圧を測定して過渡電圧波形を取得する電圧測定部と、電池の内部インピーダンスベクトル軌跡の円弧成分に対応するパラメータＲ１’とＫＷＷ関数とを含む電圧過渡応答表現式を過渡電圧波形にフィッティングさせ、パラメータＲ１’を導出する演算部と、パラメータＲ１’を用いて電池の高ＳｏＣ領域における電圧低下を推定する電圧低下推定部とを備えるのであれば、適宜構成を変更できる。
【００６８】
本発明における電圧低下推定装置は、警報を発する警報部を備えてもよい。例えば、誘導灯等の防災機器に組み込まれた電圧低下推定装置が任意の間隔で電圧低下推定方法を実行する場合、警報部は、電圧低下推定部で推定された最低電圧ＶＬＯＷが所定値よりも低くなると、視覚的および／または聴覚的な手段で警報を発する。
【００６９】
同様に、本発明に係る電圧低下推定方法は、第４ステップで推定された最低電圧ＶＬＯＷが所定値よりも低くなると、視覚的および／または聴覚的な手段で警報を発する第５ステップを含んでもよい。
【００７０】
本発明における電圧低下推定装置は、一部の構成（例えば、放電部と電圧測定部）と他の構成（例えば、演算部、記憶部および電圧低下推定部）を、それぞれ別の場所に設けることができる。例えば、一部の構成（放電部と電圧測定部）を誘導灯等の防災機器に組み込み、他の構成（演算部、記憶部、電圧低下推定部）を別の場所（例えば、ＰＣ等の中央監視部）に設けることができる。この構成によれば、誘導灯等の防災機器を遠隔で制御、監視することができる。
【００７１】
さらに、少なくとも一部の構成（例えば、放電部と電圧測定部）を複数備える場合、例えば、当該一部の構成（放電部と電圧測定部）を複数の誘導灯等の防災機器に組み込むことで、複数の誘導灯等の防災機器を統括的に遠隔で制御、監視することができる。
【００７２】
上記実施形態では、メモリー効果を有する電池として、Ｎｉ−ＭＨ電池（ニッケル水素二次電池）を例に挙げて説明したが、本発明は、Ｎｉ−ＭＨ電池以外のメモリー効果を有する電池にも適用可能である。
【符号の説明】
【００７３】
１    電圧低下推定装置
２    放電部
３    電圧測定部
４    演算部
５    記億部
６    電圧低下推定部
【図１】
2020128910.tif 000007
【図２】
2020128910.tif 000008
【図３】
2020128910.tif 000009
【図４】
2020128910.tif 000010
【図５】
2020128910.tif 000011
【図６】
2020128910.tif 000012
【図７】
2020128910.tif 000013
【図８】
2020128910.tif 000014
【図９】
2020128910.tif 000015
【図１０】
2020128910.tif 000016
【図１１】
2020128910.tif 000017
【図１２】
2020128910.tif 000018
【図１３】
2020128910.tif 000019
【図１４】
2020128910.tif 000020
【図１５】
2020128910.tif 000021
【図１６】
2020128910.tif 000022