特許 6906354 風車発電機の疲労損傷量算出装置、風力発電システム、及び風車発電機の疲労損傷量算出方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6906354

(24)【登録日】2021年7月1日

(45)【発行日】2021年7月21日

(54)【発明の名称】風車発電機の疲労損傷量算出装置、風力発電システム、及び風車発電機の疲労損傷量算出方法

(51)【国際特許分類】

   G01M 99/00 20110101AFI20210708BHJP

   F03D 17/00 20160101ALI20210708BHJP

   G01N 3/32 20060101ALI20210708BHJP

【ＦＩ】

   G01M99/00 A

   F03D17/00

   G01N3/32 Z

【請求項の数】10

【全頁数】16

(21)【出願番号】特願2017-85500(P2017-85500)

(22)【出願日】2017年4月24日

(65)【公開番号】特開2018-185171(P2018-185171A)

(43)【公開日】2018年11月22日

【審査請求日】2020年1月20日

(73)【特許権者】

【識別番号】000003078

【氏名又は名称】株式会社東芝

(73)【特許権者】

【識別番号】317015294

【氏名又は名称】東芝エネルギーシステムズ株式会社

(74)【代理人】

【識別番号】110001092

【氏名又は名称】特許業務法人サクラ国際特許事務所

(72)【発明者】

【氏名】齊藤 和宏

(72)【発明者】

【氏名】田中 明

(72)【発明者】

【氏名】奥野 研一

(72)【発明者】

【氏名】橋場 豊

(72)【発明者】

【氏名】大迫 俊樹

【審査官】 松岡 智也

(56)【参考文献】

【文献】 特開２００４−０６０４７７（ＪＰ，Ａ）

【文献】 特開２０１６−１８８６１２（ＪＰ，Ａ）

【文献】 特開２０１６−２１７１３３（ＪＰ，Ａ）

【文献】 特開２００５−２４０７８５（ＪＰ，Ａ）

【文献】 国際公開第２０１６／０４２６５２（ＷＯ，Ａ１）

【文献】 米国特許出願公開第２０１３／０１６６０８２（ＵＳ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｍ １３／００−１３／０４５、９９／００

Ｆ０３Ｄ １／００−８０／８０

(57)【特許請求の範囲】

【請求項1】

風車発電機の設置場所で一定期間中に計測された風速の平均及び風速の変化を表す風況情報を取得する風況情報取得部と、
前記取得された風況情報に基づいて、前記風車発電機のブレードの疲労損傷量を算出する疲労損傷量算出部と、
前記一定期間中に計測された前記風車発電機の所定の部位におけるひずみを表す情報を取得するひずみ情報取得部と、
前記取得されたひずみを表す情報に基づいて、前記一定期間中に前記ブレードに加わる荷重の平均及び荷重の変化を表す第１の荷重情報を算出する第１の荷重情報算出部と、
前記取得された風況情報に基づいて、前記一定期間中に前記ブレードに加わる荷重の平均及び荷重の変化を表す第２の荷重情報を算出する第２の荷重情報算出部と、
それぞれ算出された前記第１及び第２の荷重情報に基づいて、前記算出された疲労損傷量を補正する疲労損傷量補正部と、
前記一定期間の単位で前記疲労損傷量補正部によりそれぞれ補正された複数の疲労損傷量を積算する疲労損傷量積算部と、
を備える風車発電機の疲労損傷量算出装置。

【請求項2】

前記風況情報と前記ブレードの疲労損傷量との相関関係、前記所定の部位におけるひずみと前記第１の荷重情報との相関関係、前記風況情報と前記第２の荷重情報との相関関係、並びに、前記第１及び第２の荷重情報と前記ブレードの疲労損傷量との相関関係、を予め記憶する相関関係記憶部をさらに備え、
前記疲労損傷量算出部、前記第１の荷重情報算出部、前記第２の荷重情報算出部、及び、前記疲労損傷量補正部は、前記相関関係記憶部に記憶された内容を参照しつつ、前記疲労損傷量の算出、前記第１の荷重情報の算出、前記第２の荷重情報の算出、前記疲労損傷量の補正、をそれぞれ行う、
請求項１に記載の風車発電機の疲労損傷量算出装置。

【請求項3】

前記風況情報は、前記一定期間中の最大風速と最小風速との差を表す情報を含む、
請求項１又は２に記載の風車発電機の疲労損傷量算出装置。

【請求項4】

前記風況情報は、前記一定期間中の風速の変化の標準偏差を風速の平均で除した乱流強度を表す情報を含む、
請求項１から３までのいずれか１項に記載の風車発電機の疲労損傷量算出装置。

【請求項5】

前記ひずみ情報取得部は、前記風車発電機が備える支柱のひずみを表す情報を取得する、
請求項１から４までのいずれか１項に記載の風車発電機の疲労損傷量算出装置。

【請求項6】

前記積算された疲労損傷量を出力する疲労損傷量出力部、
をさらに備える請求項１から５までのいずれか１項に記載の風車発電機の疲労損傷量算出装置。

【請求項7】

前記風況情報取得部により取得された風況情報と、前記風車発電機の設計時に想定された風速の平均及び風速の変化を表す風況情報と、をあわせて出力する風況情報出力部、
をさらに備える請求項１から６までのいずれか１項に記載の風車発電機の疲労損傷量算出装置。

【請求項8】

前記第１の荷重情報及び前記第２の荷重情報と、前記第１の荷重情報と前記第２の荷重情報との差異を表す情報と、のうちの少なくともいずれかを出力する荷重情報出力部、
をさらに備える請求項１から７までのいずれか１項に記載の風車発電機の疲労損傷量算出装置。

【請求項9】

請求項１から８までのいずれか１項に記載の風車発電機の疲労損傷量算出装置と、
前記風車発電機と、
を備える風力発電システム。

【請求項10】

風車発電機の設置場所で一定期間中に計測された風速の平均及び風速の変化を表す風況情報を取得するステップと、
前記取得された風況情報に基づいて、前記風車発電機のブレードの疲労損傷量を算出するステップと、
前記一定期間中に計測された前記風車発電機の所定の部位におけるひずみを表す情報を取得するステップと、
前記取得されたひずみを表す情報に基づいて、前記一定期間中に前記ブレードに加わる荷重の平均及び荷重の変化を表す第１の荷重情報を算出するステップと、
前記取得された風況情報に基づいて、前記一定期間中に前記ブレードに加わる荷重の平均及び荷重の変化を表す第２の荷重情報を算出するステップと、
それぞれ算出された前記第１及び第２の荷重情報に基づいて、前記算出された疲労損傷量を補正するステップと、
前記一定期間の単位で、それぞれ補正された複数の疲労損傷量を積算するステップと、
を有する風車発電機の疲労損傷量算出方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明の実施形態は、風車発電機の疲労損傷量算出装置、風力発電システム、及び風車発電機の疲労損傷量算出方法に関する。

【背景技術】

【0002】

風車発電機を運用する発電事業者には、運用コストの低減が要求されている。風車発電機の運用コストでは、点検や補修、故障による停止を考慮する必要がある。この運用コストは、風車発電機導入後の時間の経過に伴い、急激に増加することが知られている。つまり、風車発電機は、絶えず変動する風の影響を受け続けるため、風車発電機本体の各部に繰り返し加わる負荷によって疲労損傷や摩耗損傷が生じ、この結果、運用コストの増加を招くことになる。

【0003】

そこで、風車発電機の設計段階では、このような風の影響を十分に考慮して、強度評価が行われる。具体的には、風車発電機の設置予定場所で事前に風況情報を収集し、想定される平均風速や風速の変動に対して十分な強度を持つ風車発電機を設計する。例えばＩＥＣの風車標準では、風の強度は、風速を表す指標である参照風速が３段階あり、一方、風の変動を表す指標である乱流強度が３段階あり、さらにこれらを組み合わせて９段階に分類されている。したがって、風車発電機は、設置予定場所の風況や実際の運用を想定して設計される。

【0004】

例えば、風車発電機を構成しているタワーやナセルなどの構造物は、適切な構造設計をすることで、風によって生じる応力を低減して十分な疲労強度を確保できる。すなわち、タワーやナセルは、経時的な疲労によって破壊する可能性が低いため、例えば風車発電機の運転中や運転停止中に、風の影響で生じる応力を実測したり、疲労損傷の状況を計算したりする必要性などがほとんどない。

【0005】

これに対して、風車発電機のブレードは、疲労強度に対する裕度を大きくとることができないため、長期の運用や実際の風速の変動によって疲労損傷が生じる可能性が高い構造物である。したがって、ブレードは、経時変化する疲労損傷の状況を正しく予測する必要性がある。

【0006】

なお、風車発電機に内蔵されるギアや軸受などは、風による変動荷重によって摩耗損傷を受ける可能性が高いものの、振動計などのセンサを設けることで重大な損傷に至る前に軽い状況の損傷を検出することが可能となる。つまり、風車発電機の実際の運用や風速変動からブレードの損傷量を推定できれば、発電事業者は、風車発電機の運転制御や、点検及び補修を合理化して、運用コストを低減することができる。

【0007】

ここで、ナセルなどに取り付けた風向風速計や荷重センサの出力に基づき算出した風車各部に働く応力から、現在の疲労損傷率若しくは劣化度を演算により求め、さらにこれらを、予め定めた疲労劣化スケジュールテーブルと対比させながら、風車の運転を制御する第１及び第２の技術が提案されている。これらの技術によれば、疲労劣化スケジュールに対して現在の疲労損傷率に余裕がある場合とそうでない場合とで風車の制御を変えることによって、変動荷重による風車ブレードなどの疲労劣化を防止し、理想的な環境下で風車を運転することができる。

【0008】

また、複数の風力発電装置を有する風力発電システムにおいて、劣化度に基づき風力発電装置をグループ化し、高劣化した風力発電装置群に含まれる各風力発電装置を優先的に電力低減制御する第３の技術が知られている。つまり、この技術は、電力系統の要求に応じた有効電力の出力制御を行いつつ、ウインドファーム（大規模集中風力発電所）としての寿命を延ばすことができる。

【先行技術文献】

【特許文献】

【0009】

【特許文献1】特許第４２４１６４４号公報

【特許文献2】特許第４９３９５０８号公報

【特許文献3】特許第５１６７３６５号公報

【発明の概要】

【発明が解決しようとする課題】

【0010】

