(57)【要約】【課題】燃料電池の電気的特性の測定を高速且つ安定に行うことができる燃料電池の特性測定装置を提供する。【解決手段】所定の直流電流に所定の交流成分を重畳した電流を燃料電池から出力させる負荷制御部３０１と、燃料電池から異なる２値の出力電流を出力させた時の各出力電圧値を計測し、計測された出力電流２値の変化量Δｉと出力電圧２値の変化量Δｖを求め、その割り算値Δｖ／Δｉを、当該燃料電池の０Ｈｚ近傍でのインピーダンス値に決定する低周波数領域インピーダンス算出部３１０とを備える。【選択図】図１
2007265895.tif 000002
【特許請求の範囲】
【請求項１】
燃料電池の特性を測定する装置において、
所定の直流電流に所定の交流成分を重畳した電流を燃料電池から出力させる負荷制御手段と、
燃料電池から異なる２値の出力電流を出力させた時の各出力電圧値を計測する計測手段と、
計測された出力電流２値の変化量Δｉと出力電圧２値の変化量Δｖを求め、その割り算値Δｖ／Δｉを、当該燃料電池の０Ｈｚ近傍でのインピーダンス値に決定する決定手段と
を備えることを特徴とする燃料電池の特性測定装置。
【請求項２】
さらに、前記重畳する交流成分の周波数を変化させながら燃料電池のインピーダンス値を計測するインピーダンス計測手段と、
計測されたインピーダンス値に基づいて、燃料電池に対応する所定の等価回路の回路定数を求める演算手段とを備え、
その演算手段において、０Ｈｚ近傍のインピーダンス値として、インピーダンスの実数部Ｚ’を前記割り算値の絶対値｜Δｖ／Δｉ｜、虚数部Ｚ”を０とする拘束条件を適用する
ことを特徴とする請求項１記載の燃料電池の特性測定装置。
【請求項３】
前記インピーダンス計測手段が、前記交流成分の周波数を高い方から低い方へと変化させながらインピーダンスを計測し、
前記演算手段が、インピーダンス値の計測毎に等価回路の回路定数を求める演算を行い、
さらに、前記演算手段が求めた等価回路定数の変化量が設定値以下となるか否かの判定を行い、設定値以下となった場合に、インピーダンス計測を終了して、回路定数を決定する判定手段を備える
ことを特徴とする請求項２記載の燃料電池の特性測定装置。
【請求項４】
前記判定手段が、所定の周波数以下になっても前記等価回路定数の変化量が設定値以下とならないことや前記等価回路定数の変化量が増加傾向にあることを検知した場合には、その旨を出力する
ことを特徴とする請求項３記載の燃料電池の特性測定装置。
【請求項５】
燃料電池の特性を測定する方法において、
所定の直流電流に所定の交流成分を重畳した電流を燃料電池から出力させる負荷制御過程と、
燃料電池から異なる２値の出力電流を出力させた時の各出力電圧値を計測する計測過程と、
計測された出力電流２値の変化量Δｉと出力電圧２値の変化量Δｖを求め、その割り算値Δｖ／Δｉを、当該燃料電池の０Ｈｚ近傍でのインピーダンス値に決定する決定過程と
を含んでいることを特徴とする燃料電池の特性測定方法。
【請求項６】
燃料電池の特性を測定する装置において、
交流インピーダンスの測定前にその試験条件で開回路電圧を測定又は演算するとともに、電流−電圧特性を測定する手段と、
開回路電圧から電流−電圧特性を引算したものを電流変化に対して微分する手段と、
その結果をインピーダンス計測を行う電流値での当該燃料電池の０Ｈｚ近傍でのインピーダンス値として決定する手段と
を備えることを特徴とする燃料電池の特性測定装置。
【請求項７】
燃料電池の特性を測定する方法において、
交流インピーダンスの測定前にその試験条件で開回路電圧を測定又は演算するとともに、電流−電圧特性を測定する過程と、
開回路電圧から電流−電圧特性を引算したものを電流変化に対して微分する過程と、
その結果をインピーダンス計測を行う電流値での当該燃料電池の０Ｈｚ近傍でのインピーダンス値として決定する過程と
を含んでいることを特徴とする燃料電池の特性測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、燃料電池の電気的特性を測定する際に用いて最適な燃料電池の特性測定装置及び方法に関する。
【背景技術】
【０００２】
燃料電池の内部損失を測定するための手法の１つとして、損失をインピーダンスに置き換えて測定する交流インピーダンス法（交流法）がある。交流インピーダンス法では、動作状態の燃料電池に電子負荷装置やバイポーラ電源を用いて微小の交流電流を負荷として重畳させ、その際の電池電圧を観測することで各周波数における内部インピーダンスが測定される。電池の負荷電流に重畳する交流電流周波数を変化させながら測定を行い、燃料電池の等価回路となる内部インピーダンスを実数部と虚数部に分けて演算し、コール・コール・プロット図（Cole-Cole Plot）を描くことで、電池内部の特性が定量的に把握される。
【０００３】
例えば特許文献１に記載されている装置では、最初に、予め定められた測定周波数Ｆ１、Ｆ２及びＦ３の交流信号を燃料電池セルに対して印加してインピーダンスを測定する。複素平面上で得られたこれらインピーダンスに対し、抵抗値が夫々取得され、更にこれら抵抗値から、劣化診断要素が生成される。劣化診断要素は、周波数Ｆ１、Ｆ２及びＦ３に対する抵抗値を夫々Ｒ１、Ｒ２及びＲ３とすると、Ｒ１、Ｒ２−Ｒ１、及びＲ３−２Ｒ２＋Ｒ１として規定されており、この劣化診断要素の値を、異なる時刻に得られた劣化診断要素の値と比較することによって、燃料電池の劣化が診断される。
【０００４】
ここで図７を参照して、従来の燃料電池のインピーダンス計測手法について具体的に説明する。従来の燃料電池のインピーダンス計測では、燃料電池に対して、振幅Ｉａｃ、周波数ｆの交流電流を重畳させた直流電流値Ｉｄｃを印加し、そのときの、燃料電池の端子電圧の周波数応答性（ゲイン、フェーズ）を測定する。周波数を変えながら順次インピーダンス計測を行い、重畳交流電流の周波数に対する燃料電池端子電圧の応答特性を測定していた。また、燃料電池端子電圧の応答特性から、燃料電池の等価回路モデルへのフィッティングを行い、燃料電池を等価回路でモデル化を行っている。
【特許文献１】特開２００５−２８５６１４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上記周波数特性を測定する際には、一般的には高周波数から低周波数方向に周波数を落としながら、燃料電池端子電圧の応答性を測定していくが、等価回路モデルにフィッティングするためには、コンデンサ成分が無視できるような、非常に遅い、例えば１Ｈｚ以下での応答性を測定する必要がある。この低い周波数領域では、燃料電池の発電で発生する生成水の排出の周波数帯域と近くなる場合があり、生成水によるガス流路詰まりや圧力損失の増加、一時的な燃料不足、酸化剤不足等により、端子電圧が不安定となり、その結果インピーダンス計測の誤差、バラツキが大きくなり、再現性も悪くなる。このため、燃料電池の等価回路において、コンデンサと並列に配置された直流抵抗のフィッティングの精度が下がる原因となっていた。
【０００６】
また、上記、生成水による発電不安定要素がない場合であっても、１Ｈｚ以下の周波数までその応答特性を取得するためには、非常に長い試験時間がかかる。
【０００７】
本発明は、上記の事情に鑑みてなされたものであり、燃料電池の電気的特性の測定を高速且つ安定に行うことができる燃料電池の特性測定装置及び方法を提供することを目的とする。より具体的には、インピーダンス計測における低周波数のインピーダンス特性、即ち、等価回路におけるコンデンサ成分が無視でき、直流抵抗のみと考えられる領域でのインピーダンス特性を素早く、かつ、安定に取得することができる燃料電池の特性測定装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、請求項１記載の発明は、燃料電池の特性を測定する装置において、所定の直流電流に所定の交流成分を重畳した電流を燃料電池から出力させる負荷制御手段と、燃料電池から異なる２値の出力電流を出力させた時の各出力電圧値を計測する計測手段と、計測された出力電流２値の変化量Δｉと出力電圧２値の変化量Δｖを求め、その割り算値Δｖ／Δｉを、当該燃料電池の０Ｈｚ近傍でのインピーダンス値に決定する決定手段とを備えることを特徴とする。
【０００９】
請求項２記載の発明は、さらに、前記重畳する交流成分の周波数を変化させながら燃料電池のインピーダンス値を計測するインピーダンス計測手段と、計測されたインピーダンス値に基づいて、燃料電池に対応する所定の等価回路の回路定数を求める演算手段とを備え、その演算手段において、０Ｈｚ近傍のインピーダンス値として、インピーダンスの実数部Ｚ’を前記割り算値の絶対値｜Δｖ／Δｉ｜、虚数部Ｚ”を０とする拘束条件を適用することを特徴とする。
【００１０】
請求項３記載の発明は、前記インピーダンス計測手段が、前記交流成分の周波数を高い方から低い方へと変化させながらインピーダンスを計測し、前記演算手段が、インピーダンス値の計測毎に等価回路の回路定数を求める演算を行い、さらに、前記演算手段が求めた等価回路定数の変化量が設定値以下となるか否かの判定を行い、設定値以下となった場合に、インピーダンス計測を終了して、回路定数を決定する判定手段を備えることを特徴とする。
【００１１】
請求項４記載の発明は、前記判定手段が、所定の周波数以下になっても前記等価回路定数の変化量が設定値以下とならないことや前記等価回路定数の変化量が増加傾向にあることを検知した場合には、その旨を出力することを特徴とする。
【００１２】
請求項５記載の発明は、燃料電池の特性を測定する方法において、所定の直流電流に所定の交流成分を重畳した電流を燃料電池から出力させる負荷制御過程と、燃料電池から異なる２値の出力電流を出力させた時の各出力電圧値を計測する計測過程と、計測された出力電流２値の変化量Δｉと出力電圧２値の変化量Δｖを求め、その割り算値Δｖ／Δｉを、当該燃料電池の０Ｈｚ近傍でのインピーダンス値に決定する決定過程とを含んでいることを特徴とする。
【００１３】
請求項６記載の発明は、燃料電池の特性を測定する装置において、交流インピーダンスの測定前にその試験条件で開回路電圧を測定又は演算するとともに、電流−電圧特性を測定する手段と、開回路電圧から電流−電圧特性を引算したものを電流変化に対して微分する手段と、その結果をインピーダンス計測を行う電流値での当該燃料電池の０Ｈｚ近傍でのインピーダンス値として決定する手段とを備えることを特徴とする。
【００１４】
請求項７記載の発明は、燃料電池の特性を測定する方法において、交流インピーダンスの測定前にその試験条件で開回路電圧を測定又は演算するとともに、電流−電圧特性を測定する過程と、開回路電圧から電流−電圧特性を引算したものを電流変化に対して微分する過程と、その結果をインピーダンス計測を行う電流値での当該燃料電池の０Ｈｚ近傍でのインピーダンス値として決定する過程とを含んでいることを特徴とする。
【発明の効果】
【００１５】
本発明によれば、周波数ｆ≒０Ｈｚ付近（低周波数領域）において、交流信号に対する応答性を測定する必要がなく、直流成分のみでインピーダンス特性を測定できるため、測定の高速化が可能となる。また、交流信号に対する応答性を測定する必要がないので、燃料電池の発電で発生する生成水による一時的な流路閉塞、あるいは圧力損失の増加が起因となる燃料ガス不足、あるいは酸化ガス不足による燃料電池端子電圧の不安定性（インピーダンス計測でのバラツキ、再現性低下の要因）を除去した測定が可能である。つまり、本発明によれば、インピーダンス計測における低周波数のインピーダンス特性、即ち、等価回路におけるコンデンサ成分が無視でき、直流抵抗のみと考えられる領域でのインピーダンス特性を素早く、かつ、安定に取得することができる。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して本発明による燃料電池の特性測定装置の実施の形態について説明する。図１は、本実施の形態の構成を示すブロック図である。なお、図１では、燃料電池１００に対する燃料ガスおよび酸化ガスの供給系および排出系等は図示を省略している。
【００１７】
本実施の形態の特性測定装置は、燃料電池１００と電子負荷装置２００とインピーダンス計測器３００からなる。インピーダンス計測器３００は、内部にＣＰＵ（中央処理装置）、ＲＡＭ、ＲＯＭ、フラッシュメモリ、ハードディスク等の記憶装置や、各種入出力インターフェースを備え、フラッシュメモリ等に記録されている所定のプログラムを実行することで各部の制御や、外部装置との信号の入出力を行う。燃料電池１００と電子負荷装置２００と電流計測部３０２は、結線１１によって電気的に直列に接続されている。電圧計測部３０３は結線１２によって燃料電池１００に並列に接続され、その両端電圧（端子電圧）を測定する。電子負荷制御部３０１は、電子負荷装置２００に対して、直流電流値と、重畳交流電流の振幅と周波数を指示・設定し、燃料電池１００に対する負荷の制御を行う。
【００１８】
インピーダンス計測器３００においてインピーダンスの計測は、低周波数領域と高周波数領域とに分けて異なる２つの手法で行われる。低周波数領域（周波数ｆ≒０Ｈｚ）では、まず、電子負荷制御部３０１によって電子負荷装置２００を制御し、所定の直流電流Ｉｄｃに所定の交流成分を重畳した電流を燃料電池１００に対して負荷電流として印加し、電流計測部３０２と電圧計測部３０３で計測された交流波形を低周波数領域インピーダンス算出部３１０に取り込む。次に、低周波数領域インピーダンス算出部３１０は、取り込んだ交流波形から、燃料電池１００から異なる２値の出力電流を出力させた時の各出力電圧値を求め、その求めた出力電流２値の変化量（偏差）Δｉと出力電圧２値の変化量（偏差）Δｖとからその割り算値Δｖ／Δｉを求め、燃料電池１００の周波数０Ｈｚ近傍でのインピーダンス値に決定する。
【００１９】
ここで、上記、低周波数領域（周波数ｆ≒０Ｈｚ）でのインピーダンス計測について図２を参照して詳細に説明する。図２は、燃料電池の電流−電圧特性（破線）を示している。本実施の形態では、低周波数領域（周波数ｆ≒０Ｈｚ）でのインピーダンス測定を、非常に遅い電流変化に対する電圧変化を測定して求めるのではなく、所定の電流Ｉｄｃの近傍でのΔｉおよびΔｖの割り算値Δｖ／Δｉから求めている。実際には、例えば、図７を参照して説明したようなインピーダンス測定における重畳交流振幅Ｉａｃを用いてＩｄｃ−Ｉａｃにおける燃料電池の端子電圧（Ｖ−）とＩｄｃ＋Ｉａｃにおける端子電圧（Ｖ＋）を測定し、下式によってインピーダンスの実数部Ｚ’と虚数部Ｚ”を決定する。
【００２０】
【数１】
2007265895.tif 000003
【００２１】
一方、高周波数領域のインピーダンスの計測は、周波数を変化させながら電子負荷装置２００によって燃料電池１００の出力電流に交流電流を重畳させた負荷を印加し、そのときに電流計測部３０２と電圧計測部３０３で計測された交流波形をインピーダンス計測部３０４に取り込み、その電流の交流分に対する電圧の交流分の応答という形で周波数ｆでのゲインとフェーズを求める、という手法で行う。高周波数領域のインピーダンス計測の際、電子負荷制御部３０１は、高い周波数から低い周波数へとあらかじめ定められた複数の値で周波数を離散的に変化させる指示を出し、インピーダンス計測部３０４は、それらの各周波数で電流及び電圧の計測波形を取り込み、フーリエ演算などを用いて演算によってインピーダンスを求める。この周波数ｆを変えて順次ゲインとフェーズを測定することによって、燃料電池１００の交流インピーダンス特性が得られる。得られたゲイン、フェーズの結果を複素平面上に表したものを、図３に計測結果の一例として示す。図３において、丸印（○）が高周波数領域のインピーダンス計測手法による計測値を示している。また、星印（☆）で示した値が上記低周波数領域のインピーダンス計測手法で算出されたものである。
【００２２】
他方、等価回路フィッティング部３０５は、高周波数領域のインピーダンスの計測において、設定された周波数が所定の周波数ｆｓ以下となった場合に、過去に計測された複数の計測結果に基づいて、あらかじめ定めた所定の等価回路の回路定数を求める演算を行う。等価回路フィッティング部３０５は、その際、周波数０Ｈｚ近傍のインピーダンス値として、インピーダンスの実数部Ｚ’を低周波数領域インピーダンス算出部３１０で決定された上記割り算値の絶対値｜Δｖ／Δｉ｜とし、虚数部Ｚ”を０とする拘束条件を適用する。すなわち、図３に示す例では、高周波数領域での複数の計測値（丸印）と、低周波数領域での算出値（星印）との両方に基づいて、等価回路の各回路定数を求める演算処理が行われることになる。
【００２３】
図４は、本実施の形態で用いる燃料電池の等価回路の一例を示している。抵抗Ｒｍは溶液抵抗に、Ｃａ，Ｃｃは電気二重層容量に、Ｒａ，Ｒｃは反応抵抗に、それぞれ対応している。ただし、図４に示す等価回路は、燃料電池の等価回路の一例であって、等価回路の構成はこれに限らず、より複雑なものや、あるいはより単純なものであってもよい。この場合、等価回路フィッティング部３０５は、インピーダンス計測結果の複数の実測値に基づき、図４に示すような等価回路を表すモデル式から求められるインピーダンスの計算結果と、インピーダンスの実測値との誤差が最も小さくなるように、等価回路の各回路定数を求める演算処理を行う。例えば最小２乗法のような既知の手法を用いることができる。
【００２４】
判定部３０６は、インピーダンス計測部３０４による複数の異なる計測結果に基づいて、等価回路フィッティング部３０５で求められた複数の回路定数の差分（回路定数の変化量）が所定の設定値以下であるか否かを判定する。例えば、各定数Ｒｍ，Ｃａ，Ｃｃ，Ｒａ，Ｒｃのそれぞれについて、それまでの複数の計測結果に基づいて推定された値と、最新の計測結果までを含めた複数の計測結果に基づいて推定された値の差を求め、それが所定の値（例えば０．０１ｍΩあるいは０．１Ｆ）以下であった否かを判定する。さらに、判定部３０６は、判定の結果が設定値以下であった場合には、インピーダンス計測を終了するための指示を出力するとともに、最新の計測結果（最も低い周波数での計測結果）を用いて求められた回路定数に、計測結果の最終値を決定する処理を行う。そして、その結果は、所定の出力装置（図示省略）に出力される。
【００２５】
図３に示すように、電子負荷制御部３０１は、周波数ｆが無限大（∞）に近い値（実際には例えば１０ｋＨｚ）から、あらかじめ定めたステップで周波数を減少させながら、燃料電池１００の直流出力電流に交流分を重畳するような電流制御を行い、インピーダンス計測部３０４が各周波数でインピーダンスを計測する。そして、交流分の周波数が設定値ｆｓ以下となった場合に、等価回路フィッティング部３０４が、等価回路の回路定数を求める演算を開始する。等価回路フィッティング部３０４が演算を開始した後、判定部３０６は、回路定数の差分の判定を開始し、差分が所定の値以下となったところで、電子負荷制御部３０１等の各部に対して計測を終了する旨の指示を出力する。
【００２６】
次に図５を参照して、図１のインピーダンス計測器３００の処理について詳細に説明する。なお、この例では、周波数ｆ（ｉ）（ｉ＝０，１，２，…）にあらかじめ所定の複数の測定周波数の値が設定されていて、燃料電池１００の等価回路として図５に示す５個の定数Ｒｍ，Ｃａ，Ｃｃ，Ｒａ，Ｒｃからなるモデルを用いるものとする。また、燃料電池１００の直流出力電流の値Ｉｄｃについては、図示してない他のルーティンにおいて設定されているものとする。
【００２７】
操作者によって図示していない所定の操作子が操作されたとすると、インピーダンス計測器３００は、直流電流Ｉｄｃに所定の振幅Ｉａｃおよび周波数の交流分を重畳して燃料電池１００の負荷とするような制御を行い、Ｉｄｃ−Ｉａｃにおける燃料電池の端子電圧（Ｖ−）とＩｄｃ＋Ｉａｃにおける端子電圧（Ｖ＋）を測定し、上記の式によってインピーダンスの実数部Ｚ’と虚数部Ｚ”を決定する（ステップＳ００）。
【００２８】
インピーダンス計測器３００は、次に、周波数設定カウンタｉとフィッティング回数カウンタｊの初期化（ｉ＝０，ｊ＝０）を行い、測定処理を開始する（ステップＳ０００）。インピーダンス計測器３００は、周波数ｆ（ｉ）において、対象燃料電池１００のインピーダンス計測を実施し、インピーダンスの実数部Ｚ’（ｉ）、虚数部Ｚ”（ｉ）を取得する（ステップＳ００１）。
【００２９】
次に、等価回路フィッティング開始判定が行われ（ステップＳ００２）、設定周波数ｆ（ｉ）が所定の周波数ｆｓ以下の場合には、ステップＳ００で算出された値と、ステップＳ００１でこれまでに計測された値とに基づいて、等価回路モデルへのフィッティングが実施され（ステップＳ００３）、そうでない場合にはカウンタｉがインクリメントされて（ｉ＝ｉ＋１）、次の設定周波数でインピーダンス計測が行われる（ステップＳ００２１からステップＳ００１）。
【００３０】
ステップＳ００３では等価回路モデルへのフィッティングの結果が、変数Ｒｅｓ（ｊ）に記憶される（ステップＳ００４）。変数Ｒｅｓ（ｊ）は、複数の計測結果から求められた回路定数の推定値Ｒｍ，Ｒａ，Ｃａ，Ｒｃ，Ｃｃを要素とする多次元の配列変数である。
【００３１】
次に、フィッティング定数差分判定を開始するか否かが判定され（ステップＳ００５）、カウンタｊが１より大きい場合にはフィッティング結果比較のための演算が行われ（ステップＳ００６）、他方、１以下の場合はカウンタｊがインクリメントされて（ステップＳ００５１）、ステップＳ００２１及びステップＳ００１〜Ｓ００５の処理が繰り返して行われる。
【００３２】
ステップＳ００６では、フィッティング結果の比較のため、カウンタｊとｊ−１におけるフィッティング結果Ｒｅｓ（ｊ）とＲｅｓ（ｊ−１）の差が求められる。次に、フィッティング結果Ｒｅｓ（ｊ）の差分判定として、ステップＳ００６で求めた各要素（各回路定数Ｒｍ，Ｒａ，Ｃａ，Ｒｃ，Ｃｃ）の差分（変化量）があらかじめ定めた設定値以下であるか否かが判定される（ステップＳ００７）。各要素の差分が設定値以下であった場合には、フィッティング結果Ｒｅｓ（ｊ）を電流値Ｉｄｃでの等価回路モデル定数に決定する処理が行われ（ステップＳ００８）、他方、いずれかの要素の差分値が設定値を超えていた場合にはステップＳ００５１でカウンタｊがインクリメントされ、ステップＳ００２１でカウンタｉがインクリメントされ、そして他の周波数ｆ（ｉ）においてインピーダンス計測ならびにフィッティングおよびその結果の判定処理が行われる（ステップＳ００１〜Ｓ００７）。
【００３３】
本実施の形態では、既定周波数ｆｓ以下の周波数となったところで、まず既定周波数ｆｓまでの交流インピーダンスデータから等価回路モデルへのフィッティングを行って回路定数を求め、その後は周波数の測定終了毎にそれ以前の交流インピーダンスデータからフィッティングを行い、その変化量が設定値以下となったところで、回路定数を決定し、測定を終了する。これによれば、高周波数領域のインピーダンス計測において、比較的低い周波領域での測定を省略することができるので、測定回数を減少させ、さらに比較的低い周波領域での測定における電気化学的反応による測定ばらつきの影響を低減することが期待できる。さらに、等価回路モデルへのフィッティング結果の差分が所定の設定以下となったところで回路定数を決定するので、測定回数を減じたことによる精度の低下も一定の範囲に抑えることができると期待できる。すなわち、本実施の形態によれば、高精度かつ短時間の測定を行うことができ、さらに生成水による測定バラツキへの影響を排除した、燃料電池の特性測定を行うことが可能となるのである。
【００３４】
以上のように、本実施の形態によれば、低周波数領域（ｆ≒０Ｈｚ付近）において、交流信号に対する応答性を測定する必要がなく、直流成分のみで、インピーダンス特性を測定できるため、測定の高速化が可能となる。また、燃料電池の発電で発生する生成水による一時的な流路閉塞、あるいは圧力損失の増加が起因となる燃料ガス不足、あるいは酸化ガス不足による燃料電池端子電圧の不安定性（インピーダンス計測でのバラツキ、再現性低下の要因）を除去した測定が可能である。
【００３５】
上述したように、本実施の形態では、燃料電池のインピーダンス特性の測定のはじめに、低周波数（ｆ≒０Ｈｚ）でのインピーダンス計測を行い、その後、周波数を変化させながらでのインピーダンス計測（高周波数から低周波数へ）を行う。また、周波数ｆｓ以下の比較的低い周波数領域（測定前に指定）に入ったら、インピーダンス計測結果が出る度に、等価回路モデルへのフィッティングを行う。その際に、低周波数（ｆ≒０Ｈｚ）でのインピーダンスの実数部Ｚ’として最初に測定した値で拘束を与えて、等価回路モデルの定数の変化が設定値以下となったら、定数の収束と判断し、インピーダンス計測を終える。これによれば、低周波数の最後までインピーダンス計測を実施し、等価回路へのフィッティングを行わなくてもよい。生成水による低周波数領域での測定バラツキ、再現性の低下を除去することが可能となる。
【００３６】
なお、上記構成において、図１の判定部３０６で、あらかじめ設定した所定の周波数以下になっても等価回路定数の変化量があらかじめ定めた設定値以下とならないことや、等価回路定数の変化量が増加傾向にあることを検知した場合には、その旨を出力するような機能を追加することも可能である。すなわち、低周波数（ｆ≒０）の拘束条件を加えたフィッティングにおいて、その回路定数が収束せずしないもしくは発散していく場合には、生成水によるフラッディング（flooding）などの不具合が発生していると判断し、その旨を表示したり、自動的に測定を停止するような処理を行うことができる。このようにすれば、フラッディングなどが発生するような環境条件でのインピーダンス測定を行うことがなくなる。
【００３７】
次に、図６を参照して、本発明の他の実施の形態について説明する。本実施の形態は、図１と同様の構成を用いて、上記の実施の形態とは異なる手法で低周波数領域（ｆ≒０）でのインピーダンス値を算出するものである。
【００３８】
本実施の形態では、まず、交流インピーダンスの測定前にその試験条件（温度、圧力、湿度などのインピーダンス測定を行いたい環境条件）における開回路電圧Vｏｃ（破線）を測定、もしくは、演算するとともに、図６に示すような電流−電圧特性（実線）を測定する。
【００３９】
次に、低周波数領域（ｆ≒０）での電流依存性を考慮した燃料電池１００の抵抗値Ｒ（ｉ）と電圧値Ｖｘとの関係を表す式（式（１））に基づいて、開回路電圧Vｏｃから、電流−電圧特性の測定で得られた各電流での電圧Ｖｘを引算したものを、電流値Ｉｘに対して微分し（式（２））、その結果Ｒ（ｉ）を低周波数（ｆ≒０）でのインピーダンスの実数部Ｚ’とする。
【００４０】
【数２】
2007265895.tif 000004
【００４１】
そして、本実施の形態においても、求めたインピーダンス値を上記の実施の形態と同様にフィッティング時の拘束条件として用いる。
【００４２】
すなわち、本実施の形態では、交流インピーダンスの測定前にその試験条件（温度、圧力、湿度など）における開回路電圧を測定もしくは演算し、また電流−電圧特性を測定する。そして、開回路電圧から電流−電圧特性を引算したものを電流変化に対して微分する。その結果を、インピーダンス計測を行う電流値Ｉｘでの当該燃料電池の０Ｈｚ近傍でのインピーダンス値として決定する。このインピーダンス値は、インピーダンス計測を行う電流値での低周波数側の拘束条件として用いることができる。
【００４３】
この実施の形態によれば、上記の実施の形態と同様に、低周波数領域（ｆ≒０）での測定を省略することができるので、計測の高速化および安定化の効果が期待できる。また、電流−電圧特性の測定において、電流ピッチを小さくしておけば、上記の実施の形態で直流電流Ｉｄｃの近傍電流での電圧変化Δｖを線形と置いた場合よりも精度が向上すると考えられる。
【００４４】
なお、本発明の実施の形態は、上記のものに限定されず、さらに他の変更を行うことが可能である。例えば、図４のステップＳ００７の判定処理においてすべての回路定数に対して比較を行うのではなく、特定の回路定数に限定することが可能である。特定の回路定数（例えばＲａのみとか、Ｒｃのみとか）の差分（変化量）が所定の範囲内にあるときに、そのときの値ですべての回路定数を決定するようにする。
【００４５】
また、本発明の構成は、図１に示す各構成を含むもの、燃料電池１００以外の構成からなるもの、インピーダンス計測器３００のみからなるもの等として考えることができる。さらに、インピーダンス計測器３００において、例えば電流計測部３０２および電圧計測部３０３を汎用の計測器から構成し、電子負荷制御部３０１、インピーダンス計測部３０４、等価回路フィッティング部３０５、判定部３０６、および低周波数領域インピーダンス算出部３１０を汎用のコンピュータを用いて構成することも可能であり、その場合には、本発明の構成を電子負荷制御部３０１、インピーダンス計測部３０４、等価回路フィッティング部３０５、判定部３０６および低周波数領域インピーダンス算出部３１０からなるものととらえることが可能であり、また、電子負荷制御部３０１、インピーダンス計測部３０４、等価回路フィッティング部３０５、判定部３０６および低周波数領域インピーダンス算出部３１０の各機能をプログラムによって実現し、それを実現するソフトウェアプログラムを製品として販売することも可能である。その際、プログラムをコンピュータ読み取り可能な記録媒体に記録して、あるいは通信回線を介して配布することが可能である。
【図面の簡単な説明】
【００４６】
【図１】本発明の燃料電池の特性計測装置の実施の形態のブロック図。
【図２】図１の構成によるインピーダンス計測手法の一例を説明するために用いる電流−電圧特性図。
【図３】図１の構成によるインピーダンス計測手法の一例を説明するための複素平面図。
【図４】燃料電池の等価回路モデルの一例を示す図。
【図５】図１の機能を説明するためのフローチャート。
【図６】本発明の燃料電池の特性計測装置の他の実施の形態を説明するために用いる電流−電圧特性図。
【図７】燃料電池の交流インピーダンス計測手法を説明するための図。
【符号の説明】
【００４７】
１００…燃料電池、２００…電子負荷装置、３００…インピーダンス計測器、３０１…電子負荷制御部（負荷制御手段）、３０２…電流計測部（計測手段）、３０３…電圧計測部（計測手段）、３０４…インピーダンス計測部（インピーダンス計測手段）、３０５…等価回路フィッティング部（演算手段）、３０６…判定部（判定手段）、３１０…低周波数領域インピーダンス算出部（計測手段，決定手段）
【図１】
2007265895.tif 000005
【図２】
2007265895.tif 000006
【図３】
2007265895.tif 000007
【図４】
2007265895.tif 000008
【図５】
2007265895.tif 000009
【図６】
2007265895.tif 000010
【図７】
2007265895.tif 000011