JP2020035081A - 異常検出装置、異常検出方法、異常検出プログラム及び異常検出システム - Google Patents

Abstract

【課題】正常時の時間波形及び周波数波形との類似度から高い精度で異常を検出する。
【解決手段】正常時周期データを、前記周期の時間幅を示すワーク毎に分割して複数の正常時ワークデータを作成する分割手段と、前記分割手段により分割された複数の正常時ワークデータのそれぞれに対して高速フーリエ変換を行って前記複数の正常時ワークデータのそれぞれを周波数領域に変換した複数の変換後正常時ワークデータを作成する変換手段と、時間領域における波形の正常モデルを正常時ワークデータから作成すると共に、周波数領域における波形の正常モデルを複数の変換後正常時ワークデータから作成するモデル作成手段と、ふたつのと、前記周期データとに基づいて、所定の指標値を算出する指標値算出手段と、前記指標値算出手段により算出された指標値と、予所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する判定手段とを有する異常検出装置。
【選択図】図２

Description

本発明は、異常検出装置、異常検出方法、異常検出プログラム及び異常検出システムに関する。

製造装置の装置パラメータの時系列データを特徴量化して、分類された不良パターンの有無による特徴量の頻度分布の有意差を検定することで、不良発生原因の装置パラメータの特定が可能な技術が従来から知られている（例えば特許文献１）。また、異常を検出した装置検査情報から、問題となる装置の特定を可能にする技術が従来から知られている（例えば特許文献２）。

特開２００５−２５１９２５号公報 特開２００４−１５３２２８号公報

しかしながら、上記の従来技術では、周期的に類似する波形が現れる時系列データ（例えば、バッチプロセスのプロセスデータ等）における波形の類似性から異常を検出することができない。これに対して、バッチプロセスにおけるバッチ毎の時間波形（すなわち、時間領域における波形）の類似性に着目して、正常時の時間波形との差を監視することで異常を検出する手法も知られているが、正常時の時間波形との部分的な差に基づく異常検出であり、検出精度が十分でない場合があった。

また、時間波形と共に、周波数波形（すなわち、周波数領域における波形）を用いた異常検出を行うことで、より高い精度での異常検出が期待できる。

本発明の実施の形態は、上記の点に鑑みてなされたもので、正常時の時間波形及び周波数波形との類似度から高い精度で異常を検出することを目的とする。

上記目的を達成するため、本発明の実施の形態は、機器の周期的な処理又は動作を表す周期データから前記機器の異常を検出する異常検出装置であって、正常時における前記機器の周期的な処理又は動作を示す正常時周期データを、前記周期の時間幅を示すワーク毎に分割して、複数の正常時ワークデータを作成する分割手段と、前記分割手段により分割された複数の正常時ワークデータのそれぞれに対して、窓関数を用いた高速フーリエ変換を行って、前記複数の正常時ワークデータのそれぞれを周波数領域に変換した複数の変換後正常時ワークデータを作成する変換手段と、時間領域における波形の正常モデルを表す第１のプロファイルを前記複数の正常時ワークデータから作成すると共に、周波数領域における波形の正常モデルを表す第２のプロファイルを前記複数の変換後正常時ワークデータから作成するモデル作成手段と、前記第１のプロファイル及び前記第２のプロファイルと、前記周期データとに基づいて、所定の指標値を算出する指標値算出手段と、前記指標値算出手段により算出された指標値と、予め設定された所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する判定手段と、を有することを特徴とする。

正常時の時間波形及び周波数波形との類似度から高い精度で異常を検出することができる。

本実施形態に係る異常検出システムの全体構成の一例を示す図である。 本実施形態に係る異常検出装置のハードウェア構成の一例を示す図である。 本実施形態に係る異常検出装置の機能構成の一例を示す図である。 本実施形態に係る正常モデル作成処理の一例を示すフローチャートである。 本実施形態に係る異常検出処理の一例を示すフローチャートである。 異常検出ワークの判定の一例を説明するための図である。 異常検出結果画面の一例を示す図である。 本実施形態に係る異常検出システムの全体構成の他の例を示す図である。

以下、本発明の実施の形態について、図面を参照しながら詳細に説明する。以降では、正常時の時間波形及び周波数波形との類似度から高い精度で異常を検出することが可能な異常検出システム１について説明する。ここで、時間波形とは、周期的に類似する波形が現れる時系列データの時間領域における波形のことである。また、周波数波形とは、当該時系列データの周波数領域における波形のことである。

なお、周期的に類似する波形が現れる時系列データとしては、例えば、バッチプロセスのプロセスデータ（すなわち、バッチプロセスの或る変数をセンシングすることで得られる時系列データ）が挙げられる。これ以外にも、例えば、周期的に何等かの処理又は動作を実行する工作機械や製造装置、産業用ロボット等の機器の振動をセンシングすることで得られる振動データ等が挙げられる。以降では、周期的に類似する波形が現れる時系列データを「周期データ」と表す。

＜全体構成＞
まず、本実施形態に係る異常検出システム１の全体構成について、図１を参照しながら説明する。図１は、本実施形態に係る異常検出システム１の全体構成の一例を示す図である。

図１に示すように、本実施形態に係る異常検出システム１は、異常検出装置１０と、センシング機器２０とが含まれる。異常検出装置１０とセンシング機器２０とは、例えばＬＡＮ（Local Area Network）等のネットワークを介して通信可能に接続されている。

