JP6641056B1 - Device abnormality diagnosis method and device abnormality diagnosis system - Google Patents
Device abnormality diagnosis method and device abnormality diagnosis system Download PDFInfo
- Publication number
- JP6641056B1 JP6641056B1 JP2019131439A JP2019131439A JP6641056B1 JP 6641056 B1 JP6641056 B1 JP 6641056B1 JP 2019131439 A JP2019131439 A JP 2019131439A JP 2019131439 A JP2019131439 A JP 2019131439A JP 6641056 B1 JP6641056 B1 JP 6641056B1
- Authority
- JP
- Japan
- Prior art keywords
- data
- abnormality
- display
- unit
- dimensions
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
【課題】エンジニアの技能に依存せずに診断を行うことができる機器の異常診断技術を提供する。【解決手段】機器の異常診断方法は、機器の状態を示すデータに含まれる診断の対象となる対象データの次元数が特定次元数になるように対象データを調整するステップと、特定次元数に調整された対象データを、特定次元数のデータを扱うアルゴリズムを用いて解析するステップと、対象データの解析に基づいて、異常の有無を識別可能な態様で対象データを表示するステップとを含む。【選択図】図12To provide a device abnormality diagnosis technique capable of performing diagnosis without depending on the skill of an engineer. An apparatus abnormality diagnosis method includes the steps of: adjusting target data so that the number of dimensions of target data to be diagnosed included in data indicating the state of the device is a specific dimension; The method includes a step of analyzing the adjusted target data using an algorithm that handles data of a specific number of dimensions, and a step of displaying the target data in a form in which the presence or absence of an abnormality can be identified based on the analysis of the target data. [Selection diagram] FIG.
Description
本発明の実施形態は、機器の異常診断方法および機器の異常診断システムに関する。 An embodiment of the present invention relates to a device abnormality diagnosis method and a device abnormality diagnosis system.
プラントに設けられる機器が計画外に停止または故障すると、プラントの運転に影響し、損害が大きくなる。そのため、機器の状態を監視し、その健全性を適切に評価することが求められている。また、評価の結果、異常の兆候がある場合には、その原因を明確にして早期に対策を施すことが求められる。 If the equipment provided in the plant stops or fails unplannedly, the operation of the plant is affected and the damage is increased. Therefore, it is required to monitor the state of the device and appropriately evaluate its soundness. Also, if there is a sign of abnormality as a result of the evaluation, it is necessary to clarify the cause and take an early action.
プラント機器としては、モータまたはポンプなどの回転機器、配管、弁、ケーブルなどがある。これらの健全性は、機器の種類に応じて、振動、温度、電流、圧力、流量、音などを測定し、その変化を監視することで評価する。しかし、プラント内の監視対象となる機器の数量または監視項目が多いため、個々の機器を人間が精密に診断することが難しい。また、機器によってはノイズまたは個体差が大きく、適切な診断には、エンジニアの経験およびスキルを要する。 Plant equipment includes rotating equipment such as a motor or a pump, piping, valves, cables, and the like. These soundnesses are evaluated by measuring vibration, temperature, current, pressure, flow rate, sound, etc. according to the type of equipment and monitoring the changes. However, it is difficult for a human to precisely diagnose each device because the number of devices to be monitored in the plant or the number of monitoring items is large. Further, depending on the device, noise or individual difference is large, and an appropriate diagnosis requires experience and skill of an engineer.
このため、技術者の熟練度に依らず効率的に診断を行う方法が考えられている。例えば、回転機器の異常時における振動または音響スペクトルの時間変化と、ピーク周波数の組を予め基準データとして登録し、測定データのスペクトルと照合することで異常の有無または原因特定する技術が知られている。また、電動弁の駆動信号をセンシングして集録し、予め設定された診断項目の許容値および診断記録を用いて異常診断を行う技術が知られている。また、近年では、AIを用いた診断手法による効率的かつ効果的な異常診断技術の開発が進められている。例えば、正常または異常データを蓄積し、機械学習によって異常の有無を判断する技術が知られている。 For this reason, a method for performing an efficient diagnosis regardless of the skill level of a technician has been considered. For example, a technique is known in which a set of a time change of a vibration or an acoustic spectrum at the time of abnormality of a rotating device and a peak frequency is registered in advance as reference data, and the presence or absence of an abnormality or a cause is identified by collating with a spectrum of measurement data. I have. Further, there is known a technique of sensing and collecting a drive signal of a motor-operated valve and performing abnormality diagnosis using a preset allowable value of a diagnostic item and a diagnostic record. In recent years, development of an efficient and effective abnormality diagnosis technique using a diagnosis method using AI has been promoted. For example, a technique is known in which normal or abnormal data is accumulated and the presence or absence of an abnormality is determined by machine learning.
前述の技術では、正常と異常の判定に、基準となる閾値が用いられている。その基準については、機器ごとにエンジニアが基準値を設ける必要がある。この基準値の設定は、熟練のエンジニアが過去の測定データまたは対象機器の仕様を考慮しながら行うため、時間と労力を要する。このため、従来の手法ではエンジニアの技能の経験およびスキルに対する依存性が残っている。 In the above-described technique, a threshold value serving as a reference is used for determining normality or abnormality. As for the standard, an engineer needs to set a standard value for each device. The setting of this reference value is performed by a skilled engineer while considering past measurement data or the specifications of the target device, and therefore requires time and effort. For this reason, the conventional technique still has a dependency on the experience and skill of the engineer.
本発明の実施形態は、このような事情を考慮してなされたもので、エンジニアの技能に依存せずに診断を行うことができる機器の異常診断技術を提供することを目的とする。 An embodiment of the present invention has been made in consideration of such circumstances, and has as its object to provide a device abnormality diagnosis technique capable of performing a diagnosis without depending on the skill of an engineer.
本発明の実施形態に係る機器の異常診断方法は、機器の状態を示すデータである信号の時間波形と周波数スペクトルのうちの少なくとも一方のデータから抽出され、前記機器の状態を示すデータのうち診断の対象となる1の対象データに含まれる個々のデータの個数である次元数を特定次元数に調整するステップと、前記特定次元数に調整された前記対象データを、前記特定次元数のデータを扱うアルゴリズムであるマッパーを用いて位相的データ解析を行うステップと、前記マッパーを用いた解析結果を異常の有無を識別可能な態様であるクラスタ表示または時系列表示で表示するステップと、を含む。 The device abnormality diagnosis method according to the embodiment of the present invention is a method for extracting, from data indicating at least one of a time waveform and a frequency spectrum of a signal which is data indicating a state of the device, a diagnosis among data indicating the state of the device. individual comprising the steps of adjustment to a specific number of dimensions the number of dimensions is the number of data, the target data adjusted to the specific number dimension, the specific number of dimensions of data included in the first target data to be Performing a topological data analysis using a mapper , which is an algorithm that handles the data processing, and displaying the analysis result using the mapper in a cluster display or a time-series display in a mode in which the presence or absence of an abnormality can be identified. .
本発明の実施形態により、エンジニアの技能に依存せずに診断を行うことができる機器の異常診断技術が提供される。 According to the embodiments of the present invention, a device abnormality diagnosis technique capable of performing diagnosis without depending on the skill of an engineer is provided.
以下、図面を参照しながら、機器の異常診断方法および機器の異常診断システムの実施形態について詳細に説明する。 Hereinafter, embodiments of a device abnormality diagnosis method and a device abnormality diagnosis system will be described in detail with reference to the drawings.
図1の符号1は、本実施形態の機器の異常診断システムである。この異常診断システム1は、発電プラントなどの施設に設けられる機器としての回転機2の異常を診断するシステムである。なお、回転機2とは、例えば、ポンプまたは送風機などに搭載されるモータである。
例えば、回転機2に故障などの異常がある場合、または異常の兆候がある場合には、その振動の加速度振幅または特定の周波数成分が、正常時と比較して異なるものとなる。正常時と異常時との差異が大きなものであれば、異常の発見がし易いが、その差異が微小なものであれば、異常の発見が困難である。また、機器としての回転機2の運転中に生じる現象に関するデータは、膨大な量となり、これら大量のデータを人手で解析することは困難である。
For example, when there is an abnormality such as a failure in the rotating
近年、機械学習および人工知能に関する技術の発達に伴い、大量のデータから有用な情報を取り出すことが可能になっている。このようなデータ解析のアルゴリズムとしては、例えば、位相的データ解析(TDA:Topological Data Analysis)またはクラスタ解析(Clustering)がある。本実施形態の異常診断システム1では、大量のデータを扱うことができるアルゴリズムを用いて、回転機2(機器)の異常の診断を行う。
In recent years, with the development of technologies related to machine learning and artificial intelligence, it has become possible to extract useful information from a large amount of data. As an algorithm for such data analysis, for example, there is a topological data analysis (TDA) or a cluster analysis (Clustering). The
従来の解析技術では、大規模、多次元のデータを扱う場合に次元の削減の処理が必要である。例えば、統計的手法による次元削減では、データ中の情報が失われてしまうといった課題がある。また、教師あり学習を用いた異常診断についても、異常となるデータラベルが少ないと学習モデルを構築することが難しいという課題がある。 In the conventional analysis technology, processing of dimension reduction is necessary when handling large-scale, multidimensional data. For example, in dimension reduction by a statistical method, there is a problem that information in data is lost. In addition, the abnormality diagnosis using supervised learning also has a problem that it is difficult to construct a learning model if there are few abnormal data labels.
そこで、本実施形態では、教師なしでデータ分析ができ、かつデータ群のトポロジカルな特徴に着目したデータ解析技術を用いて、データ群の中から異常の診断を行うようにする。 Therefore, in the present embodiment, an abnormality is diagnosed from a data group by using a data analysis technique that can analyze data without a teacher and focus on topological features of the data group.
例えば、従来の位相的データ解析では、分析のためのパラメータ設定または分析結果の解釈に人手による判断が必要である。これに対して、本実施形態では、位相的データ解析で用いる分析パラメータをスイープし、それに対する出力結果(例えば、データ数、エッジ数、ノード群数)を分析することで特徴のあるデータ(異常データ)を抽出することができる。また、どのパラメータが異常事象に寄与しているかを推定することができる。 For example, in conventional topological data analysis, manual judgment is required for setting parameters for analysis or interpreting analysis results. On the other hand, in the present embodiment, characteristic data (abnormality) is obtained by sweeping analysis parameters used in topological data analysis and analyzing output results (for example, the number of data, the number of edges, and the number of nodes). Data) can be extracted. Further, it is possible to estimate which parameter contributes to the abnormal event.
図1に示すように、異常診断システム1の使用者は、プラントに設けられた回転機2に点検用センサ3(例えば、振動計)を取り付けて、回転機2の振動データの取得を行う。なお、回転機2の運転中の状態を示すデータであるモニタリングデータとしては、振動データの他にも、動作音を収集するセンサ、サーモグラフィなどのデータがある。なお、点検用センサ3は、常設のものであっても良いし、測定の都度、設置するものでも良い。
As shown in FIG. 1, a user of the
モニタリングデータは、機器としての回転機2に取り付けられた点検用センサ3から発生する信号の時間波形と周波数スペクトルとのうちの少なくとも1つのデータである。また、モニタリングデータは、信号の時間波形または周波数スペクトルの絶対値でも良いし、信号の特徴を代表する量またはその組み合わせであっても良い。
The monitoring data is at least one of a time waveform and a frequency spectrum of a signal generated from the
なお、振動データで対象となる状態量としては、例えば、変位、速度、加速度などがある。また、サンプリング点数または周波数の範囲など、任意の測定条件を使用者が選択可能である。このようにすれば、点検用センサ3で容易に取得できる振動データに基づいて、回転機2の異常の診断を行うことができる。
It should be noted that the state quantities targeted by the vibration data include, for example, displacement, speed, acceleration, and the like. Further, the user can select an arbitrary measurement condition such as the number of sampling points or the frequency range. In this way, the abnormality of the
モニタリングデータとしては、例えば、点検用センサ3で計測された振動データを例示する。なお、モニタリングデータには、常時、回転機2の状態をモニタリングしている連続的なデータの他に、定期的に測定した断続的なデータも含まれる。以下の説明では、振動データを例に説明するが、その他のデータでも良い。
As the monitoring data, for example, vibration data measured by the
なお、モニタリングデータには、回転機2が動作することで達成される主目的の事象を示す1次データの他に、回転機2が動作することで生じる副次的な2次データが含まれる。例えば、回転機2の場合には、その回転数が1次データとなる。そして、その回転動作により生じる主目的の事象以外の振動データが2次データとなる。本実施形態では、2次データを用いて回転機2の異常の診断を行う態様を例示する。なお、1次データを用いて回転機2の異常の診断を行っても良い。
It should be noted that the monitoring data includes, in addition to the primary data indicating the main purpose event achieved by the operation of the
そして、このモニタリングデータを異常診断システム1に入力して解析を行う。例えば、モニタリングデータは、インターネット4に接続されたユーザ端末5(例えば、PC)を介して異常診断システム1に入力される。なお、モニタリングデータは、異常診断システム1に直接入力されても良い。
Then, the monitoring data is input to the
この異常診断システム1は、プラント外のクラウドサーバに実装されても良いし、プラント内のパーソナルコンピュータまたはサーバに実装されても良い。また、モニタリングデータを、サーバなどの電気機器または所定のメモリに収集した後に、これら収集したデータを異常診断システム1に転送しても良い。
This
点検用センサ3は、例えば、X軸、Y軸、Z軸の3軸方向の加速度を検出する。これにより、回転機2の加速度の振動(データ)波形を取得することができる。
The
次に、異常診断システム1のシステム構成を図1から図8に示すブロック図を参照して説明する。
Next, the system configuration of the
図1に示すように、異常診断システム1は、メイン制御部6と情報入力部7と通信部8と表示出力部9と解析処理ユニット10とを備える。これらは、メモリまたはHDDに記憶されたプログラムがCPUによって実行されることで実現される。
As shown in FIG. 1, the
さらに、異常診断システム1は、モニタリングデータベース11とプロセスデータベース12とスペックデータベース13と異常原因データベース14とを備える。これらは、メモリまたはHDDに記憶され、検索または蓄積ができるよう整理された情報の集まりである。
Further, the
本実施形態の異常診断システム1は、CPU、ROM、RAM、HDDなどのハードウェア資源を有し、CPUが各種プログラムを実行することで、ソフトウェアによる情報処理がハードウェア資源を用いて実現されるコンピュータで構成される。さらに、本実施形態の異常診断方法は、プログラムをコンピュータに実行させることで実現される。
The
メイン制御部6は、異常診断システム1を統括的に制御する。
The
情報入力部7は、異常診断システム1を使用する使用者の操作に応じて所定の情報が入力される。この情報入力部7には、マウスまたはキーボードなどの入力装置が含まれる。つまり、これら入力装置の操作に応じて所定の情報が情報入力部7に入力される。この情報入力部7には、モニタリングデータなどの所定の情報が入力される。
The
通信部8は、インターネット4などの通信回線を介してユーザ端末5などの電子機器と通信を行う。この通信部8は、無線通信装置に搭載されても良い。また、通信部8は、所定のネットワーク機器、例えば、無線LANアクセスポイントまたはアンテナに搭載されても良い。なお、通信部8は、WAN(Wide Area Network)または携帯通信網を介してユーザ端末5と通信を行っても良い。
The
表示出力部9は、解析結果を示す各種情報を表示する。例えば、クラスタ表示(図12参照)、時系列表示(図13参照)などの画面が表示される。つまり、表示出力部9は、分析結果を表示するためのグラフィックユーザインターフェース(GUI)である。
The
本実施形態の異常診断システム1には、解析結果の出力を行うディスプレイなどの表示装置が含まれる。つまり、表示出力部9は、ディスプレイに表示される画面の制御を行う。なお、ディスプレイはコンピュータ本体と別体であっても良いし、一体であっても良い。さらに、ネットワークを介して接続される他のコンピュータが備えるディスプレイに表示される画面の制御を表示出力部9が行っても良い。
The
なお、本実施形態では、表示装置としてディスプレイを例示するが、その他の態様であっても良い。例えば、ヘッドマウントディスプレイまたはプロジェクタを用いて情報の表示を行っても良い。さらに、紙媒体に情報を印字するプリンタをディスプレイの替りとして用いても良い。つまり、表示出力部9が制御する対象としてヘッドマウントディスプレイ、プロジェクタまたはプリンタが含まれても良い。
In the present embodiment, a display is exemplified as a display device, but other embodiments may be used. For example, information may be displayed using a head-mounted display or a projector. Further, a printer that prints information on a paper medium may be used instead of a display. That is, a target mounted on the
解析処理ユニット10は、特定次元数のデータを扱うアルゴリズムを用いてモニタリングデータに含まれる対象データまたは過去データの解析を行う(図2参照)。なお、本実施形態では、回転機2の運転中の状態を示すデータとして、モニタリングデータを例示しているが、モニタリングデータとプロセスデータを合せて作成されたデータセットが、回転機2の運転中の状態を示すデータであっても良い。以下の説明において、モニタリングデータという用語には、プロセスデータが含まれる場合がある。
The
ここで、本実施形態の特定次元数は、4次元以上である。なお、例えば、解析処理ユニット10は、位相的データ解析とクラスタ解析とのうち少なくともいずれか一方を用いて解析を行う。このようにすれば、特定次元数のデータを扱った解析を行うことができる。
Here, the number of specific dimensions in the present embodiment is four or more. Note that, for example, the
対象データとは、モニタリングデータに含まれる診断の対象となるデータである。なお、モニタリングデータと同時期または前後の時期に取得したプロセスデータも対象データに含まれる場合がある。 The target data is data to be diagnosed included in the monitoring data. Note that the process data acquired at the same time as or around the monitoring data may be included in the target data.
過去データとは、モニタリングデータに含まれる異常の有無が既知のデータである。例えば、現在運転中の回転機2から取得されたモニタリングデータは、異常の有無が未知である。
The past data is data in which the presence or absence of an abnormality included in the monitoring data is known. For example, monitoring data acquired from the rotating
このモニタリングデータは、対象データとして異常の有無の診断の対象となる。また、数カ月前または数年前に回転機2から取得されたモニタリングデータは、そのデータが取得された後に、実際に回転機2が正常に動作したか、それとも異常な状態となったかが把握できている。このようなモニタリングデータは、過去データとして扱われる。なお、過去のモニタリングデータと同時期または前後の時期に取得したプロセスデータも過去データに含まれる場合がある。
This monitoring data is to be diagnosed as abnormal data as target data. In addition, the monitoring data acquired from the rotating
モニタリングデータベース11は、使用者が点検用センサ3を用いて回転機2から取得したモニタリングデータが記憶される。
The
プロセスデータベース12は、回転機2に関連するプロセスを示すプロセスデータが記憶される。プロセスデータには、例えば、回転数、流量、圧力、温度、電流などの運転条件と、運転期間、使用期間などの使用条件などが含まれる。なお、このプロセスデータは、使用者が記憶させても良いし、予め記憶されているものでも良い。
The
スペックデータベース13は、回転機2のスペックを示すスペックデータが記憶される。スペックデータには、例えば、回転機2の場合は、型式、メーカー、サイズ、揚程、流量種別、軸受形式、羽数などの情報が含まれる。なお、このスペックデータは、使用者が記憶させても良いし、予め記憶されているものでも良い。
The
異常原因データベース14は、解析処理ユニット10の解析に基づいて抽出された対象データの特徴量と、回転機2の異常原因との対応関係を示す原因情報(異常原因データ)が記憶される。原因情報は、異常の特徴量と異常原因との相関関係を示すものである。
The
なお、この原因情報は、使用者が相関関係を入力、設定しても良いし、予め記憶されているものでも良い。このようにすれば、特徴量と異常原因との対応関係に基づいて、回転機2の異常原因の推定を容易に行うことができる。
The cause information may be input and set by the user, or may be stored in advance. In this manner, the cause of the abnormality of the
図2に示すように、解析処理ユニット10は、データ取得部15とデータ前処理部16と処理設定部17とデータ解析部18とパラメータチューニング部19と結果整理部20と結果表示制御部21と表示調整部22と異常原因分析部23と原因情報取得部24と原因表示制御部25とを備える。
As shown in FIG. 2, the
データ取得部15は、モニタリングデータベース11とプロセスデータベース12とスペックデータベース13とから各種データを取得する。なお、以下の説明において、モニタリングデータという用語には、データ取得部15が取得したプロセスデータまたはスペックデータが含まれる場合がある。
The
データ前処理部16は、モニタリングデータから、対象データまたは過去データなどの解析に用いるデータを抽出し、データ解析部18で解析される前に、これらデータの前処理を行う(図3参照)。ここで、対象データまたは過去データの次元数が特定次元数になるように調整される。なお、対象データと過去データとは、同一の特定次元数になるように調整される。
The
処理設定部17は、データ解析部18で解析を行うときに用いる各種設定または解析用のアルゴリズムの設定を行う(図4参照)。
The
データ解析部18は、特定次元数に調整された対象データまたは過去データを、特定次元数のデータを扱うアルゴリズムである位相的データ解析とクラスタ解析とのうち少なくともいずれか一方を用いて解析する。このデータ解析部18は、モニタリングデータに含まれる対象データの中から、異常な状態を示す対象データを検出することができる。
The
パラメータチューニング部19は、データ解析部18で解析を行うときに用いる各種設定または解析用のアルゴリズムの設定のチューニングを行う(図5参照)。
The
結果整理部20は、データ解析部18により得られたデータの解析結果が最適なものであるか否かを整理する。つまり、データ解析部18で分析した結果について、分析パラメータと分析結果との対応関係を整理する。
The
ここで、データの解析結果が最適でない場合は、パラメータチューニング部19にて設定のチューニングが行われる。つまり、データ解析部18によるデータの解析に基づいて、回転機2の異常の有無を識別可能な態様で対象データを表示可能な解析結果が得られるまで、パラメータチューニング部19によるチューニングとデータ解析部18によるデータの解析とが繰り返される。
If the data analysis result is not optimal, the
結果表示制御部21は、データ解析部18により得られたデータの解析結果を示す画面であって、回転機2の異常の有無を識別可能な態様の画面を表示出力部9に表示させる制御を行う(図6参照)。例えば、クラスタ表示(図12参照)、時系列表示(図13参照)などの画面を表示出力部9に表示させる制御を行う。
The result
ここで、本実施形態の異常の有無を識別可能な態様は、例えば、3次元以下のグラフである。例えば、クラスタ表示(図12参照)を2次元のグラフとして例示する。このクラスタ表示が3次元のグラフであっても良い。このようにすれば、通常の表示技術では可視化できない4次元以上の対象データを、3次元以下のグラフとして可視化し、回転機2の異常の有無の判定に用いることができる。なお、時系列表示(図13参照)は、異常の有無を識別可能な態様に含まれる。
Here, an aspect of the present embodiment in which the presence or absence of an abnormality can be identified is, for example, a three-dimensional or lower graph. For example, the cluster display (see FIG. 12) is illustrated as a two-dimensional graph. This cluster display may be a three-dimensional graph. With this configuration, target data of four or more dimensions that cannot be visualized by a normal display technique can be visualized as a graph of three or less dimensions, and can be used to determine whether there is an abnormality in the
位相的データ解析を用いる場合において、結果表示制御部21は、ノード数、エッジ数、入力したデータの総数、ノードの中に含まれるユニークなデータ数、ノード群数などを含む異常検出結果の整理を行う。なお、以下の説明では、ノード群をクラスタと称する場合がある。また、結果表示制御部21は、距離のメトリック、フィルター関数、インターバル、オーバーラップ、クラスタリングの手法などを含む位相的データ解析の分析パラメータの整理も行う。
When the topological data analysis is used, the result
表示調整部22は、使用者の操作に応じて表示出力部9に表示させる画面の制御を行う(図7参照)。
The
異常原因分析部23は、データ解析部18により得られたデータの解析結果に基づいて、対象データに異常がある場合に、その異常原因の分析を行う(図8参照)。
The abnormality
原因情報取得部24は、異常原因データベース14から原因情報を取得する。この原因情報に基づいて、異常原因分析部23にて分析が行われる。
The cause
原因表示制御部25は、異常原因分析部23により得られたデータの解析結果を示す画面であって、回転機2の異常原因を使用者が診断可能な態様の画面を表示出力部9に表示させる制御を行う。
The cause
図3に示すように、データ前処理部16は、対象データ調整部26と過去データ調整部27と欠損情報補間部28とデータラベル付与部29とデータスケーリング部30とスペクトル解析部31と振幅値算出部32とデータフォーマット変換部33とを備える。
As shown in FIG. 3, the
対象データ調整部26は、モニタリングデータまたはプロセスデータに含まれる診断の対象となる対象データの次元数が特定次元数になるように対象データを調整する。
The target
過去データ調整部27は、モニタリングデータまたはプロセスデータに含まれる異常の有無が既知の過去データの次元数が特定次元数になるように過去データを調整する。
The past
なお、本実施形態の「特定次元数になるようにデータを調整する」とは、任意の次元になるように、使用者がデータを間引く、または、次元圧縮することが含まれ得る。さらに、回転機2の振動に対する物理的な知識に基づいて、分析対象のデータの中から、異常に寄与し易いデータのみを抜き出すことが含まれる。例えば、2000次元の周波数のデータがある場合に、その中から異常に寄与し易い48次元のデータのみを抽出する(次元を調整する)前処理を行っても良い。
Note that “adjusting the data to have a specific number of dimensions” in the present embodiment may include the user thinning out the data or compressing the dimension to have an arbitrary dimension. Furthermore, extracting only data that easily contributes to an abnormality from the data to be analyzed based on physical knowledge of the vibration of the
欠損情報補間部28は、対象データまたは過去データを特定次元数にしたときに、それぞれの次元数の情報のうち一部の次元数の情報が欠損する場合に、その欠損した情報を補間する。このようにすれば、対象データまたは過去データのそれぞれの次元数の情報を揃えて、特定次元数のデータを扱うアルゴリズムで解析可能な対象データまたは過去データを作成することができる。
The missing
データラベル付与部29は、特定次元数が調整された対象データまたは過去データに所定のラベルを付与する。例えば、過去に分析したデータ、または、データ収集の際に有意な特徴が観測されたデータに関しては、正常のときのデータか、異常が有るときのデータかが既知ものがある。このような過去データには、異常の有無を識別可能な情報としてのラベルLが付与される(図9および図10参照)。つまり、過去データには、異常の有無を識別可能なラベルLが対応付けられる。なお、対象データは、ラベルLをブランクにする。
The data label assigning
データスケーリング部30は、モニタリングデータを正規化または標準化する処理を行う。この処理により、次元毎でスケールの異なるデータのスケールを揃えることができる。例えば、正規化処理を行うと、特徴量の値の範囲が一定の範囲に収まるようになる。主に[0,1]か、[−1,1]の範囲内に収めることが多い。また、標準化では特徴量の平均を0、分散を1にする変換を行う。
The
スペクトル解析部31は、振動データから任意の周波数成分を抽出する。この解析に基づいて、対象データまたは過去データになる振動スペクトルデータを生成することができる。このスペクトル解析部31は、例えば、高速フーリエ変換(FFT:Fast Fourier Transform)を用いて、モニタリングデータの特定の周波数成分を抽出したり、包絡線処理を行い、特定の周期的な周波数成分を抽出したりする。
The
振幅値算出部32は、モニタリングデータに含まれる振幅値を解析する。なお、この解析に基づいて、対象データまたは過去データになる振幅値データを作成することができる。
The
つまり、本実施形態では、モニタリングデータに基づいて、回転機2の振動に関する情報である振幅値データまたは振動スペクトルデータが取得可能である。
That is, in the present embodiment, amplitude value data or vibration spectrum data, which is information relating to vibration of the
データフォーマット変換部33は、データ解析部18で用いるアルゴリズムに合わせたデータフォーマットになるように、対象データまたは過去データのデータフォーマットを変換する。なお、データフォーマットの変換とは、異常分析を行うためにデータの形状を調整することである。
The data
図4に示すように、処理設定部17は、前処理条件設定部34とアルゴリズム設定部35とを備える。
As shown in FIG. 4, the
前処理条件設定部34は、前処理の条件を設定する。例えば、前述の補間または規格化の処理などの設定が変更される。
The pre-processing
アルゴリズム設定部35は、データ解析部18で扱うアルゴリズムの設定を行う。例えば、位相的データ解析を用いるか、クラスタ解析を用いるか、それとも両解析を用いるかの設定を行う。またそれらのアルゴリズムのハイパーパラメータの設定を行う。このハイパーパラメータの設定により、例えば、位相的解析を用いる場合は、フィルター関数、オーバーラップ率などが変更される。
The
図5に示すように、パラメータチューニング部19は、特定次元数設定部36とパラメータ最適化部38とを備える。
As shown in FIG. 5, the
特定次元数設定部36は、対象データまたは過去データの解析を行うときに用いる特定次元数の設定を行う。このようにすれば、異常の有無を識別可能な態様になるように特定次元数の設定を行うことができる。なお、特定次元数の設定は、前処理条件設定部34またはアルゴリズム設定部35で行っても良い。
The specific dimension
また、解析後に得られる対象データまたは過去データの表示が、異常の有無を識別可能な態様になるように特定次元数の設定が繰り返される。このようにすれば、異常の有無を識別可能な態様になるまで、特定次元数の設定を繰り返し行うことができる。この設定は、人手で行っても良いし、自動的に行っても良い。 In addition, the setting of the specific dimension number is repeated so that the display of the target data or the past data obtained after the analysis has a form in which the presence or absence of the abnormality can be identified. In this way, the setting of the specific dimension number can be repeatedly performed until the presence / absence of the abnormality can be identified. This setting may be performed manually or automatically.
パラメータ最適化部38は、結果整理部20で整理された分析パラメータと分析結果との対応関係に基づいて、データ解析部18にて分析に用いるパラメータの設定をチューニングする。このチューニングでは、解析後に得られる対象データまたは過去データの表示が、異常の有無を識別可能な態様になるように、パラメータの変更を繰り返し、パラメータの最適化を行う。
The
図6に示すように、結果表示制御部21は、クラスタ表示制御部39と時系列表示制御部40とを備える。
As shown in FIG. 6, the result
クラスタ表示制御部39は、データ解析部18により得られたデータの解析結果に基づいて、クラスタ表示(図12参照)の画面を表示出力部9に表示させる制御を行う。なお、クラスタ表示は、対象データをノードNとして表示した画面である。このようにすれば、対象データをクラスタ表示として視覚化できるため、回転機2の異常の有無の識別が行い易くなる。
The cluster
時系列表示制御部40は、データ解析部18により得られたデータの解析結果を示す時系列表示(図13参照)の画面を表示出力部9に表示させる制御を行う。なお、時系列表示は、対象データをモニタリングされた時系列順に表示し、かつ対象データが含まれるノードNを識別可能な態様で表示した画面である。このようにすれば、対象データを時系列に並べて表示できるため、回転機2の異常の有無の識別が行い易くなる。
The time-series
図7に示すように、表示調整部22は、ノード選択受付部41と対象情報表示部42とを備える。
As illustrated in FIG. 7, the
ノード選択受付部41は、クラスタ表示された画面上のノードNまたはクラスタCの選択を受け付ける処理を行う。
The node
対象情報表示部42は、ノード選択受付部41により受け付けられ、選択されたノードNに含まれる内訳情報Bを表示する(図12参照)。この内訳情報Bには、選択されたノードNに含まれる対象データの対象情報が含まれる。このようにすれば、ノードNに対応する対象データの情報を容易に得ることができる。
The target
図8に示すように、異常原因分析部23は、特徴量抽出部43と異常原因推定部44とを備える。
As illustrated in FIG. 8, the abnormality
特徴量抽出部43は、データ解析部18で解析された後の対象データの特徴量を抽出する。
The feature
異常原因推定部44は、特徴量抽出部43にて抽出された特徴量に基づいて、異常原因を推定する。つまり、特徴量抽出部43にて抽出された特徴量に基づいて、異常原因を診断可能な態様で対象データを表示可能な状態になるようにデータを整理する。このようにすれば、回転機2に異常が発見された場合に、その原因を推定することができる。
The abnormality
次に、データ解析部18で解析される対象データまたは過去データについて図9および図10を用いて説明する。
Next, target data or past data analyzed by the
図9は、データ前処理部16にて前処理された後の加速度振幅のデータの一例を含むテーブルである。このテーブルでは、モニタリングされた日付(測定日、データ収集日時)に対応付けて、特定次元数の加速度振幅が登録されている。なお、加速度振幅の特定次元数は、任意の数を選択可能である。ここで、最新の日付のデータが、異常の有無の診断の対象となる対象データTとなっている。
FIG. 9 is a table including an example of acceleration amplitude data that has been preprocessed by the
なお、特定次元数の加速度振幅の列Sが複数設けられている。これら列Sの数がデータの次元数となっている。図9のテーブルでは、所定の日付(測定日、データ収集日時)に対応付けて、5次元以上、例えば、50次元の加速度振幅が登録されている。また、過去データPに対応付けて、異常の有無を識別可能なラベルLが登録される。例えば、正常であった過去データPに「A」のラベルLが登録され、異常であった過去データPに「B」のラベルLが登録されている。 In addition, a plurality of columns S of the acceleration amplitude of a specific dimension number are provided. The number of these columns S is the number of dimensions of the data. In the table of FIG. 9, five-dimensional or more, for example, 50-dimensional acceleration amplitude is registered in association with a predetermined date (measurement date, data collection date and time). Further, a label L capable of identifying the presence or absence of an abnormality is registered in association with the past data P. For example, the label L of “A” is registered in the normal past data P, and the label L of “B” is registered in the abnormal past data P.
図10は、データ前処理部16にて前処理された後の振動スペクトルデータの一例を含むテーブルである。このテーブルでは、モニタリングされた日付に対応付けて、特定次元数ごとの振動値が登録されている。なお、周波数の特定次元数は、任意の数を選択可能である。ここで、最新の日付のデータが、異常の有無の診断の対象となる対象データTとなっている。
FIG. 10 is a table including an example of vibration spectrum data that has been preprocessed by the
なお、特定次元数の周波数の列Fが複数設けられている。これら列Fの数がデータの次元数となっている。図10のテーブルでは、所定の日付に対応付けて、5次元以上、例えば、50次元の加速度スペクトルが登録される。また、過去データPに対応付けて、異常の有無を識別可能なラベルLが登録される。例えば、正常であった過去データPに「A」のラベルLが登録され、異常であった過去データPに「B」のラベルLが登録されている。 Note that a plurality of frequency columns F having a specific number of dimensions are provided. The number of these columns F is the number of dimensions of the data. In the table of FIG. 10, five-dimensional or more, for example, 50-dimensional acceleration spectra are registered in association with a predetermined date. Further, a label L capable of identifying the presence or absence of an abnormality is registered in association with the past data P. For example, the label L of “A” is registered in the normal past data P, and the label L of “B” is registered in the abnormal past data P.
図9および図10において、対象データおよび過去データの特定次元数は統一される。ここで、それぞれの日付に対応付けて、次元数のそれぞれに当てはまるデータが抽出できない場合がある。例えば、50次元のデータのうち、数次元分のデータが欠損する場合がある。この欠損したデータは、欠損情報補間部28により補間される。なお、補間されるデータの値は、ゼロであっても良いし、欠損した項目と同じ行に含まれているデータの平均値でも良いし、この平均値の前後のデータをもとに内挿または外挿したものでも良いし、欠損した項目と同じ列に含まれているデータの平均値でも良いし、この平均値の前後のデータをもとに内挿また外挿したものでも良いし、その他の任意に定めた値であっても良い。
9 and 10, the number of specific dimensions of the target data and the past data is unified. Here, there are cases where data corresponding to each of the dimensions cannot be extracted in association with each date. For example, data of several dimensions out of 50-dimensional data may be lost. The missing data is interpolated by the missing
なお、本実施形態の過去データとしては、異常の診断対象となる回転機2から取得されたデータを用いる形態を例示するが、その他の態様であっても良い。例えば、メーカー、スペック、型式が同一構成の回転機2がある場合には、この他の回転機2から取得されたデータを、異常の診断対象となる回転機2の過去データとして用いても良い。
As the past data of the present embodiment, a mode using data acquired from the rotating
なお、特定次元数は、全ての日付で同一の次元数となるように、データ前処理部16においてデータ補完およびフォーマット変換が行われる。なお、データ前処理部16の処理は、自動的に行っても良いし、人手により任意に設定しても良い。
The
従来の機械学習を用いた統計的手法では、統計的な解析をベースにしており、取り扱うデータが正規分布などの何らかの分布に従っていることを前提としている。そのため、個別データの特徴に応じた解析手法の構築が必要ある。また、大規模、多次元のデータを扱う場合には、データの次元の削減などの処理が必要である。このような統計的手法による次元削減では、データ中の情報が失われてしまうといった課題がある。本実施形態は、このような課題を解決することができる。 A conventional statistical method using machine learning is based on a statistical analysis, and it is assumed that data to be handled follows some distribution such as a normal distribution. Therefore, it is necessary to construct an analysis method according to the characteristics of individual data. Further, when handling large-scale, multi-dimensional data, processing such as reduction of the dimension of the data is required. In the dimension reduction by such a statistical method, there is a problem that information in data is lost. The present embodiment can solve such a problem.
次に、本実施形態で用いる位相的データ解析、所謂トポロジカルデータアナリシス(以下「TDA」という。)について説明する。TDAは、大規模な多次元データに対して、統計的な特徴ではなく、データの位相幾何学的な情報(データの形)に注目して、データの特徴を分析する手法である。このTDAには、パーシステントホモロジー(Persistent homology)と、マッパー(Mapper)という2種類の手法がある。 Next, a topological data analysis, a so-called topological data analysis (hereinafter, referred to as “TDA”) used in the present embodiment will be described. TDA is a technique for analyzing data features of large-scale multi-dimensional data by focusing on topological information (data form) of the data instead of statistical features. In this TDA, there are two types of methods, a persistent homology and a mapper.
パーシステントホモロジーは、幾何学的な不変量をとらえることができる手法で、データ形状の詳細を定量的に評価することができる。また、マッパーは、幾何学的情報を把握できるシンプルなグラフで、高次元のデータを表現する手法である。 Persistent homology is a technique that can capture geometric invariants, and can quantitatively evaluate the details of data shapes. The mapper is a technique for expressing high-dimensional data with a simple graph that can grasp geometric information.
本実施形態では、TDA Mapperを用いる形態を例示する。なお、その他の解析技術を用いても良い。 In the present embodiment, a mode using TDA Mapper will be exemplified. Note that other analysis techniques may be used.
TDA Mapperにより、高次元数(特定次元数)のデータを、3次元以下のグラフとして可視化することができる。また、データを把握する上で重要な部分を、まとめて1つのノードNとして表現することができる。例えば、連続しているデータを持つノードN同士をエッジEで繋ぎ、データの集合をグラフに変換する(図12参照)。 With TDA Mapper, data with a high number of dimensions (specific number of dimensions) can be visualized as a graph of three dimensions or less. In addition, important parts for grasping data can be collectively expressed as one node N. For example, nodes N having continuous data are connected by an edge E, and a set of data is converted into a graph (see FIG. 12).
なお、本実施形態の「可視化」とは、人間が形状を認識できる態様のことを示す。例えば、ディスプレイに表示された3次元のグラフは、2次元の画像であるが、人間はその立体的な形状を認識することが可能である。 Note that “visualization” in the present embodiment indicates a mode in which a human can recognize a shape. For example, a three-dimensional graph displayed on a display is a two-dimensional image, but a human can recognize the three-dimensional shape.
なお、TDA Mapperには、対象データの座標系に解析結果が依存しにくいという特徴があり、選択した座標系に依存せず、異なる座標系のデータを比較可能であるという特徴がある。また、TDA Mapperには、小さな変位に対しては不変という特徴がある。例えば、ノイズの影響を受け難いという特徴がある。 The TDA Mapper has a feature that the analysis result hardly depends on the coordinate system of the target data, and has a feature that data in different coordinate systems can be compared without depending on the selected coordinate system. Further, the TDA Mapper has a feature that it is invariable with respect to a small displacement. For example, it has a feature that it is hardly affected by noise.
また、TDA Mapperには、多次元のデータを低次元に圧縮して、そのデータの形状を表現可能であるという特徴がある。例えば、数千点のデータのセットを入力した場合に、13の頂点と12のエッジを持つネットワークとして表現できる。 Further, the TDA Mapper has a feature that multidimensional data can be compressed into a low dimension and the shape of the data can be expressed. For example, when a set of thousands of data is input, it can be represented as a network having 13 vertices and 12 edges.
TDA Mapperでは、多次元のデータを入力し、入力したデータの距離メトリックを設定する。次に、距離を設定したデータを低次元に写像するためのフィルター関数を設定する。フィルター関数は1つまたは複数設定可能である。例えば、ユーグリッド距離をメトリックとしたデータをY軸に投影する関数をフィルターとして用いることができる。次に、投影されたデータについて、インターバルを決めてビンに分けていく。インターバル(間隔)とインターバル同士の重なり率(オーバーラップ)が、後述するノードとエッジの関係において重要となる。次に、インターバル毎に、その中に含まれているデータをクラスタリングしていく。クラスタリングのアルゴリズムとしては、例えば、k−means法がある。この場合は、クラスタ数kもパラメータとなる。インターバル毎のクラスタが、可視化した際のノードとなる。クラスタの中のデータが重複している場合、ノード同士の間に線を引く。この線をエッジと呼び、ノード間に繋がりがあることを意味している。 In TDA Mapper, multidimensional data is input, and a distance metric of the input data is set. Next, a filter function for mapping the data for which the distance has been set to a low dimension is set. One or more filter functions can be set. For example, a function that projects data using the U-grid distance as a metric on the Y axis can be used as a filter. Next, the projected data is divided into bins at determined intervals. The interval (interval) and the overlap rate (overlap) between the intervals are important in the relationship between the node and the edge described later. Next, for each interval, the data contained therein is clustered. An example of the clustering algorithm is the k-means method. In this case, the number k of clusters is also a parameter. The cluster at each interval becomes a node at the time of visualization. If data in the cluster is duplicated, draw a line between the nodes. This line is called an edge, and means that there is a connection between the nodes.
以上のような処理の過程で、使用者が設定可能な(設定する必要がある)Mapperのハイパーパラメータとしては、例えば、以下のものがある。
・メトリック:距離の定義。
・フィルター関数:密度、離心率、グラフラプラシアン、各軸への投影、SVD(特異値分解、singular value decomposition)、Isolation forest。
・インターバル:インターバル数を大きくすると、クラスタリング処理した際のクラスタ数は大きくなる。ただし、空のクラスタが多く作られる可能性がある(クラスタあたりのデータの数が少ない)。
・オーバーラップ:オーバーラップ数を大きくすると、ノード間の接続(エッジ)が増える。
・クラスタリング手法:階層型クラスタリング、k−means、その他の手法。
In the course of the above processing, the following hyper-parameters of Mapper that can be set (need to be set) by the user are as follows.
Metric: Definition of distance.
Filter functions: density, eccentricity, graph Laplacian, projection on each axis, SVD (singular value decomposition), Isolation forest.
Interval: If the number of intervals is increased, the number of clusters when performing the clustering process is increased. However, many empty clusters may be created (the number of data per cluster is small).
Overlap: When the number of overlaps is increased, connections (edges) between nodes increase.
Clustering method: Hierarchical clustering, k-means, and other methods.
パラメータチューニング部19と結果整理部20では、これらの制御パラメータと分析結果との相関を整理する(図11参照)。そして、注目しているデータ点がグラフのノードにプロットされるまでパラメータチューニング部19で分析パラメータをスイープする。
The
このパラメータチューニング操作によって、注目しているデータ点が、必ずノード上に表現されるようになる。次に、注目しているデータ点を含むノードが、別の測定日のデータ点を含むか、また別のノードとの関係があるか、すなわちエッジで結ばれているかを確認する。この際、別の測定日のデータが正常のラベルLを有していれば、注目しているデータも正常の可能性が高いと判定される。一方で正常のデータとの繋がりがない場合は、注目したデータは正常でない(異常の可能性が高い)と判定される。 This parameter tuning operation ensures that the data point of interest is represented on the node. Next, it is checked whether the node including the data point of interest includes a data point of another measurement date, and has a relationship with another node, that is, whether the node is connected by an edge. At this time, if the data on another measurement day has a normal label L, it is determined that the data of interest is also likely to be normal. On the other hand, when there is no connection with the normal data, it is determined that the focused data is not normal (the possibility of abnormality is high).
以上の手順によって、使用者が分析パラメータの設定をすることなく、異常の有無を自動で判定することが可能となる。 According to the above procedure, it is possible to automatically determine the presence or absence of the abnormality without setting the analysis parameters by the user.
図12に示すように、データ解析部18で生成され、TDA Mapperの結果であるクラスタ表示(グラフ)がディスプレイに表示される。なお、クラスタ表示の制御は、結果表示制御部21のクラスタ表示制御部39にて行われる。
As shown in FIG. 12, a cluster display (graph) generated by the
このクラスタ表示には、ノードNおよびクラスタCが表示されるグラフ表示領域Gと、ノードNに含まれる内訳情報Bが表示する内訳情報表示領域Dと、ノードNの色分けの態様を示す色ラベルコンターRとが設けられている。 The cluster display includes a graph display area G in which the node N and the cluster C are displayed, a breakdown information display area D in which the breakdown information B included in the node N is displayed, and a color label contour indicating the color classification of the node N. R is provided.
図12の例では、グラフ表示領域Gに、ノードNが丸印で表示されている。例えば、1つのノードNに含まれている対象データまたは過去データの数が多いほど、大きい丸印として表示される。また、複数のノードNがエッジEで繋がれている。このエッジEで繋がれたノード群がクラスタCである。そして、正常なデータを含むノードNのクラスタC1と、異常なデータを含むノードNのクラスタC2,C3,C4とが分けられて表示される。 In the example of FIG. 12, a node N is displayed as a circle in the graph display area G. For example, the larger the number of target data or past data included in one node N, the larger the displayed circle. Further, a plurality of nodes N are connected by an edge E. The node group connected by the edge E is a cluster C. Then, the cluster C1 of the node N containing normal data and the clusters C2, C3, C4 of the node N containing abnormal data are displayed separately.
例えば、正常なデータを含むクラスタC1が1つ表示され、異常なデータを含むクラスタC2が3つ表示されている。なお、本実施形態のクラスタCとは、原則、複数のノードNがエッジEで繋がれたものを示すが、単独のノードNであっても他のクラスタCと離間されている場合には、1つのクラスタCとして扱うものとする。 For example, one cluster C1 including normal data is displayed, and three clusters C2 including abnormal data are displayed. Note that the cluster C of the present embodiment indicates, in principle, a plurality of nodes N connected by an edge E. However, even if a single node N is separated from another cluster C, It is assumed that it is handled as one cluster C.
クラスタ表示は、対象データおよび過去データの双方に基づいて生成されたグラフが表示されている。なお、所定の1つのノードNに対象データおよび過去データの両方が含まれている場合がある。 In the cluster display, a graph generated based on both the target data and the past data is displayed. Note that there is a case where both the target data and the past data are included in one predetermined node N.
使用者は、グラフ表示領域Gに表示された所定のノードNをマウスカーソルMで選択することで、選択されたノードNに含まれる内訳情報Bを、内訳情報表示領域Dに表示させることができる。 The user can cause the breakdown information B included in the selected node N to be displayed in the breakdown information display area D by selecting a predetermined node N displayed in the graph display area G with the mouse cursor M. .
例えば、選択されたノードNに複数の対象データが含まれる場合には、それぞれの対象データが対応する日付(測定日、データ収集日時)が表示される。なお、選択されたノードNに含まれる対象データのデータ数、データラベルなどの属性情報を、内訳情報表示領域Dに表示させても良い。また、ノードNに含まれる統計量を、内訳情報表示領域Dに表示させても良い。 For example, when a plurality of target data are included in the selected node N, a date (measurement date, data collection date and time) corresponding to each target data is displayed. Note that attribute information such as the number of target data items and data labels included in the selected node N may be displayed in the breakdown information display area D. Also, the statistics included in the node N may be displayed in the breakdown information display area D.
なお、マウスカーソルMで選択するものは、クラスタCでも良い。その場合には、選択されたクラスタCに含まれる内訳情報Bを、内訳情報表示領域Dに表示させることができる。 In addition, what is selected by the mouse cursor M may be the cluster C. In that case, the breakdown information B included in the selected cluster C can be displayed in the breakdown information display area D.
色ラベルコンターRは、ノードNの表示色の設定内容が表示される。例えば、正常を示す「A」のラベルLが登録された過去データを含むノードNを黄色で表示し、異常を示す「B」のラベルLが登録された過去データを含むノードNを青色で表示する。このようにすれば、過去データの解析に基づいて、対応付けられたラベルLが識別可能な態様で、過去データを表示することができる。 In the color label contour R, the setting content of the display color of the node N is displayed. For example, a node N including past data in which a label L of “A” indicating normal is registered is displayed in yellow, and a node N including past data in which a label L of “B” indicating abnormal is registered is displayed in blue. I do. In this way, based on the analysis of the past data, the past data can be displayed in such a manner that the associated label L can be identified.
なお、使用者が確認を希望するノードNごとに色分けをしても良いし、それぞれのノードNに対応する周波数で色分けをしても良い。また、例えば、回転機2の分解点検が行われたものと、分解点検が行われていないものというラベルをデータに付けておき、色分けをしても良い。また、解析時に影響を受けたパラメータごとにノードNの色分けしても良い。
Note that the user may perform color coding for each node N for which confirmation is desired, or may perform color coding at a frequency corresponding to each node N. Further, for example, a label indicating that the
また、クラスタ表示では、注目されるデータが異常の可能性が高いと判定された場合は、特異ノード強調表示機能で、異常を有する対象データを含むノードNをGUI上で強調表示(異常可能性表示)させることができる。 Further, in the cluster display, when it is determined that the data of interest has a high possibility of abnormality, the node N including the target data having the abnormality is highlighted on the GUI by the unique node highlighting function (the possibility of abnormality is displayed). Display).
図13に示すように、データ解析部18で生成され、TDA Mapperの結果である時系列表示がディスプレイに表示される。なお、時系列表示の制御は、結果表示制御部21の時系列表示制御部40にて行われる。
As shown in FIG. 13, a time-series display generated by the
この時系列表示では、クラスタ表示にてグラフ化されていたデータが時系列順に表示される。例えば、時系列表示には、モニタリングされた日付(測定日、データ収集日時)の列K1と、それぞれ日付に対応するノードNの列K2と、それぞれノードNに対応するクラスタCの列K3と、それぞれのノードNまたはクラスタCの異常可能性を示す列K4とが設けられている。 In this time-series display, data that has been graphed in the cluster display is displayed in chronological order. For example, the time series display includes a column K1 of monitored dates (measurement date, data collection date and time), a column K2 of a node N corresponding to each date, a column K3 of a cluster C corresponding to each node N, A column K4 indicating the possibility of abnormality of each node N or cluster C is provided.
使用者は、この時系列表示を参照することで、それぞれのノードNまたはクラスタCの異常可能性の有無を表示することができる。なお、例えば、正常を示す「A」のラベルLが登録された過去データを含むノードNに対応して「正常」を表示し、異常を示す「B」のラベルLが登録された過去データを含むノードNに対応して「異常」を表示する。 By referring to the time-series display, the user can display whether or not each node N or cluster C has a possibility of abnormality. For example, “normal” is displayed corresponding to the node N including the past data in which the label “A” indicating normal is registered, and the past data in which the label “B” indicating abnormal is registered is displayed. "Abnormal" is displayed corresponding to the included node N.
つまり、結果表示制御部21は、対象データの解析に基づいて、異常の有無を識別可能な態様(例えば、時系列表示)で対象データを表示する。さらに、過去データの解析に基づいて、対応付けられたラベルLが識別可能な態様で過去データを表示する。このようにすれば、過去データの表示態様を参照して、対象データの異常の有無を識別することができる。
That is, based on the analysis of the target data, the result
なお、クラスタ表示(図12参照)と時系列表示(図13参照)とは、使用者の操作により任意に切り換えることができる。 Note that the cluster display (see FIG. 12) and the time-series display (see FIG. 13) can be arbitrarily switched by a user operation.
異常原因分析部23では、異常データの特徴量抽出機能を用いて、異常可能性があると判定されたデータの特徴量を抽出する。
The abnormality
ここで、分析対象と分析したデータを周波数成分データとした場合について説明する。例えば、異常と判定されたデータと、正常ラベルがついたデータに関して、各周波数成分の差を算出し、各周波数成分の差が大きい周波数を抽出する。この抽出した周波数と周波数成分の差を異常データの特徴量とする。差の大きい周波数成分の定義は、予め設定した閾値を超えるもの、または、上位の数%の差を抽出するなどの任意の抽出方法を選択することができる。そして、抽出した特徴量(周波数、周波数成分の差)を異常原因分析部23に入力する。
Here, a case will be described in which data to be analyzed is frequency component data. For example, a difference between each frequency component is calculated for the data determined to be abnormal and the data with a normal label, and a frequency having a large difference between the frequency components is extracted. The difference between the extracted frequency and the frequency component is used as the feature amount of the abnormal data. For the definition of the frequency component having a large difference, an extraction method that exceeds a preset threshold value or an arbitrary extraction method such as extraction of a difference of several% at the top can be selected. Then, the extracted feature amounts (the difference between the frequency and the frequency component) are input to the abnormality
本実施形態では、異常原因分析部23に特徴量が入力されることで、異常原因への特徴量の影響度が分かり、異常原因を推定することが可能となる。また、異常検知に必要なパラメータチューニングを自動化することで、異常検知を人手によるチューニングを行わないで、異常原因を推定することが可能となる。また、異常検知の結果から異常の原因を推定することが可能となる。
In the present embodiment, by inputting the characteristic amount to the abnormality
次に、回転機2の異常診断方法について図14のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。
Next, a method of diagnosing abnormality of the
図14に示すように、まず、ステップS11において、異常診断システム1の使用者は、回転機2に点検用センサ3を取り付ける。
As shown in FIG. 14, first, in step S11, the user of the
次のステップS12において、使用者は、回転機2に取り付けた点検用センサ3を用いて、回転機2の運転中の状態をモニタリングしてモニタリングデータを得る。
In the next step S12, the user monitors the operating state of the
次のステップS13において、使用者は、点検用センサ3を用いて取得したモニタリングデータを、ユーザ端末5を用いて異常診断システム1のモニタリングデータベース11に入力する。
In the next step S13, the user inputs the monitoring data acquired using the
次のステップS14において、使用者は、診断対象となる回転機2のプロセスデータを、ユーザ端末5を用いて異常診断システム1のプロセスデータベース12に入力する。
In the next step S <b> 14, the user inputs the process data of the
次のステップS15において、使用者は、診断対象となる回転機2のスペックデータまたはプロセスデータを、ユーザ端末5を用いて異常診断システム1のスペックデータベース13に入力する。
In the next step S15, the user inputs the specification data or process data of the
次のステップS16において、使用者は、診断対象となる回転機2の異常原因データを、ユーザ端末5を用いて異常診断システム1の異常原因データベース14に入力する。
In the next step S16, the user inputs the abnormality cause data of the
次のステップS17において、使用者は、異常原因データの設定であって、解析処理ユニット10の解析に基づいて抽出され得る対象データの特徴量と、回転機2の異常原因との関係性の設定を行う。つまり、異常原因データベース14には、抽出され得る対象データの特徴量と異常原因との対応関係が予め記憶される。
In the next step S17, the user sets the abnormality cause data, that is, sets the relationship between the characteristic amount of the target data that can be extracted based on the analysis of the
次のステップS18において、使用者は、異常診断システム1を用いて、モニタリングデータに含まれる対象データまたは過去データの解析を行う。
In the next step S <b> 18, the user uses the
次のステップS19において、使用者は、解析結果を異常診断システム1のディスプレイに解析結果を表示させる所定の操作を行う。使用者は、例えば、クラスタ表示の画面を表示させる操作、時系列表示の画面を表示させる操作、ノードNまたはクラスタCの選択をする操作、異常原因分析結果を示す画面を表示させる操作を行う。
In the next step S19, the user performs a predetermined operation of displaying the analysis result on the display of the
次のステップS20において、使用者は、異常診断システム1のディスプレイに表示された情報に基づいて、回転機2の異常の有無の識別を行う。
In the next step S20, the user identifies the presence or absence of an abnormality in the
次のステップS21において、使用者は、異常診断システム1のディスプレイに表示された情報に基づいて、回転機2の異常原因の特定を行う。そして、異常診断方法を終了する。
In the next step S21, the user specifies the cause of the abnormality of the
次に、解析処理ユニット10が実行する異常診断処理について図15のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。
Next, the abnormality diagnosis processing executed by the
図15に示すように、まず、ステップS31において、解析処理ユニット10のデータ取得部15は、モニタリングデータベース11に記憶されたモニタリングデータを取得する。ここで、プロセスデータベース12に記憶されたプロセスデータ、およびスペックデータベース13に記憶されたスペックデータの取得を行っても良い。
As shown in FIG. 15, first, in step S31, the
次のステップS32において、解析処理ユニット10の処理設定部17の前処理条件設定部34は、データ前処理部16にて実行される前処理条件の設定を行う。ここで、パラメータチューニング部19の特定次元数設定部36は、前処理条件設定部34の設定に基づいて、対象データまたは過去データの特定次元数の設定を行う。また、処理設定部17のアルゴリズム設定部35は、データ解析部18で用いられるアルゴリズムの設定を行う。アルゴリズムとしては、例えば、TDA Mapperを用いる。
In the next step S32, the preprocessing
次のステップS33において、解析処理ユニット10のデータ前処理部16は、データ前処理を実行する。
In the next step S33, the
次のステップS34において、解析処理ユニット10のデータ解析部18は、データ解析処理を実行する。ここで、特定次元数に調整された対象データを、特定次元数のデータを扱うアルゴリズムを用いて解析する。さらに、特定次元数に調整された過去データを、アルゴリズムを用いて解析する。
In the next step S34, the
なお、アルゴリズムとしてTDA Mapperを用いる場合には、データ解析部18のデータ解析処理において、データ入力処理とマッピング処理とカバー化処理とクラスタリング・ノード化処理とグラフ化処理が実行される。
When TDA Mapper is used as an algorithm, data input processing, mapping processing, cover processing, clustering / node conversion processing, and graph processing are executed in the data analysis processing of the
次のステップS35において、解析処理ユニット10の結果整理部20は、解析結果整理処理を実行する。ここで、データ解析部18により得られたデータの解析結果が最適なものであるか否かを整理する。
In the next step S35, the
次のステップS36において、結果整理部20は、解析結果が最適であるか否かを判定する。なお、最適であるか否かの判定は、例えば、対象データの表示が異常の有無を識別可能な態様になるか否かに基づいて設定される。ここで、解析結果が最適である場合(ステップS36がYES)は、後述のステップS38に進む。一方、解析結果が最適でない場合(ステップS36がNO)は、ステップS37に進む。
In the next step S36, the
ステップS37において、解析処理ユニット10のパラメータチューニング部19は、パラメータチューニング処理を実行する。ここで、データ解析部18で解析を行うときに用いる各種設定または解析用のアルゴリズムの設定のチューニングが行われる。そして、前述のステップS33に戻る。つまり、本実施形態では、対象データの表示が異常の有無を識別可能な態様になるように特定次元数の設定を繰り返す。
In step S37, the
ステップS38において、解析処理ユニット10の異常原因分析部23は、異常原因分析処理を実行する。ここで、特徴量抽出部43は、データ解析部18で解析された後の対象データの特徴量を抽出する。また、異常原因推定部44は、特徴量抽出部43にて抽出された特徴量に基づいて、異常原因を推定し、その異常原因を診断可能な態様で対象データを表示可能な状態になるようにデータを整理する。
In step S38, the abnormality
次のステップS39において、解析処理ユニット10の表示調整部22は、表示調整処理を実行する。ここで、表示調整部22は、使用者の操作に基づいて、ディスプレイに表示される画面の制御を行う。そして、異常診断処理を終了する。
In the next step S39, the
次に、データ前処理部16が実行するデータ前処理について図16のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。
Next, data pre-processing performed by the
図16に示すように、まず、ステップS41において、データ前処理部16の対象データ調整部26は、モニタリングデータに含まれる診断の対象となる対象データの次元数が特定次元数になるように対象データを調整する処理を行う。
As shown in FIG. 16, first, in step S41, the target
次のステップS42において、データ前処理部16の過去データ調整部27は、モニタリングデータに含まれる異常の有無が既知の過去データの次元数が特定次元数になるように過去データを調整する処理を行う。
In the next step S42, the past
次のステップS43において、データ前処理部16の欠損情報補間部28は、対象データを特定次元数にしたときに、それぞれの次元数の情報のうち一部の次元数の情報が欠損する場合に、その欠損した情報を補間する処理を行う。
In the next step S43, the missing
次のステップS44において、データ前処理部16のデータラベル付与部29は、特定次元数に調整された過去データに、異常の有無を識別可能なラベルLを対応付ける処理を行う。
In the next step S44, the data
次のステップS45において、データ前処理部16は、その他の処理を実行する。例えば、データスケーリング部30がデータを規格化する。スペクトル解析部31がスペクトル成分を算出する。振幅値算出部32が振幅値を算出する。データフォーマット変換部33がデータフォーマットを変換する。そして、データ前処理を終了する。
In the next step S45, the
次に、表示調整部22が実行する表示調整処理について図17のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。
Next, the display adjustment processing executed by the
図17に示すように、まず、ステップS51において、表示調整部22は、ディスプレイにクラスタ表示(グラフ)の画面が表示中であるか否かを判定する。ここで、グラフが表示中でない場合(ステップS51がNO)は、後述のステップS54に進む。一方、グラフが表示中である場合(ステップS51がYES)は、ステップS52に進む。
As shown in FIG. 17, first, in step S51, the
ステップS52において、表示調整部22のノード選択受付部41は、クラスタ表示された画面上のノードNまたはクラスタCの選択を受け付けたか否かを判定する。ここで、ノードNまたはクラスタCの選択を受け付けていない場合(ステップS52がNO)は、ステップS54に進む。一方、ノードNまたはクラスタCの選択を受け付けた場合(ステップS52がYES)は、ステップS53に進む。
In step S52, the node
ステップS53において、表示調整部22の対象情報表示部42は、選択されたノードNまたはクラスタCに含まれる対象データまたは過去データの内訳情報B(対象情報)を内訳情報表示領域Dに表示する処理を行う。そして、ステップS54に進む。
In step S53, the target
ステップS54において、表示調整部22は、時系列表示の画面をディスプレイに表示させる操作を受け付けたか否かを判定する。ここで、時系列表示の表示操作の受け付けがない場合(ステップS54がNO)は、後述のステップS56に進む。一方、時系列表示の表示操作の受け付けがある場合(ステップS54がYES)は、ステップS55に進む。
In step S54, the
ステップS55において、表示調整部22の処理に基づいて、結果表示制御部21の時系列表示制御部40が、時系列表示の画面をディスプレイに表示させる。この時系列表示制御部40は、対象データをモニタリングされた時系列順に表示し、かつ対象データが含まれるノードを識別可能な態様で表示する。なお、時系列表示制御部40は、過去データの解析に基づいて、対応付けられたラベルLが識別可能な態様で過去データを表示する。例えば、異常可能性を示す表示を行う(図13参照)。そして、ステップS56に進む。
In step S55, the time-series
ステップS56において、表示調整部22は、クラスタ表示の画面をディスプレイに表示させる操作を受け付けたか否かを判定する。ここで、クラスタ表示の表示操作の受け付けがない場合(ステップS56がNO)は、後述のステップS58に進む。一方、クラスタ表示の表示操作の受け付けがある場合(ステップS56がYES)は、ステップS57に進む。
In step S56, the
ステップS57において、表示調整部22の処理に基づいて、結果表示制御部21のクラスタ表示制御部39が、クラスタ表示の画面をディスプレイに表示させる。このクラスタ表示制御部39は、対象データの解析に基づいて、異常の有無を識別可能な態様で対象データを表示する。なお、クラスタ表示制御部39は、過去データの解析に基づいて、対応付けられたラベルLが識別可能な態様で過去データを表示する。例えば、異常可能性に応じてノードNを色分けして表示する(図12参照)。そして、ステップS58に進む。
In step S57, based on the processing of the
ステップS58において、表示調整部22は、異常原因分析結果を示す画面をディスプレイに表示させる操作を受け付けたか否かを判定する。ここで、異常原因分析結果の表示操作の受け付けがない場合(ステップS58がNO)は、後述のステップS60に進む。一方、異常原因分析結果の表示操作の受け付けがある場合(ステップS58がYES)は、ステップS59に進む。
In step S58, the
ステップS59において、表示調整部22の処理に基づいて、原因表示制御部25が、異常原因分析結果を示す画面をディスプレイに表示させる。そして、ステップS60に進む。
In step S59, based on the processing of the
ステップS60において、表示調整部22は、ディスプレイに表示中の解析結果を示す画面の表示を終了させる操作を受け付けたか否かを判定する。ここで、表示を終了させる操作の受け付けがない場合(ステップS60がNO)は、前述のステップS51に進む。一方、表示を終了させる操作の受け付けがある場合(ステップS60がYES)は、ステップS61に進む。
In step S60, the
ステップS61において、表示調整部22の処理に基づいて、結果表示制御部21または原因表示制御部25は、ディスプレイに表示中の解析結果を示す画面の表示を終了する。そして、表示調整処理を終了する。
In step S61, based on the processing of the
なお、本実施形態のフローチャートにおいて、各ステップが直列に実行される形態を例示しているが、必ずしも各ステップの前後関係が固定されるものでなく、一部のステップの前後関係が入れ替わっても良い。また、一部のステップが他のステップと並列に実行されても良い。 Note that, in the flowchart of the present embodiment, an example in which each step is executed in series is illustrated, but the order of each step is not necessarily fixed, and even if the order of some steps is switched. good. Also, some steps may be executed in parallel with other steps.
本実施形態のシステムは、専用のチップ、FPGA(Field Programmable Gate Array)、GPU(Graphics Processing Unit)、またはCPU(Central Processing Unit)などのプロセッサを高集積化させた制御装置と、ROM(Read Only Memory)またはRAM(Random Access Memory)などの記憶装置と、HDD(Hard Disk Drive)またはSSD(Solid State Drive)などの外部記憶装置と、ディスプレイなどの表示装置と、マウスまたはキーボードなどの入力装置と、通信インターフェースとを備える。このシステムは、通常のコンピュータを利用したハードウェア構成で実現できる。 The system of the present embodiment includes a control device in which processors such as a dedicated chip, an FPGA (Field Programmable Gate Array), a GPU (Graphics Processing Unit), or a CPU (Central Processing Unit) are highly integrated, and a ROM (Read Only). Memory) or a storage device such as a random access memory (RAM), an external storage device such as a hard disk drive (HDD) or a solid state drive (SSD), a display device such as a display, and an input device such as a mouse or a keyboard. , A communication interface. This system can be realized by a hardware configuration using an ordinary computer.
なお、本実施形態のシステムで実行されるプログラムは、ROMなどに予め組み込んで提供される。もしくは、このプログラムは、インストール可能な形式または実行可能な形式のファイルでCD−ROM、CD−R、メモリカード、DVD、フレキシブルディスク(FD)などのコンピュータで読み取り可能な非一過性の記憶媒体に記憶されて提供するようにしても良い。 The program executed by the system according to the present embodiment is provided by being incorporated in a ROM or the like in advance. Alternatively, this program is a non-transitory computer-readable storage medium such as a CD-ROM, a CD-R, a memory card, a DVD, and a flexible disk (FD) in a file in an installable or executable format. And may be provided.
また、このシステムで実行されるプログラムは、インターネットなどのネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせて提供するようにしても良い。また、このシステムは、構成要素の各機能を独立して発揮する別々のモジュールを、ネットワークまたは専用線で相互に接続し、組み合わせて構成することもできる。 Further, the program executed by this system may be stored on a computer connected to a network such as the Internet, and downloaded and provided via the network. In addition, the system can be configured by combining separate modules that independently perform the functions of the components by mutually connecting them via a network or a dedicated line.
なお、本実施形態では、異常の有無を識別可能な態様として、位相的データ解析を用いたクラスタ表示を例示しているが、その他の態様であっても良い。例えば、異常の有無を識別可能な態様としては、階層型のクラスタ解析に基づく階層型のクラスタ表示であっても良いし、非階層型のクラスタ解析に基づく非階層型のクラスタ表示であっても良い。非階層型のクラスタ解析では、例えば、K−means法などの手法を用いても良い。 Note that, in the present embodiment, cluster display using topological data analysis is exemplified as a mode in which the presence or absence of an abnormality can be identified, but other modes may be used. For example, a mode in which the presence or absence of an abnormality can be identified may be a hierarchical cluster display based on a hierarchical cluster analysis, or a non-hierarchical cluster display based on a non-hierarchical cluster analysis. good. In the non-hierarchical cluster analysis, for example, a technique such as the K-means method may be used.
以下、図18、図19、図20を用いて階層型クラスタリングを用いた場合の変形例を説明する。 Hereinafter, a modified example using hierarchical clustering will be described with reference to FIGS. 18, 19, and 20.
例えば、階層型クラスタリングがクラスタリングされていないN個の入力データについて、類似度の高い順に融合して次第に大きなクラスタを作る。そして、最終的にはN個データを1つのクラスタに統合していく手法で、データまたはクラスタを統合していく過程を、樹状図と呼ばれる木の形で表現する。 For example, N input data items that are not subjected to hierarchical clustering are merged in descending order of similarity to form a gradually larger cluster. Finally, the process of integrating data or clusters is represented by a tree called a dendrogram by a method of integrating N pieces of data into one cluster.
ここで、階層型クラスタリングのパラメータとしては、データの類似度を測る尺度と方法がある類似度を測る尺度は、ユーグリッド距離、ミンコフスキー距離などがある。また、クラスタ間の距離の測定方法としては、ウォード法または群平均法などがある。これらのパラメータは、パラメータチューニング部19で最適化される。
Here, as parameters of the hierarchical clustering, there are a scale and a method for measuring the similarity of data, and a scale for measuring the similarity includes a U-grid distance and a Minkowski distance. In addition, as a method of measuring the distance between clusters, there are a Ward method, a group average method, and the like. These parameters are optimized by the
また、パラメータチューニングを行い最適化する際の指標としては、例えば、クラスタの距離とコーフエン行列との相関係数(コーフエン相関係数)がある。 An index used for parameter tuning and optimization includes, for example, a correlation coefficient between a cluster distance and a Corfuen matrix (Corfuen correlation coefficient).
図18は、階層型クラスタによる表示結果を示している。例えば、データ前処理部16で処理した分析対象データを、データ解析部18に入力し、階層型クラスタリング手法を用いて分析し、分析結果を表示出力部9で表示させる。その結果のうち、クラスタ表示した結果である。ここでは、k次元からなる57個の振幅値データを入力している。57個の振幅値データのうち、正常の範囲にあるものと、異常を示すものとがある。
FIG. 18 shows a display result by a hierarchical cluster. For example, the analysis target data processed by the
例えば、図18において、個々の個別のデータをノードNとすると、複数のノードN、または、少なくとも1つのノードNが含まれるノード群をクラスタCと定義することができる。この階層型クラスタリングを用いた場合も、前述のTDAと同様に、正常なデータを含むノードNのクラスタC5,C6と、異常なデータを含むノードNのクラスタC7とに分けて表示させることが可能である。 For example, in FIG. 18, when each individual data is a node N, a plurality of nodes N or a node group including at least one node N can be defined as a cluster C. Also in the case of using this hierarchical clustering, similarly to the above-described TDA, it is possible to separately display the clusters C5 and C6 of the node N including normal data and the cluster C7 of the node N including abnormal data. It is.
図19に示すように、階層型クラスタリングでは、各種データ、クラスタが融合されていく過程においては、正常データは、平均的な分岐の態様となる。これに対して図20に示すように、異常データは、上層部において分岐されて表示される。なお、前述の説明は、正常データ数が異常データ数に対して多い場合を例示しているが、正常データ数が異常データ数に対して少ない場合も有り得る。 As shown in FIG. 19, in the hierarchical clustering, in the process of fusing various data and clusters, normal data has an average branching mode. On the other hand, as shown in FIG. 20, the abnormal data is branched and displayed in the upper layer. In the above description, the case where the number of normal data is larger than the number of abnormal data is exemplified. However, the number of normal data may be smaller than the number of abnormal data.
以上説明した実施形態によれば、特定次元数に調整された対象データを、特定次元数のデータを扱うアルゴリズムを用いて解析するステップを含むことにより、エンジニアの技能に依存せずに診断を行うことができる。 According to the embodiment described above, diagnosis is performed without depending on the skill of an engineer by including a step of analyzing target data adjusted to a specific number of dimensions using an algorithm that handles data of a specific number of dimensions. be able to.
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。 Although several embodiments of the present invention have been described, these embodiments are provided by way of example and are not intended to limit the scope of the invention. These embodiments can be implemented in other various forms, and various omissions, replacements, changes, and combinations can be made without departing from the gist of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are also included in the invention described in the claims and equivalents thereof.
1…異常診断システム、2…回転機、3…点検用センサ、4…インターネット、5…ユーザ端末、6…メイン制御部、7…情報入力部、8…通信部、9…表示出力部、10…解析処理ユニット、11…モニタリングデータベース、12…プロセスデータベース、13…スペックデータベース、14…異常原因データベース、15…データ取得部、16…データ前処理部、17…処理設定部、18…データ解析部、19…パラメータチューニング部、20…結果整理部、21…結果表示制御部、22…表示調整部、23…異常原因分析部、24…原因情報取得部、25…原因表示制御部、26…対象データ調整部、27…過去データ調整部、28…欠損情報補間部、29…データラベル付与部、30…データスケーリング部、31…スペクトル解析部、32…振幅値算出部、33…データフォーマット変換部、34…前処理条件設定部、35…アルゴリズム設定部、36…特定次元数設定部、38…パラメータ最適化部、39…クラスタ表示制御部、40…時系列表示制御部、41…ノード選択受付部、42…対象情報表示部、43…特徴量抽出部、44…異常原因推定部、B…内訳情報、C…クラスタ、D…内訳情報表示領域、E…エッジ、F…周波数の列、G…グラフ表示領域、K1〜K4…時系列表示の列、L…ラベル、M…マウスカーソル、N…ノード、P…過去データ、R…色ラベルコンター、S…加速度振幅の列、T…対象データ。
DESCRIPTION OF
Claims (14)
前記特定次元数に調整された前記対象データを、前記特定次元数のデータを扱うアルゴリズムであるマッパーを用いて位相的データ解析を行うステップと、
前記マッパーを用いた解析結果を異常の有無を識別可能な態様であるクラスタ表示または時系列表示で表示するステップと、
を含む、
機器の異常診断方法。 Individual data extracted from at least one of the time waveform and the frequency spectrum of a signal that is data indicating the state of the device and included in one target data to be diagnosed among the data indicating the state of the device a step of adjustment to a specific number of dimensions the number of dimensions is the number of,
The target data adjusted to the specific number of dimensions , performing a topological data analysis using a mapper that is an algorithm that handles the data of the specific number of dimensions,
A step of displaying the analysis result using the mapper in a cluster display or a time-series display in a mode in which the presence or absence of an abnormality can be identified,
including,
Device abnormality diagnosis method.
前記異常の有無を識別可能な態様は、前記クラスタ表示としての3次元以下のグラフである、
請求項1に記載の機器の異常診断方法。 The specific dimension number is four or more dimensions,
An aspect in which the presence or absence of the abnormality can be identified is a three-dimensional or less graph as the cluster display .
The method for diagnosing abnormality of a device according to claim 1.
前記特定次元数に調整された前記過去データに、前記異常の有無を識別可能なラベルを対応付けるステップと、
前記特定次元数に調整された前記過去データを、前記マッパーを用いて解析するステップと、
対応付けられた前記ラベルが識別可能な態様で前記マッパーを用いた解析結果を表示するステップと、
請求項1または請求項2に記載の機器の異常診断方法。 Adjusting the past data so that the number of dimensions of the known past data, the presence or absence of the abnormality included in the data indicating the state of the device, is the specific dimension number,
A step of associating the past data adjusted to the specific number of dimensions with a label capable of identifying the presence or absence of the abnormality,
Analyzing the past data adjusted to the specific number of dimensions using the mapper ,
And displaying the analysis results the labels applied corresponds is using the mapper in an identifiable manner,
The method for diagnosing a device abnormality according to claim 1 or 2.
請求項1から請求項3のいずれか1項に記載の機器の異常診断方法。 Including setting the number of specific dimensions,
The method for diagnosing an abnormality of an apparatus according to claim 1.
請求項4に記載の機器の異常診断方法。 Including a step of repeating the setting of the specific number of dimensions so that the display of the target data is in a form in which the presence or absence of the abnormality can be identified.
An apparatus abnormality diagnosis method according to claim 4.
請求項1から請求項5のいずれか1項に記載の機器の異常診断方法。 When the target data is the specific number of dimensions, when information of some of the number of dimensions of the information of each dimension is missing, including the step of interpolating the missing information,
The method for diagnosing an abnormality of a device according to any one of claims 1 to 5.
抽出された前記特徴量に基づいて、異常原因を診断可能な態様で前記対象データを表示するステップと、
を含む、
請求項1から請求項6のいずれか1項に記載の機器の異常診断方法。 Extracting a feature amount of the target data;
Displaying the target data in a mode in which the cause of the abnormality can be diagnosed based on the extracted feature amount;
including,
The method for diagnosing abnormality of a device according to any one of claims 1 to 6.
請求項7に記載の機器の異常診断方法。 Including a step of storing in advance a correspondence relationship between the feature amount and the cause of the abnormality,
A method for diagnosing abnormality of a device according to claim 7.
請求項1から請求項8のいずれか1項に記載の機器の異常診断方法。 An aspect in which the presence or absence of the abnormality can be identified is a cluster display in which the target data is a node.
The method for diagnosing an abnormality of a device according to any one of claims 1 to 8 .
選択された前記ノードに含まれる前記対象データの情報を表示するステップと、
を含む、
請求項9に記載の機器の異常診断方法。 Receiving a selection of the node displayed in the cluster,
Displaying information of the target data included in the selected node;
including,
The method for diagnosing abnormality of a device according to claim 9 .
請求項9または請求項10に記載の機器の異常診断方法。 Displaying the target data in a monitored chronological order, and including displaying the node including the target data in an identifiable manner,
The method for diagnosing a device abnormality according to claim 9 or claim 10 .
請求項1から請求項11のいずれか1項に記載の機器の異常診断方法。 Data indicating the state of the device is data indicating vibration during operation of the device,
An abnormality diagnosis method for a device according to any one of claims 1 to 11 .
前記特定次元数に調整された前記対象データを、前記特定次元数のデータを扱うアルゴリズムであるマッパーを用いて位相的データ解析を行うデータ解析部と、
前記マッパーを用いた解析結果を異常の有無を識別可能な態様であるクラスタ表示または時系列表示で表示する結果表示制御部と、
を備える、
機器の異常診断システム。 Individual data extracted from at least one of the time waveform and the frequency spectrum of a signal that is data indicating the state of the device and included in one target data to be diagnosed among the data indicating the state of the device a data preprocessing unit for adjustment to a particular number of dimensions the number of dimensions is the number of,
The target data adjusted to the specific dimension number, a data analysis unit that performs topological data analysis using a mapper that is an algorithm that handles the data of the specific dimension number,
A result display control unit that displays the analysis result using the mapper in a cluster display or a time-series display that is a mode in which the presence or absence of an abnormality can be identified,
Comprising,
Equipment diagnosis system.
前記機器に取り付けた点検用センサを用いて、前記機器の状態をモニタリングしてモニタリングデータを得るステップと、
前記モニタリングデータを前記異常診断システムに入力するステップと、
を含む、
機器の異常診断方法。 A device abnormality diagnosis method using the device abnormality diagnosis system according to claim 13 ,
Using an inspection sensor attached to the device, monitoring the state of the device to obtain monitoring data,
Inputting the monitoring data to the abnormality diagnosis system;
including,
Device abnormality diagnosis method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019131439A JP6641056B1 (en) | 2019-07-16 | 2019-07-16 | Device abnormality diagnosis method and device abnormality diagnosis system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019131439A JP6641056B1 (en) | 2019-07-16 | 2019-07-16 | Device abnormality diagnosis method and device abnormality diagnosis system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP6641056B1 true JP6641056B1 (en) | 2020-02-05 |
JP2021015103A JP2021015103A (en) | 2021-02-12 |
Family
ID=69320965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019131439A Active JP6641056B1 (en) | 2019-07-16 | 2019-07-16 | Device abnormality diagnosis method and device abnormality diagnosis system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6641056B1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116593145B (en) * | 2023-05-16 | 2024-01-09 | 利维智能(深圳)有限公司 | Fault detection method, device, computer equipment and storage medium |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5214577A (en) * | 1990-10-24 | 1993-05-25 | Osaka Gas Co., Ltd. | Automatic test generation for model-based real-time fault diagnostic systems |
JPH05256741A (en) * | 1992-03-11 | 1993-10-05 | Toshiba Corp | Method and apparatus for monitoring plant signal |
JP5439265B2 (en) * | 2010-04-20 | 2014-03-12 | 株式会社日立製作所 | Abnormality detection / diagnosis method, abnormality detection / diagnosis system, and abnormality detection / diagnosis program |
JP6359123B2 (en) * | 2015-01-21 | 2018-07-18 | 三菱電機株式会社 | Inspection data processing apparatus and inspection data processing method |
US10552762B2 (en) * | 2015-07-16 | 2020-02-04 | Falkonry Inc. | Machine learning of physical conditions based on abstract relations and sparse labels |
JP6863383B2 (en) * | 2016-08-31 | 2021-04-21 | 日本電気株式会社 | Defect detectors, defect detection methods, and computer-readable recording media |
JP6828626B2 (en) * | 2017-07-11 | 2021-02-10 | 富士通株式会社 | Anomaly detection program, anomaly detection method, and information processing device |
JP6950504B2 (en) * | 2017-12-08 | 2021-10-13 | 富士通株式会社 | Abnormal candidate extraction program, abnormal candidate extraction method and abnormal candidate extraction device |
-
2019
- 2019-07-16 JP JP2019131439A patent/JP6641056B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2021015103A (en) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8868985B2 (en) | Supervised fault learning using rule-generated samples for machine condition monitoring | |
KR101316486B1 (en) | Error detection method and system | |
KR102141391B1 (en) | Failure data management method based on cluster estimation | |
KR101948604B1 (en) | Method and device for equipment health monitoring based on sensor clustering | |
US20120041575A1 (en) | Anomaly Detection Method and Anomaly Detection System | |
US10496730B2 (en) | Factor analysis device, factor analysis method, and factor analysis program | |
US11138817B2 (en) | Vehicle system prognosis device and method | |
JP6828807B2 (en) | Data analysis device, data analysis method and data analysis program | |
EP2930579A2 (en) | State monitoring system, state monitoring method and state monitoring program | |
EP1958034B1 (en) | Use of sequential clustering for instance selection in machine condition monitoring | |
JP6792746B2 (en) | State identification method by segment feature analysis of frequency domain | |
US12007745B2 (en) | Apparatus for predicting equipment damage | |
KR102253230B1 (en) | Predictive diagnosis method and system of nuclear power plant equipment | |
EP4057192A1 (en) | Integrated analysis method, integrated analysis device, and integrated analysis program | |
CN116629627A (en) | Intelligent detection system of power transmission on-line monitoring device | |
JP2020071624A (en) | Abnormality diagnosing apparatus, abnormality diagnosing method and program | |
JP2022035686A (en) | Diagnostic device and parameter adjustment method | |
CN112200048A (en) | Regression model-based rotating equipment fault prediction method and system and readable storage medium | |
JP6641056B1 (en) | Device abnormality diagnosis method and device abnormality diagnosis system | |
CN117782198B (en) | Highway electromechanical equipment operation monitoring method and system based on cloud edge architecture | |
JPH0793018A (en) | Method and system for diagnosing operating state | |
TWI780434B (en) | Abnormal diagnosis device and method | |
TW202045898A (en) | Analysis apparatus, analysis method, and non-transitory computer readable storage medium | |
KR102196611B1 (en) | Predictive diagnosis method and system of nuclear power plant equipment | |
CN114660931A (en) | Method and system for diagnosing and identifying petrochemical process fault |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190716 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20190716 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20190801 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190910 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191108 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20191203 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20191226 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6641056 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |