JP2002522837A

JP2002522837A - センサの性能確認装置および方法

Info

Publication number: JP2002522837A
Application number: JP2000565441A
Authority: JP
Inventors: クイン・エス・ジョー; ガイバー・ジョン
Original assignee: Aspen Technology Inc
Current assignee: Aspen Technology Inc
Priority date: 1998-08-17
Filing date: 1998-10-27
Publication date: 2002-07-23
Anticipated expiration: 2018-10-27
Also published as: JP4528335B2; EP1112532B1; JP2008198213A; DE69813040T2; EP1112532A1; CA2340977A1; DE69813040D1; JP4308437B2; WO2000010059A1; US6356857B1; CA2340977C; ATE236418T1

Abstract

(57)【要約】【課題】【解決手段】プロセスを測定するセンサに発生する故障を検出し、識別し、分類するための装置と方法が開示される。１つまたは複数のセンサにおいて故障を識別すると、この装置と方法は故障センサに対し代替値を提供して、これらのセンサを使用するすべてのプロセスコントローラおよびプロセスモニタリングシステムが故障期間中に運転を続行できるようにする。故障センサの識別は、センサの特定サブセットを感知しないと同時に、上記サブセット外のセンサに最高の感度をもつように特有に設計される、セットの構造化残差変換によって達成される。識別された故障は、完全な故障、バイアス、ドリフト、精度劣化、または不明のいずれかに分類される。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】

本発明は、一般にプロセス変数を測定するためにセンサを使用するプロセスコ
ントローラおよびプロセス監視システムに関するものである。特に本発明は、1
つまたは複数のセンサの故障を検出・識別し、センサの故障のタイプを分類し、
エラーのあるセンサ値を正しいプロセス変数値の推定値に代替するシステムに関
するものである。

【０００２】

【発明の背景】

センサの性能確認は、プロセス工業における多くのモデルをベースとする用途
に対する重要なステップである。典型的なモデルをベースとする用途には、モデ
ル予測制御用途（ＭＰＣ）および、実験室サンプルまたはハードウェアアナライ
ザから得られるコストがかかるか頻度の少ない測定を、モデルから定期的に得ら
れる推定値に代替する推定検知利用が含まれる。

【０００３】一般的なＭＰＣ用途では、定常状態最適化を実行して、制御また操作される変
数に対する最適目標値を見出す。センサが故障している場合は、最適化した目標
値は無効になる。したがって有効なセンサは、たとえ故障しても最良の再構成に
より制御システムをオンラインに維持する。

【０００４】代表的な推定センサ用途は、予測排出監視システム（ＰＥＭＳ）の分野である
。連邦および／または州の法規は、大気汚染プラントに対し窒素酸化物（ＮＯｘ
）、酸素（Ｏ₂）および一酸化炭素（ＣＯ）などの排出を監視することを要求す
ることができる。ハードウエア連続排出監視システム（ＣＥＭＳ）は、イニシャ
ルコストとメンテナンスコストの両者において高い負担が必要である。ＰＥＭＳ
に充分に精度と信頼性があることが明らかならば、ＣＥＭＳをＰＥＭＳに置き換
えることができる。ＰＥＭＳに対する品質保証要求事項の一つは、ＰＥＭＳモデ
ルに使用している各センサの故障の有無を監視すること、およびセンサの故障を
処理する方策をもち運転の中断時間を最小にすることである。

【０００５】センサの性能確認の用語は、本特許出願中では多変量モデルベースのセンサの
性能確認を指す。これは、プラントの測定に冗長度を使用する方法である。通常
センサの測定値は、収集されたまたは過去のデータを用いる訓練手順法により確
立できる相関構造を示す。この相関構造はオンラインで監視でき、相関構造が破
壊されたときは、センサ故障が既に発生している可能性がある。しかし、この相
関構造の破壊は、プロセスの反転、プロセスの過渡現象、またはセンサの故障に
関係のない他の理由によることもある。主たる目的は、これが真にセンサの故障
であるかどうかを判断し、故障であれば異常センサを識別することである。セン
サの性能確認の各種の段階は、以下のように要約できる。

【０００６】・検出この段階は、相関構造における変化を検出する。これは、センサ
の故障によるか、またはそうでないこともある・識別この段階は、これがセンサの故障であるかどうかを判断し、特定
のセンサを識別する・推定この段階は、真の値の再構成ができ、故障値の置換ができる故障
のサイズを推定する。・分類この段階は、センサの故障のタイプ―完全な故障、バイアス、ド
リフト、または精度劣化―を分類する。特定の方法によっては、これらの段階はオーバーラップすることもある。多変量
モデルをベースとするセンサの性能確認の主題を対象とするいくつかの特許が交
付されている。その主要なものは次の通りである：・Qinその他、米国特許第５，６８０，４０９号、“Method and Apparatus fo
r detecting and identifying faulty sensors in a process（プロセス中の故
障センサ検出および識別の方法と装置）” ・Keelerその他、米国特許第５，５４８，５２８号、“Virtual Continuous E
mission Monitoring System（仮想連続排出監視システム）” ・Hopkinsその他、米国特許第５，４４２，５６２号、“Method of controlli
ng a manufacturing process using multivariate analysis（多変量解析を利用
した製造プロセス制御の方法）”

【０００７】 Qinその他は、プロセス制御におけるセンサの性能確認を扱っている。好まし
い実施形態はＰＣＡ（主要コンポーネント分析）を基本にしており、最適の再構
成手順により識別を実施する。各センサ値は、それが故障しているとの仮定に基
づいて再構成され、次に識別および分類が再構成エラーから算出した追跡インデ
ックスにより行なわれる。

【０００８】 Keelerその他は、ＰＥＭＳにおけるセンサの性能確認を扱っている。開示され
たシステムは、推定センサテクノロジーおよびＰＥＭＳに対するニューラルネッ
トワークの使用に焦点を置く。センサの性能確認テクノロジーは、識別のための
準最適再構成手順を使用し、分類を扱わず、ＰＣＡの非線型バージョンである“
エンコーダ”ニューラルネットワークを使用する。エンコーダネットワークは、
Marker Kramer“Nonlinear principal component analysis using autoassociat
ive neural networks（自動結合ニューラルネットワークを使用する非線型主要
コンポーネント解析）”, AIChE Journal, 37(2), pp. 233-243(1991)にも記載
されている。

【０００９】 Hopkinsその他は、プロセス監視（多変量統計プロセス制御）記述中でセンサ
の性能確認を扱っており、ＰＣＡまたはＰＬＳ（部分最小自乗法）を使用する。
識別は寄与率分析により行なわれる。検出は、主要コンポーネント“得点”また
は得点統計を監視し、標準信頼間隔と比較することにより行なわれる。識別は、
異常得点に対する各オリジナル測定の寄与率を調べることにより行なわれる。こ
の方法は故障タイプを分類するものではない。

【００１０】

【発明の概要】

本発明は、故障センサの検出、識別、推定、再構成および分類のための新しい
装置と方法を提供する。このアプローチは、部分最小自乗（ＰＬＳ）または主要
コンポーネント分析のような統計法を用いて、第１原則またはデータから作るこ
とのできる標準プロセスモデルを使用する。好ましい実施形態においては、プロ
セスモデルは、主要コンポーネント数を選択して故障センサ値の再構成を最適化
しているＰＣＡモデルを基本にしており、これについてはQinおよびDuniaの“De
termining the number of principal components for best reconstruction（最
適再構成のための主要コンポーネント数の決定）”, Proc. Of the 5-th IFAC S
ymposium on Dynamics and Control of Process System, 359-364, Corfu, Gree
ce, June 8-10, 1998に記載されている。

【００１１】検出段階は、モデル等式エラーに基づく検出インデックスを使用する。指数関
数的に重み付けした移動平均（ＥＷＭＡ）フィルタを、検出インデックスに適用
して一時的な過渡現象による誤ったアラームを減少させる。フィルタ処理された
検出インデックス（ＦＤＩ）は、統計的に算出したしきい値と比較され、予想さ
れる故障を検出する。予想される故障状況の検出により、本発明の識別段階を起
動する。

【００１２】本発明の主要なコンポーネントは識別段階である。検出アラームが1つまたは
複数の故障センサに起因するかどうかを判定し、さらに異常センサ（複数もあり
）を識別するために、故障の一つのサブセットに影響されないが、他に対しては
最大感度をもつ一連の検出器を構成する。これらの検出器は、最大感度をもつ構
造化残差方法（ＳＲＡＭＳ）と称される新しい方法により構成された構造化残差
（ＳＲ）を基本にする。構造化残差は、一般にGertlerとSinger,“A new struct
ural framework for parity equation based failure detection and isolation
（故障検出と分離を基にしたパリティ等式に対する新しい構造フレームワーク）
”, Automatica 26:381-388, 1990に記載されている。指数関数的に重み付けさ
れた移動平均（ＥＷＭＡ）フィルタをＳＲに適用し、一時的な過渡現象による誤
ったアラームを減少させる。またＳＲも自乗および正規化して、各ＳＲを公平に
比較する。識別は、これらの正規化した自乗フィルタ処理構造化残差（ＮＳＦＳ
Ｒ）を、統計的に推定した信頼限界と比較することにより行なわれる。ＮＳＦＳ
Ｒに加え、検出の時点から累積加算した正規化構造化残差（ＮＳＲ）に基づくイ
ンデックスを監視し、比較して、故障センサの識別に使用する。この２つの指数
が普遍化された公算比（ＧＬＲ）インデックスであり、また正規化した累積分散
（ＮＣＵＭＶＡＲ）インデックスある。ＮＣＵＭＶＡＲインデックスは、精度の
劣化したセンサを識別するために主として用いることができる。

