(57)【要約】  設備における事例ベースの異常予兆検知においては、多様な正常状態を網羅的に集める必要があるが、停止時間の長さなど運転パターンによって状態が大きく異なるために学習データ不足による誤報が発生するという問題や、学習データ不足が解消されるまで学習データ期間を長くした場合にメモリ容量と計算時間が増大して実行困難になるという問題を解決するために、設備の出力する時系列信号に基づいて設備の状態を監視する方法において、時系列信号に基づいて一定期間毎に運転パターンラベルを付与し、一定期間毎の運転パターンラベルに基づいて学習データを選定し、この選定した学習データに基づき正常モデルを作成し、時系列信号と正常モデルに基づき異常測度を算出し、この算出した異常測度に基づき設備の状態が異常か正常かの識別を行うようにした。
2013030984.tif 000002
【特許請求の範囲】
【請求項１】
設備の出力する時系列信号に基づいて設備の状態を監視する方法であって、前記時系列信号に基づいて一定期間毎に運転パターンラベルを付与し、前記一定期間毎の運転パターンラベルに基づいて学習データを選定し、該選定した学習データに基づき正常モデルを作成し、前記時系列信号と前記正常モデルに基づき異常測度を算出し、該算出した異常測度に基づき前記設備の状態が異常か正常かの識別を行うことを特徴とする設備状態監視方法。
【請求項２】
前記運転パターンラベルとして、前記設備の定常オフの状態、前記設備の定常オンの状態、前記設備の起動動作の状態、前記設備の停止動作の状態に対して異なる運転パターンラベルを付与することを特徴とする請求項１記載の設備状態監視方法。
【請求項３】
前記正常モデルを、前記運転パターンラベルごとに作成することを特徴とする請求項１記載の設備状態監視方法。
【請求項４】
設備の出力する時系列信号に基づいて設備の状態を監視する方法であって、前記時系列信号に一定期間毎に有限個にカテゴライズされた運転パターンラベルを付与する運転パターンラベル付与工程と、該運転パターンラベル付与工程により運転パターンラベルが付与された時系列信号をデータとして蓄積するデータ蓄積工程と、該データ蓄積工程で蓄積されたデータの中から前記時系列信号に付与された運転パターンラベルに基づいて所定数のデータを選択して学習データとする学習データ選択工程と、該学習データ選択工程で選択した学習データを用いて正常モデルを作成する正常モデル作成工程と、該正常モデル作成工程で作成した正常モデルとの比較に基づいて前記時系列信号の異常測度を算出する異常測度算出工程と、該異常測度算出工程で算出した異常測度に基づき異常を識別する異常識別工程とを含むことを特徴とする設備状態監視方法。
【請求項５】
前記学習データ選択工程において、前記蓄積されたデータの中から、前記時系列信号に付与された運転パターンラベルと同じ運転パターンラベルあるいは状態の近い運転パターンラベルの所定数のデータを選択して学習データとすることを特徴とする請求項４記載の設備状態監視方法。
【請求項６】
前記運転パターンラベル付与工程は、一定期間毎に前記設備の定常オフの状態、前記設備の定常オンの状態、前記設備の起動動作の状態、前記設備の停止動作の状態を含む複数の運転パターン特徴を算出する工程と、前記複数の運転パターン特徴の組合せに基づき運転パターンラベルを付与するラベル付与工程からなることを特徴とする請求項４記載の設備状態監視方法。
【請求項７】
さらに前記時系列信号の一定期間毎にマクロな変動を表す特徴を算出するマクロ特徴算出工程を備え、前記学習データ選択工程において、前記時系列信号に付与された運転パターンラベルと前記算出されたマクロ特徴に基づいて前記蓄積されたデータの中から所定数のデータを選択して学習データとすることを特徴とする請求項４記載の設備状態監視方法。
【請求項８】
前記学習データ選択工程において、前記蓄積されたデータの中から、前記時系列信号に付与された運転パターンラベルと同じ運転パターンラベルあるいは状態の近い運転パターンラベルかつ前記時系列信号に関連するマクロ特徴が類似している所定数のデータを選択して学習データとすることを特徴とする請求項７記載の設備状態監視方法。
【請求項９】
前記学習データ選択工程において、前記時系列信号に基づき顕著な状態変化の有無を判断し、前記顕著な状態変化がない場合は、前記蓄積されたデータの中から、前記時系列信号に付与された運転パターンラベルと同じ運転パターンラベルあるいは状態の近い運転パターンラベルの所定数のデータを選択して学習データとし、前記顕著な状態変化がある場合は、前記蓄積されたデータの中から、前記時系列信号に付与された運転パターンラベルと同じ運転パターンラベルあるいは状態の近い運転パターンラベルかつ前記時系列信号に関連するマクロ特徴が類似している所定数のデータを選択して学習データとすることを特徴とする請求項７記載の設備状態監視方法。
【請求項１０】
前記マクロ特徴が、センサ信号の期間全体に亘る平均、分散、最大値、最小値、又は、期間中の定常オンの時刻の平均、分散又は期間中の定常オフの時刻の平均、分散のうちの少なくとも何れか一つを含むことを特徴とする請求項７記載の設備状態監視方法。
【請求項１１】
前記顕著な状態変化の有無の判定は、前記マクロ特徴の変化に基づいて行うことを特徴とする請求項９記載の設備状態監視方法。
【請求項１２】
設備から出力される時系列信号を入力して解析するセンサ信号解析部と、該センサ信号解析部で解析した結果と前記時系列信号とを受けて前記設備の異常を診断する異常診断部と、前記センサ信号解析部及び前記異常診断部と接続してデータの入出力を行う入出力部とを備えた設備状態を監視する装置であって、
前記センサ信号解析部は、前記設備から出力された時系列信号に基づいて一定期間毎に運転パターンラベルを付与する運転パターンラベル付与手段と、該運転パターンラベル付与手段で付与された一定期間毎の運転パターンラベルに基づいて学習データを選定する学習データ作成手段と、該学習データ作成手段で作成された学習データに基づき正常モデルを作成する正常モデル作成手段と、該正常モデル作成手段で作成した正常モデルに基づいて前記設備から出力された時系列信号の異常測度を算出する異常測度算出手段と、該異常測度算出手段で算出した異常測度に基づき前記設備の状態が異常か正常かの識別を行う異常識別手段とを有することを特徴とする設備状態監視装置。
【請求項１３】
前記運転パターンラベル付与部で運転パターンラベルを付与した時系列信号を記憶するデータベース手段を更に備え、前記学習データ作成手段は、前記運転パターンラベル付与手段により前記時系列信号に付与された前記運転パターンラベルに基づいて前記データベース手段に蓄積されたデータの中から同じ運転パターンラベルあるいは状態の近い運転パターンラベルのデータを所定数選択して学習データとすることを特徴とする請求項１２記載の設備状態監視装置。
【請求項１４】
設備から出力される時系列信号を入力して設備状態を監視する装置であって、前記時系列信号を入力して外入力した時系列信号に一定期間毎に有限個にカテゴライズされた運転パターンラベルを付与する運転パターンラベル付与手段と、該運転パターンラベル付与手段により運転パターンラベルを付与された時系列信号を蓄積するデータ蓄積手段と、前記運転パターンラベル付与手段により前記時系列信号に付与された前記ラベルに基づいて前記データ蓄積手段に蓄積されたデータの中から同じ運転パターンラベルあるいは状態の近い運転パターンラベルのデータを所定数選択して学習データとする学習データ選択手段と、該学習データ選択手段で選択した学習データを用いて正常モデルを作成する正常モデル作成手段と、前記運転パターンラベル付与手段により運転パターンラベルを付与された時系列信号を前記正常モデル作成手段で作成した正常モデルと比較して前記時系列信号の異常測度を算出する異常測度算出手段と、該異常測度算出手段で算出した異常測度に基づき前記運転パターンラベルを付与された時系列信号の異常識別をする異常識別手段とを備えたことを特徴とする設備状態監視装置。  
【請求項１５】
前記運転パターンラベル付与手段に前記運転パターンラベルを付与するための条件を入力する入力手段を更に備えたことを特徴とする請求項１４記載の設備状態監視装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プラントや設備などの出力する多次元時系列データをもとに異常を早期に検知する設備状態監視方法およびその装置に関する。
【背景技術】
【０００２】
電力会社では、ガスタービンの廃熱などを利用して地域暖房用温水を供給したり、工場向けに高圧蒸気や低圧蒸気を供給したりしている。石油化学会社では、ガスタービンなどを電源設備として運転している。このようにガスタービンなどを用いた各種プラントや設備において、設備の不具合あるいはその兆候を検知する予防保全は、社会へのダメージを最小限に抑えるためにも極めて重要である。【０００３】
ガスタービンや蒸気タービンのみならず、水力発電所での水車、原子力発電所の原子炉、風力発電所の風車、航空機や重機のエンジン、鉄道車両や軌道、エスカレータ、エレベータ、機器・部品レベルでも、搭載電池の劣化・寿命など、上記のような予防保全を必要とする設備は枚挙に暇がない。【０００４】
このため、対象設備やプラントに複数のセンサを取り付け、センサ毎の監視基準に従って正常か異常かを判断することが行われている。米国特許第６,９５２,６６２号明細書（特許文献１）や米国特許第６,９７５,９６２号明細書（特許文献２）には、おもにエンジンを対象とした異常検知方法が開示されている。これは、過去のデータ例えば時系列センサ信号をデータベースとしてもっておき、観測データと過去の学習データとの類似度を独自の方法で計算し、類似度の高いデータの線形結合により推定値を算出して、推定値と観測データのはずれ度合いを出力するものである。また、特開２０１０−１９１５５６号公報（特許文献３）には、過去の正常データから、観測データに類似したコンパクトな学習データセットを抽出し、抽出した学習データを部分空間でモデル化し、観測データと部分空間の距離に基づいて異常を検知する異常検知方法が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】米国特許第６,９５２,６６２号明細書
【特許文献２】米国特許第６,９７５,９６２号明細書
【特許文献３】特開２０１０−１９１５５６号公報
【非特許文献】
【０００６】
【非特許文献１】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
【発明の概要】
【発明が解決しようとする課題】
【０００７】
特許文献１や特許文献２に記載の方法によれば、学習データとして正常時のデータを与えることにより、学習にない観測データが観察されると、これらを異常として検出することができる。しかし、その異常検知性能は学習データの質に大きく左右される。すなわち、学習データにない観測データは正常であっても異常と判断され、一方、学習データに誤って異常が混入されていると異常であっても正常と判断される可能性が生じる。そのため、学習データは正常な状態のみを正確にかつ網羅的に収集しなければ、異常検知の信頼性が著しく低下してしまうが、多様な正常状態を有する設備に対してこのような学習データ収集を行うことは非常に負荷が高く、これを人手で行うことは実質不可能である場合が多い。さらに、質の良い学習データを収集可能であっても、計算負荷の高い方法であるため、実現可能な計算時間で処理するのに許容されるデータ量が小さく、その結果網羅性が確保できなくなる場合が多い。【０００８】
特許文献３に記載の方法によれば、同様に学習データの質が重要であるが、観測データに類似したコンパクトな学習データセットを抽出するため、計算負荷を低くすることができ、網羅性を確保しやすくなる。しかし、一定期間のデータを学習データとした場合、運転と停止のパターンが複雑な設備を対象とすると、学習データの不足により類似データが見つからずに正常を異常と判定する誤報となってしまうケースが発生する。学習データ不足を解消されるまで学習データの期間を長くするとメモリ容量と計算時間が増大し、実行困難になる。また、運転パターンが変化しない場合でも、保守作業などにより状態が大きく変化した場合、一定期間のデータを学習データとしていると、その状態変化を異常として検知し誤報となってしまうケースが発生する。【０００９】
また、設備における従来の事例ベースの異常予兆検知においては、多様な正常状態を網羅的に集める必要があるが、停止時間の長さなど運転パターンによって状態が大きく異なるため、学習データ不足により誤報が発生するという課題がある。また、学習データ不足が解消されるまで学習データ期間を長くすると、メモリ容量と計算時間が増大し、実行困難になるという課題もある。【００１０】
そこで、本発明の目的は、上記課題を解決し、プラントや設備あるいは製造装置，計測装置（以下、これらを総称して設備と記す）などで運転と停止のパターンが複雑な設備に対しても、計算負荷を低く抑えたまま異常を高感度に検知することが可能な異常検知方法を備えた設備状態監視方法およびシステムを提供することにある。さらに、運転と停止のパターンの変化のほか、保守作業による状態の変化があった場合にも、適切な学習データを選択して高感度な異常検知を実現可能な設備状態監視方法およびその装置を提供することにある。
【課題を解決するための手段】
【００１１】
上記目的を達成するために、本発明では、設備の出力する時系列信号に基づいて設備の状態を監視する方法において、時系列信号に基づいて一定期間毎に運転パターンラベルを付与し、一定期間毎の運転パターンラベルに基づいて学習データを選定し、この選定した学習データに基づき正常モデルを作成し、時系列信号と正常モデルに基づき異常測度を算出し、この算出した異常測度に基づき設備の状態が異常か正常かの識別を行うようにした。【００１２】
また、本発明では、設備または装置の出力する時系列信号に基づいて設備の状態を監視する方法において、時系列信号に一定期間毎に有限個にカテゴライズされた運転パターンラベルを付与する運転パターンラベル付与工程と、この運転パターンラベル付与工程により運転パターンラベルが付与された時系列信号をデータとして蓄積するデータ蓄積工程と、このデータ蓄積工程で蓄積されたデータの中から時系列信号に付与された運転パターンラベルに基づいて所定数のデータを選択して学習データとする学習データ選択工程と、この学習データ選択工程で選択した学習データを用いて正常モデルを作成する正常モデル作成工程と、この正常モデル作成工程で作成した正常モデルとの比較に基づいて時系列信号の異常測度を算出する異常測度算出工程と、この異常測度算出工程で算出した異常測度に基づき異常を識別する異常識別工程とを含むようにした。【００１３】
また、本発明では、設備または装置から出力される時系列信号を入力して解析するセンサ信号解析部と、このセンサ信号解析部で解析した結果と時系列信号とを受けて設備の異常を診断する異常診断部と、センサ信号解析部及び異常診断部と接続してデータの入出力を行う入出力部とを備えた設備状態を監視する装置において、センサ信号解析部は、設備または装置から出力された時系列信号に基づいて一定期間毎に運転パターンラベルを付与する運転パターンラベル付与手段と、この運転パターンラベル付与手段で付与された一定期間毎の運転パターンラベルに基づいて学習データを選定する学習データ作成手段と、この学習データ作成手段で作成された学習データに基づき正常モデルを作成する正常モデル作成手段と、この正常モデル作成手段で作成した正常モデルに基づいて設備または装置から出力された時系列信号の異常測度を算出する異常測度算出手段と、この異常測度算出手段で算出した異常測度に基づき設備の状態が異常か正常かの識別を行う異常識別手段とを有するように構成した。【００１４】
更に、本発明では、設備または装置から出力される時系列信号を入力して設備状態を監視する装置において、時系列信号を入力して外入力した時系列信号に一定期間毎に有限個にカテゴライズされた運転パターンラベルを付与する運転パターンラベル付与手段と、この運転パターンラベル付与手段により運転パターンラベルを付与された時系列信号を蓄積するデータ蓄積手段と、運転パターンラベル付与手段により時系列信号に付与されたラベルに基づいてデータ蓄積手段に蓄積されたデータの中から同じ運転パターンラベルあるいは状態の近い運転パターンラベルのデータを所定数選択して学習データとする学習データ選択手段と、この学習データ選択手段で選択した学習データを用いて正常モデルを作成する正常モデル作成手段と、運転パターンラベル付与手段により運転パターンラベルを付与された時系列信号を正常モデル作成手段で作成した正常モデルと比較して時系列信号の異常測度を算出する異常測度算出手段と、この異常測度算出手段で算出した異常測度に基づき運転パターンラベルを付与された時系列信号の異常識別をする異常識別手段とを備えて構成した。
【発明の効果】
【００１５】
本発明によれば、一定期間毎にデータに運転パターンラベルを付与しておき、異常検知時には同じラベルあるいは状態の近いラベルの期間のデータを所定数集めて学習データとするため、計算負荷を増やさずに高精度な正常モデルを作成することができ、正常を異常と判定する誤報を低減できる。【００１６】
さらに本発明によれば、センサデータのマクロな変動を表すマクロ特徴を算出し、同じラベルあるいは状態の近いラベルかつマクロ特徴の類似した期間のデータを所定数集めて学習データとするため、保守作業などにより状態が大きく変化した場合にも、類似した学習データを選択して高精度な正常モデルを作成することが可能となり、正常を異常と判定する誤報を低減できる。【００１７】
以上により、ガスタービンや蒸気タービンなどの設備のみならず、水力発電所での水車、原子力発電所の原子炉、風力発電所の風車、航空機や重機のエンジン、鉄道車両や軌道、エスカレータ、エレベータ、そして機器・部品レベルでは、搭載電池の劣化・寿命など、種々の設備・部品において高感度な異常検知および容易な異常説明を両立するシステムが実現できる。
【図面の簡単な説明】
【００１８】
【図１】本発明の第１の実施例における設備状態監視システムの概略の構成を示すブロック図である。
【図２Ａ】本発明の第１の実施例におけるセンサ信号の例を示す信号リストである。
【図２Ｂ】本発明の第１の実施例におけるイベント信号の例を示す信号リストである。
【図３Ａ】本発明の第１の実施例における運転パターンラベル付与の処理の流れを示すフロー図である。
【図３Ｂ】本発明の第１の実施例における設備の可動状態を分割して４種のモードの何れかに分類した状態を示すモード分割の模式図である。
【図４Ａ】本発明の第１の実施例における運転パターンラベル付与の概念を説明する図である。
【図４Ｂ】本発明の第１の実施例におけるイベント信号から算出されたデータと運転パターンラベルの関係の例を示す図である。
【図５Ａ】本発明の第１の実施例における運転パターンラベルの管理データの例を示すデータリストである。
【図５Ｂ】本発明の第１の実施例におけるセンサ信号データの例を示すデータリストである。
【図５Ｃ】本発明の第１の実施例におけるイベント信号データの例を示すデータリストである。
【図５Ｄ】本発明の第１の実施例における運転パターンラベルの管理データの別の例を示すデータリストである
【図６Ａ】本発明の第１の実施例におけるラベル付与の条件を決めるための情報入力画面の例である。
【図６Ｂ】本発明の第１の実施例におけるラベル付与の条件を決めるための詳細情報入力画面の例である。
【図６Ｃ】本発明の第１の実施例における参照情報とラベルの関係を確認する画面の例である。
【図６Ｄ】本発明の第１の実施例におけるラベルとラベル付与参照情報とセンサ信号をグラフ表示する画面の例である。
【図７】本発明の第１の実施例における学習データ選択の処理の流れを示すフロー図である。
【図８】本発明の第１の実施例における正常モデル生成の処理の流れを示すフロー図である。
【図９】投影距離法を説明するグラフである。
【図１０】局所部分空間法を説明するアフィン部分空間図である。
【図１１】本発明の第２の実施例における設備状態監視システムの概略の構成を示すブロック図である。
【図１２】本発明の第２の実施例におけるラベル付与の条件を決めるための情報入力画面の例である。
【図１３】本発明の第２の実施例における学習データ選択の処理の流れを示すフロー図である。
【発明を実施するための形態】
【００１９】
以下、図面を用いて本発明の内容を詳細に説明する。
【実施例１】
【００２０】
図1に、本発明の設備状態監視方法を実現するシステムの一構成例を示す。本システムは、センサ信号解析部１００と異常診断部１２０、入出力部１３０とを備えて構成されている。【００２１】
本システムのセンサ信号解析部１００は、設備１０１から出力されるセンサ信号１０２とイベント信号１０３に基づき一定期間毎に運転パターンラベルを付与する運転パターンラベル付与部１０４と運転パターンラベルを付与されたセンサ信号１０２とイベント信号１０３を蓄積するデータベース１０５、センサ信号１０２とイベント信号１０３を受けて、運転パターンラベルに基づいて蓄積されたセンサ信号１０２とイベント信号１０３から学習データを選択する学習データ選択部１０６、選択した学習データを用いて正常モデルを作成する正常モデル作成部１０７、正常モデルとセンサ信号に基づき異常測度を算出する異常測度算出部１０８、算出された異常測度に基づき異常を検知する異常識別部１０９とを備えて構成される。【００２２】
異常診断部１２０は、センサ信号１０２とイベント信号１０３及び異常識別部１０９からの出力を受けて、設備１０１の異常を診断する。【００２３】
入出力部１３０は、運転パターンラベル付与部１０４、異常識別部１０９、異常診断部１２０などと接続して診断条件の入力や診断結果の出力を行う。【００２４】
状態監視の対象とする設備１０１は、ガスタービンや蒸気タービンなどの設備やプラントである。設備１０１は、その状態を表すセンサ信号１０２とイベント信号１０３を出力する。【００２５】
センサ信号１０２の例を図２Ａに示す。時刻と１０１に設けられた複数のセンサ値のデータからなり、一定間隔毎に取得される時系列信号である。センサの種類は、数百から数千と言った数になる場合もあり、例えば、シリンダ、オイル、冷却水などの温度、オイルや冷却水の圧力、軸の回転速度、室温、運転時間などである。出力や状態を表すのみならず、何かをある値に制御するための制御信号の場合もある。【００２６】
イベント信号１０３の例を図２Ｂに示す。不定期に出力される設備の操作・故障・警告を表す信号であり、時刻と操作・故障・警告を表すユニークなコードとメッセージ文字列からなる。【００２７】
次に、センサ信号解析部の構成及び動作について、図１から図１０を用いて説明する。まず、運転パターンラベル付与部１０４における処理の流れを図３Ａおよび図３Ｂを用いて説明する。まず、イベント信号１０３を入力し（Ｓ３０１）、所定の文字列またはコードの検索により起動シーケンスと停止シーケンスの切り出しを行う（Ｓ３０２）。その結果をもとに、図３Ｂに示すように、停止シーケンスの終了時刻から起動シーケンスの開始時刻までの「定常OFF」モード３１１、起動シーケンス中の「起動」モード３１２、起動シーケンスの終了時刻から停止シーケンスの開始時刻までの「定常ON」モード３１３、停止シーケンス中の「停止」モード３１４の４つの稼動状態に分割する（Ｓ３０３）。以下の説明ではこの分割をモード分割と呼び、稼動状態の種類をモードと呼ぶこととする。【００２８】
シーケンス切り出しのためには、予めシーケンスの開始イベントおよび終了イベントを指定しておき、イベント信号１０３の先頭から最後まで以下の要領でスキャンしながら切り出していく。
（１）シーケンスの途中でない場合は、開始イベントを探索する。見つかったらシーケンスの開始とする。
（２）シーケンスの途中の場合は、終了イベントを探索する。見つかったらシーケンスの終了とする。
ここで終了イベントとは、指定の終了イベントのほか、故障、警告、指定の開始イベントとする。指定の終了イベント以外で終了した場合は、異常終了として記録しておく。【００２９】
次に、予め定めた期間毎、例えば１日毎に状態遷移数をカウントする（Ｓ３０４）。状態遷移とは定常ONと定常OFFの間の遷移のことを意味する。つまり起動シーケンスと停止シーケンスの回数を数えればよい。次に、予め定めた期間の初めの状態がONかOFFのどちらであるかをチェックしその時点までのその状態の継続時間を算出する（Ｓ３０５）。その時点より以前で一番近い起動シーケンスまたは停止シーケンスを探し、それが起動シーケンスであれば状態はON、停止シーケンスであれば状態はOFFである。継続時間はそのシーケンス終了からの時間を計算すればよい。次に算出した状態遷移回数と予め定めた期間の初めの状態とその継続時間をもとに運転パターンラベルを付与する（Ｓ３０６）。【００３０】
運転パターンラベルの付与方法について図４Ａおよび図４Ｂを用いて説明する。図４Ａは設備のON/OFFの情報と１つのセンサ信号を時系列グラフとして表した図である。ONかOFFかによってセンサ値が変化することはもちろんのこと、状態の継続時間、状態遷移のパターンによってもセンサ値が大きく変化することが分かる。運転パターンラベル付与の目的は、設備状態が類似するものを学習データとして選択することであるため、例えば、図４Ａに示した実線、点線、破線、一点鎖線で囲んだデータにそれぞれ同じラベルを付与する。そのために、Ｓ３０４、Ｓ３０５で算出した状態遷移回数と予め定めた期間の初めの状態とその継続時間を利用する。【００３１】
図４Ｂに、それらの情報と運転パターンラベルの関係の例を示す。状態遷移回数は0回から1,2，…と増えるに従い異なるラベルが付与されるようにする。ただし、ラベルの種類を有限個にするため、所定の回数以上は同じものとして扱う。この例では５回以上をひとつにまとめている。また、継続時間は連続データであるため、適当な区切りを設けて分類し、それぞれが異なるラベルが付与されるようにする。さらに、初めの状態がONかOFFかによって異なるラベルが付与されるようにする。状態遷移の回数および初めの状態の継続時間の対応する欄を見、初めの状態がONならば左の番号（図４Bの例では、１〜３６）のラベル、OFFならば右の番号（図４Bの例では、３７〜７２）のラベルを付与する。【００３２】
なお、横軸、縦軸の分け方はより細かくしても粗くしてもよいし、全ての情報を用いなくてもよい。また、さらに他の情報を加えて分類に用いてもよい。また、全ての欄に異なるラベルを割り当てなくてもよい。また、この例では初めの状態がONかOFFかの２種類であるが、OFFの中でも完全停止とアイドリングの状態を区別したり、ONの中でも異なる運転モードを区別したりすることも考えられる。設備の状態が類似するか否かに関連のある情報を用いることが重要であり、ここであげた項目に限定するものではない。例えば保守作業の有無や警告の有無をラベルに反映することが考えられる。【００３３】
また、本実施例では、イベント信号を用いてモード分割する例を示したが、センサ信号に設備のON／OFFあるいはその他の運転状態を表す信号が含まれる場合も考えられる。その場合は、その信号を用いて状態遷移回数と予め定めた期間の初めの状態とその継続時間あるいは設備の状態が類似するか否かに関連のある他の情報を算出してもよい。【００３４】
運転パターンラベルを付与されたセンサ信号１０２とイベント信号１０３はデータベース１０５に蓄積される。データベース１０５における、データの管理方法の例を図５Ａないし図５Dを用いて説明する。【００３５】
図５Ａは、運転パターンラベルの管理データの例である。この例では、前述の予め定めた期間を１日の0時から24時までとしており、日付、サイト（設置場所を表すコード）、号機（同じ場所にある複数の設備を区別するＩＤ）毎に１個のデータを持つ。このデータは日付５０１、サイト５０２、号機５０３の他に運転パターンラベル５０４、センサ信号データおよびイベント信号データへのキーコード５０５と運転パターンラベルを付与するために用いた初期状態５０６、初期状態の継続時間５０７、状態遷移回数５０８と警告の有無情報５０９からなる。５０６から５０８の項目は、運転パターンラベルの付与方法に応じて変化する。また、システムに接続された設備が１台の場合はサイト５０２および号機５０３の項目は不要である。【００３６】
図５Ｂは、センサ信号データである。設備１０１からの出力１０２と同様の情報を持つが、図５Ａに示す管理データのセンサ信号データおよびイベント信号データへのキーコード５０５に紐付けられたキーコード５１１が付加されている。キーコードが決まれば日付も一意に決まるため、日付の情報は含まなくてもよい。時刻情報５１２は、イベント信号との対応付けのため必要である。ただし、一定間隔データであることが保証されている場合は、時刻情報の替わりに取得順の通し番号を用いてもよい。図５Ｃは、イベント信号データである。設備１０１からの出力１０３と同様のコード５２３とメッセージ５２４と、管理データのキーコード５０５に紐付けられたキーコード５２１と時刻情報５２２からなる。キーコードが決まれば日付も一意に決まるため、日付の情報は含まなくてもよい。【００３７】
なお、運転パターンラベル管理の期間は０時を開始とするとは限らず、任意の時刻で区切ってよい。また、期間の長さを１日とするとも限らない。ラベル付与が煩雑になるのを防ぐため、主要な繰り返しの周期にそろえるとよい。例えばガスエンジンによる発電設備などは、設備が設置された施設の操業周期に合わせて運転されるため１日とするのが適切であるが、１２時間毎、８時間毎とすることも考えられる。周期を短くすることにより、状態数は少なくなって同じ状態のデータを集めやすくなると思われる。ただし、学習するセンサデータ数を一定の量に保とうとすると、学習するセンサデータを選択する際に探索する管理データ数および集めるべき管理データ数が増えるため、短くしすぎてもいけない。逆に周期を長くすることも可能であるが、状態数が増加し同じ状態のデータを集めるのが困難になる可能性がある点に留意する必要がある。【００３８】
さらに、運転パターンラベル管理の期間を定時刻のみではなく、状態遷移のタイミングでも区切ることも考えられる。その場合の運転パターンラベル管理データの例を図５Ｄに示す。このデータは、日付５３１、サイト５３２、号機５３３、運転パターンラベル５３４、センサ信号データおよびイベント信号データへのキーコード５３５と警告の有無情報５３９を図５Ａに示す例と同様に含む。状態遷移があるときに新しくデータを作成するため、ラベル付与に用いるデータとして、その管理期間の前の状態５３６、その状態の継続時間５３７、その管理期間の状態５３８の情報を持つ。前の状態の継続時間のクラス数を上記と同じ６とすると、前の状態と今の状態のクラス数はそれぞれ２なのでラベル数は２４と少ない数で管理することが可能となる。ただし、各管理データに紐付けられるセンサ信号データの数がそれぞれ異なるため、学習データ選択の際に、管理データの個数ではなく、センサ信号の数が所定数になるようにするべきであり、そのための処理が煩雑になるというデメリットがある。【００３９】
運転パターンラベルを付与する期間の長さおよびラベル付与のルールを決めるためのＧＵＩの例を、図６Ａから図６Ｄに示す。ＧＵＩは、入出力部１３０の表示画面１３１に表示される。はじめに、図示はしないが、ラベル付与のルールを決めるために参考にするセンサ信号とイベント信号を選んでロードしておく。次に、図６Ａに示す入力画面を用いて、ラベル付与に使用する情報の分類方法を決定する。図６Ａは、運転パターンラベル付与に使用する情報を選択するための入力画面の例である。運転パターンラベル付与参照情報入力ウィンドウ６０１には、周期情報入力ウィンドウ６０２、周期単位入力ウィンドウ６０３、参照情報選択ウィンドウ６０４Ａ〜６０４Ｃ、クラス数入力ウィンドウ６０７Ａ〜６０７Ｃなどが含まれている。周期情報入力ウィンドウ６０２には、運転パターンラベルを付与する期間を数値入力する。また、周期単位入力ウィンドウ６０３には、期間の単位を時間、日、週から選択入力する。参照情報選択ウィンドウ６０４Ａ〜６０４Ｃには、ラベル付与のために使用する情報を入力する。【００４０】
使用可能な情報は予め決まっているものとし、リスト表示ボタン６０５Ａ〜６０５Ｃのクリックにより参照可能情報リスト６０６が表示されるので、リストから選択入力する。クラス数入力ウィンドウ６０７Ａ〜６０７Ｃには、参照した情報をもとに異なるラベルを付与すべきクラス数を入力する。参照可能情報には、連続データと名義データの２種類があり、名義データの場合は、名義数がクラス数としてクラス数入力ウィンドウ６０７Ａ〜６０７Ｃに入力され、ユーザ入力を受け付けないようにする。また詳細ボタン６０８Ａ〜６０８Ｃは不活性にする。【００４１】
連続データの場合は、クラス数入力ウィンドウ６０７Ａ〜６０７Ｃにクラス数を数値入力する。２以上の整数値であり、上限は特に決まりはないが、例えば１０、１００など適当な数値に定めておく。連続データの場合は、詳細ボタン６０８Ａ〜６０８Ｃの何れかの押下により、図６Ｂに示す運転パターンラベル付与参照情報詳細入力ウィンドウ６１１が表示され、クラスの境界を定義するための詳細情報を入力する。追加ボタン６０９は、参照する情報を４以上にするためのボタンであり、追加ボタン６０９の押下により情報４、情報５、・・・に対応する参照情報選択ウィンドウ６０４、リスト表示ボタン６０５、クラス数入力ウィンドウ６０７、詳細ボタン６０８が順次表示される。詳細情報も含めて入力が終了したら、ラベル定義ボタン６１０を押下する。この操作により、図６Ｃに示すラベル確認画面が表示される。【００４２】
図６Ｂは、押下された詳細ボタン６０８Ａ〜６０８Ｃの何れかに対応する参照情報のクラス境界入力画面の例である。連続データを対象とするため、クラス境界の値を入力する形とする。運転パターンラベル付与参照情報詳細入力ウィンドウ６１１には、参照情報表示ウィンドウ６１２、数値範囲表示ウィンドウ６１３、クラス境界値入力ウィンドウ６１４Ａ〜６１４Ｅが含まれる。参照情報表示ウィンドウ６１２には、押下された詳細ボタン６０８Ａ〜６０８Ｃの何れかに対応する参照情報の名前が表示される。図には、運転パターンラベル付与参照情報入力ウィンドウ６０１において、詳細ボタン６０８Ｂが押下された場合の例を示している。数値範囲表示ウィンドウ６１３には、予めロードしておいたイベント信号１０３について、入力された周期の条件のもとで算出された、運転パターンラベル付与参照情報の値の範囲を表示する。【００４３】
クラス境界値入力ウィンドウ６１４Ａ〜６１４Ｅは、クラス数入力ウィンドウ６０７Ａ〜６０７Ｃで入力された数値より１個少ない数だけ表示される。クラス数６個の場合、ｘ未満ならばクラスＮである（Ｎは１から５）という状態を決定する条件のｘを数値入力する。図の例は１未満ならばクラス１、すなわち状態遷移回数が０回の場合はクラス１である。また１以上２未満ならばクラス２、すなわち状態遷移回数が１回の場合はクラス２である。以下同様にクラス５までが定義され、クラス６の条件はその他すなわち状態遷移回数５回以上の場合である。入力終了後、終了ボタン６１５の押下により、運転パターンラベル付与参照情報入力ウィンドウ６０１に戻る。【００４４】
図６Ｃは、図６Ａの入力結果に従って定義されたラベルの確認画面の例である。運転パターンラベル確認ウィンドウ６２１は、ラベル表示ウィンドウ６２２Ａ、６２２Ｂとクラス境界条件表示ウィンドウ６２６を備えている。ラベル表示ウィンドウ６２２Ａ、６２２Ｂは２次元のマトリクスであり、条件６２３Ａ、６２３Ｂのときに、横軸項目６２４Ａ、６２４Ｂ、縦軸項目６２５Ａ、６２５Ｂに示した情報のクラスに対応するラベル番号を表示したものである。横軸項目６２４Ａ、６２４Ｂと縦軸項目６２５Ａ、６２５Ｂには、参照情報のうちクラス数の多いものから２個を選択して設定するとウィンドウ６２２の数を最小にできる。条件６２３Ａ、６２３Ｂには残りの参照情報のクラスの組合せを別々に全て記述する。【００４５】
表示されるラベルは、各情報のクラスの組合せをもとに以下に示す方法で自動的に決められる。情報１のクラスＣ１、情報２のクラスＣ２、情報３のクラスＣ３を（Ｃ１，Ｃ２，Ｃ３）と表記すると、（１，１，１）から開始して（１，１，２）、（１，１，３）、・・・（１，１，６）、（１，２，１）、・・・（２，６，６）までの全ての組合せに１から始まる通し番号をラベルとして付与する。この場合、情報１のクラス数が２、情報２のクラス数が６、情報３のクラス数が６なのでラベル総数は７２となる。情報４を参照する場合は（Ｃ１，Ｃ２，Ｃ３，Ｃ４）の全ての組合せを順に数え上げて通し番号をつける。さらに他の情報を参照する場合も同様である。クラス境界条件表示ウィンドウ６２６には参照情報の項目名とそれに対応するクラスを定義する条件が表示される。【００４６】
図６Ａに示す運転パターンラベル付与参照情報入力ウィンドウ６０１で設定された参照情報に対し、参照情報が名義データの場合はそれぞれの名義を個別にクラスに割り当てて表示し、参照情報が連続データの場合は、図６Ｂに示す運転パターンラベル付与参照情報詳細入力ウィンドウ６１１で入力されたクラス境界条件をそのまま表示する。グラフ表示ボタン６２７押下により、図６Ｄに示すグラフ表示ウィンドウ６３１が表示される。登録ボタン６２８押下によりラベル付与条件を確定し、登録する。戻るボタン６２９押下により、図６Ａに示す運転パターンラベル付与参照情報入力ウィンドウ６０１に戻る。【００４７】
図６Ｄは、入力された条件に基づき付与されたラベルをセンサデータと関連付けて確認するためのグラフ表示画面の例である。グラフ表示ウィンドウ６３１は、ラベル付与参照情報表示ウィンドウ６３２、ＯＮ／ＯＦＦ情報表示ウィンドウ６３３、センサ信号表示ウィンドウ６３４を含む。【００４８】
グラフ表示に先立ち、予めロードされたセンサ信号とイベント信号を用い、図３Ａに示す処理フローにより、運転パターンラベルを付与する。ラベル付与参照情報表示ウィンドウ６３２には、ラベルを付与するために参照された情報すなわち初期状態、状態遷移回数、初期状態継続時間のクラス番号が折れ線グラフで表示され、重ねてラベル番号が表示される。【００４９】
ＯＮ／ＯＦＦ情報表示ウィンドウ６３３には、モード分割（Ｓ３０３）によって得られるＯＮ／ＯＦＦの情報が時系列グラフとして表示される。センサ信号表示ウィンドウ６３４には、センサ選択ウィンドウ６３７により選択されたセンサ信号が時系列グラフとして表示される。スクロールバー６３５の移動によりウィンドウに表示されるデータの期間が変更される。期間表示ウィンドウ６３６には、表示中のデータの期間が表示される。【００５０】
逆にユーザによる入力に従って、スクロールバー６３５が移動し、ウィンドウ６３１，６３２，６３３の表示が変更される。戻るボタン６３８の押下により、グラフ表示ウィンドウ６３１を消去し、図６Ｃに示す運転パターンラベル確認ウィンドウ６２１に戻る。【００５１】
以上説明したＧＵＩにより、ラベル付与の条件を簡便に入力することが可能である。決定した条件を用いて付与したラベルとセンサデータを対応付けて観察することが可能であるため、同一ラベルが付与された期間の設備の状態が互いに類似しているか否かを確認でき、これによりラベル付与条件の妥当性を確認することが可能である。【００５２】
次に、学習データ選択部１０６における処理の流れを、図７を用いて説明する。まず、運転パターンラベル付与部１０４にセンサ信号１０２とイベント信号１０３とを入力し（Ｓ７０１）、予め定められた期間毎に、運転パターンラベル付与部１０４で運転パターンラベルを付与する（Ｓ７０２）。次に、学習データ選択部１０６において、以下のような手順で学習データ選択処理を行う。先ず、付与されたラベルと同一のラベルを探索対象とし（Ｓ７０３）、データベース１０５から、同一サイト同一号機の運転パターンラベル管理データをロードする（Ｓ７０４）。次に、探索対象ラベルを探索し（Ｓ７０５）、探索結果を記録する（Ｓ７０６）。ただし、警告有りのデータは探索の対象としない。また、ファイルまたはユーザインターフェースにより除外日を指定しておき、これも探索の対象からはずすことが可能なようにしておくとよい。更に、探索されたデータ数が予め定めた数に達しているかどうか確認し（Ｓ７０７）、達していなければ探索対象ラベルを設定し（Ｓ７０８）、ステップＳ７０５に戻る。【００５３】
各運転パターンラベルについて、類似した設備状態が含まれることが期待される他のラベルを優先順位つきで予め定めておき、ステップＳ７０８においては、その優先順位に従って探索対象ラベルを決める。学習データが所定数に達するまでこれを繰り返し、探索された運転パターンラベル管理データに記載されたキーコードに紐付けられるセンサ信号データとイベント信号データを時系列順にロードする（Ｓ７０９）。ただし、最後の探索で所定数を超えた分については、最後の探索対象ラベルの中で、新しいものから順に選択して合計が所定数になるようにする。【００５４】
以上の処理により、蓄積した全てのデータの中から入力データに類似した学習データを選択することができる。その結果、高精度に異常識別が可能となり、学習データ不足に起因する誤報を低減することができる。このとき、学習データ数を所定の数に制限するため、計算時間が増えることがない。【００５５】
次に、正常モデル作成部１０７における処理の流れを、図８を用いて説明する。はじめに、学習データ選択部１０６において選択されたセンサ信号データとイベント信号データを学習データとして入力する（Ｓ８０１）。次に、特徴抽出を行って特徴ベクトルを作成し（Ｓ８０２）、センサ信号データとイベント信号データのキーコードに従って別々のグループに分割し（Ｓ８０３）、そのうち１グループを除いたデータを用いて学習を行い、正常モデルを作成する（Ｓ８０４）。作成された正常モデルを用い、ステップＳ８０４で除かれた１グループのデータを入力して異常測度を算出する（Ｓ８０５）。すべてのグループのデータについて異常測度の算出が終了したかをチェックし（Ｓ８０６）、終了していなければまだ異常測度を算出していないグループについて（Ｓ８０７）、正常モデル作成（Ｓ８０４）と異常測度算出（Ｓ８０５）のステップを繰り返す。すべてのグループのデータについて異常測度の算出が終了したら（Ｓ８０６）、算出された異常測度に基づいて異常を識別するしきい値を設定する（Ｓ８０８）。最後に、すべての学習データを用いて正常モデルを作成する（Ｓ８０９）。【００５６】
次に、各ステップについて詳細に説明する。  
先ず、ステップＳ８０１において、センサ信号を入力し、次に、ステップＳ８０２において特徴抽出を行い、特徴ベクトルを得る。特徴抽出としては、センサ信号をそのまま用いることが考えられるが、ある時刻に対して±1,±2,…のウィンドウを設け，ウィンドウ幅（3,5,…）×センサ数の特徴ベクトルにより、データの時間変化を表す特徴を抽出することもできる。また、離散ウェーブレット変換(DWT: Discrete Wavelet Transform)を施して、周波数成分に分解してもよい。なお、各特徴は平均と標準偏差を用いて、平均を０、分散を１となるように変換する正準化を行うとよい。評価時に同じ変換ができるよう、各特徴の平均と標準偏差を記憶しておく。あるいは、最大値と最小値または予め設定した上限値と下限値を用いて正規化を行ってもよい。【００５７】
これらの処理は、単位およびスケールの異なるセンサ信号を同時に扱うためのものである。このとき、最低限の処理として、分散が非常に小さいセンサ信号および単調増加するセンサ信号を除く必要がある。【００５８】
また、相関解析による無効信号を削除することも考えられる。これは、多次元時系列信号に対して相関解析を行い、相関値が１に近い複数の信号があるなど、極めて類似性が高い場合に、これらは冗長だとして、この複数の信号から重複する信号を削除し、重複しないものを残す方法である。このほか、ユーザが指定するようにしてもよい。さらに、ステップＳ８０２において、主成分分析、独立成分分析、非負行列因子分解、潜在構造射影、正準相関分析など様々な特徴変換手法により、次元削減を行ってもよい。【００５９】
次に、ステップＳ８０３において学習データをグループに分割し、ステップＳ８０４において正常モデル作成を行い、ステップＳ８０５において異常測度算出部１０８で異常測度の算出を行う。【００６０】
正常モデル作成手法としては、投影距離法(ＰＤＭ: Projection Distance Method)や局所部分空間法(ＬＳＣ: Local Sub-space Classifier)が考えられる。投影距離法は、学習データに対し独自の原点をもつ部分空間すなわちアフィン部分空間（分散最大の空間）を作成する方法である。【００６１】
クラスタ毎に、図９に示すようにアフィン部分空間を作成する。図９では、３次元の特徴空間において、１次元のアフィン部分空間を作成する例を示しているが、特徴空間の次元はもっと大きくてもよく、アフィン部分空間の次元も特徴空間の次元より小さくかつ学習データ数より小さければ何次元でもかまわない。【００６２】
アフィン部分空間の算出方法について説明する。まず、学習データの平均μと共分散行列Σ を求め、次にΣの固有値問題を解いて値の大きい方から予め指定したr個の固有値に対応する固有ベクトルを並べた行列Ｕをアフィン部分空間の正規直交基底とする。この正常モデルに基づいて算出する異常測度は、各クラスタのアフィン部分空間への投影距離のdの最小値と定義する。ここで、クラスタとは、例えば図３Ｂに示すようにモードに分割された各区間のデータを集めたものとする。あるいは、ｋ平均法に代表される教師なしクラスタリング手法を利用してもよい。【００６３】
一方、局所部分空間法は、評価データqのk-近傍データを用いてｋ−１次元のアフィン部分空間を作成する方法である。【００６４】
図１０にｋ＝３の場合の例を示す。図１０に示すように、異常測度は図に示す投影距離で表されるため、評価データqに最も近いアフィン部分空間上の点bを求めればよい。【００６５】
評価データqとそのk-近傍データxi( i = 1,…,k )からbを算出するには、qをk個並べた行列Qとxiを並べた行列Xから  【００６６】
◎
【数１】
2013030984.tif 000003
【００６７】
により相関行列Cを求め、  【００６８】
◎
【数２】
2013030984.tif 000004
【００６９】
によりbを計算する。【００７０】
この方法は、評価データを入力しないとアフィン部分空間を作成できないため、ステップＳ８０４およびステップＳ８０９においては、k-近傍データを効率的に探すためのｋｄ木をモード別に構築しておく。ｋｄ木とは、k次元のユークリッド空間にある点を分類する空間分割データ構造である。座標軸の1つに垂直な平面だけを使って分割を行い、各葉ノードには1つの点が格納されるよう構成する。ステップＳ８０５では、評価データと同じモードに属するｋｄ木を利用して評価データのk-近傍データを求め、それらから前述の点ｂを求め、評価データと点ｂの距離を算出して異常測度とする。【００７１】
このほか、マハラノビスタグチ法、回帰分析法、最近傍法、類似度ベースモデル、１ク
ラスＳＶＭなど様々な方法を用いて正常モデルの作成が可能である。【００７２】
次に、ステップＳ８０８において、異常測度をもとにしきい値を設定する。具体的には、Ｓ８０５において算出した学習データに対応する異常測度を昇順にソートし、予め指定した比率に到達する値をしきい値とする。すなわち、データ数をＮ個、指定した比率をｐとするとＮｐ番目に小さい値をしきい値とする。【００７３】
学習部は正常データから構成されると考えればｐを1.0に設定すべきであり、この場合は異常測度の最大値をしきい値とする。このとき、全ての学習データに対する異常測度のを用いずに、異常識別の対象とする入力データに付与された運転パターンラベルと同じラベルを付与された学習データに対する異常測度のデータのみを用いるとなおよい。また、モード別にデータを分けて同様の処理によりしきい値設定を行うとよい。モード分割処理は正常モデル作成部１０７で前述の方法で行ってもよいし、運転パターンラベル付与部１０４で得られる結果をもとに、図５Ｂに示すセンサ信号データにモードを付加して保存しておいてもよい。【００７４】
正常モデル作成部１０７で正常モデル作成後、異常測度算出部１０８において、入力センサ信号１０２に基づき異常測度を算出する。ステップＳ８０２と同様の方法で特徴ベクトルを抽出したあと、ステップＳ８０９で作成された正常モデルに基づき、ステップＳ８０５と同様の方法で異常測度を算出する。次に、異常識別部１０９において、異常測度算出部１０８で算出された異常測度をステップＳ８０８で設定されたしきい値と比較し、しきい値以上の値の場合は異常として検知する。【００７５】
上記に説明した方法により、一定期間毎にデータに運転パターンラベルを付与しておき、異常検知時には同じラベルあるいは状態の近いラベルの期間のデータを所定数集めて学習データとするため、計算負荷を増やさずに高精度な正常モデルを作成することができ、正常を異常と判定する誤報を低減できる。【００７６】
次に、異常診断部１２０においては、設備１０１からのセンサ信号１０２とイベント信号１０３とを入力すると共に、運転パターンラベルを付与されたセンサ信号１０２とイベント信号１０３とをセンサ信号解析部１００で解析した結果として異常識別部１０９から出力される信号を受けて、異常識別部１０９において異常と判定された事象から異常の原因事象を推定し、この原因事象によってもたらされる結果事象を異常診断部１２０の表示画面（図示せず）上に表示する。この結果事象が複数あるときには、発生確率の高い順に並べて画面上に表示する。運転パターンラベルを付与されたセンサ信号１０２とイベント信号１０３とを処理することによりセンサ信号解析部１００では比較的早い処理が実現でき、異常診断部１２０での診断処理は基本的にリアルタイム処理で行うことができる。
【実施例２】
【００７７】
実施例１においては、設備の状態は、ON/OFFの切り換えパターンすなわち運転パターンに応じて変化するという仮定のもと、運転パターンラベルに基づいて学習データを選択する方法について説明してきた。しかし、同じ運転パターンラベルであっても保守作業の前後などには状態が大きく変化する場合があり、上記説明した方法では、適切な学習データを選択することができない可能性がある。【００７８】
このような場合に対応するための機能を追加した第二の実施例について図１１を用いて説明する。図１１は、本発明の設備状態監視方法を実現する第二の実施例のシステム構成を表す図である。本システムは、センサ信号解析部１１００と異常診断部１１２０、入出力部１１３０とを備えて構成されている。【００７９】
本システムのセンサ信号解析部１１００は、設備１０１から出力されるセンサ信号１０２とイベント信号１０３に基づき一定期間毎に運転パターンラベルを付与する運転パターンラベル付与部１０４と、運転パターンラベル付与部１０４でラベルが付与されたセンサ信号１０２とイベント信号１０３とからマクロ特徴を算出するマクロ特徴算出部１１０１、運転パターンラベルを付与されたセンサ信号１０２とイベント信号１０３を蓄積するデータベース１１０５、運転パターンラベルに基づいて蓄積されたセンサ信号１０２とイベント信号１０３から学習データを選択する学習データ選択部１１０６、選択した学習データを用いて正常モデルを作成する正常モデル作成部１１０７、正常モデルとセンサ信号に基づき異常測度を算出する異常測度算出部１１０８、算出された異常測度に基づき異常を検知する異常識別部１１０９とを備えて構成される。【００８０】
異常診断部１１２０は、センサ信号１０２とイベント信号１０３及び異常識別部１１０９からの出力を受けて、設備１０１の異常を診断する。【００８１】
入出力部１１３０は、運転パターンラベル付与部１０４、異常識別部１１０９、異常診断部１１２０などと接続して診断条件の入力や診断結果の出力を行う。入出力部１１３０は表示画面１１３１を備えている。【００８２】
図１に示した第一の実施例の構成との違いは、センサ信号１０２のマクロな変動を表す特徴を算出するマクロ特徴算出部１１０１を追加した点である。運転パターンラベル付与部１０４において、予め定めた期間毎、例えば１日毎に運転パターンラベルを付与した後、マクロ特徴算出部１１０１において、同じ期間毎にセンサ信号を切り出して、特徴量を算出する。特徴量としては、期間全体に亘る平均、分散、最大値、最小値、期間中の定常ONの時刻の平均、分散、期間中の定常OFFの時刻の平均、分散などが考えられる。【００８３】
図１２に、その管理データの一例を示す。この管理データは、図５Ａに示す運転パターンラベル管理データ１２０１〜１２０９（それぞれ図５Ａの５０１〜５０９に対応）に、上記した特徴量、即ち、期間全体に亘る平均、分散、最大値、最小値：１２１０、期間中の定常ONの時刻の平均、分散：１２１１、期間中の定常OFFの時刻の平均、分散：１２１２の項目を付加したものである。但し付加するのは必ずしもこれら３つの項目全てではなくても良く、何れか一つ又は二つの組合せであってもよい。【００８４】
次に、学習データ選択部１０６における、学習データ選択処理方法を説明する。基本的には第一の実施例と同様の方法で学習データを選択するが、保守作業などの状態が大きく変化する可能性がある場合には、運転パターンラベルに加え、マクロ特徴も用いて学習データを選択する。【００８５】
具体的には、図１３のフロー図に示すように、まず、運転パターンラベル付与部１０４にセンサ信号１０２とイベント信号１０３とを運転パターンラベル付与部１０４に入力し（Ｓ１３０１）、予め定められた期間毎に運転パターンラベルを付与し、マクロ特徴算出部１１０１でマクロ特徴を算出して運転パターンラベルに対応付けて保存する（Ｓ１３０２）。次に、学習データ選択部１１０６において、以下のような手順で学習データ選択処理を行う。【００８６】
先ず、実施例１の場合と同様に、図６Ａ乃至図６Ｄで説明した画面上で運転パターンラベルを探索対象として設定し（Ｓ１３０３）、データベース１１０５から、同一サイト同一号機の運転パターンラベル管理データをロードする（Ｓ１３０４）。【００８７】
次に、保存された運転パターンラベル管理データから、評価対象の期間に付与された運転パターンラベルと同一の運転パターンラベルかつマクロ特徴が類似しているデータを探索し（Ｓ１３０５）、探索結果を記録する（Ｓ１３０６）。マクロ特徴が類似しているか否かは予め設定したしきい値に基づいて判断するものとする。ただし、警告有りのデータは探索の対象としない。また、ファイルまたはユーザインターフェースにより除外日を指定しておき、これも探索の対象からはずすことが可能なようにしておくとよい。【００８８】
更に、探索されたデータ数が予め定めた数に達しているかどうか確認し（Ｓ１３０７）、探索されたデータ数が予め定めた数に達していなければ、予め設定された優先順位に従って探索対象ラベルを変更して新たな対象ラベルを設定し（Ｓ１３０８）、Ｓ１３０５に戻って、そのラベルかつマクロ特徴が類似しているデータを探索する。学習データが所定数に達するまでこれを繰り返し、探索された運転パターンラベル管理データに記載されたキーコードに紐付けられるセンサ信号データとイベント信号データを時系列順にロードする（Ｓ１３０９）。最後の探索で所定数を超えた分については、最後の探索対象ラベルの中で、特徴が類似したものから順に選択して合計が所定数になるようにする。【００８９】
マクロ特徴を用いるかどうかの処理の切り換え方法については、何通りかの方法が考えられる。第一の方法は、切り換えを行わず、常にマクロ特徴を用いて学習パターンを選択することである。この方法は、判定条件が不要のため簡便に実現可能である。【００９０】
第二の方法は、イベント信号に基づき保守作業の有無を判定し、保守作業があったならばマクロ特徴を用いて学習パターンを選択する方法である。イベント信号に保守作業の開始、終了を表す信号が含まれているならば、この方法を採用するとよい。保守作業がない場合はデータの新しい方から探索していけばよいため、所定数に達した時点で探索を打ち切ることが可能である。メインの運転パターンがあって、同じラベルが大量にある場合には、第一の方法と比較して探索の効率が良くなることが期待される。【００９１】
第三の方法は、マクロ特徴の推移を調べて保守作業の有無を判定し、保守作業があったならばマクロ特徴を用いて学習パターンを選択する方法である。例えば評価対象期間の前の数周期について各特徴の平均μと標準偏差σを算出し、評価期間の特徴値がμ±ｎσの範囲を超える特徴が一つでもあれば、保守作業など大きな状態変化があったものと判断する。【００９２】
ここで、調べる周期数とｎはパラメータである。あるいは、評価対象期間の前の数周期について１つ前の周期との特徴空間上での距離を算出し、評価対象期間の１つ前の周期との特徴空間上での距離がそれらの最大値より大きい場合、保守作業など大きな状態変化があったものと判断する。【００９３】
そのほか、評価対象期間の前の数周期についてマクロ特徴空間の重心と共分散を算出し、マハラノビス距離で判断してもよく、評価対象期間の前の数周期分を正事例の学習データとして既存の１クラス識別手法を利用してもよい。【００９４】
以上の処理により、保守作業など設備の大きな状態があった場合にも、設備状態が近い学習データを選択することが可能となり、精度の高い正常モデルを作成することが可能となり、その結果学習データ不足に起因する誤報を低減することができるようになる。【００９５】
ただし、上記方法では、運転パターンラベルが同じまたは状態が近いラベルかつセンサ信号から抽出したマクロ特徴が類似したデータ数が所定数に満たない場合が生じることが想定される。この場合に対応するため、マクロ特徴が類似したデータ数が少ない場合には、保守作業によって変化したと考えられるセンサを除外して、学習データ選択以降の処理を実行する。【００９６】
保守作業によって変化したか否かは、前述のように各特徴の平均μと標準偏差σを算出し評価期間の特徴値がμ±ｎσの範囲を超えるか否かによって判断する。あるいは、各特徴の一つ前の周期との差を算出し、評価期間の特徴の一つ前の周期との差がそれまでの最大となるか否かで判断する。
【実施例３】
【００９７】
上記に説明した実施例１及び２では、設備から出力されたセンサ信号およびイベント信号を、全て運転パターンラベルを付与して保存しておくことを想定している。しかし、この方法では、年月を経るに従いデータベースの容量が増大するという問題が生じる。この問題を解決するため、冗長なデータを削除することを考える。削除対象とするのは、同じラベルで、学習データとして選択する数より大きい予め設定された数を超えるものとし、運転パターンラベル管理データの単位で、削除してよいかどうか判断する。【００９８】
基本的には除いても影響の小さいデータを削除する。すなわち類似したデータが他にあることや実際に除いても結果が変わらないことを確認した上で削除する。【００９９】
例えば、マクロ特徴の特徴空間上の最近傍データまでの距離が最小のもののうち二番目に近いデータまでの距離が小さい方を削除する。あるいは、同じラベルのデータを用い、マクロ特徴を用いてクラスタリングを行い、所属データが多いクラスタの中から日付の古いものを削除する。あるいは、同じラベルのデータを全て用いて評価対象データを除いて学習データとする交差検証により異常測度を算出し、異常測度の最大値が最小となる期間のデータを削除する。【０１００】
また、同じラベルのデータから１期間ずつ除いて上記交差検証による異常測度算出を行い、異常測度算出結果の変化が最も小さい除外期間のデータを削除してもよい。【０１０１】
削除対象が確定したら、その管理データに紐付けられたセンサ信号データおよびイベント信号データをデータベース上から完全に削除する。管理データは削除せずに、キーコード欄にデータを削除したことを示す数値、例えば０を入力しておく。【０１０２】
上記方法により、計算負荷を増加させずに高精度な正常モデルを作成可能な上、データベース容量を削減する効果も得られる。
【産業上の利用可能性】
【０１０３】
本発明は、ガスタービンなどの電源設備、ガスタービンの廃熱などを利用して地域暖房用温水を供給する設備、工場向けに高圧蒸気や低圧蒸気を供給する設備、水力発電所での水車、原子力発電所の原子炉、風力発電所の風車、航空機や重機のエンジン、鉄道車両や軌道、エスカレータ、エレベータなどのプラントや設備などが出力する多次元時系列データをもとに異常を早期に検知する予防保全に利用可能である。
【符号の説明】
【０１０４】
１００、１１００…センサ信号解析部    １０１…設備    １０２…センサ信号    １０３…イベント信号  １０４…運転パターンラベル付与部    １０５、１１０５…データベース    １０６、１１０６…学習データ選択部    １０７、１１０７…正常モデル作成部    １０８、１１０８…異常測度算出部    １０９、１１０９…異常識別部  １２０、１１２０…異常診断部
【図１】
2013030984.tif 000005
【図２Ａ】
2013030984.tif 000006
【図２Ｂ】
2013030984.tif 000007
【図３Ａ】
2013030984.tif 000008
【図３Ｂ】
2013030984.tif 000009
【図４Ａ】
2013030984.tif 000010
【図４Ｂ】
2013030984.tif 000011
【図５Ａ】
2013030984.tif 000012
【図５Ｂ】
2013030984.tif 000013
【図５Ｃ】
2013030984.tif 000014
【図５Ｄ】
2013030984.tif 000015
【図６Ａ】
2013030984.tif 000016
【図６Ｂ】
2013030984.tif 000017
【図６Ｃ】
2013030984.tif 000018
【図６Ｄ】
2013030984.tif 000019
【図７】
2013030984.tif 000020
【図８】
2013030984.tif 000021
【図９】
2013030984.tif 000022
【図１０】
2013030984.tif 000023
【図１１】
2013030984.tif 000024
【図１２】
2013030984.tif 000025
【図１３】
2013030984.tif 000026
【国際調査報告】
2013030984.tif 000027
2013030984.tif 000028