センシング機器２０は、異常の発生を検出する対象である異常検出対象３０を計測する計測機器である。ここで、異常検出対象３０がバッチプロセスを実行する設備又は機器である場合は、センシング機器２０は、例えば、当該バッチプロセスの状態（例えば、バッチプロセスの温度や圧力、流量、濃度等）を計測する。また、異常検出対象３０が工作機械や製造装置、産業用ロボット等の機器である場合は、センシング機器２０は、例えば、当該機器の振動等を計測する。

このように、センシング機器２０は、例えばバッチプロセス等の何等かの周期的な処理（又は動作）を実行する異常検出対象３０を所定の時間毎（すなわち、サンプリング周期毎）に計測する。そして、センシング機器２０は、この計測結果を示す時系列データ（すなわち、周期データ）を作成した上で、作成した周期データを異常検出装置１０に送信する。

異常検出対象３０は、上述したように、例えば、バッチプロセスを実行する設備又は機器である。又は、異常検出対象３０は、上述したように、例えば、周期的に何等かの処理又は動作を実行する工作機械や製造装置、産業用ロボット等の機器である。なお、異常検出対象３０は、これら以外にも、周期的に何等かの処理又は動作を実行する機器であれば良く、例えば、コンベアやローラ等の産業機械、鉄道車両等であっても良い。

以降では、周期データにおける１周期分の時間区間を「ワーク」と表す。例えば、バッチプロセスを異常検出対象３０が実行する場合、１ワークは、１つのバッチ処理が開始してから終了するまでの時間区間である。同様に、例えば、周期的な何等かの処理又は動作を異常検出対象３０が実行する場合、１ワークは、１つの処理又は動作が開始してから終了するまでの時間区間である。

異常検出装置１０は、センシング機器２０から受信した周期データに基づいて、異常検出対象３０に異常が発生したことを検出するコンピュータである。なお、異常検出装置１０としては、例えば、ＰＬＣ（Programmable Logic Controller）等の制御装置が用いられても良い。

本システムの動作には、異常検出対象３０で異常が発生したことを検出するための正常モデルを作成する「モデル作成」フェーズと、異常検出対象３０の処理又は動作中に計測された周期データと正常モデルとから異常を検出する「評価」フェーズとがある。基本的に、「モデル作成」フェーズは異常検出対象３０が処理又は動作していない時に実行されるオフラインの処理であり、「評価」フェーズは異常検出対象３０の処理又は動作中に実行されるオンラインの処理である。ただし、これに限られず、「モデル作成」フェーズ及び「評価」フェーズの両方がオフラインの処理であっても良いし、「モデル作成」フェーズ及び「評価」フェーズの両方がオンラインの処理であっても良い。

異常検出装置１０は、「モデル作成」フェーズにおいて、モデル作成用の周期データから正常モデルを作成する。モデル作成用の周期データとは、例えば、正常な処理又は動作を実行している異常検出対象３０をセンシング機器２０で計測することで作成された周期データのことである。なお、モデル作成用の周期データは、異常検出対象３０の正常な処理又は動作を示すものとして、ユーザ等によって作成された周期データであっても良い。

また、異常検出装置１０は、「評価」フェーズにおいて、評価用の周期データと、正常モデルとに基づいて、異常検出対象３０で発生した異常を検出する。評価用の周期データとは、例えば、異常検出対象３０のオンライン中に、当該異常検出対象３０をセンシング機器２０が計測することによって作成された周期データのことである。

なお、本実施形態に係る異常検出システム１には、複数の種類の異常検出対象３０が含まれていても良い。この場合、本実施形態に係る異常検出装置１０は、異常検出対象３０の種類毎に正常モデルを作成すると共に、異常検出対象３０の種類毎に当該異常検出対象３０で発生した異常を検出すれば良い。

また、異常検出対象３０が複数の種類の処理又は動作を行っても良い。例えば、異常検出対象３０が複数の工程（工程Ａ〜工程Ｃ）によって製品を製造する機器である場合、当該異常検出対象３０は、工程Ａにおける動作Ａと、工程Ｂにおける動作Ｂと、工程Ｃにおける動作Ｃとを行っても良い。この場合、本実施形態に係る異常検出装置１０は、動作毎に正常モデルを作成すると共に、動作毎に異常検出対象３０で発生した異常を検出すれば良い。

また、複数のセンシング機器２０が１つの異常検出対象３０を計測しても良い。この場合、異常検出装置１０は、複数のセンシング機器２０がそれぞれ作成した複数の周期データを１つにまとめたデータにより、正常モデルの作成や異常検出を行えば良い。１つにまとめたデータとは、例えば、２台のセンシング機器２０が１台の異常検出対象３０の温度及び圧力をそれぞれ計測する場合、温度を表す時系列データと圧力を表す時系列データとが含まれる周期データのことである。

＜ハードウェア構成＞
次に、本実施形態に係る異常検出装置１０のハードウェア構成について、図２を参照しながら説明する。図２は、本実施形態に係る異常検出装置１０のハードウェア構成の一例を示す図である。

図２に示すように、本実施形態に係る異常検出装置１０は、入力装置１１と、表示装置１２と、外部Ｉ／Ｆ１３と、ＲＡＭ（Random Access Memory）１４と、ＲＯＭ（Read Only Memory）１５と、ＣＰＵ（Central Processing Unit）１６と、通信Ｉ／Ｆ１７と、補助記憶装置１８とを有する。これらの各ハードウェアは、バス１９により通信可能に接続されている。

