JP7647138B2 - Abnormality diagnosis system and abnormality diagnosis method, frequency fluctuation correction processing device and correction processing method, abnormality diagnosis program, and frequency fluctuation correction processing program - Google Patents
Abnormality diagnosis system and abnormality diagnosis method, frequency fluctuation correction processing device and correction processing method, abnormality diagnosis program, and frequency fluctuation correction processing program Download PDFInfo
- Publication number
- JP7647138B2 JP7647138B2 JP2021018234A JP2021018234A JP7647138B2 JP 7647138 B2 JP7647138 B2 JP 7647138B2 JP 2021018234 A JP2021018234 A JP 2021018234A JP 2021018234 A JP2021018234 A JP 2021018234A JP 7647138 B2 JP7647138 B2 JP 7647138B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- moving average
- learning
- acceleration data
- abnormality
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
本発明は、診断対象の稼働状態を逐次学習し、学習結果に基づき異常を診断する技術に関する。 The present invention relates to a technology that sequentially learns the operating state of a diagnostic target and diagnoses abnormalities based on the learning results.
回転機などの機器の異常診断は、診断対象の機器に設置された加速度センサの加速度データの定Q変換などに基づく診断により行われている(特許文献1,2参照)。
Abnormality diagnosis of equipment such as rotating machines is performed based on constant Q conversion of acceleration data from an acceleration sensor installed in the equipment to be diagnosed (see
基本的な診断手法として診断対象が正常に動作し続けている場合には異常時の加速度データが少ないことから、正常時の加速度データのみを学習して現在の加速度データが正常時と比較してどれだけ異なっているかを数値化することで、診断対象の正常/異常を診断している。 As a basic diagnostic technique, when the subject to be diagnosed continues to operate normally, there is little acceleration data during abnormal times, so the normal/abnormal state of the subject is diagnosed by learning only the acceleration data during normal times and quantifying how much the current acceleration data differs from normal times.
(1)従来の異常診断の手法は、以下に示す事情により診断の精度が低下するおそれがあった。 (1) Conventional methods for diagnosing abnormalities have the potential to reduce diagnostic accuracy due to the following reasons:
すなわち、前述のように診断対象が安定稼働している正常動作時の状態を学習するため、診断対象の稼働停止時(以下、「過渡期」と記述する。)の不安定な加速度が異常と判定される場合があった。 In other words, as mentioned above, because the system learns the state of the diagnostic target during normal operation when the diagnostic target is running stably, unstable acceleration when the diagnostic target is not operating (hereinafter referred to as the "transition period") may be judged to be abnormal.
また、加速度センサの設置直後や、診断対象の修繕を行った直後などは加速度センサの出力が安定しない場合もあった。例えば加速度センサは、マグネット吸着や接着剤、ボルト締めなどで診断対象に設置されるが、設置状態の馴染むまでセンサ出力が安定しない。また、診断対象の修繕を行った場合、対象機器を土台にボルト締めすることが多く、ボルト締め具合などが馴染むまでセンサ出力が安定しない。 In addition, the output of the acceleration sensor may not be stable immediately after the acceleration sensor is installed or immediately after the diagnosis target is repaired. For example, the acceleration sensor is installed on the diagnosis target using magnetic attraction, adhesives, bolts, etc., but the sensor output does not stabilize until the installation condition is settled. In addition, when the diagnosis target is repaired, the target device is often bolted to the base, and the sensor output does not stabilize until the bolt tightening is settled.
このようなセンサ出力が安定しない過渡期の加速度データは、安定稼働中に比べて加速度が大きかったり小さかったりするので、その期間のデータを学習すると診断精度が低下してしまう。 Acceleration data from such transitional periods when the sensor output is unstable can be larger or smaller than during stable operation, so learning data from those periods will result in a decrease in diagnostic accuracy.
さらに図1に示すように、ある一定の期間(例えば1週間または1か月)の診断対象の加速度データを学習すると、その期間を基準に年変動が発生する。この年変動の影響により、診断対象が正常である場合でも閾値を超えて異常が誤検出される場合がある。このような年変動は、温度により診断対象の膨張や収縮を招き、ポンプなどの場合には水に溶け込む空気量が変化することが影響しているものと考えられている。 Furthermore, as shown in Figure 1, when the acceleration data of the diagnostic target for a certain period (for example, one week or one month) is learned, annual fluctuations occur based on that period. Due to the influence of this annual fluctuation, even if the diagnostic target is normal, the threshold may be exceeded and an abnormality may be falsely detected. This annual fluctuation is thought to be caused by the expansion and contraction of the diagnostic target due to temperature, and in the case of pumps, the amount of air dissolved in water changes.
(2)そこで、加速度センサによる加速度データを学習・診断する際、加速度データを変換した周波数データをオクターブ周波数毎に区切り、区切られた区間毎に平均値などの特徴量を算出し、算出された特徴量から過渡期と判定された区間を学習・診断の対象から除外する手法が提案されている。 (2) Therefore, a method has been proposed in which, when learning and diagnosing acceleration data obtained by an acceleration sensor, the frequency data converted from the acceleration data is divided into octave frequencies, feature values such as average values are calculated for each divided section, and sections determined to be transitional periods based on the calculated feature values are excluded from the learning and diagnosing targets.
このような手法は、図2(a)に示すように、回転機1の土台3が強固であれば問題は生じない。ところが、図2(b)に示すように、不安定な土台3の上に回転機1を設置した場合、加速度センサ2による加速度データは回転機1の振動だけでなく、土台3の振動にも影響されてしまう。
As shown in Figure 2(a), this method does not cause any problems if the
その結果、図3に示すように、土台3の振動により周波数が揺らいでしまって、学習・診断の対象データ(学習データ・周波数データ)に回転機1の本来の周波数が現れないおそれがある。
As a result, as shown in Figure 3, the vibration of the
また、前記揺らぎに起因して前記区切りの位置では、周波数成分がどちらの区間に含まるかによって学習データ・診断データが大きく変わるおそれがある。例えば回転周波数の大きな成分が区切りの位置となった場合には、どちらの区間に含まるかによって前記特徴量が大きく異なってしまう。 Furthermore, due to the fluctuations, there is a risk that the learning data and diagnostic data will change significantly depending on which section the frequency component falls in at the position of the separator. For example, if a component with a large rotation frequency falls at the position of the separator, the feature amount will differ significantly depending on which section it falls in.
本発明は、このような従来の問題を解決するためになされ、加速度データの振幅の揺らぎの影響を抑えて診断精度の向上を図ることを解決課題としている。 The present invention was made to solve these conventional problems, and aims to improve diagnostic accuracy by suppressing the effects of fluctuations in the amplitude of acceleration data.
(1)本発明の一態様は、診断対象の稼働状態を逐次学習して異常を診断する異常診断システムであって、
前記診断対象の加速度データから得られた加速度振幅の移動平均に基づき揺らぎ補正を行う補正処理部と、
前記揺らぎ補正の結果を逐次学習して学習結果モデルを生成する逐次学習処理部と、
前記学習結果モデルを用いて診断対象の現在の加速度データに基づき前記異常の有無を判定する診断処理部と、を備える。
(1) One aspect of the present invention is an anomaly diagnosis system that sequentially learns an operating state of an object to be diagnosed and diagnoses an anomaly,
a correction processing unit that performs fluctuation correction based on a moving average of an acceleration amplitude obtained from the acceleration data of the diagnosis object;
a sequential learning processing unit that sequentially learns the result of the fluctuation correction and generates a learning result model;
and a diagnosis processing unit that uses the learning result model to determine the presence or absence of the abnormality based on current acceleration data of the diagnosis target.
(2)本発明の他の態様は、診断対象の稼働状態を逐次学習して異常を診断する際、診断対象の加速度データの加速度振幅の揺らぎを補正する装置であって、
前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を行ことで前記揺らぎを補正する移動平均部を備える。
(2) Another aspect of the present invention is a device for correcting fluctuations in acceleration amplitude of acceleration data of a diagnosis target when diagnosing an abnormality by sequentially learning an operating state of the diagnosis target, comprising:
The acceleration data is frequency-converted to obtain frequency data, and a moving average unit is provided for performing a moving average process of the frequency data to correct the fluctuation.
(3)本発明のさらに他の態様は、コンピュータにより診断対象の稼働状態を逐次学習して異常を診断する異常診断方法であって、
前記診断対象の加速度データから得られた加速度振幅の移動平均に基づき揺らぎ補正を行う補正処理ステップと、
前記揺らぎ補正の結果を逐次学習して学習結果モデルを生成する逐次学習処理ステップと、
前記学習結果モデルを用いて診断対象の現在の加速度データに基づき前記異常の有無を判定する診断処理ステップと、を有する。
(3) Yet another aspect of the present invention is an anomaly diagnosis method for diagnosing an anomaly by sequentially learning an operating state of an object to be diagnosed by a computer, comprising:
a correction processing step of performing fluctuation correction based on a moving average of acceleration amplitude obtained from the acceleration data of the diagnosis object;
a sequential learning processing step of sequentially learning the result of the fluctuation correction to generate a learning result model;
and a diagnosis processing step of determining the presence or absence of an abnormality based on current acceleration data of an object to be diagnosed using the learning result model.
(4)本発明のさらに他の態様は、診断対象の稼働状態を逐次学習して異常を診断する際、コンピュータにより診断対象の加速度データの加速度振幅の揺らぎを補正する方法であって、
前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を行ことで前記揺らぎを補正する移動平均ステップを有する。
(4) Yet another aspect of the present invention is a method for correcting fluctuations in acceleration amplitude of acceleration data of a diagnosis target by a computer when diagnosing an abnormality by sequentially learning an operating state of the diagnosis target, comprising:
The method includes a moving average step of performing a moving average process of a frequency based on frequency data obtained by frequency-converting the acceleration data, thereby correcting the fluctuation.
(5)なお、本発明は、コンピュータを前記異常診断システム/前記周波数の揺らぎ補正処理装置として機能させるプログラムとして構成することもできる。 (5) The present invention can also be configured as a program that causes a computer to function as the abnormality diagnosis system/frequency fluctuation correction processing device.
本発明によれば、加速度データの振幅の揺らぎの影響を抑えて異常診断の精度向上を図ることができる。 The present invention makes it possible to reduce the effects of fluctuations in the amplitude of acceleration data and improve the accuracy of abnormality diagnosis.
以下、本発明の実施形態に係る異常診断システムを説明する。この異常診断システムは、診断対象の周波数解析時に移動平均処理を施すことにより、周波数の揺らぎ補正を行って異常診断の精度向上を図っている。ここでは回転機(電動機・発電機)を診断対象の機器の一例として説明する。 The following describes an abnormality diagnosis system according to an embodiment of the present invention. This abnormality diagnosis system performs moving average processing during frequency analysis of the object to be diagnosed, thereby correcting frequency fluctuations and improving the accuracy of abnormality diagnosis. Here, a rotating machine (electric motor/generator) will be described as an example of the device to be diagnosed.
図4中の1は、ポンプ・ターボブロワ・水車などの回転体設備4に接続された診断対象の回転機を示している。この回転機1の軸受などには加速度センサ2が設置され、また加速度センサ2の出力信号を加速度データとして計測するデータ計測器5が設置されている。
In Figure 4, 1 indicates the rotating machine to be diagnosed, which is connected to rotating
このような設置状態の下、前記異常診断システム6は、データ計測器5の計測した加速度データを解析し、回転機1の稼働状態を診断する。このとき加速度センサ2・データ計測器5・前記異常診断システム6は、有線/無線によりそれぞれ接続され、データ送受信が可能となっている。以下、前記異常診断システム6を実施例1~3に基づき具体的に説明する。
In this installation state, the
図5~図7に基づき前記異常診断システム6の実施例1を説明する。ここでは前記異常診断システム6には、計測器5により計測された回転機1の加速度データが入力される。
A first embodiment of the
(1)前記異常診断システムの構成例
前記異常診断システム6は、コンピュータにより構成され、通常のコンピュータのハードウェアリソース(例えばCPU,RAM・ROMなどの一次憶装置,HDD,SSDなどの二次記憶装置)を備えている。
(1) Example of the Configuration of the Abnormality Diagnosis System The
このハードウェアリソースとソフトウェアリソース(OS,アプリケーションなど)との協働の結果、前記異常診断システム6は、図5に示すように、記憶部7,周波数の揺らぎ補正処理部8,過渡期除去部9,学習期間選定部10,逐次学習処理部11,診断処理部12,診断結果伝送部13とを実装する。
As a result of the cooperation between these hardware resources and software resources (OS, applications, etc.), the
記憶部7は、コンピュータの一次記憶装置/二次記憶装置に構築され、前記異常診断システム6に入力された加速度データが記憶され、また各部8~12の出力データが記憶される。
The
前記補正処理部8は、記憶部7の加速度データを入力とし、図6に示すように、周波数変換部15,移動平均部16を実装する。この周波数変換部15は、入力された加速度データを周波数データに変換し、変換された周波数データを移動平均部16に出力する。
The
移動平均部16には、周波数変換部15の出力した周波数データが入力され、入力された周波数データについて周波数強度の移動平均処理を行うことで周波数の揺らぎ補正を実行する。
The frequency data output by the
このとき本実施例では、事前に加速度データの計測時間に応じて定めた周波数の幅を移動平均の幅として、周波数強度の移動平均処理を行うものとする。なお、前記揺らぎ補正結果、即ち前記揺らぎ補正後の周波数データは記憶部7に出力される。
In this embodiment, the frequency width determined in advance according to the measurement time of the acceleration data is used as the width of the moving average, and moving average processing of the frequency intensity is performed. The fluctuation correction result, i.e., the frequency data after the fluctuation correction, is output to the
過渡期除去部9は、図5に示すように、記憶部7に記憶された前記揺らぎ補正結果(前記揺らぎ補正後の周波数データ)を入力とし、回転機1の稼働状態が過渡期であるか否かを判定する。このとき過渡期と判定すれば、該周波数データを逐次学習・診断の対象データ(学習データ・診断データ)から除外する。この過渡期除去結果は記憶部7に出力される。
As shown in FIG. 5, the transient
学習期間選定部10は、過渡期ではないと判定された学習データを入力として逐次学習の期間を自動的に選定する。ここで選定された学習期間データは、記憶部7に出力される。
The learning
逐次学習処理部11は、記憶部7に記憶された前記揺らぎ補正結果および学習期間データを入力とする。ここでは学習期間に応じた前記揺らぎ補正後の周波数データに基づき回転機1の正常な稼働状態を逐次的に学習し、学習結果モデルを生成する。この逐次学習の結果、生成された学習結果モデルは、記憶部7に出力される。
The sequential
診断処理部12は、記憶部7に記憶された現在の回転機1の加速度データおよび学習結果モデルを入力とする。ここでは学習結果モデルを用いて現在の加速度データから回転機1の状態を診断し、診断結果データを記憶部7に出力する。
The
なお、診断結果伝送部13は、記憶部7に記憶された診断結果データを入力とし、入力された診断結果データを監視制御システムやクラウドシステムに伝送する通信制御部として機能する。
The diagnostic
(2)前記異常診断システム6の処理内容
図7に基づき前記異常診断システム6の処理内容(S01~S09)を説明する。ここでは回転機1の正常状態を逐次学習する学習ステージ(S01~S07)と、現在の回転機1の状態が異常か否かを判定する診断ステージ(S08,S09)とが実行される。
(2) Processing Contents of the
S01:正常状態の回転機1の稼働により処理が開始され、加速度センサ2の出力信号がデータ計測器5により加速度データとして計測され、計測された加速度データが記憶部7に順次記憶される。
S01: Processing begins with the
S02:前記補正処理部8は、記憶部7に記憶された加速度について周波数の揺らぎ補正を実行する。ここでは周波数変換部15が、入力された加速度データを定Q変換(特許文献1,2参照)やフーリエ変換などの手法により周波数データに変換する。変換された周波数データに対して、移動平均部16が周波数の揺らぎ補正を実行する。
S02: The
この場合、図2(b)に示した不安定な土台3などによる振動の影響で周波数ずれが発生したとしても、「1Hz」以上のずれが発生することは殆ど無い。そのため、事前に「1Hz」の範囲の幅を移動平均の幅と定めて、周波数強度の移動平均処理を実行する。
In this case, even if a frequency shift occurs due to vibrations caused by an
このような移動平均処理による周波数の揺らぎ補正は、入力された加速度データに対する周波数解析として行われる。このとき加速度データの計測時間によって周波数の分解能が異なるので、計測時間を考慮して移動平均の幅を定める。 This type of frequency fluctuation correction using moving average processing is performed as a frequency analysis of the input acceleration data. Since the frequency resolution differs depending on the measurement time of the acceleration data, the width of the moving average is determined taking into account the measurement time.
例えば計測時間が「10秒」の場合には周波数の分解能が「0.1Hz」刻みとなる。そこで、時系列の周波数データに対して「1Hz」の範囲の移動平均を行うために移動平均の幅を「0.1Hz」として10サンプル分化し、各サンプルの平均値をグラフ化して移動平均の処理を実行する。 For example, if the measurement time is "10 seconds", the frequency resolution will be in increments of "0.1 Hz". Therefore, to perform a moving average of the time series frequency data in a range of "1 Hz", the moving average width is set to "0.1 Hz" and 10 samples are divided, and the average value of each sample is graphed to perform the moving average processing.
S03:過渡期除去部9が、前記補正後の周波数データに回転機1の稼働の過渡期のデータが含まれるか否かを確認する。例えば前記補正後の周波数データを分割し、分割された各区間の周波数ピーク値を求め、各周波数のピーク値間に閾値(1~2Hz)以上の差があれば過渡期と判定する。確認の結果、過渡期と判定されれば、学習データ・診断データから除外し、記憶部7に過渡期除去結果を出力して記憶部7の記憶データを更新する。
S03: The transient
S04:学習期間選定部10は、S03の処理後に記憶部7から学習データ(過渡期と判定されなかった学習データ)を取得し、加速度センサ2のセンサ出力が安定した期間を学習期間として選定する。
S04: The learning
例えばセンサ設置後や修繕後は実行値が高くなるため、加速度データの実効値の移動平均を計算し、該移動平均の傾きを求める。この傾きの標準偏差を計算し、該傾きが一週間「±δ」の範囲に収まった時点を学習期間の開始時期として選定することができる。 For example, since the effective value will be higher after a sensor is installed or repaired, the moving average of the effective value of the acceleration data is calculated and the slope of this moving average is obtained. The standard deviation of this slope is calculated, and the point at which the slope falls within the range of "±δ" for one week can be selected as the start of the learning period.
S05:S04で学習期間が選定された場合、逐次学習処理部11が1年間の逐次学習を行って学習結果モデルを生成し、S06に進む。この学習手法については、例えば特許文献1,2の手法を用いることができる。
S05: If a learning period is selected in S04, the sequential
このとき前記異常診断システムを構成するコンピュータのCPU性能によっては、例えば毎日一つの学習データを学習してもよく、あるいは1週間に一つの学習データを学習してよりものとする。なお、S04で学習期間が選定されなかった場合はS01の処理に戻る。 At this time, depending on the CPU performance of the computer constituting the abnormality diagnosis system, for example, one piece of learning data may be learned every day, or one piece of learning data may be learned per week. If a learning period is not selected in S04, the process returns to S01.
S06,S07:図1に示したような、年変動を抑制するためには最低1年間は逐次学習を継続する必要がある。そこで、逐次学習処理部11は、S05で1年分学習が終了しているか否かを確認する。確認の結果、1年分の学習が終了していなければS07に進んで逐次学習を再開し、学習結果モデルを生成する。一方、終了していればS08に進む。
S06, S07: As shown in FIG. 1, in order to suppress yearly fluctuations, it is necessary to continue sequential learning for at least one year. Therefore, in S05, the sequential
なお、1年分の学習終了後は、年変動の影響が無くなるので逐次学習を行う必要がないが、引き続き逐次学習を継続してもよいものとする。 After one year of study is completed, the effects of annual variation will disappear and so sequential study is no longer necessary, but sequential study may be continued.
S08,S09:診断処理部12は、学習結果モデルを用いて現在の加速度データから回転機1の状態を診断し、診断結果を記憶部7に記憶する(S08)。ただし、現在の加速度データに周波数の揺らぎ補正(S02)・過渡期除去(S03)と同様の処理を施し、かかる処理後のデータを診断データとして診断処理部12により前記診断が行われるものとする。
S08, S09: The
このとき基本的に回転機1などの設備の機器に異常が発生した場合、異常状態が継続することが多い。そこで、診断精度を向上させる観点から1度の診断結果だけではなく、連続して異常判定された場合に回転機1に異常が発生したもの診断する設定が好ましい。なお、診断結果伝送部13は、記憶部7に記憶されたS08の診断結果データを監視制御システムやクラウドなどに有線/無線により伝送する(S09)。
このように本実施例の前記異常診断システム6によれば、回転機1の加速度データを周波数データに変換し、変換された周波数データの周波数の範囲毎に周波数強度の移動平均処理を施すことで周波数の揺らぎ補正が実行される(S02)。
In this case, when an abnormality occurs in equipment of the facility such as the rotating
As described above, according to the
これにより加速度データが回転機1の振動以外の要因の影響を受ける場合に前記要因の影響を抑制することが可能となる。例えば回転機1が図2(b)の不安定な土台3上に設置されている場合であっても、該土台3の振動の影響が抑制され、正確な学習データ・診断データを得ることができる。
This makes it possible to suppress the influence of factors other than the vibration of the
その結果、S03の過渡期除去処理・S04の学習期間選定処理を経て良好な学習結果モデルを生成でき、S08の診断精度が向上する。また、正確な診断データが得られるため、回転機1の現在の稼働状況を正確に把握でき、この点でも診断精度を向上させることができる。
As a result, a good learning result model can be generated through the transient period removal process in S03 and the learning period selection process in S04, improving the diagnostic accuracy in S08. In addition, accurate diagnostic data can be obtained, so the current operating status of the
図8および図9に基づき前記異常診断システム6の実施例2を説明する。ここでは実施例1の「周波数の揺らぎ補正処理(S02)」において、周波数帯毎に移動平均の幅を変えることで実施例1よりも揺らぎの影響を抑制する。
A second embodiment of the
すなわち、図8に示すように、周波数は倍々の位置にN次の周波数が現れるため、高周波になればなるほど周波数のずれが大きくなる。そこで、本実施例では、周波数帯に応じて移動平均の幅を変えて、周波数強度の移動平均処理を実行する。 That is, as shown in FIG. 8, Nth order frequencies appear at doubled frequencies, so the higher the frequency, the greater the frequency deviation. Therefore, in this embodiment, the width of the moving average is changed according to the frequency band, and moving average processing of frequency intensity is performed.
一般的に加速度センサ2の周波数範囲は「5kHz~10kHz」なので、その周波数範囲を想定する。このとき移動平均の幅を変える周波数の区切りは任意でよいものとする。
Generally, the frequency range of the
ただし、通常、(A)「0~1kHz」は変位や速度などの低周波領域、(B)「1kHz」以上は加速度などの高周波領域とみなされるので、(A)(B)の周波数帯に分けることが好ましい。なお、さらに細分化して(a)1~3kHz、(b)3~5kHz、(c)5kH以上に分けてもよい。 However, since (A) "0 to 1 kHz" is usually considered to be a low-frequency region such as displacement and velocity, and (B) "1 kHz" or higher is considered to be a high-frequency region such as acceleration, it is preferable to divide the frequency bands into (A) and (B). It may also be further subdivided into (a) 1 to 3 kHz, (b) 3 to 5 kHz, and (c) 5 kHz or higher.
前述の(A)(B)の周波数帯に分けた場合に移動平均部16は、図9に示すように、(A)の周波数帯について「1Hz」の範囲で低周波の移動平均を行い、(B)の周波数帯について「10Hz」の範囲で高周波の移動平均を行う。この移動平均の幅「1Hz,10Hz」は一例であって任意に設定してよいものとする。
When the frequency bands are divided into (A) and (B) as described above, the moving
このように本実施例の前記異常診断システム6によれば、周波数帯に応じて移動平均の幅を変えて周波数強度の移動平均処理が実行されるため、特に回転周波数の高調波のずれの影響が抑制され、この点で診断精度をさらに向上させることができる。
In this way, according to the
図10に基づき前記異常診断システム6の実施例3を説明する。ここでは移動平均部16における移動平均の幅を周波数に応じて変化(変動)させることで周波数の揺らぎ補正が実行される。
A third embodiment of the
すなわち、回転機1の異常診断においては回転周波数が重要となる。回転機1は、回転周波数のN次毎に回転周波数が現れるため、実施例2のように移動平均の幅を固定した場合、適切な前記揺らぎ補正が困難なおそれがある。
In other words, the rotational frequency is important in diagnosing abnormalities in the
例えば回転周波数が「10Hz」の場合、「20Hz」,「30Hz」などのN次に回転周波数成分が現れ、また回転周波数が「50Hz」の場合には、「100Hz」,「150Hz」などのN次に回転周波数成分が現れる。 For example, if the rotation frequency is "10 Hz", N-th order rotation frequency components such as "20 Hz" and "30 Hz" will appear, and if the rotation frequency is "50 Hz", N-th order rotation frequency components such as "100 Hz" and "150 Hz" will appear.
そうすると図2(b)の不安定な土台3などの影響により回転周波数に「0.1Hz」のずれが発生した場合、N次毎に周波数のずれが大きくなってしまう。
As a result, if a deviation of 0.1 Hz occurs in the rotation frequency due to the influence of an
ここで「1kH」までの周波数のずれを想定すると、回転周波数が「10Hz」の場合には「1kHz」までで100次の回転周波数の高調波が発生することとなる。そのため、「1kH」までに「0.1(Hz)×100(次)=10(Hz)」の周波数ずれが発生する。 If we assume a frequency shift of up to 1 kHz, then if the rotation frequency is 10 Hz, then up to 1 kHz, 100th order rotation frequency harmonics will be generated. Therefore, up to 1 kHz, a frequency shift of 0.1 (Hz) x 100 (order) = 10 (Hz) will occur.
また、回転周波数が「50Hz」の場合には、「1kHz」までに20次の回転周波数の高調波が発生するため、「1kHz」までに「0.1(Hz)×20(次)=2(Hz)」の周波数ずれが発生する。このような事情からすれば回転周波数に応じて移動平均の幅を変更することで前記揺らぎ補正をより効果的に行える。 In addition, when the rotation frequency is "50 Hz", the 20th rotation frequency harmonic occurs by "1 kHz", so a frequency shift of "0.1 (Hz) x 20 (order) = 2 (Hz)" occurs by "1 kHz". Given this, the fluctuation correction can be performed more effectively by changing the width of the moving average according to the rotation frequency.
そこで、回転機1の回転周波数を基準にして、図10に示すように、移動平均の幅を周波数変換部15で変換された周波数データの周波数帯「1~n」に応じて自動的に変えて周波数強度の移動平均処理を行う。具体的には以下のステップ(S11~S15:図示省略)により移動平均処理を行う。
Therefore, using the rotational frequency of the
S11:回転機1の回転周波数を設定する。これは回転機1の仕様から求めて定めることができる。例えば「10Hz」などと回転周波数を設定しておくものとする。
S11: Set the rotation frequency of the
S12:回転機1の加速度データの計測時間を設定する。例えば「5秒」などと計測時間を設定しておくものとする。
S12: Set the measurement time for the acceleration data of the
S13:回転機1の周波数のずれを設定する。例えば「0.1Hz」などと周波数のずれを設定しておくものとする。
S13: Set the frequency deviation of the
S14:現在の周波数位置、即ち周波数データ中の現在位置における移動平均の幅を計算する。ここで
・現在の周波数位置:X[Hz]
・回転周波数:N[Hz]
・計測時間:T[秒]
・周波数のずれ:H[Hz]
とすれば、移動平均の幅[W]は式(1)により計算することができる。
式(1):W=X[Hz]÷N[Hz]×H[Hz]×T[秒]
S15:S15の計算により移動平均の幅を計算し、移動平均の幅を変化させながら周波数強度の移動平均処理を行う。
S14: Calculate the width of the moving average at the current frequency position, i.e., the current position in the frequency data. Current frequency position: X [Hz]
・Rotation frequency: N [Hz]
・Measurement time: T [seconds]
- Frequency deviation: H [Hz]
Then, the width of the moving average [W] can be calculated using equation (1).
Formula (1): W=X[Hz]÷N[Hz]×H[Hz]×T[seconds]
S15: The width of the moving average is calculated by the calculation in S15, and the moving average process of the frequency intensity is performed while changing the width of the moving average.
このように本実施例の前記異常診断システム6によれば、現在の周波数位置に応じた移動平均の幅が回転機1の回転周波数を基に自動的に計算されるため、低周波や高周波,高調波などを自動で考慮した移動平均処理を行うことができ、簡易に診断精度を向上させることができる。
In this way, according to the
なお、本発明は、上記実施形態に限定されるものではなく、各請求項に記載された範囲内で変形して実施することができる。例えば前記補正処理部8を前記異常診断システム6とは別のコンピュータに周波数の揺らぎ補正処理装置(周波数変換部15,移動平均部16)として構成してもよい。
The present invention is not limited to the above embodiment, and can be modified and implemented within the scope of the claims. For example, the
また、本発明は、コンピュータを前記異常判定システム6/前記周波数の揺らぎ補正装置として機能させる異常診断プログラム/周波数の揺らぎ補正処理プログラムとして構成することもできる。
The present invention can also be configured as an abnormality diagnosis program/frequency fluctuation correction processing program that causes a computer to function as the
この異常診断プログラムによれば、コンピュータによりS01~S09,S11~S15などの処理ステップ(異常診断方法)が実行される。一方、周波数の揺らぎ補正処理プログラムによれば、コンピュータによりS02,S11~S15の処理ステップ(周波数の揺らぎ補正方法)が実行される。 According to this abnormality diagnosis program, the computer executes processing steps S01 to S09, S11 to S15, etc. (an abnormality diagnosis method). On the other hand, according to the frequency fluctuation correction processing program, the computer executes processing steps S02, S11 to S15 (a frequency fluctuation correction method).
1…回転機(診断対象)
2…加速度センサ
3…土台
4…回転体設備
5…データ計測器
6…異常診断システム
7…記憶部
8…周波数の揺らぎ補正処理部(周波数の揺らぎ補正装置)
9…過渡期除去部
10…学習期間選定部
11…逐次学習処理部
12…診断処理部
13…診断結果伝送部
15…周波数変換部
16…移動平均部
1...Rotating machine (diagnosis target)
2: Acceleration sensor 3: Base 4: Rotating equipment 5: Data measuring instrument 6: Abnormality diagnosis system 7: Memory unit 8: Frequency fluctuation correction processing unit (frequency fluctuation correction device)
9: Transient period removal section 10: Learning period selection section 11: Sequential learning processing section 12: Diagnosis processing section 13: Diagnosis result transmission section 15: Frequency conversion section 16: Moving average section
Claims (10)
前記診断対象の加速度データについて周波数の揺らぎ補正を行う補正処理部と、
前記揺らぎ補正の結果を逐次学習して学習結果モデルを生成する逐次学習処理部と、
前記学習結果モデルを用いて診断対象の現在の加速度データに基づき前記異常の有無を判定する診断処理部と、
を備え、
前記補正処理部は、前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を実行する移動平均部を備え、
前記移動平均部は、周波数の幅を移動平均の幅とし、前記周波数データの周波数帯域毎に移動平均の幅を定める
ことを特徴とする異常診断システム。 An abnormality diagnosis system that sequentially learns an operating state of a diagnosis target and diagnoses an abnormality,
a correction processing unit that performs frequency fluctuation correction on the acceleration data of the diagnosis object;
a sequential learning processing unit that sequentially learns the result of the fluctuation correction and generates a learning result model;
a diagnosis processing unit that uses the learning result model to determine the presence or absence of the abnormality based on current acceleration data of the diagnosis target;
Equipped with
the correction processing unit includes a moving average unit that performs a moving average process of a frequency based on frequency data obtained by frequency-converting the acceleration data,
The moving average unit sets a frequency width as a moving average width and determines the moving average width for each frequency band of the frequency data.
An abnormality diagnosis system comprising:
前記診断対象の加速度データについて周波数の揺らぎ補正を行う補正処理部と、
前記揺らぎ補正の結果を逐次学習して学習結果モデルを生成する逐次学習処理部と、
前記学習結果モデルを用いて診断対象の現在の加速度データに基づき前記異常の有無を判定する診断処理部と、
を備え、
前記補正処理部は、前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を実行する移動平均部を備え、
前記移動平均部は、周波数の幅を移動平均の幅とし、前記加速度データを周波数変換した周波数データの周波数帯に応じて移動平均の幅を変動させる
ことを特徴とする異常診断システム。 An abnormality diagnosis system that sequentially learns an operating state of a diagnosis target and diagnoses an abnormality,
a correction processing unit that performs frequency fluctuation correction on the acceleration data of the diagnosis object;
a sequential learning processing unit that sequentially learns the result of the fluctuation correction and generates a learning result model;
a diagnosis processing unit that uses the learning result model to determine the presence or absence of the abnormality based on current acceleration data of the diagnosis target;
Equipped with
the correction processing unit includes a moving average unit that performs a moving average process of a frequency based on frequency data obtained by frequency-converting the acceleration data,
The moving average unit sets a frequency width as a moving average width, and varies the moving average width according to a frequency band of frequency data obtained by frequency-converting the acceleration data.
An abnormality diagnosis system comprising:
前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を行ことで前記揺らぎを補正する移動平均部を備え、
前記移動平均部は、周波数の幅を移動平均の幅とし、前記周波数データの周波数帯域毎に移動平均の幅を定める
ことを特徴とする補正処理装置。 A device for correcting frequency fluctuations in acceleration data of a diagnosis target when diagnosing an abnormality by sequentially learning an operating state of the diagnosis target, comprising: a moving average unit that corrects the fluctuations by performing a moving average process of the frequency based on frequency data obtained by frequency-converting the acceleration data;
The moving average unit sets a frequency width as a moving average width and determines the moving average width for each frequency band of the frequency data.
A correction processing device comprising:
前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を行ことで前記揺らぎを補正する移動平均部を備え、
前記移動平均部は、周波数の幅を移動平均の幅とし、前記加速度データを周波数変換した周波数データの周波数帯に応じて移動平均の幅を変動させる
ことを特徴とする補正処理装置。 A device for correcting frequency fluctuations in acceleration data of a diagnosis target when diagnosing an abnormality by sequentially learning an operating state of the diagnosis target, comprising: a moving average unit that corrects the fluctuations by performing a moving average process of the frequency based on frequency data obtained by frequency-converting the acceleration data;
The moving average unit sets a frequency width as a moving average width, and varies the moving average width according to a frequency band of frequency data obtained by frequency-converting the acceleration data.
A correction processing device comprising:
前記診断対象の加速度データについて周波数の揺らぎ補正を行う補正処理ステップと、
前記揺らぎ補正の結果を逐次学習して学習結果モデルを生成する逐次学習処理ステップと、
前記学習結果モデルを用いて診断対象の現在の加速度データに基づき前記異常の有無を判定する診断処理ステップと、
を有し、
前記補正処理ステップは、前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を実行する移動平均ステップを有し、
前記移動平均ステップは、周波数の幅を移動平均の幅とし、前記周波数データの周波数帯域毎に移動平均の幅を定める
ことを特徴とする異常診断方法。 A method for diagnosing an abnormality by sequentially learning an operating state of an object to be diagnosed by a computer, comprising the steps of:
a correction processing step of correcting frequency fluctuations in the acceleration data of the diagnosis object;
a sequential learning processing step of sequentially learning the result of the fluctuation correction to generate a learning result model;
a diagnosis processing step of determining the presence or absence of the abnormality based on current acceleration data of a diagnosis target using the learning result model;
having
The correction processing step includes a moving average step of performing a moving average process of a frequency based on frequency data obtained by frequency-converting the acceleration data,
The abnormality diagnosis method according to claim 1, wherein the moving average step defines a frequency width as a moving average width for each frequency band of the frequency data.
前記診断対象の加速度データについて周波数の揺らぎ補正を行う補正処理ステップと、
前記揺らぎ補正の結果を逐次学習して学習結果モデルを生成する逐次学習処理ステップと、
前記学習結果モデルを用いて診断対象の現在の加速度データに基づき前記異常の有無を判定する診断処理ステップと、
を有し、
前記補正処理ステップは、前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を実行する移動平均ステップを有し、
前記移動平均ステップは、周波数の幅を移動平均の幅とし、前記加速度データを周波数変換した周波数データの周波数帯に応じて移動平均の幅を変動させる
ことを特徴とする異常診断方法。
A method for diagnosing an abnormality by sequentially learning an operating state of an object to be diagnosed by a computer, comprising the steps of:
a correction processing step of correcting frequency fluctuations in the acceleration data of the diagnosis object;
a sequential learning processing step of sequentially learning the result of the fluctuation correction to generate a learning result model;
a diagnosis processing step of determining the presence or absence of the abnormality based on current acceleration data of a diagnosis target using the learning result model;
having
The correction processing step includes a moving average step of performing a moving average process of a frequency based on frequency data obtained by frequency-converting the acceleration data,
The abnormality diagnosis method according to claim 1, wherein the moving average step uses a frequency width as a moving average width, and varies the moving average width according to a frequency band of frequency data obtained by frequency-converting the acceleration data.
前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を行ことで前記揺らぎを補正する移動平均ステップを有し、
前記移動平均ステップは、周波数の幅を移動平均の幅とし、前記周波数データの周波数帯域毎に移動平均の幅を定める
ことを特徴とする補正処理方法 A method for correcting frequency fluctuations in acceleration data of a diagnosis target by a computer when diagnosing an abnormality by sequentially learning an operating state of the diagnosis target, comprising a moving average step of correcting the fluctuations by performing a moving average process of the frequency based on frequency data obtained by frequency-converting the acceleration data,
The moving average step determines a frequency width as a moving average width for each frequency band of the frequency data.
A correction processing method comprising:
前記加速度データを周波数変換した周波数データに基づき周波数の移動平均処理を行ことで前記揺らぎを補正する移動平均ステップを有し、
前記移動平均ステップは、周波数の幅を移動平均の幅とし、前記加速度データを周波数変換した周波数データの周波数帯に応じて移動平均の幅を変動させる
ことを特徴とする補正処理方法。 A method for correcting frequency fluctuations in acceleration data of a diagnosis target by a computer when diagnosing an abnormality by sequentially learning an operating state of the diagnosis target, comprising a moving average step of correcting the fluctuations by performing a moving average process of the frequency based on frequency data obtained by frequency-converting the acceleration data,
The moving average step uses a frequency width as a moving average width, and varies the moving average width according to a frequency band of frequency data obtained by frequency-converting the acceleration data.
A correction processing method comprising:
コンピュータを機能させることを特徴とする異常診断プログラム。 The abnormality diagnosis system according to claim 1 or 2 comprises:
An abnormality diagnosis program that causes a computer to function.
コンピュータを機能させることを特徴とする補正処理プログラム。 As the correction processing device according to claim 3 or 4 ,
A correction processing program that causes a computer to function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021018234A JP7647138B2 (en) | 2021-02-08 | 2021-02-08 | Abnormality diagnosis system and abnormality diagnosis method, frequency fluctuation correction processing device and correction processing method, abnormality diagnosis program, and frequency fluctuation correction processing program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021018234A JP7647138B2 (en) | 2021-02-08 | 2021-02-08 | Abnormality diagnosis system and abnormality diagnosis method, frequency fluctuation correction processing device and correction processing method, abnormality diagnosis program, and frequency fluctuation correction processing program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022121075A JP2022121075A (en) | 2022-08-19 |
| JP7647138B2 true JP7647138B2 (en) | 2025-03-18 |
Family
ID=82849552
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021018234A Active JP7647138B2 (en) | 2021-02-08 | 2021-02-08 | Abnormality diagnosis system and abnormality diagnosis method, frequency fluctuation correction processing device and correction processing method, abnormality diagnosis program, and frequency fluctuation correction processing program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7647138B2 (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006113002A (en) | 2004-10-18 | 2006-04-27 | Nsk Ltd | Abnormality diagnosis system for mechanical equipment |
| JP2006125976A (en) | 2004-10-28 | 2006-05-18 | Nsk Ltd | Abnormality diagnosis system for mechanical equipment |
| JP2013200144A (en) | 2012-03-23 | 2013-10-03 | Mitsubishi Electric Corp | Abnormal sound diagnosis device |
| WO2015068446A1 (en) | 2013-11-08 | 2015-05-14 | 三菱電機株式会社 | Abnormal sound diagnosis device |
| WO2017159784A1 (en) | 2016-03-17 | 2017-09-21 | Ntn株式会社 | Condition monitoring system and wind power generation device |
| JP2020030111A (en) | 2018-08-23 | 2020-02-27 | 株式会社明電舎 | Abnormality sign detection system, and abnormality sign detection method |
| JP2020046211A (en) | 2018-09-14 | 2020-03-26 | 株式会社椿本チエイン | Diagnostic device and method for diagnosis |
| US20200104200A1 (en) | 2018-09-27 | 2020-04-02 | Oracle International Corporation | Disk drive failure prediction with neural networks |
-
2021
- 2021-02-08 JP JP2021018234A patent/JP7647138B2/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006113002A (en) | 2004-10-18 | 2006-04-27 | Nsk Ltd | Abnormality diagnosis system for mechanical equipment |
| JP2006125976A (en) | 2004-10-28 | 2006-05-18 | Nsk Ltd | Abnormality diagnosis system for mechanical equipment |
| JP2013200144A (en) | 2012-03-23 | 2013-10-03 | Mitsubishi Electric Corp | Abnormal sound diagnosis device |
| WO2015068446A1 (en) | 2013-11-08 | 2015-05-14 | 三菱電機株式会社 | Abnormal sound diagnosis device |
| WO2017159784A1 (en) | 2016-03-17 | 2017-09-21 | Ntn株式会社 | Condition monitoring system and wind power generation device |
| JP2020030111A (en) | 2018-08-23 | 2020-02-27 | 株式会社明電舎 | Abnormality sign detection system, and abnormality sign detection method |
| JP2020046211A (en) | 2018-09-14 | 2020-03-26 | 株式会社椿本チエイン | Diagnostic device and method for diagnosis |
| US20200104200A1 (en) | 2018-09-27 | 2020-04-02 | Oracle International Corporation | Disk drive failure prediction with neural networks |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022121075A (en) | 2022-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3431952B1 (en) | Condition monitoring system and wind turbine generation apparatus | |
| Zhan et al. | Adaptive state detection of gearboxes under varying load conditions based on parametric modelling | |
| US11334414B2 (en) | Abnormality detecting apparatus, rotating machine, abnormality detection method, and non- transitory computer readable medium | |
| Picot et al. | Statistic-based spectral indicator for bearing fault detection in permanent-magnet synchronous machines using the stator current | |
| JP2931187B2 (en) | Pump deterioration diagnosis system | |
| CN115136487B (en) | Motor diagnostic device | |
| CN103134679B (en) | Bearing state monitors method and bearing state monitoring arrangement | |
| CN110858262B (en) | Abnormality detection device, abnormality detection method, and non-transitory computer-readable medium | |
| JP2001304954A (en) | Failure diagnosis method and device | |
| JP2001208655A (en) | Failure diagnosis method and device | |
| US12024180B2 (en) | Abnormality determination device | |
| JP3272144B2 (en) | Monitoring and diagnostic equipment for rotating equipment | |
| JP2000171291A (en) | Failure diagnosis method and failure diagnostic device | |
| JP2019045472A (en) | Status monitoring device and status monitoring method | |
| CN112149301A (en) | Model training method, fault prediction method and device and electronic equipment | |
| US20220334024A1 (en) | Vibration detection device, vibration detection method, and abnormality determination system | |
| WO2019171630A1 (en) | Diagnostic device, system, diagnostic method, and program | |
| JP7647138B2 (en) | Abnormality diagnosis system and abnormality diagnosis method, frequency fluctuation correction processing device and correction processing method, abnormality diagnosis program, and frequency fluctuation correction processing program | |
| JP7383367B1 (en) | Vibration data analysis method and analysis system for rotating equipment | |
| WO2018229836A1 (en) | Bearing diagnosis device | |
| JP7309079B2 (en) | Equipment condition monitoring device, abnormality determination system, and equipment condition monitoring method | |
| JP7552366B2 (en) | Abnormality diagnosis device, abnormality diagnosis system, and abnormality diagnosis method | |
| JPH01270623A (en) | Apparatus for diagnosing vibration of rotary machine | |
| JP7408519B2 (en) | Condition monitoring device and condition monitoring method | |
| JP2024087822A (en) | Anomaly detection device and anomaly detection method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230830 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240319 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240507 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240605 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20241001 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20241025 |
|
| 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: 20250204 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20250217 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7647138 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |