JP2016091417A - Abnormality diagnostic method, abnormality diagnostic apparatus and abnormality diagnostic program - Google Patents

Abnormality diagnostic method, abnormality diagnostic apparatus and abnormality diagnostic program Download PDF

Info

Publication number
JP2016091417A
JP2016091417A JP2014227155A JP2014227155A JP2016091417A JP 2016091417 A JP2016091417 A JP 2016091417A JP 2014227155 A JP2014227155 A JP 2014227155A JP 2014227155 A JP2014227155 A JP 2014227155A JP 2016091417 A JP2016091417 A JP 2016091417A
Authority
JP
Japan
Prior art keywords
data
unit
selection
diagnosis
abnormality diagnosis
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.)
Granted
Application number
JP2014227155A
Other languages
Japanese (ja)
Other versions
JP6492555B2 (en
Inventor
志保 袖子田
Shiho Sodekoda
志保 袖子田
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.)
IHI Corp
Original Assignee
IHI 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 IHI Corp filed Critical IHI Corp
Priority to JP2014227155A priority Critical patent/JP6492555B2/en
Publication of JP2016091417A publication Critical patent/JP2016091417A/en
Application granted granted Critical
Publication of JP6492555B2 publication Critical patent/JP6492555B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

PROBLEM TO BE SOLVED: To reduce dependence on a selector and easily select unit data.SOLUTION: An abnormality diagnostic apparatus 100 diagnoses whether a diagnosis object is abnormal on the basis of a dataset having parameter values for a plurality of items related to the diagnosis object while using a diagnostic scheme on the premise that the parameter values follow a normal distribution. The abnormality diagnostic apparatus 100 comprises a unit data selection part 21 that selects unit data D2 on the basis of which a diagnosis is conducted from selection data D1 including a plurality of datasets. The unit data selection part 21 includes a first calculation section 21a calculating a first range based on a standard deviation of the parameter values for each of the items in the datasets included in the selection data D1; and a selection section 21b selecting the datasets the parameter values of which are present in the first range for all the items out of the datasets included in the selection data D1.SELECTED DRAWING: Figure 1

Description

本発明は、異常診断方法、異常診断装置及び異常診断プログラムに関する。   The present invention relates to an abnormality diagnosis method, an abnormality diagnosis apparatus, and an abnormality diagnosis program.

従来の異常診断方法に関する技術として、例えば特許文献1には、複数の評価項目があるプラントの運転状態を監視するプラント監視方法が記載されている。特許文献1に記載された方法では、複数の評価項目毎の状態量(パラメータ値)の集まりである状態量の束(データセット)が記憶された状態量ファイルから、評価対象とする状態量の束を抽出する。そして、単位空間ファイルに記憶されたデータ(単位データ)を基準として、抽出した状態量の束のマハラノビス距離を求め、このマハラノビス距離が所定の閾値以内であるか否かに応じて、プラントの運転状態が正常か否かを判定する。   As a technique related to a conventional abnormality diagnosis method, for example, Patent Document 1 describes a plant monitoring method for monitoring an operation state of a plant having a plurality of evaluation items. In the method described in Patent Literature 1, a state quantity file that is a collection of state quantities (data sets) that is a collection of state quantities (parameter values) for each of a plurality of evaluation items is stored. Extract the bundle. Then, based on the data (unit data) stored in the unit space file, the Mahalanobis distance of the extracted bundle of state quantities is obtained, and the operation of the plant is determined depending on whether or not the Mahalanobis distance is within a predetermined threshold. Determine whether the state is normal.

特開2012−252556号公報JP 2012-252556 A

上記従来技術では、診断を行う際に基準となる単位データの選定には、明確な方法が存在していない。通常、単位データとしては、診断対象に異常が発生していない期間のデータセットや、診断対象における複数のセンサの各センサ値(パラメータ値)が正常範囲にあるデータセット等が選定される。しかし、このような単位データの選定では、選定作業が選定者の経験や感覚に依存する可能性がある。また、例えば数百もの多数のセンサが存在すると、これらのセンサのセンサ値が正常範囲にあるか否かを判別すべき場合があることから、選定作業に多大な手間及び時間がかかる可能性がある。   In the above prior art, there is no clear method for selecting unit data as a reference when making a diagnosis. Usually, as unit data, a data set in a period in which no abnormality has occurred in the diagnosis target, a data set in which the sensor values (parameter values) of a plurality of sensors in the diagnosis target are in a normal range, and the like are selected. However, in such unit data selection, there is a possibility that the selection work depends on the experience and feeling of the selector. In addition, for example, if there are hundreds of sensors, it may be necessary to determine whether or not the sensor values of these sensors are in the normal range. is there.

そこで、本発明は、選定者による依存を低減して簡易に単位データを選定することができる異常診断方法、異常診断装置及び異常診断プログラムを提供することを課題とする。   Therefore, an object of the present invention is to provide an abnormality diagnosis method, an abnormality diagnosis apparatus, and an abnormality diagnosis program that can easily select unit data while reducing the dependence on a selector.

本発明に係る異常診断方法は、診断対象に関する複数項目のパラメータ値を有するデータセットに基づき、パラメータ値が正規分布に従うことを前提にした診断手法を用いて、診断対象の異常を診断する異常診断方法であって、複数のデータセットを含む選定用データから、診断を行う際に基準となる単位データを選定するデータ選定工程を備え、データ選定工程は、選定用データに含まれる複数のデータセットにおいて、複数項目別に、パラメータ値の標準偏差に基づく第1範囲を算出する工程と、選定用データに含まれる複数のデータセットのうち、複数項目の全てにおいてパラメータ値が第1範囲内に存在するデータセットを、単位データとして選定する工程と、を含む。   The abnormality diagnosis method according to the present invention is based on a data set having a plurality of parameter values related to a diagnosis target, and uses a diagnosis method based on the assumption that the parameter values follow a normal distribution. A method comprising a data selection step of selecting unit data serving as a reference when performing diagnosis from selection data including a plurality of data sets, wherein the data selection step includes a plurality of data sets included in the selection data. In the step, the first range based on the standard deviation of the parameter value is calculated for each of the plurality of items, and the parameter values are present in the first range in all of the plurality of items among the plurality of data sets included in the selection data. Selecting a data set as unit data.

本発明に係る異常診断装置は、診断対象に関する複数項目のパラメータ値を有するデータセットに基づき、パラメータ値が正規分布に従うことを前提にした診断手法を用いて、診断対象の異常を診断する異常診断装置であって、複数のデータセットを含む選定用データから、診断を行う際に基準となる単位データを選定する単位データ選定部を備え、単位データ選定部は、選定用データに含まれる複数のデータセットにおいて、複数項目別に、パラメータ値の標準偏差に基づく第1範囲を算出する第1算出部と、選定用データに含まれる複数のデータセットのうち、複数項目の全てにおいてパラメータ値が第1範囲内に存在するデータセットを、単位データとして選定する選定部と、を含む。   The abnormality diagnosis apparatus according to the present invention is based on a data set having parameter values of a plurality of items related to a diagnosis target, and uses a diagnosis method based on the assumption that the parameter values follow a normal distribution. The apparatus includes a unit data selection unit that selects unit data serving as a reference when performing diagnosis from selection data including a plurality of data sets, and the unit data selection unit includes a plurality of data included in the selection data. In the data set, for each of a plurality of items, the first calculation unit that calculates the first range based on the standard deviation of the parameter value, and the parameter value is the first in all of the plurality of items among the plurality of data sets included in the selection data. And a selection unit that selects a data set existing within the range as unit data.

このような異常診断方法及び異常診断装置では、選定用データに含まれる複数のデータセットのうち、複数項目の全てにおいてパラメータ値が第1範囲内に存在するデータセット(つまり、第1範囲から外れたパラメータ値を有するデータセット以外)が、単位データとして選定される。よって、選定者の判断を介さずに、且つ、各パラメータ値毎に正常範囲にあるかどうかを判定することなく、選定用データから単位データを選定することができる。したがって、選定者による依存を低減して簡易に単位データを選定することが可能となる。   In such an abnormality diagnosis method and an abnormality diagnosis apparatus, among the plurality of data sets included in the selection data, a data set in which parameter values exist in the first range in all of the plurality of items (that is, out of the first range). Other than the data set having the parameter value) is selected as the unit data. Therefore, unit data can be selected from the selection data without going through the judgment of the selector and without determining whether each parameter value is in the normal range. Therefore, it is possible to easily select unit data while reducing the dependency on the selector.

本発明に係る異常診断方法は、選定用データの一部を除外するデータ除外工程をさらに備え、データ除外工程は、選定用データに含まれる複数のデータセットにおいて、複数項目別に、パラメータ値の尖度とパラメータ値の標準偏差に基づく第2範囲とを算出する工程と、尖度が所定値以上の場合、選定用データに含まれる複数のデータセットから、第2範囲から外れたパラメータ値を有するデータセットを除外する工程と、を含んでいてもよい。   The abnormality diagnosis method according to the present invention further includes a data excluding process for excluding a part of the selection data, and the data excluding process includes a parameter value for each item in a plurality of data sets included in the selection data. A step of calculating the second range based on the standard deviation of the degree and the parameter value, and when the kurtosis is equal to or greater than a predetermined value, the plurality of data sets included in the selection data have parameter values that are out of the second range And excluding the data set.

本発明に係る異常診断装置は、選定用データの一部を除外する選定用データ除外部をさらに備え、選定用データ除外部は、選定用データに含まれる複数のデータセットにおいて、複数項目別に、パラメータ値の尖度とパラメータ値の標準偏差に基づく第2範囲とを算出する第2算出部と、尖度が所定値以上の場合、選定用データに含まれる複数のデータセットから、第2範囲から外れたパラメータ値を有するデータセットを除外する除外部と、を含んでいてもよい。   The abnormality diagnosis apparatus according to the present invention further includes a selection data exclusion unit that excludes a part of the selection data, and the selection data exclusion unit is divided into a plurality of items in a plurality of data sets included in the selection data. A second calculation unit that calculates a kurtosis of the parameter value and a second range based on the standard deviation of the parameter value; and a second range from a plurality of data sets included in the selection data when the kurtosis is equal to or greater than a predetermined value And an exclusion unit that excludes a data set having a parameter value that is out of the range.

このような異常診断方法及び異常診断装置では、尖度が所定値以上の場合、第2範囲から外れたパラメータ値(いわゆる外れ値)を有するデータセットが、選定用データから除外される。これにより、選定用データの正規性を高めることができる。   In such an abnormality diagnosis method and abnormality diagnosis apparatus, when the kurtosis is equal to or greater than a predetermined value, a data set having a parameter value (so-called outlier) that is out of the second range is excluded from the selection data. Thereby, the normality of the data for selection can be improved.

本発明に係る異常診断方法は、診断対象から診断用のデータセットである診断用データを取得するデータ取得工程と、単位データ及び診断用データに基づき、診断対象が異常か或いは正常かを診断する診断工程と、単位データを選定用データとして複製し、診断工程で診断対象が正常と診断されている場合に当該選定用データに診断用データを追加するデータ更新工程と、を備えてもよい。   The abnormality diagnosis method according to the present invention diagnoses whether a diagnosis target is abnormal or normal based on a data acquisition step of acquiring diagnostic data, which is a diagnostic data set from a diagnosis target, and unit data and diagnostic data. A diagnosis process and a data update process for replicating unit data as selection data and adding diagnosis data to the selection data when the diagnosis target is diagnosed as normal in the diagnosis process may be provided.

本発明に係る異常診断装置は、診断対象から診断用のデータセットである診断用データを取得する診断用データ取得部と、単位データ及び診断用データに基づき、診断対象が異常か或いは正常かを診断する異常診断部と、単位データを選定用データとして複製し、異常診断部で診断対象が正常と診断されている場合に当該選定用データに診断用データを追加する選定用データ更新部と、を備えてもよい。   An abnormality diagnosis apparatus according to the present invention is based on a diagnosis data acquisition unit that acquires diagnosis data that is a data set for diagnosis from a diagnosis object, and based on unit data and diagnosis data, whether the diagnosis object is abnormal or normal An abnormality diagnosis unit to diagnose, a selection data update unit that duplicates unit data as selection data, and adds diagnosis data to the selection data when the diagnosis target is diagnosed as normal by the abnormality diagnosis unit; May be provided.

このような異常診断方法及び異常診断装置では、異常診断の都度、選定用データを更新することが可能となる。   In such an abnormality diagnosis method and abnormality diagnosis device, it is possible to update the selection data every time an abnormality diagnosis is made.

本発明に係る異常診断方法及び異常診断装置において、診断手法は、マハラノビスタグチメソッドであり、単位データは、マハラノビス距離を算出する際に基準となる単位空間を構成してもよい。この場合、マハラノビスタグチメソッドによる異常診断において、マハラノビス距離を算出する際に基準となる単位空間を構成する単位データを、選定者による依存を低減して簡易に選定できる。   In the abnormality diagnosis method and the abnormality diagnosis apparatus according to the present invention, the diagnosis method may be a Mahalanobis Taguchi method, and the unit data may constitute a unit space serving as a reference when calculating the Mahalanobis distance. In this case, in the abnormality diagnosis by the Mahalanobis Taguchi method, the unit data constituting the unit space used as a reference when calculating the Mahalanobis distance can be easily selected with less dependence on the selector.

本発明に係る異常診断プログラムは、コンピュータに、上記異常診断方法を実行せしめる。この異常診断プログラムにおいても、上記異常診断方法における作用効果と同様な作用効果が奏される。   The abnormality diagnosis program according to the present invention causes a computer to execute the abnormality diagnosis method. Also in this abnormality diagnosis program, the same effects as the effects in the abnormality diagnosis method are exhibited.

本発明によれば、選定者による依存を低減して簡易に単位データを選定することができる異常診断方法、異常診断装置及び異常診断プログラムを提供することが可能となる。   ADVANTAGE OF THE INVENTION According to this invention, it becomes possible to provide the abnormality diagnosis method, abnormality diagnosis apparatus, and abnormality diagnosis program which can select unit data easily, reducing the dependence by a selector.

一実施形態に係る異常診断装置の構成を示すブロック図である。It is a block diagram which shows the structure of the abnormality diagnosis apparatus which concerns on one Embodiment. MT法の概念を示す図である。It is a figure which shows the concept of MT method. 選定用データの一例を説明する図である。It is a figure explaining an example of the data for selection. パラメータ値の分布例を示すグラフである。It is a graph which shows the example of distribution of a parameter value. 異常診断処理を示すフローチャートである。It is a flowchart which shows an abnormality diagnosis process. 単位データ決定処理を示すフローチャートである。It is a flowchart which shows a unit data determination process. (a)は選定用データ除外処理を示すフローチャートである。(b)は単位データ選定処理を示すフローチャートである。(A) is a flowchart which shows the data exclusion process for selection. (B) is a flowchart showing unit data selection processing. 選定用データ更新処理を示すフローチャートである。It is a flowchart which shows the data update process for selection. 異常診断プログラムを説明するための図である。It is a figure for demonstrating an abnormality diagnosis program.

以下、添付図面を参照して、本発明の実施形態について詳細に説明する。なお、以下の説明において、同一又は相当要素には同一符号を付し、重複する説明を省略する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings. In the following description, the same or equivalent elements will be denoted by the same reference numerals, and redundant description will be omitted.

図1は、本実施形態に係る異常診断装置の構成を示すブロック図である。図2は、MT法の概念を示す図である。図3は、選定用データの一例を説明する図である。図4は、選定用データに含まれる複数のデータセットにおいて、一の項目のパラメータ値の分布例を示すグラフである。   FIG. 1 is a block diagram showing the configuration of the abnormality diagnosis apparatus according to this embodiment. FIG. 2 is a diagram illustrating the concept of the MT method. FIG. 3 is a diagram illustrating an example of selection data. FIG. 4 is a graph showing a distribution example of parameter values of one item in a plurality of data sets included in the selection data.

図1に示すように、異常診断装置100は、診断対象に関する複数項目のパラメータ値を有するデータセットに基づいて、診断対象の異常を診断する。診断対象としては、特に限定されるものではなく、例えばガスタービン発電プラントをはじめとした各種のプラント、及び、複数の機器を有する機械システム等が挙げられる。本実施形態では、ガスタービンを備えた発電用のプラントが診断対象として適用されている。   As shown in FIG. 1, the abnormality diagnosis device 100 diagnoses an abnormality of a diagnosis target based on a data set having a plurality of parameter values related to the diagnosis target. The object to be diagnosed is not particularly limited, and examples thereof include various plants including a gas turbine power plant and a mechanical system having a plurality of devices. In the present embodiment, a power generation plant including a gas turbine is applied as a diagnosis target.

データセットは、複数項目のパラメータ値が集まってなるデータ群であって、サンプルデータとも称される。複数項目のパラメータ値とは、複数種の特徴量(変数)に係るデータ値である。パラメータ値の項目としては、例えば、診断対象が備える複数種の機器(バルブ、ポンプ等)の制御値、指令値及び応答値、並びに、診断対象に設けられた複数種のセンサの各センサ値等が挙げられる。パラメータ値の項目数は、特に限定されるものではないが、ここでは数百以上とされている。   The data set is a data group in which parameter values of a plurality of items are collected, and is also referred to as sample data. The parameter values of a plurality of items are data values related to a plurality of types of feature amounts (variables). The parameter value items include, for example, control values, command values and response values of a plurality of types of devices (valves, pumps, etc.) included in the diagnosis target, and sensor values of a plurality of types of sensors provided in the diagnosis target. Is mentioned. The number of parameter value items is not particularly limited, but is several hundred or more here.

本実施形態における異常診断では、データセットの各パラメータ値が正規分布に従うことを前提にした診断手法を用いている。当該診断手法では、正常時のデータである単位データが診断の基準として用いられ、単位データの内容が診断精度に影響する。このような診断手法としては、MT(Mahalanobis Taguchi)法が挙げられる。なお、パラメータ値が正規分布に従うことを前提にした診断手法として、上記の他に、例えば回帰分析等に用いてもよい。   In the abnormality diagnosis in the present embodiment, a diagnosis method is used on the assumption that each parameter value of the data set follows a normal distribution. In the diagnostic method, unit data that is normal data is used as a reference for diagnosis, and the content of the unit data affects the diagnostic accuracy. An example of such a diagnostic technique is the MT (Mahalanobis Taguchi) method. In addition to the above, as a diagnostic method on the assumption that the parameter value follows a normal distribution, for example, it may be used for regression analysis or the like.

