(57)【要約】      （修正有）【課題】燃料電池負極側回路のパージ制御に適した監視方法及び監視装置を提供する。【解決手段】燃料電池を備える電気化学系の状態を観察する方法において、動作中の燃料電池を表すパラメーターを測定し１００、制御ベクトルを生成し２００、測定ベクトルを生成し３００、非補正時間変化と呼ばれる時間変化を計算し４００、スライディングモードにおける補正項を計算し５００、電気化学系の状態の推定値を計算し７００、測定時間を増加させながらステップを順に繰り返す。【選択図】図２Ａ
2018152326.tif 000002
【特許請求の範囲】
【請求項１】
膜によって互いに分離された負極と正極とを含む少なくとも１つの電気化学セルから形成される燃料電池（２）を備え、さらにパージ信号（Ｑａｏ）に基づいて制御されるパージバルブが設けられた負極流体排出ラインに連結された水素を分配する負極分配回路と、酸素を分配するための正極分配回路とをさらに備える、電気化学系（１）の状態（Ｘ）を観察する方法であって、以下のステップを備える：
ａ）  測定時間（ｔｋ）において、燃料電池の電流（Ｉ）、燃料電池の電圧（Ｖ）、負極分配回路の負極圧力（Ｐａ）、正極分配回路の正極圧力（Ｐｃ）と、燃料電池の温度（Ｔ）とを測定し（１００）；
ｂ）  前記パージ信号（Ｑａｏ）並びに前記電流（Ｉ）、前記温度（Ｔ）及び前記正極圧力（Ｐｃ）の前記測定値に基づいて制御ベクトル（Ｕ）を形成し（２００）；
ｃ）  前記電流（Ｉ）及び前記電圧（Ｖ）の測定値から推定された前記膜のオーム抵抗値（Ｒｍ）に基づいて、さらに負極圧力（Ｐａ）の測定値に基づいて測定ベクトル（Ｙ）を形成し（３００）；
ｄ）  非補正時間変化（
【数１】
2018152326.tif 000124
）を計算するところ（４００）、状態（Ｘ）の事前推定値（
【数２】
2018152326.tif 000125
）及び制御ベクトル（Ｕ）に基づいており、電気化学系の状態の所定の表示値（ｆ）に基づいており、状態（Ｘ）は、状態変数であって負極での窒素の量（
【数３】
2018152326.tif 000126
）及び水素の量（
【数４】
2018152326.tif 000127
）を表わすもの、並びに水の量であって前記負極でのもの（
【数５】
2018152326.tif 000128
）、前記正極でのもの（
【数６】
2018152326.tif 000129
）、及び前記膜内のもの（
【数７】
2018152326.tif 000130
）より生成され；
ｅ）  スライディングモードでの補正項（ε）を計算するところ（５００）、前記測定ベクトル（Ｙ）と、前記測定ベクトル（Ｙ）の測定推定値（
【数８】
2018152326.tif 000131
）との間の差分（Ｓ）であって前記状態（Ｘ）の事前推定値（
【数９】
2018152326.tif 000132
）に基づいて得られるものに基づいており；
ｆ）  前記測定時間（ｔｋ）における前記電気化学系（１）の前記状態（Ｘ）の推定値（
【数１０】
2018152326.tif 000133
）を計算するところ（７００）、前記非補正時間変化（
【数１１】
2018152326.tif 000134
）及び前記補正項（ε）に基づいており；
ｇ）  前記測定時間を増加させながらステップａ）からｆ）を繰り返すところ、ステップｄ）及びｅ）における前記先の推定値（
【数１２】
2018152326.tif 000135
）は、前記先行する反復で計算された前記状態推定値（
【数１３】
2018152326.tif 000136
）に対応する、
方法。
【請求項２】
前記差分（Ｓ）が前記測定ベクトル（Ｙ）と前記測定推定値（
【数１４】
2018152326.tif 000137
）との間の差分に等しい、請求項１に記載の方法。
【請求項３】
前記測定推定値（
【数１５】
2018152326.tif 000138
）は前記先行する測定時間（ｔｋ−１）で計算された状態推定値（
【数１６】
2018152326.tif 000139
）に基づいて計算されるところ、前記測定ベクトル（Ｙ）、前記状態（Ｘ）及び前記制御ベクトル（Ｕ）の間の関係を表す所定のモデル（ｈ）に基づいている、請求項１又は２に記載の方法。
【請求項４】
前記補正項（ε）は、所定の利得パラメーター（Ｋ）と前記差分（Ｓ）に適用される符号関数との積に等しい、請求項１から３のいずれかに記載の方法。
【請求項５】
前記状態ベクトル（Ｘ）は、
【数１７】
2018152326.tif 000140
のように定義され、前記制御ベクトル（Ｕ）が
【数１８】
2018152326.tif 000141
のように定義され、前記測定ベクトル（Ｙ）が
【数１９】
2018152326.tif 000142
のように定義される、請求項１から４のいずれかに記載の方法。
【請求項６】
前記状態推定値（
【数２０】
2018152326.tif 000143
）を計算するステップｆ）に状態推定変化項（
【数２１】
2018152326.tif 000144
）を計算するステップが含まれ、当該項は以下のように表される：
【数２２】
2018152326.tif 000145
ここでｆｉ＝１−５は前記状態表示値（ｆ）の項であり、σｉ＝１−５は前記利得パラメーター（Ｋ）の項であり、ｓｉ＝１，２は前記差分（Ｓ）の前記項であり、これらは
【数２３】
2018152326.tif 000146
及び
【数２４】
2018152326.tif 000147
で定義され、ここでｙ１及びｙ２は測定ベクトル（Ｙ）の項であり、
【数２５】
2018152326.tif 000148
及び
【数２６】
2018152326.tif 000149
は測定推定値（
【数２７】
2018152326.tif 000150
）の項である、
請求項１から５のいずれかに記載の方法。
【請求項７】
前記利得パラメーター（Ｋ）は、前記関係式
【数２８】
2018152326.tif 000151
及び
【数２９】
2018152326.tif 000152
が正しいと認められるように決定される、請求項６に記載の方法。
【請求項８】
膜によって互いに分離された負極及び正極を有する少なくとも１つの電気化学セルから形成される燃料電池（２）を備え、さらに水素を分配する負極分配回路であってパージ信号（Ｑａｏ）に基づいて制御されるパージバルブ（４）が設けられた負極液排出ライン（Ｌｅａ）に接続されたものと、酸素を分配する正極分配回路とを備える電気化学系（１）の状態（Ｘ）を観察するための装置（１０）であって、前記観察装置は以下を備えるもの：
・連続する測定時間（ｔｋ）の間、前記燃料電池の電流（Ｉ）の測定値、前記燃料電池の電圧（Ｖ）の測定値、前記負極分配回路の負極圧力（Ｐａ）の測定値、前記正極分配回路の正極圧力（Ｐｃ）の測定値、及び前記燃料電池の温度（Ｔ）の測定値を伝えるのに適したセンサー（２０）；
・各測定時間（ｔｋ）において、請求項１から７のいずれかに記載の観察方法を実行するように構成されたコンピューター（３０）。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は燃料電池を含む電気化学系の状態の監視方法及び監視装置に関する。またこのような監視方法及び装置を含む制御システムに関するものであり、燃料電池の負極配電回路におけるパージを制御するのに適したものである。
【背景技術】
【０００２】
燃料電池は電気化学セルのスタックによって形成されるものであり、各スタックは電解質によって互いに電気的に分離された負極及び正極をそれぞれ含むものであり、これらの電極に連続して導入される２種の反応物の間で電気化学反応が起こるものである。水素燃料電池の場合、燃料（水素）は負極と接触するように導入され、一方、酸化剤（酸素）は例えば空気中に含まれ、正極と接触するように導入される。この反応は、酸化と還元の２つの半反応に細分され、一方は負極／電解質の界面で、他方はと正極／電解質の界面で起こる。反応には２つの電極の間にイオン伝導体、すなわち電解質が存在することを必要としており、さらに外部電気回路によって形成された電子伝導体を必要とする。セルのスタックは、電気化学的反応の起こる部位であり：反応物がそこに導入されなければならず、生成物及び非反応種がそこから排出されなければならず、また生成された熱も排出されなければならない。
【０００３】
燃料電池の出口においては一般的に２つの仕組みが組み合わされる。第１には開放形態と呼ばれる形態での動作であり、そこでは負極の出口と正極の出口が流体排出管と連通する。したがって、燃料電池が理想的な化学量論比を超えて（ｓｕｐｅｒｓｔｏｉｃｈｉｏｍｅｔｒｉｃａｌｌｙ、反応物質の注入量と消費される反応物質の量との比が１より大きい状態で）供給される場合、排出ガスは、反応によって生成された水を排出するのに十分な流量を有し、相分離器を使用して回収される。消費されていない反応物を失うのを避けるために再循環装置が通常設けられ、消費されない反応物が燃料電池のインレットに再注入される。しかしながら、流体の再循環のために、正極に空気が供給される場合の負極における窒素などの非反応種は、モル濃度が徐々に増加することとなるため、定期的にセルをパージする工程を実行する必要がある。
【０００４】
デッドエンド形態と呼ばれる形態での動作に関する第２のケースでは、負極の出口が遮断され、さらにまた、必要に応じて正極の出口が遮断され、燃料電池には化学量論比が実質的に１に等しい状態で反応物が供給され、特に水素に関してはそのように行われる。したがって、負極に導入される水素の量は、反応中に消費される水素の量に対応する。この場合、非反応種及び生成された液体の水は排出されず、それらのモル濃度がセル内の負極側で徐々に増加する。したがって、定期的に実行することが必要となるのはセルをパージするステップであり、これによって非反応種及び液体の水をその場所から排出することになる。
【０００５】
いずれの場合においても、パラメーターの値を認識しておくことは有用であり、それは例えば負極における窒素及び水の量などであり、これによって負極分配回路に連結されたパージバルブの開放を制御できる。しかしながら、これらのパラメーターは、必ずしも直接測定することはできないし、リアルタイムで測定することもできない。したがって、燃料電池の電気化学系の状態が示されたところに基づいて、電気化学系の状態を推定することが可能であり、その状態の変数には負極における窒素及び／又は水の量が含まれており、制御ベクトルと測定ベクトルとに基づいており、その構成要素はシステムにおいて測定可能なパラメーターである。このような観察装置又は状態観察手段は、結果として電気化学系用のソフトウェアセンサーを形作っている。
【０００６】
Ｐｉｆｆａｒｄらによる論文、タイトル「Control by state observer of PEMFC anodic purges in dead-end operating mode, IFAC-PapersOnLine, 48-15(2015)237-243」には、観察装置であって電気化学系の状態の推定値を提供するのに適したものについて述べており、その燃料電池はデッドエンドモードで動作するものである。負極での窒素の量を表す状態変数の推定値に基づいて、負極分配回路のパージのための命令がなされる。
【０００７】
観察装置は、
【数１】
2018152326.tif 000003
のように表される制御ベクトルＵを受け取り、ここで各値はパージ弁の開放（Ｑａｏ＝１）及び閉鎖（Ｑａｏ＝０）に対応するパージ信号Ｑａｏ、冷却回路の熱伝達流体の温度の測定値Ｔ、及び燃料電池によって供給される電流の測定値Ｉである。また、
【数２】
2018152326.tif 000004
のように表される測定ベクトルＹを受け取り、ここでの単一の値は電解質膜のオーム抵抗の値Ｒｍである。制御Ｕベクトル及び測定Ｙベクトルに基づいて、
【数３】
2018152326.tif 000005
のように表される状態の推定値Ｘを提供するところ、その状態変数は、負極における窒素の量
【数４】
2018152326.tif 000006
、負極における水の量
【数５】
2018152326.tif 000007
、及び正極における水の量
【数６】
2018152326.tif 000008
、並びに膜中の水の量λｍを表すものである。状態推定値はアンセンテッドカルマンフィルタ（unscented Kalman filter, UKF）によって作成される。
【０００８】
しかしながら、観察装置の１つの目立った欠点として、燃料電池内の圧力及び温度の時間変化に関して制限したうえで仮定しなければならないことが挙げられ、これらの仮定がＵＫＦアルゴリズムの調整に必要である。
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明の目的は、先行技術の欠点を少なくとも部分的に克服することであり、より詳細には、電気化学系の状態を観察する方法を提供することであり、電気化学系は、少なくとも１つの電気化学セルであって膜によって分離された負極と正極からなるもので形成される燃料電池を有し、さらに水素を分配する負極分配回路であって、パージ信号と呼ばれる信号に基づいて制御されるパージバルブが設けられた負極液排出ラインに連結されたものと、酸素を分配するための正極分配回路とを有する。
【課題を解決するための手段】
【００１０】
この方法は、以下のステップを含む：
ａ）測定時間において以下を測定し、前記燃料電池の電流、前記燃料電池の電圧、前記負極分配回路の負極圧力、前記正極分配回路の正極圧力及び前記燃料電池の温度；
ｂ）前記パージ信号と、前記電流、前記温度、及び前記正極圧力の測定値とに基づいて制御ベクトルを形成し；
ｃ）電流及び電圧の測定値から推定された膜のオーム抵抗値に基づいて、さらに前記負極圧力の測定値に基づいて測定ベクトルを生成し；
ｄ）未補正の時間的変動を計算するところ、状態及び制御ベクトルの事前推定値に基づいており、電気化学系の状態に関して所定の通り示されたものに基づいており、前記状態とは状態変数により形成されるものであり、状態変数とは負極での窒素及び水素の量、並びに負極、正極及び膜における水の量を示すものであり；
ｅ）測定ベクトルと、状態の事前推定に基づいて得られた測定ベクトルの測定推定値との間の差分に基づいて、スライディングモードで補正項を計算し；
ｆ）前記測定時間に対する前記電気化学系の前記状態の推定値を計算するところ、計算は前記補正されていない時間的変動と前記補正項とに基づいており、
ｇ）測定時間を増加させながらステップａ）〜ｆ）を繰り返すところ、ステップｄ）及びｅ）における先の推定値が前の反復で計算された状態推定値に対応する。
【００１１】
この観察方法のいくつかの好ましい態様は以下の通りであるが、これらに限定されない。 
【００１２】
前記差分は、前記測定ベクトルと前記測定推定値との間の前記差分に等しくてもよい。
【００１３】
前記測定推定値は、前回測定時に算出された状態推定値に基づいて算出してもよく、測定ベクトルと、状態と、制御ベクトルとの関係を表す所定のモデルに基づいていてもよい。
【００１４】
補正項は所定の利得パラメーターと前記差分に適用される符号関数との積に等しくてもよい。
【００１５】
状態ベクトルは次のように定義されてもよく、
【数７】
2018152326.tif 000009
制御ベクトルは次のように定義されてもよく、
【数８】
2018152326.tif 000010
測定ベクトルが次のように定義されてもよい。
【数９】
2018152326.tif 000011
【００１６】
状態推定値を計算するステップｆ）は、状態推定変化項と呼ばれる項を計算するステップを含んでいてもよく、その項は次のように記述され：
【数１０】
2018152326.tif 000012
ここで、ｆｉ＝１−５は前記状態を示す前記項であり、σｉ＝１−５は利得パラメーターの前記項であり、ｓｉ＝１，２は前記差分の前記項であり、次のように定義され、
【数１１】
2018152326.tif 000013
及び
【数１２】
2018152326.tif 000014
ここで、ｙ１及びｙ２は前記測定ベクトルの前記項であり、さらに
【数１３】
2018152326.tif 000015
及び
【数１４】
2018152326.tif 000016
は、前記測定推定値の前記項である。
【００１７】
利得パラメーターは、次の関数を満たすように決定されてもよい。
【数１５】
2018152326.tif 000017
及び
【数１６】
2018152326.tif 000018
【００１８】
本発明はまた、膜によって互いに分離された負極及び正極を含む少なくとも１つの電気化学セルから形成される燃料電池を含む電気化学系の状態の観察装置であって、水素を分配する負極分配回路であってパージ信号に基づいて制御されるパージバルブが設けられた負極液排出ラインに連通するものと、酸素を分配する正極分配回路とを備え、
前記観察装置は以下を含む：
−  連続した測定時間に渡って、燃料電池の電流、燃料電池の電圧、負極分配回路の負極圧力、正極分配回路の正極圧力、及び燃料電池の温度の測定値を送るのに適したセンサー；
−  各測定時間において、先の特徴のいずれか１つによる観察方法を実施するように構成されたコンピューター。
【図面の簡単な説明】
【００１９】
本発明の他の態様、目的、利点及び特徴は、その好ましい実施形態の以下の詳細な説明を読むことによって、より明確に明らかになるであろう。その説明は、非限定的な例として、添付の図面を参照して提供される。
【図１】燃料電池電気化学系の部分概略図であり、一実施形態に係る電気化学系を観察するための装置に接続されており、一組の測定センサー及びセンサーからの測定値に基づいて電気化学系の状態を推定するのに適したコンピューターを含む；
【図２Ａ】一実施形態に係る観察方法によって実施される計算ステップを示すフローチャートである、また
【図２Ｂ】補正項を計算するステップを示すフローチャートである；
【図３】実施形態に係る電気化学システム及び観察装置の部分概略図である；
【図４Ａ】図４Ａから図４Ｃは、電気化学システムを動作させる状況における負極（図４Ａ）及び正極（図４Ｂ）における相対湿度の時間的変化の例であり、また負極における窒素の体積割合の時間変化の例であり（図４Ｃ）、特にその中で観察装置によって得られた状態推定値に基づいて負極分配回路のパージのための命令が発せられるものである。
【図４Ｂ】上記の通り。
【図４Ｃ】上記の通り。
【発明を実施するための形態】
【００２０】
図面及び以下の説明において、同じ参照符号は同一又は類似の要素を表す。さらに、様々な要素は、図の明瞭性を高めるために縮尺通りに示されていない。さらに、様々な実施形態及び変形形態は相互に排他的ではなく、互いに組み合わせることができる。他に示されない限り、「実質的に」、「約」及び「その程度の」という用語は、１０％以内を意味する。
【００２１】
図１は、一実施形態に係る観察装置１０の測定センサー２０を備えた燃料電池２と電気化学系１を概略的に示す。観察装置１０は、電気化学系１の状態Ｘの推定値を提供するのに適している。
【００２２】
電気化学系１は、燃料電池２と、燃料電池２の負極放出回路をパージするための装置３とを含む。ここでは特に水素燃料電池が考慮され、この燃料電池では負極側に水素が供給され、正極側に酸素を含む空気が供給される。しかしながら、本発明は、あらゆるタイプの燃料電池、特に低温、すなわち２５０℃未満の温度で動作するものに適用可能である。
【００２３】
燃料電池２は、少なくとも１つの電気化学セル、好ましくは電気化学セルのスタックを含む。各セルは、電解質膜によって互いに分離された負極及び正極を含み、このアセンブリーは膜−電極アセンブリー（ＭＥＡ）を形成している。負極、膜及び正極は、当業者に知られている従来の要素であるので、詳細には説明しない。各アセンブリーはバイポーラプレートによって隣接するセルから分離されるが、これは第１のセルの負極及び隣接するセルの正極の両方に反応種をもたらすのに適しており、また電気化学反応及び非反応種から生じる生成物を除去するのに適しており、またセル間に電流を流すのに適しているものである。また、これにより、生成された熱を除去することができるように、熱伝達流体がセル間を確実に流れるようにしてもよい。
【００２４】
燃料電池２は、インレットマニホールドとアウトレットマニホールドとの間に延在する電気化学セルの負極に水素を運ぶのに適した負極分配回路を含む。同様に、それは、インレットマニホールドとアウトレットマニホールドとの間に延びる電気化学セルの正極に酸素を運ぶのに適した正極分配回路を含む。したがって、負極及び正極のインレットマニホールドは、電気化学セルへの供給ガスの流体注入を確実にするのに適しており、アウトレットマニホールドは、電気化学セル及び燃料電池から未反応の供給ガス及び非反応種を除去するのに適している。負極及び正極の分配回路には、インレットマニホールドに連結された流体供給ラインＬａａ、Ｌａｃによって反応性ガスが供給され、流体排出ラインＬｅａ、Ｌｅｃはアウトレットマニホールドに連結される。
【００２５】
燃料電池は、電気化学セルのレベルに延在する冷却回路をさらに含む。冷却回路は、電気化学反応中に生成された熱を燃料電池から除去することを目的とする熱伝達流体を循環させるのに適している。冷却回路は、冷却回路内への／からの熱伝達流体の注入及び除去に適した流体冷却ラインＬｒに連結されている。
【００２６】
電気化学系１はまた、負極分配回路をパージするためのパージ装置３を含む。この目的のため負極排出ラインＬｅａ上に配置されたパージバルブ４を有しており、これにより負極分配回路に蓄積する可能性のある流体を短時間パージする。したがってパージバルブ４は制御開弁であり、すなわちその開動作と閉動作は制御信号Ｑａｏに基づいて制御される。この例では、燃料電池２はデッドエンドモードで動作するが、再循環モードで動作することもできる。
【００２７】
観察装置１０は、制御ベクトルＵ（ｔ）及び測定ベクトルＹ（ｔ）によって提供される情報に基づいて、電気化学系１の状態Ｘ（ｔ）のリアルタイム推定値
【数１７】
2018152326.tif 000019
を提供するのに適している。
【００２８】
この目的のために、観察装置１０は、電気化学系１の動作に関連する励起パラメーター及び／又は応答パラメーターを表す信号を提供するのに適したセンサー２０を含み、その値によって観察装置１０の制御Ｕ（ｔ）ベクトル及び測定Ｙ（ｔ）ベクトルが生成される。制御装置３にも電気的に接続されており、そこではパージバルブ４を制御することで制御信号Ｑａｏを受信する。制御信号Ｑａｏは、この実施例では、負極回路をパージする場合には１に等しい値をとり、パージバルブ４が閉じている場合にはゼロの値をとる。変形例として、信号Ｑａｏは、パージバルブの開度に関連して、０と１との間に含まれる値を取ることができる。
【００２９】
さらに、観察装置は、負極分配回路を流れるガスの負極圧力Ｐａのリアルタイムの測定値を伝えるのに適した圧力センサー２０１を含む。一例として、測定は、負極分配回路のインレットマニホールド又はアウトレットマニホールドで行ってもよい。また、燃料電池の電圧Ｖの測定値をリアルタイムで伝えるのに適した電圧センサー２０２と、電気負荷の要求に応答して燃料電池２によって供給される電流Ｉの測定値をリアルタイムで伝えるのに適した電流センサー２０３とを含む。また、観察装置は、電流Ｉと電圧Ｖの測定値に基づいて電解膜のオーム抵抗Ｒｍの値を算出する演算手段（図示せず）を備えている。オーム抵抗Ｒｍの値は電流Ｉと電圧Ｖとの測定と同様にさらに負極圧力Ｐａの測定値からもから導き出され、測定ベクトルＹ（ｔ）のパラメーターとなっている。
【００３０】
また、観察装置１０は温度センサー２０４を備えており、燃料電池の温度Ｔのリアルタイム測定を伝えるのに適しており、温度Ｔは例えば冷却回路を通って流れる熱伝達流体の温度である。例として、冷却回路の出口で温度Ｔを測定してもよい。正極分配回路を流れるガスの正極圧力Ｐｃのリアルタイム測定を伝えるのに適した圧力センサー２０５を備えている。一例として、測定は正極回路のインレットマニホールド又はアウトレットマニホールドで行うことができる。制御信号Ｑａｏ、電流Ｉの測定値、冷却温度Ｔの測定値、及び正極圧力Ｐｃの測定値は制御ベクトルＵ（ｔ）のパラメーターを構成している。制御ベクトル及び測定ベクトルのパラメーターの測定値は、時間とともに変化する。
【００３１】
観察装置は、制御ベクトルＵ（ｔ）及び測定ベクトルＹ（ｔ）に基づいて、電気化学系の状態Ｘ（ｔ）の推定値
【数１８】
2018152326.tif 000020
をリアルタイムで提供するのに適したコンピューター３０を備えている。
【００３２】
図２Ａは、本実施形態による観察方法のフローチャートを示す。この方法は、ここではスライディングモードのオブザーバークラスの状態オブザーバーとして行っている。当該観察方法は制御ベクトルＵ（ｔ）及び測定ベクトルＹ（ｔ）によって提供される情報に基づいて電気化学系１の状態Ｘ（ｔ）の推定値
【数１９】
2018152326.tif 000021
をリアルタイムで提供するのに適しており、状態表現と呼ばれる、系の動的モデルに基づいており、状態Ｘ（ｔ）の時間変化
【数２０】
2018152326.tif 000022
を表しており、スライディングモードにおける補正項が考慮されている。
【００３３】
より具体的には、観察方法における当該ステップは、連続する測定時間ｔｋにおいて実行される。時間は、時間サンプリング間隔Δｔに従って離散化される。当該方法の当該ステップは、測定時間の各インクリメントｋ＋１において繰り返される。本明細書の残りの部分を通して、記号Ｘ（ｔ）、Ｕ（ｔ）、Ｙ（ｔ）等が、考慮中の当該測定時間におけるベクトルＸ、Ｕ、Ｙ等の値であることを説明する。
【００３４】
電気化学系１の状態表示は、以下の式（１）の系によって定義される：
【数２１】
2018152326.tif 000023
（１）
ここで、ｆは状態関数であり、ｈは観察関数である。
【００３５】
電気化学系１のスライディングモードのオブザーバーは、以下の式（２）の系によって定義される。
【数２２】
2018152326.tif 000024
（２）
Ｋは利得パラメーターであり、Ｓは滑り面項と呼ばれる項であり、好ましくは推定値
【数２３】
2018152326.tif 000025
と測定ベクトルＹ：
【数２４】
2018152326.tif 000026
との間の差分として定義される。したがって、オブザーバーは、補正項と呼ばれる項、すなわちＫ．ｓｉｇｎ（Ｓ）を含み、推定値
【数２５】
2018152326.tif 000027
と測定ベクトルＹとの間の差分Ｓを最小にすることを可能にし、推定値
【数２６】
2018152326.tif 000028
は実際の状態Ｘに収束するものと推定される。
【００３６】
ステップ１００において、観察装置１０は、測定時間ｔｋにおける、パージ信号Ｑａｏの値、燃料電池の電流Ｉの値及び電圧Ｖの値、負極圧力Ｐａの値及び正極圧力Ｐｃの値、並びに熱伝達流体の温度Ｔの値を測定する。
【００３７】
ステップ２００において、パージ信号Ｑａｏ、測定電流Ｉ、測定温度Ｔ、及び測定された正極圧力Ｐｃを表すパラメーターに基づいて、測定時間ｔｋに対して、四次元の制御ベクトルＵが生成される。したがって、制御ベクトルＵは、次のように表され：
【数２７】
2018152326.tif 000029
ここで、Ｑａｏは、この例では、パージバルブが開いているか閉じているかに依拠する１又は０に等しい値である。
【００３８】
ステップ３００では、電解質膜のオーム抵抗Ｒｍ及び負極圧力Ｐａを表すパラメーターに基づいて、測定時間ｔｋについて、二次元の測定ベクトルＹが形成される。したがって、測定ベクトルＹは、以下のように表され：
【数２８】
2018152326.tif 000030
ここで、Ｒｍは、電流センサーによって測定された電流Ｉの値及び電圧センサーによって測定された電圧Ｖの値から推定される膜のオーム抵抗の値である。
【００３９】
このような制御Ｕベクトル及び測定Ｙベクトルの生成により、五次元の状態Ｘが測定時間ｔｋにおいて観察されるようになり、その中に負極での窒素の量、負極、正極及び膜での水、並びに負極での水素の量を表す状態変数が含まれる。状態ベクトルは次のように表され： 
【数２９】
2018152326.tif 000031
ここで、
【数３０】
2018152326.tif 000032
は、負極回路内の窒素のモル数であり；
【数３１】
2018152326.tif 000033
及び
【数３２】
2018152326.tif 000034
は、それぞれ負極回路及び正極回路における液相及び気相中の水のモル数であり；
【数３３】
2018152326.tif 000035
は、電解質膜中の水のモル数であり；
【数３４】
2018152326.tif 000036
は、負極回路内の水素のモル数である。
【００４０】
換言すれば、制御Ｕ（ｔ）、測定Ｙ（ｔ）、及び状態Ｘ（ｔ）のベクトルは、次のように記述される：
【数３５】
2018152326.tif 000037
（３）
【００４１】
ステップ４００では、状態Ｘの推定値
【数３６】
2018152326.tif 000038
の非補正時間変化
【数３７】
2018152326.tif 000039
と呼ばれる時間変化が算出されるが、測定時間ｔｋにおける、制御ベクトルＵ及び事前推定値
【数３８】
2018152326.tif 000040
に基づいており、上記の観察関数ｆに基づいており、ここではＲｏｂｉｎらによる発行文献Development and experimental validation of a PEM fuel cell 2D-model to study heterogeneities effects along large-area cell surface, Int. J. Hydrogen Energy, 40 (2015) 10211-10230に記載されているモデルに基づいて公式化されており、以下の式の系によって表され：
【数３９】
2018152326.tif 000041
（４）
以下のパラメーターは特にRobin 2015の発行文献で定義されている：
−  
【数４０】
2018152326.tif 000042
は膜を通る正極から負極への窒素透過流であり、特に状態変数
【数４１】
2018152326.tif 000043
及び
【数４２】
2018152326.tif 000044
に依存し、さらに制御変数Ｔ及びＰｃに依存するものである；
−  
【数４３】
2018152326.tif 000045
は、負極分配回路の出口への種Ｘのパージ流であり、特に状態変数
【数４４】
2018152326.tif 000046
、
【数４５】
2018152326.tif 000047
及び
【数４６】
2018152326.tif 000048
、並びに制御変数Ｑａｏに依存するものである；
−  
【数４７】
2018152326.tif 000049
及び
【数４８】
2018152326.tif 000050
は、それぞれ膜を通る負極から正極への水拡散の流れであり、これらは特に状態変数
【数４９】
2018152326.tif 000051
に依存し、さらに
【数５０】
2018152326.tif 000052
又は
【数５１】
2018152326.tif 000053
にそれぞれ依存し、また制御変数Ｔに依存するものである；
−  Ｆｅｏは膜を通る水の電気浸透流であり、特に状態変数
【数５２】
2018152326.tif 000054
に依存し、また制御変数Ｉに依存するものである；
−  Ｆｉｎ及びＦｏｕｔは、それぞれ正極分配回路に流入する水及び正極分配回路から流出する水の流れであり、特に制御変数Ｉ、Ｔ、Ｐｃに依存するものである；
−  
【数５３】
2018152326.tif 000055
は、正極分配回路に流入する水素の流れであり、特に、状態変数
【数５４】
2018152326.tif 000056
、
【数５５】
2018152326.tif 000057
及び
【数５６】
2018152326.tif 000058
に依存し、また制御変数Ｔに依存する。
−  Ｆはファラデー定数であり、ＥＷは膜の等価重量であり、ρｄは乾燥状態の電解膜の密度、ｅｍ及びＳは膜の厚さ及び面積である。
【００４２】
事前状態推定値
【数５７】
2018152326.tif 000059
は、第１の時間的反復ｔｋ＝１において、すなわち観察方法の初期化中において、状態ベクトルＸ（ｔｋ＝０）の所定の初期値に等しくてもよい。続く時間的反復ｔｋにおいて、事前状態推定値
【数５８】
2018152326.tif 000060
は、先行する測定時間ｔｋ−１におけるステップ７００において計算された値に対応する。
【００４３】
ステップ５００では、測定モードの時間ｔｋにおいて、スライディングモードにおける補正項εが計算され、ステップ３００で定義された測定ベクトルＹに基づいているとともに、前の測定時間で計算された同じ事前状態推定値
【数５９】
2018152326.tif 000061
に基づいており、ステップ４００と同様に先行する測定時間ｔｋ−１において計算される。これはスライディングモードでの定式を有しており、すなわち、それはＫ．ｓｉｇｎ（Ｓ）に等しく、ここでＫは所定の利得パラメーターであり、ｓｉｇｎ（．）は不連続符号関数であり、またＳは、測定ベクトルＹと、測定ベクトルＹの測定推定値
【数６０】
2018152326.tif 000062
と呼ばれる推定値との間の差分に基づいて計算されたスライディング面である。このステップは、図２Ｂを参照しつつ後述する。
【００４４】
ステップ６００では、測定時間ｔｋについて時間変化
【数６１】
2018152326.tif 000063
が計算され、非補正時間変動
【数６２】
2018152326.tif 000064
及び補正項εに基づいているが、ここでは、
【数６３】
2018152326.tif 000065
とεとの間の差分である。
【００４５】
ステップ７００では、測定時間ｔｋの状態ベクトルＸの推定値
【数６４】
2018152326.tif 000066
が計算されるが、ステップ６００で得られた時間的変動
【数６５】
2018152326.tif 000067
の時間的積分によるものである。
【００４６】
観察方法では、測定時間ｔｋを増分することによって、ステップ１００から７００の連続する時間的な反復が行われるので、更新された制御Ｕ（ｔ）及び測定Ｙ（ｔ）に基づいて状態推定値
【数６６】
2018152326.tif 000068
がリアルタイムで提供される。先行する時間的反復ｔｋにおけるステップ７００において計算された状態推定値
【数６７】
2018152326.tif 000069
の値は後の時間的反復ｔｋ＋１において、ステップ４００及び５００の「事前状態推定値」として考慮される。したがって観察方法には先行するステップ７００で計算された推定値
【数６８】
2018152326.tif 000070
を考慮する２つの閉じたループがあり、１つは非補正時間変動の計算に関連し、もう１つは補正時間に関連する。したがって、測定推定値
【数６９】
2018152326.tif 000071
と測定ベクトルＹ（ｔ）との間の差分Ｓ（ｔ）の最小化は、時間とともに徐々に行われる。この点に関して、後述するように、利得パラメーターＫは、多かれ少なかれ拘束のある有限時間で収束に達するように予め定められる。
【００４７】
したがって、事前に定義された状態Ｘ（ｔ）に対して、パージ信号及びセンサーからの測定値に基づいて制御ベクトルＵ（ｔ）及び測定Ｙ（ｔ）ベクトルを形成することにより、電気化学システム１の観察性を担保することができ、すなわち観察装置１０が、制御ベクトルＵ（ｔ）及び測定Ｙ（ｔ）ベクトルによって提供される情報に基づいて、任意の長さの時間［ｔｉ；ｔｆ］の状態Ｘ（ｔ）の推定値
【数７０】
2018152326.tif 000072
を実際に提供することができる。
【００４８】
さらに、測定ベクトルＹ（ｔ）の少なくとも１つの変数、すなわちここでは負極圧力Ｐａは、特に状態変数
【数７１】
2018152326.tif 000073
及び
【数７２】
2018152326.tif 000074
に依存し、負極回路内の窒素及び水素の量を表すことが示されている。したがって、観察装置により、状態変数
【数７３】
2018152326.tif 000075
及び
【数７４】
2018152326.tif 000076
を効率的に観察することができる。
これは、負極圧力Ｐａが特に状態変数
【数７５】
2018152326.tif 000077
に依存する限り、状態変数
【数７６】
2018152326.tif 000078
の場合にも当てはまる。「観察する」ことにより、推定された状態変数
【数７７】
2018152326.tif 000079
及び
【数７８】
2018152326.tif 000080
が実際に、実体の値と仮定される値に収束することが分かり、これは推定値
【数７９】
2018152326.tif 000081
と測定ベクトルＹとの間の差分を最小にするスライディングモード補正によるものである。
【００４９】
この観察方法はさらに、負極分配回路内及び正極分配回路内の水の量の測定値を伝えるために湿度センサーを使用する必要がないという利点がある。これはまた、燃料電池の圧力及び温度の時間的変化に対する制限を仮定することなく、状態Ｘの推定値
【数８０】
2018152326.tif 000082
をリアルタイムで提供することを可能にする。
【００５０】
図２Ｂは、図２Ａに示される観察方法のスライディングモードにおける補正項εを計算するステップ５００を詳細に示す。
【００５１】
ステップ５１０では、測定時間ｔｋについて、推定値
【数８１】
2018152326.tif 000083
が計算されるが、状態Ｘの事前推定値
【数８２】
2018152326.tif 000084
に基づいており、また以下の方程式の系によって表現される所定の観察関数ｈに基づている：
【数８３】
2018152326.tif 000085
（５）
【００５２】
上述したように、事前状態推定値
【数８４】
2018152326.tif 000086
は、一番目の時間反復ｔｋ＝１において、すなわち観察方法の開始中に、状態ベクトルＸ（ｔｋ＝０）の所定の初期値に対応していてもよい。後続の時間的反復では、事前状態推定値
【数８５】
2018152326.tif 000087
は、先行する測定時間ｔｋ−１におけるステップ７００で計算された値に対応する。
【００５３】
ステップ５２０において、推定値
【数８６】
2018152326.tif 000088
と測定値Ｙとの間の差分として定義されるスライディング面Ｓが測定時間ｔｋに対して計算され、次のようになる：
【数８７】
2018152326.tif 000089
（６）
【００５４】
ステップ５３０において、スライディング面Ｓの符号関数は、Ｓが正のときは１に、Ｓがゼロのときはゼロに、Ｓが負のときは−１に等しくなるように定義され、測定時間ｔｋに対して算出される。
【００５５】
ステップ５４０では、測定期間ｔｋに対して補正項ε＝Ｋ．ｓｉｇｎ（Ｓ）が計算され、次のようになる：
【数８８】
2018152326.tif 000090
（７）
ここで、εｉは補正項、σｉは利得パラメーターＫの項である。スライディング面の項ｓ２は項ε１、ε２及びε５のみに関連し、状態変数ｘ１、ｘ２及びｘ５も同様であるが、これらの変数ｓ２とε１、ε２及びε５とが負極に関係している限りにおいてである。
【００５６】
このようにして、スライディング面Ｓに適用される符号関数に対して利得パラメーターＫが乗算される。利得の項は、スライディング面Ｓ上での系の収束を保証するように予め定められている。これを達成するために、リアプノフ（Lyapunov）関数と呼ばれる二つの関数は次のように定義される：
【数８９】
2018152326.tif 000091
（８）
【００５７】
【数９０】
2018152326.tif 000092
を仮定することによって、すなわち、膜の水分量に関連するダイナミックレンジ（時間的変動）が膜の温度に関連するものよりも急速に変化すると仮定することによって、リアプノフ関数の時間変動
【数９１】
2018152326.tif 000093
及び
【数９２】
2018152326.tif 000094
、ここでは時間に関する項Ｖ１及びＶ２の一次導関数を、次のように表してもよい：
【数９３】
2018152326.tif 000095
（９）
【００５８】
系の収束は、時間的変動
【数９４】
2018152326.tif 000096
及び
【数９５】
2018152326.tif 000097
が負であると仮定するものであり、すなわち、
【数９６】
2018152326.tif 000098
及び
【数９７】
2018152326.tif 000099
であり、利得パラメーターＫの様々な項に関する以下の条件によって反映されるものである：
2018152326.tif 000100
（１０）
【００５９】
利得パラメーターＫの項σ１、σ２及びσ５は、その和が上記の条件を満たす限り個別に調整する必要はないことを表されている。項σ４のみは上記の基準を満たすことが好ましい。
【００６０】
変形例として、利得パラメーターＫの項の収束の条件はより限定的なものであってもよく、例えば、相対的な時間変動
【数９８】
2018152326.tif 000101
及び
【数９９】
2018152326.tif 000102
を定義することによって、厳密に負であるか、又はゼロではない負の値以下となる。
【００６１】
図３は、電気化学系１の観察装置１０の例示的な構造を示す。電気化学系１は、流体的であるとともに電気化学的である励起Ｅ（ｔ）を受け、応答Ｓ（ｔ）の形で動的に反応する。
【００６２】
観察装置１０はセンサー２０を備えており、これは各測定時間ｔｋごとに、電流Ｉ、電圧Ｖ、熱伝達流体の温度Ｔ、負極圧力Ｐａ及び正極圧力Ｐｃの測定値を伝えるのに適したものである。パージ装置３から放出される信号Ｑａｏを受け取るのにも適している。
【００６３】
コンピューター３０をさらに備えるが、プロセッサ、ランダムアクセスメモリー及び／又は読み出し専用メモリー、デジタルデータ伝送バスなどの通常のソフトウェア的手段で形成されており、状態関数ｆ及び観察関数ｈが実装されており、スライディングモード観察アルゴリズム（上述の式（２）のシステム）に沿っている。
【００６４】
コンピューターにはモジュールが含まれており、変動モジュール３１と呼ばれるものであり、状態表現の状態関数ｆが実装されている。この変動モジュール３１は各測定時間ｔｋごとの時間変動の計算に適しており、これは非補正時間変動
【数１００】
2018152326.tif 000103
というものであり、事前推定値
【数１０１】
2018152326.tif 000104
にかかるものであるが、制御ベクトルＵ及び事前推定値
【数１０２】
2018152326.tif 000105
を考慮しおり、例えば、先行する時間反復ｔｋ−１におけるステップ７００において計算された推定値である。
【００６５】
観察装置１０は、訂正モジュール３２と呼ばれるモジュールをさらに備え、これは測定時間ｔｋごとに、スライディングモードで表現された修正項ε（ｔ）を伝えることができるものであり、を測定ベクトルＹと測定推定値
【数１０３】
2018152326.tif 000106
との間の差分Ｓ（ｔ）を最小化するものである。補正項ε（ｔ）は、Ｋ．ｓｉｇｎ（Ｓ）型の定式化を有する限り、スライディングモード補正項と呼ばれ、Ｋは利得パラメーターであり、Ｓはスライディング面である。
【００６６】
補正モジュール３２は、推定値
【数１０４】
2018152326.tif 000107
を計算するためのサブモジュール３３を含み、観察関数ｈが実装される。これは
【数１０５】
2018152326.tif 000108
を計算するのに適しているが、計算は先行する推定値
【数１０６】
2018152326.tif 000109
に基づいており、これは例えば先行する時間反復ｔｋ−１におけるステップ７００で計算された推定値である。
【００６７】
比較器３４が備わっており、それは測定時間ｔｋごとにスライディング面Ｓ（ｔ）を、センサーから受信した測定ベクトルＹ（ｔ）とサブモジュール３３から受信した測定推定値
【数１０７】
2018152326.tif 000110
との間の差分として算出するのに適している。
【００６８】
これには、滑り面Ｓ（ｔ）に対して符号関数を適用するためのサブモジュール３５が含まれる。符号値（Ｓ（ｔ））は次に利得行列Ｋによって乗算されるが、これはサブモジュール３６によって最初に予め定められた項であり、滑り面Ｓ（ｔ）上の系の収束を確実なものとする。このようにして補正項ε（ｔ）が得られる。このようにして利得行列Ｋの値は時間とともに固定される。
【００６９】
コンピューター３０には比較器３７が備わっており、これにより時間変化
【数１０８】
2018152326.tif 000111
を各測定時間ｔｋごとに計算するが、これは補正されていない時間変化
【数１０９】
2018152326.tif 000112
と補正項ε（ｔ）との間の差分に当たる。その後、得られた時間変化
【数１１０】
2018152326.tif 000113
が積分器３８によって時間積分されることで測定時間ｔｋにおける電気化学系１の状態推定値
【数１１１】
2018152326.tif 000114
が生成される。
【００７０】
観測装置１０の収束は時間の経過と共に徐々に行われるが、２つの閉じたループを用いて行われ、１回の反復で計算された状態測定値
【数１１２】
2018152326.tif 000115
が補正モジュール３２の入力にも、変動モジュール３１の入力にも供給される限りは、次の時間反復のために行われる。
【００７１】
パージ装置３を観察装置１０に接続することが有利である。負極分配回路内の窒素の量
【数１１３】
2018152326.tif 000116
及び／又は水素の量
【数１１４】
2018152326.tif 000117
及び／又は水の量
【数１１５】
2018152326.tif 000118
を表す状態変数の推定値に応じて、推定値が所定の最大値に達したときは、推定値が最大値を下回るようになるまで、あるいは推定値が所定の最小値に実質的に等しくなるまで、パージ装置３がパージバルブを開いてもよい。これにより、観察装置１０による状態推定値
【数１１６】
2018152326.tif 000119
の伝達がリアルタイムに近づくにつれてますます可能になる。したがって、負極分配回路のパージを制御して、窒素の量、水素の量及び／又は水の量が燃料電池の電気化学的性能に対して悪影響を与えるものとなることを回避することが可能である。
【００７２】
図４Ａから図４Ｃは、負極ＨＲａ（図４Ａ）及び正極ＨＲｃ（図４Ｂ）における相対湿度の時間変化並びに負極での窒素飽和度の時間変化（図４Ｃ）を示しており、デッドエンドモードで動作するＰＥＭ燃料電池の負極分配回路のパージを制御する場合に関する。パージバルブは、負極での窒素の量が所定の最大値に達したときに開かれる。相対湿度及び窒素飽和パラメーターは、上記の観察方法及び装置によって得られる推定値
【数１１７】
2018152326.tif 000120
の状態値から導き出される。
【００７３】
本実施形態に係る観察方法及び装置（実線）により得られたこれらの時間的な変化を、上述のRobin 2015の刊行物に示されている物理的な基準モデル（点線）と比較する。
【００７４】
燃料電池によって供給される電流Ｉは、電気負荷によって引き出される電力から推測され、その電力は、乗用車等の国際調和排出ガス・燃費試験法（ＷＬＴＰ）によって定義される。相対湿度は、負極分配回路の入口で０％、正極分配回路の入口で５０％に固定される。酸素の化学量論比は、正極で２に固定されている。水素のそれは、燃料電池のデッドエンド動作に従って、負極において１に固定される。伝熱流体の温度は、標準的なガウスの法則に従っていてランダムであり、平均温度は３５３Ｋであり、バラつきが１０２Ｋであり、１０秒ごとにサンプリングされ、変動が毎秒毎にプラス／マイナス７Ｋに限られている。正極圧力Ｐｃ及び負極圧力Ｐａはランダムであるが、正規のガウス法（noＲｍal Gaussian law）にしたがっており、２．１０５Ｐａの平均圧力を有し、バラつきが（４．１０４）２Ｐａであり、毎秒サンプリングされ、変動はプラス／マイナス１０４Ｐａ／ｓに限られている。
【００７５】
この例では、利得パラメーターＫの項σ４が１１０に固定されることで上記決定された条件が満たされており、また項σ１、σ２及びσ５が、それぞれ０．０６、０．００及び０．０６に固定されることで、それらの和は上記で決定された条件が満たされている。項σ３は０に固定されている。
【００７６】
図４Ｃに示すように、パージバルブが閉鎖されているとき（Ｑａｏ＝０）、負極での窒素飽和は、膜を通る拡散により、所定の最大値に達するまで、ここでは０．０５％程度まで増加する傾向がある。最大値に達すると、パージ信号Ｑａｏは１に等しくなり、パージバルブは開かれる。次いで、パージバルブが閉じられるまで、負極での窒素の飽和値が低下する。このようにして、信号は時間とともに継続する。したがって、この変数は、観察装置によって得られる推定値
【数１１８】
2018152326.tif 000121
から推測され、Robin 2015のモデルによって提供される同じ変数の変化に従うと考えられる。負極及び正極における相対湿度の変化は、観察装置（図４Ａ及び４Ｂ）によって推定される状態変数
【数１１９】
2018152326.tif 000122
及び
【数１２０】
2018152326.tif 000123
から推測され、Robin 2015の参照モデルによって得られたものとよく一致している。
【００７７】
特定の実施形態について説明した。様々な改変及び変形が当業者には明らかであろう。
【００７８】
上述の観察方法及び観察装置は、デッドエンド動作モード及び燃料電池の再循環動作モードの両方に等しく適用することができる。この場合、負極アウトレットマニホールドは、再循環ループによって負極インレットマニホールドに連結される。パージバルブを開く命令は、電気化学系の状態Ｘの変数の少なくとも１つに従って発行される。
【００７９】
状態関数ｆ及び観察関数ｈは、ここでは特に、上述のRobin 2015の論文に記載されている電気化学系の動的物理モデルに基づく。それらは、当業者の知識に従って、観察中の燃料電池の特性に従って、又は他の電気化学的モデルに基づいて改変されてもよい。
【図１】
2018152326.tif 000153
【図２Ａ】
2018152326.tif 000154
【図２Ｂ】
2018152326.tif 000155
【図３】
2018152326.tif 000156
【図４Ａ】
2018152326.tif 000157
【図４Ｂ】
2018152326.tif 000158
【図４Ｃ】
2018152326.tif 000159
【外国語明細書】
2018152326000001.pdf