JP6641056B1 - Device abnormality diagnosis method and device abnormality diagnosis system - Google Patents

Device abnormality diagnosis method and device abnormality diagnosis system Download PDF

Info

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
Application number
JP2019131439A
Other languages
Japanese (ja)
Other versions
JP2021015103A (en
Inventor
啓 高倉
啓 高倉
尾崎 健司
健司 尾崎
翔一 柏瀬
翔一 柏瀬
優弥 西
優弥 西
日隈 幸治
幸治 日隈
朋美 菱沼
朋美 菱沼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Energy Systems and Solutions Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Energy Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2019131439A priority Critical patent/JP6641056B1/en
Application granted granted Critical
Publication of JP6641056B1 publication Critical patent/JP6641056B1/en
Publication of JP2021015103A publication Critical patent/JP2021015103A/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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.

特開平10−274558号公報JP-A-10-274558 特開2000−65246号公報JP 2000-65246 A 特開2019−67069号公報JP 2019-67069 A

前述の技術では、正常と異常の判定に、基準となる閾値が用いられている。その基準については、機器ごとにエンジニアが基準値を設ける必要がある。この基準値の設定は、熟練のエンジニアが過去の測定データまたは対象機器の仕様を考慮しながら行うため、時間と労力を要する。このため、従来の手法ではエンジニアの技能の経験およびスキルに対する依存性が残っている。   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.

本実施形態の機器の異常診断システムを示す構成図。FIG. 1 is a configuration diagram illustrating a device abnormality diagnosis system according to an embodiment. 解析処理ユニットを示すブロック図。FIG. 2 is a block diagram showing an analysis processing unit. データ前処理部を示すブロック図。FIG. 3 is a block diagram showing a data preprocessing unit. 処理設定部を示すブロック図。FIG. 3 is a block diagram showing a process setting unit. パラメータチューニング部を示すブロック図。FIG. 3 is a block diagram showing a parameter tuning unit. 結果表示制御部を示すブロック図。FIG. 4 is a block diagram showing a result display control unit. 表示調整部を示すブロック図。FIG. 3 is a block diagram illustrating a display adjustment unit. 異常原因分析部を示すブロック図。FIG. 3 is a block diagram illustrating an abnormality cause analysis unit. 前処理後の加速度振幅のデータの一例を示す説明図。Explanatory drawing which shows an example of the data of the acceleration amplitude after preprocessing. 前処理後の周波数のデータの一例を示す説明図。FIG. 4 is an explanatory diagram showing an example of frequency data after preprocessing. パラメータと分析結果との相関関係の一例を示すテーブル。9 is a table illustrating an example of a correlation between a parameter and an analysis result. クラスタ表示の一例を示す画面図。FIG. 4 is a screen diagram showing an example of a cluster display. 時系列表示の一例を示す画面図。FIG. 4 is a screen diagram showing an example of a time-series display. 機器の異常診断方法を示すフローチャート。5 is a flowchart illustrating a method for diagnosing device abnormality. 異常診断処理を示すフローチャート。9 is a flowchart illustrating abnormality diagnosis processing. データ前処理を示すフローチャート。9 is a flowchart showing data preprocessing. 表示調整処理を示すフローチャート。9 is a flowchart illustrating a display adjustment process. 変形例の階層型クラスタリングを用いた解析結果を示す説明図。FIG. 9 is an explanatory diagram showing an analysis result using hierarchical clustering of a modified example. 変形例の階層型クラスタリングを用いた解析結果を示す説明図。FIG. 9 is an explanatory diagram showing an analysis result using hierarchical clustering of a modified example. 変形例の階層型クラスタリングを用いた解析結果を示す説明図。FIG. 9 is an explanatory diagram showing an analysis result using hierarchical clustering of a modified example.

以下、図面を参照しながら、機器の異常診断方法および機器の異常診断システムの実施形態について詳細に説明する。   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とは、例えば、ポンプまたは送風機などに搭載されるモータである。   Reference numeral 1 in FIG. 1 denotes an apparatus abnormality diagnosis system according to the present embodiment. This abnormality diagnosis system 1 is a system that diagnoses an abnormality of a rotating machine 2 as a device provided in a facility such as a power plant. The rotating machine 2 is, for example, a motor mounted on a pump or a blower.

例えば、回転機2に故障などの異常がある場合、または異常の兆候がある場合には、その振動の加速度振幅または特定の周波数成分が、正常時と比較して異なるものとなる。正常時と異常時との差異が大きなものであれば、異常の発見がし易いが、その差異が微小なものであれば、異常の発見が困難である。また、機器としての回転機2の運転中に生じる現象に関するデータは、膨大な量となり、これら大量のデータを人手で解析することは困難である。   For example, when there is an abnormality such as a failure in the rotating machine 2 or when there is a sign of the abnormality, the acceleration amplitude or the specific frequency component of the vibration is different from that in the normal state. If the difference between the normal state and the abnormal state is large, it is easy to find the abnormality. If the difference is small, it is difficult to find the abnormality. In addition, the amount of data relating to the phenomena occurring during the operation of the rotating machine 2 as an apparatus is enormous, and it is difficult to manually analyze such a large amount of data.

近年、機械学習および人工知能に関する技術の発達に伴い、大量のデータから有用な情報を取り出すことが可能になっている。このようなデータ解析のアルゴリズムとしては、例えば、位相的データ解析(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 abnormality diagnosis system 1 of the present embodiment diagnoses an abnormality of the rotating machine 2 (device) using an algorithm that can handle a large amount of data.

従来の解析技術では、大規模、多次元のデータを扱う場合に次元の削減の処理が必要である。例えば、統計的手法による次元削減では、データ中の情報が失われてしまうといった課題がある。また、教師あり学習を用いた異常診断についても、異常となるデータラベルが少ないと学習モデルを構築することが難しいという課題がある。   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 abnormality diagnosis system 1 obtains vibration data of the rotating machine 2 by attaching an inspection sensor 3 (for example, a vibration meter) to the rotating machine 2 provided in the plant. Note that the monitoring data, which is data indicating the operating state of the rotating machine 2, includes data such as sensors for collecting operation sounds and thermography, in addition to vibration data. The inspection sensor 3 may be a permanent one, or may be installed every time a measurement is performed.

モニタリングデータは、機器としての回転機2に取り付けられた点検用センサ3から発生する信号の時間波形と周波数スペクトルとのうちの少なくとも1つのデータである。また、モニタリングデータは、信号の時間波形または周波数スペクトルの絶対値でも良いし、信号の特徴を代表する量またはその組み合わせであっても良い。   The monitoring data is at least one of a time waveform and a frequency spectrum of a signal generated from the inspection sensor 3 attached to the rotating machine 2 as an apparatus. Further, the monitoring data may be an absolute value of a time waveform or a frequency spectrum of the signal, or an amount representing a characteristic of the signal or a combination thereof.

なお、振動データで対象となる状態量としては、例えば、変位、速度、加速度などがある。また、サンプリング点数または周波数の範囲など、任意の測定条件を使用者が選択可能である。このようにすれば、点検用センサ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 rotating machine 2 can be diagnosed based on the vibration data that can be easily acquired by the inspection sensor 3.

モニタリングデータとしては、例えば、点検用センサ3で計測された振動データを例示する。なお、モニタリングデータには、常時、回転機2の状態をモニタリングしている連続的なデータの他に、定期的に測定した断続的なデータも含まれる。以下の説明では、振動データを例に説明するが、その他のデータでも良い。   As the monitoring data, for example, vibration data measured by the inspection sensor 3 is exemplified. Note that the monitoring data includes not only continuous data that constantly monitors the state of the rotating machine 2 but also intermittent data that is periodically measured. In the following description, vibration data will be described as an example, but other data may be used.

なお、モニタリングデータには、回転機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 rotating machine 2, secondary secondary data generated by the operation of the rotating machine 2. . For example, in the case of the rotating machine 2, the number of rotations is primary data. Then, vibration data other than the main purpose event generated by the rotation operation becomes secondary data. In the present embodiment, an example in which abnormality of the rotating machine 2 is diagnosed using secondary data will be described. Note that the abnormality of the rotating machine 2 may be diagnosed using the primary data.

そして、このモニタリングデータを異常診断システム1に入力して解析を行う。例えば、モニタリングデータは、インターネット4に接続されたユーザ端末5(例えば、PC)を介して異常診断システム1に入力される。なお、モニタリングデータは、異常診断システム1に直接入力されても良い。   Then, the monitoring data is input to the abnormality diagnosis system 1 for analysis. For example, the monitoring data is input to the abnormality diagnosis system 1 via a user terminal 5 (for example, a PC) connected to the Internet 4. The monitoring data may be directly input to the abnormality diagnosis system 1.

この異常診断システム1は、プラント外のクラウドサーバに実装されても良いし、プラント内のパーソナルコンピュータまたはサーバに実装されても良い。また、モニタリングデータを、サーバなどの電気機器または所定のメモリに収集した後に、これら収集したデータを異常診断システム1に転送しても良い。   This abnormality diagnosis system 1 may be mounted on a cloud server outside the plant, or may be mounted on a personal computer or server inside the plant. Moreover, after collecting the monitoring data in an electric device such as a server or a predetermined memory, the collected data may be transferred to the abnormality diagnosis system 1.

点検用センサ3は、例えば、X軸、Y軸、Z軸の3軸方向の加速度を検出する。これにより、回転機2の加速度の振動(データ)波形を取得することができる。   The inspection sensor 3 detects, for example, accelerations in three X-axis, Y-axis, and Z-axis directions. Thereby, the vibration (data) waveform of the acceleration of the rotating machine 2 can be obtained.

次に、異常診断システム1のシステム構成を図1から図8に示すブロック図を参照して説明する。   Next, the system configuration of the abnormality diagnosis system 1 will be described with reference to the block diagrams shown in FIGS.

図1に示すように、異常診断システム1は、メイン制御部6と情報入力部7と通信部8と表示出力部9と解析処理ユニット10とを備える。これらは、メモリまたはHDDに記憶されたプログラムがCPUによって実行されることで実現される。   As shown in FIG. 1, the abnormality diagnosis system 1 includes a main control unit 6, an information input unit 7, a communication unit 8, a display output unit 9, and an analysis processing unit 10. These are realized by executing a program stored in the memory or the HDD by the CPU.

さらに、異常診断システム1は、モニタリングデータベース11とプロセスデータベース12とスペックデータベース13と異常原因データベース14とを備える。これらは、メモリまたはHDDに記憶され、検索または蓄積ができるよう整理された情報の集まりである。   Further, the abnormality diagnosis system 1 includes a monitoring database 11, a process database 12, a specification database 13, and an abnormality cause database 14. These are collections of information stored in a memory or HDD and organized so that they can be searched or stored.

本実施形態の異常診断システム1は、CPU、ROM、RAM、HDDなどのハードウェア資源を有し、CPUが各種プログラムを実行することで、ソフトウェアによる情報処理がハードウェア資源を用いて実現されるコンピュータで構成される。さらに、本実施形態の異常診断方法は、プログラムをコンピュータに実行させることで実現される。   The abnormality diagnosis system 1 according to the present embodiment has hardware resources such as a CPU, a ROM, a RAM, and an HDD. When the CPU executes various programs, information processing by software is realized using the hardware resources. It is composed of a computer. Further, the abnormality diagnosis method according to the present embodiment is realized by causing a computer to execute a program.

メイン制御部6は、異常診断システム1を統括的に制御する。   The main control unit 6 controls the abnormality diagnosis system 1 overall.

情報入力部7は、異常診断システム1を使用する使用者の操作に応じて所定の情報が入力される。この情報入力部7には、マウスまたはキーボードなどの入力装置が含まれる。つまり、これら入力装置の操作に応じて所定の情報が情報入力部7に入力される。この情報入力部7には、モニタリングデータなどの所定の情報が入力される。   The information input unit 7 receives predetermined information according to an operation of a user who uses the abnormality diagnosis system 1. The information input unit 7 includes an input device such as a mouse or a keyboard. That is, predetermined information is input to the information input unit 7 in accordance with the operation of these input devices. Predetermined information such as monitoring data is input to the information input unit 7.

通信部8は、インターネット4などの通信回線を介してユーザ端末5などの電子機器と通信を行う。この通信部8は、無線通信装置に搭載されても良い。また、通信部8は、所定のネットワーク機器、例えば、無線LANアクセスポイントまたはアンテナに搭載されても良い。なお、通信部8は、WAN(Wide Area Network)または携帯通信網を介してユーザ端末5と通信を行っても良い。   The communication unit 8 communicates with an electronic device such as the user terminal 5 via a communication line such as the Internet 4. This communication unit 8 may be mounted on a wireless communication device. The communication unit 8 may be mounted on a predetermined network device, for example, a wireless LAN access point or an antenna. Note that the communication unit 8 may communicate with the user terminal 5 via a WAN (Wide Area Network) or a mobile communication network.

表示出力部9は、解析結果を示す各種情報を表示する。例えば、クラスタ表示(図12参照)、時系列表示(図13参照)などの画面が表示される。つまり、表示出力部9は、分析結果を表示するためのグラフィックユーザインターフェース(GUI)である。   The display output unit 9 displays various information indicating the analysis result. For example, screens such as a cluster display (see FIG. 12) and a time-series display (see FIG. 13) are displayed. That is, the display output unit 9 is a graphic user interface (GUI) for displaying the analysis result.

本実施形態の異常診断システム1には、解析結果の出力を行うディスプレイなどの表示装置が含まれる。つまり、表示出力部9は、ディスプレイに表示される画面の制御を行う。なお、ディスプレイはコンピュータ本体と別体であっても良いし、一体であっても良い。さらに、ネットワークを介して接続される他のコンピュータが備えるディスプレイに表示される画面の制御を表示出力部9が行っても良い。   The abnormality diagnosis system 1 according to the present embodiment includes a display device such as a display that outputs an analysis result. That is, the display output unit 9 controls a screen displayed on the display. The display may be separate from the computer main body, or may be integral with the computer. Further, the display output unit 9 may control a screen displayed on a display of another computer connected via a network.

なお、本実施形態では、表示装置としてディスプレイを例示するが、その他の態様であっても良い。例えば、ヘッドマウントディスプレイまたはプロジェクタを用いて情報の表示を行っても良い。さらに、紙媒体に情報を印字するプリンタをディスプレイの替りとして用いても良い。つまり、表示出力部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 display output unit 9 may include a head mounted display, a projector, or a printer.

解析処理ユニット10は、特定次元数のデータを扱うアルゴリズムを用いてモニタリングデータに含まれる対象データまたは過去データの解析を行う(図2参照)。なお、本実施形態では、回転機2の運転中の状態を示すデータとして、モニタリングデータを例示しているが、モニタリングデータとプロセスデータを合せて作成されたデータセットが、回転機2の運転中の状態を示すデータであっても良い。以下の説明において、モニタリングデータという用語には、プロセスデータが含まれる場合がある。   The analysis processing unit 10 analyzes target data or past data included in the monitoring data using an algorithm that handles data of a specific number of dimensions (see FIG. 2). In the present embodiment, the monitoring data is exemplified as the data indicating the operating state of the rotating machine 2. However, a data set created by combining the monitoring data and the process data is used during the operation of the rotating machine 2. May be data indicating the state of. In the following description, the term monitoring data may include process data.

ここで、本実施形態の特定次元数は、4次元以上である。なお、例えば、解析処理ユニット10は、位相的データ解析とクラスタ解析とのうち少なくともいずれか一方を用いて解析を行う。このようにすれば、特定次元数のデータを扱った解析を行うことができる。   Here, the number of specific dimensions in the present embodiment is four or more. Note that, for example, the analysis processing unit 10 performs analysis using at least one of topological data analysis and cluster analysis. In this way, it is possible to perform analysis using data of a specific number of dimensions.

対象データとは、モニタリングデータに含まれる診断の対象となるデータである。なお、モニタリングデータと同時期または前後の時期に取得したプロセスデータも対象データに含まれる場合がある。   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 machine 2 that is currently operating does not know whether there is any abnormality.

このモニタリングデータは、対象データとして異常の有無の診断の対象となる。また、数カ月前または数年前に回転機2から取得されたモニタリングデータは、そのデータが取得された後に、実際に回転機2が正常に動作したか、それとも異常な状態となったかが把握できている。このようなモニタリングデータは、過去データとして扱われる。なお、過去のモニタリングデータと同時期または前後の時期に取得したプロセスデータも過去データに含まれる場合がある。   This monitoring data is to be diagnosed as abnormal data as target data. In addition, the monitoring data acquired from the rotating machine 2 several months or several years ago can be grasped after the data is acquired, whether the rotating machine 2 actually operates normally or is in an abnormal state. I have. Such monitoring data is treated as past data. It should be noted that the process data acquired at the same time as or around the same time as the past monitoring data may be included in the past data.

モニタリングデータベース11は、使用者が点検用センサ3を用いて回転機2から取得したモニタリングデータが記憶される。   The monitoring database 11 stores monitoring data obtained from the rotating machine 2 by the user using the inspection sensor 3.

プロセスデータベース12は、回転機2に関連するプロセスを示すプロセスデータが記憶される。プロセスデータには、例えば、回転数、流量、圧力、温度、電流などの運転条件と、運転期間、使用期間などの使用条件などが含まれる。なお、このプロセスデータは、使用者が記憶させても良いし、予め記憶されているものでも良い。   The process database 12 stores process data indicating a process related to the rotating machine 2. The process data includes, for example, operating conditions such as rotation speed, flow rate, pressure, temperature, and current, and operating conditions such as an operating period and a use period. The process data may be stored by the user or may be stored in advance.

スペックデータベース13は、回転機2のスペックを示すスペックデータが記憶される。スペックデータには、例えば、回転機2の場合は、型式、メーカー、サイズ、揚程、流量種別、軸受形式、羽数などの情報が含まれる。なお、このスペックデータは、使用者が記憶させても良いし、予め記憶されているものでも良い。   The specification database 13 stores specification data indicating the specifications of the rotating machine 2. The specification data includes, for example, in the case of the rotating machine 2, information such as a model, a maker, a size, a head, a flow type, a bearing type, and the number of blades. The specification data may be stored by the user or may be stored in advance.

異常原因データベース14は、解析処理ユニット10の解析に基づいて抽出された対象データの特徴量と、回転機2の異常原因との対応関係を示す原因情報(異常原因データ)が記憶される。原因情報は、異常の特徴量と異常原因との相関関係を示すものである。   The abnormality cause database 14 stores cause information (abnormal cause data) indicating the correspondence between the feature amount of the target data extracted based on the analysis of the analysis processing unit 10 and the abnormal cause of the rotating machine 2. The cause information indicates a correlation between the feature amount of the abnormality and the cause of the abnormality.

なお、この原因情報は、使用者が相関関係を入力、設定しても良いし、予め記憶されているものでも良い。このようにすれば、特徴量と異常原因との対応関係に基づいて、回転機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 rotating machine 2 can be easily estimated based on the correspondence between the feature value and the cause of the abnormality.

図2に示すように、解析処理ユニット10は、データ取得部15とデータ前処理部16と処理設定部17とデータ解析部18とパラメータチューニング部19と結果整理部20と結果表示制御部21と表示調整部22と異常原因分析部23と原因情報取得部24と原因表示制御部25とを備える。   As shown in FIG. 2, the analysis processing unit 10 includes a data acquisition unit 15, a data preprocessing unit 16, a processing setting unit 17, a data analysis unit 18, a parameter tuning unit 19, a result organizing unit 20, a result display control unit 21, The apparatus includes a display adjustment unit 22, an abnormality cause analysis unit 23, a cause information acquisition unit 24, and a cause display control unit 25.

データ取得部15は、モニタリングデータベース11とプロセスデータベース12とスペックデータベース13とから各種データを取得する。なお、以下の説明において、モニタリングデータという用語には、データ取得部15が取得したプロセスデータまたはスペックデータが含まれる場合がある。   The data acquisition unit 15 acquires various data from the monitoring database 11, the process database 12, and the specification database 13. In the following description, the term monitoring data may include process data or specification data acquired by the data acquisition unit 15.

データ前処理部16は、モニタリングデータから、対象データまたは過去データなどの解析に用いるデータを抽出し、データ解析部18で解析される前に、これらデータの前処理を行う(図3参照)。ここで、対象データまたは過去データの次元数が特定次元数になるように調整される。なお、対象データと過去データとは、同一の特定次元数になるように調整される。   The data preprocessing unit 16 extracts data used for analysis of target data or past data from the monitoring data, and performs preprocessing of these data before being analyzed by the data analysis unit 18 (see FIG. 3). Here, the number of dimensions of the target data or the past data is adjusted to be a specific number of dimensions. Note that the target data and the past data are adjusted so as to have the same specific number of dimensions.

処理設定部17は、データ解析部18で解析を行うときに用いる各種設定または解析用のアルゴリズムの設定を行う(図4参照)。   The process setting unit 17 performs various settings or an algorithm for analysis used when the data analysis unit 18 performs analysis (see FIG. 4).

データ解析部18は、特定次元数に調整された対象データまたは過去データを、特定次元数のデータを扱うアルゴリズムである位相的データ解析とクラスタ解析とのうち少なくともいずれか一方を用いて解析する。このデータ解析部18は、モニタリングデータに含まれる対象データの中から、異常な状態を示す対象データを検出することができる。   The data analysis unit 18 analyzes the target data or the past data adjusted to the specific dimension number using at least one of a topological data analysis and a cluster analysis, which are algorithms for handling data of the specific dimension number. The data analyzer 18 can detect target data indicating an abnormal state from target data included in the monitoring data.

パラメータチューニング部19は、データ解析部18で解析を行うときに用いる各種設定または解析用のアルゴリズムの設定のチューニングを行う(図5参照)。   The parameter tuning unit 19 tunes various settings or analysis algorithm settings used when the data analysis unit 18 performs analysis (see FIG. 5).

結果整理部20は、データ解析部18により得られたデータの解析結果が最適なものであるか否かを整理する。つまり、データ解析部18で分析した結果について、分析パラメータと分析結果との対応関係を整理する。   The result organizing unit 20 organizes whether or not the analysis result of the data obtained by the data analyzing unit 18 is optimal. That is, for the results analyzed by the data analysis unit 18, the correspondence between the analysis parameters and the analysis results is arranged.

ここで、データの解析結果が最適でない場合は、パラメータチューニング部19にて設定のチューニングが行われる。つまり、データ解析部18によるデータの解析に基づいて、回転機2の異常の有無を識別可能な態様で対象データを表示可能な解析結果が得られるまで、パラメータチューニング部19によるチューニングとデータ解析部18によるデータの解析とが繰り返される。   If the data analysis result is not optimal, the parameter tuning unit 19 tunes the setting. In other words, the tuning and the data analysis by the parameter tuning unit 19 are performed until an analysis result capable of displaying the target data in a form in which the presence or absence of abnormality of the rotating machine 2 can be identified based on the data analysis by the data analysis unit 18 is obtained. The analysis of the data by 18 is repeated.

結果表示制御部21は、データ解析部18により得られたデータの解析結果を示す画面であって、回転機2の異常の有無を識別可能な態様の画面を表示出力部9に表示させる制御を行う(図6参照)。例えば、クラスタ表示(図12参照)、時系列表示(図13参照)などの画面を表示出力部9に表示させる制御を行う。   The result display control unit 21 controls the display output unit 9 to display a screen showing a result of analysis of the data obtained by the data analysis unit 18 in a form in which the presence or absence of an abnormality of the rotating machine 2 can be identified. (See FIG. 6). For example, control is performed to cause the display output unit 9 to display a screen such as a cluster display (see FIG. 12) and a time-series display (see FIG. 13).

ここで、本実施形態の異常の有無を識別可能な態様は、例えば、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 rotating machine 2. The time-series display (see FIG. 13) is included in a mode in which the presence or absence of an abnormality can be identified.

位相的データ解析を用いる場合において、結果表示制御部21は、ノード数、エッジ数、入力したデータの総数、ノードの中に含まれるユニークなデータ数、ノード群数などを含む異常検出結果の整理を行う。なお、以下の説明では、ノード群をクラスタと称する場合がある。また、結果表示制御部21は、距離のメトリック、フィルター関数、インターバル、オーバーラップ、クラスタリングの手法などを含む位相的データ解析の分析パラメータの整理も行う。   When the topological data analysis is used, the result display control unit 21 sorts the abnormality detection results including the number of nodes, the number of edges, the total number of input data, the number of unique data included in the nodes, the number of node groups, and the like. I do. In the following description, a node group may be referred to as a cluster. The result display control unit 21 also arranges analysis parameters for topological data analysis including a metric of distance, a filter function, an interval, an overlap, a clustering method, and the like.

表示調整部22は、使用者の操作に応じて表示出力部9に表示させる画面の制御を行う(図7参照)。   The display adjustment unit 22 controls a screen to be displayed on the display output unit 9 according to a user operation (see FIG. 7).

異常原因分析部23は、データ解析部18により得られたデータの解析結果に基づいて、対象データに異常がある場合に、その異常原因の分析を行う(図8参照)。   The abnormality cause analysis unit 23 analyzes the cause of the abnormality when the target data has an abnormality based on the analysis result of the data obtained by the data analysis unit 18 (see FIG. 8).

原因情報取得部24は、異常原因データベース14から原因情報を取得する。この原因情報に基づいて、異常原因分析部23にて分析が行われる。   The cause information acquisition unit 24 acquires cause information from the abnormality cause database 14. An analysis is performed by the abnormality cause analysis unit 23 based on the cause information.

原因表示制御部25は、異常原因分析部23により得られたデータの解析結果を示す画面であって、回転機2の異常原因を使用者が診断可能な態様の画面を表示出力部9に表示させる制御を行う。   The cause display control unit 25 displays on the display output unit 9 a screen showing an analysis result of the data obtained by the abnormality cause analysis unit 23, in which the user can diagnose the cause of the abnormality of the rotating machine 2. Is performed.

図3に示すように、データ前処理部16は、対象データ調整部26と過去データ調整部27と欠損情報補間部28とデータラベル付与部29とデータスケーリング部30とスペクトル解析部31と振幅値算出部32とデータフォーマット変換部33とを備える。   As shown in FIG. 3, the data pre-processing unit 16 includes a target data adjustment unit 26, a past data adjustment unit 27, a loss information interpolation unit 28, a data label assignment unit 29, a data scaling unit 30, a spectrum analysis unit 31, and an amplitude value. A calculation unit 32 and a data format conversion unit 33 are provided.

対象データ調整部26は、モニタリングデータまたはプロセスデータに含まれる診断の対象となる対象データの次元数が特定次元数になるように対象データを調整する。   The target data adjustment unit 26 adjusts the target data so that the number of dimensions of the target data to be diagnosed included in the monitoring data or the process data becomes a specific number of dimensions.

過去データ調整部27は、モニタリングデータまたはプロセスデータに含まれる異常の有無が既知の過去データの次元数が特定次元数になるように過去データを調整する。   The past data adjustment unit 27 adjusts the past data so that the number of dimensions of the past data for which the presence or absence of an abnormality included in the monitoring data or the process data is known is a specific number of dimensions.

なお、本実施形態の「特定次元数になるようにデータを調整する」とは、任意の次元になるように、使用者がデータを間引く、または、次元圧縮することが含まれ得る。さらに、回転機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 rotating machine 2 is included. For example, when there is 2000-dimensional frequency data, a pre-process of extracting (adjusting the dimensions) only 48-dimensional data that is likely to contribute to an abnormality may be performed.

欠損情報補間部28は、対象データまたは過去データを特定次元数にしたときに、それぞれの次元数の情報のうち一部の次元数の情報が欠損する場合に、その欠損した情報を補間する。このようにすれば、対象データまたは過去データのそれぞれの次元数の情報を揃えて、特定次元数のデータを扱うアルゴリズムで解析可能な対象データまたは過去データを作成することができる。   The missing information interpolating unit 28 interpolates the missing information when the target data or the past data has a specific number of dimensions and, when the information of some of the information of the respective dimensions is missing. By doing so, it is possible to prepare target data or past data that can be analyzed by an algorithm that handles data of a specific number of dimensions, by aligning the information of the respective dimensions of the target data or past data.

データラベル付与部29は、特定次元数が調整された対象データまたは過去データに所定のラベルを付与する。例えば、過去に分析したデータ、または、データ収集の際に有意な特徴が観測されたデータに関しては、正常のときのデータか、異常が有るときのデータかが既知ものがある。このような過去データには、異常の有無を識別可能な情報としてのラベルLが付与される(図9および図10参照)。つまり、過去データには、異常の有無を識別可能なラベルLが対応付けられる。なお、対象データは、ラベルLをブランクにする。   The data label assigning unit 29 assigns a predetermined label to target data or past data whose specific dimension number has been adjusted. For example, with respect to data analyzed in the past or data in which a significant feature was observed during data collection, there are known data that are normal data or abnormal data. A label L is attached to such past data as information that can identify the presence or absence of an abnormality (see FIGS. 9 and 10). That is, the label L that can identify the presence or absence of the abnormality is associated with the past data. It should be noted that the label L is blank for the target data.

データスケーリング部30は、モニタリングデータを正規化または標準化する処理を行う。この処理により、次元毎でスケールの異なるデータのスケールを揃えることができる。例えば、正規化処理を行うと、特徴量の値の範囲が一定の範囲に収まるようになる。主に[0,1]か、[−1,1]の範囲内に収めることが多い。また、標準化では特徴量の平均を0、分散を1にする変換を行う。   The data scaling unit 30 performs a process of normalizing or standardizing the monitoring data. By this processing, the scale of data having different scales for each dimension can be made uniform. For example, when the normalization processing is performed, the range of the value of the feature amount falls within a certain range. It is often within the range [0, 1] or [-1, 1]. In the standardization, conversion is performed so that the average of the feature amounts is 0 and the variance is 1.

スペクトル解析部31は、振動データから任意の周波数成分を抽出する。この解析に基づいて、対象データまたは過去データになる振動スペクトルデータを生成することができる。このスペクトル解析部31は、例えば、高速フーリエ変換(FFT:Fast Fourier Transform)を用いて、モニタリングデータの特定の周波数成分を抽出したり、包絡線処理を行い、特定の周期的な周波数成分を抽出したりする。   The spectrum analysis unit 31 extracts an arbitrary frequency component from the vibration data. Based on this analysis, vibration spectrum data that becomes target data or past data can be generated. The spectrum analysis unit 31 extracts a specific frequency component of the monitoring data, performs an envelope process, and extracts a specific periodic frequency component using, for example, a fast Fourier transform (FFT). Or

振幅値算出部32は、モニタリングデータに含まれる振幅値を解析する。なお、この解析に基づいて、対象データまたは過去データになる振幅値データを作成することができる。   The amplitude value calculator 32 analyzes the amplitude value included in the monitoring data. It should be noted that amplitude value data that becomes target data or past data can be created based on this analysis.

つまり、本実施形態では、モニタリングデータに基づいて、回転機2の振動に関する情報である振幅値データまたは振動スペクトルデータが取得可能である。   That is, in the present embodiment, amplitude value data or vibration spectrum data, which is information relating to vibration of the rotating machine 2, can be acquired based on the monitoring data.

データフォーマット変換部33は、データ解析部18で用いるアルゴリズムに合わせたデータフォーマットになるように、対象データまたは過去データのデータフォーマットを変換する。なお、データフォーマットの変換とは、異常分析を行うためにデータの形状を調整することである。   The data format conversion unit 33 converts the data format of the target data or the past data so that the data format conforms to the algorithm used in the data analysis unit 18. The conversion of the data format is to adjust the shape of the data in order to perform the abnormality analysis.

図4に示すように、処理設定部17は、前処理条件設定部34とアルゴリズム設定部35とを備える。   As shown in FIG. 4, the processing setting unit 17 includes a preprocessing condition setting unit 34 and an algorithm setting unit 35.

前処理条件設定部34は、前処理の条件を設定する。例えば、前述の補間または規格化の処理などの設定が変更される。   The pre-processing condition setting unit 34 sets pre-processing conditions. For example, settings for the above-described interpolation or normalization processing are changed.

アルゴリズム設定部35は、データ解析部18で扱うアルゴリズムの設定を行う。例えば、位相的データ解析を用いるか、クラスタ解析を用いるか、それとも両解析を用いるかの設定を行う。またそれらのアルゴリズムのハイパーパラメータの設定を行う。このハイパーパラメータの設定により、例えば、位相的解析を用いる場合は、フィルター関数、オーバーラップ率などが変更される。   The algorithm setting unit 35 sets an algorithm handled by the data analysis unit 18. For example, it is set whether to use topological data analysis, use cluster analysis, or use both analyses. It also sets hyperparameters for those algorithms. For example, when the topological analysis is used, a filter function, an overlap ratio, and the like are changed by setting the hyperparameter.

図5に示すように、パラメータチューニング部19は、特定次元数設定部36とパラメータ最適化部38とを備える。   As shown in FIG. 5, the parameter tuning unit 19 includes a specific dimension number setting unit 36 and a parameter optimizing unit 38.

特定次元数設定部36は、対象データまたは過去データの解析を行うときに用いる特定次元数の設定を行う。このようにすれば、異常の有無を識別可能な態様になるように特定次元数の設定を行うことができる。なお、特定次元数の設定は、前処理条件設定部34またはアルゴリズム設定部35で行っても良い。   The specific dimension number setting unit 36 sets a specific dimension number used when analyzing target data or past data. By doing so, it is possible to set the specific number of dimensions so that the presence or absence of the abnormality can be identified. The setting of the specific number of dimensions may be performed by the preprocessing condition setting unit 34 or the algorithm setting unit 35.

また、解析後に得られる対象データまたは過去データの表示が、異常の有無を識別可能な態様になるように特定次元数の設定が繰り返される。このようにすれば、異常の有無を識別可能な態様になるまで、特定次元数の設定を繰り返し行うことができる。この設定は、人手で行っても良いし、自動的に行っても良い。   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 parameter optimizing unit 38 tunes the setting of parameters used for analysis in the data analyzing unit 18 based on the correspondence between the analysis parameters arranged in the result organizing unit 20 and the analysis results. In this tuning, parameter change is repeated and parameter optimization is performed so that the display of target data or past data obtained after analysis is in a form in which the presence or absence of an abnormality can be identified.

図6に示すように、結果表示制御部21は、クラスタ表示制御部39と時系列表示制御部40とを備える。   As shown in FIG. 6, the result display control unit 21 includes a cluster display control unit 39 and a time-series display control unit 40.

クラスタ表示制御部39は、データ解析部18により得られたデータの解析結果に基づいて、クラスタ表示(図12参照)の画面を表示出力部9に表示させる制御を行う。なお、クラスタ表示は、対象データをノードNとして表示した画面である。このようにすれば、対象データをクラスタ表示として視覚化できるため、回転機2の異常の有無の識別が行い易くなる。   The cluster display control unit 39 controls the display output unit 9 to display a screen for cluster display (see FIG. 12) based on the analysis result of the data obtained by the data analysis unit 18. The cluster display is a screen displaying the target data as the node N. By doing so, the target data can be visualized as a cluster display, and it is easy to identify whether or not there is an abnormality in the rotating machine 2.

時系列表示制御部40は、データ解析部18により得られたデータの解析結果を示す時系列表示(図13参照)の画面を表示出力部9に表示させる制御を行う。なお、時系列表示は、対象データをモニタリングされた時系列順に表示し、かつ対象データが含まれるノードNを識別可能な態様で表示した画面である。このようにすれば、対象データを時系列に並べて表示できるため、回転機2の異常の有無の識別が行い易くなる。   The time-series display control unit 40 controls the display output unit 9 to display a screen of a time-series display (see FIG. 13) showing the analysis result of the data obtained by the data analysis unit 18. The time-series display is a screen in which the target data is displayed in the order of the monitored time series, and the nodes N including the target data are displayed in an identifiable manner. In this way, the target data can be displayed in chronological order, so that it is easy to identify whether or not there is an abnormality in the rotating machine 2.

図7に示すように、表示調整部22は、ノード選択受付部41と対象情報表示部42とを備える。   As illustrated in FIG. 7, the display adjustment unit 22 includes a node selection reception unit 41 and a target information display unit 42.

ノード選択受付部41は、クラスタ表示された画面上のノードNまたはクラスタCの選択を受け付ける処理を行う。   The node selection accepting unit 41 performs a process of accepting selection of a node N or a cluster C on a cluster-displayed screen.

対象情報表示部42は、ノード選択受付部41により受け付けられ、選択されたノードNに含まれる内訳情報Bを表示する(図12参照)。この内訳情報Bには、選択されたノードNに含まれる対象データの対象情報が含まれる。このようにすれば、ノードNに対応する対象データの情報を容易に得ることができる。   The target information display unit 42 displays the breakdown information B received by the node selection reception unit 41 and included in the selected node N (see FIG. 12). The breakdown information B includes the target information of the target data included in the selected node N. This makes it possible to easily obtain information on the target data corresponding to the node N.

図8に示すように、異常原因分析部23は、特徴量抽出部43と異常原因推定部44とを備える。   As illustrated in FIG. 8, the abnormality cause analysis unit 23 includes a feature amount extraction unit 43 and an abnormality cause estimation unit 44.

特徴量抽出部43は、データ解析部18で解析された後の対象データの特徴量を抽出する。   The feature amount extraction unit 43 extracts a feature amount of the target data analyzed by the data analysis unit 18.

異常原因推定部44は、特徴量抽出部43にて抽出された特徴量に基づいて、異常原因を推定する。つまり、特徴量抽出部43にて抽出された特徴量に基づいて、異常原因を診断可能な態様で対象データを表示可能な状態になるようにデータを整理する。このようにすれば、回転機2に異常が発見された場合に、その原因を推定することができる。   The abnormality cause estimating unit 44 estimates an abnormality cause based on the feature amount extracted by the feature amount extracting unit 43. That is, based on the feature amount extracted by the feature amount extraction unit 43, the data is arranged so that the target data can be displayed in a mode in which the cause of the abnormality can be diagnosed. In this way, when an abnormality is found in the rotating machine 2, the cause can be estimated.

次に、データ解析部18で解析される対象データまたは過去データについて図9および図10を用いて説明する。   Next, target data or past data analyzed by the data analysis unit 18 will be described with reference to FIGS.

図9は、データ前処理部16にて前処理された後の加速度振幅のデータの一例を含むテーブルである。このテーブルでは、モニタリングされた日付(測定日、データ収集日時)に対応付けて、特定次元数の加速度振幅が登録されている。なお、加速度振幅の特定次元数は、任意の数を選択可能である。ここで、最新の日付のデータが、異常の有無の診断の対象となる対象データTとなっている。   FIG. 9 is a table including an example of acceleration amplitude data that has been preprocessed by the data preprocessing unit 16. In this table, the acceleration amplitude of a specific number of dimensions is registered in association with the monitored date (measurement date, data collection date and time). It should be noted that any number of specific dimensions of the acceleration amplitude can be selected. Here, the data of the latest date is the target data T to be diagnosed as to whether there is an abnormality.

なお、特定次元数の加速度振幅の列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 data preprocessing unit 16. In this table, vibration values for each specific number of dimensions are registered in association with the monitored date. Note that an arbitrary number can be selected as the specific dimension number of the frequency. Here, the data of the latest date is the target data T to be diagnosed as to whether there is an abnormality.

なお、特定次元数の周波数の列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 information interpolation unit 28. The value of the data to be interpolated may be zero, the average value of the data included in the same row as the missing item, or the interpolation based on the data before and after this average value. Or, it may be extrapolated, the average value of the data included in the same column as the missing item, or the one interpolated or extrapolated based on the data before and after this average value, Other arbitrarily determined values may be used.

なお、本実施形態の過去データとしては、異常の診断対象となる回転機2から取得されたデータを用いる形態を例示するが、その他の態様であっても良い。例えば、メーカー、スペック、型式が同一構成の回転機2がある場合には、この他の回転機2から取得されたデータを、異常の診断対象となる回転機2の過去データとして用いても良い。   As the past data of the present embodiment, a mode using data acquired from the rotating machine 2 to be diagnosed as abnormal is exemplified, but other modes may be used. For example, when there is a rotating machine 2 having the same configuration, manufacturer, specifications, and model, data acquired from the other rotating machine 2 may be used as past data of the rotating machine 2 to be diagnosed as abnormal. .

なお、特定次元数は、全ての日付で同一の次元数となるように、データ前処理部16においてデータ補完およびフォーマット変換が行われる。なお、データ前処理部16の処理は、自動的に行っても良いし、人手により任意に設定しても良い。   The data pre-processing unit 16 performs data complement and format conversion so that the specific dimension number has the same dimension number for all dates. The processing of the data pre-processing unit 16 may be performed automatically or may be set arbitrarily by hand.

従来の機械学習を用いた統計的手法では、統計的な解析をベースにしており、取り扱うデータが正規分布などの何らかの分布に従っていることを前提としている。そのため、個別データの特徴に応じた解析手法の構築が必要ある。また、大規模、多次元のデータを扱う場合には、データの次元の削減などの処理が必要である。このような統計的手法による次元削減では、データ中の情報が失われてしまうといった課題がある。本実施形態は、このような課題を解決することができる。   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 parameter tuning unit 19 and the result sorting unit 20 sort out the correlation between these control parameters and the analysis results (see FIG. 11). Then, the parameter tuning unit 19 sweeps the analysis parameters until the data point of interest is plotted on a node of the graph.

このパラメータチューニング操作によって、注目しているデータ点が、必ずノード上に表現されるようになる。次に、注目しているデータ点を含むノードが、別の測定日のデータ点を含むか、また別のノードとの関係があるか、すなわちエッジで結ばれているかを確認する。この際、別の測定日のデータが正常のラベル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 data analysis unit 18 and resulting from the TDA Mapper is displayed on the display. The cluster display is controlled by the cluster display control unit 39 of the result display control unit 21.

このクラスタ表示には、ノード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 rotating machine 2 has been disassembled and inspected and a label that has not been disassembled and inspected may be attached to the data, and color-coded. Further, the color of the node N may be classified for each parameter affected at the time of analysis.

また、クラスタ表示では、注目されるデータが異常の可能性が高いと判定された場合は、特異ノード強調表示機能で、異常を有する対象データを含むノード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 data analysis unit 18 and being a result of the TDA Mapper is displayed on the display. The control of the time-series display is performed by the time-series display control unit 40 of the result display control unit 21.

この時系列表示では、クラスタ表示にてグラフ化されていたデータが時系列順に表示される。例えば、時系列表示には、モニタリングされた日付(測定日、データ収集日時)の列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 display control unit 21 displays the target data in a mode (for example, time-series display) in which the presence or absence of an abnormality can be identified. Further, based on the analysis of the past data, the past data is displayed in such a manner that the associated label L can be identified. By doing so, it is possible to identify the presence or absence of an abnormality in the target data by referring to the display mode of the past data.

なお、クラスタ表示(図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 cause analysis unit 23 extracts a characteristic amount of data determined to have a possibility of abnormality by using a characteristic amount extracting function of abnormal data.

ここで、分析対象と分析したデータを周波数成分データとした場合について説明する。例えば、異常と判定されたデータと、正常ラベルがついたデータに関して、各周波数成分の差を算出し、各周波数成分の差が大きい周波数を抽出する。この抽出した周波数と周波数成分の差を異常データの特徴量とする。差の大きい周波数成分の定義は、予め設定した閾値を超えるもの、または、上位の数%の差を抽出するなどの任意の抽出方法を選択することができる。そして、抽出した特徴量(周波数、周波数成分の差)を異常原因分析部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 cause analysis unit 23.

本実施形態では、異常原因分析部23に特徴量が入力されることで、異常原因への特徴量の影響度が分かり、異常原因を推定することが可能となる。また、異常検知に必要なパラメータチューニングを自動化することで、異常検知を人手によるチューニングを行わないで、異常原因を推定することが可能となる。また、異常検知の結果から異常の原因を推定することが可能となる。   In the present embodiment, by inputting the characteristic amount to the abnormality cause analysis unit 23, the degree of influence of the characteristic amount on the abnormality cause can be known, and the abnormality cause can be estimated. Further, by automating parameter tuning necessary for abnormality detection, it is possible to estimate the cause of the abnormality without manually performing the abnormality detection. Further, the cause of the abnormality can be estimated from the result of the abnormality detection.

次に、回転機2の異常診断方法について図14のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。   Next, a method of diagnosing abnormality of the rotating machine 2 will be described with reference to the flowchart of FIG. Note that FIG. 1 to FIG. 13 are appropriately referred to.

図14に示すように、まず、ステップS11において、異常診断システム1の使用者は、回転機2に点検用センサ3を取り付ける。   As shown in FIG. 14, first, in step S11, the user of the abnormality diagnosis system 1 attaches the inspection sensor 3 to the rotating machine 2.

次のステップS12において、使用者は、回転機2に取り付けた点検用センサ3を用いて、回転機2の運転中の状態をモニタリングしてモニタリングデータを得る。   In the next step S12, the user monitors the operating state of the rotating machine 2 using the inspection sensor 3 attached to the rotating machine 2 to obtain monitoring data.

次のステップS13において、使用者は、点検用センサ3を用いて取得したモニタリングデータを、ユーザ端末5を用いて異常診断システム1のモニタリングデータベース11に入力する。   In the next step S13, the user inputs the monitoring data acquired using the inspection sensor 3 to the monitoring database 11 of the abnormality diagnosis system 1 using the user terminal 5.

次のステップS14において、使用者は、診断対象となる回転機2のプロセスデータを、ユーザ端末5を用いて異常診断システム1のプロセスデータベース12に入力する。   In the next step S <b> 14, the user inputs the process data of the rotating machine 2 to be diagnosed to the process database 12 of the abnormality diagnosis system 1 using the user terminal 5.

次のステップS15において、使用者は、診断対象となる回転機2のスペックデータまたはプロセスデータを、ユーザ端末5を用いて異常診断システム1のスペックデータベース13に入力する。   In the next step S15, the user inputs the specification data or process data of the rotating machine 2 to be diagnosed to the specification database 13 of the abnormality diagnosis system 1 using the user terminal 5.

次のステップS16において、使用者は、診断対象となる回転機2の異常原因データを、ユーザ端末5を用いて異常診断システム1の異常原因データベース14に入力する。
In the next step S16, the user inputs the abnormality cause data of the rotating machine 2 to be diagnosed to the abnormality cause database 14 of the abnormality diagnosis system 1 using the user terminal 5.

次のステップ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 analysis processing unit 10 and the abnormality cause of the rotating machine 2. I do. That is, in the abnormality cause database 14, the correspondence between the feature amount of the target data that can be extracted and the abnormality cause is stored in advance.

次のステップS18において、使用者は、異常診断システム1を用いて、モニタリングデータに含まれる対象データまたは過去データの解析を行う。   In the next step S <b> 18, the user uses the abnormality diagnosis system 1 to analyze target data or past data included in the monitoring data.

次のステップ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 abnormality diagnosis system 1. The user performs, for example, an operation of displaying a cluster display screen, an operation of displaying a time-series display screen, an operation of selecting a node N or a cluster C, and an operation of displaying a screen showing a result of analyzing the cause of abnormality.

次のステップS20において、使用者は、異常診断システム1のディスプレイに表示された情報に基づいて、回転機2の異常の有無の識別を行う。   In the next step S20, the user identifies the presence or absence of an abnormality in the rotating machine 2 based on the information displayed on the display of the abnormality diagnosis system 1.

次のステップS21において、使用者は、異常診断システム1のディスプレイに表示された情報に基づいて、回転機2の異常原因の特定を行う。そして、異常診断方法を終了する。   In the next step S21, the user specifies the cause of the abnormality of the rotating machine 2 based on the information displayed on the display of the abnormality diagnosis system 1. Then, the abnormality diagnosis method ends.

次に、解析処理ユニット10が実行する異常診断処理について図15のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。   Next, the abnormality diagnosis processing executed by the analysis processing unit 10 will be described with reference to the flowchart of FIG. Note that FIG. 1 to FIG. 13 are appropriately referred to.

図15に示すように、まず、ステップS31において、解析処理ユニット10のデータ取得部15は、モニタリングデータベース11に記憶されたモニタリングデータを取得する。ここで、プロセスデータベース12に記憶されたプロセスデータ、およびスペックデータベース13に記憶されたスペックデータの取得を行っても良い。   As shown in FIG. 15, first, in step S31, the data acquisition unit 15 of the analysis processing unit 10 acquires monitoring data stored in the monitoring database 11. Here, the process data stored in the process database 12 and the specification data stored in the specification database 13 may be obtained.

次のステップS32において、解析処理ユニット10の処理設定部17の前処理条件設定部34は、データ前処理部16にて実行される前処理条件の設定を行う。ここで、パラメータチューニング部19の特定次元数設定部36は、前処理条件設定部34の設定に基づいて、対象データまたは過去データの特定次元数の設定を行う。また、処理設定部17のアルゴリズム設定部35は、データ解析部18で用いられるアルゴリズムの設定を行う。アルゴリズムとしては、例えば、TDA Mapperを用いる。   In the next step S32, the preprocessing condition setting unit 34 of the processing setting unit 17 of the analysis processing unit 10 sets preprocessing conditions to be executed by the data preprocessing unit 16. Here, the specific dimension number setting unit 36 of the parameter tuning unit 19 sets the specific dimension number of the target data or the past data based on the setting of the preprocessing condition setting unit 34. The algorithm setting unit 35 of the processing setting unit 17 sets an algorithm used in the data analysis unit 18. For example, TDA Mapper is used as the algorithm.

次のステップS33において、解析処理ユニット10のデータ前処理部16は、データ前処理を実行する。   In the next step S33, the data pre-processing unit 16 of the analysis processing unit 10 executes data pre-processing.

次のステップS34において、解析処理ユニット10のデータ解析部18は、データ解析処理を実行する。ここで、特定次元数に調整された対象データを、特定次元数のデータを扱うアルゴリズムを用いて解析する。さらに、特定次元数に調整された過去データを、アルゴリズムを用いて解析する。   In the next step S34, the data analysis unit 18 of the analysis processing unit 10 executes a data analysis process. Here, the target data adjusted to the specific number of dimensions is analyzed using an algorithm that handles data of the specific number of dimensions. Further, the past data adjusted to the specific number of dimensions is analyzed using an algorithm.

なお、アルゴリズムとして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 data analysis unit 18.

次のステップS35において、解析処理ユニット10の結果整理部20は、解析結果整理処理を実行する。ここで、データ解析部18により得られたデータの解析結果が最適なものであるか否かを整理する。   In the next step S35, the result organizing unit 20 of the analysis processing unit 10 executes an analysis result organizing process. Here, whether or not the analysis result of the data obtained by the data analysis unit 18 is optimal is arranged.

次のステップS36において、結果整理部20は、解析結果が最適であるか否かを判定する。なお、最適であるか否かの判定は、例えば、対象データの表示が異常の有無を識別可能な態様になるか否かに基づいて設定される。ここで、解析結果が最適である場合(ステップS36がYES)は、後述のステップS38に進む。一方、解析結果が最適でない場合(ステップS36がNO)は、ステップS37に進む。   In the next step S36, the result organizing unit 20 determines whether or not the analysis result is optimal. The determination as to whether or not the data is optimal is set based on, for example, whether or not the display of the target data is in a form in which the presence or absence of an abnormality can be identified. Here, when the analysis result is optimal (step S36 is YES), the process proceeds to step S38 described later. On the other hand, if the analysis result is not optimal (NO in step S36), the process proceeds to step S37.

ステップS37において、解析処理ユニット10のパラメータチューニング部19は、パラメータチューニング処理を実行する。ここで、データ解析部18で解析を行うときに用いる各種設定または解析用のアルゴリズムの設定のチューニングが行われる。そして、前述のステップS33に戻る。つまり、本実施形態では、対象データの表示が異常の有無を識別可能な態様になるように特定次元数の設定を繰り返す。   In step S37, the parameter tuning unit 19 of the analysis processing unit 10 performs a parameter tuning process. Here, tuning of various settings used for analysis by the data analysis unit 18 or setting of an algorithm for analysis is performed. Then, the process returns to step S33. That is, in the present embodiment, the setting of the specific dimension number is repeated so that the display of the target data is in a form in which the presence or absence of the abnormality can be identified.

ステップS38において、解析処理ユニット10の異常原因分析部23は、異常原因分析処理を実行する。ここで、特徴量抽出部43は、データ解析部18で解析された後の対象データの特徴量を抽出する。また、異常原因推定部44は、特徴量抽出部43にて抽出された特徴量に基づいて、異常原因を推定し、その異常原因を診断可能な態様で対象データを表示可能な状態になるようにデータを整理する。   In step S38, the abnormality cause analysis unit 23 of the analysis processing unit 10 executes an abnormality cause analysis process. Here, the feature value extraction unit 43 extracts the feature value of the target data analyzed by the data analysis unit 18. Further, the abnormality cause estimating unit 44 estimates the abnormality cause based on the characteristic amount extracted by the characteristic amount extracting unit 43, so that the target data can be displayed in such a manner that the abnormal cause can be diagnosed. Organize your data.

次のステップS39において、解析処理ユニット10の表示調整部22は、表示調整処理を実行する。ここで、表示調整部22は、使用者の操作に基づいて、ディスプレイに表示される画面の制御を行う。そして、異常診断処理を終了する。   In the next step S39, the display adjustment unit 22 of the analysis processing unit 10 executes a display adjustment process. Here, the display adjustment unit 22 controls a screen displayed on the display based on a user operation. Then, the abnormality diagnosis processing ends.

次に、データ前処理部16が実行するデータ前処理について図16のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。   Next, data pre-processing performed by the data pre-processing unit 16 will be described with reference to the flowchart in FIG. Note that FIG. 1 to FIG. 13 are appropriately referred to.

図16に示すように、まず、ステップS41において、データ前処理部16の対象データ調整部26は、モニタリングデータに含まれる診断の対象となる対象データの次元数が特定次元数になるように対象データを調整する処理を行う。   As shown in FIG. 16, first, in step S41, the target data adjustment unit 26 of the data pre-processing unit 16 sets the target data so that the number of dimensions of the target data to be diagnosed included in the monitoring data becomes the specific number of dimensions. Perform data adjustment processing.

次のステップS42において、データ前処理部16の過去データ調整部27は、モニタリングデータに含まれる異常の有無が既知の過去データの次元数が特定次元数になるように過去データを調整する処理を行う。   In the next step S42, the past data adjusting unit 27 of the data pre-processing unit 16 performs a process of adjusting the past data so that the number of dimensions of the past data for which the presence or absence of an abnormality included in the monitoring data is known is a specific number of dimensions. Do.

次のステップS43において、データ前処理部16の欠損情報補間部28は、対象データを特定次元数にしたときに、それぞれの次元数の情報のうち一部の次元数の情報が欠損する場合に、その欠損した情報を補間する処理を行う。   In the next step S43, the missing information interpolating unit 28 of the data pre-processing unit 16 determines whether or not some of the information of the respective dimensions is missing when the target data has the specified number of dimensions. , A process of interpolating the missing information.

次のステップS44において、データ前処理部16のデータラベル付与部29は、特定次元数に調整された過去データに、異常の有無を識別可能なラベルLを対応付ける処理を行う。   In the next step S44, the data label assigning unit 29 of the data pre-processing unit 16 performs a process of associating the past data adjusted to the specific number of dimensions with a label L capable of identifying the presence or absence of an abnormality.

次のステップS45において、データ前処理部16は、その他の処理を実行する。例えば、データスケーリング部30がデータを規格化する。スペクトル解析部31がスペクトル成分を算出する。振幅値算出部32が振幅値を算出する。データフォーマット変換部33がデータフォーマットを変換する。そして、データ前処理を終了する。   In the next step S45, the data preprocessing unit 16 performs other processing. For example, the data scaling unit 30 normalizes the data. The spectrum analyzer 31 calculates a spectrum component. The amplitude value calculator 32 calculates an amplitude value. The data format conversion unit 33 converts the data format. Then, the data pre-processing ends.

次に、表示調整部22が実行する表示調整処理について図17のフローチャートを用いて説明する。なお、図1から図13を適宜参照する。   Next, the display adjustment processing executed by the display adjustment unit 22 will be described with reference to the flowchart in FIG. Note that FIG. 1 to FIG. 13 are appropriately referred to.

図17に示すように、まず、ステップS51において、表示調整部22は、ディスプレイにクラスタ表示(グラフ)の画面が表示中であるか否かを判定する。ここで、グラフが表示中でない場合(ステップS51がNO)は、後述のステップS54に進む。一方、グラフが表示中である場合(ステップS51がYES)は、ステップS52に進む。   As shown in FIG. 17, first, in step S51, the display adjustment unit 22 determines whether or not a cluster display (graph) screen is being displayed on the display. If the graph is not being displayed (NO in step S51), the process proceeds to step S54 described below. On the other hand, when the graph is being displayed (YES in step S51), the process proceeds to step S52.

ステップS52において、表示調整部22のノード選択受付部41は、クラスタ表示された画面上のノードNまたはクラスタCの選択を受け付けたか否かを判定する。ここで、ノードNまたはクラスタCの選択を受け付けていない場合(ステップS52がNO)は、ステップS54に進む。一方、ノードNまたはクラスタCの選択を受け付けた場合(ステップS52がYES)は、ステップS53に進む。   In step S52, the node selection receiving unit 41 of the display adjusting unit 22 determines whether the selection of the node N or the cluster C on the cluster-displayed screen has been received. Here, when the selection of the node N or the cluster C is not received (step S52 is NO), the process proceeds to step S54. On the other hand, when the selection of the node N or the cluster C is received (YES in step S52), the process proceeds to step S53.

ステップS53において、表示調整部22の対象情報表示部42は、選択されたノードNまたはクラスタCに含まれる対象データまたは過去データの内訳情報B(対象情報)を内訳情報表示領域Dに表示する処理を行う。そして、ステップS54に進む。   In step S53, the target information display unit 42 of the display adjustment unit 22 displays the target data included in the selected node N or cluster C or the breakdown information B (target information) of past data in the breakdown information display area D. I do. Then, the process proceeds to step S54.

ステップS54において、表示調整部22は、時系列表示の画面をディスプレイに表示させる操作を受け付けたか否かを判定する。ここで、時系列表示の表示操作の受け付けがない場合(ステップS54がNO)は、後述のステップS56に進む。一方、時系列表示の表示操作の受け付けがある場合(ステップS54がYES)は、ステップS55に進む。   In step S54, the display adjustment unit 22 determines whether or not an operation for displaying a time-series display screen on the display has been received. Here, when the display operation of the time series display is not received (NO in step S54), the process proceeds to step S56 described later. On the other hand, when the display operation of the time-series display is received (YES in step S54), the process proceeds to step S55.

ステップS55において、表示調整部22の処理に基づいて、結果表示制御部21の時系列表示制御部40が、時系列表示の画面をディスプレイに表示させる。この時系列表示制御部40は、対象データをモニタリングされた時系列順に表示し、かつ対象データが含まれるノードを識別可能な態様で表示する。なお、時系列表示制御部40は、過去データの解析に基づいて、対応付けられたラベルLが識別可能な態様で過去データを表示する。例えば、異常可能性を示す表示を行う(図13参照)。そして、ステップS56に進む。   In step S55, the time-series display control unit 40 of the result display control unit 21 displays a time-series display screen on the display based on the processing of the display adjustment unit 22. The time-series display control unit 40 displays the target data in the order of the monitored time series, and displays the nodes including the target data in an identifiable manner. In addition, the time-series display control unit 40 displays the past data in a form in which the associated label L can be identified based on the analysis of the past data. For example, a display indicating the possibility of abnormality is performed (see FIG. 13). Then, the process proceeds to step S56.

ステップS56において、表示調整部22は、クラスタ表示の画面をディスプレイに表示させる操作を受け付けたか否かを判定する。ここで、クラスタ表示の表示操作の受け付けがない場合(ステップS56がNO)は、後述のステップS58に進む。一方、クラスタ表示の表示操作の受け付けがある場合(ステップS56がYES)は、ステップS57に進む。   In step S56, the display adjustment unit 22 determines whether or not an operation to display a cluster display screen on the display has been received. Here, when the display operation of the cluster display is not received (NO in step S56), the process proceeds to step S58 described later. On the other hand, if a cluster display operation has been received (YES in step S56), the process advances to step S57.

ステップS57において、表示調整部22の処理に基づいて、結果表示制御部21のクラスタ表示制御部39が、クラスタ表示の画面をディスプレイに表示させる。このクラスタ表示制御部39は、対象データの解析に基づいて、異常の有無を識別可能な態様で対象データを表示する。なお、クラスタ表示制御部39は、過去データの解析に基づいて、対応付けられたラベルLが識別可能な態様で過去データを表示する。例えば、異常可能性に応じてノードNを色分けして表示する(図12参照)。そして、ステップS58に進む。   In step S57, based on the processing of the display adjustment unit 22, the cluster display control unit 39 of the result display control unit 21 displays a cluster display screen on the display. The cluster display control unit 39 displays the target data in a manner that allows the presence or absence of an abnormality to be identified based on the analysis of the target data. The cluster display control unit 39 displays the past data in such a manner that the associated label L can be identified based on the analysis of the past data. For example, the nodes N are displayed in different colors depending on the possibility of abnormality (see FIG. 12). Then, the process proceeds to step S58.

ステップS58において、表示調整部22は、異常原因分析結果を示す画面をディスプレイに表示させる操作を受け付けたか否かを判定する。ここで、異常原因分析結果の表示操作の受け付けがない場合(ステップS58がNO)は、後述のステップS60に進む。一方、異常原因分析結果の表示操作の受け付けがある場合(ステップS58がYES)は、ステップS59に進む。   In step S58, the display adjustment unit 22 determines whether or not an operation to display a screen indicating the result of the analysis of the cause of the abnormality on the display has been received. Here, when the display operation of the abnormality cause analysis result is not received (NO in step S58), the process proceeds to step S60 described later. On the other hand, when the display operation of the abnormality cause analysis result is received (YES in step S58), the process proceeds to step S59.

ステップS59において、表示調整部22の処理に基づいて、原因表示制御部25が、異常原因分析結果を示す画面をディスプレイに表示させる。そして、ステップS60に進む。   In step S59, based on the processing of the display adjustment unit 22, the cause display control unit 25 causes the display to display a screen showing the result of the analysis of the cause of the abnormality. Then, the process proceeds to step S60.

ステップS60において、表示調整部22は、ディスプレイに表示中の解析結果を示す画面の表示を終了させる操作を受け付けたか否かを判定する。ここで、表示を終了させる操作の受け付けがない場合(ステップS60がNO)は、前述のステップS51に進む。一方、表示を終了させる操作の受け付けがある場合(ステップS60がYES)は、ステップS61に進む。   In step S60, the display adjustment unit 22 determines whether or not an operation for terminating the display of the screen showing the analysis result being displayed on the display has been received. Here, when the operation to end the display is not received (NO in step S60), the process proceeds to step S51 described above. On the other hand, when an operation to end the display is received (YES in step S60), the process proceeds to step S61.

ステップS61において、表示調整部22の処理に基づいて、結果表示制御部21または原因表示制御部25は、ディスプレイに表示中の解析結果を示す画面の表示を終了する。そして、表示調整処理を終了する。   In step S61, based on the processing of the display adjustment unit 22, the result display control unit 21 or the cause display control unit 25 ends the display of the screen showing the analysis result being displayed on the display. Then, the display adjustment processing ends.

なお、本実施形態のフローチャートにおいて、各ステップが直列に実行される形態を例示しているが、必ずしも各ステップの前後関係が固定されるものでなく、一部のステップの前後関係が入れ替わっても良い。また、一部のステップが他のステップと並列に実行されても良い。   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 parameter tuning unit 19.

また、パラメータチューニングを行い最適化する際の指標としては、例えば、クラスタの距離とコーフエン行列との相関係数(コーフエン相関係数)がある。   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 data preprocessing unit 16 is input to the data analysis unit 18, analyzed using a hierarchical clustering method, and the analysis result is displayed on the display output unit 9. Among the results, it is the result of cluster display. Here, 57 pieces of amplitude value data of k dimensions are input. Among the 57 amplitude value data, there are data in the normal range and data showing an abnormality.

例えば、図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 SYMBOLS 1 ... Abnormality diagnosis system, 2 ... Rotating machine, 3 ... Inspection sensor, 4 ... Internet, 5 ... User terminal, 6 ... Main control part, 7 ... Information input part, 8 ... Communication part, 9 ... Display output part, 10 ... Analysis processing unit, 11 ... Monitoring database, 12 ... Process database, 13 ... Spec database, 14 ... Abnormal cause database, 15 ... Data acquisition unit, 16 ... Data preprocessing unit, 17 ... Process setting unit, 18 ... Data analysis unit , 19: Parameter tuning unit, 20: Result arrangement unit, 21: Result display control unit, 22: Display adjustment unit, 23: Abnormal cause analysis unit, 24: Cause information acquisition unit, 25: Cause display control unit, 26: Target Data adjustment unit, 27: Past data adjustment unit, 28: Loss information interpolation unit, 29: Data label assignment unit, 30: Data scaling unit, 31: Spectrum Analysis unit, 32: amplitude value calculation unit, 33: data format conversion unit, 34: preprocessing condition setting unit, 35: algorithm setting unit, 36: specific dimension number setting unit, 38: parameter optimization unit, 39: cluster display Control unit, 40: time-series display control unit, 41: node selection reception unit, 42: target information display unit, 43: feature amount extraction unit, 44: abnormality cause estimation unit, B: breakdown information, C: cluster, D ... Breakdown information display area, E: edge, F: frequency column, G: graph display area, K1 to K4: time series display column, L: label, M: mouse cursor, N: node, P: past data, R ... color label contour, S: column of acceleration amplitude, T: target data.

Claims (14)

機器の状態を示すデータである信号の時間波形と周波数スペクトルのうちの少なくとも一方のデータから抽出され、前記機器の状態を示すデータのうち診断の対象となる1の対象データに含まれる個々のデータの個数である次元数特定次元数に調整するステップと、
前記特定次元数に調整された前記対象データを、前記特定次元数のデータを扱うアルゴリズムであるマッパーを用いて位相的データ解析を行うステップと、
前記マッパーを用いた解析結果を異常の有無を識別可能な態様であるクラスタ表示または時系列表示で表示するステップと、
を含む、
機器の異常診断方法。
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.
前記特定次元数は、4次元以上であり、
前記異常の有無を識別可能な態様は、前記クラスタ表示としての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 .
前記クラスタ表示された前記ノードの選択を受け付けるステップと、
選択された前記ノードに含まれる前記対象データの情報を表示するステップと、
を含む、
請求項に記載の機器の異常診断方法。
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 .
前記対象データをモニタリングされた時系列順に表示し、かつ前記対象データが含まれる前記ノードを識別可能な態様で表示するステップを含む、
請求項または請求項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 .
機器の状態を示すデータである信号の時間波形と周波数スペクトルのうちの少なくとも一方のデータから抽出され、前記機器の状態を示すデータのうち診断の対象となる1の対象データに含まれる個々のデータの個数である次元数特定次元数に調整するデータ前処理部と、
前記特定次元数に調整された前記対象データを、前記特定次元数のデータを扱うアルゴリズムであるマッパーを用いて位相的データ解析を行うデータ解析部と、
前記マッパーを用いた解析結果を異常の有無を識別可能な態様であるクラスタ表示または時系列表示で表示する結果表示制御部と、
を備える、
機器の異常診断システム。
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.
請求項13に記載の機器の異常診断システムを用いた機器の異常診断方法であって、
前記機器に取り付けた点検用センサを用いて、前記機器の状態をモニタリングしてモニタリングデータを得るステップと、
前記モニタリングデータを前記異常診断システムに入力するステップと、
を含む、
機器の異常診断方法。
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.
JP2019131439A 2019-07-16 2019-07-16 Device abnormality diagnosis method and device abnormality diagnosis system Active JP6641056B1 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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