本実施形態では、診断手法としてMT法を採用している。MT法は、品質工学において提案されているパターン認識手法の一つであって、いつもと同じ状態であるかどうかを診断する。MT法では、数百以上の特徴量(項目)を有するデータセットでも一つの判断指標に集約して扱い得る。図2に示すように、MT法による異常診断は、単位空間を単位データから生成し、この単位空間を診断の基準として、診断したい診断用データとのマハラノビス距離(Mahalanobis Distance、以下「MD値」という)を用いて異常か或いは正常かを判断する。   In the present embodiment, the MT method is adopted as a diagnostic method. The MT method is one of pattern recognition methods proposed in quality engineering, and diagnoses whether or not it is the same as usual. In the MT method, even a data set having several hundred or more feature amounts (items) can be handled as one judgment index. As shown in FIG. 2, in the abnormality diagnosis by the MT method, a unit space is generated from unit data, and the Mahalanobis distance (hereinafter referred to as “MD value”) from the diagnostic data to be diagnosed using this unit space as a reference for diagnosis. To determine whether it is abnormal or normal.

MD値が閾値よりも小さい場合、単位空間から近い状態、すなわち正常状態であるとして、正常と診断される。一方、MD値が閾値以上の場合、単位空間から遠い状態、すなわち正常状態から離れている状態であるとして、異常と診断される。単位空間は、目的に対して均質な正常な集団(ユニフォームな集団)であり、均質性を有する。この単位空間は、「定常且つ普通の状態」、「頻度の多い状態」又は「平均的な状態」ともいえる。   When the MD value is smaller than the threshold value, it is diagnosed as normal because it is close to the unit space, that is, in a normal state. On the other hand, when the MD value is equal to or greater than the threshold value, an abnormality is diagnosed as being in a state far from the unit space, that is, a state away from the normal state. The unit space is a normal group (uniform group) that is homogeneous for the purpose and has homogeneity. This unit space can also be said to be “steady and normal state”, “frequent state” or “average state”.

図1に示すように、異常診断装置100は、記憶装置10及びCPU(Central Processing Unit)20を備えている。記憶装置10は、選定用データD1及び単位データD2を記憶している。記憶装置10としては、ROM、フレキシブルディスク、CD−ROM、DVD、或いはUSBメモリ等の記録媒体、又は半導体メモリ等が例示される。   As shown in FIG. 1, the abnormality diagnosis apparatus 100 includes a storage device 10 and a CPU (Central Processing Unit) 20. The storage device 10 stores selection data D1 and unit data D2. Examples of the storage device 10 include a recording medium such as a ROM, a flexible disk, a CD-ROM, a DVD, or a USB memory, or a semiconductor memory.

選定用データD1は、単位データD2の選定に用いられるデータであって、複数のデータセットを含んでいる。ここでの選定用データD1は、過去にプラントから入力されて蓄積された複数のデータセットにより形成される。図3に示す例において、選定用データD1は、所定時間毎にプラントから入力されて蓄積された複数のデータセットdにより形成されている。なお、選定用データD1は、プラントから入力されたデータセットによらずに形成されもよく、別途に外部入力されて形成されてもよい。   The selection data D1 is data used for selection of the unit data D2, and includes a plurality of data sets. The selection data D1 here is formed by a plurality of data sets that have been input and accumulated in the past. In the example shown in FIG. 3, the selection data D1 is formed by a plurality of data sets d that are input from the plant and accumulated every predetermined time. The selection data D1 may be formed without depending on the data set input from the plant, or may be formed by external input separately.

選定用データD1に含まれるデータセットdの数は、特に限定されるものではなく、適宜設定することができる。図示する選定用データD1は、所定時間毎のデータセットdを含んでいるが、複数のデータセットを含んでいればよい。図1に戻り、単位データD2は、診断を行う際に基準となるデータであって、単位空間を構成する。単位データD2は、単位空間の生成に用いる正常な場合のデータセットを複数含んでいる。   The number of data sets d included in the selection data D1 is not particularly limited and can be set as appropriate. The selection data D1 shown in the figure includes a data set d for each predetermined time, but it is only necessary to include a plurality of data sets. Returning to FIG. 1, the unit data D2 is data used as a reference when making a diagnosis, and constitutes a unit space. The unit data D2 includes a plurality of normal data sets used for generating a unit space.

CPU20は、その機能的構成要素として、単位データ選定部21と選定用データ除外部22と診断用データ取得部23と異常診断部24と選定用データ更新部25とを有している。単位データ選定部21は、記憶装置10から選定用データD1を取得し、この選定用データD1から単位データD2を選定する。単位データ選定部21は、第1算出部21a、選定部21b及び修正部21cを含んでいる。   The CPU 20 includes a unit data selection unit 21, a selection data exclusion unit 22, a diagnosis data acquisition unit 23, an abnormality diagnosis unit 24, and a selection data update unit 25 as functional components. The unit data selection unit 21 acquires selection data D1 from the storage device 10, and selects unit data D2 from the selection data D1. The unit data selection unit 21 includes a first calculation unit 21a, a selection unit 21b, and a correction unit 21c.

第1算出部21aは、選定用データD1に含まれる複数のデータセットにおいて、その複数項目別に、パラメータ値の平均値Av、標準偏差σ及び第1範囲R1を含む統計量を算出する。第1範囲R1は、標準偏差σに基づく指定範囲である。具体的には、第1範囲R1は、指定値mと標準偏差σとの乗数からなる区間m・σに基づく範囲である。より具体的には、第1範囲R1は、平均値Avを中心とし当該平均値Avから正負側それぞれに区間m・σ拡がる信頼区間(=Av±m・σ)である(図4参照)。 The first calculation unit 21a calculates a statistic including an average value Av 1 of parameter values, a standard deviation σ 1 and a first range R1 for each of a plurality of items in a plurality of data sets included in the selection data D1. The first range R1 is a range specified based on the standard deviation sigma 1. Specifically, the first range R1 is a range based on an interval m · σ 1 that is a multiplier of a specified value m and a standard deviation σ 1 . More specifically, the first range R1 is an average value Av 1 of the center that the average value Av 1 interval m · sigma 1 spreads the positive and negative side, respectively from confidence intervals (= Av 1 ± m · σ 1) ( (See FIG. 4).

指定値mは、予め記憶された設定値である。指定値mは、不図示の入力装置を介して選定者(ユーザ)により指定されてもよいし、固定値として指定されていてもよい。指定値mは、パラメータ値の全項目で一律に同じ値とされてもよいし、少なくとも一部の項目で異なる値とされてもよい。指定値mは、例えば各項目別の正常時の分布情報に応じて、各項目で異なる値とされてもよい。また、指定値mは、診断対象の実際の運転状況に合わせて、各項目の変動範囲において適時修正してもよい。ここでの指定値mは、全ての項目で正の実数値である「2」が設定されており、その結果、複数項目の全てにおいて、「Av±2・σ」の範囲が第1範囲R1とされている。 The designated value m is a preset value stored in advance. The designated value m may be designated by a selector (user) via an input device (not shown), or may be designated as a fixed value. The designated value m may be uniformly set to the same value for all items of the parameter value, or may be set to a different value for at least some items. The designated value m may be a different value for each item, for example, according to normal distribution information for each item. In addition, the designated value m may be corrected in a timely manner within the variation range of each item in accordance with the actual driving situation to be diagnosed. The designated value m here is set to “2” which is a positive real value in all items, and as a result, the range of “Av 1 ± 2 · σ 1 ” is the first in all of the plurality of items. The range is R1.

選定部21bは、選定用データD1に含まれる複数のデータセットのうち、複数項目の全てにおいてパラメータ値が第1範囲R1内に存在するデータセット(すなわち、第1範囲R1から外れたパラメータ値を有するデータセット以外)を、単位データD2として選定する。換言すると、データセットにおける複数のパラメータ値全てが第1範囲R1内に在る場合に、そのデータセットが単位データD2に選定され、データセットにおける複数のパラメータ値の少なくとも1つでも第1範囲R1外に在る場合、そのデータセットは単位データD2に選定されない。第1範囲R1内に選定部21bは、選定した単位データD2を記憶装置10へ出力する。   The selection unit 21b includes, among the plurality of data sets included in the selection data D1, a data set in which parameter values exist in the first range R1 in all of the plurality of items (that is, parameter values that are out of the first range R1). A data set other than the data set) is selected as the unit data D2. In other words, when all of the plurality of parameter values in the data set are within the first range R1, the data set is selected as the unit data D2, and even at least one of the plurality of parameter values in the data set is the first range R1. If it is outside, the data set is not selected as unit data D2. Within the first range R1, the selection unit 21b outputs the selected unit data D2 to the storage device 10.