ところで、上記した第１〜第３の技術を適用して、風力発電装置やウインドファームの運用を最適化するためには、風力発電装置の現状の疲労損傷率を正確に計算することが前提となる。具体的には、第１の技術は、ナセルに取り付けた加速度センサから得た風速と風向きとから、各部に加わる荷重を表す荷重波形を求め、その荷重波形から応力を算出する。また、第１の技術では、風力計ではなく各種センサからの情報を基に同様に応力を算出することや、直接応力を計測することについても考えられている。一方、第２の技術では、荷重の算出について、より詳細に考えられており、風力計によって得られた情報と運転制御、運転状況に関する情報から、荷重時系列データを作成し、そこから応力時系列データ、最大主応力時系列データを計算している。

【0011】

これら第１、第２の技術を用いる疲労損傷率の計算方法は、センサによって実測した情報、例えば風速や風向などから、荷重波形を求め、そして疲労損傷の評価に必要となる応力波形を計算し、設計条件として想定した風況に対して、各部の応力波形を計算するＩＣＥの風車標準に則った方法であるといえる。

【0012】

しかしながら、この計算方法は、常に短時間で変化する風速や運転状況のデータを用いて計算が行われるため、高性能なセンサ及び計算機を用いて、非常に高速な演算が必要になる。また、風速の変化に追従したリアルタイムでの演算が難しい場合には、一時的にデータを保存しておく記憶装置なども必要になる。したがって、これらの装置を導入する必要性から、上記した方法は、風車発電機の運用コストを低減させることが難しくなっている。

【0013】

さらに、第１、第２の技術を用いる疲労損傷率の計算方法は、荷重や応力を算出する計算過程は、設計段階の計算とほぼ同様に行われる。風力発電装置の運用中には、多様な運転、機器の異常、突風や台風など、特殊な風況の組合せによって、様々な荷重、応力が風力発電装置の各部に生じることが考えられるが、これらを正確に予測することは難しい。詳述すると、設計者が設計段階で風力発電装置の各部に作用するものとして想定（計算）した荷重や応力と、実際に風力発電装置の各部に発生した荷重や応力と、の間に違いがあったとしても、発電事業者は、これらの違いを見出すことが困難である。

【0014】

また、第３の技術は、発電機出力から風力発電装置の劣化（疲労）度を推定するものであり、第１、第２の技術に較べると簡便な方法がとられている。しかしながら、発電機出力から風力発電装置の劣化度を計算する場合の評価の精度は、発電機出力と劣化度との両者の相関データ量に依存するため、発電事業者からみて十分な評価結果が得られていない可能性がある。

【0015】

このように、現状の技術は、風車発電機の実際の疲労損傷を評価できていないおそれがある。この一方で、風車発電機の疲労損傷の予測精度を向上させることは、例えば多数のセンサの高速サンプリングで得た膨大なデータの分析などを必要とし、この場合、結果的に運用コストの増加を招くことになる。

【0016】

そこで、本発明が解決しようとする課題は、風車発電機における疲労寿命の解析精度を容易に高めることができる風車発電機の疲労損傷量算出装置、風力発電システム、及び風車発電機の疲労損傷量算出方法を提供することである。

【課題を解決するための手段】

【0017】

実施の形態に係る風車発電機の疲労損傷量算出装置は、風況情報取得部、疲労損傷量算出部、ひずみ情報取得部、第１及び第２の荷重情報算出部、疲労損傷量補正部、並びに疲労損傷量積算部を備えている。風況情報取得部は、風車発電機の設置場所で一定期間中に計測された風速の平均及び風速の変化を表す風況情報を取得する。疲労損傷量算出部は、前記取得された風況情報に基づいて、前記風車発電機のブレードの疲労損傷量を算出する。ひずみ情報取得部は、前記一定期間中に計測された前記風車発電機の所定の部位におけるひずみを表す情報を取得する。第１の荷重情報算出部は、前記取得されたひずみを表す情報に基づいて、前記一定期間中に前記ブレードに加わる荷重の平均及び荷重の変化を表す第１の荷重情報を算出する。第２の荷重情報算出部は、前記取得された風況情報に基づいて、前記一定期間中に前記ブレードに加わる荷重の平均及び荷重の変化を表す第２の荷重情報を算出する。疲労損傷量補正部は、それぞれ算出された前記第１及び第２の荷重情報に基づいて、前記算出された疲労損傷量を補正する。疲労損傷量積算部は、前記一定期間の単位で前記疲労損傷量補正部によりそれぞれ補正された複数の疲労損傷量を積算する。

【発明の効果】

【0018】

本発明によれば、風車発電機における疲労寿命の解析精度を容易に高めることが可能な風車発電機の疲労損傷量算出装置、風力発電システム、及び風車発電機の疲労損傷量算出方法を提供することができる。

【図面の簡単な説明】

【0019】

【図1】第１の実施形態に係る風力発電システムの構成を概略的に示す図。

【図2】比較例の疲労損傷量算出方法を示すフローチャート。

【図3A】図１の風力発電システムにおける風車発電機の風況計測器で計測した１０分間の風速の変化を示す図。

【図3B】図３Ａの１０分間に風車発電機のブレードに加わる荷重の変化を示す図。

【図3C】図３Ａの１０分間に風車発電機のブレードに作用するモーメントを示す図。

【図4A】図３Ｂの１０分間の荷重の変化を、レインフロー法によってサイクルカウントして、ヒストグラムで整理した結果を示す図。

【図4B】図４Ａのヒストグラムの分布を関数化して確率密度にした結果を示す図。

