(57)【要約】      （修正有）【課題】  速やか且つ簡便に実際の充放電状態において電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。【解決手段】  電力貯蔵供給デバイス１００に対して充放電部１で充放電を行ない、測定部２で充放電開始後の電流値又は電圧値を所定時間間隔で測定し、計算部３で各時点における電流値又は電圧値の変化を用い電力貯蔵供給デバイス１００を周波数に応じて抵抗値が変化する抵抗成分と周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に抵抗成分の抵抗値に相当する第１パラメータ及び上記キャパシタ成分の静電容量値に相当する第２パラメータを演算し、第１のパラメータ及び第２のパラメータを用いて電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。【選択図】      図４
2006162283.tif 000002
【特許請求の範囲】
【請求項１】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価方法であって、
上記電力貯蔵供給デバイスに対して充電又は放電を行ない、上記の充電又は放電を開始した後の電流値又は電圧値を所定時間間隔で測定する第１過程と、
各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２過程と、
上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する第３過程とを備える
ことを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価方法。
【請求項２】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価方法であって、
上記電力貯蔵供給デバイスに対して充電又は放電を行ない、少なくとも１回上記の充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させた際に、電流値又は電圧値を所定時間間隔で測定する第１過程と、
電流又は電圧を変化させた際の各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２過程と、
上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する第３過程とを備える
ことを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価方法。
【請求項３】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、
上記電力貯蔵供給デバイスに対して充電又は放電を行なう充放電部と、
電流値又は電圧値を測定しうる測定部と、
上記の充電又は放電を開始した後の電流値又は電圧値を所定時間間隔で測定するように該測定部を制御する測定制御部と、
各時点における上記電流値又は電圧値に基づいて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する計算部とを備える
ことを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置。
【請求項４】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、
上記電力貯蔵供給デバイスに対して充電又は放電を行なう充放電部と、
少なくとも１回以上、該充電部が充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させるパルス印加部と、
電流値又は電圧値を測定しうる測定部と、
上記の電流又は電圧を変化させた際に、電流値又は電圧値を所定時間間隔で測定するように該測定部を制御する測定制御部と、
電流又は電圧を変化させた際の各時点における上記電流値及び電圧値に基づいて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する計算部とを備える
ことを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置。
【請求項５】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、
上記電力貯蔵供給デバイスに対して充電又は放電を行なったときの、上記の充電又は放電を開始した後の所定時間間隔の電流値又は電圧値を取り込む測定値取込部と、
取り込んだ上記電流値又は電圧値の各時点における変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータを演算する第１パラメータ演算部と、
取り込んだ上記電流値又は電圧値の各時点における変化を用い、上記電力貯蔵供給デバイスを上記等価回路で表わした場合に、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２パラメータ演算部と、
上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価部とを備える
ことを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置。
【請求項６】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、
上記電力貯蔵供給デバイスに対して充電又は放電を行ない、少なくとも１回上記の充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させた際の、所定時間間隔の電流値又は電圧値を取り込む測定値取込部と、
取り込んだ上記電流値又は電圧値の、電流又は電圧を変化させた際の各時点における変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータを演算する第１パラメータ演算部と、
取り込んだ上記電流値又は電圧値の、電流又は電圧を変化させた際の各時点における変化を用い、上記電力貯蔵供給デバイスを上記等価回路で表わした場合に、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２パラメータ演算部と、
上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価部とを備える
ことを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置。
【請求項７】
電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価プログラムであって、
コンピュータに、請求項５又は請求項６に記載の周波数特性評価装置の各機能を実現させる
ことを特徴とする、周波数特性評価プログラム。
【請求項８】
請求項５又は請求項６に記載の電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置で評価された上記周波数特性に基づき上記電力貯蔵供給デバイスの性能を評価するデバイス性能評価部と、
該デバイス性能評価部での評価結果に基づいて上記電力貯蔵供給デバイスの充電又は放電の制御を行なう充放電制御部とを備える
ことを特徴とする、保護装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価方法、周波数特性評価装置及び周波数特性評価プログラム、並びに、それを用いた保護装置に関する。
【背景技術】
【０００２】
電力貯蔵供給デバイスとしては、リチウムイオン電池、ニッケル水素電池、リチウム電池等の各種電池の１次電池や２次電池、電解コンデンサ、電気二重層キャパシタ等のキャパシタ、燃料電池等の発電型電池などが開発され、工業的に広く用いられている。
上記のような電力貯蔵供給デバイスは、一般に高容量密度、高出力密度が求められる傾向にあり、より大きい電流において電力を貯蔵、供給することが望まれ、盛んに開発が進められている。
【０００３】
ところで、電力貯蔵供給デバイスの開発を行なう際には、その内部インピーダンスを評価することがある。例えば、電力貯蔵供給デバイスの一種である電池の充放電の許容負荷電流は、電池の内部インピーダンスが大きく影響する。したがって、電池の内部インピーダンスの周波数特性を評価することは、電池の性能向上に非常に有用な情報を得ることを意味する。
【０００４】
また、電力貯蔵供給デバイスは、一般に正極と負極とが近接していて内部短絡が発生する可能性の高いデバイスである。このため、電力貯蔵供給デバイスの工業的な利用にあたっては、短絡やそれに伴うデバイス構成材料の熱分解を抑制することによりその安全性を高めることが求められている。最近では、この安全性の向上のため、内部短絡発生時の電力貯蔵供給デバイス内の温度分布のシミュレーションを行ない、その安全性の検証を行なうことが提案されている（非特許文献１，２）。
【０００５】
上記のシミュレーションを行なう場合には、電力貯蔵供給デバイスに短絡時にどのくらいの電流が流れるか評価しなくてはならないことがある。通常、短絡時に流れる電流の大きさは、電力貯蔵供給デバイスに蓄えられたエネルギー量とともに、短絡部分の抵抗や短絡の接触抵抗、そして電力貯蔵供給デバイスの内部インピーダンスに依存する。したがって、内部短絡発生時に生じるジュール熱を、シミュレーションを用いてより精密に評価するためには、内部インピーダンスの周波数特性を求め、その周波数特性に基づいた短絡電流の解析を行なうことが望ましい。
【０００６】
上述したような電力貯蔵供給デバイスの内部インピーダンスの周波数依存性を評価する場合、従来は、交流負荷に対する周波数応答による直接測定方法や、高速フーリエ変換による測定方法が良く用いられていた（非特許文献３）。
【０００７】
【非特許文献１】２００３年度電池討論会予稿集  第４０４頁〜第４０５頁
【非特許文献２】Ｊｏｕｒｎａｌ  ｏｆ  Ｐｏｗｅｒ  Ｓｏｒｃｅｓ，１１３，８１−１００（２００３）
【非特許文献３】Ａｎａｌｙｔｉｃａｌ  Ｃｈｅｍｉｓｔｒｙ  Ａ−Ｐａｇｅｓ，７５（２１），４５５Ａ−４６１Ａ（２００３）
【発明の開示】
【発明が解決しようとする課題】
【０００８】
電力貯蔵供給デバイスの内部インピーダンスは、過電圧によって変化する傾向があり、また、充電深度によって大きく変化する傾向がある。特に、電池などの電力貯蔵供給デバイスにおいて大きな内部短絡が起こった場合には、過電圧が大きい上に放電の進行とともに著しくそれが変化するため、従来、実際の短絡状態に対応した内部インピーダンスを正確に評価することが困難であった。このため、上記のシミュレーションに用いるための、短絡状態などにおける電力貯蔵供給デバイスの内部インピーダンスの周波数特性を正確に評価することは、非常に困難であった。
【０００９】
また、電力貯蔵供給デバイスにおいては、充放電の進行中には過電圧が刻々と変化するため、それに対応して内部インピーダンスが過渡的に変化する。しかし、従来の方法を用いて内部インピーダンスを評価する場合、実際の充放電時の過電圧にインピーダンス評価時の負荷信号（交流）の振幅の大きさを対応させて評価する必要があるなど、手間、時間がかかっていた。このため、過渡状態における電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することが困難であった。
【００１０】
さらに、従来の具体的な内部インピーダンスの周波数依存性の評価方法に注目すると、直接測定方法では、電力貯蔵供給デバイスに交流で充放電を行ない、その際の内部インピーダンスを通常のインピーダンスメーターなどで測定するが、１点のインピーダンスデータの測定に時間がかかったため、内部インピーダンスの周波数特性を評価するのにも時間がかかったほか、電力貯蔵供給デバイスの過渡特性、即ち、インピーダンスの周波数特性の充放電中の推移（経時変化）についての連続したデータを取ることが困難であり、また、労力がかかった。
【００１１】
また、例えば高速フーリエ変換による測定方法では、内部インピーダンスの高速な評価は可能であるが、電力貯蔵供給デバイスの実際の使用時における過電圧の経時変化にインピーダンス評価時の過電圧を対応させることが困難なため、電力貯蔵供給デバイスの過渡特性の適切な評価はできなかった。
【００１２】
また、この方法はインピーダンスが過電圧に対して線形であることを前提としているが、化学電池のように電池を得るために化学反応を使用しているデバイスの場合、インピーダンスは過電圧が大きくなるにつれて非線形性が強くなるため、過電圧の大きい大電流時や低温時のインピーダンス評価には向かなかった。
【００１３】
本発明は上記の課題に鑑みて創案されたもので、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができるようにした周波数特性評価方法、周波数特性評価装置及び周波数特性評価プログラム、並びに、それを用いた保護装置を提供することを目的とする。
【課題を解決するための手段】
【００１４】
本発明の発明者は、上記課題を解決するべく鋭意検討を行なった結果、電力貯蔵供給デバイスに対して充電又は放電を行なったときの電流値又は電圧値の変化、若しくは、電流又は電圧を変化させて充電又は放電を行なったときの電流値又は電圧値の変化から、電力貯蔵供給デバイスを抵抗成分及びキャパシタ成分を有する等価回路で表わした場合の抵抗値及び静電容量に相当する値を演算し、その値を用いることにより、上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができることを見出し、本発明を完成させた。
【００１５】
即ち、本発明の要旨は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価方法であって、上記電力貯蔵供給デバイスに対して充電又は放電を行ない、上記の充電又は放電を開始した後の電流値又は電圧値を所定時間間隔で測定する第１過程と、各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２過程と、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する第３過程とを備えることを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価方法に存する（請求項１）。これにより、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。また、この方法は、充放電開始直後に電圧値が大きく変化している場合に用いると、その間に内部インピーダンスも短時間で刻々と変化するため、充放電時の過電圧推移そのものから内部インピーダンスの周波数特性を評価できるという本法の長所を発揮できるため、好適である。
【００１６】
また、本発明の別の要旨は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価方法であって、上記電力貯蔵供給デバイスに対して充電又は放電を行ない、少なくとも１回上記の充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させた際に、電流値又は電圧値を所定時間間隔で測定する第１過程と、電流又は電圧を変化させた際の各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２過程と、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する第３過程とを備えることを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価方法に存する（請求項２）。これによっても、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。また、この方法は、特に、充放電中にわずかな電気的刺激を加えるだけで充放電を停止することなく内部インピーダンスの周波数特性を評価できるため、例えば充放電の電流が非常に大きい場合や、過充電又は過放電等の異常充放電状態の場合などのように、充放電中における内部インピーダンスの変化が大きい場合の内部インピーダンスの周波数特定を評価する場合に用いて好適である。
【００１７】
さらに、本発明の更に別の要旨は、  電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、上記電力貯蔵供給デバイスに対して充電又は放電を行なう充放電部と、電流値又は電圧値を測定しうる測定部と、上記の充電又は放電を開始した後の電流値又は電圧値を所定時間間隔で測定するように該測定部を制御する測定制御部と、各時点における上記電流値又は電圧値に基づいて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する計算部とを備えることを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置に存する（請求項３）。また、該計算部は、例えば、各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータを演算する第１パラメータ演算部と、各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを上記等価回路で表わした場合に、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２パラメータ演算部と、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価部とを備えて構成することができる。これによっても、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。また、この装置は、充放電開始直後に電圧値が大きく変化している場合に用いると、その間に内部インピーダンスも短時間で刻々と変化するため、充放電時の過電圧推移そのものから内部インピーダンスの周波数特性を評価できるという本法の長所を発揮できるため、好適である。
【００１８】
また、本発明の更に別の要旨は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、上記電力貯蔵供給デバイスに対して充電又は放電を行なう充放電部と、少なくとも１回以上、該充電部が充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させるパルス印加部と、電流値又は電圧値を測定しうる測定部と、上記の電流又は電圧を変化させた際に、電流値又は電圧値を所定時間間隔で測定するように該測定部を制御する測定制御部と、電流又は電圧を変化させた際の各時点における上記電流値及び電圧値に基づいて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する計算部とを備えることを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置に存する（請求項４）。また、該計算部は、例えば、電流又は電圧を変化させた際の各時点における上記電流値及び電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータを演算する第１パラメータ演算部と、電流又は電圧を変化させた際の各時点における上記電流値及び電圧値の変化を用い、上記電力貯蔵供給デバイスを、上記等価回路で表わした場合に上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２パラメータ演算部と、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価部とを備えて構成することができる。これによっても、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。また、この装置は、特に、充放電中にわずかな電気的刺激を加えるだけで充放電を停止することなく内部インピーダンスの周波数特性を評価できるため、例えば充放電の電流が非常に大きい場合や、過充電又は過放電等の異常充放電状態の場合などのように、充放電中における内部インピーダンスの変化が大きい場合の内部インピーダンスの周波数特定を評価する場合に用いて好適である。
【００１９】
さらに、本発明の更に別の要旨は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、上記電力貯蔵供給デバイスに対して充電又は放電を行なったときの、上記の充電又は放電を開始した後の所定時間間隔の電流値又は電圧値を取り込む測定値取込部と、取り込んだ上記電流値又は電圧値の各時点における変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータを演算する第１パラメータ演算部と、取り込んだ上記電流値又は電圧値の各時点における変化を用い、上記電力貯蔵供給デバイスを上記等価回路で表わした場合に、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２パラメータ演算部と、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価部とを備えることを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置に存する（請求項５）。これによっても、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。また、この装置は、充放電開始直後に電圧値が大きく変化している場合に用いると、その間に内部インピーダンスも短時間で刻々と変化するため、充放電時の過電圧推移そのものから内部インピーダンスの周波数特性を評価できるという本法の長所を発揮できるため、好適である。
【００２０】
また、本発明の更に別の要旨は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価装置であって、上記電力貯蔵供給デバイスに対して充電又は放電を行ない、少なくとも１回上記の充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させた際の、所定時間間隔の電流値又は電圧値を取り込む測定値取込部と、取り込んだ上記電流値又は電圧値の、電流又は電圧を変化させた際の各時点における変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータを演算する第１パラメータ演算部と、取り込んだ上記電流値又は電圧値の、電流又は電圧を変化させた際の各時点における変化を用い、上記電力貯蔵供給デバイスを上記等価回路で表わした場合に、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する第２パラメータ演算部と、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価部とを備えることを特徴とする、電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置に存する（請求項６）。これによっても、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。また、この装置は、特に、充放電中にわずかな電気的刺激を加えるだけで充放電を停止することなく内部インピーダンスの周波数特性を評価できるため、例えば充放電の電流が非常に大きい場合や、過充電又は過放電等の異常充放電状態の場合などのように、充放電中における内部インピーダンスの変化が大きい場合の内部インピーダンスの周波数特定を評価する場合に用いて好適である。
【００２１】
さらに、本発明の更に別の要旨は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価プログラムであって、コンピュータに、上記の周波数特性評価装置の各機能を実現させることを特徴とする、周波数特性評価プログラムに存する（請求項７）。これによっても、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。
【００２２】
また、本発明の別の要旨は、上記の電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価装置で評価された上記周波数特性に基づき上記電力貯蔵供給デバイスの性能を評価するデバイス性能評価部と、該デバイス性能評価部での評価結果に基づいて上記電力貯蔵供給デバイスの充電又は放電の制御を行なう充放電制御部とを備えることを特徴とする、保護装置に存する（請求項８）。これにより、実際の充放電状態において、内部インピーダンスの周波数特性を基に性能を評価して電力貯蔵供給デバイスを制御することができるようになる。
【発明の効果】
【００２３】
本発明の電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する周波数特性評価方法、周波数特性評価装置及び周波数特性評価プログラムによれば、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。
また、本発明の保護装置によれば、実際の充放電状態において、内部インピーダンスの周波数特性を基に性能を評価して、例えば過充電などから電力貯蔵供給デバイスを保護する等のために、電力貯蔵供給デバイスを制御することができるようになる。
【発明を実施するための最良の形態】
【００２４】
以下、本発明について、その実施の形態を示して説明するが、本発明は以下の実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において、任意に変形して実施することができる。
［Ｉ．第１実施形態］
本発明は、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価するためのものである。電力貯蔵供給デバイスとは、電力を貯蔵又は供給することが可能なデバイスを広く意味する。電力貯蔵供給デバイスの具体例としては、アルカリ電池、リチウム電池等の１次電池、鉛酸電池、ニカド電池、ニッケル水素電池、リチウムイオン電池等の２次電池、燃料電池等の発電型電池、電解コンデンサ、電気二重層キャパシタ等のキャパシタなどが挙げられる。
【００２５】
また、本発明の第１実施形態としての周波数特性評価方法では、以下に説明する第１過程、第２過程及び第３過程を経て、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。
即ち、第１過程では、電力貯蔵供給デバイスに対して充電又は放電を行ない、上記の充電又は放電を開始した後の電流値又は電圧値を所定時間間隔で測定する。
また、第２過程では、測定した各時点における電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する。
さらに、第３過程では、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。
【００２６】
なお、本実施形態の周波数特性評価方法はどのような場合に適用してもよい。例えば、充放電開始後の電圧値（過電圧）の変化が大きいときに適用することができる。即ち、充放電開始後の電圧値（過電圧）の変化が大きい場合には、電力貯蔵供給デバイス内の電荷移動反応速度の過電圧に対する非線形性より、内部インピーダンスが充放電の進行とともに大きく変化している。従来の方法では、このような条件での充放電中の内部インピーダンスの周波数特性を評価することは困難が伴なっていた。しかし、本実施形態の周波数特性評価方法では、充放電中の電圧そのものを数値処理することによってインピーダンスの周波数特性を評価することができるので、このような条件下でも内部インピーダンスの周波数特性を簡便に評価することができる。
以下、各過程毎に詳細に説明する。
【００２７】
［１．第１過程］
第１過程では、上述したように、電力貯蔵供給デバイスに対して直流で充電又は放電を行ない、上記の充電又は放電を開始した後の電流値又は電圧値を所定時間間隔で測定する。
充電と放電とのいずれを行なうかは任意であり、充電のみを行なうようにしてもよく、放電のみを行なうようにしてもよく、充電及び放電の両方を組み合わせて行なうようにしてもよい。通常は、充電時の内部インピーダンスの周波数特性を評価したい場合には充電を行なうようにし、放電時の内部インピーダンスの周波数特性を評価したい場合には放電を行なうようにする。なお、本明細書において、充電及び放電を区別せずに述べる場合、単に「充放電」という。
【００２８】
また、充放電を行なう前に、電力貯蔵供給デバイスは、予め所望の充電深度及び所望の温度に調整しておくことが好ましい。この際、充電深度や温度は、評価しようとする実際の環境に応じて任意に設定すればよい。例えば、通常の使用状態における周波数特性を評価する場合には、充電深度及び温度は、実際に使用する場合の充電深度及び使用温度に調整しておけばよい。さらに、例えば、シミュレーションに用いる周波数特性を評価する場合には、充電深度及び温度は、シミュレーションにおいて想定する充電深度及び想定温度に調整しておけばよい。
さらに、充電速度や放電速度、充放電を行なう時間などの充放電条件も任意であり、周波数特性を評価しようとする実際の環境等に応じて任意に設定すればよい。
【００２９】
また、上記のように、第１過程においては、充放電の開始後、電流値又は電圧値を所定時間間隔で測定する。このとき、測定される電圧値（即ち、電力貯蔵供給デバイスの正極及び負極の端子間の電圧）や電流値は、充放電開始後の電力貯蔵供給デバイスの内部インピーダンスの周波数特性を反映したものとなるため、これらを測定することにより、上記周波数特性の評価を行なうことができる。
【００３０】
第１過程において、電流値と電圧値とのいずれを測定するかは、第２過程で行なう演算に応じて適当に選択すればよい。例えば、電力貯蔵供給デバイスに定電流で充放電を行なう場合には、電流値は一定の定数値になることが分かっているため、周波数特性の評価を行なうには電圧値を測定するようにすれば良い。また、例えば電力貯蔵供給デバイスに定電圧で充放電を行なう場合には、電圧値は一定の定数値になることが分かっているため、周波数特性の評価を行なうには電流値を測定するようにすれば良い。さらに、電圧値及び電流値の両方が変化しうるように充放電を行なう場合には、周波数特性の評価を行なうには、電流値及び電圧値の両方を測定するようにすることが好ましい。もちろん、定電圧や定電流で充放電を行なう場合に、電流値及び電圧値の両方を測定するようにしてもよい。この場合、周波数特性の評価結果の精度や信頼性をより高めることができる。
【００３１】
また、電流値や電圧値は、所定時間間隔で測定する。所定時間は、周波数特性の評価ができる限り任意であり、例えば測定開始から終了に向けて徐々に長時間となるように測定間隔を広げるようにしてもよいが、簡単のため、一定の時間とすることが好ましい。さらに、この所定時間（測定間隔）を小さくすれば、高周波数相当の内部インピーダンスの周波数特性まで評価することができる。
【００３２】
上記所定時間の具体的な範囲としては、その下限は、通常１μｓ以上、好ましくは１０μｓ以上、より好ましくは１００μｓ以上が望ましい。この下限を下回ると、データ点数が多くなりすぎ、内部インピーダンスの評価にあたり計算機の負荷が大きくなる虞があるためである。一方、上限については制限は無いが、集積回路に組み込む微小デバイス等の場合には、通常１０００ｍｓ以下、好ましくは１００ｍｓ以下、より好ましくは１０ｍｓ以下が望ましい。この上限を上回ると、高周波数領域での内部インピーダンス評価ができず、電力貯蔵供給デバイス性能向上に必要な情報が得られない虞があるためである。
【００３３】
また、測定開始から測定終了までの間に測定する測定データ点数に制限はなく、周波数特性の評価ができる限り任意であるが、通常１０点以上、好ましくは２０点以上、より好ましくは５０点以上、また、通常１０００００００点以下、好ましくは１００００００点以下、より好ましくは１０００００点以下が望ましい。この範囲の下限を下回ると、内部インピーダンスの周波数応答曲線が滑らかな形状に得られないため、即ち、周波数特性の信頼性が低下する虞があるため、その結果を用いて得られる熱シミュレーションの結果の精度が低下する虞がある。また、上限を上回ると、データ点数が多くなりすぎ、インピーダンスの算出にあたり計算機の負荷が大きくなる虞がある。
【００３４】
さらに、測定開始から測定終了までの総測定時間も制限なく任意であるが、通常０．０１秒以上、好ましくは０．１秒以上、より好ましくは０．５秒以上、また、通常６０秒以下、好ましくは３０秒以下、より好ましくは１０秒以下である。この範囲の下限を下回ると、電力貯蔵供給デバイスのインピーダンスの低周波数特性が充分に得られない虞がある。また、上限を上回ると、上記の測定間隔（即ち、所定時間）を微少な時間とすることの意味がなくなり、周波数特性の評価に要する時間を短縮するという利点を得られなくなる虞がある。
【００３５】
また、測定する電流値や電圧値の精度は、高いほど好ましい。具体的には、電流値の場合、その最小測定単位が通常１００ｍＡ以下、好ましくは１０ｍＡ以下、より好ましくは１ｍＡ以下が望ましい。一方、電圧値の場合、その最小測定単位が通常１０ｍＶ以下、好ましくは１ｍＶ以下、より好ましくは０．１ｍＶ以下が望ましい。これらの範囲の上限を上回った場合、得られる内部インピーダンスの周波数特性の精度が低下する虞がある。
【００３６】
［２．第２過程］
第２過程では、第１過程で測定した各時点における電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する。なお、第１過程において測定した時点のうちの一部の時点の測定結果についてだけ上記の演算を行なうようにしてもよいが、周波数特性評価の精度を高める観点からは、測定した全ての時点の測定結果から上記演算を行なうことが好ましい。
【００３７】
充放電開始後の電流と電圧との関係は、電力貯蔵供給デバイスの内部インピーダンスを反映した挙動を示す。したがって、第１過程で測定した電流値及び電圧値の時間挙動を解析すれば、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。
具体的には、例えば、以下に示す式を用いて演算を行なうことができる。なお、以下に示す式の一部は、「“Ｅｌｅｃｔｒｏｃｈｅｍｉｃａｌ  Ｍｅｔｈｏｄｓ，ｓｅｃｏｎｄ  ｅｄｉｔｉｏｎ”  Ａｌｌｅｎ  Ｊ．  Ｂａｒｄ，Ｌａｒｒｙ  Ｒ．  Ｆａｕｌｋｎｅｒ（Ｊｏｈｎ  Ｗｉｌｅｙ＆Ｓｏｎｓ）、２００１」に記載のものを用いた。
【００３８】
一般に、電力貯蔵供給デバイスの内部インピーダンスＺは、式（１）に示すように、実数成分ＺREAL及び虚数成分Ｚimaginaryよりなる。
【数１】
2006162283.tif 000003
なお、上記式（１）において、ｊは虚数単位を表わす。
【００３９】
ここで、電力貯蔵供給デバイスの内部インピーダンスが、図１に示すような、抵抗（周波数に応じて抵抗値が変化する抵抗成分）とキャパシタ（周波数に応じて静電容量値が変化するキャパシタ成分）との直列接続により構成される等価回路で表わされると仮定すると、以下の式（２）、式（３）が成り立つ。
【数２】
2006162283.tif 000004
ここで、Ｒは上記等価回路の抵抗の抵抗値（第１パラメータ）を表わし、Ｃは上記等価回路のキャパシタの静電容量（第２パラメータ）を表わし、ωは負荷された交流の角速度を表わす。なお、本実施形態の周波数特性評価方法では、直流時の電流又は電圧の時間変化を測定すること、並びに、抵抗値Ｒ及び静電容量Ｃが時間とともに変化すると想定したことにより、この角速度ωを変化させる手間を省きながらも電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができることを、利点の一つとしている。
【００４０】
上記の式（１）〜式（３）を変形すると、以下の式（４）及び式（５）が導かれる。
【数３】
2006162283.tif 000005
【００４１】
一方、電力貯蔵供給デバイスに充放電を開始した後の電流又は電圧が変化している期間や、電流又は電圧の印加を開始した直後から印加終了までの期間（即ち、第２実施形態におけるパルスの印加時等）の電位変化は下記式（６）で与えられる。
【数４】
2006162283.tif 000006
ここで、ΔＥは電流又は電圧印加後の電圧変化、具体的には、充放電時における充放電開始前からの電圧変化を表わし、ΔＩは電流又は電圧印加後の電流変化、具体的には、充放電時における充放電開始前からの電流変化を表わし、Δqは電流又は電圧印加後のキャパシタ成分に蓄えられた電気量をそれぞれ示す。
【００４２】
上記式（６）は、さらに変形して下記式（７）のようになる。
【数５】
2006162283.tif 000007
ここで、ｔは時間を表わし、具体的には、充放電を開始してからの経過時間を表わす。
【００４３】
なお、電力貯蔵供給デバイスを流れる電流が定電流である場合、例えば、定電流で充放電を行なう場合や、後述する第２実施形態のようにパルスを印加して電流の波形を変化させた際に電流の波形が矩形波や階段状波になった場合などには、上記式（７）は、さらに下記式（８）のように変形することができる。また、連続する電流及び電圧の変化を各時点における階段状波の集合として取り扱うこともでき、その集合を構成する階段状波それぞれについても、下記式（８）を適用することができる。
【数６】
2006162283.tif 000008
【００４４】
また、上記式（８）を微分形で表すと、下記式（９）のようになる。
【数７】
2006162283.tif 000009
【００４５】
以上の式（１）〜式（９）を用いれば、電力貯蔵供給デバイスに電流又は電圧の印加を開始した後の電流値又は電圧値の変化を用いて、印加開始後の時間に対する抵抗値Ｒ及び静電容量Ｃの変化が得られ、さらに、この抵抗値Ｒ及び静電容量Ｃの値より印加開始後の時間に対するインピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryの推移が得られる。
【００４６】
なお、本実施形態においては、第１過程で充放電を行なう際に、電力貯蔵供給デバイスに直流の電流又は電圧を印加することにより充放電を行なっている。このように直流の電流又は電圧を印加した場合に、印加開始後の時間ｔに対する実数成分ＺREAL及び虚数成分Ｚimaginaryの変化の関係を、交流の電流又は電圧を印加した際の交流の角速度ωに対する実数成分ＺREAL及び虚数成分Ｚimaginaryの変化の関係に変換するには、印加開始後の時間ｔの逆数１／ｔを角速度ωと見なせば良い。また、角速度ωは、さらに必要に応じて周波数に変えることもできる。
【００４７】
また、直流で充放電を行なう場合（第２実施形態のように電流又は電圧の波形を矩形波状又は階段状に変化させる場合を含む）には、インピーダンスの絶対値、電圧変化及び電流変化の関係として、以下の式を用いることができる。
【数８】
2006162283.tif 000010
【００４８】
第２過程においては、例えば上述した式（１）〜式（１０）を用いて、具体的に充放電開始後の時間に対する上記抵抗値Ｒ及び上記静電容量Ｃの推移を得る。即ち、第１過程において電流又は電圧を測定した時点などの、サンプリングした各時点における上記抵抗値Ｒ及び上記静電容量Ｃを求める。このとき、上記の式（１）〜式（１０）の使い方には種々の組み合わせが考えられるが、各時点における上記抵抗値Ｒ及び上記静電容量Ｃを求めることができる限りその使い方は任意である。以下に、使い方の例を示す。なお、図２（ａ）〜図２（ｃ）は、説明のために第１過程において定電流で放電を行なった場合の電圧の測定値をプロットした曲線の一例として示すグラフである。
【００４９】
（算出法１）
第１過程において、定電流で充放電を行ない、充放電開始後の電圧を測定する。
また、第２過程では、第１過程で測定した電圧と時間ｔとの関係｛（電圧変化）対（時間）のプロットで得られる曲線｝に対し、上記式（８）を用いる。具体的には、例えば、充放電時電流をΔＩ、充放電開始直前の電圧と充放電開始後の電圧の差をΔＥとし、図２（ａ）に示すように、ある時点に対して曲線の接線を得て、その切片すなわちΔI×Ｒより上記抵抗値Ｒを求め、また、その傾きすなわちΔＩ／Ｃより上記静電容量Ｃを求める。なお、ΔＩは、充電側の値を正とし、放電側の値を負とする。また、上記のように、印加開始後の時間ｔの逆数１／ｔを角速度ωとする。このようにして求めた上記抵抗値Ｒ及び上記静電容量Ｃの値により、上記式（２）及び上記式（３）を用いて各時点ｔにおける内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryが得られる。
【００５０】
（算出法２）
第１過程において、定電流で充放電を行ない、充放電開始後の電圧を測定する。
また、第２過程では、第１過程で測定した電圧と時間ｔとの関係｛（電圧変化）対（時間）のプロットで得られる曲線｝に対し、上記式８を用いる。具体的には、例えば、充放電時電流をΔＩ、充放電開始直前の電圧と充放電開始後の電圧の差をΔＥとし、図２（ｂ）に示すように、ある時点に対して曲線の接線を得て、その切片すなわちΔI×Ｒより上記抵抗値Ｒを求める。また、別途、対応する時点について上記式（１０）から内部インピーダンスの絶対値｜Ｚ｜を算出し、このようにして求めた上記抵抗値Ｒ及び絶対値｜Ｚ｜の値を用い、上記式（５）から上記静電容量Ｃが得られる。なお、ΔＩは、充電側の値を正とし、放電側の値を負とする。また、上記のように、印加開始後の時間ｔの逆数１／ｔを角速度ωとする。このようにして求めた上記抵抗値Ｒ及び上記静電容量Ｃの値により、上記式（２）及び上記式（３）を用いて各時点ｔにおける内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryが得られる。
【００５１】
（算出法３）
第１過程において、定電流で充放電を行ない、充放電開始後の電圧を測定する。
また、第２過程では、第１過程で測定した電圧と時間ｔとの関係｛（電圧変化）対（時間）のプロットで得られる曲線｝に対し、上記式８を用いる。具体的には、例えば、充放電時電流をΔＩ、充放電開始直前の電圧と充放電開始後の電圧の差をΔＥとし、図２（ｃ）に示すように、ある時点に対して曲線の接線を得て、その傾きすなわちΔＩ／Ｃより上記静電容量Ｃを求める。また、別途、対応する時点について上記式（１０）から内部インピーダンスの絶対値｜Ｚ｜を算出し、このようにして求めた上記静電容量Ｃ及び絶対値｜Ｚ｜の値を用い、上記式（５）から上記抵抗値Ｒが得られる。なお、ΔＩは、充電側の値を正とし、放電側の値を負とする。また、上記のように、印加開始後の時間ｔの逆数１／ｔを角速度ωとする。このようにして求めた上記抵抗値Ｒ及び上記静電容量Ｃの値により、上記式（２）及び上記式（３）を用いて各時点ｔにおける内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryが得られる。
【００５２】
（算出法４）
第１過程において充放電を行ない、充放電開始後の電圧を測定する。
また、第２過程では、式（７）を用いて各時点における上記抵抗値Ｒを上記静電容量Ｃの関数として得る。また、別途、対応する時点について上記式（１０）から内部インピーダンスの絶対値｜Ｚ｜を算出し、このようにして求めた上記抵抗値Ｒの関数及び絶対値｜Ｚ｜の値を用い、上記式（５）から上記静電容量Ｃを得る。ここで、充放電時電流をΔＩ、充放電開始直前の電圧と充放電開始後の電圧の差をΔＥとする。ΔＩは、充電側の値を正とし、放電側の値を負とする。また、上記のように、印加開始後の時間ｔの逆数１／ｔを角速度ωとする。さらに、この静電容量Ｃと先の式（７）とから抵抗値Ｒを得る。このようにして求めた上記抵抗値Ｒ及び上記静電容量Ｃの値により、上記式（２）及び上記式（３）を用いて各時点ｔにおける内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryが得られる。
なお、この算出法４は、定電流以外の充放電時、例えば、内部インピーダンスを求めるための電流の負荷の矩形波状（あるいは階段状波状）の制御が不安定であったり、電流ではなく電圧側を制御するような場合で、それぞれ電流値が負荷中に大きく変化してしまうような場合などに用いて好適である。
【００５３】
［３．第３過程］
第３過程では、上記第１のパラメータ（即ち、抵抗値Ｒ）及び第２のパラメータ（即ち、静電容量Ｃ）を用いて、上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。具体的には、第２過程で求めた上記抵抗値Ｒ及び上記静電容量Ｃの値により、上記式（２）及び上記式（３）を用いて各時点ｔにおける内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryを得て、その実数成分ＺREAL及び虚数成分Ｚimaginaryに対してナイキストプロットを行なう。
【００５４】
各時点における実数成分ＺREAL及び虚数成分Ｚimaginaryに対してナイキストプロットを行なうことにより、例えば図３に示すような曲線が得られる。なお、図３は横軸を実数成分ＺREALとし、縦軸を虚数成分Ｚimaginaryとしたナイキストプロットの一例であり、白抜きの円が各時点の測定プロットを表わし、それを結ぶ実線が周波数特性を示す曲線を表わす。
以上のようにして、本実施形態の周波数特性評価方法が目的としていた電力貯蔵供給デバイスの内部インピーダンスの周波数特性を、上記のような曲線として得ることができる。なお、通常は、上記周波数特性を表わす曲線はアーク状の曲線として得られることが多い。
【００５５】
また、得られた周波数特性を表わす曲線から、常法を用いて、上記曲線に対する上記抵抗成分の抵抗値Ｒや上記キャパシタの静電容量Ｃを求めることも可能である。このような処理により、例えば、電力貯蔵供給デバイスがリチウム電池である場合を例に挙げると、得られた曲線の横軸との交点（直列抵抗成分）からはセパレータ部や端子部などの抵抗値が得られ、上記曲線の円弧状部分からは電荷移動反応の抵抗値や固液界面での電気二重層容量が得られる。
【００５６】
［４．測定装置］
上述したような本実施形態の周波数特性評価方法は、全てを人の手により行なうようにしてもよいが、その操作の一部又は全部を、何らかの装置を用いて行なうようにすることが好ましい。例えば、図４に示すような周波数特性評価装置によって行なうこともできる。
図４に示す周波数特性評価装置は、充放電部１と、測定部２と、計算部３とを備えていて、さらに、計算部３は、記憶部４と、測定値取込部５と、演算部６と、周波数特性評価部７と、制御部８とを備えている。なお、図４においては、データの動きは実線の矢印で示し、制御部８の制御信号は破線の矢印で示してある。
【００５７】
充放電部１は、評価対象である電力貯蔵供給デバイス１００に充電又は放電を行なわせるものである。したがって、充放電部１により、例えば定電流や定電圧などの任意の条件において電力貯蔵供給デバイス１００に充放電を行なわせることができるようになっている。
【００５８】
また、測定部２は、充放電部１により電力貯蔵供給デバイス１００が充放電を開始させられた後、上記充放電開始後の充放電の電流又は電圧を測定するものである。電流値又は電圧値を常時測定するようにしても良いが、ここでは、測定部２は測定制御部（図示省略）を備えていて、この測定制御部の制御によって上記の電流値又は電圧値を所定時間間隔で測定するようになっているものとする。詳しくは、測定部２は、電流値を測定する電流値測定部２１と、電圧値を測定する電圧値測定部２２とを備えていて、測定制御部の制御に従って、電力貯蔵供給デバイス１００の充放電の電流値及び電圧値を両方とも測定できるようになっている。また、電流値測定部２１で測定された各時点の電流値及び電圧値は、記憶部４に送られるようになっている。
【００５９】
さらに、記録部４は、測定部２（即ち、電流値測定部２１及び電圧値測定部２２）から送られてきた各時点における電流値及び電圧値を記録する部分である。記録された電流値及び電圧値は、測定値取込部５によって読み出され、第１及び第２のパラメータの演算に用いられることになる。
【００６０】
また、測定値取込部５は、制御部８の制御により、記録部４に記録された各時点の電流値及び電圧値を取り込んで演算部６に送る部分である。具体的には、測定部取込部５は、記録部４から電流値を取り込んで演算部６に送る電流値取込部５１と、記録部４から電圧値を取り込んで演算部６に送る電圧値取込部５２とを備えていて、この電流値取込部５１及び電圧値取込部５２が、制御部８の制御にしたがって、所定時間間隔で測定された各時点の電流値及び電圧値の取り込みをそれぞれ行なうようになっている。なお、測定部２が電流値又は電圧値を常時測定するようになっている場合には、測定した電流値又は電圧値の中から所定時間間隔の各時点の測定データを取り込むようにすれば良い。
【００６１】
さらに、演算部６は、制御部８の制御により、測定値取込部５（即ち、電流値取込部５１及び電圧値取込部５２）を介して記録部４から取り込まれた電流値及び電圧値から、第１パラメータ（即ち、上記抵抗値Ｒ）と第２パラメータ（即ち、上記静電容量Ｃ）とを演算して求める部分である。具体的には、演算部６は、第１パラメータを演算する第１パラメータ演算部６１と、第２パラメータを演算する第２パラメータ演算部６２とを備えていて、第１パラメータ演算部６１及び第２パラメータ演算部６２が、上述した式（１）〜式（１０）などを用いて、制御部８の制御により、それぞれ第１及び第２のパラメータを求めるようになっている。また、求められた第１及び第２のパラメータは、周波数特性評価部７に送られるようになっている。
【００６２】
また、周波数特性評価部７は、制御部８の制御により、演算部６（即ち、第１パラメータ演算部６１及び第２パラメータ演算部６２）から送られた第１及び第２のパラメータから電力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価する部分である。具体的には、送られた各時点の第１及び第２のパラメータから内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryを得て、これらをナイキストプロットするようになっている。
【００６３】
さらに、制御部８は、測定値取込部５、演算部６及び周波数特性評価部７の動作の制御を行なうものである。
本例の周波数特性評価装置は以上のように構成されているので、使用時には、まず、充放電部１及び測定部２を電力貯蔵供給デバイス１００に接続し、充放電部１により充放電を開始して、測定部２により充放電開始後の電流値及び電圧値を所定時間間隔で測定する（第１過程）。測定された各時点の電流値及び電圧値は、記憶部４に記憶される。
【００６４】
次に、制御部８は、測定値取込部５を制御して、記憶部４に記憶された各時点における電流値及び電圧値を取り込ませ、演算部６に送らせる。さらに、制御部８は、演算部６を制御して、取り込んだ電流値及び電圧値から各時点それぞれの第１及び第２のパラメータを演算させる（第２過程）。演算結果は、周波数特性評価部７に送られる。
【００６５】
また、制御部８は、周波数特性評価部７を制御して、送られた各時点の第１及び第２のパラメータそれぞれについて内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryを計算させ、それらをナイキストプロットして、電周波数特性を評価させる（第３過程）。これにより、本例の周波数特性評価装置は、力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価することができる。
なお、得られた評価結果は図示しない出力部から出力するようにしてもよいが、ここでは、記憶部４に保存しておくものとする。
【００６６】
以上のように、ここで例示した本例の周波数特性評価装置を用いれば、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価することができる。従来は、充電あるいは放電の進行中において、充放電を停止することなく、矩形波状又はステップ状の電流や電圧を印加しながら、電力貯蔵供給デバイスの内部インピーダンスの過渡的な変化の評価することは行われていた。しかし、従来の技術はある特定の周波数におけるインピーダンス評価に留まっており、現実の充放電状態において任意の周波数における周波数特性を、短時間かつ簡便に評価することはできていなかった。具体的には、例えば直接測定法では、印加する交流の周波数を振りながら、各周波数でのインピーダンスを評価するという原理上、短時間で評価できず、また、例えばフーリエ変換法では、過電圧に対してインピーダンスが線形であるという前提が必要であったため、現実の充放電状態においてインピーダンスの周波数特性を評価することはできなかった。したがって、上記のように速やか且つ簡便に実際の充放電状態において上記周波数特性を評価できることは、電力貯蔵供給デバイスの評価などを行なう上で、極めて大きな利点である。
【００６７】
また、この周波数特性評価装置は、例えば、過電圧が大きい場合の充放電（高電流や低温化）での内部インピーダンスの周波数特性評価に用いることができる。充放電中の電流電圧そのもので評価を行なうことができ、内部インピーダンスの過電圧に対する非線形性を気にせず評価することができるためである。
【００６８】
ただし、ここで例示した周波数特性評価装置は、あくまで本実施形態の周波数特性評価方法を実現するための装置の一例であり、この装置は、本発明の要旨を逸脱しない範囲において任意に変形して実施することができる。
例えば、計算部３は、充放電部１及び測定部２と別体に形成し、例えば充放電部１及び測定部２を計算部３に対して着脱可能に形成してもよい。
また、例えば、充放電部１は、本例のように測定部２と別体とせず、一体に構成しても良い。
【００６９】
さらに、例えば、本例では電流値測定部２１と電圧値測定部２２とを測定部２に一体に構成したが、両者は別体に構成しても良い。また、上述したように、定電流または定電圧で充放電を行なう場合などにおいては、電流及び電圧のうち一方のみを測定するようにしても良い。
【００７０】
また、例えば、本例では充放電時の電流値及び電圧値の両方を測定するようにしたが、記録部に設定値や推定値などを記録させておき、これを測定値取込部５に取り込ませるようにしてもよい。例えば、定電流で充放電を行なう場合には設定した電流値を記録部４に記録させ、また、定電圧で充放電を行なう場合には設定した電圧値を記録部４に記録させ、これらの設定した電流値又は電圧値を測定値取込部５が取り込むように構成してもよい。
【００７１】
さらに、例えば、本例では電流値及び電圧値の両方を取込部５が記憶部４から取り込むものとして説明したが、定電流や定電圧で充放電を行なっている場合などのように、電流値及び電圧値のうちのいずれか一方が定数となる場合には、その定数となる値については取込部５による取り込みを行なわず、予め演算部６内で用いる演算の式に定数として盛り込んで第１及び第２のパラメータの演算を行なうようにしてもよい。
【００７２】
また、充放電時の電力貯蔵供給デバイス１００の温度を所定の温度に設定するため、恒温槽などの温度調整装置を取り付けて併用するようにしても良い。もちろん、上記周波数特性評価装置に温度調整部を設け、周波数特性評価装置自体が電力貯蔵供給デバイス１００の温度を制御できるようにしてもよい。
【００７３】
さらに、制御部８に、充放電部１の充放電や測定部２の測定時点を制御する機能を備えさせても良い。この場合、充放電部１の充放電を制御するには、制御部８には電力貯蔵供給デバイス１００に印加する電流又は電圧の大きさなどを制御する機能を備えさせれば良い。また、測定部２の測定時点を制御するために、制御部８に測定制御部の機能を備えさせて、制御部８が電力貯蔵供給デバイス１００の端子間の電流値又は電圧値を測定する上記の所定間隔を制御できるようにしても良い。
【００７４】
ところで、上述の周波数特性評価装置はコンピュータを用いて実現することも可能である。具体的には、計算部３や、計算部３の主要計算部３１、或いは、充放電部１の充放電や測定部２の測定時点を制御する制御装置（図示せず）などをコンピュータを用いて実現することができる。このように、コンピュータに上記周波数特性評価装置の各機能を実現させるためのプログラム（周波数特性評価プログラム）は、コンピュータ読取可能な記録媒体に格納することができ、記録媒体に格納した状態で頒布し、市場を流通させ、また売買することができる。
【００７５】
そして、このような記録媒体に格納されている周波数特性評価プログラムをコンピュータにインストールし、コンピュータがこの周波数特性評価プログラムを実行することで、上述した周波数特性評価装置を実現することができるため、記録媒体に格納されている周波数特性評価プログラムによって上述した周波数特性評価装置によって得られる効果と同様の効果が得られることになる。
【００７６】
ここで、記録媒体には、例えば半導体メモリなどのメモリ、磁気ディスク、光ディスク（例えばＣＤ−ＲＯＭ等）、光磁気ディスク（ＭＯ）、磁気テープ、ハードディスク、フレキシブルディスク、ＩＣカード、ＲＯＭカートリッジ、パンチカード、コンピュータ内部の記憶装置（ＲＡＭやＲＯＭなどのメモリ）、外部記憶装置等のプログラムを記録することができるものが含まれる。また、バーコードなどの符号が印刷された印刷物等のコンピュータ読取可能な種々の媒体を利用することもできる。なお、ＣＤ−ＲＯＭ、ＭＯ、磁気テープ、ＩＣカード等を可搬型記録媒体ともいう。
【００７７】
なお、ここでは、周波数特性評価プログラムをコンピュータ読取可能な記録媒体に格納しているが、記録媒体に格納しなくても良い。また、周波数特性評価プログラムは、上述のような記録媒体に格納せずに、例えば伝送媒体としてのネットワーク（通信ネットワーク）を介して伝送（送受信）することもでき、これにより周波数特性評価プログラムを頒布し、市場を流通させ、また売買することもできる。例えば、周波数特性評価プログラムをウェブサーバ等にアップロードしたり、インターネット等のネットワークを経由して例えばウェブサーバ等からダウンロードしたりすることもできる。
【００７８】
［II．第２実施形態］
本発明の第２実施形態の周波数特性評価方法も、上記第１実施形態と同様、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価するためのものである。
また、第２実施形態としての周波数特性評価方法では、以下に説明する第１過程、第２過程及び第３過程を経て、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。
【００７９】
即ち、第１過程では、電力貯蔵供給デバイスに対して充電又は放電を行ない、少なくとも１回上記の充電又は放電を開始した後の電流又は電圧を階段状又は矩形波状に変化させた際に、電流値又は電圧値を所定時間間隔で測定する。
また、第２過程では、電流又は電圧を変化させた際の各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する。
さらに、第３過程では、上記第１のパラメータ及び第２のパラメータを用いて上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。
【００８０】
なお、本実施形態の周波数特性評価方法はどのような場合に適用してもよいが、例えば、充放電の電流値が大きい場合や、過充電時のように電池内部が変化しているような場合などに用いることができる。即ち、充放電の電流値が大きい場合や、過充電時のように電池内部が変化しているような場合には、内部インピーダンスの過電圧に対する非線形性が強く、また、充放電の進行とともに内部インピーダンスは大きく変化している。従来の方法では、このような条件での充放電中の内部インピーダンスの周波数特性を評価することは困難が伴なっていた。しかし、本実施形態の周波数特性評価方法では、充放電中のにわずかな電気的な刺激を加えるだけで内部インピーダンスの周波数特性の評価することができるので、このような条件下でも内部インピーダンスの周波数特性を簡便に評価することができる。
以下、各過程毎に詳細に説明する。
【００８１】
［１．第１過程］
第１過程では、上述したように、電力貯蔵供給デバイスに対して直流で充電又は放電を行ない、少なくとも１回、上記の充電又は放電を開始した後の電流又は電圧を、電流又は電圧を印加することにより階段状又は矩形波状に変化させた際に、電流値又は電圧値を所定時間間隔で測定する。
【００８２】
充電と放電とのいずれを行なうかは任意であり、充電のみを行なうようにしてもよく、放電のみを行なうようにしてもよく、充電及び放電の両方を組み合わせて行なうようにしてもよい。通常は、充電時の内部インピーダンスの周波数特性を評価したい場合には充電を行なうようにし、放電時の内部インピーダンスの周波数特性を評価したい場合には放電を行なうようにする。
【００８３】
また、第１実施形態で説明したのと同様に、充放電を行なう前に、電力貯蔵供給デバイスは、予め所望の充電深度及び所望の温度に調整しておくことが好ましい。
さらに、充電速度や放電速度、充放電を行なう時間などの充放電条件が任意であり、周波数特性を評価しようとする実際の環境等に応じて任意に設定すればよい点も、第１実施形態と同様である。
【００８４】
ただし、本実施形態では、電力貯蔵供給デバイスに対して直流で充放電を行ない、この充放電の開始後に、少なくとも１回、電流及び電圧の一方を階段状又は矩形波状に変化させ（以下適宜、これを「波形変化をさせる」と表現する）、上記の変化させた期間において、電流及び電圧の他方を所定時間間隔で測定する。電力貯蔵供給デバイスに充放電を行なっているときに電流及び電圧の一方を変化させると、その一方の変化に応答して電流及び電圧の他方が変動することになるが、この変動は、第１実施形態で説明した充放電開始後の電流又は電圧の変化と同様に、電力貯蔵供給デバイスの内部インピーダンスを反映した挙動を示すものである。したがって、上記変動を測定すれば、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。
【００８５】
即ち、例えば図５に示すように、パルス電圧を期間Ｔ1及び期間Ｔ2の２回だけ印加して電圧の波形を矩形波状に変化させたときには、その期間Ｔ1及び期間Ｔ2それぞれの電流値を所定時間間隔で測定する。また、例えば図６に示すように、パルス電流を期間Ｔ3及び期間Ｔ4の２回だけ印加して電流の波形を矩形波状に変化させたときには、その期間Ｔ3及び期間Ｔ4の電圧値を所定時間間隔で測定する。なお、図５及び図６はそれぞれ、本実施形態の第１過程における電流値及び電圧値を説明するためのグラフであり、電圧及び電流それぞれの横軸は、同じ時間軸を表わすものとする。
【００８６】
さらに、例えば図７に示すように、パルス電圧を印加して期間Ｔ5及び期間Ｔ6に分けて階段状に電圧の波形を変化させたときには、その期間Ｔ5及びＴ6それぞれの電流値を所定時間間隔で測定する。また、例えば図８に示すように、パルス電流を印加して期間Ｔ7及び期間Ｔ8に分けて階段状に電流の波形を変化させたときには、その期間Ｔ7及びＴ8それぞれの電圧値を所定時間間隔で測定する。なお、図７及び図８はそれぞれ、本実施形態の第１過程における電流値及び電圧値を説明するためのグラフであり、電圧及び電流それぞれの横軸は、同じ時間軸を表わすものとする。
【００８７】
波形変化をさせる場合の波形は、階段状及び矩形波状のいずれであってもよいが、充放電の進行に伴う内部インピーダンスの変化に対し、その周波数特性評価のための操作の影響を小さくするためには、矩形波状に変化させることが好ましい。なお、評価のための一連の操作において、波形変化を２回以上行なう場合には、変化させる波形は階段状波及び矩形波状のいずれか一方のみを行うようにしても良く、両者を組み合わせて（例えば交互に）行なうようにしてもよい。
【００８８】
また、波形変化をさせる際に、電流と電圧とのいずれを変化させるかは任意である。通常、電力貯蔵供給デバイスは定電流で充放電を行なうことが多いが、このように定電流で操作する場合は、波形変化として電流を変化させることが一般的に望ましい。一方、定電圧で操作する場合は、波形変化として電圧を変化させることが一般的に望ましい。また、評価中の電流及び電圧の測定において、通常は、電流の測定精度が高い場合は波形変化として電圧を変化させるほうが得られる内部インピーダンスの周波数特性の精度が上がり、電圧の測定精度が高い場合は波形変化として電流を変化させるほうが得られる内部インピーダンスの周波数特性の精度が上がる。なお、一連の評価において、電流及び電圧のいずれか一方のみを変化させるようにしても良く、両者を組み合わせて（例えば交互に）変化させるようにしても良い。
【００８９】
また、正及び負（即ち、電流及び電圧を加えるか、差し引くか）のどちらに変化させるかも任意である。なお、一連の評価において、正及び負のいずれか一方の方向のみに変化させるようにしても良く、正及び負の方向を組み合わせて（例えば交互に）変化させるようにしても良い。
【００９０】
さらに、電流又は電圧を変化させる場合、変化の回数に制限はなく１以上であれば任意の回数だけ変化さえればよいが、周波数特性評価結果の信頼性を高める観点から、通常は、２回以上が望ましい。また、それに加えて周波数特性の時間推移を評価するためには、適当な時間間隔を置いて波形を変化させることが望ましい。
【００９１】
また、変化の大きさにも制限はなく、周波数特性を評価しようとする実際の環境等に応じて任意に選択すればよい。
ただし、矩形波状に変化をさせる場合には、基底とする充放電（即ち、波形を変化させない状態の充放電）によって電力貯蔵供給デバイスに出入りする電流容量に対し、矩形波状に波形を変化させたことによる電流容量の変化の割合が、通常１０％以下、好ましくは５％以下、より好ましくは１％以下とすることが望ましい。これは、充放電中の内部インピーダンスの周波数特性の推移を評価するような場合、できるだけ、評価しようとしている電力貯蔵供給デバイスの状態を変化させないようにすることが好ましいためである。
【００９２】
一方、階段状に変化させる場合には、基底とする充放電によって電力貯蔵供給デバイスに出入りする電流容量に対し、階段状に波形を変化させたことによる、階段形状の全段の電流容量の変化の割合が、通常１０％以下、好ましくは５％以下、より好ましくは１％以下とすることが望ましい。また、基底とする充放電によって電力貯蔵供給デバイスに出入りする電流容量に対し、階段状に波形を変化させたことによる、階段形状の各段それぞれの電流容量の変化の割合は、通常５％以下、好ましくは２．５％以下、より好ましくは０．５％以下とすることが望ましい。これは、矩形波の場合と同様、充放電中の内部インピーダンスの周波数特性の推移を評価するような場合、できるだけ、評価しようとしている電力貯蔵供給デバイスの状態を変化させないようにすることが好ましいためである。
【００９３】
さらに、１回の変化の期間（即ち、図５〜図８の期間Ｔ1〜Ｔ8など）それぞれの長さも任意であり、電力貯蔵供給デバイスの種類などに応じて適宜設定することができる。
ただし、矩形波状に変化させる場合、変化の期間は、通常１０μｓ以上、好ましくは１００μｓ以上、より好ましくは１０００μｓ以上、また、通常１０００ｓ以下、好ましくは１００ｓ以下、より好ましくは１０ｓ以下が望ましい。これは、通常の化学電池等の電力貯蔵供給デバイスでは１ｋＨｚ〜０．００１Ｈｚあたりに特徴的な内部インピーダンスの周波数特性が出やすいことと、１回の変化の期間が長すぎると電力貯蔵供給デバイス内部の状態（リチウム濃度、温度など）が基底の充放電だけの状態（波形を変化させていない状態）から乖離してしまう虞があることとの兼ね合いのためである。
【００９４】
一方、階段状に変化させる場合、変化の期間は、階段状の各段の期間として、通常１０μｓ以上、好ましくは１００μｓ以上、より好ましくは１０００μｓ以上、また、通常５００ｓ以下、好ましくは５０ｓ以下、より好ましくは５ｓ以下が望ましい。これも、通常の化学電池等の電力貯蔵供給デバイスでは１ｋＨｚ〜０．００１Ｈｚあたりに特徴的な内部インピーダンスの周波数特性が出やすいことと、１回の変化の期間が長すぎると電力貯蔵供給デバイス内部の状態（リチウム濃度、温度など）が基底の充放電だけの状態（波形を変化させていない状態）から乖離してしまう虞があることとの兼ね合いのためである。
【００９５】
なお、電流又は電圧を２回以上変化させる場合には、各変化の期間は一定であっても良く、異なっていても良い。
また、特に階段状に変化させる場合、ひとつの階段状波内の各段の大きさも、一定であっても良く、異なっていても良い。
【００９６】
さらに、階段状に波形を変化させる場合、波形の変化を順に大きくするように変化させても良く、順に小さくするように変化させても良い。
また、階段状に波形を変化させる場合、ひとつの階段状波内の段の数は任意であるが、通常２以上、好ましくは３以上、また、通常５０以下、好ましくは２０以下が望ましい。
【００９７】
さらに、本実施形態の第１過程において、電流値と電圧値とのいずれを測定するかは、第２過程で行なう演算に応じて適当に選択すればよい。したがって、電流を矩形波状又は階段状に変化させた場合には少なくとも電圧を測定するようにすれば良く、電圧を矩形波状又は階段状に変化させた場合には少なくとも電流を測定するようにすれば良い。これは、電流や電圧を矩形波状又は階段状に変化させる際には、変化させたときの電流や電圧の設定値を周波数特性の評価に用いることができるためである。もちろん、電流値及び電圧値の両方を測定するようにしてもよい。この場合、周波数特性の評価結果の精度や信頼性をより高めることができる。
【００９８】
また、電流値や電圧値は、所定時間間隔で測定するが、この所定時間は、第１実施形態と同様に、周波数特性の評価ができる限り任意である。
さらに、測定開始から測定終了までの間に測定する測定データ点数に制限はなく、第１実施形態と同様に、周波数特性の評価ができる限り任意である。
また、測定開始から測定終了までの総測定時間も、第１実施形態と同様に、制限なく任意である。
さらに、測定する電流値や電圧値の精度は、第１実施形態と同様、高いほど好ましい。
【００９９】
［２．第２過程］
第２過程では、電流又は電圧を変化させた際の各時点における上記電流値又は電圧値の変化を用い、上記電力貯蔵供給デバイスを、周波数に応じて抵抗値が変化する抵抗成分と、周波数に応じて静電容量値が変化するキャパシタ成分とからなる等価回路で表わした場合に、上記抵抗成分の抵抗値に相当する第１パラメータ、及び、上記キャパシタ成分の静電容量値に相当する第２パラメータを演算する。なお、第１実施形態と同様に、第１過程において測定した一部の時点の測定結果についてだけ上記の演算を行なうようにしてもよいが、周波数特性評価の精度を高める観点からは、測定した全ての時点の測定結果から上記演算を行なうことが好ましい。
【０１００】
上記のように、充放電開始後に電流及び電圧の一方を変化させた際の各時点における、電流及び電圧の他方の大きさ（電流値又は電圧値）は、電力貯蔵供給デバイスの内部インピーダンスを反映した挙動を示す。したがって、第１過程で測定した電流値及び電圧値の時間挙動を解析すれば、電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価することができる。
【０１０１】
具体的には、例えば、第１実施形態で説明した式（１）〜式（１０）を用い、第１実施形態と同様に演算を行なうことができる。即ち、式（１）〜式（１０）を用いれば、電力貯蔵供給デバイスに電流又は電圧の印加を開始した後に電流及び電圧の一方を変化させ、この変化させた際の他方の値（電流値又は電圧値）の変化を用いて、波形変化を開始した時点からの経過時間に対する抵抗値Ｒ及び静電容量Ｃの変化が得られ、さらに、この抵抗値Ｒ及び静電容量Ｃの値より波形変化を開始した時点からの経過時間に対するインピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryの変化が得られる。
【０１０２】
ただし、本実施形態では、式（６）〜式（１０）において、電流又は電圧印加後の電圧変化ΔＥ及び電流変化ΔＩは、具体的には、基底となっている電圧及び電流からの電圧変化及び電流変化を表わす。例えば、電流の波形を矩形波状又は階段状に変化させた場合には、電圧変化ΔＥは電流を矩形状又は階段状に変化させる直前からの電圧変化を表わし、電流変化ΔＩは、充放電中に波形を変化させた期間の、当該矩形波状又は階段状の電流値（つまり、基底となっている電流と、矩形波状又は階段状に変化させられたパルス負荷時の電流値との差）を表わす。また、例えば、電圧の波形を矩形波状又は階段状に変化させた場合には、電流変化ΔＩは電圧を矩形状又は階段状に変化させる直前からの電流変化を表わし、電圧変化ΔＥは、充放電中に波形を変化させた期間の、当該矩形波状又は階段状の電圧値（つまり、基底となっている電圧と、矩形波状又は階段状に変化させられたパルス負荷時の電圧値との差）を表わす。
さらに、本実施形態では、式（６）〜式（９）において、ｔは電流又は電圧の波形を矩形波状又は階段状に変化させた開始時点からの経過時間を表わす。
【０１０３】
また、本実施形態においても、第１過程で充放電を行なう際に、電力貯蔵供給デバイスに直流の電流又は電圧を印加することにより充放電を行なっているので、第１実施形態と同様に、波形変化を開始した時点、即ち、波形を変化させ始めた時点からの経過時間ｔに対する実数成分ＺREAL及び虚数成分Ｚimaginaryの変化の関係を、交流の電流又は電圧を印加した際の交流の角速度ωに対する実数成分ＺREAL及び虚数成分Ｚimaginaryの変化の関係に変換するには、印加開始後の時間ｔの逆数１／ｔを角速度ωと見なせば良い。また、角速度ωは、さらに必要に応じて周波数に変えることもできる。
【０１０４】
このように、第２過程においては、例えば上述した式（１）〜式（１０）を用いて、具体的に波形変化を開始した時点からの経過時間に対する上記抵抗値Ｒ及び上記静電容量Ｃの推移を得る。即ち、第１過程において電流又は電圧を測定した時点等の、サンプリングした各時点における上記抵抗値Ｒ及び上記静電容量Ｃを求める。このとき、上記の式（１）〜式（１０）の使い方には種々の組み合わせが考えられるが、各時点における上記抵抗値Ｒ及び上記静電容量Ｃを求めることができる限りその使い方は任意である。具体的には、例えば、ΔＩを波形変化を開始する直前からの電流変化とし、ΔＥを波形変化を開始する直前からの電圧変化として、第１実施形態で上述した算出法１〜算出法４と同様の操作により波形変化をさせている期間の各時点における上記の抵抗値Ｒ（第１パラメータ）及び静電容量Ｃ（第２パラメータ）を求めることにより、内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryが得られる。
【０１０５】
即ち、算出法１では、第１過程において、電流を矩形波状又は階段状に変化させ、変化させた期間の電圧値を測定する。そして、測定した電圧値に対して式（８）を用いて上記抵抗値Ｒ及び上記静電容量Ｃを求める。
また、算出法２では、第１過程において、電流を矩形波状又は階段状に変化させ、変化させた期間の電圧値を測定する。そして、測定した電圧値に対して式（８）を用いて上記抵抗値Ｒを求め、さらに、式（５）及び式（１０）を用いて上記静電容量Ｃを求める。
【０１０６】
さらに、算出法３では、第１過程において、電流を矩形波状又は階段状に変化させ、変化させた期間の電圧値を測定する。そして、測定した電圧値に対して式（８）を用いて上記静電容量Ｃを求め、さらに、式（５）及び式（１０）を用いて上記抵抗値Ｒを求める。
また、算出法４では、第１過程において、電流及び電圧の一方を矩形波状又は階段状に変化させ、変化させた期間の電流及び電圧の他方の値を測定する。そして、式（７）を用いて抵抗値Ｒを静電容量Ｃの関数として得た後、式（５）及び式（１０）を用いて静電容量Ｃを求めて、更に式（７）を用いて抵抗値Ｒを求める。
【０１０７】
さらに、上記のように、算出法１〜算出法４では、波形変化を開始した時点からの経過時間ｔの逆数１／ｔを角速度ωとする。
ただし、第１実施形態とは異なり、本実施形態の周波数評価方法ではパルス電流又はパルス電圧を印加してベースライン（即ち、波形を変化させない状態の電流及び電圧の波形）から波形を変化させたときの応答を用いるため、電流や電圧のベースラインが変化することを考慮して、電圧変化ΔＥ及び電流変化ΔＩを算出する際に、そのベースラインの変化を考慮した補正を行なうと、周波数特性評価の精度をより高めることができ、好ましい。具体的には、電圧変化ΔＥ及び電流変化ΔＩを算出する際に、上記のベースラインの変化量を差し引くことが好ましい。
【０１０８】
［３．第３過程］
第３過程では、第１実施形態と同様に、上記第１のパラメータ（即ち、抵抗値Ｒ）及び第２のパラメータ（即ち、静電容量Ｃ）を用いて、上記電力貯蔵供給デバイスの内部インピーダンスの周波数特性を評価する。具体的には、第２過程で求めた上記抵抗値Ｒ及び上記静電容量Ｃの値により、上記式（２）及び上記式（３）を用いて各時点ｔにおける内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryを得て、その実数成分ＺREAL及び虚数成分Ｚimaginaryに対してナイキストプロットを行なう。
【０１０９】
以上のようにして、本実施形態の周波数特性評価方法が目的としていた電力貯蔵供給デバイスの内部インピーダンスの周波数特性を、上記のような曲線として得ることができる。
また、第１実施形態と同様に、得られた周波数特性を表わす曲線から、常法を用いて、上記曲線に対する上記抵抗成分の抵抗値Ｒや上記キャパシタの静電容量Ｃを求めることも可能であり、リチウム電池のセパレータ部や端子部などの抵抗値、並びに、電荷移動反応の抵抗値や固液界面での電気二重層容量などが得られることも同様である。
【０１１０】
［４．測定装置］
上述したような本実施形態の周波数特性評価方法は、全てを人の手により行なうようにしてもよいが、その操作の一部又は全部を、何らかの装置を用いて行なうようにすることが好ましい。例えば、図９に示すような周波数特性評価装置によって行なうこともできる。
図９に示す周波数特性評価装置は、第１実施形態で説明した周波数特性評価装置にパルス印加部９を備えさせたもので、充放電部１と、測定部２と、計算部３と、パルス印加部９とを備えていて、さらに、計算部３は、記憶部４と、測定値取込部５と、演算部６と、周波数特性評価部７と、制御部８とを備えている。なお、図９において、図４と同じ符号で示す部分は図４と同様のものを表わし、さらに、データの動きは実線の矢印で示し、また、制御部８の制御信号は破線の矢印で示してある。
【０１１１】
充放電部１は、第１実施形態の周波数特性評価装置と同様のものである。
また、パルス印加部９は、充放電部１が充放電を開始した後、電流又は電圧にパルス電流又はパルス電圧を印加して、電流及び電圧の一方の波形を階段状又は矩形波状に変化させるものである。印加するパルスの回数、大きさ、印加する時間、変化させる波形、電流及び電圧のどちらに印加するかなどは、適宜設定できるようになっている。
【０１１２】
また、測定部２は、第１実施形態の周波数特性評価装置と同様のものであり、充放電部１により電力貯蔵供給デバイス１００が充放電を開始させられた後、測定制御部（図示省略）の制御に従って、上記充放電開始後に波形を変化させた期間（即ち、電流及び電圧の一方の波形を変化させている期間）の電流又は電圧を所定時間間隔で測定するものである。したがって、パルス印加部９によって電流又は電圧の波形が変化させられた期間の電流値及び電圧値についても、所定時間間隔で測定できるようになっている。なお、本実施形態の測定部２も、第１実施形態と同様に、電流値又は電圧値を常時測定するようにしても良く、その場合、測定値取込部５が測定した電流値又は電圧値の中から波形を変化させた期間の所定時間間隔の各時点の測定データを取り込むようにすれば良い。
【０１１３】
さらに、記録部４、測定値取込部５、演算部６、周波数特性評価部７及び制御部８は、周波数特性評価のために充放電開始後にパルス印加部９が電流又は電圧の波形を変化させた期間の測定結果を用いるようになっている他は、いずれも第１実施形態同様になっている。したがって、測定部２（即ち、電流値測定部２１及び電圧値測定部２２）から送られてきた、パルス印加部９が電流又は電圧の波形を変化させた期間の各時点における電流値及び電圧値は、記録部４に記録され、その電流値及び電圧値は、制御部８の制御により、測定値取込部５（即ち、電流値取込部５１及び電圧値取込部５２）に取り込まれて演算部６に送られるようになっている。また、演算部６（即ち、第１パラメータ演算部６１及び第２パラメータ演算部６２）は、制御部８の制御により、送られてきた電流値及び電圧値から第１パラメータ（即ち、上記抵抗値Ｒ）と第２パラメータ（即ち、上記静電容量Ｃ）とを演算し、周波数特性評価部７に送るようになっている。さらに、周波数特性評価部７は、制御部８の制御により、送られてきた第１及び第２のパラメータから内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryを得て、これらをナイキストプロットして、電力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価するようになっている。
【０１１４】
本例の周波数特性評価装置は以上のように構成されているので、使用時には、まず、充放電部１及び測定部２を電力貯蔵供給デバイス１００に接続し、充放電部１により充放電を開始する。その後、パルス印加部９から電流又は電圧にパルスを印加して、電流及び電圧の一方の波形を階段状又は矩形波状に変化させる。この間、測定部２により充放電開始後の電流値及び電圧値を所定時間間隔で測定するようにする（第１過程）。測定された各時点の電流値及び電圧値は、記憶部４に記憶される。
【０１１５】
次に、制御部８は、測定値取込部５を制御して、記憶部４に記憶された各時点における電流値及び電圧値のうち、波形変化をさせている期間の電流値及び電圧値を取り込ませ、演算部６に送らせる。さらに、制御部８は、演算部６を制御して、取り込んだ電流値及び電圧値から各時点それぞれの第１及び第２のパラメータを演算させる（第２過程）。演算結果は、周波数特性評価部７に送られる。
【０１１６】
また、制御部８は、周波数特性評価部７を制御して、送られた各時点の第１及び第２のパラメータそれぞれについて内部インピーダンスＺの実数成分ＺREAL及び虚数成分Ｚimaginaryを計算させ、それらをナイキストプロットして、電周波数特性を評価させる（第３過程）。これにより、本例の周波数特性評価装置は、力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価することができる。
なお、得られた評価結果は図示しない出力部から出力するようにしてもよいが、ここでは、記憶部４に保存しておくものとする。
【０１１７】
以上のように、ここで例示した本例の周波数特性評価装置を用いれば、第１実施形態の周波数特性評価装置と同様に、速やか且つ簡便に、実際の充放電状態において、電力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価することができる。
また、この周波数特性評価装置は、充放電中の電力貯蔵供給デバイスの内部インピーダンスの周波数特性の時間推移を簡単に評価することもできる。例えば、電力貯蔵供給デバイスに高電流で過充電をしているような場合、具体例としては、過電圧が大きく、且つ、電力貯蔵供給デバイスが徐々に壊れていくような場合に、上記周波数特性の時間推移を評価しようとする際には、本実施形態の周波数特性評価装置は、非常に有効である。なお、充放電開始直後に電流値又は電圧値が比較的変化が小さい状態の電力貯蔵供給デバイス１００の内部インピーダンスの周波数特性を評価する場合にも本例の周波数特性評価装置は用いて好適である。
【０１１８】
ただし、ここで例示した周波数特性評価装置は、あくまで本実施形態の周波数特性評価方法を実現するための装置の一例であり、この装置は、第１実施形態で説明した周波数特性評価装置と同様、本発明の要旨を逸脱しない範囲において任意に変形して実施することができる。
また、さらに、例えばパルス印加部９は、充放電部１、測定部２、計算部３などと任意の組み合わせで一体に構成しても良く、また、着脱可能に別体に構成しても良い。なお、パルス印加部９は、例えば、充放電部１と並列に配設された充放電器などにより構成することができる。
【０１１９】
さらに、本実施形態の周波数特性評価装置も、第１実施形態の周波数特性評価装置と同様、コンピュータを用いて実現することが可能であり、また、コンピュータに上記周波数特性評価装置の各機能を実現させるためのプログラム（周波数特性評価プログラム）についても、第１実施形態と同様に実現することができる。
【０１２０】
［III．第３実施形態］
第１及び第２実施形態で説明したような電力貯蔵供給デバイスの内部インピーダンスの周波数特性評価方法を利用して、電力貯蔵供給デバイスの保護装置を作製することができる。例えば、第１及び第２実施形態で説明した方法により評価した周波数特性に基づき、評価対象である電力貯蔵供給デバイスの性能を評価するデバイス性能評価部と、このデバイス性能評価部での評価結果に基づいて上記電力貯蔵供給デバイスの充放電の制御を行なう充放電制御部とを備えた保護装置によれば、上記のように実際の使用状態に基づいて評価した電力貯蔵供給デバイスの周波数特性により、その充放電を確実に制御することが可能となる。また、上記のデバイス性能評価部は、具体的には、例えば第１又は第２実施形態で説明した周波数特性評価装置で評価された周波数特性（上記の例では、記憶部４に記憶されたもの）に基づいて電力貯蔵供給デバイスの性能を評価するように構成することができる。
【０１２１】
保護装置の具体例を挙げると、図１０のようなものが挙げられる。図１０は、第１実施形態で説明した周波数特性評価装置に、保護装置１０を備えさせた構成を表わす。なお、図１０において、図４と同様の符号で示す部分は、図４と同様のものを表わす。
保護装置１０は、図１０に示すように、デバイス性能評価部１１と、充放電制御部１２とを備えている。なお、保護装置は保護回路を備えていてもよい。
【０１２２】
デバイス性能評価部１１は、記憶部４に記憶された電力貯蔵供給デバイス１００の内部インピーダンスの評価結果を読み込み、この評価結果に基づいて電力貯蔵供給デバイスの性能を評価する部分である。本例では、記憶部４から読み込んだ評価結果から得られる直列抵抗成分Ｒが、電力貯蔵供給デバイス１００の初期の状態或いは正常の状態に対して、所定の倍率の大きさ以上である場合に過充電と判定するようになっている。上記の所定の倍率は、例えば上記第２実施形態の方法で予め過充電時の内部インピーダンスの周波数特性の時間推移を評価しておき、それに基づいて決定することができる。また、直列抵抗成分をある特定の周波数での内部インピーダンスを用いて求めることができる場合は、その周波数における内部インピーダンスを用いて直列抵抗成分を求めるようにしてもよい。ここでは、周波数１０００Ｈｚでの内部インピーダンスから得られる直列抵抗成分の値が性常時の３倍以上となった場合に、電力貯蔵供給デバイス１００が過充電状態であると判定するようになっているものとするが、この設定は適宜変更しても良い。
【０１２３】
また、充放電制御部１２は、デバイス性能評価部１１の評価結果に基づき、充放電部１の動作を制御する部分である。本例では、デバイス性能評価部１１において電力貯蔵供給デバイス１００が過充電状態であると判定された場合に、充放電部１を制御して充放電操作を停止するようになっている。
【０１２４】
以上のように、周波数特性評価装置に保護装置１０を備えるようにすれば、電力貯蔵供給デバイス１００が過充電状態になった場合でも、内部インピーダンスの周波数特性に基づいてデバイス性能評価部１１が上記の過充電状態を感知し、充放電制御部１２が充放電を停止するので、過充電により電力貯蔵供給デバイス１００に問題が生じることを防止することができる。
【０１２５】
なお、図１０の例では周波数特性評価装置として第１実施形態で説明したものを用いたが、第２実施形態で説明したような周波数特性評価装置に上記の保護装置１０を適用するようにしても良い。第２実施形態で例示した周波数特性評価装置に上記の保護装置１０を適用した場合には、その構成は例えば図１１のようになる。なお、図１１において図９で用いたのと同様の符号で示す部分は、図９と同様のものを表わす。さらに、保護装置１０、デバイス性能評価部１１及び充放電制御部１２は、図１０で説明したものと同様のものである。
【実施例】
【０１２６】
以下、実施例を用いて本発明について更に具体的に説明するが、本発明は以下の実施例に限定されるものではなく、本発明の要旨を逸脱しない範囲において任意に変形して実施することができる。
【０１２７】
［実施例１］
（１）電流及び電圧の測定
電力貯蔵供給デバイスとして、外装がラミネート樹脂で形成された容量６ｍＡｈのリチウムイオン電池（縦８５ｍｍ、横４５ｍｍ、厚み０．４ｍｍ）を用意した。なお、用いたリチウムイオン電池の組成は、以下の表１のとおりである。
【０１２８】
【表１】
2006162283.tif 000011
【０１２９】
このリチウムイオン電池の電圧を、満充電である４．１Vに調整した。
次に、このリチウムイオン電池に対し、充放電機（ＫＥＩＴＨＬＥＹ社  ソースメータ２４４０）を接続し、この充放電機はパーソナルコンピュータに接続した。充放電機の動作の制御並びに電流および電圧の記録は、このパーソナルコンピュータで実施できるよう、プログラムを組んだ。
【０１３０】
次に、リチウムイオン電池を恒温槽に入れ、温度３０℃に保持した。
続いて、充放電器の放電電流を１２ｍＡに設定後、リチウムイオン電池の放電を定電流で開始した。この際、放電の開始から１秒後まで、測定時間間隔０．１ｍｓｅｃで電流とリチウムイオン電池の電圧とを測定した（第１過程）。
また、同様の測定を、放電電流を１８ｍＡとして同様の測定を行なった（第１過程）。
電圧の測定結果を図１２に示す。
【０１３１】
（２）内部インピーダンスの解析
図１２の測定結果より、前述の算出法１を用いて抵抗値Ｒ及び静電容量Ｃを求め（第２過程）、それからＺrealとＺimaginaryとを求めてナイキストプロットを行ない、リチウムイオン電池の内部インピーダンスの周波数特性を得た（第３過程）。得られたナイキストプロットを図１３に示す。
図１３には、実際に高電流で放電を開始した際の内部インピーダンスの示す曲線は、放電の進行にともなって過電圧が増加することの影響を受け、半円状から大きくひずんでいることが示されている。また、電流が大きいほうが、過電圧が大きいため、その影響で内部インピーダンスがより小さいことが示されている。
【０１３２】
［実施例２］
（１）電流及び電圧の測定
電力貯蔵供給デバイスとして、実施例１と同様のリチウムイオン電池を用意し、電圧を、満充電である４．１Vに調整した。
このリチウムイオン電池に対し、充放電機（ＫＥＩＴＨＬＥＹ社  ソースメータ２４４０）２台を並列に接続し、また、この充放電機をパーソナルコンピュータに接続した。充放電機の動作の制御並びに電流および電圧の記録は、このパーソナルコンピュータで実施できるよう、プログラムを組んだ。ここで、パーソナルコンピュータは、充放電機の１台（充放電機Ａ）が４２ｍＡの定電流でリチウムイオン電池に充電を行ない、さらに、もう一台（充放電機Ｂ：パルス印加部）が約２０秒ごとに１２ｍＡで１秒間幅の充電を行なってリチウムイオン電池を流れる電流の波形を矩形波状に変化させるように制御を行なうようにした。
【０１３３】
次いで、リチウムイオン電池を恒温槽に入れ、温度３０℃に保持した。
続いて、充放電機Ａの電流の設定を４２ｍＡに設定後、リチウムイオン電池の過充電を開始した。充放電機Ａによる充電を継続した状態で、約２０秒おきに１秒間だけ充放電機Ｂによって１２ｍＡの充電電流を加え、電流の波形を矩形波状に変化させた。この際、測定時間間隔０．１ｍｓｅｃで電流とリチウムイオン電池の電圧とを測定した（第１過程）。充放電機Ｂにより電流の波形を矩形波状に変化させた期間に測定された電圧の測定結果のうち、充放電機Ａによる充放電開始後２０秒、７４秒及び６３４秒後から波形を変化させた期間での測定結果を図１４に示す。
【０１３４】
（２）内部インピーダンスの解析
図１４の測定結果より、前述の算出法１を用いて抵抗値Ｒ及び静電容量Ｃを求め（第２過程）、それからＺrealとＺimaginaryとを求めてナイキストプロットを行ない、リチウムイオン電池の内部インピーダンスの周波数特性を得た（第３過程）。得られたナイキストプロットを図１５に示す。
過充電時においてはリチウムイオン電池内部は不安定で時間とともに変化していくが、本実施例の方法では従来方法のように過充電を停止することなく、過充電中開始後の内部インピーダンスの周波数特性の時間推移を評価することができる。例えば、図１５を見てわかるように、過充電開始直後は内部インピーダンスの変化は小さいが、過充電開始後に時間が長く経過すると、著しく内部インピーダンスが大きくなっていることがわかる。
【産業上の利用可能性】
【０１３５】
本発明は、電力貯蔵供給デバイスの内部インピーダンスの評価方法に関する。こうした内部インピーダンスの評価は、例えば、リチウムイオン電池などにおいて充電や放電における電流負荷の限界を評価するための手法や、その内部や外部に短絡が生じた際の、電池内の温度分布を、特定の過程を経て計算機シミュレーションし、精度よく評価する手法などに応用できる。また、リチウムイオン電池の過充電時や高温下への暴露時などのように、電池の内部インピーダンスが安定せず時間とともに急速に変化するような場合の内部インピーダンスの評価にも適する。
【図面の簡単な説明】
【０１３６】
【図１】本発明の第１及び第２実施形態について説明するもので、電力貯蔵供給デバイスの内部インピーダンスを表わす、抵抗とキャパシタとの直列接続により構成される等価回路の図である。
【図２】本発明の第１実施形態について説明するもので、（ａ）〜（ｃ）はいずれも、第１過程において定電流で放電を行なった場合の電圧の測定値をプロットした曲線の一例として示すグラフである。
【図３】本発明の第１及び第２実施形態について説明するもので、横軸を実数成分ＺREALとし、縦軸を虚数成分Ｚimaginaryとしたナイキストプロットの一例である。
【図４】本発明の第１実施形態としての周波数特性評価装置を表わすブロック図である。
【図５】本発明の第２実施形態の第１過程における電流値及び電圧値を説明するためのグラフである。
【図６】本発明の第２実施形態の第１過程における電流値及び電圧値を説明するためのグラフである。
【図７】本発明の第２実施形態の第１過程における電流値及び電圧値を説明するためのグラフである。
【図８】本発明の第２実施形態の第１過程における電流値及び電圧値を説明するためのグラフである。
【図９】本発明の第２実施形態としての周波数特性評価装置を表わすブロック図である。
【図１０】本発明の第３実施形態としての保護装置について説明するためのブロック図である。
【図１１】本発明の第３実施形態としての保護装置の変形例について説明するためのブロック図である。
【図１２】本発明の実施例１で測定したリチウムイオン電池の電圧の時間変化を表わすグラフである。
【図１３】本発明の実施例１で得られたリチウムイオン電池の内部インピーダンスの周波数特性を表わすナイキストプロットである。
【図１４】本発明の実施例２で測定したリチウムイオン電池の電圧の時間変化を表わすグラフである。
【図１５】本発明の実施例２で得られたリチウムイオン電池の内部インピーダンスの周波数特性を表わすナイキストプロットである。
【符号の説明】
【０１３７】
１  充放電部
２  測定部
３  計算部
４  記憶部
５  測定値取込部
６  演算部
７  周波数特性評価部
８  制御部
９  パルス印加部
１０  保護装置
１１  デバイス性能評価部
１２  充放電制御部
２１  電流値測定部
２２  電圧値測定部
３１  主要計算部
５１  電流値取込部
５２  電圧値取込部
６１  第１パラメータ演算部
６２  第２パラメータ演算部
【図１】
2006162283.tif 000012
【図２】
2006162283.tif 000013
【図３】
2006162283.tif 000014
【図４】
2006162283.tif 000015
【図５】
2006162283.tif 000016
【図６】
2006162283.tif 000017
【図７】
2006162283.tif 000018
【図８】
2006162283.tif 000019
【図９】
2006162283.tif 000020
【図１０】
2006162283.tif 000021
【図１１】
2006162283.tif 000022
【図１２】
2006162283.tif 000023
【図１３】
2006162283.tif 000024
【図１４】
2006162283.tif 000025
【図１５】
2006162283.tif 000026