JP6619908B1 - 診断装置、診断方法および診断プログラム - Google Patents

Abstract

【課題】回転機械を有する装置の測定電流に基づく異常検知精度が向上された診断装置を提供する。【解決手段】診断装置は、回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断装置であって、前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するよう構成された実効値取得部と、取得された複数の前記実効値の分布状態を表す対象分布情報を算出するよう構成された分布情報算出部と、算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するよう構成された検知部と、を備える。【選択図】 図２

Description

本開示は、回転機械を有する装置の診断技術に関する。

例えば特許文献１〜２には、電動機などの回転機械の異常検知を回転機械の電流値に基づいて行う手法が開示されている。具体的には、特許文献１は、ファンモータを流れる電流の実効値と、この実効値が得られた際の電源周波数および電源電圧に応じて設定された閾値とを比較して、ファンモータのモータ異常の判定（検知）を行うことを開示する。また、特許文献２は、設備状態が変化すると、その電流波形の振幅確率密度関数も変化することから、電動機の定格電流の基準の正弦波から求めた参照振幅確率密度関数と、電動機の稼働時に計測した電流波形から求めた点検時振幅確率密度とのカルバック・ライブラー距離（情報距離）を算出することで、電動機の異常の検知を行うことを開示する。

特開２０１３−０５０２９４号公報 特開２０１１−２５７３６２号公報

電動機の電流値は負荷やモータの特性により大きく異なる。このため、電動機の異常検知を実効値の閾値監視で実施する場合、負荷やモータの特性ごとに閾値を設定する必要があるが、実際にはその数が膨大となるなど、個別に設定することは困難である。また、定格電流の正弦波を基準として異常度合いの計算を行う方法についても同様に、電流値が負荷やモータによって大きく異なるため、実際には正常と判定すべきところを異常と判定知してしまう可能性がある。

上述の事情に鑑みて、本発明の少なくとも一実施形態は、回転機械を有する装置の測定電流に基づく異常検知精度が向上された診断装置を提供することを目的とする。

本発明の少なくとも一実施形態に係る診断装置は、
回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断装置であって、
前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するよう構成された実効値取得部と、
取得された複数の前記実効値の分布状態を表す対象分布情報を算出するよう構成された分布情報算出部と、
算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するよう構成された検知部と、を備える。

また、本発明の少なくとも一実施形態に係る診断方法は、
回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断方法であって、
前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するステップと、
取得された複数の前記実効値の分布状態を表す対象分布情報を算出するステップと、
算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するステップと、を備える。

また、本発明の少なくとも一実施形態に係る診断プログラムは、
回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断プログラムであって、
コンピュータに、
前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するよう構成された実効値取得部と、
取得された複数の前記実効値の分布状態を表す対象分布情報を算出するよう構成された分布情報算出部と、
算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するよう構成された検知部と、を実現させるためのプログラムである。

本発明の少なくとも一実施形態によれば、回転機械を有する装置の測定電流に基づく異常検知精度が向上された診断装置が提供される。

本発明の一実施形態に係る診断装置および診断対象機械を含む診断系を概略的に示す図である。 本発明の一実施形態に係る診断装置の機能を概略的に示す図である。 本発明の一実施形態に係るある種類の異常についての電流の実効値の正常時と異常時とを対比して示す図である。 本発明の一実施形態に係る他の種類の異常についての電流の実効値の正常時と異常時とを対比して示す図である。 本発明の一実施形態に係る回転機械が異常時の実効値の確率密度分布を示す図である。 本発明の一実施形態に係る回転機械が正常時の実効値の確率密度分布を示す図である。 本発明の一実施形態に係る複数の診断対象機械の正常電流波形を対比して示す図である。 本発明の一実施形態に係る診断の流れを説明するための図である。 本発明の一実施形態に係る診断方法を示す図である。

以下、添付図面を参照して本発明の幾つかの実施形態について説明する。ただし、実施形態として記載されている又は図面に示されている構成部品の寸法、材質、形状、その相対的配置等は、本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
例えば、「ある方向に」、「ある方向に沿って」、「平行」、「直交」、「中心」、「同心」或いは「同軸」等の相対的或いは絶対的な配置を表す表現は、厳密にそのような配置を表すのみならず、公差、若しくは、同じ機能が得られる程度の角度や距離をもって相対的に変位している状態も表すものとする。
例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
例えば、四角形状や円筒形状等の形状を表す表現は、幾何学的に厳密な意味での四角形状や円筒形状等の形状を表すのみならず、同じ効果が得られる範囲で、凹凸部や面取り部等を含む形状も表すものとする。
一方、一の構成要素を「備える」、「具える」、「具備する」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。

図１は、本発明の一実施形態に係る診断装置１および診断対象機械９を含む診断系６を概略的に示す図である。図２は、本発明の一実施形態に係る診断装置１の機能を概略的に示す図である。

図１に示すように、診断装置１は、回転機械９１を有する診断対象機械９を、その回転機械９１の回転時の測定電流に基づいて診断する装置である。回転機械９１は、例えば、供給された電気エネルギーを回転力に変換するモータやポンプなどの電動機、力学的エネルギー（機械エネルギー）により回転子を回転させることで、力学的エネルギーを電気エネルギーに変換する発電機、このような発電機を回転駆動させるガスタービンや蒸気タービンなどである。また、診断対象機械９は、例えば、回転機械９１を主要なコンポーネントとする電動機や発電機自体、このような電動機や発電機を備える例えば車両などの機械（装置）、蒸気タービンやガスタービンなどの回転軸に発電機の駆動軸が連結されて発電を行う発電装置などである。なお、発電装置は、例えば火力発電プラントやガスタービン・コンバインドサイクル（ＧＴＣＣ）発電プラントなどの発電プラントを構成するものであっても良い。

