(57)【要約】【課題】ウェーブレット変換により歪を含む振動音響解析と機器異常部位推定とを精度良く且つ容易に行い、異常部位早期発見による二次損傷防止並びに余寿命予測による計画的保守を実現する。【解決手段】特徴量抽出処理工程で、運転中機器の振動音響の時系列データから特徴量を抽出し、機械学習診断工程で、機器正常作動状態での特徴量と機器の一般計測データを蓄積し、機械学習により構築した機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルに一般計測データと特徴量抽出処理工程で抽出された特徴量とを与え、確率分布モデルに基づいて異常度を算出し、異常部位推定工程で、予め算出された異常度の損傷時データセットと、機器の一般計測データから算出される異常度の実計測データセットを比較して機器の各部位毎の類似度を求め、類似度が高い順番に損傷可能性の高い部位として出力する。【選択図】図１
2019027814.tif 000002
【特許請求の範囲】
【請求項１】
運転中の機器において計測された振動音響の時系列データから特徴量を抽出する特徴量抽出処理工程と、
該特徴量抽出処理工程で抽出され且つ前記機器が正常に作動している状態での特徴量と前記機器の一般計測データとを蓄積し、機械学習により前記機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルを構築し、該確率分布モデルに前記機器の一般計測データと前記特徴量抽出処理工程で抽出された特徴量とを与え、前記確率分布モデルに基づいて異常度を算出する機械学習診断工程と
を有する振動音響解析方法。
【請求項２】
前記特徴量抽出処理工程は、
運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値、ピーク値、最大ピーク差、歪度、尖度、平均波高率、絶対平均振幅を抽出する基本特徴量抽出処理工程と、
前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するスペクトル解析特徴量抽出処理工程と、
前記時系列データをウェーブレット変換し複数の周波数帯に相当するスケールレベルに分解し、該スケールレベル毎のパワー時系列データを求める連続ウェーブレット変換処理工程と、
該連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値、最大値、最大波高率、絶対平均値を抽出するウェーブレット変換特徴量抽出処理工程と、
前記連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するウェーブレット変換スペクトル解析特徴量抽出処理工程と
を有する請求項１記載の振動音響解析方法。
【請求項３】
前記機器の一般計測データは、回転数、出力、仕事量、圧力、温度、速力の少なくとも一つである請求項１又は２記載の振動音響解析方法。
【請求項４】
運転中の機器において計測された振動音響の時系列データから特徴量を抽出する特徴量抽出処理部と、
該特徴量抽出処理部で抽出され且つ前記機器が正常に作動している状態での特徴量と前記機器の一般計測データとを蓄積し、機械学習により前記機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルを構築し、該確率分布モデルに前記機器の一般計測データと前記特徴量抽出処理部で抽出された特徴量とを与え、前記確率分布モデルに基づいて異常度を算出する機械学習診断部と
を備えた振動音響解析装置。
【請求項５】
前記特徴量抽出処理部は、
運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値、ピーク値、最大ピーク差、歪度、尖度、平均波高率、絶対平均振幅を抽出する基本特徴量抽出処理部と、
前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するスペクトル解析特徴量抽出処理部と、
前記時系列データをウェーブレット変換し複数の周波数帯に相当するスケールレベルに分解し、該スケールレベル毎のパワー時系列データを求める連続ウェーブレット変換処理部と、
該連続ウェーブレット変換処理部で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値、最大値、最大波高率、絶対平均値を抽出するウェーブレット変換特徴量抽出処理部と、
前記連続ウェーブレット変換処理部で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するウェーブレット変換スペクトル解析特徴量抽出処理部と
を備えた請求項４記載の振動音響解析装置。
【請求項６】
前記機器の一般計測データは、回転数、出力、仕事量、圧力、温度、速力の少なくとも一つである請求項４又は５記載の振動音響解析装置。
【請求項７】
請求項１〜３の何れか一項に記載の振動音響解析方法の前記機械学習診断工程において、機器の各部位毎に損傷時に発生する異常度の損傷時データセットを予め算出し、
前記機械学習診断工程で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとを比較して前記機器の各部位毎の類似度を求め、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力する異常部位推定工程を有する機器異常部位推定方法。
【請求項８】
請求項４〜６の何れか一項に記載の振動音響解析装置の前記機械学習診断部は、機器の各部位毎に損傷時に発生する異常度の損傷時データセットを予め算出するよう構成され、
前記機械学習診断部で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとを比較して前記機器の各部位毎の類似度を求め、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力する異常部位推定部を備えた機器異常部位推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、振動音響解析方法及び装置と機器異常部位推定方法及び装置に関するものである。
【背景技術】
【０００２】
一般に、エンジン等の回転体を含む機器において、異常を検知するために、歪を含む振動や音響を計測し時系列データを取ることが行われている。
【０００３】
前記時系列データの解析には、フーリエ変換等を用いたスペクトル解析が行われることが多い。単純なスペクトル解析では、時間領域の情報が欠損してしまうことがある。又、単純なスペクトル解析では、計測データに含まれる周波数分布が広い場合には、適切なスペクトル解析が得にくくなることがある。
【０００４】
これに対して、ウェーブレット変換を用いた手法が提案されている。該ウェーブレット変換では、前記フーリエ変換によって周波数特性を求める際に失われる時間領域の情報を残すことができ、該時間領域と周波数領域とに分解して解析が行われる。
【０００５】
尚、前記ウェーブレット変換を用いた手法と関連する一般的技術水準を示すものとしては、例えば、特許文献１がある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１０−４８６８４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、エンジン等の回転体を含む機器は、軸受部、ギア部、動力発生部といった多数の部品により構成されており、計測した振動音響データには各部からの情報が混在しているため、不適切な解析結果が示される場合があった。又、前記機器は、使用状態に応じて出力が増減され、振動音響の大きさが著しく変化するため、計測される振動音響データの大小のみで異常部位の推定を行うことは精度の面で適切な対応であるとは言えなかった。
【０００８】
因みに、回転速度が数百［ｍｉｎ−１］程度となる比較的低速の軸受等において、その振動音響データには、数［ｋＨｚ］〜数十［ｋＨｚ］の高周波帯に大きな成分が現れることが多い。これは、傷等による振動音響信号はパルス状の起振源により軸受部周囲がその固有振動数で振動し、主に高周波成分で構成されることが多いためである。このため、計測信号をスペクトル解析しても、回転速度付近の低周波領域のピーク周波数と振幅とを適切に抽出できないことがある。この対策としては、例えば、ローパスフィルタ等を用いてエンベロープ処理を行った後にスペクトル解析を行うことで低周期のピーク周波数と振幅とを抽出することが一般的である。この手法は、低周波領域の解析には適している反面、高周波領域の信号を解析することができない。又、ローパスフィルタの性能や特性により、ピーク周波数と振幅とが正しい値ではなくなり、評価に利用することが困難となる場合がある。
【０００９】
これに対し、本発明者等の研究により、以下のことが判明した。即ち、計測された時系列データをウェーブレット変換にて時間領域と周波数領域との複数のベクトルに分解し、続いて、各周波数帯の時系列ベクトルに対してスペクトル解析を行うことで、低周波領域から高周波領域までの広い範囲でピーク周波数と振幅とを適切に抽出することが可能となる。
【００１０】
但し、ウェーブレット変換による解析を行うことで一つの時系列データが多数の時系列データに分解されるため、この解析結果を人が判定することは多大な労力を必要とし、実用的であるとは言えなかった。
【００１１】
本発明は、上記従来の問題点に鑑みてなしたもので、ウェーブレット変換により歪を含む振動音響解析と機器の異常部位の推定とを精度良く且つ容易に行うことができ、異常部位の早期発見による二次損傷の防止、並びに余寿命予測による計画的保守を実現し得る振動音響解析方法及び装置と機器異常部位推定方法及び装置を提供しようとするものである。
【課題を解決するための手段】
【００１２】
上記目的を達成するために、本発明の振動音響解析方法は、運転中の機器において計測された振動音響の時系列データから特徴量を抽出する特徴量抽出処理工程と、
該特徴量抽出処理工程で抽出され且つ前記機器が正常に作動している状態での特徴量と前記機器の一般計測データとを蓄積し、機械学習により前記機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルを構築し、該確率分布モデルに前記機器の一般計測データと前記特徴量抽出処理工程で抽出された特徴量とを与え、前記確率分布モデルに基づいて異常度を算出する機械学習診断工程と
を有する。
【００１３】
前記振動音響解析方法において、前記特徴量抽出処理工程は、
運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値、ピーク値、最大ピーク差、歪度、尖度、平均波高率、絶対平均振幅を抽出する基本特徴量抽出処理工程と、
前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するスペクトル解析特徴量抽出処理工程と、
前記時系列データをウェーブレット変換し複数の周波数帯に相当するスケールレベルに分解し、該スケールレベル毎のパワー時系列データを求める連続ウェーブレット変換処理工程と、
該連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値、最大値、最大波高率、絶対平均値を抽出するウェーブレット変換特徴量抽出処理工程と、
前記連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するウェーブレット変換スペクトル解析特徴量抽出処理工程と
を有することができる。
【００１４】
前記振動音響解析方法において、前記機器の一般計測データは、回転数、出力、仕事量、圧力、温度、速力の少なくとも一つとすることができる。
【００１５】
一方、本発明の振動音響解析装置は、運転中の機器において計測された振動音響の時系列データから特徴量を抽出する特徴量抽出処理部と、
該特徴量抽出処理部で抽出され且つ前記機器が正常に作動している状態での特徴量と前記機器の一般計測データとを蓄積し、機械学習により前記機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルを構築し、該確率分布モデルに前記機器の一般計測データと前記特徴量抽出処理部で抽出された特徴量とを与え、前記確率分布モデルに基づいて異常度を算出する機械学習診断部と
を備えることができる。
【００１６】
前記振動音響解析装置において、前記特徴量抽出処理部は、
運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値、ピーク値、最大ピーク差、歪度、尖度、平均波高率、絶対平均振幅を抽出する基本特徴量抽出処理部と、
前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するスペクトル解析特徴量抽出処理部と、
前記時系列データをウェーブレット変換し複数の周波数帯に相当するスケールレベルに分解し、該スケールレベル毎のパワー時系列データを求める連続ウェーブレット変換処理部と、
該連続ウェーブレット変換処理部で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値、最大値、最大波高率、絶対平均値を抽出するウェーブレット変換特徴量抽出処理部と、
前記連続ウェーブレット変換処理部で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するウェーブレット変換スペクトル解析特徴量抽出処理部と
を備えることができる。
【００１７】
前記振動音響解析装置において、前記機器の一般計測データは、回転数、出力、仕事量、圧力、温度、速力の少なくとも一つとすることができる。
【００１８】
更に、本発明は、前記振動音響解析方法の前記機械学習診断工程において、機器の各部位毎に損傷時に発生する異常度の損傷時データセットを予め算出し、
前記機械学習診断工程で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとを比較して前記機器の各部位毎の類似度を求め、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力する異常部位推定工程を有する機器異常部位推定方法とすることができる。
【００１９】
又、本発明は、前記振動音響解析装置の前記機械学習診断部が、機器の各部位毎に損傷時に発生する異常度の損傷時データセットを予め算出するよう構成され、
前記機械学習診断部で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとを比較して前記機器の各部位毎の類似度を求め、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力する異常部位推定部を備える機器異常部位推定装置とすることができる。
【発明の効果】
【００２０】
本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置によれば、ウェーブレット変換により振動音響解析と機器の異常部位の推定とを精度良く且つ容易に行うことができ、異常部位の早期発見による二次損傷の防止、並びに余寿命予測による計画的保守を実現し得るという優れた効果を奏し得る。
【図面の簡単な説明】
【００２１】
【図１】本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置の実施例を示すフローチャートである。
【図２】本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置の実施例における特徴量抽出処理の詳細を示すフローチャートである。
【図３】本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置の実施例におけるブロック図である。
【図４】本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置の実施例における計測信号の波形例及びウェーブレット変換により解析した波形例を示す線図である。
【図５】本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置の実施例における部位の異常時のパターンを示す線図である。
【発明を実施するための形態】
【００２２】
以下、本発明の実施の形態を添付図面を参照して説明する。
【００２３】
図１〜図５は本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置の実施例である。
【００２４】
本実施例の振動音響解析方法の場合、図１に示す如く、特徴量抽出処理工程と、機械学習診断工程とを有している。
【００２５】
前記特徴量抽出処理工程は、運転中の機器において計測された振動音響の時系列データから特徴量を抽出するようになっている。
【００２６】
前記機械学習診断工程は、前記特徴量抽出処理工程で抽出され且つ前記機器が正常に作動している状態での特徴量と前記機器の一般計測データ（例えば、回転数、出力、仕事量、圧力、温度、速力）とを蓄積し、機械学習により前記機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルを構築し、該確率分布モデルに前記機器の一般計測データと前記特徴量抽出処理工程で抽出された特徴量とを与え、前記確率分布モデルに基づいて異常度を算出するようになっている。
【００２７】
図１の例では、ステップＳ１０として計測を開始し、多点計測の場合は同時に計測を行うようになっている。ステップＳ１１として一定時間の時系列データの入力を行い、ステップＳ１２として計測を停止するようになっている。
【００２８】
この後、ステップＳ１３，Ｓ１４として前記特徴量抽出処理工程と機械学習診断工程とを行うようになっている。
【００２９】
前記機械学習診断工程で算出された異常度をステップＳ１５として保存し蓄積するようになっている。
【００３０】
前記ステップＳ１５で異常度を保存し蓄積した後、ステップＳ１６として設定計測間隔時間が経過したか否かを判断し、設定計測間隔時間が経過していなければ、ステップＳ１７として一定時間待機してステップＳ１６に戻り、前記設定計測間隔時間が経過していれば、ステップＳ１０へ戻って計測を繰り返し行うようになっている。
【００３１】
前記ステップＳ１４の機械学習診断工程で算出された異常度を、ステップＳ１８として定義済みの正常範囲と比較し、前記異常度が正常範囲外である場合には、ステップＳ１９として警報を発令すると共に、ステップＳ２０として異常部位推定工程を行うようになっている。該異常部位推定工程は、機器異常部位推定方法において行われる工程となっている。
【００３２】
前記特徴量抽出処理工程は、詳述すると、前記ステップＳ１１において計測データとして入力された一定時間の時系列データを利用し、図２に示す如く、ステップＳ１００としての基本特徴量抽出処理工程と、ステップＳ２００としてのスペクトル解析特徴量抽出処理工程と、ステップＳ３００としての連続ウェーブレット変換処理工程と、ステップＳ４００としてのウェーブレット変換特徴量抽出処理工程と、ステップＳ５００としてのウェーブレット変換スペクトル解析特徴量抽出処理工程とを有している。
【００３３】
因みに、前記ステップＳ１１において計測データとして入力された一定時間の時系列データは、機器の計測部位の数がｎ個である場合、該機器の計測部位の数に相当するｎ個だけ存在する。
【００３４】
前記基本特徴量抽出処理工程は、運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値ＭＰ１、
【数１】
2019027814.tif 000003
ピーク平均値ＭＰ２、
【数２】
2019027814.tif 000004
（但し、ピークは、極大値、極小値の両方）
最大ピーク差ＭＰ３、
【数３】
2019027814.tif 000005
（但し、x+p：サンプル期間中の最大値、x-p：サンプル期間中の最小値）
歪度ＭＰ４、
【数４】
2019027814.tif 000006
但し、σ：標準偏差
【数５】
2019027814.tif 000007
尖度ＭＰ５、
【数６】
2019027814.tif 000008
平均波高率ＭＰ６、
【数７】
2019027814.tif 000009
（但し、［波高率(Crest factor)］＝［ピーク値］／［実効値］であるが、ピーク値はピーク平均値ＭＰ２を取っている。）
絶対平均振幅ＭＰ７
【数８】
2019027814.tif 000010
を抽出するようになっている（ステップＳ１１０参照）。尚、前記ステップＳ１００の基本特徴量抽出処理工程においては、ステップＳ１１０で抽出された特徴量を、ステップＳ１２０として定義済みの正常範囲と比較し、前記特徴量が正常範囲外である場合には、ステップＳ１３０として警報を発令するようになっている。
【００３５】
前記スペクトル解析特徴量抽出処理工程は、前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するようになっている（ステップＳ２１０参照）。尚、前記ステップＳ２００のスペクトル解析特徴量抽出処理工程においては、ステップＳ２１０で抽出された特徴量を、ステップＳ２２０として定義済みの正常範囲と比較し、前記特徴量が正常範囲外である場合には、ステップＳ２３０として警報を発令するようになっている。
【００３６】
前記連続ウェーブレット変換処理工程は、前記時系列データをウェーブレット変換し複数（例えば、ｍ個）の周波数帯に相当するスケールレベルに分解し、該スケールレベル毎のパワー時系列データを求めるようになっている（ステップＳ３１０参照）。
【００３７】
前記ウェーブレット変換特徴量抽出処理工程は、連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値ＷＰ１、
【数９】
2019027814.tif 000011
最大値ＷＰ２、
【数１０】
2019027814.tif 000012
最大波高率ＷＰ３、
【数１１】
2019027814.tif 000013
絶対平均値ＷＰ４
【数１２】
2019027814.tif 000014
を抽出するようになっている（ステップＳ４１０参照）。尚、前記ステップＳ４００のウェーブレット変換特徴量抽出処理工程においては、ステップＳ４１０で抽出された特徴量を、ステップＳ４２０として定義済みの正常範囲と比較し、前記特徴量が正常範囲外である場合には、ステップＳ４３０として警報を発令するようになっている。
【００３８】
前記ウェーブレット変換スペクトル解析特徴量抽出処理工程は、前記連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するようになっている（ステップＳ５１０参照）。尚、前記ステップＳ５００のウェーブレット変換スペクトル解析特徴量抽出処理工程においては、ステップＳ５１０で抽出された特徴量を、ステップＳ５２０として定義済みの正常範囲と比較し、前記特徴量が正常範囲外である場合には、ステップＳ５３０として警報を発令するようになっている。
【００３９】
前記基本特徴量抽出処理工程で抽出された特徴量（ステップＳ１１０参照）と、前記スペクトル解析特徴量抽出処理工程で抽出された特徴量（ステップＳ２１０参照）と、前記ウェーブレット変換特徴量抽出処理工程で抽出された特徴量（ステップＳ４１０参照）と、前記ウェーブレット変換スペクトル解析特徴量抽出処理工程で抽出された特徴量（ステップＳ５１０参照）と、ステップＳ６００として計測される機器の一般計測データ（例えば、回転数、出力、仕事量、圧力、温度、速力）とを、ステップＳ７００として保存し蓄積するようになっている。
【００４０】
前記異常部位推定工程は、前記機械学習診断工程で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとを比較して前記機器の各部位毎の類似度を求め、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力するようになっている。
【００４１】
一方、図３は本実施例におけるブロック図であり、振動音響解析装置として、特徴量抽出処理部１０と、機械学習診断部２０とを備え、機器異常部位推定装置として、前記振動音響解析装置の特徴量抽出処理部１０及び機械学習診断部２０と、異常部位推定部３０とを備えている。
【００４２】
前記特徴量抽出処理部１０は、運転中の機器において計測された振動音響の時系列データから特徴量を抽出するようになっている。
【００４３】
前記機械学習診断部２０は、前記特徴量抽出処理部１０で抽出され且つ前記機器が正常に作動している状態での特徴量と前記機器の一般計測データとを蓄積し、機械学習により前記機器正常作動状態での特徴量と一般計測データとの相関関係に関する確率分布モデルを構築し、該確率分布モデルに前記機器の一般計測データと前記特徴量抽出処理部１０で抽出された特徴量とを与え、前記確率分布モデルに基づいて異常度を算出するようになっている。
【００４４】
前記特徴量抽出処理部１０は、基本特徴量抽出処理部１１と、スペクトル解析特徴量抽出処理部１２と、連続ウェーブレット変換処理部１３と、ウェーブレット変換特徴量抽出処理部１４と、ウェーブレット変換スペクトル解析特徴量抽出処理部１５とを備えている。
【００４５】
前記基本特徴量抽出処理部１１は、運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値、ピーク値、最大ピーク差、歪度、尖度、平均波高率、絶対平均振幅を抽出するようになっている。
【００４６】
前記スペクトル解析特徴量抽出処理部１２は、前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するようになっている。
【００４７】
前記連続ウェーブレット変換処理部１３は、前記時系列データをウェーブレット変換し複数の周波数帯に相当するスケールレベルに分解し、該スケールレベル毎のパワー時系列データを求めるようになっている。
【００４８】
前記ウェーブレット変換特徴量抽出処理部１４は、前記連続ウェーブレット変換処理部１３で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値、最大値、最大波高率、絶対平均値を抽出するようになっている。
【００４９】
前記ウェーブレット変換スペクトル解析特徴量抽出処理部１５は、前記連続ウェーブレット変換処理部１３で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値を抽出するようになっている。
【００５０】
前記異常部位推定部３０は、前記機械学習診断部２０で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとを比較して前記機器の各部位毎の類似度を求め、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力するようになっている。
【００５１】
次に、上記実施例の作用を説明する。
【００５２】
図１において、先ず、機器の各部位の計測が開始され（ステップＳ１０参照）、一定時間の時系列データの入力が行われ（ステップＳ１１参照）、計測が停止された後（ステップＳ１２参照）、特徴量抽出処理工程が行われ（ステップＳ１３参照）、続いて、機械学習診断工程が行われる（ステップＳ１４参照）。
【００５３】
前記機械学習診断工程で算出された異常度は保存されて蓄積され（ステップＳ１５参照）、この後、設定計測間隔時間が経過したか否かが判断され（ステップＳ１６参照）、設定計測間隔時間が経過していなければ、一定時間待機（ステップＳ１７参照）してステップＳ１６に戻り、前記設定計測間隔時間が経過していれば、ステップＳ１０へ戻って計測が繰り返し行われる。
【００５４】
前記機械学習診断工程で算出された異常度は、定義済みの正常範囲と比較され（ステップＳ１８参照）、前記異常度が正常範囲外である場合には、警報が発令されると共に（ステップＳ１９参照）、異常部位推定工程が行われる（ステップＳ２０参照）。
【００５５】
前記特徴量抽出処理工程においては、計測データとして入力された一定時間の時系列データが利用され（ステップＳ１１参照）、図２に示す如く、基本特徴量抽出処理工程（ステップＳ１００参照）と、スペクトル解析特徴量抽出処理工程（ステップＳ２００参照）と、連続ウェーブレット変換処理工程（ステップＳ３００参照）と、ウェーブレット変換特徴量抽出処理工程（ステップＳ４００参照）と、ウェーブレット変換スペクトル解析特徴量抽出処理工程（ステップＳ５００参照）とが行われる。
【００５６】
前記基本特徴量抽出処理工程では、運転中の機器において計測された振動音響の時系列データの計測加速度から特徴量として実効値ＭＰ１（[数１］参照）、ピーク平均値ＭＰ２（[数２］参照）、最大ピーク差ＭＰ３（[数３］参照）、歪度ＭＰ４（[数４］参照）、尖度ＭＰ５（[数６］参照）、平均波高率ＭＰ６（[数７］参照）、絶対平均振幅ＭＰ７（[数８］参照）が抽出される（ステップＳ１１０参照）。尚、前記基本特徴量抽出処理工程において抽出された特徴量は、定義済みの正常範囲と比較され（ステップＳ１２０参照）、前記特徴量が正常範囲外である場合には、警報が発令される（ステップＳ１３０参照）。
【００５７】
前記スペクトル解析特徴量抽出処理工程では、前記時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値が抽出される（ステップＳ２１０参照）。尚、前記スペクトル解析特徴量抽出処理工程において抽出された特徴量は、定義済みの正常範囲と比較され（ステップＳ２２０参照）、前記特徴量が正常範囲外である場合には、警報が発令される（ステップＳ２３０参照）。
【００５８】
前記連続ウェーブレット変換処理工程では、前記時系列データがウェーブレット変換され複数（例えば、ｍ個）の周波数帯に相当するスケールレベルに分解され、該スケールレベル毎のパワー時系列データが求められる（ステップＳ３１０参照）。
【００５９】
ここで、本実施例における計測信号の波形及びウェーブレット変換により解析した波形は、例えば、図４に示すようになる。図４の例では、スケールパラメータを２ｓ（ｓ＝０〜７即ちｍ＝８）として１オクターブ毎に分解しており、１オクターブ毎に８スケールまで解析している。但し、各波形の縦軸は同一ではなく、プロット幅に合わせて適宜調整されている。尚、分解方法は自由であり、より高周波の成分を見たい場合には、１以下のスケールパラメータも適用できる。
【００６０】
前記ウェーブレット変換特徴量抽出処理工程では、連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データから各周波数帯の特徴量として実効値ＷＰ１（[数９］参照）、最大値ＷＰ２（[数１０］参照）、最大波高率ＷＰ３（[数１１］参照）、絶対平均値ＷＰ４（[数１２］参照）が抽出される（ステップＳ４１０参照）。尚、前記ウェーブレット変換特徴量抽出処理工程において抽出された特徴量は、定義済みの正常範囲と比較され（ステップＳ４２０参照）、前記特徴量が正常範囲外である場合には、警報が発令される（ステップＳ４３０参照）。
【００６１】
前記ウェーブレット変換スペクトル解析特徴量抽出処理工程では、前記連続ウェーブレット変換処理工程で求められたスケールレベル毎のパワー時系列データをスペクトル解析した後の特徴量として最大ピーク次数及びそのピーク値、第二ピーク次数及びそのピーク値、第三ピーク次数及びそのピーク値、第四ピーク次数及びそのピーク値、・・・第ｎピーク次数及びそのピーク値が抽出される（ステップＳ５１０参照）。尚、前記ステップＳ５００のウェーブレット変換スペクトル解析特徴量抽出処理工程において抽出された特徴量は、定義済みの正常範囲と比較され（ステップＳ５２０参照）、前記特徴量が正常範囲外である場合には、警報が発令される（ステップＳ５３０参照）。
【００６２】
前記基本特徴量抽出処理工程で抽出された特徴量（ステップＳ１１０参照）と、前記スペクトル解析特徴量抽出処理工程で抽出された特徴量（ステップＳ２１０参照）と、前記ウェーブレット変換特徴量抽出処理工程で抽出された特徴量（ステップＳ４１０参照）と、前記ウェーブレット変換スペクトル解析特徴量抽出処理工程で抽出された特徴量（ステップＳ５１０参照）と、例えば、回転数、出力、仕事量、圧力、温度、速力等の機器の一般計測データ（ステップＳ６００参照）は、保存され蓄積される（ステップＳ７００参照）。ここで、機器の計測部位の数はｎ個であるが、前記一般計測データは、診断対象となる機器（１ユニット）に対して１セットとなる。
【００６３】
前記異常部位推定工程においては、前記機械学習診断工程で予め算出された異常度の損傷時データセットと、前記機器の一般計測データから算出される異常度の実計測データセットとが比較されて前記機器の各部位毎の類似度が求められ、該機器の各部位毎の類似度が高い順番に損傷可能性の高い部位として出力される。
【００６４】
図５は本実施例における部位の異常時のパターンを示す線図である。図５に示す線図は、横軸を周波数レベル、縦軸を異常度とし、計測位置「イ」「ロ」・・・それぞれについて異常無の状態のデータに対し、例えば、部位Ａの異常時のパターンや部位Ｂの異常時のパターンが表示され、どの周波数レベルでの異常度が高くなっているかにより、損傷可能性の高い部位の特定が可能となる。
【００６５】
これにより、本実施例では、計測された時系列データをウェーブレット変換にて時間領域と周波数領域との複数のベクトルに分解し、続いて、各周波数帯の時系列ベクトルに対してスペクトル解析を行うことで、低周波領域から高周波領域までの広い範囲でピーク周波数と振幅とを適切に抽出することが可能となる。
【００６６】
しかも、一つの時系列データからウェーブレット変換により分解される多数の時系列データの解析結果を人が判定しなくて済み、多大な労力を必要とせず、実用的となる。
【００６７】
こうして、ウェーブレット変換により歪を含む振動音響解析と機器の異常部位の推定とを精度良く且つ容易に行うことができ、異常部位の早期発見による二次損傷の防止、並びに余寿命予測による計画的保守を実現し得る。
【００６８】
尚、本発明の振動音響解析方法及び装置と機器異常部位推定方法及び装置は、上述の実施例にのみ限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。
【符号の説明】
【００６９】
１０        特徴量抽出処理部
１１        基本特徴量抽出処理部
１２        スペクトル解析特徴量抽出処理部
１３        連続ウェーブレット変換処理部
１４        ウェーブレット変換特徴量抽出処理部
１５        ウェーブレット変換スペクトル解析特徴量抽出処理部
２０        機械学習診断部
３０        異常部位推定部
【図１】
2019027814.tif 000015
【図２】
2019027814.tif 000016
【図３】
2019027814.tif 000017
【図４】
2019027814.tif 000018
【図５】
2019027814.tif 000019