JPWO2016103462A1 - 燃料電池の状態検出装置及び方法 - Google Patents

Abstract

アノードガス及びカソードガスの供給を受けて発電する燃料電池の状態検出装置であって、少なくともアノード極の状態量に対して応答性を示す周波数帯を含む高周波数帯から選択された周波数に基づく高周波数インピーダンス、及び少なくともカソード極の状態量に対して応答性を示す周波数帯を含む低周波数帯から選択された周波数に基づく低周波数インピーダンスを取得するインピーダンス取得手段と、取得された前記高周波数インピーダンス及び前記低周波数インピーダンスを組み合わせて、前記燃料電池の内部状態としての前記アノード極の状態量と前記カソード極の状態量をそれぞれ推定する内部状態量推定手段と、を備えた燃料電池の状態検出装置。

Description

この発明は、燃料電池の状態検出装置及び方法に関する。

燃料電池の電圧値やインピーダンス値を測定し、これらの値に基づいて燃料電池の内部状態を検出する燃料電池の状態検出装置が知られている。

例えば、特許第４６４０６６１号には、電解質膜抵抗に対応する第１の周波数領域における第１のインピーダンス、及び電解質膜抵抗と触媒層抵抗との合算値に対応する第１の周波数領域よりも低い第２の周波数領域における第２のインピーダンスをそれぞれ算出し、第２のインピーダンスと第１のインピーダンスとの差である差分インピーダンスに基づいて触媒層の含水量を算出することが提案されている。

また、特開２００５−２８５６１４号公報には、燃料電池の複素インピーダンス曲線（コール・コールプロット）の実軸との交点における周波数Ｆ₁、酸素が反応する際の反応抵抗（カソード極の反応抵抗）を表す第１の領域内の周波数Ｆ₂、及び酸素の拡散に関する抵抗を表す第２の領域内の周波数Ｆ₃に対応する複素インピーダンスを取得し、取得した複素インピーダンスから内部抵抗値を求めることが記載されている。

しかしながら、特許第４６４０６６１号では、アノード極の状態量とカソード極の状態量をそれぞれ把握することはできない。また、特開２００５−２８５６１４号公報においても、インピーダンス曲線にはアノード極の状態及びカソード極の状態が混在しており、アノード極の状態量とカソード極の状態量を個別に把握することは難しい。

本発明は、このような問題点に着目してなされたものであり、燃料電池におけるアノード極の状態量やカソード極の状態量等の内部状態量を個別に検出し得る燃料電池の状態検出装置及び方法を提供することを目的とする。

本発明のある態様によれば、アノードガス及びカソードガスの供給を受けて発電する燃料電池の状態検出装置が提供される。より詳細には、この状態検出装置は、少なくともアノード極の状態量に対して応答性を示す周波数帯を含む高周波数帯から選択された周波数に基づく高周波数インピーダンス、及び少なくともカソード極の状態量に対して応答性を示す周波数帯を含む低周波数帯から選択された周波数に基づく低周波数インピーダンスを取得するインピーダンス取得手段を備える。また、状態検出装置は、取得された高周波数インピーダンス及び低周波数インピーダンスを組み合わせて、燃料電池の内部状態としてのアノード極の状態量とカソード極の状態量をそれぞれ推定する内部状態量推定手段を備える。

図１は、本発明の実施形態による燃料電池セルの斜視図である。 図２は、図１の燃料電池のII−II断面図である。 図３は、本発明の実施形態による燃料電池システムの概略構成図である。 図４Ａは、低周波数帯の交流電圧を入力した場合において燃料電池の簡易等価回路モデルに流れる電流の経路を示す図である。 図４Ｂは、図４Ａの場合よりも高い周波数帯の交流電圧を入力した場合において燃料電池の簡易等価回路モデルに流れる電流の経路を示す図である。 図４Ｃは、図４Ｂの場合よりも高い周波数帯の交流電圧を入力した場合において燃料電池簡易等価回路モデルに流れる電流の経路を示す図である。 図４Ｄは、高周波数帯の交流電圧を入力した場合において燃料電池の簡易等価回路モデルに流れる電流の経路を示す図である。 図５は、一実施形態に係る状態量の推定の流れを示すフローチャートである。 図６は、一実施形態に係る状態量の推定の流れを示すフローチャートである。 図７は、定常時及び非定常時における燃料電池のＩ−Ｖ特性線をそれぞれ示した図である。 図８は、一実施形態に係る状態量の推定の流れを示すフローチャートである。 図９は、カソード極の電気二重層容量の候補の周波数応答を示している。 図１０Ａは、アノード極の電気二重層容量の候補の周波数応答を示している。 図１０Ｂは、アノード極１１２の反応抵抗値の候補の周波数応答を示している。 図１１は、一実施形態に係る状態量の推定の流れを示すフローチャートである。 図１２は、定常時における燃料電池１のＩ−Ｖ特性線を示している。 図１３は、Ｉ−Ｖ特性線における傾きΔＶ／ΔＩの算出を行うための１組の電流及び電圧の設定方法の一例を説明する図である。 図１４は、一実施形態に係る燃料電池システムにおいて、インピーダンス計測に係る要部を概略的に示したブロック図である。

以下、図面等を参照して本発明の実施形態について説明する。

燃料電池セルは、燃料極としてのアノード極と酸化剤極としてのカソード極とによって電解質膜を挟んで構成されている。燃料電池のセルでは、水素を含有するアノードガスがアノード極に供給される一方で、酸素を含有するカソードガスがカソード極に供給されて、これらガスを用いることで発電が行われる。アノード極及びカソード極の両電極において進行する電極反応は、以下の通りである。

アノード極：２Ｈ₂ → ４Ｈ⁺＋４ｅ^-
カソード極：４Ｈ⁺＋４ｅ^-＋Ｏ₂ → ２Ｈ₂Ｏ
図１及び図２は、本発明の一実施形態による燃料電池セル１０の構成を説明するための図である。図１は燃料電池セル１０の斜視図であり、図２は図１の燃料電池セル１０のII−II断面図である。

図１及び図２に示すように、燃料電池セル１０は、膜電極接合体（ＭＥＡ）１１と、ＭＥＡ１１を挟むように配置されるアノードセパレータ１２及びカソードセパレータ１３と、を備える。

ＭＥＡ１１は、電解質膜１１１と、アノード極１１２と、カソード極１１３とから構成されている。ＭＥＡ１１は、電解質膜１１１の一方の面側にアノード極１１２を有しており、他方の面側にカソード極１１３を有している。

電解質膜１１１は、フッ素系樹脂により形成されたプロトン伝導性のイオン交換膜である。電解質膜１１１は、湿潤状態で良好な電気伝導性を示す。なお、電解質膜１１１としては、想定される燃料電池の対応に応じて、例えばリン酸（Ｈ₃ＰＯ₄）を所定のマトリックスに含浸させたものなどの他の材料を用いるようにしても良い。

アノード極１１２は、触媒層１１２Ａとガス拡散層１１２Ｂとを備える。触媒層１１２Ａは、白金又は白金等が担持されたカーボンブラック粒子により形成された部材であって、電解質膜１１１と接するように設けられる。ガス拡散層１１２Ｂは、触媒層１１２Ａの外側に配置される。ガス拡散層１１２Ｂは、ガス拡散性及び導電性を有するカーボンクロスで形成された部材であって、触媒層１１２Ａ及びアノードセパレータ１２と接するように設けられる。

アノード極１１２と同様に、カソード極１１３も触媒層１１３Ａとガス拡散層１１３Ｂとを備える。触媒層１１３Ａは電解質膜１１１とガス拡散層１１３Ｂとの間に配置され、ガス拡散層１１３Ｂは触媒層１１３Ａとカソードセパレータ１３との間に配置される。

アノードセパレータ１２は、ガス拡散層１１２Ｂの外側に配置される。アノードセパレータ１２は、アノード極１１２にアノードガス（水素ガス）を供給するための複数のアノードガス流路１２１を備えている。アノードガス流路１２１は、溝状通路として形成されている。

カソードセパレータ１３は、ガス拡散層１１３Ｂの外側に配置される。カソードセパレータ１３は、カソード極１１３にカソードガス（空気）を供給するための複数のカソードガス流路１３１を備えている。カソードガス流路１３１は、溝状通路として形成されている。

アノードセパレータ１２及びカソードセパレータ１３は、アノードガス流路１２１を流れるアノードガスの流れ方向とカソードガス流路１３１を流れるカソードガスの流れ方向とが互いに逆向きとなるように構成されている。なお、アノードセパレータ１２及びカソードセパレータ１３は、これらガスの流れ方向が同じ向きに流れるように構成されてもよい。

このような燃料電池セル１０を自動車用電源として使用する場合には、要求される電力が大きいため、数百枚の燃料電池セル１０を積層した燃料電池スタックとして使用する。そして、燃料電池スタックにアノードガス及びカソードガスを供給する燃料電池システムを構成して、車両を駆動させるための電力を取り出す。なお、本実施形態では、燃料電池セル１０が積層された燃料電池スタックの単位で後述するインピーダンス測定を行うが、燃料電池セル１０一枚の単位や燃料電池スタックの一部（例えば数十枚のセル）の単位でインピーダンス測定を行うようにしても良い。

また、燃料電池スタックにおいては、一枚の燃料電池セル１０におけるアノード極１１２、カソード極１１３、及び電解質膜１１１が、複数枚直列に配置されることによって総和としてのアノード極、カソード極、及び電解質膜が構成されることとなる。しかしながら、以下では説明の便宜上、この総和としてのアノード極、カソード極、及び電解質膜についても、セル単体のアノード極１１２、カソード極１１３、及び電解質膜１１１と同一の符号を付す。

図３は、本発明の一実施形態による燃料電池システム１００の概略図である。

燃料電池システム１００は、燃料電池１と、カソードガス給排装置２と、アノードガス給排装置３と、電力システム５と、コントローラ６と、を備える。