例えば、図１に示すように、診断装置１は、例えば診断対象機械９が備える回転機械（電動機）に供給される交流電流や、回転機械（発電機）から出力される交流電流など、診断対象機械９における電流Ｉを測定可能な電流測定装置７に接続されることで、この電流Ｉをリアルタイムに監視するように構成されても良い。図１に示す実施形態では、診断対象機械９は、回転機械９１である電動機を備える機械である。そして、この電動機と、電動機を回転駆動するための電力を入力（供給）する例えば配電盤、分電盤あるいは制御盤などとなる電気盤８とを接続する配線（電源線８１）を流れる交流の電流Ｉを、小型クランプを用いて電流測定装置７で測定するようになっている。また、電流測定装置７は、測定した電流Ｉ（測定電流）の測定値Ｉｍを、有線または無線などの通信媒体７１（図１では通信線）を介して診断装置１に送信（入力）するようになっている。

例えば診断対象機械９が、回転機械９１として、発電機、および、この発電機を回転駆動するガスタービンや蒸気タービンを備える場合には、発電機から電気盤８などに出力される電流Ｉを測定電流としても良い。つまり、診断対象機械９が、関連して回転する２以上の回転機械９１を備える場合には、そのうちの一つの回転機械９１から出力される電流Ｉ、あるいは、そのうちの一つの回転機械９１に入力（供給）される電流Ｉを測定電流としても良い。例えば、蒸気タービン等による発電により生じる電流Ｉを測定することで、蒸気タービン等および発電機の両方が正常であるか否かを診断することが可能となり、少なくとも一方に異常が生じている場合の検知が可能となる。

ただし、図１に示すような実施形態に本発明は限定されない。他の幾つかの実施形態では、診断対象機械９は、電流測定装置７に接続されていなくても良い。この場合、電流測定装置７による上記の測定によって、その測定値Ｉｍの時間推移のデータを記憶しておき、このデータを、例えばＵＳＢメモリ、ＳＤメモリ、コンパクトフラッシュ（登録商標）などの持ち運び可能な記憶媒体を介して診断装置１に入力しても良い。あるいは、この測定値Ｉｍの時間推移のデータを、通信ネットワークを介して診断装置１が受信するように構成しても良い。

以下、診断装置１について、図２を用いて詳細に説明する。
図２に示すように、診断装置１は、実効値取得部２と、分布情報算出部３と、検知部４と、を備える。この診断装置１は、コンピュータで構成されている。例えば、診断装置１は、図１に示すようなラップトップ型のコンピュータであっても良い。診断装置１が備える各機能部について、それぞれ説明する。

なお、コンピュータは、図示しないＣＰＵなどのプロセッサ１１や、ＲＯＭやＲＡＭといったメモリ（記憶装置１２）などを備えている。そして、主記憶装置にロードされたプログラム（診断プログラム１０）の命令に従ってプロセッサ１１が動作（データの演算など）することで、上記の各機能部を実現する。換言すれば、上記の診断プログラム１０は、コンピュータに上記の各機能部を実現させるためのソフトウェアであり、一時的な信号ではなく、コンピュータによる読み込みが可能で持ち運び可能な上記のような記憶媒体に記憶されても良い。

実効値取得部２は、上述したような診断対象機械９の測定電流の所定期間（以下、第１所定期間）での時間推移である電流波形Ｗ（以下、対象電流波形Ｗｔ）における規定周期数毎の実効値Ｉｅを取得するよう構成された機能部である。上記の規定周期数は１以上（１周期数以上）であり、第１所定期間は測定電流の規定周期数分の電流波形Ｗの複数分の時間以上を有する期間である。つまり、対象電流波形Ｗｔは、複数周期に渡る測定電流の時間推移で形成（構成）される。よって、対象電流波形Ｗｔにおいて時間軸に沿って並ぶ規定周期数分ごとの電流波形Ｗ（以下、単位電流波形Ｗｕ）に基づいて実効値Ｉｅがそれぞれ算出されることで、対象電流波形Ｗｔから複数の実効値Ｉｅが得られる。

例えば、上記の複数の実効値Ｉｅの算出は、上記の実効値取得部２が行っても良い。この場合には、電流測定装置７から診断装置１に入力される電流の測定値Ｉｍは電流の瞬時値である。そして、実効値取得部２は、対象電流波形Ｗｔの全体を取得した後、取得した対象電流波形Ｗｔを形成する複数の単位電流波形Ｗｕを、時間軸に沿って重複する部分がないように識別した上で、その各々の単位電流波形Ｗｕの実効値Ｉｅを算出しても良い。または、実効値取得部２は、対象電流波形Ｗｔの各部分を時間の経過に従って順次取得しつつ、得られた単位電流波形Ｗｕについての実効値Ｉｅの算出を順次行っても良い。

あるいは、複数の実効値Ｉｅの算出は、上述した電流測定装置７などの他の装置が行っても良い。この場合には、電流測定装置７から診断装置１に入力される電流の測定値Ｉｍは実効値Ｉｅである。そして、実効値取得部２は、第１所定期間で得られることが期待される数の電流の実効値Ｉｅが得られた場合に、後述する処理を行っても良い。

より詳細には、上記の実効値Ｉｅは、対象電流波形Ｗｔを形成する各々の単位電流波形Ｗｕからのサンプリングにより得られた複数の電流の瞬時値に基づいて算出される。つまり、１つの単位電流波形Ｗｕから、相互に異なる時間における電流の瞬時値が例えば等間隔でサンプリングされる。より具体的には、上記のサンプリングによるサンプリング数は、単位電流波形Ｗｕを構成する１以上の周期の１周期あたり９００以上１２００以下であるのが良い。これによって、必要な精度を得つつ、処理負荷が過大とならないようにできる。図１〜図２に示す実施形態では、上記のサンプリング数は、例えば１０９０など、１０７０〜１１１０の範囲としている。なお、上記の規定周期数は１周期としている。このサンプリング数を適切に設定することで、後述する対象分布情報に、回転機械の稼働状態（正常または異常）が適切に反映されるように図っている。

なお、既に説明したように、上記の測定電流は、診断対象機械９（回転機械９１）に供給される電力の電流を測定したものであっても良いし、診断対象機械９（回転機械９１）から出力される電力の電流を測定したものであっても良い。また、電流測定装置７を用いてリアルタイムに測定される電流の瞬時値を第１所定期間の間測定することで、対象電流波形Ｗｔを取得しても良いし、既に測定済みの対象電流波形Ｗｔのデータを取得しても良い。

分布情報算出部３は、上記の実効値取得部２によって取得された複数の実効値Ｉｅの分布状態を表す分布情報（以下、対象分布情報Ｄｔ）を算出するよう構成された機能部である。この分布情報は、確率密度関数などの確率分布であっても良いし、分布の形や、この分布の形を数値化した標準偏差や分散値などであっても良い。図１〜図２に示す実施形態では、対象分布情報Ｄｔの作成にあたって取得する実効値Ｉｅの数は４００以上としている。これによって、回転機械の稼働状態（正常または異常）が対象分布情報Ｄｔに適切に反映されるようにすることが可能である。

検知部４は、上記の分布情報算出部３によって算出された対象分布情報Ｄｔに基づいて、診断対象機械９の異常検知を実行するよう構成された機能部である。この異常検知（異常検知処理）は、診断対象機械９が正常か異常かを判定することであり、測定電流の実効値Ｉｅに影響を与えるような異常の発生を判定できる。例えば、ミスアライメント（芯ずれ）、ポンプにおいてパイプ内の圧力低下により液体が気化した気泡を生じるキャビテーション、ディスク円盤接触、ベルトの緩み、地絡などの異常の検知が可能である。

図１〜図２に示す実施形態では、図２に示すように、診断装置１は、電流測定装置７に接続されており、電流測定装置７からリアルタイムに電流の測定値Ｉｍが入力されるようになっており、これを実効値取得部２が取得するようになっている。また、診断装置１においては、実効値取得部２と分布情報算出部３とが接続されており、分布情報算出部３は、実効値取得部２から入力された複数の実効値Ｉｅに基づいて、対象分布情報Ｄｔを算出するようになっている。そして、検知部４は、分布情報算出部３に接続されており、分布情報算出部３から対象分布情報Ｄｔが入力されると、上記の異常検知を実行するようになっている。また、この検知部４は、出力装置１４に異常検知の結果を出力するようになっている。本実施形態では、出力装置１４はディスプレイなどの表示装置であるが、他の実施形態では、ユーザに異常検知の結果を知らせることができれば、出力装置１４は、表示装置以外であっても良く、表示装置、プリンタ、スピーカー、あるいはＬＥＤなどの発光装置のうちの少なくとも１つを含んでいても良い。

なお、診断対象機械９が蒸気タービン等と発電機とを備える場合には、蒸気タービン等が発電機を駆動させる。このような関係にあるため、例えば発電機の発電による電流の時間推移を対象電流波形Ｗｔとして診断を実行する場合に異常が判定されると、蒸気タービン等または発電機の少なくとも一方の異常が疑われることになる。よって、検知部４による診断を１次診断として、異常が検知された場合には詳細な調査を行っても良い。

上記の構成によれば、診断対象機械９が備える例えばタービンや発電機、モータなどの回転機械９１の回転時に、回転機械９１（例えば電動機）に供給される交流電流や、回転機械９１（例えば発電機）から出力される交流電流などの電流Ｉの瞬時値（測定電流）の時間推移である電流波形Ｗ（対象電流波形Ｗｔ）における規定周期数毎の電流波形Ｗ（後述する単位電流波形Ｗｕ）の実効値Ｉｅを取得する。そして、こうして取得された複数の実効値Ｉｅの、例えば確率密度関数や標準偏差などの分布状態を表す情報（対象分布情報Ｄｔ）に基づいて、診断対象機械９の異常検知を実行する。

本発明者らは、回転機械９１の異常の種類によっては、異常発生時において電流の瞬時値の分布状態に基づいて異常検知を実行した場合（特許文献２参照）よりも、電流の実効値の分布状態に基づいて異常検知を実行した場合の方が、異常判定の精度が高いことを見出している。よって、上記の対象電流波形Ｗｔに基づいて得られる複数の電流の実効値Ｉｅの分布状態（ばらつき）に基づいて回転機械９１の診断を実行することで、回転機械９１が正常であるか、異常であるかをより適切に診断することができる。

