(57)【要約】  本発明は、異常なオンライン解析またはラボラトリ計測を検出、および異常プロセス状態に起因する異常品質エクスカーションを検出するための方法およびシステムである。【選択図】図２
2010519636.tif 000002
【特許請求の範囲】
【請求項１】
オンラインアナライザまたはラボラトリ解析によって計測される運転パラメータの異常検出のための推論推定を計算するオンラインシステムであって、
ａ）リアルタイムデータコレクタ、
ｂ）リアルタイムデータプリプロセッサ、
ｃ）推論モデルを含む、モデルおよび統計的有意性計算器、
ｄ）複数のモデル結果を要約する手段、および
ｅ）事象解析用の診断ディスプレイ
を含むことを特徴とするシステム。
【請求項２】
前記データコレクタが、前記モデルへの入力値としてのプロセス計測値を計測することを特徴とする請求項１に記載のシステム。
【請求項３】
前記プロセス計測値に温度、圧力および流量が含まれることを特徴とする請求項２に記載のシステム。
【請求項４】
前記データプリプロセッサが、データの全体エラー検出を含むことを特徴とする請求項１に記載のシステム。
【請求項５】
前記データプリプロセッサが、データの数学的変換を含むことを特徴とする請求項１に記載のシステム。
【請求項６】
前記データプリプロセッサが、ノイズを除去するためのフィルタを含むことを特徴とする請求項１に記載のシステム。
【請求項７】
前記モデルおよび統計的有意性計算器が、現在の運転パラメータ（Ｙｅｓｔｉｍａｔｅ）、将来の予測運転パラメータ（Ｙｐｒｅｄｉｃｔｅｄ）並びに前記モデルおよび／または生の計測値における長期ドリフト問題を判定することを特徴とする請求項１に記載のシステム。
【請求項８】
前記モデルおよび統計的有意性計算器が、前記運転パラメータの値の推定値を求めることを特徴とする請求項１に記載のシステム。
【請求項９】
前記入力値の全てが現在のレベルのまま維持される場合に、前記モデルが前記運転パラメータの最終的な定常状態値の予測将来値（Ｙｐｒｅｄｉｃｔ）を決定することを特徴とする請求項２に記載のシステム。
【請求項１０】
前記モデルおよび統計的有意性計算器が、プロセス制御アプリケーション用に将来の運転パラメータ（Ｙｃｏｎｔｒｏｌ）を更に決定することを特徴とする請求項７に記載のシステム。
【請求項１１】
前記Ｙｅｓｔｉｍａｔｅのためのモデルへの入力値が、対応するラボラトリ解析またはオンライン解析と時間同期されることを特徴とする請求項７に記載のシステム。
【請求項１２】
前記推論モデルが、カルマンフィルタ、部分最小二乗（ＰＬＳ）回帰、ニューラルネットおよび段階的回帰を使用して構築可能であることを特徴とする請求項１に記載のシステム。
【請求項１３】
前記推論モデルが、部分最小二乗（ＰＬＳ）モデルであることを特徴とする請求項１１に記載のシステム。
【請求項１４】
前記モデルがバイアス項を含むことを特徴とする請求項１１に記載のシステム。
【請求項１５】
運転パラメータのオンライン解析またはラボラトリ解析が提供されるたびに、前記バイアス項が更新されることを特徴とする請求項１４に記載のシステム。
【請求項１６】
前記統計的有意性計算器が、前記計測された運転パラメータと前記運転パラメータの推定値の間に統計的に有意な差が存在するかどうかを判定することを特徴とする請求項８に記載のシステム。
【請求項１７】
前記統計的有意性計算器が、前記計測された運転パラメータと前記運転パラメータのモデル推定値の差に基づいて更新されるバイアス項を含むことを特徴とする請求項１６に記載のシステム。
【請求項１８】
前記モデルが、前記運転パラメータの予測値と前記運転限界を比較することにより、将来の運転限界エクスカーションを決定することを特徴とする請求項９に記載のシステム。
【請求項１９】
前記推論モデルが、ラボラトリ解析のためのサロゲート計測値の識別を含むことを特徴とする請求項１に記載のシステム。
【請求項２０】
前記推論モデルが、定常運転の期間の識別と、大部分の定常状態運転の削除を含むことを特徴とする請求項１に記載のシステム。
【請求項２１】
前記推論モデルが、入力および出力計測のための正常範囲の識別を含むことを特徴とする請求項１に記載のシステム。
【請求項２２】
前記推論モデルが、極めて異常な運転期間中のデータを除外するためのデータフィルタを含むことを特徴とする請求項１に記載のシステム。
【請求項２３】
前記推論モデルが、単入力単出力モデルを含むことを特徴とする請求項１に記載のシステム。
【請求項２４】
前記推論モデルが、前記データからの低速トレンドの除去を含むことを特徴とする請求項１に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概して連続工業プロセスにおける異常事象の早期検出に関し、より詳しくは、主要な運転パラメータの不正な計測値の検出および主要な運転パラメータの将来的な異常エクスカーションの予測に関する。
【背景技術】
【０００２】
連続工業プロセスが経済的にほぼ最適な条件で運転されている場合、製品品質仕様等の主要な運転パラメータは最大または最小限界で運転されている。従って、連続工業プロセス（製油所や化学プラント等）の効率的な運転および異常事象の防止のどちらにとっても、これらのパラメータの現在値および予測将来値を把握することが非常に重要である。例えば、異常な品質エクスカーションは、仕様範囲からの製品の逸脱、プロセス装置の突然の誤動作（蒸気形成によるポンプのキャビテーション等）、およびプロセス性能の劣化（触媒上へのコークスの蓄積により反応が失われたり煙管内でのコークス生成により熱伝達が失われたり等）を引き起こす恐れがある。
【０００３】
プロセスストリームの品質や他の主要な運転パラメータを直接計測すると、多大な費用がかかり、且つ問題が発生しやすくなる恐れがある。オンライン解析の場合は、高い初期設置コストと高い保守コストの両方がかかる。オンライン解析では、多くの場合、専用のプロセスサンプリングシステムと、解析装置用に環境から保護されたフィールドシェルタとを必要とする。この装置の保守には、特別に訓練された要員と高度な予防的保守の取り組みとが必要な場合がある。しかし、保守はほとんどの場合、把握されているオンラインアナライザ問題に対して行われるだけである。最近のオンラインアナライザシステムには、テストおよびキャリブレーションのための標準的なサンプルと、連続的な装置診断を実行するマイクロコンピュータとが組み込まれている。
【０００４】
多くの場合、現場ではオンライン解析とラボラトリ解析を併用して品質計測を実施することを選択するか、オンライン解析ではなくラボラトリ解析を利用して品質計測を実施することを選択する。現場サンプルの採取およびそれらのサンプルの解析には広範囲にわたり人為的介入が必要であるため、このラボラトリ解析は、通常、頻繁的でなく（毎日１回〜毎週１回の頻度で）、正常変動性が大きく、エラーレートが高い。
【０００５】
オンライン解析およびラボラトリ解析の手法を補足するために、より簡単に入手可能なプロセス計測値（主に、温度、圧力、流量）を基に品質パラメータの推論推定値を作成することができる。推論計測値には従来の用途として２つの用途があり、１つは、より低速でサンプリングされ、閉ループプロセス制御アプリケーション内で使用されるアナライザ値の連続的な推定値を作成することであり、もう１つはアナライザ値およびラボラトリ値を検証することである。これらの用途については、実際のオンラインアナライザ値またはラボラトリ計測値をモデルに迅速に更新することにより、十分な機能を持たないモデルであってもそれなりに十分な性能を実現することが可能である。モデルが、次のアナライザサンプルを推定するためにある程度の能力を備えていれば、モデルの動作は、最後のアナライザサンプルを次のアナライザサンプルの推定値として使用する場合と比較してほとんど変わらない。しかし、新しいアナライザサンプル値が最小変化限度および最大変化限度内に確実に収まるようにする場合を除き、迅速な更新を使用するモデルは、異常なアナライザサンプル値を検出する場合または異常プロセス事象による異常な品質エクスカーションを予測する場合には適切でない。
【０００６】
これらの用途については、モデル内の欠陥と実際の異常事象を区別するような問題は発生し得ない。このため、モデルを構築するにあたっては、最高品質のトレーニングデータだけしか使用できない。
【０００７】
連続プロセス業界における推論計測の大部分は、ニューラルネット、段階的回帰、部分最小二乗（ＰＬＳ）等のプロセスデータ駆動方法（モデル構造とモデルパラメータは運転データから決定される）を使用して開発されている。或いは、運転データからモデルパラメータのみが決定される第１原理エンジニアリングモデルに基づいて推論計測を行うこともできる。これらの手法を使用して開発されたモデルの品質は、モデルを構築するために選択されたデータの品質、または第１原理エンジニアリングモデル内のパラメータにフィッティングさせるために選択されたデータの品質に大きく影響される。データ選択の方法、データ解析の方法、およびデータの条件付けの方法は、簡単な統計的仮定に基づく一般的な手法に基づくのではなく、データの特性に合わせて調整される必要がある。連続プロセス産業で高品質な推論推定を開発することができなかった場合に、プロセスデータの特性に適合しない無効なデータ選択方法やデータ条件付けの方法まで原因をたどることが可能な場合が多い。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
異常検出用の方法またはシステム（ＡＥＤ）（例えば、米国特許出願公開第２００６／００５８８９８号明細書を参照）の目的は、プロセスおよび装置の問題が深刻な出来事に拡大するのを防止することにある。本方法またはシステムでは、最初に、プロセスまたは装置の現れつつある問題の早期警告を、警報システムが始動される前にオペレータに提供することによって、次に、問題の根本的原因を局所化して診断するための主要な情報を、オペレータに提供することによって目的を達成する。
【課題を解決するための手段】
【０００９】
本発明では、主要な運転パラメータ値（特にオンラインアナライザまたはラボラトリ解析からの）の異常値を検出、および最近の異常プロセス状態に起因する、製品品質等の主要プロセスパラメータの異常エクスカーションを予測するための方法およびシステムを、異常事象検出を使用して説明する。要約すると、本発明は、時間動特性を明らかにする入力値の前処理とＰＬＳ等の従来技術アルゴリズムとを組み合わせた、運転パラメータ尺度の異常事象を検出するための推論モデルと、運転パラメータの高品質推論モデルを構築して異常値の検出および異常エクスカーションの予測を行うための方法（図１）と、モデルをオンラインで実施するためのシステム（図２）と、を含む。
【００１０】
オンラインで実施するためのシステムは、以下のステップを含む。
１．リアルタイムデータを前処理／時間同期する
２．現在の品質の推定値、将来の予測品質値、およびプロセス制御アプリケーションで使用する品質の推定値を算出する
３．実際の計測値に基づいてモデルをオンラインで更新する
４．モデルの推定結果を解釈する
５．複数の運転パラメータの計測値の正常／異常ステータスを集約および要約する
６．異常事象診断情報をコンソールのオペレータに提供する
７．モデル診断情報を保守技術者に提供する
【００１１】
本発明の中心は、オンラインアナライザによってまたはラボラトリ解析によって計測される主要な運転パラメータの値を推定する推論モデルである。このモデルでは、温度、圧力、流量等の簡単に入手可能なプロセス計測値を使用する。ニューラルネットワーク（ＮＮ）、部分最小二乗（ＰＬＳ）、線形回帰等の推論計測モデルを計算するための従来技術アルゴリズムが多数存在する。しかし、これらの方法では、連続工業プロセスの特性である時間動特性を組み込むのが難しい。
【００１２】
推論モデルについて本発明では、従来技術の方法、特にＰＬＳ（しかし、これに限定されるものではない）を、モデル入力データの３種類の前処理と組み合わせることにより、連続工業プロセスに固有の時間動特性を処理する（時間同期として知られているステップ）。時間同期の各形態は、モデルが適用される特定の使用に合わせて構成されている。これにより、２つの部分からなるモデル計算がもたらされる。この計算では、最初に時間同期された各入力が計算され、次にこれらの時間同期された入力を集約して出力値の推定値とする。入力は、アナライザ／ラボラトリ計測の現在の値を推定するため、アナライザ計測の将来の値を予測するため、モデルベースの予測制御アルゴリズムへの入力として使用されるため、または比例・積分・微分（ＰＩＤ）アルゴリズム等の標準的な制御アルゴリズムの入力として使用されるために、種々の方法で時間同期される。
【００１３】
推論モデルをトレーニングするために、本発明は、連続工業プロセスによって生成されるデータの特性に対して構成されるモデルトレーニングデータセットを作成するモデル開発手法を含む。この手法は、以下のステップを含む。
・設計された実験データではなく正常な運転データを使用する。
・少数の運転点で長期化する定常状態運転を明らかにする。
・プロセス時間動特性を明らかにする。
・モデル入力間の相互相関を明らかにする。
・解析／運転パラメータへの計測されない影響を明らかにする。
【００１４】
オンラインシステムでは、データの前処理を行い、オフラインモデル開発と一貫性のある方法で推論モデルを計算する。更に、オンラインシステムは、モデルをリアルタイムで適応させる方法と、異常事象が存在するかどうかについて推論モデル計算およびアナライザ／ラボラトリ計測を解釈する方法とを含む。オペレータは異常事象の指標を受信すると、本システムを利用して基礎となる異常指標および支援データにドリルダウンすることができる。オペレータは、異常指標および支援データを自身で分析して、修正処置を決定する必要がある。
【００１５】
典型的な異常事象検出アプリケーションの場合は、個々のコンソールオペレータの処理範囲内で、１０〜２０の主要なプロセスパラメータ／品質値を有する。ファジーペトリネットを使用して、これらの主要なプロセスパラメータモデルを、以下のいずれかの方法で処理する。
・マージして、主要な全てのプロセスパラメータの正常／異常ステータスに関する単一の要約トレンドを提供する。
・同じプロセスサブエリア内の他のモデルとマージして、そのプロセスサブエリアの正常／異常ステータスに関する要約トレンドを提供する。
【００１６】
このようにして、オンラインシステムは、要約されたプロセス条件の正常／異常ステータスをプロセスオペレータに提供することができる。
【図面の簡単な説明】
【００１７】
【図１】推論モデルを開発するための方法を示す。本発明の方法部分の主要なステップを示す。実線のボックスは従来技術のステップを示し、点線のボックスは本発明でカバーされる領域を示す。
【図２】異常事象検出のために推論推定値を計算するオンラインシステムを示す。この図は本発明のオンラインシステム部分についてデータおよび計算の流れを示す。
【図３】異常事象の可能性への統計的有意性のファジー変換を示す。この図は、統計的有意性の概念を使用して、ファジーネットによる異常事象の可能性の項への全てのモデルの残差を正規化する方法を示す。
【図４】解析異常性の可能性を要約するためのファジーネットを示す。この図は、あらゆる解析の中での異常性の複合的な可能性を個々の可能性から算出する方法を示す。
【図５】Ｙｐｒｅｄｉｃｔｅｄが品質上限を超える９０％可能性の例を示す。この図では、Ｙｐｒｅｄｉｃｔｅｄのデータ分布を使用して、異常解析エクスカーションを予測する際にファジーネットワークにおける限度を決定する方法を示す。
【図６】品質限度の超過の可能性を要約するファジーネットを示す。この図は、異常な限度超過の複数の予測をマージして、異常な限界超過の可能性を要約する単一の共通トレンド表示とする方法を示す。
【図７】オペレータ用ディスプレイの視覚化コンポーネントを示す。この図は、潜在的な異常事象をオペレータが診断するのを支援する、オンラインシステムの主要なオペレータ用ディスプレイを示す。
【図８】推論モデルを開発するための典型的な従来技術の方法を示す。この図は、推論モデルを開発するための典型的な従来技術の方法における主要なステップを示す。
【図９】トレーニングデータから定常状態運転を除外しない場合の結果の例を示す。この図では、定常状態運転からデータを除外しないと、推論モデルによる推定結果にどのような影響を及ぼすかを例示する。
【図１０】ステップの変更に対する一次むだ時間プロセスの動的応答を示す。この図では、連続工業プロセスで生じる時間動特性に典型的なものである動的応答を示す。
【図１１】一組の蒸留塔を示す。
【図１２】アナライザへのサンプルラインの部分的なプラギングと関連付けられた設定時間トレンドを示す。た設定時間トレンドを示す。
【図１３】アナライザのドリフティングキャリブレーションと関連付けられた設定時間トレンドを示す。
【発明を実施するための形態】
【００１８】
本発明では、以下の処理を実行する推論モデルを使用して製品品質等の主要な運転パラメータでの異常事象を検出するシステムを提供する。
・主要な運転パラメータの現在値を推定する
・主要な運転パラメータの無効な計測値を早期に検出する
・主要な運転パラメータの将来的な異常エクスカーションを予測する
【００１９】
更に、本発明では、オンラインシステム用に必要な推論モデルおよび解析メカニズムを作成する方法についても提供する。
【００２０】
Ｉ．推論計測による異常事象検出の目的
異常事象や、製品品質等の主要なパラメータ計測値に関しては、以下の内容を把握することが重要である。
・オンラインアナライザまたはラボラトリによって正しい値が提供されているか？
・閉ループプロセス制御で使用されているモデルは品質の値を正しく予測しているか？
・プロセスの変化または異常事象に起因する重大な品質エクスカーションが存在する可能性があるか？
・アナライザまたはラボラトリの値の予測値はどうなっているか？
【００２１】
ＩＩ．推論モデルを使用したオンライン異常事象検出
図２に示すモデル計算プロセスは、以下の４つの主要部分に分割することができる。
Ａ．リアルタイムのデータ収集および前処理
Ｂ．モデルおよび統計的有意性の計算
Ｃ．プロセスオペレータに対する複数のモデル結果の要約
Ｄ．事象解析に関する診断表示の提供
【００２２】
Ａ．リアルタイムデータ前処理
データは様々な周期でシステムに取り込まれる。標準的なプロセス計測値（温度、圧力、流量等）は、通常、非常に頻繁に利用可能となり、少なくとも１分あたり一度の頻度でシステムに取り込まれる。特定の解析に左右されるが、オンラインアナライザからの値は、通常、１０〜３０分毎に利用可能となる。ラボラトリの結果は、通常、１日に一度〜１週間に一度の頻度で利用可能である。ラボラトリの結果を除く他の計測値は全て規則的なサンプルスケジュールでシステムに取り込まれる。ラボラトリの結果については、現場サンプルの収集およびラボラトリ解析の実行において手作業が必要となるため、サンプルサイクルに大きなばらつきが発生し得る。
【００２３】
次に行うべきステップは、データの全体エラー検出である。現場装置によって異なるが、ＨＡＲＴ（登録商標）プロトコルまたはＦｏｕｎｄａｔｉｏｎ  Ｆｉｅｌｄバスプロトコル（Ｈａｒｔ  Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ  Ｆｏｕｎｄａｔｉｏｎの商標であり「ＨＡＲＴ（Ｈｉｇｈｗａｙ  Ａｄｄｒｅｓｓｉｂｌｅ  Ｒｅｍｏｔｅ  Ｔｒａｎｓｄｕｃｅｒ）」と表現する）をサポートする現場装置から付随診断情報が提供される場合がある。生のモデル入力信号を送信機のゼロおよびスパンと照らし合わせて、信号がアクティブであり送信機の範囲内に入っていることを確認する。このチェックでは次の２種類のケースを扱う。ＤＣＳシステムには、４〜２０ミリアンペアの信号範囲から外れている入力信号に対して「不正な計測値」フラグを使用するものもあれば、入力信号が再度４〜２０ミリアンペアの範囲内に収まるまで、この範囲内にある最後の有効な値を維持するものもある。更に、システムでは、ＤＣＳで使用されているフラグを認識し、例えば、ＤＣＳ計測ステータスを「走査停止」または「手動」に設定することにより、現場装置が現場保守またはキャリブレーションを受けていることを示す。
【００２４】
最小および最大運転範囲チェックに加えて、現場装置の機能に関する診断情報を提供するコンピュータによってオンラインアナライザを実行することが可能である。オンラインアナライザの場合は、現在のサンプル値と前の結果を比較することが可能であり、これにより、新しい値の変化量は新しい値が伝送されたことを確認するのに十分であったかどうかがチェックされる共に、新しい値の変化量が不適切な量でなかったかどうかがチェックされる。この不適切な変化に関するチェックは、新しい値を推定値と比較することにより、および生のバイアス（Ｒａｗ  Ｂｉａｓ）の変化量が不適切でないことをチェックすることにより、補完することができる。
【００２５】
モデル入力、オンラインアナライザ値、またはラボラトリ値が、任意の全体エラーチェックで不良と判定された場合、オペレータにはその旨が直ちに通知される。障害が継続する場合、オペレータは、閉ループ制御アプリケーションを無効にする等のフォローアップ操作を行う、アナライザの現場チェックを要求する、入力の計測器チェックを要求する、またはラボラトリによる新規現場サンプルの採取およびその解析を要求することが必要となる場合がある。
【００２６】
生のモデル入力／出力信号が全体エラーチェックに合格すると、モデル入力として必要な数学的変換が計算される。アナライザまたはラボラトリの結果においても、推論モデルの形態と一貫性を持たせるために数学的変換が必要な場合がある。
【００２７】
最後に、モデル入力はフィルタリングされ、過大なノイズが除去される。スパイク挙動を伴う計測値は、ユニオンフィルタ等のスパイク除去フィルタで処理することができる。一方、高周波ノイズは高周波除去フィルタ（通常は指数フィルタ）を使用して低減することが可能である。
【００２８】
Ｂ．モデルおよび統計的有意性の計算
オンラインアナライザおよびラボラトリ計測に対する推論モデルは、通常、部分最小二乗（ＰＬＳ）回帰、ニューラルネット、または段階的線形回帰を使用してフィッティングされる。図２では次の３種類の異常事象が識別される。第１は品質計測値自体における差し迫った問題（Ｙｅｓｔｉｍａｔｅ値を使用して識別）、第２は異常事象によって引き起こされる重大な品質エクスカーションの予測（Ｙｐｒｅｄｉｃｔｅｄ値を使用して識別）、第３はモデルまたは生の計測のいずれかにおける長期的なドリフト問題（生のバイアス（Ｒａｗ  Ｂｉａｓ）を使用して識別）である。モデルの別の用途として、プロセス制御アプリケーションで使用する推定値（Ｙｃｏｎｔｒｏｌ）を算出することができる。
【００２９】
推論品質尺度の計算（Ｙｅｓｔｉｍａｔｅ、Ｙｐｒｅｄｉｃｔｅｄ、Ｙｃｏｎｔｒｏｌ）
この３種類の推論品質尺度を計算するためのモデルは、以下の３つの主成分から構成されている。
・モデル入力と品質計測との時間同期
・推論モデル計算
・バイアス計算と更新
【００３０】
バイアス計算の副次的な結果として、モデル計算の問題またはアナライザのキャリブレーションドリフトを識別するために生のバイアス（Ｒａｗ  Ｂｉａｓ）が解析される。
【００３１】
時間同期
品質計測値は、通常、プロセスの変化に対して応答するのが最も遅い。従って、モデルへの入力は、品質計測と時間同期される必要がある。この時間同期は、図２の中で用語「Ｇ」で示している。これは、微分差分方程式またはラプラス変換伝達関数によって非常に簡単に表現される動的モデルであり、以下に示すような一次むだ時間モデル等がある。
【数１】
2010519636.tif 000003
式中、
Ｘ（ｓ）  −  入力信号のラプラス変換
Ｙ（ｓ）  −  出力信号のラプラス変換
Ｇ（ｓ）  −  伝達関数の記号
ｅ−ＴＤ＊ｓ  −  純粋なむだ時間ＴＤに関するラプラス変換
Ｔａｕ  −  一次常微分方程式の時定数
ｓ  −  ラプラス変換パラメータ
【００３２】
このモデルには定常状態ゲイン項がないことに留意されたい。定常状態の関係は全て、推論モデルのパラメータに組み込まれている。モデルの開発段階を通じて、この時間同期モデルＧ（ｓ）は、各候補入力と品質計測値との間の最良な単入力単出力動的モデルとして開発されている。多入力単出力（ＭＩＳＯ）動的モデルを開発する必要はない。
【００３３】
Ｙの３つの推定のそれぞれで使用される時間同期は、異なるものであるが、関連付けられている。
・Ｙｅｓｔｉｍａｔｅ計算では、種々の入力と品質計測を時間同期するために動的モデル「Ｇ」が使用される。
・Ｙｃｏｎｔｒｏｌ計算では、動的モデル「Ｇ」の代わりに、むだ時間項ｅ−ＴＤ＊ｓを持たない動的モデル「Ｇ＊」が使用される。尚、他のモデルパラメータは全て「Ｇ」の場合と同じである。
・Ｙｐｒｅｄｉｃｔｅｄ計算では、入力は時間同期されない。時間同期が存在しないので、全ての入力が現在値のまま維持される場合、Ｙｐｒｅｄｉｃｔｅｄ値は「Ｙ」の最終的な定常状態値の推定値となる。
【００３４】
推論モデル計算
３種類の推論推定の全てについて、推論モデルは同じである。モデルの計算には、カルマンフィルタ、部分最小二乗（ＰＬＳ）回帰、ニューラルネット、または段階的回帰等の多種多様なモデルを使用できる。しかし、好ましいのはＰＬＳモデルである。入力が、通常、相互相関しており、計算が非常に簡単であり、モデル結果への各入力の寄与（ｃｏｎｔｒｉｂｕｔｉｏｎ）が容易に特定できるためである。変換によって容易に補償できないほど強力な非直線性をプロセスが有する場合は、ニュートラルネットモデルがふさわしいかもしれない。
【００３５】
ＰＬＳモデルの場合は、オンライン計算のために複雑な潜在変数構造を維持する必要がない。元々の潜在変数の形態では、潜在変数Ｌｋは、以下のように平均中心の、ユニットバリアンスでスケーリングされた入力変数
【数２】
2010519636.tif 000004
の一次結合で形成されている。
【数３】
2010519636.tif 000005
式中、
ａｉ，ｊ  −  入力「ｉ」と潜在変数「ｊ」の潜在変数ローディング（係数）
Ｂｎ  −  潜在変数「ｎ」の回帰係数
Ｌｋ  −  ＰＬＳモデルにおける重要な潜在変数
ＹＮ  −  時間＝ＮでのＰＬＳモデル出力
Ｘｊ，Ｎ  −  時間同期した生のまたは変換された入力
ＬＶ＿Ｂｉａｓ  −  モデル開発に基づくモデル定数
【００３６】
入力ごとにファクタおよび係数の全てを単一の係数にまとめることにより、一般的なＰＬＳモデル計算をかなり簡単に行うことができる。
【数４】
2010519636.tif 000006
式中、
σｊ  −  Ｘｊの標準偏差である。
【００３７】
バイアス計算と更新
品質を推定するモデルは、モデルの誤差を補償するため、および実際の計測値に影響を及ぼす計測されないまたはモデル化されない結果を補償するために、継続的に更新される必要がある。モデルの付加的なバイアス項を更新するのが最も簡単な手法である。バイアスの初期値は、モデル開発段階で決定されるモデル定数である。これは、モデル推定値と実際の値（ここでは「生のバイアス（Ｒａｗ  Ｂｉａｓ）」と称する）との残差に基づいて更新される。
生のバイアスＮ＝品質計測値Ｎ−品質推定値Ｎ      式８
【００３８】
この計算が実行されるのは、新しいオンライン解析値またはラボラトリ解析値が提供され、全てのモデル入力が全体エラーチェックに合格し、且つ新しい品質計測値（Ｑｕａｌｉｔｙ  Ｍｅａｓｕｒｅｍｅｎｔ）がその全体エラーチェックに合格した場合にのみ実行される。オンラインアナライザの多くは、低頻度（通常は１０〜３０分毎であり、解析に左右される）でしか新しい解析値を提供できない。解析値の変化は、新しい値が提供されたことを示すトリガとして使用可能である。
【００３９】
オンラインアナライザは、解析装置の診断チェックを実施するまたはキャリブレーションサンプルを時々実行するコンピュータによってしばしば実行される。バイアスの更新は次の場合には実行されるべきでない。診断チェックで問題が示された場合、プロセスサンプルではなくキャリブレーションサンプルが実行された場合、サンプル結果が不当である場合（サンプル間で値の変化がない、またはサンプル間の変化が大きすぎる場合）、またはプロセスが異常状態（例えば、シャットダウン状態）にある場合である。通常、プロセスストリームサンプリングからオンライン解析結果の提供までの時間は、サンプル間で一貫性があるので、このむだ時間は時間同期化モデルの一部となる。
【００４０】
品質計測用の新しい値がラボラトリ解析からもたらされる場合は、現場サンプルが採取された時間を使用して、生のバイアス（Ｒａｗ  Ｂｉａｓ）の計算のために正しい比較時間を決定する必要がある。手動による現場サンプルからラボラトリ結果の提供まで時間は、大幅に変動する可能性がある。従って、ラボラトリの結果の準備が整うと、このむだ時間が計上される。このため、現場のオペレータが現場サンプルを時間「Ｋ」で採取し、ラボラトリ解析値が時間「Ｎ」に届く場合、以下の式が成立する。
生のバイアスＮ＝品質計測値Ｎ−品質推定値Ｋ      式９
【００４１】
「生のバイアス（Ｒａｗ  Ｂｉａｓ）」から「バイアス（Ｂｉａｓ）」を更新する従来アルゴリズムとして複数のものが存在する。中には、指数フィルタ等のローパスフィルタを使用して生のバイアス（Ｒａｗ  Ｂｉａｓ）を処理するものがある。
バイアスＮ＝ａ＊バイアスＮ−１＋（１−ａ）＊生のバイアスＮ      式１０
ここで
【数５】
2010519636.tif 000007
サンプル期間  −  ある解析結果から次の解析結果までの時間
フィルタ時間  −  指数フィルタの時定数
【００４２】
プロセス制御アプリケーションに使用する場合とは対照的に、推論モデルを異常事象検出に使用する場合は、フィルタ時間をサンプル時間に比べて長い時間間隔に設定する必要がある。通常は、サンプル時間の２０倍を超える長さに設定する。例えば、オンラインアナライザのサンプル時間が１０分である場合は、フィルタ時間を３時間より長く設定する。これにより、「ａ」の値は０．９５より大きくなる。
【００４３】
更新されるバイアスをどのアルゴリズムで計算するかに関わらず、異常事象検出の主な要件は、異常事象の影響がバイアス更新計算に吸収されることによって見えなくならないように長いフィルタ時間を指定することである。異常事象の影響が品質計測値に現れるのが、指定されたフィルタ時間より遅くなればなるほど、異常事象の検出も益々難しくなる。従って、検出するのが最も難しい異常事象は計測値の低速ドリフトであり、この場合はフィルタ時間も同等に遅くする必要がある。
【００４４】
Ｃ．ファジーネットの使用によるモデル結果の正規化およびモデルの結合
異常事象検出の場合、異常事象の可能性があることをいつ結果は示しているかをコンソールオペレータが把握できるように、数値モデルの結果にはコンソールオペレータ用の解釈が必要である。この解釈を作成するためのメカニズムでは、簡単な「ファジファイア（ｆｕｚｚｉｆｉｅｒ）」およびモデル残差の統計的変動（標準偏差で表現されるような）を使用して、異常事象の可能性の推定値を作成する必要がある。
【００４５】
図３に、実施方法の例を示す。実線は「ファジファイア」として知られる数学的関数である。この関数では正規化されたモデル残差（残差の標準偏差に基づいて正規化されている）を、異常事象が存在する可能性に変換する。システムユーザは、モデル残差結果が明らかに正常である領域を選択する。この例では、モデル残差が標準偏差−０．５と０．５の間にある場合は明らかに正常であり、モデル残差が標準偏差５．５より大きいか、標準偏差−５．５より小さい場合、モデル残差は明らかに異常性を示している。これらの領域間における異常性の可能性は、１／２標準偏差での可能性ゼロから標準偏差５．５での可能性１００％までの線形補間となる。点線および破線は可能性の代表値（６０％と９０％）であり、ここでコンソールオペレータは、異常状態の可能性を通知されることになる。
【００４６】
異常性の具体的な確率は、残差のデータ分布関数によって決定される。表１は、平均値の種々の標準偏差内に収まる正常値のパーセンテージを示す（データ分布関数はガウス分布または正規分布であると想定する）。
【００４７】
【表１】
2010519636.tif 000008
【００４８】
しかし、残差のデータ分布は、通常、ガウス分布または正規分布に適合しない。むしろ、正常／定常状態運転中の残差分布とプロセスの軽度の不調時の残差分布という２つの確率分布が複合化したものである。このような事情にも関わらず、上表は適切な設定を概念的に案内する際には依然として有用である。
【００４９】
モデル残差が「異常事象の存在の可能性」という一般的な尺度に変換されると、複数のモデル残差を結合してプロセス運転の広範囲を要約することができる。これを行うには、多数のモデルの中から最も高い異常の可能性を選出する。
【００５０】
品質解析計測の問題点の検出
品質計測での問題を検出するには、実際の品質計測値Ｙと推定された品質計測値Ｙｅｓｔｉｍａｔｅとの間の統計的に大きな差を探索する。この差は、品質推定残差である。この残差の統計的な正常域は、軽度のプロセス不調時に残差の標準偏差から算出する必要がある。図３の例に示すように、「正常域」である可能性が高い領域を実際に定義すると、データの９９．９％が存在する範囲ということになる。この一組の限界点は、ゼロ（この残差の平均）からおよそ３．５標準偏差離れた点に対応する。これらの限界点を変更して、感度を上げたり下げたり調整することができる。
【００５１】
典型的な異常事象検出（ＡＥＤ）のオンライン実装の場合、全ての解析間での複合的な異常の可能性は、個々の可能性から計算される（図４を参照）。全ての解析を対象とした正規化された品質推定残差（残差を正常運転中のその標準偏差によって割り算することにより正規化される）は、ファジーペトリネットに入力される。ファジファイアでは、正規化された各残差を異常の可能性に変換する。ここで、−０．５ＳＤと０．５ＳＤとの間にある値は１００％正常であり、５．５ＳＤより大きいか、または５．５ＳＤより小さい値は１００％異常である。異常の最大の可能性を選択することにより、全ての解析問題が要約される。可能性が６０％（３．５標準偏差に対応する）を超えている場合、オペレータには、品質計測問題が予想される旨の第１の通知が与えられる。
【００５２】
重大な品質エクスカーションの予測
将来的に重大な品質エクスカーションが発生する可能性は、Ｙｐｒｅｄｉｃｔｅｄ．の値と品質限度とを比較することにより推定される。モデルの品質によっては、Ｙｐｒｅｄｉｃｔｅｄ．の値とＹの実際の将来値との間には若干の誤差が発生する。この誤差の推定には、品質推定残差（Ｙｅｓｔｉｍａｔｅ−Ｙ）の確率分布が最適である。ガウス分布を想定した場合、図５はＹの将来値とＹｐｒｅｄｉｃｔｅｄとの間のばらつきを示す。この例では、Ｙｐｒｅｄｉｃｔｅｄは１．３標準偏差であり、品質上限を超えている。このＹｐｒｅｄｉｃｔｅｄ値については、Ｙの将来値の９０％は品質上限を上回り、１０％が品質上限を下回る。特定の可能性推定値について、この表２では、特定の可能性値について必要とされる上限を超えるオフセットを示している。
【００５３】
【表２】
2010519636.tif 000009
【００５４】
図６に示すファジーネットでは、この方法がオンラインシステムでどのように実施されるかを示している。ファジーネットへの入力は、正規化された、Ｙｐｒｅｄｉｃｔｅｄ．と品質上限との差である。
【数６】
2010519636.tif 000010
【００５５】
ファジファイアでは、この値を可能性テーブルに基づいて将来的な限度違反の可能性に変換する。Ｙｐｒｅｄｉｃｔｅｄ．が、品質上限より下で１．８３標準偏差より小さい場合、運転が品質限度内に収まる可能性の推定値は１００％となる。Ｙｐｒｅｄｉｃｔｅｄ．が、品質上限より上で１．６４標準偏差より大きい場合、品質上限を超える可能性の推定値は１００％となる。このようなコーナーポイント値は厳密には正しくないが、これらの値により、ファジファイアでは、２．５５標準偏差（オペレータが第１の通知を取得するポイント）で６０％の可能性、そして１．３標準偏差で９０％の可能性を正確に算出することができる。
【００５６】
生のバイアス（Ｒａｗ  Ｂｉａｓ）を調べてモデルの問題／キャリブレーションのドリフトを確認
モデルの問題または長期のキャリブレーションドリフトは、生のバイアス（Ｒａｗ  Ｂｉａｓ）の値に現れる。この２つの問題は、互いに区別するのが困難ないし不可能な場合がある。良好な推定機能を使用してモデルが稼動状態にされていると仮定すると、モデルの問題は、品質値に実際に影響するのにモデル内に含まれない何らかの要因における変化、相関構造の変化、またはプロセス運転条件の変化に起因する入力と出力との間の時間動特性の関係の変化に起因する。キャリブレーションのドリフトは、ゼロと計測値の範囲の両方またはいずれかに影響を及ぼす場合がある。
【００５７】
これらの問題は、生のバイアス（Ｒａｗ  Ｂｉａｓ）の統計的に大きな変化であって、恒久的な変化または統計的有意性の限界の上下を循環する変化をチェックすることにより検出される。品質解析計測問題を検出する場合と同様に、これらの問題のチェックでも、いつ、生のバイアス（Ｒａｗ  Ｂｉａｓ）の値がその正常範囲から外れ、その範囲外に複数時間とどまったかを調べる。循環チェックでは、生のバイアス（Ｒａｗ  Ｂｉａｓ）が統計的有意性の限界を超えた回数をカウントする。生のバイアス（Ｒａｗ  Ｂｉａｓ）の統計的に正常な範囲は、正常運転中のその標準偏差から算出可能である。「正常範囲」の一般的な定義を示すならば、データの９９．９％がある範囲となろう。この一組の限界点は、生のバイアス（Ｒａｗ  Ｂｉａｓ）の予測平均であるゼロからおよそ３．５標準偏差離れた点に対応する。これらの限界点を変更して、感度を上げたり下げたり調整することができる。
【００５８】
統計的限界から逸脱した状態が続いた場合、キャリブレーションサンプルを使用したアナライザのチェック、更に必要があればアナライザの再キャリブレーションを計測器の技術者に行わせるように、オペレータは指示される。
【００５９】
Ｄ．オペレータ用インターフェイスディスプレイ
オペレータ用ディスプレイの目的は、潜在的な異常事象が検出されたことをオペレータに通知し、オペレータがその状況を診断するのに必要な支援情報を迅速に提供することにある。この情報は、図７に示すように、一連のトレンドおよび分布棒グラフによって提供される。
【００６０】
トレンド１には、将来的な異常品質エクスカーションの可能性をオペレータに知らせるための情報が表示される。このトレンド例では、時間４１４５あたりでエクスカーションの発生が予測されている。点線はエクスカーションの６０％の可能性に対応し、一方、実線はエクスカーションの９０％の可能性に対応している。この例では、エクスカーションは、予測したように発生する（しかし、必ずしも完全に予測された範囲ではない）。
【００６１】
トレンド２では、推定されたアナライザ値と実際のアナライザ値との比較を示している。このトレンドは、オペレータに推定値の品質を植え込むことに加えて、実際のアナライザで障害が発生した場合にオペレータが推定値を把握する必要がある場合にも有用である。
【００６２】
トレンド３では、推定された値と実際の値との差が統計的に大きいかどうかをオペレータに示す。点線はアナライザ問題の６０％の可能性に対応し、一方、実線はアナライザ問題の９０％の可能性に対応している。
【００６３】
棒グラフ４では、どの計測において、アナライザ値の平均値からのアナライザ値の変化が示されているかを表示している。これらの入力は、アナライザ値の変化の推定値への寄与が高いものから順に配置されている。算出されるアナライザの変化への各モデル入力の寄与は、式５における当該変動から得られる。
（ＹＮ−ＹＮ，ｆｉｌｔ）＝Ａ１＊（Ｘ１，Ｎ−Ｘ１，Ｎ，ｆｉｌｔ）＋Ａ２＊（Ｘ２，Ｎ−Ｘ２，Ｎ，ｆｉｌｔ）＋．．．Ａｋ＊（Ｘｋ，Ｎ−Ｘｋ，Ｎ，ｆｉｌｔ）      式１３
入力Ｘｊと関連付けられた実線バーのサイズは次のとおりである。
実線バーのサイズｊ、Ｎ＝Ａｊ＊（Ｘｊ，Ｎ−Ｘｊ，Ｎ，ｆｉｌｔ）      式１４
点線の参照ボックスのサイズ（予想される寄与を示す）は、トレーニングデータセットからのこの項の標準偏差である。
点線ボックスのサイズｊ＝標準偏差（Ａｊ＊（Ｘｊ，Ｎ−Ｘｊ，Ｎ，ｆｉｌｔ））トレーニングデータ      式１５
【００６４】
ＩＩＩ．推論品質モデルの開発
正確なモデルの開発で鍵となるのは、高品質のトレーニングデータセットを備えることである。高品質のトレーニングデータセットを決定する指標は次のとおりである。
・トレーニングデータのＳ／Ｎ比が良好である
・計測可能な、解析への全ての影響を含む。
・生データの適切な数学的変換（特に時間同期化で必要）を含む。
・異常なプロセス運転に関する注釈を含む。
・プロセス／制御システム構成の補助的指標を含む。
・対象とする運転範囲に及ぶ。
・種々の運転点のバランスの取れた表現を含む。
・実際のプロセス状況を表現していない不良データ（アウトライアー）を除外する。
【００６５】
図８では推論モデルを開発する典型的な従来技術の方法を示し、一方、図１では異常事象検出のため、並びにプロセス時間動特性と相互関係したモデル入力とを補償するために必要な標準的な方法に対する拡張機能を示す。図１に示す方法を構成する各ステップの内容を以下に説明する。以下のパラグラフに振られた番号は、図１の振番されたボックスに対応している。
【００６６】
１．推論計測運転領域の識別
このステップでは、履歴プロセス運転が、計画された推論モデル開発と関連していた時期、およびモデルが適用される運転モードを識別することが必要である。これらのモデルは経験的なものであるため、データ内の相関構造に依存する。その相関構造は、運転モードに応じて変化する可能性がある。この段階では、運転内に存在する運転モードと、この運転モードでプロセスが運転されていた時期を識別する必要がある。その後、モデル構築段階では、これらの運転モードのうち、どの運転モードを結合すれば、運転をカバーするのに必要な合計のモデル数を削減できるかを決定する必要がある。更に、品質の尺度に影響を与えると考えられる装置構成モードが存在する場合もある。
【００６７】
連続工業プロセスでは、必要な保守のために、まれに周期的な変更を経験する。変更中、現場では、ボトルネックの解消または他の理由でプロセス装置のアップグレードを選択する場合がある。この種のプロセス変更により、変更前の運転からのデータが、モデル開発に対して無効にされる場合がある。従って、プロセスに加えた重要な機械的変更の日付、またはプロセス制御方策に加えた主要な変更の日付は、モデル開発のためにデータを収集可能な最も早い時期として留意する必要がある。
【００６８】
２．ラボラトリ解析におけるサロゲート計測値の識別
このステップでは、時間同期化の計算を目的としてオフライン品質計測値の「代わりに使用できる」オンライン計測値を選択することが必要である。ラボラトリサンプルはまれにしか取得されないので、ラボラトリ計測値と同時に変化するとみられる類似のオンライン計測値に依存する必要がある。例えば、９５％沸点は、最終沸点ＦＢＰと強く関連付けられているとする。ＦＢＰ推論モデルでは９５％沸点と強く関連付けられた計測値を使用するとみられ、この強く関連付けられた計測値と９５％沸点との時間同期は、ＦＢＰとの時間同期と同じになると想定される。
【００６９】
サロゲートオンライン計測値を選択する場合、オンラインアナライザの自動サンプル点がラボラトリ解析の手動サンプル点と物理的に近接しているかどうかに注意することが重要である。このような２つのサンプル点の配置間隔が大きくなるほど、時間同期で発生する誤差も大きくなる。物理的分離が管の長さだけである場合は、介在する管の容量を正常プロセスフローで除算した結果に相当するむだ時間を、時間同期モデルに加算（手動サンプリング点が自動サンプリング点の下流にある）または時間同期モデルから減算（手動サンプリング点が自動サンプリング点の上流にある）することにより、これを補正することができる。
【００７０】
サロゲートオンライン計測もラボラトリ解析の一部である場合は、オンライン計測値とラボラトリ計測値とを比較することにより、サンプル位置の違いに起因する時間シフトを算出することができる。
【００７１】
３．プロセスデータの収集
この方法では、正常運転のデータを使用して推論モデルを作成することが可能なため、プロセス用に計画された実験的試験にかかるコストを回避できる。計画された実験によってプロセスが誘起されない場合、このモデルを構築するのに必要な情報はプロセスの不調および運転点の変化からもたらされる。
【００７２】
連続工業プロセスは、主に、限られた数の定常状態運転点でしか運転されない。運転のこれらの定常期間は、プロセス計測間の関係に関する情報にほとんど寄与しない。従って、プロセス事象を十分に捕捉してこれらの計測関係をモデル化するために、長期にわたる運転データの収集が必要である（例えば、６ヶ月〜２年）。プロセス履歴システムは長年にわたるデータを格納する（しばしば高分解能で）ので、データのソースとなる可能性が高い。
【００７３】
通常、生のオンライン計測データは、データ圧縮を行わずに、およそ１分間に１度の割合でスポット値として収集する必要がある。モデルに含める候補である計測値に加えて、プロセス性能、プロセス運転モード、装置構成、およびプロセス制御構成の指標である変数も収集する必要がある。
【００７４】
異常プロセス事象または大きな外乱に関する注釈もプロセス計測データと共に収集する必要がある。
【００７５】
４．異常データサンプルの識別
データ内に不正な値が多数含まれることがよくある。一般に識別される内容は、飽和／オフスケール計測、計測器のキャリブレーション、および履歴システムデータ収集障害である。通常、中度から重度のプロセス外乱は、プロセス運転のシャットダウンまたは重大な変化を引き起こさない限り、異常データとして分類されない。
【００７６】
５．評価：サンプルはアウトライアーであるか
データサンプルは、その値が計測される実際のプロセス条件を表していない場合に限り（例えば、送信機の範囲外でプロセスが運転されていた）、アウトライアーとして判断されるべきである。品質計測値が不良である期間はいずれも、破棄する必要がある。但し、最終的には少数の候補モデル入力しかモデルに取り込まれないので、モデルトレーニングデータに取り込み可能なサンプル時間は評価対象の特定の入力セットに応じて変化する。
【００７７】
６．定常運転の期間の識別
定常運転期間からの限られた量のデータだけがトレーニングデータセットに含められる。共通の定常運転点が複数存在する場合は、それぞれの定常運転点から同様のデータ量を選択する必要がある。定常運転からのデータ量の合計は、通常、合計データセットの２５％程度とするべきである。
【００７８】
７．大部分の定常状態運転の削除
データフィルタを構築することにより、アルゴリズムで定常運転を識別し、削除することができる。このフィルタに基づいて、品質値の変化率（ＲＯＣ）または他の重要な運転性能指標を判定するのが最も一般的である。このようなデータフィルタは、以下の計算を行うことにより構築できる。
ＲＯＣＮ＝ＹＮ−ＹＮ，ｆｉｌｔｅｒｅｄ      式１６
ここで
ＹＮ，ｆｉｌｔｅｒｅｄ  −  ローパスフィルタにかけられたＹＮの値
データフィルタ  −  下限≦ＲＯＣＮ≦上限を除外
【００７９】
フィルタリングされる項ＹＮ，ｆｉｌｔｅｒｅｄは、現在の平均運転点の推定値として働く。現在の値ＹＮからＹＮ，ｆｉｌｔｅｒｅｄを引き算すると、このデータ級数は、平均値がゼロに近い静止級数（ＲＯＣＮ）に変換される。ゼロに近いＲＯＣＮの値は定常運転を示し、この値に情報量はほとんどない。下限と上限の一般的な選択肢は次のとおりである。
下限＝−標準偏差ＲＯＣ
上限＝＋標準偏差ＲＯＣ
このように、定常運転の大部分は削除されている。
【００８０】
図９では、定常運転を削除しない場合の結果を示している。結果として得られた推論モデルでは、変化が小さい場合は良好な推定特性を示すが、変化が大きい場合に示す推定特性は不十分なものとなろう。
【００８１】
８．入力および出力における正常な範囲の識別
これらの範囲を識別する目的は、計測値が飽和している可能性がある（例えば、送信機の範囲外で）データ、または非常に非線形な運転の領域でプロセスが運転されている（例えば、高いタワーデルタ圧力によって示されるタワーフラッディング）データをアルゴリズムによって除外することにある。この情報はまた、オンラインシステムの全体エラー検出計算でも使用される。
【００８２】
９．データフィルタの作成
前述のフィルタに加えて、追加のデータフィルタを作成することにより、極めて異常な運転の時間中のデータを自動的に除外することができる。
【００８３】
１０．評価：不十分なＳ／Ｎ比
Ｓ／Ｎ比が不十分となる原因は、データ内に強力な信号が不足している（計測値をその設定点の近くに維持する良好なプロセス制御が原因となる場合が多い）ことにあると考えられる。または、信号内の非常に大きなノイズが原因となる場合もある。強力な信号の不足が原因である場合、これらの尺度はモデルから効果的に外される。
【００８４】
１１．平滑化フィルタ／スパイクフィルタを使用したデータのフィルタリング
最も一般的に必要とされるフィルタは、スパイクフィルタ、ミッドパスフィルタ、およびローパスフィルタである。スパイクフィルタは、計測値の変化を、予め設定された最大変化までと制限する場合がある。データに適用されるフィルタはいずれもオンラインシステム内で再生される必要があるため、これらは物理的に実現可能とする必要がある。フィルタは、位相遅れを最小限に抑えるように選択する。
【００８５】
１２．生データの数学的変換の作成
モデルのフィッティングを向上させるために、複数の周知の変換が採用される。例えば、組成アナライザおよびタワーオーバーヘッド圧力のログ変換、無次元比への流量の変換、トレイ温度への圧力補償の適用等がある。
【００８６】
１３．ＳＩＳＯ動特性モデルの構築
データ内の動特性を処理するために、各候補入力とアナライザ値との間の動的モデルが開発される。これらは、信号入力・信号出力モデル（ＳＩＳＯ）であり、通常は、単なる一次モデルにむだ時間モデルを加えたものである。
【数７】
2010519636.tif 000011
式中、
Ｘ（ｓ）  −  入力信号のラプラス変換
Ｙ（ｓ）  −  出力信号のラプラス変換
Ｇ（ｓ）  −  伝達関数の記号
ｅ−ＴＤ＊ｓ  −  純粋なむだ時間ＴＤに関するラプラス変換
Ｔａｕ  −  一次常微分方程式の時定数
ｓ  −  ラプラス変換パラメータ
【００８７】
多入力単出力（ＭＩＳＯ）動的モデルの開発は必要でも望ましくもない。入力は、通常、相互相関しているので、ＭＩＳＯモデリング結果はモデリングケースの運用のために選択された特定の入力セットに大きく影響される。
【００８８】
１４．データからの低速トレンドおよびバイアスの除去
計測されない恒久的な負荷外乱がデータ内に存在しない場合、微分データのモデリングからの係数は、十分に評価されたデータのモデリングからの定数と類似すべきである。計測されない恒久的な負荷外乱は、モデルに極度にバイアスをかける可能性がある。このバイアスに関する問題を回避するために、モデルは十分に評価されたデータではなく微分データの形式を使用して開発される。
【００８９】
この変換を実行する特定の方法では、指数フィルタ等のローパスフィルタを使用して、フィルタリングされた値を現在の値から引き算する。
【数８】
2010519636.tif 000012
式中、
Ｙ  −  品質計測値
【数９】
2010519636.tif 000013
−  フィルタリングされた値が差し引かれた品質尺度
Ｘｉ  −  候補入力ｉ
【数１０】
2010519636.tif 000014
−  フィルタリングされた値が差し引かれた候補入力
ａ  −  指数フィルタ定数（ＸおよびＹについて同じ値）
【００９０】
指数フィルタ定数は、品質に影響するプロセスの支配的な時定数の１〜３倍となるように選択すべきである。これはステップ１３で構築されたＳＩＳＯモデルから推定可能である。図１０に示すように、品質計測値は、プロセス変化に対する最終的な応答の６４％〜９５％で、１〜３倍の時定数の後にそれぞれ応答する。これより短い時定数を選択すると、高周波ノイズが増幅され、信号強度が低下する。これより長い時定数を選択すると、低周波の恒久的な外乱が増幅される。フィルタ定数が１〜３倍の時定数の範囲内にあれば、Ｓ／Ｎ比は最も改善される。
【００９１】
１５．推論モデルの構築
前のステップで作成したトレーニングデータセットを使用すれば、ツールおよび方法を構築する標準的なモデルによって推論モデルを構築できる。この目的は、良好な推論モデルをもたらす最小数の入力を使用することにある。式１７を参照すると、ほぼ同等の推定強度を有する相関した代替入力の中から選択する場合は、ＴａｕおよびＴＤに対して最も大きな値を有する入力を選択することにより、推論モデルは異常エクスカーションをより早く予測し、プロセスを場合によってはより適切に制御することができるようになる。
【００９２】
例
図１１は一組の蒸留塔を示す。ここで、第１の塔Ｔ−１０１はエタン、メタン、水素等の軽質炭化水素を除去するように設計され（脱エタン塔（ｄｅ−ｅｔｈａｎｉｚｅｒ）として知られている）、第２の塔Ｔ−２０１はＣ３炭化水素（主にプロパン）を抽出するように設計されている。
【００９３】
検出される必要がある異常事象が複数ある。タワーＴ−２０１の運転にとって重要なのは、Ｔ２０１から取り出される蒸留液の組成である。これはアナライザＡ２０１で計測され、推論コントローラ「Ａ−２０１コントローラ」によって制御される。
【００９４】
このアプリケーション用の推論推定器を設計する第１のステップは、検出する異常事象のタイプを把握することである。この場合、推定器は、アナライザサンプルループのプラギングやアナライザキャリブレーションでのドリフティング等のアナライザ運転における問題を識別し、ポンプＰ−２０１の吸引時の過剰な軽質炭化水素を予測することによりポンプＰ−２０１のキャビテーションを防止する。
【００９５】
アナライザの問題を検出する場合は、現行の工業的手法に従って、Ｔ−２０１の温度、圧力、および流量の計測値を基に推論を構築する。アナライザの問題は、「生のバイアス」、アナライザ信号、およびモデル推定の時間トレンドを組み合わせることにより明らかにされる。
【００９６】
図１２は、アナライザへのサンプルラインの部分的なプラギングと関連付けられた設定時間トレンドを示す。この例では、プラギングはサンプル時間５００付近（グラフ上では破線の垂直線で示す）で生じるが、データではわずかな兆候（アナライザ信号の正常な変動性より低い（図１２Ｂ））しか見えない。工業的手法ではサンプル値の間の最小の変化によってアナライザ信号を検証することができるが、この例に示すような小さな変動の場合は、この種の検証チェックはトリガされる可能性が低い。
【００９７】
サンプル時間５００とサンプル時間８００との間では、プラギングされたサンプル示度と実際の組成とがかなり厳密に一致している。これは単に、プラギングが生じて以降、プロセス組成が大して変化しておらず、組成の設定点付近で運転されているからである。これは運転に潜む落とし穴である。
【００９８】
プラギングが唯一影響するのは、サンプル時間８００で組成が大きく動き始めたときである。時間８００の後、モデル推定アルゴリズムのバイアス更新部のために、アナライザ値とモデル推定値はかなり厳密に一致しているが、両方とも正しくない。アナライザとモデル推定は互いに確認し合うため、閉ループ制御またはオペレータのいずれの操作もこの不正な情報に基づくことになる。プラギングが部分的であるために、プロセス流体の一部はアナライザに漏出する。これにより、アナライザは実際の組成の方向へゆっくりと移動することができる。
【００９９】
「生のバイアス」信号を解析すると状況は直ちに明らかになる。統計的に大きな変化は、この場合、３つの標準偏差よりも大きな「生のバイアス」変化として求められ、図１２Ｄにおいて水平に延びる点線で表される。この設定の厳密な値は、不正な指標となるのを避けるために「生のバイアス」の長期履歴を調べて調整され、通常は、履歴信号特性に応じて３つの標準偏差と５つの標準偏差との間に分布し得る。
【０１００】
図１３では、アナライザのドリフティングキャリブレーションと関連付けられた設定時間トレンドを示す。この例では、サンプル時間５００付近で最終的な１０％ドリフトが生じる（グラフでは点線の垂直線で示される）が、アナライザ（図１３Ｂ）またはモデル推定（図１３Ｃ）データでは明らかな兆候は見られない。
【０１０１】
「生のバイアス」信号を解析すると、状況は徐々に明らかになる。サンプル時間７２０までには統計的に大きな変化に達し、この場合、図１３Ｄにおいて水平に延びる点線で表されるように、３つの標準偏差よりも大きな「生のバイアス」変化として求められる。この設定の厳密な値は、不正な指標となるのを避けるために「生のバイアス」の長期履歴を調べて調整され、通常は、履歴信号特性に応じて３つの標準偏差と５つの標準偏差との間に分布し得る。
【０１０２】
アナライザ問題の検出とは対照的に、Ｐ−２０１の吸引時の組成の予測では、ポンプのキャビテーションを防止するのに十分なリードタイムを確保するために異なるプロセス計測の選択が必要である。Ｔ−２０１の計測においてＴ−２０１原料内での軽質炭化水素の増大が示された場合、軽質炭化水素はＴ−２０１内にすでに存在しており、タワーからそれらを取得するにはタワーのオーバーヘッドを通すしかないので、ポンプのキャビテーションを防ぐのは難しい。
【０１０３】
代わって、推論推定器の設計は、Ｔ−２０１への原料内の過度の軽質炭化水素を示す追加のリードタイムを提供するために、Ｔ−１０１での計測に応じて行う必要がある。入力計測の選択は更に上流へと移動するので、推論推定器の精度は劣化する可能性がある。このため、このような入力計測の選択は、上記したようなアナライザ問題の検出には望ましくない。
【０１０４】
特許文献引用
【０１０５】
【表３】
2010519636.tif 000015
【０１０６】
【表４】
2010519636.tif 000016
【０１０７】
非特許文献引用
（１）Ｃａｒｄｏｓｏ，Ｊ．  ｅｔ  ａｌ  “Ｆｕｚｚｙ  Ｐｅｔｒｉ  Ｎｅｔｓ：  Ａｎ  Ｏｖｅｒｖｉｅｗ”，１３ｔｈ  Ｗｏｒｄ  Ｃｏｎｇｒｅｓｓ  ｏｆ  ＩＦＡＣ，Ｖｏｌ．Ｉ：Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  ＩＩ，Ｄｉｓｃｒｅｔｅ  Ｅｖｅｎｔ  Ｓｙｓｔｅｍｓ，Ｓａｎ  Ｆｒａｎｃｉｓｃｏ，ＣＡ，ＵＳＡ，Ｊｕｎｅ  ３０−Ｊｕｌｙ  ５，１９９６，ｐｐ．４４３−４４８
（２）Ｅｒｂａｙ，Ａ．Ｓ．  ＆  Ｕｐａｄｈｙａｙａ，Ｂ．Ｒ．  ”Ａ  Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ  Ｂａｓｅｄ  Ｏｎ−Ｌｉｎｅ  Ｓｉｇｎａｌ  Ｖａｌｉｄａｔｉｏｎ  Ｓｙｓｔｅｍ  ｆｏｒ  Ｎｕｃｌｅａｒ  Ｐｏｗｅｒ  Ｐｌａｎｔｓ”，Ｎｕｃｌｅａｒ  Ｔｅｃｈｎｏｌｏｇｙ，Ｖｏｌ  １１９，ｐｐ．６３−７５，Ｊｕｌｙ  １９９７
（３）Ｇｅｒｔｌｅｒ，Ｊ．，”Ｓｕｒｖｅｙ  ｏｆ  Ｍｏｄｅｌ−Ｂａｓｅｄ  Ｆａｉｌｕｒｅ  Ｄｅｔｅｃｔｉｏｎ  ａｎｄ  Ｉｓｏｌａｔｉｏｎ  ｉｎ  Ｃｏｍｐｌｅｘ  Ｐｌａｎｔｓ”，ＩＥＥＥ  Ｃｏｎｔｒｏｌ  Ｓｙｓｔｅｍｓ  Ｍａｇａｚｉｎｅ，ｐｐ  ３−１１  （Ｄｅｃ．  １９８８）
（４）Ｇｅｒｔｌｅｒ，Ｊ．  ｅｔ  ａｌ．，‘Ｉｓｏｌａｔｉｏｎ  Ｅｎｈａｎｃｅｄ  Ｐｒｉｎｃｉｐａｌ  Ｃｏｍｐｏｎｅｎｔ  Ａｎａｌｙｓｉｓ’，ＡＩＣｈＥ  Ｊｏｕｒｎａｌ，Ｖｏｌ  ４５  Ｉｓｓｕｅ  ２，ｐｐ  ３２３−３３４，１９９９
（５）Ｇｒｏｓｓ，Ｋ．Ｃ．  ｅｔ  ａｌ，”Ａｐｐｌｉｃａｔｉｏｎ  ｏｆ  ａ  Ｍｏｄｅｌ−Ｂａｓｅｄ  Ｆａｕｌｔ  Ｄｅｔｅｃｔｉｏｎ  Ｓｙｓｔｅｍ  ｔｏ  Ｎｕｃｌｅａｒ  Ｐｌａｎｔ  Ｓｉｇｎａｌｓ”，Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｉｎｔｅｌｌｉｇｅｎｔ  Ｓｙｓｔｅｍ  Ａｐｐｌｉｃａｔｉｏｎ  ｔｏ  Ｐｏｗｅｒ  Ｓｙｓｔｅｍｓ，Ｊｕｌ．６−１０，１９９７，Ｓｅｏｕｌ，Ｋｏｒｅａ  ｐｐ．  ６６−７０
（６）Ｈｏｂｅｒｔ，Ｋ．Ｅ．  ＆  Ｕｐａｄｈｙａｙａ，Ｂ．Ｒ．  ”Ｅｍｐｉｒｉｃａｌ  Ｐｒｏｃｅｓｓ  Ｍｏｄｅｌｉｎｇ  Ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  Ｓｉｇｎａｌ  Ｖａｌｉｄａｔｉｏｎ”，Ａｎｎａｌｓ  ｏｆ  Ｎｕｃｌｅａｒ  Ｅｎｅｒｇｙ，Ｖｏｌ．２１，Ｎｏ．７，ｐｐ  ３８７−４０３，１９９４
（７）Ｌｏｎｇ，Ｔ；ｅｔ  ａｌ；”Ｓｅｎｓｏｒ  Ｆｕｓｉｏｎ  ａｎｄ  Ｆａｉｌｕｒｅ  Ｄｅｔｅｃｔｉｏｎ  Ｕｓｉｎｇ  Ｖｉｒｔｕａｌ  Ｓｅｎｓｏｒｓ”；Ｐｒｏｃｅｅｄｉｎｇｓ  ｏｆ  ｔｈｅ  １９９９  Ａｍｅｒｉｃａｎ  Ｃｏｎｔｒｏｌ  Ｃｏｎｆｅｒｅｎｃｅ；ｖｏｌ．４；Ｊｕｎ．１９９９；ｐｐ  ２４１７−２４２１
（８）Ｌｏｒｂｅｒ，Ａ，ｅｔ  ａｌ，”Ａ  Ｔｈｅｏｒｅｔｉｃａｌ  Ｆｏｕｎｄａｔｉｏｎ  ｆｏｒ  ｔｈｅ  ＰＬＳ  Ａｌｇｏｒｉｔｈｍ”，Ｊｏｕｒｎａｌ  ｏｆ  Ｃｈｅｍｏｍｅｔｒｉｃｓ，Ｖｏｌ．１，ｐｐ  １９−３１，１９８７
（９）Ｍａｎｕｓ，Ｈ．，”Ｖａｌｉｄａｔｉｎｇ  Ｄａｔａ  ｆｒｏｍ  Ｓｍａｒｔ  Ｓｅｎｓｏｒｓ”  Ｃｏｎｔｒｏｌ  Ｅｎｇｉｎｅｅｒｉｎｇ，ｐｐ．６３−６６，Ａｕｇ．１９９４
（１０）Ｍａｒｔｅｎｓ，Ｈ．，  ＆  Ｎａｅｓ，Ｔ．，”Ｍｕｌｔｉｖａｒｉａｔｅ  Ｃａｌｉｂｒａｔｉｏｎ”，Ｊｏｈｎ  Ｗｉｌｅｙ  ＆  Ｓｏｎｓ，１９８９
（１１）Ｐｉｏｖｏｓｏ，Ｍ．Ｊ．，  ｅｔ  ａｌ．  ”Ｐｒｏｃｅｓｓ  Ｄａｔａ  Ｃｈｅｍｏｍｅｔｒｉｃｓ”，ＩＥＥＥ  Ｔｒａｎｓ  ｏｎ  Ｉｎｓｔｒｕｍｅｎｔａｔｉｏｎ  ａｎｄ  Ｍｅａｓｕｒｅｍｅｎｔ，Ｖｏｌ．４１，Ｎｏ．２，  Ａｐｒｉｌ  １９９２，ｐｐ．２６２−２６８
（１２）Ｑｉｎ，Ｊ．Ｓ，ａｎｄ  Ｌｉ，Ｗ．，Ｄｅｔｅｃｔｉｏｎ，  ｉｄｅｎｔｉｆｉｃａｔｉｏｎ，ａｎｄ  ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ  ｏｆ  ｆａｕｌｔｙ  ｓｅｎｓｏｒｓ  ｗｉｔｈ  ｍａｘｉｍｉｚｅｄ  ｓｅｎｓｉｔｉｖｉｔｙ，ＡＩＣｈＥ  Ｊｏｕｒｎａｌ，Ｖｏｌ  ４５  Ｉｓｓｕｅ  ９，ｐｐ  １９６３−１９７６，１９９９
（１３）Ｓｏｗｉｚｒａｌ，Ｈ；”Ｖｉｒｔｕａｌ  Ｓｅｎｓｏｒｓ”；Ｐｒｏｃｅｅｄｉｎｇｓ  ＳＰＩＥ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｓｏｃｉｅｔｙ  ｆｏｒ  Ｏｐｔｉｃａｌ  Ｅｎｇｉｎｅｅｒｉｎｇ；ｖｏｌ．２４０９；Ｆｅｂ．７−９，１９９５；ｐｐ  ２４６−２５４
（１４）Ｓｐｏｅｌｄｅｒ，Ｈ；”Ｖｉｒｔｕａｌ  Ｉｎｓｔｒｕｍｅｎｔａｔｉｏｎ  ａｎｄ  Ｖｉｒｔｕａｌ  Ｅｎｖｉｒｏｎｍｅｎｔｓ”  ＩＥＥＥ  Ｉｎｓｔｒｕｍｅｎｔａｔｉｏｎ  ａｎｄ  Ｍｅａｓｕｒｅｍｅｎｔ  Ｍａｇａｚｉｎｅ；ｖｏｌ．２，Ｉｓｓｕｅ  ３；Ｓｅｐ．１９９９；ｐｐ  １４−１９
（１５）Ｔｓｏｕｋａｌａｓ，Ｌ．Ｓ．，  ＆  Ｕｈｒｉｇ，Ｒ．Ｅ．”Ｎｅｕｒａｌ  ａｎｄ  Ｆｕｚｚｙ  Ｓｙｓｔｅｍｓ  ｉｎ  Ｅｎｇｉｎｅｅｒｉｎｇ”，Ｊｏｈｎ  Ｗｉｌｅｙ，Ｎｅｗ  Ｙｏｒｋ，１９９６
（１６）Ｕｐａｄｈｙａｙａ，Ｂ．Ｒ．，ｅｔ  ａｌ  ”Ｍｕｌｔｉｖａｒｉａｔｅ  Ｓｔａｔｉｓｔｉｃａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｔｅｃｈｎｉｑｕｅｓ  ｆｏｒ  Ｆａｕｌｔ  Ｄｅｔｅｃｔｉｏｎ  ａｎｄ  Ｄｉａｇｎｏｓｉｓ”，ＩＳＡ  Ｔｒａｎｓａｃｔｉｏｎｓ，Ｖｏｌ．２９，Ｎｏ．４，ｐｐ  ７９−９５，１９９０
（１７）Ｖｅｎｋａｔａｓｕｂｒａｍａｎｉａｎ，Ｖ，ｅｔ  ａｌ  ”Ａ  Ｒｅｖｉｅｗ  ｏｆ  Ｐｒｏｃｅｓｓ  Ｆａｕｌｔ  Ｄｅｔｅｃｔｉｏｎ  ａｎｄ  Ｄｉａｇｎｏｓｉｓ”，Ｐａｒｔｓ  １，２，＆  ３，Ｃｏｍｐｕｔｅｒｓ  ａｎｄ  Ｃｈｅｍｉｃａｌ  Ｅｎｇｉｎｅｅｒｉｎｇ，ｖｏｌ  ２７，２００３
【図１】
2010519636.tif 000017
【図２】
2010519636.tif 000018
【図３】
2010519636.tif 000019
【図４】
2010519636.tif 000020
【図５】
2010519636.tif 000021
【図６】
2010519636.tif 000022
【図７】
2010519636.tif 000023
【図８】
2010519636.tif 000024
【図９】
2010519636.tif 000025
【図１０】
2010519636.tif 000026
【図１１】
2010519636.tif 000027
【図１２】
2010519636.tif 000028
【図１３】
2010519636.tif 000029
【国際調査報告】
2010519636.tif 000030