燃料電池１は、上述のように複数枚の燃料電池セル１０（単位セル）を積層した積層電池である。燃料電池１は、アノードガス及びカソードガスの供給を受けて、車両の走行に必要な電力を発電する。燃料電池１は、電力を取り出す出力端子として、アノード極側端子１Ａと、カソード極側端子１Ｂと、を有している。

カソードガス給排装置２は、燃料電池１にカソードガスを供給するとともに、燃料電池１から排出されるカソードオフガスを外部に排出する。カソードガス給排装置２は、カソードガス供給通路２１と、カソードガス排出通路２２と、フィルタ２３と、エアフローセンサ２４と、カソードコンプレッサ２５と、カソード圧力センサ２６と、水分回収装置（ＷＲＤ；Water Recovery Device）２７と、カソード調圧弁２８と、を備える。

カソードガス供給通路２１は、燃料電池１に供給されるカソードガスが流れる通路である。カソードガス供給通路２１の一端はフィルタ２３に接続され、他端は燃料電池１のカソードガス入口部に接続される。

カソードガス排出通路２２は、燃料電池１から排出されるカソードオフガスが流れる通路である。カソードガス排出通路２２の一端は燃料電池１のカソードガス出口部に接続され、他端は開口端として形成される。カソードオフガスは、カソードガスや電極反応によって生じた水蒸気等を含む混合ガスである。

フィルタ２３は、カソードガス供給通路２１に取り込まれるカソードガスに含まれる塵や埃等を除去する部材である。

カソードコンプレッサ２５は、フィルタ２３よりも下流側のカソードガス供給通路２１に設けられる。カソードコンプレッサ２５は、カソードガス供給通路２１内のカソードガスを圧送して燃料電池１に供給する。

エアフローセンサ２４は、フィルタ２３とカソードコンプレッサ２５との間のカソードガス供給通路２１に設けられる。エアフローセンサ２４は、燃料電池１に供給されるカソードガスの流量を検出する。

カソード圧力センサ２６は、カソードコンプレッサ２５とＷＲＤ２７との間のカソードガス供給通路２１に設けられる。カソード圧力センサ２６は、燃料電池１に供給されるカソードガスの圧力を検出する。カソード圧力センサ２６で検出されたカソードガス圧力は、燃料電池１のカソードガス流路等を含むカソード系全体の圧力を代表する。

ＷＲＤ２７は、カソードガス供給通路２１とカソードガス排出通路２２とに跨って接続される。ＷＲＤ２７は、カソードガス排出通路２２を流れるカソードオフガス中の水分を回収し、その回収した水分を用いてカソードガス供給通路２１を流れるカソードガスを加湿する装置である。

カソード調圧弁２８は、ＷＲＤ２７よりも下流のカソードガス排出通路２２に設けられる。カソード調圧弁２８は、コントローラ６によって開閉制御され、燃料電池１に供給されるカソードガスの圧力を調整する。

次に、アノードガス給排装置３について説明する。

アノードガス給排装置３は、燃料電池１にアノードガスを供給するとともに、燃料電池１から排出されるアノードオフガスをカソードガス排出通路２２に排出する。アノードガス給排装置３は、高圧タンク３１と、アノードガス供給通路３２と、アノード調圧弁３３と、アノード圧力センサ３４と、アノードガス排出通路３５と、バッファタンク３６と、パージ通路３７と、パージ弁３８と、を備える。

高圧タンク３１は、燃料電池１に供給するアノードガスを高圧状態に保って貯蔵する容器である。

アノードガス供給通路３２は、高圧タンク３１から排出されるアノードガスを燃料電池１に供給する通路である。アノードガス供給通路３２の一端は高圧タンク３１に接続され、他端は燃料電池１のアノードガス入口部に接続される。

アノード調圧弁３３は、高圧タンク３１よりも下流のアノードガス供給通路３２に設けられる。アノード調圧弁３３は、コントローラ６によって開閉制御され、燃料電池１に供給されるアノードガスの圧力を調整する。

アノード圧力センサ３４は、アノード調圧弁３３よりも下流のアノードガス供給通路３２に設けられる。アノード圧力センサ３４は、燃料電池１に供給されるアノードガスの圧力を検出する。アノード圧力センサ３４で検出されたアノードガス圧力は、バッファタンク３６や燃料電池１のアノードガス流路等を含むアノード系全体の圧力を代表する。

アノードガス排出通路３５は、燃料電池１から排出されたアノードオフガスを流す通路である。アノードガス排出通路３５の一端は燃料電池１のアノードガス出口部に接続され、他端はバッファタンク３６に接続される。アノードオフガスには、電極反応で使用されなかったアノードガスや、カソードガス流路１３１からアノードガス流路１２１へとリークしてきた窒素等の不純物ガスや水分等が含まれる。

バッファタンク３６は、アノードガス排出通路３５を流れてきたアノードオフガスを一時的に蓄える容器である。バッファタンク３６に溜められたアノードオフガスは、パージ弁３８が開かれる時に、パージ通路３７を通ってカソードガス排出通路２２に排出される。

パージ通路３７は、アノードオフガスを排出するための通路である。パージ通路３７の一端はアノードガス排出通路３５に接続され、他端はカソード調圧弁２８よりも下流のカソードガス排出通路２２に接続される。

パージ弁３８は、パージ通路３７に設けられる。パージ弁３８は、コントローラ６によって開閉制御され、アノードガス排出通路３５からカソードガス排出通路２２に排出するアノードオフガスのパージ流量を制御する。

パージ弁３８が開弁状態となるパージ制御が実行されると、アノードオフガスは、パージ通路３７及びカソードガス排出通路２２を通じて外部に排出される。この時、アノードオフガスは、カソードガス排出通路２２内でカソードオフガスと混合される。このようにアノードオフガスとカソードオフガスとを混合させて外部に排出することで、混合ガス中のアノードガス濃度（水素濃度）が排出許容濃度以下の値に設定される。

電力システム５は、電流センサ５１と、電圧センサ５２と、走行モータ５３と、インバータ５４と、バッテリ５５と、ＤＣ／ＤＣコンバータ５６と、を備える。

電流センサ５１は、燃料電池１から取り出される出力電流を検出する。電圧センサ５２は、燃料電池１の出力電圧、つまりアノード極側端子１Ａとカソード極側端子１Ｂの間の端子間電圧を検出する。電圧センサ５２は、燃料電池セル１０の１枚ごとの電圧を検出するように構成されてもよいし、燃料電池セル１０の複数枚ごとの電圧を検出するように構成されてもよい。

走行モータ５３は、三相交流同期モータであって、車輪を駆動するため駆動源である。走行モータ５３は、燃料電池１及びバッテリ５５から電力の供給を受けて回転駆動する電動機としての機能と、外力によって回転駆動されることで発電する発電機としての機能と、を有する。

インバータ５４は、ＩＧＢＴ等の複数の半導体スイッチから構成される。インバータ５４の半導体スイッチは、コントローラ６によってスイッチング制御され、これにより直流電力が交流電力に、又は交流電力が直流電力に変換される。走行モータ５３を電動機として機能させる場合、インバータ５４は、燃料電池１の出力電力とバッテリ５５の出力電力との合成直流電力を三相交流電力に変換し、走行モータ５３に供給する。これに対して、走行モータ５３を発電機として機能させる場合、インバータ５４は、走行モータ５３の回生電力（三相交流電力）を直流電力に変換し、バッテリ５５に供給する。

バッテリ５５は、燃料電池１の出力電力の余剰分及び走行モータ５３の回生電力が充電されるように構成されている。バッテリ５５に充電された電力は、必要に応じてカソードコンプレッサ２５等の補機類や走行モータ５３に供給される。

ＤＣ／ＤＣコンバータ５６は、燃料電池１の出力電圧を昇降圧させる双方向性の電圧変換機である。ＤＣ／ＤＣコンバータ５６によって燃料電池１の出力電圧を制御することで、燃料電池１の出力電流等が調整される。

コントローラ６は、中央演算装置（ＣＰＵ）、読み出し専用メモリ（ＲＯＭ）、ランダムアクセスメモリ（ＲＡＭ）及び入出力インタフェース（Ｉ／Ｏインタフェース）を備えたマイクロコンピュータで構成される。コントローラ６には、電流センサ５１や電圧センサ５２等の各種センサからの信号の他、アクセルペダルの踏み込み量を検出するアクセルストロークセンサ（図示せず）等のセンサからの信号が入力される。

コントローラ６は、燃料電池システム１００の運転状態に応じて、アノード調圧弁３３やカソード調圧弁２８、カソードコンプレッサ２５等を制御し、燃料電池１に供給されるアノードガスやカソードガスの圧力や流量を調整する。

また、コントローラ６は、走行モータ５３の要求電力やカソードコンプレッサ２５等の補機類の要求電力、バッテリ５５の充放電要求等に基づいて、目標出力電力を算出する。コントローラ６は、目標出力電力に基づいて、予め定められた燃料電池１のＩＶ特性（電流電圧特性）を参照して燃料電池１の目標出力電流を算出する。そして、コントローラ６は、燃料電池１の出力電流が目標出力電流となるように、ＤＣ／ＤＣコンバータ５６によって燃料電池１の出力電圧を制御し、走行モータ５３や補機類に必要な電流を供給する制御を行う。

さらに、コントローラ６は、燃料電池１の各電解質膜１１１の湿潤度（含水量）が発電に適した状態となるように、カソードコンプレッサ２５等を制御する。

また、コントローラ６は、後述する第１〜第６実施形態において、燃料電池１の出力電圧に所定周波数の交流信号が重畳された電圧値の振幅値を、同じく交流信号が重畳された電流値の振幅値を除して、所定周波数における燃料電池１のインピーダンスＺを算出する。

上述のように説明した燃料電池システム１００において、コントローラ６、電流センサ５１、電圧センサ５２、及びＤＣ／ＤＣコンバータ５６により燃料電池１の状態検出装置が構成される。

