[go: up one dir, main page]

JP2018124117A - 状態監視システムおよび風力発電装置 - Google Patents

状態監視システムおよび風力発電装置 Download PDF

Info

Publication number
JP2018124117A
JP2018124117A JP2017015297A JP2017015297A JP2018124117A JP 2018124117 A JP2018124117 A JP 2018124117A JP 2017015297 A JP2017015297 A JP 2017015297A JP 2017015297 A JP2017015297 A JP 2017015297A JP 2018124117 A JP2018124117 A JP 2018124117A
Authority
JP
Japan
Prior art keywords
value
change
evaluation value
vibration waveform
data
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.)
Pending
Application number
JP2017015297A
Other languages
English (en)
Inventor
隆 長谷場
Takashi Haseba
隆 長谷場
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.)
NTN Corp
Original Assignee
NTN Corp
NTN Toyo Bearing Co Ltd
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 NTN Corp, NTN Toyo Bearing Co Ltd filed Critical NTN Corp
Priority to JP2017015297A priority Critical patent/JP2018124117A/ja
Priority to PCT/JP2018/001794 priority patent/WO2018142986A1/ja
Priority to CN201880009272.4A priority patent/CN110234972B/zh
Priority to US16/481,797 priority patent/US20200025648A1/en
Publication of JP2018124117A publication Critical patent/JP2018124117A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0224Process history based detection method, e.g. whereby history implies the availability of large amounts of data
    • G05B23/0227Qualitative history assessment, whereby the type of data acted upon, e.g. waveforms, images or patterns, is not relevant, e.g. rule based assessment; if-then decisions
    • G05B23/0232Qualitative history assessment, whereby the type of data acted upon, e.g. waveforms, images or patterns, is not relevant, e.g. rule based assessment; if-then decisions based on qualitative trend analysis, e.g. system evolution
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00Measuring characteristics of vibrations by using a detector in a fluid
    • G01H3/04Frequency
    • G01H3/08Analysing frequencies present in complex vibrations, e.g. comparing harmonics present
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • G01M13/045Acoustic or vibration analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05BINDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00Control
    • F05B2270/30Control parameters, e.g. input parameters
    • F05B2270/334Vibration measurements
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/20Pc systems
    • G05B2219/26Pc applications
    • G05B2219/2619Wind turbines

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Mechanical Engineering (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Sustainable Development (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Wind Motors (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

【課題】振動波形のトレンドデータの変化の検知感度の向上を簡便に実現する状態監視システムおよび風力発電装置を提供する。
【解決手段】
振動センサは、機械要素の振動波形を計測する。処理装置は、振動波形の変化を検知する。また、処理装置は、評価値演算部と検知部とを含む。評価値演算部は、一定時間内に前記振動センサから出力される振動波形データの実効値を特徴付ける評価値を、時間的に連続して演算する。検知部は、評価値に基づいて振動波形の変化を検知する。また、評価値演算部は、評価値として、一定時間内における実効値の分布の尖度および歪度に基づいた値を演算する。
【選択図】図6

Description

この発明は、装置を構成する機械要素の状態を監視する状態監視システムに関し、特に、風力発電装置を構成する機械要素の状態を監視する状態監視システムに関する。
風力発電装置においては、風力を受けるブレードに接続される主軸を回転させ、増速機により主軸の回転を増速させた上で発電機のロータを回転させることによって発電が行なわれる。主軸、増速機の回転軸、発電機の回転軸の各々は、転がり軸受によって回転自在に支持されており、そのような軸受の異常を診断する状態監視システム(CMS:Condition Monitoring System)が知られている。このような状態監視システムにおいては、軸受に固設された振動センサにより測定される振動波形データを用いて、軸受に損傷が発生しているか否かが診断される。
このような振動波形データの変化を検知する方法としては、振動波形データの実効値を算出し、算出した実効値のトレンドデータの変化を検知する方法が知られている。トレンドデータの変化を検知したことをトリガとして振動波形データの計測を開始することにより、計測された振動波形データを用いて機械要素の異常を診断することができる。
そのような方法の1つとして、前回周期で取得された実効値と今回周期で取得された実効値との差分値を計算し、その差分値が閾値を超えたときに、振動波形データの変化を検知する方法がある。たとえば、特開2012−252651号公報(特許文献1)には、発電プラントから送信されるプロセスデータの前回周期と今回周期の差分を抽出するように構成された監視装置が開示されている。
特開2012−252651号公報
しかしながら、上述した、差分値が閾値を超えるか否かに基づいてトレンドデータの変化を検知する方法では、差分値が大きくなるまで該変化を検知することが難しいという問題があった。これは、主軸ならびに増速機および発電機の回転速度によって各々の軸受の振動の大きさが異なるため、振動波形データに重畳するノイズの影響も回転速度によって異なってくることに起因する。そのため、トレンドデータの変化を検知するためには、ノイズに起因する差分値よりも大きい値に閾値を設定する必要がある。しかしながら、閾値を大きい値に設定すると、トレンドデータに変化が起きていても、該変化に起因する差分値が閾値を超えるまでは該変化を検知することができない場合が起こり得る。このことにより、たとえば、軸受の損傷に起因してトレンドデータが変化した場合、重大な故障に発展するまで該変化を検知できない可能性がある。その結果、故障の兆候である軸受の損傷を早期に検出することが困難となる。
また、実効値のトレンドデータの分布範囲(拡がり)を表す数値レンジは、主軸等の回転速度やノイズの影響度などによって異なるため、差分値の数値レンジもトレンドデータ間で異なったものとなる。その結果、トレンドデータの有意な変化を検知するためには、トレンドデータの数値レンジに応じて閾値を設定し直す必要があった。すなわち、トレンドデータの数値レンジが小さければ閾値を小さい値に設定する一方で、トレンドデータの数値レンジが大きければ閾値を大きい値に設定することが求められていた。このように、トレンドデータの変化の検知感度を確保するためには、トレンドデータごとにその数値レンジに応じた閾値を個別に設定しなければならないという問題があった。
そこで、この発明は、かかる課題を解決するためになされたものであり、その目的は、振動波形のトレンドデータの変化の検知感度の向上を簡便に実現する状態監視システムおよび風力発電装置を提供することである。
この発明のある局面に従えば、状態監視システムは、振動センサと、処理装置とを備える。状態監視システムは、装置を構成する機械要素の状態を監視する。振動センサは、機械要素の振動波形を計測する。処理装置は、評価値演算部と診断部とを含み、振動波形の変化を検知する。評価値演算部は、一定時間内に振動センサから出力される振動波形データの実効値を特徴付ける評価値を、時間的に連続して演算する。検知部は、評価値の推移に基づいて振動波形の変化を検知する。さらに、評価値演算部は、評価値として、一定時間内における実効値の分布の尖度および歪度に基づいた値を演算する。
この発明によれば、振動波形のトレンドデータの変化の検知感度の向上を簡便に実現する状態監視システムおよび風力発電装置を提供することができる。
本発明の実施の形態に係る状態監視システムが適用された風力発電装置の構成を概略的に示した図である。 図1に示したデータ処理装置の構成を機能的に示す機能ブロック図である。 振動波形データの差分値の時間的変化の一例を示した図である。 尖度と歪度の定義を説明する図である。 データのトレンド変化時の分布の概念図である。 図3に示した振動波形データ例の評価値の時間的変化を示した図である。 本発明の実施の形態に係る状態監視システムにおける振動波形データの変化を検知するための制御処理を説明するフローチャートである。
以下、本発明の実施の形態について、図面を参照しながら詳細に説明する。なお、図中の同一または相当部分には同一符号を付してその説明は繰返さない。
図1は、この発明に係る状態監視システムが適用された風力発電装置の構成を概略的に示した図である。図1を参照して、風力発電装置10は、主軸20と、ブレード30と、増速機40と、発電機50と、主軸用軸受(以下、単に「軸受」と称する。)60と、振動センサ70と、データ処理装置80とを備える。増速機40、発電機50、軸受60、振動センサ70およびデータ処理装置80は、ナセル90に格納される。ナセル90は、タワー100によって支持される。
主軸20は、ナセル90内に進入して増速機40の入力軸に接続され、軸受60によって回転自在に支持される。そして、主軸20は、風力を受けたブレード30により発生する回転トルクを増速機40の入力軸へ伝達する。ブレード30は、主軸20の先端に設けられ、風力を回転トルクに変換して主軸20に伝達する。
軸受60は、ナセル90内において固設され、主軸20を回転自在に支持する。軸受60は、転がり軸受によって構成され、たとえば、自動調芯ころ軸受や円すいころ軸受、円筒ころ軸受、玉軸受等によって構成される。なお、これらの軸受は、単列のものでも複列のものでもよい。
振動センサ70は、軸受60に固設される。振動センサ70は、軸受60の振動波形を計測し、計測した振動波形データをデータ処理装置80へ出力する。振動センサ70は、たとえば、圧電素子を用いた加速度センサによって構成される。
増速機40は、主軸20と発電機50との間に設けられ、主軸20の回転速度を増速して発電機50へ出力する。一例として、増速機40は、遊星ギヤ、中間軸および高速軸等を含む歯車増速機構によって構成される。なお、特に図示はしないが、増速機40内にも、複数の軸を回転自在に支持する複数の軸受が設けられている。
発電機50は、増速機40の出力軸に接続され、増速機40から受ける回転トルクによって発電する。発電機50は、たとえば、誘導発電機によって構成される。なお、発電機50内にも、ロータを回転自在に支持する軸受が設けられている。
データ処理装置80は、ナセル90内に設けられ、軸受60の振動波形データを振動センサ70から受ける。データ処理装置80は、予め設定されたプログラムに従って、軸受60の振動波形データの変化を検知する。また、風力発電装置10外の解析部180および通知部170に振動波形データを送信する(図2参照)。
図2は、図1に示したデータ処理装置80の構成を機能的に示す機能ブロック図である。図2を参照して、データ処理装置80は、ローパスフィルタ(以下、「LPF(Low Pass Filter)」と称する。)110と、実効値演算部120と、記憶部130と、評価値演算部140と、検知部150と、閾値設定部160とを含む。
LPF110は、軸受60の振動波形データを振動センサ70から受ける。LPF110は、その受けた振動波形データにつき、予め定められた周波数(たとえば、400Hz)よりも低い信号成分を通過させ、高周波成分を遮断する。
実効値演算部120は、軸受60の振動波形データをLPF110から受ける。実効値演算部120は、軸受60の振動波形データの実効値(「RMS(Root Mean Square)値」とも称される。)を算出し、その算出された振動波形データの実効値を記憶部130へ出力する。
記憶部130は、実効値演算部120により演算された軸受60の振動波形データの実効値を時々刻々と格納する。記憶部130は、たとえば、読み書き可能な不揮発性のメモリ等によって構成される。
記憶部130は、少なくとも一定時間内(たとえば7日間)における軸受60の振動波形データの実効値を格納するように構成される。記憶部130は、たとえば、所定の時間間隔(たとえば2時間)で、軸受60の振動波形データを実効値演算部120から受けると、一定時間内における振動波形データの実効値のうちの最も古い振動波形データの実効値を消去するとともに、新たに入力された振動波形データの実効値を追加するように構成される。
すなわち、記憶部130は、所定の時間間隔で、一定時間内における軸受60の振動波形データの実効値を更新する。後述するように、記憶部130に格納された一定時間内における軸受60の振動波形データの実効値が読み出され、その読み出された実効値を用いて振動波形データの変化が検知される。また、記憶部130は、後述する解析部180に振動波形データの実効値を出力する。
評価値演算部140は、一定時間内における軸受60の振動波形データの実効値を記憶部130から読み出すと、読み出した一定時間内における振動波形データの実効値を特徴付ける評価値を演算する。評価値演算部140は、評価値を時間的に連続して演算するように構成される。すなわち、評価値演算部140は、所定の時間間隔で、評価値を更新する。評価値演算部140における評価値の演算の詳細については後述する。
閾値設定部160は、検知部150において振動波形データの変化を検知するために用いられる閾値を設定するために用いられる。閾値設定部160は、設定された閾値を検知部150へ出力する。なお、閾値設定部160における閾値の設定はユーザが任意に決定してもよいし、振動波形データに基づいて決定されるように構成してもよい。
検知部150は、評価値を評価値演算部140から受け、閾値を閾値設定部160から受ける。検知部150は、評価値と閾値とを比較することにより、振動波形データの変化を検知する。具体的には、評価値が閾値より大きい場合、検知部150は振動波形データの変化を検知する。一方、評価値が閾値以下である場合、検知部150は振動波形データの変化を検知しない。また、検知部150は、検知結果を解析部180および通知部170に出力する。
通知部170は、遠隔地にいるユーザに検知結果を、たとえば視覚的手段や音声などの方法で通知する。
解析部180は、検知部150から振動波形データの変化が検知したという情報が入力された場合、当該検知をトリガとして振動波形データの計測を開始する。具体的には、解析部180は、トリガ発生時点以降記憶部130に格納される振動波形データの実効値を読み出す。解析部180は読み出した振動波形データの実効値を解析することにより、軸受60の異常を診断する。該振動波形データの解析により、風力発電装置10の振動波形データの変化の原因(たとえば、軸受60の損傷および環境の変化など)等をより詳細に調べることができる。解析部180による振動波形データの解析は、自動解析用のプログラムを用いてもよいし、ユーザが手動で行ってもよい。
以下、検知部150における振動波形データの変化を検知する方法を説明する。最初に、図3を参照して、比較例として、実効値の差分値を用いて振動波形データの変化を検知する方法を説明する。
図3は、記憶部130に格納される軸受60の振動波形データの実効値の時間的変化の一例、および、該実効値の差分値の時間的変化を示した図である。なお、本願明細書において、実効値の差分値とは、今回格納された実効値から前回格納された実効値を減算した値である。
図3を参照して、実効値は時間的に変動している。実効値の時系列変化の傾向(以下、トレンドとも称する)としては、時刻t1よりも前の期間では、実効値の数値レンジが一定範囲に収まっている。これに対して、時刻t1よりも後の期間では、実効値の変動が大きくなっている。このときの実効値の数値レンジは上限側が高くなるように広がっている。その結果、時刻t1よりも前に比べて、数値レンジの中央が上昇している。
このように図3の例では、図中に丸印で囲んだ領域42に示されるように、時刻t1付近で実効値のトレンドに変化が生じている。なお、この時刻t1付近でのトレンド変化は、たとえば、時刻t1以降における大きな実効値の変化で表される計測対象の状態変化、または、風力発電装置10の設置場所における風の吹き方を示す風況などの環境の変化を示す。したがって、このトレンドデータの変化を検知することが求められる。
ここで、時刻t1におけるような、上記計測対象の状態変化、もしくは、環境の変化を示すトレンドデータの変化を、差分値に基づいて検知することを考える。図3では、閾値Tdは、時刻t1よりも前の期間における差分値の数値レンジよりも高い値に設定されている。
図3に示すように、時刻t1では差分値は閾値Tdを下回っている。そのため、差分値を用いては時刻t1におけるトレンドデータの変化を検知することができない。なお、時刻t1よりも遅い時刻t2にて、差分値が閾値Tdを超えている。よって、時刻t1よりも遅い時刻t2にてトレンドデータの変化が検知されることとなる。
図3において、差分値に基づいてトレンドデータの変化が検知される時刻t2と、実際に計測対象の状態変化もしくは環境の変化に対応してトレンドデータが変化した時刻t1とのずれを小さくするためには、閾値Tdを低下させることが必要となる。しかしながら、時刻t1における差分値は、時刻t1よりも前の時刻における差分値と同程度の値を示すことから、閾値Tdを下げることで、時刻t1よりも前のトレンドデータのトレンド変化(以下、トレンドデータの変化とも称する)が起こっていない期間に、誤ってトレンドデータの変化が検知されることとなる。
このように実効値の差分値を用いてトレンドデータの変化を検知する方法では、誤検知を防ぐ観点から、差分値の数値レンジによって閾値Tdが制約を受ける。その結果、当該方法では、差分値の数値レンジが十分に大きくなるまでトレンドデータの変化を捉えることができないという問題がある。
また、実効値の数値レンジによって差分値の数値レンジも異なるため、トレンドデータの変化を検知できる感度(以下、検知感度とも称する)を確保するためには、トレンドデータごとにその数値レンジに応じた閾値を個別に設定しなければならないという問題がある。
そこで、本実施の形態では、一定時間内における振動波形データの実効値を特徴付ける評価値を時間的に連続して演算し、その演算された評価値の推移に基づいて振動波形の変化を検知するように構成する。上記構成において、評価値は、一定時間内における実効値の分布の尖度(Kurtosis)および歪度(Skewness)に基づいた値とする。
尖度および歪度は、分布の形を表す統計値であり、差分値とは異なり、無次元化された値である。そのため、実効値の数値レンジとは関係なく、一定時間内における実効値の分布の特徴を表すことができる。したがって、様々な実効値の数値レンジに対して、様々な閾値を設定する必要がなくなり、閾値を統一することができる。これにより、トレンドデータの変化の検知感度の向上を簡便に実現することができる。
以下、図4および5を参照して、尖度および歪度の定義を説明する。
図4(A)は、尖度の定義を説明する図である。一定時間内における実効値の分布の形を表す統計値として、尖度は、該分布の尖り具合を示す。一般に、尖度は、正規分布のときに0となり(グラフ32参照)、正規分布に比べて厚い裾を持つときに正の値となり(グラフ33参照)、正規分布に比べて薄い裾を持つときに負の値となる(グラフ31参照)傾向がある。なお、本実施の形態で扱うデータにおいては、その分布の尖度はほぼ正となる。すなわち、本実施の形態においては、尖度の絶対値が小さくなるほど、平均値の周りにデータが集中していることを表す。
より詳細には、分布の裾の厚さとは、分布の平均値の周りにデータが集中している度合いを表す。以下の説明では、一定時間内における実効値のデータの個数をnとして、それらをx,x,・・・xと表すこととする。実効値のデータx,x,・・・xの分布において、平均値をμ、標準偏差をσ、尖度をKとすると、μ、σおよびKはそれぞれ、以下の式(1)〜(3)で与えられる。
Figure 2018124117
Figure 2018124117
Figure 2018124117
図4(B)は、歪度の定義を説明する図である。歪度は、分布の左右対称性(歪み)を表す。歪度は、分布が左右対称のときに0となり(グラフ35参照)、分布が左右対称のときに比べ、分布が負側(左側)に偏っているときに正の値となり(グラフ34参照)、分布が左右対称のときに比べて分布が正側(右側)に偏っているときに負の値となる(グラフ36参照)。すなわち、歪度の絶対値が大きくなるほど、データの分布が正または負に偏っていることを表す。
より詳細には、歪度をSとすると、Sは以下の式(4)で与えられる。
Figure 2018124117
次に、図5を参照して、本明細書の振動波形データのような時系列データにおいて、トレンド変化が起こったときに出現し得る尖度および歪度の変化を説明する。図5は、データのトレンド変化時の分布の概念図である。
まず、データの上昇トレンドが起こっているときを考える。上昇トレンドが発生すると、発生前と比較して、正側(右側)に外れ値が出始めると考えられる。そのため、上昇トレンドが起こったときのデータの分布は、それより前の時間の分布(グラフ38参照)に比べ、より正側に範囲が拡大する(裾が厚くなる)と考えられる(グラフ39参照)。すなわち、尖度の値は正に大きくなり、歪度の値は正に大きくなると考えられる。
次に、データの下降トレンドが起こっているときを考える。下降トレンドが発生すると、発生前と比較して負側(左側)に外れ値が出始めると考えられる。そのため、下降トレンドが起こったときのデータの分布は、それより前の時間の分布(グラフ38参照)に比べ、より負側に範囲が拡大する(裾が厚くなる)と考えられる(グラフ37参照)。すなわち、尖度の値は正に大きくなり、歪度の値は負に大きくなると考えられる。
すなわち、本明細書の振動波形データのような時系列データにおいて、上昇または下降トレンドが起こったときには、尖度の値は正に大きくなり、歪度の値は正または負に大きくなると考えられる。換言すると、尖度および歪度の絶対値の値は、共に大きくなる。そこで、本実施の形態では、振動波形データの変化を検知するための評価値として、一定時間内における実効値の分布の尖度および歪度に基づいた値を演算するように構成される。より好ましくは、評価値として、尖度と歪度の積の絶対値を演算するように構成される。
本実施の形態では、上述の通り、評価値として、一定時間内における実効値の分布の尖度Kおよび歪度Sの積の絶対値を演算する。評価値をPとすると、評価値Pは以下の式(5)で与えられる。
Figure 2018124117
式(5)から分かるように、評価値Pは、尖度Kが大きくなるほど大きな値となる。また、評価値Pは、歪度Sの絶対値が大きくなるほど大きな値となる。したがって、一定時間内における実効値の分布において、負側(左側)にデータの裾が厚くなる場合(図5のグラフ37参照)、もしくは、正側(右側)にデータの裾が厚くなる場合(図5のグラフ39参照)、評価値Pは値が大きくなる。
図6は、図3に示した実効値の時間的変化に対する、評価値Pの時間的変化を示す図である。
図6を参照して、評価値Pは、時刻t1付近において急激に増加している。これは、時刻t1付近で、トレンドデータの変化に対応して、一定時間内における実効値の分布にも変化が生じたことを示している。詳細には、上述のように、一定時間内における実効値の分布において、負側または正側にデータが集中する歪みが生じたことを示している。
図6に示すように、時刻t1での評価値Pの値と同程度の値の閾値Tpを設定することで、時刻t1付近でのトレンドデータの変化を検知することができる。評価値Pは尖度Kと歪度Sの積の絶対値であるため、尖度Kおよび歪度Sと同じく、無次元化された値である。すなわち、様々な実効値の数値レンジに対して、統一した閾値Tpを設定することができる。これによれば、差分値によっては検知することが困難な変化も検知することができる。この結果、トレンドデータの変化の検知感度を向上させることが可能となる。
図7は、本実施の形態に係る状態監視システムにおける振動波形の変化を検知するための制御処理を説明するフローチャートである。図7に示される制御処理は、データ処理装置80により、所定の時間間隔で繰り返し実行される。
図7を参照して、データ処理装置80は、ステップS01において、軸受60の振動波形データを振動センサ70から受ける。続いて、ステップS02において、LPF110は、軸受60の振動波形データに対してフィルタ処理を実行する。
次に、ステップS03において、フィルタ処理が施された軸受60の振動波形データをLPF110から受けると、データ処理装置80は実効値演算部120において、軸受60の振動波形データの実効値を算出する。ステップS04において、データ処理装置80は、実効値演算部120において算出された振動波形データの実効値を記憶部130に格納する。
次に、ステップS05において、データ処理装置80は実効値演算部120において、所定の条件を満たす実効値を全実効値データから抽出する。具体的には、データ処理装置80は記憶部130に記憶された実効値のうち、所定の期間分の最新のデータにおいて、発電機出力が規定値以上であり、回転速度が規定値以上であるという条件を満たすデータのみを抽出する。
データ処理装置80の評価値演算部140は、ステップS06において、ステップS05において抽出された実効値のデータ数が規定個数以上であるか否かを判定する。ステップS05において抽出された振動波形データの実効値のデータ数が規定個数未満である場合(S06にてNO)、以降の処理S07〜S09はスキップされて、処理がメインルーチンに戻される。
一方、ステップS05において抽出されたデータ数が規定個数以上である場合(S06にてYES)、ステップS07に処理が進められ、データ処理装置80は評価値演算部140において、抽出された規定個数の振動波形データの実効値の評価値Pを演算する。ここで、評価値Pは、上述の通り、実効値の尖度Kと歪度Sとの積の絶対値である。
ステップS08において、データ処理装置80は、検知部150により、演算された評価値Pと閾値Tpとを比較する。評価値Pが閾値Tp未満である場合(S08にてNO)、データ処理装置80は、以降の処理S09をスキップして処理をメインルーチンに戻す。一方、評価値Pが閾値Tp以上である場合(S08にてYES)、データ処理装置80は、ステップS09において、検知部150は通知部170および解析部180に検知結果を出力する(図2参照)。そして、通知部170は、ユーザに振動波形の変化の検知を通知する。解析部180は、該検知後に記憶部130に格納される振動波形データの実効値を解析することにより、風力発電装置10の異常を診断する。この結果、振動波形の変化の原因となる事象(たとえば、重大な故障の予兆)を、早期に把握することが可能である。
以上のように、本実施の形態によれば、一定時間内における軸受60の振動波形データの実効値を特徴付ける評価値として、一定時間内における実効値の分布の尖度および歪度に基づいて算出する。このようにすると、トレンドデータの数値レンジを考慮して閾値を設定することが不要となるため、差分値によっては検知することが困難な変化も検知することができる。この結果、トレンドデータの変化の検知感度を向上させることが可能となる。具体的には、差分値では検知することが困難であった、重大な故障の予兆である機械要素の損傷などを検知することができる。
好ましくは、評価値として、一定時間内における実効値の分布の尖度および歪度の積の絶対値を用いる。このようにすると、一定時間内における実効値の分布において、正側または負側に裾が厚くなる変化が生じた場合には、当該変化を反映するように評価値も変化する。よって、この評価値の変化を捉えることで、トレンドデータの変化を検知することができる。
今回開示された実施の形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施の形態の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
10 風力発電装置、20 主軸、30 ブレード、40 増速機、42 実効値変化、50 発電機、60 軸受、70 振動センサ、80 データ処理装置、90 ナセル、100 タワー、120 実効値演算部、130 記憶部、140 評価値演算部、150 検知部、160 閾値設定部、170 通知部、180 解析部、P 評価値、Td,Tp 閾値。

Claims (4)

  1. 装置を構成する機械要素の状態を監視する状態監視システムであって、
    前記機械要素の振動波形を計測するための振動センサと、
    前記振動波形の変化を検知するための処理装置とを備え、
    前記処理装置は、
    一定時間内に前記振動センサから出力される振動波形データの実効値を特徴付ける評価値を、時間的に連続して演算する評価値演算部と、
    前記評価値の推移に基づいて、前記振動波形の変化を検知する検知部とを含み、
    前記評価値演算部は、前記評価値として、前記一定時間内における前記実効値の分布の尖度および歪度に基づいた値を演算するように構成される、状態監視システム。
  2. 前記評価値演算部は、前記評価値として、前記一定時間内における前記実効値の分布の尖度と歪度の積の絶対値を演算するように構成される、請求項1に記載の状態監視システム。
  3. 前記検知部は、前記評価値が閾値を超えたときに、前記振動波形の変化を検知するように構成される、請求項1または2に記載の状態監視システム。
  4. 請求項1〜3のいずれか1項に記載の状態監視システムを備える、風力発電装置。
JP2017015297A 2017-01-31 2017-01-31 状態監視システムおよび風力発電装置 Pending JP2018124117A (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2017015297A JP2018124117A (ja) 2017-01-31 2017-01-31 状態監視システムおよび風力発電装置
PCT/JP2018/001794 WO2018142986A1 (ja) 2017-01-31 2018-01-22 状態監視システムおよび風力発電装置
CN201880009272.4A CN110234972B (zh) 2017-01-31 2018-01-22 状态监视系统及风力涡轮机
US16/481,797 US20200025648A1 (en) 2017-01-31 2018-01-22 Condition monitoring system and wind turbine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017015297A JP2018124117A (ja) 2017-01-31 2017-01-31 状態監視システムおよび風力発電装置

Publications (1)

Publication Number Publication Date
JP2018124117A true JP2018124117A (ja) 2018-08-09

Family

ID=63039579

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017015297A Pending JP2018124117A (ja) 2017-01-31 2017-01-31 状態監視システムおよび風力発電装置

Country Status (4)

Country Link
US (1) US20200025648A1 (ja)
JP (1) JP2018124117A (ja)
CN (1) CN110234972B (ja)
WO (1) WO2018142986A1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2887298T3 (es) * 2019-01-08 2021-12-22 Nordex Energy Se & Co Kg Método para el funcionamiento de una turbina eólica
US11939955B2 (en) 2019-03-28 2024-03-26 Ntn Corporation Condition monitoring system
GB2588975B (en) * 2019-11-18 2022-06-15 Octonion Sa A method of determining changes in stationary states of a signal
JP7401327B2 (ja) * 2020-01-31 2023-12-19 高砂熱学工業株式会社 診断装置、診断方法、診断プログラムおよび診断システム
JP7379241B2 (ja) * 2020-03-25 2023-11-14 三菱重工業株式会社 回転機の診断監視装置及び方法
WO2022036710A1 (en) * 2020-08-21 2022-02-24 Abb Schweiz Ag Method of indication selection for a fault diagnosis of a gearbox
US11401918B2 (en) * 2020-12-02 2022-08-02 General Electric Renovables Espana, S.L. System and method for controlling a wind turbine
US11708815B2 (en) 2021-02-08 2023-07-25 General Electronic Company System and method for controlling a wind turbine
US11774324B2 (en) 2021-03-12 2023-10-03 General Electric Renovables Espana, S.L. System and method for detecting actual slip in a coupling of a rotary shaft
US11913429B2 (en) 2021-04-29 2024-02-27 General Electric Renovables Espana, S.L. System and method for slip detection and surface health monitoring in a slip coupling of a rotary shaft

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5943634A (en) * 1996-05-14 1999-08-24 Csi Technology, Inc. Vibration data analysis based on time waveform parameters
JP4232162B2 (ja) * 2004-12-07 2009-03-04 三菱電機株式会社 圧縮機検査装置
AU2010220551B2 (en) * 2009-03-02 2013-10-10 Suzlon Energy Gmbh Method for monitoring wind turbines
JP5725833B2 (ja) * 2010-01-04 2015-05-27 Ntn株式会社 転がり軸受の異常診断装置、風力発電装置および異常診断システム
CN101858778A (zh) * 2010-05-28 2010-10-13 浙江大学 基于振动监测的风力发电机组故障自动诊断方法
JP2011252762A (ja) * 2010-06-01 2011-12-15 Jfe Advantech Co Ltd 軸受状態監視方法及び軸受状態監視装置
JP2012008030A (ja) * 2010-06-25 2012-01-12 Toshiba Plant Systems & Services Corp 回転機軸受診断装置
GB201020381D0 (en) * 2010-12-02 2011-01-12 Crane John Uk Ltd Component failure detection system
WO2013133002A1 (ja) * 2012-03-08 2013-09-12 Ntn株式会社 状態監視システム
CN106471247B (zh) * 2014-06-24 2019-06-28 Ntn株式会社 状态监视系统和使用该系统的风力发电系统
EP2988187B1 (en) * 2014-08-22 2017-03-29 ABB Schweiz AG A method for assessing the condition of rotating machinery connected to an electric motor
CN104748992A (zh) * 2015-03-13 2015-07-01 阿海珐有限公司 用于监测和评估旋转的机器的状态的装置和方法

Also Published As

Publication number Publication date
CN110234972A (zh) 2019-09-13
CN110234972B (zh) 2021-09-21
WO2018142986A1 (ja) 2018-08-09
US20200025648A1 (en) 2020-01-23

Similar Documents

Publication Publication Date Title
JP2018124117A (ja) 状態監視システムおよび風力発電装置
CN108780025B (zh) 状态监视系统和风力发电装置
JP5561835B2 (ja) 物品監視のためのシステムおよび方法
JP6820771B2 (ja) 状態監視システムおよび風力発電装置
WO2017159784A9 (ja) 状態監視システムおよび風力発電装置
US11047768B2 (en) Abnormality diagnosis device and abnormality diagnosis method
JP6158100B2 (ja) 風力発電機の異常程度判定システムまたは風力発電機の異常程度判定方法
JP2015527586A (ja) 機械システムの変化を検出し、前記機械システムに関連する制限値を、機械システムの現在の条件を反映するように適合させるためのモニタリングシステム及び方法
JP6824076B2 (ja) 状態監視システムおよび風力発電装置
US10895873B2 (en) Machine health monitoring of rotating machinery
JP2019052964A (ja) 状態監視システムおよびデータ処理装置
CN113383215A (zh) 用于模式匹配轴承振动诊断的系统和过程
EP3349003B1 (en) Rotating machine abnormality detection device, rotating machine abnormality detection method, and rotating machine
JP2010175446A (ja) 状態診断装置
JP2016170085A (ja) 異常診断装置及び異常診断方法
US11746752B2 (en) Method and system for detecting machine defects
JP6869156B2 (ja) 状態監視装置および状態監視方法
JP5751606B2 (ja) 機械設備における異常診断システム
JP6808405B2 (ja) 回転機器の監視装置
JP6639288B2 (ja) 状態監視システムおよび風力発電装置
CN119768673A (zh) 监视装置及监视方法
JP6736987B2 (ja) 風力発電装置の回転部品の状態監視装置
US20260043712A1 (en) Monitoring device and monitoring method
JP7600857B2 (ja) 軸受の異常診断装置、異常診断方法、およびプログラム
WO2015151237A1 (ja) 風車の診断装置