【００１３】この時、故障の大きさは、モデル、故障データ、および故障センサが正しく識
別されたとの仮定に基づいて最適の形で推定される。これは、例えばMartensお
よびNaes “Multivariate Calibration”, John Wiley and Sons, New York, 19
89に記載された公開されている既存の技術を使用する。故障方向の知識（故障セ
ンサの識別から知る）および推定される故障の大きさを使用し、正しいセンサ値
を再構成する。

【００１４】故障分類段階は、センサ故障のタイプに関する診断情報を提供する。特に４つ
のタイプの故障が考えられる、つまり、完全な故障、バイアス、ドリフトおよび
精度の劣化である。完全な故障は、識別した故障センサ測定値に関する回帰分析
を行うことにより求められ、回帰線がゼロの傾きをもつことを統計的に推定して
表わされる。他の３つの故障タイプを分類するには、識別時点以降の推定される
故障サイズに関する回帰分析を行う。バイアスは、推定される故障サイズ回帰線
がゼロの傾きおよびゼロでないオフセットをもつとする統計推定により示され、
わずかな残差エラーをもつ。ドリフトは、推定される故障サイズ回帰線がゼロで
ない傾きをもつことを統計的に推定して示され、わずかな残差エラーをもつ。精
度劣化は、推定された故障サイズ回帰線がゼロの傾き、ゼロのオフセットおよび
有意な残差エラーをもつことを統計的に推定して示される。また精度の劣化は、
故障がＮＣＵＭＶＡＲインデックスによってのみ識別可能な場合にも、示される
。

【００１５】

【発明の実施の形態】

本発明の上記およびその他の目的、特徴および利点は、添付の図面に図示され
た本発明の好ましい実施形態のさらに詳しい記述で明らかとなるであろう。図中
では一貫して同一部分に対しては同一符号を用いている。図面は必ずしも縮尺通
りでなく、本発明の原理を示するために強調が加えられている。

【００１６】図１は、本発明のセンサの性能確認装置を示している。センサ確認装置の外部
入力はセンサ測定値１０８から成り、その測定値には、例えば製造プラントにお
ける分散制御システム（ＤＣＳ）５０から得られる状態情報を含む。本発明の目
的は、センサ測定値の冗長情報を用いて、センサがいつ故障したかを検出し、故
障したセンサに選択的に代替値を提供することである。センサ測定値を使用する
プロセスコントローラ６０または推定センサ７０のような他の装置は、この時必
要に応じ代替値を使用できる。これにより、これらの装置は、故障したセンサの
存在下で引き続き動作を続行できる。

【００１７】再び図１によれば、センサの性能確認装置１００は、センサ確認の各種の段階
を実行する８つの演算ユニットから構成される。マイクロプロセッサ、または他
の所望のタイプのコンピュータを用いて実行できる演算ユニットは、次の通りで
ある。

【００１８】・前処理ユニット１１０、これは拡大縮小とオフセット情報１１２を維持して
、ＤＣＳ５０からのセンサ値１０８のベクトルを正規化した値１１８のベクトル
に変換する。これらの正規化した値は、標準化ユニットをもち、好ましい実施形
態においては、そのユニットに対し、その値を平均値０と標準偏差１をもつよう
に設計する。これは、故障センサを識別するために、有効なモデル行列１２２、
検出インデックス１４２および構造化残差変換１３２を形成するのに必要である
。拡大縮小およびオフセット情報は、好ましい実施形態においてはセンサの過去
のデータセットの統計分析から求められる。

【００１９】・モデルユニット１２０、これはモデル行列１２２を維持して、所定の時点で
前処理されたセンサ値１１８のベクトルを等式エラー１２８と称されるエラーベ
クトルに変換する。等式エラーは、構造化残差ユニット１３０で使用して構造化
残差値１３８および故障サイズ１３６を算定し、また検出ユニット１４０で使用
して検出インデックス１４２を計算する。また構造化残差ユニットもモデル行列
を使用し、構造化残差変換１３２およびその対応する故障推定変換１３４を設計
する。好ましい実施形態においては、モデル行列は、主要コンポーネント分析を
使用してセンサの過去のデータセットから形成される。

【００２０】・構造化残差ユニット１３０、これは構造化残差変換１３２の使用セットを維
持し、現在の等式エラー１２８を１セットの構造化残差値１３８に変換する。識
別ユニット１５０はこれらの構造化残差値を用いて、故障センサの同一性を判定
する。構造化残差ユニットもまた推定ユニット１６０により使用される１セット
の故障推定変換１３４を維持し、識別した故障に対する正規化した故障サイズ１
３６を得る。

【００２１】・検出ユニット１４０、これは等式エラー１２８から検出インデックス１４２
を計算し、検出インデックスを検出しきい値１４４と比較する。検出インデック
スが検出しきい値を上回る場合は、検出イベント１４８を発生する。この検出イ
ベントを使用して識別ユニット１５０を起動し、また、例えばオペレータアラー
ム発生、またはオペレータコンソール９０もしくは格納装置９２へのイベント記
録に、センサ確認装置の外部で利用できる。

【００２２】・識別ユニット１５０、これは構造化残差値１３８を使用して、識別インデッ
クスのいくつかのセット１５２を計算する。検出イベント１４８が発生した場合
は、識別ユニットは、後述の特定のロジックに従って、識別インデックスをしき
い値１５４の対応するセットと比較して故障特性を識別する。識別に成功すると
、識別イベント１５８を発生し、それを使用して推定ユニットを起動する。この
識別イベントは、どのセンサが故障しているかに関する情報を含み、例えばオペ
レータアラーム発生、またはオペレータコンソール９０もしくは格納装置９２へ
のイベント記録に、センサ確認装置の外部で利用できる。

【００２３】・推定ユニット１６０、これは識別イベント１５８により起動される。推定ユ
ニットは、構造化残差ユニット１３０に、等式エラー１２８に適正な故障推定変
換１３４（識別情報に基づいて選択される）を適用するよう要求し、その結果正
規化した故障サイズ情報１３６を生成する。次に推定ユニットは、適切な拡大縮
小１１６を１３６に適用して、オリジナルセンサ値の単位に対し拡大縮小された
故障サイズ推定１６８を生成する。

【００２４】・代替ユニット１７０、これは、対応する故障センサ測定値１０８から推定さ
れた故障サイズ１６８を差し引き、１セットの代替値を得る。これらの代替値を
センサ確認装置の外部で利用できる。例えばそれらをＤＣＳ５０に書込み、他の
装置で利用する。代替ユニットは、推定ユニット１６０が起動しているときだけ
機能する。代替ユニットは、オペレータ１７２の要求によっても機能を停止する
ことができる。この場合には、代替ユニットの出力１７８はセンサ確認装置に対
するオリジナル入力１０８と同一である。

【００２５】・分類ユニット１８０、これはセンサ故障のタイプに関する分析情報を計算す
る。特に４種の故障−完全な故障、バイアス、ドリフトおよび精度劣化−が考え
られる。分類ユニットは、推定ユニット１６０が起動している時だけ機能し、ま
た識別イベント１５８が発生したときだけ機能する。分類は、識別時点以降のセ
ンサデータ１０８の分析および故障サイズ推定１６８に基づく。分類ユニットに
より生成される分類情報１８８は、センサ確認装置の外部で利用できる。例えば
オペレータアラーム発生、またはオペレータコンソール９０もしくは格納装置９
２へのイベント記録に使用できる。これらの演算ユニットのそれぞれの好ましい実施形態を以下に詳しく述べる。

【００２６】１．前処理ユニット前処理ユニット１１０は、拡大縮小およびオフセット情報を計算および維持し
て、ＤＣＳからの測定されたセンサ値１０８と正規化された値１１８との間で変
換を行い、それをモデルユニットおよび構造化残差ユニットに入力する。ｉ番目
のセンサに対する拡大縮小値Ｓ_iおよびオフセット値Ｏ_iは、好ましい実施形態に
おいては、ｉ番目のセンサの平均値μ_iおよび標準偏差σ_iから下記のように計算
される。

【００２７】

【数１】

【００２８】次にｉ番目のセンサの測定された値ｕ_iは、次の式により標準化された値ｘ_iに
変換される。

【００２９】

【数２】

【００３０】 μ_iおよびσ_iは、代表的なサンプルデータから、当業者にとって公知の標準法
により推定される。

