(57)【特許請求の範囲】
【請求項１】
風力発電設備の健全性を評価するためのシステムであって、
前記風力発電設備の運転状態と相関がある複数の特性項目の指標である特性値の組合せからなる複数のデータセットを単位期間ごとに取得するように構成されたデータセット取得部と、
前記単位期間内の前記複数のデータセットの少なくとも一部を、前記風力発電設備の出力を示すパラメータのクラス毎に分類するためのクラス分類部と、
前記クラス分類部によって各クラスに分類された前記データセットを平均化し、平均データセットをクラス毎に算出する平均処理部と、
前記風力発電設備の正常時における前記特性値の組合せからなる基準データセットと前記平均データセットとの偏差を示す偏差指標値をクラス毎に算出するための偏差指標値算出部と、
少なくとも一つのクラスにおいて前記偏差指標値が閾値を上回った場合、前記単位期間において前記風力発電設備の異常が生じたと判定するように構成された異常判定部と、を備えることを特徴とする風力発電設備の診断システム。
【請求項２】
前記基準データセットは、前記パラメータのクラス毎に分類され、前記パラメータのクラス毎に平均化されており、
前記偏差指標値算出部は、同一クラスに属する前記平均データセットと前記基準データセットとの偏差を示す前記偏差指標値をクラス毎に算出するように構成されたことを特徴とする請求項１に記載の風力発電設備の診断システム。
【請求項３】
条件（Ａ）又は条件（Ｂ）の少なくとも一方からなる第１判定条件を満たさないデータセットを含む非対象単位期間に属する前記複数のデータセットを排除し、前記第１判定条件を全てのデータセットが満たす対象単位期間に属する対象データセットを選別するように構成された対象データセット選別部をさらに備え、
前記クラス分類部は、前記複数のデータセットのうち、前記対象データセットについてのみ前記パラメータのクラス毎に分類するように構成されており、
前記条件（Ａ）は、前記風力発電設備のロータ回転数ωｒが、カットイン風速時におけるロータ回転数ωｒ０以上かつ定格回転数ωｒ＿ｒａｔｅｄ未満の閾値ωｒ＿ｔｈ以上であるという条件であり、
前記条件（Ｂ）は、前記風力発電設備の発電機回転数ωｇが、電力系統への発電機の併入時における発電機回転数ωｇ０以下の閾値ωｇ＿ｔｈ以上であるという条件であることを特徴とする請求項１又は２に記載の風力発電設備の診断システム。
【請求項４】
前記基準データセットは、前記条件（Ａ）又は前記条件（Ｂ）の前記少なくとも一方を満たすことを特徴とする請求項３に記載の風力発電設備の診断システム。
【請求項５】
前記対象データセット選別部は、前記第１判定条件に加えて、条件（Ｃ）又は条件（Ｄ）の少なくとも一方からなる第２判定条件を全てのデータセットが満たす前記対象単位期間におけるデータセットを前記対象データセットとして選別するように構成され、
前記条件（Ｃ）は、前記風力発電設備に一時的な出力制限が課せられていないという条件であり、
前記条件（Ｄ）は、前記風力発電設備の運転状態に応じた温度変化を示す前記風力発電設備の部位の温度変動率が閾値以下であるという条件であることを特徴とする請求項３又は４に記載の風力発電設備の診断システム。
【請求項６】
前記偏差指標値算出部は、前記基準データセットにより構成される単位空間に対する、前記対象データセットにより構成される信号空間のマハラノビス距離を前記偏差指標値として算出するように構成され、
前記異常判定部は、前記マハラノビス距離に基づいて、前記対象単位期間における前記風力発電設備の異常の有無を判定するように構成されたことを特徴とする請求項１乃至５の何れか一項に記載の風力発電設備の診断システム。
【請求項７】
前記対象データセットを構成する前記特性値のうち、１以上の着目特性項目以外の特性項目の特性値の組み合わせからなる分析用対象データセットと、前記基準データセットを構成する前記特性値のうち、前記着目特性項目以外の特性項目の特性値の組み合わせからなる分析用基準データセットと、の偏差を示す分析用偏差指標値を算出し、前記偏差指標値に対して該分析用偏差指標値が占める割合に基づき、各々の前記特性項目の前記偏差指標値への寄与度を算出するように構成された寄与度算出部と、
前記異常判定部により前記風力発電設備に異常があると判定された場合に、前記寄与度算出部により算出された前記寄与度に基づいて、前記複数の特性項目のうち、前記異常の要因となっている特性項目を特定するように構成された異常要因特定部と、をさらに備えることを特徴とする請求項１乃至６の何れか一項に記載の風力発電設備の診断システム。
【請求項８】
前記偏差指標値算出部によって算出された前記マハラノビス距離と、前記寄与度算出部によって算出された前記特性項目の各々の前記寄与度とに基づいて、前記特性項目の各々のマハラノビス距離成分を算出するＭＤ成分算出部と、
前記ＭＤ成分算出部によって算出された前記特性項目の各々の前記マハラノビス距離成分の時系列変化を表示するための表示部と、をさらに備えることを特徴とする請求項７に記載の風力発電設備の診断システム。
【請求項９】
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、前記ロータの回転を前記発電機に伝えるためのドライブトレインと、をさらに備え、
前記特性項目は、風速、前記発電機の送電端出力、前記ロータの回転数、及び前記複数の風車翼のピッチ角を含み、
前記異常判定部は、前記ロータ、前記ドライブトレイン、または、前記発電機の少なくとも一つの異常を判定するように構成されたことを特徴とする請求項１乃至８の何れか一項に記載の風力発電設備の診断システム。
【請求項１０】
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、前記風車翼のピッチ角を制御するためのアクチュエータと、をさらに備え、
前記特性項目は、前記アクチュエータの温度、前記発電機の送電端出力及び外気温を含み、
前記異常判定部は、前記アクチュエータの異常を判定するように構成されたことを特徴とする請求項１乃至９の何れか一項に記載の風力発電設備の診断システム。
【請求項１１】
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、前記ロータに接続され、前記ロータの回転を前記発電機に伝えるように構成された主軸と、前記主軸を回転可能に支持するように構成された主軸受と、をさらに備え、
前記特性項目は、前記主軸受の温度、前記主軸受の振動情報および外気温を含み、
前記振動情報は、
前記主軸受の加速度信号を周波数解析して得られる第１信号のうち、前記主軸の回転数の整数倍に対応する周波数の信号強度として定義される前記主軸受の加速度レベル、
前記主軸受の加速度信号のピークを結ぶ包絡線信号を周波数解析して得られる第２信号のうち、前記主軸の回転数の整数倍に対応する周波数の信号強度として定義される前記主軸受の包絡線レベル、
前記主軸受の半径方向速度信号を周波数解析して得られる第３信号のオーバーオール値、または、
前記主軸受の半径方向加速度信号を周波数解析して得られる第４信号のオーバーオール値、
の少なくとも一つを含み、
前記異常判定部は、前記主軸受の異常を判定するように構成されたことを特徴とする請求項１乃至１０の何れか一項に記載の風力発電設備の診断システム。
【請求項１２】
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された同期発電機と、前記同期発電機の界磁電流を調節して前記同期発電機の端子電圧を制御するための自動電圧調節器と、をさらに備え、
前記特性項目は、前記同期発電機の巻線温度、前記自動電圧調節器の電圧、前記自動電圧調節器の電流、前記同期発電機の界磁電圧、および、外気温を含み、
前記異常判定部は、前記同期発電機の異常を判定するように構成されたことを特徴とする請求項１乃至１１の何れか一項に記載の風力発電設備の診断システム。
【請求項１３】
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、をさらに備え、
前記特性項目は、前記発電機の回転軸を軸支する発電機軸受の温度、前記発電機軸受の振動情報および外気温を含み、
前記振動情報は、
前記発電機軸受の加速度信号を周波数解析して得られる第５信号のうち、前記発電機の回転軸の回転数の整数倍に対応する周波数の信号強度として定義される前記発電機軸受の加速度レベル、
前記発電機軸受の加速度信号のピークを結ぶ包絡線信号を周波数解析して得られる第６信号のうち、前記発電機の回転軸の回転数の整数倍に対応する周波数の信号強度として定義される前記発電機軸受の包絡線レベル、
前記発電機軸受の半径方向速度信号を周波数解析して得られる第７信号のオーバーオール値、または、
前記発電機軸受の半径方向加速度信号を周波数解析して得られる第８信号のオーバーオール値、
の少なくとも一つを含み、
前記異常判定部は、前記発電機軸受の異常を判定するように構成されたことを特徴とする請求項１乃至１２の何れか一項に記載の風力発電設備の診断システム。
【請求項１４】
前記データセット取得部は、前記単位期間が１分〜６０分の範囲で、好ましくは１０分程度でデータセットを取得するように構成され、
前記偏差指標値算出部は、前記データセット取得部において前記単位期間にわたって取得された前記複数のデータセットを用いて、３０分〜６００分の範囲で好ましくは６０分の期間で前記偏差指標値を含む分析値を算出するように構成されたことを特徴とする請求項１乃至１３の何れか一項に記載の風力発電設備の診断システム。
【請求項１５】
風力発電設備の健全性を評価するための方法であって、
前記風力発電設備の運転状態と相関がある複数の特性項目の指標である特性値の組合せからなる複数のデータセットを単位期間ごとに複数取得するデータセット取得ステップと、
前記単位期間内の前記複数のデータセットの少なくとも一部を、前記風力発電設備の出力を示すパラメータのクラス毎に分類するクラス分類ステップと、
前記クラス分類ステップで各クラスに分類された前記データセットを平均化し、平均データセットをクラス毎に算出する平均処理ステップと、
前記風力発電設備の正常時における前記特性値の組合せからなる基準データセットと前記平均データセットとの偏差を示す偏差指標値をクラス毎に算出する偏差指標値算出ステップと、
少なくとも一つのクラスにおいて前記偏差指標値が閾値を上回った場合、前記単位期間において前記風力発電設備の異常が生じたと判定する異常判定ステップと、を備えることを特徴とする風力発電設備の診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は風力発電設備の診断システム及び診断方法に関する。
【背景技術】
【０００２】
風力発電装置の状態監視において、多次元情報の処理を行う品質工学手法を用いることが知られている。
例えば、特許文献１には、マハラノビス・タグチメソッド（以下「ＭＴ法」という。）を用いて診断対象の健全性を評価する方法及び該方法を風車に適用した風車の健全性診断装置が開示されている。この健全性評価方法では、計測時間に関連付けられた複数の特性値が特性項目別に格納された被診断データファイルから被診断データセットを抽出するとともに、特定の特性項目の特性値が予め定義されている所定の基準範囲（正常範囲）に属している比較用データファイルから基準データセットを抽出する。そして、被診断データセットと基準データセットとを用いて、風車の状態を表す状態指標値であるマハラノビス距離を算出し、該マハラノビス距離に基づいて診断対象である風車の異常を判定するようになっている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第５１０１３９６号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１の方法では、各特性項目について、単位期間毎（例えば１分間隔）の計測データ（特性値）からなる被診断データセット及び基準データセットを用いてマハラノビス距離を算出する。ここで、各計測データのサンプリング周期が単位期間よりも十分短い場合には、該単位期間に取得された全ての計測データの平均値を各特性項目の特性値としてマハラノビス距離が算出される。
このような、ＭＴ法をはじめとした多次元情報を扱う品質工学手法を用いた風車の健全性診断において、より精度良く異常判定を行うことが望まれる。
【０００５】
上述の事情に鑑みて、本発明の少なくとも一実施形態は、精度の良好な異常判定が可能な風力発電設備の診断システムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
（１）本発明の少なくとも一実施形態に係る風力発電設備の診断システムは、
風力発電設備の健全性を評価するためのシステムであって、
前記風力発電設備の運転状態と相関がある複数の特性項目の指標である特性値の組合せからなる複数のデータセットを単位期間ごとに取得するように構成されたデータセット取得部と、
前記単位期間内の前記複数のデータセットの少なくとも一部を、前記風力発電設備の出力を示すパラメータのクラス毎に分類するためのクラス分類部と、
前記クラス分類部によって各クラスに分類された前記データセットを平均化し、平均データセットをクラス毎に算出する平均処理部と、
前記風力発電設備の正常時における前記特性値の組合せからなる基準データセットと前記平均データセットとの偏差を示す偏差指標値をクラス毎に算出するための偏差指標値算出部と、
少なくとも一つのクラスにおいて前記偏差指標値が閾値を上回った場合、前記単位期間において前記風力発電設備の異常が生じたと判定するように構成された異常判定部と、を備える。
【０００７】
例えば特許文献１に記載のように、データセットの単位期間内における時間平均と、基準データセットとの偏差を示す偏差指標値に基づいて風力発電設備の異常判定を行う場合、データセットに含まれる特性値（例えば出力）に大きなばらつきがあったとしても、特性値のばらつきは平均化により均される。このように平均化されたデータセットから算出された偏差指標値に基づく風力発電設備の異常判定は、単位期間内における比較的大きい又は小さい特性値の影響が低減されたものとなる。
この点、上記（１）の構成では、風力発電設備の出力を示すパラメータ（例えば出力）のクラス毎にデータセットを平均化して、クラス毎に偏差指標値を算出する。このため、単位期間内においてパラメータのばらつきが大きい場合であってもばらつきを丸めることなく、クラス毎に算出された偏差指標値をそれぞれ閾値と比較することで、精度の良好な異常判定が可能となる。
【０００８】
（２）幾つかの実施形態では、上記（１）の構成において、
前記基準データセットは、前記パラメータのクラス毎に分類され、前記パラメータのクラス毎に平均化されており、
前記偏差指標値算出部は、同一クラスに属する前記平均データセットと前記基準データセットとの偏差を示す前記偏差指標値をクラス毎に算出するように構成される。
上記（２）の構成では、クラス毎に平均化された平均データセットと、クラス毎に平均化された基準データセットのうち同一クラスに属する平均データセットと基準データセットを用いて偏差指標値を算出するので、同じ状況下で取得されたデータ同士を比較することができるため、より的確な異常判定が可能となる。
【０００９】
（３）幾つかの実施形態では、上記（１）又は（２）の構成において、
条件（Ａ）又は条件（Ｂ）の少なくとも一方からなる第１判定条件を満たさないデータセットを含む非対象単位期間に属する前記複数のデータセットを排除し、前記第１判定条件を全てのデータセットが満たす対象単位期間に属する対象データセットを選別するように構成された対象データセット選別部をさらに備え、
前記クラス分類部は、前記複数のデータセットのうち、前記対象データセットについてのみ前記パラメータのクラス毎に分類するように構成されており、
前記条件（Ａ）は、前記風力発電設備のロータ回転数ωｒが、カットイン風速時におけるロータ回転数ωｒ０以上かつ定格回転数ωｒ＿ｒａｔｅｄ未満の閾値ωｒ＿ｔｈ以上であるという条件であり、
前記条件（Ｂ）は、前記風力発電設備の発電機回転数ωｇが、電力系統への発電機の併入時における発電機回転数ωｇ０以下の閾値ωｇ＿ｔｈ以上であるという条件である。
風力発電設備が運転中であるか否かにより、特性項目の指標である特性値が異なる場合がある。例えば、風車が運転中であれば、ロータ回転数が所定値以上であるのに対し、風車が停止中であればロータ回転数は所定値以下である。また、風力発電設備の発電機が運転中であれば、発電機回転数が所定値以上であるのに対し、発電機が停止中であれば、発電機回転数は所定値以下である。このように、風力発電設備の運転状態の違いが偏差指標値に少なからず影響するので、この影響を排除した適切な異常判定が望まれる。この点、上記（３）の構成では、風力発電設備の運転状態が所定の場合のデータセットのみを選別して偏差指標値を算出するので適切に異常判定を行うことができる。
【００１０】
（４）幾つかの実施形態では、上記（３）の構成において、前記基準データセットは、前記条件（Ａ）又は前記条件（Ｂ）の前記少なくとも一方を満たす。
上記（４）の構成では、条件（Ａ）又は条件（Ｂ）を満たす対象データセットと基準データセットを用いて偏差指標値を算出するので、風力発電設備の運転状態が所定の場合のデータセットのみを選別して偏差指標値を算出するので適切に異常判定を行うことができる。
【００１１】
（５）幾つかの実施形態では、上記（３）又は（４）の構成において、
前記対象データセット選択部は、前記第１判定条件に加えて、条件（Ｃ）又は条件（Ｄ）の少なくとも一方からなる第２判定条件を全てのデータセットが満たす前記対象単位期間におけるデータセットを前記対象データセットとして選別するように構成され、
前記条件（Ｃ）は、前記風力発電設備に一時的な出力制限が課せられていないという条件であり、
前記条件（Ｄ）は、前記風力発電設備の運転状態に応じた温度変化を示す前記風力発電設備の部位の温度変動率が閾値以下であるという条件である。
風力発電設備の運転状態によって、特性項目間（例えば発電機出力と、風車翼のピッチ角と、ロータ回転数との間）の相関関係が異なる場合がある。例えば、風力発電設備に一時的な出力制限が課せられていれば、特性項目間の相関関係が出力制限されていない場合と比べて変化する。また、風力発電設備の運転開始直後には、特性項目の指標である特性値（例えば、風力発電設備を構成する部品の温度）の変動率が通常運転時に比べて小さい。このように、風力発電設備の運転状態の違いが偏差指標値に少なからず影響するので、この影響を排除した適切な異常判定が望まれる。この点、上記（５）の構成では、風力発電設備の運転状態が所定の場合のデータセットのみを選別して偏差指標値を算出するので適切に異常判定を行うことができる。
【００１２】
（６）幾つかの実施形態では、上記（１）〜（５）の何れかの構成において、
前記偏差指標値算出部は、前記基準データセットにより構成される単位空間に対する、前記対象データセットにより構成される信号空間のマハラノビス距離を前記偏差指標値として算出するように構成され、
前記異常判定部は、前記マハラノビス距離に基づいて、前記対象単位期間における前記風力発電設備の異常の有無を判定するように構成される。
上記（６）の構成では、基準データセットおよび対象データセットに基づいてマハラノビス距離（以下において「ＭＤ値」ともいう。）を算出し、マハラノビス距離に基づいて風力発電設備の異常判定を行う。このため、例えば、ユークリッド距離を用いて異常判定を行う場合に比べて、複数の特性項目間の相関関係を考慮した評価が可能である。
【００１３】
（７）幾つかの実施形態では、上記（１）〜（６）の何れかの構成において、
前記対象データセットを構成する前記特性値のうち、１以上の着目特性項目以外の特性項目の特性値の組み合わせからなる分析用対象データセットと、前記基準データセットを構成する前記特性値のうち、前記着目特性項目以外の特性項目の特性値の組み合わせからなる分析用基準データセットと、の偏差を示す分析用偏差指標値を算出し、前記偏差指標値に対して該分析用偏差指標値が占める割合に基づき、各々の前記特性項目の前記偏差指標値への寄与度を算出するように構成された寄与度算出部と、
前記異常判定部により前記風力発電設備に異常があると判定された場合に、前記寄与度算出部により算出された前記寄与度に基づいて、前記複数の特性項目のうち、前記異常の要因となっている特性項目を特定するように構成された異常要因特定部と、をさらに備える。
上記（７）の構成によれば、各特性項目の偏差指標値への寄与度を算出して、該寄与度に基づいて異常の要因となっている特性項目を特定することができ、異常要因に対する適切な対応が可能となる。
【００１４】
（８）幾つかの実施形態では、上記（７）の構成において、
前記偏差指標値算出部によって算出された前記マハラノビス距離と、前記寄与度算出部によって算出された前記特性項目の各々の前記寄与度とに基づいて、前記特性項目の各々のマハラノビス距離成分を算出するＭＤ成分算出部と、
前記ＭＤ成分算出部によって算出された前記特性項目の各々の前記マハラノビス距離成分の時系列変化を表示するための表示部と、をさらに備える。
上記（８）の構成によれば、マハラノビス距離と各特性項目の寄与度とに基づいて算出される各特性項目のマハラノビス距離成分の時系列変化を表示することにより、各特性項目のマハラノビス距離成分及びマハラノビス距離成分の和であるマハラノビス距離を視覚的に認知することができる。このため、風車の異常を視覚的に判断することができる。
【００１５】
（９）幾つかの実施形態では、上記（１）〜（８）の何れかの構成において、
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、前記ロータの回転を前記発電機に伝えるためのドライブトレインと、をさらに備え、
前記特性項目は、風速、前記発電機の送電端出力、前記ロータの回転数、及び前記複数の風車翼のピッチ角を含み、
前記異常判定部は、前記ロータ、前記ドライブトレイン、または、前記発電機の少なくとも一つの異常を判定するように構成される。
上記（９）の構成によれば、風速、発電機の送電端出力、ロータの回転数、及び複数の風車翼のピッチ角を特性項目として算出された偏差指標値に基づいて、風力発電設備のロータ、ドライブトレイン、または、発電機の少なくとも一つの異常を判定することができる。
【００１６】
（１０）幾つかの実施形態では、上記（１）〜（９）の何れかの構成において、
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、前記風車翼のピッチ角を制御するためのアクチュエータと、をさらに備え、
前記特性項目は、前記アクチュエータの温度、前記発電機の送電端出力及び外気温を含み、
前記異常判定部は、前記アクチュエータの異常を判定するように構成される。
上記（１０）の構成によれば、風車翼のピッチ角を制御するためのアクチュエータの温度、発電機の送電端出力及び外気温を特性項目として算出された偏差指標値に基づいて、風力発電設備の風車翼のピッチ角を制御するためのアクチュエータの異常を判定することができる。
【００１７】
（１１）幾つかの実施形態では、上記（１）〜（１０）の何れかの構成において、
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、前記ロータに接続され、前記ロータの回転を前記発電機に伝えるように構成された主軸と、前記主軸を回転可能に支持するように構成された主軸受と、をさらに備え、
前記特性項目は、前記主軸受の温度、前記主軸受の振動情報および外気温を含み、
前記振動情報は、
前記主軸受の加速度信号を周波数解析して得られる第１信号のうち、前記主軸の回転数の整数倍に対応する周波数の信号強度として定義される前記主軸受の加速度レベル、
前記主軸受の加速度信号のピークを結ぶ包絡線信号を周波数解析して得られる第２信号のうち、前記主軸の回転数の整数倍に対応する周波数の信号強度として定義される前記主軸受の包絡線レベル、
前記主軸受の半径方向速度信号を周波数解析して得られる第３信号のオーバーオール値、または、
前記主軸受の半径方向加速度信号を周波数解析して得られる第４信号のオーバーオール値、
の少なくとも一つを含み、
前記異常判定部は、前記主軸受の異常を判定するように構成される。
上記（１１）の構成によれば、主軸受の温度、主軸受の振動情報および外気温を特性項目として算出された偏差指標値に基づいて、風力発電設備の主軸受の異常を判定することができる。
【００１８】
（１２）幾つかの実施形態では、上記（１）〜（１１）の何れかの構成において、
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された同期発電機と、前記同期発電機の界磁電流を調節して前記同期発電機の端子電圧を制御するための自動電圧調節器と、をさらに備え、
前記特性項目は、前記同期発電機の巻線温度、前記自動電圧調節器の電圧、前記自動電圧調節器の電流、前記同期発電機の界磁電圧、および、外気温を含み、
前記異常判定部は、前記同期発電機の異常を判定するように構成される。
上記（１２）の構成によれば、同期発電機の巻線温度、自動電圧調節器の電圧、自動電圧調節器の電流、同期発電機の界磁電圧、および、外気温を特性項目として算出された偏差指標値に基づいて、風力発電設備の同期発電機の異常を判定することができる。
【００１９】
（１３）幾つかの実施形態では、上記（１）〜（１２）の何れかの構成において、
前記風力発電設備は、複数の風車翼を備えるロータと、前記ロータの回転により駆動されるように構成された発電機と、をさらに備え、
前記特性項目は、前記発電機の回転軸を軸支する発電機軸受の温度、前記発電機軸受の振動情報および外気温を含み、
前記振動情報は、
前記発電機軸受の加速度信号を周波数解析して得られる第５信号のうち、前記発電機の回転軸の回転数の整数倍に対応する周波数の信号強度として定義される前記発電機軸受の加速度レベル、
前記発電機軸受の加速度信号のピークを結ぶ包絡線信号を周波数解析して得られる第６信号のうち、前記発電機の回転軸の回転数の整数倍に対応する周波数の信号強度として定義される前記発電機軸受の包絡線レベル、
前記発電機軸受の半径方向速度信号を周波数解析して得られる第７信号のオーバーオール値、または、
前記発電機軸受の半径方向加速度信号を周波数解析して得られる第８信号のオーバーオール値、
の少なくとも一つを含み、
前記異常判定部は、前記発電機軸受の異常を判定するように構成される。
上記（１３）の構成によれば、発電機の回転軸を軸支する発電機軸受の温度、発電機軸受の振動情報および外気温を特性項目として算出された偏差指標値に基づいて、風力発電設備の発電機軸受の異常を判定することができる。
【００２０】
（１４）幾つかの実施形態では、上記（１）〜（１４）の何れかの構成において、
前記データセット取得部は、前記単位期間が１分〜６０分の範囲で、好ましくは１０分程度でデータセットを取得するように構成され、
前記偏差指標値算出部は、前記対象データセット選別部において前記データセット取得単位期間にわたって取得された前記複数のデータセットから選別された対象データセットを用いて、３０分〜６００分の範囲で好ましくは６０分の期間で前記偏差指標値を含む分析値を算出するように構成される。
【００２１】
（１５）本発明の少なくとも一実施形態に係る風力発電設備の診断方法は、
風力発電設備の健全性を評価するための方法であって、
前記風力発電設備の運転状態と相関がある複数の特性項目の指標である特性値の組合せからなる複数のデータセットを単位期間ごとに複数取得するデータセット取得ステップと、
前記単位期間内の前記複数のデータセットの少なくとも一部を、前記風力発電設備の出力を示すパラメータのクラス毎に分類するクラス分類ステップと、
前記クラス分類ステップで各クラスに分類された前記データセットを平均化し、平均データセットをクラス毎に算出する平均処理ステップと、
前記風力発電設備の正常時における前記特性値の組合せからなる基準データセットと前記平均データセットとの偏差を示す偏差指標値をクラス毎に算出する偏差指標値算出ステップと、
少なくとも一つのクラスにおいて前記偏差指標値が閾値を上回った場合、前記単位期間において前記風力発電設備の異常が生じたと判定する異常判定ステップと、を備えることを特徴とする風力発電設備の診断方法。
【００２２】
上記（１５）の方法では、風力発電設備の出力を示すパラメータ（例えば出力）のクラス毎にデータセットを平均化して、クラス毎に偏差指標値を算出する。このため、単位期間内においてパラメータのばらつきが大きい場合であってもばらつきを丸めることなく、クラス毎に算出された偏差指標値をそれぞれ閾値と比較することで、精度の良好な異常判定（健全性診断）が可能となる。
【発明の効果】
【００２３】
本発明の少なくとも一実施形態によれば、精度の良好な健全性の診断が可能な風力発電設備の診断システムが提供される。
【図面の簡単な説明】
【００２４】
【図１】一実施形態に係る風力発電設備の概略構成図である。
【図２】一実施形態に係る風力発電設備の診断システムの構成の概略を示すブロック図である。
【図３】一実施形態に係るデータセットを構成するひと組の特性項目を示す表である。
【図４】一実施形態に係るデータセットのクラスの分類を示す図である。
【図５】クラス分類及び平均化を行った平均化データセットの一例を示す表である。
【図６】マハラノビス距離算出対象の単位空間データ（基準データセットの集合）の要素の一例を示す図である。
【図７】マハラノビス距離算出対象の信号空間データ（平均データセットの集合）の要素の一例を示す図である。
【図８】図５に示されるデータセットを用いてマハラノビス距離を算出した結果を示す図である。
【図９】寄与度の算出方法を説明するための表である。
【図１０】一実施形態に係る表示部により表示される各特性項目のマハラノビス距離成分の時系列変化の一例である。
【発明を実施するための形態】
【００２５】
以下、添付図面を参照して本発明の幾つかの実施形態について説明する。ただし、実施形態として記載されている又は図面に示されている構成部品の寸法、材質、形状、その相対的配置等は、本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
【００２６】
まず、一実施形態に係る風力発電設備の診断システムが診断対象とする風力発電設備について説明する。図１は、一実施形態に係る風力発電設備の概略構成図である。図１に示すように、風力発電設備１は、複数の風車翼２と、風車翼２が取り付けられたハブ４とを備えるロータ３と、ロータ３の回転により駆動されるように構成された発電機１０と、ロータ３の回転を発電機１０に伝えるためのドライブトレイン９と、を備える。風車翼２が風を受けると、風車翼２を含むロータ３が回転する。ロータ３の回転エネルギーは、ドライブトレイン９を介して発電機１０に伝達され、発電機１０によって電気エネルギーに変換されるようになっている。
なお、ドライブトレイン９を含む種々の機器は、水上又は地上に立設されたタワー１４の上に設置されたナセル１２に収納されていてもよい。また、ハブ４は、ハブカバー５によって覆われていてもよい。また、発電機１０は電力系統に連系される同期発電機であってもよい。
【００２７】
ドライブトレイン９は、ハブ４（ロータ３）に接続される主軸６と、主軸６を回転可能に支持する主軸受７（７ａ，７ｂ）と、ロータ３の回転を変速して発電機１０に伝えるためのトランスミッション８と、を含む。トランスミッション８は、油圧ポンプ及び油圧モータを含む油圧トランスミッションであってもよい。あるいは、トランスミッション８は、主軸６の回転を増速して発電機１０に出力するように構成された機械式（ギア式）増速機であってもよい。
なお、一実施形態では、風力発電設備１は、ドライブトレインを介さずに主軸６の回転により発電機１０を駆動するダイレクトドライブ型の風力発電装置であってもよい。
【００２８】
一実施形態では、風力発電設備１は、風車翼２のピッチ角を制御するためのアクチュエータ（不図示）をさらに備える。該アクチュエータは、電動ピッチモータであってもよく、又は、油圧シリンダを用いた油圧アクチュエータであってもよい。
【００２９】
一実施形態では、発電機１０は同期発電機であり、風力発電設備１は、同期発電機（発電機１０）の界磁電流を調節して同期発電機（発電機１０）の端子電圧を制御するための自動電圧調節器（不図示）をさらに備える。
【００３０】
次に、一実施形態に係る風力発電設備の診断システムについて説明する。図２は、一実施形態に係る風力発電設備の診断システムの構成の概略を示すブロック図である。図２に示すように、風力発電設備の診断システム１００は、風力発電設備１の健全性を評価するための偏差指標値の算出に用いるデータセットを生成するためのデータセット生成部２０及び基準データセット生成部３０と、算出された偏差指標値に基づいて風力発電設備１の異常の有無を判定するための健全性評価部４０と、を備える。
【００３１】
データセット生成部２０は、データセット取得部２２と、対象データセット選別部２４と、クラス分類部２６と、平均処理部２８と、を備える。
基準データセット生成部３０は、基準データセット取得部３２と、基準データセット選別部３４と、基準データセットクラス分類部３６と、基準データセット平均処理部３８と、を備える。
健全性評価部４０は、偏差指標値算出部４２と、異常判定部４３と、寄与度算出部４４と、異常要因特定部４５と、ＭＤ成分算出部４６と、表示部４７と、を備える。
【００３２】
以下に、診断システム１００が備える各部について、偏差指標値としてマハラノビス距離を算出する実施形態に付いて、図を参照して説明する。
【００３３】
データセット生成部２０は、風力発電設備１の健全性を評価するためのマハラノビス距離（偏差指標値）の算出において信号空間を構成するデータセットを生成するように構成される。
データセット取得部２２は、風力発電設備１の運転状態と相関がある複数の特性項目の指標である特性値の組合せからなる複数のデータセットを単位期間Ｔｉ（ｉ＝１，２，…）ごとに取得するように構成される。
図３は、一実施形態に係るデータセットを構成するひと組の特性項目を示す表である。同図に示すように、一実施形態では、発電機回転数、発電機出力（発電機１０の送電端出力）、風速、及び、３枚の風車翼２（第１翼〜第３翼）それぞれのピッチ角、の６種類の特性項目により、データセットを構成する。
各特性項目の指標である特性値とは、典型的には各特性項目についての計測値である。
各特性項目によってサンプリング周期が異なる場合、全ての特性項目について単位時間あたりのデータ数を揃えるようにしてもよい。例えば、複数の特性項目のうち最も長いサンプリング周期が１分である場合、サンプリング周期が１分よりも短い他の特性項目については、１分間に計測された計測値の平均を算出し、算出された１分間平均値によりデータセットを構成するようにしてもよい。
【００３４】
複数の特性項目は、診断対象に応じて適宜選択できる。
例えば、ロータ３、ドライブトレイン９、または、発電機１０を異常判定の対象とする場合、風速、発電機１０の送電端出力、ロータ３の回転数、及び複数の風車翼２（第１翼〜第３翼）のピッチ角を特性項目としてもよい。
【００３５】
また、例えば、風車翼２のピッチ角を制御するためのアクチュエータを異常判定の対象とする場合、該アクチュエータの温度、発電機１０の送電端出力及び外気温を特性項目としてもよい。
【００３６】
また、例えば、風力発電設備１の主軸受７（７ａ，７ｂ）を異常判定の対象とする場合、主軸受７（７ａ，７ｂ）の温度、主軸受７（７ａ，７ｂ）の振動情報および外気温を特性項目としてもよい。
主軸受７（７ａ，７ｂ）の振動情報は、主軸受前記の加速度信号を周波数解析して得られる第１信号のうち、主軸６の回転数の整数倍に対応する周波数の信号強度として定義される主軸受前記の加速度レベルであってもよく、主軸受前記の加速度信号のピークを結ぶ包絡線信号を周波数解析して得られる第２信号のうち、主軸６の回転数の整数倍に対応する周波数の信号強度として定義される前記主軸受の包絡線レベルであってもよく、主軸受前記の半径方向速度信号を周波数解析して得られる第３信号のオーバーオール値であってもよく、主軸受前記の半径方向加速度信号を周波数解析して得られる第４信号のオーバーオール値であってもよく、又は、これらのうち２つ以上を含んでいてもよい。
【００３７】
また、例えば、同期発電機（発電機１０）を異常判定の対象とする場合、同期発電機（発電機１０）の巻線温度、自動電圧調節器の電圧、自動電圧調節器の電流、同期発電機（発電機１０）の界磁電圧、および、外気温を特性項目としてもよい。
【００３８】
また、例えば、発電機１０の回転軸を軸支する発電機軸受を異常判定の対象とする場合、発電機軸受の温度、発電機軸受の振動情報および外気温を特性項目としてもよい。
発電機軸受の振動情報は、発電機軸受の加速度信号を周波数解析して得られる第５信号のうち、発電機１０の回転軸の回転数の整数倍に対応する周波数の信号強度として定義される前記発電機軸受の加速度レベルであってもよく、発電機軸受の加速度信号のピークを結ぶ包絡線信号を周波数解析して得られる第６信号のうち、発電機１０の回転軸の回転数の整数倍に対応する周波数の信号強度として定義される発電機軸受の包絡線レベルであってもよく、発電機軸受の半径方向速度信号を周波数解析して得られる第７信号のオーバーオール値であってもよく、発電機軸受の半径方向加速度信号を周波数解析して得られる第８信号のオーバーオール値であってもよく、又は、これらのうち２つ以上を含んでいてもよい。
【００３９】
複数の特性値の組合せからなるデータセットは、単位期間Ｔｉ毎に取得する。単位期間Ｔｉの長さは、例えば１分〜６０分の範囲の期間である。各単位期間Ｔｉの長さは、例えば１０分間である。
【００４０】
対象データセット選別部２４は、データセット取得部２２において単位期間Ｔｉ毎に取得された複数のデータセットの中に、対象データセット判定条件を満たさないデータセットが含まれる時には、そのデータセットが含まれる単位期間Ｔｉを非対象単位期間とし、該非対象単位期間に属するデータセット全てを排除する。一方、対象データセット選別部２４は、データセット取得部２２において単位期間Ｔｉ毎に取得された複数のデータセットの全てが対象データセット判定条件を満たすデータセットである場合には、その単位期間Ｔｉを対象単位期間とし、該対象単位期間に含まれる複数のデータセットを対象データセットとして選別する。
風力発電設備１の運転状態の違いは偏差指標値に少なからず影響するので、この影響を排除した適切な異常判定が望まれる。そこで、上述のように対象データセット選別部２４によって風力発電設備１の運転状態が所定の場合のデータセットのみを選別して偏差指標値を算出することで、適切に異常判定を行うことができる。
【００４１】
対象データセット選別部２４が対象データセットを選別するための対象データセット判定条件は、下記条件（Ａ）又は下記条件（Ｂ）の少なくとも一方からなる第１判定条件である。
条件（Ａ）：風力発電設備１のロータ回転数ωｒが、カットイン風速時におけるロータ回転数ωｒ０以上かつ定格回転数ωｒ＿ｒａｔｅｄ未満の閾値ωｒ＿ｔｈ以上である。
条件（Ｂ）：風力発電設備１の発電機回転数ωｇが、電力系統への発電機の併入時における発電機回転数ωｇ０以下の閾値ωｇ＿ｔｈ以上である。
即ち、条件（Ａ）により、風力発電設備１の運転中以外のデータが対象データセットから排除される。また、条件（Ｂ）により、発電機１０の運転中以外のデータが対象データセットから排除される。
【００４２】
一実施形態では、対象データセット選別部２４が対象データセットを選別するための対象データセット判定条件は、第１判定条件に加えて、下記条件（Ｃ）又は下記条件（Ｄ）の少なくとも一方からなる第２判定条件をさらに含む。他の一実施形態では、対象データセット選別部２４が対象データセットを選別するための対象データセット判定条件は、下記条件（Ｃ）又は下記条件（Ｄ）の少なくとも一方からなる第２判定条件である。
条件（Ｃ）：風力発電設備１に一時的な出力制限が課せられていない。
条件（Ｄ）：風力発電設備１の運転状態に応じた温度変化を示す風力発電設備１の部位の温度変動率が閾値以下である。
すなわち、条件（Ｃ）により、風力発電設備１に一時的な出力制限が課せられ、例えば、発電機出力、風車翼のピッチ角及びロータ回転数といった特性項目間の相関関係が、出力制限前から変化したデータが対象データセットから排除される。また、条件（Ｄ）により、風力発電設備１の運転開始直後のみに起こる、風力発電設備１の部位（例えば、発電機の巻線や、発電機軸受）の特異的な温度変動が生じる期間のデータが対象データセットから排除される。
【００４３】
クラス分類部２６は、単位期間Ｔｉ内の複数のデータセットのうち、対象データセット選別部２４にて対象データセットとして選別されたデータセットを、風力発電設備１の出力を示すパラメータのクラス毎に分類するように構成される。クラス分類を行うことで、同一のクラスに属する対象データセット（信号空間を構成する）と基準データセット（単位空間を構成する）とからマハラノビス距離を算出することができ、これにより同じ状況下で取得されたデータ同士を比較することができるため、より的確な異常判定が可能となる。
クラス分類部２６におけるクラス分類は、例えば「風速」、「発電機出力」又は、「ロータ回転数」等、風力発電設備１の出力を示す、いずれかの特性項目に着目して行うことができる。図３に示す複数の特性項目を用いて風力発電設備１の異常判定を行う場合、例えば「発電機出力」に着目してクラス分類を行うことができる。
【００４４】
図４は、一実施形態に係るデータセットのクラスの分類を示す図であり、図３に示す特性項目を用いて評価を行う際に用いることができるクラス分類である。
一実施形態では、図４に示すように、発電機出力（横軸）について一定の幅の区間に区切られたＢＩＮｊ（ｊ＝１〜ｎ）が設定される。図４に示す例では、１００ｋＷの区間ごとにＢＩＮが設定され、ＢＩＮ１（６５０〜７５０ｋＷ）、ＢＩＮ２（７５０〜８５０ｋＷ）、ＢＩＮ３（８５０〜９５０ｋＷ）…等、各ＢＩＮ（クラス）毎に発電機出力の区間が割り当てられる。
各対象データセットは、各データセットに含まれる特性項目「発電機出力」の特性値（計測値）に応じて、いずれかのＢＩＮｊに振り分けられる。（すなわち、発電機出力に着目したクラス分類が行われる。）
【００４５】
平均処理部２８は、単位期間Ｔｉ毎に、クラス分類部２６によって各クラスに分類された対象データセットを平均化し、クラス（ＢＩＮｊ）毎に平均データセットＤｉｊを算出する。すなわち、各クラス毎に、対象データセットに含まれる特性項目の指標である特性値の各々について平均値を求め、各特性値の平均値の組合せからなる平均データセットＤｉｊを、単位期間Ｔｉ毎、かつクラス（ＢＩＮｊ）毎に算出する。
【００４６】
図５は、データセット取得部２２、対象データセット選別部２４、クラス分類部２６及び平均処理部２８の処理を繰り返すことで得られる複数の平均データセットＤｉｊの一例を示す表である。
図５に示す表は、以下のように取得したものである。
まず、データセット取得部２２にて、図３に示す６つの特性項目（発電機回転数、発電機出力、風速及び第１翼〜第３翼のピッチ角）からなるデータセットを、単位期間Ｔｉ（Ｔ１〜Ｔ５）にわたって取得した。なお、各単位期間Ｔ１〜Ｔ５の長さは１０分間である。
次に、取得したデータセットについて、クラス分類部２６にて、図４に示すＢＩＮｊ毎にクラス分類を行った。ＢＩＮｊは発電機出力によって分類されたクラスであり、各ＢＩＮｊの幅は１００ｋＷである。
平均処理部２８にて、ＢＩＮｊにクラス分類されたデータセットについて、単位期間Ｔｉ毎の平均値（１０分間平均）を算出し、平均データセットＤｉｊ（Ｄ１１〜Ｄ５５）を得た。各平均データセットＤｉｊ（Ｄ１１〜Ｄ５５）は、６つの特性項目の特性値の組合せ（１行ごとに示される）により構成される。
【００４７】
このようにして得られた平均データセットＤｉｊの各々は、後で説明する偏差指標値算出部４２にて、信号空間を構成するデータセットとしてマハラノビス距離の算出に用いられる。なお、偏差指標値算出部４２では、上述のようにして得られた平均データセットＤｉｊ毎に、マハラノビス距離が算出される。
【００４８】
基準データセット生成部３０は、風力発電設備１の健全性を評価するためのマハラノビス距離（偏差指標値）の算出において単位空間を構成するデータセット（基準データセット）を生成するように構成される。基準データセットは、上述した平均データセットＤｉｊと同じ特性項目についての特性値の組合せのうち、風力発電設備１の正常時における前記特性値の組合せである。
【００４９】
基準データセット取得部３２は、データセット取得部２２にて取得されるデータセットを構成する複数の特性項目の指標である特性値の組合せからなる複数のデータセットを単位期間Ｔｉ毎に取得するように構成される。
データセット取得部２２にて取得されるデータセットが図３に示される特性項目により構成されるとき、基準データセット取得部３２にて取得される基準データセットもまた、図３に示される特性項目により構成される。
【００５０】
基準データセット選別部３４は、基準データセット取得部３２において単位期間Ｔｉ毎に取得された複数のデータセットの中から、風力発電設備１の正常時における特性値の組合せであるデータセットを選別する。すなわち、各々の特性項目の特性値が、風力発電設備１の正常時であることを示す範囲の中に存在するデータセットのみを選別し、それ以外のデータセットは排除される。
【００５１】
さらに、基準データセット選別部３４は、対象データセット選別部２４と同様に、基準データセット取得部３２において単位期間Ｔｉ毎に取得された複数のデータセットの中に、基準データセット判定条件を満たさないデータセットが含まれる時には、そのデータセットが含まれる単位期間Ｔｉを非対象単位期間とし、該非対象単位期間に属するデータセット全てを排除する。一方、基準データセット選別部３４は、基準データセット取得部３２において単位期間Ｔｉ毎に取得された複数のデータセットの全てが基準データセット判定条件を満たすデータセットである場合には、その単位期間Ｔｉを対象単位期間とし、該対象単位期間に含まれる複数のデータセットを選別する。
基準データセット選別部３４による選別を行うことで、同一の条件を満たす対象データセットと基準データセットを用いて偏差指標値を算出することができる。よって、同じ条件下で取得されたデータ同士を比較することができるため、より的確な異常判定が可能となる。
【００５２】
基準データセット選別部３４がデータセットを選別するための基準データセット判定条件は、上記対象データセット判定条件と同じ条件であり、本実施形態において、下記条件（Ａ）又は下記条件（Ｂ）の少なくとも一方からなる第１判定条件である。
条件（Ａ）：風力発電設備１のロータ回転数ωｒが、カットイン風速時におけるロータ回転数ωｒ０以上かつ定格回転数ωｒ＿ｒａｔｅｄ未満の閾値ωｒ＿ｔｈ以上である。
条件（Ｂ）：風力発電設備１の発電機回転数ωｇが、電力系統への発電機の併入時における発電機回転数ωｇ０以下の閾値ωｇ＿ｔｈ以上である。
一実施形態では、基準データセット選別部３４が対象データセットを選別するための基準データセット判定条件は、第１判定条件に加えて、下記条件（Ｃ）又は下記条件（Ｄ）の少なくとも一方からなる第２判定条件をさらに含む。他の一実施形態では、基準データセット選別部３４が対象データセットを選別するための基準データセット判定条件は、下記条件（Ｃ）又は下記条件（Ｄ）の少なくとも一方からなる第２判定条件である。
条件（Ｃ）：風力発電設備１に一時的な出力制限が課せられていない。
条件（Ｄ）：風力発電設備１の運転状態に応じた温度変化を示す風力発電設備１の部位の温度変動率が閾値以下である。
【００５３】
基準データセットクラス分類部３６は、単位期間Ｔｉ内の複数のデータセットのうち、基準データセット選別部３４にて選別されたデータセットを、風力発電設備１の出力を示すパラメータのクラス毎に分類するように構成される。クラス分類の際のクラスの設定は、クラス分類部２６でのクラスの設定と同様である。
すなわち、本実施形態では、クラス分類部２６にけるクラス分類と同様に、図４に示すように、発電機出力（横軸）について一定の幅（１００ｋＷ）の区間に区切られたＢＩＮｊ（ｊ＝１〜ｎ）が設定される。そして、基準データセット選別部３４にて選別されたデータセットは、各データセットに含まれる特性項目「発電機出力」の特性値（計測値）に応じて、いずれかのＢＩＮｊに振り分けられる。（すなわち、発電機出力に着目したクラス分類が行われる。）
【００５４】
基準データセット平均処理部３８は、単位期間Ｔｉ毎に、基準データセットクラス分類部３６によって各クラスに分類されたデータセットを平均化し、クラス（ＢＩＮｊ）毎に平均基準データセットＳＤｉｊを算出する。すなわち、各クラス毎に、対象データセットに含まれる特性項目の指標である特性値の各々について平均値を求め、各特性値の平均値の組合せからなる平均基準データセットＳＤｉｊを、単位期間Ｔｉ毎、かつクラス（ＢＩＮｊ）毎に算出する。
【００５５】
以上に説明した基準データセット取得部３２、基準データセット選別部３４、基準データセットクラス分類部３６及び基準データセット平均処理部３８の処理を繰り返すことで、平均基準データセットＳＤｉｊを複数の単位期間Ｔｉにわたって集積することで、ＢＩＮｊ毎の基準データセットＤｒｅｆ＿ｊが得られる。
なお、図５に示すデータの集合を基準データセットとして用いることができる。この場合、平均データセットＤｉｊとして算出したデータセットを、平均基準データセットＳＤｉｊとすることができ、平均基準データセットＳＤｉｊを単位期間Ｔ１〜Ｔ５にわたって集積したものを、各ＢＩＮｊにおける基準データセットＤｒｅｆ＿ｊとする。
なお、後に説明する偏差指標値算出部４２では、ＢＩＮｊ毎の基準データセットＤｒｅｆ＿ｊを用いて、ＢＩＮｊ毎にマハラノビス距離が算出される。
【００５６】
基準データセットは、診断対象の風力発電設備１の運転開始前に予め取得されたものであってもよい。この場合、基準データセットは、診断対象の風力発電設備１と同一機種の風力発電設備を用いて、同一サイトにて取得された基準データセットであってもよい。
あるいは、基準データセットは、診断対象の風力発電設備１の運転を開始してから、当該風力発電設備１の運転中に取得されるようになっていてもよい。
【００５７】
健全性評価部４０は、データセット生成部２０及び基準データセット生成部３０で得られた平均データセットＤｉｊ及び基準データセットＤｒｅｆ＿ｊに基づいて、風力発電設備１の健全性の評価を行う。
【００５８】
偏差指標値算出部４２は、風力発電設備１の正常時における特性値の組合せからなる基準データセットＤｒｅｆ＿ｊと平均データセットＤｉｊとの偏差を示す偏差指標値として、マハラノビス距離ＭＤｉｊを算出する。すなわち、基準データセットＤｒｅｆ＿ｊにより構成される単位空間に対する、平均データセットＤｉｊにより構成される信号空間のマハラノビス距離ＭＤｉｊを偏差指標値として算出する。ここでは、同一のＢＩＮｊ（クラス）に属する基準データセットＤｒｅｆ＿ｊと平均データセットＤｉｊとの偏差を示すマハラノビス距離ＭＤｉｊが、各単位期間ＴｉについてＢＩＮｊ毎（クラス毎）に算出される。
【００５９】
なお、ＭＴ法では、目的に対して均質な集団（正常データの集合）を単位空間とし、対象データの単位空間の中心からの距離をマハラノビス距離（ＭＤ値）として求める。ＭＤ値が小さければその対象データは正常である可能性が大きく、ＭＤ値が大きければその対象データは異常である可能性が大きい。
【００６０】
ここで、図６及び図７を参照して、マハラノビス距離の算出方法について説明する。図６は、マハラノビス距離算出対象の単位空間データ（基準データセットの集合）の要素の一例を示す図であり、図７は、マハラノビス距離算出対象の信号空間データ（平均データセットの集合）の要素の一例を示す図である。
図６及び図７の「項目１」〜「項目ｎ」の各々は特性項目を示す。
また、「data(k)」及び「data(l)」におけるｋ及びｌは、それぞれ基準データセット及び平均データセットの番号を示す。すなわち、図６の単位空間データは、data(1)〜data(k)のｋ個の基準データセットを含んでいる。また、図７の信号空間データは、data(1)〜data(l)のｌ個の平均データセットを含んでいる。マハラノビス距離は平均データセット（信号空間データ）毎に計算されるので、図７の信号空間データを用いれば、data(1)〜data(l)の平均データセット毎に１つずつ、合計ｌ個のマハラノビス距離が算出される。
また、Ｘｎｋ及びＹｎｌは、それぞれ、各データセットにおける各項目（特性項目）の特性値を示す。
なお、マハラノビス距離を算出するための前提条件として、単位空間を構成するデータセット数ｋは、特性項目の数ｎよりも大きい必要がある（ｋ＞ｎ）。
【００６１】
0006320307.tif 000002
【数１】
0006320307.tif 000003
次に、上記式（Ａ）で算出した各項目毎の平均を用いて、下記式（Ｂ）により単位空間データについて共分散行列ＣＯＶ（ｎ×ｎ行列）を求める。
【数２】
0006320307.tif 000004
そして、信号空間データdata(1) 〜data(l)と、上記式（Ａ）により求めた平均及び上記式（Ｂ）により求めた共分散行列の逆行列を用いて、下記式（Ｃ）によりマハラノビス距離の２乗値Ｄ２を算出する。
【数３】
0006320307.tif 000005
さらに、単位空間データのＤ２の平均を１に正規化するためにＤ２をｎで除算してＭＤ２値を算出する。このＭＤ２値をマハラノビス距離（評価値）とする（ＭＤ２値＝Ｄ２／ｎ）。
【００６２】
偏差指標値算出部４２は、対象データセット選別部２４においてデータセット取得単位期間にわたって取得された複数のデータセットから選別された対象データセットを用いて、３０分〜６００分の範囲で、好ましくは６０分の期間で、マハラノビス距離を算出するようになっていてもよい。
【００６３】
異常判定部４３は、少なくとも一つのクラス（ＢＩＮｊ）においてマハラノビス距離ＭＤｉｊ（偏差指標値）が閾値を上回った場合、単位期間Ｔｉにおいて風力発電設備１の異常が生じたと判定するように構成される。
異常判定部４３は、予めマハラノビス距離ＭＤｉｊ（偏差指標値）の閾値が記憶された記憶部を有していてもよい。異常判定部４３は、偏差指標値算出部４２により算出されたマハラノビス距離を含む情報を受け取って、該マハラノビス距離と、記憶部に記憶された閾値とを比較して、算出されたマハラノビス距離が閾値よりも大きい場合に、風力発電設備１の異常が生じたと判定するようになっていてもよい。
【００６４】
図８は、図５に示されるデータセットを用いてマハラノビス距離を算出した結果を示す図である。図８において、横軸はデータセット番号であり、縦軸はマハラノビス距離（ＭＤ２値）である。
ここで、マハラノビス距離（ＭＤ２値）の算出対象は、図５の表中のデータセットＮｏ．１〜２５の各データセットである。
例えば、ＢＩＮ１に属するデータセットＮｏ．１〜５については、ＢＩＮ１の基準データセットＤｒｅｆ＿１により構成される単位空間データを用いてマハラノビス距離（ＭＤ２値）が算出される。同様に、ＢＩＮ２に属するデータセットＮｏ．６〜１０については、ＢＩＮ２の基準データセットＤｒｅｆ＿２により構成される単位空間データを用いてマハラノビス距離（ＭＤ２値）が算出される。他のデータセットについても、同様にマハラノビス距離（ＭＤ２値）が算出される。
【００６５】
なお、図５に記載される各特性項目（発電機回転数等）についての数値データは、実測値ではなく、検証用に作成したデータである。また、データセットＮｏ．８及びＮｏ．１８の発電機回転数（図５の表においては単に「回転数」と記載される）の欄にはカッコなしの数値とカッコ付きの数値の２種類の数値が記載されているが、単位空間データを構成する基準データセットの要素としてはカッコなしの数値を用い、信号空間データを構成する平均データセットの要素としては、カッコ内の数値を用いて各マハラノビス距離（ＭＤ２値）を算出した。
【００６６】
図８に示すグラフでは、Ｎｏ．８とＮｏ．１８以外のＭＤ２値は、０付近で小さい。これらのＭＤ値は、風力発電設備１が正常であることを示している。
これに対し、データセットＮｏ．８とＮｏ．１８から算出されたＭＤ２値は、他のデータセットから算出されたマハラノビス距離に比べて大きい。このことから、Ｎｏ．８及びＮｏ．１８のデータセットが取得された単位期間Ｔ３において、風力発電設備１に何らかの異常が発生したものと判定することができる。
なお、図８に示すようなマハラノビス距離の計算結果が得られる場合、異常判定部４３においてマハラノビス距離（ＭＤ２値）が５０〜１００の範囲内に閾値を設けておくことで、適切に風力発電設備１の異常を検出することができる。
【００６７】
寄与度算出部４４は、偏差指標値算出部４２により算出されたマハラノビス距離に対して、どの特性項目が寄与しているかを示す指標（寄与度）を算出する。また、異常要因特定部４５は、異常判定部４３により風力発電設備１に異常があると判定された場合に、寄与度算出部４４により算出された寄与度に基づいて、複数の特性項目のうち、異常の要因となっている特性項目を特定するように構成される。
ここで、項目１〜項目３からなる３つの特性項目を用いてマハラノビス距離を算出した場合の寄与度を算出し、異常要因を特定する方法の一例について図９を用いて説明する。図９は、寄与度の算出方法を説明するための表である。なお、図９において、各検証の行の項目１〜項目３の欄に記載されている○は、各検証において当該項目がマハラノビス距離又は以下に説明する分析用マハラノビス距離の算出に用いられることを示し、×は、各検証において当該項目がマハラノビス距離又は以下に説明する分析用マハラノビス距離の算出に用いられないことを示す。
【００６８】
まず、上述で説明したマハラノビス距離算出方法により、対象データセット及び基準データセットを構成する項目１〜項目３の全ての特性項目の特性値を用いて、マハラノビス距離（ＭＤ２値）を算出する（図９の表の検証１）。この際、特性項目数が３であるので、ＭＤ２値はＤ２／３である。
また、特性項目の中から、１つの着目特性項目を設定し、全特性項目から着目特性項目を除いた特性項目の特性値の組合せからなる分析用対象データセット及び分析用基準データセットを用いて、分析用マハラノビス距離（分析用偏差指標値）を算出する（図９の表の検証２〜４）。例えば、検証２では、項目１を着目特性項目と設定し、それ以外の項目２及び項目３を用いて分析用マハラノビス距離を算出する。この際、特性項目数が２であるので、ＭＤ２値はＤ２／２である。検証３及び検証４においては、それぞれ、項目２及び項目３を着目特性項目として、それ以外の２つの項目を用いて分析用マハラノビス距離を算出する。
【００６９】
そして、検証１で算出したマハラノビス距離に対して、検証２〜検証４で算出した分析用マハラノビス距離が占める割合に基づいて、各々の特性項目のマハラノビス距離への寄与度を算出する。
例えば、全ての特性項目（項目１〜項目３）を用いて算出したマハラノビス距離に対して、項目１を着目特性項目として検証２で算出したマハラノビス距離が占める割合が大きければ、項目１の寄与度は小さいと判断できる。逆に、全ての特性項目（項目１〜項目３）を用いて算出したマハラノビス距離に対して、項目１を着目特性項目として検証２で算出したマハラノビス距離が占める割合が小さければ、項目１の寄与度は大きいと判断できる。
【００７０】
あるいは、各着目特性項目についてＳＮ比利得率を求めることによって寄与度を算出するようにしてもよい。この場合、算出された各々のマハラノビス距離（ＭＤ２値）について、下記式（Ｄ）よりＳＮ比を求める。
【数４】
0006320307.tif 000006
表１０に示すように検証１〜検証４のそれぞれについて、求められるＳＮ比をη、η’１、η’２、η’３とする。そして、着目特性項目毎に、ＳＮ比利得率を算出する。項目１、項目２、項目３がそれぞれ着目特性項目である場合、ＳＮ比利得は、それぞれ、η−η’１、η−η’２、η−η’３である。
このようにして求めたＳＮ比利得は、マハラノビス距離における各特性項目の寄与度を示す。例えば、着目特性項目を項目１としたときのＳＮ比利得が大きい場合には、項目１はマハラノビス距離を大きくすることに寄与しており、風力発電設備１の異常要因となっている可能性があることを示す。また、例えば、着目特性項目を項目２としたときのＳＮ比利得が小さい場合には、項目２はマハラノビス距離の大きさに寄与しておらず、風力発電設備１の異常要因とはなっていないと判定される。
【００７１】
ＭＤ成分算出部４６は、偏差指標値算出部４２によって算出されたマハラノビス距離と、寄与度算出部４４によって算出された特性項目の各々の寄与度とに基づいて、特性項目の各々のマハラノビス距離成分を算出するように構成される。また、表示部４７は、ＭＤ成分算出部４６によって算出された特性項目の各々のマハラノビス距離成分の時系列変化を表示するように構成される。
図１０は、表示部４７により表示される各特性項目のマハラノビス距離成分の時系列変化の一例である。ここでは、項目１〜項目４の４つの特性項目によりマハラノビス距離が算出される場合の例を示している。また、図１０においては、予め設定されたマハラノビス距離の閾値が合わせて表示されている。
例えば、図１０においては、時刻がｔ１〜ｔ２の間の期間にマハラノビス距離が閾値を超えているため、当該期間において風力発電設備１に異常が発生していたと判定することができる。また、風力発電設備１に異常が発生していたと判定されたｔ１〜ｔ２の期間において項目４の寄与度が大きくなっているため、風力発電設備１の時刻ｔ１〜ｔ２における異常の要因は、項目４であると判定することができる。
【００７２】
以上、本発明の実施形態について説明したが、本発明は上述した実施形態に限定されることはなく、上述した実施形態に変形を加えた形態や、これらの形態を適宜組み合わせた形態も含む。
【００７３】
本明細書において、「ある方向に」、「ある方向に沿って」、「平行」、「直交」、「中心」、「同心」或いは「同軸」等の相対的或いは絶対的な配置を表す表現は、厳密にそのような配置を表すのみならず、公差、若しくは、同じ機能が得られる程度の角度や距離をもって相対的に変位している状態も表すものとする。
例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
また、本明細書において、一の構成要素を「備える」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。
【符号の説明】
【００７４】
１          風力発電設備
２          風車翼
３          ロータ
４          ハブ
５          ハブカバー
６          主軸
７          主軸受
８          トランスミッション
９          ドライブトレイン
１０        発電機
１２        ナセル
１４        タワー
２０        データセット生成部
２２        データセット取得部
２４        対象データセット選別部
２６        クラス分類部
２８        平均処理部
３０        基準データセット生成部
３２        基準データセット取得部
３４        基準データセット選別部
３６        基準データセットクラス分類部
３８        基準データセット平均処理部
４０        健全性評価部
４２        偏差指標値算出部
４３        異常判定部
４４        寄与度算出部
４５        異常要因特定部
４６        ＭＤ成分算出部
４７        表示部
１００      診断システム
【図１】
0006320307.tif 000007
【図２】
0006320307.tif 000008
【図３】
0006320307.tif 000009
【図４】
0006320307.tif 000010
【図５】
0006320307.tif 000011
【図６】
0006320307.tif 000012
【図７】
0006320307.tif 000013
【図８】
0006320307.tif 000014
【図９】
0006320307.tif 000015
【図１０】
0006320307.tif 000016