(57)【特許請求の範囲】
【請求項１】
電流供給源からバッテリに供給されるか、前記バッテリから負荷に供給される電流の量を検出する電流計と、
前記バッテリの両端電圧を測定する電圧計と、
充電制御信号に応答して、前記電流供給源から供給される電流を、一定の周波数を有するパルス状のパルス電流として前記バッテリに伝達するか、放電制御信号に応答して、前記バッテリから出力される電流を、一定の前記周波数を有する前記パルス電流として前記負荷に伝達する制御スイッチ部と、
一定の前記周波数を有する前記充電制御信号と前記放電制御信号を生成し、前記バッテリと等価のバッテリ等価モデルに前記電流計および前記電圧計から測定された前記バッテリの電流および電圧を適用し、前記バッテリ等価モデルを構成する素子のパラメータ値を予測するプロセッサとを備え、
前記バッテリ等価モデルは、
前記バッテリの内部抵抗と、
前記バッテリの開放電圧を示す内部キャパシタと、
前記バッテリの内部の電気化学的反応による非線形的な動作特性を示すダイナミック抵抗とダイナミックキャパシタの並列回路とが直列に接続され、
前記充電制御信号および前記放電制御信号は、前記パラメータ値である前記内部抵抗、前記内部キャパシタ、前記ダイナミック抵抗、および前記ダイナミックキャパシタ毎の予測のための他の一定の周波数を有し、
前記パルス電流の任意の１周期内の電流供給期間が前記パルス電流に対応して変化する前記バッテリの電圧が飽和状態になるまで活性化されるように設定されることを特徴とするバッテリパラメータ管理システム。
【請求項２】
前記制御スイッチ部は、
前記充電制御信号に応答して、前記電流供給源から前記バッテリに供給される電流をスイッチングする充電制御スイッチと、
前記放電制御信号に応答して、前記バッテリから前記負荷に供給される電流をスイッチングする放電制御スイッチとを備えることを特徴とする請求項１記載のバッテリパラメータ管理システム。
【請求項３】
前記プロセッサは、
前記電流供給源から前記バッテリに電流が供給される充電期間中の一部期間において、前記パラメータである前記内部抵抗、前記内部キャパシタ、前記ダイナミック抵抗、および前記ダイナミックキャパシタ毎の予測のための他の一定の前記周波数を有する前記パルス電流が前記電流供給源から前記バッテリに供給されるように前記充電制御信号を設定し、
前記バッテリから前記負荷に電流が供給される放電期間中の一部期間において、前記パラメータ値である前記内部抵抗、前記内部キャパシタ、前記ダイナミック抵抗、および前記ダイナミックキャパシタ毎の予測のための前記他の一定の周波数を有する前記パルス電流が前記バッテリから前記負荷に供給されるように前記放電制御信号を制御することを特徴とする請求項２記載のバッテリパラメータ管理システム。
【請求項４】
前記バッテリの温度を測定する温度計をさらに備え、
前記プロセッサは、前記バッテリの温度、前記バッテリの電流および電圧を用いて前記バッテリの電荷量を予測することを特徴とする請求項１記載のバッテリパラメータ管理システム。
【請求項５】
前記バッテリ等価モデルは、
前記ダイナミック抵抗と前記ダイナミックキャパシタの並列回路の一端子が前記電圧計の一端子に接続され、
前記バッテリの前記内部抵抗および前記内部キャパシタが前記並列回路の他の一端子および前記電圧計の他の一端子に直列に接続されたことを特徴とする請求項１記載のバッテリパラメータ管理システム。
【請求項６】
請求項１記載のバッテリパラメータ管理システムで適用するバッテリ等価モデルのバッテリパラメータ推定方法において、
前記バッテリを充電または放電するためのパルス電流供給期間を、前記バッテリ等価モデルのパラメータである内部抵抗、内部キャパシタ、ダイナミック抵抗、およびダイナミックキャパシタ毎の予測のための第１期間ないし第３期間に区分し、前記第１期間ないし前記第３期間に対応してそれぞれ異なる一定の周波数を有するパルス電流を供給するパルス電流供給ステップと、
前記第１期間において、前記バッテリの充電電流の変化量、前記バッテリの放電電流の変化量、および前記バッテリの電圧の変化量を用いて前記内部抵抗の抵抗値を推定する内部抵抗の抵抗値推定ステップと、
前記第２期間において、前記バッテリに供給された総電荷量、前記バッテリの放電電荷量、および前記バッテリの電圧の変化量を用いて前記内部キャパシタのキャパシタンスを推定する内部キャパシタのキャパシタンス推定ステップと、
前記第３期間において、前記バッテリの充電電流の変化量、前記バッテリの放電電流の変化量、前記バッテリの電圧の変化量、および前記ダイナミック抵抗と前記ダイナミックキャパシタの時定数を用いて前記ダイナミック抵抗の抵抗値および前記ダイナミックキャパシタのキャパシタンスを推定するダイナミック素子のパラメータ推定ステップとを含むことを特徴とするバッテリパラメータ推定方法。
【請求項７】
前記第１期間ないし前記第３期間は、前記パルス電流の周波数によって決定され、
前記第１期間において、前記バッテリのインピーダンスが前記内部抵抗の抵抗値と等しくなり、
前記第２期間において、前記バッテリのインピーダンスが前記内部キャパシタのキャパシタンスと周波数との積に反比例する値を有し、
前記第３期間において、前記バッテリのインピーダンスが前記ダイナミック抵抗の抵抗値と前記内部抵抗の抵抗値との和と等しい値を有することを特徴とする請求項６記載のバッテリパラメータ推定方法。
【請求項８】
前記内部抵抗の抵抗値推定ステップにおける前記内部抵抗の抵抗値Ｒｓは、数式
【数１】
0006317541.tif 000015
を用いて推定し、
ここで、ΔＶは、前記第１期間の前記バッテリの電圧の変化量であり、ΔＩは、前記第１期間の前記バッテリに供給される電流の変化量であることを特徴とする請求項７記載のバッテリパラメータ推定方法。
【請求項９】
前記内部キャパシタのキャパシタンス推定ステップにおける前記内部キャパシタのキャパシタンスＣｉは、数式
【数２】
0006317541.tif 000016
を用いて推定し、
ここで、ΔＱは、前記第２期間の予め設定された時間差Δｔ＝ｔ０−ｔ１の間に前記バッテリに供給された総電荷量であり、ｉは、前記バッテリに供給される電流であり、ΔＶＯＣは、前記総電荷量によって変化した前記バッテリの開放電圧の最大値の変化量であることを特徴とする請求項７記載のバッテリパラメータ推定方法。
【請求項１０】
前記ダイナミック素子のパラメータ推定ステップは、
前記バッテリに供給される電流の変化量および前記バッテリの電圧の変化量を用いて前記ダイナミック抵抗の抵抗値を推定するダイナミック抵抗の抵抗値推定ステップと、
前記ダイナミック抵抗の抵抗値と前記時定数を用いて前記ダイナミックキャパシタのキャパシタンスを推定するダイナミックキャパシタ推定ステップとを含むことを特徴とする請求項７記載のバッテリパラメータ推定方法。
【請求項１１】
前記ダイナミック抵抗の抵抗値推定ステップは、前記パルス電流の１周期中に前記電流の大きさが急激に増加して一定の電流値を維持する期間を含み、
前記ダイナミックキャパシタ推定ステップは、前記ダイナミック抵抗の抵抗値推定ステップの後、前記電流の大きさが再び急激に減少して一定の電流値を維持する期間を含み、
前記ダイナミックキャパシタ推定ステップの前記電流の大きさが再び急激に減少する前に、バッテリ電圧の大きさが飽和することを特徴とする請求項１０記載のバッテリパラメータ推定方法。
【請求項１２】
前記ダイナミック抵抗の抵抗値Ｒｄは、数式
【数３】
0006317541.tif 000017
を用いて推定し、
ここで、ｔ２は、バッテリに供給される電流の大きさが急に増加する時点であり、ｔ３は、前記バッテリ電圧の大きさが飽和する時点であり、ΔＩは、バッテリに供給される電流の変化量であることを特徴とする請求項１１記載のバッテリパラメータ推定方法。
【請求項１３】
前記ダイナミックキャパシタのキャパシタンスＣｄは、数式
【数４】
0006317541.tif 000018
を用いて推定し、
ここで、Ｒｄは、ダイナミック抵抗の抵抗値、ｔ４は、電流が急に減少する時点、ｔ５は、時定数の値が１の時点、ｔ６は、時定数の値が２の時点、そして、ｔ７は、時定数の値が３の時点であることを特徴とする請求項１２記載のバッテリパラメータ推定方法。
【請求項１４】
前記数式に含まれた各時点は、３つの数式
【数５】
0006317541.tif 000019
を満足し、
ここで、ΔＶｄは、前記ｔ４の時点の前記バッテリの電圧値であることを特徴とする請求項１３記載のバッテリパラメータ推定方法。
【請求項１５】
前記パルス電流供給ステップは、
前記パルス電流に対応して変化する前記バッテリの電圧が飽和状態になるまで、前記パルス電流の任意の１周期内での電流供給期間が活性化されているように設定されることを特徴とする請求項６記載のバッテリパラメータ推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、バッテリに関するものであり、特に、リチウムイオンのような二次電池の固有特性および動作特性に影響を与える要素を考慮したバッテリモデルを構成する素子のパラメータを推定するバッテリパラメータ管理システムおよびバッテリパラメータ推定方法に関するものである。
【背景技術】
【０００２】
ハイブリッド電気自動車の電気駆動モータの動力源として用いられるリチウムイオンバッテリのような二次電池は、過放電される場合には永久的にバッテリの特性を失い、過充電される場合には爆発および火災が発生するため、充電および放電に注意を払わなければならない。バッテリの安定した動作、すなわち、充電および放電のために、バッテリの現在状態、つまり、バッテリの残存電荷量（Ｓｔａｔｅ  ｏｆ  Ｃｈａｒｇｅ、ＳＯＣ）を推定することは必須である。
【０００３】
バッテリの状態を正確に予測するためには、まず、バッテリに適したバッテリ等価モデル（Ｅｑｕｉｖａｌｅｎｔ  Ｍｏｄｅｌ）を選定することが一般的である。バッテリの外部ではバッテリの両端の電圧とバッテリに流れる電流のみを測定することができるが、このように測定された電圧および電流を用いてバッテリ等価モデルに用いられた電気素子のパラメータを推定する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１２-１５９４１４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
バッテリ等価モデルの構造が過度に複雑に設定された場合、素子の電気的、そして物理的特性はもちろんのこと、素子間の電気的関連性によって、バッテリ等価モデルに用いられた電気素子のパラメータを推定するための計算式が複雑になることは当然である。バッテリ等価モデルの構造が複雑な場合、複雑な計算式を演算する機能ブロックが非常に多く要求されるため、バッテリ等価モデルを構成する素子のパラメータの推定をハードウェアで実現する時、ハードウェアの占める面積は相当なものになるはずである。素子のパラメータの推定をソフトウェアで実現する場合にも、高性能の高価なプロセッサが必要になる。
【０００６】
したがって、ハードウェアまたはソフトウェアで推定する両方の場合とも、消費する面積および費用の面において、１つの集積回路で実現することが容易でない問題があった。
【０００７】
本発明が解決しようとする技術的課題は、単純な構造を有するバッテリ等価モデルを構成する素子のパラメータを簡単に推定することができるバッテリパラメータ管理システムを提供することである。
【０００８】
本発明が解決しようとする他の技術的課題は、単純な構造を有するバッテリ等価モデルを構成する素子のパラメータを簡単に推定することができるバッテリパラメータ推定方法を提供することである。
【課題を解決するための手段】
【０００９】
上記の技術的課題を達成するための、本発明の一実施形態にかかるバッテリパラメータ管理システムは、電流計と、電圧計と、制御スイッチ部と、プロセッサとを備える。前記電流計は、電流供給源から前記バッテリに供給されるか、前記バッテリから負荷に供給される電流の量を検出する。前記電圧計は、前記バッテリの両端電圧を測定する。前記制御スイッチ部は、充電制御信号に応答して、前記電流供給源から供給される電流を前記バッテリにスイッチングするか、放電制御信号に応答して、前記バッテリから供給される電流を前記負荷にスイッチングする。前記プロセッサは、前記充電制御信号および前記放電制御信号を生成し、前記バッテリと等価のバッテリ等価モデルに前記電流計および前記電圧計から測定されたバッテリの電流および電圧を適用し、前記バッテリ等価モデルを構成する素子のパラメータ値を予測する。前記バッテリ等価モデルは、バッテリの内部抵抗と、バッテリの開放電圧を示す内部キャパシタと、バッテリ内部の電気化学的反応による非線形的な動作特性を示すダイナミック抵抗とダイナミックキャパシタの並列回路とが直列に接続される。
【００１０】
上記の他の技術的課題を達成するための、本発明の一実施形態にかかるバッテリパラメータ推定方法は、請求項１記載のバッテリパラメータ管理システムで適用され、パルス電流供給ステップと、内部抵抗の抵抗値推定ステップと、内部キャパシタのキャパシタンス推定ステップと、ダイナミック素子のパラメータ推定ステップとを含む。前記パルス電流供給ステップは、前記バッテリを充電または放電する時間期間中のパルス電流供給期間の間に一定の周波数を有するパルス電流を供給する。前記内部抵抗の抵抗値推定ステップは、前記パルス電流供給期間中の第１区間において、前記バッテリの充電電流の変化量、前記バッテリの放電電流の変化量、および前記第１区間でのバッテリ電圧の変化量を用いて前記内部抵抗の抵抗値を推定する。前記内部キャパシタのキャパシタンス推定ステップは、前記パルス電流供給期間中の第２区間において、前記バッテリに供給された総電荷量、前記バッテリの放電電荷量、および前記第２区間での前記バッテリ電圧の変化量を用いて前記内部キャパシタのキャパシタンスを推定する。前記ダイナミック素子のパラメータ推定ステップは、前記パルス電流供給期間中の第３区間において、前記バッテリの充電電流の変化量、前記バッテリの放電電流の変化量、前記バッテリ電圧の変化量、および前記ダイナミック抵抗と前記ダイナミックキャパシタの時定数を用いて前記ダイナミック抵抗の抵抗値および前記ダイナミックキャパシタのキャパシタンスを推定する。
【発明の効果】
【００１１】
本発明にかかるバッテリパラメータ管理システムおよびバッテリパラメータ推定方法は、単純な構造を有するバッテリ等価モデルを用い、等価モデルに含まれた素子のパラメータをリアルタイムで極めて簡単に推定することができるという利点がある。
【図面の簡単な説明】
【００１２】
【図１】本発明の一実施形態にかかるバッテリ管理システムを示す。
【図２】本発明の他の実施形態にかかるバッテリパラメータ推定方法に用いられたバッテリ等価モデルを示す。
【図３】バッテリに供給されるパルス電流の周波数に応じたバッテリのインピーダンスを示す。
【図４】本発明にかかるバッテリ管理システムに適用される充電および放電動作の一実施形態を示す。
【図５】内部抵抗の抵抗値推定ステップに適用される第１区間での電圧および電流波形を示す。
【図６】内部キャパシタのキャパシタンス推定ステップに適用される第２区間での電圧および電流波形を示す。
【図７】ダイナミック素子のパラメータ推定ステップに適用される第３区間での電圧および電流波形を示す。
【図８】第４期間での電圧波形を拡大したものである。
【図９】ダイナミック素子のパラメータ推定ステップに適したパルス電流の周波数を決定する過程を示す。
【図１０】本発明にかかるバッテリパラメータ推定方法を適用したコンピュータシミュレーションの結果を示す。
【図１１】本発明によって推定されたインピーダンス値を基に推定した開放電圧と内部キャパシタに降下した電圧との関係を示す。
【発明を実施するための形態】
【００１３】
本発明と本発明の動作上の利点および本発明の実施によって達成される目的を十分に理解するためには、本発明の例示的な実施形態を説明する添付図面および添付図面に記載された内容を参照しなければならない。
【００１４】
以下、添付した図面を参照して本発明の好ましい実施形態を説明することにより、本発明を詳細に説明する。各図面に提示された同一の参照符号は同一の部材を示す。
【００１５】
本発明は、自動車などに使用している実際のバッテリの残存電荷量を直接測定しなくても予測できるように、実際のバッテリに対応するバッテリ等価モデルを設定し、実際に測定可能な実際のバッテリに供給される電流および電圧を用いてバッテリ等価モデルを構成する電気素子のパラメータの値を予測することができるバッテリパラメータ推定方法および当該推定方法が実際に実行されるバッテリパラメータ管理システムを提案する。
【００１６】
図１は、本発明の一実施形態にかかるバッテリ管理システムを示す。
【００１７】
図１を参照すれば、バッテリ管理システム１００は、温度計１０１と、電流計１０２と、電圧計１０３と、バッテリ１１０と、センシング抵抗１２０と、制御スイッチ部１３０と、電流供給源１４０と、負荷１５０と、プロセッサ１６０と、ユーザインタフェース装置１７０とを備える。
【００１８】
温度計１０１は、バッテリ１１０の温度を測定し、電流計１０２は、バッテリ１１０に直列に接続されたセンシング抵抗１２０に並列に接続され、バッテリ１１０に入出力される電流の量を間接的に検出し、電圧計１０３は、バッテリ１１０の両端電圧、すなわち、開放電圧を測定する。
【００１９】
制御スイッチ部１３０は、プロセッサ１６０で生成された充電制御信号および放電制御信号に応答して、電流供給源１４０から供給される電流をバッテリ１１０に伝達するか、バッテリ１１０から出力される電流を負荷１５０に伝達する。制御スイッチ部１３０は、充電制御信号に応答して、電流供給源１４０から供給される電流をバッテリにスイッチングする充電制御スイッチ１３１と、放電制御信号に応答して、バッテリ１１０から出力される電流を負荷１５０にスイッチングする放電制御スイッチ１３２とを備える。
【００２０】
プロセッサ１６０は、第一として、充電制御信号および放電制御信号を生成し、制御スイッチ部１３０に含まれた２つのスイッチ１３１、１３２のスイッチング動作を制御する。
【００２１】
第二として、バッテリ１１０の温度、バッテリ１１０に供給されるかバッテリ１１０から出力される電流、およびバッテリ１１０の開放電圧をモニタし、モニタされたバッテリの温度、電流および開放電圧を用いてバッテリ等価モデルに適用された素子のパラメータを計算する。この時、バッテリ等価モデルに適用された素子のパラメータを計算するために、バッテリが充電される期間中の一部期間において、バッテリに供給される電流が一定の周波数を有するパルス状となるように充電制御信号を制御し、バッテリが放電される期間中の一部期間において、バッテリから負荷に供給される電流が一定の周波数を有するパルス状となるように前記放電制御信号を制御する。
【００２２】
パーソナルコンピュータのようなユーザインタフェース装置１７０は、プロセッサ１６０で果たす機能の一部を担当するものであり、必要に応じてシステムに追加して用いることができる。
【００２３】
図２は、本発明にかかるバッテリパラメータ推定方法に用いられたバッテリ等価モデルを示す。
【００２４】
図２を参照すれば、本発明にかかるバッテリパラメータ推定方法に用いられるバッテリ等価モデル１１０は、並列に接続されたダイナミック抵抗２０４およびダイナミックキャパシタ２０３と、内部抵抗２０２と、内部キャパシタ２０１とが直列に接続された構造を有する。ここで、内部抵抗２０２は、バッテリ固有の抵抗成分を代表し、内部キャパシタ２０１は、バッテリに負荷が接続されていない状態でのバッテリの開放電圧ＶＯＣを代表する。並列接続されたダイナミック抵抗２０４とダイナミックキャパシタ２０３は寄生成分によるもので、バッテリ内部の電気化学的反応による非線形的な動作特性を代表する。
【００２５】
本発明では、バッテリ等価モデル１１０を構成する素子のパラメータを推定する時、処理速度が速いため、高価なプロセッサを必要としないようにするために、図２に示されたような簡単なバッテリ等価モデルに適用されるバッテリパラメータ推定方法を提案する。
【００２６】
以下の説明において、バッテリの電圧は、電圧計１０３で測定された電圧を意味し、バッテリに供給される電流、すなわち、バッテリに流入または流出する電流は、電流計１０２で測定された電流を意味する。
【００２７】
図３は、バッテリに供給されるパルス電流の周波数に応じたバッテリのインピーダンスを示す。
【００２８】
図３を参照すれば、図１に示されたバッテリ管理システム１００のバッテリ１１０に供給されるパルス電流の周波数ｓが０から第１周波数ｆ１に到逹するまでのバッテリのインピーダンスＺ（ｓ）は数式１、そして、第１周波数ｆ１は数式２のようにそれぞれ表現することができる。
【００２９】
【数１】
0006317541.tif 000002
【００３０】
【数２】
0006317541.tif 000003
【００３１】
ここで、ｓは、パルス電流の周波数であり、Ｃｉは、内部キャパシタ２０１のキャパシタンスである。数式１によれば、バッテリ１１０に印加されるパルス電流の周波数が第１周波数ｆ１より低い場合、内部キャパシタＣｉが開放（ｏｐｅｎ）されるため、バッテリのインピーダンスＺ（ｓ）は、内部キャパシタＣｉと周波数ｓとの積に反比例し、つまり、周波数ｓが増加するほど減少する。数式２によれば、第１周波数ｆ１は、内部抵抗２０２の抵抗値Ｒｓおよびダイナミック抵抗２０４の抵抗値Ｒｄの和とダイナミックキャパシタ２０３のキャパシタンスＣｄとの積に反比例する。
【００３２】
第１周波数ｆ１と第２周波数ｆ２との間におけるバッテリのインピーダンスＺ（ｓ）は数式３、そして、第２周波数ｆ２は数式４のようにそれぞれ表現することができる。
【００３３】
【数３】
0006317541.tif 000004
【００３４】
【数４】
0006317541.tif 000005
【００３５】
数式３によれば、パルス電流が第１周波数ｆ１と第２周波数ｆ２との間の周波数を有する時には、内部キャパシタＣｉが短絡（ｓｈｏｒｔ）し、ダイナミックキャパシタＣｄは開放される。数式４によれば、第２周波数は、ダイナミック抵抗ＲｄとダイナミックキャパシタＣｄとの積に反比例する。
【００３６】
パルス電流の周波数が第３周波数ｆ３以上の場合、バッテリのインピーダンスは数式５、そして、第３周波数ｆ３は数式６のようにそれぞれ表現することができる。
【００３７】
【数５】
0006317541.tif 000006
【００３８】
【数６】
0006317541.tif 000007
【００３９】
数式５によれば、パルス電流の周波数が高周波である第３周波数ｆ３以上の時、内部キャパシタ２０１とダイナミックキャパシタ２０３とがすべて短絡する。すなわち、すぐ直前の電圧センシング値とそれに続く瞬間の電圧センシング値が一定の値以上に変化する場合、このような電圧の変化は、バッテリ等価モデルの内部抵抗２０２に起因するため、本実施形態にかかる方法を通じて内部抵抗２０２の抵抗値を推定することができる。キャパシタは急激な電流の変化には応答せずに短絡するため、この時のキャパシタのインピーダンスは０になる。
【００４０】
数式６によれば、第３周波数ｆ３は、内部抵抗２０２の抵抗値Ｒｓと内部キャパシタ２０１のキャパシタンスＣｉとの積に反比例する。
【００４１】
以下、図３に示されたパルス電流の周波数の変化に応じたバッテリのインピーダンスの変化に基づき、バッテリ等価モデルを構成する内部抵抗２０２の抵抗値Ｒｓ、内部キャパシタ２０１のキャパシタンスＣｉ、ダイナミック抵抗２０４の抵抗値Ｒｄ、およびダイナミックキャパシタ２０３のキャパシタンスＣｄを推定する方法について説明する。
【００４２】
以下のすべての説明は、バッテリを充電する一定の時間期間中の一部時間期間の間に一定の周波数を有するパルス電流をバッテリに供給することを前提とし、パルス電流の最適な周波数を決定する方法については後述する。
【００４３】
図４は、本発明にかかるバッテリ管理システムに適用される充電および放電動作の一実施形態を示す。
【００４４】
図４を参照すれば、充電制御信号を用いて充電制御スイッチ１３１の開閉（Ｓｗｉｔｃｈ  Ｏｎ／Ｏｆｆ）を調整することにより、バッテリ１１０を定電流（Ｃｏｎｓｔａｎｔ  Ｃｕｒｒｅｎｔ）で充電する定電流チャージ（ＣＣ充電：ＣＣ  Ｃｈａｒｇｅ）期間中の一部期間４１０、４２０において、電流供給源１４０からバッテリ１１０に供給される定電流をパルス状とすることができる。同様に、放電制御信号を用いて充電制御スイッチ１３１の開閉を調整することにより、バッテリ１１０から定電圧（Ｃｏｎｓｔａｎｔ  Ｖｏｌｔａｇｅ）放電される定電圧チャージ（ＣＶ充電：ＣＶ  Ｃｈａｒｇｅ）期間中の一部期間４３０において、バッテリ１１０から負荷１５０に供給する電流をパルス状とすることができる。以下の説明におけるパルス電流の波形は、図４に示された方式で生成されたものである。
【００４５】
図３、数式５および数式６を参照すれば、第１区間、すなわち、パルス電流の周波数が高周波である第３周波数ｆ３以上の区間では、内部キャパシタ２０１とダイナミックキャパシタ２０３とがすべて短絡し、第１区間でのインピーダンスＺ（ｓ）は内部抵抗２０２の抵抗値Ｒｓと等しい。
【００４６】
本発明では、第１区間に対応する第３周波数ｆ３以上の高周波のパルス電流が必要であるが、このために、充電制御スイッチ１３１の開閉速度を高くし、高周波のパルス電流を生成させることもでき、パルス電流の値が瞬間的に遷移する狭い時間期間を用いることもできる。高周波のパルス電流を生成してバッテリに提供する代わりに、低周波のパルス電流をバッテリに提供すれば、パルス電流の周波数に応じた消費電力を低減することができるという利点がある。
【００４７】
図５は、内部抵抗の抵抗値推定ステップに適用される第１区間での電圧および電流波形を示す。
【００４８】
図５を参照すれば、内部抵抗２０２の抵抗値Ｒｓは、下記の数式７に表現したように、パルス電流供給期間中の第１区間において、バッテリに供給される電流の変化量ΔＩおよびバッテリ電圧の変化量ΔＶを用いて推定する。
【００４９】
【数７】
0006317541.tif 000008
【００５０】
以下の説明では、低周波のパルス電流の遷移期間を用い、高周波のパルス電流の遷移期間でのバッテリの電圧および電流を用いて第１区間でのインピーダンスを予測する。
【００５１】
図５を参照すれば、本発明では、パルス電流の大きさが低い電流レベル（ｃｕｒｒｅｎｔ  ｌｅｖｅｌ）から高い電流レベルに遷移（ｔｒａｎｓｉｅｎｔ）する第１センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ１）でインピーダンスを予測するが、第１センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ１）では、バッテリの電圧も低い電圧レベルから高い電圧レベルに遷移する。第１センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ１）において、パルス電流の電流レベル差である電流の変化量ΔＩとバッテリ１１０の両端の電圧レベル差である電圧の変化量ΔＶを数式７に代入すると、内部抵抗２０２の抵抗値Ｒｓを容易に推定することができる。
【００５２】
図３、数式１および数式２を参照すれば、第２区間、すなわち、パルス電流の周波数が第１周波数ｆ１未満、すなわち、低周波の場合、内部キャパシタ２０１が開放（ｏｐｅｎ）されるため、バッテリ１１０の電圧はバッテリの開放電圧ＶＯＣと等しくなり、第２区間でのバッテリ１１０のインピーダンスＺ（ｓ）は、内部キャパシタ２０１のキャパシタンスＣｉと周波数ｓとの積に反比例する。第２区間では、パルス電流の変化に対して内部キャパシタ２０１は非常に遅い反応を示すため、以下で説明する第２センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ２）の間に電流の値を積分して総電荷量ΔＱを計算し、総電荷量を第２センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ２）での開放電圧の変化量ΔＶＯＣで除することにより、内部キャパシタ２０１のキャパシタンスＣｉを推定することができる。
【００５３】
図５の右下の表示は、図５のｘ軸が時間（Ｔｉｍｅ）であり、ｙ軸は大きさ（Ａｍｐｌｉｔｕｄｅ）であることを意味し、以下の図においても同様に適用される。
【００５４】
図６は、内部キャパシタのキャパシタンス推定ステップに適用される第２区間での電圧および電流波形を示す。
【００５５】
図６を参照すれば、開放電圧の変化量ΔＶＯＣは、第２センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ２）の開始地点ｔ０での開放電圧と第２センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ２）の終了時点ｔ１での開放電圧との差を意味する。内部キャパシタ２０１のキャパシタンスＣｉは、下記の数式８に表現したように、パルス電流供給期間中の第２センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ２）において、バッテリに供給される総電荷量ΔＱおよびバッテリ開放電圧の変化量ΔＶＯＣを用いて推定する。
【００５６】
【数８】
0006317541.tif 000009
【００５７】
ここで、総電荷量ΔＱは、Δｔ（＝ｔ１−ｔ０）時間の間にバッテリ１１０に供給された総電荷量であり、ｉは、バッテリ１１０に供給される電流であり、ΔＶＯＣは、総電荷量ΔＱによって増加したバッテリ開放電圧の変化量である。
【００５８】
第１周波数ｆ１と第２周波数ｆ２との間に定義される第３区間では、以下に説明する方式でダイナミック素子のパラメータを推定する。
【００５９】
図７は、ダイナミック素子のパラメータ推定ステップに適用される第３区間での電圧および電流波形を示す。
【００６０】
図７を参照すれば、ダイナミック素子のパラメータ推定ステップは、ダイナミック抵抗２０４の抵抗値Ｒｄとダイナミックキャパシタ２０３のキャパシタンスＣｄをそれぞれ求める２つの過程に区分することができ、パルス電流の１周期中に２つのセンシング過程（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ３、Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ４）が行われる。パルス電流の１周期中に電流の大きさが急激に増加した後、一定の電流値を維持する第３センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ３）ｔ２〜ｔ３の間には、ダイナミック抵抗２０４の抵抗値Ｒｄを推定し、第３センシング期間ｔ２〜ｔ３の経過後電流の大きさが急激に減少した後、一定の電流値を維持する第４センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ４）ｔ４〜ｔ７の間には、ダイナミックキャパシタ２０３のキャパシタンスＣｄを推定する。
【００６１】
まず、第３センシング期間ｔ２〜ｔ３で行われるダイナミック抵抗２０４の抵抗値Ｒｄを推定する過程について説明する。
【００６２】
第３センシング期間ｔ２〜ｔ３では、パルス電流の大きさが変化（ΔＩＢ）することによって変化するバッテリの電圧変化ΔＶＢは、数式９のように表示することができる。
【００６３】
【数９】
0006317541.tif 000010
【００６４】
バッテリ１１０の内部キャパシタ２０１のキャパシタンスＣｉがダイナミックキャパシタ２０３のキャパシタンスＣｄに比べて非常に大きいため、パルス電流の周波数を適切に調整すると、内部キャパシタ２０１のインピーダンスを、並列に接続されたダイナミック抵抗２０４とダイナミックキャパシタ２０３の並列インピーダンスに比べて大きくすることができる。したがって、パルス電流が急に低い電流レベルから高い電流レベルに遷移する瞬間ｔ２のパルス電流の変化は、内部キャパシタ２０１に降下する電圧ΔＶＯＳで示される。この瞬間ｔ２の後は、高い電流レベルを有するパルス電流は、内部キャパシタ２０１にはそれ以上影響を与えず、ダイナミック抵抗２０４およびダイナミックキャパシタ２０３の並列回路にすべて降下ΔＶｄする。
【００６５】
数式９に表示されたダイナミック抵抗２０４の抵抗値Ｒｄは、数式１０のようにまとめることができる。
【００６６】
【数１０】
0006317541.tif 000011
【００６７】
ここで、ｔ２は、バッテリ１１０に供給される電流の大きさ、すなわち、電流レベルが急に増加する時点であり、ｔ３は、第３センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ３）内でバッテリ１１０の電圧レベルが飽和（ｓａｔｕｒａｔｉｏｎ）する時点である。バッテリ２０１に供給される電流の大きさが急に減少する時点ｔ４の前にバッテリ電圧の大きさが飽和しなければならないが、パルス電流の周波数を適切に選択すればこのような条件を満足させることができる。これについては後述する。
【００６８】
第４センシング期間ｔ４〜ｔ７が開始される時点のｔ４は、バッテリ１１０に供給される電流の大きさが急に減少する瞬間であり、この時、内部キャパシタ２０１に降下するバッテリ１１０の開放電圧ΔＶＯＳは瞬時に減少するが、並列回路に降下した電圧ΔＶｄは指数関数的に減少する。これは、ダイナミックキャパシタ２０３にダイナミック抵抗２０４が並列に接続されているからであり、ｔ４の瞬間に、バッテリ電圧は並列回路に降下した電圧ΔＶｄになる。第４センシング期間ｔ４〜ｔ７の最後になる時点のｔ７は、ダイナミック抵抗２０４およびダイナミックキャパシタ２０３の時定数（ｔｉｍｅ  ｃｏｎｓｔａｎｔ；τ）の３倍（３τ）になる瞬間である。
【００６９】
説明の便宜のために、第４センシング期間ｔ４〜ｔ７をダイナミック抵抗２０４およびダイナミックキャパシタ２０３の時定数（ｔｉｍｅ  ｃｏｎｓｔａｎｔ；τ）の３倍（３τ）になる瞬間に定めたが、これより短い時間期間または長い時間期間に設定することも可能である。
【００７０】
以下、ダイナミックキャパシタ２０３のキャパシタンスＣｄを推定する過程について説明する。
【００７１】
図８は、第４センシング期間での電圧波形を拡大したものである。
【００７２】
図８を参照すれば、第４センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ４）の開始時点ｔ４でのバッテリの電圧のレベルΔＶｄは、時間が経過するほど指数関数的に減少し、第４センシング期間（Ｓｅｎｓｉｎｇ  Ｐｅｒｉｏｄ４）は、ダイナミック抵抗２０４の抵抗値Ｒｄおよびダイナミックキャパシタ２０３のキャパシタンスＣｄによって決定される時定数を用いて３つのサブ期間に区別する。
【００７３】
３つのサブ期間は、初期ｔ４の電圧レベルΔＶｄの約３６．２％の大きさになる瞬間ｔ５、初期ｔ４の電圧レベルΔＶｄの１３．５％の大きさになる瞬間ｔ６、および初期ｔ４のレベル電圧ΔＶｄの０．０５％の大きさになる瞬間ｔ７にそれぞれ区別される。
【００７４】
これらは、数式１１のように表現することができる。
【００７５】
【数１１】
0006317541.tif 000012
【００７６】
数式１１および時定数の概念を用いると、ダイナミックキャパシタ２０３のキャパシタンスＣｄを数式１２のように推定することができる。
【００７７】
【数１２】
0006317541.tif 000013
【００７８】
数式１２は、３つのサブ期間を平均してダイナミックキャパシタ２０３のキャパシタンスＣｄを求めるものと表現されているが、３つのサブ期間のうち、選択された１つのサブ期間に対して推定するか、または２つの選択されたサブ期間を平均してダイナミックキャパシタのキャパシタンスＣｄを推定することも可能である。
【００７９】
図９は、ダイナミック素子のパラメータ推定ステップに適したパルス電流の周波数を決定する過程を示す。
【００８０】
図９を参照すれば、ダイナミック素子のパラメータ推定ステップにおいて、バッテリに供給される電流の大きさが急に減少する時点ｔ４の前にバッテリ電圧の大きさが飽和すべき条件を満足させるパルス電流の周波数を決定するために、パルス電流の周波数を変更させていきながら、この時のバッテリの電圧（ΔＶｄ１、ΔＶｄ２、ΔＶｄ３、ΔＶｄ４）の差を一定の閾値Ｖｔｈと比較し、これは、数式１３のように表現することができる。
【００８１】
【数１３】
0006317541.tif 000014
【００８２】
ここで、ｎは、２以上の自然数である。
【００８３】
ここで、２つの互いに異なる周波数を有するパルス電流に対応するバッテリ電圧の差が閾値Ｖｔｈより小さい場合、２つの周波数のうちの１つの周波数をパルス電流の周波数に決定すると、バッテリに供給される電流の大きさが急に減少する時点ｔ４の前にバッテリ電圧の大きさが飽和すべき条件を満足する。
【００８４】
図１０は、本発明にかかるバッテリパラメータ推定方法を適用したコンピュータシミュレーションの結果を示す。
【００８５】
図１０を参照すれば、任意にパラメータ値を設定した後、本発明にかかるバッテリパラメータ推定方法を適用した時、一定時間の経過後任意に設定したパラメータ値を正確に推定することが分かる。前記コンピュータシミュレーションにおいて、内部抵抗２０２の抵抗値Ｒｓは０．１５Ω、内部キャパシタ２０１のキャパシタンスＣｉは２０００ｐＦ、ダイナミック抵抗２０４の抵抗値Ｒｄは０．０４Ω、そして、ダイナミックキャパシタ２０３のキャパシタンスＣｄは３ｐＦにそれぞれ設定した。シミュレーションの開始後３０秒前に、内部抵抗２０２の推定抵抗値（Ｅｓｔｉｍａｔｅｄ  Ｒｓ）、内部キャパシタ２０１の推定キャパシタンス（Ｅｓｔｉｍａｔｅｄ  Ｃｉ）、ダイナミック抵抗２０４の推定抵抗値（Ｅｓｔｉｍａｔｅｄ  Ｒｄ）、そして、ダイナミックキャパシタ２０３の推定キャパシタンスＣｄはそれぞれ設定した値と等しい値と推定され、内部抵抗２０２の抵抗値Ｒｉが最も遅く推定されたことが分かる。
【００８６】
図１０の最上部の実験結果を参照すれば、シミュレーションにおいて、内部キャパシタ２０１の両端電圧ＶＣは、時間に応じて変化しない一定の値を有するが、パルス状のバッテリ１１０の推定された両端電圧ＶＯＣは、シミュレーションが開始された時から約３０秒後に内部キャパシタ２０１の両端電圧ＶＣに収斂することが分かる。
【００８７】
また、図１０に示された実験結果の上から２番目の実験結果を参照すれば、時間に応じて変化しない一定の値を有する内部キャパシタ２０１の両端電圧ＶＣとバッテリ１１０に供給される電圧（Ｖｓｅｎｓｉｎｇ）との差（Ｖｓｅｎｓｉｎｇ−ＶＣ；Ｓｅｎｓｉｎｇ  ｅｒｒｏｒ）は一定の値に収斂しないため、単にバッテリ１１０に供給される電圧（Ｖｓｅｎｓｉｎｇ）と内部キャパシタ２０１の両端電圧ＶＣを有するだけではいかなる値も類推することができない。しかし、時間に応じて変化しない一定の値を有する内部キャパシタ２０１の両端電圧ＶＣとバッテリ１１０の推定された両端電圧ＶＯＣとの差（ＶＯＣ−ＶＣ；Ｅｓｔｉｍａｔｅｄ  ｅｒｒｏｒ）は、シミュレーションが開始されて約３０秒後に０（ｚｅｒｏ）に収斂しているため、本発明にかかるバッテリパラメータ推定方法を用いる場合、シミュレーションが開始された時から約３０秒後は、内部キャパシタ２０１の両端電圧ＶＣを正確に類推できることが分かる。
【００８８】
図１１は、本発明によって推定されたインピーダンス値を基に推定した開放電圧と内部キャパシタに降下した電圧との関係を示す。
【００８９】
図１１の上部に示された推定された開放電圧（細い実線）と内部キャパシタ２０１の実際のキャパシタンスＶＣの電圧（太い実線）との関係をみると、本発明によって推定された開放電圧（細い実線）と内部キャパシタ２０１の実際のキャパシタンス（太い実線）時間に対する増加傾向はほぼ類似することが分かる。また、図１１の下部に示された開放電圧ＶＯＣと内部キャパシタ２０１のキャパシタンスＶＣとの差ＶＯＣ−ＶＣをみると、最大誤差が０．６ｍＶになることが分かり、これは無視すべき誤差である。
【００９０】
図１１を参照すれば、本発明にかかる推定したバッテリモデルのパラメータを用いる場合、バッテリに頻繁な電流の出入りがあっても、小さい推定誤差（０．６ｍＶ以下）を有しながら開放電圧ＶＯＣを推定することが可能であることが分かる。
【００９１】
以上、本発明に関する技術思想を添付図面とともに述べたが、これは、本発明の好ましい実施形態を例示的に説明したものであって、本発明を限定するものではない。また、本発明の属する技術分野における通常の知識を有する者であれば、誰でも本発明の技術的思想の範疇を離脱しない範囲内で多様な変形および模倣可能であることが自明である。
【符号の説明】
【００９２】
１００：バッテリパラメータ管理システム
１１０：バッテリ
１０１：温度計
１０２：電流計
１０３：電圧計
１２０：センシング抵抗
１３０：制御スイッチ部
１４０：電流供給源
１５０：負荷
１６０：プロセッサ
１７０：ユーザインタフェース装置
【図１】
0006317541.tif 000020
【図２】
0006317541.tif 000021
【図３】
0006317541.tif 000022
【図４】
0006317541.tif 000023
【図５】
0006317541.tif 000024
【図６】
0006317541.tif 000025
【図７】
0006317541.tif 000026
【図８】
0006317541.tif 000027
【図９】
0006317541.tif 000028
【図１０】
0006317541.tif 000029
【図１１】
0006317541.tif 000030