本実施形態では、燃料電池１においてアノード極１１２の状態量である反応抵抗値Ｒ_a及び電気二重層容量Ｃ_a、カソード極１１３の状態量である反応抵抗値Ｒ_c及び電気二重層容量Ｃ_c、並びに電解質膜１１１の状態量である電解質膜抵抗値Ｒ_mを考慮した簡易的な等価回路モデルを設定し、該簡易等価回路モデルに基づいて燃料電池１の状態推定を行う。

なお、電解質膜抵抗値Ｒ_mは、電解質膜１１１の湿潤度に応じてその値が定まる状態量である。通常、電解質膜１１１が乾燥するにつれて電解質膜抵抗値Ｒ_mが高くなる傾向にある。

また、アノード極１１２の反応抵抗値Ｒ_aは、アノード極１１２におけるアノードガスの反応に応じて増減し、例えばアノードガスが不足している等の当該反応の進行が円滑に行われない要因があるとこれに応じて反応抵抗値Ｒ_aは上昇する。

さらに、アノード極１１２の電気二重層容量Ｃ_aは、燃料電池１においてアノード極１１２が有する電気容量を表すようにモデル化したものである。従って、電気二重層容量Ｃ_aはアノード極１１２を構成する材料や大きさ等の種々の要素に基づいて決定されることとなる。

また、カソード極１１３の反応抵抗値Ｒ_cは、カソード極１１３におけるカソードガスの反応に応じて増減し、例えばカソードガスが不足している等の当該反応の進行が円滑に行われない要因があるとこれに応じて反応抵抗値Ｒ_cは上昇する。

さらに、カソード極１１３の電気二重層容量Ｃ_cは、カソード極１１３が有する電気容量を表すようにモデル化したものである。従って、電気二重層容量成分の値Ｃ_cはカソード極１１３を構成する材料や大きさ等の種々の要素に基づいて決定されることとなる。

ここで、本発明者らは、燃料電池１の簡易等価回路モデルにおいて、燃料電池１の出力電流に重畳された交流信号（交流電流）の燃料電池内部において流れる経路には、周波数依存特性が存在することを見出している。以下では交流電流が燃料電池内部において流れる経路における周波数依存特性について説明する。

図４Ａ〜図４Ｄは、本実施形態に係る燃料電池１の簡易等価回路モデルにおいて、燃料電池１の出力電流に重畳された交流電流が燃料電池内部において流れる経路を、交流電流の周波数帯ごとに模式的に示した図である。

図４Ａにおいては、例えば０Ｈｚ近傍の低周波数帯（以下では、第１周波数帯とも記載する）に属する周波数の交流電流の経路を示している。また、図４Ｂにおいては、例えば、数Ｈｚ程度の第１周波数帯と比較して若干高い周波数帯（以下では、第２周波数帯とも記載する）に属する周波数の交流電流の経路を示している。さらに、図４Ｃにおいては、例えば、数十Ｈｚ〜数ＫＨｚの第２周波数帯よりも若干高い周波数帯（以下、第３周波数帯とも記載する）に属する周波数の交流電流の経路を示している。また、図４Ｄでは、例えば、数十ＫＨｚ以上の最も高い周波数帯（以下、第４周波数帯とも記載する）に属する周波数の交流電流の経路を示している。なお、図４Ａ〜図４Ｄにおいては交流電流の経路は太字で示している。

先ず、図４Ａに示した第１周波数帯に属する周波数の交流電流では、周波数が低いために値の変動が緩やかであり、電流の値が一定値である直流に近い性質を有することとなる。従って、このように直流に近い性質の交流電流は、アノード極１１２の電気二重層容量及びカソード極１１３の電気二重層容量側部分には流れ込まないか、又は流れ込んだとしてもその大きさが無視できる程度に小さくなる。すなわち、図で示すように、交流電流は、実質的に、アノード極１１２の反応抵抗、電解質膜抵抗、及びカソード極１１３の反応抵抗の部分にのみ流れることとなる。

次に、図４Ｂに示す第２周波数帯に属する周波数の交流電流では、上記第１周波数帯に属する周波数の交流電流と比較すると値の変動が大きくなり、より交流としての性質が強くなる。したがって、図に示すように、カソード極１１３の電気二重層容量側にも交流電流が流れはじめると考えられる。

一方で、アノード極１１２の反応抵抗値Ｒａは、カソード極１１３の反応抵抗値Ｒ_cと比較してかなり小さい値をとることが知られているので、アノード極１１２の反応抵抗側には比較的電流が流れやすい。したがって、第２周波数帯における周波数程度の交流電流では、依然として、アノード極１１２の電気二重層容量側部分には流れ込まないか、又は流れ込んだとしてもその大きさが無視できる程度に小さくなる。と考えられる。

さらに、図４Ｃに示す第３周波数帯に属する周波数の交流電流では、上記第２周波数帯に属する周波数の交流電流と比較して値の変動がより大きくなるので、さらに交流としての性質が強くなる。したがって、アノード極１１２の電気二重層容量の影響も無視できなくなり、アノード極１１２の電気二重層容量部分にも電流が流れると考えられる。

一方で、この第３周波数帯では、上記交流電流の値の変動速度に対して、カソード極１１３における酸化還元反応が追随できなくなり、当該酸化還元反応が見かけ上生じなくなるという状態が生じる。

したがって、実質的にカソード極１１３におけるカソードガスの反応が生じないこととなるので、上記酸化還元反応に起因するカソード極１１３の反応抵抗の影響を無視することができる。

すなわち、第３周波数帯においては、図４Ｃに示すように、交流電流がカソード極１１３の反応抵抗に流れ込まないか、又は流れ込んだとしてもその大きさが無視できる程度に小さくなるので、実質的に電気二重層容量成分にのみ流れると考えられる。

なお、アノード極１１２においては交流電流の値の変動に対する酸化還元反応の追従性能が比較的高く、当該酸化還元反応は、第３周波数帯ではまだ交流電流の値の変動に追従することができる。したがって、図に示すように、第３周波数帯に属する周波数の交流電流は、依然としてアノード極１１２の反応抵抗を流れると考えられる。

そして、図４Ｄに示す第４周波数帯に属する周波数の交流電流では、上記第３周波数帯に属する周波数の交流電流と比較して値の変動がさらに大きくなるので、カソード極１１３だけではなくアノード極１１２における酸化還元反応も、当該交流電流の値の変動に追随できなくなる。

従って、カソード極１１３に加えてアノード極１１２における反応も実質的に生じないこととなり、カソード極１１３の反応抵抗及びアノード極１１２の反応抵抗の双方の影響を無視することができる。

すなわち、第４周波数帯においては、交流電流がカソード極１１３及びアノード極１１２の双方の反応抵抗に流れ込まないか、又は流れ込んだとしてもその大きさが無視できる程度に小さくなる。したがって、図に示すように、第４周波数帯に属する周波数の交流電流は、カソード極１１３及びアノード極１１２のそれぞれの電気二重層容量側にのみ流れることとなる。

以上の説明から理解されるように、上述の第１周波数帯から選択される周波数の交流電流、第２周波数帯から選択される周波数の交流電流、第３周波数帯から選択される周波数の交流電流、及び第４周波数帯から選択される周波数の交流電流に対して、燃料電池の簡易等価回路における各要素に流れる経路が異なる。

したがって、本発明者らは、このように周波数に応じた交流電流の経路の違いを利用して、簡易等価回路に基づいて得られるインピーダンスの式

（ただし、ｊは虚数単位を意味する。）
を参照し、各周波数帯に属する周波数に基づくインピーダンスから種々の状態量を個々に推定することに想到した。

例えば、上記第４周波数帯（以下、「電解質膜応答周波数帯」とも記載する）から選択される周波数の交流電流は、電解質膜抵抗、アノード極１１２の電気二重層容量、及びカソード極１１３の電気二重層容量の部分に流れるので、当該電解質膜応答周波数帯から選択される周波数に基づくインピーダンス（以下、「電解質膜応答インピーダンス」とも記載する）には、電解質膜抵抗値Ｒ_mの情報が含まれている。

なお、この電解質膜応答周波数帯はいわゆるＨＦＲ（High Frequency Resistance）に用いられる周波数帯である。したがって、インピーダンスの式（１）においてω→∞とすると、インピーダンスＺは、電解質膜抵抗値Ｒ_mにほぼ一致するとみなすことができる。

また、第３周波数帯（以下では「アノード極応答周波数帯」とも記載する）から選択される周波数の交流電流は、電解質膜抵抗、アノード極１１２の反応抵抗、アノード極１１２の電気二重層容量、及びカソード極１１３の電気二重層容量の部分に流れるので、当該アノード極応答周波数帯から選択される周波数に基づくインピーダンス（以下、「アノード極応答インピーダンス」と記載する）には、少なくともアノード極１１２の反応抵抗値Ｒ_a及びアノード極１１２の電気二重層容量値Ｃ_aの情報が含まれている。

特にこの場合、図４Ｃで示した等価回路においてカソード極１１３の反応抵抗を無視することができるので、インピーダンスの式は、

で与えられる。

さらに、第２周波数帯から選択される周波数の交流電流は、電解質膜抵抗、アノード極１１２の反応抵抗、カソード極１１３の反応抵抗、及びカソード極１１３の電気二重層容量の部分に流れるので、当該第２周波数帯から選択される周波数に基づくインピーダンスには、状態量として電解質膜抵抗値、アノード極１１２の反応抵抗値、カソード極１１３の反応抵抗値Ｒ_c、及びカソード極１１３の電気二重層容量値Ｃ_cの情報が含まれている。

また、最も低い周波数帯である第１周波数帯（以下では「低周波数帯」とも記載する）から選択される周波数の交流電流は、電解質膜抵抗、アノード極１１２の反応抵抗、及びカソード極１１３の反応抵抗の部分に流れるので、当該低周波数帯から選択される周波数に基づくインピーダンス（以下、低周波数インピーダンスと記載する）には、少なくともカソード極１１３の反応抵抗値Ｒ_cの情報が含まれている。