入力装置１１は、例えばキーボードやマウス、タッチパネル等であり、ユーザが各種操作を入力するのに用いられる。表示装置１２は、例えばＬＣＤ（Liquid Crystal Display）等であり、異常検出装置１０の処理結果を表示する。なお、異常検出装置１０は、入力装置１１及び表示装置１２のうちの少なくとも一方を有していなくても良い。

外部Ｉ／Ｆ１３は、外部装置とのインタフェースである。外部装置には、記録媒体１３ａ等がある。異常検出装置１０は、外部Ｉ／Ｆ１３を介して、記録媒体１３ａ等の読み取りや書き込みを行うことができる。記録媒体１３ａには、例えば、フレキシブルディスク、ＣＤ（Compact Disc）、ＤＶＤ（Digital Versatile Disc）、ＳＤメモリカード、ＵＳＢメモリ等がある。なお、記録媒体１３ａには、本実施形態に係る異常検出装置１０の各機能を実現するプログラムが格納されていても良い。

ＲＡＭ１４は、プログラムやデータを一時保持する揮発性の半導体メモリである。ＲＯＭ１５は、電源を切ってもプログラムやデータを保持することができる不揮発性の半導体メモリである。ＲＯＭ１５には、例えば、異常検出装置１０の起動時に実行されるＢＩＯＳ（Basic Input/Output System）、ＯＳ（Operating System）設定やネットワーク設定等が格納されている。

ＣＰＵ１６は、ＲＯＭ１５や補助記憶装置１８等からプログラムやデータをＲＡＭ１４上に読み出し、処理を実行することで、異常検出装置１０全体の制御や機能を実現する演算装置である。

通信Ｉ／Ｆ１７は、異常検出装置１０が他の機器等と通信を行うためのインタフェースである。異常検出装置１０は、通信Ｉ／Ｆ１７を介して、センシング機器２０から周期データを受信することができる。

補助記憶装置１８は、プログラムやデータを格納している不揮発性のメモリであり、例えば、ＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）等である。補助記憶装置１８に格納されるプログラムやデータには、本実施形態に係る異常検出装置１０の各機能を実現するプログラム、異常検出装置１０全体を制御する基本ソフトウェアであるＯＳ、及びＯＳ上において各種機能を提供するアプリケーションソフトウェア等がある。なお、補助記憶装置１８は、格納しているプログラムやデータを所定のファイルシステムやＤＢ（データベース）等により管理している。

本実施形態に係る異常検出装置１０は、上記のハードウェア構成を有することにより、後述するような各種処理を実現できる。

＜機能構成＞
次に、本実施形態に係る異常検出装置１０の機能構成について、図３を参照しながら説明する。図３は、本実施形態に係る異常検出装置１０の機能構成の一例を示す図である。

図３に示すように、本実施形態に係る異常検出装置１０は、データ取得部１０１と、分割部１０２と、周波数変換部１０３と、モデル作成部１０４と、指標値算出部１０５と、異常判定部１０６と、出力部１０７とを有する。これら各機能部は、異常検出装置１０にインストールされた１以上のプログラムがＣＰＵ１６に実行させる処理により実現される。

また、本実施形態に係る異常検出装置１０は、周期データ記憶部１１０と、モデル記憶部１２０とを有する。これら各記憶部は、例えば補助記憶装置１８を用いて実現可能である。なお、これら各記憶部のうちの少なくとも１つの記憶部が、異常検出装置１０とネットワークを介して接続される記憶装置等を用いて実現されていても良い。

周期データ記憶部１１０は、モデル作成用の周期データと、評価用の周期データとを記憶する。これらのモデル作成用の周期データ及び評価用の周期データは、サンプリング周期毎に、センシング機器２０により計測された計測値が含まれる時系列データである。言い換えれば、周期データは、例えば、横軸を時間、縦軸を計測値とした時間領域で表される。なお、計測値は、センシング機器２０のセンシング対象がバッチプロセスの或る変数である場合、当該変数の値（例えば、温度や圧力、流量等）である。又は、計測値は、センシング機器２０のセンシング対象が振動である場合、加速度等である。

なお、モデル作成用の周期データと、評価用の周期データとを異なるデータとして区別して周期データ記憶部１１０に記憶されている必要はない。例えば、１つの周期データにおいて、或る時間幅のデータ（例えば、時刻ｔ＝ｔ_１〜時刻ｔ＝ｔ_２までの間に計測された計測値が含まれる周期データ）をモデル作成用の周期データとし、別の或る時間幅（例えば、時刻ｔ＝ｔ_３〜時刻ｔ＝ｔ_４まで間に計測された計測値が含まれる周期データ）を評価用の周期データとしても良い。

データ取得部１０１は、「モデル作成」フェーズにおいて、モデル作成用の周期データを周期データ記憶部１１０から取得する。また、データ取得部１０１は、「評価」フェーズにおいて、評価用の周期データを周期データ記憶部１１０から取得する。

分割部１０２は、「モデル作成」フェーズにおいて、データ取得部１０１により取得された周期データ（モデル作成用の周期データ）をワーク毎に分割する。以降では、ワーク単位のモデル作成用の周期データを「モデル作成用ワークデータ」と表す。モデル作成用の周期データにＮ個のワークが含まれる場合、このモデル作成用の周期データは、Ｎ個のモデル作成用ワークデータに分割される。