【００３１】２．モデルユニットモデルユニット１２０は、下記の式で表わされる標準プロセスモデルを維持す
る。

【００３２】

【数３】

【００３３】ここでｘ^*∈Ｒⁿは、前処理されたセンサ値のベクトルであり、Ｂ∈Ｒ^mxnはモ
デル行列１２２であり、ｅ^*∈Ｒ^mは等式エラー１２８である。行列Ｂは、プロセ
スの質量バランスもしくはエネルギーバランスから算出するか、または主要コン
ポーネント解析（ＰＣＡ）もしくは部分最小自乗（ＰＬＳ）のような統計法を用
いてプロセスデータから算出できる。

【００３４】２．１ＰＣＡに基づくモデルユニット好ましい実施形態においては、標準プロセスモデルはＰＣＡを用いてプロセス
データから求められる。ＰＣＡは、一般に例えばJ. Edward Jackson“A User's
Guide to Principal Components（主要コンポーネントに対するユーザーズガイ
ド）”、John Wiley 1991に記載されている。ＰＣＡにおいては、センサ値ｘ^*(t
)∈Ｒⁿは、次の形に分解される。

【００３５】

【数４】

【００３６】ここでＰ∈Ｒ^nx(n-m)は、ｘ^*(t)の相関行列の主たる固有値λ₁≧λ₂≧Λλ_n-m に関連する固有ベクトルである。この場合には、モデル等式エラーｅ^*(t)は、残
差コンポーネント（λ_n-m+1, Λλ_n）に対応し、またGertler, Li, HuangとMcAv
oy“Isolation enhanced principal component analysis（分離拡張主要コンポ
ーネント解析）”Proc. IFAC Workshop on On-line Fault Detection and Super
vision in the Chemical Process Industry, 4-5 June, 1998に記載されている
ように次の式が成り立つ。

【００３７】

【数５】

【００３８】好ましい実施形態においては、最適数の主要コンポーネントを選択して最適再
構成機能を達成する、これについてはQinおよびDuniaの“Determining the numb
er of principal components for best reconstruction（最適再構成のための主
要コンポーネント数の決定）”、Proc. Of the 5-th IFAC Symposium on Dynami
cs and Control of Process System, 359-364, Corfu, Greece, June 8-10, 199
8に詳述されている。

【００３９】３．構造化残差ユニット構造化残差ユニット１３０は、特定方向の故障を感知するかまたは感知しない
変換を維持する。これらの変換は、以下の故障表現の中で説明する。

【００４０】３．１故障の表現正常な条件下では、等式エラーｅ^*(t)１２８は主とし測定雑音を含み、この雑
音はゼロの平均ガウス雑音と仮定できる。センサ故障が発生すると、センサ測定
はプロセス変数の正常値と故障値を含み、次式のようになる。

【００４１】

【数６】

【００４２】ここでｆ_i(t)∈Ｒ^liは故障の大きさのベクトルであり、Ξ_i∈Ｒ^nxliは故障方
向の行列である。ｉ番目のセンサの単一センサ故障を表わすと次式になる

【００４３】

【数７】

【００４４】上記の式は、単位行列のｉ番目の列である。複数センサ故障を同時に表わすに
は、Ξ_iは単位行列の対応する列を単に含むだけでよい。センサ故障を含む式６
の測定値関係を用いると、等式エラー１２８は下記のように表わすことができる
。

【００４５】

【数８】

【００４６】故障は、一般に等式エラーｅ(t)を増大させる。等式エラーの大きさは、後述
のようなセンサ故障を検出するのに用いられる。個別のベクトルまたは行列Ξ_i
が各センサ故障を表わすため、故障の識別は方向ベクトルを用いることにより達
成できる。

【００４７】３．２単一故障に対する構造化残差変換構造化残差ユニット１３０は、構造化残差と呼ばれる１セットの変換を維持す
る。後述する識別ユニット１５０は、故障センサを識別するために構造化残差変
換を用いる。

【００４８】構造化残差は、一般にGertlerとSinger、“A new structural framework for
parity equation based failure detection and isolation（故障検出と分離を
基にしたパリティ等式に対する新しい構造フレームワーク）”、Automatica 26:
381-388, 1990に記載されている。本発明においては、構造化残差に対する新し
い設計法を開示しており、“Structured Residual Approach with Maximized Se
nsitivity（最大感度をもつ構造化残差アプローチ）”またはＳＲＡＭＳと呼ば
れる。ｊ番目のセンサにおける単一センサ故障の場合に対しては、次の式が成り立つ
。

【００４９】

【数９】

【００５０】ここでｂ_jは行列Ｂのｊ列である。したがってモデル等式エラーは次式となる
。

【００５１】

【数１０】

【００５２】選択された残差構造を用いると、１セットの構造化残差ｒ(t)は、等式エラー
ベクトルｅ(t)に変換行列Wをあらかじめ乗じることにより生成できる。

【００５３】

【数１１】

【００５４】ｒ(t)の各エレメントを一つの特定のセンサ故障に対して関与するが、他の故
障に対して関与しないように設計することが望ましい。ｊ番目センサが故障して
いるｒ(t)のｉ番目のエレメントは、次式で表わされる。

【００５５】

【数１２】

【００５６】ここでｗ_i ^Tは行列Ｗ∈Ｒ^nxmのｉ行である。ｅ(t)がランダム雑音であるため、残
差ｒ_ij(t)を故障方向ベクトルｂ_jに対して感知するかまたはしないようにだけ設
計できる。ＳＲＡＭＳデザイン基準は次のように記述される。ｗ_iをｒ_ij(t)がｉ番目センサ故障を感知しないが、他のものは最大感度で感知
するように選ぶ。数学的には次の形と同一である。

【００５７】

【数１３】

【００５８】ただし、

【００５９】

【数１４】

【００６０】幾何学的には、ｗ_iはｂ_iに対して直角であるように選択し、他の故障方向ｂ_j
ｊ≠ｉに対するその角度は最小に選択する。この問題は、次のようにも表現する
ことができる。

【００６１】

【数１５】

【００６２】ただし、

【００６３】

【数１６】

【００６４】が条件である。ここで、

【００６５】

【数１７】

【００６６】制約ｗ_i ^Tｂ_i ^o＝０を満足するために、ベクトルｗ_iは下記のように選択する必
要がある。

【００６７】

【数１８】

【００６８】ここでＳ_wiはｂ_iの９０°直交補空間である。したがってｗ_iの解はｚ_iの解に
変換される。ラグランジュ乗数λを用いて、ｚ_iの解は下記の値を最大とするこ
とにより見出せる。

【００６９】

【数１９】

【００７０】（Ｉ−ｂ_i ^oｂ_i ^oT）ｂ_i ^o＝０であるから上の式は成り立つ。次のように表わす
。

【００７１】

【数２０】

【００７２】および

【００７３】

【数２１】

【００７４】上式はＢ^oのＳ_wiへの投影であり、最大化の問題は次のように書き換えできる
。

【００７５】

【数２２】

【００７６】Ｊをｚ_iに関して微分すると次の値が得られる。

【００７７】

【数２３】

【００７８】ただし、次の関係に注意が必要である

【００７９】

【数２４】

【００８０】式２０は次のように変形できる

【００８１】

【数２５】

【００８２】または

【００８３】

【数２６】

【００８４】したがってｚ_ijは、最大の固有ベクトルに相当するＢ_i ^oｂ_i ^oTの固有ベクトル
である。ｚ_iが判明すると、ｗ_iは式１５から容易に計算できる。

【００８５】３．３複数故障に対する構造化残差変換ＳＲＡＭＳ設計法は、複数センサが同時に故障するケースに応用することがで
きる。複数センサ故障に対する等式エラーベクトルは、次の式で表わされる。

【００８６】

【数２７】

【００８７】ここでＢΞ_i∈Ｒ^mxqは、故障したセンサに対応するBのｑ個の列を含み、また
ｆ_i(t) ∈Ｒ^qは故障の大きさのベクトルである。ＢΞ_iの列は線型関係をもつこ
とに留意が必要である。ＢΞ_iに単一値分解（当業者によく知られている行列分
解）を実施し、ゼロでない単独値のみを維持することにより次の式が得られる

【００８８】

【数２８】

【００８９】ここで、

【００９０】

【数２９】

【００９１】複数センサ故障に対してＷ_qを設計するために、Ｗ_qはＵ_qの直交補空間内に選
択する。即ち、

【００９２】

【数３０】

【００９３】単一センサ故障ケースに類似の次式となる。さらに、次式が最大とされる。

【００９４】

【数３１】

【００９５】ただし、

【００９６】

【数３２】

【００９７】ここで、

【００９８】

【数３３】

【００９９】ラグランジュ乗数を用いた後、ｚ_qは次式の解であることが明らかになる。

【０１００】

【数３４】

【０１０１】式２９から次の関係式が得られる。

【０１０２】

【数３５】

【０１０３】または

【０１０４】

【数３６】