以下では、各実施形態において、上記電解質膜応答周波数帯、アノード極応答周波数帯、及び低周波数帯燃料のうちの少なくとも２つを用いた各状態量の推定の詳細について説明する。

なお、一般に「周波数ｆ」と「角周波数ω」との間にはω＝２πｆの関係があることは知られており、これらの間には無次元の定数２πを乗じた差異しかないため、各実施形態においては説明の簡略化のため、「周波数」と「角周波数」を同一視し、いずれを表す場合にも「ω」の記号を用いる。

（第１の実施形態）
以下、第１の実施形態について説明する。

図５は、本実施形態に係る状態量の推定の流れを示すフローチャートである。

図示のように、先ず、ステップＳ１０１において、電解質膜応答周波数帯における１点の周波数ω_Hを選択し、周波数ω_Hに基づくインピーダンスＺ（ω_H）を求める。

具体的には、コントローラ６が、インピーダンス計測タイミングにおいて、燃料電池１から出力される出力電圧及び出力電流に、電解質膜応答周波数帯の周波数ω_Hの交流信号が重畳されるようにＤＣ／ＤＣコンバータ５６を制御する。

さらに、コントローラ６は、電圧センサ５２で測定された出力電圧の値Ｖにフーリエ変換を施し電圧振幅値Ｖ（ω_H）を得て、電流センサ５１で測定された出力電流の値Ｉにフーリエ変換処理を施し電圧振幅値Ｉ（ω_H）を得て、これらの比Ｖ（ω_H）／Ｉ（ω_H）をインピーダンスＺ（ω_H）として求める。なお、インピーダンスＺ（ω_H）を計測する手法は、電解質膜応答周波数帯以外のアノード極応答周波数帯や低周波数帯から選択した周波数に対して行う場合も同様であるので、以降はその詳細な説明は省略する。

次に、ステップＳ１０２において、コントローラ６は、得られたインピーダンスＺ（ω_H）から電解質膜抵抗値Ｒ_mを推定する。具体的には、上述のように、電解質膜応答周波数帯はいわゆるＨＦＲ計測にて用いられる周波数帯であるところ、当該高周波数帯から選択された周波数ω_Hに基づくインピーダンスＺ（ω_H）又はその実部成分Ｚ_r（ω_H）は、ほぼ電解質膜抵抗値Ｒ_mに一致することとなる。すなわち、インピーダンスＺ（ω_H）又はその実部成分Ｚ_r（ω_H）の値をそのまま電解質膜抵抗値Ｒ_mと推定する。

ステップＳ１０３において、コントローラ６は、アノード極応答周波数帯における２点の周波数ω₁、ω₂を選択し、当該周波数ω₁、ω₂に基づくアノード極応答インピーダンスＺ（ω₁）、Ｚ（ω₂）を求める。

ステップＳ１０４において、コントローラ６は、推定された電解質膜抵抗値Ｒ_m及び得られた２つのインピーダンスＺ（ω₁）、Ｚ（ω₂）から、アノード極１１２の反応抵抗値Ｒ_a及びアノード極１１２の電気二重層容量値Ｃ_aを推定する。

この推定の態様を具体的に説明する。先ず、アノード極応答周波数帯における２点の周波数ω₁、ω₂を選択する場合には、上述のようにカソード極１１３の反応抵抗を無視することができ、したがって、インピーダンスの式として、簡易等価回路に基づくインピーダンスの式（１）からカソード極１１３の反応抵抗値Ｒ_aを除いた式（２）を用いることができる。

ここで、式（２）において、既知の値である２点の周波数ω₁、ω₂及びこれらに基づくインピーダンスＺ（ω₁）、及びＺ（ω₂）の組み合わせを代入し、その実部Ｚ_r（ω₁）及びＺ_r（ω_m2）をとる。そして、推定された電解質膜抵抗値Ｒ_mが既知であることを考慮すると、未知数をＲ_a及びＣ_aとする２つの方程式が得られることとなる。したがって、得られた２つの方程式を解けばＲ_a及びＣ_aを求めることができる。

未知数Ｒ_a、及びＣ_aを求める方法の一例を示す。先ず、式（２）の実部をとって変形すると、

となる。横軸ω²、縦軸１／Ｚ_rとする平面を考慮すると、当該平面において式（３）は直線を表し、その傾きｍ_rは

で与えられる。ここで、２点の周波数ω₁、ω₂が既知であるので、この２点の周波数ω₁、ω₂及びこれらに対応するインピーダンス計測値の実部Ｚ_r（ω₁）及びＺ_r（ω₂）を上記平面にプロットすると、これらを結ぶ直線が定まり、傾きｍ_rの値が決まる。すなわち、式（４）の未知数はＲ_a及びＣ_aである。

次に、式（３）が表す直線の切片ａは、

で与えられる。ここで、傾きｍ_rの値と同様に点の周波数ω₁、ω₂及びこれらに対応するインピーダンス計測値の実部Ｚ_r1及びＺ_r2により切片ａの値も決まる。そして、Ｚ_rはインピーダンス計測値の実部Ｚ_r1及びＺ_r2に相当するので、式（５）の未知数はＲ_aのみである。

したがって、式（５）によれば、アノード極１１２の反応抵抗値Ｒ_aを、

として求めることができる。

また、式（６）により定まったＲ_aを式（４）に代入することで、アノード極１１２の電気二重層容量値Ｃ_aを、

として求めることができる。

なお、Ｒ_a、及びＣ_aを求める計算方法は、上記した計算方法に限られるものではなく、種々の適切な計算方法を用いることができる。

次に、ステップＳ１０５において、コントローラ６は、低周波数帯における１点の周波数ω_Lを選択し、当該周波数ω_Lに基づくインピーダンスＺ（ω_L）を計測する。

ステップＳ１０６において、コントローラ６は、既に推定されている電解質膜抵抗値Ｒ_m、アノード極１１２の反応抵抗値Ｒ_a、及びアノード極１１２の電気二重層容量値Ｃ_a、並びに計測されたインピーダンスＺ（ω_L）を用いて、カソード極１１３の電気二重層容量値Ｃ_cを推定する。

この推定の態様を具体的に説明する。上述のように低周波数帯における周波数ω_Lの交流電流は、燃料電池１の簡易等価回路における全ての回路要素、すなわちアノード極１１２の反応抵抗及び電気二重層容量、電解質膜抵抗、並びにカソード極１１３の反応抵抗及び電気二重層容量の部分を流れる。したがって、周波数ω_Lに基づき得られる低周波数インピーダンスＺ（ω_L）では、アノード極１１２の反応抵抗値Ｒ_a及び電気二重層容量値Ｃ_a、電解質膜抵抗値Ｒ_m、並びにカソード極１１３の反応抵抗値Ｒ_c及び電気二重層容量値Ｃ_cの情報が含まれる。したがって、インピーダンスの式としては、上記全ての回路要素が考慮されている式（１）を用いる必要がある。

式（１）において、既知の値である周波数ω_L及びこれに基づくインピーダンスＺ（ω_L）を代入し、その実部Ｚ_r（ω_L）及び虚部Ｚ_i（ω_L）をとる。そして、推定された電解質膜抵抗値Ｒ_m、アノード極１１２の反応抵抗値Ｒ_a、及びアノード極１１２の電気二重層容量値Ｃ_aが既知であることを考慮すると、未知数をＲ_c及びＣ_cとする２つの方程式が得られることとなる。したがって、この２つの方程式を解けば未知数をＲ_c及びＣ_cを求めることができる。

このように未知数Ｒ_c及びＣ_cを求める方法の一例を示す。先ず、式（１）の実部をとって変形すると、

となる。

また、式（１）の虚部をとって変形すると、

となる。

ここで、周波数ω_L、周波数ω_Lに対応するインピーダンス計測値の実部Ｚ_r（ω_L）及び虚部Ｚ_i（ω_L）、並びにＲ_a及びＣ_aが既知であるので、これらを式（８）及び式（９）に代入して変形すると、カソード極１１３の電気二重層容量値Ｃ_cは、

ただし、式（１０）中、ωはω_Lであり、Ａは、下記の式（１１）のように定義される。

さらに、カソード極１１３の反応抵抗値Ｒ_cは、

として求められる。

ただし、式（１２）中のＡは上記式（１１）のように定義され、式（１２）中のＢは下記式（１３）のように定義される。

以上のように、ステップＳ１０１〜ステップＳ１０６の工程により、燃料電池１の状態量として、電解質膜抵抗値Ｒ_m、アノード極１１２の反応抵抗値Ｒ_a、アノード極１１２の電気二重層容量値Ｃ_a、カソード極１１３の反応抵抗値Ｒ_c、及びカソード極１１３の電気二重層容量値Ｃ_cが推定されることとなる。

上記した本実施形態によれば、以下の効果を得ることができる。本実施の形態では、コントローラ６、電流センサ５１、電圧センサ５２、及びＤＣ／ＤＣコンバータ５６により状態検出装置が構成される。また、インピーダンス取得手段及び内部状態量推定手段は、コントローラ６により構成される。

本実施形態によれば、アノードガス及びカソードガスの供給を受けて発電する燃料電池１の状態検出装置のインピーダンス取得手段は、少なくともアノード極１１２の状態量Ｒ_a、Ｃ_aに対して応答性を示す周波数帯を含む高周波数帯（アノード極応答周波数帯及び電解質膜応答周波数帯）から選択された周波数ω_H、ω₁、ω₂に基づく高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）、及び少なくともカソード極の状態量Ｒ_c、Ｃ_cに対して応答性を示す周波数帯を含む低周波数帯から選択された周波数ω_Lに基づく低周波数インピーダンスＺ（ω_L）を取得する（ステップＳ１０１、ステップＳ１０３、ステップＳ１０５）。