【図5A】平均風速を基準にした荷重変動分布の形状パラメータの分析例を示す図。

【図5B】乱流強度を基準にした荷重変動分布の形状パラメータの分析例を示す図。

【図6】疲労評価の対象部位の疲労寿命線図を例示した図。

【図7】平均風速と乱流強度との関係を示すマトリックス。

【図8A】乱流強度及び平均風速と疲労損傷量との関係を示す図。

【図8B】乱流強度と疲労損傷量との関係を示す図。

【図9A】風車発電機の運転中の平均風速と風速変化との関係を示すマトリックス。

【図9B】風車発電機の停止中の平均風速と風速変化との関係を示すマトリックス。

【図10】風車発電機の運転中の平均風速と疲労損傷量との関係を示す図。

【図11】風車発電機の運転中の平均風速と荷重変化との関係を示す図。

【図12】ひずみを表す情報又は風況情報に基づくブレードへの荷重変化と、平均風速と、の関係を示す図。

【図13】図１の風車発電機の疲労損傷量算出装置による疲労損傷量算出方法を示すフローチャート。

【図14】第２の実施形態に係る風力発電システムの構成を概略的に示す図。

【発明を実施するための形態】

【0020】

以下、実施の形態を図面に基づき説明する。
＜第１の実施の形態＞
図１に示すように、本実施形態に係る風力発電システム１０は、風車発電機（風力発電装置）２０及び風車発電機の疲労損傷量算出装置３０を備えている。風車発電機２０は、筐体としてのナセル１５と、ナセル１５を下方から支持する支柱であるタワー１４と、複数のブレード（風車ブレード）１２と、ハブを介して複数のブレード１２を支持するボス１９と、を主に備えている。

【0021】

ナセル１５は、回転軸１７、変速機構１６、発電機本体１８を内蔵する。回転軸１７の一端部は、ボス１９に固定され、回転軸１７の他端部は、変速機構１６と連結されている。変速機構１６は、適宜のカップリング機構などを介して発電機本体１８と連結されている。例えば３つ設けられたブレード１２は、それぞれの基端部がボス１９（回転軸１７）の軸まわりに１２０度（ｄｅｇ）の間隔をおいて各々固定されている。

【0022】

ここで、風車発電機２０の運転中（発電中）には、各ブレード１２は、風力を受けて回転するように、各ブレード本体の圧力面（正圧面）及び負圧面の向きが調整される。なお、風車発電機２０の運転停止中（非発電中）には、各ブレード１２は、風力を受けても例えば回転しないように、各ブレード本体の圧力面及び負圧面の向きが調整される。このような複数のブレード１２、並びにボス１９及び回転軸１７は、ナセル１５に対して回転可能に支持されている。

【0023】

つまり、風車発電機２０の運転中において、ボス１９及び回転軸１７と共に一体となって回転する複数のブレード１２は、風力から得た流体エネルギを回転エネルギに変換する。また、回転軸１７の回転エネルギ（駆動力）は、変速機構１６により減速又は増速されて発電機本体１８へ伝達される。発電機本体１８は、伝達されたこの回転エネルギを用いて発電を行う。

【0024】

一方、ナセル１５には、風速計や風向計などの機能を有する風況計測器（風況センサ）２３が設けられている。風向や風速を計測可能な風況計測器２３は、風車発電機２０の設置された設置場所で、例えば風速の平均及び風速の変化を表す風況情報を計測し、計測したこの風況情報（計測結果）を出力する。なお、風況計測器２３は、単に、風速や風向を計測し、この計測結果を、計測した時刻情報と共に出力するものであってもよい。一方、タワー１４には、ひずみゲージ２１、２２や、作業者が出入りする出入口２４が設けられている。

【0025】

ひずみゲージ２１、２２は、タワー１４の胴部（本体部分）に取り付けられており、タワー１４の例えば軸方向のひずみを各々計測する。なお、タワー１４の例えば径方向のひずみを計測させるように、ひずみゲージ２１、２２を取り付けてもよい。また、ひずみゲージ２１とひずみゲージ２２とは、図１に示すように、タワー１４における、例えば高さの異なる位置にそれぞれ取り付けられている。ひずみゲージ２１、２２の計測結果を基づき、風力によりタワー１４やブレード１２に加わる荷重やその荷重の変化を算出することが可能となる。さらに、ひずみゲージ２１、２２の計測結果と、タワー１４の各部の寸法値と、ひずみゲージ２１、２２の取り付け位置と、に基づいて、タワー１４に加わる曲げモーメントなどを算出することも可能となる。

【0026】

ここで、タワー１４に対して曲げモーメントが加わる方向は、ナセル１５の取り付けの向きや位置によって影響を受ける。しかしながら、例えば２つのひずみゲージを、タワー１４の高さ方向において同一の高さ位置で、かつタワー１４を３６０度（ｄｅｇ）、周回する周回角度方向において互いを９０度（ｄｅｇ）ずらした位置に、それぞれ設置して、タワー１４の軸方向のひずみを計測させることで、タワー１４において、曲げ応力が最大となる方向を求めることができる。また、９０度取り付け位置をずらしたこのような２つのひずみゲージに加え、上述したように、取り付け高さをずらしたひずみゲージをさらに追加してもよい。

【0027】

