(57)【要約】【課題】学習データの完全性、異常の混入を許容でき、プラントなどの設備において異常の早期・高精度な発見を可能とする方法及びそのシステムを提供する。【解決手段】上記目的を達成するために、本発明は、（１）時間的なデータの振舞いに着目し、時間を追って軌跡をクラスタに分割する。（２）分割したクラスタ郡に対して、部分空間でモデルし、はずれ値を異常候補として算出する。（３）学習データをリファレンスとして活用（比較・参照など）し、経時変化、環境変動、保守（部品交換）、稼動状態による状態遷移を把握する。（４）モデル化は、データのN個抜き（N=0,1,2,・・・）の回帰分析法や投影距離法等の部分空間法（例えばN=1の場合は、異常データが１個混入していると考え、これを除いてモデル化する）、或いは局所部分空間法によるものとする。なお、回帰分析法における直線の当てはめは、最も低次の回帰分析に相当する。【選択図】図１
2010092355.tif 000002
【特許請求の範囲】
【請求項１】
プラントまたは設備の異常を早期に検知する異常検知方法であって、
複数のセンサからデータを取得し、前記データの時間的な変化に基づいて、データ空間の軌跡を複数のクラスタに分割し、
前記分割した前記クラスタ群に対して、部分空間法でモデル化し、
はずれ値を異常候補として算出することを特徴とする異常検知方法
【請求項２】
前記モデル化は、前記複数のセンサから過去に得られたデータからなる学習データを対象にモデル化したリファレンスモデルを併用して、はずれ値を異常候補として算出することを特徴とする請求項１に記載の異常検知方法。
【請求項３】
前記モデル化は、データのN個抜き（N=0,1,2,・・・）の回帰分析法、または、部分空間法を用いて行うことを特徴とする請求項１記載の異常検知方法。
【請求項４】
前記部分空間法は、投影距離法、ＣＬＡＦＩC法、または、観測データの近傍を対象とする局所部分空間法であることを特徴とする請求項３記載の異常検知方法。
【請求項５】
前記はずれ値の検出は、前記回帰分析法、前記投影距離法、前記ＣＬＡＦＩC法、または、観測データの近傍を対象とする前記局所部分空間法の出力を統合して行うものであることを特徴とする請求項４記載の異常検知方法。
【請求項６】
プラントまたは設備の異常を早期に検知する方法において、
複数のセンサからデータを取得し、
時間的なデータの時間的変化に基づいて、データ空間の軌跡を複数のクラスタに分割し、
分割したクラスタが分かるように、データを表示することにより、状態を可視化することを特徴とする異常検知方法。
【請求項７】
前記複数のセンサからの過去のデータに基づくクラスタ群に対する異常検知により、第１報を出力し、前記複数のセンサの現在のデータに基づくクラスタ群に対する異常検知により第２報を出力することを特徴とする請求項１記載の異常検知方法。
【請求項８】
プラントまたは設備の異常を早期に検知する異常検知システムにおいて、
複数のセンサからデータを取得するデータ取得部と、
前記データの時間的な変化に基づいて、データ空間の軌跡を複数のクラスタ群に分割する軌跡分割部と、
分割した前記クラスタ群に対して、部分空間でモデル化するモデル化部と、
はずれ値を異常候補として算出する異常検出部を有し、
前記軌跡分割部においては前記分割したクラスタ群に対して、部分空間でモデル化し、はずれ値を異常候補として算出することを特徴とする異常検知システム
【請求項９】
前記データ取得部において取得した前記複数のセンサからの過去のデータを蓄積し、かつ、モデル化してリファレンスモデルを作成する学習データベース部を有し、前記学習データベース部における前記リファレンスモデルを併用して、はずれ値を異常候補として算出することを特徴とする請求項８記載の異常検知システム。
【請求項１０】
前記モデル化は、データのN個抜き（N=0,1,2,・・・）の回帰分析法、または、部分空間法を用いて行うことを特徴とする請求項８記載の異常検知システム。
【請求項１１】
前記部分空間法は、投影距離法、ＣＬＡＦＩＳ法、または、観測データの近傍を対象とする局所部分空間法であることを特徴とする請求項８記載の異常検知システム。
【請求項１２】
前記はずれ値の検出は、前記回帰分析法、前記投影距離法、前記ＣＬＡＦＩＳ法、または、観測データの近傍を対象とする前記局所部分空間法の出力を統合して行うものであることを特徴とする請求項８記載の異常検知システム。
【請求項１３】
プラントまたは設備の異常を早期に検知する異常システムにおいて、
複数のセンサからデータを取得するデータ取得部と、
前記データの時間的な変化に基づいて、データ空間の軌跡を複数のクラスタに分割する軌跡分割部と、
前記分割したクラスタ群に対して、部分空間でモデル化するモデル化部と、
はずれ値を異常候補として算出する異常検出部と、
データ表示部とを有し、
前記分割した複数のクラスタが分かるように、データを表示することにより、
状態を可視化することを特徴とする異常検知システム。
【請求項１４】
前記複数のセンサからの過去のデータに基づくクラスタ群に対する異常検知により、第１報を出力し、前記複数のセンサの現在のデータに基づくクラスタ群に対する異常検知により第２法を出力することを特徴とする請求項１３記載の異常検知システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プラントや設備などの異常を早期に検知する方法に関する。
【背景技術】
【０００２】
電力会社では、ガスタービンの廃熱などを利用して地域暖房用温水を供給したり、工場向けに高圧蒸気や低圧蒸気を供給したりしている。石油化学会社では、ガスタービンなどを電源設備として運転している。このようにガスタービンなどを用いた各種プラントや設備において、その異常を早期に発見することは、社会へのダメージを最小限に抑えることができ、極めて重用である。
【０００３】
ガスタービンや蒸気タービンのみならず、水力発電所での水車、原子力発電所の原子炉、風力発電所の風車、航空機や重機のエンジン、鉄道車両、エスカレータ、エレベータ、機器・部品レベルでも、搭載電池の劣化・寿命など、早期に異常を発見しなければならない設備は枚挙に暇がない。最近では、健康管理のため、脳波測定・診断に見られるように、人体に対する異常（各種症状）検知も重要になりつつある。
【０００４】
このため、例えば米国のSmartSignal社では、US6,952,662やUS6,975,962に記載のように、おもにエンジンを対象に、異常検知の業務をサービスしている。そこでは、過去のデータをデータベース（DB）としてもっておき、観測データと過去の学習データとの類似度を独自の方法で計算し、類似度の高いデータの線形結合により推定値を算出して、推定値と観測データのはずれ度合いを出力する。General Electric社のように、US6,216,066出願内容を見ると、異常検知をk-meansクラスタリングにより検出している例もある。
【０００５】
【特許文献１】US6,952,662
【特許文献２】US6,975,962
【特許文献３】US6,216,066
【非特許文献１】Stephan W. Wegerich；Nonparametric modeling of vibration signal features for equipment health monitoring、Aerospace Conference, 2003. Proceedings. 2003 IEEE，Volume 7, Issue, 2003 Page(s):3113 - 3121
【発明の開示】
【発明が解決しようとする課題】
【０００６】
SmartSignal社が用いている手法では、データベースに格納する過去の学習データは様々な状態を網羅的に含む必要がある。もし、学習データにない観測データが観察されると、これらはすべて学習データには包含されないものとして扱われ、はずれ値と判断され、正常信号においても異常と判定されてしまい、検査の信頼性が著しく低下してしまう。このため、ユーザは、過去のあらゆる状態のデータすべてをDBとして蓄えることが必須になる。
【０００７】
一方で、学習データに異常が混入されている場合、異常を表している観測データとの乖離度が低くなり、これを見逃してしまう。このため、学習データに異常が含まれないよう、十分なるチェックが必要である。
【０００８】
このように、SmartSignal社が提案する学習データに基づく方法では、ユーザは、網羅的なデータ収集と異常の排除という負荷を背負うことになる。特に、経時変化や、周囲の環境変動、部品交換などの保守作業の有無などに、きめこまかく対応する必要がある。このような対応を行うことは、実質難しく、不可能な場合が多い。
【０００９】
General Electric社の方法では、k-meansクラスタリングのため、信号の挙動を見ておらず、その点で本質的な異常検知にはなっていない。
【００１０】
そこで、本発明の目的は、上記課題を解決し、学習データの不完全性や異常の混入を許容し、ユーザ負荷を軽減することを図り、さらに早期に異常を高感度に検知することが可能な異常検知方法およびシステムを提供することにある。
【課題を解決するための手段】
【００１１】
上記目的を達成するために、本発明は、（１）時間的なデータの振舞いに着目し、時間を追って軌跡をクラスタに分割する、（２）分割したクラスタ群に対して、部分空間でモデルし、はずれ値を異常候補として算出する。（３）学習データをリファレンスとして活用（比較・参照など）し、経時変化、環境変動、保守（部品交換）、稼動状態による状態遷移を把握する。（４）モデル化は、データのN個抜き（N=0,1,2,・・・）の回帰分析法や投影距離法等の部分空間法（例えばN=1の場合は、異常データが１個混入していると考え、これを除いてモデル化する）、或いは局所部分空間法によるものとする。
【００１２】
さらに、（５）部分空間法など、複数の識別器の出力を統合して。異常判断を行う。なお、回帰分析法における直線の当てはめは、最も低次の回帰分析に相当する。
【発明の効果】
【００１３】
本発明によれば、学習データが完全でなくとも、異常の混入を許容でき、プラントなどの設備において異常の早期・高精度な発見を可能とする。即ち、SmartSignal社のように正常部の学習、複数の状態に応じた網羅的データ収集の必要性がない。
【発明を実施するための最良の形態】
【００１４】
以下実施例に従って、本発明の内容を詳細に説明する。
【実施例１】
【００１５】
説明を簡略化するために、本願においては、同一の構成要素には同一の符号を付加している。本発明の実施の一形態を図1に、対象とする信号の例を図２に示す。対象信号は、図２に示すような複数の時系列信号である。ここでは、系列１，２，３、４という４種類の信号を表している。実際には、信号は４種類に限るものではなく、数百から数千と言った数になる場合もある。
【００１６】
各信号が、対象プラントや設備に設けられた複数のセンサからの出力に相当する。例えば、シリンダ、オイル、冷却水などの温度、オイルや冷却水の圧力、軸の回転速度、室温、運転時間などが、いろいろなセンサから、１日に数回ないしリアルタイム等、定めた間隔で観測される。出力や状態を表すのみならず、何かを制御するための制御信号（入力）の場合もある。ON/OFF制御のこともあれば、一定値になるように制御する場合もある。これらのデータは、お互いに相関が高いものや低いものがある。これら全ての信号が、対象になり得る。これらのデータを見て、異常の有無を判断する。ここでは、多次元時系列信号として扱う。
【００１７】
図１に示した異常検知方法について説明する。まず、多次元時系列信号を多次元信号取得部１にて取得する。次に、取得した多次元時系列信号が、欠損している場合があるため、欠損値の修正・削除を欠損値の修正・削除部２にて行う。欠損の修正は、例えば、前後のデータの置換えや移動平均の置換えが一般的である。削除は、多くのデータが同時に０にリセットされた時など、データとしての異常を排除するものである。欠損値の修正・削除は、設備の状態やエンジニアの知識を状態データ／知識３というDBに蓄えておき、これをもとに行う場合もある。次に、修正・削除された多次元時系列信号に関して、相関解析による無効信号の削除を、相関解析による無効信号の削除部４にて行う。これは、図３にて相関行列の例を示すように、多次元時系列信号に対して相関解析を行い、相関値が１に近い複数の信号がある場合など、極めて類似性が高い場合に、これらは冗長だとして、この複数の信号から重複する信号を削除し、重複しないものを残すものである。この場合も、状態データ／知識３に格納した情報に基づいて、削除する。
【００１８】
次に、データの次元削減を主成分分析部５にて行う。ここでは、主成分分析により、Ｍ次元の多次元時系列信号を、次元数ｒのｒ次元多次元時系列信号に線形変換する。主成分分析は、ばらつき最大となる軸を生成するものである。KL変換でも構わない。次元数ｒは、主成分分析により求めた固有値を降順に並べ、大きい方から加算した固有値を全固有値の和で割り算した累積寄与率なる値に基づいて決める。
【００１９】
次に、ｒ次元の多次元時系列信号に対して、軌跡分割によるクラスタリングを軌跡分割によるクラスタリング部６にて行う。図４にその模様を示す。図４の左上の３次元表示（特徴空間と呼ぶ）は、主成分分析後のｒ次元の多次元時系列信号を寄与率の高い３次元で表示したものである。この状態では、まだ対象設備の状態が複雑なものとして観測されていることが分かる。図４の残りの八つの３次元表示は、軌跡を時間に沿って追跡し、クラスタリングを実施したものであり、各クラスタを表している。クラスタリングは、時間に沿ってデータ間の距離が定めたしきい値を超えれば、別のクラスタとし、しきい値を超えなければ、同じクラスタとして扱うものである。これにより、クラスタ１、３、９、１０、１７は、運転ＯＮの状態のクラスタであり、クラスタ６、１４、２０は、運転ＯＦＦの状態のクラスタに別れていることが分かる。クラスタ２など、図示していないクラスタは、過渡期のものである。これらのクラスタを分析すると、運転ＯＮの状態では、線状に軌跡移動していることが伺え、運転ＯＦＦの状態では、不安定な軌跡移動していることが伺える。このように、軌跡分割によるクラスタリングは、いくつかのメリットがあることが分かる。
運転ＯＮの状態、運転ＯＦＦの状態など、複数の状態に分類できる。
運転ＯＮの状態に見られるように、これらのクラスタは、例えば線形など、低次元のモデルで表現できる。
これらのクラスタリングは、設備のアラーム信号や保全情報を加味して、これらのひも付きとして、実施してもよい。具体的には、各クラスタに、アラーム信号などの情報が属性として付加される。
【００２０】
図５に、特徴空間にて、クラスタリングによりラベルを付けた状態のほかの例を表す。図６は、一つの時系列信号の上に、クラスタリングの結果を表示したものである。この場合、クラスタは１６個生成でき、時系列信号が、１６個のクラスタに分割されている様子が分かる。運転時間（累積時間）も重ねて表示した。水平になっている部分は、運転ＯＦＦである。運転ＯＮと運転ＯＦＦが、精度良く分離できていることが分かる。
【００２１】
次に、クラスタリングされた各クラスタを対象に、低次元部分空間でのモデル化を、各クラスタを対象にモデル化部８にて行う。正常部限定である必要はなく、異常が混入することは問題ない。ここでは、例えば、回帰分析によりモデル化を行う。回帰分析の一般式は、下記のとおりである。「ｙ」が、クラスタごとのｒ次元の多次元時系列信号に相当する。「X」が、ｙを説明する変数である。「y~」がモデルになる。「ｅ」は偏差である。
y:目的変数（ｒ列）
ｂ：回帰係数（1+p列）
Ｘ：説明変数行列(ｒ行、1+p列)
||y-Xb||⇒min
b=(X'X)-1X'y（'は転置を表す）
y~=Xb=X(X'X) -1X'y（説明変数の影響を表す部分）
e=y-y~（y~で近似できない部分。説明変数の影響を取り除いた部分）
ただし、rank X=p+1
ここでは、各クラスタのｒ次元の多次元時系列信号に対し、データのN個抜き（N=0,1,2,・・・）の回帰分析を行う。例えばN=１の場合は、異常信号が１種混入していると考え、これを除いた信号を「X」としてモデル化する。Ｎ＝０の場合は、全ｒ次元多次元時系列信号を扱うことになる。
【００２２】
回帰分析以外に、CLAFIC法や投影距離法等の部分空間法を適用してもよい。そして、モデルからの偏差をモデルからの偏差算出部９にて求める。図２４に一般的なCLAFIC法を図解して示す。２クラス、2次元パターンの場合を示す。各クラスの部分空間、即ち、ここでは１次元の直線として表される部分空間を求める。
【００２３】
一般的には、各クラスのデータの自己相関行列を固有値分解して、固有ベクトルを基底として求める。値が大きい、上位何個かの固有値に対応する固有ベクトルを用いる。未知パターンｑ（最新の観測パターン）が入力されると、部分空間への正射影の長さ、或いは部分空間への投影距離を求める。そして、正射影の長さが最大、或いは投影距離が短いクラスに、未知パターン（最新の観測パターン）ｑを分類する。
【００２４】
図２４では未知パターンｑ（最新の観測パターン）は、クラスＡに分類される。図２に示した多次元時系列信号では、基本的に正常部を対象とするため、１クラス識別の問題になるので（図１にて図示）、クラスＡを正常部とし、未知パターンｑ（最新の観測パターン）からクラスＡまでの距離を求めて、これを偏差とする。そして、偏差が大きいと、はずれ値と判断する。このような部分空間法では、異常値が若干混ざっていても、次元削減し、部分空間にした時点で、その影響が緩和される。部分空間法適用のメリットである。
【００２５】
なお、投影距離法では、各クラスの重心を原点とする。各クラスの共分散行列にＫＬ展開を適用して得られた固有ベクトルを基底として用いる。いろいろな部分空間法が立案されているが、距離尺度を有するものならば、はずれ度合いが算出可能である。なお、密度の場合も、その大小により、はずれ度合いを判断可能である。CLAFIC法は、正射影の長さを求めることから、類似度尺度である。
【００２６】
このように、部分空間にて距離や類似度を計算し、はずれ度合いを評価することになる。投影距離法などの部分空間法は、距離に基づく識別器のため、異常データが利用できる場合の学習法として、辞書パターンを更新するベクトル量子化や距離関数を学習するメトリック学習を使うことができる。
【００２７】
また、局所部分空間法と呼ぶ、未知パターンｑ（最新の観測パターン）に近いｋ個の多次元時系列信号を求め、各クラスの最近傍パターンが原点となるような線形多様体を生成し、その線形多様体への投影距離が最小となるクラスに未知パターンを分類する方法も適用できる（図３１の局所部分空間法の枠内参照）。局所部分空間法も部分空間法の一種である。
【００２８】
既に説明したクラスタリング後の各クラスタに対して、局所部分空間法を適用することになる。ｋは、パラメータである。異常検知では、先ほどと同様に、１クラス識別の問題になるので、大多数のデータが属するクラスＡを正常部とし、未知パターンｑ（最新の観測パターン）からクラスＡまでの距離を求めて、これを偏差とする。
【００２９】
この手法では、例えば、ｋ個の多次元時系列信号を用いて形成される部分空間への、未知パターンｑ（最新の観測パターン）からの正射影した点を推定値として算出することもできる（図３１の局所部分空間法の枠内に説明した推定値なるデータ）。また、ｋ個の多次元時系列信号を、未知パターンｑ（最新の観測パターン）に近い順に並べ替え、その距離に反比例した重み付けを行って、各信号の推定値を算出することもできる。投影距離法などでも、同様に推定値を算出できる。
【００３０】
パラメータｋは通常は１種類に定めるが、パラメータｋをいくつか変えて実行すると、類似度に応じて対象データを選択することになり、それらの結果から総合的に判断できるため、一層効果的である。局所部分空間法では、クラスタ内の選択したデータを対象とするため、異常値が若干混ざっていても、局所部分空間にした時点で、その影響が大きく緩和される。
【００３１】
局所部分空間法の「局所」という概念は、回帰分析にも適用可能である。即ち、「ｙ」として、観測未知パターンｑに近いｋ個の多次元時系列信号を求め、このｙのモデルとして、「y~」を求め、偏差「e」を算出する。
【００３２】
なお、単に、１クラス識別の問題と考えれば、１クラスサポートベクターマシンなどの識別器も適用可能である。この場合、高次空間に写像する、radial basis functionなどのカーネル化が使えることになる。１クラスサポートベクターマシンでは、原点に近い側が、はずれ値、即ち異常になる。ただし、サポートベクターマシンは、特徴量の次元は大きくても対応できるが、学習データ数が増えると計算量が膨大となるという欠点もある。
【００３３】
このため、MIRU２００７（画像の認識・理解シンポジウム、Meeting on Image Recognition and Understanding 2007）にて発表されている、「IS-2-10  加藤丈和，野口真身，和田俊和（和歌山大），酒井薫，前田俊二（日立）；パターンの近接性に基づく１クラス識別器」などの手法も適用可能であり、この場合、学習データ数が増えても、計算量は膨大なものとならないというメリットがある。
【００３４】
次に、回帰分析を例にとり、実験例を説明する。図７に、Ｎ＝０とし、線形回帰分析により、ｒ次元多次元時系列信号をモデル化し、そのモデルと実測値の偏差を図示した例を示す。図８は、参考に、軌跡分割によるクラスタリングを実施しない場合である。図７の場合、偏差が大きいのは、運転ＯＦＦの区間、および、運転ＯＮの区間で時系列信号が振動的な振舞いをしている場合である。最後に、はずれ値検出部１０にてはずれ値を求める。ここでは、しきい値との大小をチェックする。検出された異常信号は、主成分分析後であるため、これを逆変換し、元の信号がどのような割合で合成されて、異常と判断されたかを確認することもできる。
【００３５】
このように、軌跡分割によるクラスタリングを中心に、低次元モデルで多次元時系列信号を表現することにより、複雑な状態を分解でき、簡単なモデルで表現できるため、現象を理解しやすいという利点がある。また、モデルを設定するため、SmartSignal社の方法のように完全に、データを完備する必要はない。データの欠落があっても良いというメリットがある。
【００３６】
次に、局所部分空間法の適用例を図２８に示す。信号を前半と後半に分け（交差確認と呼ぶ検証方法に従った）、それぞれを学習データとして、残りのデータまでの距離を求めた例である。パラメータｋは１０とした。ｋをいくつか変えて、それらの多数決をとれば、安定な結果が得られる（後述するバギングという手法と同様の考えに基づく）。この局所部分空間法では、データのN個抜きが、自動で行われるというメリットがある。同図適用例では、運転OFF中の不規則な挙動が、検知されている。
【００３７】
上記例では、クラスタリングの必要性も緩和されるが、観測データが属しているクラスタ以外を学習データとし、このデータと観測データに局所部分空間法を適用してもよい。この方法によれば、別クラスタからの乖離度を評価できる。投影距離法も同様である。図２９にそれらの例を示す。観測データが属しているクラスタ以外を学習データとした。この考えは、時系列データのように類似データが連続する場合に、最も類似したデータを「局所」領域から排除できるため、効果的である。なお、データのN個抜きは、特徴量（センサ信号）として説明したが、時間軸方向のデータであってもよい。
【００３８】
次に、データの表現形態を、いくつか図を用いて説明する。図９に幾つかの例を示す。図９の左側の図は、主成分分析後のｒ次元時系列信号を2次元表示したものである。データの振舞いを可視化した例となる。図９の右側の図は、軌跡分割によるクラスタリングを実施し、クラスタを図示したものである。クラスタ毎に、簡単な低次のモデル（ここでは直線）で表現した例である。図１０の左側の図は、データの動く速さが分かるように図示した例である。後述のWavelet解析を適用すれば、速さ、即ち周波数も分析可能、そして多変量として扱うことが可能となる。図１０の右側の図は、図９の右側の図で示したモデルからの偏差が分かるように表示した例である。
【００３９】
図１１の左側の図は、別の例である。距離基準などに基づき類似と判定したクラスタをマージし（同図では、隣接したクラスタのマージを示した）、マージ後のモデルを示し、かつモデルからの偏差を図示した例である。図１１の右側の図は、状態を表現したものである。状態Ａ、Ｂ、Ｃという３種類の状態を分けて表示した。状態を分けて考えると、図１２の左側の図に示すように、状態Ａの変化などを図示可能となる。
【００４０】
図６の例で考えると、同じ運転ＯＮの状態でも、運転ＯＦＦの前後で、異なる挙動を示しており、これらを特徴空間の中で表現できることなる。図１２の右側の図は、過去の学習データから得られたモデル（低次の部分空間）からの変化を示したものであり、状態変化を観察できる。このように、データを加工し、加工したデータをユーザに示し、現在の状況を可視化することにより、より良い理解を促すことができる。
【００４１】
次に、別の実施例を説明する。既に説明したブロックは省略する。図１３に異常検知方法を示す。ここでは、各クラスタに対し、ｒ次元多次元時系列信をランダムに定めた個数分選択する。ランダム選択により、
全ての信号を使う場合に見えなかった特性が現れる
有効ではない信号を除く
全ての組み合わせより、短時間で計算できる
という利点がある。
さらに、時間軸方向にｒ次元多次元時系列信号をランダムに定めた個数分選択するという選択も考えられる。ここでは、クラスタを単位とすることもあるが、クラスタ内を区分し、これらを定めた個数分、ランダムに選択する。
【００４２】
図１４に別の実施例を示す。アラーム信号を処理する部分１２を付加した。図１５（ａ）に示すようにアラーム信号の発生履歴を取得する。そして、そのヒストグラムを表示する。頻度が高い区間は、異常の度合いが高いことは容易に想像ができる。そこで、図１５（ｂ）に示すように、ヒストグラムの頻度も考慮して、図１に示す異常特定部１３にて、アラーム信号とはずれ値を組み合わせて異常の度合いや信頼度を付加したり、異常判定を行う。
【００４３】
図１６に別の実施例を示す。Wavelet（変換）解析を付加した例である。Wavelet解析信号付与部１４において、Ｍ次元の多次元時系列信号を対象に、図１７に示すWavelet解析を行い、これらの信号をＭ次元の多次元時系列信号に加える。Ｍ次元の多次元時系列信号と置換することも可能である。このような新規に加えられたり、置換されたりした多次元時系列信号を対象に、局所部分空間法などの識別器により、異常を検知する。なお、図１７（ａ）は後で説明する図２５におけるスケール１の信号に該当し、図１７（ｂ）は後で説明する図２５におけるスケール８の変動に該当し、図１７（ｃ）は図２５におけるスケール４の変動に該当し、図１７（ｄ）は図２５におけるスケール２の変動に該当する。
【００４４】
Wavelet解析は、多重解像度表現を与えるものである。図２５にWavelet解析を図解して示す。スケール１の信号がもとの信号である。これを順次隣と加算してスケール２の信号を作り、元の信号との差を算出して、スケール２の変動信号を生成する。順次、これを繰り返すと、最後に、スケール８の一定値の信号とその変動信号が得られ、結局、元の信号は、スケール２、４、８の各変動信号とスケール８の直流信号に分解できる。従って、このようなスケール２、４、８の各変動信号を、新たな特徴的信号とみなして、多次元時系列信号に追加して扱う。
【００４５】
パルスやインパルスなどの非定常信号では、フーリエ変換を行って得られる周波数スペクトルは全域に広がり、個々の信号について特徴を抽出することが困難である。時間に局在したスペクトルを得ることができるWavelet変換は、化学プロセスのような、パルスやインパルスなどを含む、非定常な信号を多く含むデータを対象にした場合に都合がよい。
【００４６】
また、一次遅れをもつ系では、時系列の状態だけでは、そのパターンを観測困難であるが、時間・周波数領域上では、識別可能な特徴が現れることがあり、Wavelet変換が有効な場合が多い。
【００４７】
なお、Wavelet解析の応用は、電気学会編集で、新誠一著の「ウェーブレット解析の産業応用」2005年、朝倉出版にくわしい。化学プラントの制御系診断、空調プラント制御における異常検知、セメントの焼成プロセスの異常監視、ガラス溶解炉制御など、多くの対象に適用されている。
【００４８】
本実施例における、従来技術との相違点は、Wavelet解析を多重解像度表現として扱い、もとの多次元時系列信号の情報を、Wavelet変換により顕在化した点にある。その上で、これらを多変量として扱うことにより、異常が微弱な段階から早期に検知できる。即ち、予兆として早期に検知可能となる。
【００４９】
Wavelet解析を分類に適用したほかの例を説明する。異常データが少なからずあり、教示できる場合の例である。使用する記号を以下に示す。
【００５０】
【数１】
2010092355.tif 000003
まず、以下のモデルを仮定する。
【００５１】
【数２】
2010092355.tif 000004
次に、判別式としてつぎの式を考える。
【００５２】
【数３】
2010092355.tif 000005
また、パラメータ推定に関しては、次式で定義される尤度を最大にする最尤推定法を用いることにする。
【００５３】
【数４】
2010092355.tif 000006
[判別式]
【００５４】
【数５】
2010092355.tif 000007
上式に基づいて、分類が可能になる。特に、上記はマルチクラスへの分類の例であるが、クラス数Ｋを２とすれば、２分類問題、すなわち異常検知の問題となる。そして、上記観測値にWavelet解析の結果を適用すれば、時間に局在した異常をも対象とした異常検知を行うことができる。このように、異常データが少なからずあり、教示できる場合には、統計的な意味で、より精度の高い異常検知を行うことができる。
【００５５】
図１８に別の実施例を示す。散布図・相関解析部１５を付加した例である。ｒ次元の多次元時系列信号を対象に、散布図・相関解析を行った例を図１９に示す。図１９の相互相関解析では、遅れのラグ（lag）を考慮している。通常、相互相関関数の最大値の位置をラグと呼ぶ。この定義に従えば、二つの現象に関する時間のずれは、相互相関関数のラグに等しい。
【００５６】
ラグの正負は、二つの現象のどちらが早く起こるかにより決まる。このような散布図解析や相互相関解析の結果は、時系列信号間の相関を表すものであるが、各クラスタの特徴付けにも活用でき、クラスタ間の類似性の判断指標となり得るものである。例えば、ラグの量の一致度により、クラスタ間の類似性を判断する。これにより、図１１に示した類似したクラスタのマージ等が可能となる。マージ後のデータを用いてモデル化する。なお、マージの方法は、他の手法であってもよい。
【００５７】
図２０に別の実施例を示す。時間・周波数解析部１６を付加した例である。ｒ次元の多次元時系列信号を対象に、時間・周波数解析を行った例を図２１に示す。時間・周波数解析、或いは散布図・相関解析を行い、これらの信号をＭ次元の多次元時系列信号に加える、またはＭ次元の多次元時系列信号と置換することも可能である。
【００５８】
図２２に別の実施例を示す。学習データのＤＢ１７とモデル化（１）１８を付加した例である。図２３に、その詳細を示す。モデル化（１）により、学習データを対象にそれを複数のモデルとしてモデル化し、観測データとの類似性を判断して該当モデルを当てはめ、観測データとの偏差を算出する。モデル化（２）は、図１と同様の部分であり、これから、観測データから得られるモデルとの偏差を算出する。
【００５９】
そして、モデル化（１）（２）のそれぞれの偏差から、状態変化を算出し、総合偏差を算出する。ここで、モデル化（１）（２）を均等に扱うことも可能であるが、重み付けをしてもよい。即ち、学習データを基本と考えれば、モデル（１）の重みを大きくし、観測データを基本と考えれば、モデル（２）の重みを大きくする。
【００６０】
図１２に示した表現に従えば、モデル（１）で構成した部分空間モデルをクラスタ間で比較すれば、そしてそれらが本来同一状態のクラスタならば、その状態変化を知ることができる。そして、観測データの部分空間モデルがそれから移動していれば、状態変化を読み取ることができる。状態変化が、部品交換など意図するものであれば、即ち、設計側が知っており、それによる変化を許容すべきであれば、モデル（１）の重みを小さくし、モデル（２）の重みを大きくする。状態変化が、意図しないものであれば、モデル（１）の重みを大きくする。例えば、パラメータαを、モデル（１）の重みとして使えば、
α×モデル（１）＋（１−α）×モデル（２）
として定式化できる。
【００６１】
モデル（１）の重みを、古いものほど小さくするという忘却形のものでもよい。この場合、最近のデータに基づくモデルが重視されることになる。
【００６２】
図２３において、物理モデル２２は、シミュレーションにより対象エンジンなどを模擬したモデルである。対象知識が十分にある場合は、対象エンジンなどを離散時間（非）線形状態空間モデル（状態方程式などで表現）で表現できるので、その中間値や出力などを推定可能となる。従って、この物理モデルに従えば、このモデルからの偏差に基づいて、異常検知を行うことも可能となる。
【００６３】
もちろん、物理モデルに従って、学習データのモデル（１）を修正することも可能である。或いは、逆に、学習データのモデル（１）に従って、物理モデルを修正することも可能である。物理モデルの変形として、過去の実績としての知見を、物理モデルとして組み込むことも可能である。アラームの発生や部品交換に伴うデータの遷移を、物理モデルに組み込むことも可能である。或いは、学習データ（個々のデータや重心位置など）を、アラームの発生や部品交換に伴い、移動させてもよい。
【００６４】
なお、図２３に対して、図１から図２２までに示したように、物理モデルに対して、統計モデルを主に使うのは、データを生み出すプロセスについての理解が少ない場合に、統計モデルが有効であることによる。距離や類似性は、データの生成過程が不明瞭でも、定義できる。画像が対象である場合も、画像生成過程が不明瞭なとき、統計モデルは有効である。対象に関する知識が、すこしでも利用できる場合には、物理モデル２２が活用できることになる。
【００６５】
図２６に、異常検知、及び異常検知後の診断を示す。同図において、設備からの時系列信号から、今までに説明した方法により、異常を検知する。同時に、保守のイベント（アラームや作業実績など）などの付帯情報を取り込み、異常を高感度に検知する。同図に示すように、早期に予兆として発見できれば、故障となって稼動停止となる前に、何らかの対策がうてることになる。そこで、この予兆に基づき、異常診断を行う。
【００６６】
異常診断は、予兆を内包しているセンサを特定する現象診断と、故障を引き起こす可能性のあるパーツを特定する原因診断に分けると考えやすい。異常検知部では、異常診断部に対して、異常の有無という信号のほか、特徴量に関する情報を出力する。異常診断部は、これらの情報をもとに診断を行う。
【００６７】
図２７に、ハードウェア構成を示す。異常検知を実行するプロセッサに、対象とするエンジンなどのセンサデータを入力し、欠損値の修復などを行って、データベースＤＢに格納する。プロセッサは、観測データ、ＤＢデータを用いて、異常検知を行い、異常検知信号を出力する。表示部では、実施例において述べたような各種表示を行う。トレンドを表示することも可能とする。
【００６８】
データベースＤＢは、熟練エンジニアがＤＢを操作できる。特に、異常事例や対策事例を教示でき、格納できる。（１）学習データ（正常）、（２）異常データ、（３）対策内容が、格納される。データベースＤＢを、熟練エンジニアが手を加えられる構造にすることにより、洗練された、有用なデータベースができあがることになる。データ操作は、学習データ（個々のデータや重心位置など）を、アラームの発生や部品交換に伴い、移動させることにより行う。また、取得データを追加することも可能である。異常データがあれば、データの移動に、一般化ベクトル量子化などの手法も適用できる。
【００６９】
上記実施例では、エンジンなどの設備を対象に説明したが、時系列信号のたぐいならば、対象にこだわらない。人体の測定データへも適用可能である。本実施例に従えば、状態の数、遷移の回数が多くても、対応可能である。
【００７０】
また、実施例にて説明した各機能、例えばクラスタリング、主成分分析や、Ｗａｖｅｌｅｔ解析などは、必ず実施するというものでもなく、対象の信号の性質に応じて、適宜実施すればよい。
【００７１】
クラスタリングも、時間軌跡のみならず、混合分布に対するEM（Expectation-Maximization ）アルゴリズムやｋ-meansクラスタリングなどを含め、データマイニング分野における手法が使えることは言うまでもない。得られたクラスタは、これを対象に識別器を適用してもよいが、クラスタをグルーピングして、これを対象に識別器を適用してもよい。最も簡単な例は、日々の観測データが属するクラスタと、それが属するクラスタ以外に分けることである（図１０の右側の特徴空間にて図示した着目データなる現データと、それより時間的に前の過去データに該当）。また、センサ信号（特徴量）の選択は、ラッパー法(例えば、backward stepwise selectionにより、全ての特徴量がある状態から、最も不要な特徴を一つずつ取り除く)など、既存の手法が適用可能である。
【００７２】
さらに、識別器は、いくつかの識別器を準備し、それらの多数決をとることも可能である。即ち、異なる識別器群を用いたアンサンブル（集団）学習が適用できる。図３０にその構成例を示す。例えば、識別器１は投影距離法、識別器２は局所部分空間法、識別器３は線形回帰法と言ったものである。事例データに基づくものならば、任意の識別器が適用可能である。
【００７３】
複数の識別器を使う理由は、識別器はそれぞれ異なる基準、異なる対象データ範囲（セグメント分けやその統合に依存）で、はずれ具合いを求めるため、その結果に微妙な違いが生まれるためである。このため、多数決をとって安定化するか、OR（はずれ値そのもの、すなわち多値の場合は最大値検出）論理で、いずれかの識別器で異常が検知されれば異常が発生したとして出力し、異常をのこらず検知しようとするか、AND（多値の場合は最小値検出）論理で、いずれの識別器でも同時に異常が検知されれば異常が発生したとして出力し、誤検知を最小限にするといった、上位基準で識別器を構成する。もちろん、アラーム信号、部品交換などの保全情報などの情報を加味して、上記統合を行うことも可能であることは言うまでもない。
【００７４】
識別器h1、h2、・・・をすべて同じ識別器とし、対象データ範囲（セグメント分けやその統合に依存）を変えて学習することも可能である。例えば、パターン認識の代表的手法であるバギングやブースティングなどの手法も適用できる。この手法の適用により、異常検知に関して、より高い正解率が確保できる。ここで、バギングは、N個のデータから重複を許してデータをK個取ってきて（復元抽出）、このK個を元に1個目の識別器h1を作り、またN個のデータから重複を許してデータをK個取ってきて、このK個(1個目の識別器と中身が違う)を元に2個目の学習器h2を作ることを続け、異なるデータからいくつかの識別器を作り、実際に判別器として使う時は多数決を取る方法である。
【００７５】
ブースティング（Adaboostと呼ぶ手法）は、N個のデータに均等な重み 1/N をまずは割り当て、1個目の識別器h1はN個全てのデータを用いて学習し、学習した後、データN個に関して正解率を調べ、それを元に信頼度β1(>0)を求める。1個目の識別器が正解したデータの重みをexp(-β1)倍して重みを減らし、正解できなかったデータの重みをexp(β1)倍して重みを増やす。
【００７６】
2個目の識別器h2はN個全てのデータを用いて重み付きの学習を行い、信頼度β2(>0)を求め、データの重みを更新する。2個とも正解したデータの重みは軽く、2個とも間違えたデータの重みは重くなる。以降、これを繰返してM個の識別器を作り、実際に判別器として使う時は信頼度付き多数決を取るものである。これらの手法を、クラスタ群を対象に適用することにより、性能向上が期待できる。
【００７７】
図３１に、図３０に示した識別器を含んだ、異常検知全体の構成例の一例を示す。軌跡クラスタリング、特徴選択などを経て、アンサンブル学習を行い、高い識別率を達成するものである。識別器出力の統合の方法は上述した通りであるが、どのクラスタに、どの識別器を適用するかという組み合わせは、いくつか存在する。例えば、観測データとは異なるクラスタに対して局所部分空間法を適用して、異なるクラスタからのはずれ具合いを把握し（推定値も算出する）、観測データと同じクラスタに対しては回帰分析法を適用して、自クラスタからのはずれ具合いを把握する。
【００７８】
そして、それらの識別器出力を統合して、異常判定を行うことができる。他のクラスタからのはずれ具合いを、投影距離法や回帰分析法により行うことも可能である。自クラスタからのはずれ具合いを投影距離法により行うこともできる。クラスタは、アラーム信号が活用できる場合には、アラーム信号の重症度レベルに応じて、重症アラーム信号が付加されていないクラスタを対象にすることもできる。
【００７９】
クラスタ間の類似性を判断し、類似クラスタを統合して、これを対象にすることもできる。識別器出力の統合は、はずれ値の加算、最大・最小、OR/ANDなどのスカラー変換処理でもよいし、識別器の出力をベクトル的に、多次元として扱うこともできる。もちろん、識別器出力のスケールは、極力一致させることとする。
【００８０】
上述したクラスタとの関連の持たせ方に関して、さらに、他クラスタを対象にして第１報の異常検知を行い、自クラスタのデータが収集された時点で自クラスタを対象にして第２報の異常検知を行ってもよい。このようにして、顧客への注意喚起を促すことが可能になる。このように、本実施例は、対象クラスタ群との関係において、信号の振舞い、behaviorに、より着目した実施例と言える。
【００８１】
上述したいくつかの実施例に関する総合的効果をさらに補足する。えば、発電設備を所有している会社では、機器の保守費用削減を希望しており、保証期間中に機器を点検、部品交換を実施している。これは時間ベースの設備保全と言われている。しかし、最近は機器の状態を見て、部品交換を実施する状態ベースの保全に移行しつつある。状態保全を実施するには、機器の正常・異常データを収集する必要があり、このデータの量、質が状態保全の品質を決めてしまう。しかし、異常データの収集は、まれなケースも多く、大型の設備になるほど、異常データを収集することは困難である。従って、正常データから、はずれ値を検出することが重要となる。上述したいくつかの実施例によれば、
正常データから、異常を検知できる、
データ収集が不完全でも精度の高い異常検知が可能となる、
異常データが包含されていても、この影響を許容できる、
といった直接的効果に加え、
ユーザにとって、現象を理解しやすい、
エンジニアの知識を活用できる
物理モデルも併用できる、
と言った副次的な効果がある。
【産業上の利用可能性】
【００８２】
プラント、設備の異常検知として利用することが出来る。
【図面の簡単な説明】
【００８３】
【図１】本発明の第１の実施例を示すチャートである。
【図２】多次元時系列信号の例である。
【図３】相関行列の例である。
【図４】軌跡分割クラスタリングの適用例である。
【図５】軌跡分割クラスタリングの適用例である。
【図６】軌跡分割クラスタリングの適用例である。
【図７】軌跡分割クラスタリング実施時のモデルとの偏差の例である。
【図８】軌跡分割クラスタリング未実施時のモデルとの偏差の例である。
【図９】特徴空間にてデータ表示をした一例である。
【図１０】特徴空間にてデータ表示をした他の例である。
【図１１】特徴空間にてデータ表示をしたさらに他の例である。
【図１２】特徴空間にてデータ表示をあいたさらに他の例である。
【図１３】本発明の第２の実施例を示すチャートである。
【図１４】本発明の第３の実施例を示すチャートである。
【図１５】アラーム信号のヒストグラム例である。
【図１６】本発明の第４の実施例を示すチャートである。
【図１７】Ｗａｖｅｌｅｔ（変換）解析の例である。
【図１８】本発明の第５の実施例を示すチャートである。
【図１９】散布図解析および相互相関解析の例である。
【図２０】本発明の第６の実施例を示すチャートである。
【図２１】時間・周波数解析の例である。
【図２２】本発明の第７の実施例を示すチャートである。
【図２３】本発明の第７の実施例の詳細を支援すチャートである。
【図２４】部分空間法の一例である。
【図２５】Ｗａｖｅｌｅｔ変換の説明図である。
【図２６】異常診断との関係を説明する図である。
【図２７】ハードウェア構成図
【図２８】局所部分空間法の適用例である。
【図２９】投影距離法、局所部分空間法の適用例である。
【図３０】複数識別器の統合による識別例である。
【図３１】複数識別器の統合による異常検知例である。
【符号の説明】
【００８４】
1・・・多次元信号取得部、
2・・・欠損値の修正・削除部、
3・・・状態データ・知識データベース、
4・・・相関解析による無効信号の削除部、
6・・・軌跡分割クラスタリング、
7・・・アラーム信号／保全情報、
8・・・各クラスタを対象にモデル化部、
9・・・モデルからの偏差算出部、
10・・・はずれ値検出部、
11・・・各クラスタを特徴選択でモデル化部、
12・・・アラーム信号などの一定区間累積ヒストグラム、
13・・・異常特定部、
14・・・Ｗａｖｅｌｅｔ（変換）解析部、
15・・・各クラスタ軌跡散布図・相関解析部、
16・・・各クラスタ毎時間・周波数解析部、
17・・・学習データ、
18・・・モデル化（１）部、
19・・・プロセッサ、
20・・・表示器、
21・・・データベース
22・・・物理モデル。
【図１】
2010092355.tif 000008
【図２】
2010092355.tif 000009
【図３】
2010092355.tif 000010
【図４】
2010092355.tif 000011
【図５】
2010092355.tif 000012
【図６】
2010092355.tif 000013
【図７】
2010092355.tif 000014
【図８】
2010092355.tif 000015
【図９】
2010092355.tif 000016
【図１０】
2010092355.tif 000017
【図１１】
2010092355.tif 000018
【図１２】
2010092355.tif 000019
【図１３】
2010092355.tif 000020
【図１４】
2010092355.tif 000021
【図１５】
2010092355.tif 000022
【図１６】
2010092355.tif 000023
【図１７】
2010092355.tif 000024
【図１８】
2010092355.tif 000025
【図１９】
2010092355.tif 000026
【図２０】
2010092355.tif 000027
【図２１】
2010092355.tif 000028
【図２２】
2010092355.tif 000029
【図２３】
2010092355.tif 000030
【図２４】
2010092355.tif 000031
【図２５】
2010092355.tif 000032
【図２６】
2010092355.tif 000033
【図２７】
2010092355.tif 000034
【図２８】
2010092355.tif 000035
【図２９】
2010092355.tif 000036
【図３０】
2010092355.tif 000037
【図３１】
2010092355.tif 000038