【０１０５】したがって、ｚ_qは、単に、最大固有値に関連するＢ_q ^oＢ_q ^oTの固有ベクトルで
ある。

【０１０６】３．４従来の構造化残差方法との比較従来の構造化残差方法（一般にGertlerとSinger、“A new structural framew
ork for parity equation based failure detection and isolation（故障検出
と分離を基にしたパリティ等式に対する新しい構造フレームワーク）” , Autom
atica 26:381-388, 1990、およびGerterとSinger、“Augmented models for sta
tistical fault isolation in complex dynamic systems（複雑動的システムに
おける統計的故障分離に対する拡張モデル）” , Proceeding of the American
Control Conference, pp317-322, 1985）は、ｗ_iを対象とする特定の故障に対し
て感知しないものを選択するが、他の故障に対する感度を最大にすることはない
。代表的な構造化残差設計では、ｗ_iの選択は、特有のものではなく、ある程度
任意である。この設計おける任意性により、故障を分離する可能性を最大にしな
い準最適解を導く。本発明のＳＲＡＭＳ方法は、ｗ_iの特有の設計を実現して、
他の故障に対する感度を最大にする反面、対象の故障を感知しない。

【０１０７】３．５構造化残差ユニットにより維持される変換本発明の好ましい実施形態において、構造化残差ユニットは、先のセクション
３．２に記載の単一センサ故障に対応する構造化残差変換の中核セットを維持す
る。これらの中核構造化残差変換は、構造化残差変換１３２の初期動作セットお
よび関連の故障推定変換１３４を形成する。

【０１０８】ＤＳＣからの１つまたは複数のセンサ測定値１０８が、センサ確認装置への入
力において不良と表示されるときは、構造化残差変換および関連の故障推定変換
の新しい動作セットが、構造化残差ユニットにより形成される。これらは次の通
りである。ｂ個のセンサが入力において不良と表示されると、ｎ−ｂ＋１の新し
い構造化残差変換が設計される。第１構造化残差変換が、不良のｂ個のセンサに
対応する恒式行列の列を選択するΞを用いて設計される。この構造化残差変換は
、これらの不良センサ値を感知しないが、他のセンサのいずれの故障も感知する
。残りのｎ−ｂの構造化残差変換は、不良のｂ個のセンサに対応する恒式行列の
列に加え、相互のセンサに対応するさらにもう１つの列を選択するΞにより設計
される。これにより入力で不良と表示されるセンサの外にも、別の故障センサを
識別することが可能になる。

【０１０９】その他の構造化残差変換は、選択的に設計して、動作セットに加えることがで
きる。例えばΞに基づく１セットの変換を、恒式行列からのすべてのペアの列か
ら選択する。これにより、同時に故障しているペアを識別できる。しかし、複数
故障を検出する好ましい方法は、いくつかのセンサ確認装置を用いることである
。例えば図１において、プロセスコントローラ６０が入力として５つのセンサ値
を使用する場合は、最高５つのセンサ確認装置１００を使用できる。各センサ確
認装置は、入力の１つおよびプロセスコントローラに入力されない他の冗長セン
サ測定値を使用する。

【０１１０】別の方法は、いくつかのセンサ確認装置をもつものであり、そのそれぞれは特
定のセンサの故障の検出および識別に最適化され、再構成のために最適化される
別のセンサ確認装置に信号供給する。

【０１１１】４．検出ユニット検出ユニット１４０は、現在の等式エラーに基づいて検出インデックス１４２
を計算することにより、前処理されたセンサ値の相関構造を監視する。この検出
インデックスが、統計的に推定された信頼限界に基づくしきい値１４４を越える
場合には、検出イベント１４８が発生し、識別ユニットが起動される。好ましい
実施形態では、正規化され、フィルタ処理された故障検出インデックス（ＮＦＦ
ＤＩ）は、以下のように使用される。

【０１１２】４．１故障検出インデックス故障検出インデックス（ＦＤＩ）は、下記のように定義される。

【０１１３】

【数３７】

【０１１４】故障が存在しないときは、ＦＤＩ＝‖ｅ^*(t)‖²＝‖Ｂｘ^*(t)‖²である。した
がって、ＦＤＩに対する検出アラームしきい値は、実際には‖ｅ^*(t)‖²に対す
る信頼限界である。

【０１１５】４．２正規化故障検出インデックス実際の工業プロセスにおいて、プロセス変数ｘ^*(t)は、一般には正規分布して
いない。しかし、等式エラーｅ^*(t)１２８は、主として測定雑音を含み、したが
ってゼロの平均ガウス雑音と考えることができる。したがって、ｅ^*(t)^TＲ_e ^-1ｅ ^* (t)は、ｍの自由度をもつχ²分布を満足する、ここで、Ｒ_e＝Ｅ[ｅ(t)ｅ^T(t)]
は等式エラー共分散行列である。このように正規化された故障検出インデックス
を、

【０１１６】

【数３８】

【０１１７】と定義することにより、ＮＦＤＩに対するｄ_aの検出アラームしきい値は、当業
者には既知の方法で、χ²分布を用い任意の有意水準αに対し求めることができ
る。

【０１１８】ＮＦＤＩ≡‖η(t)‖²、ここでη(t)＝Ｒ_e ^-1/2ｅ(t)およびＲ_e ^-1/2は逆等式エ
ラー共分散行列の行列平方根であることに留意すること。Ｒ_e ^-1/2は等式エラー
の相関性解除の効果をもつ。

【０１１９】４．３正規化され、フィルタ処理された故障検出インデックス好ましい実施形態においては、指数関数的に重み付けされた移動平均（ＥＷＭ
Ａ）フィルタを等式エラーに適用して、測定データ中の過渡現象および雑音の効
果を減少させる。

【０１２０】

【数３９】

【０１２１】したがって次式が成り立つ。

【０１２２】

【数４０】

【０１２３】および

【０１２４】

【数４１】

【０１２５】ここで、

【０１２６】

【数４２】

【０１２７】正規化され、フィルタ処理された故障検出インデックス１４２は、次のように
定義される。

【０１２８】

【数４３】

【０１２９】故障がなければ、η(t)は正規分布するため、^-η(t)はまた正規分布し、Ｅ[^-
η(t)]＝Ｅ[(t)η]＝０である。しかし、ＥＷＭＡフィルタは、^-η(t)の分散を
大幅に減少させるため、信頼限界ｄ_aを直接 ^-ｄ_aに適用することはできない。^-
η(t)の分散は次の式により表わされる。

【０１３０】

【数４４】

【０１３１】ここで、

【０１３２】

【数４５】

【０１３３】は自己相関関数である。η(t)が独立的かつ同一分布している場合は、ρ(k)は任
意のｋ＞０に対して０である。しかし、実際のデータに対しては、動的な過渡現
象のためにこの条件を満足することは不可能である。式３６の右辺のゼロでない
ρ(k)は、 ^-ｄ_aに対する検出アラームしきい値の決定に寄与する。‖η(t)‖²に
対するしきい値がｄ_aであるから、式３６から、‖^-η(t)‖²に対する検出アラー
ムしきい値１４４が下記の式であることは、容易に知ることができる：

【０１３４】

【数４６】

【０１３５】 │ρ(k)│≦１、したがって次式が成り立つことに留意すること。

【０１３６】

【数４７】

【０１３７】したがって、ＮＦＦＤＩに対するしきい値は、通常ＮＦＤＩに対するしきい値
よりも小さい。実際に、ｐ(k)は、一定の時間ｋ後にはゼロに減少し、検出ユニ
ットだけが固定数の項を維持して ^-ｄ_aを計算する。有意性αは、検出ユニット
の調節可能なパラメータである。図４はＮＦＦＤＩＮが故障発生を急速に増加している様子を示す。

【０１３８】５．識別ユニット識別ユニット１５０は、検出ユニットが検出イベントを発生する場合に、起動
される。識別ユニットは、インデックス１５２のいくつかのセットを計算し、そ
れを用いて、検出された故障を１つまたは複数のセンサの故障に識別できるかど
うかを判定する。各セットのインデックスは、構造化残差ユニット１３０により
維持される有効セットの構造化残差から算出される。一般に１つのセットにおけ
る単一インデックスがしきい値よりも小さく、一方でそのセットの中のすべての
他のインデックスがしきい値よりも大きい場合には、既知の方向の故障が識別さ
れる。この故障の方向は、インデックスに関連する構造化残差を定義する方向に
一致する。最も一般的な場合は、この故障は単一センサの故障に相当する。

【０１３９】本発明の好ましい実施形態において、識別ユニットは、３セットのインデック
ス１５２およびそれらの関連しきい値１５４を維持する。これらのインデックス
のセットは、正規化され、自乗され、フィルタ処理された構造化残差（ＮＳＦＳ
Ｒ）、一般化公算比（ＧＬＲ）インデックス、および正規化累積加算分散（ＮＣ
ＵＭＶＡＲ）インデックスである。関連しきい値は、ｇ_NSFSR、ｇ_GLRおよびｇ_NC _UMVAR により表わされる。図２に、故障識別プロセスを要約したフローダイアグ
ラムを示す。

