(57)【要約】  本発明は、流動接触分解ユニットのプロセス・ユニットにおける異常事象を検知する方法である。本方法は、プロセス・ユニットの動作を統計および工学モデルと比較する。統計モデルは、これらの装置の正常動作の主成分分析により構築される。また、当該工学モデルは、部分最小二乗法分析および変数間の相関分析に基づいている。プロセス・ユニットと正常なモデル結果の動作の差違が異常状態を示す場合、異常状態の原因を特定して是正する。【選択図】図１
2008512800.tif 000002
【特許請求の範囲】
【請求項１】
流動接触分解ユニット（ＦＣＣＵ）のいくつかのプロセス・ユニットにおける異常事象検知（ＡＥＤ）を行なう方法であって、
（ａ）前記プロセス・ユニットからのオンライン測定値を、対応するプロセス・ユニットの正常動作のモデルの組と比較する工程；
（ｂ）現在の動作が、期待される正常動作とは異なっていて、プロセス・ユニットにおける異常状態の存在を示すか否かを判定する工程；
（ｃ）前記プロセスのオペレータが前記ＦＣＣＵにおける異常状態の根本原因を特定することを支援する工程；および
（ｄ）前記ユニットを正常動作に戻すべく是正処置を実行する工程
を含むことを特徴とする方法。
【請求項２】
前記モデルの組は、設備グループおよびオペレーティング・モードに対応し、１個以上のオペレーティング・モードを含みうる各グループ毎に１個のモデルが対応することを特徴とする請求項１に記載の方法。
【請求項３】
前記モデルの組は、設備グループおよびオペレーティング・モードに対応し、各グループおよび各モードに１個のモデルが対応することを特徴とする請求項１に記載の方法。
【請求項４】
前記設備グループは、同一グループ内に全ての主要な物質およびエネルギー相互作用を含むことを特徴とする請求項２に記載の方法。
【請求項５】
前記設備グループは、同一グループ内に高速リサイクルを含むことを特徴とする請求項４に記載の方法。
【請求項６】
正常動作の前記モデルの組は、主成分モデルを含むことを特徴とする請求項５に記載の方法。
【請求項７】
正常動作の前記モデルの組は、工学モデルを含むことを特徴とする請求項６に記載の方法。
【請求項８】
各プロセス・ユニットにおける正常動作の前記モデルの組は、主成分モデルまたは工学モデルのいずれかであることを特徴とする請求項１に記載の方法。
【請求項９】
各プロセス・ユニットにおける正常動作の前記モデルは、主成分分析（ＰＣＡ）、部分最小二乗法に基づく推論および相関に基づく工学モデルを用いて決定されることを特徴とする請求項１に記載の方法。
【請求項１０】
前記流動接触分解ユニットおよび下流塔は、１２個の異常モニターに分割されていることを特徴とする請求項９に記載の方法。
【請求項１１】
前記プロセス・ユニットは、ＦＣＣＵシステムの動作区画に分割されていることを特徴とする請求項１に記載の方法。
【請求項１２】
１０個の動作区画が存在することを特徴とする請求項１１に記載の方法。
【請求項１３】
前記１０個の動作区画は、反応炉−再生装置、ライトエンド塔、触媒循環、スタックバルブ、サイクロン、送風機、炭素収支、主精留装置への触媒キャリーオーバー、湿式ガス圧縮機、バルブ流量モデルを含むことを特徴とする請求項１１に記載の方法。
【請求項１４】
前記異常事象モニターは各々、当該領域における異常状態の確率を示す連続信号を生成することを特徴とする請求項９に記載の方法。
【請求項１５】
前記モデルは、センサーにより測定されたプロセス変数値を含むことを特徴とする請求項９に記載の方法。
【請求項１６】
異なるプロセス・ユニットに対するモデルの前記主成分は、同一センサーで測定された、いくつかのプロセス変数値を含むことを特徴とする請求項９に記載の方法。
【請求項１７】
前記モデルは、特定ユニット、送風機、再生装置サイクロン、バルブ／流量および湿式ガス圧縮機周辺のタグ間の整合性を更に識別して、関係パターンに何らかの初期破綻があれば表示することを特徴とする請求項１３に記載の方法。
【請求項１８】
前記モデルは、特殊原因による動作における偽陽性を除去すべくモデル計算を抑制する工程を更に含むことを特徴とする請求項１７に記載の方法。
【請求項１９】
（ａ）前記モデルを決定する工程が、疑わしいデータに基づく粗いモデルから出発し、
（ｂ）前記粗いモデルを用いて高品質の訓練データを収集してモデルを改良し、
（ｃ）工程（ｂ）を繰り返して更にモデルを改良する
ことを特徴とする請求項９に記載の方法。
【請求項２０】
前記訓練データは、前記プロセス・ユニットのモデルの履歴データを含むことを特徴とする請求項１９に記載の方法。
【請求項２１】
前記モデルは、変換された変数を含むことを特徴とする請求項２０に記載の方法。
【請求項２２】
前記変換された変数は、蒸留塔における還流対供給速度、高純度蒸留における組成の対数、圧力補償された温度測定値、副流歩留まり、流量対バルブ位置および反応速度対ｅｘｐ（温度）を含むことを特徴とする請求項２１に記載の方法。
【請求項２３】
２個の変数の測定値のいくつかの対は、変数の１個により、動的伝達関数を用いて時間同期化されることを特徴とする請求項１９に記載のモデル。
【請求項２４】
プロセス動作における動作点の変化に影響を受けるプロセス測定値の変数は、変化する平均を減算することにより偏差変数に変換されることを特徴とする請求項２０に記載のモデル。
【請求項２５】
前記モデルは、ノイズについて修正されることを特徴とする請求項２０に記載の方法。
【請求項２６】
前記モデルは、変数のノイズが多い測定値のフィルタリングまたは削除により修正されることを特徴とする請求項２５に記載の方法。
【請求項２７】
変数の測定値は、スケーリングされることを特徴とする請求項２０に記載の方法。
【請求項２８】
前記測定値は、当該変数の期待正常範囲にスケーリングされることを特徴とする請求項２７に記載の方法。
【請求項２９】
異常現象モニターのリストがオペレータ向けに自動的に識別、分離、ランク付けおよび表示されることを特徴とする請求項４に記載の方法。
【請求項３０】
オペレータに対し、事象の調査を支援すべく診断情報が異なるレベルで詳細に提示されることを特徴とする請求項９に記載の方法。
【請求項３１】
主成分の個数は、主成分の係数同士の大きさが等しくなるように選択されていることを特徴とする請求項２０に記載の方法。
【請求項３２】
前記主成分は、オンライン測定値により提供されるプロセス変数を含むことを特徴とする請求項４に記載の方法。
【請求項３３】
測定値のいくつかの対は、動的フィルタを用いて、前記変数の１個と時間同期化されることを特徴とする請求項３２に記載の方法。
【請求項３４】
プロセス動作の動作点変化に影響を受けるプロセス測定値変数は、偏差変数に変換されることを特徴とする請求項３２に記載の方法。
【請求項３５】
主成分の個数は、連続する成分により表わされる全体プロセス変動の大きさにより選択されることを特徴とする請求項３２に記載の方法。
【請求項３６】
石油精製所のＦＣＣＵのプロセス・ユニットのいくつかにおける異常事象検知（ＡＥＤ）を行なうシステムであって、
（ａ）プロセス・ユニットの動作を記述するプロセス・ユニットのモデルの組；
（ｂ）前記プロセス・ユニットに異常状態が存在することを示すべく現在の動作が期待正常動作と異なるか否かを示すディスプレイ；および
（ｃ）ＨＤＣプロセス・ユニットにおける異常状態の根本原因を示すディスプレイ
を含むシステム。
【請求項３７】
各プロセス・ユニットの前記モデルは、主成分モデルまたは工学モデルのいずれかであることを特徴とする請求項３６に記載のシステム。
【請求項３８】
前記ＦＣＣＵは、３個の動作区画に分割されていて、各々の区画に主成分モデルを有することを特徴とする請求項３７に記載のシステム。
【請求項３９】
前記主成分は、オンライン測定により提供されるプロセス変数を含むことを特徴とする請求項３８に記載のシステム。
【請求項４０】
前記モデルは、オペレータにより誘発された通知および偽陽性を除去すべくモデル計算を抑制する工程を更に含むことを特徴とする請求項３８に記載のシステム。
【請求項４１】
（ａ）前記モデルを導く工程が、疑わしいデータに基づく初期モデルの取得から出発し、
（ｂ）前記初期モデルを用いてデータを整備してモデルを改良し、
（ｃ）工程（ｂ）を繰り返してモデルを改良する
ことを特徴とする請求項３７に記載のシステム。
【請求項４２】
前記訓練データ集合は、モデル構築用に前記プロセス・ユニットの履歴データを含むことを特徴とする請求項４１に記載のシステム。
【請求項４３】
前記モデルは、変換された変数を含むことを特徴とする請求項４２に記載のシステム。
【請求項４４】
前記変換された変数は、蒸留塔における還流対合計生成物流量、蒸留塔における組成およびオーバーヘッド圧力の対数、圧力補償温度測定値、流量対バルブ位置および反応床の温度差および圧力差を含むことを特徴とする請求項４３に記載のシステム。
【請求項４５】
測定値のいくつかの対は、動的フィルタを用いて、前記変数の１個に時間同期化されることを特徴とする請求項４２に記載のシステム。
【請求項４６】
プロセス動作における動作点の変化に影響を受けるプロセス測定値の変数は、偏差変数に変換されることを特徴とする請求項４２のシステム。
【請求項４７】
モデルの識別に先立って、変数の測定値がスケーリングされることを特徴とする請求項４２に記載のシステム。
【請求項４８】
前記測定値は、当該変数の期待正常範囲によりスケーリングされることを特徴とする請求項４７に記載のシステム。
【請求項４９】
主成分の個数は、連続する成分により表わされる全般的プロセス変動の大きさにより選択されることを特徴とする請求項４２に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、供給予熱装置、反応炉、再生装置、湿式ガス圧縮機、主精留装置および下流ライトエンド処理塔を含む流動接触分解ユニット（ＦＣＣＵ）の動作に関する。特に、本発明は、プロセスが正常動作から逸脱している場合の検知およびプロセスの異常部分を分離する通知の自動生成に関する。
【背景技術】
【０００２】
接触分解は、重質油をより有用なガソリンおよびより軽量の生成物に変換する最重要且つ広範に利用されている精練プロセスの一つである。当該プロセスは、現代的な精練所の中心であるＦＣＣＵにおいて実行される。ＦＣＣＵは、反応炉および再生装置を含む、複雑且つ緊密に統合されたシステムである。図２３に、典型的なＦＣＣＵのレイアウトを示す。新規供給およびリサイクルの流れは、熱交換器により予熱され、フィードライザ底部で装置に投入されて高温再生触媒に混ぜ合わせられる。ＦＣＣプロセスは触媒を、蒸気に通された際に流体として振舞う極めて微細な粒子（約７０ミクロン）の形状で使用する。平均ライザ反応炉温度は９００〜１０００°Ｆ、油供給温度は５００〜８００°Ｆ、触媒の再生装置出口温度は１２００〜１５００°Ｆの範囲にある。当該プロセスは、フィードライザ線内で高温の供給油に触媒を接触させる工程を含む。触媒からの熱で供給油を蒸発させて所望の反応温度まで上昇させる。分解反応は、供給油がライザ内の高温触媒に接触した時に始まり、反応炉内で油蒸気が触媒から分離されるまで続く。分解反応が進行するにつれて、触媒の表面におけるコークスの形成により、触媒が次第に不活性化する。使用済み触媒は再生装置内に流入し、コークス堆積物を空気で焼き払うことにより再生される。煙道ガスおよび触媒は、サイクロン集塵器において分離される。流動化された触媒は、反応ゾーンと再生ゾーンの間を連続的に循環して、再生装置から供給油と反応炉へ熱を伝導する媒体として作用する。触媒および炭化水素蒸気は機械的に分離され、触媒に残留している油は、触媒が再生装置に入る前に水蒸気蒸留により除去される。いくつかの装置の触媒は、再生装置を出る際に水蒸気蒸留されて、触媒が油供給と接触する前に吸着酸素を除去する。炭化水素蒸気は、合成精留装置へ送られて液体およびガス状生成物に分離される。これらは次いで、下流のライトエンド塔で更に精製されて、ガソリンその他の商品生成物が製造される。ＦＣＣＵおよび下流装置の完全な模式図を図２４に示す。
【０００３】
ＦＣＣＵの複雑な動的特性に起因して、各種の根本原因により異常なプロセス動作が生じ易く、深刻な問題に発展して工場の操業停止にまで至る恐れがある。これらの動作は、生産の機会損失から、設備の損傷、環境汚染、死傷事故の発生に到る深刻な安全および経済的な打撃を与えかねない。オペレータの主要な仕事は、異常な状況の原因を突き止めて、適時的且つ効率的に補償または是正措置を行なうことである。
【０００４】
【特許文献１】米国特許第５，８５９，９６４号明細書
【特許文献２】米国特許第５，９４９，６７８号明細書
【特許文献３】米国特許第６，５２２，９７８号明細書
【特許文献４】米国特許第６，３６８，９７５号明細書
【特許文献５】米国特許第６，４６６，８７７号明細書
【特許文献６】米国特許第６，５２１，０８０号明細書
【特許文献７】米国特許第６，５６４，１１９号明細書
【特許文献８】米国特許第６，６３６，８４２号明細書
【非特許文献１】Ｃａｒｄｏｓｏ，Ｊ．ら著「Ｆｕｚｚｙ  Ｐｅｔｒｉ  Ｎｅｔｓ：Ａｎ  Ｏｖｅｒｖｉｅｗ」（１３ｔｈ  Ｗｏｒｄ  Ｃｏｎｇｒｅｓｓ  ｏｆ  ＩＦＡＣ、Ｖｏｌ．Ｉ：Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  ＩＩ、Ｄｉｓｃｒｅｔｅ  Ｅｖｅｎｔ  Ｓｙｓｔｅｍｓ、Ｓａｎ  Ｆｒａｎｃｉｓｃｏ、ＣＡ、ＵＳＡ、１９９６年６月３０日〜７月５日、４４３〜４４８ページ）
【非特許文献２】Ｊａｃｋｓｏｎ，Ｅ．著「Ａ  Ｕｓｅｒ’ｓ  Ｇｕｉｄｅ  ｔｏ  Ｐｒｉｎｃｉｐａｌ  Ｃｏｍｐｏｎｅｎｔ  Ａｎａｌｙｓｉｓ」（Ｊｏｈｎ  Ｗｉｌｅｙ  ＆  Ｓｏｎｓ、１９９１）
【非特許文献３】Ｋｏｕｒｔｉ，Ｔ．著「Ｐｒｏｃｅｓｓ  Ａｎａｌｙｓｉｓ  ａｎｄ  Ａｂｎｏｒｍａｌ  Ｓｉｔｕａｔｉｏｎ  Ｄｅｔｅｃｔｉｏｎ：Ｆｒｏｍ  Ｔｈｅｏｒｙ  ｔｏ  Ｐｒａｃｔｉｃｅ」（ＩＥＥＥ  Ｃｏｎｔｒｏｌ  Ｓｙｓｔｅｍｓ  Ｍａｇａｚｉｎｅ、２００２年１０月、１０〜２５ページ）
【非特許文献４】Ｋｕ，Ｗ．著「Ｄｉｓｔｕｒｂａｎｃｅ  Ｄｅｔｅｃｔｉｏｎ  ａｎｄ  Ｉｓｏｌａｔｉｏｎ  ｆｏｒ  Ｓｔａｔｉｓｔｉｃａｌ  Ｐｒｏｃｅｓｓ  Ｃｏｎｔｒｏｌ  ｉｎ  Ｃｈｅｍｉｃａｌ  Ｐｒｏｃｅｓｓｅｓ」（ＰｈＤ論文、Ｌｅｈｉｇｈ  Ｕｎｉｖｅｒｓｉｔｙ、１９９４年８月１７日）
【非特許文献５】Ｍａｒｔｅｎｓ，Ｈ．およびＮａｅｓ，Ｔ．著「Ｍｕｌｔｉｖａｒｉａｔｅ  Ｃａｌｉｂｒａｔｉｏｎ」（Ｊｏｈｎ  Ｗｉｌｅｙ＆Ｓｏｎｓ、１９８９）
【非特許文献６】Ｐｉｏｖｏｓｏ，Ｍ．Ｊ．ら著「Ｐｒｏｃｅｓｓ  Ｄａｔａ  Ｃｈｅｍｏｍｅｔｒｉｃｓ」（ＩＥＥＥ  Ｔｒａｎｓ  ｏｎ  Ｉｎｓｔｒｕｍｅｎｔａｔｉｏｎ  ａｎｄ  Ｍｅａｓｕｒｅｍｅｎｔ、Ｖｏｌ．４１、Ｎｏ．２、１９９２年４月、２６２〜２６８ページ）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
現在実用化されているのは、高度なプロセス制御装置アプリケーションを用いて軽微なプロセス撹乱に応答して自動的にプロセスを調整し、中程度〜重度の異常動作の場合はプロセスへの人の介入に依存して、極めて重大な異常動作の場合は自動非常プロセス停止システムを用いることである。コンソールのオペレータに異常なプロセス動作が始まったことを知らせる通常の方式は、プロセス警報を介するものである。これらの警報は、主要なプロセス測定値（温度、圧力、流量、レベルおよび組成）が動作範囲の所定の静的な組から逸脱した場合に作動する。この通知技術は、主要な測定値同士がＦＣＣＵ等の複雑なプロセスにおける相関が求められた場合、偽陽性の発生率を低く抑えながらタイムリな警報を発することは困難である。
【０００６】
典型的なＦＣＣＵの動作に対応する６００種を超える主要なプロセス測定値が存在する。従来型の分散制御システム（ＤＣＳ）配下で、オペレータはセンサーのリストおよびその傾向を調べ、正常ＦＣＣＵ動作について記憶している知識と比較して、自身のスキルを用いて潜在的な問題を検知しなければならない。動作中のＦＣＣＵには極めて多くセンサーがあるため、異常を見過ごしやすい場合がある。現在のＤＣＳに基づく監視技術では、オペレータが有する唯一の自動化された検知支援機能は、所定の限度を超えた場合に各センサーに警報を発することに基づくＤＣＳ警報システムである。ＦＣＣＵのように大規模且つ複雑なあらゆるプロセスにおいて、この種の通知は往々にして到着が遅過ぎるため、オペレータが問題を抑制すべく予防的処置を講じるための充分な時間が取れない。本発明は、ＦＣＣＵのオペレータにより効果的に通知を提供する。
【課題を解決するための手段】
【０００７】
本発明は、ＦＣＣＵのプロセス・ユニット用に異常事象を検知する方法である。異常事象検知（ＡＥＤ）システムは、多数の高度に統合された動的プロセス・ユニットを含んでいる。本方法は、現在の動作を、対応ユニットの正常動作の各種のモデルと比較する。ユニットの動作と正常動作の差違がプロセス・ユニットにおける異常状態を示す場合、異常状態の原因が特定されて、オペレータが是正処置を講じるべく関連する情報を効率的に提示する。
【発明を実施するための最良の形態】
【０００８】
本発明は、異常事象検知（ＡＥＤ）技術を用いて、オペレータに対し、ＦＣＣＵ区画における異常状態の早期通知を提供する方法である。
【０００９】
スナップショットに基づいてオン／オフの表示だけを行なう警報技術とは対照的に、本方法はファジー論理を用いて、動作上の問題に寄与する異常性を立証する複数の証拠を組み合わせて、その確率をリアルタイムに推定する。この確率を、連続信号としてオペレータに提示することにより、現行の単一センサーの警報に基づくオン／オフ方法に付随するチャタリングが除去される。オペレータには、問題の完全な調査および根本原因まで掘り下げて的を射た処置を講じることを可能にする一組のツールが提供される。この方式は、オペレータに対し、従来の警報システムより数分〜数時間も早い、異常動作の高度な警告を与えることが示されてきた。このような早期の通知により、オペレータは情報に裏付けられた決断を行ない、問題の拡大や災害を回避するための是正処置を講じることができる。本方法はＦＣＣＵに成功裏に適用されている。例えば、図２７に、触媒循環問題（下位問題の詳細は後述する）に対する完全な掘り下げを示す。
【００１０】
ＦＣＣＵアプリケーションは、特定の動作に関する知識の各種の出典を用いて、主成分分析法（ＰＣＡ）、部分最小二乗法（ＰＬＳ）に基づく推論モデル、相関に基づく工学モデル、および関連するセンサー変換からの表示を組み合わせて、いくつかのファジー論理ネットワークに統合する。このファジー論理ネットワークは、証拠を集約して、潜在的な問題の信頼度を示す。従って、ネットワークは問題の初期拡大段階において、より高い信頼度を以ってこれを検知することができ、オペレータに対し、重大な事故を回避すべく補償的または是正的処置を講じるために必須なリードタイムを与えることができる。これは、ＤＣＳシステムからの単一センサー警報に基づいてＦＣＣＵを監視する現行の実用方式に対する主要な優位性である。ＦＣＣＵの複雑且つ高速の動的特性に起因して、オペレータが動作上の問題を回避するには警報が届くのが遅すぎるか、或いは、膨大な警報がオペレータに殺到して当人を混乱させ、対応を支援するどころか却って阻害してしまう場合が極めて多い。
【００１１】
接触分解ユニットは、設備グループ（主要な機能区画または動作区画と称する）に分けられる。これらの設備グループは、その設計に応じて、異なる接触分解ユニット毎に異なっていてよい。接触分解ユニットの特定のプロセス・ユニットを含む設備グループを選択する手順を付録１に記述している。
【００１２】
好適な実施形態において、本発明は流動接触分解ユニット（ＦＣＣＵ）動作を以下の全般的モニター即ち、
１．全般的ＦＣＣＵ例外動作
２．全般的ＦＣＣＵ極限動作
３．全般的触媒ライトエンド例外動作
４．全般的触媒ライトエンド極限動作
およびこれらの特別要注意モニターに分類する。
１．反応炉−再生装置触媒循環
２．再生装置スタックバルブ動作
３．サイクロン動作
４．送風機動作
５．炭素収支チェック
６．主精留装置への触媒キャリーオーバー
７．湿式ガス圧縮機
８．バルブ−流量整合性モデル
【００１３】
全般的モニターは、全般的動作に何らかの逸脱があれば検知して多数のセンサーをカバーすべく「グロスモデル検査」を実行する。特別要注意モニターは、潜在的に重大な要注意領域をカバーして、早期検知用の的を絞ったモデルで構成されている。本アプリケーションは、これらのモニターの全てに加え、正常／通常の動作関連事象から生じる通知を抑制して特別な原因の動作による偽陽性を除去する等を行なう、いくつかの実用的なツールを提供する。
【００１４】
Ａ．オペレータ・インターフェース
オペレータ・ユーザー・インターフェースは、オペレータにプロセスの全体像を示すため、システムに必須な構成要素である。ディスプレイは、オペレータにＦＣＣＵ動作の全体像を素早く提示して、拡大しつつある異常があればその確率を示すことを意図している。
【００１５】
図２５に、本システムのオペレータ・インターフェースを示す。オペレータ・インターフェースの設計検討に関する詳細な説明を、付録１のＩＶ節内の「ＡＥＤ用のＰＣＡモデルおよび簡単な工学モデルの配備」節のサブセクション「オペレータ対話環境およびインターフェースの設計」に記述している。本インターフェースは上述の異常事象モニターで構成されている。これは、各動作域において重要な異常兆候のリストを表わすために開発された。モデルの結果を主要なセンサーの状態と比較することにより異常兆候が生成される。ファジー論理を用いて異常な徴候を集約し、ある問題生じる単一の確率を評価する。各区画の正常動作に関する固有の知識に基づいて、ある問題の確率を評価すべくセンサーおよびモデル残差からの入力を取得するためにファジー論理ネットワークを開発した。図２６は、図２８に示す対応ファジー論理ネットワークを用いる触媒循環問題の確率を示す。図２７は、触媒循環問題の完全な掘り下げを示す。図２８の緑色のノードは、下位問題を組み合わせて「触媒循環問題」の最終的な確実性を決定する様子を示す。異常状態の推定された確率を連続的な傾向として操作チームに提示して状態の進行状況を示す。図２９に、全ての領域について連続的な信号徴候を提示している、当該問題のオペレータ向け表示を示す。これはオペレータにとって、各々のセンサーの状態を個別に検査する必要なしにプロセスの健全性の全体像を得るために重要な利点をもたらす。更に重要なのは、このカバー範囲が広く、事象を見逃す可能性が低くなるため、オペレータに「安心感」を与える点である。従って、これは正常性インジケータとしても利用可能である。確率が０．６に達した場合、問題インジケータは黄色（注意）に変わり、確率が０．９に達した場合、インジケータは赤色（警告）になる。
【００１６】
本発明は、触媒循環（ＣＣＲ）、反応炉−再生装置動作（ＦＣＣ）、および触媒ライトエンド（ＣＬＥ）動作の領域をカバーすべく３種類の主成分分析（ＰＣＡ）を含んでいる。ＰＣＡモデルのカバー範囲は、異なるプロセス・ユニットの相互作用に基づいて決定されており、各モデルはこれを考慮するために重なり合うセンサーを備えている。反応炉−再生装置領域には顕著な相互作用があるため、１種類のＰＣＡモデルがそれら両方の動作をカバーすべく設計されている。触媒循環ＰＣＡは、反応炉−再生装置間の触媒流を特に監視することを狙いとして的を絞ったモデルである。ＦＣＣＵからの生成物を処理する触媒ライトエンド（ＣＬＥ）塔は、別々のＰＣＡに含まれている。また、重大な事象に発展する恐れがある状況を監視することを意図した多くの特別要注意モニターが存在する。その目的は、オペレータが処置を行なうために充分なリードタイムを持てるよう、早めに問題を検知することである。
【００１７】
正常動作時において、オペレータは供給速度の変更やセットポイントの移動等、ＰＣＡモデル内のいくつかのセンサーに持続期間が短い高残差を生じ得るいくつかの定常的な処置を実行する。そのような通知は冗長であって新たな情報を与えないため、本発明ではそれらの開始を検知して通知を抑制するための仕組みを組み込んでいる。
【００１８】
オペレータは、色が緑色から黄色に、次いで赤色に変化する警報三角形により切迫しつつある問題について知らされる。本アプリケーションは、優先順位付けされた下位問題のリストを閲覧することにより更に問題を調べる掘り下げ能力をオペレータに提供する。この新規な方法は、下位問題まで深堀する能力をオペレータに。これによりオペレータは、根本原因の探索を絞り込むことができる。図２９に、触媒循環、ＦＣＣ例外およびＦＣＣ極限動作が注意喚起する警告を有する様子を示す。これは、オペレータが補償または是正処置が講じることができるよう、状態の根本原因の分離と診断を支援する。図３０に、黄色の三角形をクリックした結果、即ち逸脱しているセンサーを逸脱の度合いでソートした残差を示すパレート図を示す。
【００１９】
本アプリケーションは、パレート図方式を用いてオペレータに対し極めて広範囲に情報を示す。一連の提示は、正常動作からの個々の逸脱の降順になされている。これにより少数の重大な不良要因まで絞り込まれたプロセスを簡潔且つ手短に見ることが可能になるため、コンソール・オペレータは講じるべき処置について情報に裏付けられた決定を行なうことができる。図３０に、パレート図として構成されたセンサーのリストを介してこの機能を示す。個々のバーをクリックすれば、センサーのタグ傾向対モデル予測を示すカスタマイズされたプロット図が作成される。オペレータはまた、「マルチ傾向ビュー」を用いて問題含みのセンサー群の傾向を見ることができる。例えば、図３１に、図３０のパレート図のセンサーの値の傾向とモデル予測を示す。図３２は同じ概念を、今回は正規化された射影偏差誤差に基づいて、バルブ流量・モデルのランキングに適用したものを示す。この場合、バーをクリックすれば、正常動作の境界（図３３）に関する現在の動作点を示すＸ−Ｙ散布図を生成する。当該アプリケーションの別の例を、図３４にパレート図、および送風機モニター用のＸ−Ｙプロット図に示す。
【００２０】
ＰＣＡモデルに加え、工学関係および部分最小二乗法に基づく推論を用いて構築された多くの特別要注意モニターがある。これらは、送風機および湿式ガス圧縮機等、反応炉−再生装置領域内の重要な設備に適用される。これらのモニターの基礎として、単一の異常信号を生成するファジー論理ネットワークがある。
【００２１】
要約すれば、本発明の利点には以下のものが含まれる。
１．監視のため、全般的ＦＣＣＵ動作を１０個動作領域、即ち反応炉−再生装置、触媒ライトエンド塔、触媒循環、スタックバルブ、サイクロン、送風機、炭素収支、主精留装置への触媒キャリーオーバー、湿式ガス圧縮機、バルブ−流量モデルへ分解すること。
２．全般的ＦＣＣＵの動作状態を１２個の単一警報に要約すること。
３．ＰＣＡモデルが、モデルによりカバーする６００個以上のセンサーのモデル予測を提供すること。
４．これら６００個以上のセンサーの異常な逸脱が、３個のＰＣＡモデルの二乗誤差和に基づいて５種の警報に要約されること。
５．再生装置スタックバルブ、再生装置サイクロン、送風機、炭素収支、触媒キャリーオーバーおよび湿式ガス圧縮機に対する工学モデル推論により、的確且つ早期な検知能力が強化されること。
６．バルブ−流量モデルが、制御処置に影響するため不調の源であるかまたはその影響を受ける制御ループを監視する強力な方法を提供すること。
７．特別な原因／定常動作から生じる事象が抑制されて、偽陽性が除去される。６００個以上の個別タグから単に１２個の信号へ大幅に規模を縮小することにより、偽陽性率が大幅に減少する。ＰＣＡモデリング方式は本質的に、単一センサー警報問題を洗練された方法で解決する。
【００２２】
Ｂ．ＦＣＣＵ用のＡＥＤモデルの構築および配備
本アプリケーションは、ＦＣＣＵの異常動作を検知すべく、ＰＣＡモデル、工学モデルおよび発見的方法を有する。第１工程は、動作上の履歴的な問題の要注意ユニットを分析することを含む。この問題識別工程は、アプリケーションの範囲を定めるために重要である。
【００２３】
これらのモデルの構築を、付録１に一般的に記述する。流動接触分解ユニットのこれらのモデルを構築に関するいくつかの特定の要注意事項を以下に述べる。
【００２４】
問題の識別
アプリケーション構築の第１工程は、プロセス動作の利点となる重要な問題を識別することである。異常事象検知アプリケーションは一般に、２種類の異なる問題に適用することができる。第一は、何らかの異常事象を求めてプロセス領域全体を監視する汎用的な異常事象アプリケーションである。このタイプは、数百個の測定値を使用するが、特定の異常動作の履歴的な記録を一切必要としない。本アプリケーションは、異常事象を検知してプロセスの部分（タグ）に紐付けるだけである。問題の診断には、オペレータまたは技術者のスキルを必要とする。
【００２５】
第２のタイプは、特定の異常動作に着目する。このタイプは、異常が検知された場合に特定の診断を提供する。通常、少数（５〜２０個）の測定値しか含まないが、事象の履歴データレコードを必要とする。このモデルは、ＰＣＡ／ＰＬＳまたは簡単な工学相関（質量／エネルギー収支ベース）に基づいている。本明細書は、広い範囲をカバーすべく両方の種類のアプリケーションを記述している。オペレータまたは技術者は次いで、自身のプロセス知識／専門性に頼って正確に原因を診断する。通常、大部分の事象は、主に計測器やバルブに関する問題の結果であると思われる。
【００２６】
異常事象検知の問題を選択する際に、以下の問題含みの特性を考慮しなければならない。異常が稀（３〜４ヶ月毎）であるため、異常事象検知器を作成する努力が無駄に思えるかもしれない。また、特定の異常が３〜４ヶ月毎にしか発生しない場合、個々のオペレータは数年間にわたり当該事象を見ない可能性もある。その結果、一旦事象が起きた際に、オペレータはどう対処すべきかわからない恐れがある。従って問題の識別を充分広範に行なって、オペレータが定期的にアプリケーションと対話するようにしなければならない。
【００２７】
問題を特定しようとする場合に、異常事象検知のアプリケーションを正当化するのに充分な数の異常事象が存在しないとの誤った印象を現場担当者から受けるのが一般的である。一般に、異常事象がプロセスに対して影響を及ぼす頻度が過小評価されがちである。その理由は以下の通りである。
往々にして、異常事象の記録および解析が行なわれない。重大な損失を与えたものだけが追跡および解析される。
オペレータは毎日これらを取扱うため、異常事象も正常動作の一部と見なしがちである。
定期的に繰り返し発生する異常事象が存在しない限り、本アプリケーションは異常事象を定期的に（例えば週５回超）「見る」ためにプロセスの充分広範な部分を包含していなければならない。
【００２８】
１．ＰＣＡモデル
ＰＣＡモデルは、ＦＣＣＵ  ＡＥＤの心臓部である。ＰＣＡは、実際のプロセス変数を、元の変数の線形結合である主成分（ＰＣ）と呼ばれる一組の「直交」または独立変数に変換する。基礎となるプロセスには、プロセスに影響を及ぼす特定の独立した事象を表わす多くの自由度がある。これらの異なる独立した事象は、プロセス変動とし当該プロセスのデータに出現する。プロセス変動は、供給速度の変更等の意図的な変更、または周囲温度の変化等の意図しない撹乱に起因する場合がある。
【００２９】
各々の主成分は、当該プロセスに対するこれらの異なる独立事象により生じたプロセス変動の一意な部分を捕捉する。主成分は、変動が減少する順序で抽出される、後続する各々の主成分は、全体のプロセス変動性のより小さい部分を捕捉する。主要な主成分は、プロセス変動の重要な基礎となる発生源を表わしていなければならない。例えば、供給速度の変更が最も大きいプロセス変動の単一発生源であるため、第１の主成分はしばしば供給速度変更の効果を表わす。
【００３０】
本アプリケーションは、プロセスの主成分分析法（ＰＣＡ）に基づいており、「正常動作」の経験的モデルを精製する。ＰＣＡモデルの構築プロセスは、付録１の「ＡＥＤ用ＰＣＡモデル構築」節に詳述している。以下に、ＦＣＣＵ用の異常事象検知アプリケーション構築へ向けてＰＣＡを適用する場合に必要な特別の考慮について述べる。
【００３１】
ＦＣＣＵ  ＰＣＡモデルの構築
本アプリケーションは、反応炉−再生装置領域（ＦＣＣ−ＰＣＡ）、触媒循環（ＣＣＲ−ＰＣＡ）、および触媒ライトエンド塔（ＣＬＥ−ＰＣＡ）をカバーするＰＣＡモデルを有している。これにより、ＦＣＣ動作全般および初期警報を広範囲にカバーすることが可能になる。
【００３２】
ＰＣＡモデルの構築は以下の工程を含む。
１）入力データおよび動作範囲の選択
２）履歴データの収集および前処理
３）データおよびプロセスの分析
４）初期モデルの生成
５）モデルのテストおよびチューニング
６）モデルの配備
【００３３】
ＰＣＡモデルの構築に関する一般的な原理を、付録１の「ＡＥＤ用ＰＣＡモデルの構築」節内のサブセクションＩ「概念的ＰＣＡモデルの設計」に記述している。これらの工程はモデル構築における努力の主要部をなす。ＰＣＡモデルはデータ駆動であるため、正常動作を表わす訓練データの高い品質および量が極めて重要である。基本的な構築戦略は、極めて粗いモデルから出発し、次いで当該モデルの忠実度を逐次改良してくことである。これには、モデルがどの程度実際のプロセス動作に匹敵するかを観察し、および観察に基づいてモデルを再訓練することが必要である。そのような工程について、次に簡潔に述べる。
【００３４】
入力データおよび動作範囲の選択
ＰＣＡモデル内のタグのリストがカバー範囲を表わすため、要注意領域内の全てのタグの総合的なリストから出発する。測定値および変数の選択プロセスを、付録１の「ＡＥＤ用ＰＣＡモデルの構築」節内のサブセクションＩＩ「入力データおよび動作範囲の選択」で概説している。信頼できない、または誤った振る舞いを示すことがわかっている測定値は全てリストから除外すべきである。初期ＰＣＡモデルが得られたならば、反復的な手順を用いて更に測定値の削減が実行される。
【００３５】
履歴データの収集および前処理
正常動作の良いモデルを構築するには正常動作の訓練データ集合が必要である。このデータ集合は以下のようでなければならない。
・正常動作範囲にわたる
・正常動作データだけを含む
【００３６】
ある拠点で異常事象の履歴の完全な記録を保存していることは極めて稀であるため、履歴データは訓練データ集合を生成する出発点としてしか利用できない。オペレータログ、オペレータ変更ジャーナル、警報ジャーナル、計器保守記録等の操作記録は、異常プロセス履歴の部分的な記録を提供する。データ収集のプロセスを、付録Ｉ．の「ＡＥＤ用ＰＣＡ  モデルの構築」節内のサブセクションＩＩＩ「履歴データの収集」に詳述している。
【００３７】
ＦＣＣＵの場合、履歴データは、夏期と冬期の両期間をカバーすべく１．５年の動作にわたっている。１分毎にデータ平均を取った結果、各々のタグについて時点の数が約７００，０００個以上になることがわかった。基礎となる情報を保持したまま、データ集合をより扱い易くするために工学的判定が適用されて、６番目毎の点を保持した結果、各センサーについて約１００，０００個以上の点が得られた。これにより、代表的な挙動をＰＣＡモデルが補足することができる。
【００３８】
全てのタグが内部に含まれる変動範囲／情報を判定すべく、平均、最小／最大および標準偏差等の基本的統計量を計算した。また、動作ログを調べて、既知のユニット停止または異常動作を有するウインドウ内に含まれるデータを除外した。各々の候補測定値を精査して、訓練データ集合に含める際の適正性を判定した。
【００３９】
バランスのよい訓練データ集合の生成
動作ログを用いて、履歴データを既知の異常動作が生じた期間と、識別された異常動作が生じていない期間に分ける。識別された異常動作が存在しないデータが予備訓練データ集合となる。
【００４０】
これらの除外がなされたならば、第１の粗いＰＣＡモデルを構築することができる。これは極めて粗いモデルになるため、保持される主成分の正確な数は重要ではない。これは、モデルに含まれる測定値の個数の５％以下であるべきである。ＰＣの個数は、最終的にプロセスの自由度の数に一致しなければならないが、これはプロセス撹乱の異なる全ての発生源を含んでいるため、通常は未知である。何個の主成分を含めるべきかを決定する標準的な方法がいくつか存在する。またこの段階で、可変スケーリングへの統計的なアプローチを用いるべきである。即ち、全ての変数をユニット分散に応じてスケーリングする。
【００４１】
ここで、訓練データ集合をこの予備モデルに通して、データがモデルに一致しない期間を識別しなければならない。これらの期間を調べて、その時点で異常事象が生じていたか否かを知る必要がある。そうであると判定された場合、これらの期間には、既知の異常事象が生じている時間としてフラグを立てなければならない。これらの期間を訓練データ集合から除外して、修正されたデータでモデルを再構築しなければならない。データおよびプロセス分析を用いてバランスの取れた訓練データ集合を生成するプロセスを、付録１の「ＡＥＤ用ＰＣＡモデルの構築」節の第ＩＶ節「データおよびプロセス分析」に概説している。
【００４２】
初期モデルの構築
モデル構築の戦略は、極めて粗いモデル（疑わしい訓練データ集合の結果）から出発し、次いで当該モデルを用いて高品質訓練データ集合の収集することである。次いでこのデータを用いてモデルを改良し、次いで当該モデルを用いて品質が向上した訓練データを収集し続ける。このプロセスは、モデルが満足すべき状態になるまで繰り返される。
【００４３】
特定の測定値が選択されて、訓練データ集合が構築されたならば、標準的な統計ツールを用いて短時間でモデルを構築することができる。各々の主成分により捕捉された百分率分散を表わすそのようなプログラムの例を図４５に示す。モデル構築プロセスについて、付録１の「ＡＥＤ用ＰＣＡモデルの構築」節内の第Ｖ節「モデルの生成」に記述している。
【００４４】
モデルのテストおよびチューニング
初期モデルが構築されたならば、新規の訓練データ集合を生成することにより強化する必要がある。これは、当該モデルを用いてプロセスを監視することにより行なわれる。モデルが潜在的な異常状況を示したならば、技術者はプロセス状況を調べて分類しなければならない。技術者は３種類の異なる状況、即ち、何らかの特別なプロセス動作が生じている、実際に異常な状況が生じている、或いは、プロセスは正常で誤った徴候を示している、のいずれかを発見するであろう。
【００４５】
プロセスデータは、ガウスまたは正規分布をなさないであろう。従って、残余誤差の変動性から異常事象を検知すべくトリガを設定する標準的な統計的方法を用いるべきではない。その代わりに、トリガ位置は、モデルの使用経験に基づいて経験的に設定する必要がある。付録１の「ＡＥＤ用ＰＣＡモデルの構築」節内の第ＶＩ節「モデルのテストおよびチューニング」にモデルのテストおよび強化手順を記述している。
【００４６】
ＰＣＡモデルの配備
プロセス・ユニットにＡＥＤをうまく配備するには、精密なモデル、うまく設計されたユーザー・インターフェース、および適切なトリガ位置の組合せを必要とする。ＰＣＡモデルを配備する詳細な手順を、付録１の「ＡＥＤ用ＰＣＡモデルおよび簡単な工学モデルの配備」に記述している。
【００４７】
時間の経過につれて、構築担当者または拠点技術者は、モデルの１個を改良する必要がると判断する場合がある。プロセスの状態が変化したか、モデルが誤った徴候を提示しているかのいずれかである。この場合、訓練データ集合は、追加的なプロセスデータにより拡張されて、改良されたモデル係数を得ることができる。トリガ位置は、先に述べたのと同じ経験則を用いて再計算することができる。
【００４８】
もはや適切に現在のプロセス動作を表わさない古いデータは、訓練データ集合から削除しなければならない。特定の種類の動作がもはや実行されていない場合、当該動作からの全てのデータを削除しなければならない。大規模なプロセス変更の後では、訓練データおよびＡＥＤモデルを最初から作り直す必要があるかもしれない。
【００４９】
ＦＣＣＵ  ＰＣＡモデルは３８８個のタグの初期セットから出発し、次いで２２８個のタグへ修正された。触媒循環ＰＣＡ（ＣＣＲ−ＰＣＡ）モデルは２４個のタグを含んでいて、必須の触媒循環機能を監視する。触媒ライトエンドＰＣＡ（ＣＬＥ−ＰＣＡ）のタグは３６６個から２５６個まで絞り込まれ、回収に関わる下流区画、即ち主精留装置、脱エタン吸収装置、脱ブタン装置、スポンジ吸収装置、ＬＰＧスクラバーおよびナフサ・スプリッタ（図２４）をカバーする。ＦＣＣ−ＰＣＡモデルの詳細を付録２Ａに示し、触媒循環ＰＣＡモデルを付録２Ｂに記述し、ＣＬＥ−ＰＣＡモデルを付録２Ｃに記述している。
【００５０】
ＩＩ．ＡＥＤ工学モデル
工学モデルの構築
工学モデルは、異常状態の特定の検知に的を絞った推論および相関に基づくモデルで構成されている。工学モデルの構築に関する詳細な説明は、付録１の「ＡＥＤ用の簡単な工学モデル」節に見られる。
【００５１】
ＦＣＣＵアプリケーションの工学モデル要件は、履歴プロセスデータの工学的評価と、コンソール・オペレータおよび設備専門化へのインタビューを実施することにより決定された。工学的評価は、ＦＣＣＵ動作の重大な要注意領域および最悪ケースのシナリオを含んでいた。工学的評価の結論について述べるために、ＦＣＣＵ  ＡＥＤアプリケーション用に以下の工学モデルが構築された。
・触媒循環モニター
・スタックバルブ・モニター
・再生装置サイクロン動作モニター
・送風機動作モニター
・炭素収支モニター
・主精留装置モニターへの触媒キャリーオーバー
・湿式ガス圧縮機モニター
・バルブ−流量整合性モニター
【００５２】
推論を構築する手順は、先に議論したＰＣＡモデルの場合と全く同様である。しかし、特定の出力が予測されていない（全てのデータは入力である）ＰＣＡモデルの場合とは異なり、推論については予測のため望ましい変数がある。特定の入力に基づいて出力タグをモデル化するために部分最小二乗法（ＰＬＳ）を使用する。ＰＣＡの場合のように、これは測定値の選択とデータの前処理を必要とする。しかし、この場合は、測定値の選択は、出力タグにとって最適の予測手段である変数の観点から行なわれる。これは、プロセスの専門家と対話して数回の繰り返しを行ない、入力リストを最適な組に絞り込むことを含む。
【００５３】
触媒循環モニターは６個の下位問題領域、即ち（ａ）触媒循環の動作範囲、（ｂ）触媒循環ＰＣＡモデル残差、（ｃ）Ｒｘ−Ｒｇ  Ｊ屈曲密度、（ｄ）Ｒｘ−Ｒｇ触媒レベル、（ｅ）異常ＲｘＲｇデルタＰ制御、（ｆ）エネルギーと圧力収支の触媒循環計算との間の整合性、を用いて触媒循環の健全性を監視する。触媒循環は効率的なＦＣＣ動作の主要な要素であり、問題の早期検知は大幅な節約につながる。問題の完全な切り分けを図２７および図２８の対応するファジー・ネットワークに示す。
【００５４】
再生装置スタックバルブは、反応炉−再生装置の圧力差を維持するために必須である。これは、反応炉温度がスタックバルブ開口部を操作することにより反応炉−再生装置圧力差を調整する反応炉カスケード温度制御チェーンへの重要なリンクである。バルブを監視するために、（ａ）スタックバルブの正常動作範囲が調べられて、（ｂ）スタックバルブ開口部と差動圧力コントローラ出力との間の整合性が調べられる。図３５に、再生装置スタックバルブ・モニターの掘り下げを示す。付録３Ａにこのモニターの詳細を記述している。
【００５５】
再生装置サイクロンを用いて煙道ガスから触媒微粉を凝結させて触媒の損失を防止する。触媒はＦＣＣＵでリサイクルすべく触媒ホッパーで集められる。このモニターは、主要なモデルパラメータ、即ち煙道ガス温度、再生装置上部圧力、煙道ガスＯ２モデル、微粉ホッパー重量変化率、およびサイクロン差圧を検査する。付録３Ｂに本モニターの詳細を述べ、図３６にディスプレイを示す。
【００５６】
送風機は再生装置に空気を供給し、この空気を用いて反応炉から出た使用済み触媒に堆積したコークスを焼き払う。送風機はこのように、安定したＦＣＣ動作を維持する設備の必須部分である。送風機モニターは、タービン速度、デルタ空気温度、蒸気圧供給、空気流、蒸気流対タービン、排気温度を検査する。この場合の推論モデルは、（ａ）送風機への空気流（ｂ）タービンへの蒸気流、（ｃ）再生装置温度、および（ｄ）送風機吐出であった。推理における予測タグの詳細を付録３Ｃに示す。図３７にモニターの掘り下げを示す。
【００５７】
炭素収支モニターは、反応炉−再生装置の炭素収支を検査する。使用される証拠は、触媒循環ＰＣＡモデルのＴ統計量、煙道ガスＣＯ濃度、煙道ガスＯ２レベル、および他のいくつかのセンサーである。このモニターを図３８に示し、付録３Ｄにその詳細を述べている。
【００５８】
主精留装置への触媒キャリーオーバーは、反応炉ストリッパ・レベル、反応炉差圧、主精留装置へのスラリーのポンプ注入、およびストレーナー差圧を監視する。図３９にモニターを示す。付録は３Ｅにモニターの詳細を述べている。
【００５９】
湿式ガス圧縮機は、主精留装置のオーバーヘッド生成物を取り込み、下流のライトエンド塔で更に処理すべく圧縮する。ＷＧＣもまた塔圧を維持するため、モニターすべき別の必須の要注意領域である。当該モニターは、第２段の吸入流、タービンへの蒸気、第１段吐出流、触媒ガス出口温度を検査する。当該モニターの推論モデルは、（ａ）第２段圧縮機吸入流、（ｂ）タービンへの蒸気流（ｃ）第１段圧縮機吐出流、および（ｄ）触媒ガス吐出である。これらの推論の詳細を付録３．Ｆに記述し、図４０にモニターを示す。
【００６０】
流量−バルブ位置整合性モニターは、測定された流量（バルブ全体の圧力降下に関して補償済み）を流量のモデル推定と比較して導かれた。これらは、プロセスにおいて制御ループの状態が直接監視されているため強力な検査方法である。流量のモデル推定は、係数をバルブ曲線の式（線形または放物線と仮定される）にフィッティングさせることにより履歴データから得られた。初期アプリケーションにおいて、１２個の流量／バルブ位置整合性モデルが構築された。再生装置使用済みエアレーション蒸気バルブの例を図３３に示す。履歴的に信頼できない性能を示した制御ループについても、いくつかのモデルが構築された。バルブ流量モデルの詳細を付録３Ｇに記述している。
【００６１】
バルブ−流量モデルの不一致に加え、値超過用いて、制御バルブが制御可能な範囲を超えている旨をオペレータに通知する追加的な検査方法がある。図４１にファジー・ネットの両方の構成要素を示し、値超過の例を図４２に示す。
【００６２】
長期センサードリフトを補償すべく、時間依存ドリフト項がモデル推定に加えられた。オペレータはまた、センサーを再較正した後、または手動でバイパスバルブが変更された場合にドリフト項のリセットを要求することができる。流量推定装置に対するこの変更により、オンライン検知アルゴリズム内の実装の堅牢性が大幅に向上した。
【００６３】
工学モデル配備
ＡＥＤ内で工学モデルを実装するための手順は相当直接的である。ユニット（例：再生装置サイクロン、スタックバルブ、送風機、湿式ガス圧縮機動作等）内における特定の既知の種類の挙動を識別するモデルの場合、通知のトリガ位置は、コンソール・オペレータの入力と組み合わせた履歴データの分析から決定された。計算モデル（例：流量／バルブ位置モデル）の場合、通知のトリガ位置は、最初にモデル残差の標準偏差から導かれた。動作の最初の数ヶ月間、既知のＡＥＤ徴候はオペレータにより確認され、トリガ位置が適切であることが保障され、必要に応じて変更された。付録１の「ＡＥＤ用のＰＣＡモデルおよび簡単な工学モデルの配備」節に工学モデル配備の詳細を記述している。
【００６４】
ある状況の下では、バルブ／流量診断により、オペレータに冗長な通知が行く可能性がある。バルブ／流量の対に生じた問題に対する単一の警報をオペレータに与えるべく、バルブ／流量診断にモデル抑制を適用した。
【００６５】
Ｃ．追加的なＡＥＤツール
円滑な日常のＡＥＤ動作を実現すべく、ＡＥＤモデルを保守して実際の要注意事項への対応を支援するための各種のツールが提供される。
【００６６】
事象抑制／タグ無効化
定常的な毎日の運用において、オペレータは通常、多くの操作（例：セットポイントの変更、タグの保守、デコーキング、乾燥機の交換、再生）その他のプロセス変更を行なう。そのような既知の事象を予め抑制するために、システムは事象抑制を可能にする。セットポイント動作を実装する場合は常に、対応するＰＶその他の関連タグにおける工程の変更が通知を生成する可能性がある。実際には、ＡＥＤモデルが既にそのような変更を知らない場合、結果的に異常信号が生じ得る。これを抑制するために、ファジー・ネットは状態の検査および抑制すべきタグのリストを使用する。他の状況において、ＰＣＡモデル、バルブ流量モデル、およびファジー・ネットにおけるタグを指定された期間一時的に無効化することができる。大抵の場合、オペレータが再起動するのを忘れないように、再起動を１２時間後に設定している。タグが長期間サービスから除外されている場合、無効化することができる。現在抑制されている事象のリストを図４３に示す。
【００６７】
事象の詳細のログ取得
上述のようなシステムから最大の利点を引き出すには、オペレータを訓練して、ＡＥＤシステムを日常的な業務プロセスに組み込むことが必要である。是正処置を講じる最終的な権限は依然としてオペレータにあるため、ＡＥＤの性能発揮および強化のためその入力が重要である。ＡＥＤ事象詳細を体系的に捕捉して精査してフィードバックを与えるために、オペレータに対しＡＥＤ事象書式が提示される。これらは、事象の記録を維持して、ＡＥＤの利点を評価するのに役立った。ＡＥＤが稼動した時点から、いくつかの重大な事象が捕捉されて、操作担当者向けに文書化された。サンプル書式を図４４に示す。
【００６８】
代替的なソリューションの方が良い場合−反復事象のための是正処置
特定の反復する問題が識別された場合、構築担当者は問題を解決するための、より良いプロセスが存在しないことを確認しなければならない。特に、構築担当者は異常事象検知アプリケーションの構築を試みる前に以下の点を確認しなければならない。
・  問題を永久に解決することができるか？往々にして、現場担当者が問題を調査して永久に解決する充分な時間が無かったために問題が存在する。当該問題に組織的な関心が向けられたならば永久的な解決策が見つかる場合が多い。これが最善のアプローチである。
・問題を直接測定することができるか？問題を検知するためのより信頼性が高い方法は、当該プロセスにおける問題を直接測定できるセンサーを設置することである。これを用いて、プロセス制御アプリケーションを介して問題を防止することも可能である。これが次善のアプローチである。
・異常動作に対するアプローチを測定する推論的測定法を構築することができるか？推論的測定法は、ＰＣＡ異常事象モデルに極めて近い。信頼性を以って問題状態へのアプローチを測定するために利用可能なデータ（例：デルタ圧力を用いるタワーフラッディング）が存在する場合、これを状態の存在を検知するためだけでなく、状態が生じるのを防止する制御アプリケーションの基礎としても用いることができる。これが三番目に良いアプローチである。
【００６９】
異常事象検知アプリケーションは警報システムを代替しない
プロセス問題が急に生じる場合は常に、警報システムが異常事象検知アプリケーションと同程度に速く問題を識別する。オペレータが原因を診断するのを支援するために、事象のシーケンス（例：測定値が例外的である順序）が警報の順序よりも役に立つ場合がある。アプリケーションがオンラインになったならば、この可能性を調べるべきである。
【００７０】
しかし、異常事象検知のアプリケーションは、異常事象が緩慢に（１５分超）拡大する場合、オペレータに高機能の注意喚起を与えることができる。これらのアプリケーションはプロセスデータのパターンの変化に感応するものであり、単一の変数の大きな変動を必要としない。従って、警報を回避することができる。警報システムが、プロセスが狭い動作領域から逸脱した場合にオペレータに対して警告する（真の安全警報ではない）ように構成されている場合、本アプリケーションはこれらの警報を代替することが可能である。
【００７１】
単に異常事象の存在を検知することに加え、ＡＥＤシステムはまた、オペレータが事象を調査できるよう逸脱したセンサーを切り離す。現代的な工場は数千個のセンサーを備えており、人間がその全てをオンラインで監視することは不可能なことを考慮すればこれは重要な利点である。ＡＥＤシステムはこのように、異常な状況を効率的且つ効果的に取扱うべくオペレータ・ツールキットへの別の強力な付加機能とみなすことができる。
【００７２】
付録１
様々な規模の事象や撹乱が常にプロセス動作に影響を及ぼしている。ほとんどの場合、これらの事象や撹乱はプロセス制御システムにより処置される。しかし、プロセス制御システムがプロセス事象を適切に処置できない場合は常にオペレータがプロセス動作へ想定外の介入をする必要がある。発明者は、この状況を異常動作と定義し、その原因を異常事象と定義する。
【００７３】
異常動作を検知して、オペレータが根本原因の箇所を切り離すのを支援するために用いるモデルの組を生成してオンライン配備するための方法論およびシステムが開発されている。好適な実施形態において、モデルは主成分分析（ＰＣＡ）を用いる。これらのモデルの組は、既知の工学関係を表わす簡単なモデルおよび履歴データベースに存在する正常データパターンを表わす主成分分析法（ＰＣＡ）モデルの両方で構成されている。これら多くのモデル計算の結果が組み合わされて、プロセスが異常動作に入りつつあるか否かをプロセスオペレータが容易に監視できるようにする少数の要約時間傾向となる。
【００７４】
図１に、オンライン・システム内の情報が各種の変換、モデル計算、ファジー・ペトリネット、および圧密処理を経由して流れ、プロセス領域の正常性／異常性を示す概略傾向に到る様子を示す。本システムの核心は、プロセス動作の正常性を監視するために用いる各種のモデルである。
【００７５】
本発明に記載するＰＣＡモデルは、連続的な精製および化学プロセスを広範に監視し、設備およびプロセスにおいて生じつつある問題を迅速に検知することを目的とする。その意図は、特定のコンソールのオペレータ・ポストの責任権限の下で、全てのプロセス設備およびプロセス動作の包括的な監視機能を提供することである。これには、数百〜数千のプロセス測定値を有する多くの大規模な精製または化学プロセス動作ユニット（例：蒸留塔、反応炉、圧縮機、熱交換トレイン等）が含まれていてよい。監視機能は、長期にわたる性能低下ではなく、数分〜数時間の時間尺度で広がる問題を検知するように設計されている。プロセスおよび設備の問題は事前に特定されている必要がない。これは、文献に引用されている、特定の重要なプロセス問題を検知してプロセス動作のはるかに小さい部分を包含すべく構築されたＰＣＡモデルを利用するのとは対照的である。
【００７６】
この目的を達成すべく、ＰＣＡモデルを構築して配備する方法は、以下を含む連続的な精製および化学プロセスへのそれらのアプリケーションのために必要とされる多くの新規な機能拡張を含んでいる。
・測定値入力を選択、解析、変換するためのＰＣＡモデル基準および方法の設備範囲を確立するための基準
・主成分モデル、ＰＣＡの変動に基づく多変数統計モデルの開発
・付随する統計指標を再構成する簡単な工学関係に基づくモデルの構築
・例外計算および連続的なオンラインモデル更新を提供するオンラインデータの前処理
・ファジー・ペトリネットを用いた、モデル指標が正常か異常かの解釈
・ファジー・ペトリネットを用いた、多数のモデル出力の、プロセス領域の正常性／異常性を示す単一の連続的な概略傾向への組合せ
・操作および保守作業を反映させるため、モデルおよびファジー・ペトリネットとのオペレータの対話環境設計。
【００７７】
これらの拡張は、ＰＣＡおよび簡単な工学モデルを効果的に用いることができるたよう、連続的な精製および化学工場の操業の特徴、および対応するデータ特徴を扱うために必要である。これらの拡張は、第Ｉ種および第ＩＩ種の誤りの多くを防止し、且つ異常事象をより速く示すという利点をもたらす。
【００７８】
本節ではＰＣＡに一般的背景を記述しない。そのためには、非特許文献２などの標準的な教科書を参照されたい。
【００７９】
古典的なＰＣＡ技術では、以下の統計的な仮定をするが、その多くは正常連続的な精製および化学工場のプロセス動作から発生するデータによりある程度破られる。
１．プロセスは静的、即ちその平均および分散は時間に対して不変である。
２．変数同士の相互相関は、正常プロセス動作にわたり線形である
３．プロセスノイズ確率変数は相互に独立している。
４．プロセス変数の共分散行列は縮退していない（即ち半正定値）。
５．データは、「適切に」スケーリングされる（標準的な統計手法は単位分散へのスケーリング）。
６．（未補償）プロセスダイナミクス（このための標準的な部分補償は当該モデルの遅延変数の包含すること）が存在しない。
７．全ての変数がある程度の相互相関を有する。
８．データは多変数正規分布を示す。
【００８０】
その結果、ＰＣＡモデルを構築する際の入力および後続段階の選択、解析および変換において、違反の程度を評価して補償するため各種の調整が実行される。
【００８１】
これらのＰＣＡモデルがオンライン配備されたならば、モデル計算は、既知の動作の効果および保守作業の影響を除去し、故障したかまたは「不正な動作をしている」入力を無効にし、プロセスを通じた事象の伝播をオペレータが観察して承認することを可能にし、プロセスが正常に戻ったならば自動的に計算を再開すべく、特定の例外処理を必要とする。
【００８２】
ＰＣＡモデルの利用は、正常動作の間は真でなければならない既知の工学関係に基づく単純な冗長性検査により補われる。これらは、物理的に冗長な測定値を検査するような単純なものであってもよく、或いは物質収支と工学収支を検査するような複雑なものであってもよい。
【００８３】
冗長性検査の最も簡単な形式は、例えば以下のような簡単な２×２検査である。
・温度１＝温度２
・流量１＝バルブ特性曲線１（バルブ１の位置）
・処理ユニット１への物質流量＝処理ユニット１からの物質流量
【００８４】
これらは、図２のバルブ流量プロット等の簡単なｘ−ｙプロットとしてオペレータに提示される。各々のプロットには正常動作の領域があり、このプロット線では灰色の部分で示す。この領域の外側の動作は異常であるとして信号が送られる。
【００８５】
多重冗長性もまた、単一の多次元モデルにより調べることができる。多次元冗長性の例として以下のものがある。
・圧力１＝圧力２＝．．．．＝圧力ｎ
・処理ユニット１への物質流量＝処理ユニット１からの物質流量＝．．．＝処理ユニット２への物質流量
【００８６】
多次元検査は、「ＰＣＡ的」モデルにより表わされる。図３において、３種の独立且つ冗長な尺度Ｘ１、Ｘ２、およびＸ３がある。Ｘ３が１ずつ変化する都度、Ｘ１はａ１３ずつ変化し、またＸ２はａ２３ずつ変化する。この関係の組は、単一の主成分方向（Ｐ）を有するＰＣＡモデルとして表わされる。この種のモデルは、広範なＰＣＡモデルと同様の方法でオペレータに提示される。２次元の冗長性検査と同様に、灰色の部分は正常動作の領域を示す。Ｐの主成分負荷量は、変動性が最大の方向からＰを判定する従来の方法ではなく、工学方程式から直接計算される。
【００８７】
プロセス動作の特徴は、例外動作がこれらの関係を、プロセス動作の正常な範囲、正常領域での装置の変化、および保守作業を通じて正確に保つことを必要とする。例外動作の例として以下のものがある。
・流量計周辺のバイパス弁の開口
・上流の補償／下流の圧力変化
・領域測定値の再カリブレーション
・動作モードに基づくプロセス流の方向変換
【００８８】
ファジー・ペトリネットを用いてＰＣＡモデルおよび工学冗長性検査を組み合わせて、プロセスオペレータに対し、自分の制御下にあるプロセス動作が正常であることを連続的に要約提示する（図４）。
【００８９】
各々のＰＣＡモデルから複数の統計指標が作成され、これら指標をプロセスオペレータが扱うプロセス装置の構成および階層に対応付ける。従来の２乗予測誤差ＳＰＥの和による指標の感度は、ＰＣＡモデルが対応する完全なプロセス領域の指定された部分からの入力に対するＳＰＥ指標への寄与分だけを含むサブセット指標を生成することにより向上する。ＰＣＡモデルからの各々の統計指標はファジー・ペトリネットへ送られて当該指標を０〜１の縮尺に変換され、正常動作（値ゼロ）から異常動作（値１）までの範囲を連続的に示す。
【００９０】
各々の冗長性検査もまた、ファジー・ネットを用いて、連続的な正常−異常指標に変換される。これらのモデルが異常性を表示するために用いる２種の異なる指標がある。即ちモデルからの偏差および動作範囲外の偏差である（図３に示す）。これらの偏差は、誤差の２乗とＰＣＡモデルのホテリングＴ２乗指標との和に等しい。２より大きい次元を検査する場合、どの入力に問題があるかを識別することが可能である。図３において、Ｘ３−Ｘ２関係が依然として正常エンベロープ内に含まれるため、問題は入力Ｘ１にある。各々の偏差尺度は、ファジー・ペトリネットにより０〜１の縮尺に変換されて、正常動作（値ゼロ）から異常動作（値１）までの範囲を連続的に示す。
【００９１】
オペレータの権限下にある各々のプロセス領域の場合、適用可能な正常−異常指標の組を単一の正常−異常指標に組み合わせる。これは、ファジーペトリ論理を用いて異常動作の最悪ケースの表示を選択して行なわれる。このようにして、オペレータはプロセス領域内の全ての検査の高水準な要約を得られる。本節はファジー・ペトリネットの一般的背景を記述しない。そのためには、非特許文献１を参照されたい。
【００９２】
異常事象へのアプリケーションを構築するための全体プロセスを図５に示す。基本的な構築戦略は反復的であり、構築担当者は粗いモデルから出発し、次いで、正常動作および異常動作の両方の間、当該モデルがどの程度実際のプロセス動作を表わすかの観察に基づいて、当該モデルの能力を逐次向上させる。次いで、モデルはこれらの観察に基づいて再構築および再訓練される。
【００９３】
異常事象検知用のＰＣＡモデルの構築
Ｉ．概念的ＰＣＡモデル設計
全体的な設計目的は以下の通りである。
・コンソールのオペレータが自分の操作権限下にある全ての処理ユニットについてプロセス動作の連続的な状態（正常対異常）を把握できるようにする。
・オペレータが自身の操作権限内で急速（数分から数時間）に拡大しつつある異常事象の早期発見できるようにする。
・オペレータに対し、異常事象の根本原因を診断するために必要な主要プロセスの情報だけを提供する。
【００９４】
実際の根本原因診断は本発明の範囲外である。コンソールのオペレータは、自身のプロセス知識および訓練に基づいてプロセス問題が診断することが求められる。
【００９５】
広範なプロセス範囲を有することが、異常動作監視の全体的な成功にとって重要である。オペレータがシステムを学んで自身の技術を維持するために、定期的にシステムを用いることが必要である。特定の異常事象が稀にしか生じないため、プロセスの小さい部分の異常動作監視は稀にしかオペレータが使わず、恐らくは遂に異常事象を検知した際にオペレータにシステムを無視させる。この広範な範囲は、相当に経済的な関心事である特定のプロセス問題の検知に基づいてモデルを設計する、公開されているモデリングの目的とは対照的である（非特許文献３を参照）。
【００９６】
数千のプロセス測定値が、単一コンソールのオペレータの操作権限下の処理ユニット内にある。連続的な精製および化学プロセスは、これらの測定値中に顕著な時間ダイナミクスを示し、データ同士の相互相関を破る。このことは、相互相関を維持できる個別のＰＣＡモデルにプロセス装置を分割することを必要とする。
【００９７】
概念上のモデル設計は、以下の４種の主な決定からなる。
・プロセス装置を、対応するＰＣＡモデルを有する設備グループに再分割する。
・プロセス動作期間を、異なるＰＣＡモデルを必要とするプロセス動作モードに再分割する。
・設備グループ内のどの測定値を各ＰＣＡモデルへの入力として指定すべきかを識別する。
・設備グループ内のどの測定値が既知の事象または他の例外動作を抑制するフラグとして機能すべきかを識別する。
【００９８】
Ａ．処理ユニットの対応範囲
最初の意思決定は、単一ＰＣＡモデルで対応する設備のグループを形成することである。これに含まれる特定の処理ユニットは、プロセス統合／相互作用を理解している必要がある。多変数拘束コントローラの設計と同様に、ＰＣＡモデルの境界は、全ての重要なプロセス相互作用およびプロセスの変動や撹乱の主要な上流および下流における兆候を包含していなければならない。
【００９９】
以下の規則を用いてこれらの設備グループを判定する。
【０１００】
設備グループは、同一設備グループ内で全ての主要な物質とエネルギーの融合、およびの迅速な再利用を含むことにより定義される。プロセスが多変数拘束コントローラを用いる場合、制御モデルは処理ユニット同士の相互作用位置を明示的に識別する。さもなければ、相互作用は当該プロセスの工学解析を通じて識別する必要がある。
【０１０１】
プロセスグループは、プロセス装置グループ同士の相互作用が最小である点で分割されなければならない。最も明確な分割点が生じるのは、相互作用が次の下流ユニットへの供給を含む単一のパイプを介してのみ生じる場合である。この場合、温度、圧力、流量、および供給の組成が下流の設備グループに対する一次影響であり、直ぐ下流にあるユニットの圧力は上流の設備グループへの一次影響である。これらの一次影響の測定値は、上流および下流の両方の設備グループのＰＣＡモデルに含まれているべきである。
【０１０２】
上流と下流の設備グループの間にプロセス制御アプリケーションの影響を含める。プロセス制御アプリケーションは、上流で下流の設備グループの間に追加的な影響経路を提供する。フィードフォワードおよびフィードバック経路の両方が存在し得る。このような経路が存在する場合、これらの経路を駆動する測定値は両方の設備グループに含まれていなければならない。プロセス制御アプリケーションの解析により、処理ユニット同士の主要な相互作用が示される。
【０１０３】
長時間動特製（ｓｉｇｎｉｆｉｃａｎｔ  ｔｉｍｅ  ｄｙｎａｍｉｃｓ）が存在する場所（例：貯蔵タンク、長いパイプライン等）は設備グループを分割する。ＰＣＡモデルは基本的に、急速なプロセス変動（例：数分〜数時間の期間にわたり生じるもの）を扱う。プロセスに影響を及ぼすのに数時間、数日、または数週間もかかる影響はＰＣＡモデルに適していない。これらの影響が正常なデータパターンにとって重要な場合、これらの影響の測定値を動的に補償して、それらの有効時間を他のプロセス測定値と同期化させる必要がある（動的補償についての説明を参照）。
【０１０４】
Ｂ．プロセス動作モード
プロセス動作モードは、プロセスの挙動が大幅に異なる特定の期間として定義される。これらの例として、異なる等級の製品（例：ポリマー生産）の生産、顕著なプロセス遷移（例：始動、停止、原料の交換）、激変した原料（例：オレフィン生産においてエタンに代えてナフサを分解）の処理、または異なる構成のプロセス設備（異なる組の処理ユニットが稼動）がある。
【０１０５】
これらの顕著な動作モードが存在する場合、各々の主要動作モードに対して個別のＰＣＡモデルが構築すれ必要があり得る。必要とされるモデルが少ない方が良い。構築担当者は、特定のＰＣＡモデルが類似の動作モードに対応可能であると仮定すべきである。この仮定は、各動作モードからの新規データをモデルに通して当該モデルが正しく動作するか否かを確認することにより検証されなければならない。
【０１０６】
Ｃ．履歴的プロセス問題
異常事象検知システムの構築に組織的な関心が持たれるのは、従来から経済的に重大な影響を及ぼすプロセス問題が存在しているためである。しかし、これらの重要な問題を分析して、当該問題に取り組む最適なアプローチを識別する必要がある。特に、構築担当者は、異常事象検知アプリケーションの構築を試みる前に以下の点を確認しなければならない。
１．問題を恒久的に解決することが可能か？往々にして、現場担当者が問題を調査して恒久的に解決する充分な時間が無かったために問題が存在する。当該問題に組織的な関心が向けられたならば恒久的な解決策が見つかる場合が多い。これが最善のアプローチである。
２．問題を直接測定することが可能か？問題を検知するためのより信頼性が高い方法は、当該プロセスにおける問題を直接測定できるセンサーを設置することである。これを用いて、プロセス制御アプリケーションを介して問題を防止することも可能である。これが次善のアプローチである。
３．異常動作に対するアプローチを測定する推論的測定法を構築することが可能か？推論的測定法は通常、部分最小二乗法、ＰＬＳ、ＰＣＡ異常事象モデルに極めて近いモデルを用いて構築される。推論的測定法を構築する他の一般的な代替案として、ニューラルネットワークおよび線形回帰モデルが含まれる。信頼性を以って問題状態へのアプローチを測定するために利用可能なデータ（例：デルタ圧力を用いるタワーフラッディング）が存在する場合、これを状態の存在を検知するためだけでなく、状態が生じるのを防止する制御アプリケーションの基礎ベースとしても用いることができる。これが三番目に良いアプローチである。
【０１０７】
問題状態の直接測定、およびこれらの状態の推論的測定値は、異常検知モデルの全体的なネットワークに容易に一体化することが可能である。
【０１０８】
ＩＩ．入力データおよび動作範囲の選択
設備グループ内に何千ものプロセス測定値がある。予備設計では以下を行なう。
・全てのカスケード２次コントローラ測定値、および特にこれらのユニットへの最終２次出力（フィールド制御バルブへの信号）を選択する
・コンソールのオペレータがプロセスを監視するのに用いる主要な測定値（例：自身の操作系統に現われるもの）を選択する
・担当技術者がプロセスの性能を測定するために用いる何らかの測定値を選択する
・供給速度、供給温度、または供給品質の何らかの上流測定値を選択する
・プロセスの動作領域、特に圧力に影響を及ぼす下流状態の測定値を選択する
・重要な測定のために追加的な冗長測定値を選択する
・非線形変換の計算に必要となるかもしれない測定値を選択する
・撹乱（例：周囲温度）の何らかの外部測定値を選択する
・プロセス状態の重要な尺度であるとプロセスの専門家が考える他のあらゆる測定値を選択する
【０１０９】
上記のリストから以下の特性を有する測定値のみを含める。
・測定値に誤りまたは問題がある性能の履歴が含まれていない
・測定値は満足すべき信号対ノイズ比を有する
・測定値はデータ集合内の他の測定値と相互相関を求められている
・測定値は正常動作期間の１０％を超えて飽和していない
・測定値は、稀にしか変化しない固定セットポイントに密接には制御されていない（制御階層の最終一次側）
・測定値は、長期にわたる「不良値」動作、またはトランスミッタの限界まで飽和していない
・測定値は、極めて非線形であることが知られている値の範囲を超えない
・測定値は未加工測定値からの冗長計算ではない
・フィールド制御バルブへの信号は、時間の１０％を超えて飽和していない
【０１１０】
Ａ．モデル入力を選択するための評価
ＰＣＡ異常検知モデル、信号対ノイズ比、および相互相関への潜在的入力に優先順位付けを行なうための２種の統計的な基準がある。
【０１１１】
１）信号対ノイズテスト
信号対ノイズ比は、入力信号における情報内容の尺度である。
【０１１２】
信号対ノイズ比は、以下のように計算される。
１．未加工信号は、プロセスのものに等しい近似動的時定数を有する指数フィルタを用いてフィルタリングされる。連続的な精製および化学プロセスの場合、この時定数は通常、３０分〜２時間の範囲にある。他の低域通過フィルタを用いてもよい。指数フィルタの場合の式は以下の通りである。
Ｙｎ＝Ｐ＊Ｙｎ−１＋（１−Ｐ）＊Ｘｎ  指標フィルタの式    式１
Ｐ＝Ｅｘｐ（−Ｔｓ／Ｔｆ）    フィルタ定数の計算    式２
ここで、
Ｙｎ  現在フィルタリングされる値
Ｙｎ−１  以前にフィルタリングされた値
Ｘｎ  現在の未加工値
Ｐ  指数フィルタ定数
Ｔｓ  測定値のサンプル時間
Ｔｆ  フィルタ時定数
２．残留信号は、未加工信号からフィルタリングされた信号を減算することにより得られる。
Ｒｎ＝Ｘｎ−Ｙｎ    式３
３．信号対ノイズ比は、フィルタリングされた信号の標準偏差を残留信号の標準偏差で除算した比である。
Ｓ／Ｎ＝σＹ／σＲ    式４
【０１１３】
全ての入力が、所定の最小値（例：４）より大きいＳ／Ｎを示すことが好ましい。Ｓ／Ｎが当該最小値を下回る入力については、モデルに含めるべきか否かを判定するために個別に調べる必要がある。
【０１１４】
Ｓ／Ｎを計算するために用いるデータ集合は、ノイズ内容の推定を極端に長引かせるため、長期間にわたる定常状態動作は全て除外しなければならない。
【０１１５】
２）相互相関テスト
相互相関は、入力データ組の情報冗長性の尺度である。任意の２つの信号間の相互相関は次式により計算される。
１．各入力組ｉ、ｋ間の共分散Ｓｉｋを計算する。
【数１】
2008512800.tif 000003
２．共分散からの入力の各組について相関係数を計算する。
ＣＣｉｋ＝Ｓｉｋ／（Ｓｉｉ＊Ｓｋｋ）１／２    式６
【０１１６】
入力がモデルに含まれていてはならないことをフラグする二つの状況がある。第１の状況は生じるのは、特定の入力と残りの入力データ組との間に顕著な相互相関が存在しない場合である。各々の入力に対して、当該データ組において有意な相関係数、例えば０．４、を有する少なくとも１個の他の入力が存在しなければならない。
【０１１７】
第２の状況が生じるのは、異なる識別器を有する何らかの計算に際して同一入力情報が（偶然に）２回含まれた場合である。１に近い（例えば０．９５超）相関係数を示す一対の入力は全て、両方の入力をモデルに含めるべきか否か判定すべく個別の検査が必要である。入力が物理的には独立しているが論理的には冗長である（即ち、２個の独立熱電対が別々に同じプロセス温度を測定している）場合、これら両方の入力はモデルに含まれるべきである。
【０１１８】
２つの入力がお互いの変換（即ち、温度および圧力補償された温度）である場合、これらの測定値の１個と残りのデータ組との間で顕著に向上した相互相関が存在しない限り、オペレータが馴染んでいる測定値を含めることが好ましい。次いで、より高い相互相関を有する方を含めるべきである。
【０１１９】
３）飽和した変数の識別および扱い
精製および化学プロセスは往々にして厳しいか緩い拘束に遭遇して、結果的にモデル入力に対する飽和値および「不良値」が生じる。一般的な拘束は、計器トランスミッタの上下範囲、アナライザ範囲、最大および最小制御バルブ位置、およびプロセス制御アプリケーション出力限度である。入力は、モデル構築およびこれらのモデルのオンライン利用の両方において入力を事前処理する際に特別な扱いを要する飽和に関していくつかのカテゴリに分類することができる。
【０１２０】
標準的なアナログ計器（例：４〜２０ミリアンペアの電子トランスミッタ）において、不良値は以下の２種の別々の原因で生じる恐れがある。
・実際のプロセス状態がフィールドトランスミッタの範囲外にある
・フィールドとの接続が絶たれた
【０１２１】
これらの状況のいずれかが生じた場合、プロセス制御システムを個々の測定値に基づいて設定して、当該測定値が不良値であることを示す特別のコードを当該測定値に割り当てるか、或いは測定の最終良好値を維持する。これらの値は次いで、プロセス制御システムに対して実行されるあらゆる計算全体にわたり伝播する。「最終良好値」オプションが設定された場合、これは検知および除外が困難な誤った計算に至る恐れがある。「不良値」コードがシステム全体に伝播した場合に典型的に、不良測定値に依存する全ての計算も同様に不良としてフラグが立てられる。
【０１２２】
プロセス制御システムに設定されたオプションとは無関係に、不良値を含むそれらの期間は、訓練またはテストデータ集合に含まれてはならない。構築担当者は、プロセス制御システムにおいてどのオプションが設定されたかを識別し、次いで、不良値であるサンプルを除外するためのデータフィルタを設定する必要がある。オンライン実装の場合、プロセス制御システムでどのオプションが選択されたかとは無関係に、入力は事前処理して不良値を欠落値としてフラグ付けしなければならない。
【０１２３】
通常は長期間にわたり不良値であるこれらの入力をモデルから除外すべきである。
【０１２４】
拘束変数とは、測定値がある限界にあって、当該測定値が実際のプロセス状態（値がトランスミッタ範囲の最大または最小限にデフォルト化されている場合とは逆。「不良値」節で記述）に一致する。このプロセス状況は、以下のいくつかの理由により生じる。
・プロセスの各部分は、特に優先的な状況、例えばフレアシステムへの圧力解放流がない限り、通常は不活性である。これらの優先的状況が有効である時間は、データフィルタを設定することにより訓練および検証データ集合から除外しなければならない。オンライン実装の場合、これらの優先事象は、選択されたモデル統計量を自動抑制するトリガ事象である
・プロセス制御システムは、プロセス動作制限、例えば製品仕様限界、に反してプロセスを駆動すべく設計されている。これらの拘束は通常、２種のカテゴリに分類される。即ち時々飽和するものと、通常は飽和しているものである。通常は飽和している入力はモデルから除外しなければならない。稀に（例えば全時間の１０％未満）飽和するだけの入力はモデルに含まれてよいが、飽和していない場合には期間に基づいてスケーリングしなければならない。
【０１２５】
Ｂ．プロセス制御アプリケーションからの入力
プロセス制御アプリケーションは、プロセスデータの相互相関構造に極めて重大な影響を及ぼす。特に以下の通りである。
・被制御変数の変動が大幅に減少するため、プロセスに顕著なプロセス撹乱が出現したか、または、オペレータが主要なセットポイントを変えて意図的に動作点を動かした場合の短期間を除いて、被制御変数内の動きは基本的にノイズである。
・被制御変数における通常の変動は、制御システムにより被操作変数（最終的にはフィールド内の制御バルブへ送られた信号）へ転送される。
【０１２６】
精製および化学プロセスの正常動作は通常、２種類の異なる制御機構により制御される。即ち古典的制御カスケード（図６に示す）および、より最近の多変数拘束コントローラＭＶＣＣ（図７に示す）である。
【０１２７】
１）カスケード構造からモデル入力を選択
図６に、典型的な「カスケード」プロセス制御アプリケーションを示す。これは精製および化学プロセス用に極めて一般的な制御機構である。そのようなアプリケーションからの多くの潜在的なモデル入力が存在するが、当該モデルの候補である唯一のものは未加工プロセス測定値（本図の「ＰＶ」）およびフィールドバルブへの最終出力である。
【０１２８】
カスケード制御機構の最終一次側のＰＶは、極めて重要な測定値であるにもかかわらず、当該モデルへ含めるには適していない候補である。この測定値は通常、制御機構の目的がこの測定値をセットポイントに保つことであるため、動きが極めて限られる。自身のセットポイントが変えられた場合に最終一次側のＰＶ内に動きがあってもよいが、通常これは稀である。一次セットポイントの不定期的移動から生じるデータパターンは通常、モデルがデータパターンを特徴付けるに足る影響力を訓練データパターンに有していない。
【０１２９】
最終一次側のセットポイントの変化から生じるデータパターンを特徴付けることがこのように困難なため、オペレータが当該セットポイントを動かしたならば、当該モデルの指標である二乗予測誤差ＳＰＥの和が顕著に増加しそうである。従って、最終一次側のセットポイントがわずかに変化しても「既知事象抑制」を起こす候補トリガである。オペレータが最終一次側セットポイントを変える都度、「既知事象抑制」論理がＳＰＥ計算からその影響を自動的に除去する。
【０１３０】
構築担当者が最終一次側のＰＶをモデルに含めるのであれば、当該測定値は、オペレータがセットポイントを変えた時点からプロセスが新規セットポイントの値の近くへ移動するまでの短い期間に基づいてスケーリングすべきである（例えば、セットポイントが１０から１１まで変化した場合、ＰＶが１０．９５に達したときに新規セットポイントの９５％以内に）。
【０１３１】
最終一次側のＰＶと極めて強い相関（例えば相関係数が０．９５を超える）を有する測定値もまたあり得る。例えば、最終一次側のＰＶとして用いられる、温度測定の近傍に置かれた冗長な熱電対である。これらの冗長な測定値は、最終一次側のＰＶ用に選択されたのと同じ方法で扱われなければならない。
【０１３２】
カスケード構造は、各々の二次側に設定ポイント制限を持つことができ、フィールド制御バルブへの信号に対して出力制限を設けることができる。これらの潜在的に拘束された動作の状態を調べることは、セットポイントに関連付けられた測定値が拘束された方法で操作されたか否か、または、フィールドバルブに対する信号が拘束されたか否か見るために重要である。これらの拘束された動作の間の日付を用いてはならない。
【０１３３】
２）多変数拘束コントローラ（ＭＶＣＣ）からのモデル入力の選択／計算
図７に、精製および化学プロセス向けの極めて一般的な制御機構である典型的ＭＶＣＣプロセス制御アプリケーションを示す。ＭＶＣＣは動的数学モデルを用いて、被操作変数ＭＶ（通常は調整制御ループのバルブ位置またはセットポイント）の変化がどのように制御変数ＣＶ（プロセス状態を測定する従属温度、圧力、組成および流量）を変えるかを予測する。ＭＶＣＣは、プロセス動作を動作限度まで押し上げようと試みる。これらの限度は、ＭＶ限度またはＣＶ限度のいずれかであってよく、外部オプティマイザにより決定される。プロセス動作の限度の数は、コントローラが操作可能なＭＶの数から、制御されている物質収支の数を差し引いたものに等しい。従って、ＭＶＣＣが１２個のＭＶ、３０個のＣＶ、および２個のレベルを有する場合、プロセスは１０個の限度に向けて動作される。ＭＶＣＣはまた、プロセスに対する測定された負荷変動の影響を予測して、これらの負荷変動（フィードフォワード変数ＦＦとして知られる）を補償する。
【０１３４】
未加工ＭＶまたはＣＶが、ＰＣＡモデルに含めるべき良い候補であるか否かは、ＭＶまたはＣＶがＭＶＣＣによる動作限度に対して保持されている時間の割合に依存する。「拘束変数」の節で議論したように、時間の１０％を超えて拘束される未加工変数は、ＰＣＡモデルへ含める候補としては適していない。通常、非拘束変数は「拘束変数」節の議論に従って扱われなければならない。
【０１３５】
無拘束ＭＶが調整制御ループに対するセットポイントである場合、セットポイントを含めるべきではなく、代わりに当該調整制御ループの測定値を含めるべきである。また、当該調整制御ループからのフィールドバルブに対する信号も含めるべきである。
【０１３６】
無拘束ＭＶがフィールドバルブ位置に対する信号である場合、これをモデルに含めるべきである。
【０１３７】
Ｃ．冗長測定値
プロセス制御システムデータベースは、ＰＣＡモデルへの候補入力の中で顕著な冗長性を有する可能性がある。冗長性の一種として「物理的冗長性」があり、これはプロセス設備内で互いに物理的に近接して配置された多数のセンサー（熱電対等）が存在する場合である。別の種類の冗長性として「計算的冗長性」があり、これは未加工センサーが新たな変数（例：圧力補償された温度、または容積測定流量測定値から計算されて質量流量）に数学的に組み合わされた場合である。
【０１３８】
一般的なルールとして、未加工測定値および当該測定値から計算された入力は共に当該モデルに含まれてはならない。一般に好適なのは、プロセスオペレータが最も馴染んでいる測定値のバージョンを含めることである。このルールにおける例外は、当該モデル用にデータの相互相関構造を向上させるために未加工入力を数学的に変換しなければならない場合である。その場合、未加工測定ではなく、変換された変数を当該モデルに含めるべきである。
【０１３９】
物理的冗長性は、モデルの相互検証情報を提供するために極めて重要である。一般的なルールとして、物理的に冗長な未加工測定値はモデルに含めるべきである。多数の物理的に冗長な測定値が存在する場合、これらの測定値は、主成分（可変スケーリングに関する節を参照）の選択を阻害するのを防止すべく特別にスケーリングされなければならない。一般的なプロセスの例は、原子炉暴走を捕捉すべく原子炉に設置された多数の熱電対から生じる。
【０１４０】
極めて巨大なデータベースをマイニングする場合、構築担当者は全ての候補入力の中で相互相関計算を行なうことにより冗長な測定値を識別することができる。相互相関が極めて高い（例えば０．９５超）測定値ペアを個別に調べて各々の組を物理的に冗長または計算的に冗長のいずれかに分類しなければならない。
【０１４１】
ＩＩＩ．履歴データの収集
構築に要する努力の多くは、正常なプロセス動作の全てのモードを含むことがわかっている良好な訓練データ集合の生成にある。このデータ集合は以下を満たさなければならない。
【０１４２】
正常動作範囲にわたること。動作範囲の小さい部分にわたるデータ集合は殆どがノイズからなる。定常状態動作の間のデータの範囲に比べたデータの範囲は、データ集合の情報の品質を良く示している。
【０１４３】
全ての正常動作モード（季節的なモード変動を含む）を含んでいること。各々の動作モードは、異なる相互相関構造を有していてよい。動作モードを特徴付けるパターンがモデルにより捕捉されない限り、これらの非モデル化動作モードは異常動作として現われる。
【０１４４】
正常動作データだけを含んでいること。強い異常動作データが訓練データに含まれている場合、モデルは誤ってこれらの異常動作を正常動作としてモデル化してしまう。その結果、後で当該モデルを異常動作と比較した際に、異常動作を検知することができない。
【０１４５】
履歴は、オンライン・システムで用いるデータになるべく類似しているべきである。オンライン・システムは、異常事象を検知するのに充分速い周期でスポット値を提供する。連続的な精製および化学的操作の場合、このサンプリング周期は約１分である。データ履歴管理機能の限度内で、訓練データは可能な限り１分のスポット値に等価でなければならない。
【０１４６】
データ収集の戦略は、長期間（通常９〜１８ヶ月間の範囲）の動作履歴から開始し、これらの期間から明白な、または文書化された異常事象を除外するように試みることである。そのような長期間を用いることにより以下が実現する。
・比較的小さい異常事象は、モデルパラメータに大きく影響する程度には充分な強さを以って訓練データ集合に現われない
・大部分の動作モードが発生して、データ内に現われているはずである。
【０１４７】
Ａ．履歴データ収集に関する事項
１）データ圧縮
多くの履歴データベースはデータ圧縮を用いて、必要なデータ記憶容量を最小化する。不都合なことに、これを実行すればデータの相互相関構造を乱す恐れがある。プロジェクトの開始時点において、データベースのデータ圧縮を無効にしてデータのスポット値の履歴を取っておくべきである。可能な場合は常に非圧縮データを用いて最終的なモデルを構築しなければならない。平均化された値は、それらが利用できる唯一のデータであって、利用できる最短のデータ平均でない限り、用いるべきではない。
【０１４８】
２）データ履歴の長さ
モデルが正常なプロセスパターンを適切に表わすために、訓練データ集合は、全ての正常動作モード、正常な動作変動、およびプロセスに生じる変化や正常な微小攪乱の例を含んでいなければならない。これは、長期間（例：９〜１８ヶ月）にわたるプロセス動作のデータを用いることにより実現される。特に、精製および化学プロセスにおいて季節的（春夏秋冬）な動作の違いは極めて重要であり得る。
【０１４９】
時折、これらの長い範囲のデータが未だ利用できない場合（例：プロセス設備の切り替えその他の重要な再構成の後）がある。これらの場合、モデルは訓練データの短い（例：６週間）初期集合から始めて、更なるデータが収集されるにつれて訓練データ集合が拡張されて、モデルが安定するまで当該モデルが毎月更新される（例：モデル係数は新規データを追加しても変化しない）。
【０１５０】
３）補助履歴データ
この期間に対する各種の操作ジャーナルもまた収集しなければならない。これを用いて、動作期間を異常とみなすが、または訓練データ集合から除外する必要がある何らかの特別なモードにおいて動作している。特に、重要な履歴的異常事象をこれらのログから選択してモデル用のテストケースとすることができる。
【０１５１】
４）特定の測定履歴の欠如
往々にして、セットポイントおよびコントローラ出力は、工場プロセスデータの履歴管理機能で履歴化されないことが多い。これらの値の履歴化は、プロジェクトの開始時点で直ちに開始しなければならない。
【０１５２】
５）動作モード
もはや適切に現在のプロセス動作を表わさない旧データは、訓練データ集合から削除しなければならない。大規模なプロセス変更の後で、訓練データおよびＰＣＡモデルを最初から作り直すことが必要な場合がある。特定の種類の動作がもはや実行されていない場合、当該動作からの全てのデータを訓練データ集合から削除しなければならない。
【０１５３】
動作ログを用いて、異なる動作モードの下でいつプロセスが稼動したかを識別しなければならない。これらの異なるモードは別々のモデルを必要とする場合がある。モデルがいくつかの動作モードに対応することを意図している場合、各々の動作モデルからの訓練データ集合内のサンプルの数はほぼ等しくなければならない。
【０１５４】
６）サンプリングレート
構築担当者は、拠点のプロセス履歴管理機能を使用して数ヶ月分のプロセスデータ、好ましくは１分間隔のスポット値、を収集しなければならない。これが入手可能でない場合、なるべく平均化されておらず、且つ解像度が最も高いデータを使用すべきである。
【０１５５】
７）稀にサンプリングされる測定値
品質測定値（アナライザおよび研究室サンプル）は他のプロセス測定値よりも極めて長いサンプリング周期を有し、数十分毎〜１日１回の範囲である。これらの測定値をモデルに含めるには、これらの品質測定値の連続的な推定値を生成する必要がある。図８に、連続的な品質推定値のオンライン計算を示す。これと同じモデル構造を構築して履歴データに適用しなければならない。この品質推定値は次いで、ＰＣＡモデルへの入力になる。
【０１５６】
８）モデルにより起動されるデータ注釈
極めて明らかな異常を除いて、履歴データの品質を判定するのは困難である。異常動作データを含めることによりモデルが偏る恐れがある。大量の履歴データを用いる戦略は、訓練データ集合における異常な動作により生じるモデルのバイアスをある程度補償する。プロジェクトの開始に先立って履歴データから構築されたモデルは、品質に関しては疑いを持たなければならない。初期訓練データ集合は、プロジェクトが継続している間に生じるプロセス状況の高い品質注釈を含むデータ集合により代替されなければならない。
【０１５７】
モデル構築の戦略は、初期の「粗い」モデル（疑わしい訓練データ集合の結果）から出発し、次いでモデルを用いて高品質訓練データ集合の収集を起動することである。モデルを用いてプロセスを監視するに従い、正常動作、特別動作、および異常動作に関する注釈およびデータが集められる。モデルが異常動作のフラグを立てるか、または異常事象をモデルが見逃した都度、事象の原因および持続期間に注釈が付けられる。このように、プロセス動作を監視するモデルの能力に対するフィードバックを訓練データに取り込むことができる。次いでこのデータを用いてモデルを改良し、次いで当該モデルを用いて品質が向上した訓練データを収集し続ける。このプロセスは、モデルが満足すべき状態になるまで繰り返される。
【０１５８】
ＩＶ．データ及びプロセス解析
Ａ．初期の粗いデータ解析
動作ログを用いて、プロセスの主要な性能指標を調べることにより、履歴データを、異常動作の存在が判明している期間と、識別された異常動作が存在しない期間に分ける。識別された異常動作が存在しないデータが訓練データ集合となる。
【０１５９】
ここで、各測定値は、訓練データ集合の候補であるか否かを判定すべく、自身の全履歴を調べられる必要がある。除外しなければならない測定値は以下の通りである。
・「不良値」である期間が長いもの
・自身のトランスミッタの上限または下限に固定された期間が長いもの
・極めて僅かな変動性しか示さないもの（自身のセットポイントに密に制御されているものを除く）
・自身の動作範囲に関して極めて大きい変動性を連続的に示すもの
・データ集合内の他のどの測定値とも殆ど相互相関を示さないもの
・信号対ノイズ比が低いもの
【０１６０】
データを調べる間、測定値が短期的に「不良値」を示すか、或いは自身のトランスミッタ上限または下限に短期的に固定されている期間もまた除外されなければならない。
【０１６１】
これらの除外がなされたならば、第１の粗いＰＣＡモデルが構築される筈である。これは極めて粗いモデルになるため、保持される主成分の正確な個数は重要でない。これは通常、モデルに含まれる測定値の個数の約５％である。ＰＣの数は最終的に、プロセス内の自由度の個数に合致する筈であるが、これはプロセス撹乱の全ての異なる発生源を含んでいるため、通常は未知である。含めるべき主成分の個数を判定するいくつかの標準的な方法がある。また、この段階では、可変スケーリングへの統計的アプローチを用いるべきである。即ち全ての変数を単位分散にスケーリングする。
Ｘ’＝（Ｘ−Ｘａｖｇ）／σ    式７
【０１６２】
ここで訓練データ集合をこの事前モデルに通して、データがモデルに合致しない期間を識別しなければならない。これらの期間を調べて、その時点で異常事象が生じていたか否かを判定しなければならない。そうであると判定されたならば、これらの期間はまた、既知の異常事象が生じている時間としてフラグを立てられなければならない。これらの期間を訓練データ集合から除外し、修正されたデータによりモデルを再構築しなければならない。
【０１６３】
Ｂ．外れ値および異常動作期間の除外
以下を通じて、明らかな異常事象を除去する。
【０１６４】
文書化された事象の除外。拠点において異常事象の履歴の完全な記録が残されていることは極めて珍しい。しかし、顕著な動作問題は、オペレータのログ、オペレータの変更ジャーナル、警報ジャーナル、および計器保守記録等の動作記録に文書化されなければならない。これらは、異常事象の履歴の部分的な記録しか提供しない。
【０１６５】
主要性能指標ＫＰＩが異常である期間の除外。供給レート、製品レート、製品品質等の測定値が共通の主要性能指標である。各々のプロセス動作には、装置に固有の追加的なＫＰＩがあってよい。この限られた測定値の集合を注意深く調べることにより、異常動作の期間が明確に示される。図９に、ＫＰＩのヒストグラムを示す。このＫＰＩの動作目的はこれを最大化することであるため、このＫＰＩが低い場合の動作期間は異常動作である可能性がある。プロセス品質は、最適動作が仕様限界以内であって、正常な供給速度の変動に対する感度が低いため、往々にして最も解析が容易なＫＰＩである。
【０１６６】
Ｃ．ノイズの補償
ノイズとは、プロセスに関して有用な情報を含んでいない測定信号の高周波内容を指す。ノイズは、オリフィス・プレートを横断する二相流またはレベルの撹乱等の特定のプロセス状況により生じ得る。ノイズは、電気インダクタンスにより生じ得る。しかし、恐らくはプロセス撹乱により生じた顕著なプロセス変動性は有用な情報にあって、フィルタリングで除去してはならない。
【０１６７】
精製および化学的プロセスの測定で生じる２種の一次ノイズがある。即ち測定スパイクおよび指数相関を有する連続ノイズである。測定スパイクにより、信号は、自身の以前の値に近い値に戻る前に、サンプル数が少ない割に不合理なほど大幅に飛ぶ。ノイズ・スパイクは、ユニオン・フィルタ等の従来型スパイク除外フィルタを用いて除外される。
【０１６８】
信号におけるノイズの量は、信号対ノイズ比として知られる尺度により定量化できる（図１０参照）。これは、高周波ノイズに起因する信号変動性の程度に対するプロセス変動に起因する信号変動性の程度の比として定義される。４未満の値が、信号が顕著なノイズを含んでいてモデルの効果を阻害し得ることを示す典型的な値である。
【０１６９】
構築担当者が顕著なノイズを含む信号に遭遇するたびに、三種の選択の一つを行なう必要がある。好ましい順に以下の通りである。
・ノイズの発生源を除去することにより信号を固定する（最良の対応策）
・フィルタリング技術を用いてノイズを除去／最小化する
・信号を当該モデルから除去する
【０１７０】
信号対ノイズ比が２〜４の間にある信号では通常、指数相関を有する連続ノイズは、指数フィルタ等の従来型低域通過フィルタにより除去することができる。指数フィルタの式は以下の通りである。
Ｙｎ＝Ｐ＊Ｙｎ−１＋（１−Ｐ）＊Ｘｎ    指数フィルタの式    式８
Ｐ＝Ｅｘｐ（−Ｔｓ／Ｔｆ）    フィルタ定数の計算    式９
Ｙｎは現在のフィルタリングされた値
Ｙｎ−１は以前にフィルタリングされた値
Ｘｎは現在の未加工値
Ｐは指数フィルタ定数
ＴＳは測定値のサンプル時間
Ｔｆフィルタが時定数
【０１７１】
信号対ノイズ比が極めて低い（例えば２未満の）信号は、直接モデルに含まれる程度にはフィルタリング技術によっても充分に改善できないであろう。入力が重要と見なされる場合、変数のスケーリングはスケーリング係数のサイズを大幅に増す（通常は２〜１０の範囲）ことにより、モデルの感度を下げるべく設定されなければならない。
【０１７２】
Ｄ．変換された変数
変換された変数は、二つの異なる理由によりモデルに含まれなければならない。
【０１７３】
第一に、特定の設備およびプロセス化学の工学的解析に基づいて、プロセス内で既知の非線形性を変換して当該モデルに含めなければならない。ＰＣＡの仮定の一つが、モデルの変数が線形に相関していることであるため、プロセスまたは設備の顕著な非線形性はこの相互相関構造を壊して、モデルからの偏差として現われる。これは、モデルの使用可能範囲に影響を及ぼす。
【０１７４】
既知の非線形変換の例として以下のものがある。
・蒸留塔における還流対供給比
・高純度蒸留の組成のログ
・圧力補償された温度測定
・副流の生成
・バルブ位置への流れ（図２）
・指標的温度変化への反応率
【０１７５】
第二に、過去に生じたプロセス問題からのデータもまた、これらの問題が当該プロセスの測定値にどのように現われるかを理解するために調べなければならない。例えば、装置のデルタ圧力と供給速度の関係は、フラッディング点に達するまでは比較的線形であるが、そこからデルタ圧力は指数的に増加する。タワーフラッディングは、この線形相関が破れたことで検出されるため、デルタ圧力および供給速度の両方が含まれなければならない。別の例として、触媒流問題が移送ラインのデルタ圧力にしばしば見られる。従って、絶対圧力測定値をモデルに含めるのではなく、デルタ圧力を計算して含めなければならない。
【０１７６】
Ｅ．動的変換
図１１に、プロセスダイナミクスが２個の測定値の現在値同士の相互相関をどのように乱すかを示す。遷移時間中、一方の値は常に変化しているが他方はそうでなく、従って遷移時間中は現在値同士の相互相関は存在しない。しかし、これら２個の測定値は、動的伝達関数を用いて主要な変数を変換することにより、時間同期に戻すことができる。データの時間同期化には通常、無駄時間動的モデル（式９にラプラス変換形式で示す）を有する第１オーダーで充分である。
【数２】
2008512800.tif 000004
Ｙ：未加工データ
Ｙ’：時間同期化データ
Ｔ：時定数
Θ：無駄時間
Ｓ：ラプラス変換パラメータ
【０１７７】
この技術が必要とされるのは、モデルで用いる変数同士に顕著な動的分離がある場合のみである。通常、変数の１〜２％だけがこの処理を必要とする。これは、オペレータによりしばしば大きな刻みで変えられるセットポイント等の独立変数、およびモデル化されている主処理ユニットのかなり上流にある測定値に当てはまる。
【０１７８】
Ｆ．平均的な動作点の除去
連続的な精製および化学プロセスは、ある動作点から別の動作点へ常に移動されている。これらは、オペレータまたは最適化プログラムが主要なセットポイントを変更した意図的なものであったてもよく、或いは、熱交換器の汚れや触媒の非活性化等の遅いプロセス変動により生じる場合がある。その結果、未加工データは静止していない。これらの動作点の変更は、静止データを生成するために除去する必要がある。さもなければ、これらの変更は異常事象として誤って現われる。
【０１７９】
プロセスの測定値は偏差変数、即ち移動平均動作点からの偏差に変換される。異常事象検知用のＰＣＡモデルを生成する際に、平均動作点を除去する当該変換が必要である。これは、自身の未加工値から測定値の指数的にフィルタリングされた値（指数フィルタ式については式８、９を参照）を減算し、当該モデルでこの差を用いることにより行なわれる。
Ｘ’＝Ｘ−Ｘｆｉｌｔｅｒｅｄ    式１０
Ｘ’：動作点変化を除去すべく変換された測定値
Ｘ：元の未加工測定値
Ｘｆｉｌｔｅｒｅｄ：指標的にフィルタリングされた未加工測定値
【０１８０】
指数フィルタ用の時定数は、プロセスの主要な時定数とほぼ同じサイズでなければならない。多くの場合、時定数は約４０分で充分である。この変換の結果、ＰＣＡモデルへの入力が、移動平均動作点からのプロセスの最近の変化の測定値となる。
【０１８１】
この変換を正確に実行すべく、多くの場合毎分またはそれ以上速くオンライン・システムと合致するサンプリング周期でデータを集めなければならない。これは結果的に、１年分の動作データに対応すべく各々の測定について５２５，６００個のサンプルを集めることになる。この変換を計算したならば、データ集合を再サンプリングして、より管理しやすいサンプル数、通常は３０，０００〜５０，０００サンプルの範囲内に減らす。
【０１８２】
Ｖ．モデル構築
特定の測定値が選択されて、訓練データ集合が構築されたならば、標準的なツールを用いてモデルは素早く構築することができる。
【０１８３】
Ａ．モデル入力のスケーリング
ＰＣＡモデルの性能は、入力のスケーリングに依存している。スケーリングへの従来型のアプローチは、訓練データ集合内で各々の入力を自身の標準偏差σで分割することである。
Ｘｉ’＝Ｘｉ／σｉ    式１１
【０１８４】
多数のほぼ同一の測定値（固定された触媒反応床の多数の温度測定値）を含む入力集合の場合、このアプローチは更に、測定値をほぼ同一の測定値の個数の二乗根で除算することで修正される。
冗長なデータ群の場合
Ｘｉ’＝Ｘｉ／（σｉ＊ｓｑｒｔ（Ｎ））    式１２
ここで、Ｎ＝冗長なデータ群における入力の個数
【０１８５】
これらの従来方式のアプローチは、連続的な精製および化学プロセスからの測定値に対しては不適当である。プロセスは通常、指定された動作点でうまく制御されているため、データの分布は定常状態動作からのデータと、「乱された」および動作点が変化した動作からのデータとの組合せである。これらのデータは、優位な定常状態動作データからの過度に小さい標準偏差を有する。その結果生じるＰＣＡモデルは、プロセス測定値における微小又は適度な偏差に極端に影響される。
【０１８６】
連続的な精製および化学プロセスの場合、スケーリングは、正常プロセス撹乱の最中、または連続定常状態動作中に生じる変動性の程度に含まれない動作点変化の最中に生じる変動性の程度に基づくべきである。通常の無拘束変数の場合、スケーリング係数を決定する二つの異なる方式がある。
【０１８７】
第一の方式は、プロセスが定常状態で動作しておらず、且つ顕著な異常事象が生じていない期間を識別するものである。限られた個数の測定値が定常状態動作の主要な指標の役割を果たす。これらは典型的なプロセスの主要性能指標であって、通常はプロセスの供給速度、製品生産速度、および製品品質を含んでいる。これらの主要な尺度を用いて、動作を正常な定常状態動作、正常な撹乱された動作、および異常動作の期間に区分する。正常な撹乱された動作期間からの標準偏差が大部分の測定値に対する良好なスケーリング係数を与える。
【０１８８】
撹乱された動作に基づいてスケーリングを明示的に計算する代替的な方式は、以下のように訓練データ集合全体を用いるものである。スケーリング係数は、平均から３標準偏差外れたデータ分布に注目することにより近似することができる。例えば、データの９９．７％が平均の３標準偏差以内に存在し、データの９９．９９％が平均の４標準偏差以内に存在する筈である。平均から９９．７％と９９．９９％の間にあるデータ値の拡がりは、データ集合内の「撹乱された」データの標準偏差の近似としての役割を果たすことができる。図１２を参照されたい。
【０１８９】
最後に、測定値がしばしば拘束される（飽和変数に関する議論を参照）場合、スケーリング係数として用いる標準偏差の計算には変数が拘束されない期間だけを用いるべきである。
【０１９０】
Ｂ．主成分の個数を選択する
ＰＣＡは実際のプロセス変数を主成分ＰＣと呼ばれる一組の独立変数に変換する。これは、元の変数の線形結合である（式１３）。
ＰＣｉ＝Ａｉ，１＊Ｘ１＋Ａｉ，２＊Ｘ２＋Ａｉ，３＊Ｘ３＋．．．      式１３
【０１９１】
プロセスには、プロセスに影響を及ぼす特定の独立した事象を表わす多くの自由度がある。これらの異なる独立した事象は、プロセス変動とし当該プロセスのデータに出現する。プロセス変動は、供給速度の変更等の意図的な変更、または周囲温度の変化等の意図しない撹乱に起因する場合がある。
【０１９２】
各々の主成分は、プロセスに対するこれらの異なる独立した効果により生じるプロセスの変動性の一部をモデル化する。主成分は、データ集合における変動が減少する方向で抽出され、後続する各々の主成分は、次第にプロセス変動性をモデル化しなくなっていく。顕著な主成分はプロセス変動の顕著な発生源を表わす。例えば、供給速度の変更が最も大きいプロセス変動の発生源であるため、第１の主成分は通常、供給速度変更の効果を表わす。ある時点で構築担当者は、主成分によりモデル化されたプロセス変動が、いつ独立発生源を表わさなくなったかを決定しなければならない。
【０１９３】
主成分の正確な個数を選択する工学的アプローチは、主成分の一次寄与分である変数のグループがもはや工学的に意味を持たなくなった時点で止めることである。ＰＣによりモデル化されたプロセス変動の一次要因は、元の変数（それは、負荷と呼ばれる）の係数Ａｉ，ｎに着目することにより識別される。規模が比較的大きいそれらの係数は、特定のＰＣの主要な寄与分である。プロセスを良く理解する者であれば、ＰＣの主要な寄与分である変数の群に着目し、当該ＰＣに名前（例：供給速度効果）を割り当てることができる筈である。次第に多くのＰＣがデータから抽出されるにつれて、係数のサイズが更に等しくなる。この時点で、特定のＰＣによりモデル化されている変動は基本的にノイズである。
【０１９４】
いつＰＣが単にノイズをモデル化しているかを判定する従来の統計的方法は、各々の新規ＰＣによりモデル化されているプロセス変動がいつ一定になるかを識別することである。これはＰＲＥＳＳ統計で測定され、各々の連続するＰＣ（図１３）によりモデル化された変化の量をプロットする。不都合なことに、このテストは、精製および化学プロセス用に構築されたＰＣＡモデルには往々にして曖昧である。
【０１９５】
ＶＩ．モデルのテスト及びチューニング
プロセスデータは、ガウス即ち正規分布を示さない。その結果、残留誤りの３標準偏差における異常事象を検知するトリガを設定する標準的な統計的方法を用いてはならない。その代わり、トリガ位置は、モデルを用いて経験に基づいて経験的に設定されなければならない。
【０１９６】
最初に、拠点技術者が受容できる頻度、通常は毎日５、６回異常事象が通知されるようにトリガ・レベルを設定しなければならない。これは、訓練データ集合のＳＰＥＸ（これはＱ統計またはＤＭＯＤＸ統計とも呼ばれる）統計に着目して決定することができる。このレベルは、実際の異常事象を見逃すことなく、しかも誤警報に拠点技術者が振り回されないように設定されている。
【０１９７】
Ａ．モデルの強化
初期モデルが生成されたならば、新たな訓練データ集合を生成することにより強化する必要がある。これはモデルを用いてプロセスを監視することにより行なわれる。モデルが潜在的に異常な状況を示す場合、技術者はプロセス状況を調べて分類しなければならない。技術者は以下の３種の異なる状況を見出す。即ち何らかの特殊なプロセス動作が生じている、実際に異常な状況が生じている、或いは、プロセスは正常であって異常は虚偽の通知である。
【０１９８】
新たな訓練データ集合は、特殊動作および正常動作からのデータで構成されている。初期モデルを生成するために実行されたのと同じ分析をデータに対して実行氏、且つモデルを再計算することが必要である。この新たなモデルでは、依然としてトリガレバーが経験的に設定されるが、今回はより良い注釈が付けられたデータを伴うため、このトリガ点を調整して、真の異常事象が生じた場合のみ通知を与えるようにすることができる。
【０１９９】
異常事象検知のための単純な工学モデル
プロセス設備の物理的、化学的および機械的の設計、並びに多くの類似測定の挿入により、連続的な精製および化学プロセスからのデータにかなりの程度の冗長性をもたらす。この冗長性は、同一の測定が存在する場合は物理的冗長性と呼ばれ、物理的、化学的または機械的関係を用いてプロセス状態の独立であるが等価な推定を実行する場合は計算的冗長性と呼ばれる。このクラスのモデルを工学的冗長性モデルと呼ぶ。
【０２００】
Ｉ．２次元工学的冗長性モデル
これはモデルの最も簡単な形式であり、以下のような一般形式を有する。
Ｆ（ｙｉ）＝Ｇ（ｘｉ）＋フィルタリング済みバイアス（ｆｉｌｔｅｒｅｄ  ｂｉａｓ）ｉ＋オペレータ・バイアス（ｏｐｅｒａｔｏｒ  ｂｉａｓ）＋誤差ｉ    式１４
未加工バイアスｉ＝Ｆ（ｙｉ）−｛Ｇ（ｘｉ）＋フィルタリング済みバイアスｉ＋オペレータ・バイアス｝＝誤差ｉ    式１５
フィルタリング済みバイアスｉ＝フィルタリング済みバイアスｉ−１＋Ｎ＊未加工バイアスｉ−１    式１６
Ｎ−収束係数（例：０．０００１）
正常動作範囲：ｘｍｉｎ＜ｘ＜ｘｍａｘ
正常モデル偏差：−（最大誤差）＜誤差＜（最大誤差）
【０２０１】
オペレータが、モデルのシフトを必要とする何らかのフィールド事象（例：バイパスフロー）があったと判定する都度、「オペレータ・バイアス」項は更新される。オペレータの命令があれば、式１４が正確に満足される（誤差ｉ＝０）ように、オペレータ・バイアス項が更新される。
【０２０２】
工学的冗長性モデルにバイアスを与える持続的な未測定プロセス変動に対処すべく「フィルタリングされたバイアス」項は連続的に更新される。収束係数「Ｎ」は、ユーザーが指定した、通常は数日間にわたる期間の後で、あらゆる持続的な変化を除去すべく設定される。
【０２０３】
「正常動作範囲」および「正常モデル偏差」は、工学的冗長性モデルの履歴データから決定される。大多数の場合、最大誤差の値は単一の値であるが、これはまた、ｘ軸位置に依存する値のベクトルであり得る。
【０２０４】
物質収支、エネルギー収支、推定アナライザ表示対実際のアナライザ表示、圧縮機曲線等、任意の２次元方程式をこのように表わすことができる。図１４に、２次元のエネルギー収支を示す。
【０２０５】
適用例として、流量対バルブ位置モデルについてより詳細に述べる。
【０２０６】
Ａ．流量対バルブ位置モデル
特に有益な工学的冗長性モデルは、流量対バルブ位置モデルである。図２にこのモデルをグラフ的に示す。このモデルの特定の形式は以下の通りである。
【数３】
2008512800.tif 000005
ここで、
Ｆｌｏｗ：制御バルブを通る測定された流量
Ｄｅｌｔａ＿Ｐｒｅｓｓｕｒｅ＝測定された最も近い上流圧力−測定された最も近い下流圧力
Ｄｅｌｔａ＿Ｐｒｅｓｓｕｒｅｒｅｆｅｒｅｎｃｅ：正常動作の間の平均デルタ圧力
ａ：履歴データにフィッティングされたモデルパラメータ
Ｃｖ：履歴データから経験的に決定されたバルブ特性曲線
ＶＰ：制御バルブへの信号（実際の制御バルブ位置でない）
このモデルの目的は以下の通りである。
・詰まりつつある／詰まった制御バルブの検知
・凍結／故障した流量測定の検知
・制御システムが流量の制御を喪失した箇所の制御バルブ動作の検知
【０２０７】
この流量対バルブ式の特定の構成は、人的要因の理由で選択されている。この形式で方程式をｘ−ｙプロットするのがオペレータにとって最も容易に理解される。これらのモデルのいずれも、オペレータが最も理解しやすそうな方法で構成することが重要である。
【０２０８】
Ｂ．流量対バルブ位置モデルの構築
連続的な精製および化学プロセスが受ける長い期間にわたる定常状態動作ため、制御バルブの動作全体にわたる充分なデータを得るには長い履歴記録（１〜２年）が必要である。図１５に、定常動作が長期間にわたる流量、バルブ位置、およびデルタ圧力データの典型的な範囲を示す。最初の段階は、図に示すように、動作に何らかの重要な変動がある所で短い期間を切り離すことである。これは次いで、履歴内で様々な期間から取り出された正常動作の期間と混合されなければならない。
【０２０９】
往々にして、上流圧力（多くの場合ポンプ吐出）または下流圧力を利用できない。そのような場合、欠落している測定値が当該モデルにおける固定されたモデルパラメータになる。両方の圧力が欠落している場合、当該モデルに圧力効果を含めることは不可能である。
【０２１０】
バルブ特性曲線は、線形バルブ曲線、二次バルブ曲線、または区分的線形関数のいずれかにフィットする。区分的線形関数は、最も柔軟であって、任意の形状のバルブ特性曲線にもフィットする。
【０２１１】
バルブを直接横断して測定値が得られた場合、「ａ」の理論値は１／２である。測定値はそこには配置されない。「ａ」は、圧力測定値の実際の位置決めに対応すべく経験的に決定されたパラメータとなる。
【０２１２】
往々にして、デルタ圧力が変化する期間がほとんど無い場合もある。正常動作期間中のデルタ圧力のノイズがモデルフィッティング・プログラムを混乱させる恐れがある。これを克服するために、モデルを二相で構築する。即ち、最初は、デルタ圧力が変化する期間だけを含む小さいデータ集合を用いてモデルのフィッティングを行なう。次いで、決定された値で圧力従属パラメータ（「ａ」と、恐らくは欠落している上流または下流圧力）を固定して、大きい方のデータを用いてモデルを再構築する。
【０２１３】
Ｃ．流量対バルブ異常表示のファジー・ネット処理
任意の２次元工学的冗長性モデルと同様に、二種の異常性の尺度、即ち「正常動作範囲」および「正常モデル偏差」がある。「正常モデル偏差」は、正規化された指標即ち誤差／最大誤差に基づいている。これは、４型ファジー識別器（図１６）へ送られる。構築担当者は、正規化された指標を用いて、標準的な方法で正常（値ゼロ）から異常（値１）への遷移を検出することができる。
【０２１４】
「正常動作範囲」指標は、正常領域からのバルブ位置の距離である。これは通常、バルブ位置の変化の結果バルブ内の流量が殆どまたは全く変化しないバルブの動作の領域を表わす。構築担当者は再度、４型ファジー識別器を用いて、正常動作範囲の上下端、および正常動作から異常動作への遷移の両方をカバーすることができる。
【０２１５】
Ｄ．複数の流量／バルブモデルのグループ化
オペレータが好む流量／バルブモデルのグループ化の一般的な方法は、これらのモデルの全てを単一のファジー・ネットワークに入れて、傾向インジケータがオペレータに対し、それらの重要な流量コントローラの全てが機能していることを教える。その場合、ファジー・ネットワーク（図４）へのモデル表示は、各々の流量／バルブモデルに対して「正常動作範囲」および「正常モデル偏差」の表示を含んでいる。この傾向は、最悪モデル表示からの識別器の結果を含んでいよう。
【０２１６】
共通の設備タイプがグループ化された場合、このグループに注目するオペレータが好む別の方法として、流量／バルブのパレート図を用いるものがある（図１７）。このチャートにおいて、上位１０個の異常なバルブは、左端の最も異常なものから右端の最異常でないものまで動的に配置されている。各々のパレート・バーもまた、正常範囲内にあるモデル異常性表示の変化の程度を示す参照ボックスを備えている。図１７の図は、「バルブ１０」が実質的に正常ボックスの外側にあるが、他のものは全て正常に振舞っていることを示す。オペレータは次に、図２と同様に「バルブ１０」用のプロットを調べて、流量制御ループに関する問題を診断する。
【０２１７】
ＩＩ．多次元工学的冗長性モデル
次元が２より大きくなったならば、高次元工学的冗長性検査を扱うべく単一の「ＰＣＡ状」モデルが構築される。
【０２１８】
多次元冗長性の例として以下のものがある。
・圧力１＝圧力２＝．．．．＝圧力ｎ
・処理ユニット１へ入る物質流量＝処理ユニット１から出る物質流量＝．．．処理ユニット２へ入る物質流量
【０２１９】
測定値の較正誤差のため、これらの式は各々で係数を補償する必要がある。従って、最初に構築すべきモデル集合は以下の通りである。
Ｆ１（ｙｉ）＝ａ１Ｇ１（ｘｉ）＋フィルタリング済みバイアス１、ｉ＋オペレータ・バイアス１＋誤差１、ｉ
Ｆ２（ｙｉ）＝ａｎＧ２（ｘｉ）＋フィルタリング済みバイアス２、ｉ＋オペレータ・バイアス２＋誤差２、ｉ
Ｆｎ（ｙｉ）＝ａｎＧｎ（ｘｉ）＋フィルタリング済みバイアスｎ、ｉ＋オペレータ・バイアスｎ＋誤差ｎ、ｉ    式１８
【０２２０】
これらのモデルは、２次元の工学的冗長性モデルが構築されたのと同一の方法で構築される。
【０２２１】
この組の多次元検査はここで「ＰＣＡ状」モデルに変換される。この変換は、ＰＣＡモデルにおける主成分を、主成分係数（負荷）が、この独立した効果に起因する測定値の比例的変化を表わす、プロセスに対する独立した効果のモデルとしての解釈に依存する。図３において、３個の独立且つ冗長な処置、Ｘ１、Ｘ２、およびＸ３がある。Ｘ３が１ずつ変化する都度、Ｘ１はａ１だけ変化し、またＸ２はａ２だけ変化する。この関係の組は、係数がスケーリングされてない工学単位である、単一の主成分モデル（Ｐ）として表わされる。
Ｐ＝ａ１Ｘ１＋ａ２Ｘ２＋ａ３Ｘ３    式１９
ここにａ３＝１
【０２２２】
当該モデルのこの工学単位のバージョンは、以下のように標準的なＰＣＡモデル形式に変換することができる。
【０２２３】
標準的な統計概念との類似性を示すならば、各次元Ｘの変換係数は、正常動作範囲に基づいていてよい。例えば、平均の周囲の３σを用いて正常動作範囲を定義すれば、スケーリングされた変数は以下のように定義される。
Ｘｓｃａｌｅ＝Ｘｎｏｒｍａｌ  ｏｐｅｒａｔｉｎｇ  ｒａｎｇｅ／６σ    式２０
（正常動作データの９９．７％は平均から３σ以内に入る筈である）
Ｘｍｉｄ＝Ｘｍｉｄ  ｐｏｉｎｔ  ｏｆ  ｏｐｅｒａｔｉｎｇ  ｒａｎｇｅ    式２１
（「平均」を正常動作範囲の中点として明示的に定義）
Ｘ’＝（Ｘ−Ｘｍｉｄ）／Ｘｓｃａｌｅ    式２２
（平均とσが決定された際の標準的なＰＣＡスケーリング）
Ｘｉ用のＰ’負荷は以下の通りである。
【数４】
2008512800.tif 000006
（負荷ベクトルを正規化する旨の要件）
これはＰを以下のように変換する
Ｐ’＝ｂ１＊Ｘ１＋ｂ２＊Ｘ２＋・・・＋ｂｎ＊ＸＮ    式２４
Ｐ’「標準偏差」＝ｂ１＋ｂ２＋・・・＋ｂｎ    式２５
【０２２４】
この変換により、多次元工学的冗長性モデルはここで、標準的なＰＣＡ構造を計算、例外処理、オペレータへの提示および対話に用いて取り扱うことができる。
【０２２５】
異常事象検知のためにＰＣＡモデルおよび簡単な工学モデルの配備
Ｉ．オペレータおよび既知の事象抑制
抑制論理は、以下の目的に必要とされる。
・測定可能な異常事象から虚偽の表示を除外する方法を提供する
・オペレータが調べた異常な表示を消去する方法を提供する
・保守のため一時的にモデルまたは測定値を無効にする方法を提供する
・不良動作するモデルを、それらが再調整されるまで使用不可にする方法を提供する
・不良動作する計器を永久に使用不可にする方法を提供する
【０２２６】
抑制には２種類ある。即ち外部の測定可能な事象により自動的に起動される抑制と、オペレータにより起動される抑制である。これらの２種類の抑制の背後にある論理を図１８、１９に示す。これらの図はファジー化モデル指標に生じている抑制を示しているが、抑制は、特定の測定値に対して、特定のモデル指標に対して、モデル全体に対して、或いはプロセス領域内のモデルの組合せに対しても生じ得る。
【０２２７】
オペレータが起動した抑制には、抑制が終了した時点を判定する２個のタイマーがある。１個のタイマーは、抑制された情報が正常状態を戻ってそのままであることを確かめる。このタイマーの典型的な値は１５〜３０分である。第２のタイマーは、正常状態に戻ったか否かに拘らず、異常事象の検査を再起動させる。このタイマーの典型的な値は、オペレータの勤務シフト（８〜１２時間）に等しいか、または半永久的な抑制の場合には極めて長時間となる。
【０２２８】
事象に基づく抑制の場合、測定可能なトリガが必要である。これは、オペレータによるセットポイントの変更、突然の測定値変化、またはデジタル信号であり得る。この信号は、図２０に示すタイミング信号に変換される。このタイミング信号は、以下の式を用いてトリガ信号から生成される。
Ｙｎ＝Ｐ＊Ｙｎ−１＋（１−Ｐ）＊Ｘｎ    指数フィルタの式    式２６
Ｐ＝Ｅｘｐ（−Ｔｓ／Ｔｆ）    フィルタ定数の計算    式２７
Ｚｎ＝Ｘｎ−Ｙｎ    タイミング信号の計算  式２８
ここで、
Ｙｎ  トリガ信号の現在のフィルタリング済み値
Ｙｎ−１  トリガ信号の以前のフィルタリング済み値
Ｘｎ  トリガ信号の現在の値
Ｚｎ  図２０に示すタイミング信号
Ｐ  指数フィルタ定数
Ｔｓ  測定値のサンプル時間
Ｔｆ  フィルタ時定数
【０２２９】
タイミング信号が閾値（図２０に０．０５として示す）を超える限り、事象は抑制されたままである。構築担当者は、フィルタ時定数Ｔｆを変えることにより抑制の長さを設定する。この機能のために簡単なタイマーを用いてもよいが、このタイミング信号は異なるサイズのトリガ信号に対応しており、変化が大きい場合には抑制をより長く、変化が小さい場合には抑制をより短く生成する。
【０２３０】
図２１に事象抑制とオペレータ抑制がＰＣＡモデル内の所定の入力組を無効にする様子を示す。自動的に抑制される入力組がオンラインモデル性能から決定される。オペレータが見たくない表示をＰＣＡモデルが出す場合は常に、この表示を誤差の２乗和指標に対する少数の個々の寄与分まで辿ることができる。これらの個々の寄与分を抑制すべく、この指標の計算は以下のように修正される。
【数５】
2008512800.tif 000007
ｗｉ：入力ｉ（通常は１に等しい）に対する寄与分の重み
ｅｉ：入力ｉからの２乗和誤差の寄与分
【０２３１】
トリガ事象が生じた場合、抑制対象である入力の各々について寄与分の重みがゼロに設定される。これらの入力が再有効化される場合、寄与分の重みは段階的に値１に戻される。
【０２３２】
ＩＩ．ＰＣＡモデル分解
ＰＣＡモデルは広範なプロセス設備範囲を用いて構築されているが、モデル指標を分離して、オペレータのプロセスに対する認識によりよく合致して、異常事象の兆候に対する感度を向上させるグループに分けることができる。
【０２３３】
再び式２９を参照すると、誤差の２乗和のいくつかのグループ分けが可能である。
【数６】
2008512800.tif 000008
【０２３４】
通常、これらのグループ分けは設備（例：装置の再沸器区画）のより小さいサブユニット周辺、または設備の機能に関連する（例：製品の品質）サブグループ分けに基づいている。
【０２３５】
各々の寄与分ｅｉがプロセスノイズに基づいて常に誤差の２乗和に追加されるため、ノイズに起因する指標のサイズは当該指標に寄与する入力の個数に線形に増加する。誤差の２乗和計算への寄与がより小さいため、当該指標の信号対ノイズ比が向上して、異常事象に対する指標の応答性が向上する。
【０２３６】
同様の方法で、各々の主成分を設備のグループ分けに一致するように再分割することができ、各々のサブグループ用にホテリングＴ２指標に類似した指標を生成することができる。
【数７】
2008512800.tif 000009
【０２３７】
これらの指標の閾値は、テストデータをモデルに通して、テストデータに対するそれらの性能に基づいて閾値の感度を設定することにより計算される。
【０２３８】
これらの新たな指標は、通常のＰＣＡモデルが扱われるのと同一の方法でオペレータ向けに解釈される。元の入力に基づくパレート図を、誤差の２乗和の指標への最大寄与分、およびＴ２計算における最大Ｐへの最大寄与分について示す。
【０２３９】
ＩＩＩ．重なり合うＰＣＡモデル
入力がいくつかのＰＣＡモデルに現われることにより、モデルに影響を及ぼす全ての相互作用がモデルの範囲内で含まれるようになる。これにより、これらの入力が誤差の２乗和指標の主要な寄与分である場合、オペレータ向けに複数の表示を与えることがあり得る。
【０２４０】
この問題を回避するために、複数のＰＣＡモデルに現われるあらゆる入力の対し、それらのＰＣＡモデルの１個をその一次モデルとして割り当てる。一次ＰＣＡモデル向けの式２９における寄与分の重みは１に保持される一方、非一次ＰＣＡモデルの重みはゼロに設定される。
【０２４１】
ＩＶ．オペレータ対話環境およびインターフェース設計
オペレータ・インターフェースの一次目的は以下の通りである。
・オペレータの権限下にある主要プロセス領域が正常である旨の表示を連続的に提供する
・背景にあるモデル情報への迅速な（マウス・クリック１、２回）ナビゲーションを提供する
・オペレータに対し、モデルを使用可能にする制御を提供する
図２２は、オペレータが用いる一次インターフェースにおいてこれらの設計目的がどのように表わされるかを示す。
【０２４２】
ファジー・ペトリネットからの最終出力は、図４に示すように正常性傾向である。この傾向は、ファジー識別関数に定義された異常の最大尤度を示すモデル指標を表わす。要約に示される傾向の個数には自由度があり、オペレータとの議論で決定される。この傾向の上に、オペレータが処置を講じるべき時にヘルプ信号への２本の基準線があり、黄色の線は通常は０．６の値で設定され、赤色の線は通常は０．９の値で設定される。これらの線は、オペレータに対し、いつ処置を講じるべきであるかの指針を与える。傾向が黄色の線を越えた場合、図４における緑色の三角形が黄色に変わり、傾向が赤色の線を越えた場合、緑色の三角形は赤色になる。当該三角形にはまた、最も異常な表示を与えているモデルに関連付けられたディスプレイにオペレータを案内する機能を備えている。
【０２４３】
当該モデルがＰＣＡモデルであるか、または設備グループ（例：全ての制御バルブ）の一部である場合、緑色の三角形を選択することによりパレート図が作成される。ＰＣＡモデルの場合、これはモデル指標への多数の最大寄与分のうち、これは最も異常（左側）から最も異常でない（右側）までを示す。大抵は、主要な異常事象インジケータは最初の２、３個の測定値の中にある。パレート図は、異常の兆候であると見なされる前にどの程度測定が通常とは異なるかをオペレータに示すべく各バーの周辺に赤い囲み枠を含んでいる。
【０２４４】
ＰＣＡモデルの場合、オペレータには、棒グラフパレート内における順序に合致する傾向パレートが与えられる。傾向パレートにより、各プロットは２個の傾向、即ち実際の測定（シアンで）および、すべてが正常（褐色）な場合に測定値が取るべきであった値のＰＣＡモデルからの推定を有する。
【０２４５】
バルブ／流量モデルの場合、パレートの下での詳細は、２次元の流量対バルブ位置モデル・プロットである。このプロットから、オペレータはモデルに対しオペレータ・バイアスを適用することができる。
【０２４６】
設備がグループ分けされていない場合、緑色の三角形を選択することにより、オペレータは概略傾向の下での最悪２次元モデルへ直接誘導する。
【０２４７】
各々の棒の下でオン／オフ動作のボタンを選択することにより、オペレータによる抑制をパレート図レベルで行なうことができる。
【０２４８】
【表１】
2008512800.tif 000010
【０２４９】
【表２】
2008512800.tif 000011
【０２５０】
【表３】
2008512800.tif 000012
【０２５１】
【表４】
2008512800.tif 000013
【０２５２】
【表５】
2008512800.tif 000014
【０２５３】
【表６】
2008512800.tif 000015
【０２５４】
【表７】
2008512800.tif 000016
【０２５５】
【表８】
2008512800.tif 000017
【０２５６】
付録３
工学モデル／推論
Ａ．再生装置スタックバルブ・モニター
再生装置スタックバルブＡ、Ｂの値を差圧コントローラ出力に対して相互検証する。正常な状態の下では、それらは全て合致しなければならない。
【０２５７】
【表９】
2008512800.tif 000018
【０２５８】
【表１０】
2008512800.tif 000019
【０２５９】
Ｄ．炭素収支：
本モニターは、触媒循環ＣＣＲ−ＰＣＡモデルの第４主成分のＴ統計量に的を絞る。
【０２６０】
Ｅ．主精留装置への触媒キャリーオーバー
本モニターは、以下の変数が限度内にあるか否かを調べる
（ａ）反応炉ストリッパ・レベル
（ｂ）反応炉の差圧
（ｃ）主精留装置底部ストレーナー差圧、および
（ｄ）主精留装置からのスラリーポンプ動作圧
【０２６１】
【表１１】
2008512800.tif 000020
【０２６２】
Ｇ．バルブ−流量モデル
ＡＥＤアプリケーション用に構築された合計１２個のバルブモデルがある。全てのバルブモデルにはバイアス更新が実装されている。流量は、以下のようにデルタ圧を補償される。
補償流量＝ＦＬ／（ＤＰ／ＳｔｄＤＰ）＾ａ、ここに、
ＦＬ＝実際の流量、ＤＰ＝上流圧−下流圧、ＳｔｄＤＰ＝標準デルタ圧、ａはパラメータである。次いで、モデルの整合性を検証すべく推定補償流量と実際の補償流量の間のプロット図を作成する（Ｘ−Ｙプロット）。以下は、１２個のバルブ流量モデルのリストである。以下のモデルの変数の順序は（ＯＰ、ＦＬ、ＵｐＰ−ＤｎＰ、ＳｔｄＤＰ、ａ、Ｂｏｕｎｄ）である。
【０２６３】
【表１２】
2008512800.tif 000021
【図面の簡単な説明】
【０２６４】
【図１】オンライン・システムの情報が、各種の変換、モデル計算、ファジー・ペトリネット、および圧密処理を経由して流れ、プロセス領域の正常性／異常性を示す概略傾向に到達する様子を示す。
【図２】オペレータにバルブ流量のプロットを簡単なｘ−ｙプロット図として示す。
【図３】ＰＣＡモデルとして表わされた３次元冗長性を示す。
【図４】ファジー・ネットワーク機構の模式図を示す。
【図５】異常事象アプリケーション構築の全般的プロセスの模式図を示す。
【図６】プロセス制御カスケードの内部構造の模式図を示す。
【図７】多変数拘束コントローラＭＶＣＣの解剖の模式図を示す。
【図８】現在の品質のオンライン推論的評価の模式図を示す。
【図９】履歴データのＫＰＩ解析を示す。
【図１０】信号対ノイズ比を示す。
【図１１】プロセスダイナミクスが２種の測定値の現在値同士の相互相関を乱し得る様子を示す。
【図１２】プロセスデータの確率分布を示す。
【図１３】押圧統計量を示す。
【図１４】２次元エネルギー収支モデルを示す。
【図１５】長期間定常動作した場合の流量、バルブ位置、およびデルタ圧力データの典型的な動きを示す。
【図１６】第４種のファジー識別器を示す。
【図１７】流量対バルブのパレート図を示す。
【図１８】オペレータ抑制ロジックの模式図を示す。
【図１９】事象抑制ロジックの模式図を示す。
【図２０】事象抑制の持続期間の設定を示す。
【図２１】ＰＣＡモデルにおいて所定の組の入力を無効にする事象抑制およびオペレータ抑制を示す。
【図２２】オペレータが用いる一次インターフェースで設計目的がどのように表わされるかを示す。
【図２３】ＦＣＣＵの模式的レイアウトを示す。
【図２４】ＦＣＣＵおよびライトエンド塔の全般的模式図を示す。
【図２５】ＦＣＣＵ動作における全ての問題含みのモニターのオペレータ向け表示を示す。
【図２６】触媒循環問題に対するファジー論理に基づく連続異常性インジケータを示す。
【図２７】触媒循環問題の完全な掘り下げを、全ての肯定的な証拠と共に示す。
【図２８】触媒循環問題のファジー論理ネットワークを示す。
【図２９】触媒循環、ＦＣＣ例外およびＦＣＣ極限異常事象モニターにおける警報を示す。
【図３０】図２９のＦＣＣ例外シナリオに関係するタグのパレート図を示す。
【図３１】図３０のタグのマルチ傾向を示す。これは、タグ値およびモデル予測も示す。
【図３２】逸脱するバルブ流量モデルのランク付けされたリスト（パレート図）を示す。
【図３３】バルブ流量モデル、バルブ開口対流量のＸ−Ｙプロット図を示す。
【図３４】送風機モニターにおけるパレート図およびＸ−Ｙプロット線を示す。
【図３５】再生装置スタックバルブ・モニターにおける掘り下げを示す。
【図３６】再生装置サイクロン・モニターにおける掘り下げを示す。
【図３７】送風機モニターにおける掘り下げを示す。
【図３８】炭素収支モニターにおける掘り下げを示す。
【図３９】主精留装置への触媒キャリーオーバーにおける掘り下げを示す。
【図４０】湿式ガス圧縮機における掘り下げを示す。
【図４１】バルブ流量モニターのファジー・ネットワークを示す。
【図４２】制御可能な範囲から外れたバルブの例を示す。
【図４３】事象抑制ディスプレイを示す。
【図４４】ＡＥＤ事象供給バック形式を示す。
【図４５】各々の連続的なＰＣによりモデル化された変動量をプロットする標準的な統計プログラムを示す。
【図１】
2008512800.tif 000022
【図２】
2008512800.tif 000023
【図３】
2008512800.tif 000024
【図４】
2008512800.tif 000025
【図５】
2008512800.tif 000026
【図６】
2008512800.tif 000027
【図７】
2008512800.tif 000028
【図８】
2008512800.tif 000029
【図９】
2008512800.tif 000030
【図１０】
2008512800.tif 000031
【図１１】
2008512800.tif 000032
【図１２】
2008512800.tif 000033
【図１３】
2008512800.tif 000034
【図１４】
2008512800.tif 000035
【図１５】
2008512800.tif 000036
【図１６】
2008512800.tif 000037
【図１７】
2008512800.tif 000038
【図１８】
2008512800.tif 000039
【図１９】
2008512800.tif 000040
【図２０】
2008512800.tif 000041
【図２１】
2008512800.tif 000042
【図２２】
2008512800.tif 000043
【図２３】
2008512800.tif 000044
【図２４】
2008512800.tif 000045
【図２５】
2008512800.tif 000046
【図２６】
2008512800.tif 000047
【図２７】
2008512800.tif 000048
【図２８】
2008512800.tif 000049
【図２９】
2008512800.tif 000050
【図３０】
2008512800.tif 000051
【図３１】
2008512800.tif 000052
【図３２】
2008512800.tif 000053
【図３３】
2008512800.tif 000054
【図３４】
2008512800.tif 000055
【図３５】
2008512800.tif 000056
【図３６】
2008512800.tif 000057
【図３７】
2008512800.tif 000058
【図３８】
2008512800.tif 000059
【図３９】
2008512800.tif 000060
【図４０】
2008512800.tif 000061
【図４１】
2008512800.tif 000062
【図４２】
2008512800.tif 000063
【図４３】
2008512800.tif 000064
【図４４】
2008512800.tif 000065
【図４５】
2008512800.tif 000066
【国際調査報告】
2008512800.tif 000067