なお、ひずみゲージの取り付け（貼り付け）位置は、ブレード１２に加わる力に応じて曲げモーメントが大きくなるタワー１４の低層階が好ましいものの、タワー１４の剛性がタワー１４の周回方向に応じて変わらないことがより好ましいので、図１に示すように、出入口２４が設けられた高さ位置を避けたほうがよい。

【0028】

ところで、図２は、上述した第１の技術における疲労寿命解析プロセスを比較例として示すフローチャートである。この比較例には、風車の劣化度を計算して、現在の劣化度と、長期間にわたる劣化度スケジュールと、を対比させた情報を用いて、風車を制御する方法が開示されている。この比較例では、現在の劣化度を計算する方法として、風速計などのセンサから取得した情報を用いて（Ｓ２１）、風車に加わる荷重波形を計算し（Ｓ２２）、この荷重波形を基に座標変換を行って特定の係数を算出し（Ｓ２３）、この算出結果から風車に生じる応力波形を求める（Ｓ２４）。さらに、この応力波形からレインフローカウント法を用いて応力頻度を計算し（Ｓ２５）、計算した応力頻度から、個々の繰返し応力に対する疲労損傷率を線形累積損傷則から求める（Ｓ２６）。一連のこの疲労損傷量算出方法は、ＩＥＣ（International Electrotechnical Commission：国際電気標準会議）の風車設計標準などにおける一般的な疲労設計法に準拠した公知の技術であるといえる。

【0029】

図３Ａは、ある一定期間の計測例として、風況計測器２３で計測した１０分間の風速の変化を示している。また、図３Ｂは、上記の１０分間にブレード１２に加わる荷重の変化を示し、さらに、図３Ｃは、上記の１０分間にブレード１２に作用するモーメントの変化を示している。このように風速は絶えず変化し、それに伴い、ブレード１２に作用する荷重、モーメントも変動する。このような荷重及びモーメントの変動から、疲労評価の対象
部位の応力変動とそれらの回数（繰り返し回数）とを計算すれば、疲労損傷評価を行うことが可能となる。しかしながら、このような高速に変化する事象を逐一計算して累積する疲労損傷評価は、高性能な計算機やセンサが必要になる。

【0030】

また、風速の変化から風車のブレードに加わる繰返し応力を算出する場合、いわゆるタワーシャドウ効果（回転中のブレードがタワーの前を通過するときに、ブレードの回転が減速する現象）などの影響、つまりブレードの回転に伴う応力変化、が考慮されない算出結果となる。各ブレードの例えばルート部（根元部）に、ひずみゲージなどのセンサを取り付ければ、上記のタワーシャドウ効果の影響なども、繰返し応力の算出結果に加味されることになるが、センサの使用個数が増え、これに伴い計算負荷も増加することになる。また、この場合、複雑な計算を要することになるので、発電事業者にとっては、疲労損傷評価の妥当性が判断し難いものとなる。

【0031】

そこで、本実施形態に係る風車発電機の疲労損傷量算出装置３０は、例えば、高速サンプリングで得た膨大なデータを分析しなくても、発電事業者からみて妥当性のある風車発電機の疲労損傷評価を行えるものである。本実施形態の疲労損傷量算出装置３０は、風速の変化に応じて、疲労損傷率を逐一計算して累積するのではなく、ある一定期間（所定期間）毎に疲労損傷率を求め、最終的に一定期間毎の疲労損傷率を積算する。

【0032】

具体的には、本実施形態の風車発電機の疲労損傷量算出装置３０は、図１に示すように、風況情報取得部３１、疲労損傷量算出部３２、ひずみ情報取得部３３、第１の荷重情報算出部３４、第２の荷重情報算出部３５、疲労損傷量補正部３６、疲労損傷量積算部３７、疲労損傷量出力部３８、並びに相関関係記憶部３９を備えている。これら風況情報取得部３１や疲労損傷量算出部３２を始めとする各構成要素（図１中に機能ブロックで示された構成要素）は、各種の電子部品を組み合わせて構成されるハードウェアによって実現されていてもよいし、所定のプログラムを実行して得られるソフトウェアによって実現されていてもよい。

【0033】

風況情報取得部３１は、風車発電機２０の設置場所で一定期間中に計測された風速の平均（平均風速）及び風速の変化（風速の変化幅）を表す風況情報を取得する。つまり、風況情報取得部３１は、風車発電機２０のナセル１５に設けられた風況計測器２３より例えば出力される上記の風況情報を取得する構成であってもよいし、風況計測器２３が計測した例えば風速や風向及びその計測時刻を風況計測器２３から受け取り、受け取ったこれらの情報を基に風況情報を算出して当該風況情報を取得する構成であってもよい。

【0034】

ここで、上記した風況情報は、一定期間中の最大風速と最小風速との差を表す情報や、一定期間中の風速の変化の標準偏差（変動風速の標準偏差）を風速の平均で除した乱流強度を表す情報を含む概念である。本実施形態の風況情報取得部３１は、例えば１０分間の平均風速と乱流強度とを代表値（指標値）として取得する。

【0035】

一方、疲労損傷量算出部３２は、風況情報取得部３１により取得された風況情報に基づいて、風車発電機２０のブレード１２の疲労損傷量を算出する。疲労損傷量は、疲労損傷率、寿命消費量、寿命消費率なとどと同じ意味を表す用語であり、例えば全く損傷していないノーダメージの部品が数値「０」又は「０％」などで表現され、例えば完全破壊の状態の部品が数値「１」又は「１００％」などで表現される。本実施形態の疲労損傷量算出部３２は、風況情報取得部３１により取得された平均風速と乱流強度とを用いて、疲労損傷量（疲労損傷量の基準値）を算出する。