【０１４０】５．１正規化され、自乗され、フィルタ処理された構造化残差理論的には、各構造化残差ｒ_i(t)は、センサが故障しないときはゼロであるか
、またはゼロに近いと想定される。しかし、モデル化エラー、測定値ノイズおよ
び他の不確定性のために、ｒ_i(t)はゼロに等しくない。したがって、各ｒ_i(t)ｉ
,ｊ＝１，２，Λ，ｎに対するゼロでない信頼限界は、統計技法を用いて求めら
れる。

【０１４１】故障がないか、またはｉ番目の構造化残差が故障のセンサ（複数もある）では
、ｉ番目の構造化残差は正規分布に従う。

【０１４２】

【数４８】

【０１４３】ここでＲ_e＝Ｅ｛ｅ(t)ｅ^T(t)｝は等式エラー共分散行列である。したがって、

【０１４４】

【数４９】

【０１４５】次に、ｉ番目の正規化構造化残差（ＮＳＲ）を、

【０１４６】

【数５０】

【０１４７】として定義し、正規化され自乗された構造化残差（ＮＳＳＲ）をｖ_i ²(t)と定義
することにより、ＮＳＳＲに対する信頼限界はχ_a ²(1)になる、ここで、αは有
意水準である。

【０１４８】ＥＷＭＡフィルタをＮＳＲに適用すると、正規化されフィルタ処理された構造
化残差（ＮＦＳＲ）は、次式で表わされる。

【０１４９】

【数５１】

【０１５０】故障がないか、またはｉ番目の構造化残差が故障のセンサ（複数もある）では
、ＮＦＳＲも次式による正規分布に従う。

【０１５１】

【数５２】

【０１５２】および

【０１５３】

【数５３】

【０１５４】ここで

【０１５５】

【数５４】

【０１５６】は自己相関関数である。

【０１５７】したがって、 ^-ｖ_i ²(t)により定義された、正規化され自乗されフィルタ処理
された構造化残差（ＮＳＦＳＲ）、に対するアラームしきい値は、次の式により
示されるアラームしきい値をもつ。

【０１５８】

【数５５】

【０１５９】すべてのセンサが正常であれば、 ^-ｖ_i ²(t)は、次の式により示されるアラー
ムしきい値を通常下回る。

【０１６０】

【数５６】

【０１６１】ｉ番目の構造化残差に関連するセンサ（複数もあり）がＳＲＡＭＳ設計により
故障の場合には、^-ｖ_i ²(t)は故障の影響を受けることはない。しかし他の残差^-
ｖ_i ²(t)ｊ≠ｉは大幅に増大する、なぜなら、それらの他のセンサの故障に対す
る感度は最大であるからである。

【０１６２】したがって、選択された有意水準αおよび時点ｔを用いると、ＮＳＦＳＲsの
中で、そのしきい値ｇⁱ _NSFSRより小さい正確な1つのＮＳＦＳＲが存在すること
により、故障センサを識別し、識別イベント１５８を発生することになる。この
例は図７に示され、その中で上のウインドウは故障センサ５に対応しないＮＳＦ
ＳＲは極めて大きくなるのに対し、下のウインドウはセンサ５に対するＮＳＦＳ
Ｒはそのしきい値を下回ったままであることを示す。

【０１６３】多数のセンサの性能を確認することが多いため、たとえ各インデックスがその
しきい値を上回る確率が低くても、最低限１つの残差がそのしきい値を上回る確
率をもつことは有意である。例えば２０の残差があると仮定すると、ｒ_i(t)ｉ＝
１，２，Λ，２０である。少なくとも一つのＮＳＦＳＲがそのしきい値を上回る
確率は（α＝０．０５と仮定して）、すべてのセンサが正常であってもＰ＝１−(１−α)²⁰＝１−０．９５²⁰＝０．６４２である。したがって、ＮＳＦＳＲは、故障検出に対して使用されない。代わりに
、検出ユニットは、ＮＦＦＤＩを用いて検出イベントをトリガーし、その後に識
別ユニットが、各種セットのインデックスを使用して、ＮＳＦＳＲから始めて、
故障センサを識別する。

【０１６４】検出ユニットが検出イベントを発生しても、いなくても、選択的にＮＳＦＳＲ
インデックスを各サイクルにおいて計算できる。しかし、次の２つのセクション
に記載したインデックスセットは、故障検出の時点以降、構造化残差情報を累積
加算する。したがって、それらは検出イベントの発生以降においてのみ有効であ
る。

【０１６５】５．２一般化された公算比インデックスセンサ故障が平均して有意な変化、例えば完全故障、バイアス、ドリフトを生
じる場合は、一般化された公算比（ＧＬＲ）は、通常これを検出するのに適正な
検定法である（Benveniste, BassevilleおよびMoustakides“The asymptotic lo
cal approach to change detection and model validation（検出とモデル確認
変化に対する漸近軌跡方法）”、IEEE Trans. Auto. Cont. 32(7):538-592, Jul
y 1987）。

【０１６６】以下の分析を明瞭に示すために、各構造化残差変換は、単一センサにおける故
障に相当すると仮定する。しかし、一般には、ＧＬＲインデックスは、構造化残
差設計には関係なく構成し、監視することができる。

【０１６７】式４２に示すように、センサが故障でなければ、ｖ_i(t)はゼロの平均であり、
単位分散をもつ正規分布を示す。しかしｊ番目のセンサがｔ_fの時点で故障を発
生すると、次式

【０１６８】

【数５７】

【０１６９】となり、ゼロの平均ではなく、分散は変化しないことが明らかである。言い換え
れば次式が成り立つ。

【０１７０】

【数５８】

【０１７１】ここで

【０１７２】

【数５９】

【０１７３】ＳＲＡＭＳ設計より、μ_jj＝０である。

【０１７４】ｖ_i(t)の平均変化を検出するために、下記のＧＬＲ関数を定義する。

【０１７５】

【数６０】

【０１７６】 μ_ijに関してＳ_tf ^t（μ_ij）を微分することにより、μ_ijの最適の推定値が求
められる。

【０１７７】

【数６１】

【０１７８】 ^<μ_ijを用いると、Ｓ_tf ^t（μ_ij）は最大となり次式が成り立つ。

【０１７９】

【数６２】

【０１８０】ｊ番目のセンサが故障であれば、μ_jj＝０、およびｖ_j(t)〜Ｘ(０，１)である
。式５０から次式が成り立つ。

【０１８１】

【数６３】

【０１８２】または

【０１８３】

【数６４】

【０１８４】したがって、選択された有意水準αおよび時点ｔを用いると、次のｎＧＬＲs
の中で、

【０１８５】

【数６５】

【０１８６】しきい値ｇ_GLR＝χ_a ²(１)よりも小さい確実に１つのＧＬＲが存在することによ
り、故障したセンサを特定し、識別イベント１５８を発生する。

【０１８７】５．３正規化された累積加算分散インデックスセンサが精度劣化故障を生じた場合、ＧＬＲインデックスセットはそれを識別
するのが困難である。精度劣化故障の位置を識別するため、ＮＳＲの平均および
分散が、故障の生じた時点以降計算される。

【０１８８】

【数６６】

【０１８９】説明を判りやすくするために、各構造化残差変換が単一センサの故障に相当す
ると仮定すれば、故障のセンサがないか、またはｉ番目のセンサが故障のときは
、次式が成り立つ。

【０１９０】

【数６７】

【０１９１】したがって、選択された有意水準αおよび時点ｔを用いると、ｎ個のＮＣＵＭ
ＶＡＲs間で、しきい値ｇ_NCUMVAR＝[１／(ｔ−ｔ_f) ]χ_a ²（ｔ−ｔ_f−１）よりも小さい確実に1つのＮＣＵＭＶＡＲが存在すると、精度故障をもつセ
ンサを識別定し、識別イベント１５８を発生する。

【０１９２】６．推定ユニット推定ユニット１６０は、識別ユニット１５０が識別イベント１５８を発生した
場合だけ、起動される。この時、推定ユニットは、識別イベントを発生するのに
機能した構造化残差変換に対応するモデルに基づく故障の大きさ、故障データお
よび行列Ξ_iを最適の形で推定する。この最適推定は、例えばMartensとNaes “M
ultivariate Calibration（多変量キャリブレーション）”、John Wiley and So
ns, New York, 1989に記載された公開されている既知の方法に基づくものである
。

【０１９３】最適推定は、方向Ｕ_qでの時間関数ｆ_q(t)を選択することにより得られる。こ
こでは、式２４と同様に、Ｕ_qが非単一値だけをその中に維持するＢΞ_iにおける
単一値分解の直交左側行列である。言い換えれば、次式を最小にすることである
。

【０１９４】

【数６８】

【０１９５】この問題の少なくとも１つの自乗解は次式になる。

【０１９６】

【数６９】

【０１９７】最初の故障の大きさｆ_i(t)は、式２５を用いて推定できる。

【０１９８】

【数７０】

