(57)【特許請求の範囲】
【請求項１】
燃料電池において発電に使用されなかった水素を、前記燃料電池へ水素を供給する通路に還流させる水素還流経路を備えた燃料電池システムであって、
前記水素還流経路に設けられるガスパージ弁と、
前記水素還流経路を含む水素循環系内における水素質量流量を推定する水素質量流量推定手段と、
発電電流に応じて前記水素質量流量の閾値を設定する閾値設定手段と、
前記水素質量流量が前記閾値以下になったときに、前記ガスパージ弁を開放するパージ制御手段と、
を含んで構成されたことを特徴とする燃料電池システム。
【請求項２】
前記水素循環系内の混合ガス中の不純物ガス存在量を推定する不純物ガス存在量推定手段を備え、
前記パージ制御手段が、前記水素質量流量が前記閾値以下で、かつ、前記不純物ガス存在量が閾値以上になったとき、前記ガスパージ弁を開放することを特徴とする請求項１記載の燃料電池システム。
【請求項３】
前記水素質量流量推定手段及び／又は不純物ガス存在量推定手段が、
前記水素循環系内の混合ガス中の音速を求めると共に、該音速に基づいて前記混合ガスの平均密度を求め、該平均密度に基づいて前記混合ガス中の水素濃度及び／又は不純物ガス濃度を推定する濃度推定手段を含んで構成され、
該濃度推定手段によって求められる濃度に基づいて水素質量流量及び／又は不純物ガス存在量を推定することを特徴とする請求項１又は２記載の燃料電池システム。
【請求項４】
前記水素質量流量推定手段が、
前記水素循環系内の混合ガス中における超音波の伝播時間に基づいて前記混合ガスの体積流量を検出する体積流量検出手段を含んで構成され、
前記濃度推定手段が、前記体積流量検出手段による伝播時間の計測結果に基づいて音速を求めて、前記水素濃度及び／又は不純物ガス濃度を推定することを特徴とする請求項３記載の燃料電池システム。
【請求項５】
前記水素循環系内の混合ガス中における超音波の伝播時間に基づいて前記混合ガスの体積流量を検出する体積流量検出手段と、
前記混合ガスの圧力を検出する圧力検出手段と、
前記混合ガスの温度を検出する温度検出手段と、
前記超音波の伝播時間に基づいて前記混合ガス中の音速を求めると共に、該音速に基づいて前記混合ガスの平均密度を求め、該平均密度に基づいて前記混合ガス中の水素ガス濃度及び不純物ガス濃度を推定する濃度推定手段と、
を備え、
前記水素質量流量推定手段が、前記体積流量，圧力，温度及び水素ガス濃度に基づいて水素質量流量を推定し、
前記不純物ガス存在量推定手段が、前記圧力，温度及び不純物ガス濃度に基づいて前記不純物ガス存在量を推定することを特徴とする請求項２記載の燃料電池システム。
【請求項６】
前記水素還流経路に除湿手段を介装し、
前記濃度推定手段が、前記除湿手段の下流側における超音波の伝播時間に基づいて水素濃度及び／又は不純物ガス濃度を推定することを特徴とする請求項３〜５のいずれか１つに記載の燃料電池システム。
【請求項７】
前記超音波の伝播時間の計測部位近傍に、前記混合ガスの湿度を検出する湿度検出手段を設け、
前記濃度推定手段が、前記湿度に基づいて求められる水蒸気濃度と前記平均密度とに基づいて、前記水素ガス濃度及び／又は不純物ガス濃度を推定することを特徴とする請求項３〜５のいずれか１つに記載の燃料電池システム。
【請求項８】
前記超音波の伝播時間の計測部位近傍に、前記混合ガスの圧力を検出する圧力検出手段及び混合ガスの温度を検出する温度検出手段を設け、
前記濃度推定手段が、飽和水蒸気状態での水蒸気濃度を前記圧力及び温度から求め、該水蒸気濃度と前記平均密度とに基づいて、前記水素ガス濃度及び／又は不純物ガス濃度を推定することを特徴とする請求項３〜５のいずれか１つに記載の燃料電池システム。
【請求項９】
前記濃度推定手段が、前記混合ガスの水蒸気濃度を予め定められた一定値とし、該水蒸気濃度と前記平均密度とに基づいて、前記水素ガス濃度及び／又は不純物ガス濃度を推定することを特徴とする請求項３〜５のいずれか１つに記載の燃料電池システム。
【請求項１０】
燃料電池において発電に使用されなかった水素を、前記燃料電池へ水素を供給する通路に還流させる水素還流経路を備えた燃料電池システムにおいて、
前記水素還流経路を含む水素循環系内における水素質量流量を推定し、
前記水素質量流量が発電電流に応じて設定される閾値以下になったときに、前記水素還流経路から混合ガスをパージさせることを特徴とする燃料電池システムの制御方法。
【請求項１１】
燃料電池において発電に使用されなかった水素を、前記燃料電池へ水素を供給する通路に還流させる水素還流経路を備えた燃料電池システムにおいて、
前記水素還流経路を含む水素循環系内における水素質量流量を推定し、
前記水素循環系内の混合ガス中の不純物ガス存在量を推定し、
前記水素質量流量が発電電流に応じて設定される閾値以下で、かつ、前記不純物ガス存在量が閾値以上になったとき、前記水素還流経路から混合ガスをパージさせることを特徴とする燃料電池システムの制御方法。
【請求項１２】
前記水素循環系内の混合ガス中における超音波の伝播時間に基づいて、前記混合ガス中の音速を求め、
該音速に基づいて前記混合ガスの平均密度を求め、
該平均密度に基づいて前記混合ガス中の水素濃度及び／又は不純物ガス濃度を推定し、
前記濃度に基づいて水素質量流量及び／又は不純物ガス存在量を推定することを特徴とする請求項１０又は１１記載の燃料電池システムの制御方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、燃料電池システム及び制御方法に関し、詳しくは、燃料ガスを循環して燃料電池に還流させて再利用するタイプの燃料電池システムにおいて、還流系からのパージを制御する技術に関する。
【０００２】
【従来の技術】
従来から、燃料電池システムにおいて、燃料電池で消費されなかった水素を、燃料電池に還流して再利用を図ることが行われている（特開平１０−０５５８１４号公報及び特開２０００−５８０９２号公報参照）。
また、上記のように燃料ガスを還流させて再利用するシステムの場合、燃料ガスを循環し続けるにしたがって循環ライン内に不純物ガス（例えば燃料電池の発電作用には不必要な窒素ガス）が徐々に蓄積・増加していくことが知られている。
【０００３】
このような不必要なガスを循環させ続けることは循環動力を無駄に使うばかりでなく、燃料電池に供給される燃料中に占める不純物ガスの割合が増加することで燃料電池の出力低下を招くことになる。
そこで、前記特開２０００−０５８０９２号に開示される構成では、燃料ガスの循環ラインから不純物ガスを含む燃料ガスをパージさせるパージ弁を備え、循環ライン中の不純物ガス濃度が一定値を超えたときに前記パージ弁を開放させるようになっていた。
【０００４】
【発明が解決しようとする課題】
ところで、燃料ガスを還流させて再利用する燃料電池システムにおいては、燃料ガス供給源からの水素供給流量をＱ（Ｌ／ｍｉｎ）とすると、このＱ（Ｌ／ｍｉｎ）に対してある比率αの流量αＱ（Ｌ／ｍｉｎ）がエジェクタによって吸込まれて還流され、燃料電池には合計流量である（１＋α）Ｑ（Ｌ／ｍｉｎ）が流入する。
【０００５】
そして、燃料電池にてＱ（Ｌ／ｍｉｎ）の水素が消費され、消費されずに残ったαＱ（Ｌ／ｍｉｎ）が未利用ガスとして燃料電池から排出されて循環するという関係にある。
すなわち、燃料ガス供給源から流入する水素流量Ｑ（Ｌ／ｍｉｎ）が燃料電池で消費され、この流量Ｑ（Ｌ／ｍｉｎ）に依存した量αＱ（Ｌ／ｍｉｎ）が循環している。
【０００６】
一方、特開２０００−０１２０５９号公報に開示されるように、燃料電池に流入する水素質量流量が燃料電池の発電効率に関係しており、燃料電池で消費される水素質量流量以上の量を燃料電池に流入させることが求められ、また、燃料電池の急激な負荷増大に対応するための余裕分も必要なことから、循環される水素質量流量が所定以上であることが要求される。
【０００７】
換言すれば、循環される水素質量流量（余裕分）が所定以上であれば、燃料電池の発電に影響がないため、循環ラインからのパージは必要ないことになるが、従来のように、不純物ガス濃度に基づいてパージを制御する構成では、不純物ガス濃度が同じでも燃料電池における消費量によって循環ラインにおける水素質量流量が変化することに対応できないため、適切なタイミングでパージを行わせることができないという問題があった。
【０００８】
また、従来のように不純物ガス濃度に基づいてパージを制御する構成では、過渡的な不純物ガス濃度の変化に影響されてパージが行われ、無用な燃料ガスの損失を増大させる可能性があった。
例えば、循環ラインの不純物ガスがある量蓄積されている状態で、燃料電池の負荷の増大変化があると、一時的に循環ラインの水素質量が低下し、それによって相対的に不純物ガス濃度が増大するため、従来制御では、係る一時的な不純物ガス濃度の増大に基づいて、パージが行われる可能性がある。
【０００９】
しかし、有効なパージを実行するという観点からすると、不純物ガスの蓄積存在量が所定量以下の場合はパージを控えることが望まれるため、従来制御では、頻繁なパージ制御によって無用な燃料ガスの損失を増大させる可能性があった。本発明は上記問題点に鑑みなされたものであり、燃料ガスの無用な損失を低減させて燃料電池システムのエネルギ効率を高めることが可能な燃料電池システム及び制御方法を提供することを目的とする。
【００１０】
【課題を解決するための手段】
上記目的を達成するため、請求項１，１０記載の発明では、水素還流経路を含む水素循環系内における水素質量流量を推定し、前記水素質量流量が発電電流に応じて設定される閾値以下になったときに、水素還流経路から混合ガスをパージさせる構成とした。上記構成によると、水素循環系内の混合ガス（水素ガス，不純物ガス，水蒸気）の流量変化に応じて変化する水素質量流量が閾値以下となり、発電効率を維持できなくなると、パージを行って燃料ガスと共に不純物ガスをパージさせる。ここで、発電電流に応じて発電効率の維持のために必要とされる水素質量流量が変化することに対応して、水素質量流量の閾値が変更される。
【００１１】
請求項２，１１記載の発明では、水素循環系内における水素質量流量及び不純物ガス存在量をそれぞれ推定し、前記水素質量流量が前記閾値以下で、かつ、前記不純物ガス存在量が閾値以上になったときに、水素還流経路から混合ガスをパージさせる構成とした。上記構成によると、たとえ燃料電池の負荷変化等によって水素質量流量が発電電流に応じて設定される閾値以下になっても、不純物ガス存在量が閾値を下回るときには、パージが行われない。
【００１３】
請求項３，１２記載の発明では、水素質量流量及び／又は不純物ガス存在量の推定に用いる水素濃度及び／又は不純物ガス濃度を、水素循環系内の混合ガス中の音速に基づいて求められる混合ガスの平均密度に基づいて推定させる構成とした。
【００１４】
上記構成によると、混合ガス中における音速から混合ガスの平均密度が分かり、例えば、混合ガスが水素と不純物ガスとからなる場合には、水素濃度と不純物ガス濃度との合計が１で、各ガスの濃度と密度との乗算値を加算した結果が平均密度となるから、これらの式に基づいて水素濃度及び不純物ガス濃度を求めることが可能である。
【００１５】
請求項４記載の発明では、水素質量流量の推定に、水素循環系内の混合ガス中における超音波の伝播時間に基づいて体積流量を検出する手段を用いる構成とし、濃度推定に前記体積流量の検出のための伝播時間の計測結果を用いる構成とした。上記構成によると、体積流量の検出のために計測される伝播時間に基づいて、同時に、水素濃度及び／又は不純物ガス濃度の推定が行われる。
【００１６】
請求項５記載の発明では、超音波の伝播時間に基づいて体積流量を検出する手段と、混合ガスの圧力を検出する手段と、混合ガスの温度を検出する手段と、を備え、前記伝播時間に基づいて混合ガス中の音速を求め、該音速に基づいて求められる平均密度に基づいて水素ガス濃度及び不純物ガス濃度を推定し、前記体積流量，圧力，温度及び水素ガス濃度に基づいて水素質量流量を推定すると共に、前記圧力，温度及び不純物ガス濃度に基づいて不純物ガス存在量を推定する構成とした。
【００１７】
上記構成によると、体積流量の検出に伴って濃度推定が行われ、また、圧力，温度の変化に対応して、水素質量流量及び不純物ガス存在量が推定される。請求項６記載の発明では、除湿手段の下流側での伝播時間に基づいて、水素濃度及び／又は不純物ガス濃度を推定させる構成とした。上記構成によると、水素循環系の混合ガスが、水素，不純物ガス，水蒸気からなる場合に、伝播時間の計測部位の上流側で除湿されることで、混合ガスが水素と不純物ガスとの２成分から構成されることになり、伝播時間から水素の密度と不純物ガスの密度との平均が求められることになる。従って、未知数が水素濃度及び不純物ガス濃度のみとなり、濃度の方程式と密度の方程式とから、水素濃度及び不純物ガス濃度が求められることになる。
【００１８】
請求項７記載の発明では、伝播時間の計測部位近傍で湿度を検出し、該湿度に基づいて求められる水蒸気濃度と伝播時間から求められる平均密度とに基づいて、水素ガス濃度及び／又は不純物ガス濃度を推定する構成とした。上記構成によると、平均密度が、水素，不純物ガス，水蒸気の平均密度として求められる構成であっても、水蒸気の濃度を既知とすることで、水素濃度及び不純物ガス濃度が求められることになる。
【００１９】
請求項８記載の発明では、伝播時間の計測部位近傍で、圧力及び温度を検出し、飽和水蒸気状態での水蒸気濃度を前記圧力及び温度から求め、該水蒸気濃度と伝播時間から求められる平均密度とに基づいて、水素ガス濃度及び／又は不純物ガス濃度を推定する構成とした。上記構成によると、水素循環系内が略飽和水蒸気状態であれば、圧力及び温度から水蒸気濃度が求められ、水蒸気の濃度を既知とすることで、水素濃度及び不純物ガス濃度が求められることになる。
【００２０】
請求項９記載の発明では、混合ガスの水蒸気濃度を予め定められた一定値とし、該水蒸気濃度と伝播時間から求められる平均密度とに基づいて、水素ガス濃度及び／又は不純物ガス濃度を推定する構成とした。上記構成によると、水蒸気の濃度を一定値と仮定することで、水素濃度及び不純物ガス濃度が求められる。
【００２１】
【発明の効果】
請求項１，１０記載の発明によると、発電効率を維持でき、また、負荷変動に対応できるだけの余裕分を確保できる水素質量流量になっているか否かに基づいてパージを行わせることができ、燃料電池の負荷の変化に応じて常に最適なタイミングでパージを行わせることができるという効果がある。
【００２２】
請求項２，１１記載の発明によると、水素質量流量が低下しても、不純物ガス存在量が少ない場合にはパージが行われないので、過渡的に水素質量流量が低下したときに、パージすべき不純物ガス存在量が少ないのにパージが行われてしまうことを回避でき、有効なパージのみを行わせて、より一層無用なパージによる水素ガスの損失を低減させることができるという効果がある。
【００２３】
請求項３，１２記載の発明によると、混合ガス中の音速から水素濃度及び／又は不純物ガス濃度を推定するので、水素質量流量・不純物ガス存在量の推定に必要な濃度情報を音速に基づき得ることができるという効果がある。
【００２４】
請求項４記載の発明によると、水素質量流量を推定するための体積流量の検出に伴って同時に濃度情報を得ることができ、簡便な構成で水素質量流量・不純物ガス存在量の推定を行わせることができるという効果がある。請求項５記載の発明によると、水素質量流量を推定するための体積流量の検出に伴って同時に濃度情報を得ることができ、かつ、水素循環系内の圧力，温度の変化に対応して水素質量流量及び不純物ガス存在量を高精度に推定することができるという効果がある。
【００２５】
請求項６記載の発明によると、混合ガス中の超音波の伝播時間に基づき水素濃度及び／又は不純物ガス濃度を推定させる構成において、除湿後の混合ガスにおける伝播時間を計測させるので、水蒸気ガスの濃度測定を行うことなく、水素濃度及び不純物ガス濃度を推定することができるという効果がある。請求項７記載の発明によると、湿度の検出結果から水蒸気濃度を求めることで、水蒸気を含む混合ガス中における超音波の伝播時間から水素濃度・不純物ガス濃度を推定させることができるという効果がある。
【００２６】
請求項８記載の発明によると、水素循環系内が飽和水蒸気状態であるときに、圧力，温度から簡便に水蒸気濃度を求めて、水蒸気を含む混合ガス中における超音波の伝播時間から水素濃度・不純物ガス濃度を推定させることができるという効果がある。請求項９記載の発明によると、水素循環系内の湿度，圧力，温度が略一定と見なせる条件下であるときに、水蒸気を含む混合ガス中における超音波の伝播時間から水素濃度・不純物ガス濃度を簡便に推定させることができるという効果がある。
【００２７】
【発明の実施の形態】
以下に本発明の実施の形態を図に基づいて説明する。
図１は、実施形態における燃料電池システムを示す図である。
この図１において、水素ガス供給源１からの水素ガスは、供給調圧バルブ２，水素ガス供給ライン３，エジェクタ４，水素ガス供給ライン５を介して、燃料電池スタック２０の燃料極（図示省略）に導入され、燃料電池スタック２０にて消費されなかった燃料ガスは、水素ガス還流路７，９，１０（水素還流経路）を介し、前記エジェクタ４によって水素ガス供給ライン５に還流される。
【００２８】
上記構成によって、燃料電池スタック２０において発電に使用されなかった水素は、再び燃料電池２０に還流され、水素の循環系が水素ガス還流路７，９，１０及び水素ガス供給ライン５によって構成される。
前記水素ガス還流路７から分岐延設される排出ライン１２には、パージ弁１１が介装されている。
【００２９】
前記排出ライン１２の分岐部よりも下流側の水素ガス還流路９には、上流側から順に、除湿器８（除湿手段），温度計３７（温度検出手段），流量計３０（体積流量検出手段），圧力計３６（圧力検出手段）が介装されている。
また、前記水素ガス供給ライン５には、加湿器６が介装されている。
信号演算処理部４０は、前記温度計３７，流量計３０，圧力計３６からの検出信号を入力し、後述する演算処理に基づいて前記パージ弁１１の開放を制御する。
【００３０】
図２は、前記流量計３０の詳細を示す。
前記流量計３０は、本実施形態において超音波式流量計であり、超音波受発信器３１，３２を備え、上流側の超音波受発信器３１から下流側に向けて斜めに発した音波が、超音波反射部３３で反射して下流側の超音波受発信器３２に受信され、下流側の超音波受発信器３２から上流側に向けて斜めに発した音波が、超音波反射部３３で反射して上流側の超音波受発信器３１に受信されるようになっており、下流側に向かう超音波の伝播時間ｔ１と、上流側に向かう超音波の伝播時間ｔ２とが計測されるようになっている。
【００３１】
ここで、水素ガス還流路９の軸線と超音波との角度をθ、ガスの音速をＣ、超音波の伝播経路３４の長さをＬ、水素ガス還流路９での平均流速をＶとすると、前記伝播時間ｔ１，ｔ２（sec）は、
ｔ１＝Ｌ／（Ｃ＋Ｖcosθ）  ……（１）
ｔ２＝Ｌ／（Ｃ−Ｖcosθ）  ……（２）
で表され、上記の２式から、平均流速Ｖ（ｍ/ｓ），音速Ｃ（ｍ/ｓ）を求めると、
Ｖ＝Ｌ／（２cosθ）＊（１／ｔ１−１／ｔ２）  ……（３）
Ｃ＝Ｌ／２＊（１／ｔ１＋１／ｔ２）  ……（４）
となる。
【００３２】
従って、水素ガス還流路９内における平均体積流量Ｑ（ｍ3/s）は、管路断面積をＳ、流量補正係数をＫとすると、
Ｑ＝Ｖ・Ｓ・Ｋ  ……（５）
で求められるが、前記平均体積流量Ｑは、基準温度及び基準圧力状態での体積流量であるため、そのときの圧力ｐ及び温度Ｔに応じた補正を施して、Ｑ’（ノルマル立方メートル／秒）を算出する。
【００３３】
Ｑ’（Ｎｍ3/s）＝Ｑ×（ｐ／ｐ０）×（Ｔ０／Ｔ）  ……（６）
尚、上式で、ｐ０は基準圧力を，Ｔ０は基準温度を示す。
本実施形態では、水素質量流量dmH2/dtに基づいてパージを制御する構成となっており、水素ガス濃度をｘ，水素密度をρｘとすると、水素質量流量dmH2/dt（kg/s）は、
dmH2/dt＝x・ρｘ・Ｑ’  ……（７）
として算出される（水素質量流量推定手段）。
【００３４】
ここで、循環系には、水素，水蒸気，不純物ガス（窒素）からなる混合ガスが流れるが、前記超音波式の流量計３０は除湿器８の下流側に配設されるから、前記音速Ｃは、水蒸気を除く水素と不純物ガス（窒素）との２成分からなる混合ガス中における音速Ｃとなる。
また、音速Ｃを求めた混合ガスの成分が既知である場合には、音速Ｃは、各成分の混合比（濃度）の関数で表され（特開２０００−３０４７３２号参照）、水素と不純物ガス（窒素）の２成分の場合、係る混合ガスの平均密度ρavと音速Ｃとは、図３に示すような相関になる。
【００３５】
従って、混合ガスでの音速Ｃが分かれば、図３に従って混合ガスの平均密度ρavが求まり、平均密度ρavから各成分の濃度を求めることができる。
即ち、水素ガスの濃度をｘ，水素ガスの密度をρｘとし、不純物ガスの濃度をｙ，不純物ガスの密度をρｙとすると、混合ガスは水素ガスと不純物ガスとの２成分からなるから、
ｘ＋ｙ＝１  ……（８）
ｘ・ρｘ＋ｙ・ρｙ＝ρav  ……（９）
であるから、
ｘ＝（ρｙ−ρav）／（ρｙ−ρｘ）  ……（１０）
ｙ＝（ρav−ρｘ）／（ρｙ−ρｘ）  ……（１１）
となり、水素ガスの密度ρｘ及び不純物ガスの密度ρｙと、音速Ｃから求められる平均密度ρavとから、水素ガス濃度ｘ及び不純物ガス濃度ｙを算出することができる（濃度推定手段）。
【００３６】
尚、上記水素ガスの濃度ｘ及び不純物ガス濃度ｙを求めるときに用いる水素ガスの密度ρｘ及び不純物ガスの密度ρｙは、そのときの圧力・温度から設定される。
上記のようにして水素ガス濃度ｘ及び不純物ガス濃度ｙが求められると、係る濃度ｘ，ｙに前記体積流量Ｑ’を乗算することで、水素ガス体積流量Ｑｘ’及び不純物ガス体積流量Ｑｙ’が求められ、該体積流量Ｑｘ’，Ｑｙ’にガス密度ρｘ，ρｙを乗算することで、水素質量流量と不純物ガス質量流量が求められる。
【００３７】
燃料電池に供給する水素質量流量は燃料電池の発電効率に影響を与えるため、燃料電池で消費される水素質量流量以上の量を燃料電池に流入させることが求められ、また、燃料電池の急激な負荷増大に対応するための余裕分も必要なことから、循環ラインにおける水素質量流量、即ち、消費量分に付加される余裕分が所定以上であることが要求される。
【００３８】
一方、還流ガス流量は、エジェクタ４によってある比率の流量が自動的に還流されるようになっているため、還流ガスの中の水素質量流量が上記要求に見合った量よりも少ないと、燃料電池では、必要な水素ガス量が不足し、逆に、不必要なガス成分（不純物ガス）が充満することになる。
そこで、前記信号演算処理部４０では、上記のようにして算出した水素質量流量と閾値とを比較し、水素質量流量が閾値以下になると、燃料電池で必要な水素量が確保できる状態となるまで（例えば一定時間）、パージ弁１１を開いて循環系内に溜まった不純物ガスをパージする（パージ制御手段）。
【００３９】
ここで、燃料電池にとって必要な水素ガス質量流量とは、燃料電池の運転状態に依存しており、燃料電池を変換器と見なすと、その出力である発電電流値に対応した量に設定するのが合理的であり、前記閾値を発電電流に応じて変更すると良い（閾値設定手段）。
図４に示すように、前記水素ガス濃度ｘと体積流量Ｑ’とを変数とすると、水素質量流量一定の状態は双曲線群で表され、ある発電電流値に対応した水素質量流量ａで定常運転する場合を想定すると、図５に示すように、その双曲線以下の領域（斜線領域）は水素質量流量が不足した状態を表すことになり、係る領域に入るとパージ弁１１を開くようにする。
【００４０】
図６のフローチャートは、上記パージ制御の流れを示すものであり、まず、ステップＳ１では、超音波の伝播時間ｔ１，ｔ２の計測を行う。
そして、ステップＳ２，３では、前記伝播時間ｔ１，ｔ２に基づき体積流量Ｑの演算を行い、同時並行で、ステップＳ４〜６では、水素ガス濃度ｘの演算を行う。
【００４１】
ステップＳ２では、前記伝播時間ｔ１，ｔ２に基づき平均流速Ｖを演算し（式（３））、該平均流速Ｖに基づいて体積流量Ｑを演算する（式（５））。
ステップＳ３では、前記体積流量Ｑに圧力，温度に応じた補正を施して、体積流量Ｑ’を求める。
ステップＳ４では、前記伝播時間ｔ１，ｔ２に基づき音速Ｃを演算する（式（４））。
【００４２】
ステップＳ５では、前記音速Ｃを平均密度ρavに変換する（図３参照）。
ステップＳ６（濃度推定手段）では、前記平均密度ρavに基づいて水素ガス濃度ｘ（及び不純物ガス濃度ｙ）を演算する（式（１０））。
ステップＳ７（水素質量流量推定手段）では、前記体積流量Ｑ’，水素ガス濃度ｘ及び水素ガス密度ρｘに基づいて水素ガス質量流量を演算する。
【００４３】
ステップＳ８（閾値設定手段）では、そのときの発電電流に応じて前記水素ガス質量流量の閾値を設定する。
ステップＳ９では、ステップＳ７で求めた水素ガス質量流量が、ステップＳ８で設定した閾値以下であるか否かを判別し、水素ガス質量流量が閾値以下であれば、ステップＳ１０（パージ制御手段）へ進んで、パージ弁１１を開放してパージを行わせる。
【００４４】
ところで、定常運転状態では、上記のように、水素ガス質量流量が閾値以下になったときにパージを行わせることで、無用なパージを回避することができるが、特に車両用燃料電池システムでは、燃料電池の負荷変動が大きく、かつ、変動頻度が多い。
そして、負荷変動時には、燃料電池の発電電流に対応して設定される閾値の値にもよるが、水素質量流量が不足した状態に入り込みやすくなり、負荷に連動して頻繁にパージ弁１１を開放されると、パージ弁開放で排出される水素燃料ガスも無用に増加することになってしまう。
【００４５】
上記のような過渡的な水素質量流量の不足に対しては、循環系における不純物ガス存在量が少ない状態では、たとえ上記の水素質量流量不足領域に一時的に在ったとしても、パージ弁を開放しないようにするとよい。
前記不純物ガス存在量は、還流系の管路等で構成される容積Ｖを一定と見なせば、式（１１）で示される不純物ガス濃度ｙと、圧力ｐと、温度Ｔとから不純物ガス存在量に対応した量を推定することができる（不純物ガス存在量推定手段）。
【００４６】
ここで、上記不純物ガス存在量は、流量（Ｌ／ｍｉｎ）ではなく、また、不純物ガス濃度ｙそのものと対応しているのではなくて、不純物ガス濃度ｙを圧力ｐと温度Ｔで補正した値が存在量に対応することになる。
【００４７】
【数１】
0003882667.tif 000002
【００４８】
水素質量流量が閾値以下になったときに、上記のようにして求められる不純物ガス存在量に対応する値が閾値以上であるときに、パージを実行させるようにすれば、負荷変動時に無用にパージが行われることが回避される。
尚、水素質量流量のみに基づくパージ制御での図５に示されるパージ領域に対して、圧力、温度を一定と仮定すると、不純物ガス濃度ｙ＝１−水素ガス濃度であるから、上記の水素質量流量及び不純物ガス存在量に基づくパージでは、図７のように不純物ガス濃度の閾値に相当するｂによっても囲まれる斜線領域に示されるより限定した領域でパージが行われることになる。
【００４９】
図８のフローチャートは、上記の水素質量流量及び不純物ガス存在量に応じたパージ制御を示すものであり、ステップＳ２１〜Ｓ２９の処理は、前記図６のフローチャートのステップＳ１〜Ｓ９と同様に行われる。
ステップＳ２９（不純物ガス存在量推定手段）で、水素質量流量が閾値以下であると判別されると、ステップＳ３０へ進んで、不純物ガス存在量を演算する。
【００５０】
ステップＳ３１では、前記不純物ガス存在量が予め記憶された閾値以上であるか否かを判別し、前記不純物ガス存在量が閾値以上であれば、ステップＳ３２へ進んで、パージを行わせる。
ところで、図１に示されるシステム構成では、除湿器８の下流側で、水素及び不純物ガスの混合ガスの流量を計測させる構成としたが、循環系の混合ガスが除湿されず、水素，水蒸気，不純物ガスの３成分になっている場合、或いは、図９又は図１０に示すように、除湿器８の上流側位置でガス流量を検知する場合でも、水蒸気濃度ｚを既知とすることで、前記伝播時間ｔ１，ｔ２に基づき求められる平均密度ρavから水素ガス濃度ｘ及び不純物ガス濃度ｙを求めることができる。
【００５１】
図９に示す実施形態では、流量計３０の上流側近傍に湿度計３５を介装させてあり、該湿度計３５（湿度検出手段）を備える場合には、該湿度の検出結果と圧力ｐとから、水蒸気濃度ｚを求めることができる。
また、図１０に示す実施形態では、湿度計３５を備えないが、循環系内が飽和水蒸気状態である場合には、圧力ｐと温度Ｔとに基づいて温度Ｔでの飽和水蒸気圧を求めて、水蒸気分圧（濃度）を推定することができる。
【００５２】
ここで、流量計３０は、水素，不純物ガス，水蒸気の混合ガスにおける超音波の伝播時間を計測することになるから、伝播時間ｔ１，ｔ２に基づく平均密度ρavは、３成分混合ガスの平均密度ρavであり、前記式（８），（９）は、水蒸気の密度をρｚとすると、下記のように変更される。
ｘ＋ｙ＋ｚ＝１  ……（１３）
ｘ・ρｘ＋ｙ・ρｙ＋ｚ・ρｚ＝ρav  ……（１４）
ここで、上式を変形して、
ｘ＋ｙ＝１−ｚ  ……（１５）
ｘ・ρｘ＋ｙ・ρｙ＝ρav−ｚ・ρｚ  ……（１６）
として、これを未知数ｘ，ｙについて解くと、
ｘ＝[（１−ｚ）ρｙ＋ｚ・ρｚ−ρav]／（ρｙ−ρｘ）……（１７）
ｙ＝[ρav−ｚ・ρｚ−（１−ｚ）ρｙ]／（ρｙ−ρｘ）……（１８）
となり、水蒸気濃度ｚが既知であれば、上式から水素ガス濃度ｘ及び不純物ガス濃度ｙが求められることになる。
【００５３】
図１１は、超音波流量計３０が水素，不純物ガス，水蒸気の混合ガスにおける伝播時間を計測する構成において、水蒸気濃度ｚを既知として、該水蒸気濃度ｚから水素濃度ｘを演算するパージ制御を示す。
図１１に示すフローチャートは、ステップＳ２５で平均密度ρavを演算すると、次のステップＳ２６Ａ（濃度推定手段）では、湿度計３５に基づき、或いは、循環系内が飽和水蒸気状態である場合に圧力・温度から、水蒸気濃度ｚを演算する。
【００５４】
そして、ステップＳ２６Ｂでは、前記水蒸気濃度ｚを用いて水素ガス濃度ｘを演算する。
上記ステップＳ２６Ａ，ステップＳ２６Ｂ以外の各ステップでは、図８と同様な処理が行われる。
尚、循環系内における湿度，圧力，温度が略一定である場合には、湿度計３５，圧力計３６，温度計３７を設けなくても、水蒸気濃度ｚを予め一定値と見なして、水素ガス濃度ｘ及び不純物ガス濃度ｙを求めることが可能である。
【００５５】
また、流量計３０を図１２に示すように、水素ガス循環系を構成する供給ライン５に介装させる構成としても良い。
この場合、計測される水素質量流量は、エジェクタ４によって還流される分と、燃料電池スタック２０における消費分に対応する供給ライン３からの供給分との合計になるから、発電電流に応じて設定する水素質量流量の閾値を、供給ライン３からの供給分だけ嵩上げすれば、前記実施形態と同様にしてパージ制御を行わせることができる。
【図面の簡単な説明】
【図１】除湿器の下流側に流量計を備える実施形態のシステム構成図。
【図２】超音波流量計の伝播時間と諸物理量の関係を説明するための模式図。
【図３】音速演算値と混合ガスの平均密度との関係を示す線図。
【図４】水素質量流量と濃度と混合ガス流量の関係を表す線図。
【図５】定常運転状態での水素質量流量不足領域を表す線図。
【図６】水素質量流量に基づくパージ制御を示すフローチャート。
【図７】水素質量流量と不純物ガス存在量とに基づくパージ領域を表す線図。
【図８】水素質量流量と不純物ガス存在量とに基づくパージ制御を示すフローチャート。
【図９】湿度計を備える実施形態のシステム構成図。
【図１０】湿度計を備えない実施形態のシステム構成図。
【図１１】水蒸気濃度から水素ガス濃度を求めるパージ制御を示すフローチャート。
【図１２】流量計を燃料ガス供給系に配置した実施形態を示すシステム構成図。
【符号の説明】
１…水素ガス供給源
２…供給調圧バルブ
３…水素ガス供給ライン
４…エジェクタ
５…水素ガス供給ライン
６…加湿器
７，９，１０…水素ガス還流路
８…除湿器
１１…パージ弁
１２…排出ライン
３０…流量計
３５…湿度計
３６…圧力計
３７…温度計
４０…信号演算処理部
【図１】
0003882667.tif 000003
【図２】
0003882667.tif 000004
【図３】
0003882667.tif 000005
【図４】
0003882667.tif 000006
【図５】
0003882667.tif 000007
【図６】
0003882667.tif 000008
【図７】
0003882667.tif 000009
【図８】
0003882667.tif 000010
【図９】
0003882667.tif 000011
【図１０】
0003882667.tif 000012
【図１１】
0003882667.tif 000013
【図１２】
0003882667.tif 000014