(57)【要約】      （修正有）【課題】変化の激しい起動や停止シーケンスに対しても、その激しい変化およびその変化の統計確率的特性を監視しながら、異常を検知可能な異常検知方法を提供する。【解決手段】設備１０１における取得できないセンサデータを推定することにより、センサのサンプリング時刻が、シーケンスの開始と同期していない場合に発生する異なる時間で取得した同じシーケンスのセンサデータ間にある時間ずれ、および多次元センサを同期してデータを取得できない場合に発生する異なるセンサのセンサデータ間にある時間ずれを解消する。【選択図】図１Ａ
2014048697.tif 000002
【特許請求の範囲】
【請求項１】
プラント又は設備の異常を検知する方法であって、
プラント又は設備に取り付けたセンサから間欠的に出力されたセンサ信号及び該センサ信号と同じ期間における前記プラント又は設備の起動シーケンス又は停止シーケンスの開始および終了に対応するイベント信号を入力し、
該入力したセンサ信号のうち前記起動シーケンス又は終了シーケンスの開始のイベント信号と終了のイベント信号との間の区間に対応するセンサ信号を切り出し、
該切り出したセンサ信号のある時刻における信号値とその確率分布を推定し、
該推定した確率分布に基づいて特徴量を抽出し、
該抽出した特徴量に基づいて前記プラント又は設備の異常を検知する
ことを特徴とする設備状態監視方法。
【請求項２】
前記切り出したセンサ信号のある時刻における信号値とその確率分布を推定することを、
前記切り出したセンサ信号を前記起動シーケンス又は終了シーケンスの開始のイベント信号を起点とした時間に合わせ、
該時間を合わせたセンサ信号のデータを推定する時刻を決定し、
該決定した時刻におけるセンサデータを推定し、
該推定したセンサデータの確率分布推定する
ことにより行うことを特徴とする請求項１記載の設備状態監視方法。
【請求項３】
前記センサデータを推定する手法を画面上に表示された複数の手法の中から選択し、該選択した手法に基づいて前記センサデータを推定することを特徴とする請求項２に記載の設備状態監視方法。
【請求項４】
前記推定したセンサデータに関する情報を画面上に表示することを特徴とする請求項２記載の設備状態監視方法。
【請求項５】
前記抽出した特徴量に基づいて前記プラント又は設備の異常を検知することを、前記プラント又は設備が正常に作動しているときの前記センサ信号を用いて前記センサ信号の正常空間或いは識別境界を決定し、前記抽出した特徴量が前記決定した正常空間或いは識別境界の内部にあるか否かを判断して前記プラント又は設備の異常を検知することを特徴とする請求項１記載の設備状態監視方法。
【請求項６】
プラント又は設備の異常を検知する装置であって、
プラント又は設備に取り付けたセンサから間欠的に出力されたセンサ信号及び該センサ信号と同じ期間における前記プラント又は設備の起動シーケンス又は停止シーケンスの開始および終了に対応するイベント信号を入力し、該入力したセンサ信号のうち前記起動シーケンス又は終了シーケンスの開始のイベント信号と終了のイベント信号との間の区間に対応するセンサ信号を切り出し、該切り出したセンサ信号を前記起動シーケンス又は終了シーケンスの開始のイベント信号を起点とした時間に合わせる処理を行うデータ前処理部と、
該データ前処理部で処理したセンサ信号のある時刻における信号値とその確率分布を推定する確率分布推定部と、
該確率分布推定部で推定した確率分布に基づいて特徴量を抽出する特徴量抽出部と、
該特徴量抽出部で抽出した特徴量に基づいて前記プラント又は設備の異常を検出する異常検出部と
入出力する情報を表示する画面を有して前記異常検出部で検出した前記プラント又は設備の異常に関する情報を前記画面上に表示する入出力部と
を備えたことを特徴とする設備状態監視装置。
【請求項７】
前記確率分布推定部は、
前記切り出したセンサ信号を前記起動シーケンス又は終了シーケンスの開始のイベント信号を起点とした時間に合わせたセンサ信号のデータを推定する時刻を決定する推定時刻決定部と、
該推定時刻決定部で決定した時刻におけるセンサデータを推定するセンサデータ推定部と、
該センサデータ推定部で推定したセンサデータの統計確率分布推定する統計確率分布推定部と、
を備えることを特徴とする請求項６記載の設備状態監視装置。
【請求項８】
前記入出力部は前記センサデータ推定部でセンサデータを推定するための複数の手法を前記画面上に表示し、前記センサデータ推定部は、前記表示された複数の手法の中から前記画面上で選択された手法に基づいて前記センサデータを推定することを特徴とする請求項７に記載の設備状態監視装置。
【請求項９】
前記入出力部は、前記センサデータ推定部で推定したセンサデータに関する情報を前記画面上に表示することを特徴とする請求項７記載の設備状態監視装置。
【請求項１０】
前記異常検出部は、前記プラント又は設備が正常に作動しているときの前記センサ信号を用いて前記センサ信号の正常空間或いは識別境界を決定する学習部と、前記特徴量抽出部で抽出した特徴量が前記決定した正常空間或いは識別境界の内部にあるか否かを判断して前記プラント又は設備の異常を検知する異常検知部とを備えていることを特徴とする請求項６記載の設備状態監視装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プラントや設備などから出力される多次元時系列データをもとに、変化の激しい起動や停止シーケンス時の設備の不具合あるいはその兆候を早期に検知する、また、コストダウンのため、長くしたサンプリング間隔で取得できない激しい変化を復元してその変化の統計確率的特性を監視する設備状態監視方法及び設備状態監視装置に関する。
【背景技術】
【０００２】
電力会社では、ガスタービンの廃熱などを利用して地域暖房用温水を供給したり、工場向けに高圧蒸気や低圧蒸気を供給したりしている。石油化学会社では、ガスタービンなどを電源設備として運転している。このようにガスタービンなどを用いた各種プラントや設備において、設備の不具合あるいはその兆候を検知する予防保全は、社会へのダメージを最小限に抑えるためにも極めて重要である。特に、起動や停止のような変化の激しいシーケンスにおいては、故障が多発し易いため、その期間中の異常を早期に検知することが重要である。
【０００３】
ガスタービンや蒸気タービンのみならず、水力発電所での水車、原子力発電所の原子炉、風力発電所の風車、航空機や重機のエンジン、鉄道車両や軌道、エスカレータ、エレベータ、また、切削、穴開けなどの加工装置は、不良発生を防ぐため、そのパフォーマンスに異常があれば、ただちに検知することが要求されている。
【０００４】
このため、対象設備やプラントに複数のセンサを取り付け、センサ毎の監視基準に従って正常か異常かの自動判別が行われている。
【０００５】
設備あるいは製造装置や計測装置を対象とした通常運転時の異常検知に有効である一例は特許文献１（特開２０１１−０７０６３５号）に開示されている。特許文献１の開示例においては、設備の多次元センサデータを特徴空間へ写像し、特徴空間で正常モデルを作成してから、新たに入力されたセンサデータが正常モデルへの投影距離を異常測度とし、その異常測度の所定閾値を越えたかどうかによって異常検知を行っている。
【０００６】
また、時系列センサデータの統計確率的特性を同時に監視可能な統計確率的特性を表すパラメータを算出しながら、異常検知する代表的な手法として、非特許文献１および非特許文献２に開示されている方法がある。この方法では、直接にセンサ波形から時刻毎に算出した統計確率的なパラメータを用い、正常モデルを作成し、そのモデルからの離れ度合いを用いて異常を検知する。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１１−０７０６３５号公報
【非特許文献】
【０００８】
【非特許文献１】センサモデルと投票に基づく発電機の異常検知の検討、情報科学技術フォーラム講演論文集８（３）、１３９−１４２、２００９
【非特許文献２】尤度ヒストグラムに基づく異常検出、パターン認識・メディア理解研究会（ＰＲＭＵ）、２０１１
【発明の概要】
【発明が解決しようとする課題】
【０００９】
特許文献１に開示されている技術は、変化の激しい起動や停止シーケンス時、あるいは、負荷変動の大きい加工装置に発生する異常を予兆検知するのが困難である。図１２に運転モードの違いによる特徴空間上における特許文献１に記載の方法による局所空間の変化を示す。図からわかるように、（ａ）の定常運転時の場合、取得した正常センサデータのデータ値は同じレベルであるため、特徴空間上においては、正常なセンサデータから作られた正常な局所空間が小さく、異常が発生する場合に、異常データが正常局所空間から大きく離れ、識別されやすい。
【００１０】
一方、変化の激しいシーケンス、例えば、（ｂ）起動シーケンスの場合に、取得した正常センサデータのデータ値の変化は大きいため、特徴空間上において、正常なセンサデータから作られた正常な局所空間が定常運転時より大きく広がっており、シーケンス期間中、異常が発生する場合に、特徴空間上において異常データが正常な局所空間にあり、異常として検知されることは困難である。
【００１１】
また、非特許文献１および非特許文献２に開示された方法によれば、時刻毎に異常度を算出するため、変化が激しい場合のシーケンスの異常を検知できる。しかし、プラントなどの設備において、コストダウンのため、長いサンプリング間隔のセンサデータしか取得していない場合に、その激しい変化を十分にキャッチできない。また、センサのサンプリング時刻が、シーケンスの開始と同期していない場合、異なる時間で取得した同じシーケンスのセンサデータに時間ずれが発生する。また、多次元センサを同期してデータを取得できない場合に、センサ間のセンサデータにも時間ずれが発生する。このため、上記に示した非特許文献１および非特許文献２に開示されている技術は、時刻毎での統計確率的なパラメータを算出できず、異常検知を行うことができない。
【００１２】
そこで、本発明では、上記した従来技術の課題を解決し、変化の激しい起動や停止シーケンスに対しても、その激しい変化およびその変化の統計確率的特性を監視しながら、異常を検知可能な異常検知方法を備えた設備状態監視方法および設備状態監視装置を提供するものである。
【課題を解決するための手段】
【００１３】
上記課題を解決するために、本発明は、プラント又は設備の異常を検知する方法を、プラント又は設備に取り付けたセンサから間欠的に出力されたセンサ信号及び該センサ信号と同じ期間における前記プラント又は設備の起動シーケンス又は停止シーケンスの開始および終了に対応するイベント信号を入力し、該入力したセンサ信号のうち前記起動シーケンス又は終了シーケンスの開始のイベント信号と終了のイベント信号との間の区間に対応するセンサ信号を切り出し、該切り出したセンサ信号のある時刻における信号値とその確率分布を推定し、該推定した確率分布に基づいて特徴量を抽出し、該抽出した特徴量に基づいて前記プラント又は設備の異常を検知するようにした。
【００１４】
又、上記課題を解決するために、本発明では、プラント又は設備の異常を検知する装置を、プラント又は設備に取り付けたセンサから間欠的に出力されたセンサ信号及び該センサ信号と同じ期間における前記プラント又は設備の起動シーケンス又は停止シーケンスの開始および終了に対応するイベント信号を入力し、該入力したセンサ信号のうち前記起動シーケンス又は終了シーケンスの開始のイベント信号と終了のイベント信号との間の区間に対応するセンサ信号を切り出し、該切り出したセンサ信号を前記起動シーケンス又は終了シーケンスの開始のイベント信号を起点とした時間に合わせる処理を行うデータ前処理部と、該データ前処理部で処理したセンサ信号のある時刻における信号値とその確率分布を推定する確率分布推定部と、該確率分布推定部で推定した確率分布に基づいて特徴量を抽出する特徴量抽出部と、該特徴量抽出部で抽出した特徴量に基づいて前記プラント又は設備の異常を検出する異常検出部と、入出力する情報を表示する画面を有して前記異常検出部で検出した前記プラント又は設備の異常に関する情報を前記画面上に表示する入出力部とを備えて構成した。
【発明の効果】
【００１５】
本発明によれば、変化の激しいところで機器制限のため取得できないセンサデータを密に推定することによって、そこに発生する異常を捕らえることを可能にしたため、変化の激しいシーケンス時に発生する異常を検知できるようになった。
【００１６】
さらに本発明によれば、取得できないセンサデータを推定することにより、センサのサンプリング時刻が、シーケンスの開始と同期していない場合に発生する異なる時間で取得した同じシーケンスのセンサデータ間にある時間ずれ、および多次元センサを同期してデータを取得できない場合に発生する異なるセンサのセンサデータ間にある時間ずれを解消することを可能にしたため、シーケンス期間中の任意の時刻におけるセンサ波形の統計確率的特性を監視することができるようになった。
【００１７】
以上により、ガスタービンや蒸気タービンなどの設備のみならず、水力発電所での水車、原子力発電所の原子炉、風力発電所の風車、航空機や重機のエンジン、鉄道車両や軌道、エスカレータ、エレベータ、そして機器・部品レベルでは、搭載電池の劣化・寿命など、種々の設備・部品において高感度な異常検知および容易な異常説明を両立するシステムが実現できるようになった。
【図面の簡単な説明】
【００１８】
【図１Ａ】本発明の設備状態監視システムの概略の構成を示すブロック図である。
【図１Ｂ】学習時の処理の流れを示すフロー図である（１クラス）。
【図１Ｃ】学習時の処理の流れを示すフロー図である（多クラス）。
【図１Ｄ】異常検知時の処理の流れを示すフロー図である。
【図１Ｅ】センサデータ推定時刻決定処理の流れを示すフロー図である。
【図２Ａ】センサ波形における開始・終了時のイメージ図である。
【図２Ｂ】センサデータ切出し開始・終了時刻の決定処理の流れを示すフロー図である。
【図２Ｃ】センサデータ切出し開始・終了時刻の決定用指標を示す図である。
【図３Ａ】イベントデータの例を示すフローである。
【図３Ｂ】イベントデータを受けて時刻合わせを行い処理のイメージ図である。
【図３Ｃ】イベントデータを受けて時刻合わせを行い別処理のイメージ図である。
【図４】センサデータ推定時刻決定部における処理の流れを示すフローである。
【図５Ａ】センサデータ推定処理を説明する図である。
【図５Ｂ】センサデータ推定の別処理を説明する図である。
【図５Ｃ】サンプリング点における補正処理を説明する図である。
【図６Ａ】確率分布推定処理を説明する図である。
【図６Ｂ】確率分布推定の別処理を説明する図である。
【図７Ａ】特徴量抽出処理の流れを示すフローである。
【図７Ｂ】特徴量を示す図である。
【図８Ａ】センサデータ収束判定処理の流れを示すフローである。
【図８Ｂ】センサデータ収束判定指標を示す図である（一定値に収束）。
【図８Ｃ】センサデータ収束判定指標を示す図である（一定値を中心に振動）。
【図９Ａ】異常検知結果を表示するＧＵＩ図である（２Ｄ表示）。
【図９Ｂ】異常検知結果を表示するＧＵＩ図である（３Ｄ表示）。
【図１０】シーケンスの切出し・センサデータの推定手法等設定のためのＧＵＩである。
【図１１Ａ】センサデータ推定前後計測曲線をチェックするＧＵＩ図である。
【図１１Ｂ】センサデータ推定後曲線・センサモデル・統計確率分布をチェックするＧＵＩ図である。
【図１２】運転モードの違いによって特徴空間上における局所空間の変化を示す図である。
【発明を実施するための形態】
【００１９】
本発明では、プラントなどの設備において、変化の激しい起動時や停止時のシーケンスを実行しているときに発生する設備の不具合或いはその兆候を検知する設備状態監視方法及び設備状態監視装置に関するもので、シーケンスの開始時刻を合わせ、間欠的に出力されたセンサデータの推定時刻を決定し、その時刻におけるセンサデータを推定することにより、時系列推移を考慮した時刻ごとの確率分布に基づいて異常検知を行うものである。
以下に、本発明の実施例を、図を用いて説明する。
【実施例１】
【００２０】
図１Ａに、本実施例における設備状態監視方法を実現するシステムの一構成例を示す。
本システムは、設備１０１から、または、データベース１１１から出力されるサンプリングセンサデータ１００２とイベントデータ１００１およびユーザからのユーザ指示１００３を受けて異常を検知する異常検知システム１０、途中結果や異常検知結果を記憶する記憶媒体１１、および途中結果や異常検知結果を表示する表示装置１２を備えて構成される。
【００２１】
また、異常検知システム１０は、データを処理するデータ前処理部１０２、データベース１１１からのセンサデータ１００２およびイベントデータ１００１をデータ前処理部１０２にて処理した後、センサデータ推定時刻を決定する推定時刻決定部１１２、そして、設備１０１からのセンサデータ１００２およびイベントデータ１００１をデータ前処理部１０２にて処理した後、センサデータ推定時刻決定部１１２にて決定された時刻におけるセンサデータを推定するセンサデータ推定部１０３、さらにその時刻において統計確率分布推定を行う統計確率分布推定部１０４、統計確率分布を用いて特徴量を抽出する特徴量抽出部１０５、特徴量抽出部１０５で抽出された特徴量を用いて学習する学習部１１３、学習後に学習部１１３から出力される正常空間あるいは識別境界１００４を用いて異常を検知する異常検知部１０６を備えて構成される。
【００２２】
さらに、データ前処理部１０２は、イベントデータ１００１からユーザが指定したシーケンスの開始時刻を探索するイベントデータ解析部１０２１、指定シーケンスの開始時刻の情報を用いて受け取ったセンサデータ１００２からセンササンプリングデータを切出すための開始・終了時刻を算出してセンサデータ１００２の切出しを行うセンサデータ切出し部１０２２、切出したセンサデータの時刻を合わせるセンサデータ時刻合わせ部１０２３を備えている。
【００２３】
また、学習部１１３、識別境界１００４、及び異常検知部１０６で、識別部１０７（１０７´）を構成している。
【００２４】
本システムの動作には、データベース１１１に蓄積されたデータを用いてセンサデータ推定時刻を決定する「推定時刻決定」フェーズ、蓄積されたデータを用いて異常検知に用いる正常空間あるいは識別境界１００４を決定する「学習」フェーズ、推定時刻において補正を実施した後の入力センサデータと正常空間あるいは識別境界に基づき実際に異常検知を行う「異常検知」フェーズの三つのフェーズがある。基本的に、前の二つのフェーズ「推定時刻決定」フェーズと「学習」フェーズとはオフラインの処理、最後の三つ目のフェーズ「異常検知」はオンラインの処理である。ただし、「異常検知」をオフラインの処理とすることも可能である。以下の説明では、それらを「推定時刻決定時」、「学習時」、「異常検知時」という言葉で区別する。
【００２５】
また、図１Ａにおける実線の矢印１００は「異常検知」フェーズにおけるデータの流れを示す「異常検知」パス、点線の矢印１００’は「学習」フェーズにおけるデータの流れを示す「学習」パス、破線の矢印１００”は「推定時刻決定」フェーズにおけるデータの流れを示す「推定時刻決定」パスである。
【００２６】
状態監視の対象とする設備１０１は、ガスタービンや蒸気タービンなどの設備やプラントである。設備１０１は、その状態を表すセンサデータ１００２とイベントデータ１００１を出力する。
【００２７】
本実施例においては、先ずオフラインで「推定時刻決定」フェーズの処理を行い、次に、「推定時刻決定」フェーズの処理の結果を用いてオフラインで「学習」フェーズの処理を行う。そして、「推定時刻決定」フェーズの処理の結果と「学習」フェーズの結果を用いて、オンライン処理による「異常検知」フェーズを実行する。
【００２８】
センサデータ１００２は、設備１０１に取り付けられた複数のセンサから、一定間隔毎に取得される多次元時系列データである。センサの種類は、設備やプラントの大きさにもよるが、数百から数千と言った数になる場合もあり、例えば、シリンダ、オイル、冷却水などの温度、オイルや冷却水の圧力、軸の回転速度、室温、運転時間などである。出力や状態を表すのみならず、何かをある値に制御するための制御データの場合もある。
【００２９】
「推定時刻決定時」における処理の流れについて、図１Ｅを用いて説明する。本処理は、「推定時刻決定」パス１００”に沿ってデータベース１１１から抽出されたイベントデータ１００１とセンサデータ１００２とを用いて実施する。
【００３０】
具体的には、データ前処理部１０２のイベントデータ解析部１０２１にて、データベース１１１から出力されたイベントデータ１００１およびユーザ指示１００３を入力し(Ｓ１３１)、ユーザ指示１００３により指定されたシーケンスの開始時刻を入力したイベントデータ１００１の中から探索する(Ｓ１３２)。
【００３１】
一方、センサデータ切出し部１０２２にて、データベース１１１から出力されたセンサデータ１００２を入力し(Ｓ１３４)、イベントデータ解析部１０２１から得られたシーケンス開始時刻に対応するセンサデータ切出し開始・終了時刻を算出し、データベース１１１から入力したセンサデータ１００２からセンサデータを切出す(Ｓ１３５)。
【００３２】
次に、切り出されたセンサデータは、センサデータ時刻合わせ部１０２３に送られ、センサデータ時刻合わせ部１０２３において時刻合わせが行われ(Ｓ１３６)、推定時刻決定部１１２に送られてセンサデータ推定時刻が決定され(Ｓ１３７)、この決定された推定時刻が保存または出力される(Ｓ１３８)。
【００３３】
「学習時」の処理の流れについて、図１Ｂと図１Ｃを用いて説明する。本処理は「推定時刻決定」パス１００’に沿ってデータベース１１１から抽出されたイベントデータ１００１とセンサデータ１００２とを用いて実施する。
【００３４】
図１Ｂに示しているのは１クラス識別器１０７を利用する場合の学習、図１Ｃに示しているのは多クラス識別器１０７´を利用する場合の学習である。
【００３５】
図１Ｂまたは図１Ｃにおいて、まず、イベントデータ解析部１０２１は、データベース１１１から出力されたイベントデータ１００１およびユーザ指示１００３を入力し(Ｓ１０１)、ユーザ指示１００３により指定されたシーケンスの開始時刻を入力したイベントデータ１００１の中から探索する(Ｓ１０２)。
【００３６】
一方、センサデータ切出し部１０２２は、データベース１１１から出力されたセンサデータ１００２を入力し(Ｓ１０４)、イベントデータ解析部１０２１から得られたシーケンス開始時刻に対応するセンサデータ切出し開始・終了時刻を算出し、データベース１１１から入力したセンサデータ１００２からセンサデータを切出す(Ｓ１０５)。この切り出されたセンサデータは、センサデータ時刻合わせ部１０２３にて、時刻合わせが行われる(Ｓ１０６)。
【００３７】
次に、得られた時刻合わせ後のセンサデータを用いて学習を行う。推定時刻決定部１１２から出力されたセンサデータ推定時刻をセンサデータ推定部１０３へ入力し(Ｓ１０３)、センサデータ推定部１０３にて入力したセンサデータ推定時刻の情報に基づいてセンサデータの時刻を推定する(Ｓ１０７)。次に、統計確率分布推定部１０４にて、時刻が推定されたセンサデータの統計確率分布を推定し(Ｓ１０８)、推定された統計確率分布に基づき特徴量抽出部１０５にて推定したセンサデータの特徴量を抽出する(Ｓ１０９)。
【００３８】
流れの最後に、図１Ｂの１クラスの識別器１０７を利用する場合には、特徴量抽出部１０５で抽出したセンサデータの特徴量を用い、識別器１０７の学習部１１３で学習をすることにより正常空間を作り(Ｓ１１０)、この作った正常空間を出力する(Ｓ１１１)。
【００３９】
一方、図１Ｃの多クラスの識別器１０７´を利用する場合、ユーザ指示１００３によってデータベース１１１から読み込んだセンサデータの正常か異常かのインディックスがあるファイルを入力してセンサデータが正常か異常かを教示し(Ｓ１１２)、その後、特徴量抽出部１０５で抽出した特徴量を用い、識別器１０７´の学習部１１３で学習をすることにより正常か異常かを識別する識別境界１００４を決定し(Ｓ１１０’)、この決定した識別境界１００４を出力する(Ｓ１１１’)。
【００４０】
次に、新たに観測されたセンサデータの「異常検知時」の処理の流れについて、図１Ｄを用いて説明する。本処理は「異常検知」パス１００に沿って設備１０１から抽出されたイベントデータ１００１とセンサデータ１００２とを用いて実施する。まず初めに、イベントデータ解析部１０２１は、設備１０１から出力されたイベントデータ１００１およびユーザ指示１００３を入力し(Ｓ１２１)、ユーザにより指定されたシーケンスの開始時刻を探索する(Ｓ１２２)。
【００４１】
一方、センサデータ切出し部１０２２は、設備１０１から出力されたセンサデータ１００２を入力し(Ｓ１２４)、イベントデータ解析部１０２１から得られたシーケンス開始時刻に対応するセンサデータ切出し開始終了時刻を算出し、センサデータを切出し(Ｓ１２５)、センサデータ時刻合わせ部１０２３にて、時刻合わせを行う(Ｓ１２６)。
【００４２】
次に、学習時予め推定時刻決定部１１２にて決定して保存しておいたセンサデータ推定時刻をセンサデータ推定部１０３に入力し(Ｓ１２３)、センサデータ推定部１０３にて推定時刻決定部１１２から入力したセンサデータ推定時刻においてセンサデータ時刻合わせ部１０２３から入力した時刻合わせしたセンサデータについてセンサデータ推定を行い(Ｓ１２７)、統計確率分布推定部１０４にて推定したセンサデータの統計確率分布を推定し(Ｓ１２８)、推定された統計確率分布に基づき特徴量抽出部１０５にて特徴量を抽出する(Ｓ１２９)。
【００４３】
最後に、特徴量抽出部１０５で抽出した特徴量を用い、識別器１０７（１０７´）において学習部１１３で作成した正常空間又は識別境界１００４を用いて異常検知部１０６で異常識別を行い (Ｓ１３０)、検知結果を出力・表示する(Ｓ１３１)。
【００４４】
次に、図１Ａに示した各部の動作について順に詳細に説明する。すなわち、センサデータ切出し部１０２２における切出し開始・終了時刻の決定、センサデータ時刻合わせ部１０２３におけるセンサデータの時刻合わせ、推定時刻決定部１１２におけるセンサデータ推定時刻の決定、センサデータ推定部１０３におけるセンサデータの推定、統計確率分布推定部１０４における確率分布の推定、特徴量抽出部１０５における特徴量抽出について、図２Ａないし図８Ｃを用いて、詳しく説明する。
【００４５】
［切出し開始・終了時刻の決定］  
センサデータ切出し部１０２２において、まず、センサデータ切出し開始・終了時刻を算出し、そして、その切出し開始・終了時刻を用いて、その間にあるセンサデータを切出す。
【００４６】
図２Ａは取得したセンサ波形における切出し開始・終了時のイメージを示す図である。図２Ａの（ａ）と（ｂ）の例は、切出しを開始してから終了するまでの間にセンサ波形の立ち上がりと下がりとの両方を含み、開始・終了時におけるセンサデータ値は同じレベルとなるものである。（ａ）の例では、立ち上がりと下がりの間でセンサ波形は滑らかに変化しているのに対し、（ｂ）の例では、立ち上がりと下がりの間で波形がギザギザとなっている。また（ｃ）の起動シーケンスのみがある例、あるいる（ｄ）の停止シーケンスのみがある例は、切出し開始・終了時のセンサデータ値が別々のレベルになるものである。
【００４７】
次に、図２Ｂおよび図２Ｃを用いてセンサデータ切出しを行うため切出し開始・終了時刻を算出する処理の流れおよび切出し開始・終了判定指標について説明する。
【００４８】
図２Ｂはセンサデータ切出し開始・終了判定の流れを示す図である。センサデータ切出し部１０２２は、まず、ユーザ指示１００３を入力し(Ｓ２０１)、ユーザ指示により、モード（開始or終了）の算出は自動か非自動かを決定する(Ｓ２０２)。次に、イベントデータ解析部１０２１から得られた指定シーケンスの開始時刻を入力する(Ｓ２０３)。また、設備１０１、あるいは、データベース１１１から出力されたセンサデータ１００２を入力する(Ｓ２０４)。Ｓ２０３からの指定シーケンスの開始時刻、Ｓ２０４からのセンサデータ、Ｓ２０２からの開始モードの自動か非自動かの決定結果を受け、切出し開始時刻の計算を開始する(Ｓ２０５)。
【００４９】
切出し開始時刻の計算は自動的に実施する場合に、窓で一部センサデータを切出し(Ｓ２０６)、開始判定指標を算出し(Ｓ２０７)、開始判定を行い(Ｓ２０８)、「Ｎｏ」と判定された場合、時刻が増大方向へ窓を移動し(Ｓ２０９)、開始判定(Ｓ２０６〜Ｓ２０８)を繰り返す。「Ｙｅｓ」と判定された場合、センサデータ切出し開始時刻を出力・保存する(Ｓ２１１)。
【００５０】
一方、自動算出を実施しない場合には、指定シーケンスの開始時刻をセンサデータ切出し開始時刻とし(Ｓ２１０)、センサデータ切出し開始時刻を出力する(Ｓ２１１)。
【００５１】
センサデータ切出し開始時刻を算出した後、センサデータ切出し終了時刻の計算を実施する。センサデータ切出し終了時刻の計算は、Ｓ２１１からの切出し開始時刻、Ｓ２０２から終了モードの自動か非自動かの決定結果を受け、切出し終了時刻の計算を開始する(Ｓ２１２)。
【００５２】
切出し終了時刻の計算を自動的に実施する場合には、切出し開始時刻以後のセンサデータを対象とし、窓で一部切出し(Ｓ２１３)、終了判定指標を算出し(Ｓ２１４)、終了判定を行い(Ｓ２１５)、「Ｎｏ」と判定された場合、時刻が増大方向へ窓を移動し(Ｓ２１６)、終了判定(Ｓ２１３〜Ｓ２１５)を繰り返す。「Ｙｅｓ」と判定された場合、センサデータ切出し終了時刻を出力・保存する(Ｓ２１８)。
【００５３】
一方、自動算出を実施しない場合には、センサデータ切出し開始時刻以降の所定センサデータ数の時刻をセンサデータ切出し終了時刻とし(Ｓ２１７)、センサデータ切出し終了時刻を出力する(Ｓ２１８)。
【００５４】
開始・終了判定指標について、図２Ｃに一例を示す。この例では、前後の２個のセンサデータを直線で繋ぎ、その直線の傾きを開始・終了判定指標とする。指標が所定閾値より大きくなったところの時刻をセンサデータ切出し開始時刻とする。また、指標が所定閾値より小さくなったところの時刻をセンサデータ切出し終了時刻とする。
【００５５】
［センサデータの時刻合わせ］  
センサデータの時刻合わせ部１０２３での処理は、センサデータ切出し部１０２２から得られた切出し開始時刻を用いて行われる。
【００５６】
イベントデータ１００１の例を図３Ａに示す。イベントデータ１００１は、不定期に出力される設備の操作・故障・警告を表す信号であり、時刻と操作・故障・警告を表すユニークなコードとメッセージ文字列からなる。例えば、起動シーケンスの開始時および停止シーケンスの開始時に対応しているメッセージ文字列はそれぞれ「Request module on」および「Request module off」である。また、イベントデータ１００１において、同じ指定シーケンスが異なる時間で行われるため、その開始時刻は複数個ある。
【００５７】
センサデータ時刻合わせ部１０２３で実行するイベントデータ１００１を用いてセンサデータ１００２の時刻合わせ処理の第一の例を図３Ｂに示す。図には、(a)時刻合わせ前のセンサデータと(b)時刻合わせ後のセンサデータを示す。図示の（a）のように、異なる時間で取得した同指定シーケンスのセンサデータを算出したそれぞれの切出し開始時刻で経過時刻を算出し、（ｂ）に示すように、ゼロ時刻を合わせるように切出したセンサデータを同一時間軸上に時刻を合わせる。開始時刻からの経過時間の時間間隔は一定間隔ではなくてもよい。または、一番短い時間間隔に合わせて一定間隔にしてもよい。（ｂ）の表において、時刻合わせ後のセンサデータの一覧に、表示されている数字は取得したセンサデータ、空白となっているのは取得できないセンサデータである。
【００５８】
センサデータ時刻合わせ部１０２３で実行するイベントデータ１００１を用いたセンサデータの時刻合わせ処理の第二の例を図３Ｃに示す。本例では、図示のように、（ｂ）の修正センサデータ列の切出し開始時刻ts、correctと切出し終了時刻te、correctを（ａ）の参照センサデータ列の切出し開始時刻ts、refと切出し終了時刻te、ref合わせるように、参照センサデータ列の時間間隔Δtrefから、（ｂ）の修正センサデータ列の時間間隔Δt’correctを
【００５９】
【数１】
2014048697.tif 000003
【００６０】
により変更して、（ｃ）の時刻合わせ後の修正センサデータ列を得る。
【００６１】
［センサデータ推定時刻の決定］  
推定時刻決定部（１１２）で行われるセンサデータ推定時刻の決定処理の流れを、図４を用いて説明する。データベース１１１からの学習用正常なセンサデータをデータ前処理部１０２で処理して得られた時刻合わせた後のセンサデータ列をセンサデータの推定時刻決定部１１２に入力し(Ｓ４０１)、窓（ウィンドウ）で一部センサデータを切出し(Ｓ４０２)、激しさ評価指標を算出し(Ｓ４０３)、激しさ評価指標とサンプリング間隔との関係式を用いて激しさ評価指標からサンプリング間隔を算出し(Ｓ４０５)、終了するかとの判定を行う(Ｓ４０６)。「Ｎｏ」と判定された場合、時刻が増大方向へ窓を移動し(Ｓ４０７)、センサデータの激しさ算出からサンプリング間隔を算出する処理 (Ｓ４０２〜Ｓ４０５)を繰り返す。「Ｙｅｓ」と判定された場合、サンプリング間隔を用いて窓内の推定時刻を算出し(Ｓ４０８)、推定時刻を保存・出力する(Ｓ４０９)。
【００６２】
本発明では、時系列データの激しさ評価指標は、時系列波形の周波数の高低、および、時系列波形の上昇・降下の落差の大きさにより定量化されるものとして定義する。すなわち、時系列波形の周波数が高い、または、時系列波形の上昇・降下の落差が大きい場合、激しさが大きい、逆に、時系列波形の周波数が低い、または、時系列波形の上昇・降下の落差が小さい場合、激しさが小さい。
【００６３】
具体的に、例えば、窓で切出した一部のデータに対し、フーリエ解析を行い、パワースペクトルを計算し、最大パワースペクトルの周波数値をそのデータ列の周波数とし、ある最大周波数で正規化されたデータ列の周波数を周波数による激しさIfreq（freq）とする。同じデータに対し、前後データの差分の最大値に対し、ある最大データの差分で正規化してデータの差分による激しI｜Δｙ｜（｜Δｙ｜）とする。ある周波数やあるデータの差分は、例えば、すべてのセンサデータに対して統計的に計算して得られた最大値を利用しても良い。これに限らない。また、次に、このデータ列の激しさは、下式で計算される。
【００６４】
【数２】
2014048697.tif 000004
【００６５】
激しさ評価指標はその他の定義を使っても良い。  
激しさ評価指標とサンプリング間隔との関係式は別途に予め実験やシミュレーションを通して求めておく(Ｓ４０４)。図示のように、サンプリング間隔の最大値はデータ取得のサンプリング間隔であり、最小値は1秒である。激しさ評価指標とサンプリング間隔との間には反比例の関係を持つ。
【００６６】
センサデータ推定時刻の決定は予め決められた一定間隔でも良く、指定した個数推定するよう等間隔に決めても良い。
【００６７】
以上のように、センサデータ推定時刻の決定を行うことにより、処理コストをダウンし、高効率に処理を実施することが可能となる。
【００６８】
［センサデータの推定］  
センサデータ推定部でのセンサデータの推定（推定センサデータの算出）について、図５Ａ、図５Ｂと図５Ｃを用いて説明する。推定センサデータは同じ取得センサデータ列の取得センサデータおよび推定センサデータの近傍における同指定シーケンスの異なる時間で取得した他の取得センサデータから重み付き加算の形で算出される。
【００６９】
図５Ａはセンサデータ推定の第一の例を示す。第一の例では、各取得センサデータ間のセンサデータ推定値はその両側にある取得センサデータから線形的に算出される。y(x)は取得できないデータの推定値、yxは取得したセンサデータ値、jは時刻合わせた後データが0秒から取得時のサンプリング間隔でカウントしたサンプリング番号(j=1〜n)、iは異なる時間で取得した同じ指定シーケンスの番号(i=1〜m)とすると、y(x)は
【００７０】
【数３】
2014048697.tif 000005
【００７１】
により算出される。
【００７２】
また、図５Ｂで示す第二の例では、推定センサデータと同じ取得センサデータ列のすべての取得センサデータを用いて非線形的に推定センサデータを算出する。センサデータの推定値y(x)は、
【００７３】
【数４】
2014048697.tif 000006
【００７４】
と表わされる。ここで、αは重み係数であり、αは下記のように高次元写像されたxにより求められる。
【００７５】
【数５】
2014048697.tif 000007
【００７６】
高次元写像関数は下式
【００７７】
【数６】
2014048697.tif 000008
【００７８】
を用いる。λは経験により決めた係数である。
【００７９】
また、β1iとβ2iは重み係数であり、これらは周辺取得センサデータの分散から計算される。
【００８０】
センサデータの推定には、スプライン法やバイキュービック法等もある。いくつかの手法のどちらか一方を使っても良いし、切り替えても良い。切り替えには、例えば激しさ指標を使う。
【００８１】
取得センサデータで区切られたそれぞれの区間におけるセンサデータを推定する際に、用いられる手法が異なる場合に、手法の切り替わるところにおいて、垂直な段差のようなデータのずれが発生する。
【００８２】
図５Ｃはサンプリング点における段差を補正する一例を示す。この例では、推定line1と推定line2の点xj における段差に対し、サンプリング点の前後に補正スペース2l (１≦データ取得時のサンプリング間隔)を設け、補正スペースxj-l からxj+l 内で、補正曲線ｙ´（ｘ）ように二つのセンサデータを線形的につなげる。すなわち、推定した二つのセンサデータの境界（繋ぎ目）に発生する垂直な段差を斜めの段差に変更することによって、段差のような不連続な繋がりから、滑らかな繋がりにすることである。
【００８３】
このように、補正スペース内での補正センサデータy’ (x) は下式で求められる。
【００８４】
【数７】
2014048697.tif 000009
【００８５】
重み係数w(x) は
【００８６】
【数８】
2014048697.tif 000010
【００８７】
により計算される。
【００８８】
以上のように、センサデータ推定を行うことにより、機器制限のために取得できないデータを推定し、特に、実際にデータを取得できないシーケンスの激しい変化を再現することが可能となる。
【００８９】
［統計確率分布の推定］  
統計確率分布推定部１０４における統計確率分布の推定について、各推定時刻において、複数の同指定シーケンスの異なる時間で取得したセンサデータの推定後の値を用い、各推定時刻における確率分布の推定方法について、図６Ａと図６Ｂを用いて説明する。
【００９０】
図６Ａに示している例は、通常に各推定時刻におけるセンサデータが正規分布に従う場合の確率分布Ｇである。この場合に、確率分布Ｇはガウス関数で表され、その時刻におけるセンサデータの平均μと標準偏差σを用いて下式のようになる。
【００９１】
【数９】
2014048697.tif 000011
【００９２】
一方、図６Ｂに示した例は、各推定時刻におけるセンサデータが正規分布に従わない場合の確率分布Ｇの一例である。この場合に、例えば、マルチガウス関数で近似しても良い。その他の関数で近似しても良い。マルチガウス関数で近似する場合には、下式のようになる。
【００９３】
【数１０】
2014048697.tif 000012
【００９４】
以上の統計確率分布Ｇの推定により、各時刻におけるセンサデータの分布状況を把握することが可能となる。また、各時刻における新たに観測されたセンサデータについて、正常・異常の割合がどのぐらいあるかを見て確認することが可能となる。
【００９５】
［特徴量抽出］  
特徴量抽出部１０５で実施する特徴量抽出処理の流れを図７Ａを用いて説明する。まず、統計確率分布推定部１０４からの各推定時刻における統計確率分布Ｇを入力する(Ｓ７０１)。
【００９６】
次に、各推定時刻における統計確率分布Ｇを用いてその時刻の異常度ｖ（ｔ）を(数１０)を用いて算出する(Ｓ７０２)。
【００９７】
【数１１】
2014048697.tif 000013
【００９８】
そして、後で説明する特徴量抽出部１０５において実施するセンサ波形の収束判定から得られたシーケンス収束時刻を入力し(Ｓ７０３)、センサデータの開始時刻からシーケンスの収束時刻まで(数１１)で累積することで特徴量である尤度を算出する(Ｓ７０４)。
【００９９】
【数１２】
2014048697.tif 000014
【０１００】
すべてのセンサに対し、Ｓ７０１〜Ｓ７０３を実施し、最後に、すべてのセンサの尤度を統合すれば、(数１２)で表す尤度ヒストグラムが得られる。添字s1〜snはセンサ番号である。
【０１０１】
【数１３】
2014048697.tif 000015
【０１０２】
図７Ｂには抽出した特徴量である尤度ヒストグラムを示す。  
図７Ａの処理Ｓ７０３にて入力したシーケンス収束時刻は、センサデータを切出す時に、ユーザ指示により「自動算出を実施しない」ことを選択した際、デフォルトの所定数まで切出したセンサデータを対象として、センサ波形の収束判定を行うことにより求めたものである。
【０１０３】
図８Ａには、その収束時刻を求める処理の流れを示す。まず、切出した正常サンプリングセンサデータを入力する(Ｓ８０１)。次に、窓で一部センサデータを切出し(Ｓ８０２)、収束判定指標を算出し(Ｓ８０３)、収束判定を行う(Ｓ８０４)。「Ｎｏ」と判定された場合、時刻が増大方向へ窓を移動し(Ｓ８０５)、収束判定(Ｓ８０２〜Ｓ８０４)を繰り返し、「Ｙｅｓ」と判定された場合、センサデータ収束時刻を出力する(Ｓ８０６)。
【０１０４】
シーケンス収束時とは、シーケンスが開始してから、該シーケンスのセンサデータが一定値に収束し始める時刻、あるいは、一定値を中心に振動し始める時刻のことである。前者の収束判定指標を示すイメージは図８Ｂとなり、後者の収束判定指標を示すイメージは図８Ｃとなる。
【０１０５】
図８Ｂに一定値に収束するケースの収束判定指標を示す。収束判定指標は窓で切出したサンプリングセンサデータを対象とする主成分分析により得られた第一主軸、あるいは線形回帰により得られた回帰直線の傾きである。また、最後の一部のサンプリングデータの［最小値、最大値］範囲内に到達し、最大値と最小値との差分が所定閾値以下との制約条件を加え、収束判定指標が所定の閾値より小さいいくつかの時刻のうち、最初の時刻はシーケンス収束時刻とする。
【０１０６】
図８Ｃに一定値を中心に振動するケースの収束判定指標を示す。この場合の収束判定指標も窓で切出したサンプリングセンサデータを対象とする主成分分析により得られた第一主軸の傾き（図８Ｃに示したように、第一主軸が水平軸とのなす角度）である。また、最後の一部のサンプリングデータの頂点に余弦波で合わせた後、類似度を算出し、類似度が所定閾値以上に達することを制約条件にさらに加え、収束判定指標が所定の閾値より小さくなったいくつかの時刻のうち、最初の時刻をシーケンス収束時刻とする。
【０１０７】
［ＧＵＩ］  
次に、各処理をすることにあたって使用するＧＵＩについて、図９Ａないし図１１Ｂを用いて説明する。
【０１０８】
図９Ａ及び図９Ｂは、図１Ｂ〜Ｄに示したフロー図におけるイベントデータおよびユーザ指示を入力する処理ステップＳ１０１(Ｓ１２１)、学習センサデータを入力する処理ステップＳ１０４(Ｓ１２４)、図１Ｃにおける正常・異常指示を入力する処理ステップＳ１１２、図１Ｂ又は図１Ｃの出力処理ステップである正常空間あるいは正常・異常識別境界を出力する処理ステップＳ１１１(Ｓ１１１’)、図１Ｅに示したフロー図における推定時刻を出力する処理ステップＳ１３８、図２Ｂに示したフロー図における切出し開始・終了時刻を出力する処理ステップＳ２１１、Ｓ２１８、図４に示したフロー図における推定時刻を出力する処理ステップＳ４０９、図７Ａに示したフロー図における特徴量を出力する処理ステップＳ７０５、図８Ａに示したフロー図におけるセンサデータ収束時刻を出力する処理ステップＳ８０６、図１Ｄに示したフロー図における異常検知結果を出力する処理ステップＳ１３１に係わるＧＵＩ、また、異常検知テスト結果や異常検知結果の表示に係わるＧＵＩである。
【０１０９】
本ＧＵＩは、特徴量を表示するパネル９００、センサデータ、センサデータの正常異常を示すインディクス、イベントデータ、パラメータ等を保存するファイル一式を格納するフォルダを選択するための「参照」ボタン９０１２、選択したフォルダを表示するための「入力フォルダ」ボックス９０１１、Ｓ１１１(Ｓ１１１’)処理ステップから受けた正常空間(正常・異常識別境界)、Ｓ１３８とＳ４０９処理ステップから受けた決定した推定時刻、Ｓ２１１とＳ２１８処理ステップから受けた切出し開始・終了時刻、Ｓ７０５処理ステップから受けた特徴量である尤度ヒストグラム、Ｓ８０６処理ステップから受けたセンサ波形の収束時刻、また、図示していないが、Ｓ１０７(Ｓ１２７)処理ステップから受けた推定後のセンサデータ、Ｓ１０８(Ｓ１２８)処理ステップから受けた抽出した統計確率分布等の途中結果およびＳ１３１処理ステップから受けた異常検知結果のファイル一式を格納するフォルダを選択するための「参照」ボタン９０２２、選択したフォルダを表示するための「出力フォルダ」ボックス９０２１、今回の学習および異常検知テストに係わるデータを登録するための「データ期間登録」登録ボックス９０３、検知手法を選択するための「異常検知手法」選択ボックス９０４、異常検知のための詳細な設定に係わる「その他の設定」ボタン９０５、データベース１１１からのデータを用いて学習および異常検知テストを実行する「学習および異常検知テスト実行」ボタン９０６、設備１０１からのデータを対象として異常検知をする「異常検知実行」ボタン９０７、異常検知結果を表示するデータの「表示期間」ボックス９０８、特徴量や異常結果の表示等の表示項目を選択する「表示項目」ボックス９０９、２次元表示をするか３次元表示をするかを選択するための「表示形式」ボックス９１０、以上の表示の設定に基づき異常検知を実施し、センサデータを推定処理ステップＳ１０７(Ｓ１２７)から受けた異常検知結果を表示する「異常検知結果表示」ボタン９１１、および、統計確率分布ステップＳ１０８(Ｓ１２８)から受けた途中結果である推定後センサデータおよび統計確率分布を受け、表示するための「途中結果表示」ボタン９１２を備えて構成される。
【０１１０】
「入力フォルダ」ボックス９０１１、「出力フォルダ」ボックス９０２１でそれぞれフォルダを選択し、「データ期間登録」ボックス９０８でデータ期間を登録し、「異常検知手法」選択ボックス９０４で異常検知手法を選択し、「その他の設定」ボタン９０５でその他の設定が終わった後、「学習および異常検知テスト実行」ボタン９０６を押すことにより、図１Ｂ、あるいは、図１Ｃに示した学習処理、およびデータベース１１１からのデータを用いて図１Ｄに示した異常検知流れに沿った異常検知テスト処理が実行される。ボタン９０６が一旦押されると、異常検知テスト処理が完了するまでは、「異常検知実行」ボタン９０７、「異常検知結果表示」ボタン９１１、「途中結果表示」ボタン９１２が押せない状態になる。
【０１１１】
学習処理および異常検知テスト処理が完了したら、「異常検知実行」ボタン９０７が押せる状態に切り替わる。また、「異常検知結果表示」ボタン９１１、「途中結果表示」ボタン９１２も押せる状態に切り替わる。この場合に、「表示期間」ボックス９０８に学習データ、あるいは、異常検知テストデータの期間を登録し、「表示項目」ボックス９０９および「表示形式」ボックス９１０を選択して「異常検知結果表示」ボタン９１１、あるいは「途中結果表示」ボタン９１２を押すと、表示期間にある途中結果、あるいは異常検知結果が「表示」パネル９００に表示される。
【０１１２】
次に、「異常検知実行」ボタン９０７を押す。これにより、図示していない設備１０１に繋がっている一時データ保存用の記憶媒体から「データ期間登録」ボックス９０３に登録されている期間のデータを読み込んで異常検知を行う。異常検知実行が完了したら、「表示期間」ボックス９０８に異常検知データの期間を登録し、「表示項目」ボックス９０９および「表示形式」ボックス９１０を選択して「異常検知結果表示」ボタン９１１、「途中結果表示」ボタン９１２を押すと、表示期間にある異常検知データの途中結果、あるいは異常検知結果が「表示」パネル９００に表示される。
【０１１３】
表示に関連するボタンを押す前には、「表示」パネル９００には実行の進捗が表示されている。たとえば、最初に、“設定してください”と表示する。設定がいったん始まったら、すぐに、“設定実施中”と切り替わる。そして、入力フォルダ、出力フォルダ、データ期間登録、異常検知手法の設定が完了し、「学習および異常検知テスト実行」ボタン９０６を押すと、“学習および異常検知テスト実行中”と表示する。
【０１１４】
「学習および異常検知テスト実行」が完了したら、“学習および異常検知テスト実行が済み、「異常検知実行」ボタン９０７を押して異常検知を実施するか、表示設定をし、表示ボタンを押して、表示してください”と表示する。「異常検知実行」ボタン９０７を押さず、「表示期間」、「表示項目」、「表示形式」のいずれかの表示設定をいったん開始したら、“表示設定中”と切り替える。表示設定が完了したら、“表示設定が完了し、表示ボタンを押して、表示してください”と表示する。
【０１１５】
「異常検知結果表示」ボタン９１１か「途中結果表示」ボタン９１２を押したら、設定に従った学習および異常検知テストの結果が表示される。一方、「異常検知実行」ボタン９０７を押したら、“異常検知テスト実行中”と表示し、「異常検知テスト実行」が完了したら、“異常検知テスト実行が済み、表示設定をしてください”と表示する。表示設定をいったん開始したら、“表示設定中”と切り替える。表示設定が完了したら、“表示設定が完了し、表示ボタンを押して、表示してください”と表示し、「異常検知結果表示」ボタン９１１か「途中結果表示」ボタン９１２を押したら、設定に従った異常検知の結果が表示される。
【０１１６】
図９Ａに本実施例によるＧＵＩ表示の一例を示す。この例では、表示パネル９００において、特徴量９００１、異常バー９００２、および表示に係わる表示関連項目９００３が表示されている。表示関連項目９００３には、表示データの種類（データベースのデータを利用した異常検知テスト結果か、設備からのデータを利用した異常検知の結果かのこと）、表示期間、本結果が得られるための学習期間、評価期間を表示する。異常バー９００２は、異常のある特徴量の位置を黒く表示する。
【０１１７】
図９Ａに示した例では、データベース１１１からのデータを利用した異常検知テスト結果を示す一例である。また、図示をしていないが、設備１０１からのデータを利用した異常検知結果を示すことも可能である。また、「表示形式」に３Ｄと選択すれば、図９Ｂの「表示」パネル９００に示される３Ｄの特徴量９００１’が表示される。
【０１１８】
図９Ａ又は図９Ｂに示したようなＧＵＩを表示することにより、特徴量である尤度ヒストグラムの状況、異常検知の結果を視覚的に確認することができるため、ユーザにとって判断しやすくなる。
【０１１９】
図１０は、異常検知のための詳細な設定に係わり、図９にある「その他の設定」ボタン９０５を押すことで呼び出されるＧＵＩである。本ＧＵＩは、図１Ｂ〜図１Ｄにおける、センサデータ切出し開始・終了時刻算出し、センサデータを切出し処理ステップＳ１０５(Ｓ１２５)、センサデータ推定ステップＳ１０７(Ｓ１２７)の設定にも係わる。
【０１２０】
本ＧＵＩは、「シーケンス設定」領域１００１、「センサデータ推定設定」領域１００２、「データ設定」領域１００３、「識別器設定」領域１００４、「設定状況一覧」表示パネル１００５、および「保存」ボタン１００６を備えて構成される。
【０１２１】
「シーケンス設定」領域１００１において、「編集」ボタン１００１６を押すことで、すべての項目を編集することが可能となる。編集項目には、「シーケンス種類」と「シーケンス切出し」がある。それぞれ、「シーケンス種類」にはシーケンス種類を選択するためのボックス１００１１、「シーケンス切出し」には、シーケンス切出し開始・終了時刻自動算出項目に、「あり」を表示するチェックボックス１００１２１と１００１２３、およびそれぞれの「あり」の後ろに自動算出時使用する指標を選択するためのボックス１００１２２と１００１２４がある。「シーケンス種類」選択ボックス１００１１にて、起動や停止等異常検知したいシーケンスの種類を選択することができる。「シーケンス切出し」にて、開始・終了時刻を自動算出するかを決めることができる。
【０１２２】
自動算出を実施する場合に、それぞれの「あり」チェックボックス１００１２１と１００１２３をチェックし、指標ボックス１００１２２と１００１２４にて、使用指標を選択する。自動算出を実施しない場合に、「あり」チェックボックスをチェックせずに、使用指標選択ボックスを空白にする。この場合には、デフォルトのシーケンス切出し開始・終了時刻を利用する。
【０１２３】
図１０に示した例では、シーケンス種類を選択するためのボックス１００１１に起動シーケンスを入力した例を示している。また、この例では、シーケンス切出し開始・終了時刻の自動算出を行わないため、それぞれの「あり」をチェックせず、指標選択用ボックスを空白にしている。「決定」ボタン１００１７を押すことにより、「シーケンス設定」領域１００１で設定された内容が登録される。
【０１２４】
「センサデータ推定設定」領域１００２において、「編集」ボタン１００２６を押すことで、すべての項目を編集することが可能となる。編集項目には、「推定手法」、「パラメータ」、「推定間隔」がある。それぞれ、「推定手法」には線形法・非線形法・混合法を選択するためのチェックボックス１００２１１、１００２１３、１００２１５、各分類方法に対応する更なる詳細な方法を選択するためのボックス１００２１２、１００２１４、１００２１６がある。
【０１２５】
推定手法を選択する際に、「推定手法」の線形法・非線形法・混合法を選択するためのチェックボックス１００２１１、１００２１３、１００２１５の何れかをチェックし、その後ろに対応する手法選択のボックス１００２１２、１００２１４、１００２１６によって推定手法を決めることができる。「パラメータ」にはパラメータの種類を選択するための選択ボックス１００２２１、選択されたパラメータの具体的な数字を入力するためのボックス１００２２２、および１つのパラメータの種類選択と数値入力が終わった後、次のパラメータの種類選択と数値入力を実施するための「追加」ボタン１００２２３がある。
【０１２６】
「推定間隔」には、「指定あり」の場合に表示するチェックボックス１００２３２、および「指定あり」がチェックされたら推定間隔を入力するためのボックス１００２３３がある。推定間隔を指定しない場合に、「指定あり」チェックボックをチェックせず、その後ろの秒数も入力しない。この場合に、自動的に正常な学習データを用いてそれぞれのセンサ波形の激しさによって推定時刻の決定が行われる。指定する場合には、「指定あり」チェックボックスをチェックし、その後ろの秒数も入力する。これにより、指定された秒数間隔に推定時刻が設定される。
【０１２７】
図１０に示した例では、非線形法のチェックボックス１００２１３がチェックされ、対応する手法選択ボックス１００２１４にカーネルを用いた推定法が選択される。また、「パラメータ」の「種類」選択ボックス１００２２１によってパラメータ１／パラメータ２が選択され、数値が５．０と１０．０と入力されている。推定間隔は指定ありで、推定間隔は１秒と設定される。「決定」ボタン１００２７を押すことにより、「センサデータ推定設定」領域１００２で設定された内容が登録される。
【０１２８】
「データ設定」領域１００３において、「編集」ボタン１００３６を押すことで、すべての項目を編集することが可能となる。編集項目には、「学習／評価データ切わけ指定」および「除外データ」がある。さらに、「学習／評価データ切わけ指定」には、「あり」のチェックボックス１００３１１、指定「あり」の場合に、それぞれ学習データ期間を入力するボックス１００３１２、評価データ期間を入力するボックス１００３１３、指定「なし」の場合に、「なし」のチェックボックス１００３２１、自動的に学習データ・評価データを分けることに係わる評価手法のfold数を入力するボックス１００３２２がある。「除外データ」には、「あり」のチェックボックス１００３３１、と「あり」の場合に、データ登録するための「データ登録」ボックス１００３３２がある。
【０１２９】
図１０に示した例では、学習／評価データ切わけ指定の「あり」のチェックボックス１００３１１がチェックされ、それぞれ学習データボックスと評価データボックス１００３１２と１００３１３に指定期間を入力している。また、除外データがないため、「除外データ」の「あり」のチェックボックス１００３３１をチェックせず、「データ登録」ボックス１００３３２にも空白している。「決定」ボタン１００３７を押すことにより、「データ設定」領域１００３で設定された内容が登録される。
【０１３０】
「識別器設定」領域１００４において、「編集」ボタン１００４６を押すことで、すべての項目を編集することが可能となる。編集項目には、「識別器種類」および「詳細項目」がある。それぞれの内容に対し、「識別器種類」ボックス１００４１、および「詳細項目」ボックス１００４２がある。「識別器種類」ボックス１００４１では、識別器の種類選択ができる。例えば、サポートベクターマシン、ベイズ識別器、ｋ近傍識別器、ニューラルネットワーク等がある。
【０１３１】
「詳細項目」ボックス１００４２では、「識別器種類」ボックス１００４１にて選んだ識別器に対応する詳細な項目の選択ができる。例えば、識別器のクラス数について、１クラスや多クラスを選択できる。１クラスを選択した場合、学習時図１Ｂの処理流れに沿って学習を行い、正常空間を求める。一方、多クラスを選択した場合、学習時図１Ｃの処理流れに沿って学習を行い、正常・異常識別境界を求める。本実施例では、「識別器種類」ボックス１００４１には識別器１と選択し、「詳細項目」ボックス１００４２に多クラスと選択している。「決定」ボタン１００４７を押すことにより、「識別器設定」領域１００４で設定された内容が登録される。
【０１３２】
上記入力した内容が入力された時点で、自動的に「設定状況一覧」１００５に表示される。それぞれの設定項目を編集している場合に、それぞれの項目名の後ろに「編集中」と表示される。決定が終わった時点、それぞれの「決定」ボタン１００１６、１００２６、１００３６、１００４６を押したら、それぞれの項目名の後ろに「編集中」が「決定済」と変更される。また修正したい項目があったら、修正したい設定項目がある領域の「編集」ボタンを押して、編集を行い、編集をし終わった後、再度該当する領域の「決定」ボタン１００１７、１００２７、１００３７、１００４７の何れかを押すことにより修正を完了する。
【０１３３】
「設定状況一覧」１００５に表示された内容を確認し終わった後、「保存」ボタン１００６を押すと、図１０で設定した内容が保存され、図１０のＧＵＩが消える。
【０１３４】
図１０のＧＵＩを用いてセンサデータ推定等を行う詳細な設定を実施した後、図９Ａ又は図９Ｂの「学習＆異常検知テスト実行」ボタン９０６を押して、学習および異常検知テストを実施した後、「異常検知実行」ボタン９０７を押し、異常検知を実施する。異常検知を実施する際に用いられる正常空間あるいは識別境界は「学習時」求められたものを利用する。
【０１３５】
学習、異常検知テスト、異常検知を行った後、得られた途中結果である推定後計測曲線、センサモデル、センサモデルにおけるある時刻の統計確率分布をチェックすることに係わるＧＵＩについて、図１１Ａ及び図１１Ｂを用いて説明する。
【０１３６】
図１１Ａ及び図１１Ｂに示したＧＵＩは、「センサ設定」領域」１１０１、「表示設定」領域１１０２、および表示実行のための「表示」ボタン１１０３、「表示パネル」領域１１０４を備えて構成される。「センサ設定」１１０１には「センサ種類」項目があり、種類選択は選択ボックス１１０１１によって実施される。「表示設定」１１０２には「表示データ日付」、「表示内容」、「確率分布表示」がある。
【０１３７】
表示データ日付は、「表示データ日付」ボックス１１０２１によって入力される。表字内容は「表示内容」の選択ボックス１１０２２１によって選択される。また選択ボックス下にある「設定プロパティ」１１０２２２によって、表示する内容のプロパティを選択する。「確率分布表示」には、「あり」の場合にチェックするチェックボックス１１０２３１、および「あり」の場合に、設定用の「設定プロパティ」１１０２３２がある。
【０１３８】
図１１Ａは、推定前後計測曲線を表示する一例を示す。図１１Ａに示した例では、「表示内容」の選択ボックス１１０２２１に「推定前＆推定後計測曲線」と選択し、またその他の項目に適切な選択肢を選び、「表示」ボタン１１０３を押すと、「表示パネル」領域１１０４に時間とセンサ値との関係を示すグラフ１１０５が表示される。推定前のデータ列１１０５１は離散的であり、推定後のセンサ波形１１０５２は連続的である。
【０１３９】
また、本グラフに係わる設定項目は領域１１０５３に表示される。表示項目内容には、図１０で説明したＧＵＩ上で設定された内容も含み、多次元センサのセンサ番号、計測内容、データ取得時間、データ種類(学習データか評価データか、および、データベースからか設備からか、との内容を含む)、収束時刻、センサデータの推定手法、推定時に用いたパラメータの値、推定時間間隔種類、推定時間間隔、推定前のマーカーの種類、推定後曲線の種類等がある。領域１１０５３の表示はマウスの右クリックにて表示されないように選択できるようにしても良い。
【０１４０】
図１１Ａで、「表示内容」選択ボックス１１０２２２で「センサモデル＆推定後計測曲線」を選択し、またその他の項目に適切な選択肢を選び、「表示」ボタン１１０３を押すと、領域１１０４には図１１Ｂに示すような、センサモデル毎の推定計測曲線を示すグラフ１１０６が表示される。グラフ１１０６において、一点鎖線１１０６１はセンサモデルの平均値曲線（μ）、点線１１０６２、１１０６３はセンサモデルの平均値曲線プラスマイナス標準偏差の３倍（μ±３σ）、実線１１０６４は推定後の計測曲線を表す。また、本図形に係わる設定項目は領域１１０６５に表示される。表示項目内容には、それぞれの線形はどんな曲線であるかを示す。領域１１０６５の表示はマウスの右クリックにて表示されないように選択しても良い。
【０１４１】
さらに、統計確率分布を見たい時刻にマウスを合わせ(グラフ１１０６の矢印で示した個所)、右クリックにて表示するように選択すれば、領域１１０４においてグラフ１１０６の下に、図１１Ｂに示すように、ある指定時刻における統計確率分布１１０７が表示される。ある時刻における統計確率分布１１０７にはガウス曲線１１０７１、観測データ１１０７２、および領域１１０７３に本統計確率分布に関連する項目が表示される。領域１１０７３に表示される表示項目内容には、センサ番号、計測内容、本統計確率分布の経過時刻、平均値の数値、標準偏差の数値、観測データの推定値の該統計確率分布における確率値と異常度がある。領域１１０７３の表示はマウスの右クリックにて表示されないように選択しても良い。
【０１４２】
図１１Ａ又は図１１Ｂで説明したようなＧＵＩにより、センサデータ推定手法の選択やパラメータの設定等を簡便に入力することが可能である。また、センサデータ推定前後の結果を確認することが可能であるため、選択した手法や設定したパラメータの妥当性を確認することが可能である。さらに、確率分布および新たに観測されたデータの所在位置や異常度を観察し、シーケンスの進行状況を確認することが可能である。
【０１４３】
以上、本発明者によってなされた発明を実施例に基づき具体的に説明したが、本発明は前記実施例に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能であることは言うまでもない。
【符号の説明】
【０１４４】
１０１・・・設備    １００・・・異常検知パス    １００’・・・学習パス    １００”・・・推定時刻決定パス    １００１・・・イベントデータ    １００２・・・センサデータ    １００３・・・ユーザ指示    １００４・・・識別境界or正常空間    １０２・・・データ前処理部    １０２１・・・イベントデータ解析部    １０２２・・・センサデータ切出し部    １０２３・・・センサデータ時刻合わせ部    １０３・・・センサデータ推定部    １０４・・・統計確率分布推定部    １０５・・・特徴量抽出部    １０６・・・異常検知部    １１１・・・データベース    １１２・・・推定時刻決定部    １１３・・・学習部。
【図１Ａ】
2014048697.tif 000016
【図１Ｂ】
2014048697.tif 000017
【図１Ｃ】
2014048697.tif 000018
【図１Ｄ】
2014048697.tif 000019
【図１Ｅ】
2014048697.tif 000020
【図２Ａ】
2014048697.tif 000021
【図２Ｂ】
2014048697.tif 000022
【図２Ｃ】
2014048697.tif 000023
【図３Ａ】
2014048697.tif 000024
【図３Ｂ】
2014048697.tif 000025
【図３Ｃ】
2014048697.tif 000026
【図４】
2014048697.tif 000027
【図５Ａ】
2014048697.tif 000028
【図５Ｂ】
2014048697.tif 000029
【図５Ｃ】
2014048697.tif 000030
【図６Ａ】
2014048697.tif 000031
【図６Ｂ】
2014048697.tif 000032
【図７Ａ】
2014048697.tif 000033
【図７Ｂ】
2014048697.tif 000034
【図８Ａ】
2014048697.tif 000035
【図８Ｂ】
2014048697.tif 000036
【図８Ｃ】
2014048697.tif 000037
【図９Ａ】
2014048697.tif 000038
【図９Ｂ】
2014048697.tif 000039
【図１０】
2014048697.tif 000040
【図１１Ａ】
2014048697.tif 000041
【図１１Ｂ】
2014048697.tif 000042
【図１２】
2014048697.tif 000043