次に、上述した検知部４についての幾つかの実施形態について、詳細に説明する。
図３Ａは、本発明の一実施形態に係るある種類の異常についての電流の実効値の正常時と異常時とを対比して示す図である。図３Ｂは、本発明の一実施形態に係る他の種類の異常についての電流の実効値の正常時と異常時とを対比して示す図である。図４Ａは、本発明の一実施形態に係る回転機械が異常時の実効値Ｉｅの確率密度分布を示す図である。また、図４Ｂは、本発明の一実施形態に係る回転機械が正常時の実効値Ｉｅの確率密度分布を示す図である。なお、図３Ａおよび図３Ｂの縦軸および横軸のスケールは同じである。同様に、図４Ａおよび図４Ｂの縦軸および横軸のスケールは同じである。

幾つかの実施形態では、図２に示すように、診断装置１は、診断対象機械９が正常状態にある際の上記の測定電流の所定期間（以下、第２所定期間）での時間推移である電流波形Ｗ（以下、正常電流波形Ｗｂ）における規定周期数毎の実効値Ｉｅの分布状態を表す正常分布情報Ｄｂを記憶するよう構成された記憶部５を、さらに備えても良い。そして、検知部４は、対象分布情報Ｄｔと正常分布情報Ｄｂとに基づいて、上記の異常検知を実行する。正常分布情報Ｄｂの作成にあたって取得する実効値Ｉｅの数は、対象分布情報Ｄｔのものと同数であっても良いし、異なっていても良く、また、４００以上であっても良い。また、正常電流波形Ｗｂからの電流の瞬時値のサンプリング数についても、対象電流波形Ｗｔのものと同じであっても良いし、異なっていても良い。

例えば、図３Ａ〜図３Ｂに示すように、上記の複数の実効値Ｉｅは、回転機械９１に異常が生じている状態である異常状態と、異常が生じていない状態である正常状態とで、挙動が異なる。具体的には、図３Ａ〜図３Ｂに示すように、横軸を、電流の瞬時値の１周期分の波形の数（波形数）、縦軸を実効値Ｉｅとしたグラフに異常状態および正常状態の各々における上記の複数の実効値Ｉｅをプロットすると、異常状態の時（異常時）の実効値Ｉｅの振幅は、正常状態の時（正常時）の実効値Ｉｅの振幅よりも大きく変化する。なお、図３Ａおよび図３Ｂには、回転機械９１に相互に異なる種類の異常が生じた場合が示されている。

このため、診断対象機械９の異常時および正常時の両方における複数の実効値Ｉｅの確率密度分布をそれぞれ算出すると、例えば図４Ａ〜図４Ｂのようになる。図４Ａ〜図４Ｂのグラフは、横軸が実効値Ｉｅであり、縦軸が確率密度である。また、図４Ａが異常時を示し、図４Ｂが正常時を示すが、異常時の方が正常時よりも分布が広がっており、ばらつき度合いが大きくなっている。すなわち、正常時と異常時とで実効値Ｉｅの分布状態に違いが出てくることから、検知部４は、対象分布情報Ｄｔと正常分布情報Ｄｂとの分布状態の違いに基づいて、異常検知を実行する。

なお、上記の第１所定期間および上記の第２所定期間は、互いに重なる部分のない期間であり、これらの期間の長さは通常は同じとするが、相互に異なる長さの期間であっても良い。また、記憶装置１２に正常電流波形Ｗｂ、あるいは正常電流波形Ｗｂを形成する単位電流波形Ｗｕに基づいて算出された実効値Ｉｅが記憶されていても良い。この場合には、診断装置１は、検知部４による異常検知の前に記憶装置１２の記憶情報に基づいて正常分布情報Ｄｂを算出し、記憶部５に記憶しても良い。この記憶部５は、診断装置１が備える記憶装置１２の所定の記憶領域に形成される。

上記の構成によれば、対象分布情報Ｄｔと、診断対象機械９（回転機械９１）の正常時の電流波形Ｗ（正常電流波形Ｗｂ）に基づいて得られる正常分布情報Ｄｂとの比較などに基づいて、診断対象機械９の異常検知を実行する。これによって、診断対象機械９（回転機械９１）が正常であるか異常であるかを適切に判定することができる。

図５は、本発明の一実施形態に係る複数の診断対象機械９の正常電流波形Ｗｂを対比して示す図である。また、図６は、本発明の一実施形態に係る診断の流れを説明するための図である。

幾つかの実施形態では、上記の正常分布情報Ｄｂは、診断対象機械９の個体毎に用意しても良い。本発明者らは、図５に示すように、同じ種類の診断対象機械９であったとしても、個体間で、正常電流波形Ｗｂの振幅の大きさの程度が同じとはならず、実効値Ｉｅの変化の程度が同じとはならない場合があることを見出している。

