(57)【要約】  プロセスプラントにおけるプロセスデータの監視のためのシステムおよび方法であって、多変量プロセス変数データはそれぞれ複数の観測値を有する複数のプロセス変数を含み、各プロセス変数を観測値成分のセットを有するプロセス変数のベクトルとして定義することを含み、観測値成分のセットは、プロセス変数の観測値に対応する時間依存のプロセスデータを含んでいる。また、このシステムおよび方法は、多変量変換をそれぞれプロセス変数の一つに対応する複数のプロセス変数のベクトル変換の関数として計算される。各プロセス変数のベクトル変換は、プロセス変数を単一化する単変量変数の関数であり、多変量変換に基づいてプロセスの動作が表現され、このプロセスの動作の表現は、プロセス変数のそれぞれに対する単変量変数によるプロセスデータの多変量射影を指定する。【選択図】  なし【特許請求の範囲】
【請求項１】
プロセスプラントにおける多変量プロセスデータの監視を行う方法であって、前記多変量プロセスデータがそれぞれ複数の観測値を有する複数のプロセス変数を有しており、該方法は、
前記プロセス変数の観測値に対応する時間依存のプロセスデータを有する観測値成分のセットを有するプロセス変数ベクトルとして各プロセス変数を定義するステップと、
多変量変換をそれぞれプロセス変数の一つに対応する複数のプロセス変数のベクトル変換の関数として計算するステップであって、前記プロセス変数のベクトル変換のそれぞれは、前記プロセス変数を単一化する単変量変数の関数である、ステップと、
前記プロセスの動作を前記多変量変換に基づいて表現するステップであって、前記プロセスの動作の表現は、前記プロセス変数のそれぞれに対する前記単変量変数によって前記プロセスデータの多変量射影を指定する、ステップと
を包含する方法。
【請求項２】
請求項１記載の方法であって、前記複数のプロセス変数のそれぞれに対する観測値のセットを転置することをさらに包含し、プロセス変数ベクトルとして各プロセス変数を定義するステップは、前記観測値の転置されたセットを有するプロセス変数ベクトルとして各プロセス変数を定義することを包含している方法。
【請求項３】
請求項１記載の方法であって、各プロセス変数のベクトル変換を前記プロセス変数を単一化する単変量変数の関数として計算することをさらに包含している方法。
【請求項４】
請求項３記載の方法であって、各プロセス変数のベクトル変換を計算することは、各プロセス変数のベクトル変換を前記プロセス変数を単一化する単変量変数の関数として、および前記プロセス変数ベクトルに対して一意的な位置パラメータの関数として計算することを包含している方法。
【請求項５】
請求項３記載の方法であって、各プロセス変数のベクトル変換を計算することは、各プロセス変数のベクトル変換を前記プロセス変数を単一化する単変量変数の関数として、および前記観測値成分のセット内の各観測値成分に適用されるスケーリングパラメータの関数として計算することを包含している方法。
【請求項６】
請求項１記載の方法であって、前記複数のプロセス変数のベクトル変換の関数としての多変量変換の計算は、多変量変換を前記複数のプロセス変数のベクトル変換の積として計算することを包含している方法。
【請求項７】
請求項１記載の方法であって、前記多変量変換に基づいて前記プロセスの動作を表現することは、前記多変量変換に適用されたスケーリングパラメータに基づいてプロセスの動作を表現することを包含している方法。
【請求項８】
請求項１記載の方法であって、前記観測値成分のセット内の各観測値成分は、ベクトルの次元を有している方法。
【請求項９】
請求項８記載の方法であって、
所定の期間にわたってプロセスプラント内のプロセスからプロセスデータを収集すること、および
各ベクトルの次元を、前記観測値のサンプリングレートの関数として、およびデータ収集の時間期間の関数として決定すること
をさらに包含している方法。
【請求項１０】
請求項１記載の方法であって、前記単変量変数は、前記プロセス変数を束ねるように指定された変数を包含している方法。
【請求項１１】
請求項１記載の方法であって、前記多変量変換に基づいて前記プロセスの動作を表現することは、前記多変量変換の視覚化を前記単変量変数の関数として生成することを包含している方法。
【請求項１２】
請求項１１記載の方法であって、前記多変量変換の前記単変量変数の関数として視覚化は、前記プロセス変数ベクトルについて一意的な位置パラメータによって決定される単変量変数に対して一意的な位置を有する各プロセス変数を表示する視覚化を生成することを包含している方法。
【請求項１３】
請求項１記載の方法であって、各プロセス変数は、既知のプロセス変数制御境界を有し、前記方法は、各プロセス変数に一意的なスケーリングパラメータを使用して各プロセス変数のベクトル変換を同一のプロセス変数制御境界に正規化して前記プロセス変数に共通のプロセス変数の制御境界を生成することをさらに包含し、前記多変量変換に基づいて前記プロセスの動作を表現することは、前記共通のプロセス変数の境界に関連して各プロセス変数を監視することを包含している方法。
【請求項１４】
請求項１記載の方法であって、少なくとも一つのプロセス変数が未知のプロセス変数制御境界を含み、前記方法が、
前記少なくとも一つのプロセス変数に対するデータのセットを収集するステップであって、前記収集されたプロセスデータのセットは、前記プロセスが正常に動作しているときに前記プロセスの正常動作を表し、前記少なくとも一つの多変量プロセスデータの複数の観測値を含む、ステップと、
前記プロセスの正常動作の下に未知のプロセス変数制御境界を決定するステップと、
各プロセス変数に一意的なスケーリングパラメータを使用して、同一のプロセス変数制御境界に各プロセス変数のベクトル変換を正規化して、前記プロセス変数に対する共通のプロセス変数制御境界を生成するステップと
をさらに包含し、
前記多変量変換に基づいて前記プロセスの動作を表現することは、前記決定されたプロセス変数の境界に関連して各プロセス変数を監視することを包含している方法。
【請求項１５】
請求項１記載の方法であって、前記プロセスのリアルタイム動作の測度を有するプロセスプラント内のプロセスからプロセスデータを収集することをさらに包含している方法。
【請求項１６】
請求項１５記載の方法であって、前記プロセスプラント内のプロセスからプロセスデータを収集することは、前記プロセスがオンラインのときにプロセスのリアルタイム動作の測度を有するプロセスプラント内のプロセスからの監視されたオンラインプロセスデータを収集することを包含している方法。
【請求項１７】
請求項１５記載の方法であって、前記多変量変換に基づいて前記プロセスの動作を表現することは、前記多変量変換の視覚化を前記単変量変数のリアルタイムの関数として生成することを包含している方法。
【請求項１８】
プロセスプラントにおけるプロセスの監視を行うためのシステムであって、前記プロセスがそれぞれ複数の観測値を有する複数のプロセス変数を有する多変量プロセスデータを生成し、各プロセス変数は、前記プロセス変数の観測値に対応する時間依存のプロセスデータを有する観測値成分のセットを有するプロセス変数のベクトルとして定義され、前記システムは、
多変量変換を、それぞれ前記プロセス変数の一つに対応する複数のプロセス変数のベクトル変換の関数として計算するように構成されたプロセス変数ツールであって、各プロセス変数のベクトル変換は前記プロセス変数を単一化する単変量変数の関数である、プロセス変数ツールと、
前記多変量変換に基づいて前記プロセスの動作を表現するように構成された監視ツールであって、前記プロセスの動作の表現は、各プロセス変数に対する前記単変量変数によって前記プロセスデータの多変量射影を指定する、監視ツールと
を備えているシステム。
【請求項１９】
請求項１８記載のシステムであって、各プロセス変数は、前記複数のプロセス変数のそれぞれに対する観測値の転置されたセットを有するプロセス変数ベクトルとして定義されるシステム。
【請求項２０】
請求項１８記載のシステムであって、前記プロセス変数ツールは、各プロセス変数のベクトルの変換を前記プロセス変数を単一化する単変量変数の関数として計算するように構成されているシステム。
【請求項２１】
請求項２０記載のシステムであって、前記プロセス変数ツールは、各プロセス変数のベクトルの変換を前記プロセス変数を単一化する単変量変数の関数として、および前記プロセス変数ベクトルに対して一意的な位置パラメータの関数として計算するように構成されているシステム。
【請求項２２】
請求項２０記載のシステムであって、前記プロセス変数ツールは、各プロセス変数のベクトルの変換を前記プロセス変数を単一化する単変量変数の関数として、および前記観測値成分のセット内の各観測値成分に適用されるスケーリングパラメータの関数として計算するように構成されているシステム。
【請求項２３】
請求項１８記載のシステムであって、前記プロセス変数ツールは、前記多変量変換を前記複数のプロセス変数のベクトル変換の積として計算するように構成されているシステム。
【請求項２４】
請求項１８記載のシステムであって、前記プロセス変数ツールは、前記多変量変換に適用されたスケーリングパラメータに基づいてプロセスの動作を表現するように構成されているシステム。
【請求項２５】
請求項１８記載のシステムであって、前記観測値成分のセット内の各観測値成分は、ベクトルの次元を有しているシステム。
【請求項２６】
請求項２５記載のシステムであって、
所定の期間にわたってプロセスプラント内のプロセスからプロセスデータを収集するように構成されたデータ収集ツールをさらに備え、各ベクトルの次元が、前記観測値のサンプリングレートの関数として、およびデータ収集の時間期間の関数として決定されるシステム。
【請求項２７】
請求項１８記載のシステムであって、前記単変量変数は、前記プロセス変数を束ねるように指定された変数を包含しているシステム。
【請求項２８】
請求項１８記載のシステムであって、前記監視ツールは、前記多変量変換の視覚化を前記単変量変数の関数として生成するように構成されているシステム。
【請求項２９】
請求項２８記載のシステムであって、前記監視ツールは、前記プロセス変数ベクトルについて一意的な位置パラメータによって決定される単変量変数に対して一意的な位置を有する各プロセス変数を表示する視覚化を生成するように構成されているシステム。
【請求項３０】
請求項１８記載のシステムであって、各プロセス変数は、既知のプロセス変数制御境界を有し、前記システムは、各プロセス変数に一意的なスケーリングパラメータを使用して各プロセス変数のベクトル変換を同一のプロセス変数制御境界に正規化して前記プロセス変数に対する共通のプロセス変数の制御境界を生成するように構成された正規化ツールをさらに有し、前記監視ツールは、前記共通のプロセス変数の境界に関連して各プロセス変数を監視するように構成されているシステム。
【請求項３１】
請求項１８記載のシステムであって、少なくとも一つのプロセス変数が未知のプロセス変数制御境界を含み、前記システムが、
前記少なくとも一つのプロセス変数に対するデータのセットを収集するように構成されたデータ収集ツールであって、前記収集されたプロセスデータのセットは、前記プロセスが正常に動作しているときに前記プロセスの正常動作を表し、前記少なくとも一つの多変量プロセスデータの複数の観測値を含む、データ収集ツールと、
前記プロセスの正常動作の下に未知のプロセス変数制御境界を決定するように構成された解析ツールと、
各プロセス変数に一意的なスケーリングパラメータを使用して、同一のプロセス変数制御境界に各プロセス変数のベクトル変換を正規化して、前記プロセス変数に対する共通のプロセス変数制御境界を生成する正規化ツールと
をさらに包含し、
前記監視ツールは、前記共通のプロセス変数の境界に関連して各プロセス変数を監視するように構成されているシステム。
【請求項３２】
請求項１８記載のシステムであって、前記プロセスのリアルタイム動作の測度を有するプロセスプラント内のプロセスからプロセスデータを収集するように構成されたデータ収集ツールをさらに備えているシステム。
【請求項３３】
請求項３２記載のシステムであって、前記データ収集ツールは、前記プロセスがオンラインのときにプロセスのリアルタイム動作の測度を有するプロセスプラント内のプロセスからの監視されたオンラインプロセスデータを収集するように構成されているシステム。
【請求項３４】
請求項３２記載のシステムであって、前記監視ツールは、前記多変量変換の視覚化を前記単変量変数のリアルタイムの関数として生成するように構成されているシステム。
【請求項３５】
プロセスプラントにおけるプロセス監視を行うためのシステムであって、該システムは、
多変量プロセスデータを単一化する単変量変数に基づく複数の観測値をそれぞれ有する複数の多変量プロセス変数のそれぞれに対して単変量統計表現を生成するように構成された第１の解析ツールと、
前記複数の多変量プロセス変数の複数の単変量統計表現に基づいて前記プロセスの動作の単変量統計表現を生成するように構成された第２の解析ツールと、
前記プロセスの動作の単変量統計表現に基づいて、単変量変数の関数として前記プロセスの動作を監視するように構成された第３の解析ツールと
を備えたシステム。
【請求項３６】
請求項３５記載のシステムであって、各プロセス変数を観測値成分のセットを有するプロセス変数のベクトルとして定義するように構成された第４の解析ツールをさらに備え、前記観測値成分のセットは、多変量プロセス変数の観測値に対応する時間依存のプロセスデータを有しているシステム。
【請求項３７】
請求項３５記載のシステムであって、前記複数の多変量プロセス変数の複数の単変量統計表現の組合せに基づいて、前記複数の多変量プロセス変数を前記複数の多変量プロセス変数の単変量統計表現に変換するように構成された第４の解析ツールをさらに備えたシステム。
【請求項３８】
請求項３７記載のシステムであって、前記第２の解析ツールは、前記複数の多変量プロセス変数の単変量表現に基づいて前記プロセスの動作の単変量統計表現を生成するように構成されているシステム。
【請求項３９】
請求項３５記載のシステムであって、前記第３の解析ツールは、前記単変量変数の関数としての前記プロセスの動作の単変量統計表現に基づいて前記プロセスのリアルタイムの動作を監視するように構成されているシステム。
【請求項４０】
請求項３５記載のシステムであって、前記第３の解析ツールは、前記単変量変数の関数としての前記プロセスの動作の単変量統計表現の視覚化を生成するように構成されているシステム。
【請求項４１】
請求項４０記載のシステムであって、前記第３の解析ツールは、前記単変量変数に対して一意的な位置を有する各プロセス変数を表示する視覚化を生成するように構成されているシステム。
【請求項４２】
請求項３５記載のシステムであって、各多変量プロセス変数は既知のプロセス変数制御境界を有し、前記システムは、多変量プロセス変数の各単変量統計表現を共通のプロセス変数の制御境界に正規化するように構成された第４の解析ツールをさらに備え、前記第３の解析ツールは、前記共通のプロセス変数制御境界に関連して前記プロセスの動作の単変量統計表現によって表現されるように各多変量プロセス変数を監視するように構成されているシステム。
【請求項４３】
請求項３５記載のシステムであって、前記多変量プロセス変数の少なくとも一つが未知のプロセス変数制御境界を有し、前記システムは、
前記少なくとも一つの多変量プロセス変数に対するデータのセットを収集するように構成された第４の解析ツールであって、前記収集された多変量プロセス変数データのセットは、前記プロセスが正常に動作しているときに前記プロセスの正常動作を表し、前記少なくとも一つの多変量プロセス変数の複数の観測値を有している、第４の解析ツールと、
前記プロセスの正常動作の下に前記未知のプロセス変数制御境界を決定するように構成された第５の解析ツールと、
多変量プロセス変数の各単変量統計表現を共通のプロセス変数の制御境界に正規化するように構成された第６の解析ツールと
をさらに備え、
前記第３の解析ツールは、前記共通のプロセス変数の制御境界に関連して前記プロセスの動作の単変量統計表現によって表されるように各多変量プロセス変数を監視するように構成されているシステム。
【請求項４４】
請求項３５記載のシステムであって、複数の多変量プロセス変数のそれぞれに対する単変量統計表現は、前記第１の解析ツールによって、以下の式のようにして計算されるシステム。
ｐi＝［（ｙ−αi）2＋ｋiｘi,j2］
ここで、ｙ＝単変量変数、
αi＝各多変量プロセス変数のプロセス変数のベクトルに一意的な位置パラメータ
ｋi＝スケーリングパラメータ
Ｘij＝プロセス変数のベクトルの観測値成分のセット内の観測値成分
ｉ＝１，...，ｎ、ここで、ｎは、多変量プロセス変数の数に相当、および
ｊ＝１，...，ｍ、ここで、ｍは、観測値の収集におけるサンプリングキーボードと観測値の収集におけるトータル時間とによって決定される各観測値成分に対するベクトルの次元に対応している。
【請求項４５】
請求項４４記載のシステムであって、多変量プロセス変数のそれぞれのプロセス変数のベクトルが、前記第４の解析ツールによって、以下の式のようにして計算されるシステム。
Ｘi＝［ｘi,1，ｘi,2，．．．，ｘi,m］T
ここで、Ｔは転置演算を含む。
【請求項４６】
請求項３５記載のシステムであって、前記プロセスの動作の単変量統計表現は、前記第２の解析ツールによって、以下の式のようにして計算されるシステム。
Ｍv（ｙ）＝Ｋp／Ｐ（ｙ）
ここで、Ｋp＝スケーリングファクター
Ｐ（ｙ）＝ｐ1（ｙ）＊ｐ2（ｙ）．．＊ｐn（ｙ）
ｙ＝単変量変数
ｐn（ｙ）＝ｎ番目の多変量プロセス変数に対する単変量統計表現、および
ｎ＝多変量プロセス変数の数である。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、米国仮出願第６０／８２７６３１号に基づく優先権を主張し、その内容は、参照によってここに取り入れられる。
【０００２】
本開示は、概略的にはプロセスプラントの診断に関し、より詳細には、多変量統計技術を含む監視及び診断システムに関する。
【背景技術】
【０００３】
化学、石油又は他のプロセスなどのプロセス制御システムは、典型的には、少なくとも一つのホスト又はオペレータワークステーションと、例えばフィールドデバイス等の一つ又はそれ以上のプロセス制御および計装機器に、アナログ、デジタル又はアナログ／デジタル結合バスを介して通信可能に接続された一つ又はそれ以上の集中型又は非集中型のプロセスコントローラを有している。例えば、バルブ、バルブポジショナ、スイッチ、トランスミッタおよびセンサ（即ち、温度、圧力、流速センサ）等のフィールドデバイスは、プロセスプラント環境に配置され、バルブの開閉、プロセスパラメータの測定、流量の増減などのプロセス内の機能を実行する。周知のＦＵＮＤＡＴＩＯＮ（登録商標）Ｆｉｅｌｄｂｕｓ（以下「Ｆｉｅｌｄｂｕｓ」と称する。）プロトコル又はＨＡＲＴ（登録商標）プロトコルに従うスマートフィールドデバイスは、制御計算、警告機能やプロセスコントローラ内で共通的に実現されるたの制御機能も果たし得る。
【０００４】
典型的には、プロセスプラント環境内に位置するプロセスコントローラは、フィールドデバイス及び／又はフィールドデバイスに関連する多の情報によって生成され又はこれに関連するプロセス測定値又はプロセス変数を表す信号を受信する。コントローラアプリケーションは、例えば、プロセス制御の決定、受信した情報に基づく制御信号の生成、およびＨＡＲＴおよびＦｉｅｌｄｂｕｓフィールドデバイス等のフィールドデバイスにおいて行われている制御モジュール又はブロックの調整を行う異なる制御モジュールを実現する。プロセスコントローラ内の制御モジュールは、通信ライン又は信号パス上にフィールドデバイスに向けて制御信号を送出する。
【０００５】
フィールドデバイスおよびプロセスコントローラからの情報は、典型的には、オペレータワークステーション、保守ワークステーション、パーソナルコンピュータ、ハンドヘルドデバイス、データヒストリアン、レポート生成、集中型データベースなどの一つ又はそれ以上の他のハードウエアデバイスに利用可能にし、オペレータ又は保守要員が、例えばプロセス制御ルーチンの設定を変更し、プロセスコントローラ又はスマートフィールドデバイス内の制御モジュールの動作を変更し、プロセス又はプロセスプラント内の特定のデバイスの状態を閲覧し、人の訓練又はプロセス制御ソフトウエアのテストの目的でプロセスの動作をシミュレートし、およびプロセスプラント内の問題又はハードウエア欠陥を診断するのを可能にする。
【０００６】
知られているように、プロセスプラント環境、特に多数のフィールドデバイスおよび支援装置を有するプロセスプラント内では、頻繁に問題が生ずる。これらの問題は、壊れた又は欠陥のあるデバイス、不正なモードでの存在、不適切に調整されたプロセス制御ループ、プロセスプラント内のデバイス間の通信における一つ又はそれ以上の欠陥などの形態を採る。これらおよび他の問題は、実際は非常に多く、一般的には、異常な状態でのプロセスの運転という結果（即ち、異常な状態にあるプロセスプラント）をもたらし、通常これは、プロセスプラントの準最適性能に関連している。
【０００７】
多くの診断ツールおよびアプリケーションは、プロセスプラント内の問題の原因を検出し決定し、一旦問題が発生しおよび検出されたならオペレータ又は保守要員がその問題を診断し修正するように開発されている。例えば、直接又はワイヤレスバス、イーサネット（Ｒ）、モデム、電話回線などの通信接続を介してプロセスコントローラに典型的に接続されるオペレータワークステーションは、エマーソンプロセスマネジメント社によって販売されているＤｅｌｔａＶ（登録商標）およびＯｖａｔｉｏｎ（登録商標）などのソフトウエアを実行するようにされたプロセッサおよびメモリを有している。これらの制御システムは、多くの制御モジュールおよび制御ループ診断ツールを有している。同様に、コントローラアプリケーションと同様の通信接続を介して又はプロセス制御（ＯＰＣ）のためのオブジェクトリンキングアンドエンベッディング（ＯＬＥ）やハンドヘルド接続のような異なる通信接続を介して、フィールドデバイスなどのプロセス制御デバイスに接続され得る保守ワークステーションは、典型的には、プロセスプラント内のフィールドデバイスによって生成される保守アラーム又は警報を閲覧し、プロセスプラント内のデバイスをテストし、プロセスプラント内のフィールドデバイスおよび他のデバイス上で保守アクティビティを実行するように設計されている。同様の診断アプリケーションが、プロセスプラント内のサポート装置内の診断の問題に向けて開発されている。
【０００８】
このように、ＡＭＳ（登録商標）一式の商業的に入手が可能なソフトウエア：Intelligent Device Manager from Emerson Process Managementは、フィールドデバイスの動作状態を究明し追跡するために、フィールドデバイスと通信しこれに関連するデータを格納する。「フィールドデバイス管理システムにおいて使用する統合された通信ネットワーク」と題された米国特許第５，９６０，２１４号も参照されたい。いくつかの例では、ＡＭＳソフトウエアは、フィールドデバイス内のパラメータを変更し、例えば自己較正ルーチン又は自己診断ルーチンなどのアプリケーションをフィールドデバイスにそれ自身上で実行させ、フィールドデバイスの状態又は健全性についての情報を得るために使用される。この情報は、例えば、状態情報（例えば、アラーム又は他の同様なイベントが起こったかどうか）、デバイス構成情報（例えば、フィールドデバイスが現在又は後に構成される方式およびフィールドデバイスによって使用される測定ユニットのタイプ）、デバイスパラメータ（例えば、フィールドデバイスの範囲値及び他のパラメータ）などを含み得る。もちろん、この情報は、保守要員によってフィールドデバイスの問題を監視し、保守し及び／又は診断するのに使用され得る。
【０００９】
同様に、多くのプロセスプラントは、例えば、種々の回転する装置の動作状態のモニター、診断および最適化するためのＣＳＩシステムズによって供給されているＲＢＭアプリケーションのようなソフトウエアアプリケーションを有している。保守要員は、通常、プラント内で回転する装置の性能を保守および監視し、その回転装置の修理又は交換の時期および必要性を決定するためにこれらのアプリケーションを使用する。同様に、多くのプロセスプラントは、発電および配電の制御と維持とを行うために、Liebert and ASCO社によって提供されるような電源制御および診断アプリケーションを含み得る。例えば、リアルタイムオプティマイザ（ＲＴＯ＋）のような制御オプティマイザのような制御最適化アプリケーションを、プロセスプラントの制御アクティビティを最適化するために実行することもよく知られている。このような最適化アプリケーションは、例えば利益のようないくつかの所望の最適化変数に関連してプロセスプラントの運転を最適化するために入力がどのように変更され得るかを予測するために、典型的には、プロセスプラントの複雑なアルゴリズム及び／又はモデルを使用する。
【００１０】
これらおよび他の診断および最適化アプリケーションは、一つ又はそれ以上のオペレータ又はワークステーションにおけるシステムワイドな基準で典型的に実現され、プロセスプラントの動作状態に関してオペレータ又は保守要員に予め構成されたディスプレイ又はプロセスプラント内のデバイス若しくは装置を提供する。典型的な表示には、プロセスプラントなどの内部のプロセスコントローラ又は他のデバイスによって生成されたアラームを受け取る警告表示、プロセスプラント内のプロセスコントローラ又は他のデバイスの動作状態を示す制御表示、プロセスプラント内のデバイスの動作状態を示す保守表示などが含まれる。同様に、これらおよび他の診断アプリケーションは、オペレータ又は保守要員が制御ループを再調整し又は他の制御パラメータをリセットするのを可能とし、それらのフィールドデバイスの現在の状態を決定するために一つ又はそれ以上のフィールドデバイス上でテストするのを可能とし、又はフィールドデバイス又は他の装置を較正するのを可能とする。
【００１１】
これらの種々のアプリケーションおよびツールは、プロセスプラント内の問題を同定しおよび修正するのに非常に役に立つが、これらの診断アプリケーションは一般的にプロセスプラント内に問題が既に出現した後にのみ、従って、異常な状態がプロセスプラント内に既に存在した後に使用されるように構成されている。不運にも、異常な状態はそれが発見され同定されこれらのツールを使用して修正される前のしばらくの間存在し、結果的に問題が発見され同定され修正される間、プロセスプラントの準最適性能がもたらされる。多くの場合、制御オペレータは、アラーム、警報又はプロセスプラントの低い性能に基づいていくつかの問題が存在することを最初に検出するであろう。オペレータは、次に、可能性のある問題を保守要員に通知するであろう。保守要員は、実際の問題を検出し又は検出せず、運転試験又は他の診断アプリケーションをさらに促し、又は実際の問題を同定するのに必要な他のアクティビティを実行することが必要である。一旦問題が同定されると、保守要員は部品を注文し修理のスケジュールを立て、これらの全ては、問題の出現とその問題の修正との間の時間の重大な期間をもたらし、その間、プロセスプラントは一般的にプラントの準最適動作に関連する異常な状態で動作する。
【００１２】
加えて、多くのプロセスプラントは、比較的短い時間でプロセスプラント内に多大なコスト又はダメージを結果としてもたらす異常な状態を経験することができる。例えば、いくつかの異常状態は、これらの異常状態がたとえ短時間であっても、装置への重大なダメージ、原料の損失、又はプロセスプラント内の予期しない休止時間を生じる。従って、問題が出現した後のプラント内の問題の単なる検出は、どんなに素早く問題修正が行われても、重大なロス又はダメージをプロセスプラント内に結果としてもたらす。結果的に、異常状態が初めの段階で発生するのを妨げようとするのは、それに対して単に反応し異常状態が生じた後にプロセスプラント内の問題を修正しようと試みることに代えて好ましいことである。
【００１３】
プロセスプラント内の何れかの重大なロスが発生する前に予想された異常状態を妨げるためのステップを採るという目的で、これらの異常状態が実際に生じる前にユーザがプロセスプラント内の或る異常状態の出現を予想するのを可能とするデータを収集するのに使用され得る一つの技術がある。この手法は、「Root Cause Diagnostics」と題された米国特許出願第０９／９７２，０７８号、米国特許出願第７，０８５，６１０号（米国特許出願第０８／６２３，５６９号、米国特許出願第６，０１７，１４３号を基礎とする）に開示されている。これらの出願の両方の全体が参照によってここに組み込まれる。一般的に言って、この技術は、プロセスプラント内のフィールドデバイスのような多くのデバイスのそれぞれに、統計的データ収集および処理ブロック又は統計的処理モニタリング（ＳＰＭ）ブロックを置いている。統計的データ収集および処理ブロックは、プロセス変数データを収集し、平均、メジアン、標準偏差などのその収集されたデータに関連する或る統計的な測度を決定する。これらの統計的な測度は、次にユーザに送られ、既知の異常状態の将来的な出現を提示するパターンを認識するように分析される。特定の疑わしい将来の異常状態が検出されると、内在する問題を修正するステップが採られ、それによって初期の段階で異常状態が避けられる。
【００１４】
主成分分析（ＰＣＡ）は、プロセス制御産業で多次元データセットの解析について使用されている多変量解析技術である。ＰＣＡ技術は、一般には、高次元空間から低次元空間へのデータの減少を含み、これはオリジナルのデータにおける重要な変動の殆どを説明する。例えば、ＰＣＡは、多変量データ空間の次元（即ち、多次元データ）を小さい次元（即ち、２又は３次元）に減少させて、多次元データが低次元の主要成分空間に射影される（スコアされる）ように使用され得る。産業プロセスにおける欠陥検出のためのＰＣＡ技術の典型的な実現に関する更なる詳細は、Ｌ．Ｈ．Ｃｈｉａｎｇ他の「産業システムおける欠陥の検出と診断」（Springer-Verlag London Limited, pp. 35-54 (2001)）と、Ｅ．Ｌ．Ｒｕｓｓｅｌｌ他の「化学プロセスにおける欠陥検出および診断のためのデータ駆動技術」（pringer-Verlag London Limited, (2000)）とにおいて見ることができる。
【００１５】
多くのソフトウエアパッケージがＰＣＡ機能を提供しているが、しかし、これらのパッケージの殆どはデータベースからのオフラインプロセスデータを使用している。換言すれば、ＰＣＡはオフライン多変量統計解析ツールを有していた。このことは、履歴データ又は格納された多次元データのようなプロセスデータにおける異常状態の事後分析には有用であるが、それはリアルタイム（例えば、オンライン、アットライン又はインラインのデータ分析）の異常状態の検出に使用することはできない。いくつかの現存するＰＣＡソフトウエアパッケージは、データのリアルタイム分析が可能であるかも知れないが、しかし、そのソフトウエアがしばしばＯＰＣサーバーなどのインターフェイスを介して制御システムからのプロセスデータにアクセスするであろう場合に限られる。その結果、プラントの人員は、ソフトウエアパッケージの制御システムへの及び制御システムとのインターフェイスを含めて、制御システムとは別のソフトウエアパッケージの保守の負担を不幸にも負うであろう。更に、異常状態の検出は、異常状態阻止ツールが種々の異常状態に対して油断がなく正確であることを示すためには、全ての時間と種々の異常状態（例えば、連続プロセス対一回だけの解）にとって有用であるべきである。
【００１６】
また、従来のＰＣＡ技術は、プロセスの動力学を表す情報を容易に搬送できない。ＰＣＡ分析の結果、即ち、スコアは、内在するプロセスデータの統計的表現である散乱したプロットによりしばしば表示される。スコアデータが散乱したプロットに記録され又は説明されるような方法にも拘わらず、不正確に特徴化されたデータに関する小さな不確実性が一般的に存在する。以下に記述するように、更なる情報無しに、プロセスオペレータは、与えられたスコアが正常か異常かをしばしば確信することができない。例えば、ＰＣＡ分析の最初の２つの主成分が多次元データにおける最大の変動を捕らえ、それらはプロセスにおける最大の変動を含む多次元データの挙動についての最も重大な情報を表しているので、ＰＣＡ分析からのその最初の２つのスコアがプロットされ得る。
【００１７】
これらの理由から、従来のＰＣＡ技術は、異常状態の検出においては頼りにならない性能しか示さない。ＰＣＡスコアは、人が存在しないとき（即ち、誤認警報）に異常状態を提示することができ、反対に、プロセスが異常に振る舞う（即ち、誤った診断）ときに正常運転を提示し得る。更に、ＰＣＡを実現するために、ＰＣＡモデルがトレーニングデータ（例えば、「正常な」プロセスの動作からのデータ）を使用して開発され、結果として生まれるモデルは異常なプロセス状態を検出するために適用される。トレーニングデータは、主成分に関してモデルを引き出すのに使用される。そのようなものとして、ＰＣＡスコアは実際のプロセス変数に直接関係しているものではなく、ＰＣＡを使用して監視されている実際のプロセス変数の物理的プロセスの制限は分析において失われ、それによって統計的解釈（例えば、「データの良好性」の定量化）を要求する。
【発明の概要】
【課題を解決するための手段】
【００１８】
本開示の特定の特徴によれば、プロセス制御システムおよびその何れかの要素の実装の監視および診断を容易にするための多くの技術が開示されている。このような監視および診断は、一般的に、主成分分析（ＰＣＡ）および主成分回帰、部分最小二乗法、Ｆｉｓｈｅｒ判別分析（ＦＤＡ）又は正準変量分析などの他の多変量統計解析技術をプロセス制御システムに統合し又は包含している。次に、監視の多変量的性質は、診断、異常状態防止および他の故障検出をサポートする。また、このような監視および診断は、一般的に、多変量プロセスデータを単独の監視変数を用いて監視する単変量解析を統合し又は包含している。多変量又は単変量解析のそれぞれは、異常動作の検出および防止のような多くの故障監視および検出技術を容易にし得る。
【００１９】
特に、多変量プロセス変数を束ねて一つにするように指定された、単一の監視変数で多変量プロセスデータを監視するための方法が開示されている。多変量プロセスデータはそれぞれ複数の観測値を有する複数のプロセス変数を含んでいる。各プロセス変数は、時間依存のデータを含むプロセス変数ベクトルとして定義される。各プロセス変数ベクトルは、単変量変数、プロセス変数ベクトルに対する一意的な位置インディケータ、およびスケーリングファクターの関数として変換される。単変量変数は種々のプロセス変数を単一化する。多変量変換は、種々のプロセス変数のベクトル変換によって計算され、多変量変換は多変量射影を指定し、多変量データは単変量変数によって射影される。
【００２０】
多変量変換は単変量変数の関数として動作し、プロセス変数のそれぞれは、同じスケール上で監視され得る。プロセス変数のそれぞれは、単変量解析に従って記述されたプロセス変数のそれぞれに共通のプロセスの境界を決定するために、単変量解析（例えば、同じ応答に対して調整された）によって正規化される。従って、プロセスのモデルは多変量変換に基づいており、モデルは、プロセス変数のそれぞれに対する単変量変数によって、プロセスの多変量射影を指定する。プロセス変数は、範囲外（異常）イベントに遭遇したプロセス変数を特定するために、同じスケール上同時に監視され得る。
【００２１】
単変量方法は、訓練データを使用する多変量技術に比較して、もし全てのプロセスの境界が既知であるなら訓練データを使用することなく完了するが、もし境界の全てが既知ではないならプロセス変数の境界を決定するための方法も開示される。もしプロセス変数の境界が未知なら、正常なプロセスの動作状態の下で集められた訓練データが、プロセス変数を正規化しおよび共通のプロセス変数の境界を展開するのに使用され得る。
【図面の簡単な説明】
【００２２】
【図１】図１は、分散型プロセス制御システムと、一つ又はそれ以上のオペレータ及び保守ワークステーション、コントローラフィールドデバイス及び支援装置を含むネットワークとを有するプロセスプラントの典型的なブロックダイアグラムである。
【００２３】
【図２】図２は、図１のプロセスプラントの部分の典型的なブロックダイアグラムであり、プロセスプラントの異なる要素内に位置する異常状態防止システムの種々の構成要素間の通信相互接続を例示している。
【００２４】
【図３】図３は、プロセスプラントの監視動作に対する動作上の要求又は制御境界に伴うプロセス変数の図式プロットである。
【００２５】
【図４】図４は、動作上の要求又は制御境界に比較して多くのプロセス変数の図式の多変量視覚化又はプロットであり、視覚化又はプロットは、特定の実施形態及びオンラインプロセス監視についての開示の特徴に従って生成されるユーザインターフェイスに組み込むのに適している。
【００２６】
【図５】図５は、プロセス変数の一つが動作上の要求又は制御境界を超えた後の図４の視覚化プロットの他の図示であり、異常状態は故障状態を示している。
【００２７】
【図６】図６は、主成分分析（ＰＣＡ）を実現するように構成され、本開示の一つの特徴に従ったトレーニング及び運転のためのプロセス制御システム内にインスタンスが作成された典型的な機能ブロックを表している。
【００２８】
【図７】図７は、本開示の他の特徴に従ったＰＣＡ技術の実現の間を表すもう一つの典型的なＰＣＡ機能ブロックの表現である。
【００２９】
【図８】図８は、本開示のもう一つの特徴に従う多変量監視及び診断モジュール（ＭＭＤＭ）のブロックダイアグラムであり、いくつかの実施形態におけるその成分は、図６及び図７の機能ブロックにおいて実現され得る。
【００３０】
【図９】図９−１２は、本開示のもう一つの特徴に従って生成された例示のユーザインターフェイスディスプレイであり、これは、いくつかの実施形態において、多変量監視及び故障検出のために、図６及び図７の機能ブロックによって、又は図８のＭＭＤＭツールによって生成又は支援され得る。
【図１０】図９−１２は、本開示のもう一つの特徴に従って生成された例示のユーザインターフェイスディスプレイであり、これは、いくつかの実施形態において、多変量監視及び故障検出のために、図６及び図７の機能ブロックによって、又は図８のＭＭＤＭツールによって生成又は支援され得る。
【図１１】図９−１２は、本開示のもう一つの特徴に従って生成された例示のユーザインターフェイスディスプレイであり、これは、いくつかの実施形態において、多変量監視及び故障検出のために、図６及び図７の機能ブロックによって、又は図８のＭＭＤＭツールによって生成又は支援され得る。
【図１２】図９−１２は、本開示のもう一つの特徴に従って生成された例示のユーザインターフェイスディスプレイであり、これは、いくつかの実施形態において、多変量監視及び故障検出のために、図６及び図７の機能ブロックによって、又は図８のＭＭＤＭツールによって生成又は支援され得る。
【００３１】
【図１３】図１３−１４は、本開示のもう一つの特徴に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【図１４】図１３−１４は、本開示のもう一つの特徴に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【００３２】
【図１５】図１５−１９は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【図１６】図１５−１９は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【図１７】図１５−１９は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【図１８】図１５−１９は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【図１９】図１５−１９は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【００３３】
【図２０】図２０−２１は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【図２１】図２０−２１は、更にもう一つの実施形態に従って生成され、多変量監視及び故障検出のために図９−１２のユーザインターフェイスディスプレイと同様の方法で生成され又は支援される更なる例示のユーザインターフェイスディスプレイである。
【００３４】
【図２２】図２２−２３は、図６及び７の機能ブロックで使用される、本開示のもう一つの特徴に従うＰＣＡベースの定常状態検出システム及び技術のブロックダイアグラムである。
【図２３】図２２−２３は、図６及び７の機能ブロックで使用される、本開示のもう一つの特徴に従うＰＣＡベースの定常状態検出システム及び技術のブロックダイアグラムである。
【００３５】
【図２４】図２４は、加熱炉の図式ダイアグラムであり、本開示の技術及びＭＭＤＭツールが監視、診断及び故障検出のために適用される例示のプロセスである。
【００３６】
【図２５】図２５は、主成分ｔ１及びｔ２に基づいて２次元空間にマッピングされたトレーニングデータスコアのグラフ表示であり、トレーニングデータは、図２４の加熱炉に関連してＰＣＡベースの監視、診断及び故障検出に使用されている。
【００３７】
【図２６】図２６−２７は、図２４の加熱炉についてＣＡベースの監視、診断及び故障検出に関連してＭＭＤＭツールによって生成される例示のユーザインターフェイスディスプレイである。
【図２７】図２６−２７は、図２４の加熱炉についてＣＡベースの監視、診断及び故障検出に関連してＭＭＤＭツールによって生成される例示のユーザインターフェイスディスプレイである。
【００３８】
【図２８】図２８は、本開示のもう一つの実施形態に従う動態的作用の検出及び決定のための多変量統計に関連するデータ妥当性検証に使用される原油流量データの図式プロットである。
【００３９】
【図２９】図２９は、図２８に表される流量に関連する実際の出力温度データを、動態的作用の検出又は決定のための多変量統計技術のＰＣＲモデルによって予想される出力温度と比較した図式プロットである。
【００４０】
【図３０】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【図３１】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【図３２】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【図３３】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【図３４】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【図３５】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【図３６】図３０−３６は、動態的作用の検出又は決定のための多変量統計技術の実現に関連して生成される例示のユーザインターフェイスディスプレイである。
【００４１】
【図３７】図３７は、一実施形態に従う動態的作用の検出又は決定のための多変量統計技術のフローダイアグラムである。
【００４２】
【図３８】図３８−３９は、例示の加熱炉についての実際の出力温度を、本開示のもう一つの特徴に従う多変量統計モデルによって予想されるデータと比較した図式プロットである。
【図３９】図３８−３９は、例示の加熱炉についての実際の出力温度を、本開示のもう一つの特徴に従う多変量統計モデルによって予想されるデータと比較した図式プロットである。
【００４３】
【図４０】図４０−４１は、多変量統計モデルに基づくコークス化検出の一実施形態に従う実際と予想との出力温度データを比較する図式プロットである。
【図４１】図４０−４１は、多変量統計モデルに基づくコークス化検出の一実施形態に従う実際と予想との出力温度データを比較する図式プロットである。
【００４４】
【図４２】図４２は、本開示のもう一つの特徴に従うコークス化についてのもう一つの非回帰ベースの技術を実証する加熱炉の出力温度についてのＰＣＡモデルパラメータの図式プロットである。
【００４５】
【図４３】図４３は、図４２の図式プロットに表される非回帰ベースのコークス化検出技術の一実施形態のフローダイアグラムである。
【発明を実施するための形態】
【００４６】
図１を参照すれば、異常得状態防止システムが実現されている例示のプロセスプラント１０は、一つ又はそれ以上の通信ネットワークを介して指示装置に互いに接続された多くの制御および保守システムを含んでいる。特に、図１のプロセスプラント１０は、一つ又はそれ以上のプロセス制御システム１２および１４を含んでいる。プロセス制御システム１２は、ＰＲＯＶＯＸ若しくはＲＳ３システムのような従来のプロセス制御システム又は他のどのような制御システムでもあり得、コントローラ１２Ｂと入力／出力（Ｉ／Ｏ）カード１２Ｃとに接続されたオペレータインターフェイス１２Ａを含み、入力／出力カード１２Ｃは次にアナログおよびHighway Addressable Remote Transmitter（ＨＡＲＴ）フィールドデバイス１５のような種々のフィールドデバイスに接続されている。プロセス制御システム１４は、分散型プロセス制御システムであり得、イーサネット（Ｒ）のようなバスを解して一つ又はそれ以上の分散コントローラ１４Ｂに接続された一つ又はそれ以上のオペレータインターフェイス１４Ａを含んでいる。分散コントローラ１４Ｂは、例えば、Emerson Process Management of Austin、テキサス、によって販売されているＤｅｌｔａＶ（登録商標）コントローラ又は他の所望のタイプのコントローラであり得る。分散コントローラ１４Ｂは、Ｉ／Ｏデバイスを介して、一つ又はそれ以上の、例えばＨＡＲＴ又はＦｉｅｌｄｂｕｓフィールドデバイス又は他のスマート若しくは非スマートタイプのフィールドデバイスのようなフィールドデバイス１６に接続され、これには、例えば、ＰＲＯＦＩＢＵＳ、ＷＯＲＬＤＦＩＰ、Ｄｅｖｉｃｅ−Ｎｅｔ、ＡＳ−ＩｎｔｅｒｆａｃｅａｎｄＣＡＮｐｒｏｔｏｃｏｌｓの何れかを使用するものが含まれる。既知であるが、フィールドデバイス１６は、アナログ又はデジタル情報を他のデバイス情報と同様にプロセス変数に関連するコントローラ１４Ｂに提供する。オペレータインターフェイス１４Ａは、例えば制御最適化、診断エキスパート、ニューラルネットワーク、チューナー等を含む、プロセスの制御のためにプロセス制御オペレータにとって利用可能なツール１７，１９を格納し及び実行する。
【００４７】
また更に、ＡＭＳ（登録商標）アプリケーション及び／又は以下に記述される監視、診断及び通信アプリケーションを実行するコンピュータのような保守ツールは、保守、監視及び診断アクティビティを実行するために、プロセス制御システム１２および１４に、又はその中の個々のデバイスに接続され得る。例えば、保守コンピュータ１８は、コントローラ１２Ｂ及び／又はデバイス１５に所望の通信ライン又はネットワーク（ワイヤレス、ハンドヘルドデバイスネットワークを含む）を介して接続され得、デバイス１５上の他の保守アクティビティと通信し、いくつかの例では、再構築し又はそれを実行する。同様に、ＡＭＳＴＭアプリケーションのような保守アプリケーションは、分散型のプロセス制御システム１４に関連する一つ又はそれ以上のユーザインターフェイス１４Ａにインストールされ、それによって実行されて、デバイス１６の運転状態に関するデータ収集を含む保守及び監視機能を果たす。
【００４８】
また、プロセスプラント１０は、タービン、モーター等の種々の回転（及び他の）装置を有し、これらは、保守コンピュータ２２に常設又は一時的な通信リンク（読み取り次に取り外される装置２０に接続されたバス、ワイヤレス通信システム又はハンドヘルドデバイス）を介して接続されている。以下のアプリケーション、モジュール及びツールと同様に、回転装置２０及びプラント内の他の装置の動作状態を診断し、監視し及び最適化するために、保守コンピュータ２２は、ＣＳＩ（Emerson Process Management Company）によって提供されているような商業的に入手可能なアプリケーションを含む多くの監視及び診断アプリケーションを格納し実行し得る。保守要員は、プラント１０内の回転装置２０の保守及び性能を監視し、回転装置２０に関する問題を決定し、及び回転装置２０を修理し又は取り替えるかどうか及びその時期を決定するために、アプリケーション２３を通常使用する。いくつかのケースでは、外部コンサルタント又はサービス組織が、装置２０に関連するデータを一時的に入手し又は測定し得、回転装置２０の分析を行って、このデータを、問題、低い性能又は回転装置２０に影響を与える他の問題を検出するのに使用する。これらのケースでは、分析を行っているコンピュータは、通信ラインを介してシステム１０の他の部分に接続されているとは限らず、一時的にのみ接続されているかも知れない。
【００４９】
同様に、システム１０に関連する電力発生分配装置２５を有する電力発生分配システム２４は、例えばバスを介して実行されシステム１０内の電力発生分配装置２５の動作を監視するもう一つのコンピュータ２６に接続されている。コンピュータ２６は、例えばLiebert and ASCO又は他の会社によって提供されているような公知の電力制御及び診断アプリケーション２７を、電力発生分配装置２５を制御し及び保守するために実行し得る。ここでも、多くのケースでは、外部コンサルタント又はサービス組織が、電力発生分配装置２５に関連するデータを一時的に入手し又は測定し得、装置２５の分析を行って、このデータを、問題、低い性能又は装置２５に影響を与える他の問題を検出するのに使用する。これらのケースでは、分析を行っている（コンピュータ２６のような）コンピュータは通信ラインを介してシステム１０の他の部分には接続されているとは限らず、一時的にのみ接続されているかも知れない。
【００５０】
図１に図示したように、コンピュータシステム３０は、異常状態防止システム３５の少なくとも一部を実現し、特に、コンピュータシステム３０は構成アプリケーション３８及び任意には異常動作検出システム４２を格納し実現し、その多くの実施形態が以下でより詳細に説明されるであろう。加えて、コンピュータシステム３０は、警報／警告システム４３を実現し得る。
【００５１】
一般的に言えば、異常状態防止システム３５は、異常動作検出システム、フィールドデバイス１５、１６に任意に配置されているモジュール又はツール（図１には示されていない）、コントローラ１２Ｂ、１４Ｂ、回転装置２０又はその支援コンピュータ２２、電力発生装置２５又はその支援コンピュータ２６、及びプロセスプラント１０内の他の所望のデバイス及び装置、並びに／又はコンピュータシステム３０内の異常動作検出システム４２と通信し得、これらの異常動作検出システムのそれぞれを構成し、デバイス又はそれらが監視を行っているサブシステムの動作に関する情報を受け取る。異常状態防止システム３５は、有線接続バス４５を介してプラント１０内の少なくともいくつかのコンピュータ又はデバイスに接続され、あるいは、例えばワイヤレス接続、ＯＰＣ（プロセス制御のためのＯＬＥ）を使用する専用接続、データを集めるためのハンドヘルドデバイスに依存するような間欠的な接続などを含む他の所望の通信接続を介して接続され得る。同様に、異常状態防止システム３５は、プロセスプラント１０内のフィールドデバイス及び装置に関するデータをＬＡＮ又はインターネット、電話接続など（図１ではインターネット接続４６として示されている）の公衆接続を介して取得し得、このようなデータは、例えば第三者のサービスプロバイダによって収集される。更に、異常状態防止システム３５は、プラント１０内のコンピュータ／デバイスに、、例えばイーサネット（Ｒ）、Ｍｏｄｂｕｓ，ＨＴＭＬ，ＸＭＬ，専有技術等の種々の技術及び／又はプロトコルを介して通信可能に接続される。従って、ＯＰＣを使用する異常状態防止システム３５をプラント１０内のコンピュータ／デバイスに接続する特定の例がここに記載されているけれども、当業者は異常状態防止システム３５をプラント１０内のコンピュータ／デバイスに接続する種々の他の方法を同様に使用し得ることを認識するであろう。何れの場合にも、異常状態防止システム３５は、コンピュータ、デバイス又は他のプラント１０の他の特徴と通信し及びこれらからプロセス変数データを受け取り、これらには、プロセス制御システム（例えば、ＤｅｌｔａＶ（登録商標）、Ｏｖａｔｉｏｎ又は他の分散制御システムシステム）、デバイス、種々の標準プロトコル（例えば、ＦＯＵＮＤＡＴＩＯＮ（登録商標）  Ｆｉｅｌｄｂｕｓ，  ＨＡＲＴ，ＯＰＣ，Ｍｏｄｂｕｓ，ワイヤレス）に従うコンピュータ、及び種々のトランスミッタ、センサ、デバイスと共に実装され又はプロセスプラント１０内に分散しているアクチュエータなどが含まれるが、これらに限定されるものではない。更に以下に論ずるように、異常状態防止システム３５によって受信され使用されたデータは、データヒストリアンからのデータのような履歴データであるが、オンラインデータ（例えば、プロセスがオンラインのときに収集されるデータ）であり、これは、必ずしもこれに限定されないが、リアルタイムで収集されたデータと同様に、プロセスが動作しているときにオンライン、アットライン又はインラインで収集されたデータを含んでいる。
【００５２】
バックグラウンドとして、ＯＰＣは、プラント又はプロセス制御システムからプロセスデータにアクセスするメカニズムを確立する標準である。典型的には、ＯＰＣサーバは、例えばフィールドデバイスからのプロセス情報を顕在化させ又は提供するプロセス制御システムにおいて実現される。ＯＰＣクライアントは、ＯＰＣサーバへの接続を生成し、フィールドデバイスへのプロセス情報の書き込み又はフィールドデバイスからのプロセス情報の読み出しを生成する。ＯＰＣサーバは、ＯＬＥ技術（即ち、コンポーネントオブジェクトモデル又はＣＯＭ）を使用して、クライアントによって実装されたソフトウエアアプリケーションがフィールドデバイス又は他のプラント装置からデータにアクセスし得るように、このようなクライアントと通信を行う。
【００５３】
図２は、図１の例示のプロセスプラント１０の部分５０を表し、異常状態防止システム３５及び／又は警報／警告システム４３が例示のプロセスプラント１０の部分５０における種々デバイスと通信し得る一つの方法を記述している。一方、図２は、異常状態防止システム３５と、ＨＡＲＴおよびＦｉｅｌｄｂｕｓフィールドデバイス内の一つ又はそれ以上の異常動作検出システムとの間の接続を例示し、同様の通信が、異常状態防止システム３５と、図１に示されている他のデバイスおよび装置の何れかを含むプロセスプラント１０内の他のデバイス若しくは装置との間に出現し得ることが理解されるであろう。
【００５４】
図２に図示されているプロセスプラント１０の部分５０は、所望のどのような通信若しくはコントローラプロトコルにも適合する所望のどのようなタイプのＩ／Ｏデバイスであってもよい一つ又はそれ以上のフィールドデバイス６４および６６に入力／出力（Ｉ／Ｏ）カード又はデバイス６８および７０を介して接続される一つ又はそれ以上のプロセスコントローラ６０を有する分散型のプロセス制御システム５４を含んでいる。プロセス制御システム５４およびフィールドデバイス６６は、所望の他のどのような通信プロトコルも使用することができるが、ＨＡＲＴフィールドデバイスとして図示され、フィールドデバイス６６はＦｉｅｌｄｂｕｓフィールドデバイスとして図示されている。加えて、フィールドデバイス６４および６６のそれぞれは、例えば、センサ、バルブ、トランスミッタ、ポジショナ等のようなどのようなタイプであってもよく、所望のどのような解放された、占有された、又は他の通信若しくはプログラミングプロトコルであってもよく、Ｉ／Ｏデバイス６６および７０は、フィールドデバイス６４および６６によって使用されている所望のプロトコルと互換性を有していなければならないことが理解される。
【００５５】
何れにしても、構成エンジニア、プロセス制御オペレータ、保守要員、プラントマネージャー、スーパーバイザなどのプラントの人間によってアクセス可能な一つ又はそれ以上のユーザインターフェイス又はコンピュータ７２および７４（これは、どのようなタイプのパーソナルコンピュータ、ワークステーションなどであり得る）が、所望のどのようなハードウエア又はワイヤレス通信構造を使用しても、および例えばイーサネット（Ｒ）プロトコルのような所望のどのような若しくは適切な通信プロトコルを使用しても、通信ライン又はバス７６を介して、プロセスコントローラ６０に接続され得る。加えて、データベース７８が通信バス７６に接続されて、オンラインプロセス変数データ、パラメータデータ、状態データおよびプロセスプラント１０内のプロセスコントローラ６０およびフィールドデバイス６４および６６に関連する他のデータと同様に、構成情報を収集しおよび格納するデータヒストリアンとして動作する。このように、データベース７８は、プロセスコントローラ６０およびフィールドデバイス６４および６６内にダウンロードされ格納されたときに、プロセス制御システム５４に対する制御構成情報と同様にプロセス構成モジュールを含む現在の構成を格納する構成データベースとして動作する。同様に、データベース７８は、プロセスプラント１０内のフィールドデバイス６４および６６によって収集された統計データ、フィールドデバイス６４および６６によって収集されたプロセス変数から決定された統計データ、および以下に記述する他のタイプのデータを含む、履歴異常状態防止データを格納する。
【００５６】
プロセスコントローラ６０、Ｉ／Ｏデバイス６８および７０、並びにフィールドデバイス６４および６６は、典型的には、過酷なプラント環境内におよびこれに分散して配置されるが、ワークステーション７２および７４は、通常、オペレータ、保守要員等が容易にアクセスし得る制御室、保守室又は他のあまり過酷でない環境に位置している。
【００５７】
一般的に言えば、プロセスコントローラ６０は、多くの異なる独立して実行される制御モジュール又はブロックを使用して、制御戦略を実現する一つ又はそれ以上のコントローラアプリケーションを格納しおよび実行する。制御モジュールのそれぞれは、一般に機能ブロックと称されるものによって構成され、各機能ブロックは、全体に亘る制御ルーチンの一部又はサブルーチンであり、他の機能ブロックと（リンクと称される通信を介して）結合して動作してプロセスプラント１０内のプロセス制御ループを実現する。よく知られているように、オブジェクト指向プログラミングプロトコルにおけるオブジェクトである機能ブロックは、トランスミッタ、センサ又は他のプロセスパラメータ測定デバイスに関連する入力機能、ＰＩＤ、ファジー論理などを実行する制御ルーチンに関連する制御機能、バルブのようないくつかのデバイスの動作を制御してプロセスプラント１０内のいくつかの物理的機能を実行する制御若しくは出力機能のうちの一つを典型的に実行する。もちろん、モデル予測制御（ＭＰＣｓ）、最適化などのハイブリッドおよび他のタイプの複雑な機能ブロックも存在する。ＦｉｅｌｄｂｕｓプロトコルおよびＤｅｌｔａＶ（登録商標）システムプロトコルは、オブジェクト指向プログラミングプロトコルにおいて設計され実現された制御モジュールおよび機能ブロックを使用するが、制御モジュールは、例えば、連続機能ブロック、ラダーロジックなどを含む所望の制御プログラミングスキームを使用して設計され得、機能ブロック又は他の何れか特定のプログラミング技術を使用して設計されるように限定されることはないことを理解すべきである。
【００５８】
図２に示すように、保守ワークステーション７４は、プロセッサ７４Ａ、メモリ７４Ｂおよびディスプレイデバイス７４Ｃを有している。メモリ７４Ｂは、図１に関連して説明した異常状態防止アプリケーション３５および警報／警告システム４３を格納して、これらのアプリケーションがプロセッサ７４Ａ上で実行され得てディスプレイデバイス７４Ｃ（又はプリンタのような他の何れかの表示デバイス）を介して情報がユーザに提供されるようにする。
【００５９】
一つ又はそれ以上のフィールドデバイス６４および６６のそれぞれは、以下に述べるように、異常動作検出のための検知デバイス又はルーチンによって送られる一つ又はそれ以上のプロセス変数に関連する統計的データの集合を実装するためのルーチンのような格納ルーチンのためのメモリ（図示せず）を含んでいる。また、一つ又はそれ以上のフィールドデバイス６４および６６のそれぞれは、統計的データの集合を実装するためのルーチン及び／又は異常動作検出のためのルーチンのようなルーチンを実行するプロセッサ（図示せず）を有している。統計的データの集合及び／又は異常動作検出は、必ずしもソフトウエアによって実現されるとは限らない。むしろ、当業者は、このようなシステムは、一つ又はそれ以上のフィールドデバイス及び／又は他のデバイス内のソフトウエア、ファームウエア及び／又はハードウエアによって実現され得ると認識するであろう。
【００６０】
図２に示すように、フィールドデバイス６４および６６のいくつか（潜在的には全て）は、異常動作検出（即ち、異常状態防止）ブロック８０および８２を有し、これは以下により詳細に記述されるであろう。図２のブロック８０および８２は、デバイス６４の一つおよびデバイス６６の一つに位置しているものとして示されているが、これら又は同様のブロックは、フィールドデバイス６４および６６のいくつかにおいて配置され得、コントローラ６０、Ｉ／Ｏデバイス６８，７０又は図１に示されているデバイスの何れかのような他のデバイスに配置され得る。加えて、ブロック８０および８２は、デバイス６４および６６の何れかのサブセットであり得る。
【００６１】
一般的に言って、ブロック８０および８２又はこれらのブロックのサブエレメントは、プロセス変数データのようなデータを、それらが配されているデバイス及び／又は他のデバイスから集める。加えて、ブロック８０および８２又はこれらのブロックのサブエレメントは、変数データを処理し得、いくつもの理由のためにそのデータの分析を行う。例えば、バルブに関連しているとして示されているブロック８０は、バルブが固着状態にあるかどうかを決定するためにバルブプロセス変数データを分析する固着バルブ検出ルーチンを有し得る。加えて、ブロック８０は、バルブ内のプロセス変数又は他のデータを集めるブロックＳＰＭ１−ＳＰＭ４などの一つ又はそれ以上の統計的プロセス監視（ＳＰＭ）ブロック又はユニットのセットを含み得、例えば、収集されたデータの平均、メジアン、標準偏差、二乗平均平方根（ＲＭＳ）、変化率、範囲、最小値、最大値などを決定するために、及び／又は収集されたデータにおけるドリフト、バイアス、ノイズ、スパイク（急上昇）などを検出するために、集められたデータ上で一つ又はそれ以上の統計的計算を実行する。生成される特定の統計データも、それが生成される方法も、重要ではない。従って、異なるタイプの統計データが、上記の特定のタイプに加えて又はその代わりに生成され得る。加えて、公知の技術を含む種々の技術がこのようなデータを生成するために使用され得る。統計的プロセス監視（ＳＰＭ）ブロックという用語は、少なくとも一つのプロセス変数又は他のプロセスパラメータ上の統計的プロセス監視を実行し、そしてデバイス又はそのためにデータが収集される外部のデバイス内の所望のソフトウエア、ファームウエア又はハードウエアによって実行され得る機能性を表すために、ここでは使用される。ＳＰＭは一般的にデバイスデータが収集されるデバイスに配置されるので、ＳＰＭはより正確なプロセス変数データを定量的および定性的に取得することができることを理解すべきである。結果として、ＳＰＭブロックは、一般的に、プロセス変数データが収集されるデバイスの外部に位置するブロックより収集されたプロセス変数データに関してより良い統計的計算を行うことができる。
【００６２】
図２ではブロック８０および８２はＳＰＭを含むように示されているけれども、ＳＰＭブロックは、その代わりにブロック８０および８２から分離した独立型のブロックであり得、対応するブロック８０および８２と同じデバイスに又は異なるデバイスに配置され得ることを理解すべきである。ここで論じたＳＰＭブロックは、公知のＦＯＵＮＤＡＴＩＯＮ（登録商標）ＦｉｅｌｄｂｕｓのＳＰＭブロック、又は公知のＦＯＵＮＤＡＴＩＯＮ（登録商標）ＦｉｅｌｄｂｕｓのＳＰＭブロックと比較して異なる若しくは追加の能力を有するＳＰＭブロックを含み得る。統計的プロセス監視（ＳＰＭ）ブロックという用語は、プロセス変数のようなデータを収集し、そして平均、標準偏差などの統計的な測度を決定するためにこのデータ上でいくつかの統計的処理を実行するどのようなタイプのブロック又は要素にも言及するようにここでは使用される。結果として、この用語は、これらの要素が機能ブロックであろうとが他のタイプのブロック、プログラム、ルーチン若しくは要素であろうが、ＦＯＵＮＤＡＴＩＯＮ（登録商標）Ｆｉｅｌｄｂｕｓプロトコル若しくはＰｒｏｆｉｂｕｓ、ＨＡＲＴ、ＣＡＮ等のプロトコルに準拠するか否かに関係なく、この機能を果たすソフトウエア、ファームウエア又はハードウエア及び／又は他の要素をカバーすることを意図している。所望なら、ブロック８０および８２の基本的な動作は、参照によってここに取り入れられる米国特許第６０１７１４３号に記載されているように少なくとも部分的に実行され又は実現され得る。
【００６３】
ブロック８０および８２は、図２ではＳＰＭブロックを含むように示されているけれども、ＳＰＭブロックはロック８０および８２に必須ではないことを理解すべきである。例えば、ブロック８０および８２の異常動作検出ルーチンはＳＰＭブロックによって処理されていないプロセス変数を使用して動作することができる。もう一つの例として、ブロック８０および８２のそれぞれは、他のデバイスに位置する一つ又はそれ以上のＳＰＭブロックによって提供されるデータを受け取り、そレに基づいて動作することができる。さらにもう一つの例として、プロセス変数データは、多くの典型的なＳＰＭブロックによって提供されていない方法で処理され得る。丁度その一つの例として、プロセス変数データは、バンドパスフィルタやいくつかの他のタイプのフィルタのような有限インパルス応答（ＦＩＲ）又は無限インパルス応答（ＩＩＲ）によってフィルタされ得る。もう一つの例として、プロセス変数データは、それが特定の範囲に止まるように切り揃えられることができる。もちろん、公知のＳＰＭブロックは、このような異なる又は追加の処理の可能性を提供す量に改変されることができる。
【００６４】
トランスミッタに関連しているとして示されている図２のブロック８２は、トランスミッタによって収集されたプロセス変数データを分析してプラント内のラインが詰まっているか否かを決定する詰まりライン検出ユニットを有し得る。加えて、ブロック８２は、トランスミッタ内のプロセス変数又は他のデータを収集し、例えば収集されたデータの平均、メジアン、標準偏差などを決定するために一つ又はそれ以上の統計的計算をその収集されたデータについて行うブロックＳＰＭ１−ＳＰＭ４のような一つ又はそれ以上のＳＰＭブロック又はユニットを含み得る。ブロック８０および８２は、それぞれ４つのＳＰＭブロックを含むように示されているが、ブロック８０および８２は統計的データを収集しおよび決定するための他の数のＳＰＭブロックをその中に有し得る。
【００６５】
異常状態防止システムの実装および構成とその構成要素に関する詳細は、米国特許公開公報第２００５／０１９７８０３号、現在は米国特許第７０７９９８４号（「プロセスプラントにおける異常状態防止」）、米国特許公開公報第２００５／０１９７８０６号（「プロセスプラントにおける異常状態防止のための構成システムおよび方法」）および米国特許公開公報第２００５／０１９７８０５号（「プロセスプラントにおける異常状態防止のためのデータ表示システム」）に見られ、これらのそれぞれは、全ての目的のために参照によってここに取り入れられる。
【００６６】
上記および参照文献に記載されている異常状態防止システムおよび技術においては、ＳＰＭ（又は異常状態防止）ブロック８０，８２は、一つ又はそれ以上の異常状態防止モジュールに関連し又はその構成要素と考え得る。異常状態防止ブロックは、フィールドデバイス内にあり、そこではより早いサンプルデータが利用可能であり、異常状態防止モジュールはホストシステム又はコントローラ内に存在し得る。異常状態防止モジュールは、一つ又はそれ以上の異常状態防止ブロックからデータを取得し、このデータをより大きなシステムについての決定を行うのに使用される。より一般的には、異常状態防止モジュールは、フィールドデバイス、設備又は他の装置（例えばバルブ、ポンプなど）のタイプそれぞれについての診断をサポートするために、一つ又はそれ以上の機能ブロック（例えば異常状態防止ブロック）からデータを受け取るように改良されおよび構成され得る。とはいえ、異常状態防止モジュールに関連する機能ブロックは、それのために改良される特定の装置以外のデバイスに備えられ実装され得る。このような場合、異常状態防止モジュールは、分散特性を有している。他の異常状態防止モジュールは、特定のフィールドデバイスの診断専用であるにも拘わらず、プロセスコントローラ６０のような一つのデバイス内に全体が実装され得る。何れにしても、診断ルーチン又は技術は、装置（又はプロセス）の異常状態又は動作の検出、予測および防止のために各タイプの装置に対して開発され得る。記述を容易にする目的のみのために、「異常状態防止モジュール」という用語は、ここではこのようなルーチン又は技術に言及して使用されるであろう。従って、異常状態防止モジュールは、診断を行うのに必要な測定のセットに応答し、さらに、(ｉ)そのモジュールによって検出されるべき異常状態のセット、および(ii)規則のセット、を含み、測定値における変化を対応する異常状態に結びつける。更に、以下に開示されている技術の記述における異常状態防止モジュールへの参照は、その技術が同様に異常状態防止ブロックと関連して利用され得ることの理解とともに説明される。
【００６７】
いくつかの場合には、構成アプリケーション３８又は異常状態防止システム３５の他の構成要素は、各異常状態防止モジュールについてのテンプレートの改良又は生成をサポートし得る。例えば、ＤｅｌｔａＶ（登録商標）制御システムによって提供される構成および開発プラットフォームは、対応する混合テンプレートブロックからの異常状態防止モジュールの特定のインスタンス、又はインスタンス化を生成するのに使用され得る。
【００６８】
異常状態防止機能として図２に関連して示しおよび記述したけれども、上記モジュールおよびブロックは、以下に記述するプロセス監視、診断および故障検出のために構成された多変量統計解析技術の実装をより一般的に指向している。いくつかの場合には、以下の技術は、異常状態防止モジュール又はブロックを含み又はこれらに統合され得る。何れの場合にも、システムおよび技術（および何れのモジュール、機能ブロック、アプリケーション、ソフトウエア又は他の構成要素若しくはその特徴）への以下の参照は、上述のワークステーションツール１７、１９、オペレータインターフェイス１２Ａ、１４Ａ、アプリケーション２３、異常状態防止システム２５およびインターフェイス７２、７４で使用され、含まれ、統合されさもなければ関連づけられ得る。
【００６９】
図３に移行、ここに記述されている多くの多変量監視及び診断技術が、視覚化および多変量および単変量の診断ツールの特徴とともに構築され組み込まれている。例えば、プロセス変数（ＰＶ）がその対応するプロセス要求とともに、図３におけるしばしばシューハートチャートとして言及されるチャートに示されている。各プロセス変数は、上部制御境界（ＵＣＬ）１００、下部制御境界（ＬＣＬ）１０２、および目標値１０４をプロットしたシューハートチャートの観点から記述されている。ＵＣＬ１００およびＬＣＬ１０２は真の物理的境界であり、統計的な境界ではない。動作中、各プロセス変数はそのＵＣｌおよびＬＣＬの間に束縛され、プロセス変数が目標値である理想的な値であることが好ましい。これらの境界内では、プロセス変数は制御されていると言え、そうでない場合は制御不能と言われる。上記のように、ＵＣｌおよびＬＣＬは物理的境界を表し、プロセス動作の要件によって提供され得る。この視覚化の方法は、プロセス変数データの多面的見解に基づいているけれども、プロセス変数のこの記述は、以下に記述するプロセス変数の正規化に基づくさらなる視覚化技術との結合において有用であり得る。
【００７０】
しかしながら、ＰＣＡの基本のアルゴリズムの要約は、以下にさらに記述される種々の特徴への参照に提供される。ＰＣＡはここでは、この開示の種々の特徴とともに使用される多変量統計解析として開示されているけれども、主成分回帰（ＰＣＲ）、部分最小二乗法（ＰＬＳ）、Ｆｉｓｈｅｒ判別分析（ＦＤＡ）および正準変量分析（ＣＶＡ）等を含みしかしこれらに限定されない他の多変量統計解析が、ＰＣＡに代えて使用され得ることが理解されるべきである。
【００７１】
与えられたプロセス、プロセスユニット、プロセスデバイスなどに対して、多くの測定されたプロセス変数が存在する。これらのそれぞれは、ＰＣＡアルゴリズムへの入力変数として使用することができる。プロセス変数のそれぞれからのデータは、ある時間の長さについて、あるいはある数のポイントについて収集され、ここではプロセス変数の観測結果として言及される。一般的には、ＰＣＡにおいて使用されるプロセス変数の全ては、同時にサンプリングされる。
【００７２】
もし、ｍが入力変数の数であり、ｎが各入力変数の観測結果の数であるなら、行列Ｘは、全ての入力変数に対する観測値の全てを含むｎ×ｍの行列である。典型的なプロセスにおいては、いくつかの変数が他より著しく大きな重要性を有している。全てのプロセス変数はＰＣＡモデル上で等しい影響を有していることを保証するために、Ｘデータは自動スケール化（各入力変数に対して、平均を差し引き、そして標準偏差で除算する）され得る。入力変数が自動スケール化された後、標本共分散行列が、Ｓ＝ＸT・Ｘ／（ｎ−１）によって計算され、ここでＴは行列Ｘの転置演算である。
【００７３】
固有値−固有ベクトル分解が標本共分散行列：Ｓ＝Ｖ・Ｄ・ＶTの上で行われ、ここでＤは、ｍ個の固有値を含む対角行列であり、最大から最小まで分類されている。Ｖの列は固有ベクトルに対応し、Ｔは行列Ｖの転置演算である。最大の固有値およびそれに対応する固有ベクトルは、オリジナルのプロセスデータにおける変動の最大量を含むオリジナルの入力空間における方向に相当している。これは第１の主成分であると考えられる。２番目に大きい固有値およびそれに対応する固有ベクトルは、変動の次の最大値を含み、第１の主成分に直交する方向に相当している。このことが、全てのオリジナルのデータを表す新たな直交ベクトル空間が生成されるまで続けられる。
【００７４】
固有値および固有ベクトルが計算された後、ＰＣＡは、固有値および固有ベクトルの何れの一つが、通常プロセス内での重要な変動に対応する重要な主成分であるか、何れが通常はノイズに相当する重要でないものであるかを決定する。Ｋａｉｓｅｒの規則（即ち、１より大きい固有値を有する主成分を選択する）、パーセント分散テスト（即ち、データ中の全変動の所定のパーセント（例えば８０％又は９０％）を説明する主成分の数を選択する）、平行分析、相互検証などを含む、重要な主成分の数を選択するための多くの異なる方法が提案されている。
【００７５】
何れの方法が使用されるにせよ、最終的な結果は、重要な主成分の数、ａ、を決定し、これはＰＣＡモデルに対して保持されるであろう。次に、負荷行列Ｐ∈Ｒm x aが、Ｖから最初のａ列をとることにより生成される。オリジナルの観測値の負荷行列によって定義される部分空間への射影はスコア行列と呼ばれ、Ｔ＝Ｘ・Ｐによって表記される。Ｔの列はＰＣＡモデルのスコアと呼ばれ、ｔ1、ｔ2...ｔaとして言及される。負荷行列Ｐは、自動スケール化で使用される平均および標準偏差に沿って、まとめてＰＣＡモデルであると考えられている。このＰＣＡモデルは、次に、今後のデータセットに適用される。
【００７６】
ＰＣＡモデルが生成された後、統計に基づいた境界が、オリジナルのプロセスデータに基づいて、オリジナルのデータによって定義されるプロセスの正常な動作を定義するために生成される。使用される一つの方法は、ＨｏｔｅｌｌｉｎｇのＴ2テストであるが、統計に基づいた境界を誘導するための他の方法も使用し得ることが理解されるべきである。スコアが境界を超えたとき、何かがプロセス内で異常であることを示し得る。また、変動の量を制限する２−σおよび３−σの境界も、代わりに適用可能である。
【００７７】
ＰＣＡモデルからのスコアデータに基づいて、いくつかのプロットを行うことができる。最も単純なプロットは、スコア対時間のプロットである。このプロットでは、上限および下限の境界が統計に基づいて生成され、何れかの境界が超えられたときに、アラームが誘起される。第２の共通のプロットは、スコアの２つのプロットである。最もしばしば最初の２つのスコアｔ1、ｔ2がプロットされ、その理由は、これらはプロセスにおける最大の変動量を記述するからである。多くの例示の２−スコア散布図が以下に記述されるユーザインターフェイスに示されている。これらおよび他の場合において、この統計に基づいた境界は、正常なプロセスデータの周囲の円又は楕円を含み又は伴っている。再び、これらの境界が超えられると、アラームが誘起される。
【００７８】
図４および図５に移行して、一般的に多変量実時間プロセスデータのオンライン監視及び分析のサポートを指向する開示の特徴に従って、プロセス視覚化技術は、上記シューハートチャートに関連する制御境界のような、プロセス変数データの多次元（又は多変量）のセット内の境界の制御に使用されている。多変量データセットにもかかわらず、境界は以下に記述するような変換を使用して、単変量の形式で表現され得る。一般的に言えば、この変換は、多変量データの監視のための単変量技術に基づいている。この技術を使用すれば、プロセス変数のいかなる数のプロセスデータおよび制御境界も、単一の射影パラメータを使用することにより、正規化されて同じ基準をもたらすことができる。
【００７９】
これに対して、オフラインアプリケーションに対する異なる統計的ツールが存在し、それらのいくつかはオンライン監視に使用することができる。例えば、上記のように、主成分分析（ＰＣＡ）は多変量データ空間の次元を数次元（しばしば２−３次元）に減少させるのにしばしば使用される技術である。しかしながら、ＰＣＡを実装するために、ＰＣＡモデルは訓練データ又は「正常な」プロセス動作からのデータを使用して改良され、次に、そのモデルはプロセスのさらなる監視に適用される。上記で論じたように、誘導されたモデルは主成分に関するものであり、これは本質的に訓練データによって決定された共分散行列の支配的な固有ベクトルである。そのようなＰＣＡの結果も、スコアとして言及され、真のプロセス変数に直接関連づけられておらず、そのため、監視されている変数の物理的プロセス境界は解析において失われている。この理由により、ＰＣＡツールを使用するときは、いくつかの統計的解釈が、「データの良好性」を定量化するためにしばしば使用される（即ち、９５％又は９９％の信頼区間がＨｏｔｅｌｌｉｎｇのＴ2統計を使用して計算される）。
【００８０】
単変量解析を使用して、決定論的な方法が、多変量プロセスデータの監視のために、単一監視変数とともに提案されている。この方法は、オフラインおよびオンラインデータ監視および分析アプリケーションの両方に容易に適用され得る。この方法およびモデルの展開は、以下に記述される。
【００８１】
はじめに、Ｘ1、Ｘ2．．．Ｘnはプロセス変数を表しているとする。各変数Ｘi、ｉ＝１、．．．、ｎは、時間に依存するデータｘi,jを含むベクトルであり、ｊはＸiのｊ番目のサンプル又は成分である。例えば、Ｘi＝［ｘi,1，ｘi,2，．．．，ｘi,m］Tで、ＴはＸi上の転置演算を表し、ｍは、データサンプリング速度およびデータ収集のトータル時間によって決定されるベクトルの次元である。このような用語により、各プロセス変数は変換、ｐi＝［（ｙ−αi）2＋ｋiｘi,j2］に従って記述され、ｉ＝１，．．．，ｎであり、ｊ＝１，．．．，ｍであり、ｙは多変量プロセス変数を束ねるように指定された変数であり、αiはＸiに割り当てられた一意的な位置であり、ｋiはスケーリングパラメータであり、ｘi,jはＸiのｊ番目の成分である。変換ｐiは各変数Ｘi、ｉ＝１，．．．，ｎに対して一意的に定義され得るので、多変量変換は、Ｐ（ｙ）＝ｐ1（ｙ）＊ｐ2（ｙ）．．＊ｐn（ｙ）のように定義され得る。
【００８２】
次のＰが、射影モデルＭv（ｙ）＝Ｋp／Ｐ（ｙ）を誘導するのに使用され、ここで、Ｋpはスケーリングファクターであり、Ｍv（ｙ）は多変量射影を指定する。これらの帰属を使用して、多変量データが単一の変数ｙによって射影される。
【００８３】
図４および図５は、変数に対するそれぞれの物理的動作境界に関連する単一プロット上の多重プロセス変数（ＰＶ１、ＰＶ２等）の視覚化の結果を示している。図４は、各プロセス変数が動作境界内に留まっているデータポイントに対するプロット１０６を表している。図５は、プロセス変数の一つがその制御境界を超えた場合のもう一つのデータポイントに対するプロット１０８を表している。
【００８４】
この方法（制御境界が超えられたことを視覚化がどのようにして表示するかを含めて）をさらに示すために、６つのプロセス変数Ｘ1，Ｘ2，Ｘ3，Ｘ4，Ｘ5およびＸ6を使用して例が提供されている。正常な動作条件の下では、これらの変数のそれぞれは、｜ＸiL｜（データポイントは中心に置かれていると想定する）によって定義された動作境界を有している。この情報の下に、各変換ｐｉは、係数ｋiと同じレベルに調整され、図４および図５の視覚化が得られる。これは、多変量監視の設定のグループを効率的に定義する。以上のように、６つのプロセス変数のそれぞれは、変数ｙ（水平軸）に関連して、αiの＝−１０，−５，０，５，１０，１５によって与えられる一意的な位置を有している。更に、全てのプロセス変数は、ｋiによって調整されて、ＸiLによって与えられる最大動作範囲に対する等しい最大応答又は境界を達成する。
【００８５】
さて、もし監視ステージの間に６つのピークの何れかがプロセス変数インディケータを貫いたら、特定の変数が範囲イベントから外れることに遭遇するであろう。これは図５に示してある。これから分かるように、プロセス変数＃２が許容範囲の境界に近く、プロセス変数＃４は明らかに境界を超えている。それでもなお、図５に示す状態は、６つのプロセス変数に対して正常な動作を明確に示している。
【００８６】
監視ステージの間、図４および図５のイベントは孤立していない又は単独ではないイベントである。それらは実際には、連続して変化する変数の視覚の断片又はフレームである。このことは、新たなプロセスデータがモデルに（与えられたサンプリング速度で）入力されるに伴って、応答は即座に変化する。従って、監視イメージは、０．１〜１秒又はそれより長いサンプリング速度に従ってリフレッシュされる。
【００８７】
この監視アプローチでは、全てのプロセス変数に対する境界が知られており、開示されている方法は訓練データを必要としない。例えば、多変量監視を生成するために、ある人は上記設計手法を適用し、そして単変量モデルが完成し、これは他のアプローチでは不可能である。他の場合には、もしプロセス変数の境界の全てが分かっているわけではないなら、予め決められた境界のないプロセス変数の訓練データが、正常なプロセス動作状態の下でプロセス変数の境界を決定するために使用され得、そのモデルが次に計算され得る。境界が既知又は未知の何れかの場合において、プロセス変数のそれぞれは、共通のプロセス変数の境界を決めるために、単変量の解析を使用して正規化され得（例えば、図４参照）、そして、プロセス変数のそれぞれが共通のプロセス変数の境界に対して監視され得る。正規化のさらなる例は、既知および未知の両方のプロセス変数の境界を使用して、以下にさらに論ずるように、単変量解析法により実装され得る。
【００８８】
開示されている単変量解析法は、プロセス変数の特定のセット又はタイプに限定されるものではない。それは、プロセス変数のどのような混合（例えば、プロセス入力、プロセス出力、および外部変数）にも適用することができる。加えて、単変量法は、プロセスデータのパターンを生成しおよび解析するユニークな方法を提供する。これらのパターンは、さらに異常状態防止アプリケーションを含むプロセス診断に使用され得る。例えば、ほとんどの時間、２又はそれ以上のピークが同期して上昇又は下降する傾向を有しているなら、その場合は対応するプロセス変数が関連づけられる。そうでなければ、これらの変数は無関係である。同様に、単変量技術は、プロセス変数が不良であるかどうかを評価するのに使用され得る。手短言えば、開示されている単変量データモデリングおよび可視化方法は、単一の変数を使用してオンラインおよびオフラインの状況の両方に対して、多変量監視の実装技術を提供する。全てのプロセス変数の境界が既知であるなら、この方法は、訓練データを要求せず、プロセスオペレータによって容易に理解される監視結果を提供する。
【００８９】
図４およびず５に示されている方法で表示するためのプロセス変数を正規化するのに使用され得る３つの代替の方法に関連するさらなる情報が提供される。それらを記述するために、図３に示されている典型的なプロセス変数の記述（シューハートチャート）が考慮される。この図では、プロセス変数は、ＵＣＬ１００とＬＣＬ１０２との間に束縛され、プロセス変数に対する理想的な値は、目標値である。これらの境界内で、プロセス変数は制御されていると言え、そうでなければ制御不能であると言える。上記で論じたように、ＵＣＬとＬＣＬは、プロセス動作の要求により提供される物理的な境界を表している。これに基づいて、以下のデータ正規化方法が使用され得る。
【００９０】
第１の方法は、プロセス変数の制御境界が各プロセス変数について既知であるときに使用され得る。特に、プロセス変数ＰＶは、制御境界および目標値の関数として正規化され得る。
【数１】
2010506256.tif 000002
ここで、目標値＝Ｔ，ＣＬ＝｜ＵＣＬ−Ｔ｜＝｜ＬＣＬ−Ｔ｜である。
【００９１】
第２の方法は、プロセス変数の制御境界が全てのプロセス変数については知られていないときに使用され得る。特に、プロセス変数ＰＶは、プロセスが正常に動作しているときのプロセス変数について収集された訓練データのセットの関数として正規化され、訓練データのセット内のプロセス変数の観測値を含んでいる。
【数２】
2010506256.tif 000003
ここで、正常なプロセスの動作は訓練データのセットＴ｛ＰＶ｝によって表され、ＭPVはＴ｛ＰＶ｝の平均であり、Ｔ｛ＰＶ｝においてＣＬ＝｜ＰＶ−ＭPV｜max である。
【００９２】
第３の方法は、統計的自動スケーリングアプローチを使用し、これも、プロセス変数制御境界が全てのプロセス変数について未知であるときに使用される。特に、プロセス変数ＰＶは、プロセスが正常に動作しているときのプロセス変数について収集された訓練データのセットの関数として正規化され、訓練データのセット内のプロセス変数の観測値を含んでいる。
【数３】
2010506256.tif 000004
ここで、正常なプロセスの動作は訓練データのセットＴ｛ＰＶ｝によって表され、σ2はＴ｛ＰＶ｝の分散であり、ＭPVはＴ｛ＰＶ｝の平均であり、ｋσ、ｋ＝１，２，．．．、ｎは統計的な境界（例えば、３σ）を決定する。
【００９３】
これらの方法の変形と他のユーザ定義制御境界（およびそれらの決定方法）は、要望に応じて記載されている技術に関連して利用され得る。
【００９４】
上述のデータ監視及び分析技術は、高速のオンライン視覚化及びプロセスデータの挙動の理解にとって理想的である。それは、プロセス診断および異常状態防止アプリケーションで使用され得る。例えば、Ｆｉｅｌｄｂｕｓ上を拡散する多数のインテリジェントデバイスからのデータは、素早く分析され監視されることができる。同様のことがセンサアレイおよびバイオマイクロアレイからのデータについても言える。監視結果は、プロセス最適化および質の制御にも使用され得る。更に、この方法は、データ分析を高めるために統計的データモデリング方法とともに使用され得る。例えば、スピードアップとスコアデータの釈明を高めるために、それはＰＣＡ法と並行して使用され得る。
【００９５】
本開示の特定の特徴に従って、上記視覚化技術は、以下に記述するように、他の多変量統計技術を介する後の処理のような他の方法でプロセスデータを表示するオペレータインターフェイスに統合され得る。例えば、それらの物理的な境界に関連するプロセス変数の視覚化は、プロセスの現在の状態についてのより正確な決定を可能とするために、ＰＣＡスコアプロットに関連して使用され得る。
【００９６】
上記単変量技術は、多変量実時間プロセスデータを監視する方法を提示する確定的な多変量データ解析ツール（その例は以下に記述される）に組み込まれ得る。このように、単変数監視モデルは、多変量データの挙動の視覚化に使用され得る。
【００９７】
より一般的には、このアプローチの利点は、全てのプロセス変数が単一の監視変数に束ねられることである。開示されている単変量の方法は、オンラインおよびオフラインアプリケーションの両方に適している。上述のように、開示されている方法は、全てのプロセス変数の動作境界が既知であるとき、訓練データを必要としない。プロセス変数に関連する動作境界は、単変量監視モデルを展開するのに使用される。いくつかのケースでは、訓練データは、その動作境界が既知でないときのプロセス変数に対して使用される。このモデルは、例えば、ＤｅｌｔａＶ（登録商標）を含む多くのソフトウエアプラットフォームで実装され得る。
【００９８】
図６および図７を参照して、本開示の他の特徴は、主成分分析（ＰＣＡ）および他の多変量統計技術によって提供される機能を指向している。以下に記述されるＰＣＡおよび他の多変量統計技術は、分散型制御システムのＤｅｌｔａＶ（登録商標）およびＯｖａｔｉｏｎなどのプロセス制御システムの統合された構成要素であり得る。このような統合は、オンライン監視アプリケーションに対するＰＣＡ技術のプラント要員の利用をサポートする。上述のように、ＰＣＡおよび他の技術によってサポートされる視覚化及びアラームは、図４および図５に関連する上記記載のもののように、他の視覚化スキームに沿って制御システムに完全に統合され得る。これらのおよび他の理由により、プラント要員は、別個のアプリケーションを監視し、保守し、サポートする必要性無しに、簡素化され又は単一のユーザインターフェイスを利用し得る。また、プロセス制御システムの一部としてＰＣＡを実装することは、以下に記述するように、オンラインプロセスデータを使用する監視及び診断を実装することをより実用的にする。
【００９９】
図６および図７は、ＰＣＡベースの技術を、以下で論ずるように、プロセス制御システムに実装するための例示の機能ブロック１１０、１１２を示している。一般的に言えば、各ＰＣＡ機能ブロック１０、１１２は、訓練期間にデータを収集することができ、訓練機関の終わりにＰＣＡモデルを発展させ、次にそのＰＣＡモデルを全ての将来のデータに適用する。機能ブロック１１０，１１２の左側には、未処理のプロセス変数であるＰＣＡへの入力がある。機能ブロックの出力は、最重要の固有値に関連するスコアである。
【０１００】
例えば、機能ブロック１１０，１１２は、テキサス州オースチンのＥｍｅｒｓｏｎ  Ｐｒｏｃｅｓｓ  Ｍａｎａｇｅｍｅｎｔによって販売されているＤｅｌｔａＶ（登録商標）のＣｏｎｔｒｏｌＳｔｕｄｉｏを介して利用可能にされた特別注文の機能ブロックであり得る。この例では、２０個までのプロセス変数がＰＣＡへの入力として提供され得る。もちろん、一般的に入力がいくつであってもＰＣＡに使用することができる。ＰＣＡ機能ブロック１１０，１１２への入力は、ＤＣＳで標準的に利用可能なプロセス変数である。従属のプロセス変数（ＹＩＮ）に対する入力も、主成分回帰（ＰＣＲ）を行うのに使用することができ、これはここに記載されている開示された一つ又はそれ以上の方法において利用可能な統計的技術である。開示されている例は、単一の変数であるけれども、一般的には、多数従属変数でもあり得ることが理解されるべきである。
【０１０１】
ＰＣＡ機能ブロック１１０，１１２は、ＰＣＡモデルの構成およびアプリケーションの両方に対する多くの動作モードを含んでいる。説明している例示の実施形態では、ＰＣＡ機能ブロック１１０，１１２は、（ＬＥＡＲＮ機能ブロック上の真／偽によって決定されるように）学習と監視の２つの異なるモードで実行される。学習モードの間、この機能ブロックは入力変数のそれぞれからデータを収集する。例えば、ユーザが学習コマンドを与えると、機能ブロック１１０，１１２は、全ての入力変数についてのプロセスデータの収集を開始する。プロセスデータは、コンピュータメモリに格納され、データの収集は、ユーザが監視コマンドを与えるまで独立して継続する。
【０１０２】
十分な量のデータが収集された後、ユーザはＰＣＡブロックに監視コマンドを与える。このコマンドが与えられた後、機能ブロック１１０，１１２は、収集されたデータに基づいてＰＣＡを展開するためにＰＣＡアルゴリズムを実行する。モデルの展開は、以下のステップを含んでいる：入力データを自動スケーリングする（各入力変数について、その平均を減算する、その標準偏差で除算する）、スケーリングされたデータから共分散行列を計算する、固有値および固有ベクトルを計算する、ＰＣＡ負荷行列を保持するために重要な固有ベクトルの数を決定する、およびＰＣＡに対してσ−ベースの境界又はＴ2のＨｏｔｅｌｌｉｎｇの統計を使用して統計的境界を計算する（例えば、９５％又は９９％）。
【０１０３】
この計算が完了したとき、ＰＣＡ機能ブロック１１０，１１２は、監視モードに移行し、そこでそのブロックに入力された新たなプロセスデータに基づいてスコアを計算する。特に、新たなプロセスデータは、リアルタイムで生成されたオンラインプロセスデータである。もし何れかのスコアが境界を超えたなら、これは異常状態を示し、および機能ブロック１１０，１１２上のアラームパラメータがセットされる。アラームパラメータは、プロセス制御システム又はネットワークの他の部分の何れかと連携し得る。
【０１０４】
ＰＣＡの結果の視覚化プロットのそれぞれは、制御システムのオペレータインターフェイスの一部として提供され又は生成され得る。例えば、インターフェイスは、時間と共にスコアをプロットするＤｅｌｔａＶ  Ｐｒｏｃｅｓｓ  Ｈｉｓｔｏｒｙ  Ｖｉｅｗ  にＰＣＡスコアプロットを生成するのに使用され得る。これに代えて又はこれに加えて、２つのＰＣＡスコアに対するデータが、上記および以下で示しおよび言及するように、２次元のグラフにプロットされ得る。計算された統計的な境界は、オペレータインターフェイスを介して提供される円形又は楕円の何れかの形状を結果として生じる。このスコアが境界を超えると、アラームが誘起される。
【０１０５】
いくつかのケースでは、ここに記載した多変量監視及び診断技術は、異常状態防止を指向するシステムのモジュール（又は他の要素若しくはシステム構成要素）に実装され得る。従って、他の異常状態防止アルゴリズムと同様に、この技術は異常なプロセスの状態をそれが発生する前に検出することを指向している。これらのアルゴリズムは、多くの異常状態に対して「油断が無く」、正確で連続的であることを一般的に誇示している。従って、異常状態防止モードは、一時的ではなく、非リアルタイムの解決策ではなく、連続プロセスであり得る。
【０１０６】
ここに記載されている多変量監視及び診断技術は、バッチ処理においても使用することができ、そこでは、機能ブロック１１０，１１２は入力変数に対するプロセスデータの異なるセットを収集し、各セットは、プロセスの異なるオンライン状態に対応している。例えば、プロセスが種々の能力で開始されたときに、ユーザは、ＰＣＡ機能ブロックがプロセスが正常に動作しているときに特定の能力レベル（例えば、特定のオンライン状態）に対する入力変数についてのプロセスデータを収集するように学習コマンドを与える。その後、ユーザは、ＰＣＡアルゴリズムを実行して特定のオンライン状態（例えば、能力レベル）について収集されたデータに基づいてＰＣＡモデルを展開するように、ＰＣＡ機能ブロック１１０，１１２に対して監視コマンドを与える。その結果、プロセスが正常に動作しているときのそのプロセスに対するそれぞれ異なるオンライン状態に対応する種々のＰＣＡモデルが展開され得る。その後、プロセスが特定のオンライン状態（例えば、５０％の能力）で実行されているとき、ＰＣＡ機能ブロック１１０，１１２は、対応するＰＣＡモードを使用して監視モードに移行し、そこでそれはスコアをブロックに入力された新たなプロセスデータに基づいて計算する。従って、ユーザは、オンライン、プロセスを監視するプロセスのリアルタイム動作、異常状態の検出などを分析するために、プロセスのオンライン状態に最も近い一つを選択するように複数の多変量統計モデルから選択し得る。一実施例では、ユーザは、特定の訓練セッションを選択し、プロセスの特定のオンライン状態にカスタマイズされた新たなモデルを構築するように選択し得る。
【０１０７】
本開示のこの特徴は、オンライン状態（即ち、リアルタイム）の多変量監視および診断アプリケーションに対する一般的な（オープンソースの）データ駆動アプローチを提示している。更には、このアプローチは、オンラインモデルの展開のためのデータの生成および操作と、診断、故障検出などで使用されるモデルパラメータの比較とを一般的に提示している。開示されているアプローチは、プロセスの監視および診断の統合プラットフォームの定義を含んでおり、究極的には、テキサス州オースチンのＥｍｅｒｓｏｎ  Ｐｒｏｃｅｓｓ  Ｍａｎａｇｅｍｅｎｔによって販売されているＰｌａｎｔＷｅｂ構造のようなプラントワイド又はネットワーク化された構造を使用することにより実装され得る。提示されている統合プラットフォームは、多変量アプローチに基づいており、そして、以下に記述する検出技術で使用されているものを含む種々のコンピュータアルゴリズムに適合することができる。データ駆動プラットフォームのオープンソースの性質に鑑みれば、開示されている方法およびシステムは、異なる入力領域（例えば、ＦＯＵＮＤＡＴＩＯＮ（登録商標）  Ｆｉｅｌｄｂｕｓ、ＨＡＲＴ（登録商標）、Ｍｏｄｂｕｓ、ＯＰＣ、無線等）からのデータを受け取ることができる。受け取ったデータは、何れかの数の異なるトランスミッタ、センサおよびアクチュエータから収集されたオンラインプロセスデータであり得、多変量プロセス領域を定義するのに使用される。同様のデータが、プロセス制御システムによって使用され、データ履歴に格納され得る。
【０１０８】
開示されているシステムおよびモジュールを介するデータの流れは、図８に示されている。プロセスからのデータが一度開示されているモジュールに移行すると、そのデータは、同じ時間期間の間、オンライン訓練、オンライン監視又はその両方に使用される。２つのスイッチＳＷ１およびＳＷ２は、要望に応じて異なるデータパスをサポートする方法で実装される。もし訓練が望まれるなら、データは、ユーザ／オペレータが蓄積プロセスを停止するまで、蓄積される。このプロセスは、ＳＷ１をオンライン訓練状態にすることにより遂行される。ファイル１２０に格納されている蓄積データの何れは、次にモデルビルダー又はジェネレータ１２２において使用される。生成されたそれぞれのモデルは、示されているように格納され、監視モードに関係している構成要素によって、例えばＳＷ２を介して利用可能にされる。モデル構築の動作が終了すると、訓練データとともにモデルは、典型的にはｔｘｔファイルとして、データベース又はメモリ（ドライブＣ：／）にあるモデルホルダに格納される。
【０１０９】
このデータによって構築されたモデルは、これらに限定されるものではないが、主成分分析（ＰＣＡ）、人工神経ネットワーク（ＡＮＮ）、ファジー論理およびベイズ意思決定木などを含む。アルゴリズムのそれぞれのためのツールは、例えば異常状態防止モジュールにより、ドロップダウンウィンドウの様式で提供され得る。これらのモジュールは、以下に記述するように、予測の展開に使用され得る。同様のモデルが、オンライン監視にも使用され得る。モデルがオンライン監視モードにあるとき、異常状態防止モジュールはモデル出力（Ｏ／Ｐ）を使用してツリータイプの出力：統計的に正規化されたプロセス変数（ＰＶ）およびプロセス変数（ＰＶ）の評価を生成する。図８に示すように、出力のそれぞれ（ＰＶの評価、正規化されたＰＶ＆物理的境界、およびスコア＆統計的境界）は、アラーム状態（アラーム診断）及び／又はモードの検証を生成するためにさらに使用される。
【０１１０】
統計的出力は、多変量データポイントの低次元の表現であるスコアにより構成されている。低次元空間が、多変量領域における最大のデータ変動を有する方向に基づいて得られる。ＰＣＡモデリングは、上記で論じたように、スコアを展開するのに使用されるが、他のモデリング技術も同様にスコアを生成する。また、統計的出力は、９５％（９９％）の信頼境界又はσ−ベースの境界を含み得る。
【０１１１】
正規化されたプロセス変数の概念は、既知および未知の両方のプロセス変数の境界に関連して記載された。この概念に基づき、バーチャート（又は他の二次元プロット又はグラフ）が正規化されたデータを表すのに使用され得る。これらのチャートは、また、アラームの特定に有用な物理的又はプロセスの境界を含んでいる。本開示の更なる特徴に関連して以下に記述するように、単一のユーザインターフェイスにおける統計的境界と物理的境界の組合せは、異常状態防止アラーム状態を正確に特定する有用な方法を提供し、異常状態防止構想の支えとなる技術を提示する。
【０１１２】
プロセス変数の評価は、基本的にプロセス変数の可変性に対して与えられる評価であるので、異常状態防止に関して有用である。プロセス変数がより大きい変動を示せば、その評価は高くなる。これは、同じプロセス状態の異なる時点での監視に対して得られる評価が容易に交叉して確認され得るという点で有用である。例えば、もし、全てのプロセス変数が２つ又はそれ以上のモニタリングモデルについて同じ評価であるなら、モニタリングモデルの信頼性は増大する。プロセス変数の評価を得るために、以下の計算が使用され得る：
【数４】
2010506256.tif 000005
ここで、ＰＶjratingはｊ番目のプロセス変数の評価であり、Ｓはｎ×ｎの自動スケールされた共分散行列であり、ｎは多変量データ空間の次元であり、ｔｒＳはＳのトレースであり、ｋ＜ｎは主成分空間であり、ｉ＝１，．．．，ｎに対するσiおよびｐiはそれぞれ主固有値およびＳの固有ベクトルであり、ｐi,jは固有ベクトルｐiのｊ番目の成分であり、ｐi,jｓは負荷として言及される。評価は展開されたＰＣＡのそれぞれについて計算され、他のモデルパラメータとともにモデルファイルに格納され得る。
【０１１３】
上記の式は、モデルパラメータからなり、このようにして得られたプロセス変数評価は、そのモデルに密接に関係している。従って、２つのモデルが正確に同じプロセスイベントを記述しているなら、それらは同じオーダーであるとしてプロセス変数を評価する。その場合、モデルは、それぞれの評価に基づいて検証され得る。このことはまた、モデルが異なる場合にも一般的に当てはまる。例えば、一方のモデルがＰＣＡベースであり、他方がＡＮＮベースであり、それらのプロセス変数評価が大体同じオーダーである。プロセスデータが非線型である（即ち、正規分布から大きく外れている）場合は、例外が含まれ得る。
【０１１４】
他のモデルパラメータもモデルの検証に使用され得ることがさらに注目される。固有ベクトルの負荷ｐi,jおよび固有値σiのようなＰＣＡパラメータ、又は隠れ層の数および出力の数のようなＡＮＮパラメータもまた、有用なパラメータである。全てのモデルは、検証タスクに使用され得るパラメータの特定のセットを示す。
【０１１５】
本開示のもう一つの特徴に従えば、モデル出力データがアラーム状態、より一般的には、混合の、オペレータに対する多変量グラフィカルインターフェイスを生成するのに使用され得る。図９−図１３は、正規化されたプロセス変数、プロセス変数評価（図１３参照）、スコア、およびアラーム状態を表示する例示のグラフィカルインターフェイスを表している。グラフィカルインターフェイスはまた、他のモデル情報およびグラフがアクセスされるタブを含んでいる。図１３に示すように、プロセス変数評価は、図１３の上側左の変数の記述に隣接して、丸括弧内に与えられる（例えば(2)）。この特定の例では、１０のプロセス変数があり、１から１０の評価範囲である。プロセス変数Ｓ２１５３１は評価１を有し、最大の変動を有する変数であることをほのめかしている。
【０１１６】
図１４を参照すれば、ＢｉＰｌｏｔタブ１３０へのアクセスにより、オペレータに、スコアプロットパネル１３２に表示されているような２つのスコアの提示を見る機会が与えられる。プロセス変数とともに同時にスコアを表示するのにバイプロットが使用され、これはそれらの負荷ｐi,jによって定義されるグラフである。ラインは、どれだけ多くの主成分がプロセス変数内に存在するかを示している。例えば、評価１のプロセス変数Ｓ２１５３１は、最大の第１の主成分ｔ１（水平軸）を有している。これに対して、評価８のプロセス変数Ｓ２１００８は、最大の第３の主成分ｔ３（垂直軸）を有している。これは、ｔ１はデータ空間における最大の変動の方向を定義するので完全に筋が通っているが、ｔ３は著しく低い。プロセス変数の他の主成分との関係を、グラフの上側右の部分におけるウィンドウにおける表示を単に変更することにより、調査することができる（例えば、ｔ１−ｔ２に対するスコアが得られる）。
【０１１７】
このバイプロットは、オンラインスコア（黒い点）をプロセス変数に関連づける容易な方法を提供する。そのスコアに最も近いラインは、そのスコアに対する主要なプロセス変数を決定する。ラインは常に同じ位置にあり、スコアはグラフ上で一意的な動的成分なので、オペレータは、どのプロセス変数又はプロセス変数ぼグループが、与えられた時間に射影されたスコアに最も影響があるかを素早く理解することができる。従って、オペレータは、単一のディスプレイを介して全てのプロセス変数に素早く精通するようになることができる。
【０１１８】
実装のオプションに関して、開示されているユーザインターフェイスディスプレイ（およびそれらの生成の原因であるシステム又はモジュール）は、ＤｅｌｔａＶ（登録商標）、Ｏｖａｔｉｏｎ、又は他の制御システムプラットフォームに実装され得る。例えば、ＰＣＡベースの機能ブロックはこの目的で使用される。あるいは、．ＮＥＴ技術を有するビジュアルベーシック（ＶＢ）も、実装プラットフォームとして使用され得る。一つのケースでは、例示の実装方策はＶＢ．ＮＥＴ、ＤｅｌｔａＶ（登録商標）、Ｏｖａｔｉｏｎのいくつかの組み合わせを含み得る。上記にも拘わらず、実装は、異なるソフトウエアプラットフォーム（即ち、Ｊａｖａ（Ｒ），Ｄｅｌｐｈｉ，Ｃ＋＋）上でも達成され得る。図９−１３のオペレータディスプレイインターフェイスに関する更なる詳細は、以下に提供される。
【０１１９】
ＰＣＡの有用性にも拘わらず、ＰＣＡ技術は、プロセスの動的特性についての情報を容易に伝達しないということが注目される。これは、一般的に、ＰＣＡ分析の結果（スコア）を表示するために、散布図を使用したことに起因する。散布図は有用ではあるが、一般的には、プロセスデータの静的表現である。上記で論じたように、多次元の訓練又は監視データセットは、より低い次元の主成分空間上に射影される（スコアされる）。典型的に、最初の２つの主成分はデータの中の最大の変動を捕らえ、そして結果として、それらはデータの挙動に関して最も重要な情報により散布図を生成する。ＨｏｔｅｌｌｉｎｇのＴ2分布は、９５％又は９９％の信頼域を計算し、又はσ−ベースの境界を適用するのにも使用され得る。この領域の知識は、異常状態の検出にとって重要であり得る。何故なら、もし訓練データが正常なプロセスの性能を表しているなら、その領域内のどのような新たなスコアポイントに対しても、データポイントが正常であることは９５％（９９％）確かであり得るからである。これに対して、境界領域の外側のスコアポイントについては、データポイントが正常ではない（異常）ことは９５％（９９％）確かであり得るからである。
【０１２０】
上記から明らかなように、散布図からのスコアの記録と解釈は、データが不正確に解釈されることに関連する不正確さの程度を結果として生ずる。不幸にも、上記で論じたように、経験は、このアプローチがしばしば誤ったアラームに導くことを示している。この主な理由は、統計的に計算された境界は、プロセス変数の実際の物理的境界を考慮していないからである。各スコアは平均を中心とした入力変数の線型結合である。例えば、信頼限界領域の内側と外側で、与えられたスコアが正常又は異常であるという１００％の確実さは存在せず、従来のＰＣＡは異常状態の検出において必ずしも信頼できないという、５％（１％）の曖昧さが存在する。例えば、何も存在ないとき（誤ったアラーム）にＰＣＡスコアが異常状態を誘起し、そして反対にプロセスが実際に異常に振る舞っているときにＰＣＡスコアが正常な振る舞いを示し得る。異常状態であると決定し、異常アラームを発行するとき、これは一般的に許容できないことは明らかである。疑問なのは、物理的境界の同じ線型結合をとるときに何が起こるのかである。プロセス変数は異なるスケールを有しているので、これらの境界がどのようにＰＣＡスコアプロットに変換されるかは、必ずしも明らかではない。
【０１２１】
図９−１３のインターフェイスディスプレイは、オペレータインターフェイスが、先行する懸念に対処している間に、現在のプロセスデータに適用されるＰＣＡモデルを視覚化するのに使用されている、例示の方法を表している。手短に言えば、ＰＣＡ又は他の多変量解析から生成されたＰＣＡスコアプロット１４６と、単変量解析から生成されたＰＣＡスコアプロット１４８の両方が含まれる。これらの２つのプロットの組み合わせは、プロセスの現在の状態の指標を与える。図９の例では、ＰＣＡスコアは物理的境界内にあり、全てのプロセス変数がそれらの物理的境界内にあるので、このことは、プロセスが正常な動作状態にあることを示し、インディケータブロック１５０はその旨を強調している。
【０１２２】
図１０は、ＰＣＡスコアプロット１４６に従って、どのようにして全てが正常であるように見えるかの例を表している。しかしながら、物理的境界プロット１４８に示すように、プロセス変数の一つがその動作境界を超えている。従って、インディケータブロック１５２はアラーム見逃しを強調している。
【０１２３】
反対に、プロセス変数の全てがまだそれらの物理的動作境界内にある間、ＰＣＡスコアプロットが問題のないことを示すことも可能である。図１１はこのケースを表しており、インディケータブロック１５４によって示されているように、我々は偽アラームと呼んでいる。従って、真のアラーム状態は、ＰＣＡプロット上の統計的境界と物理的動作境界の両方が超えられたときに検出される。このことは、図１２においてインディケータブロック１５６で表されている。
【０１２４】
図１３を参照すれば、図示表現が、スライドバーとして示されている収集されたデータの先のセットに対するユーザ選択可能な図と、ドロップダウンメニューとして示されている主成分のユーザ選択可能な図とを伴っている。スライドバー又は他の図示表現は、ユーザが先に収集されたデータを先のプロセスの動作に基づいて選択し、先に収集されたデータに基づいてＰＣＡスコアを生成することを可能とする。この図示表現は、次に、ユーザが異常挙動に対する先のどのようなプロセス動作にもアクセスし解析することを容易に行えるように、ＰＣＡスコアで更新される。またさらに、プロセス変数は上記で論じたそれらの共通のプロセス境界を含むプロセスの境界の関数として表現され得ることは注目される。同様に、異なる主成分の選択は、ユーザが、プロセス上で異なる重要性を有する異なるプロセス変数の動作と、従って更新された図示表現とを見ることを可能とする。
【０１２５】
図１５−１９は、図９−１４の実施形態に関連して示したのと同様の方法で、制御境界のＰＣＡスコアとともに動的プロセスデータ視覚化を表す例示のユーザインターフェイスの異なる状態を示している。上述のように、７つのプロセス変数の制御境界は固有の値に正規化され、そしてこの場合に対するＰＣＡスコアが９５％信頼域の外側にあるなら、警告が発行される。これらの２つの図示は、入ってくるデータに対して同時にオンラインで実行される。加えて、ＰＣＡスコアプロットは、プロセス変数におけるピークはデータ入力に従う変化を監視している間、動的な一つのポイントを含んでいる。全てのピークが制御境界１６０の下にあるとき、プロセスは制御下にあり動作は正常である。対応する統計的スコアも、境界円１６２の内側にあるはずである。その状態は図１６に示されている。同様に、一つ又はそれ以上のピークが境界１６０を横切ったとき、制御不能又は異常状態が出現し、図１７に示すように、ＰＣＡスコアは境界円１６２の外側にあるはずである。
【０１２６】
偽の診断呼出もまた、経験され得る。図１８では、スコアは境界円１６２の外側にあり（異常）、一方、全てのプロセス変数は境界１６０（正常）の下にある。このケースは明らかに無視されるべき偽の警報である。しかしながら、その反対も起こり得る。図１９に示されているケースは、スコアが境界円１６２の内側にあり（正常）、一方、プロセス変数のインディケータはプロセス変数＃２が境界を破っている（異常）を示している。このケースは、検出の見逃し又は警報の見逃しである。
【０１２７】
２つのオンライン図示表現は、データポイントが正常又は異常のプロセス状態に相当しているかどうかを正確に診断することができる。更に、異常状態が発生したなら、正常状態を破っているプロセス変数がすぐに特定される。プロセスの動作に対するこれらの変数の重要性は、同じではない。特定の監視の要求に対するプロセス変数の重要性は、負荷行列ＰにおけるＰＣＡ負荷値から決定される。もし、警報を伝達しているプロセス変数が小さなＰＣＡ負荷値に関係しているなら、それはプロセス動作に対する非常に小さな重要性しか有しておらず、その警報は、単なる警告に分類され得る。これに対して、もし、負荷値が高ければ、プロセスの動作にとって変数は重要であり得、高い警報が忠告される。
【０１２８】
先の視覚化技術および添付のオペレータインターフェイスは、プロセス制御システムアーキテクチャ又はプラットフォーム内で実現され得、何れの場合にも、上述のＰＣＡブロックによって機能的に実装されることによってサポートされ得る。
【０１２９】
図２０および２１は、開示技術に従って生成されたさらにもう一つのオペレータインターフェイスを表しており、両方とも、上述のようにオンライン監視、診断および偽アラーム検出機能を有している。特に、図２０は、プロセス変数の境界のパーセンテージとして、バーグラフ内のプロセス変数を表し、上述のように共通のプロセス変数境界を含んでいる。
【０１３０】
多変量統計技術を用いた定常状態の動作の検出
【０１３１】
図２２および図２３は、ＰＣＡベースの統計的定数を使用した定常状態動作の認識と検出に関連する開示のもう一つの特徴を指向している。
【０１３２】
多くの異常状態防止アルゴリズムは、一つ又はそれ以上のプロセス変数の初期状態を学習し、次にプロセス変数がある様式で変化したときにアラームを誘起することによる検出に依存している。しかしながら、これらの異常状態防止アルゴリズムの多くは、初期条件のための訓練前および検出前の両方においてプロセスが定常状態にあることをも要求する。オペレータが多くの一つ又はそれ以上のプロセス変数を見、それらが定常状態であるかどうかを言うことは一般的には容易である。しかしながら、同様の決定を容易に行い得るアルゴリズムを生成することはより困難である。
【０１３３】
定常状態を検出することを指向する技術がここで開示されている。この技術は、複数のプロセス変数の監視に基づいて決定を行うことができる。
【０１３４】
プ過去においては、ロセス変数ｘが与えられると、平均と標準偏差は、与えられた長さ（例えば、５分）の非重複サンプリングウィンドウ上で計算され得る。平均（バーｘ）は以下の式により計算され：
【数５】
2010506256.tif 000006
標準偏差は以下の式により計算される：
【数６】
2010506256.tif 000007
ここで、ｎはサンプル数であり、ｘ1、ｘ2、．．．ｘnはサンプリングウィンドウ上でとったプロセス変数ｘのサンプルである。
【０１３５】
バーｘ1およびｓ1は一つのサンプリングウィンドウで計算されるそれぞれ平均および標準偏差とし、バーｘ2およびｓ2は、次のサンプリングウィンドウ上で計算されるそれぞれ平均および標準偏差とする。次に、もし、｜バーｘ1−バーｘ2｜≦３・ｓなら、プロセスは定常状態にあると記述される。反対に、もし、｜バーｘ1−バーｘ2｜＞３・ｓなら、プロセスは定常状態にはない。
【０１３６】
この計算で使用されたパラメータｓがどのように構成されるかが取り組まれる。もし、プロセスが定常状態にあり、サンプリングウィンドウが正しく選択されているなら、ｓ1≒ｓ2である。この場合には、ｓ＝ｓ1又はｓ＝ｓ2かどうかは問題ではない。しかしながら、定常状態検出のためのより小さい境界を結果的に生ずるより頑健なアルゴリズム（従って、異常状態防止モジュールが実行を開始したときに実際にプロセスが定常状態にあるというより大きい確実性を提供する）は、ｓ1およびｓ2の最小値又はｓ＝ｍｉｎ（ｓ1，ｓ2）を使用する。
【０１３７】
このような場合では、定常状態を宣言する基準は、｜バーｘ1−バーｘ2｜＞３・ｍｉｎ（ｓ1，ｓ2）である。
【０１３８】
もちろん、「３」という複数は、ユーザが設定可能なパラメータである。しかしながら、これは、特に多くのプロセス変数が存在する場合に、異常状態防止システムの構成を複雑にする。
【０１３９】
プロセス変数の平均と標準偏差を計算し、定常状態を決定する機能は、記載していない他の機能と同様に、上述のように統計的プロセス監視（ＳＰＭ）ブロックにカプセル化される。
【０１４０】
本開示のもう一つの特徴に従えば、定常状態検出技術は、ＰＣＡモデルによって生成されたスコア（プロセス変数データから直接というよりむしろ）を使用することを含み、複数のプロセス変数を有するプロセスが定常状態にあるかどうかを決定する。図２２は、ＰＣＡブロック１７２に結合したＳＰＭブロック１７０を表しており、複数のプロセス変数（ＩＮ1，ＩＮ2，ＩＮ3，．．．）を有するシステムが現在定常状態にあるかどうかを決定するのに使用される。ＰＣＡブロック１７２は、上述の機能ブロックに対応している。
【０１４１】
図２２に示されている例は、ＰＣＡモデルの第１のスコアを使用して定常状態の検出を実装する。ＰＣＡモデルに関して上記したように、最大の固有値の値に対応し第１の負荷ベクトルから計算される第１のスコアは、プロセスにおける変動の最大量を表している。例えば、いくつかのプロセスでは、これは、他のプロセス変数の多くを介して増殖する負荷の変化に対応している。このような場合には、ただ、第１の主成分が定常状態の基準に合致しさえすればよく、その場合はシステム全体が定常状態にあると決定され得る。
【０１４２】
第１の成分が最も統計的に重要である一方、別の実施形態は、図２３に示すように、ＰＣＡモデルの複数のスコアを使用し得る。この場合には、ＰＣＡ機能ブロック１７２がＳＰＭブロック１７０に結合され、各ＳＰＭブロック１７０はＰＣＡモデルのスコアの一つを監視する。この場合では、選択方法（例えば、Ｋａｉｓｅｒの規則、パーセント分散、平行分析など）に従って重要であると決定された各主要成分に対する一つのＳＰＭブロック１７０が存在するであろう。論理ブロック１７４は、システム全体についての定常状態の検出を行うのに使用され得る。もし、全てのＳＰＭブロック１７０が定常状態を示したら、定常状態はシステム全体について決定され得る。もし、ブロック１７０の何れも非定常プロセスを示したら、システム全体が定常状態にはないと決定され得る。
【０１４３】
この技術を実装する方法は、プロセスのＰＣＡモデルを生成し、第１のＰＣＡスコアを統計的プロセス監視（ＳＰＭ）ブロックに結合し、およびもしＳＰＭブロックが第１のＰＣＡスコアが定常状態にあることを宣言してシステムが定常状態であることを宣言する全てのプロセス変数からのプロセスデータを収集することを含み得る。ＳＰＭブロックが連続する非重複サンプリングウィンドウ上で入力変数の平均（バーｘ1およびバーｘ2）並びに標準偏差（ｓ1およびｓ2）を計算し、もし、｜バーｘ1−バーｘ2｜≦ｎ・ｍｉｎ（ｓ1，ｓ2）であるなら、ここでｎは実数であるが、入力変数は定常状態であることを宣言する。
【０１４４】
いくつかのケースでは、定常状態は、｜バーｘ1−バーｘ2｜≦ｎ・ｓ1によって決定される。あるいは、これに代えて、定常状態は、｜バーｘ1−バーｘ2｜≦ｎ・ｓ2によって決定される。例示の実施形態では、ｎ＝３である。
【０１４５】
上述の技術および方法は、第２およびより多くのＰＣＡスコアに結合した一つ又はそれ以上の追加のＳＰＭブロックと、ＳＰＭブロックのそれぞれの定常状態の検出を入力としてとり、もし全てのＳＰＭブロックが定常状態を示すならシステム全体が定常状態であることを示す出力信号を生成する論理モジュールとにより実装され得る。論理モジュールは、もしＳＰＭブロックの何れかが非定常プロセスであると決定したなら、システム全体が非定常プロセスであると追加的に宣言し得る。
【０１４６】
これに代えて又はこれに加えて、他の多変量技術が定常状態の検出のためのＳＰＭブロックに使用され得る。このような技術には、ＰＬＳ，ＰＣＲ，ＦＤＡ、ＣＶＡなどが含まれる。
【０１４７】
図２２および２３では、複数の別々のユニットを含むとして示されているけれども、開示されている技術および方法は、単一の、統合された機能ブロック（例えば、Ｆｉｅｌｄｂｕｓ機能ブロック）、フィールドデバイスインターフェイスモジュール、制御システム、単体のソフトウエア、アプリケーションにおいて実装され得る。
【０１４８】
過渡的ダイナミクスの検出
【０１４９】
図２７ー３７は、他の多変量技術を使用した過渡的ダイナミックスの検出に関連する本開示のもう一つの特徴を指向している。
【０１５０】
動的挙動は一般的にはプロセスの動作中に遭遇するものである。例えば、プロセスは、流速のような入力又は負荷変数の変化による過渡的な挙動を示す。これは加熱炉では一般的なことであり、なぜならそれらは通常、様々な原油の流量の下に動作しているからである。
【０１５１】
故障又は異常状態検出が行われているとき、正常動作状態と過渡的データとを見分けることは重要である。これは、プロセスは過渡的な領域では正常な方法で振る舞うことはなく、誤った診断が容易に構築されるからである。加熱炉では、原油の流量が増加したとき、原油の出力温度が上昇し、正常なレベルに下げるのに時間を要する。これは、流量における変化が操作された変数（燃料流量）における変化より早いからである。換言すれば、一次動的システムにおいて、新たな定常点に移行するには、最終温度は、時定数プラス時間のズレに等しくなるような時間を必要とする。それ故に、動的過渡的領域と定常状態又は設定点、故障若しくは異常状態が検出されたときの動作領域との間の見分けは、有用である。
【０１５２】
本開示のもう一つの特徴に従えば、過渡的動作が、例えば主成分分析（ＰＣＡ）に基づく多変量統計技術を使用して検出される。開示されている技術は、加熱炉における動的過渡に関連して記述されているが、開示の技術は、他の制御事情に関連して実装するのにも適している。また、この検出方法は、加熱炉における異常状態と動的過渡との間の違いを識別するために、ＰＣＡベースの技術を主成分回帰（ＰＣＲ）に統合することもできる。加熱管路内のコークス化（又は閉塞）の検出の例は、提案されている方法の頑健さを表すのに使用されるであろう。検出アルゴリズムにおける誤ったアラームを回避するために、プロセスの過渡的挙動を識別することは重要である。過渡的モードは誤った結果および誤ったアラームに導く。開示されている過渡検出技術は、蒸留カラム、熱交換などの他のプロセス又はユニットの動作における過渡モードに使用することができることも注目される。
【０１５３】
図２４に移って、加熱炉は精油所および石油化学プラントでは重油なユニットである。それらは原油又は重質炭化水素の温度をある程度の高温に上げるのに使用される。１８０で概略的に示されている加熱炉は、一つ又はそれ以上の管状通路１８２を含んでいる。供給路は等分されて管状通路１８２に別れている。各管状通路１８２は、管状通路１８２の内部を流れる供給物が燃焼する燃料によって加熱される熱交換器として機能している。燃料の流量は、通常、原油（供給の流れ）が所望の目標最終温度を達成するように操作される。
【０１５４】
加熱炉における一つの大きな問題は、コークス化である。加熱炉における高温により、原油は分解してコークスと呼ばれる残留炭素を形成する。コークスは時間と共に管の内部表面に堆積し、一般的に管又は加熱炉の性能を低下させる。このプロセスは閉塞と呼ばれている。時間が経過すると、加熱炉の性能は低くなり、加熱炉をきれいにするために、突然のシャットダウンが起こる。いくつかの例では、いくつかの管が重大な閉塞状態となる。加熱炉の管の内部のコークスの形成は、ユニットクリーニングをスケジュールするために監視されることが強く望まれる。
【０１５５】
加熱炉は、多変量監視および診断技術を使用する分析に適した他の多変量データ構造、および時折ＭＭＤＭ（即ち、多変量監視及び診断モジュール）ツールと称されるここに記載のツールを提供する。
【０１５６】
表１に説明するプロセス変数（ＰＶｓ）は、ここに記載の技術に従ってＭＭＤＭにより監視され処理され得る。これらのプロセス変数は、図２４にも示されている。プロセス変数は、上述のように、プロセス制御システム又はネットワークを介して所望の様式で、そしてＰｌａｎｔＷｅｂデジタルプラントアーキテクチャのような種々の異なるアーキテクチャおよびプラットフォームに従って、例えばＥｍｅｒｓｏｎのＦＯＵＮＤＡＴＩＯＮ（登録商標）  Ｆｉｅｌｄｂｕｓデバイス（例えば、Ｍｉｃｒｏ  Ｍｏｔｉｏｎのコリオリフローメーター、Ｒｏｓｅｍｏｕｎｔ温度デバイスなど）を使用することにより、供給され得ることが注目される。
【表１】
2010506256.tif 000008
【０１５７】
図２５は、主成分ｔ1およびｔ2によって定義される空間における訓練データスコアを表している。潜在する訓練データが種々の流量について、制御された変数の単一の設定点で集められる。約３０，０００のリアルタイムプロセスデータポイントが収集された。訓練データは、ＭＭＤＭツールのモデルを展開するのに使用され、ＭＭＤＭツールによって生成された散布図の一つが、図２５に示されているプロットと一致している。散布図の８つの集団は、訓練データにおける８つの異なる流量領域と一致している。
【０１５８】
図２５のプロットに示されている破線１９０の内側の領域は、図２６および２７の楕円に一致しおり、これは、同様に、９９％信頼性楕円を表す。つまり、楕円の内側の領域は、楕円の境界内のどのようなスコアも、加熱炉ユニットの正常動作に帰属する９９％の確率を与える。
【０１５９】
図２６および２７に移行して、管状通路（Ｔout）の出口温度は、本開示のこの特徴の従属変数として使用され、他の６つの変数は独立変数として取り扱われる。独立変数は、監視のために選択され、スコアプロットを有するＰＣＡモデルと同様に、正規化されたプロセス変数データモデルを生成するために使用される。プロセス変数データは、モデルの展開のための訓練データを提供する。このモデルは、次に、図２６に示すように、オンライン監視のためのオペレータユーザインターフェイスを定義するのに使用される。上述のように、このインターフェイスは、プロセスアラーム情報インディケータと同様に、独立プロセス変数を正規化およびスコアプロットの両方の様式で表示する。一般的に言えば、このインターフェイスは、図８に関連して上記したようにして生成され得、図９−２１に関連して記述した例示のインターフェイスの一つ又はそれ以上の何れかの特徴を包含し得る。
【０１６０】
従属プロセス変数（Ｔout）に関しては、我々は、訓練データを回帰タイプのＰＣＲモデルを生成するのに使用した。これは、モデル構築セッションの間にＭＭＤＭによって達成される。ＰＣＲモデルは、次に、オペレータインターフェイスウィンドウの「ＰＣＲ」タブを押し、これによって図２７に示すディスプレイウィンドウが現れることによりアクセスされ得る。ディスプレイウィンドウは、スコア成分Ｔ1，Ｔ2，．．．，Ｔnについて計算された回帰係数を示している。
【０１６１】
展開されたＰＣＲモデルを使用することによる過渡的な挙動の検出を有効にするために、検証データセット（４３００データポイント）が様々な原油流量について生成された。図２８は、検証データのための、原油流量に対する７つの異なる領域を示している。従って、６つの異なる過渡（動的）領域が検証データに存在し、ＰＣＲモデルは検証データからこれらの領域を予測するのに使用される。図２９は、予測されたＴout（Ｙ）と実際のＴoutを表している。図２９に示される重複するプロットは、従属変数（Ｙ）の定常状態および過渡的領域の両方における予測におけるＰＣＲモデルの成功を表している。過渡的領域における予測は、僅かに遅れており、期待とは一致していない。それにも拘わらず、検証データセットに存在する６つの過渡的領域は、図２９において明らかに特定可能である。
【０１６２】
図３０−３４は、Ｈｏｔｅｌｌｉｎｇ（Ｔ2）計算の使用と、過渡的領域検出に対するオペレータインターフェイスの対応する部分とを表している。特に、図３０は、過渡的領域が開始したときのＴ2プロットを示している。明らかに、Ｔ2は、原油の流量が変化し始めるに伴って急上昇している。加えて、図３０は、動作のポイントがスコアプロットｔ2−ｔ3の楕円の外にあることを示している。図３１は、たとえ独立変数の何れもが境界の外になくても、動的（過渡的）領域はＰＣＡによって検出されていることを示している。動的過渡的な進展と共に、Ｔ2は許容し得る境界の外に留まり、スコアプロットは楕円の外にある（図３２）。一度プロセスが定常状態の領域に落ち着くと、Ｔ2応答は境界の下にあり、スコアプロットは楕円領域の内側にある。以上のように、ＰＣＡスコアプロットとＴ2表現は、過渡の時定数を見積もるのと同様に、プロセス変数における過渡を検出するのに使用される。図３４は、全体の検証データについてＴ2により検出される６つの過渡的領域を示している。各過渡的領域の幅は、過渡の持続期間を表している。
【０１６３】
図３５および３６は、ＰＣＲ分析において制御された変数を使用しないアプローチの結果を表している。その際、制御された変数の組み込みの有用性と、それが提案されている方法が検出の頑健性にどのように影響するかとが記述される。この研究における加熱炉の例では、原油の最終温度（ＴＦ）は制御された変数である。原油を所望の設定値、最終温度まで加熱することが望まれている。既に述べたように、モデルは、同じ訓練データについての６つの独立変数のうちの５つを使用して展開され、原油の最終温度（ＴＦ）はリストから除かれる。６つの変数分析に使用したのと同様に、同じ検証データを使用することにより、新たなモデルが検証された。図３５は、新たなモデルはプロセスの動的挙動を検出することができないことを示している。
【０１６４】
Ｔ2も新たなモデルを使用して検証データについて得られた。図３６は、６つの変数モデルで得られた急上昇はもはや存在しないことを示している。
【０１６５】
先の記述に基づいて、この方法における制御された変数を含めることは、過渡的又は動的挙動を検出するのに有用である。このことは、制御された変数が負荷又は入力変数の変化に直接影響されるので、予想されることである。
【０１６６】
図３７は、コークス化の検出のための手法に向けられている。
【０１６７】
提示した議論から、過渡的期間の間、プロセスは動的に振る舞い、予測されるＹは実際のＹとは異なるであろうことは明らかである。従って、過渡的な間、コークス化検出は中止される。そうするための方法はＴ2の挙動を使用することである。検出アルゴリズムは、Ｔ2が境界の外にあるなら、コークス化検出を動作しないようにする。もし、Ｔ2が許容される境界内に戻ってきたら、コークス化検出は開始される。この手法は、図３７に示されている。開示されている検出方法および技術は、例えば、蒸留カラム、熱交換、タービン等の、また、化学、製油、電力、ライフサイエンスなどの産業の他のプロセスおよびユニット動作にも適用され得る。
【０１６８】
本開示の先の特徴に従って、ＰＣＡベースの監視技術がプロセス動作における過渡的挙動のオンライン検出に対して実装された。ＰＣＲ回帰技術も、プロセス動作における過渡的挙動のオンライン検出対して実装され得る。ＰＣＡおよびＰＣＲ技術の両方の組合せは、プロセスの警報の特徴付けをサポートし、開示されている技術は、異常状態防止タイプの警報とプロセス過渡によって引き起こされるものとを区別することができる。代替の設計が、多変量回帰、人工神経ネットワーク、およびベイズ意思決定木などの技術を使用することにより可能である。
【０１６９】
例示のアプリケーション：加熱炉におけるＰＣＡおよびＰＣＲを用いたコークス検出
【０１７０】
図３１−４１に移行して、本開示の他の特徴が、加熱炉管状通路におけるコークス化（又は閉塞）を検出するＰＣＡおよびＰＣＲ技術に向けられている。以下に記述される実施形態では、ＰＣＡおよびＰＣＲ技術が、オンラインデータを収集しおよび分析し、並びにオンラインでコークス化を検出するためにプロセスを監視するのに使用される。開示されている方法および技術は、どのようなコントローラ、フィールドデバイス、およびここに記述されている他のプロセス制御システムの構成要素とも統合され得る。
【０１７１】
以下に記述するように、開示されているコークス化検出方法は、一般的に３つのステップを含んでおり、即ち：訓練データを集めること、ＰＣＡおよびＰＣＲ技術を使用してモデルを展開すること、および目標となるプロセス変数の実際および予測された値を使用してコークス化障害を検出すること、である。
【０１７２】
収集された訓練データは、一般的にはリアルタイムのデータであり、これは故障のない、そして負荷変数の期待される動作範囲を表している。例えば、訓練データは、動作している加熱炉で遭遇されると期待される原油流量の最小および最大を含み得る。展開されるモデルは、最小および最大流量の間の幾らかの流量を有していることから恩恵を受け得る。いくつかのケースでは、訓練データは、負荷変数の各範囲をカバーするように収集され、変化する負荷変数に対する過渡的（動的）データをも含み得る。一般的には、制御された変数に対するサンプリング時間および設定点は、訓練データのセットを通じて定数のままであり得る。
【０１７３】
先の特性を有するデータが与えられると、訓練データは種々の方法で収集され得る。例えば、データはオンラインで又はプロセスデータヒストリアンから抽出され得る。何れの場合にも、開示されている方法の実践は、訓練データが収集される方法に限定されない。
【０１７４】
複数のプロセス変数が、管状通路における頑健なコークス化検出をサポートするのに使用される。この例のケースでは、開示されている方法は、表１に述べられ図２４に示されているプロセス変数のセットを使用する。これらのうち、管状通路の出口温度（Ｔout）が監視された従属プロセス変数であり、管状通路内でコークス形成に伴って減少する。コークス化ヒーターについては、管の内側にコークス（故障）が付着し始めるにつれて、高温側から低温側へ管を横切る全体の熱移動が低下する。正常動作（故障のない）に比較して、管状通路の出口温度（Ｔout）もまた減少する。表１の他の６つの変数は、独立であると考えられる。
【０１７５】
加熱炉は、上述の多変量監視及び診断モジュール（ＭＭＤＭ）を使用して、その分析を説明するのに適した多変量データ構造を提供する。監視のために一度独立変数が選択されると、選択された変数は、スコアプロットを有するＰＣＡモデルと同様に、正規化されたプロセス変数データモデルを生成するのに使用される。収集されたプロセス変数データは、モデルの展開のための訓練データを提供する。展開されたモデルは、次に、上述のように、オンライン監視のためのオペレータインターフェイスを生成するのに使用される。例えば、図２５に示すように、６つの独立変数が、４つの異常アラーム状態インディケータに関連して、正規化およびスコアプロットの両方の様式で表示される。このモデルおよびオペレータインターフェイスの詳細は、上記で述べられている。
【０１７６】
従属プロセス変数（Ｔout）に関しては、訓練データは回帰タイプのＰＣＲモデルを生成するのにも用いられる。上述のように、ＰＣＲモデルの生成は、モデル構築手法の実装の間に、ＭＭＤＭツールによって自動的に達成される。ＰＣＲモデルは、オペレータインターフェイスにおけるＰＣＲタブを選択することによってアクセスされ、図２７に示す例示の表示ウィンドウが生成される。上述のように、このディスプレイウィンドウは、スコア成分Ｔ1，Ｔ2，．．．，Ｔnについて計算された回帰係数を示している。
【０１７７】
開示されているコークス検出に適した加熱炉監視システムおよび方法の実施形態が、記述される。
【０１７８】
Ｔout（Ｙ）の実際の値が予測される値ハットＴout（ハットＹout）より一貫して低ければ、コークスの形成が示唆される。このことは、実際の管出口の温度が一貫して予測された温度より低いことを意味している。このケースでは、故障アラーム又は警報の何れも生成されない。Ｙの絶対値又は平均値は、この目的のために使用され得る。以下にさらに示すように、予測されたハットＹからの実際のＹの偏差は、コークス化の重大さの優れた示唆を与える。これは、実際のユニット性能に基づいてユニット洗浄又はシャットダウンを計画するために、プラント技術者又はオペレータによって使用され得る優れたツールである。比Ｙ／ハットＹが時間の関数としてプロットされることが推奨される。この方法により、プラント技術者又はオペレータがユニットの性能を継続的に監視しおよびユニットの低下又はコークス化を検出する。これに代えて、２つの値の差が時間の関数として監視され得る。ハットＹに対するＹの比に基づいて、プラント要員はどのような行動を採るべきかを決定することができる。コークス化は遅いプロセスであるため、プラント要員はアラーム又は警報を実行するように選択しないかも知れない。しかしながら、警報又はアラームを閾値の境界に基づいて活性化することは容易である。例えば、その比がある境界（λ）まで下がったら、警報が発行される。
【０１７９】
例えば、もし、
【数７】
2010506256.tif 000009
なら、警報が発行される。
【０１８０】
種々の警報が、コークス化の重大さを示すために発行される。例えば、λ1は低度のコークス化を表し、λ2は中程度のコークス化を示すなどである。
【０１８１】
より一般的には、開示されているコークス検出方法は、以下のステップを含む。即ち：十分な訓練データを得ること、統計的モデルを展開すること、およびユニット性能を監視するためのモデルを使用すること、である。
【０１８２】
開示されている技術のための訓練および検証
加熱炉に対する訓練および検証データは、ＧＳＥ  Ｓｙｓｔｅｍｓから商業的に入手可能な高忠実性シュミレータＳＳＰＲＯにより、又はこれから得られた。サンプリングレートは１秒である。燃料の流量は、原油の目標最終温度が達成されるように操作された。シュミレーションは、種々の重大さのレベルのコークス化を導入し得るように柔軟性を有していた。シュミレーションは、訓練データおよびシュミレーションデータを得るために使用された。７つのプロセス変数が記録され、ＭＭＤＭツールによって提供される多変量解析において使用された。
【０１８３】
訓練は、種々の流量で実施された。原油の流量と入口温度により、現実のプラントの動作を開始するのに関して、一律のノイズを有することが可能となった。約３０，０００のリアルタイムのプロセスデータポイントが収集された。訓練データはＭＭＤＭモデルを展開するのに使用され、ＭＭＤＭツールによって生成された散布図は、図２５に示されており、８つの流れが訓練データに留まっていることに相当する８つの群を表している。図２５は、９９％信頼性の楕円を表示している。楕円の内部の領域は、楕円の境界内の何れのスコアも、加熱炉の正常動作に帰属するという９９％の信頼性を与える。
【０１８４】
提案されている異常状態検出の頑健さをテストするために、いくつかのリアルタイムシュミレーションデータファイルが種々のコークス化レベルで収集された。シュミレーションデータポイントは、良いデータと誤ったデータとの両方を含んでいる。コークス化を管状通路に導入するために、全体の熱伝導係数（Ｕ）が故意に下げられた。いくつかのデータセットが以下のように収集された。セット１では、係数Ｕは１０００から９９５に、次に９９０に、９８０に、９５０に下げられた。セット２では、係数Ｕは１０００から５００まで下げられた。
【０１８５】
この例示の適用では、この特定の場合に対する展開されたＰＣＲモデルは、Ｙ＝Ｂ0＋Ｂ1＊Ｔ1＋Ｂ2＊Ｔ2＋Ｂ3＊Ｔ3，ここで、Ｂ0＝３１６．９２２７８９５，Ｂ1＝−０．２１５１５４４８５，Ｂ2＝０．２７４１２０９８０，およびＢ3＝−０．２７６０８３０３９である。Ｔ1，Ｔ2，Ｔ3は、変数の寄与に基づくオンライン（検証データ）で決定されたスコア成分である。
【０１８６】
障害（コークス化）の検出のために新たな監視にこのモデルを適用する前に、検証データの新たなセットを使用することにより、このモデルがテストされた。このステップの目的は、この方法の頑健さを確実にすることである。図３８は、Ｙ（Ｔout）の予測においてこのＰＣＲモデルがどれくらい正確かを示している。このモデルは、９９．９％より高い正確さを有しており、データが種々の動作条件の中で散乱している（図２５）という優れた正確性である。もし、動作条件の範囲が狭ければ、正確さはさらに高くなるであろう。
【０１８７】
また、このモデルは、訓練データで使用されているより１２％高い流量についても、Ｙ（Ｔout）の予測に使用された。換言すれば、このモデルは、訓練データで使用された流量の範囲の外の流量に対しても、プロセスの挙動を予測するのに使用された。図３９は、予測は高い正確性を以てここでも動作した。それにもかかわらず、どのモデルも、訓練範囲の外側のデータに対して使用されることはあまり推奨されることではない。
【０１８８】
上述のモデルは、検証データにおけるコークス化の検出に使用された。このモデルは、全てのシュミレーションされた検証データについてコークス化を検出することができる。図４０は、実際のＹが予測されたＹよりどれくらい低いかを示している。Ｙに対する実際のＹの比は、コークス化が増大する（全体の熱伝導が低下する）に伴って低くなる。予測されたＹは、正常な状態（故障／コークス化のない）でのＹの値である。予測されたＹに関する実際のＹの挙動から、プラント要員は管状通路内のコークス化の重大性を決定することができる。
【０１８９】
データの他のセットが提案されている検出の方法論を検証するために使用される。このセットは、高い値の突然の大きな異常状態（コークス化）が導入された良好なデータを含んでいる。図４１は、予測されたＹに対する実際のＹの比がどのように変化するかを示している。図４１の最初の部分は、良好なモデルが正常の挙動をどのように予測するかを示している。
【０１９０】
要約すれば、上述の技術（および内在するモデル）は、ＰＣＡ／ＰＣＲを使用して加熱炉管状通路におけるコークス、閉塞および他の障害状態を検出することを指向している。そのためには、原油の入口温度、加熱炉管状通路の原油流量、トータル原油流量、管状通路出口温度、最終温度、燃料流量、および燃料出口温度が監視され、管状通路出口温度が加熱炉の性能インディケータとして使用され得る。いくつかの例では、性能インディケータは、ＰＣＡを用いてモデル化され、そのモデルについての訓練データの範囲内で動作条件に対するＰＣＲを使用することによって予測される。それにもかかわらず、ＰＣＲは、訓練データの範囲外での動作条件に対する加熱炉のキーとなる性能インディケータの予測に使用され得る。
【０１９１】
図４２および４３において取り組まれている開示の特徴に移行すれば、加熱炉におけるコークス化の多変量監視及び診断に対するＰＣＡパラメータベースの技術が開示されている。このアプローチは、複数の主成分分析（ＰＣＡ）モデルの展開と、異なるプロセス動作条件に対するパラメータの比較とを基礎としている。コークス化の検出はさておき、開示されている技術、方法は、種々の異なるプロセス装置（例えば、加熱炉におけるススの堆積、反応装置、コンプレッサ、熱交換器、蒸留カラム、膜などの閉塞および他の故障、）を巻き込んだ種々のプロセスの遅い変化の決定又は検出に対して有効である。
【０１９２】
上述のように、ＰＣＡベースの技術は、オンライン監視及び診断に関連して実装するための分散型プロセス制御システム（ＤＣＳ）に統合され得る。また、開示されている技術は、オンラインモデルの展開およびモデルパラメータの比較のためのデータの生成と操作を可能とする。また、ＰＣＡベースの技術は、主成分回帰（ＰＣＲ）分析を使用して加熱炉におけるコークス化レベルを決定するという挑戦によって表される具体的な詳細に適用される。
【０１９３】
本開示のこの特徴は、コークス化のような遅い種々のプロセス変化を、回帰ベースのアプローチを使用することなく検出するための方法を提供している。開示されている方法は、その代わりに、プロセス動作の間に連続して繰り返し又はさもなければ絶えずＰＣＡモデルを展開することと、および展開されたモデルの正常なプロセス動作状態の下で展開されたベースラインモデルに対して比較することとに基づいている。開示されている技術は、以下に述べる理由により回帰を使用する他のものに比較して、特定のアプリケーション（例えば、異常状態防止アプリケーション）におけるより頑健な解を提供し得る。一つの主たる理由は、回帰ベースの技術は、データポイントがフィールドデバイス通信リンク（例えば、Ｆｉｅｌｄｂｕｓ、Ｍｏｄｂｕｓ、無線など）から読み取られる毎に予測を生成するからである。結果として、回帰ベースの技術は、プロセス状態の現在の予測であり、これは次に、診断状態を決定するためにロジックを介して処理され又は評価される。このプロセスは、データ列における異常ノイズと同様に、設定点および定常状態変化のような動的プロセスによる誤った警報／アラームをしばしば生成する。このタイプの警報は、オペレータにとって面倒であり、従って、もし全く除去されないなら最小化されなければならない。
【０１９４】
これに対して、開示されているアプローチは、連続又はオンデマンドの何れかで、プロセス動作のバックグラウンドにおいて予め決められた次元のデータファイルを生成する。これらのデータファイルは、次に、プロセスの性能を評価するために使用されるＰＣＡモデル最新情報を生成するのに使用される。結果として、開示されている方法は、現在のプロセスデータの値には全く依存しないので、より頑健である。むしろ、開示されている方法は、例えば、データファイルのサイズによって決められる長期のデータ挙動に基づくプロセス性能診断に基づいて警報レポートを生成し得る。開示されている方法の更なる詳細は、以下に記述される。
【０１９５】
最初に、ＰＣＡを展開するために訓練データファイルから始められる。訓練データは、正常なプロセス動作条件の間に収集されたデータを表している。上記のように、データファイルはデータの行（観測値）と列とを含み、それぞれの列は配置されたプロセス変数（ＰＶ）に対応している。プロセス変数は、行データ又は行データから得られる統計的な符号定数（即ち、異常状態防止ブロックデータ）の何れかである。従って、データファイルは、ｍ＞＞ｎという通常の仮定を伴うｎ列×ｍ行を含んでいる。このような配置で、このモデルは、ｍ個の観測値を伴うｎプロセス変数に対処し、ｍ×ｎのデータ行列Ｘを含み、これは上述のように自動スケールされている。実際のモデルは、次にｎ×ｎの共分散行列Ｓのスペクトル成分を計算することにより構築され、ここでＳ＝ＸT・Ｘ／（ｍ−１）であり、これはＳをＳ＝Ｖ・Ｄ・ＶTと書き換え、ここで、
【数８】
2010506256.tif 000010
であり、ｓi2，ｉ＝１，．．，ｎは、最大から最小に並び替えられた固有値（分散）であり、Ｖの列は対応する固有ベクトルである。最大の固有値およびその対応する固有ベクトルは、オリジナルのプロセスデータにおける最大の変動量を含むオリジナルの入力空間における方向に相当している。これは、第１の主成分であると考えられる。第２番目に大きい固有値とその対応する固有ベクトルは、第１の主成分に直角の方向に相当し、次に大きい変動量を含んでいる。これが、全てのオリジナルのプロセスデータを説明する新たな直交ベクトル空間が生成するまで続けられる。
【０１９６】
固有値および固有ベクトルが計算された後、何れが重要な主成分であるかを決定することが必要であり、それは、通常、プロセスにおける重要な変動に対応し、重要でないものはノイズに対応している。重要および重要でない主要成分を決定するために、この問題を取り扱う多くの方法が存在する。
【０１９７】
どの方法が選択されるにしても、最終の結果は、重要な主要成分ａの数の決定であり、それはＰＣＡモデルについて保持されるであろう。次に、負荷行列Ｐ∈Ｒn×aがＶから最初のａ列（固有ベクトル）をとることにより生成される。オリジナルの観測値の負荷行列による部分空間への射影は、スコア行列と称され、Ｔ＝Ｘ・Ｐと表記される。Ｔの列はＰＣＡモデルのスコアと称され、そしてｔ1，ｔ2，．．．ｔnとして言及され、ｎ×ａ行列ＰはＰＣＡモデルとして言及される。
【０１９８】
訓練データファイルＸ。本開示のこの特徴の目的のために、プロセスの挙動が決定されるのに影響を与えＸのうちの最初のｋ＜ｎを表す独立のプロセス変数を選択し、同様にプロセスの挙動を記述しＸのうちの最後のｎ−ｋ＞０によって表される従属変数を選択することにより、訓練データファイルＸが選択される。もし、回帰技術が実装されると、Ｘの中の最初のｋ個のプロセス変数が、最後のｎ−ｋ個のプロセス変数の予測に使用される。しかしながら、ＰＣＡモデルを展開するとき、Ｘはｋ個の独立プロセス変数だけを有し、そしてｎ＝ｋであることに注意すべきである。それらは、次に従属変数を予測するのに採用されるスコアｔ1，ｔ2，．．．ｔaを生成するのに使用される主成分ベクトル（負荷ベクトル）を有する負荷行列Ｐを生成するのに使用される。開示されているアプローチでは、独立および従属変数は組み合わせられて訓練データファイルＸを生成する。
【０１９９】
一度ＰＣＡモデルが特定されると、負荷行列Ｐ＝［ｐ1，ｐ2，．．．，ｐa］が定義されるようになり、ここでｐi＝［ｐi,1，．．，ｐi,k，ｐi,k+1，．．．，ｐi,n］Tはｉ＝１，．．，ａに対する負荷ベクトルである。ｐiの最初のｐi,1，．．，ｐi,k成分は、Ｘにおける最初のｋ個の独立変数に対応する負荷値である。同様に、ｐi,k+1，．．，ｐi,n成分は、Ｘにおける最後のｎ−ｋ個の従属変数に相当する負荷値である。結果として、Ｐはオリジナルのプロセス変数に直接関連づけられた相関のある負荷値の空間を定義する。Ｘにおける各従属変数、それらのｎ−ｋ個は、以下の式よって与えられるこの空間における負荷値距離（metric）を有し：
【数９】
2010506256.tif 000011
ここで、 ｊ＝ｋ＋１，…，ｎであり、σiはＤにおけるｉ番目の固有値である。
【０２００】
負荷値距離（metric）ｚjは以下の性質を有するｊ番目の従属プロセス変数のモデルの結果（outcome）である。ここで、ＸおよびＹは、同じプロセス動作から異なる時間に集められたデータファイルであり、ｚjおよびｖjは同じ従属プロセス変数のＸおよびＹにそれぞれ対応するモデルの結果（outcome）であり、もし、ｚj≒Ｖjならｊ番目のプロセス変数によって記述されるプロセスの挙動は、ＸおよびＹに対して同様である。これとは逆に、もし、ｚj≠ｖjなら、ｊ番目のプロセス変数によって記述されるプロセスの挙動は、ＸおよびＹに対して異なる。
【０２０１】
この性質はプロセスの診断には有用であり、特に変化の遅いするプロセス挙動の診断には有用である。その理由は、負荷値距離（metric）ｚj はオペレータによって特定される時間期間に亘るプロセスの挙動を反映しており、従って、回帰アプローチによって与えられる最も最近のアップデートよりも完全な診断を提供するからである。
【０２０２】
加熱炉アプリケーション。先の技術は、加熱炉におけるコークス化のレベルを決定するのに使用され得る。上述のように、データファイルＸは表１および図２４によって定義されるプロセス変数により構成されている。その実施例では、Ｔoutは一意的な従属変数であり、一方、他の全ての変数は独立である。しかしながら、本開示のこの特徴に従う技術では、全ての変数がＰＣＡモデルの生成に使用され、ＰＣＲは必要ではない。
【０２０３】
このような割り当てで、ベースラインモデルが正常動作の条件（即ち、０％コークス化）の下で生成され、負荷値距離（metric）ＺがＴoutについて計算された。続いて、同様のことが、加熱炉におけるコークス化の異なる％に対して行われ、図４２は得られた結果を示している。図４２から明らかなように、コークス化の％が上昇するに伴うＺにおける単調な減少が存在し、これは、コークス化の％が変化することにより、負荷空間における構造的な変化が存在することを示唆している。変化は合作的であり、その時に一つの変数を監視することにより容易に決定し得ないことは注目される。
【０２０４】
上記に基づいて、Ｚ上の異なる境界を置くことにより、診断／故障検出方法又はシステムにおいて使用され得る。各モデルはＺの挙動をモデルデータを集めるのに使用される時間の長さに基づいて表現するので、この技術を使用して発行される警報およびアラームは、より信頼できる。さらに、実行の間、この検出方法は、オペレータに対して完全にトランスペアレントである。これは、この検出方法が監視コンソールのバックグラウンドで実行され、設定された境界が超えられ又は違反されたときにのみ異常状態のオペレータに警報を出すように構成され得るからである。図４３のフローチャートは、この方法およびシステムが一実施形態に従ってどのようにして実装されるかを示している。このアプローチで、誤ったアラームおよび警報が避けられ、にもかかわらずプロセスパラメータの性能が経時的に厳格に監視され得ることは注目される。
【０２０５】
本開示のこの特徴により、ＰＣＡモデルパラメータは、オンライン診断について定義されおよび適用され、加熱炉におけるコークス化および種々のプロセス装置の故障および異常状態に有用である。このモデルパラメータは、ＰＣＡ負荷および分散情報を使用して導かれ得、負荷値距離（metric）は従属又は独立のプロセス変数について定義され得る。いくつかのケースでは、開示されている方法は、瞬時の又は最新の変化よりもむしろ長期間のコークス化ヒーターの観察に使用され得る。例えば開示されている方法は、オンラインの長期に亘る協調的な診断に使用され得る。これに代えて又はこれに加えて、開示されている方法は、回帰分析へのアプローチの代替を提供し得る。開示されている方法は、例えば、ＤｅｌｔａＶ（登録商標）、Ｏｖａｔｉｏｎなどを含む多くの制御システムプラットフォームに関連して、そしてＲｏｓｅｍｏｕｎｔ３４２０ＦＦインターフェイスモジュールのような種々のプロセス装置およびデバイスとともに実装され得る。あるいは、開示されている方法およびシステムは、単独型の異常状態検出アプリケーションとして実装され得る。何れの場合も、開示されている方法およびシステムは、警報を生成するように構成され得、さもなければ加熱炉におけるコークス化のレベルの基準をサポートし得る。
【０２０６】
当業者は、上記の例示のシステムおよび方法が種々の方法で改変され得ることを認識するであろう。例えば、ブロックは省略し又は並び替えられ、追加のブロックが加えられ得る。例えば、図７に関しては、ブロック１４６はフローの異なる点に実装され得る。同様に、ブロック１４８は、別のルーチンの一部として実装され得、従って、それは、適切なコマンドが受信されて別のルーチンの実装を開始するときに依存して、図７のフローの種々のポイントに実際に出現し得る。
【０２０７】
異常状態防止モジュールおよび異常状態防止ブロックを含む上記の例は、開示されているシステム、方法、および技術の実践がこのような状況に限定されないということの理解と共に開示されている。むしろ、開示されているシステム、方法、および技術は、異なる構成の構造、コンポーネント配置、又は他の別々の部品、ユニット、構成要素、若しくは監視、データ集合の選択が可能な事項の集合など、どのような診断システム、アプリケーション、ルーチン、技術又は手法とも使用するのに適している。診断において使用されているプロセスパラメータを特定する他の診断システム、アプリケーションなどは、ここに記載されているシステム、方法および技術から開発されさもなければ恩恵を受け得る。パラメータの特定の仕様は、ここに関連するプロセスデータを配置し、監視しおよび格納するに使用され得る。更に、開示されているシステム、方法および技術は、プロセス制御システムの診断の特徴に関連して単独で使用される必要はなく、特に、このような特徴が開発段階又は開発の早期の段階にあるときは、単独で使用される必要はない。むしろ、開示されているシステム、方法、および技術は、プロセス制御システム、プロセスプラント又はプロセス制御ネットワークなどとの使用に適している。
【０２０８】
ここに開示されている方法、プロセス、手法および技術は、ハードウエア、ファームウエア、およびソフトウエアのどのような組合せをも使用して実装され得る。従って、ここに記載されているシステムおよび技術は、標準の多目的プロセッサにおいて又は記載したように特別に設計されたハードウエア又はファームウエアを使用して実装され得る。ソフトウエアに実装されたとき、そのソフトウエアは、磁気ディスク、レーザディスク、若しくは他の媒体、コンピュータ、プロセッサ、Ｉ／Ｏデバイス、フィールドデバイス、インターフェイスデバイスなどのラム若しくはロム若しくはフラッシュメモリなどに格納され得る。同様に、そのソフトウエアは、例えば、コンピュータ読み取り可能なディスク又は他の可搬性のコンピュータ記憶メカニズム又は通信媒体を含む既知の又は所望の配送方法を介してユーザ又はプロセス制御システムに配送される。通信媒体は、コンピュータ読み取り可能な命令、データ構造、プログラムモジュール、又は搬送波若しくは他の搬送メカニズムなのような変調されたデータ信号における他のデータを具体化する。「変調されたデータ信号」とは、一つ又はそれ以上のその特性のセットを有する信号又はその信号のエンコード情報に変更された信号をいう。限定されない一例として、通信媒体は、有線ネットワーク又は直接有線接続のような有線媒体、音響、ラジオ周波数、赤外のようなワイヤレス線媒体および他の無線媒体を含む。従って、そのソフトウエアは、電話回線、インターネットなど（これらは、可搬性の記憶媒体を介してこのようなソフトウエアを提供するものと同じ又は互いに交換可能であるように見える）の通信チャンネルを介してユーザ又はプロセス制御システムに配送され得る。
【０２０９】
本発明は、特定の実施形態を参照して記述されているけれども、これらは例示のみを意図し本発明を限定するものではなく、当業者には、開示されている実施例に、本発明の精神および範囲から逸脱することなく、変更、付加、削除が為され得ることが明らかであろう。
【符号の説明】
【０２１０】
１０  プロセスプラント
１２  プロセス制御システム
１２Ａ  オペレータインターフェイス
１２Ｂ  コントローラ
１２Ｃ  入力／出力カード
１４  プロセス制御システム
１４Ａ  オペレータインターフェイス
１４Ｂ  分散コントローラ
１５  フィールドデバイス
１６  フィールドデバイス
１８  保守コンピュータ
１９  ワークステーションツール
２０  回転装置
２２  保守コンピュータ
２３  アプリケーション
２４  電力発生分配システム
２５  電力発生分配装置
２６  支援コンピュータ
２７  診断アプリケーション
３０  コンピュータシステム
３５  異常状態防止システム
３８  構成アプリケーション
４２  異常動作検出システム
４３  警報／警告システム
４５  有線接続バス
４６  インターネット接続
５４  プロセス制御システム
６０  プロセスコントローラ
６４  フィールドデバイス
６６  フィールドデバイス
６８  Ｉ／Ｏデバイス
７０  Ｉ／Ｏデバイス
７２  ワークステーション
７４  保守ワークステーション
７４Ａ  プロセッサ
７４Ｂ  メモリ
７４Ｃ  ディスプレイデバイス
７６  通信バス
７８  データベース
【図１】
2010506256.tif 000012
【図２】
2010506256.tif 000013
【図３】
2010506256.tif 000014
【図４】
2010506256.tif 000015
【図５】
2010506256.tif 000016
【図６】
2010506256.tif 000017
【図７】
2010506256.tif 000018
【図８】
2010506256.tif 000019
【図９】
2010506256.tif 000020
【図１０】
2010506256.tif 000021
【図１１】
2010506256.tif 000022
【図１２】
2010506256.tif 000023
【図１３】
2010506256.tif 000024
【図１４】
2010506256.tif 000025
【図１５】
2010506256.tif 000026
【図１６】
2010506256.tif 000027
【図１７】
2010506256.tif 000028
【図１８】
2010506256.tif 000029
【図１９】
2010506256.tif 000030
【図２０】
2010506256.tif 000031
【図２１】
2010506256.tif 000032
【図２２】
2010506256.tif 000033
【図２３】
2010506256.tif 000034
【図２４】
2010506256.tif 000035
【図２５】
2010506256.tif 000036
【図２６】
2010506256.tif 000037
【図２７】
2010506256.tif 000038
【図２８】
2010506256.tif 000039
【図２９】
2010506256.tif 000040
【図３０】
2010506256.tif 000041
【図３１】
2010506256.tif 000042
【図３２】
2010506256.tif 000043
【図３３】
2010506256.tif 000044
【図３４】
2010506256.tif 000045
【図３５】
2010506256.tif 000046
【図３６】
2010506256.tif 000047
【図３７】
2010506256.tif 000048
【図３８】
2010506256.tif 000049
【図３９】
2010506256.tif 000050
【図４０】
2010506256.tif 000051
【図４１】
2010506256.tif 000052
【図４２】
2010506256.tif 000053
【図４３】
2010506256.tif 000054
【国際調査報告】
2010506256.tif 000055
2010506256.tif 000056
2010506256.tif 000057
2010506256.tif 000058
2010506256.tif 000059