そして、燃料電池１の状態検出装置の内部状態量推定手段は、取得された高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）及び低周波数インピーダンスＺ（ω_L）を組み合わせて燃料電池１の内部状態としてのアノード極１１２の状態量Ｒ_a、Ｃ_aとカソード極１１３の状態量Ｒ_c、Ｃ_cをそれぞれ推定する。

これによれば、取得された高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）及び低周波数インピーダンスＺ（ω_L）という異なる周波数帯から得られるインピーダンス情報に基づき、周波数の大小に応じた電流変動に対するアノード極１１２の反応とカソード極１１３の反応の追従速度差を利用して、少なくともアノード極１１２の状態量Ｒ_a、Ｃ_aとカソード極１１３の状態量Ｒ_c、Ｃ_cをそれぞれ個別に検知することができる。したがって、高精度なアノード極１１２の状態量Ｒ_a、Ｃ_aとカソード極１１３の状態量（Ｒ_c、Ｃ_c）の情報を得ることができ、結果としてこれらの状態量を利用して行われる燃料電池１の動作制御をより適切なものとすることができる。

さらに、本実施形態によれば、内部状態量推定手段は、高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）に基づいてある内部状態量Ｒ_m、Ｒ_a、Ｃ_aを推定し、推定された該内部状態量Ｒ_m、Ｒ_a、Ｃ_a及び低周波数インピーダンスＺ（ω_L）に基づいて他の内部状態量Ｒ_c、Ｃ_cを推定している。

これにより、一つの周波数帯である低周波数帯の低周波数インピーダンスＺ（ω_L）だけでは確定できない内部状態量Ｒ_c、Ｃ_cを、他の周波数帯である高周波数帯の高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）により推定された内部状態量Ｒ_m、Ｒ_a、Ｃ_aに基づいて確定させることができる。すなわち、複数種類の内部状態量Ｒ_m、Ｒ_a、Ｃ_a、Ｒ_c、Ｃ_cにおけるそれぞれの区別をより確実に行うことができる。

なお、逆に、内部状態量推定手段は、低周波数インピーダンスＺ（ω_L）に基づいてある内部状態量を推定し、推定された内部状態量及び高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）に基づいて他の内部状態量を推定するようにしても良い。

また、本実施形態によれば、上記高周波数帯（アノード極応答周波数帯及び電解質膜応答周波数帯）は、燃料電池１のアノード極１１２の状態量Ｒ_a、Ｃ_aに対して応答性を示す周波数帯であるアノード極応答周波数帯と、アノード極応答周波数帯よりも高い周波数帯であり燃料電池１の電解質膜の状態量Ｒ_mに対して応答性を示す電解質膜応答周波数帯と、を含む。そして、インピーダンス取得手段は、高周波数インピーダンスＺ（ω_H）、Ｚ（ω₁）、Ｚ（ω₂）として、アノード極応答周波数帯から選択された周波数に基づくアノード極応答インピーダンスＺ（ω₁）、Ｚ（ω₂）及び電解質膜応答周波数帯から選択された周波数に基づく電解質膜応答インピーダンスＺ（ω_H）の双方を取得している（ステップＳ１０１、ステップＳ１０３）。

これにより、燃料電池１の電解質膜１１１の状態量Ｒ_m及びアノード極１１２の状態量Ｒ_a、Ｃ_aをそれぞれ、電解質膜応答インピーダンスＺ（ω_H）及びアノード極応答インピーダンスＺ（ω₁）、Ｚ（ω₂）に基づいて推定することができる。

さらに、本実施形態によれば、内部状態量推定手段は、電解質膜応答インピーダンスＺ（ω_H）に基づいて電解質膜１１１の状態量Ｒ_mを推定し（ステップＳ１０２）、推定された該電解質膜１１１の状態量Ｒ_m、及びアノード極応答周波数帯インピーダンスＺ（ω₁）、Ｚ（ω₂）に基づいてアノード極１１２の状態量Ｒ_a、Ｃ_aを推定する（ステップＳ１０４）。

これにより、推定された該電解質膜１１１の状態量Ｒ_m、及びアノード極応答周波数帯インピーダンスＺ（ω₁）、Ｚ（ω₂）に基づいて、アノード極１１２の状態量Ｒ_a、Ｃ_aを他の状態量からより確実に切り分けて推定することができる。

特に、本実施形態では、アノード極１１２の状態量Ｒ_a、Ｃ_aには、該アノード極１１２の反応抵抗値Ｒ_a及び電気二重層容量値Ｃ_aが含まれており、カソード極１１３の状態量Ｒ_c、Ｃ_cには、該カソード極１１３の反応抵抗値Ｒ_c及び電気二重層容量値Ｃ_cが含まれている。そして、内部状態量推定手段は、アノード極応答インピーダンスＺ（ω₁）、Ｚ（ω₂）に基づいて、アノード極１１２の反応抵抗値Ｒ_a、及びアノード極１１２の電気二重層容量値Ｃ_aを推定する（ステップＳ１０４）。また、内部状態量推定手段は、推定された電解質膜１１１の状態量Ｒ_m、アノード極１１２の反応抵抗値Ｒ_a、及びアノード極１１２の電気二重層容量値Ｃ_a、並びに低周波数インピーダンスＺ（ω_L）に基づいてカソード極１１３の反応抵抗値Ｒ_cを推定する（ステップＳ１０６）。

これによれば、カソード極１１３の反応抵抗値Ｒ_c以外の情報が全て含まれる低周波数帯の低周波数インピーダンスＺ（ω_L）に対して、アノード極応答インピーダンス（Ｚ（ω₁）、Ｚ（ω₂））に基づいて推定されたアノード極１１２の反応抵抗値Ｒ_a及び電気二重層容量値Ｃ_a、電解質膜応答インピーダンスＺ（ω_H）に基づいて推定された電解質膜１１１の状態量Ｒ_mを適用することができる。

したがって、目的の状態量Ｒ_c以外の情報を含む低周波数帯の低周波数インピーダンスＺ（ω_L）から、当該目的の状態量Ｒ_cを好適に切り分けて推定することができる。

（第２実施形態）
以下では、第２実施形態について説明する。なお、既に説明した第１の実施形態の要素と同様の要素には同一の符号を付す。

図６は、第２実施形態に係る状態量の推定の流れを示すフローチャートである。図６におけるステップＳ１０１〜ステップＳ１０４については、図５におけるＳ１０１〜ステップＳ１０４と同様であるので、その詳細な説明を省略する。第２実施形態では、低周波数帯の周波数で低周波数インピーダンスを計測することに代えて、予め設定された燃料電池１のＩ−Ｖ特性線図（Ｉ−Ｖ特性図）における特性線の直線部分の傾きを低周波数インピーダンスとみなして取得するようにしている。

図示のように、ステップＳ１０１〜ステップＳ１０４を経た後、すなわち、アノード極１１２の反応抵抗値Ｒ_a及び電気二重層容量値Ｃ_aの推定値を取得した後、ステップＳ２０５において、燃料電池１のＩ−Ｖ特性図における特性線の直線部分の傾きΔＶ／ΔＩを、低周波数インピーダンスＺ（ω_L）とみなして取得する。

図７には、定常時及び非定常時における燃料電池１のＩ−Ｖ特性線をそれぞれ示す。なお、この燃料電池１のＩ−Ｖ特性線は予め実験等に基づいて定められるものである。特性線Ｃｖ１は、定常時におけるＩ−Ｖ特性を示しており、特性線Ｃｖ２は、非定常時におけるＩ−Ｖ特性を示している。ここで、定常時とは、発進時や停車時等の急加速状態ではない安定走行時における燃料電池１が出力特性を意味する。

特に、図から理解されるように、定常時の特性線Ｃｖ１の定常領域Ｐでは、傾きΔＶ／ΔＩの変動が小さくほぼ一定値をとっており、直線形態をとる。したがって、定常領域Ｐでは出力電流Ｉにかかわらず、傾きΔＶ／ΔＩを一定値とみなすことができる。

このように、ΔＶ／ΔＩの値が一定の定常領域Ｐは、定常時の特性線Ｃｖ１のΔＶ／ΔＩの値が、所定値以下となるように横軸（出力電流Ｉ）の区間である。

本実施の形態では、コントローラ６は、この定常領域ＰにおけるΔＶ／ΔＩの値を予め図示しないメモリ等に記憶させて置き、低周波数インピーダンスＺ（ω_L）を取得するタイミングで当該メモリからΔＶ／ΔＩの値を読み出して、低周波数インピーダンスＺ（ω_L）とみなす。このようにして得られた低周波数インピーダンスＺ（ω_L）は現実の値に良く整合する。

そして、ステップＳ２０６において、低周波数インピーダンスＺ（ω_L）として取得したΔＶ／ΔＩの値を用いてカソード極１１３の反応抵抗値Ｒ_cの推定を行う。

具体的に説明する。上述した式（１）においてωが低周波数（ω→０）であると仮定すると、

が成り立つと考えられる。したがって、式（１４）においてインピーダンスＺをΔＶ／ΔＩに置き換えると、

となる。

これにより、ステップＳ１０１〜ステップＳ１０４の過程で推定された電解質膜抵抗値Ｒ_m、及びアノード極１１２の反応抵抗値Ｒ_aを式（１５）に代入することで、カソード極１１３の反応抵抗値Ｒ_cを算出することができる。

以上、説明した本実施形態にかかる燃料電池１の状態検出装置によれば、インピーダンス取得手段としてのコントローラ６は、燃料電池１のＩ−Ｖ特性線における傾きΔＶ／ΔＩを低周波数インピーダンスＺ（ω₁）として取得する。すなわち、低周波数インピーダンスＺ（ω₁）を直接計測することなく取得することができる。