修正部21cは、選定部21bで選定され記憶装置10へ出力される単位データD2を、データ修正ユニット30を介して入力されたユーザによる操作入力に応じて修正する。データ修正ユニット30は、キーボードや操作ボタン等の入力部30aと、ディスプレイ等の表示部30bと、を有している。データ修正ユニット30では、選定部21bで選定され記憶装置10へ出力される単位データD2を表示部30bに表示すると共に、当該単位データD2を修正すべくユーザにより入力部30aを介して操作入力がなされた場合、その操作入力を修正部21cへ出力する。   The correction unit 21c corrects the unit data D2 selected by the selection unit 21b and output to the storage device 10 according to the operation input by the user input via the data correction unit 30. The data correction unit 30 includes an input unit 30a such as a keyboard and operation buttons, and a display unit 30b such as a display. In the data correction unit 30, the unit data D2 selected by the selection unit 21b and output to the storage device 10 is displayed on the display unit 30b, and an operation input is made by the user via the input unit 30a to correct the unit data D2. If so, the operation input is output to the correction unit 21c.

選定用データ除外部22は、診断開始時から未だ単位空間が生成されていない(異常診断部24において単位空間を生成する際、当該生成が初回となる)場合、記憶装置10から選定用データD1を取得し、この選定用データD1の一部を除外する。選定用データ除外部22は、第2算出部22a及び除外部22bを含んでいる。   If the unit space has not yet been generated from the start of diagnosis (when the unit space is generated in the abnormality diagnosis unit 24, the selection data excluding unit 22 is the first generation), the selection data exclusion unit 22 selects the selection data D1 from the storage device 10. And a part of the selection data D1 is excluded. The selection data excluding unit 22 includes a second calculating unit 22a and an excluding unit 22b.