図５のグラフは、横軸が波形数、縦軸が実効値Ｉｅである。そして、このグラフに、相互に異なる個体である複数（図５では４台）の診断対象機械９（９ａ〜９ｄ）の正常電流波形Ｗｂをそれぞれプロットすると、図５に示すように実効値Ｉｅの振幅の変化が相互に異なっている。実効値Ｉｅの振幅の変化の程度などが異なれば、正常分布情報Ｄｂも相互に異なり得る。よって、診断対象機械９の個体毎の正常分布情報Ｄｂを予め測定などを通して取得（保存）し、対象分布情報Ｄｔと、この対象分布情報Ｄｔを取得した診断対象機械９の正常分布情報Ｄｂとの比較などに基づいて異常検知を行うことで、誤検知を防止することが可能となる。

図２に示す実施形態では、上述した実効値取得部２には、診断対象機械９に接続された電流測定装置７から測定電流の瞬時値または実効値Ｉｅである測定値Ｉｍが入力されるように構成されている。そして、図６に示すように、実効値取得部２は、回転機械９１の回転開始から回転停止までの期間Ｔ内に設定される第１稼働期間Ｔａに入力される測定値Ｉｍに基づいて、正常電流波形Ｗｂに関する複数の実効値Ｉｅを取得する第１取得部２１と、この第１稼働期間Ｔａの経過後で、かつ、回転機械９１の回転停止までの間に設定される第２稼働期間Ｔｂに入力される測定値Ｉｍに基づいて、対象電流波形Ｗｔに関する複数の実効値Ｉｅを取得する第２取得部２２と、を有している。よって、回転機械９１は、第１稼働期間Ｔａと第２稼働期間Ｔｂとの間は継続して回転している状態にある。

図６に示す実施形態では、回転機械９１の回転開始後、回転状態が安定したとみなせる時期以降において、回転機械９１の正常性を確認した後、第１稼働期間Ｔａを設定している。その後、定期的に第２稼働期間Ｔｂを設定し、診断対象機械９の回転機械９１の正常性を監視している。また、回転の停止後に再度回転を開始した後も、同様に、第１稼働期間Ｔａの後に第２稼働期間Ｔｂを設定している。

上記の構成によれば、回転機械９１が停止することなく回転を継続している状態において、正常分布情報Ｄｂの算出の元になる測定値Ｉｍ（測定電流の瞬時値またはその実効値Ｉｅ）が取得された後に、対象分布情報Ｄｔの算出の元になる測定値Ｉｍが取得される。換言すれば、対象電流波形Ｗｔを取得する前における回転機械９１の回転時の診断対象機械９の状態を正常と定義しつつ、その後において回転機械が回転を継続している状態で得られた測定値Ｉｍを診断のために監視する。

本発明者らは、回転機械９１の異常の種類によっては、異常発生時において電流の瞬時値の分布状態に基づいて異常検知を実行した場合（特許文献２参照）よりも、電流の実効値Ｉｅの分布状態に基づいて異常検知を実行した場合の方が、異常判定の精度が高いことを見出している。これによって、正常時の診断対象機械９の個体間に、電流波形Ｗに基づいて得られる複数の実効値Ｉｅのばらつき度合いに違いがあるような場合であっても、診断対象機械毎に正常状態を適切に定義することが可能となる。また、正常状態に対する変化を回転機械９１の継続した回転時に定期的などで監視することで、その傾向に基づいて異常の予兆診断を行うこともできる。

また、上述した実施形態において、幾つかの実施形態では、上述した対象分布情報Ｄｔおよび正常分布情報Ｄｂは、確率密度関数などの確率分布であり、検知部４は、各々が確立分布である対象分布情報Ｄｔと正常分布情報Ｄｂとの距離に基づいて、異常検知を実行しても良い。この距離は、例えば周知なカルバック・ライブラー距離や相対ピアソン距離など、２つの確率分布（確率密度関数）の違いを定量化することが可能な指標値であれば良い。なお、正常分布情報Ｄｂである確率分布をｐ（ｘ）、対象分布情報Ｄｔである確率分布をｐ’（ｘ）とすると、相対ピアソン距離は、例えば、∫ｑ_α（ｘ）［｛ｐ（ｘ）／ｑ_α（ｘ）｝−１］^２ｄｘで算出可能であり、この際、ｑ_α＝αｐ＋（１−α）ｐ’ ０≦α＜１の関係を有する。

図２に示す実施形態では、検知部４は、対象分布情報Ｄｔである確率分布（例えば確率密度関数など）と、正常分布情報Ｄｂである確率分布との間の相対ピアソン距離に基づいて、診断対象機械９の異常検知を行うようになっている。実効値Ｉｅの確率分布が正規分布に近い場合には、確率分布の端部においては確率密度がゼロとなるが、このような場合であっても相対ピアソン距離を用いることによって、ノイズに対してロバストな異常検知を行うことが可能となる。

以下、上述した診断装置１が行う処理に対応する診断方法について、図７を用いて説明する。図７は、本発明の一実施形態に係る診断方法を示す図である。

この診断方法は、回転機械９１を有する診断対象機械９を、その回転機械９１の回転時の測定電流に基づいて診断する方法である。図６に示すように、診断方法は、実効値取得ステップ（Ｓ１）と、分布情報算出ステップ（Ｓ２）と、検知ステップ（Ｓ３）と、を備える。これらの各ステップについて説明する。
これらのステップについて、図１に示す診断系６において実行する場合を例に説明する。

実効値取得ステップ（Ｓ１）は、対象電流波形Ｗｔにおける規定周期数毎の実効値Ｉｅを取得するステップである。実効値取得ステップ（Ｓ１）は、既に説明した実効値取得部２が実行する処理内容と同様であるため、詳細は省略する。