なお、例えば、Ｉ−Ｖ特性線における傾きΔＶ／ΔＩの値として低周波数インピーダンスＺ（ω₁）を取得すること、及び計測により低周波数インピーダンスＺ（ω₁）を取得することの双方の方法により、低周波数インピーダンスＺ（ω₁）を取得し、これら双方の方法により得られた低周波数インピーダンスＺ（ω₁）相互を比較・補正する等して取得されたより高精度の低周波数インピーダンスＺ（ω₁）をカソード極１１３の反応抵抗値Ｒ_cの推定に用いるようにしても良い。

さらに、本実施の形態では、インピーダンス取得手段としてのコントローラ６は、燃料電池１のＩ−Ｖ特性線Ｃｖ１における傾きの値の変動が所定値以下となる定常領域Ｐにおいて、傾きΔＶ／ΔＩを低周波数インピーダンスＺ（ω₁）として取得している。

このように傾きΔＶ／ΔＩの変動が比較的小さい定常領域Ｐでは、出力電流Ｉの計測値にかかわらず、傾きΔＶ／ΔＩの値を一定とみなして差し支えないので、出力電圧Ｖや出力電流Ｉの計測値ごとに傾きΔＶ／ΔＩの値を算出する必要が無くなり、演算量を低減させることができる。

（第３の実施形態）
第３の実施形態について説明する。なお、既に説明した実施形態の要素と同様の要素には同一の符号を付す。

図８は、本実施形態に係る状態量の推定の流れを示すフローチャートである。図示のように、本実施形態では、図５に示すステップＳ１０１及びステップＳ１０２に相当する電解質膜応答周波数帯の周波数を用いた電解質膜抵抗値Ｒ_mの推定が省略される。

特に本実施の形態では、特有のステップＳ３０４においてアノード極応答周波数帯の２点の周波数ω₁、ω₂において取得したアノード極応答インピーダンスＺ（ω₁）及びＺ（ω₂）を用いて、状態量としてのアノード極１１２の反応抵抗値Ｒ_a、及びアノード極１１２の電気二重層容量値Ｃ_a、カソード極１１３の電気二重層容量値Ｃ_c、及び電解質膜抵抗値Ｒ_mが推定される（ステップＳ３０４）。

以下、ステップＳ３０４における状態量推定の一態様を説明する。

本実施の形態においても、上述したインピーダンスの式（２）に基づいて計算が行われる。式（２）の実部をとって式（３）を得て、式（３）に基づき式（４）を得る工程については、第１の実施形態に係るアノード極１１２の反応抵抗値Ｒ_a及びアノード極１１２の電気二重層容量値Ｃ_aの推定の場合と同様である。

そして、式（４）を変形すれば、

が得られる。なお、ｍ_rは上述のように、２つのインピーダンスＺ（ω₁）及びＺ（ω₂）を結ぶ直線の傾きであって既知の値である。

一方、式（２）の虚部をとると、

ここで、上記式（１７）に式（１６）のＲ_aを代入して両辺にωを乗ずると、

そして、上記既知の周波数ω₁及びω₂、並びにこれに対応するインピーダンス計測値の虚数成分Ｚ_i1及びＺ_i2をそれぞれ、式（１８）に代入して２つの式を得、この２つの式の差をとってカソードの電気二重層容量Ｃ_cを消去すると、未知数であるアノードの電気二重層容量Ｃ_aに関する４次方程式、

が得られる。

式（１９）の４次方程式を解き、Ｃ_aが虚数値を取りえないことを考慮すると、アノードの電気二重層容量Ｃ_aの候補として２つの解、

が得られる。なお、式（１９）の４次方程式は、当業者により知られている種々の解法を用いることができる。

ただし、ｔ₁は下記のように定義される定数である。

以上、本発明の実施形態について説明したが、上記実施形態は本発明の適用例の一部を示したに過ぎず、本発明の技術的範囲を上記実施形態の具体的構成に限定する趣旨ではない。

さらに、式中のＡ₂、Ａ₁、及びＡ₀は、それぞれ、

さらに、Ｃ_a1及びＣ_a2のそれぞれを、上記式（１６）に代入することで、当該Ｃ_a1及びＣ_a2に対応して反応抵抗の推定値の候補としてＲ_a1及びＲ_a2が定まる。推定値の候補Ｒ_a1及びＲ_a2は下記のようになる。

ここで、上述のアノード極１１２の電気二重層容量値の候補Ｃ_a1及びＣ_a2、並びに反応抵抗値の候補Ｒ_a1及びＲ_a2から、現実の特性に適合する真の推定値を決定する必要がある。その方法の一例を説明する。

本実施の形態では、この真の推定値の決定を、Ｃ_a1、Ｒ_a1、Ｃ_a2、及びＲ_a2の値のみから判断するのではなく、上記式（１７）におけるインピーダンス虚部の式を変形して得られるカソード極１１３の電気二重層容量値Ｃ_cの式、

図９には、カソード極１１３の電気二重層容量値の候補Ｃ_c1、Ｃ_c2の周波数応答を示している。なお、当グラフは、予め実験等により算出された周波数ω₁及びω₂をアノード極応答周波数帯の範囲で連続的に変化させて得られた電気二重層容量値の候補Ｃ_c1、Ｃ_c2のデータに基づくものである。

なお、当グラフでは、Ｃ_c1が描く線を破線、Ｃ_c2が描く線を実線で表す。また、周波数ω_dは、アノード極１１２の反応抵抗値及び電気二重層容量値の候補の組（Ｃ_a1、Ｒ_a1）及び（Ｃ_a2、Ｒ_a2）について（Ｃ_a1、Ｒ_a1）＝（Ｃ_a2、Ｒ_a2）となる周波数である。すなわち、周波数ω_dでは、Ｃ_a1、Ｒ_a1、Ｃ_a2、及びＲ_a2における上記式（２０）、（２１）、（２４）、及び（２５）の根号内部が０となる。

図示のように、周波数ω＜ω_dとなる領域においては、電気二重層容量値の推定値候補Ｃ_c2が基本的に０以下の値をとり、ω_dの直前ではＣ_c2の値が周波数変化に対して極端に敏感であることから、周波数ω＜ω_dとなる領域ではＣ_c1が現実的に採用されるべき真の推定値である。

したがって、カソード極１１３の電気二重層容量値及び反応抵抗値についても、周波数ω＜ω_dとなる領域においては、上記Ｃ_c1に対応するＣ_a1及びＲ_a1がそれぞれ採用される。

一方で、ω＞ω_dとなる領域においては、カソード極１１３の電気二重層容量値の候補（Ｃ_c1、Ｃ_c2）の変化を見ただけではＣ_c1及びＣ_c2のいずれを採用すべきか判断が難しい。そこで、この判断は、アノード極１１２の反応抵抗値及び電気二重層容量値の候補の組（Ｃ_a1、Ｒ_a1）及び（Ｃ_a2、Ｒ_a2）を直接検討することにより行う。

図１０Ａは、アノード極１１２の電気二重層容量の候補Ｃ_a1、Ｃ_a2の周波数応答を示している。また、図１０Ｂは、アノード極１１２の反応抵抗値の候補Ｒ_a1、Ｒ_a2の周波数応答を示している。なお、これらグラフも、予め実験等により算出された周波数ω₁及びω₂をアノード極応答周波数帯の範囲で連続的に変化させて得られた候補の組（Ｃ_a1、Ｒ_a1）及び（Ｃ_a2、Ｒ_a2）のデータに基づくものである。

図１０Ａを参照すると、ω＞ω_dとなる領域においては、アノード極１１２の電気二重層容量値の候補Ｃ_a1の周波数に極端に敏感になっている。したがって、ω＞ω_dとなる領域においては、アノード極１１２の電気二重層容量値の真の推定値としてＣ_a2が現実的に採用されるべき値である。したがって、周波数ω＞ω_dとなる領域においては、Ｃ_a2とこれに対応するＲ_a1がそれぞれ採用されるべきである。

なお、図１０Ｂを参照すると理解されるように、周波数ω_dはより小さい、ω＜ω_dの領域では、反応抵抗値の候補Ｒ_a2が周波数変化に対して極端に敏感であることから、反応抵抗値の候補Ｒ_a1が現実的に採用されるべき真の推定値である判断される。したがって、この周波数ω＜ω_dとなる領域においては、Ｒ_a1に対応するＣ_a1及びＲ_a1がそれぞれ採用されるべきであり、この点はカソード極１１３の電気二重層容量値の周波数応答に基づく考察と整合することがわかる。

また、周波数ω＝ω_dのときは、（Ｃ_a1、Ｒ_a1）＝（Ｃ_a2、Ｒ_a2）となるので、これらの候補の組のどちらを真の候補の組として採用しても良い。

以上の考察に基づけば、真の推定値の決定においては、候補の組（Ｃ_a1、Ｒ_a1）及び（Ｃ_a2、Ｒ_a2）の中から決定すべき対象が周波数に応じて変化することがわかる。具体的には、アノード極応答周波数帯の２点の周波数ω₁、ω₂と、周波数ω_dと、の大小に応じて、候補の組（Ｃ_a1、Ｒ_a1）及び（Ｃ_a2、Ｒ_a2）から適切な方を決定する。さらに、決定されたアノード極１１２の電気二重層容量値Ｃ_a及び反応抵抗値Ｒ_aの推定値を式（３）に代入すれば、周波数ωとインピーダンス計測値の実部Ｚ_rが既知であることから、電解質膜抵抗値Ｒ_mが求まる。

このように求められたアノード極１１２の電気二重層容量値Ｃ_a、反応抵抗値Ｒ_a、及び電解質膜抵抗値Ｒ_mの推定値を用いて、以降のステップＳ１０５及びステップＳ１０６が第１実施形態と同様に行われ、カソード極１１３の反応抵抗値Ｒ_cも推定される。

以上、説明した本実施形態に係る燃料電池１の状態判定によれば、インピーダンス取得手段及び内部状態量推定手段としてのコントローラ６により、高周波数インピーダンスとして、アノード極応答インピーダンスＺ（ω₁）及びＺ（ω₂）のみが取得され、アノード極応答周波数帯インピーダンスＺ（ω₁）及びＺ（ω₂）に基づいてアノード極１１２の状態量Ｃ_a及びＲ_aが推定される。