【0036】

図４Ａは、図３Ｂで示した１０分間の荷重変化（ブレード１２に加わる荷重変動）を、レインフロー法によってサイクルカウントし、荷重変動Ｆａと、各荷重変動Ｆａの頻度を足し合わせた合計回数Ｎとを求め、これをヒストグラムで整理した結果である。上記のレインフロー法は、レインフロー計数法（rainflow counting method）などともいい、不規則な繰り返し変動荷重を受ける構造体において、疲労寿命を予測するための応力頻度、若しくは、ひずみ頻度の計数法の一つである。図４Ａに示すように、ブレード１２に作用する荷重変動Ｆａは、小さいＦａの回数（頻度）が多く、まれに大きなＦａがある。

【0037】

一方、図４Ｂは、ヒストグラムの分布を関数化して確率密度にした結果を示している。荷重変動Ｆａの分布をモデル化するのに適した確率関数としては、ワイブル分布などが挙げられる。ワイブル分布は、２つのパラメータ、分布の形を決める形状パラメータ、分布の幅（横軸の大きさ）を決める寸法パラメータによって定義できる。これにより、図３Ｃに示した１０分間のモーメントの変化も、荷重変動と同様に、モーメント変動Ｍａの大きさとその頻度とを示す分布をモデル化することができる。

【0038】

すなわち、図３Ａ〜図３Ｃに示した風速、荷重、モーメントは、それぞれ１０分間の平均風速、１０分間の乱流強度、１０分間の荷重変動の合計回数、１０分間の荷重変動分布の形状パラメータ、１０分間の荷重変動分布の尺度パラメータ、１０分間のモーメント変動の合計回数、１０分間のモーメント変動分布の形状パラメータ、１０分間のモーメント変動分布の尺度パラメータ、という代表値に変換できる。

【0039】

このような処理を、別の１０分間、つまり、平均風速、乱流強度が異なる場合についても行って分析すると、平均風速や乱流強度と合計回数、各形状パラメータ、各尺度パラメータの関係を分析できる。この結果、平均風速、乱流強度とパラメータとの明確な相関が明らかになれば、風況から荷重変動やモーメント変動が求められることになる。例えば、尺度パラメータは、平均速度と乱流強度とを掛け合わせたパラメータで整理すると、良好な相関が得られる。一方、明確な相関が認められない場合もある。

【0040】

図５Ａ、図５Ｂは、荷重変動分布の形状パラメータの分析例を示している。図５Ａ、図５Ｂに示すように、平均風速や乱流強度との間に、全てのプロットで成立する相関が認められない場合、平均風速や乱流強度の範囲を区切ってマトリックス化し、風況のパラメータを用いて関数化する。なお、疲労損傷評価に際して要求される精度は、パラメータ毎に異なる。また、要求される疲労損傷評価の精度によっても異なる。

【0041】

続いて、風況情報（平均風速及び乱流強度）から、１０分間の荷重変動の合計回数、１０分間の荷重変動分布の形状パラメータ、１０分間の荷重変動分布の尺度パラメータを求めることができるようになった以降の疲労損傷評価プロセスについて説明する。ワイブル分布の形状パラメータと尺度パラメータとを用いると、荷重変動の値がある範囲をとる場合の確率を計算できる。１０分間の合計回数を掛ければ、ある値の荷重変動の１０分間の回数となる。モーメントの場合も同様である。

【0042】

図６は、疲労評価の対象部位の疲労寿命線図を例示したものである。つまり、図６は、疲労評価の対象部位における材料の設計Ｓ−Ｎ線図（応力振幅と破断回数の関係）と荷重が加わった場合に疲労評価対象部位に生じる応力との関係を用いて、荷重変動Ｆａと疲労評価対象部位（本実施形態ではブレード１２）の破断回数Ｎｆとの関係に直したものである。図６に示すように、荷重変動Ｆａの大きさが決まれば、破断回数Ｎｆが決定される。したがって、ある値の荷重変動に対応した破断回数がわかれば、累積疲労損傷則にて疲労損傷量を求めることが可能となる。すなわち、ある１０分間の風況下で取り得る、荷重変動の範囲とそれらの各回数（破断回数）を推定することができ、これにより、１０分間の疲労損傷量（疲労損傷量の予測値）を算出できる。

【0043】

様々な条件（個々のケース）について整理すれば、図７に示すように、平均風速及び乱流強度（風況情報）から疲労損傷量を求める例えば風車発電機運転中のマトリックス（風況情報とブレード１２の疲労損傷量との相関関係）を作成することが可能となる。各マトリックスの内容は、図８Ａ、図８Ｂに示すように、ケース毎に区分された、平均風速及び乱流強度を表す風況情報から疲労損傷量を算出するための評価式で構成される。

【0044】