分布情報算出ステップ（Ｓ２）は、上記の実効値取得ステップによって取得された複数の実効値Ｉｅの対象分布情報Ｄｔを算出するステップである。分布情報算出ステップ（Ｓ２）は、既に説明した分布情報算出部３が実行する処理内容と同様であるため、詳細は省略する。

検知ステップ（Ｓ３）は、上記の分布情報算出ステップによって算出された対象分布情報Ｄｔに基づいて、診断対象機械９の異常検知を実行するステップである。この検知ステップ（Ｓ３）は、既に説明した検知部４が実行する処理内容と同様であるため、詳細は省略する。

図７に示す実施形態では、ステップＳ０において、上述した正常分布情報Ｄｂを取得する。具体的には、診断対象機械９が備える回転機械９１の回転開始から回転停止までの期間Ｔ内に設定される第１稼働期間Ｔａにおいて入力される測定値Ｉｍに基づいて、正常電流波形Ｗｂに関する複数の実効値Ｉｅを取得すると共に、取得された複数の実効値Ｉｅに基づいて正常分布情報Ｄｂを算出し、記憶部５に記憶している。なお、測定値Ｉｍが電流の実効値Ｉｅである場合には、測定値Ｉｍが入力されることで複数の実効値Ｉｅが取得されるが、測定値Ｉｍが測定電流の瞬時値である場合には、既に説明したように第２所定期間の間の測定電流の瞬時値に基づいて複数の実効値Ｉｅを算出する。こうして取得された複数の実効値Ｉｅに基づいて、正常分布情報Ｄｂをする。

その後、ステップＳ１において、上記の実効値取得ステップを実行し、対象電流波形Ｗｔから得られる複数の実効値Ｉｅを取得する。ステップＳ２において、分布情報算出ステップを実行し、対象分布情報Ｄｔを算出する。そして、ステップＳ３において、検知ステップを実行し、診断対象機械９の異常検知を実行する。図６に示す実施形態では、対象分布情報Ｄｔと正常分布情報Ｄｂとの相対ピアソン距離を算出し、この相対ピアソン距離に基づいて異常検知を実行している。また、ステップＳ４において、異常検知の結果を例えばディスプレイなどの出力装置１４に出力する。

本発明は上述した実施形態に限定されることはなく、上述した実施形態に変形を加えた形態や、これらの形態を適宜組み合わせた形態も含む。
＜付記＞

（１）本発明の少なくとも一実施形態に係る診断装置（１）は、
回転機械（９１）を有する診断対象機械（９）を前記回転機械（９１）の回転時の測定電流に基づいて診断する診断装置（１）であって、
前記測定電流の第１所定期間での時間推移である対象電流波形（Ｗｔ）における規定周期数毎の実効値（Ｉｅ）を取得するよう構成された実効値取得部（２）と、
取得された複数の前記実効値（Ｉｅ）の分布状態を表す対象分布情報（Ｄｔ）を算出するよう構成された分布情報算出部（３）と、
算出された前記対象分布情報（Ｄｔ）に基づいて、前記診断対象機械（９）の異常検知を実行するよう構成された検知部（４）と、を備える。

上記（１）の構成によれば、診断対象機械（９）が備える例えばタービンや発電機、モータなどの回転機械（９１）の回転時に、回転機械（９１）（電動機）に供給される交流電流や、回転機械（９１）（発電機）から出力される交流電流などの電流の瞬時値（測定電流）の時間推移である電流波形（Ｗｔ）における規定周期数毎の電流波形（後述する単位電流波形）の実効値（Ｉｅ）を取得する。そして、こうして取得された複数の実効値（Ｉｅ）の、例えば確率密度関数や標準偏差などの分布状態を表す情報（対象分布情報（Ｄｔ））に基づいて、診断対象機械（９）の異常検知を実行する。これによって、上記の対象電流波形（Ｗｔ）に基づいて得られる複数の電流の実効値（Ｉｅ）の分布状態（ばらつき）に基づいて回転機械（９１）の診断を実行することで、回転機械（９１）が正常であるか、異常であるかをより適切に診断することができる。

（２）幾つかの実施形態では、上記（１）の構成において、
前記診断対象機械（９）が正常状態にある際の前記測定電流の第２所定期間での時間推移である正常電流波形（Ｗｂ）における前記規定周期数毎の前記実効値（Ｉｅ）の分布状態を表す正常分布情報（Ｄｂ）を記憶するよう構成された記憶部（５）を、さらに備え、
前記検知部（４）は、前記対象分布情報（Ｄｔ）と前記正常分布情報（Ｄｂ）とに基づいて、前記異常検知を実行する。

上記（２）の構成によれば、対象分布情報（Ｄｔ）と、診断対象機械（９）（回転機械（９１））の正常時の電流波形（Ｗｂ）に基づいて得られる正常分布情報（Ｄｂ）との比較などに基づいて、診断対象機械（９）の異常検知を実行する。これによって、個体診断対象機械（９）（回転機械（９１））が正常であるか異常であるかを適切に判定することができる。

（３）幾つかの実施形態では、上記（２）の構成において、
前記対象分布情報（Ｄｔ）および前記正常分布情報（Ｄｂ）は確率分布であり、
前記検知部（４）は、前記対象分布情報（Ｄｔ）と前記正常分布情報（Ｄｂ）との距離に基づいて、前記異常検知を実行する。