また、分割部１０２は、「評価」フェーズにおいて、データ取得部１０１により取得された周期データ（評価用の周期データ）をワーク毎に分割する。以降では、ワーク単位の評価用の周期データを「評価用ワークデータ」と表す。評価用の周期データにＭ個のワークが含まれる場合、この評価用の周期データは、Ｍ個の評価用ワークデータに分割される。

周波数変換部１０３は、「モデル作成」フェーズにおいて、モデル作成用ワークデータ毎に、当該モデル作成用ワークデータに対して窓関数を用いた高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）を行って、窓（ウインドウ）毎に周波数領域に変換する。これにより、各モデル作成用ワークデータにおいて、窓毎に、縦軸をスペクトル強度、横軸を周波数としたパワースペクトルが得られる。例えば、各モデル作成用ワークデータ内にＬ個の窓（ウインドウ）がそれぞれ含まれる場合、各モデル作成用ワークデータにおいて、Ｌ個のパワースペクトルがそれぞれ得られる。以降では、周波数変換部１０３により変換されたモデル作成用ワークデータを「変換後モデル作成用ワークデータ」と表す。

同様に、周波数変換部１０３は、「評価」フェーズにおいて、評価用ワークデータ毎に、当該評価用ワークデータに対して窓関数を用いた高速フーリエ変換を行って、窓（ウインドウ）毎に周波数領域に変換する。これにより、各評価用ワークデータにおいて、窓毎に、縦軸をスペクトル強度、横軸を周波数としたパワースペクトルが得られる。例えば、各評価用ワークデータ内にＬ個の窓（ウインドウ）がそれぞれ含まれる場合、各評価用ワークデータにおいて、Ｌ個のパワースペクトルがそれぞれ得られる。以降では、周波数変換部１０３により変換された評価用ワークデータを「変換後評価用ワークデータ」と表す。

モデル作成部１０４は、「モデル作成」フェーズにおいて、各モデル作成用ワークデータの平均を算出することで、時間領域における波形（つまり、時間波形）の正常モデルである時間波形プロファイルを作成する。そして、モデル作成部１０４は、時間波形プロファイルをモデル記憶部１２０に保存する。なお、各モデル作成用ワークデータの平均とは、各ワーク間で対応する時刻における計測値の平均のことである。また、各ワーク間で対応する時刻とは、各ワークの開始時刻を基準として、各ワークで相対時刻が同一である各時刻のことである。

また、モデル作成部１０４は、「モデル作成」フェーズにおいて、各変換後モデル作成用ワークデータの平均を算出することで、周波数領域における波形（つまり、周波数波形）の正常モデルである周波数波形プロファイルを作成する。そして、モデル作成部１０４は、周波数波形プロファイルをモデル記憶部１２０に保存する。なお、各変換後モデル作成用ワークデータの平均とは、各ワーク間で同一の周波数におけるパワースペクトル値（又は「スペクトル強度」とも称される。）の平均のことである。

指標値算出部１０５は、「評価」フェーズにおいて、モデル記憶部１２０に記憶されている時間波形プロファイルと、各評価用ワークデータのそれぞれとの類似度を表す第１の指標値を算出する。また、指標値算出部１０５は、「評価」フェーズにおいて、モデル記憶部１２０に記憶されている周波数波形プロファイルと、各変換後評価用ワークデータのそれぞれとの類似度を表す第２の指標値を算出する。なお、評価用ワークデータ数がＭ個である場合、Ｍ個の第１の指標値とＭ個の第２の指標値とが算出される。

ここで、第１の指標値としては、例えば、時間波形プロファイルと評価用ワークデータとの相関係数が挙げられる。同様に、第２の指標値としては、例えば、周波数波形プロファイルと変換後評価用ワークデータとの相関係数が挙げられる。以降では、第１の指標値及び第２の指標値は相関係数であるものとして説明する。

ただし、第１の指標値及び第２の指標値は相関係数に限られない。第１の指標値としては、例えば、時間波形プロファイルと評価用ワークデータとの平均値や中央値、マハラノビス距離等を用いることができる。同様に、第２の指標値としては、例えば、周波数波形プロファイルと変換後評価用ワークデータとの平均値や中央値、マハラノビス距離等を用いることができる。

このように、第１の指標値としては、時間波形プロファイルと評価用ワークデータとの類似度を表す任意の指標値を用いることができる。同様に、第２の指標値としては、周波数波形プロファイルと評価用ワークデータとの類似度を表す任意の指標値を用いることができる。

異常判定部１０６は、「評価」フェーズにおいて、各第１の指標値と、予め設定された第１の閾値との大小関係を比較判定する。また、異常判定部１０６は、各第２の指標値と、予め設定された第２の閾値との大小関係を比較判定する。

そして、異常判定部１０６は、例えば、第１の指標値が第１の閾値以下である評価用ワークデータに対応するワークを異常検出ワークと判定する。同様に、異常判定部１０６は、第２の指標値が第２の閾値以下である評価用ワークデータに対応するワークを異常検出ワークと判定する。

出力部１０７は、「評価」フェーズにおいて、例えば、異常判定部１０６により異常検出ワークと判定された評価用ワークデータに関する情報が含まれる画面（後述する異常検出結果画面）を表示装置１２に出力する。

ここで、評価用ワークデータに関する情報としては、例えば、当該評価用ワークデータを識別するワークＩＤ、ワークの年月日、ワークを代表する時刻（例えば、ワークの開始時刻、終了時刻、又は開始時刻と終了時刻との中間を示す時刻等）、及び指標値（第１の指標値又は第２の指標値）等が挙げられる。