ここで、同じ風況でもブレード１２に加わる繰返し応力は、風車発電機２０の運転状況によって変わる。そこで、上述したマトリックスは、図９Ａ、図９Ｂに示すように、風車発電機２０の運転中及び運転停止中といった運転状況に応じて変更を加えたものが予め用意されている。つまり、図９Ａ、図９Ｂは、風車発電機２０の運転中、運転停止中において、平均風速と風速変化（乱流強度など）との相対的な関係を個々のケース毎に区分したマトリックスを例示している。さらに、図１０は、風車発電機２０の運転中の状況において、ケース毎の平均風速と疲労損傷量との相関関係を表している。図１０の相関関係では、実線に比べて、破線は、風速の変化が大きく、疲労損傷量も大きくなっている。図１０に示すような情報を利用することで、ケース及び平均風速から疲労損傷量を求めることが可能となる。

【0045】

すなわち、図１に示すように、相関関係記憶部３９は、平均風速及び乱流強度を表す風況情報とブレード１２の疲労損傷量との相関関係（例えば図８Ａ、図８Ｂ、図１０などに例示した情報）を予め記憶している。さらに、疲労損傷量算出部３２は、相関関係記憶部３９に記憶された内容を参照しつつ、風車発電機２０におけるブレード１２の疲労損傷量を算出する。

【0046】

また、本実施形態に係る風車発電機の疲労損傷量算出装置３０は、疲労評価の信頼性を向上させるために、風況情報を基に算出したブレード１２に加わり得る荷重（荷重の基準値）とひずみの計測値を基に算出したブレード１２に加わり得る荷重とを比較する。つまり、図１に示すように、ひずみ情報取得部３３は、一定期間中に計測された風車発電機２０の所定の部位（本実施形態では風車発電機２０が備えた支柱であるタワー１４）におけるひずみを表す情報を取得する。具体的には、ひずみ情報取得部３３は、タワー１４に取り付けられたひずみゲージ２１、２２を少なくとも含む２以上のひずみゲージの計測結果（例えば上記した９０度取り付け位置をずらした２つのひずみゲージの計測結果）を取得する。

【0047】

また、第１の荷重情報算出部３４は、ひずみ情報取得部３３により取得されたひずみを表す情報に基づいて、一定期間中にブレード１２に加わる荷重の平均及び荷重の変化を表す第１の荷重情報を算出する。さらに、第２の荷重情報算出部３５は、風況情報取得部３１により取得された風況情報に基づいて、一定期間中にブレード１２に加わる荷重の平均及び荷重の変化を表す第２の荷重情報を算出（推定）する。また、疲労損傷量補正部３６は、第１及び第２の荷重情報算出部３４、３５によりそれぞれ算出された第１及び第２の荷重情報に基づいて、疲労損傷量算出部３２により算出された疲労損傷量を補正する。

【0048】

ここで、図１１は、前記した図１０に示した例と同じケースにおいて、平均風速と、ブレード１２に加わり得る荷重変化（荷重変化の推定値）と、の相関関係を表している。一方、図１２は、計測されたひずみを表す情報に基づくブレード１２への荷重変化ΔＦ２及び風況情報を基に推定したブレード１２への荷重変化ΔＦ１と、平均風速と、の相関関係を示している。図１２に示すように、風況情報を基に推定した荷重変化ΔＦ１は、ひずみを表す情報に基づくブレード１２への荷重変化ΔＦ２と異なる場合がある。この理由は、設計時に用いる荷重計算式や適宜のツールの計測精度などが不十分な場合も考えられるが、他にもブレードのピッチ制御、発電機出力、制御不具合による想定外の荷重の影響など、多様な場合が考えられる。

【0049】

そこで、図１に示すように、相関関係記憶部３９は、風況情報とブレード１２の疲労損傷量との相関関係の他に、前記所定の部位（タワー１４）におけるひずみと前記した第１の荷重情報との相関関係、風況情報と前記した第２の荷重情報との相関関係、並びに、第１及び第２の荷重情報とブレード１２の疲労損傷量との相関関係、を予め記憶している。一方、第１の荷重情報算出部３４、第２の荷重情報算出部３５及び疲労損傷量補正部３６は、相関関係記憶部３９に記憶された内容を参照しつつ、第１の荷重情報の算出、第２の荷重情報の算出、疲労損傷量の補正、をそれぞれ行う。これにより、ブレード１２に加わる応力を適切に評価することが可能となる。

【0050】

また、疲労損傷量積算部３７は、一定期間の単位で疲労損傷量補正部３６によりそれぞれ補正された複数の疲労損傷量を積算する。ここで、既述してきた一定期間とは、例えば５分、１０分、２０分などといった、所定の時間間隔であって、疲労損傷量積算部３７は、一定期間を５分間として求めた疲労損傷量と、一定期間を１０分間として求めた疲労損傷量などと、を積算するものであってもよい。一方、疲労損傷量出力部３８は、疲労損傷量積算部３７により積算された疲労損傷量を出力する。疲労損傷量出力部３８は、例えば表やグラフなどで疲労損傷量を可視的に表示する表示装置として構成されていてもよいし、風車発電機の疲労損傷量算出装置３０とは別体の表示装置へ疲労損傷量を表示出力するためのインターフェースとして構成されていてもよい。

【0051】

次に、風車発電機の疲労損傷量算出装置３０による疲労損傷量算出方法を図１３に示すフローチャートに基づき説明する。図１３に示すように、まず、風況情報取得部３１は、風車発電機２０の設置場所で（ナセル１５に設けられた風況計測器２３により）、一定期間中に計測された風速の平均及び風速の変化を表す風況情報を取得する（Ｓ１）。次に、疲労損傷量算出部３２は、風況情報取得部３１により取得された風況情報に基づいて、風車発電機２０におけるブレード１２の疲労損傷量を算出する（Ｓ２）。

