(57)【特許請求の範囲】
【請求項１】
生産システムにおけるシステム出力を最適化するための方法であって、
プロセッサによって、１つまたは複数の観察された初期入力変数を、環境変数とシステム応答変数とを含む複数の出力変数に分離し、
前記プロセッサを使用して、１つまたは複数の初期制御変数と前記システム応答変数との間の関係を決定するノンパラメトリック推定を構築し、
前記決定された関係を使用して、大域の入力／出力マッピング関数を推定し、前記環境変数の範囲を推定し、
前記入力／出力マッピング関数を最大化することによって、前記初期制御変数から１つまたは複数の最適制御変数を特定し、さらに前記環境変数の値を特定し、
前記最適制御変数のうちの１つまたは複数を生産システムに組み込んで、前記生産システムの生産出力を増大させ、
前記最適制御変数は、前記初期制御変数の状態に関連した複数のＫＰＩ（key performance indicator）値を用いて、いずれの初期制御変数の値が、前記大域の入力／出力マッピング関数の推定の結果として最大ＫＰＩ値になるかを検査することによって決定される、方法。
【請求項２】
相関分析を実行して、前記出力変数のそれぞれについての相関スコアを決定することをさらに含み、前記相関分析は、
前記出力変数のそれぞれと前記初期制御変数との間でピアソンの積率相関を計算し、
前記出力変数のそれぞれについての前記相関スコアを閾値と比較し、
前記閾値を超える相関スコアを有する各出力変数をシステム応答変数に指定し、前記閾値以下の相関スコアを有する各出力変数を環境変数に指定することを含む、請求項１に記載の方法。
【請求項３】
前記１つまたは複数の観察された初期入力変数を分離することは、因果分析を実行して前記出力変数を識別することを含み、前記因果分析は、
制御変数の数を増大させ、
前記出力変数に予測分析を実行することを含み、
前記制御変数の数を増やすことによって出力変数における予測誤差が減る場合は、前記出力変数がシステム応答変数として指定され、前記制御変数の数を増やすことによって前記出力変数における前記予測誤差がわずかしか変化しない場合には、所定の閾値に従って、前記出力変数が環境変数として指定される、請求項１に記載の方法。
【請求項４】
グラフィカルユーザインターフェースを使用して、前記１つまたは複数の最適制御変数をユーザに推奨することをさらに含む、請求項１に記載の方法。
【請求項５】
前記環境変数の範囲が、履歴記録、予測アルゴリズム、および領域知識からなるグループに基づいている、請求項１に記載の方法。
【請求項６】
前記生産システムが火力発電所であり、前記最適制御変数のうちの１つまたは複数を組み込むことが、前記火力発電所での投入燃料の成分構成を変え、前記火力発電所の出力効率を上げる、請求項１に記載の方法。
【請求項７】
前記環境変数が、外部の刺激を表し、前記初期入力変数による影響を受けることができない、請求項１に記載の方法。
【請求項８】
生産システムにおけるシステム出力を最適化するシステムであって、
１つまたは複数の観察された初期入力変数を、環境変数とシステム応答変数とを含む複数の出力変数に分離し、
初期制御変数と前記システム応答変数との間の関係を決定するノンパラメトリック推定を構築し、
前記決定された関係を使用して、大域の入力／出力マッピング関数を推定し、前記環境変数の範囲を推定し、
前記入力／出力マッピング関数を最大化することによって、前記初期制御変数から１つまたは複数の最適制御変数を特定し、さらに前記環境変数の値を特定し、
前記最適制御変数のうちの１つまたは複数を生産システムに組み込んで、前記生産システムの生産出力を増大させる
ように構成されたプロセッサと、
前記１つまたは複数の観察された初期入力変数、前記１つまたは複数の初期制御変数、前記複数の出力変数、および前記１つまたは複数の最適制御変数を記憶するためのメモリと
を有し、
前記最適制御変数は、前記初期制御変数の状態に関連した複数のＫＰＩ（key performance indicator）値を用いて、いずれの初期制御変数の値が、前記大域の入力／出力マッピング関数の推定の結果として最大ＫＰＩ値になるかを検査することによって決定される、システム。
【請求項９】
前記プロセッサは、相関分析を実行して、前記出力変数のそれぞれについての相関スコアを決定するようにさらに構成され、
前記相関分析が、
前記出力変数のそれぞれと前記初期制御変数との間でピアソンの積率相関を計算し、
前記出力変数のそれぞれについて前記相関スコアを閾値と比較し、
前記閾値を超える相関スコアを有する各出力変数をシステム応答変数に指定し、前記閾値以下の相関スコアを有する各出力変数を環境変数に指定することを含む、請求項８に記載のシステム。
【請求項１０】
前記プロセッサは、因果分析を実行して前記出力変数を識別するようにさらに構成され、前記因果分析が、
制御変数の数を増大させ、
前記出力変数に予測分析を実行することを含み、
前記制御変数の数を増やすことによって出力変数における予測誤差が減る場合は、前記出力変数がシステム応答変数として指定され、前記制御変数の数を増やすことによって前記出力変数における前記予測誤差がわずかしか変化しない場合には、所定の閾値に従って、前記出力変数が環境変数として指定される、請求項８に記載のシステム。
【請求項１１】
前記１つまたは複数の最適制御変数をユーザに表示するように構成された、グラフィカルユーザインターフェースをさらに有する、請求項８に記載のシステム。
【請求項１２】
前記環境変数の範囲が、履歴記録、予測アルゴリズム、および領域知識からなるグループに基づいている、請求項８に記載のシステム。
【請求項１３】
前記生産システムが火力発電所であり、前記プロセッサが、前記最適制御変数のうちの１つまたは複数を組み込むときに、前記火力発電所での投入燃料の成分構成を変え、前記火力発電所の出力効率を上げるようにさらに構成された、請求項８に記載のシステム。
【請求項１４】
前記環境変数が、外部の刺激を表し、前記初期入力変数による影響を受けることができない、請求項８に記載のシステム。
【請求項１５】
前記システム応答変数が、人間の入力および前記システムの構造によって決定される、請求項８に記載のシステム。
【請求項１６】
生産システムにおけるシステム出力を最適化するためのコンピュータ読取り可能なプログラムを含む非一時的なコンピュータ読取り可能な記憶媒体であって、前記コンピュータ読取り可能なプログラムは、コンピュータ上で実行されると、
プロセッサにより、１つまたは複数の観察された初期入力変数を、環境変数とシステム応答変数とを含む複数の出力変数に分離するステップと、
前記プロセッサを使用して、１つまたは複数の初期制御変数と前記システム応答変数との間の関係を決定するノンパラメトリック推定を構築するステップと、
前記決定された関係を使用して、大域の入力／出力マッピング関数を推定し、前記環境変数の範囲を推定するステップと、
前記入力／出力マッピング関数を最大化することによって、前記初期制御変数から１つまたは複数の最適制御変数を特定し、さらに前記環境変数の値を特定するステップと、
前記最適制御変数のうちの１つまたは複数を生産システムに組み込んで、前記生産システムの生産出力を増大させるステップと
を、前記コンピュータに実行させ、
前記最適制御変数は、前記初期制御変数の状態に関連した複数のＫＰＩ（key performance indicator）値を用いて、いずれの初期制御変数の値が、前記大域の入力／出力マッピング関数の推定の結果として最大ＫＰＩ値になるかを検査することによって決定される、コンピュータ読取り可能な記憶媒体。
【請求項１７】
前記コンピュータ読取り可能なプログラムは、コンピュータ上で実行されると、相関分析を実行するステップを前記コンピュータにさらに実行させ、前記相関分析が、
前記出力変数のそれぞれと前記初期制御変数との間でピアソンの積率相関を計算することと、
前記出力変数のそれぞれについての前記積率相関に基づく相関スコアを閾値と比較することと、
前記閾値を超える相関スコアを有する各出力変数をシステム応答変数に指定し、前記閾値以下の相関スコアを有する各出力変数を環境変数に指定することと
を含む、請求項１６に記載のコンピュータ読取り可能な記憶媒体。
【請求項１８】
前記１つまたは複数の観察された初期入力変数を前記分離することが、因果分析を実行して前記出力変数を識別することを含み、前記因果分析が、
制御変数の数を増大させることと、
前記出力変数に予測分析を実行することと
を含み、
前記制御変数の数を増やすことによって出力変数における予測誤差が減る場合は、前記出力変数がシステム応答変数として指定され、前記制御変数の数を増やすことによって前記出力変数における前記予測誤差がわずかしか変化しない場合には、所定の閾値に従って、前記出力変数が環境変数として指定される、請求項１６に記載のコンピュータ読取り可能な記憶媒体。
【請求項１９】
前記環境変数の範囲が、履歴記録、予測アルゴリズム、および領域知識からなるグループに基づいている、請求項１６に記載のコンピュータ読取り可能な記憶媒体。
【請求項２０】
前記生産システムが火力発電所であり、前記最適制御変数のうちの１つまたは複数を組み込むことが、前記火力発電所での投入燃料の成分構成を変え、前記火力発電所の出力効率を上げる、請求項１６に記載のコンピュータ読取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願情報
本出願は、２０１５年７月１７日出願の仮米国特許出願第62/193,629号の優先権を主張し、参照によりその全体をここに組み込む。
【０００２】
本発明は、システム出力の最適化に関し、より詳細には、生産システムにおける主要業績評価指標（key performance indicator）の最適化に関する。
【背景技術】
【０００３】
関連技術の説明
多くの生産システム（たとえば、発電所、化学プラントなど）では、システム出力は、制御可能な１つまたは複数の入力変数によって決定される。たとえば、火力発電所では、投入燃料（たとえば石炭など）の原料構成が、火力発電所の出力効率に大きく影響を及ぼすことがある。この出力効率は、通常、発電される電力と消費される燃料の量との比として計算される。したがって、このようなプラントでの重要な目標は、システム効率を改善するために、所望の特性を有する投入燃料をどのように選択するかを決定することである。
【０００４】
このような静的入力に加えて、生産プロセス中に何らかの変数を調整する動的制御方式もある。たとえば、生産システムの多くでは、オペレータがパラメータ（たとえば、電圧、電流の回転角、またはある種の成分など）を調整する必要がある。したがって、出力品質（たとえば、主要業績評価指標（ＫＰＩ）の値など）を最大化するために、こうした制御変数の組合せをどのように決定するかは、生産の管理および最適化にとって大きい関心事である。
【発明の概要】
【０００５】
本発明の原理の一態様によれば、生産システムにおけるシステム出力を最適化するための方法が提供される。この方法は、プロセッサにより、１つまたは複数の初期入力変数を複数の出力変数に分離することを含んでおり、この出力変数は、環境変数およびシステム応答変数を含む。この方法はまた、プロセッサを使用して、１つまたは複数の初期制御変数とシステム応答変数との間の関係を決定するノンパラメトリック推定（nonparametric estimation）を構築することと、その決定した関係を使用して、大域の入力／出力マッピング関数を推定し、環境変数の範囲を推定することとを含む。この方法はさらに、入力／出力マッピング関数と環境変数の範囲とを最大化することによって、初期制御変数から１つまたは複数の最適制御変数を生成することを含む。この方法はさらに、最適制御変数の１つまたは複数を生産システムに組み込んで、この生産システムの生産出力を増大させることを含む。
【０００６】
本発明の原理の別の態様によれば、生産システムにおけるシステム出力を最適化するためのシステムが提供される。このシステムは、１つまたは複数の初期入力変数を、環境変数とシステム応答変数とを含む複数の出力変数に分離し、初期制御変数とシステム応答変数との間の関係を決定するノンパラメトリック推定を構築し、この決定された関係を使用して、大域の入力／出力マッピング関数を推定し、環境変数の範囲を推定し、入力／出力マッピング関数と環境変数の範囲とを最大化することによって、これら初期制御変数から１つまたは複数の最適制御変数を生成し、該最適制御変数の１つまたは複数を生産システムに組み込んで、この生産システムの生産出力を増大させるように構成されたプロセッサを備えている。このシステムはまた、１つまたは複数の該初期入力変数、１つまたは複数の該初期制御変数、複数の該出力変数、および１つまたは複数の該最適制御変数を記憶するためのメモリを備えている。
【０００７】
本発明の原理のさらに別の態様によれば、生産システムにおけるシステム出力を最適化するためのコンピュータ読取り可能なプログラムを含んでいる非一時的なコンピュータ読取り可能な記憶媒体が提供され、このコンピュータ読取り可能なプログラムは、コンピュータ上で実行されると、プロセッサにより、１つまたは複数の初期入力変数を、環境変数とシステム応答変数とを含む複数の出力変数に分離するステップと、このプロセッサを使用して、１つまたは複数の初期制御変数とシステム応答変数との間の関係を決定するノンパラメトリック推定を構築するステップと、この決定された関係を使用して、大域の入力／出力マッピング関数を推定し、環境変数の範囲を推定するステップと、入力／出力マッピング関数と環境変数の範囲とを最大化することによって、これら初期制御変数から１つまたは複数の最適制御変数を生成するステップと、最適制御変数の１つまたは複数を生産システムに組み込んで、この生産システムの生産出力を増大させるステップとを、コンピュータに実行させる。
【０００８】
これらのおよび他の特徴および利点は、添付の図面に関連して読まれるはずである、それらの例示的な実施態様の以下の詳細な説明から明らかとなるであろう。
【図面の簡単な説明】
【０００９】
本開示は、以下の図を参照して好ましい実施形態の以下の説明において詳細を提供する。
【００１０】
【図１】図１は、本発明の原理による、主要業績評価指標（ＫＰＩ）最適化ワークフローを実行するためのシステム／方法を示す、ハイレベルのブロック／流れ図である。
【図２】図２は、本発明の原理による、ノンパラメトリックフレームワークによってシステム用のＫＰＩ値を最適化するための方法を示す流れ図である。
【図３】図３は、本発明の原理による、ノンパラメトリックフレームワークによってシステム用のＫＰＩ値を最適化するための方法を示す流れ図である。
【図４】図４は、本発明の原理の一実施形態による、本発明の原理を適用する例示的な処理システムを示すブロック図である。
【図５】図５は、本発明の原理の一実施態様による、ＫＰＩ最適化ワークフローを実行する例示的なシステムを示すブロック図である。
【発明を実施するための形態】
【００１１】
本発明の原理によれば、変数依存性と多様体正則化（manifold regularization）を両方とも使用するノンパラメトリックフレームワークを組み込むことによって、生産システムにおいて主要業績評価指標（ＫＰＩ）値を最適化するためのシステムおよび方法が提供される。これらのシステムおよび方法は、生産システムにおける出力の量または質を改善することができる。目標を明確に定義および測定することができる限り、提供されるこのシステムおよび方法を使用し、制御変数を調整して、該出力ＫＰＩ値を最大化することができる。
【００１２】
従来の方法は、通常、ニューラルネットワークを使用して、入力および出力のマッピング関数を推定し、これらの関数は、数多くのパラメータを調整することを必要とし、局所最適解に陥る場合もある。本発明の原理の一実施形態によれば、大域的最適解を提供できる、多様体正則化されたカーネル回帰（kernel regression）のフレームワークが使用される。この方法により、制御パラメータを最適なＫＰＩ値に位置付けるのが、さらに便利になる。さらに、本発明の原理の一実施形態によれば、データ駆動型の数学的な方法を使用して、入力変数によってあまり良好に説明できない変数を分離し、それによってオンライン最適化を改善することができる。
【００１３】
一実施形態では、入力変数と出力変数とを互いに関連付ける、ノンパラメトリック推定が構築される。入力変数と出力変数との間の関係は、事前発見プロセス（pre-discovery process）によって、または領域知識（domain knowledge）を用いて、該入力変数の依存構造を考慮する。一実施形態では、観察された変数は出力変数に分離され、この出力変数は、システム状態変数および環境変数を含んでいる。システム状態変数は、制御変数に依存しており、環境変数は、この制御仕組みに影響されないが、代わりに外界からの刺激に依存している。
【００１４】
一実施形態では、多様体正則化が損失関数（loss function）とともに追加されて、入力／出力マップを推定する。
【００１５】
アプリケーションにおいて、本発明の原理の一実施形態によれば、環境変数のオンライン推定が作成される。この入力は、推定された関係関数に挿入され、最適制御変数が、この条件下に置かれる。
【００１６】
本明細書に記載の実施形態は、全体的にハードウェアでもよく、全体的にソフトウェアでもよく、またはハードウェア要素とソフトウェア要素を両方含んでいてもよい。好ましい実施形態では、本発明は、ソフトウェアで実装されており、このソフトウェアはファームウェア、常駐ソフトウェア、マイクロコードなどを含むが、これらに限定されない。
【００１７】
実施態様は、コンピュータもしくは任意の命令実行システムが使用する、またはそれに関連するプログラムコードを提供する、コンピュータ使用可能またはコンピュータ読取り可能な媒体からアクセスできるコンピュータプログラム製品を含んでいてもよい。コンピュータ使用可能またはコンピュータ読取り可能な媒体は、命令実行システム、機器、もしくは装置が使用するか、またはそれらと関連するプログラムを記憶、伝達、伝搬、または移送する任意の装置を備えていてもよい。媒体は、磁気、光、電子、電磁、赤外線、もしくは半導体のシステム（もしくは機器、もしくは装置）、または伝搬媒体であってもよい。媒体は、半導体あるいは固体の記憶装置、磁気テープ、取外し可能なコンピュータディスケット、ランダムアクセスメモリ（ＲＡＭ）、リードオンリメモリ（ＲＯＭ）、硬質磁気ディスク、光ディスクなど、コンピュータ読取り可能な記憶媒体を含んでいてもよい。
【００１８】
記憶媒体または記憶装置がコンピュータに読み込まれて、本明細書に記載の手順を実行する場合、コンピュータの動作を構成および制御するために、コンピュータプログラムはそれぞれ、一般目的または特殊目的のプログラム可能なコンピュータによって読み取ることが可能な、機械読取り可能な記憶媒体または記憶装置（たとえば、プログラムメモリもしくは磁気ディスク）に具体的に記憶してもよい。また、この発明のシステムは、コンピュータプログラムで構成された、コンピュータ読取り可能な記憶媒体において実施されることを考慮してもよく、そのように構成された記憶媒体は、コンピュータに、特定の事前定義された方式で動作して本明細書に記載の機能を実行させる。
【００１９】
プログラムコードを記憶および／または実行するのに適したデータ処理システムは、システムバスを介して記憶素子に直接または間接に結合された、少なくとも１つのプロセッサを備えていてもよい。記憶素子は、プログラムコードを実際に実行する際に利用されるローカルメモリ、大容量記憶装置、および、少なくとも何らかのプログラムコードを一時的に記憶して、実行中に大容量記憶装置からコードを取り出す回数を減らすキャッシュメモリを含むことができる。入力／出力装置すなわちＩ／Ｏ装置（キーボード、表示装置、ポインティング装置などを含むが、これらに限定されない）は、直接または介在するＩ／Ｏ制御装置を介して、このシステムに結合していてもよい。
【００２０】
また、ネットワークアダプタは、データ処理システムが、介在する専用ネットワークまたは公衆ネットワークを介して、他のデータ処理システム、またはリモートプリンタもしくは記憶装置に結合できるように、システムに結合されてもよい。モデム、ケーブルモデム、およびイーサネットカードは、ネットワークアダプタの現在利用可能なタイプのうちの、ほんのいくつかに過ぎない。
【００２１】
次に、図を詳細に参照すると、同じ番号は、同じまたは同様の要素を表しており、最初に図１を参照すると、主要業績評価指標（ＫＰＩ）最適化ワークフローを実行するためのシステム／方法１００を示すハイレベルのブロック図／流れ図が、本発明の原理の一実施形態に従って例示的に示してある。
【００２２】
システム１００は、要素１１０、１１５、１２０、１２５、１３０、１３５、１４０、１４５、１５０、１５５、および１６０から構成されている。一実施形態では、システム１００の主要な５つの構成要素、すなわち、変数セパレータ１１０、第１の関数推定器１２５、第２の関数推定器１３０，オンライン推定器１４５、およびオンラインオプティマイザ１５５からなる。
【００２３】
一実施形態では、変数セパレータ１１０は、観察された変数１０５（たとえば時系列など）を出力変数に分離する。出力変数は、２組の変数、すなわち環境変数１３５およびシステム応答変数１１５にグループ分けされる。環境変数１３５は、外界からの何らかの形の刺激を表す変数であり、制御変数１２０による（直接または間接の）影響を受けることができない。システム応答変数１１５は、制御変数１２０の下でのシステム１００の応答を表す変数である。一実施形態では、システム応答変数１１５は、人間の入力およびシステム１００の構造によって決定される変数である。
【００２４】
一実施形態では、変数セパレータ１１０は、観察された変数１０５と制御された変数１２０との間の相関関係を決定することによって、２組の変数１１５、１３５を分離する。一実施形態では、制御変数１２０とシステム応答変数１１５とを互いに関連付ける変数セパレータ１１０によって、ノンパラメトリック推定が構築される。制御変数１２０とシステム応答変数１１５の間の関係は、事前発見プロセスによって、または領域知識を介して、制御変数１２０の依存構造を考慮に入れる。一実施形態では、制御変数１２０とシステム応答変数１１５との間の関係を決定するために、カーネル回帰技法が使用される。
【００２５】
一実施形態では、相関分析が実行される。全ての出力変数と制御変数１２０との間でピアソンの積率相関（Pearson product-moment correlation）が計算され、その絶対値がスコアとして計算される。この相関は、以下のように表すことができる。
【００２６】
【数１】
0006718500.tif 000002
ここで、「ｘ」は該出力変数であり、「ｙ」は制御変数１２０である。出力変数のそれぞれの相関スコアが閾値と比較される。閾値を超える相関スコアを有する各出力変数が、システム応答変数として指定され、閾値以下の相関スコアを有する各出力変数が、環境変数として指定される。相関分析が実行された後、環境変数１３５は、制御変数１２０との相関が低くなり、システム応答変数１１５は、制御変数１２０との相関が高くなる。１つに複数の制御変数が存在する場合、各出力変数のスコアは、全ての入力変数との相関を加算したものになることに留意されたい。
【００２７】
一実施形態では、因果分析が実行されて、２組の変数１１５、１３５を識別する。入力変数を使用して観察された変数を予測する性能（予測誤差）が検査される。予測のために、さらに多くの制御変数１２０が追加される。制御変数をさらに多くした結果、出力変数の予測中の誤りが少なくなる場合、該出力変数は、システム応答変数１１５として識別される。そうでない場合、さらに多くの制御変数１２０を追加しても、予測誤差に著しい変化がない場合、既定の閾値に従って、対応する出力変数が、環境変数１３５として識別される。
【００２８】
因果分析中には、制御変数１２０は、システム応答変数１１５上での予測確度を高くするが、環境変数１３５を予測する際には著しい差をもたらさない。
【００２９】
一実施形態では、第１の関数推定器１２５が、制御変数１２０とシステム応答変数１１５との間の関係を決定する。一実施形態では、その関係は、関数Ｓ＝ｇ（Ｃ）として表され、ここで、Ｓはシステム変数であり、ｇ（Ｃ）は制御変数Ｃの関数である。
【００３０】
一実施形態では、制御変数１２０、すなわちＣｉ、およびシステム応答変数１１５、すなわちＳｉが、ｉ＝１、２、．．．、Ｎについて提供される。これらの変数を使用して、以下の関数が推定される。
【００３１】
【数２】
0006718500.tif 000003
ここで、ｈは帯域幅であり、各αｉはモデルの係数であり、以下の目的関数を最小化することによって推定することができる。
【００３２】
【数３】
0006718500.tif 000004
ここで、
【００３３】
【数４】
0006718500.tif 000005
は、正則化パラメータであり、ＬはＣｉのペアの類似度行列上のラプラシアン行列（Laplacian matrix）である。
【００３４】
一実施形態では、第２の関数推定器１３０が、大域の入力／出力マッピング関数、
【００３５】
【数５】
0006718500.tif 000006
を推定し、第１の関数推定器１２５において回復された関数を考慮に入れ、ここで、
【００３６】
【数６】
0006718500.tif 000007
はＫＰＩ値１４０であり、Ｅは環境変数１３５である。
【００３７】
一実施形態では、ｉ＝１、２、．．．、Ｎについて、制御変数１２０であるＣｉ、第１の推定モジュール１２５によって提供される関係（すなわち、
【００３８】
【数７】
0006718500.tif 000008
）、環境変数１３５であるＥｉ、およびＫＰＩ値１４０である
【００３９】
【数８】
0006718500.tif 000009
が提供される。これらの変数を使用して、以下の関数が推定される。
【００４０】
【数９】
0006718500.tif 000010
ここで、
【００４１】
【数１０】
0006718500.tif 000011
であり、μは正則化パラメータであり、ＬはＣｉおよびＥｉのペアの類似度行列上のラプラシアン行列である。
【００４２】
一実施形態では、オンライン推定器１４５は、環境変数１３５の範囲を推定する。
【００４３】
一実施形態では、入力／出力マッピング関数
【００４４】
【数１１】
0006718500.tif 000012
が推定された後、最急上昇法（gradient ascent method）を使用して、ＫＰＩ値１４０の局所最適を決定する。このシステムは、ＫＰＩ値１４０のこの決定された局所最適を使用して、最適制御変数
【００４５】
【数１２】
0006718500.tif 000013
を決定する。次いで、システム１００は、たとえばグラフィカルユーザインターフェース１６０を介して、推奨制御変数１２０をユーザに提供する。一実施形態によれば、グラフィカルユーザインターフェース１６０により、ユーザは、最適制御変数
【００４６】
【数１３】
0006718500.tif 000014
の１つまたは複数を生産システム１６５の動作に組み込むことができる。
【００４７】
一実施形態では、オンラインオプティマイザ１５５は、第２の関数推定器１３０において推定される関数を最大化することによって、最適制御変数
【００４８】
【数１４】
0006718500.tif 000015
を特定し、更にオンライン推定器１４５の環境変数１５０の値を特定する。各ＫＰＩ値１４０は、制御変数１２０の状態に関連している。すなわち、各ＫＰＩ値１４０は、１組の制御変数１２０のある値／構成に関数推定器１３０を適用した結果である。最適制御変数は、いずれの制御変数１２０の値が、結果として最大ＫＰＩ値１４０になるのかを検査することによって決定される。
【００４９】
一実施形態では、生産システム１６５の入力と出力の比を増大するために、最適制御変数のうちの１つまたは複数が生産システム１６５に組み込まれる。生産システム１６５は、たとえば、発電所、化学プラント、生産工場などでもよい。
【００５０】
たとえば、火力発電所では、たとえば、投入燃料（たとえば石炭など）の成分構成を変えて、消費される燃料の量に対する発電される電気の出力効率を高めるために、１つまたは複数の制御変数１２０がシステム１６５に組み込まれる。
【００５１】
次に、図２〜図３を参照すると、ノンパラメトリックフレームワークによって、システムのＫＰＩ値を最適化するための方法２００の流れ図が、本発明の原理の一実施形態に従って例示的に示してある。
【００５２】
２１０において、観察された変数１０５および制御変数１２０が、システム１００によって受信される。
【００５３】
２１５において、変数セパレータ１１０は、観察された変数１０５を分離する。一実施形態では、観察された変数１０５は、分離された後、２組の変数、すなわち環境変数１３５およびシステム応答変数１１５にグループ分けされる。環境変数１３５は、外界からの何らかの形の刺激を表す変数であり、制御変数１２０による（直接または間接に）影響を受けることができない。システム応答変数１１５は、制御変数１２０の下でのシステム１００の応答を表す変数である。一実施形態では、システム応答変数１１５は、人間の入力およびシステム１００の構造によって決定される変数である。一実施形態では、環境変数は制御可能ではないが、システム応答変数は制御可能である。一実施形態では、制御不可能な変数１３５の分離は、関数推定の不確実性を低減し、制御変数１２０と出力ＫＰＩ値１４０との間のマッピング関数のモデリング力を改善することができる。
【００５４】
一実施形態では、観察された変数１０５と制御された変数１２０との間の相関関係を決定することによって、２組の変数１１５、１３５が分離される。
【００５５】
２２０において、制御変数１２０とシステム応答変数１１５とを互いに関連付ける変数セパレータ１１０によって、ノンパラメトリック推定が構築される。この関連付けは、事前発見プロセスによって、または領域知識を介して、入力変数の依存構造を考慮に入れる。
【００５６】
２２５において、相関分析が実行される。環境変数１３５は、制御変数１２０との相関が低いが、システム応答変数１１５は、制御変数１２０との相関が高い。
【００５７】
２３０において、２組の変数１１５、１３５を識別するために、因果分析が実行される。因果分析中、制御変数１２０は、システム応答変数１１５上での予測確度を高くするが、環境変数１３０を予測する際には著しい差をもたらさない。ターゲットについての関連した予測変数を識別し、関係のない特徴／変数を除去することによって、該入力／出力マッピング関数の推定がより正確になる。
【００５８】
２３５において、第１の関数推定器１２５は、制御変数１２０とシステム応答変数１１５との間の関係を回復する。一実施形態では、その関係は、関数
【００５９】
【数１５】
0006718500.tif 000016
として表され、ここで、Ｓは該システム変数であり、ｇ（Ｃ）は制御変数Ｃの関数である。
【００６０】
２４０において、第２の関数推定器１３０は、第１の関数推定器１２５において回復された該関数を考慮し、大域の入力／出力マッピング関数
【００６１】
【数１６】
0006718500.tif 000017
を回復し、ここで、
【００６２】
【数１７】
0006718500.tif 000018
はＫＰＩ値１４０であり、Ｅは環境変数１３５である。入力／出力マッピングにおいて、いずれの制御変数１２０が最適な（最高の）ＫＰＩ値１４０を提供できるかを決定することによって、最適制御変数値を決定することができる。一実施形態では、このプロセスは、「ヒルクライミング」、すなわち制御変数１２０の複数の初期選択から開始し、勾配方向（この勾配は、極大点に達するまで、制御変数１２０に関するＫＰＩ値１４０の最大増加の方向として定義される）に沿って繰り返すことに依存する。次いで、極大値の全てを比較することによって、最終的な最適制御変数として、大域的最適制御変数１２０の値を選択することができる。
【００６３】
一実施形態では、ｉ＝１、２、．．．、Ｎについて、制御変数１２０（Ｃｉ）、第１の推定モジュール１２５によって提供される関係（すなわち、
【００６４】
【数１８】
0006718500.tif 000019
）、環境変数１３５（Ｅｉ）、およびＫＰＩ値１４０（
【００６５】
【数１９】
0006718500.tif 000020
）が提供される。これらの変数を使用して、以下の関数が推定される。
【００６６】
【数２０】
0006718500.tif 000021
ここで、
【００６７】
【数２１】
0006718500.tif 000022
であり、μは正則化パラメータであり、ＬはＣｉおよびＥｉのペアの類似度行列上のラプラシアン行列である。
【００６８】
２４５において、オンライン推定器１４５は、たとえば環境変数１３５の履歴値のうちの最大値および最小値に基づいて、環境変数１３５の範囲を推定する。環境変数の範囲は、たとえば、履歴記録、予測アルゴリズム、領域知識などに基づくことができる。
【００６９】
一実施形態では、環境変数１３５は、以下の形式、すなわちリアルタイム予測、ならびに各環境変数１３５の下限および上限（範囲）とすることができる。
【００７０】
２５０において、入力／出力マッピング関数
【００７１】
【数２２】
0006718500.tif 000023
が推定された後、最急上昇法を使用して、ＫＰＩ値１４０の局所最適を決定する。２５５において、システムは、ＫＰＩ値１４０の決定された局所最適を使用して、最適制御変数を決定する。
【００７２】
２６０において、たとえばグラフィカルユーザインターフェース１６０を介して、最適制御変数がユーザに推奨される。
【００７３】
一実施形態では、オンラインオプティマイザ１５５は、第２の関数推定器１３０において推定される関数を最大化することによって、最適制御変数
【００７４】
【数２３】
0006718500.tif 000024
を特定し、更にオンライン推定器１４５の環境変数１５０の値を特定する。
【００７５】
反復性の最急上昇は、
【００７６】
【数２４】
0006718500.tif 000025
として定義することができ、ここで、
【００７７】
【数２５】
0006718500.tif 000026
は、ポイントｘでの関数
【００７８】
【数２６】
0006718500.tif 000027
の勾配である。一実施形態では、局所最適解を避けるために、関数は、無作為に選択された異なる初期解から初期化される。
【００７９】
２６５において、推奨最適制御変数が、たとえば、発電所、化学プラント、生産工場などの生産システム１６５に組み込まれて、たとえば、投入燃料に対して出力電力を増大させるなど、生産システム１６５の出力を増大させる。
【００８０】
ほとんどの生産システム１６５では、何らかの形の制御変数１２０によってシステム出力が決定される。本発明の原理の一実施形態では、推奨制御変数がそのような生産システム１６５に組み込まれて、生産システム１６５の入力と出力の比を増大させる。たとえば、火力発電所では、推奨制御変数がシステム１６５に組み込まれて、たとえば、消費される燃料の量に対して発電される電気の出力効率を高めるように、投入燃料（たとえば石炭など）の成分構成を変える。
【００８１】
次に、図４を参照すると、本発明の原理を適用してもよい例示的な処理システム４００が、本発明の原理の一実施形態に従って例示的に示してある。処理システム４００は、１つまたは複数のネットワーク、バス、または伝送媒体４０２を介して、他の構成要素に動作可能に結合された少なくとも１つのプロセッサ（ＣＰＵ）４０４を備えており、システムバス４０２は、たとえば、ＷｉＦｉ、Ｂｌｕｅｔｏｏｔｈなど１つまたは複数の有線または無線のネットワークを介して任意の構成要素に接続されてもよく、１つまたは複数の伝送媒体は、たとえば同軸ケーブル、光ファイバなど閉じた伝送媒体、および、たとえば無線伝送、マイクロ波伝送など閉じていない伝送媒体を含んでいてもよい。キャッシュ４０６、リードオンリメモリ（ＲＯＭ）４０８、ランダムアクセスメモリ（ＲＡＭ）４１０、入力／出力（Ｉ／Ｏ）アダプタ４２０、サウンドアダプタ４３０、ネットワークアダプタ４４０、ユーザインターフェースアダプタ４５０、および表示アダプタ４６０が、動作可能にシステムバス１０２に結合されている。
【００８２】
第１の記憶装置４２２および第２の記憶装置４２４が、Ｉ／Ｏアダプタ４２０によって、システムバス４０２に動作可能に結合されている。記憶装置４２２および４２４は、ディスク記憶装置（たとえば、磁気ディスク記憶装置または光ディスク記憶装置）、固体磁気装置などのいずれとすることもできる。記憶装置４２２および４２４は、同じタイプの記憶装置または互いに異なるタイプの記憶装置とすることができる。
【００８３】
スピーカ４３２は、サウンドアダプタ４３０によって、システムバス４０２に動作可能に結合されている。送受信機４４２は、ネットワークアダプタ４４０によって、システムバス４０２に動作可能結合されている。表示装置４６２は、表示アダプタ４６０によって、システムバス４０２に動作可能に結合されている。
【００８４】
第１のユーザ入力装置４５２、第２のユーザ入力装置４５４、および第３のユーザ入力装置４５６が、ユーザインターフェースアダプタ４５０によって、システムバス４０２に動作可能に結合されている。ユーザ入力装置４５２、４５４、および４５６は、キーボード、マウス、キーパッド、画像キャプチャ装置、動作検知装置、マイクロホン、こうした装置のうち少なくとも２つの装置の機能を組み込んだ装置など、任意のものとすることができる。ユーザは、ユーザ入力装置４５２、４５４、および４５６を使用して、推奨制御変数を生産システム１６５の機能に組み込んでもよい。もちろん、本発明の原理の精神を維持しながら、他のタイプの入力装置を使用することもできる。ユーザ入力装置４５２、４５４、および４５６は、同じタイプのユーザ入力装置または互いに異なるタイプのユーザ入力装置とすることができる。ユーザ入力装置４５２、４５４、および４５６を使用して、システム４００に対して情報を入力／出力する。
【００８５】
もちろん、処理システム４００は、当業者なら容易に思いつくように、他の要素（図示せず）を備えていてもよく、またある種の要素を省略してもよい。たとえば、当業者に容易に理解できるように、他の様々な入力装置および／または出力装置は、それらの具体的な実装形態に応じて、処理システム４００内に備えることができる。たとえば、様々なタイプの無線および／または有線の、入力装置および／または出力装置を使用することができる。さらに、当業者には容易に理解できるように、様々な構成での追加のプロセッサ、制御装置、メモリなどを利用することもできる。処理システム４００のこれらおよび他の変形形態は、本明細書において提示された本発明の原理の教示が与えられ、当業者が容易に思いつく。
【００８６】
さらに、図１を参照して説明されるシステム１００は、本発明の原理のそれぞれの実施態様を実装するためのシステム１００および方法を示すハイレベルのブロック図／流れ図であることを理解されたい。処理システム４００の一部またはその全ては、図５のシステム５００の要素のうち１つまたは複数の要素で実装してもよい。
【００８７】
さらに、処理システム４００は、たとえば図２〜図３の方法２００の少なくとも一部を含め、本明細書に記載の方法の少なくとも一部を実行してもよいことを理解されたい。同様に、図５のシステム５００の一部またはその全てを使用して、たとえば、図１のシステム１００、図２〜図３の方法２００の少なくとも一部を含め、本明細書に記載の方法の少なくとも一部を実行してもよい。
【００８８】
次に、図５を参照すると、ＫＰＩ最適化ワークフローを実行する例示的なシステム５００が、本発明の原理の一実施態様に従って例示的に示してある。
【００８９】
説明して明らかにするため、システム５００の多くの態様が単数形で説明されているが、これらの態様は、システム５００の説明に関して述べた項目のうち複数の項目に適用することができる。たとえば、システム５００に関しては、単一の表示装置５１５に言及してもよいが、本発明の原理の精神を維持しながら、本発明の原理の各技法に従って２つ以上の表示装置５１５を使用することができる。さらに、表示装置５１５は、本発明の原理の精神を維持しながら複数形に拡張することができる、システム５００に関係のある一態様に過ぎないことが理解される。
【００９０】
一実施形態では、システム５００は、複数の構成要素を含んでもよく、これらの構成要素には、１つまたは複数の回路５０２、制御装置５０４、変数セパレータ５１０、表示装置５１５、第１の関数推定器５２０、第２の関数推定器５３０、オンライン推定器５４０、オンラインオプティマイザ５５０、および／または記憶装置５６０（たとえば、コンピュータ読取り可能な記憶媒体）が含まれてもよい。前述の構成要素は、たとえば、１つまたは複数のネットワーク、バス、または伝送媒体５０１によって接続されていてもよく、１つまたは複数のネットワークは、たとえば、ＷｉＦｉ、Ｂｌｕｅｔｏｏｔｈなど１つまたは複数の有線または無線のネットワークを含んでいてもよく、また、この１つまたは複数の伝送媒体は、たとえば同軸ケーブル、光ファイバなど閉じた伝送媒体、および、たとえば無線伝送、マイクロ波伝送など閉じていない伝送媒体を含んでいてもよい。さらに、データ（たとえば、システム状態のデータ、温度、電圧など）は、本発明の原理の様々な実施形態に従って収集してもよい。バスは、有線接続または無線接続を介して、任意の構成要素に接続されていてもよい。
【００９１】
一実施形態では、変数セパレータ５１０を利用して、観察された変数１０５（たとえば時系列など）を、２組の変数、すなわち環境変数およびシステム応答変数に分離およびグループ分けしてもよい。一実施形態では、変数セパレータは、観察された変数と制御された変数との間の相関関係を決定することによって、２組の変数を分離する。一実施形態では、制御変数１２０とシステム応答変数１１５とを互いに関連付ける変数セパレータ５１０によって、ノンパラメトリック推定が構築される。この関連付けは、事前発見プロセスによって、または領域知識を介して、制御変数１２０の依存構造を考慮に入れる。
【００９２】
少なくとも１つの記憶装置５６０を利用して、たとえば、観察された変数１０５、制御変数１２０、環境変数１３５、システム応答変数１２５、ＫＰＩ値１４０などを記憶してもよい。
【００９３】
一実施形態では、第１の関数推定器５２０が、制御変数とシステム応答変数との間の関係を回復する。一実施形態では、この関係は、関数
【００９４】
【数２７】
0006718500.tif 000028
として表され、ここで、
【００９５】
【数２８】
0006718500.tif 000029
はシステム変数であり、
【００９６】
【数２９】
0006718500.tif 000030
は制御変数Ｃの関数である。
【００９７】
一実施形態では、第２の関数推定器５３０が、大域の入力／出力マッピング関数
【００９８】
【数３０】
0006718500.tif 000031
を回復し、第１の関数推定器５２０において回復された関数を考慮に入れ、ここで、
【００９９】
【数３１】
0006718500.tif 000032
はＫＰＩ値であり、Ｅは環境変数である。
【０１００】
一実施形態では、オンライン推定器５４０は、環境変数の範囲を推定する。環境変数の範囲は、たとえば、履歴記録、予測アルゴリズム、領域知識などに基づくことができる。
【０１０１】
一実施形態では、オンラインオプティマイザ５５０は、第２の関数推定器５３０において推定される関数を最大化することによって、最適制御変数
【０１０２】
【数３２】
0006718500.tif 000033
を特定し、更にオンライン推定器５４０の環境変数の値を特定する。
【０１０３】
前述の説明は、あらゆる点で説明的かつ例示的なものであるが、限定的なものではないと理解すべきであり、本明細書に開示された本発明の範囲は、詳細な説明から判断すべきではなく、むしろ特許請求の範囲から、特許法によって許された範囲全体に従って解釈されるものと判断すべきである。図に示し、本明細書に記載した実施形態は、本発明の原理を例示するものに過ぎず、本発明の範囲および精神から逸脱することなく、様々な改善を当業者が実施してもよいことを理解されたい。本発明の範囲および精神から逸脱することなく、当業者は、他の様々な特徴の組合せを実施することもできる。特許法によって要求される詳細および独自性とともに本発明の態様をこうして説明してきたので、特許証によって保護された、特許請求されかつ望まれているものは、添付の特許請求の範囲に記載されている。
【図１】
0006718500.tif 000034
【図２】
0006718500.tif 000035
【図３】
0006718500.tif 000036
【図４】
0006718500.tif 000037
【図５】
0006718500.tif 000038