＜正常モデル作成処理＞
次に、「モデル作成」フェーズにおいて、正常モデルを作成する処理（正常モデル作成処理）について、図４を参照しながら説明する。図４は、本実施形態に係る正常モデル作成処理の一例を示すフローチャートである。

まず、データ取得部１０１は、モデル作成用の周期データを周期データ記憶部１１０から取得する（ステップＳ１０１）。

次に、分割部１０２は、上記のステップＳ１０１で取得したモデル作成用の周期データをワーク毎に分割して、モデル作成用ワークデータを作成する（ステップＳ１０２）。以降では、モデル作成用の周期データがＮ個のモデル作成用ワークデータに分割されたものとして説明する。

以降のステップＳ１０３の処理と、ステップＳ１０４〜ステップＳ１０５の処理とは並列に実行される。ただし、ステップＳ１０３の処理が実行された後、ステップＳ１０４〜ステップＳ１０５の処理が実行されても良いし、ステップＳ１０４〜ステップＳ１０５の処理が実行された後、ステップＳ１０３の処理が実行されても良い。又は、例えば、ステップＳ１０４の処理、ステップＳ１０３の処理、ステップＳ１０５の処理の順で実行されても良い。

モデル作成部１０４は、Ｎ個のモデル作成用ワークデータの平均を算出して、時間波形プロファイルを作成する。すなわち、例えば、１ワーク内の時間幅をＴ、各ワークの開始時刻を基準としたワーク内の相対時刻をｔ（１≦ｔ≦Ｔ）、ワークｎ（１≦ｎ≦Ｎ）のモデル作成用ワークデータの相対時刻ｔにおける計測値をｕ_ｎ（ｔ）とした場合、モデル作成部１０４は、相対時刻ｔ毎に、Ｐ（ｔ）＝（ｕ_１（ｔ）＋・・・＋ｕ_Ｎ（ｔ））／Ｎを算出する。このＰ（ｔ），１≦ｔ≦Ｔが時間波形プロファイルである。

そして、モデル作成部１０４は、作成した時間波形プロファイルをモデル記憶部１２０に保存する（ステップＳ１０３）。

周波数変換部１０３は、モデル作成用ワークデータ毎に、窓関数を用いた高速フーリエ変換を行って、窓（ウインドウ）毎に周波数領域に変換する（ステップＳ１０４）。

より具体的には、周波数変換部１０３は、例えば、２０４８個のデータ値が含まれる時間幅をウインドウ幅、オーバラップ率を５０％として、各モデル作成用ワークデータに対して、ウインドウ毎に高速フーリエ変換を行ってパワースペクトルを算出する。これにより、例えば、１つのモデル作成用ワークデータに対して、Ｋ＝６５５３６／（２０４８／２）＝６４個のパワースペクトルが算出される。すなわち、Ｎ個のモデル作成用ワークデータのそれぞれに対して、Ｋ個のパワースペクトルが算出される。このＫ個のパワースペクトルが変換後モデル作成用ワークデータである。

次に、モデル作成部１０４は、Ｎ個の変換後モデル作成用ワークデータの平均を算出して、周波数波形プロファイルを作成する。すなわち、各変換後モデル作成用ワークデータに含まれる周波数を表す番号（周波数番号）をｋ（１≦ｋ≦Ｋ）、ワークｎ（１≦ｎ≦Ｎ）の変換後モデル作成用ワークデータの周波数番号ｋにおけるパワースペクトル値をｖ_ｎ（ｋ）とした場合、モデル作成部１０４は、周波数番号ｋ毎に、Ｑ（ｋ）＝（ｖ_１（ｋ）＋・・・＋ｖ_Ｎ（ｋ））／Ｎを算出する。このＱ（ｋ），１≦ｋ≦Ｋが周波数波形プロファイルである。

そして、モデル作成部１０４は、作成した周波数波形プロファイルをモデル記憶部１２０に保存する（ステップＳ１０５）。

以上のように、本実施形態に係る異常検出装置１０は、時間領域における波形の正常モデルを表す時間波形プロファイルと、周波数領域における波形の正常モデルを表す周波数波形プロファイルとを作成する。後述するように、これらの正常モデルとの類似度を評価することで、異常検出対象３０で異常が発生したか否かが判定される。

＜異常検出処理＞
次に、「評価」フェーズにおいて、正常モデルを用いて、異常検出対象３０で異常が発生したことを検出する処理（異常検出処理）について、図５を参照しながら説明する。図５は、本実施形態に係る異常検出処理の一例を示すフローチャートである。

まず、データ取得部１０１は、評価用の周期データを周期データ記憶部１１０から取得する（ステップＳ２０１）。

次に、分割部１０２は、上記のステップＳ２０１で取得した評価用の周期データをワーク毎に分割して、評価用ワークデータを作成する（ステップＳ２０２）。以降では、評価用の周期データがＭ個の評価用ワークデータに分割されたものとして説明する。

以降のステップＳ２０３の処理と、ステップＳ２０４〜ステップＳ２０５の処理とは並列に実行される。ただし、ステップＳ２０３の処理が実行された後、ステップＳ２０４〜ステップＳ２０５の処理が実行されても良いし、ステップＳ２０４〜ステップＳ２０５の処理が実行された後、ステップＳ２０３の処理が実行されても良い。又は、例えば、ステップＳ２０４の処理、ステップＳ２０３の処理、ステップＳ２０５の処理の順で実行されても良い。