これにより、電解質膜応答インピーダンスの計測に基づく電解質膜抵抗値Ｒ_mの推定を省略してコントローラ６に対する負荷を軽減しつつも、アノード極１１２の状態量Ｃ_a及びＲ_aを推定することができ、最終的にカソード極１１３の状態量である反応抵抗値Ｒ_cも推定することができる。

（第４の実施形態）
第４の実施形態について説明する。なお、既に説明した実施形態の要素と同様の要素には同一の符号を付す。

図１１は、本実施形態に係る状態量の推定の流れを示すフローチャートである。図示のように、本実施形態では、第３実施形態と同様に、ステップＳ１０３においてアノード極応答インピーダンスＺ（ω₁）、Ｚ（ω₂）を求め、ステップＳ３０４においてアノード極１１２の反応抵抗値Ｒ_a及び電気二重層容量値Ｃ_a、カソード１１３の電気二重層容量値Ｃ_c、及び電解質膜抵抗値Ｒ_mの推定値を求める。

その後、第２実施形態の場合と同様に、ステップＳ２０５において燃料電池１のＩ−Ｖ特性に基づいて低周波数インピーダンスΔＶ／ΔＩを取得して、ステップＳ２０６においてこのように取得した低周波数インピーダンスΔＶ／ΔＩ及び電解質膜抵抗値Ｒ_mの推定値からカソード極１１３の反応抵抗値Ｒ_cが推定される。

したがって、本実施形態にかかる燃料電池１の状態判定によれば、低周波数インピーダンスＺ（ω₁）を直接計測することなく推定することができるとともに、電解質膜応答インピーダンスの計測に基づく電解質膜抵抗値Ｒ_mの推定を省略することができるので、コントローラ６に対する負荷をより一層軽減することができる。

（第５の実施形態）
第５の実施形態について説明する。なお、既に説明した実施形態の要素と同様の要素には同一の符号を付す。

本実施の形態では、第２実施形態及び第４実施形態に係るステップＳ２０５において、図７の定常時の特性線Ｃｖ１の定常領域ＰにおけるΔＶ／ΔＩの値をメモリしておく態様に代えて、ΔＶ／ΔＩの値を算出するために実際の出力電圧Ｖ及び出力電流Ｉの計測値を用いる。

図１２は、定常時における燃料電池１のＩ−Ｖ特性線を示している。特に本実施の形態では、所定の計測タイミングにおいて電流センサ５１で測定された出力電流Ｉ₁、Ｉ₂と、同タイミングにおいて電圧センサ５２で測定された出力電圧Ｖ₁、Ｖ₂に対して、−（Ｖ₁−Ｖ₂）／（Ｉ₁−Ｉ₂）を計算することで傾きΔＶ／ΔＩが計算される。

すなわち、出力電圧及び出力電流の計測値に応じて、低周波数インピーダンスとみなす傾きΔＶ／ΔＩが定められることとなる。

本実施の形態では、このように燃料電池１のＩ−Ｖ特性線における傾きΔＶ／ΔＩが、電流及び電圧の２組の計測値（Ｉ₁、Ｖ₁）、（Ｉ₂、Ｖ₂）に基づいて算出される。これにより、定常領域Ｐにおいて一定値とみなして定めた傾きΔＶ／ΔＩを用いる場合と比較して、実際の特性をより高精度に反映したΔＶ／ΔＩの値を得ることができる。結果として、このΔＶ／ΔＩの値を低周波数インピーダンスとみなして算出されるカソード極１１３の反応抵抗値Ｒ_cの推定値の精度も向上することとなる。

（第６の実施形態）
第６の実施形態について説明する。なお、既に説明した実施形態の要素と同様の要素には同一の符号を付す。

本実施の形態では、Ｉ−Ｖ特性線における傾きΔＶ／ΔＩを求めるために、第５実施形態のように出力電流及び出力電圧の２組の計測値（Ｉ₁、Ｖ₁）、（Ｉ₂、Ｖ₂）を計測することに代えて、出力電流及び出力電圧の一つの計測値（Ｉ₃、Ｖ₃）と、事前に設定した１点（Ｉ_set、Ｖ_set）を用いてＩ−Ｖ特性線における傾きΔＶ／ΔＩの算出を行う。

図１３は、Ｉ−Ｖ特性線における傾きΔＶ／ΔＩの算出を行うための１組の電流及び電圧の設定方法の一例を説明する図である。なお、当図では図面の明確化のため、定常時の特性線Ｃｖ１を破線で示している。図示のように、本実施形態では、図の黒塗り四角形で示した点が、上記（Ｉ_set、Ｖ_set）に該当する。特に、Ｉ_set＝０である。

したがって、上記計測値（Ｉ₃、Ｖ₃）及び事前設定値（Ｉ_set、Ｖ_set）により、−（Ｖ_set−Ｖ₃）／（Ｉ_set−Ｉ₃）を計算することで傾きΔＶ／ΔＩの値が算出される。

上述のように本実施の形態によれば、Ｉ−Ｖ特性線における傾きΔＶ／ΔＩの値は、電流及び電圧の１組の計測値（Ｉ₃、Ｖ₃）及び事前に設定した１組の電流及び電圧の値（Ｉ_set、Ｖ_set）に基づいて算出される。

したがって、燃料電池１のＩ−Ｖ特性線における傾きΔＶ／ΔＩにあたり、当該傾きの値を算出するのに用いるＩ−Ｖ特性線上の２点のうち、一点を事前に設定された（Ｉ_set、Ｖ_set）を用いて演算量を抑えつつ、もう一点に計測値（Ｉ₃、Ｖ₃）を用いることで計算の精度も一定以上に確保することができる。

（第７の実施形態）
以下、第７の実施形態について説明する。なお、既に説明した実施形態の要素と同様の要素には同一の符号を付す。

本実施形態では、第１実施形態等において行われる燃料電池１のインピーダンスの計測にあたり、交流信号を重畳した出力電流Ｉ及び出力電圧Ｖを測定する構成に代えて、燃料電池１に所定の測定用電流源から電流Ｉを供給し、当該供給電流Ｉと出力される電圧Ｖとに基づいてインピーダンスＺ＝Ｖ／Ｉを算出するいわゆる励起電流印加法が行われる。

図１４は、本実施形態の燃料電池システム１００において、インピーダンス計測に係る要部を概略的に示したブロック図である。

図示のように、本実施形態に係る燃料電池システム１００では、燃料電池１に交流電流を調整しつつ印加する印加交流電流調整部２００が設けられている。

印加交流電流調整部２００は、スタックとして構成された燃料電池１の正極端子（カソード極側端子）１Ｂ及び負極端子（アノード極側端子）１Ａの他に、中途端子１Ｃに接続されている。なお、中途端子１Ｃに接続された部分は図に示すようにアースされている。

そして、印加交流電流調整部２００は、中途端子１Ｃに対する正極端子１Ｂの正極側交流電位差Ｖ１を測定する正極側電圧測定センサ２１０と、中途端子１Ｃに対する負極端子１Ａの負極側交流電位差Ｖ２を測定する負極側電圧測定センサ２１２と、を有している。

さらに、印加交流電流調整部２００は、正極端子１Ｂと中途端子１Ｃからなる回路に交流電流Ｉ１を印加する正極側交流電源部２１４と、負極端子１Ａと中途端子１Ｃからなる回路に交流電流Ｉ２を印加する負極側交流電源部２１６と、これら交流電流Ｉ１及び交流電流Ｉ２の振幅や位相を調整するコントローラ２１８と、正極側交流電位差Ｖ１、Ｖ２及び交流電流Ｉ１、Ｉ２に基づいて燃料電池１のインピーダンスＺの演算を行う演算部２２０と、を有している。

本実施形態では、コントローラ２１８は、正極側交流電位差Ｖ１と負極側交流電位差Ｖ２が等しくなるように、交流電流Ｉ１と交流電流Ｉ２の振幅及び位相を調節する。なお、このコントローラ２１８は、図３に示すコントローラ６により構成されても良い。

また、演算部２２０は、図示しないＡＤ変換器やマイコンチップ等のハードウェア、及びインピーダンスを算出するプログラム等のソフトウェア構成を含み、正極側交流電位差Ｖ１を交流電流Ｉ１で除して、中途端子１Ｃから正極端子１ＢまでのインピーダンスＺ１を算出し、負極側交流電位差Ｖ２を交流電流Ｉ２で除して、中途端子１Ｃから負極端子１ＡまでのインピーダンスＺ２を算出する。さらに、演算部２２０は、インピーダンスＺ１とインピーダンスＺ２の和をとることで、燃料電池１の全インピーダンスＺを算出する。

上記した本実施形態に係る燃料電池の状態推定装置によれば、以下の効果を得ることができる。

本実施形態に係る燃料電池の状態推定装置は、燃料電池１に接続されて、該燃料電池１に交流電流Ｉ１，Ｉ２を出力する交流電源部２１４，２１６と、燃料電池１の正極側１Ｂの電位から中途部分１Ｃの電位を引いて求めた電位差である正極側交流電位差Ｖ１と、燃料電池１の負極側１Ａの電位から中途部分１Ｃの電位を引いて求めた電位差である負極側交流電位差Ｖ２と、に基づいて交流電流Ｉ１，Ｉ２を調整する交流調整部としてのコントローラ２１８と、調整された交流電流Ｉ１，Ｉ２並びに正極側交流電位差Ｖ１及び負極側交流電位差Ｖ２に基づいて燃料電池１のインピーダンスＺを演算するインピーダンス演算部２２０と、を有する。