【0052】

さらに、ひずみ情報取得部３３は、ひずみゲージ２１、２２により一定期間中に計測された風車発電機２０の所定の部位（タワー１４の軸方向）におけるひずみを表す情報を取得する（Ｓ３）。次いで、第１の荷重情報算出部３４は、ひずみ情報取得部３３により取得されたひずみを表す情報に基づいて、一定期間中にブレード１２に加わる荷重の平均及び荷重の変化を表す第１の荷重情報を算出する（Ｓ４）。一方、第２の荷重情報算出部３５は、風況情報取得部３１により取得された風況情報に基づいて、一定期間中にブレード１２に加わる荷重の平均及び荷重の変化を表す第２の荷重情報を算出する（Ｓ５）。

【0053】

続いて、疲労損傷量補正部３６は、第１及び第２の荷重情報算出部３４、３５によりそれぞれ算出された第１及び第２の荷重情報に基づいて、疲労損傷量算出部３２により算出された疲労損傷量を補正する（Ｓ６）。次に、疲労損傷量積算部３７は、補正された疲労損傷量の数（個数）が、規定した数を超えている場合（Ｓ７のＹＥＳ）、前記一定期間の単位で疲労損傷量補正部３６によりそれぞれ補正された複数の疲労損傷量を積算する（Ｓ８）。さらに、疲労損傷量出力部３８は、疲労損傷量積算部３７により積算された疲労損傷量を例えば表示装置に表示出力する（Ｓ９）。

【0054】

既述したように、本実施形態に係る風車発電機の疲労損傷量算出装置３０を含む風力発電システム１０及び風車発電機の疲労損傷量算出方法によれば、センサ類の追加、高機能なサンプリング装置の設置、膨大なデータを分析処理する例えば高性能な演算装置の追加によるモニタリングコストの増加、などを抑えつつ、ブレード１２についての信頼度の高い疲労損傷量を容易に求めることができる。つまり、本実施形態によれば、風車発電機２０における疲労寿命の解析精度を容易に高めることができるので、発電事業者では、解析精度が高められた疲労寿命の解析結果を基に、風車発電機２０に対する適切な点検、補修、更新計画などを立案することができる。

【0055】

＜第２の実施の形態＞
次に、第２の実施形態を図１４に基づき説明する。なお、図１４中において、図１に示した第１の実施形態中の構成要素と同一の構成要素については、同一の符号を付与し重複する説明を省略する。

【0056】

図１４に示すように、第２の実施形態に係る風力発電システム５０は、第１の実施形態の風力発電システム１０が備えていた風車発電機の疲労損傷量算出装置３０に代えて、風車発電機の疲労損傷量算出装置７０を備えている。風車発電機の疲労損傷量算出装置７０は、風車発電機の疲労損傷量算出装置３０の構成に加え、風況情報出力部７１及び荷重情報出力部７２をさらに備えている。

【0057】

風況情報出力部７１は、風況情報取得部３１により取得された風況情報と、風車発電機２０の設計時に想定された風速の平均（平均風速）及び風速の変化（乱流強度）を表す風況情報と、をあわせて出力する。一方、荷重情報出力部７２は、第１及び第２の荷重情報算出部３４、３５によりそれぞれ算出された第１の荷重情報及び第２の荷重情報と、前記第１の荷重情報と前記第２の荷重情報との差異を表す情報と、のうちの少なくともいずれかを出力する。なお、風況情報出力部７１及び荷重情報出力部７２は、例えば表やグラフなどで上記の各風況情報や上記の差異を表す情報を可視的に表示する表示装置として構成されていてもよいし、風車発電機の疲労損傷量算出装置３０とは別体の表示装置へ各風況情報や前記差異を表す情報を表示出力するためのインターフェースとして構成されていてもよい。

【0058】

したがって、本実施形態に係る風車発電機の疲労損傷量算出装置７０を含む風力発電システム５０によれば、設計時に想定した風況に対する実際の風況や、実際の風況による風車発電機の応答性がどのような状況にあるのかを、例えば発電事業者が容易に把握することが可能となる。これにより、発電事業者では、風車発電機の有益な運用を実現することができる。

【0059】

以上、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これらの実施形態やその変形例は、発明の範囲や要旨に含まれると共に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

【符号の説明】

【0060】

１０，５０…風車発電システム、１２…ブレード、１４…タワー、２０…風車発電機、２１，２２…ひずみゲージ、２３…風況計測器、３０，７０…風車発電機の疲労損傷量算出装置、２３…風況計測器、３１…風況情報取得部、３２…疲労損傷量算出部、３３…ひずみ情報取得部、３４…第１の荷重情報算出部、３５…第２の荷重情報算出部、３６…疲労損傷量補正部、３７…疲労損傷量積算部、３８…疲労損傷量出力部、３９…相関関係記憶部、７１…風況情報出力部、７２…荷重情報出力部。

【図1】

【図2】

【図3A】

【図3B】

【図3C】

【図4A】

【図4B】

【図5A】

【図5B】

【図6】

【図7】

【図8A】

【図8B】

【図9A】

【図9B】

【図10】

【図11】

【図12】

【図13】

【図14】