指標値算出部１０５は、モデル記憶部１２０に記憶されている時間波形プロファイルと、Ｍ個の評価用ワークデータのそれぞれとの類似度を表す相関係数（第１の指標値）を算出する（ステップＳ２０３）。すなわち、Ｍ個の評価用ワークデータのそれぞれをｐ_ｍ（ｔ），１≦ｔ≦Ｔとした場合、指標値算出部１０５は、各ｍ（１≦ｍ≦Ｍ）に対して、時間波形プロファイルＰ（ｔ）と、ｍ番目の評価用ワークデータｐ_ｍ（ｔ）との相関係数ｒ_ｍを算出する。なお、ｔは、各ワークの開始時刻を基準としたワーク内の相対時刻である。

これにより、各ｍ（１≦ｍ≦Ｍ）に対して、時間波形プロファイルとｍ番目の評価用ワークデータとの類似度を表す第１の指標値として相関係数ｒ_ｍが得られる。

周波数変換部１０３は、評価用ワークデータ毎に、窓関数を用いた高速フーリエ変換を行って、窓（ウインドウ）毎に周波数領域に変換する（ステップＳ２０４）。

より具体的には、周波数変換部１０３は、上記のステップＳ１０４と同様に、例えば、２０４８個のデータ値が含まれる時間幅をウインドウ幅、オーバラップ率を５０％として、各評価用ワークデータに対して、ウインドウ毎に高速フーリエ変換を行ってパワースペクトルを算出する。これにより、例えば、１つの評価用ワークデータに対して、Ｋ＝６５５３６／（２０４８／２）＝６４個のパワースペクトルが算出される。すなわち、Ｍ個の評価用ワークデータのそれぞれに対して、Ｋ個のパワースペクトルが算出される。このＫ個のパワースペクトルが変換後評価用ワークデータである。

次に、指標値算出部１０５は、モデル記憶部１２０に記憶されている周波数波形プロファイルと、Ｍ個の評価用ワークデータのそれぞれとの類似度を表す相関係数（第２の指標値）を算出する（ステップＳ２０５）。すなわち、Ｍ個の変換後評価用ワークデータのそれぞれをｑ_ｍ（ｋ），１≦ｋ≦Ｋとした場合、指標値算出部１０５は、各ｍ（１≦ｍ≦Ｍ）に対して、周波数波形プロファイルＱ（ｋ）と、ｍ番目の変換後評価用ワークデータｑ_ｍ（ｋ）との相関係数ｓ_ｍを算出する。なお、ｋは、周波数番号である。

これにより、各ｍ（１≦ｍ≦Ｍ）に対して、周波数波形プロファイルとｍ番目の評価用ワークデータとの類似度を表す第２の指標値として相関係数ｓ_ｍが得られる。

異常判定部１０６は、各相関係数ｒ_ｍ（１≦ｍ≦Ｍ）と第１の閾値ｔｈ_ｔとの大小関係を比較判定すると共に、各相関係数ｓ_ｍ（１≦ｍ≦Ｍ）と第２の閾値ｔｈ_ｆとの大小関係を比較判定する（ステップＳ２０６）。

そして、異常判定部１０６は、各相関係数ｒ_ｍ（１≦ｍ≦Ｍ）のうち、第１の閾値ｔｈ_ｔ以下と判定された相関係数ｒ_ｍに対応するワークを異常検出ワークと判定する。同様に、異常判定部１０６は、各相関係数ｓ_ｍ（１≦ｍ≦Ｍ）のうち、第２の閾値ｔｈ_ｆ以下と判定された相関係数ｓ_ｍに対応するワークを異常検出ワークと判定する。

ここで、異常検出ワークと判定される場合の一例を図６（ａ）及び図６（ｂ）に示す。図６（ａ）は、第１の閾値ｔｈ_ｔ＝０．７として、相関係数ｒ_ｍがｔｈ_ｔ以下であるワークを異常検出ワークと判定する場合を示している。図６（ａ）に示す例では、ワークＩＤ「５５」、「８０」、「８１」のワークが、時間領域における異常検出ワークと判定されている。また、図６（ｂ）は、第２の閾値ｔｈ_ｆ＝０．５として、相関係数ｓ_ｍがｔｈ_ｆ以下であるワークを異常検出ワークと判定する場合を示している。図６（ｂ）に示す例では、ワークＩＤ「８０」のワークが、周波数領域における異常検出ワークと判定されている。

このように、異常判定部１０６は、第１の指標値と第１の閾値との比較結果に応じて異常検出ワークを判定すると共に、第２の指標値と第２の閾値との比較結果に応じて異常検出ワークを判定する。

次に、出力部１０７は、上記のステップＳ２０６で異常検出ワークと判定された評価用ワークデータに関する情報が含まれる異常検出結果画面を表示装置１２に出力する（ステップＳ２０７）。ここで、異常検出結果画面の一例を図７に示す。図７に示す異常検出結果画面１０００は、上記のステップＳ２０６でワークＩＤ「５５」、「８０」、「８１」のワークが時間領域における異常検出ワークと判定されたと共に、ワークＩＤ「８０」のワークが周波数領域における異常検出ワークと判定された場合を示している。

図７に示す異常検出結果画面１０００には、時間領域における異常検出ワークの評価用ワークデータに関する情報が表示される第１の表示欄１１００と、周波数領域における異常検出ワークの評価用ワークデータに関する情報が表示される第２の表示欄１２００とが含まれる。これにより、ユーザは、異常検出対象３０で異常が検出されたワークの情報を知ることができる。