コントローラ２１８は、燃料電池１の正極側の正極側交流電位差Ｖ１が負極側の負極側交流電位差Ｖ２と実質的に一致するように、正極側交流電源部２１４により印加される交流電流Ｉ１及び負極側交流電源部２１６により印加される交流電流Ｉ２の振幅及び位相を調節する。これにより、正極側交流電位差Ｖ１の振幅と負極側交流電位差Ｖ２の振幅とが等しくなるので、正極端子１Ｂと負極端子１Ａが実質的に等電位となる。したがって、インピーダンス計測のための交流電流Ｉ１、Ｉ２が負荷５３に流れることが防止されるので、燃料電池１による発電に影響を与えることが防止される。

また、燃料電池１が発電状態の場合に上記インピーダンス計測を実行する場合、当該発電により生じた電圧に計測用交流電位が重畳されることとなるので、正極側交流電位差Ｖ１及び負極側交流電位差Ｖ２の値自体が大きくなるが、正極側交流電位差Ｖ１及び負極側交流電位差Ｖ２の位相や振幅自体が変わるわけではないので、燃料電池１が発電状態ではない場合と同様に高精度なインピーダンス計測を実行することができる。

以上、本発明の実施形態について説明したが、上記実施形態は本発明の適用例の一部を示したに過ぎず、本発明の技術的範囲を上記実施形態の具体的構成に限定する趣旨ではない。例えば、各実施形態における、アノード極応答インピーダンス、電解質膜応答インピーダンス、及び低周波数インピーダンスを取得する工程（ステップＳ１０１、ステップＳ１０３、及びステップＳ１０５）等は、各実施形態において説明した工程順に限定されることなく、任意に変更することが可能である。

例えば、アノード極応答インピーダンス、電解質膜応答インピーダンス、及び低周波数インピーダンスを取得する工程を全て行った後に、各状態量の推定を行うようにしても良い。

また、燃料電池１において複数の内部状態量を推定する態様は、上記各実施の形態に説明した態様のみに限られない。

例えば、第１実施形態や第３実施形態におけるステップＳ１０５における低周波数帯から一つの周波数ω_Lを選択する態様に代えて、低周数帯で２つの周波数ω_L1、ω_L2を選択して低周波数インピーダンスＺ（ω_L1）及びＺ（ω_L2）を求めるようにしても良い。これにより、最終的にカソード極１１３の反応抵抗Ｒ_cだけでなく、カソード極１１３の電気二重層容量Ｃ_aの推定値も求めることができる。

また、燃料電池１の簡易等価回路の態様も、上記各実施の形態で用いたものに限定されない。例えば、上記各実施形態において説明した各極の反応抵抗や電気二重層容量等の回路素子以外にも、拡散抵抗、電子輸送抵抗、及びアイオノマ抵抗等の他の要素を含む等価回路を設定し、これら他の要素に基づく内部状態量としての拡散抵抗値、電子輸送抵抗値、及びアイオノマ抵抗値等を推定の対象とするようにしても良い。

以下では、各実施形態において、上記電解質膜応答周波数帯、アノード極応答周波数帯、及び低周波数帯のうちの少なくとも２つを用いた各状態量の推定の詳細について説明する。

この推定の態様を具体的に説明する。先ず、アノード極応答周波数帯における２点の周波数ω_１、ω_２を選択する場合には、上述のようにカソード極１１３の反応抵抗を無視することができ、したがって、インピーダンスの式として、簡易等価回路に基づくインピーダンスの式（１）からカソード極１１３の反応抵抗値Ｒ _ｃ を除いた式（２）を用いることができる。

ここで、式（２）において、既知の値である２点の周波数ω_１、ω_２及びこれらに基づくインピーダンスＺ（ω_１）、及びＺ（ω_２）の組み合わせを代入し、その実部Ｚ_ｒ（ω_１）及びＺ_ｒ（ω _２ ）をとる。そして、推定された電解質膜抵抗値Ｒ_ｍが既知であることを考慮すると、未知数をＲ_ａ及びＣ_ａとする２つの方程式が得られることとなる。したがって、得られた２つの方程式を解けばＲ_ａ及びＣ_ａを求めることができる。

Claims (12)

1. アノードガス及びカソードガスの供給を受けて発電する燃料電池の状態検出装置であって、
   少なくともアノード極の状態量に対して応答性を示す周波数帯を含む高周波数帯から選択された周波数に基づく高周波数インピーダンス、及び少なくともカソード極の状態量に対して応答性を示す周波数帯を含む低周波数帯から選択された周波数に基づく低周波数インピーダンスを取得するインピーダンス取得手段と、
   取得された前記高周波数インピーダンス及び前記低周波数インピーダンスを組み合わせて、前記燃料電池の内部状態としての前記アノード極の状態量と前記カソード極の状態量をそれぞれ推定する内部状態量推定手段と、
   を備えた燃料電池の状態検出装置。
2. 請求項１に記載の燃料電池の状態検出装置であって、
   前記内部状態量推定手段は、
   前記高周波数インピーダンスに基づいてある内部状態量を推定し、推定された該内部状態量及び前記低周波数インピーダンスに基づいて他の内部状態量を推定するか、又は
   前記低周波数インピーダンスに基づいてある内部状態量を推定し、推定された該内部状態量及び前記高周波数インピーダンスに基づいて他の内部状態量を推定する燃料電池の状態検出装置。
3. 請求項１又は請求項２に記載の燃料電池の状態検出装置であって、
   前記高周波数帯は、前記燃料電池のアノード極の状態量に対して応答性を示す周波数帯であるアノード極応答周波数帯と、前記アノード極応答周波数帯よりも高い周波数帯であり前記燃料電池の電解質膜の状態量に対して応答性を示す電解質膜応答周波数帯と、を含み、
   前記インピーダンス取得手段は、前記高周波数インピーダンスとして、前記アノード極応答周波数帯から選択された周波数に基づくアノード極応答インピーダンス及び前記電解質膜応答周波数帯から選択された周波数に基づく電解質膜応答インピーダンスの少なくともいずれか一方を取得する燃料電池の状態検出装置。
4. 請求項３に記載の燃料電池の状態検出装置であって、
   前記インピーダンス取得手段は、前記高周波数インピーダンスとして、前記アノード極応答インピーダンス及び前記電解質膜応答インピーダンスの双方を取得し、
   前記内部状態量推定手段は、前記電解質膜応答インピーダンスに基づいて前記電解質膜の状態量を推定し、推定された該電解質膜の状態量、及び前記アノード極応答周波数帯インピーダンスに基づいて前記アノード極の状態量を推定する燃料電池の状態検出装置。
5. 請求項３に記載の燃料電池の状態検出装置であって、
   前記インピーダンス取得手段は、前記高周波数インピーダンスとして、前記アノード極応答インピーダンスのみを取得し、
   前記内部状態量推定手段は、前記アノード極応答周波数帯インピーダンスに基づいて前記アノード極の状態量を推定する燃料電池の状態検出装置。
6. 請求項４に記載の燃料電池の状態検出装置であって、
   前記アノード極の状態量には、該アノード極の反応抵抗値及び電気二重層容量値が含まれ、
   前記カソード極の状態量には、該カソード極の反応抵抗値及び電気二重層容量値が含まれ、
   前記内部状態量推定手段は、
   前記アノード極応答インピーダンスに基づいて、前記アノード極の反応抵抗値及び前記アノード極の電気二重層容量値を推定し、
   推定された前記電解質膜の状態量、前記アノード極の反応抵抗値、前記アノード極の電気二重層容量値、及び前記低周波数インピーダンスに基づいて前記カソード極の反応抵抗値及び電気二重層容量値の少なくともいずれか一方を推定する燃料電池の状態検出装置。
7. 請求項１〜６のいずれか１項に記載の燃料電池の状態検出装置であって、
   前記インピーダンス取得手段は、
   前記低周波数インピーダンスとして前記燃料電池のＩ−Ｖ特性線における傾きの値を取得する燃料電池の状態検出装置。
8. 請求項７に記載の燃料電池の状態検出装置であって、
   前記インピーダンス取得手段は、
   前記燃料電池のＩ−Ｖ特性線における傾きの値の変動が所定値以下となる定常時において、前記傾きの値を前記低周波数インピーダンスとして取得する燃料電池の状態検出装置。
9. 請求項７又は請求項８に記載の燃料電池の状態検出装置であって、
   前記Ｉ−Ｖ特性線における傾きは、電流及び電圧の２組の計測値に基づいて算出される燃料電池の状態検出装置。
10. 請求項７又は請求項８に記載の燃料電池の状態検出装置であって、
    前記Ｉ−Ｖ特性線における傾きは、電流及び電圧の１組の計測値及び事前に設定した１組の電流及び電圧の値に基づいて算出される燃料電池の状態検出装置。
11. 請求項１〜請求項１０のいずれか１項に記載の燃料電池の状態検出装置において、
    前記燃料電池が積層電池として構成され、
    前記積層電池に接続されて該積層電池に交流電流を出力する交流電源部と、
    前記積層電池の正極側の電位から該積層電池の中途部分の電位を引いて求めた電位差である正極側交流電位差と、前記燃料電池の負極側の電位から前記中途部分の電位を引いて求めた電位差である負極側交流電位差と、に基づいて交流電流を調整する交流調整部と、
    前記調整された交流電流並びに前記正極側交流電位差及び前記負極側交流電位差に基づいて前記燃料電池の前記インピーダンス計測値を演算するインピーダンス演算部と、
    を有する燃料電池の状態検出装置。
12. アノードガス及びカソードガスの供給を受けて発電する燃料電池の状態検出方法であって、
    少なくともアノード極の状態量に対して応答性を示す周波数帯を含む高周波数帯から選択された周波数に基づく高周波数インピーダンス、及び少なくともカソード極の状態量に対して応答性を示す周波数帯を含む低周波数帯から選択された周波数に基づく低周波数インピーダンスを取得するインピーダンス取得手段と、
    取得された前記高周波数インピーダンス及び前記低周波数インピーダンスを組み合わせて、前記燃料電池の内部状態としての前記アノード極の状態量と前記カソード極の状態量をそれぞれ推定する内部状態量推定手段と、
    を備えた燃料電池の状態検出方法。

Images