上記（３）の構成によれば、対象分布情報（Ｄｔ）である確率分布（例えば確率密度関数など）と、正常分布情報（Ｄｂ）である確率分布との距離に基づいて、診断対象機械（９）の異常検知を実行する。これによって、上記の距離と閾値との比較に基づいて正常、異常を判定するなど、診断対象機械（９）の異常検知を容易に行うことができる。

（４）幾つかの実施形態では、上記（３）の構成において、
前記距離は相対ピアソン距離である。

上記（４）の構成によれば、検知部（４）は、対象分布情報（Ｄｔ）である確率分布と、正常分布情報（Ｄｂ）である確率分布との間の相対ピアソン距離に基づいて、診断対象機械（９）の異常検知を行う。実効値（Ｉｅ）の確率分布が正規分布に近い場合には、確率分布の端部においては確率密度がゼロとなるが、このような場合であっても相対ピアソン距離を用いることによって、ノイズに対してロバストな異常検知を行うことができる。

（５）幾つかの実施形態では、上記（２）〜（４）の構成において、
前記実効値取得部（２）には、前記診断対象機械（９）に接続された電流測定装置（７）から前記測定電流の瞬時値または前記実効値（Ｉｅ）である測定値（Ｉｍ）が入力されるように構成されており、
前記実効値取得部（２）は、
前記回転機械（９１）の回転開始から回転停止までの期間内に設定される第１稼働期間（Ｔａ）に入力される前記測定値（Ｉｍ）に基づいて、前記正常電流波形（Ｗｂ）に関する前記複数の実効値（Ｉｅ）を取得する第１取得部（２１）と、
前記第１稼働期間（Ｔａ）の経過後で、かつ、前記回転停止までの間に設定される第２稼働期間（Ｔｂ）に入力される前記測定値（Ｉｍ）に基づいて、前記対象電流波形（Ｗｔ）に関する前記複数の実効値（Ｉｅ）を取得する第２取得部（２２）と、を有する。

上記（５）の構成によれば、回転機械（９１）が停止することなく回転を継続している状態において、正常分布情報（Ｄｂ）の算出の元になる測定値（Ｉｍ）（測定電流の瞬時値または実効値（Ｉｅ））が取得された後に、対象分布情報（Ｄｔ）の算出の元になる測定値（Ｉｍ）が取得される。換言すれば、対象電流波形（Ｗｔ）を取得する前における回転機械（９１）の回転時の診断対象機械（９）の状態を正常と定義しつつ、その後において回転機械（９１）が回転を継続している状態で得られた測定値（Ｉｍ）を診断のために監視する。これによって、正常時の診断対象機械（９）の個体間に、電流波形に基づいて得られる複数の実効値（Ｉｅ）のばらつき度合いに違いがあるような場合であっても、診断対象機械（９）毎に正常状態を適切に定義することができる。また、正常状態に対する変化を回転機械（９１）の継続した回転時に定期的などで監視することで、その傾向に基づいて異常の予兆診断を行うこともできる。

（６）幾つかの実施形態では、上記（５）の構成において、
前記回転機械（９１）は、前記第１稼働期間（Ｔａ）と前記第２稼働期間（Ｔｂ）との間は継続して回転している状態にある。

上記（６）の構成によれば、上記（５）と同様の効果を奏する。なお、同一の回転機械（９１）である場合など、回転機械（９１）が変わる等の条件変更がなければ、回転機械（９１）は、第１稼働期間（Ｔａ）と第２稼働期間（Ｔｂ）との間は継続して回転している状態にある必要はない。

（７）幾つかの実施形態では、上記（１）〜（６）の構成において、
前記実効値（Ｉｅ）は、前記規定周期数の電流波形からのサンプリングにより得られた複数の前記測定電流の瞬時値に基づいて算出され、
前記サンプリングによるサンプリング数は、前記規定周期数を構成する１周期あたり９００以上１２００以下である。

上記（７）の構成によれば、サンプリング数を９００以上１２００以下とすることにより、回転機械（９１）の稼働状態（正常または異常）が適切に反映された対象分布情報（Ｄｔ）を得ることができ、適切な診断精度での診断を行うことができる。

（８）幾つかの実施形態では、上記（１）〜（７）の構成において、
前記規定周期数は１である。
上記（８）の構成によれば、規定周期数を１周期とすることにより、回転機械（９１）の稼働状態（正常または異常）が適切に反映された対象分布情報（Ｄｔ）を得ることができ、適切な診断精度での診断を行うことができる。

（９）幾つかの実施形態では、上記（１）〜（８）の構成において、
前記分布情報算出部（３）は、前記実効値取得部（２）によって取得された４００以上の前記実効値（Ｉｅ）に基づいて、前記対象分布情報（Ｄｔ）を算出する。

上記（９）の構成によれば、４００以上の実効値（Ｉｅ）に基づいて対象分布情報（Ｄｔ）を算出することで、回転機械（９１）の稼働状態（正常または異常）が対象分布情報（Ｄｔ）に適切に反映されるようにすることができる。