このとき、第１の表示欄１１００及び第２の表示欄１２００では、時間領域における異常検出ワークと判定されると共に周波数領域における異常検出ワークと判定されたワークの評価用ワークデータに関する情報を所定の表示態様で表示する。図７に示す例では、ワークＩＤ「８０」の評価用ワークデータに関する情報が網掛けされ、強調表示されている。これにより、ユーザは、時間領域と周波数領域との両方で異常と判定されたワークの情報を知ることができる。

なお、異常検出結果画面を出力する場合にも、出力部１０７は、上記のステップＳ２０６における判定結果を任意の出力先に出力しても良い。このような出力先としては、例えば、補助記憶装置１８、ネットワークを介して接続される所定のサーバ装置等が挙げられる。

以上のように、本実施形態に係る異常検出装置１０は、ワーク単位の評価用ワークデータと時間波形プロファイルとの類似度を示す第１の指標値を算出すると共に、ワーク単位の評価用ワークデータと周波数波形プロファイルとの類似度を示す第２の指標値を算出する。そして、本実施形態に係る異常検出装置１０は、これらの第１の指標値及び第２の指標値を用いて、時間領域と周波数領域との両方で異常の有無を判定する。

このように、本実施形態に係る異常検出装置１０では、時間領域と周波数領域との両方で正常モデルとの類似度を判定する。これにより、異常検出対象３０で発生した異常を高い精度で検出することができるようになる。

なお、本実施形態では、上記のステップＳ２０６で相関係数ｒ_ｍ（１≦ｍ≦Ｍ）と第１の閾値ｔｈ_ｔとを比較判定すると共に、相関係数ｓ_ｍ（１≦ｍ≦Ｍ）と第２の閾値ｔｈ_ｆとを比較判定したが、これに限られない。例えば、予め設定された重みをα、β（ただし、α＋β＝１，α≧０，β≧０）、Ｒ_ｍ＝α×ｒ_ｍ＋β×ｓ_ｍとして、各ｍ（１≦ｍ≦Ｍ）に対して、このＲ_ｍと、予め設定された閾値ｔｈとを比較判定しても良い。この場合、異常判定部１０６は、各Ｒ_ｍ（１≦ｍ≦Ｍ）のうち、閾値ｔｈ以下と判定されたＲ_ｍに対応するワークを異常検出ワークと判定すれば良い。

＜異常検出システム１の他の例＞
ここで、本実施形態に係る異常検出システム１の全体構成の他の例について、図８を参照しながら説明する。図８は、本実施形態に係る異常検出システム１の全体構成の他の例を示す図である。

図８に示すように、本実施形態に係る異常検出システム１は、異常検出装置１０と、センシング機器２０と、表示端末４０とを有し、例えばインターネット等のネットワークＮを介して通信可能に接続される構成であっても良い。言い換えれば、表示端末４０のユーザに対して、異常検出装置１０による正常モデル作成処理及び異常検出処理が、クラウド型のサービスとして提供されても良い。

図８に示す異常検出システム１では、異常検出装置１０は、出力部１０７によって異常検出結果画面を表示端末４０に送信する。これにより、表示端末４０には、例えば図７に示す異常検出結果画面１０００が表示される。なお、表示端末４０としては、例えば、ＰＣ（パーソナルコンピュータ）やスマートフォン、タブレット端末等を用いることができる。

＜まとめ＞
以上のように、本実施形態に係る異常検出システム１は、例えばオフラインにおいて、異常検出対象３０の正常な処理又は動作を示す周期データを用いて、時間領域における波形の正常モデルを表す時間波形プロファイルと、周波数領域における波形の正常モデルを表す周波数波形プロファイルとを作成する。

そして、本実施形態に係る異常検出システム１は、オンラインで動作している異常検出対象３０の処理又は動作をセンシング機器２０で計測することで取得された周期データと、時間波形プロファイルと、周波数波形プロファイルとを用いて、当該異常検出対象３０の異常を検出する。このように、本実施形態に係る異常検出システム１では、時間領域と周波数領域との両方で正常モデルとの類似度を判定することで、異常検出対象３０の異常を高い精度で検出することができるようになる。

本発明は、具体的に開示された上記の実施形態に限定されるものではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。

１ 異常検出システム
１０ 異常検出装置
２０ センシング機器
３０ 異常検出対象
１０１ データ取得部
１０２ 分割部
１０３ 周波数変換部
１０４ モデル作成部
１０５ 指標値算出部
１０６ 異常判定部
１０７ 出力部
１１０ 周期データ記憶部
１２０ モデル記憶部

Claims (8)

1. 機器の周期的な処理又は動作を表す周期データから前記機器の異常を検出する異常検出装置であって、
   正常時における前記機器の周期的な処理又は動作を示す正常時周期データを、前記周期の時間幅を示すワーク毎に分割して、複数の正常時ワークデータを作成する分割手段と、
   前記分割手段により分割された複数の正常時ワークデータのそれぞれに対して、窓関数を用いた高速フーリエ変換を行って、前記複数の正常時ワークデータのそれぞれを周波数領域に変換した複数の変換後正常時ワークデータを作成する変換手段と、
   時間領域における波形の正常モデルを表す第１のプロファイルを前記複数の正常時ワークデータから作成すると共に、周波数領域における波形の正常モデルを表す第２のプロファイルを前記複数の変換後正常時ワークデータから作成するモデル作成手段と、
   前記第１のプロファイル及び前記第２のプロファイルと、前記周期データとに基づいて、所定の指標値を算出する指標値算出手段と、
   前記指標値算出手段により算出された指標値と、予め設定された所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する判定手段と、
   を有することを特徴とする異常検出装置。