【０１９９】ここで（）⁺は、当業者には既知のMoore-Penrose疑似逆関数である。構造化
残差ユニット１３０は、それが維持する各構造化残差に対して推定行列１３４

【０２００】

【数７１】

【０２０１】を維持する。故障を識別した構造化残残差に対応する推定行列を使用して、等式
エラー１２８を、前処理されたセンサ値スペース１３６内の故障ベクトルの推定
に写像する。次に、推定ユニットは、前処理ユニットからの逆拡大縮小化１１６
を適用して、最初のセンサスペース１６８内に故障ベクトルの推定を得る。

【０２０２】

【数７２】

【０２０３】図５は、センサ２のバイアス故障に対する推定ユニットの出力を示す。図６は
、センサ５のドリフト故障に対する推定ユニットの出力を示す。両者の場合、故
障サイズは正確に推定される。

【０２０４】７．代替ユニット代替ユニット１７０は、推定ユニットにより推定された（式５６）最初のセン
サスペース１６８における故障ベクトルの推定を単純に行い、センサ値１０８の
測定されたベクトルから、推定された故障を差し引いて代替値１７８のベクトル
を生成する。

【０２０５】

【数７３】

【０２０６】代替ユニットは、推定ユニット１６０が機能するときだけ機能し、また識別ユ
ニット１５０が識別イベント１５８を発生した場合だけ機能する。さらに安全機
能として、装置のオペレータは、代替ユニットの機能を確実に無効にできる。

【０２０７】８．分類ユニット分類ユニット１８０は、センサの故障のタイプに関する診断情報を計算する。
特に４タイプの故障、すなわち完全故障、バイアス、ドリフトおよび精度劣化が
考慮される。

【０２０８】分類ユニットは、分類イベント１５８が発生しているときだけ機能する。セン
サが最初に故障であると識別された時点から、測定されたセンサ値のバッファを
維持し、当該センサに対する推定故障サイズのバッファもまた維持する。これら
の２つのバッファにおける情報を分析し、それを用いて該当センサに対する故障
のタイプを分類する。分析は主として統計的推定を使用し、したがって分類ユニ
ットの調節可能なパラメータである有意水準αを必要とする。

【０２０９】説明を判りやすくするために、単一センサ故障が、識別ユニット１５０により
識別されていることを前提とする。分類ロジックのフローダイアグラムは、図３
に示されている。

【０２１０】８．１分類タイプ：完全故障完全故障は、故障センサの測定値１０８に関する回帰分析を実施することによ
り求められ、回帰線が良く一致し、ゼロの傾きをもつ統計推定により表示される
。

【０２１１】｛ｋ,ｕ_i(k)│ｋ＝ｔ_id,Λ,ｔ｝を回帰線が計算される点のシーケンスとする
。ここで、ｔ_idはセンサが故障と識別された点であり、またtは現在の時点であ
る。

【０２１２】標準最小自乗回帰を用いて、点と最良の一致線との間の残差値を最小にするこ
とにより、傾きβおよびｙ−切片αが次のように推定される。

【０２１３】

【数７４】

【０２１４】ここで

【０２１５】

【数７５】

【０２１６】傾きが統計的にゼロに等しいかどうかを調べるために、下記の統計値を作成す
る。

【０２１７】

【数７６】

【０２１８】ここで ^<βの標準偏差は次式で求められる。

【０２１９】

【数７７】

【０２２０】そして残差の標準偏差σは次式により推定される。

【０２２１】

【数７８】

【０２２２】ここで、

【０２２３】

【数７９】

【０２２４】は残差エラーである。式６０により得られる統計値は、ｔ−ｔ_id−１の自由度をもつスチューデントｔ
分布を満足する（例えばBirkesおよびDodge“Alternative Method of Regressio
n（回帰の代替方法）”、John Wiley & Sons 1993を参照のこと）。回帰線の一致度は、特定のしきい値を上回る線型相関係数により求められる。

【０２２５】８．２分類タイプ：バイアスバイアスは、推定ユニット１６０により計算された故障センサの故障サイズ推
定１６８に関する回帰分析を実施することにより求められ、回帰線が良く一致し
、ゼロの傾きをもち、ゼロでないバイアスをもつことを統計的に推定して示され
る。｛ｋ,ｆ_i(k)│ｋ＝ｔ_id,Λ,ｔ｝を、回帰線を計算する一連の点とする。ここで
、ｔ_idは、センサが故障と識別される点であり、tは現時点である。ゼロの傾き
に対する検定は、故障サイズシーケンスｆ_i(k)がセンサ値シーケンスｕ_i(k)に替
わる以外は、完全欠陥の分類と同じである。回帰線の一致度は、特定のしきい値
を越える線型相関係数により求められる。ゼロでないバイアスを検定するために次の統計値を用いる。

【０２２６】

【数８０】

【０２２７】これは、ｔ−ｔ_id−１の自由度をもつスチューデントｔ分布をも満足する。

【０２２８】８．３分類タイプ：ドリフトドリフトは、推定ユニット１６０により計算された故障センサの故障サイズ推
定１６８に関する回帰分析を行うことにより求められ、回帰線が良く適合し、ゼ
ロでない傾きをもつことを統計的に推定して示される。

【０２２９】８．４分類タイプ：精度劣化精度劣化は、推定ユニット１６０により計算された故障センサの故障サイズ推
定１６８に関する回帰分析を実施することにより求められ、回帰線が良く適合せ
ず、ゼロの傾きをもち、またゼロのバイアスをもつことを統計的に推定して示さ
れる。

【０２３０】精度劣化は、識別イベント１５８がセクション５．３のＮＣＵＭＶＡＲ検定に
より発生した場合にも表示される。

【０２３１】９．均等物上記の詳細は、本発明の好ましい実施形態に関するものである。本発明は、こ
れらの詳細には限定されるものではないと理解されるべきである。当業者によっ
て行われる代替および変更は、本発明の精神から逸脱することなく細部の実施形
態に対して実施できる。これらの改変は、本発明の範囲内にあるものとする。

【図面の簡単な説明】

【図１】本発明の好ましい実施形態の８つの演算ユニットおよびこれらのユニット間の
関係を示す。

【図２】図１の実施形態における識別プロセスのフローダイアグラムを示す。

【図３Ａ】図１の実施形態における分類プロセスのフローダイアグラムを示す。

【図３Ｂ】図１の実施形態における分類プロセスのフローダイアグラムを示す。

【図４】ボイラプロセスでの故障の発生時に、急速に増大し、しきい値を上回る本発明
の検出インデックスを示す。

【図５】ボイラプロセスのセンサ２におけるバイアス故障の発生時の推定ユニットの出
力を示す。

【図６】ボイラプロセスのセンサ５におけるドリフト故障の発生時の推定ユニットの出
力を示す。

【図７Ａ】ボイラプロセスのセンサ５におけるドリフト故障の発生時に応答する本発明の
ＮＳＦＳＲインデックスを示す。

【図７Ｂ】ボイラプロセスのセンサ５におけるドリフト故障の発生時に応答する本発明の
ＮＳＦＳＲインデックスを示す。