（１０）本発明の少なくとも一実施形態に係る診断方法は、
回転機械（９１）を有する診断対象機械（９）を前記回転機械（９１）の回転時の測定電流に基づいて診断する診断方法であって、
前記測定電流の第１所定期間での時間推移である対象電流波形（Ｗｔ）における規定周期数毎の実効値（Ｉｅ）を取得するステップ（例えば図６のＳ１）と、
取得された複数の前記実効値（Ｉｅ）の分布状態を表す対象分布情報（Ｄｔ）を算出するステップ（例えば図６のＳ２）と、
算出された前記対象分布情報（Ｄｔ）に基づいて、前記診断対象機械（９）の異常検知を実行するステップ（例えば図６のＳ３）と、を備える。
上記（１０）の構成によれば、上記（１）と同様の効果を奏する。

（１１）本発明の少なくとも一実施形態に係る診断プログラム（１０）は、
回転機械（９１）を有する診断対象機械（９）を前記回転機械（９１）の回転時の測定電流に基づいて診断する診断プログラム（１０）であって、
コンピュータに、
前記測定電流の第１所定期間での時間推移である対象電流波形（Ｗｔ）における規定周期数毎の実効値（Ｉｅ）を取得するよう構成された実効値取得部（２）と、
取得された複数の前記実効値（Ｉｅ）の分布状態を表す対象分布情報（Ｄｔ）を算出するよう構成された分布情報算出部（３）と、
算出された前記対象分布情報（Ｄｔ）に基づいて、前記診断対象機械（９）の異常検知を実行するよう構成された検知部（４）と、を実現させるためのプログラムである。
上記（１１）の構成によれば、上記（１）と同様の効果を奏する。

１ 診断装置
１０ 診断プログラム
１１ プロセッサ
１２ 記憶装置
１４ 出力装置
２ 実効値取得部
２１ 第１取得部
２２ 第２取得部
３ 分布情報算出部
４ 検知部
５ 記憶部
６ 診断システム
７ 電流測定装置
７１ 通信媒体
８ 電気盤
９ 診断対象機械
９１ 回転機械
Ｉ 電流
Ｉｅ 実効値
Ｉｍ 測定値
Ｗ 電流波形
Ｗｂ 正常電流波形
Ｗｔ 対象電流波形
Ｗｕ 単位電流波形
Ｄｂ 正常分布情報
Ｄｔ 対象分布情報
Ｔ 回転機械の回転開始から回転停止までの期間
Ｔａ 第１稼働期間
Ｔｂ 第２稼働期間

Claims (11)

1. 回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断装置であって、
   前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するよう構成された実効値取得部と、
   取得された複数の前記実効値の分布状態を表す対象分布情報を算出するよう構成された分布情報算出部と、
   算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するよう構成された検知部と、を備える診断装置。
2. 前記診断対象機械が正常状態にある際の前記測定電流の第２所定期間での時間推移である正常電流波形における前記規定周期数毎の前記実効値の分布状態を表す正常分布情報を記憶するよう構成された記憶部を、さらに備え、
   前記検知部は、前記対象分布情報と前記正常分布情報とに基づいて、前記異常検知を実行する請求項１に記載の診断装置。
3. 前記対象分布情報および前記正常分布情報は確率分布であり、
   前記検知部は、前記対象分布情報と前記正常分布情報との距離に基づいて、前記異常検知を実行する請求項２に記載の診断装置。
4. 前記距離は相対ピアソン距離である請求項３に記載の診断装置。
5. 前記実効値取得部には、前記診断対象機械に接続された電流測定装置から前記測定電流の瞬時値または前記実効値である測定値が入力されるように構成されており、
   前記実効値取得部は、
   前記回転機械の回転開始から回転停止までの期間内に設定される第１稼働期間に入力される前記測定値に基づいて、前記正常電流波形に関する前記複数の実効値を取得する第１取得部と、
   前記第１稼働期間の経過後で、かつ、前記回転停止までの間に設定される第２稼働期間に入力される前記測定値に基づいて、前記対象電流波形に関する前記複数の実効値を取得する第２取得部と、を有する請求項２〜４のいずれか１項に記載の診断装置。
6. 前記回転機械は、前記第１稼働期間と前記第２稼働期間との間は継続して回転している状態にある請求項５に記載の診断装置。
7. 前記実効値は、前記規定周期数の電流波形からのサンプリングにより得られた複数の前記測定電流の瞬時値に基づいて算出され、
   前記サンプリングによるサンプリング数は、前記規定周期数を構成する１周期あたり９００以上１２００以下である請求項１〜６のいずれか１項に記載の診断装置。
8. 前記規定周期数は１である請求項１〜７のいずれか１項に記載の診断装置。
9. 前記分布情報算出部は、前記実効値取得部によって取得された４００以上の前記実効値に基づいて、前記対象分布情報を算出することを特徴とする請求項１〜８のいずれか１項に記載の診断装置。
10. 回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断方法であって、
    前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するステップと、
    取得された複数の前記実効値の分布状態を表す対象分布情報を算出するステップと、
    算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するステップと、を備える診断方法。
11. 回転機械を有する診断対象機械を前記回転機械の回転時の測定電流に基づいて診断する診断プログラムであって、
    コンピュータに、
    前記測定電流の第１所定期間での時間推移である対象電流波形における規定周期数毎の実効値を取得するよう構成された実効値取得部と、
    取得された複数の前記実効値の分布状態を表す対象分布情報を算出するよう構成された分布情報算出部と、
    算出された前記対象分布情報に基づいて、前記診断対象機械の異常検知を実行するよう構成された検知部と、を実現させるための診断プログラム。