2. 前記分割手段は、
   前記周期データを、前記ワーク毎に分割して、複数のワークデータを作成し、
   前記変換手段は、
   前記分割手段により分割された複数のワークデータのそれぞれに対して、窓関数を用いた高速フーリエ変換を行って、前記ワークデータのそれぞれを周波数領域に変換した複数の変換後ワークデータを作成し、
   前記指標値算出手段は、
   前記第１のプロファイルと、前記複数のワークデータとから複数の第１の指標値を算出すると共に、前記第２のプロファイルと、前記複数の変換後ワークデータとから複数の第２の指標値とを算出し、
   前記判定手段は、
   前記複数の第１の指標値と、前記複数の第２の指標値と、予め設定された所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する、ことを特徴とする請求項１に記載の異常検出装置。
3. 前記第１の指標値は、前記第１のプロファイルと前記ワークデータとの相関係数、平均値、中央値、又はマハラノビス距離のいずれかであり、
   前記第２の指標値は、前記第２のプロファイルと前記ワークデータとの相関係数、平均値、中央値、又はマハラノビス距離のいずれかである、ことを特徴とする請求項２に記載の異常検出装置。
4. 前記判定手段は、
   前記複数の第１の指標値のそれぞれと、予め設定された所定の第１の閾値とを比較判定すると共に、前記複数の第２の指標値のそれぞれと、予め設定された所定の第２の閾値とを比較判定することで、前記機器で異常が発生したか否かを判定する、ことを特徴とする請求項２又は３に記載の異常検出装置。
5. 前記判定手段は、
   前記第１の指標値と前記第２の指標値との重み付け和と、予め設定された所定の閾値とを比較判定することで、前記機器で異常が発生したか否かを判定する、ことを特徴とする請求項２又は３に記載の異常検出装置。
6. 機器の周期的な処理又は動作を表す周期データから前記機器の異常を検出する異常検出装置が、
   正常時における前記機器の周期的な処理又は動作を示す正常時周期データを、前記周期の時間幅を示すワーク毎に分割して、複数の正常時ワークデータを作成する分割手順と、
   前記分割手順により分割された複数の正常時ワークデータのそれぞれに対して、窓関数を用いた高速フーリエ変換を行って、前記複数の正常時ワークデータのそれぞれを周波数領域に変換した複数の変換後正常時ワークデータを作成する変換手順と、
   時間領域における波形の正常モデルを表す第１のプロファイルを前記複数の正常時ワークデータから作成すると共に、周波数領域における波形の正常モデルを表す第２のプロファイルを前記複数の変換後正常時ワークデータから作成するモデル作成手順と、
   前記第１のプロファイル及び前記第２のプロファイルと、前記周期データとに基づいて、所定の指標値を算出する指標値算出手順と、
   前記指標値算出手順により算出された指標値と、予め設定された所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する判定手順と、
   を実行することを特徴とする異常検出方法。
7. 機器の周期的な処理又は動作を表す周期データから前記機器の異常を検出する異常検出装置に、
   正常時における前記機器の周期的な処理又は動作を示す正常時周期データを、前記周期の時間幅を示すワーク毎に分割して、複数の正常時ワークデータを作成する分割手順と、
   前記分割手順により分割された複数の正常時ワークデータのそれぞれに対して、窓関数を用いた高速フーリエ変換を行って、前記複数の正常時ワークデータのそれぞれを周波数領域に変換した複数の変換後正常時ワークデータを作成する変換手順と、
   時間領域における波形の正常モデルを表す第１のプロファイルを前記複数の正常時ワークデータから作成すると共に、周波数領域における波形の正常モデルを表す第２のプロファイルを前記複数の変換後正常時ワークデータから作成するモデル作成手順と、
   前記第１のプロファイル及び前記第２のプロファイルと、前記周期データとに基づいて、所定の指標値を算出する指標値算出手順と、
   前記指標値算出手順により算出された指標値と、予め設定された所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する判定手順と、
   を実行させることを特徴とする異常検出プログラム。
8. センシング対象の機器の周期的な処理又は動作を表す周期データをセンシングするセンシング機器と、該センシング機器から送信された周期データから前記機器の異常を検出する異常検出装置とが含まれる異常検出システムであって、
   正常時における前記機器の周期的な処理又は動作を示す正常時周期データを、前記周期の時間幅を示すワーク毎に分割して、複数の正常時ワークデータを作成する分割手段と、
   前記分割手段により分割された複数の正常時ワークデータのそれぞれに対して、窓関数を用いた高速フーリエ変換を行って、前記複数の正常時ワークデータのそれぞれを周波数領域に変換した複数の変換後正常時ワークデータを作成する変換手段と、
   時間領域における波形の正常モデルを表す第１のプロファイルを前記複数の正常時ワークデータから作成すると共に、周波数領域における波形の正常モデルを表す第２のプロファイルを前記複数の変換後正常時ワークデータから作成するモデル作成手段と、
   前記第１のプロファイル及び前記第２のプロファイルと、前記周期データとに基づいて、所定の指標値を算出する指標値算出手段と、
   前記指標値算出手段により算出された指標値と、予め設定された所定の閾値とに基づいて、前記機器で異常が発生したか否かを判定する判定手段と、
   を有することを特徴とする異常検出システム。

Images