───────────────────────────────────────────────────── フロントページの続きＦターム(参考） 2F076 BA13 BA17 BA18 BD13 BE05 BE08 5H223 AA02 AA20 BB01 CC09 DD05 DD09 EE02 FF05

Claims

【特許請求の範囲】

【請求項１】測定されるプロセス中の1つまたは複数のセンサの故障を検
出する装置であって、測定されたセンサ値を含む信号の動作ベクトルを受け取り、前記測定されたセ
ンサ値を正規化して、前処理されたセンサ値を生成する前処理ユニットと、その前処理ユニットに結合され、前処理されたセンサ値を、主として測定雑音
を含む等式エラー値に変換するモデルユニットと、そのモデルユニットに結合され、構造化残差変換と称する複数の変換を含む構
造化残差ユニットであって、各変換が等式エラー値を構造化残差値に変換し、各
変換がセンサのサブセットの故障を感知しないように設計されている、構造化残
差ユニットと、前記モデルユニットに結合され、前記等式エラー値の間の関係を監視し、予測
された関係からの前記関係の有意偏差の発生により検出イベントを生成する検出
ユニットと、前記構造化残差ユニットおよび検出ユニットに結合され、検出イベントにより
起動されて、前記構造化残差値を用いて1つまたは複数のセンサの故障の有無を
判定し、その判定により識別イベントを生成する識別ユニットと、前記前処理ユニット、構造化残差ユニットおよび識別ユニットに結合され、識
別イベントにより起動されて、前記識別された故障センサのそれぞれについて故
障サイズを推定する推定ユニットと、その推定ユニットに結合され、識別されたすべての故障について、対応する前
記測定されたセンサ値から前記推定された故障サイズを差し引くことにより、前
記動作信号値中の前記故障の測定されたセンサ値について代替値を計算する代替
ユニットと、前記推定ユニットに結合され、前記推定ユニットが起動しているときに起動し
て、前記識別されたセンサ故障を故障タイプの固定セットに分類する分類ユニッ
トとを備えた装置。
【請求項２】請求項１において、前記前処理ユニットにおける正規化が、
前記測定されたセンサ値を拡大縮小およびオフセットすることにより達成される
装置。
【請求項３】請求項１において、前記動作信号が、前記測定されたセンサ
値に関する既知の状態の情報を含む装置。
【請求項４】請求項１において、前記モデルユニットにおける変換が、前
記前処理されたセンサ値に行列を乗じて前記等式エラーを生成することにより達
成される装置。
【請求項５】請求項４において、前記行列を主要コンポーネント分析の残
差部分から算出する装置。
【請求項６】請求項４において、前記行列を、各センサ値に１つの、１セ
ットの部分最小自乗モデルの残差部分から算出する装置。
【請求項７】請求項４において、前記行列を、前記測定されるプロセスの
質量バランスまたはエネルギーバランスから算出する装置。
【請求項８】請求項１において、前記構造化残差ユニットの各変換が前記
等式エラーとのベクトル点乗積の形をもつ装置。
【請求項９】請求項１において、前記構造化残差ユニットの各変換が、セ
ンサのサブセットの故障を感知せず、サブセット外の他のすべてのセンサの故障
を最大限感知するように設計されている装置。
【請求項１０】請求項１において、前記構造化残差ユニットの各変換が、
単一センサの故障を感知せず、他のすべてのセンサの故障を最大限感知するよう
に設計されている装置。
【請求項１１】請求項１において、各変換を定義する前記センサのサブセ
ットが、既知の不具合状態のすべてのセンサを少なくとも含む装置。
【請求項１２】請求項１において、前記検出ユニットが前記等式エラー値
の間の関係を、前記等式エラーの関数である検出インデックスを計算し、前記検出インデックスをしきい値と比較して、予測された関係からの前記関係
の有意偏差の発生を検出することにより監視する装置。
【請求項１３】請求項１２において、前記検出インデックスが、前記等式
エラー値を自乗して加算することにより得られる装置。
【請求項１４】請求項１２において、前記等式エラーに行列を乗じ前記等
式エラー値の相関性が解除されて、相関性を解除された等式エラー値が生成され
、前記検出インデックスが、前記相関性を解除された等式エラー値を自乗して加
算することにより得られる装置。
【請求項１５】請求項１３または１４において、前記しきい値が統計手法
を用いて決定される装置。
【請求項１６】請求項１３または１４において、前記検出インデックスが
少なくとも時間に関してフィルタ処理され、過渡現象および雑音の影響を平滑化
して、フィルタ処理された検出インデックスが生成され、前記フィルタ処理された検出インデックスを前記検出インデックスの代わりに
用いて、前記等式エラー値の間の関係を監視し、前記フィルタ処理された検出インデックスをしきい値と比較して、予測された
関係からの前記関係の有意偏差の発生を検出する装置。
【請求項１７】請求項１６において、前記検出インデックスが、指数関数
的に重み付けされた移動平均値フィルタの適用によりフィルタ処理される装置。
【請求項１８】請求項１７において、前記フィルタ処理された検出インデ
ックスについての前記しきい値が、自己相関関数を用いて前記検出インデックス
しきい値から計算される装置。
【請求項１９】請求項１において、前記識別ユニットが、各構造化残差値
を対応するしきい値と比較し、前記識別された構造化残差値と称する１つの構造化残差値を除くすべての前記
構造化残差値が、対応するしきい値を上回る場合には、識別イベントが発生し、前記判定された故障センサが前記センサのサブセットであり、設計により、そ
のサブセットについて、前記識別された構造化残差値に対応する前記構造化残差
変換が前記サブセットの故障の存在を感知しない装置。
【請求項２０】請求項１において、前記識別ユニットが、前記構造化残差
値に変換操作を適用して、変換された構造化残差インデックスを生成し、前記識別ユニットが、各変換された構造化残差インデックスを対応するしきい
値と比較し、前記変換され識別された構造化残差インデックスと称する1つの構造化残差を
除くすべての前記変換された構造化残差インデックスが、対応するしきい値を上
回る場合には、識別イベントが発生し、前記判定された故障センサが前記センサのサブセットであり、設計により、そ
のサブセットについて、前記変換された構造化残差インデックスに対応する前記
構造化残差変換が前記サブセットの故障の存在を感知しない装置。
【請求項２１】請求項２０において、前記変換操作が自乗および拡大縮小
の操作を含む装置。
【請求項２２】請求項２０において、前記変換操作が、指数関数的に重み
付けした、時間における移動平均フィルタを適用することを含む装置。
【請求項２３】請求項２２において、各変換された構造化残差インデック
スについての前記しきい値が、自己相関関数を用いて前記対応する構造化残差し
きい値から計算される装置。
【請求項２４】請求項１において、前記識別ユニットが、前記構造化残差
値に、前記故障検出イベントの時点から累積する累積操作を適用して、累積され
た構造化残差インデックスを生成し、前記識別ユニットが、各累積された構造化残差インデックスを対応するしきい
値と比較し、前記累積され識別された構造化残差インデックスと称する1つの構造化残差を
除くすべての前記累積された構造化残差インデックスが、対応するしきい値を上
回る場合には、識別イベントが発生し、前記判定された故障センサが前記センサのサブセットであり、設計により、そ
のサブセットについて、前記累積され識別された構造化残差インデックスに対応
する前記構造化残差変換が前記サブセットの故障の存在を感知しない装置。
【請求項２５】請求項２４において、前記累積操作で、前記故障検出イベ
ントの時点からの前記構造化残差の統計的平均における有意変化を検出するよう
設計された、一般化公算比インデックスを計算する装置。
【請求項２６】請求項２４において、前記累積操作で、前記故障検出イベ
ントの時点からの前記構造化残差の統計的平均を計算し、累積された平均インデ
ックスを生成する装置。
【請求項２７】請求項２４において、前記累積操作で、前記故障検出イベ
ントの時点からの構造化残差の統計的分散を計算し、累積された分散インデック
スを生成する装置。
【請求項２８】請求項１において、前記識別されたセンサ故障のそれぞれ
について推定される故障サイズを、前記等式エラー値への前記識別されたセンサ
故障の影響を最小にすることにより計算する装置。
【請求項２９】請求項１において、前記分類ユニットが、識別された各故
障センサに対し１つずつ、複数の回帰線を計算し、この各回帰線は前記識別の時
点からの前記対応する故障センサについて測定されたセンサ値を用いて計算し、
複数のセンサ値回帰線を生成し、前記分類ユニットが、識別された各故障センサに対し１つずつ、複数の回帰線
を計算し、この各回帰線は前記識別の時点からの前記対応する故障センサについ
て推定される故障サイズを用いて計算し、複数の故障サイズ回帰線を生成し、前記分類ユニットが、センサ値回帰線ならびに故障サイズ回帰線からの傾き、
バイアスおよび残差情報を用いて、前記識別されたセンサ故障を故障タイプの固
定セットに分類する装置。
【請求項３０】請求項２９において、前記対応するセンサ値回帰線の残差
がしきい値を下回り、前記対応するセンサ値回帰線の傾きが統計的にゼロに等し
い場合には、故障センサを“完全故障”タイプとして分類する装置。
【請求項３１】請求項２９において、前記対応する故障サイズ回帰線の残
差がしきい値を下回り、前記対応する故障サイズ回帰線の傾きが統計的にゼロに
等しく、前記対応する故障サイズ回帰線のバイアスが統計的にゼロに等しくない
場合には、故障センサを “バイアス” タイプとして分類する装置。
【請求項３２】請求項２９において、前記対応する故障サイズ回帰線の残
差がしきい値を下回り、前記対応する故障サイズ回帰線の傾きが統計的にゼロに
等しくない場合には、故障センサを “ドリフト” タイプとして分類する装置。
【請求項３３】請求項２９において、前記対応する故障サイズ回帰線の残
差がしきい値を上回り、前記対応する故障サイズ回帰線の傾きおよびバイアスが
ともに統計的にゼロに等しい場合には、故障センサを “精度劣化” タイプとし
て分類する装置。
【請求項３４】請求項２６において、前記識別イベントが、対応するしき
い値を上回る累積された分散インデックスにより発生する場合には、故障センサ
を “精度劣化” タイプとして分類する装置。
【請求項３５】コンピュータを用いて、測定されるプロセスにおける１つ
または複数のセンサ故障を検出する方法であって、測定されたセンサ値を含む信号の動作ベクトルを受け取り、前記測定されたセ
ンサ値を前処理して、前処理されたセンサ値を生成するステップと、前記前処理されたセンサ値を、主として測定雑音を含む等式エラー値に変換す
るステップと、複数の変換を前記等式エラー値に適用して複数の構造化残差値を生成するステ
ップであって、前記変換が、構造化残差変換と称され、センサのサブセットの故
障を感知しないように設計されているステップと、前記等式エラー値の間の関係を監視し、予測された関係からの前記関係の有意
偏差の発生により検出イベントを生成するステップと、検出イベントが発生した場合に、前記構造化残差値を用いて1つまたは複数の
センサの故障の有無を判定し、その判定により識別イベントを生成するステップ
と、識別イベントが発生した場合に、前記識別された故障センサのそれぞれについ
て故障サイズを推定するステップと、識別イベントが発生した場合、識別されたすべての故障について、対応する前
記測定されたセンサ値から前記推定された故障サイズを差し引くことにより、前
記動作信号中の故障の測定されたセンサ値を修正された値に置き換えるステップ
と、識別イベントが発生した場合、前記識別されたセンサ故障を故障タイプの固定
セットに分類するステップとを備えた方法。
【請求項３６】請求項３５において、前記測定されたセンサ値を前処理す
るステップが、拡大縮小およびオフセットすることにより達成される方法。
【請求項３７】請求項３５において、前記動作信号が、前記測定されたセ
ンサ値に関する既知の状態の情報を含む方法。
【請求項３８】請求項３５において、前記前処理されたセンサ値を等式エ
ラー値に変換するステップが、前記前処理されたセンサ値に行列を乗じることに
より達成される方法。
【請求項３９】請求項３８において、前記行列を主要コンポーネント分析
の残差部分から算出する方法。
【請求項４０】請求項３８において、前記行列を、各センサ値に１つの、
１セットの部分最小自乗モデルの残差部分から算出する方法。
【請求項４１】請求項３８において、前記行列を、前記測定されるプロセ
スの質量バランスまたはエネルギーバランスから算出する方法。
【請求項４２】請求項３５において、各構造化残差変換が、前記等式エラ
ーへの点乗積として適用される係数のベクトルからなる方法。
【請求項４３】請求項３５において、各構造化残差変換が、センサのサブ
セットの故障を感知せず、サブセット外の他のすべてのセンサの故障を最大限感
知するように設計されている方法。
【請求項４４】請求項３５において、各構造化残差変換が、単一センサの
故障を感知せず、他のすべてのセンサの故障を最大限感知するように設計されて
いる方法。
【請求項４５】請求項３５において、各変換を定義する前記センサのサブ
セットが、既知の不具合状態のすべてのセンサを少なくとも含む方法。
【請求項４６】請求項３５において、前記等式エラー値の間の関係を監視
するステップが、前記等式エラーの関数である検出インデックスを計算することと、前記検出インデックスをしきい値と比較して、予測された関係からの前記関係
の有意偏差の発生を検出することとからなる方法。
【請求項４７】請求項４６において、前記検出インデックスが、前記等式
エラー値の自乗値を加算することにより得られる方法。
【請求項４８】請求項４６において、前記等式エラーに行列を乗じ前記等
式エラー値の相関性が解除されて、相関性を解除された等式エラー値が生成され
、前記検出インデックスが、前記相関性を解除された等式エラー値を自乗して加
算することにより得られる方法。
【請求項４９】請求項４７または４８において、前記しきい値が統計技法
を用いて決定される方法。
【請求項５０】請求項４７または４８において、前記検出インデックスが
少なくとも時間に関してフィルタ処理され、過渡現象および雑音の影響を平滑化
して、フィルタ処理された検出インデックスが生成され、前記フィルタ処理された検出インデックスを前記検出インデックスの代わりに
用いて、前記等式エラー値の間の関係を監視し、前記フィルタ処理された検出インデックスをしきい値と比較して、予測された
関係からの前記関係の有意偏差の発生を検出する方法。
【請求項５１】請求項５０において、前記検出インデックスが、指数関数
的に重み付けされた移動平均フィルタの適用によりフィルタ処理される方法。
【請求項５２】請求項５１において、前記フィルタ処理された検出インデ
ックスについての前記しきい値が、自己相関関数を用いて前記検出インデックス
しきい値から計算される方法。
【請求項５３】請求項３５において、１つまたは複数のセンサの故障の有
無を判定するステップが、各構造化残差値を対応するしきい値と比較することと、前記識別された構造化残差値と称する１つの構造化残差値を除くすべての前記
構造化残差値が、対応するしきい値を上回る場合には、識別イベントを発生する
ことと、前記故障センサが前記センサのサブセットであると判定することとからなり、
そのサブセットについて、設計により、前記識別された構造化残差値に対応する
構造化残差変換が前記サブセットの故障の存在を感知しない方法。
【請求項５４】請求項３５において、前記構造化残差値に変換操作を適用
して、変換された構造化残差インデックスを生成する方法であって、１つまたは
複数のセンサの故障の有無を判定するステップが、各変換された構造化残差インデックスを対応するしきい値と比較することと、前記変換され識別された構造化残差インデックスと称する構造化残差の1つを
除くすべての前記変換された構造化残差インデックスが、対応するしきい値を上
回る場合には、識別イベントを発生することと、前記故障センサが前記センサのサブセットであると判定することとからなり、
そのサブセットについて、設計により、前記変換された構造化残差インデックス
に対応する構造化残差変換が前記サブセットの故障の存在を感知しない方法。
【請求項５５】請求項５４において、前記変換操作が自乗および拡大縮小
の操作を含む方法。
【請求項５６】請求項５４において、前記変換操作が、指数関数的に重み
付けした、時間における移動平均フィルタを適用することを含む方法。
【請求項５７】請求項５６において、各変換された構造化残差インデック
スについての前記しきい値が、自己相関関数を用いて前記対応する構造化残差し
きい値から計算される方法。
【請求項５８】請求項３５において、前記構造化残差値に、前記故障検出
イベントの時点から累積する累積操作を適用して、累積された構造化残差インデ
ックスを生成する方法であって、1つまたは複数のセンサの故障の有無を判定す
るステップが、各累積された各構造化残差インデックスを対応するしきい値と比較することと
、前記累積され識別された構造化残差インデックスと称する１つの構造化残差を
除くすべての前記累積された構造化残差インデックスが、対応するしきい値を上
回る場合には、識別イベントを発生することと、前記故障センサが前記センサのサブセットであると判定することとからなり、
そのサブセットについて、設計により、前記累積され識別された構造化残差イン
デックスに対応する前記構造化残差変換が前記サブセットの故障の存在を感知し
ない方法。
【請求項５９】請求項５８において、前記累積操作で、前記故障検出イベ
ントの時点からの前記構造化残差の統計的平均における有意変化を検出するよう
設計された、一般化公算比インデックスを計算する方法。
【請求項６０】請求項５８において、前記累積操作で、前記故障検出イベ
ントの時点からの前記構造化残差の統計的平均を計算し、累積された平均インデ
ックスを生成する方法。
【請求項６１】請求項５８において、前記累積操作で、前記故障検出イベ
ントの時点からの構造化残差の統計的分散を計算し、累積された分散インデック
スを生成する方法。
【請求項６２】請求項３５において、前記識別されたセンサ故障のそれぞ
れについて推定される故障サイズを、前記等式エラー値への前記識別されたセン
サ故障の影響を最小にすることにより計算する方法。
【請求項６３】請求項３５において、前記識別されたセンサ故障を故障タ
イプの固定セットに分類するステップが、識別された各故障センサに対し１つずつ、複数の回帰線を計算し、この各回帰
線は前記識別の時点からの前記対応する故障センサについて測定されたセンサ値
を用いて計算し、複数のセンサ値回帰線を生成することと、識別された各故障センサに対し１つずつ、複数の回帰線を計算し、この各回帰
線は前記識別の時点からの前記対応する故障センサについて推定される故障サイ
ズを用いて計算し、複数の故障サイズ回帰線を生成することと、センサ値回帰線ならびに故障サイズ回帰線からの傾き、バイアスおよび残差情
報を用いて、前記識別されたセンサ故障を故障タイプの固定セットに分類するこ
ととにより達成される方法。
【請求項６４】請求項６３において、前記対応するセンサ値回帰線の残差
がしきい値を下回り、前記対応するセンサ値回帰線の傾きが統計的にゼロに等し
い場合には、故障センサを “完全故障” タイプとして分類する方法。
【請求項６５】請求項６３において、前記対応する故障サイズ回帰線の残
差がしきい値を下回り、前記対応する故障サイズ回帰線の傾きが統計的にゼロに
等しく、前記対応する故障サイズ回帰線のバイアスが統計的にゼロに等しくない
場合には、故障センサを “バイアス” タイプとして分類する方法。
【請求項６６】請求項６３において、前記対応する故障サイズ回帰線の残
差がしきい値を下回り、前記対応する故障サイズ回帰線の傾きが統計的にゼロに
等しくない場合には、故障センサを “ドリフト” タイプとして分類する方法。
【請求項６７】請求項６３において、前記対応する故障サイズ回帰線の残
差がしきい値を上回り、前記対応する故障サイズ回帰線の傾きおよびバイアスが
ともに統計的にゼロに等しい場合には、故障センサを “精度劣化”タイプとし
て分類する方法。
【請求項６８】請求項６０において、前記識別イベントが、対応するしき
い値を上回る累積された分散インデックスにより発生する場合には、故障センサ
を “精度劣化”タイプとして分類する方法。