第2算出部22aは、選定用データD1に含まれる複数のデータセットにおいて、その複数項目別に、パラメータ値の平均値Av、標準偏差σ、尖度γ及び第2範囲R2を含む統計量を算出する。尖度γは、パラメータ値の分布がどれほど平均値Av周辺に集中しているかを表すものである。尖度γは、他の値から大きく外れた値である外れ値がパラメータ値として含まれると、大きな値を示す傾向を有する。尖度γは、下式により求めることができる。
γ=((a−Av+(a−Av
+…(a−Av)/(n・σ
〜a:パラメータ値
The second calculation unit 22a includes, in a plurality of data sets included in the selection data D1, a statistic including an average value Av 2 , a standard deviation σ 2 , a kurtosis γ, and a second range R2 for each of the plurality of items. Is calculated. The kurtosis γ represents how much the parameter value distribution is concentrated around the average value Av 2 . The kurtosis γ tends to show a large value when an outlier, which is a value greatly deviating from other values, is included as a parameter value. The kurtosis γ can be obtained by the following equation.
γ = ((a 1 −Av 2 ) 2 + (a 2 −Av 2 ) 2
+ ... (a n −Av 2 ) 2 ) / (n · σ 2 2 )
a 1 to a n : Parameter value

第2範囲R2は、標準偏差σに基づく指定範囲である。具体的には、第2範囲R2は、指定値kと標準偏差σとの乗数からなる区間k・σに基づく範囲である。より具体的には、第2範囲R2は、平均値Avを中心とし当該平均値Avから正負側それぞれに区間k・σ拡がる信頼区間(=Av±k・σ)である。 The second range R2 is a specified range based on the standard deviation sigma 2. Specifically, the second range R2 is a range based on a section k · σ 2 that is a multiplier of the specified value k and the standard deviation σ 2 . More specifically, the second range R2 is the average value Av 2 was mainly the average value Av 2 interval k · sigma 2 spreads in the positive and negative side, respectively from confidence intervals (= Av 2 ± k · σ 2).

指定値kは、予め記憶された指定値である。指定値kは、不図示の入力装置を介してユーザにより指定されてもよいし、固定値として指定されていてもよい。指定値kは、パラメータ値の全項目で一律に同じ値とされてもよいし、少なくとも一部の項目で異なる値とされてもよい。指定値kは、例えば各項目別の正常時の分布情報に応じて、各項目で異なる値とされてもよい。また、指定値kは、診断対象の実際の運転状況に合わせて、各項目の変動範囲において適時修正してもよい。ここでの指定値kは、パラメータ値の全項目で、指定値mよりも大きい正の実数値である「3」が設定されており、その結果、複数項目の全てにおいて、「Av±3・σ」の範囲が第2範囲R2とされている。 The designated value k is a designated value stored in advance. The designated value k may be designated by the user via an input device (not shown) or may be designated as a fixed value. The designated value k may be set to the same value for all items of the parameter value, or may be set to a different value for at least some of the items. The designated value k may be a different value for each item according to, for example, normal distribution information for each item. Further, the designated value k may be corrected in a timely manner within the variation range of each item in accordance with the actual driving situation to be diagnosed. Here, the designated value k is set to “3” which is a positive real value larger than the designated value m for all items of the parameter value, and as a result, “Av 2 ± 3” is set for all of the plurality of items. The range of “σ 2 ” is the second range R2.

除外部22bは、第2算出部22aで算出した尖度γが所定値以上の場合、選定用データD1に外れ値が含まれていると判定し、第2範囲R2から外れたパラメータ値を外れ値とし、当該外れ値を有するデータセットを選定用データD1から除外する。換言すると、データセットにおける複数のパラメータ値全てが第2範囲R2内に在る場合に、そのデータセットは選定用データD1から除外されない一方、データセットにおける複数のパラメータ値の少なくとも1つでも第2範囲R2外に在る場合、そのデータセットは選定用データD1から除外される。除外部22bは、外れ値を有するデータセットを除外した選定用データD1を、記憶装置10へ出力する。所定値は、外れ値の有無を判定する上で一般的又は経験的に予め設定される値である。所定値は、例えば5以上の値とされ、より確実に判定する場合には10以上の値とされる。   When the kurtosis γ calculated by the second calculation unit 22a is equal to or greater than a predetermined value, the excluding unit 22b determines that the selection data D1 includes an outlier, and deviates from the parameter value outside the second range R2. The data set having the outlier is excluded from the selection data D1. In other words, when all of the plurality of parameter values in the data set are within the second range R2, the data set is not excluded from the selection data D1, while at least one of the plurality of parameter values in the data set is the second value. If it is outside the range R2, the data set is excluded from the selection data D1. The excluding unit 22b outputs the selection data D1 excluding the data set having an outlier to the storage device 10. The predetermined value is a value generally or empirically set in advance when determining the presence or absence of an outlier. The predetermined value is, for example, a value of 5 or more, and is a value of 10 or more in the case of more reliable determination.

診断用データ取得部23は、診断用のデータセットである診断用データを、診断対象から取得する。診断用データ取得部23は、取得した診断用データを異常診断部24及び選定用データ更新部25へ出力する。   The diagnostic data acquisition unit 23 acquires diagnostic data, which is a diagnostic data set, from the diagnosis target. The diagnostic data acquisition unit 23 outputs the acquired diagnostic data to the abnormality diagnosis unit 24 and the selection data update unit 25.

異常診断部24は、記憶装置10から単位データD2を取得すると共に、診断用データ取得部23から診断用データを取得し、取得した単位データD2及び診断用データに基づいて診断対象が異常か或いは正常かを診断する。   The abnormality diagnosis unit 24 acquires the unit data D2 from the storage device 10, acquires the diagnosis data from the diagnosis data acquisition unit 23, and the diagnosis target is abnormal based on the acquired unit data D2 and the diagnosis data. Diagnose whether it is normal.

具体的には、異常診断部24は、次の処理を実行する。すなわち、記憶装置10から単位データD2を取得し、この単位データD2から単位空間を生成する。診断用データ取得部23から診断用データを取得し、生成した単位空間を利用してMD値を演算する。MD値が閾値よりも小さい場合、プラントの状態が正常であると診断する一方、MD値が当該閾値以上の場合、診断対象の状態が異常であると診断する。閾値は、予め設定されて記憶された値である。なお、このような単位データD2及び診断用データに基づく異常診断の手法としては、公知手法を用いることができる。   Specifically, the abnormality diagnosis unit 24 executes the following process. That is, the unit data D2 is acquired from the storage device 10, and a unit space is generated from the unit data D2. The diagnostic data is acquired from the diagnostic data acquisition unit 23, and the MD value is calculated using the generated unit space. When the MD value is smaller than the threshold value, the state of the plant is diagnosed as normal. On the other hand, when the MD value is equal to or larger than the threshold value, the diagnosis target state is diagnosed as abnormal. The threshold value is a value set and stored in advance. In addition, a well-known method can be used as a method of abnormality diagnosis based on such unit data D2 and diagnostic data.

異常診断部24は、診断結果を選定用データ更新部25へ出力する。また、異常診断部24は、診断結果を外部のディスプレイやスピーカ等の報知装置(不図示)へ出力して報知する。   The abnormality diagnosis unit 24 outputs the diagnosis result to the selection data update unit 25. In addition, the abnormality diagnosis unit 24 outputs and notifies the diagnosis result to a notification device (not shown) such as an external display or a speaker.

選定用データ更新部25は、診断開始時から少なくとも1度は単位空間が生成されている(異常診断部24において単位空間を生成する際、当該生成が初回以外となる)場合、選定用データD1を更新する。具体的には、記憶装置10から単位データD2を取得し、この単位データD2を選定用データD1として複製する。そして、異常診断部24で診断対象が正常と診断されている場合、複製してなる選定用データD1に診断用データを追加して記憶装置10へ出力する一方、異常診断部24で診断対象が正常と診断されていない場合、複製してなる選定用データD1をそのままを記憶装置10へ出力する。   When the unit space is generated at least once from the start of diagnosis (when the unit space is generated by the abnormality diagnosis unit 24, the selection data update unit 25 generates data other than the first time). Update. Specifically, the unit data D2 is acquired from the storage device 10, and this unit data D2 is copied as selection data D1. When the abnormality diagnosis unit 24 determines that the diagnosis target is normal, the diagnosis data is added to the duplicated selection data D1 and output to the storage device 10, while the abnormality diagnosis unit 24 selects the diagnosis target. If it is not diagnosed as normal, the selection data D1 that is duplicated is output to the storage device 10 as it is.

次に、異常診断装置100により実施される異常診断方法について説明する。   Next, an abnormality diagnosis method performed by the abnormality diagnosis apparatus 100 will be described.

図5は、異常診断処理を示すフローチャートである。図5に示すように、異常診断装置100では、CPU20において次の異常診断処理を所定時間周期で繰り返し実行し、診断対象の異常を所定時間周期で診断する。すなわち、まず、単位データ決定処理を実行して単位データD2を決定し、当該単位データD2を異常診断部24で取得する(S1)。診断対象から出力された診断用データを、診断用データ取得部23を介して異常診断部24で取得する(S2:データ取得工程)。   FIG. 5 is a flowchart showing the abnormality diagnosis process. As shown in FIG. 5, in the abnormality diagnosis apparatus 100, the CPU 20 repeatedly performs the next abnormality diagnosis process at a predetermined time period, and diagnoses the abnormality to be diagnosed at a predetermined time period. That is, first, unit data determination processing is executed to determine unit data D2, and the unit diagnosis D24 is acquired by the abnormality diagnosis unit 24 (S1). The diagnosis data output from the diagnosis target is acquired by the abnormality diagnosis unit 24 via the diagnosis data acquisition unit 23 (S2: data acquisition step).

続いて、異常診断部24により、単位データD2に基づいて単位空間を生成する(S3)。具体的には、単位データD2の複数項目それぞれにおいて、平均値及び標準偏差を求めると共に、パラメータ値と平均値と標準偏差とから規準化値を求めて、相関行列を生成する。生成した相関行列の逆行列を求めて、単位空間を生成する。そして、異常診断部24により、単位空間を利用して、診断用データに基づきMD値を算出する(S4)。   Subsequently, the abnormality diagnosis unit 24 generates a unit space based on the unit data D2 (S3). Specifically, in each of the plurality of items of the unit data D2, an average value and a standard deviation are obtained, and a normalized value is obtained from the parameter value, the average value, and the standard deviation, and a correlation matrix is generated. A unit space is generated by obtaining an inverse matrix of the generated correlation matrix. Then, the abnormality diagnosis unit 24 calculates the MD value based on the diagnostic data using the unit space (S4).

続いて、算出したMD値が閾値よりも小さいか否かを判定する(S5)。MD値が閾値よりも小さい場合、診断用データは正常時データであり、診断対象の状態が正常であると診断する(S6)。一方、演算したMD値が閾値以上の場合、診断用データは異常時データであり、診断対象の状態が異常であると診断する(S7)。なお、上記S3〜上記S7の処理が診断工程を構成する。   Subsequently, it is determined whether or not the calculated MD value is smaller than a threshold value (S5). When the MD value is smaller than the threshold value, the diagnosis data is normal data, and it is diagnosed that the state of the diagnosis target is normal (S6). On the other hand, when the calculated MD value is greater than or equal to the threshold value, the diagnosis data is abnormal data, and the diagnosis target state is diagnosed as abnormal (S7). In addition, the process of said S3-said S7 comprises a diagnostic process.

ここで、本実施形態の異常診断装置100においては、上述したように、診断対象の異常を診断するに当たり、単位空間の生成前に単位データ決定処理を実行する。以下、単位データ決定処理について詳説する。   Here, in the abnormality diagnosis apparatus 100 of the present embodiment, as described above, when diagnosing abnormality of the diagnosis target, unit data determination processing is executed before generation of the unit space. Hereinafter, the unit data determination process will be described in detail.

図6は、単位データ決定処理を示すフローチャートである。図6に示すように、単位データ決定処理では、まず、単位空間の生成が初回か否かを判定する(S11)。ここでは、例えば、上記S3の処理の実行に初回判定フラグを持たせておき、この初回判定フラグが立っているか否かによって判定することができる。   FIG. 6 is a flowchart showing the unit data determination process. As shown in FIG. 6, in the unit data determination process, it is first determined whether or not the unit space is generated for the first time (S11). Here, for example, an initial determination flag is provided in the execution of the process of S3, and the determination can be made based on whether or not the initial determination flag is set.

上記S11でYESの場合、選定用データ除外部22により、選定用データD1から外れ値を除外する選定用データ除外処理を実行する(S12:データ除外工程)。続いて、単位データ選定部21により、選定用データD1から単位データD2を選定する単位データ選定処理を実行する(S13:データ選定工程)。   In the case of YES in S11, the selection data exclusion unit 22 executes selection data exclusion processing for excluding outliers from the selection data D1 (S12: data exclusion process). Subsequently, the unit data selection unit 21 executes unit data selection processing for selecting the unit data D2 from the selection data D1 (S13: data selection step).

上記S11でNOの場合、選定用データ更新部25により、選定用データD1を更新する選定用データ更新処理を実行する(S14:データ更新工程)。続いて、上記S14と同様に、単位データ選定部21により、選定用データD1から単位データD2を選定する単位データ選定処理を実行する(S15:データ選定工程)。   In the case of NO in S11, the selection data update unit 25 executes selection data update processing for updating the selection data D1 (S14: data update process). Subsequently, similarly to S14, the unit data selection unit 21 executes unit data selection processing for selecting the unit data D2 from the selection data D1 (S15: data selection step).

上記S13又はS15の後、選定した単位データD2に関する情報を、データ修正ユニット30の表示部30bに表示する。そして、例えば単位データD2の項目別の分散や過去の経験等から、選定者によってデータ修正ユニット30の入力部30aを介して操作入力がなされた場合、その操作入力に応じて単位データD2を修正する(S16)。以上により、単位データD2が決定され、この単位データD2が記憶装置10に出力されて記憶されることとなる。   After S13 or S15, information on the selected unit data D2 is displayed on the display unit 30b of the data correction unit 30. Then, for example, when an operation input is made by the selector through the input unit 30a of the data correction unit 30 based on, for example, the distribution of the unit data D2 by item or past experience, the unit data D2 is corrected according to the operation input. (S16). Thus, the unit data D2 is determined, and this unit data D2 is output to the storage device 10 and stored.

図7(a)は、選定用データ除外処理を示すフローチャートである。図7(a)に示すように、選定用データ除外処理においては、記憶装置10から選定用データD1を取得し、選定用データD1に含まれる複数のデータセットの複数項目別に、パラメータ値の平均値Av、標準偏差σ、尖度γ及び第2範囲R2を含む統計量を算出する(S21)。何れかの項目で尖度γが所定値以上であるか否かに基づいて、選定用データD1に含まれる複数のデータセットのパラメータ値に外れ値が存在する否かを判定する(S22)。 FIG. 7A is a flowchart showing the selection data excluding process. As shown in FIG. 7A, in the selection data excluding process, the selection data D1 is acquired from the storage device 10, and the average of the parameter values is obtained for a plurality of items of a plurality of data sets included in the selection data D1. Statistics including value Av 2 , standard deviation σ 2 , kurtosis γ, and second range R 2 are calculated (S 21). Based on whether or not the kurtosis γ is greater than or equal to a predetermined value in any item, it is determined whether or not there is an outlier in the parameter values of a plurality of data sets included in the selection data D1 (S22).

尖度γが所定値以上ではない場合、選定用データD1に含まれる複数のデータセットのパラメータ値に外れ値が無いとして、そのまま処理を終了する。一方、尖度γが所定値以上の場合、外れ値を有するデータセットを選定用データD1から除外する(S23)。具体的には、第2範囲R2から外れたパラメータ値を外れ値として、選定用データD1に含まれているデータセットのうち当該外れ値を有するデータセットを、選定用データD1から除外する。外れ値を除外した後、再度統計量を算出し(上記S21)、尖度γに基づき外れ値が存在するか否かを判定し(上記S22)、外れ値を有するデータセットを選定用データD1から除外する(上記S23)という処理を、外れ値が存在しなくなるまで(上記S22の判定処理がNOとなるまで)繰り返し行う。そして、選定用データD1を記憶装置10へ出力する。   If the kurtosis γ is not greater than or equal to the predetermined value, it is determined that there are no outliers in the parameter values of the plurality of data sets included in the selection data D1, and the process is terminated. On the other hand, when the kurtosis γ is equal to or larger than the predetermined value, the data set having an outlier is excluded from the selection data D1 (S23). Specifically, a parameter value out of the second range R2 is regarded as an outlier, and a data set having the outlier among the data sets included in the selection data D1 is excluded from the selection data D1. After the outlier is excluded, the statistic is calculated again (S21), it is determined whether or not the outlier exists based on the kurtosis γ (S22), and the data set having the outlier is selected as the selection data D1. Is repeatedly performed until there is no outlier (the determination process in S22 is NO). Then, the selection data D1 is output to the storage device 10.

図7(b)は、単位データ選定処理を示すフローチャートである。図7(b)に示すように、単位データ選定処理においては、記憶装置10から選定用データD1を取得し、選定用データD1に含まれる複数のデータセットの複数項目別に、パラメータ値の平均値Av、標準偏差σ、第1範囲R1を含む統計量を算出する(S31)。選定用データD1に含まれる複数のデータセットのうち、複数項目の全てにおいてパラメータ値が第1範囲R1内に存在するデータセット(すなわち、第1範囲R1から外れたパラメータ値を有するデータセット以外)を、単位データD2として選定する(S32)。 FIG. 7B is a flowchart showing unit data selection processing. As shown in FIG. 7B, in the unit data selection process, the selection data D1 is acquired from the storage device 10, and the average value of the parameter values for each of the plurality of items of the plurality of data sets included in the selection data D1. A statistic including Av 1 , standard deviation σ 1 , and first range R1 is calculated (S31). Among the plurality of data sets included in the selection data D1, the data sets in which the parameter values exist in the first range R1 in all of the plurality of items (that is, other than the data sets having parameter values outside the first range R1). Is selected as unit data D2 (S32).

図8は、選定用データ更新処理を示すフローチャートである。図8に示すように、選定用データ更新処理においては、記憶装置10から単位データD2を取得し、この単位データD2を選定用データD1として複製する(S41)。続いて、異常診断部24による診断結果が正常か否か(上記S5の判定処理がYESか否か)を判定する(S42)。診断結果が正常の場合、当該診断に用いられた診断用データ(上記S2で取得した診断用データ)を選定用データD1に追加し、この選定用データD1を記憶装置10へ出力する(S43)。一方、診断結果が正常ではない場合、診断用データを選定用データD1に追加せずに、選定用データD1をそのまま記憶装置10へ出力する。   FIG. 8 is a flowchart showing the selection data update process. As shown in FIG. 8, in the selection data update process, the unit data D2 is acquired from the storage device 10, and this unit data D2 is duplicated as the selection data D1 (S41). Subsequently, it is determined whether or not the diagnosis result by the abnormality diagnosis unit 24 is normal (whether the determination process of S5 is YES) (S42). When the diagnosis result is normal, the diagnosis data (diagnosis data acquired in S2) used for the diagnosis is added to the selection data D1, and the selection data D1 is output to the storage device 10 (S43). . On the other hand, when the diagnosis result is not normal, the selection data D1 is output to the storage device 10 as it is without adding the diagnosis data to the selection data D1.

次に、上述した一連の異常診断に係る処理をCPU(コンピュータ)20に実行せしめる異常診断プログラムを説明する。   Next, an abnormality diagnosis program that causes the CPU (computer) 20 to execute the above-described series of abnormality diagnosis processes will be described.

図9は、異常診断プログラムを説明するための図である。図9に示すように、異常診断プログラムPは、CPU20に挿入されてアクセスされるものであり、記憶装置10に形成されたプログラム格納領域12内に格納されている。異常診断プログラムPは、単位データ選定モジュール51と、選定用データ除外モジュール52と、異常診断モジュール53と、診断用データ取得モジュール54と、選定用データ更新モジュール55と、を備えて構成される。   FIG. 9 is a diagram for explaining the abnormality diagnosis program. As shown in FIG. 9, the abnormality diagnosis program P is inserted into the CPU 20 and accessed, and is stored in the program storage area 12 formed in the storage device 10. The abnormality diagnosis program P includes a unit data selection module 51, a selection data exclusion module 52, an abnormality diagnosis module 53, a diagnosis data acquisition module 54, and a selection data update module 55.

単位データ選定モジュール51により実現される機能は、上述した単位データ選定部21の機能と同様である。選定用データ除外モジュール52により実現される機能は、上述した選定用データ除外部22の機能と同様である。異常診断モジュール53により実現される機能は、上述した異常診断部24の機能と同様である。診断用データ取得モジュール54により実現される機能は、上述した診断用データ取得部23の機能と同様である。選定用データ更新モジュール55により実現される機能は、上述した選定用データ更新部25の機能と同様である。   The functions realized by the unit data selection module 51 are the same as the functions of the unit data selection unit 21 described above. The function realized by the selection data exclusion module 52 is the same as the function of the selection data exclusion unit 22 described above. The function realized by the abnormality diagnosis module 53 is the same as the function of the abnormality diagnosis unit 24 described above. The function realized by the diagnostic data acquisition module 54 is the same as the function of the diagnostic data acquisition unit 23 described above. The function realized by the selection data update module 55 is the same as the function of the selection data update unit 25 described above.

以上、本実施形態に係る異常診断装置100、異常診断方法、及び異常診断プログラムPでは、選定用データD1に含まれる複数のデータセットにおいて、複数項目別に、パラメータ値の標準偏差σに基づく第1範囲R1が算出される。そして、選定用データD1に含まれる複数のデータセットのうち、複数項目の全てにおいてパラメータ値が第1範囲R1内に存在するデータセットが、単位データD2として選定される。 As described above, in the abnormality diagnosis device 100, the abnormality diagnosis method, and the abnormality diagnosis program P according to the present embodiment, the first based on the standard deviation σ 1 of the parameter value for each of a plurality of items in a plurality of data sets included in the selection data D1. One range R1 is calculated. Then, among the plurality of data sets included in the selection data D1, a data set having parameter values within the first range R1 in all of the plurality of items is selected as the unit data D2.

これにより、単位データD2を選定者(人)が選定する必要が無くなり、選定者による判断の介入を低減して、且つ、各パラメータ値毎に正常範囲にあるかどうかを判定することもなく、単位データD2を選定することができる。したがって、単位データD2の選定を自動的に選定可能なレベルまで定式化し、当該選定の安定化を図ることができ、選定者による依存を低減して簡易に単位データD2を選定することが可能となる。その結果、人の経験によるところが大きくて手間がかかるという問題(例えば、選定者次第で異常データが単位データD2に混ざってしまうことや、項目が多いと各項目別にパラメータ値を確認するのは難しい等)を解消することが可能となる。   This eliminates the need for the selector (person) to select the unit data D2, reduces the judgment intervention by the selector, and determines whether each parameter value is in the normal range, Unit data D2 can be selected. Therefore, the selection of the unit data D2 can be automatically formulated to a selectable level, the selection can be stabilized, and the unit data D2 can be easily selected with less dependence on the selector. Become. As a result, there is a problem that it depends on human experience and is troublesome (for example, depending on the selected person, abnormal data is mixed with the unit data D2, or if there are many items, it is difficult to check the parameter value for each item. Etc.) can be eliminated.

本実施形態では、選定用データD1に含まれる複数のデータセットにおいて、複数項目別に、パラメータ値の尖度γとパラメータ値の標準偏差σに基づく第2範囲R2とが算出される。尖度γが所定値以上の場合、選定用データD1に含まれる複数のデータセットから、第2範囲R2から外れたパラメータ値を有するデータセットが除外される。これにより、外れ値を有するデータセットを選定用データD1から除外でき、選定用データD1の正規性を高めることができる。 In the present embodiment, in a plurality of data sets included in the selection data D1, the kurtosis γ of the parameter value and the second range R2 based on the standard deviation σ 2 of the parameter value are calculated for each of a plurality of items. When the kurtosis γ is equal to or greater than a predetermined value, a data set having a parameter value outside the second range R2 is excluded from a plurality of data sets included in the selection data D1. Thereby, a data set having an outlier can be excluded from the selection data D1, and the normality of the selection data D1 can be improved.

特に本実施形態では、選定用データD1に含まれる数百以上の項目のパラメータ値の各分布において、平均値Avからm・σに亘る第1範囲R1に入るものが「正常」と判断され、自動的に単位データD2として選定される。また、平均値Avからk・σに亘る第2範囲R2に入らない外れ値を含むデータセットが、選定用データD1から除外される。よって、数百以上の項目数を含む選定用データD1からでも、指定値m,kを設定するだけで、選定者による依存を低減して単位データD2を自動的に得ることができる。 In particular, in the present embodiment, in each distribution of parameter values of several hundred or more items included in the selection data D1, those that fall within the first range R1 ranging from the average value Av 1 to m · σ 1 are determined as “normal”. And automatically selected as unit data D2. The data set containing outliers not enter the second range R2 ranging from the average value Av 2 to k · sigma 2 is excluded from the selection data D1. Therefore, even from the selection data D1 including the number of items of several hundred or more, it is possible to automatically obtain the unit data D2 by reducing the dependence by the selector only by setting the designated values m and k.

本実施形態では、単位データD2が選定用データD1として複製され、診断対象が正常と診断されている場合、当該選定用データD1に診断用データが追加される。これにより、異常診断の都度、選定用データD1を自動更新して単位データD2を選定することができる。その結果、例えば、正常時データとみなされるパラメータ値が時間的に変化する場合や、単位空間の初回生成時に選定用データD1にデータセットが十分に含まれていない場合でも、異常診断の都度、単位データD2を適正化することができる。   In this embodiment, when the unit data D2 is replicated as the selection data D1, and the diagnosis target is diagnosed as normal, the diagnosis data is added to the selection data D1. As a result, the unit data D2 can be selected by automatically updating the selection data D1 at each abnormality diagnosis. As a result, for example, even when a parameter value regarded as normal data changes with time, or even when a data set is not sufficiently included in the selection data D1 when the unit space is initially generated, The unit data D2 can be optimized.

本実施形態では、データ修正ユニット30の入力部30aを介して修正の操作入力がなされた場合、その操作入力に応じて、上記S32で選定した単位データD2が修正される。これにより、選定者が実用に合わせて個別にチューニングする(項目別の分散や人の経験を考慮してデータを追加/削除修正する)ことができ、診断精度を向上させることが可能となる。   In this embodiment, when a correction operation input is made via the input unit 30a of the data correction unit 30, the unit data D2 selected in S32 is corrected according to the operation input. As a result, the selected person can individually tune according to practical use (adding / deleting / correcting data in consideration of item-by-item distribution and human experience), and the diagnostic accuracy can be improved.

以上、本発明の好適な実施形態について説明したが、本発明は上記実施形態に限られるものではなく、各請求項に記載した要旨を変更しない範囲で変形し、又は他のものに適用したものであってもよい。   The preferred embodiments of the present invention have been described above. However, the present invention is not limited to the above-described embodiments. The present invention is modified without departing from the scope described in the claims or applied to others. It may be.

例えば、上記実施形態において、第1範囲R1及び第2範囲R2は、パラメータ値の平均値Av,Avを中心に分布しているが、パラメータ値の最頻値を中心に分布してもよい。上記実施形態において、外れ値を含むデータセットを選定用データD1から除外する工程(選定用データ除外処理)は、人により行われてもよい。 For example, in the above-described embodiment, the first range R1 and the second range R2 are distributed around the average values Av 1 and Av 2 of the parameter values, but may be distributed around the mode value of the parameter values. Good. In the above embodiment, the step of excluding the data set including outliers from the selection data D1 (selection data excluding process) may be performed by a person.

上記実施形態において、単位データ選定部21、選定用データ除外部22、診断用データ取得部23、異常診断部24及び選定用データ更新部25は、1つのCPU20により構成されているが、これらの少なくとも何れかが別の情報処理装置により構成されていてもよいし、それぞれ別の情報処理装置により構成されてもよい。同様に、上記実施形態において、選定用データD1及び単位データD2は同じ記憶装置10に記憶されているが、互いに別の記憶装置に記憶されてもよい。   In the above embodiment, the unit data selection unit 21, the selection data exclusion unit 22, the diagnosis data acquisition unit 23, the abnormality diagnosis unit 24, and the selection data update unit 25 are configured by one CPU 20, but these At least one of them may be configured by different information processing apparatuses, or may be configured by different information processing apparatuses. Similarly, in the above embodiment, the selection data D1 and the unit data D2 are stored in the same storage device 10, but may be stored in different storage devices.

上記実施形態において、単位データ選定部21、選定用データ除外部22、診断用データ取得部23、異常診断部24、選定用データ更新部25、データ修正ユニット30、選定用データD1及び単位データD2における各間の入力又は出力は、有線で行ってもよいし無線で行ってもよい。   In the above embodiment, the unit data selection unit 21, the selection data exclusion unit 22, the diagnosis data acquisition unit 23, the abnormality diagnosis unit 24, the selection data update unit 25, the data correction unit 30, the selection data D1, and the unit data D2 The input or output between each may be performed by wire or wirelessly.

21 単位データ選定部
21a 第1算出部
21b 選定部
22 選定用データ除外部
22a 第2算出部
22b 除外部
23 診断用データ取得部
24 異常診断部
25 選定用データ更新部
100 異常診断装置
d データセット
D1 選定用データ
D2 単位データ
P 異常診断プログラム
R1 第1範囲
R2 第2範囲
21 unit data selection unit 21a first calculation unit 21b selection unit 22 selection data exclusion unit 22a second calculation unit 22b exclusion unit 23 diagnosis data acquisition unit 24 abnormality diagnosis unit 25 selection data update unit 100 abnormality diagnosis device d data set D1 Selection data D2 Unit data P Abnormality diagnosis program R1 First range R2 Second range

Claims (9)

診断対象に関する複数項目のパラメータ値を有するデータセットに基づき、前記パラメータ値が正規分布に従うことを前提にした診断手法を用いて、前記診断対象の異常を診断する異常診断方法であって、
複数の前記データセットを含む選定用データから、診断を行う際に基準となる単位データを選定するデータ選定工程を備え、
前記データ選定工程は、
前記選定用データに含まれる複数の前記データセットにおいて、前記複数項目別に、前記パラメータ値の標準偏差に基づく第1範囲を算出する工程と、
前記選定用データに含まれる複数の前記データセットのうち、前記複数項目の全てにおいて前記パラメータ値が前記第1範囲内に存在するデータセットを、前記単位データとして選定する工程と、を含む、異常診断方法。
Based on a data set having parameter values of a plurality of items related to a diagnosis object, using a diagnosis method on the assumption that the parameter value follows a normal distribution, an abnormality diagnosis method for diagnosing an abnormality of the diagnosis object,
From a selection data including a plurality of the data sets, comprising a data selection step of selecting unit data serving as a reference when performing diagnosis,
The data selection process includes
Calculating a first range based on a standard deviation of the parameter values for each of the plurality of items in the plurality of data sets included in the selection data;
Selecting, as the unit data, a data set in which the parameter value exists in the first range in all of the plurality of items among the plurality of data sets included in the selection data. Diagnosis method.
前記選定用データの一部を除外するデータ除外工程をさらに備え、
前記データ除外工程は、
前記選定用データに含まれる複数の前記データセットにおいて、前記複数項目別に、前記パラメータ値の尖度と前記パラメータ値の標準偏差に基づく第2範囲とを算出する工程と、
前記尖度が前記所定値以上の場合、前記選定用データに含まれる複数の前記データセットから、前記第2範囲から外れた前記パラメータ値を有するデータセットを除外する工程と、を含む、請求項1に記載の異常診断方法。
A data excluding step of excluding a part of the selection data;
The data exclusion step includes
Calculating a kurtosis of the parameter value and a second range based on a standard deviation of the parameter value for each of the plurality of items in the plurality of data sets included in the selection data;
Removing the data set having the parameter value out of the second range from the plurality of data sets included in the selection data when the kurtosis is equal to or greater than the predetermined value. The abnormality diagnosis method according to 1.
前記診断対象から診断用の前記データセットである診断用データを取得するデータ取得工程と、
前記単位データ及び前記診断用データに基づき、前記診断対象が異常か或いは正常かを診断する診断工程と、
前記単位データを前記選定用データとして複製し、前記診断工程で前記診断対象が正常と診断されている場合に当該選定用データに前記診断用データを追加するデータ更新工程と、を備える、請求項1又は2に記載の異常診断方法。
A data acquisition step of acquiring diagnostic data which is the data set for diagnosis from the diagnosis target;
Based on the unit data and the diagnostic data, a diagnostic step of diagnosing whether the diagnostic object is abnormal or normal;
A data update step of replicating the unit data as the selection data and adding the diagnostic data to the selection data when the diagnosis target is diagnosed as normal in the diagnosis step. 3. The abnormality diagnosis method according to 1 or 2.
前記診断手法は、マハラノビスタグチメソッドであり、
前記単位データは、マハラノビス距離を算出する際に基準となる単位空間を構成する、請求項1〜3の何れか一項に記載の異常診断方法。
The diagnostic method is a Mahalanobis Taguchi method,
The abnormality diagnosis method according to claim 1, wherein the unit data constitutes a unit space that serves as a reference when calculating the Mahalanobis distance.
診断対象に関する複数項目のパラメータ値を有するデータセットに基づき、前記パラメータ値が正規分布に従うことを前提にした診断手法を用いて、前記診断対象の異常を診断する異常診断装置であって、
複数の前記データセットを含む選定用データから、診断を行う際に基準となる単位データを選定する単位データ選定部を備え、
前記単位データ選定部は、
前記選定用データに含まれる複数の前記データセットにおいて、前記複数項目別に、前記パラメータ値の標準偏差に基づく第1範囲を算出する第1算出部と、
前記選定用データに含まれる複数の前記データセットのうち、前記複数項目の全てにおいて前記パラメータ値が前記第1範囲内に存在するデータセットを、前記単位データとして選定する選定部と、を含む、異常診断装置。
Based on a data set having a plurality of parameter values related to a diagnosis target, using a diagnostic method based on the assumption that the parameter values follow a normal distribution, an abnormality diagnosis apparatus for diagnosing the abnormality of the diagnosis target,
A unit data selection unit that selects unit data serving as a reference when performing diagnosis from data for selection including a plurality of the data sets,
The unit data selection unit
A first calculator that calculates a first range based on a standard deviation of the parameter value for each of the plurality of items in the plurality of data sets included in the selection data;
A selection unit that selects, as the unit data, a data set in which the parameter value exists in the first range in all of the plurality of data sets included in the selection data; Abnormality diagnosis device.
前記選定用データの一部を除外する選定用データ除外部をさらに備え、
前記選定用データ除外部は、
前記選定用データに含まれる複数の前記データセットにおいて、前記複数項目別に、前記パラメータ値の尖度と前記パラメータ値の標準偏差に基づく第2範囲とを算出する第2算出部と、
前記尖度が前記所定値以上の場合、前記選定用データに含まれる複数の前記データセットから、前記第2範囲から外れた前記パラメータ値を有するデータセットを除外する除外部と、を含む、請求項5に記載の異常診断装置。
A selection data exclusion unit for excluding a part of the selection data;
The selection data excluding unit is
A second calculator that calculates a kurtosis of the parameter value and a second range based on a standard deviation of the parameter value for each of the plurality of data sets included in the selection data;
An exclusion unit that excludes a data set having the parameter value out of the second range from the plurality of data sets included in the selection data when the kurtosis is equal to or greater than the predetermined value. Item 6. The abnormality diagnosis device according to Item 5.
前記診断対象から診断用の前記データセットである診断用データを取得する診断用データ取得部と、
前記単位データ及び前記診断用データに基づき、前記診断対象が異常か或いは正常かを診断する異常診断部と、
前記単位データを前記選定用データとして複製し、前記異常診断部で前記診断対象が正常と診断されている場合に当該選定用データに前記診断用データを追加する選定用データ更新部と、を備える、請求項5又は6に記載の異常診断装置。
A diagnostic data acquisition unit for acquiring diagnostic data which is the data set for diagnosis from the diagnostic target;
Based on the unit data and the diagnostic data, an abnormality diagnosis unit that diagnoses whether the diagnosis target is abnormal or normal;
A selection data update unit that replicates the unit data as the selection data and adds the diagnosis data to the selection data when the abnormality diagnosis unit diagnoses the diagnosis target as normal. The abnormality diagnosis device according to claim 5 or 6.
前記診断手法は、マハラノビスタグチメソッドであり、
前記単位データは、マハラノビス距離を算出する際に基準となる単位空間を構成する、請求項5〜7の何れか一項に記載の異常診断装置。
The diagnostic method is a Mahalanobis Taguchi method,
The abnormality diagnosis apparatus according to claim 5, wherein the unit data constitutes a unit space that serves as a reference when calculating the Mahalanobis distance.
コンピュータに、請求項1〜4の何れか一項記載の異常診断方法を実行せしめる、異常診断プログラム。   An abnormality diagnosis program for causing a computer to execute the abnormality diagnosis method according to claim 1.
JP2014227155A 2014-11-07 2014-11-07 Abnormality diagnosis method, abnormality diagnosis device, and abnormality diagnosis program Active JP6492555B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014227155A JP6492555B2 (en) 2014-11-07 2014-11-07 Abnormality diagnosis method, abnormality diagnosis device, and abnormality diagnosis program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014227155A JP6492555B2 (en) 2014-11-07 2014-11-07 Abnormality diagnosis method, abnormality diagnosis device, and abnormality diagnosis program

Publications (2)

Publication Number Publication Date
JP2016091417A true JP2016091417A (en) 2016-05-23
JP6492555B2 JP6492555B2 (en) 2019-04-03

Family

ID=56016204

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014227155A Active JP6492555B2 (en) 2014-11-07 2014-11-07 Abnormality diagnosis method, abnormality diagnosis device, and abnormality diagnosis program

Country Status (1)

Country Link
JP (1) JP6492555B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017217069A1 (en) * 2016-06-13 2017-12-21 株式会社日立製作所 Rotor diagnostic apparatus, rotor diagnostic method, and rotor diagnostic program
JP2019178625A (en) * 2018-03-30 2019-10-17 三菱重工業株式会社 Abnormality diagnosis system of pump equipment and abnormality diagnosis method of pump equipment
WO2020194852A1 (en) * 2019-03-27 2020-10-01 株式会社島津製作所 Pump monitoring device, vacuum pump, and product-accumulation diagnosis data processing program
DE112020003659T5 (en) 2019-08-01 2022-04-21 Mitsubishi Electric Corporation METHOD FOR DIAGNOSIS OF ANOMALIES, DEVICE FOR DIAGNOSIS OF ANOMALIES, AND PROGRAM FOR DIAGNOSIS OF ANOMALIES
JP7367366B2 (en) 2019-07-23 2023-10-24 オムロン株式会社 Anomaly detection device, anomaly detection method, and anomaly detection program

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011106467A (en) * 2011-02-28 2011-06-02 Hitachi Ltd Method for displaying display screen in system for diagnosing performance of gas turbine
JP2012008030A (en) * 2010-06-25 2012-01-12 Toshiba Plant Systems & Services Corp Rotator bearing diagnostic device
JP2012252556A (en) * 2011-06-03 2012-12-20 Mitsubishi Heavy Ind Ltd Plant monitoring device, plant monitoring program and plant monitoring method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012008030A (en) * 2010-06-25 2012-01-12 Toshiba Plant Systems & Services Corp Rotator bearing diagnostic device
JP2011106467A (en) * 2011-02-28 2011-06-02 Hitachi Ltd Method for displaying display screen in system for diagnosing performance of gas turbine
JP2012252556A (en) * 2011-06-03 2012-12-20 Mitsubishi Heavy Ind Ltd Plant monitoring device, plant monitoring program and plant monitoring method

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017217069A1 (en) * 2016-06-13 2017-12-21 株式会社日立製作所 Rotor diagnostic apparatus, rotor diagnostic method, and rotor diagnostic program
JPWO2017217069A1 (en) * 2016-06-13 2019-03-14 株式会社日立製作所 Rotating machine diagnostic device, rotating machine diagnostic method, and rotating machine diagnostic program
US11209487B2 (en) 2016-06-13 2021-12-28 Hitachi, Ltd. Rotor diagnostic apparatus, rotor diagnostic method, and rotor diagnostic program
JP2019178625A (en) * 2018-03-30 2019-10-17 三菱重工業株式会社 Abnormality diagnosis system of pump equipment and abnormality diagnosis method of pump equipment
WO2020194852A1 (en) * 2019-03-27 2020-10-01 株式会社島津製作所 Pump monitoring device, vacuum pump, and product-accumulation diagnosis data processing program
CN113631817A (en) * 2019-03-27 2021-11-09 株式会社岛津制作所 Pump monitoring device, vacuum pump, and data processing program for diagnosing product accumulation
JPWO2020194852A1 (en) * 2019-03-27 2021-12-02 株式会社島津製作所 Data processing program for pump monitoring equipment, vacuum pumps and product deposition diagnostics
JP7188560B2 (en) 2019-03-27 2022-12-13 株式会社島津製作所 Data processing program for pump monitor, vacuum pump and product deposition diagnostics
CN113631817B (en) * 2019-03-27 2024-03-08 株式会社岛津制作所 Pump monitoring device, vacuum pump, and recording medium
JP7367366B2 (en) 2019-07-23 2023-10-24 オムロン株式会社 Anomaly detection device, anomaly detection method, and anomaly detection program
DE112020003659T5 (en) 2019-08-01 2022-04-21 Mitsubishi Electric Corporation METHOD FOR DIAGNOSIS OF ANOMALIES, DEVICE FOR DIAGNOSIS OF ANOMALIES, AND PROGRAM FOR DIAGNOSIS OF ANOMALIES

Also Published As

Publication number Publication date
JP6492555B2 (en) 2019-04-03

Similar Documents

Publication Publication Date Title
JP6492555B2 (en) Abnormality diagnosis method, abnormality diagnosis device, and abnormality diagnosis program
JP5765336B2 (en) Fault analysis apparatus, fault analysis method and program
US20170261403A1 (en) Abnormality detection procedure development apparatus and abnormality detection procedure development method
US20110288846A1 (en) Technique and tool for efficient testing of controllers in development (h-act project)
US11892819B2 (en) Control device, control system, control method, and computer-readable storage medium
US11748321B2 (en) Time-series data condensation and graphical signature analysis
JP2013041491A (en) Abnormality diagnostic device
JP2016091043A (en) Diagnostic job generating system, diagnostic job generating method, and diagnostic job generation displaying method
JPWO2020148838A1 (en) Estimator, estimation method, and program
CN111143101B (en) Method, device, storage medium and electronic equipment for determining fault source
JP2019125206A (en) Information processing device, information processing method, and program
JP6310865B2 (en) Source code evaluation system and method
JP5948998B2 (en) Abnormality diagnosis device
JP6625183B1 (en) Remodeling support device and remodeling support method
JP6320269B2 (en) Software test support apparatus and software test support program
JP6347771B2 (en) Abnormality diagnosis apparatus, abnormality diagnosis method, and abnormality diagnosis program
JP2018018267A (en) Software development support apparatus, software development support method, and program
JP2013041492A (en) Abnormality diagnostic device
Staron et al. Industrial self-healing measurement systems
JP2015228269A (en) Abnormality diagnostic device
JP2020098538A (en) Control device, control method and control program
JP7073432B2 (en) Impact explanation variable identification device, impact explanation variable identification method and impact explanation variable identification program
JP7391765B2 (en) Plant monitoring support equipment, methods and programs
WO2021124861A1 (en) Information processing device, information processing method, and information processing program
TWI837371B (en) Support device, display device, support method and support program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170926

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180706

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180724

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180921

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: 20190205

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190218

R151 Written notification of patent or utility model registration

Ref document number: 6492555

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151