(57)【要約】  本発明では、２つ以上の経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）を備えた燃焼プロセス（ＣＰ）から生じるガスの量（Ｇ）の推定のための経験的アンサンブルに基づく仮想センサシステム（ＶＳ）を提供する。ガスの量（Ｇ）は各経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の中で推定され、結合関数（ｆ）は、個々の経験的モデル（ｙ１、ｙ２、．．．、ｙｍ）の１つずつからのガスの推定値より正確であるガスの量（Ｇ）の結合推定値を提供するため、経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）からの結果を結合する。仮想センサシステム（ＶＳ）の全性能は、経験的モデル（ｙ１、ｙ２、．．．、ｙｍ）の個数を増加することにより高められる。
2010537192.tif 000002
【特許請求の範囲】
【請求項１】
経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の１つずつが、前記プロセス（ＣＰ）からの経験的データ（ＥＤ）を使用して訓練されるように用意され、前記プロセス（ＣＰ）の１つ以上のセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）から１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を受け取り、前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて、前記ガスの量（Ｇ）を表現する信号出力値（ｙ１、ｙ２、．．．、ｙｍ）を計算するためさらに用意されている、２つ以上の経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）と、
前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）を受け取り、前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の関数として、前記ガスの量（Ｇ）を表現する仮想センサ出力値（ｙＲ）を連続的に計算するため用意されている結合関数（ｆ）と
を備える、燃焼プロセス（ＣＰ）から生じるガスの量（Ｇ）の推定のためのアンサンブルに基づく仮想センサシステム（ＶＳ）。
【請求項２】
前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）のすべてが同一の構造を有する請求項１に記載の仮想センサシステム（ＶＳ）。
【請求項３】
前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）のすべてが同じ信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）のセットを受け取るため用意されている請求項１に記載の仮想センサシステム（ＶＳ）。
【請求項４】
前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）がニューラルネットワークである請求項１に記載の仮想センサシステム（ＶＳ）。
【請求項５】
前記結合関数（ｆ）が前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の平均値として前記仮想センサ出力値（ｙＲ）を連続的に計算するため用意されている請求項１に記載の仮想センサシステム（ＶＳ）。
【請求項６】
前記結合関数（ｆ）が、１つ以上の前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を受け取り、前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）が前記１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて動的に重み付けされ、仮想センサ出力値（ｙＲ）を計算するため用意されている請求項１に記載の仮想センサシステム（ＶＳ）。
【請求項７】
前記結合関数（ｆ）が、１つ以上の前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を受け取り、前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）と前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）と前記経験的モデル（ＮＮＲ）の構造とに基づいて仮想センサ出力値（ｙＲ）を計算するため用意されている経験的モデル（ＮＮＲ）である請求項１に記載の仮想センサシステム（ＶＳ）。
【請求項８】
前記センサが前記仮想センサ出力値（ｙＲ）の所定の性能要件を実現するため、ある程度の個数の前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）をインスタンス化できるように用意されている請求項１に記載の仮想センサシステム（ＶＲ）。
【請求項９】
１つ以上の前記センサ（Ｓ１、Ｓ２、．．．、Ｓｍ）がガスの量（Ｇ）の推定のためのアンサンブルに基づく仮想センサシステム（ＶＳ）であり、連結できるように用意されている請求項１に記載の仮想センサシステム（ＶＲ）。
【請求項１０】
前記プロセス（ＣＰ）からの経験的データを使って経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）のアンサンブルを訓練するステップと、
前記訓練された経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）に前記プロセス（ＣＰ）の１つ以上のセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）からの前記１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を供給するステップと、
前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）において、前記ガスの量（Ｇ）を表現する信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の計算を実行するステップと、
前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）を連続的に合成し、前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の関数として、前記ガスの量（Ｇ）を表現する仮想センサ出力値（ｙＲ）を計算するステップと
を含む、１つ以上のセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）からの１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）から、燃焼プロセス（ＣＰ）から生じるガスの量（Ｇ）を推定する方法。
【請求項１１】
前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）のすべてが同一の構造を有する請求項１０に記載の方法。
【請求項１２】
前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）のすべてに同じ信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）のセットを供給するステップを含む請求項１０に記載の方法。
【請求項１３】
前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）がニューラルネットワークである請求項１０に記載の方法。
【請求項１４】
前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の平均値として、ガスの量（Ｇ）を表現する前記仮想センサ出力値（ｙＲ）を連続的に計算するステップを含む請求項１０に記載の方法。
【請求項１５】
１つ以上の前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を連続的に受け取り、前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）が前記１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて動的に重み付けされ、仮想センサ出力値（ｙＲ）を計算するステップを含む請求項１０に記載の方法。
【請求項１６】
１つ以上の前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を受け取り、前記信号出力値（ｙ１、ｙ２、．．．、ｙｍ）と前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）と前記経験的モデル（ＮＮＲ）の構造とに基づいて仮想センサ出力値（ｙＲ）を計算するステップを含む請求項１０に記載の方法。
【請求項１７】
前記仮想センサ出力値（ｙＲ）の所定の性能要件に基づいて要求された個数の前記経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）を計算するステップを含む請求項１０に記載の方法。
【請求項１８】
１つ以上の前記信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）がそれ自体で請求項１１に記載された方法による仮想センサ出力値（ｙＲ）であることによって反復的である請求項１０に記載された方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、経験的アンサンブルに基づく仮想センシングのための方法及びシステムに関する。特に、燃焼プロセスからのＮＯｘ、ＣＯ２などの排出量を測定する仮想ガスセンサのための方法及びシステムに関する。
【背景技術】
【０００２】
ＮＯｘは、燃焼中に生成される一窒素酸化物（ＮＯ及びＮＯ２）の総称である。ＮＯｘは、燃焼空気にある二原子窒素の高温酸化を介して形成される。さらに、ある種の石炭及び石油のような窒素添加燃料の燃焼は、燃料に結合した窒素をＮＯｘに変換する。
【０００３】
大気ＮＯｘは、最終的に酸性雨の一因となる硝酸を形成する。１９９７年に５４カ国によって批准された京都議定書では、ＮＯ２を温室ガスとして分類し、長距離越境大気汚染条約のいわゆるＧｏｔｈｅｎｂｕｒｇ議定書が行っているように、ＮＯ２排出の全世界的な削減を求めている。
【０００４】
その結果、ＮＯｘ排出は数カ国において規制され、例えば、スウェーデンでは、１９９２年以降、燃焼プラントからのＮＯｘ排出に対して課徴金が課された。また、ノルウェーでは、２００７年以降、すべてのＮＯｘ排出に対して一般経費が課された。同様に、フランス及びイタリアでも課徴金が課されている。さらに、米国では、ＮＯｘ排出量許可証の制度がある。
【０００５】
従って、特定のプラント又は燃焼プロセスから排出されるＮＯｘの量を測定することがある。
しかしながら、優れたセンサを開発することは、例えば、高温及び煤煙を伴う過酷な動作環境のために困難である。必要とされる感度は高く、典型的に、ＮＯのレベルは約１００〜２０００ｐｐｍであり、ＮＯ２の濃度は２０〜２００ｐｐｍである。また、ガス流からの冷却のような様々な誤差が発生する。
【０００６】
類似した考慮に基づいて、炭素酸化物及び硫黄酸化物のようなその他のガスを測定する必要がある。
【０００７】
一般に、利用可能な計器が測定のため適切でない様々な状況があり、以下のリストは、最も一般的な状況を指定する。（ＢｉｏＣｏｍｐ  Ｓｙｓｔｅｍｓ社によって、彼らのウェブページ、ｈｔｔｐ：／／ｗｗｗ．ｂｉｏｃｏｍｐｓｙｓｔｅｍｓ．ｃｏｍ／ｔｅｃｈｎｏｌｏｇｙ／ｖｉｒｔｕａｌｓｅｎｓｏｒｓ／ｉｎｄｅｘ．ｈｔｍ（２００８年７月２５日）に当初提案された）
【０００８】
１．着目している物理量はオンラインで測定されない。典型的な事例は、サンプルが分析のため定期的に研究室に送られる場合である。サンプルは、環境排出、生産量又はプロセス条件を制御するために分析される大気サンプル、水サンプル、石油サンプル又は材料サンプルである。
【０００９】
２．利用可能な物理センサは、特に、自動制御で用いるには過度に低速である。
【００１０】
３．物理センサは極端に下流であり、例えば、最終製品が製造ばらつきを検出するため連続的に監視される。しかし、是正の判断を行うには、情報の送信が遅すぎる。
【００１１】
４．物理センサは過度に高価である。
【００１２】
５．物理センサを設置する手段がなく、例えば、物理的空間がない。
【００１３】
６．センサ環境が過度に不良である。
【００１４】
７．物理センサが不正確である。利用可能な物理センサは内在する不正確さ又は劣化のいずれかの影響を受ける場合がある。ベンチュリ流量計におけるスケーリングは典型的な例である。
【００１５】
８．物理センサは維持に費用がかかる。
【００１６】
ソフトセンシング又はプロキシセンシングとしても知られている仮想センシング技術は、費用がかかる、又は非実用的である物理測定装置及びセンサシステムに対する実現可能であり、かつ、経済的な代替案を提供するため使用されるソフトウェアベースの技術である。仮想センシングシステムは、着目している量の推定値を計算するため、他のオンライン測定量及びプロセスパラメータから入手できる情報を使用する。
【００１７】
多種多様な仮想センシング技術が利用可能であり、分析技術と経験技術との２つの主要なカテゴリーに分類されることができる。
【００１８】
分析技術は、測定推定値の計算の基礎を、着目している量とその他の利用可能な測定量及びパラメータとの関係を規定する物理法則の近似による。
【００１９】
「第１原理」モデルに基づく分析的技術を使用する重大な利点は、物理的に測定できない量が関係する物理モデル方程式から導出可能であるとき、物理的に測定できない量の計算を可能にすることである。
【００２０】
分析的アプローチの主な欠点は、分析的アプローチが効果的であるために正確な定量的数学モデルを必要とすることである。大規模システムに対し、このような情報は入手できない場合があり、又は、このような情報は編集するため過度に費用がかかり、かつ、時間がかかる場合がある。同様に、プラント又はプロセスに変更が加えられる場合、エンジニアリング作業が物理モデルを更新し修正するために必要になる。モデリングツールがこのようなモデル構築及び維持管理を支援するため利用できるが、プロセス専門家がモデルを最新の状態に保つために必要とされる。
【００２１】
経験的技術は、測定推定値の計算の基礎を、同量の利用可能な履歴測定データと、この履歴測定データとその他の利用可能な測定量及びパラメータとの相関とに置く。非測定量の履歴データは、一時的に設置されたセンサシステムを用いる実際の測定キャンペーン、研究室分析の記録、又はオンラインで実行するには計算のために過度に費用がかかる複雑な分析モデルを用いる詳細な推定のいずれかから導出されることができる。詳細な推定は、履歴データを利用できないことが明かである測定できない量を推定するために経験的仮想センサを開発することを望む場合、唯一の可能な選択肢である。
【００２２】
経験的仮想センシングは、線形回帰（Ｎ．Ｒ．Ｄｒａｐｅｒ  ａｎｄ  Ｈ．Ｓｍｉｔｈ，１９９８．Ａｐｐｌｉｅｄ  Ｒｅｇｒｅｓｓｉｏｎ  Ａｎａｌｙｓｉｓ，Ｗｉｌｅｙ  Ｓｅｒｉｅｓ  ｉｎ  Ｐｒｏｂａｂｉｌｉｔｙ  ａｎｄ  Ｓｔａｔｉｓｔｉｃｓを参照）、重み付き最小二乗回帰（Ａ．Ｂｊｏｒｃｋ，１９９６．Ｎｕｍｅｒｉｃａｌ  Ｍｅｔｈｏｄｓ  ｆｏｒ  Ｌｅａｓｔ  Ｓｑｕａｒｅ  Ｐｒｏｂｌｅｍｓ，Ｃａｍｂｒｉｄｇｅ．を参照）、カーネル回帰（Ｊ．Ｓ．Ｓｉｍｏｎｏｆｆ，１９９６．Ｓｍｏｏｔｈｉｎｇ  Ｍｅｔｈｏｄｓ  ｉｎ  Ｓｔａｔｉｓｔｉｃｓ，Ｓｐｒｉｎｇｅｒ．を参照）、回帰ツリー（Ｌ．Ｂｒｅｉｍａｎ，Ｊ．Ｆｒｉｅｄｍａｎ，Ｒ．Ａ．Ｏｌｓｈｅｎ  ａｎｄ  Ｃ．Ｊ．Ｓｔｏｎｅ，１９８４．Ｃｌａｓｓｉｆｉｃａｔｉｏｎ  ａｎｄ  ｒｅｇｒｅｓｓｉｏｎ  ｔｒｅｅｓ．Ｗａｄｓｗｏｒｔｈ．を参照）、サポートベクトル回帰（Ｈ．Ｄｒｕｃｋｅｒ，Ｃ．Ｊ．Ｃ．Ｂｕｒｇｅｓ，Ｌ．Ｋａｕｆｍａｎ，Ａ．Ｓｍｏｌａ  ａｎｄ  Ｖ．Ｖａｐｎｉｋ，１９９７．Ｓｕｐｐｏｒｔ  Ｖｅｃｔｏｒ  Ｒｅｇｒｅｓｓｉｏｎ  Ｍａｃｈｉｎｅｓ．Ａｄｖａｎｃｅ  ｉｎ  Ｎｅｕｒａｌ  Ｉｎｆｏｒｍａｔｉｏｎ  Ｐｒｏｃｅｓｓｉｎｇ  Ｓｙｓｔｅｍ  ９，ＮＩＰＳ  １９９６，１５５−１６１，ＭＩＴ  Ｐｒｅｓｓ．を参照）、ニューラルネットワーク回帰（Ｊ．Ｈｅｒｔｚ，Ａ．Ｋｒｏｇｈ，ａｎｄ  Ｒ．Ｐａｌｍｅｒ，１９９１．Ｉｎｔｒｏｄｕｃｔｉｏｎ  ｔｏ  ｔｈｅ  Ｔｈｅｏｒｙ  ｏｆ  Ｎｅｕｒａｌ  Ｃｏｍｐｕｔａｔｉｏｎ，Ａｄｄｉｓｏｎ−Ｗｅｓｌｅｙ：Ｒｅｄｗｏｏｄ  Ｃｉｔｙ，Ｃａｌｉｆｏｒｎｉａ．を参照）などの多種多様の統計モデリング法又は機械学習モデリング法を使用して実施されることができる関数近似及び関数回帰技術に基づく。
【００２３】
データ駆動型モデリングとしても知られている経験的モデリングは、条件を分析し、オペレーショナルデータからプロセスの進展を予測するため使用される技術の集合を対象とする。経験的モデリングは、現実の実際的な事例では、どちらも多くの場合に欠けているが、プロセスの詳細な物理的理解も、プラントとプラントのコンポーネントの材料特性、幾何学的構成、及びその他の特徴の知識も必要としないという利点がある。
【００２４】
基礎となるプロセスモデルは、多くの場合に「学習」と呼ばれる手続を通じて、測定又はシミュレーションされたプラントデータを一般的な線形又は非線形モデルに当てはめることにより特定される。この学習プロセスは、能動的でも受動的でもよく、プロセス変数間の関係の特定とモデルへの組み込みを伴う。能動的学習プロセスは、勾配に基づくパラメータ調整を通じて誤差関数を最小化する反復プロセスを伴う。受動的学習プロセスは、数学的反復を必要とせず、代表的なデータベクトルを訓練マトリクスに編集することだけで構成される。
【００２５】
経験的モデルを設計する際の重要な検討事項は、訓練データが正確な予測が問い合わせされる条件の例を提供しなければならないことである。これは、すべての可能な条件が訓練データの中に存在しなければならないということではなく、訓練データはこれらの条件を適切に対象とすべきであるということである。経験的モデルは、内挿予測を提供するが、訓練データは、この予測が十分に正確であるため、内挿場所の上下を適切に対象としなければならない。正確な外挿、すなわち、訓練データの外側にあるデータの推定値を提供することは、殆どの経験的モデルに対し、可能性がなく、また信頼性もない。
【００２６】
経験的モデルは、モデルを開発するため使用されたデータが収集された運転条件と同一又は類似の運転条件に適用されたときに限り確実に正確である。プラント条件又は運転が著しく変化するとき、モデルは学習された空間の外側で外挿するように強制され、その結果は低信頼性である。この観察は、非線形モデルが、既知の線形形式で外挿する線形モデルと異なり、未知の形式で外挿するので、非線形経験的モデルに対して、特に、真実である。人工ニューラルネットワーク及び局所多項式回帰モデルはどちらも非線形であり、一方、主成分分析及び部分最小二乗のような変換に基づく技術は線形技術である。外挿は、線形モデルを使用するとしても、測定されたプロセス変数の間の純粋な線形関係の存在が期待されないので、経験的モデルのため推奨されない。さらに、プロセスに対する線形近似は、これらの極端な領域における訓練データの密度が非常に低いか、又は非存在であるため外挿中に妥当性が低い。
【００２７】
人工ニューラルネットワークモデル（Ｊ．Ｈｅｒｔｚ，Ａ．Ｋｒｏｇｈ，ａｎｄ  Ｒ．Ｐａｌｍｅｒ，１９９１．Ｉｎｔｒｏｄｕｃｔｉｏｎ  ｔｏ  Ｔｈｅ  Ｔｈｅｏｒｙ  ｏｆ  Ｎｅｕｒａｌ  Ｃｏｍｐｕｔａｔｉｏｎ，Ａｄｄｉｓｏｎ−Ｗｅｓｌｅｙ：Ｒｅｄｗｏｏｄ  Ｃｉｔｙ，Ｃａｌｉｆｏｒｎｉａを参照）は、非線形加算装置として動作する単純な計算ノードの層を含む。これらのノードは、重み付き結合線と高度に相互結合され、これらの重みは、訓練データが訓練プロセス中にニューラルネットワークへ与えられたときに調整される。うまく訓練されたニューラルネットワークは種々のタスクを実行可能であり、その中で最も一般的なタスクは、出力値の予測、分類、関数近似、及びパターン認識である。
【００２８】
関連付けられた結合重みのセットを有するニューラルネットワークの層だけが正当な処理層として認識される。ニューラルネットワークの入力層は、関連付けられた重みのセットを有するので、真の処理層ではない。それに対して、出力層は関連付けられた重みのセットを有する。従って、ニューラルネットワークの中で層の個数を記述する最も効率的な用語法は、中間層という用語の使用による。中間層は、出力層を除いた正当な層である。
【００２９】
ニューラルネットワーク構造は、ある程度の個数の中間層と、出力層とにより構成される。ニューラルネットワークの計算能力は、単一の非線形中間層を備えるニューラルネットワークは、十分な個数の中間ノードが与えられる場合、任意の非線形関数を近似できるという一般的な関数近似理論によって証明された。
【００３０】
ニューラルネットワーク訓練プロセスは、ニューラルネットワークの重みを小さい乱数に初期化することから始まる。ネットワークは、その後に、入力ベクトルと、ターゲットと呼ばれることがよくある対応する所望の出力とのセットからなる訓練データが与えられる。ニューラルネットワーク訓練プロセスは、特定の入力ベクトル／ターゲットのペアが与えられた場合に、ネットワークの出力を所望の値に近づけるための内部重みの反復調整である。重みはネットワークが所望の出力を計算する可能性を高めるように調整される。訓練プロセスは、ネットワークの出力値と所望の出力値との間の平均二乗誤差（ＭＳＥ）を最小化しようとする。ＭＳＥ関数の最小化は圧倒的に最も一般的なアプローチであるが、その他の誤差関数が利用可能である。
【００３１】
ニューラルネットワークは、工業設備からのプロセスデータを監視するパターン認識問題に適用可能である強力なツールである。ニューラルネットワークは、非線形システムを監視し、複雑なデータセットの中の欠陥パターンを認識するため適している。反復訓練プロセスであるため、ニューラルネットワークモデルを作り上げるために必要とされる計算量は、その他のタイプの経験的モデルの計算量より多い。従って、計算要件は、典型的にその他の経験的モデルタイプより制限的であるモデルサイズの上限の原因となる。
【００３２】
コミッティモデリングとしても知られているアンサンブルモデリング（Ｔ．Ｇ．Ｄｉｅｔｔｅｒｉｃｈ（Ｅｄ．），２０００．Ｅｎｓｅｍｂｌｅ  Ｍｅｔｈｏｄｓ  ｉｎ  Ｍａｃｈｉｎｅ  Ｌｅａｒｎｉｎｇ，Ｌｅｃｔｕｒｅ  Ｎｏｔｅｓ  ｉｎ  Ｃｏｍｐｕｔｅｒ  Ｓｃｉｅｎｃｅ；Ｖｏｌ．１８５７．Ｓｐｒｉｎｇｅｒ−Ｖｅｒｌａｇ，Ｌｏｎｄｏｎ，ＵＫを参照）は、単一の予測モデルを構築する代わりに、コンポーネントモデルのセットが作り上げられ、それらの独立した予測が単一の集計された予測を生成するため結合される技術である。結果として得られるコンパウンドモデル（アンサンブルと呼ばれる）は、単一のコンポーネントモデルより一般的に正確であり、過剰適合現象に対してよりロバストであるという傾向があり、非常に低減された変動を有し、準最適モデル訓練手続と関連付けられることがある不安定性問題を回避する。
【００３３】
アンサンブルでは、各モデルは、一般的に別々に訓練され、各コンポーネントモデルの予測出力がその後にアンサンブルの出力を生成するために結合される。しかしながら、いくつかのモデルの出力の結合は、それらのモデルの予測の間に何らかの形式の「不一致」が存在する場合に限り役に立つ（Ｍ．Ｐ．Ｐｅｒｒｏｎｅ  ａｎｄ  Ｌ．Ｎ．Ｃｏｏｐｅｒ，１９９２．Ｗｈｅｎ  ｎｅｔｗｏｒｋｓ  ｄｉｓａｇｒｅｅ：ｅｎｓｅｍｂｌｅ  ｍｅｔｈｏｄｓ  ｆｏｒ  ｈｙｂｒｉｄ  ｎｅｕｒａｌ  ｎｅｔｗｏｒｋｓ，Ｎａｔｉｏｎａｌ  Ｓｃｉｅｎｃｅ  Ｆｕｎｄａｔｉｏｎ，ＵＳＡを参照）。当然、同一のモデルの結合は性能改善を生じない。一般的に採用される一つの方法は、訓練中に各モデルが調べる訓練セットを変更することにより、モデル間の不一致を生成することを試みる、いわゆるバギング法である（Ｌ．Ｂｒｅｉｍａｎ，１９９６．Ｂａｇｇｉｎｇ  Ｐｒｅｄｉｃｔｏｒｓ，Ｍａｃｈｉｎｅ  Ｌｅａｒｎｉｎｇ，２４（２），ｐｐ．１２３−１４０を参照）。バギングは、訓練セットのランダムサンプリングに関して各モデルを訓練することにより、アンサンブルのための個体を作成し、最終的な予測を形成する際に、各コンポーネントモデルに等しい重みを与えるアンサンブル法である。アンサンブル生成及びコンポーネントモデル集計のためのその他のより精緻なスキームが存在し、新しいスキームが考え出される。
【００３４】
全体的なモデル変動を低減するためのアンサンブルの使用は、悪条件問題を回避し、過剰なモデル変動の類似した制御を実現するためにニューラルネットワークモデルの訓練及びアーキテクチャを制約する正則化法と密接な関係がある（Ａ．Ｖ．Ｇｒｉｂｏｋ，Ｊ．Ｗ．Ｈｉｎｅｓ，Ａ．Ｕｒｍａｎｏｖ，ａｎｄ  Ｒ．Ｅ．Ｕｈｒｉｇ．２００２．Ｈｅｕｒｉｓｔｉｃ，Ｓｙｓｔｅｍａｔｉｃ，ａｎｄ  Ｉｎｆｏｒｍａｔｉｏｎａｌ  Ｒｅｇｕｌａｒｉｚａｔｉｏｎ  ｆｏｒ  Ｐｒｏｃｅｓｓ  Ｍｏｎｉｔｏｒｉｎｇ．Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｊｏｕｒｎａｌ  ｏｆ  Ｉｎｔｅｌｌｉｇｅｎｔ  Ｓｙｓｔｅｍｓ，１７（８），ｐｐ  ７２３−７５０，Ｗｉｌｅｙを参照）。
【００３５】
米国特許第５３８６３７３号明細書の「Ｖｉｒｔｕａｌ  ｃｏｎｔｉｎｕｏｕｓ  ｅｍｉｓｓｉｏｎ  ｍｏｎｉｔｏｒｉｎｇ  ｓｙｓｔｅｍ  ｗｉｔｈ  ｓｅｎｓｏｒ  ｖａｌｉｄａｔｉｏｎ」は、プラントの運転を制御するために、ニューラルネットワークに基づく排出のための仮想センサの使用を開示する。
【００３６】
米国特許第６８８２９２９号明細書の「ＮＯｘ  ｅｍｉｓｓｉｏｎ−ｃｏｎｔｒｏｌ  ｓｙｓｔｅｍ  ｕｓｉｎｇ  ａ  ｖｉｒｔｕａｌ  ｓｅｎｓｏｒ」は、エンジンの運転を制御するために、ニューラルネットワークに基づく、排出のための仮想センサの使用を開示する。
【００３７】
米国特許第７２８０９８７号明細書の「Ｇｅｎｅｔｉｃ  ａｌｇｏｒｉｔｈｍ  ｂａｓｅｄ  ｓｅｌｅｃｔｉｏｎ  ｏｆ  ｎｅｕｒａｌ  ｎｅｔｗｏｒｋ  ｅｎｓｅｍｂｌｅ  ｆｏｒ  ｐｒｏｃｅｓｓｉｎｇ  ｗｅｌｌ  ｌｏｇｇｉｎｇ  ｄａｔａ」は、望ましい当てはめ関数値を有するアンサンブルを選択するために多目的当てはめ関数を用いるアルゴリズムを使用して、地球物理データを処理するニューラルネットワークアンサンブルを生成する方法を開示する。
【００３８】
仮想センシングは、ＮＯｘ及びその他のガスを測定するための魅力的な解決法であるが、上記の参照されたシステムより簡単に実施でき、高精度であり、ロバスト性があり、安定性がある仮想センシングのためのシステムが必要である。
【発明の概要】
【発明が解決しようとする課題】
【００３９】
本発明は、経験的モデリングのアンサンブルモデリングとの組み合わせによってガス検知に適した仮想センサの精度、ロバスト性、安定性及び簡易性の問題を解決する。
【課題を解決するための手段】
【００４０】
実施形態では、本発明は、燃焼プロセスから生じるガスの量の推定のためのアンサンブルに基づく仮想センサシステムであって、各経験的モデルがプロセスからの経験的データを使用して訓練されるように用意され、プロセスの１つ以上のセンサから１つ以上の信号入力値を受け取り、信号入力値に基づいて、ガスの量を表現する信号出力値を計算するためさらに用意されている、２つ以上の経験的モデルと、信号出力値を受け取り、信号出力値の関数として、ガスの量を表現する仮想センサ出力値を連続的に計算するため用意されている結合関数とを備える仮想センサシステムである。
【００４１】
実施形態では、本発明は、１つ以上のセンサからの１つ以上の信号入力値から、燃焼プロセスから生じるガスの量を推定する方法であって、経験的モデルのアンサンブルをプロセスからの経験的データを用いて訓練するステップと、プロセスの１つ以上のセンサからの１つ以上の信号入力値を訓練された経験的モデルに供給するステップと、信号入力値に基づいて経験的モデルの中の信号出力値の計算を実行するステップと、信号出力値を連続的に結合し、信号出力値の関数として、ガスの量を表現する仮想センサ出力値を計算するステップとを備える方法である。
【００４２】
本発明の実施形態では、結合関数（ｆ）は、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の平均値として仮想センサ出力値（ｙＲ）を連続的に計算するため用意されている。平均値は、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の幾何平均もしくは算術平均値として、又は、中央値として計算されることができる。
【００４３】
平均計算は、実施が容易であることに加えて、単一ノードの仮想センサを使うと不可能の場合のある所要の精度を達成することも可能にすることが示される。
【００４４】
本発明の実施形態では、経験的モデル又は内部ノードのすべてが同一構造でもよい。この設定は、所要の個数の内部ノードがテンプレートノードに基づいて仮想センサシステムの中で簡単にインスタンス化できるという利点がある。さらに、ノードはすべてが燃焼プロセスのセンサから同じ信号入力値のセットを受け取るため用意されてもよい。センサからの信号はすべてのノードに配布され、特別な場合を取り扱う余計な作業が回避される。
【００４５】
実施形態では、本発明による仮想センサシステムの精度は、より多数の経験的モデルをインスタンス化することによって高められてもよい。従って、精度を高めるためにシステムの複雑さを増大させることは不要である。単にアンサンブルのサイズを増大することによってより良い結果を達成するこの方法は、例えば、アンサンブルの選択に重点を置くその他の方法とは異なる。
【図面の簡単な説明】
【００４６】
【図１】本発明による仮想センサシステムの実施形態を示すブロック図である。
【図２】５０個の個別推定値（細線）と、実際の値（太い破線）と、アンサンブル出力（太い連続線）との間の比較を示すグラフである。
【図３】右に向かってアンサンブルサイズが増加する状態で、ＮＯｘを測定する本発明による仮想センサシステムの実施形態の性能をｐｐｍ単位で示す図である。
【図４】設備キャリブレーションを示す図である。
【図５】本発明の実施形態によるＮＯｘ測定の入力パラメータ及び値を示す図である。
【図６】１０入力のテストデータに関するＰＥＭＳ（予測排出量監視システム）性能を示す図である。
【図７】８入力のテストデータに関するＰＥＭＳ性能を示す図である。
【図８】７２８個の個別の出力（赤）、実際の値（緑）、及び、アンサンブル出力（青）の間の比較を示す図である。
【図９】本発明による仮想センサシステムの実施形態におけるアンサンブルの平均絶対誤差（ＭＡＥ）を示す図である。
【図１０】本発明の実施形態により仮想センサシステムが連結されることが可能な方法を示す図である。
【発明を実施するための形態】
【００４７】
図１は本発明による燃焼プロセス（ＣＰ）から生じるガスの量（Ｇ）を測定するため使用される仮想センサシステムの実施形態のブロック図である。
【００４８】
本発明の実施形態では、燃焼プロセス（ＣＰ）から生じるガスの量（Ｇ）の推定のためのアンサンブルに基づく仮想センサシステム（ＶＳ）は、２つ以上の経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）を備え、経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の１つずつがガスの量（Ｇ）を推定するため用意され、結合関数（ｆ）が、個別の経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の１つずつからの信号出力値（ｙ１、ｙ２、．．．、ｙｍ）より正確であるガスの量（Ｇ）の推定値を提供するために経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）からの結果を結合するため用意されている。ガスの量（Ｇ）は、当業者によって理解されるように濃度又は質量排出量として与えられることができる。測定されるべき燃焼プロセスの中で生成されるガスの例は、ＮＯｘ、ＣＯ２、Ｏ２などである。しかしながら、本発明は、当業者によって理解されるように、その他のプロセスからのその他のガスの量を測定するために使用されてもよい。
【００４９】
より詳細には、本発明の本実施形態では、経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の１つずつは、燃焼プロセス（ＣＰ）からの経験的データ（ＥＤ）を使用して訓練されるように用意されている。本発明の実施形態では、経験的データは、仮想センサシステム（ＶＳ）が配置された燃焼プロセス（ＣＰ）からの履歴測定データである。非測定量からの経験的データ（ＥＤ）は、図１に示されるように、センサ値（ＩＡ及びＩＢ）を用いる一時的に設置されたセンサシステム（ＳＡ及びＳＢ）を使って、ならびに、固定センサ（Ｓ１、Ｓ２、．．．、Ｓｍ）と組み合わされた実際の測定キャンペーンから、研究室分析の記録から、又は、オンラインで実行するには計算のために過度に費用がかかる複雑な分析モデルを使った詳細な推定値から導出されることができる。しかしながら、訓練データは、当業者によって理解されるように、その他の類似したプロセスに由来することも可能である。訓練データは、すべての経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）に対し同一でもよく、相違してもよく、例えば、すべてのプロセス測定量が経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の１つずつ訓練データのため含まれるとは限らない。これは、経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の間に多様性を設ける１つの方法である。当業者によって理解されるように、異なる初期化パラメータを設定することによってそれぞれに異なって初期化されても構わない。
【００５０】
各経験的モデルは、プロセス（ＣＰ）の１つ以上のセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）から１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を受け取り、信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて、該ガスの量（Ｇ）を表現する各経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）からの信号出力値（ｙ１、ｙ２、．．．、ｙｍ）を計算するためさらに用意されている。その上、仮想センサシステム（ＶＳ）は、各経験的モデルから信号出力値（ｙ１、ｙ２、．．．、ｙｍ）を受け取り、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の関数として、ガスの量（Ｇ）を表現する仮想センサ出力値（ｙＲ）を連続的に計算するため用意された結合関数（ｆ）を備える。
【００５１】
実施形態では、本発明は、１つ以上のセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）からの１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）から、燃焼プロセス（ＣＰ）から生じるガスの量（Ｇ）を推定する方法である。この方法は、プロセス（ＣＰ）からの経験的データを使って経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）のアンサンブルを訓練するステップと、訓練された経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）にプロセス（ＣＰ）の１つ以上のセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）からの１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を供給するステップと、信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）において、ガスの量（Ｇ）を表現する信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の計算を実行するステップと、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）を連続的に結合し、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の関数として、ガスの量（Ｇ）を表現する仮想センサ出力値（ｙＲ）を計算するステップとを備える。
【００５２】
本発明の実施形態では、経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）又は内部ノードのすべてが同一構造を有することがある。この設定は、所要の個数の内部ノードがテンプレートノードに基づいて仮想センサシステムの中で簡単にインスタンス化できるという利点がある。本実施形態では、経験的モデルの対応する入力及び出力のフォーマットも同様に同一でもよく、すなわち、経験的モデルＮＮ１の入力１のフォーマットは、経験的モデルＮＮ２〜ＮＮｎなどの入力１のフォーマットと同じである。
【００５３】
ノードはすべてが燃焼プロセスのセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）から同じ信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）のセットを受け取るため用意されてもよい。センサからの信号はすべてのノードに配布され、特別な場合を取り扱う余計な作業が回避される。
【００５４】
経験的モデリングは本文献中で上述され、異なる技術を用いて実施されることができる。本発明の実施形態では、経験的モデルはニューラルネットワークである。
【００５５】
仮想センサシステムの結合関数（ｆ）は、異なる基準に基づいて出力値（ｙＲ）を計算するため用意されてもよい。本発明の実施形態では、結合関数（ｆ）は、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の平均値として、仮想センサ出力値（ｙＲ）を連続的に計算するため用意される。平均値は、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）の幾何もしくは算術平均値として、中央値として、又は、２つの中間値の平均のような平均と中央値との組み合わせとして計算されることができる。中央値計算を用いる本発明による仮想センサシステムの性能は、中央値計算が使用されるとき、出力が一般的に個別のノイズ又は不規則性による影響を受けないので、殆どの場合に、平均値計算より優れていることがわかる。
【００５６】
このアプローチは、ニューラルネットワークのような経験的回帰モデルの性能の中に予想することができる内在する変動を打ち消す。この変動の原因は、様々の程度の訓練データの過剰適合（すなわち、データの中でノイズをモデル化するという結果になる）、訓練前のニューラルネットワーク・パラメータの典型的にランダムなインスタンス化、及び、ニューラルネットワークモデルをデータに当てはめるため使用される非決定論的な勾配降下技術から生じることがある。
【００５７】
図２は３つの要因の組み合わせから生じることがある一種の変動を示し、一組のニューラルネットワーク仮想センサモデルが沖合石油プラットフォームから排出された水中の残留石油濃度を推定するため開発された。同図は、５０個のモデルの個別出力と、推定された実際の期待値と、５０個の個別推定値のアンサンブル結合とを表している。
【００５８】
本発明の実施形態では、結合関数（ｆ）は、経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）からの信号出力値（ｙ１、ｙ２、．．．、ｙｍ）に加えて、プロセスセンサ（Ｓ１、Ｓ２、．．．、Ｓｍ）から直接的に該信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）のうちの１つ以上を受け取り、仮想センサ出力値（ｙＲ）を計算するため用意されている。本発明の本実施形態では、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）は、１つ以上の信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）に基づいて、個別、動的に重み付けされる。動的重み付けは、センサ又はセンサからの伝送線のうちの１つ以上に関連したノイズ及び外乱からの仮想センサ出力値への影響を低減し得る。本発明の関連した実施形態では、結合関数（ｆ）は、信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）を受け取り、信号出力値（ｙ１、ｙ２、．．．、ｙｍ）と信号入力値（Ｉ１、Ｉ２、．．．、Ｉｍ）と経験的モデル（ＮＮＲ）の構造とに基づいて、仮想センサ出力値（ｙＲ）を計算するため用意された経験的モデル（ＮＮＲ）である。
【００５９】
図３は、本発明による仮想センサシステム（ＶＳ）の実施形態の精度又は性能がノードの個数に伴って高まる状況を表している。所与のアプリケーションにおける仮想センサシステムに対する性能要件は変化することがあり、不必要に多数のノードは仮想センサシステム（ＶＳ）の初期化プロセスを鈍化することがある。本発明の実施形態では、仮想センサシステム（ＶＳ）は、特定の性能基準を提供するため、ある程度の個数の該経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）をインスタンス化することが可能であるように用意されている。本発明の実施形態では、仮想センサシステム（ＶＳ）は、ガスの量（Ｇ）を表現する仮想センサ出力値（ｙＲ）の所定の性能要件を達成するため、必要とされる該経験的モデル（ＮＮ１、ＮＮ２、．．．、ＮＮｎ）の個数を動的に割り付けるように用意されている。性能要件は、例えば、ｐｐｍ（パーツ・パー・ミリオン）単位で与えられることがある。
【００６０】
本発明の実施形態では、仮想センサシステム（ＶＳ）は図１０から理解できるように連結されてもよい。ここでは、燃焼プロセスからのＯ２が本発明による仮想センサシステムの実施形態において推定される方法が明らかにされている。Ｏ２濃度は、燃焼室コンフィギュレーション、第８段抽気流、ブリード弁気流、燃料流、及び、軸流コンプレッサ気流に基づいて推定される。推定されたＯ２濃度は、付加的なプロセス測定値、すなわち、火炎温度、気圧、周囲湿度、及び、周囲温度と共に、ＮＯｘ仮想センサシステムへの入力として使用される。仮想センサシステムの連結は、システムの性能を高める上に、経験的モデルの構造及びシステムの訓練を簡素化することがある。
【００６１】
異なるアンサンブルサイズを使用する本発明のテストは、アンサンブル性能がアンサンブルサイズの増加に伴って高まることを明らかにしてきた。単にアンサンブルのサイズを増加させることによってより優れた結果を達成するこの方法は、例えば、アンサンブルの選択に重点を置くその他の方法と相違している。これらのテストでは、アンサンブルサイズは、最小値である２個のコンポーネントモデルから最大値である５９個のコンポーネントモデルまで変化させられた。各アンサンブルサイズに対し、１００回の個別トライアルが実施され、（平均絶対誤差として表現された）結果として生じる性能が計算された。収集された結果は図３に要約され、値が約２０〜３０個の個体数のアンサンブルサイズで漸減することを表している。図８は７００個より多くの出力を備えた極端な事例を示している。
【００６２】
ＰＥＭＳ（パラメトリック排出量監視システム）テクノロジーは、当初、ガスタービンの窒素酸化物（ＮＯｘ）排出量を監視するＣＥＭＳ（連続排出量監視システム）に対し、より高い費用対効果の代替案を得るために開発された。ＣＥＭＳは、物理的な汚染物質分析装置測定量を使用して、ガスもしくは粒状物質濃度又は排出率の決定のため必要な全設備である。ＮＯｘ排出量を直接的に測定する代わりに、ＰＥＭＳは、燃焼温度、圧力、及び、燃料消費量のような主要な運転パラメータからＮＯｘ排出量を計算するので、従ってすべての観点において仮想センサであるとみなすことができる。
【００６３】
本発明の実施形態では、ノルウェー大陸棚の沖合石油プラットフォームで動作するＧＥ  ＬＭ２５００  ＤＬＥガスタービンは、排出量を最小限に抑えるため、最適パラメータ設定を特定するためにマッピングされた。マッピングを実行するため、物理的な排出量監視設備が設置され、タービンは最適パラメータ設定が特定される負荷の範囲で駆動される。結果はタービン負荷をパラメータ設定にマッピングする表と考えることができる。
【００６４】
マッピング中にタービンは連続的に回転させられるので、取得されたデータはＰＥＭＳの構築のため理想的ではない。推奨される手続は、タービンのマッピングが終了した後、異なるタービン負荷で付加的なデータを収集することである。しかしながら、この手続は、余計なダウンタイムコストを発生させる可能性がある。
【００６５】
取得されたデータは図４に示され、１秒間隔でサンプリングされ、％単位のＣＯ２、％単位のＯ２、ｐｐｍ単位のＣＯ、ｐｐｍ単位のＴＨＣ、ｐｐｍ単位のＮＯｘ、及び１５％のＯ２に対し補正されたｐｐｍ単位のＮＯｘの値である。
【００６６】
ＰＥＭＳモデリングのため使用されたデータは、測定設備の２回の強調表示されたキャリブレーション間の約５時間のデータであった。
【００６７】
本実施形態では、選択されたタービンからのプロセスデータが２つの異なるタービン制御システム（ＡＢＢ及びＷｏｏｄｗａｒｄ）から利用できた。このデータは、陸上のヒストリアン・データシステムへ一部しか反映されず、すなわち、タービンに関連付けられたすべての測定量が陸上で利用可能であるとは限らなかった。
【００６８】
ＡＢＢシステムからの殆どの測定量はデータヒストリアンに反映されたが、Ｗｏｏｄｗａｒｄシステムからの測定量は、制御システムを停止しプログラム変更を行うことなしに反映させることができなかったので使用されなかった。着目しているタービンに対し、４０個の測定量が陸上のプロセスヒストリアンにおいて最終的に利用可能であった。
【００６９】
排出量データは、異なるクロック、制御システムのタイムスタンプ及び陸上データヒストリアンのタイムスタンプと一致しないタイムスタンプを備えた携帯型コンピュータシステムで取得された。排出量データとプロセスデータとを同期させるため、２つのデータ系列が、図４に示されているように、キャリブレーション点を表し、プロセス時系列と排出量時系列の両方において整合性を示した有意な変化を視覚的にマッチングすることにより手動で同期させられた。タービンマッピングはデータ内に確実なパターンを作り出すため、この場合に、この方法は可能だった。その他の場合、この手動式同期は実行することが非常に困難であり得るので、従って、使用されたすべてのデータロギング設備のクロックの正確な同期が必要とされる。
【００７０】
上述の排出量データ及びプロセスデータが与えられると、ある程度の個数の試行的なＰＥＭＳモデルが代替的なＰＥＭＳ設計及びコンフィギュレーションを模索するために開発された。選択されたタービンのため利用可能であるすべてのプロセス測定量の中から、１０個の測定量からなるサブセットがＰＥＭＳへの入力に使用されるように選択された。
【００７１】
選択された入力は次のとおりである。
燃料ガス供給圧
ガスジェネレータ・コンプレッサ・吐き出し圧−ＰＳ３
ガスジェネレータ・排気温度−Ｔ５４
パワータービン排気温度
燃料ガスレギュレータ位置（インナーリング）
燃料ガスレギュレータ位置（パイロットリング）
燃料ガスレギュレータ位置（アウターリング）
第８段ブリード弁位置
ＣＤＰブリード弁位置
ガスジェネレータ空気取り入れ口温度
【００７２】
図５は、着目している５時間において、これらの１０個の測定量に対応する時系列の概要を示す。
【００７３】
これらの入力が与えられると、ＰＥＭＳが、多数のモデルが個別に構築され、集計されたアンサンブルモデルの中に結合される。この場合、アンサンブルＰＥＭＳモデルは、２０個の個別ＰＥＭＳモデルの結合であった。
【００７４】
これらのモデルを訓練してテストするため、プロセスデータ及び排出量データの５時間の原データセットは、訓練セット、検証セット、及びテストセットに分割され、訓練セットはモデルを構築するため使用され、検証セットはモデリングを制御するため（すなわち、モデルを訓練データに過剰適合させることを回避するため）使用され、テストセットはモデル性能を評価するため使用された。
【００７５】
原データセットを分割するため、データの４０％が訓練のため無作為に選択され、３０％が検証のため無作為に選択され、残りの３０％がテストのため保たれた。
【００７６】
テストデータセット（すなわち、モデルを構築するため訓練中に使用されなかったデータ）に基づくＰＥＭＳ性能の結果は図６にグラフで示され、０．２８４７２ｐｐｍの平均絶対誤差を与える。ここで、
【数１】
2010537192.tif 000003
であり、ｙｉは期待値であり、
【数２】
2010537192.tif 000004
はモデル推定値である。
【００７７】
このＰＥＭＳアプローチのＳＡＣ（非ＤＬＥ）タービンへの用途の実現可能性を調査するため、選択された測定量のうちの２つ（すなわち、旧式の標準的な燃焼器ＳＡＣタービンでは利用できない２つのブリード弁位置）が除外され、８個の測定量だけが図７に示されるように入力に取り込まれる付加的なテストが実行された。
【００７８】
この場合のテストデータセットに関するＰＥＭＳ性能の結果は図９にグラフで示され、０．３７４５３ｐｐｍのＭＡＥを与える。
【００７９】
８入力を用いるＰＥＭＳの平均誤差は、１０入力全部を用いるＰＥＭＳの平均誤差より約３０％高い。しかしながら、絶対値において、８入力ＰＥＭＳの誤差は、３ｐｐｍ未満という（ＧＥ  ＬＭ２５００  ＤＬＥのような）低ＮＯｘタービンに対する現在の精度要件と比較すると依然として低い。
【００８０】
本実施形態では、訓練データとテストデータとの間に高い類似性が存在する。たとえ、訓練データ及びテストデータが完全に分離されたデータセットであるとしても（原データセットから非復元無作為抽出されているとしても）、訓練データ及びテストデータは依然として同じ時系列から取得され、テストセット中の点が訓練セットの中に非常に類似した点を有する可能性は非常に高い。これにもかかわらず、取得された０．２８ｐｐｍと要求された３ｐｐｍとの間の精度の「マージン」は、本実施形態にある程度の信頼を与えるために十分に大きい。
【００８１】
別の実施形態では、複数のモデルが生成され、アンサンブルの一部になる特有のモデルを選択する仕組みが使用される。このことは、静的に、すなわち、訓練フェーズの後に１回に限り、始めに望ましくないモデルを廃棄することにより行われるか、又は、動的に、すなわち、現在の運転状態が与えられると、その運転状態の中又は付近においてより優れた性能を実証したコンポーネントモデルを優先する重み付けスキームを導入することにより行われる。
【００８２】
さらに別の実施形態では、ハイブリッド型アンサンブルモデル、すなわち、コンポーネントモデルが必ずしも同じタイプではないが、例えば、ニューラルネットワークと同様に、その他の回帰モデル、又は経験的モデル及び分析的モデルの組み合わせとで構成されるアンサンブルが使用される。
【図１】
2010537192.tif 000005
【図２】
2010537192.tif 000006
【図３】
2010537192.tif 000007
【図４】
2010537192.tif 000008
【図５】
2010537192.tif 000009
【図６】
2010537192.tif 000010
【図７】
2010537192.tif 000011
【図８】
2010537192.tif 000012
【図９】
2010537192.tif 000013
【図１０】
2010537192.tif 000014
【国際調査報告】
2010537192.tif 000015
2010537192.tif 000016
2010537192.tif 000017
2010537192.tif 000018