(57)【要約】【課題】  候補地点に風車等を設置してよいか否かを、乱流強度等とは異なる新たな評価基準によって判定することに適した乱流評価装置等を提供する。【解決手段】  乱流評価装置１は、候補地点における自然エネルギーの乱流を評価するためのものである。境界条件記憶部５は、候補地点を含む候補空間における自然エネルギーの境界入力条件を記憶する。シミュレーショ部１１は、境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算する。乱流評価指標演算部１５は、シミュレーション部１１が演算した候補空間における自然エネルギーの変化と、境界入力条件記憶部５が記憶する境界入力条件とを用いて、候補地点における乱流評価指標を演算する。判定部１７は、乱流評価指標が評価基準値を超える場合に、候補地点の自然エネルギーの乱流が大きいと判定する。【選択図】  図１
2018091328.tif 000002
【特許請求の範囲】
【請求項１】
候補地点における自然エネルギーの乱流を評価するための乱流評価装置であって、
前記候補地点を含む候補空間における前記自然エネルギーの境界入力条件を記憶する境界入力条件記憶手段と、
前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算するシミュレーション手段と、
前記シミュレーション手段が演算した前記候補空間における前記自然エネルギーの変化と、前記境界入力条件記憶手段が記憶する前記境界入力条件とを用いて、前記候補地点における乱流評価指標を演算する乱流評価指標演算手段と、
乱流評価指標が評価基準値を超える場合に、前記候補地点の前記自然エネルギーの乱流が大きいと判定する判定手段を備える乱流評価装置。
【請求項２】
前記自然エネルギーの変化は、前記自然エネルギーが変化する速さであり、
前記境界入力条件は、前記候補空間の境界の一部又は全部における前記自然エネルギーの変化の最大の速さを含み、
前記候補空間に含まれる評価点における前記自然エネルギーの変化の速さの標準偏差を演算する標準偏差演算手段を備え、
前記乱数評価指標演算手段は、前記乱流評価指標を、前記評価点における前記標準偏差の、前記最大の速さに対する比として演算する、請求項１記載の乱流評価装置。
【請求項３】
実測して得られた、又は、前記シミュレーション手段により得られた、前記評価点における前記自然エネルギーにおいて、単位時間当たりの変化の平均が判定基準値を超えた方向がある場合に、
前記境界入力条件は、前記評価点において前記判定基準値に対応する前記自然エネルギーの変化を生じ得るものであり、
前記シミュレーション手段は、前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算し、
前記乱流評価指標演算手段は、前記判定基準値を超えた方向における前記自然エネルギーの変化に対して乱流評価指標を演算する、請求項１又は２に記載の乱流評価装置。
【請求項４】
候補地点における自然エネルギーの乱流を評価するための乱流評価装置であって、
前記候補地点を含む候補空間における前記自然エネルギーの境界入力条件を記憶する境界入力条件記憶手段と、
前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算するシミュレーション手段と、
前記シミュレーション手段が演算した前記候補空間における前記自然エネルギーの変化を用いて、前記候補地点における乱流評価指標を演算する乱流評価指標演算手段と、
前記候補地点において前記自然エネルギーを利用して発電する発電装置の設計疲労線を記憶する設計ＤＥＬ記憶手段と、
複数の疲労荷重式を用いて前記発電装置の蓄積疲労を計算する蓄積疲労演算手段と、
前記蓄積疲労が前記設計疲労線を越える時点を演算する寿命評価手段を備え、
前記乱流評価指標演算手段は、複数の方位のそれぞれからの前記自然エネルギーに対して、前記候補地点における前記乱流評価指標を演算し、
前記疲労荷重式は、少なくとも、
前記乱流評価指標が前記評価基準値を超える方位からの前記自然エネルギーによる蓄積疲労を計算するための第１疲労荷重式と、
前記乱流評価指標が前記評価基準値を超えない方位からの前記自然エネルギーによる蓄積疲労を計算するための第２疲労荷重式を含み、
前記第１疲労荷重式と前記第２疲労荷重式とは異なり、
前記蓄積疲労演算手段は、複数の方位のそれぞれに対して、前記乱流評価指標の値に応じて前記複数の疲労荷重式のうちの一つを選択して前記蓄積疲労を演算する、乱流評価装置。
【請求項５】
前記疲労荷重式は、前記発電装置における前記自然エネルギーの変化の速さ及び前記自然エネルギーの変化の速さの標準偏差に対して直線近似することにより得られたものである、請求項４記載の乱流評価装置。
【請求項６】
候補地点における自然エネルギーの乱流を評価するための乱流評価装置における乱流評価方法であって、
前記乱流評価装置は、前記候補地点を含む候補空間における前記自然エネルギーの境界入力条件を記憶する境界入力条件記憶手段を備え、
前記乱流評価装置が備えるシミュレーション手段が、前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算するシミュレーションステップと、
前記乱流評価装置が備える乱流評価指標演算手段が、前記シミュレーション手段が演算した前記候補空間における前記自然エネルギーの変化と、前記境界入力条件記憶手段が記憶する前記境界入力条件とを用いて、前記候補地点における乱流評価指標を演算する乱流評価指標演算ステップと、
前記乱流評価装置が備える判定手段が、乱流評価指標が評価基準値を超える場合に、前記候補地点の前記自然エネルギーの乱流が大きいと判定する判定ステップを含む乱流評価方法。
【請求項７】
コンピュータにおいて、請求項６記載の乱流評価方法を実現するためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本願発明は、乱流評価装置、乱流評価方法及びプログラムに関し、特に、候補地点における自然エネルギーの乱流を評価するための乱流評価装置等に関するものである。
【背景技術】
【０００２】
近年、風車事故が多発している。事故調査の結果、「金属疲労」が事故の主原因と報告されている。
【０００３】
発明者らは、風力発電について研究・開発し、実用化を行ってきた（例えば特許文献１参照）。風車は、一般に、ＩＥＣ規格に基づき設計される（例えば非特許文献１参照）。また、風車の乱流（現地風条件）を評価する値として、乱流強度が知られている。乱流強度は、一般的に、風車の標準偏差の平均風速に対する比として定義されている（非特許文献２の１９５頁参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１６−１８６２５３号公報
【非特許文献】
【０００５】
【非特許文献１】独立行政法人新エネルギー・産業技術総合開発機構，付属書Ａ（参考）「ＩＥＣ規格による風力発電機の設計条件」，[online]， [平成28年11月7日検索]、インターネット＜ＵＲＬ：http://www.nedo.go.jp/content/100107253.pdf＞
【非特許文献２】独立行政法人新エネルギー・産業技術総合開発機構，風力発電導入ガイドブック，[online]，2008年2月，改訂第9版、[平成28年11月7日検索]、インターネット＜ＵＲＬ：http://www.nedo.go.jp/content/100079735.pdf＞
【発明の概要】
【発明が解決しようとする課題】
【０００６】
現時点では、乱流が金属疲労の原因になるとの報告はない。そのため、現在の風車は、ＩＥＣ規格等に従っていれば、いずれの場所に設置しても問題がないとされている。
【０００７】
しかしながら、発明者らは、風車の事故や故障に対して乱流が影響している可能性に着目した。風車がＩＥＣ規格等を満足していても、設置候補地点における乱流の状況によっては、その設置候補地点に風車を設置してはいけない可能性がある。
【０００８】
ここで、乱流強度を基準にして、候補地点に風車を設置してもよいか否かを判定することも考えられる。しかしながら、乱流強度は、例えば風車のハブ高さのように、一定の高さにおける風の標準偏差と平均風速を用いて評価するものである。そのため、候補地点における風条件を十分に評価することができない可能性があった。
【０００９】
そこで、本願発明は、候補地点に風車等を設置してよいか否かを、乱流強度等とは異なる新たな評価基準によって判定することに適した乱流評価装置等を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本願発明の第１の観点は、候補地点における自然エネルギーの乱流を評価するための乱流評価装置であって、前記候補地点を含む候補空間における前記自然エネルギーの境界入力条件を記憶する境界入力条件記憶手段と、前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算するシミュレーション手段と、前記シミュレーション手段が演算した前記候補空間における前記自然エネルギーの変化と、前記境界入力条件記憶手段が記憶する前記境界入力条件とを用いて、前記候補地点における乱流評価指標を演算する乱流評価指標演算手段と、乱流評価指標が評価基準値を超える場合に、前記候補地点の前記自然エネルギーの乱流が大きいと判定する判定手段を備えるものである。
【００１１】
本願発明の第２の観点は、第１の観点の乱流評価装置であって、前記自然エネルギーの変化は、前記自然エネルギーが変化する速さであり、前記境界入力条件は、前記候補空間の境界の一部又は全部における前記自然エネルギーの変化の最大の速さを含み、前記候補空間に含まれる評価点における前記自然エネルギーの変化の速さの標準偏差を演算する標準偏差演算手段を備え、前記乱数評価指標演算手段は、前記乱流評価指標を、前記評価点における前記標準偏差の、前記最大の速さに対する比として演算するものである。
【００１２】
本願発明の第３の観点は、第１又は第２の観点の乱流評価装置であって、実測して得られた、又は、前記シミュレーション手段により得られた、前記評価点における前記自然エネルギーにおいて、単位時間当たりの変化の平均が判定基準値を超えた方向がある場合に、前記境界入力条件は、前記評価点において前記判定基準値に対応する前記自然エネルギーの変化を生じ得るものであり、前記シミュレーション手段は、前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算し、前記乱流評価指標演算手段は、前記判定基準値を超えた方向における前記自然エネルギーの変化に対して乱流評価指標を演算するものである。
【００１３】
本願発明の第４の観点は、候補地点における自然エネルギーの乱流を評価するための乱流評価装置であって、前記候補地点を含む候補空間における前記自然エネルギーの境界入力条件を記憶する境界入力条件記憶手段と、前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算するシミュレーション手段と、前記シミュレーション手段が演算した前記候補空間における前記自然エネルギーの変化と、前記境界入力条件記憶手段が記憶する前記境界入力条件とを用いて、前記候補地点における乱流評価指標を演算する乱流評価指標演算手段と、前記候補地点において前記自然エネルギーを利用して発電する発電装置の設計疲労線を記憶する設計ＤＥＬ記憶手段と、複数の疲労荷重式を用いて前記発電装置の蓄積疲労を計算する蓄積疲労演算手段と、前記蓄積疲労が前記設計疲労線を越える時点を得る寿命評価手段を備え、前記乱流評価指標演算手段は、複数の方位のそれぞれからの前記自然エネルギーに対して、前記候補地点における前記乱流評価指標を演算し、前記疲労荷重式は、少なくとも、前記乱流評価指標が前記評価基準値を超える方位からの前記自然エネルギーによる蓄積疲労を計算するための第１疲労荷重式と、前記乱流評価指標が前記評価基準値を超えない方位からの前記自然エネルギーによる蓄積疲労を計算するための第２疲労荷重式を含み、前記第１疲労荷重式と前記第２疲労荷重式とは異なり、前記蓄積疲労演算手段は、複数の方位のそれぞれに対して、前記乱流評価指標の値に応じて前記複数の疲労荷重式の一つを選択して前記蓄積疲労を演算する。
【００１４】
本願発明の第５の観点は、第４の観点の乱流評価装置であって、前記疲労荷重式は、前記発電装置における前記自然エネルギーの変化の速さ及び前記自然エネルギーの変化の速さの標準偏差に対して直線近似することにより得られたものである。
【００１５】
本願発明の第６の観点は、候補地点における自然エネルギーの乱流を評価するための乱流評価装置における乱流評価方法であって、前記乱流評価装置は、前記候補地点を含む候補空間における前記自然エネルギーの境界入力条件を記憶する境界入力条件記憶手段を備え、前記乱流評価装置が備えるシミュレーション手段が、前記境界入力条件の下で、前記候補空間における前記自然エネルギーの変化を演算するシミュレーションステップと、前記乱流評価装置が備える乱流評価指標演算手段が、前記シミュレーション手段が演算した前記候補空間における前記自然エネルギーの変化と、前記境界入力条件記憶手段が記憶する前記境界入力条件とを用いて、前記候補地点における乱流評価指標を演算する乱流評価指標演算ステップと、前記乱流評価装置が備える判定手段が、乱流評価指標が評価基準値を超える場合に、前記候補地点の前記自然エネルギーの乱流が大きいと判定する判定ステップを含むものである。
【００１６】
本願発明の第７の観点は、コンピュータにおいて、第６の観点の乱流評価方法を実現するためのプログラムである。なお、本願発明を、第７の観点のプログラムを記録するコンピュータ読み取り可能な記録媒体として捉えてもよい。
【発明の効果】
【００１７】
本願発明の各観点によれば、シミュレーションにより得られた自然エネルギーのデータを利用して、乱流評価指標という新たな基準を用いて、候補地点における自然エネルギーの乱流を評価できる。ここで、乱流評価指標は、従来の乱流強度とは異なる。乱流強度は、ハブ高さ等における固定点における風速の標準偏差と平均速度によって得られる値である。他方、乱流評価指標は、境界入力条件換算値であり、シミュレーションにおいて用いられる境界入力条件と、ハブ高さ等における標準偏差等を用いて得られる。そして、後に示すように、発明者らは、この乱流評価指標を使用して、候補地点における乱流が自然エネルギーを利用した発電装置の故障の要因となりうる程度か否かを原因となるか否かの判断できることを実験により明らかにした。
【００１８】
さらに、本願発明の第２の観点によれば、例えば、評価点が風車のハブ高さであり、境界入力条件において風車の頂上の標高における最大風速Ｕin（境界風速）である場合には、「乱流評価指標＝σu（ハブ高さ風速標準偏差）／Ｕin（境界風速）」という計算式により得ることができる。
【００１９】
さらに、本願発明の第３の観点によれば、例えば風速が大きい時間帯などのように、変化が大きい時間帯において乱流評価指標を演算して乱流を評価する。このように、危険性のある方向及びその時間帯を特定して、乱流評価指標により評価することができる。
【００２０】
さらに、本願発明の第４の観点によれば、乱流評価指標に対応した疲労荷重式を用いることにより、発電装置の寿命予測をすることができる。特に、第５の観点によれば、例えば風力発電であれば、第２の観点にあるように乱流指標が風速と風速標準偏差と風速の比から得られること、さらには、発明者らが見出した風速と風速標準偏差が蓄積疲労と直線近似できることを利用して、評価基準値を超えるものと超えないものとに対し、風速と風速標準偏差により直線近似することにより得られた疲労荷重式を利用することにより、精度の高い寿命予測が可能となる。
【図面の簡単な説明】
【００２１】
【図１】本願発明の実施の形態の一例に係る乱流評価装置の構成の一例を示すブロック図である。
【図２】図１の乱流評価装置１の動作の一例を示すフロー図である。
【図３】本願で新たに導入する乱流評価指標を説明するための図である。
【図４】（ａ）実際の風車での歪みセンサを設置した状況と、（ｂ）〜（ｄ）歪み計測データを示す図である。
【図５】東風のときの風況及び歪み計測データの分析と、風車疲労荷重評価を示す図である。
【図６】風況シミュレーションにより得られたデータを示す図である。
【図７】本願発明の実施の形態の他の一例に係る乱流評価装置の構成の一例を示すブロック図である。
【図８】図７の乱流評価装置２１の動作の一例を示すフロー図である。
【図９】東方位と北方位を対象とした風速と風速標準偏差と疲労（ＤＥＬ）の分析結果を示す。
【図１０】過去の実測データを用いて行ったブレードの寿命評価を示す。
【図１１】過去の実測データに加えてシミュレーションにより得られた乱流評価指標を用いて行ったブレードの寿命評価を示す。
【発明を実施するための形態】
【００２２】
以下、図面を参照して本願発明の実施例について説明する。なお、本願発明は、以下の実施例に限定されるものではない。
【実施例１】
【００２３】
図１は、本願発明の実施の形態の一例に係る乱流評価装置の構成の一例を示すブロック図である。図２は、図１の乱流評価装置１の動作の一例を示すフロー図である。図１及び図２を参照して、乱流評価装置１の構成及び動作の一例を説明する。
【００２４】
ここで、乱流評価装置１は、自然エネルギーとして風を使って発電する風力発電装置（風車）を候補地点に設置してよいか否かを判定するものであるとする。また、候補空間は、風車を設置する空間であるとする。候補空間の最も高い場所は、風車の最も高い箇所（通常、ブレード先端高さである。）の標高となる。評価点は、風車のハブ高さであるとする。
【００２５】
従来知られている乱流強度は、「乱流強度＝σu（ハブ高さ風速標準偏差）／Ｕ（ハブ高さ平均風速）」と表現されるように、風車によって定まる一定の高さにおける風のデータのみを使用する。これは、風車で実測できるため、広く利用されている。しかしながら、乱流強度は、いわば、ハブ高さにおける風という一部のサンプルによって候補地点の評価をしようとしている。そのため、風力発電の事業者は、乱流強度によっては、風車を、全体として評価することはできない。
【００２６】
他方、発明者らは、自ら研究・開発した高精度なシミュレーション装置によって、任意の条件の下で、風のデータを得ることができる。そして、新たに導入する乱流評価指標は、「乱流評価指標＝σu（ハブ高さ風速標準偏差）／Ｕin（境界風速）」であり、シミュレーションにおいて用いる境界入力条件を参照して、境界風速換算値として得られるものである。そのため、乱流評価装置１は、境界入力条件の下でのシミュレーションにより得られた風のデータを利用して、境界入力条件を加味した乱流評価指標により、候補地点における乱流の評価をすることができる。よって、風力発電の事業者は、シミュレーションを利用して、候補地点に風車を設置した状態を想定して、候補地点に風車を設置してよいか否かを判断することができる。
【００２７】
このように、本願により新たに導入する乱流評価指標は、候補地点に風車を設置する空間を評価することができる。さらに、各風向によっても高精度に評価できるため、風力発電の事業者にとって、例えば、故障の危険性がある風向を事前に知り、その風向のときには発電を止めて故障を回避し、他の風向のときは発電することなどの適切な発電制御を行うことにより、故障による停止を回避して、効率よく電気を得ることができる。このように、候補地点の評価は、故障を予防する対応策が可能か否かによっても評価が変わる。本願発明によれば、方向に応じた制御などの使用態様も考慮することができる。
【００２８】
図１を参照して、乱流評価装置１は、地形データ記憶部３（本願請求項の「地形データ記憶手段」の一例）と、境界入力条件記憶部５（本願請求項の「境界入力条件記憶手段」の一例）と、基準値記憶部７と、乱流評価指標記憶部９と、シミュレーション部１１（本願請求項の「シミュレーション手段」の一例）と、制御部１２と、標準偏差演算部１３（本願請求項の「標準偏差演算手段」の一例）と、乱流評価指標演算部１５（本願請求項の「乱流評価指標演算手段」の一例）と、判定部１７（本願請求項の「判定手段」の一例）と、表示部１９（本願請求項の「表示手段」の一例）を備える。
【００２９】
地形データ記憶部３は、候補地点を含む地域の地形を特定するための地形データを記憶する。境界入力条件記憶部５は、境界入力条件を記憶する。本実施例では、境界入力条件には、境界風速が含まれるとする。ここで、境界風速は、風車のブレード先端高さの標高での最大風速であるとする。以下では、境界風速は１０ｍ／ｓとする。基準値記憶部７は、判断処理の基準値を記憶する。本実施例では、判定基準値と評価基準値を記憶する。乱流評価指標記憶部９は、乱流評価指標演算部１５が演算した乱流評価指標を記憶する。
【００３０】
図２を参照して、図１の乱流評価装置１の各構成の動作の一例を説明する。
【００３１】
まず、シミュレーション部１１は、境界入力条件記憶部５が記憶する境界入力条件の下で、地形データ記憶部３が記憶する地形データを参照して、評価期間での評価点（ハブ高さ）における風のデータを得る（ステップＳＴ１）。
【００３２】
次に、制御部１２は、例えば４５°で分割して、東（東を中心に４５°の範囲から吹く風。以下同様。）、南東、南、南西、西、北西、北、北東の８つに分割するように、風向を複数に分割し、その一つを選択する（ステップＳＴ２）。なお、上下を加えて３次元的に分割してもよい。
【００３３】
制御部１２は、シミュレーション部１１が得た風のデータにおける風向を利用して、ステップＳＴ２で設定された風向の風の風速の経時的データを得る（ステップＳＴ３）。ステップＳＴ３で得た風速の経時的データは、一つ又は複数の期間のものとする。連続する期間に属する風のデータを、「風速群」という。期間を分ける基準は、例えば、他の風向から設定された風向になり他の風向に変わるまでを一つの風速群としたり、同じ風向であっても風速が所定値未満から所定値以上となって所定値未満になるまでを一つの風速群としたり、それらを組み合わせて風速群としたりしてよい。
【００３４】
制御部１２は、一つの風速群を選択する（ステップＳＴ４）。そして、選択した風速群での単位時間の平均風速が、判定基準値を超えるか否かを判断する（ステップＳＴ５）。本実施例では、単位時間は１０分間であり、判定基準値は９．０ｍ／ｓである。単位時間の平均風速が判定基準値を超えるならば、ステップＳＴ６に進む。単位時間の平均風速が判定基準値を超えないならば、ステップＳＴ１０に進む。
【００３５】
ステップＳＴ６において、標準偏差演算部１３は、ハブ高さにおける風速の標準偏差を演算する。そして、乱流評価指標演算部１５は、標準偏差演算部１３が演算した風速の標準偏差と、境界入力条件記憶部５が記憶する境界速度を用いて、候補地点における乱流評価指標を演算する（ステップＳＴ７）。
【００３６】
判定部１７は、乱流評価指標が評価基準値を超えるか否かを判定する（ステップＳＴ９）。ここで、評価基準値は、０．２である。判定部１７は、乱流評価指標が評価基準値を超えないならば、ステップＳＴ１０に進む。判定部１７は、乱流評価指標が評価基準値を超えるならば、候補地点における風の乱流が大きいと判定し、表示部１９に対して、風向と、乱流評価指標が評価基準値を超えることを表示させ（ステップＳＴ９）、ステップＳＴ１０に進む。
【００３７】
制御部１２は、全ての風速群に対して処理を行ったか否かを判断する（ステップＳＴ１０）。処理を行っていない風速群があるならば、ステップＳＴ４に戻り、処理を行っていない一つの風速群を選択する。全ての風速群に対して処理を行ったならば、全ての風向に対して処理を行ったか否かを判断する（ステップＳＴ１１）。全ての風向に対して処理を行っていないならば、ステップＳＴ２に戻り、処理を行っていない一つの風向を選択する。全ての風向に対して処理を行ったならば、図２の処理を終了する。
【００３８】
図３〜図６を参照して、発明者らが行った実験を説明する。
【００３９】
図３は、本願で新たに導入する乱流評価指標を説明するための図である。
【００４０】
図３（ａ）は、風車の鉛直方向の理想的な層流の風速分布を示す。地表から滑らかに増加し、羽根（ブレード）の部分では均等に分布することにより、各羽根には均等に力が働き、回転スピードはほぼ同一となる。そのため、３枚の羽根の根元にかかるストレスもほぼ同じとなる。風車性能が規定される風条件は、図３（ａ）のような状態である。
【００４１】
しかしながら、例えば風上に山があるような地点では、図３（ｂ）にあるように、地形の変化により、地表近くの風速や風の周期が乱される。図３（ｃ）にあるように、風の乱れは、図３（ａ）の理想的な風速分布から離れることとなる。図３（ｃ）では、受風面の下部の風速が、乱流化の影響で減速されると、３枚の羽根の回転スピードや、根元にかかるストレスが異なることになる。特に、乱流化の影響で風の周期も小さくなるために、羽根の根元には、サイクリックにストレスが発生し、そのストレスが蓄積し、最終的には危機損傷につながる可能性がある。そして、これによって金属疲労等が進行し、故障や事故が発生する可能性が高まる。図３（ｄ）は、後のシミュレーションで得られたデータによる乱流化の状況を示すものである。縦軸は地表からの高さを示し、横軸は速度分布を示す。矢印は、ハブ高さにおいて、理想的な風速分布から減速した幅を示す。
【００４２】
ここで、複雑地形特有の風条件で金属疲労が生じても、風車の構造強度で十分に耐えることができるのであれば、その地形に風車を設置することができる。そこで、実験では、実際に測定した風車機器に入力される風速のデータや、風車ブレード歪みのデータを測定し、複雑地形特有の風条件と風車の構造強度への影響を特定したうえで、複雑地形での風条件の再現結果（リアムコンパクト社による風況解析結果）に基づく、数値風況面から、風車事故を提言する最適配置基準の定量化を行う。この最適配置基準を利用することにより、最適配置計画を支援することができる。
【００４３】
実験では、まず、実際の風車において、風況と風車挙動の同期計測データを比較し、歪みデータ変動量（振幅）の大きい風向を特定した。ここで、風況は、ナセル風速計・風向計を使用して風速・風向変動を解析した。風車挙動は、３枚のブレード根元ルート部に歪みケージを設置し、歪み変動を解析した。
【００４４】
そして、特定した風向を対象に、ブレード曲げ疲労荷重評価を実施して、複雑地形での風車疲労への影響を把握した。歪みデータより２０年間のブレード曲げ疲労荷重を算定し、ブレード曲げ疲労ダメージを評価した（ブレード曲げ疲労ダメージ＝歪みデータより算定した疲労荷重（２０年）／設計許容疲労荷重（２０年））。ブレード曲げ疲労ダメージが１を超えた場合、耐用年数２０年経過前に疲労破壊に至ることとなる。
【００４５】
図４（ａ）は、実際の風車での歪みセンサを設置した状況を示す。この風車は、ブレード先端高さ１００ｍであり、ブレードの直径が８０ｍ、ハブ高さは６０ｍである。ブレード３本の根元（ルート部：ハブ接合面から約１．３ｍ）は、既設設備に後付けで取り付け可能な場所であり、歪みセンサを設置した。ブレード根元の強度評価（ブレードの曲げ荷重評価）は、ブレードを通して風力エネルギーが入力されるため、重要である。
【００４６】
図４（ｂ）〜（ｄ）は、歪み計測データを示す。（ｂ）東風では、11/13の8:00〜10:00に、大きな歪みが計測されている。（ｃ）北風では、11/9の17:00〜22:00に、大きな歪みが計測されている。（ｄ）南風では、11/14の2:00〜5:00に、大きな歪みが計測されている。
【００４７】
風向別に風速及び歪み計測データを分析すると、東風、北風及び南風で、最大９／ｓ程度まで風速が上昇した際に、歪み変動（変動）が最大になっている。東風、北風及び南風は、大きな歪みが出た時間帯で、ほぼ同程度の風速まで上昇したが、歪み変動は、東風が最も大きかった。疲労ダメージは、東風が2.03、北風は0.95、南風が0.93であった。東風は１を超えるため、耐用年数２０年経過前に疲労破壊に至ることが予想される。
【００４８】
図５を参照して、東風のときの風況及び歪み計測データの分析と、風車疲労荷重評価を具体的に説明する。図５（ａ）を参照して、歪み変動が大きな11/13の8:00〜13:00を具体的に分析する。ここで、サンプリング周期は、５０Ｈｚ（０．０２ｓ）である。この時間帯の風速及び風向の時系列データを分析すると、連続的に東方向の風が発生し、乱流化の影響により風の周期も短くなっている。図５（ｂ）を参照して、風速９．７ｍ／ｓ（１０分間平均）のとき、疲労ダメージの最大が２．０３になっている。風速の標準偏差は、約２．４ｍ／ｓである。東風が連続して発生した8:00〜10:00にかけて疲労ダメージが２程度となっている。その後、南東に移行し１未満となっている。そのため、例えば東風のときには発電を停止し、南東の風に移行したときに発電を再開するなどように、風向に応じた制御を行うことにより、故障を回避しつつ、継続した発電を実現することができる。
【００４９】
続いて、図６を参照して、風況シミュレーションにより得られたデータについて説明する。ここで、地形による乱流の影響を評価するために、東風のときの風上の山があるとき（標高修正なし）と、山がないとき（標高修正あり）を比較した。
【００５０】
図６（ａ）は、東風について、標高修正のある場合と無い場合を比較したものである。標高修正が無い場合は、標高修正がある場合に比較して、風速の変化が大きかった。標高修正がある場合の風速の変化は、風上に山がない北風や南風での風速の変化と同様であった。
【００５１】
実測風況データと風況シミュレーション結果を比較すると、東風（風向９０deg付近）の場合、風向変動周期は、両者が同一の周期的変動（周期６〜７ｓ程度）が見られ、風速変動が大きかった。また、風速の標準偏差は、両者ともに約２．４ｍ／ｓ（シミュレーションではＵ成分）であり、ほぼ一致する傾向が見られた。そのため、実機では、シミュレーション結果で得られたように、乱流が発生したと考えられる。
【００５２】
図３（ｄ）で示したように、境界入力条件として与えた風速に対して、実機では、ブレード面で減速（乱流化）している。シミュレーション結果では、地形の影響を受けて、受風面下部に向けて平均風速が小さくなっており乱流が発生していると考えられる。
【００５３】
乱流評価について、実測の場合、東風風速９．７ｍ／ｓ、標準偏差２．４で、乱流強度２５％となった。東風において疲労ダメージ２．０３で、５．８８年で設計疲労荷重の計画値に到達する結果となった。
【００５４】
シミュレーション結果では、東風において、本願が新たに提案する乱流評価指標（境界風速（１０ｍ／ｓ）換算値）は、0.236（＝ハブ高さ標準偏差2.36m/s／境界風速10m/s）となった。なお、乱流強度は、64%（ハブ高さ標準偏差2.36m/s／平均風速3.71m/s）であった。他方、北風及び南風では、それぞれ、乱流評価指標は0.099及び0.111、乱流強度は0.11及び0.09である。これより、乱流評価指標が０．２を超えないという配置基準（すなわち、乱流評価指標の上限を０．２とする基準）は妥当であり、今回の実験で、乱流評価指標0.236で、疲労ダメージ2.03となり、定量的に評価することができた。
【００５５】
以上をまとめると、風車事故を低減できる配置は、以下のように考えることができる。まず、例えば風況観測データをもとに風速が９ｍ／ｓを超えた場合において、その風向に乱流を引き起こす障害（地形・地物）の有無を確認する。次に、障害が引き起こす風向をシミュレーションし、乱流評価指標が0.2以内であることを確認する。すなわち、今回の実験で、乱流評価指標が0.2を超えたとき、疲労ダメージが１以上となった。
【実施例２】
【００５６】
図７は、本願発明の実施の形態の他の一例に係る乱流評価装置の構成の一例を示すブロック図である。図８は、図７の乱流評価装置２１の動作の一例を示すフロー図である。図７及び図８を参照して、乱流評価装置２１の構成及び動作の一例を説明する。
【００５７】
乱流評価装置２１は、自然エネルギーとして風を使って発電する風力発電装置（風車）を候補地点に設置してよいか否かを判定するものであるとする。
【００５８】
図７を参照して、乱流評価装置２１は、地形データ記憶部３（本願請求項の「地形データ記憶手段」の一例）と、境界入力条件記憶部５（本願請求項の「境界入力条件記憶手段」の一例）と、基準値記憶部７と、乱流評価指標記憶部９と、シミュレーション部１１（本願請求項の「シミュレーション手段」の一例）と、制御部１２と、標準偏差演算部１３（本願請求項の「標準偏差演算手段」の一例）と、乱流評価指標演算部１５（本願請求項の「乱流評価指標演算手段」の一例）を備える。これらの各部は、図１の乱流評価装置１における同じ符号を付したものと同様に動作し、乱流評価指標演算部１５は、各方向の乱流評価指標を演算する。
【００５９】
乱流評価装置２１は、さらに、実測値記憶部２３と、設計ＤＥＬ記憶部２５（本願請求項の「設計ＤＥＬ記憶手段」の一例）と、疲労荷重式演算部２７と、蓄積疲労演算部２９（本願請求項の「蓄積疲労演算手段」の一例）と、寿命評価部３１（本願請求項の「寿命評価手段」の一例）を備える。
【００６０】
図８を参照して、シミュレーション部１１は、境界入力条件記憶部５が記憶する境界入力条件の下で、地形データ記憶部３が記憶する地形データを参照して、評価期間での評価点（ハブ高さ）における風のデータを得る（ステップＳＴＪ１）。
【００６１】
次に、乱流評価指標演算部１５は、１６方位に分割された各方位に対して、標準偏差演算部１３が演算した風速の標準偏差と、境界入力条件記憶部５が記憶する境界速度を用いて、候補地点における乱流評価指標を演算する（ステップＳＴＪ２）。
【００６２】
ステップＳＴＪ１及びＳＴＪ２は、具体的には、図１においてステップＳＴ８及びＳＴ９を省略したものと同様にして演算することができる。
【００６３】
疲労荷重式演算部２７は、各方位の一部又は全部に対して、直線近似を行い、疲労荷重式を演算する（ステップＳＴＪ３）。
【００６４】
蓄積疲労演算部２９は、各方位の乱流評価指標の値に応じて、疲労荷重式を選択し（ステップＳＴＪ４）、実測値記憶部２３が記憶する一定期間の実測データを用いて蓄積疲労を演算する（ステップＳＴＪ５）。なお、実測データがない場合には、シミュレーション等により得られた値により代用してもよい。
【００６５】
疲労荷重式演算部２７が、全部の方位に対して疲労荷重式を演算するならば、蓄積疲労演算部２９は、各方位に対応する疲労荷重式を選択して蓄積疲労を演算すればよい。
【００６６】
疲労荷重式演算部２７が一部の方位に対して疲労荷重式を演算する場合には、乱流評価指標が評価基準値を超える方位の少なくとも一つに対して疲労荷重式を演算し、蓄積疲労演算部２９は、乱流評価指標が評価基準値を超える方位に対しては、乱流評価指標が評価基準値を超える方位から得られた疲労荷重式の一つを選択する。疲労荷重式演算部２７は、乱流評価指標が評価基準値を超えない方位の少なくとも一つに対して疲労荷重式を演算し、蓄積疲労演算部２９は、乱流評価指標が評価基準値を超えない方位に対しては、乱流評価指標が評価基準値を超えない方位から得られた疲労荷重式の一つを選択する。
【００６７】
図９〜図１１では、疲労荷重式演算部２７は、乱流評価指標が評価基準値０．２を超えるもの、評価基準値０．２を超えず０．１を超えるもの、０．１を超えないものから、それぞれ一つの方位を選択し、疲労荷重式を演算し、蓄積疲労演算部２９は、乱流評価指標が評価基準値０．２を超える方位、評価基準値０．２を超えず０．１を超える方位、０．１を超えない方位に対して、それぞれ、対応する疲労荷重式を選択して蓄積疲労を計算している。
【００６８】
寿命評価部３１は、蓄積疲労演算部２９が演算した蓄積疲労が、時の経過に応じて増加して、設計ＤＥＬ記憶部２５が記憶する風車の設計疲労線を超える時期を求めて、風車の寿命を評価する（ステップＳＴＪ６）。
【００６９】
図９〜図１１を参照して、具体的な例について説明する。
【００７０】
図９は、東方位（乱流評価指標０．２３６）と北方位（乱流評価指標０．０９９）を対象とした分析結果を示す。東方位及び北方位を対象に、風速、風速標準偏差、乱流強度及び疲労（ＤＥＬ）の実測データを用いて相互にデータの関係性を分析した結果、以下の２点が判明した。すなわち、疲労（ＤＥＬ）と風速の関係は、直線近似が可能である。疲労（ＤＥＬ）と風速標準偏差（σ）の関係は、直線近似が可能である。
【００７１】
図９（ａ）及び（ｂ）は、それぞれ、東方位及び北方位の１０分の風速（ｍ／ｓ）と風速標準偏差（ｍ／ｓ）の関係を示すグラフである。横軸は風速（ｍ／ｓ）であり、縦軸は風速標準偏差（ｍ／ｓ）である。直線近似は、東方位が傾きａ＝０．２３であり、北方位が傾きａ＝０．１４である。東方位の場合、北方位に比べて風速標準偏差の傾きが大きく、北方位と比較して風速が大きくなるほど風速標準偏差が大きくなる。
【００７２】
図９（ｃ）は、風速と疲労との関係を示す。横軸は風速（ｍ／ｓ）であり、縦軸は疲労（ＤＥＬ）である。直線近似は、東方位が傾きａ＝０．１５であり、北方位が傾きａ＝０．０８である。図９（ｄ）は、風速標準偏差と疲労との関係を示す。横軸は風速標準偏差（ｍ／ｓ）であり、縦軸は疲労（ＤＥＬ）である。直線近似は、東方位が傾きａ＝０．４１であり、北方位が傾きａ＝０．４６である。よって、風速標準偏差（σ）が大きいほど、風速と疲労（ＤＥＬ）の傾きも大きくなることが判明した。
【００７３】
図１０は、過去の風速・風向の実測データ（１年分）を用いて20年間のストレス（蓄積疲労）を積算し、耐用年数である２０年の設計疲労線と比較してブレードの寿命評価を行ったものである。
【００７４】
図１０（ａ）は、１２方位のそれぞれの方向角を示す。図１０（ｂ）は、各方位の実測疲労（ＤＥＬ）線を示すグラフである。横軸は風速（ｍ／ｓ）であり、縦軸は疲労（ＤＥＬ）である。図１０（ｃ）は、実測データにおいて、各方位の風向で、風速に応じた時間を積算することにより得られた、実測ＤＥＬと設計ＤＥＬを示す。１年間の実測ＤＥＬ及び設計ＤＥＬの蓄積疲労は、それぞれ、13,976及び20,425である。20年間の蓄積疲労は、それぞれ、279,524及び408,496である。ストレス（蓄積疲労）割合を、実測ＤＥＬ／設計ＤＥＬとする。寿命（年）は、２０（年、耐用年数）／ストレス割合であり、２０／０．６８＝２９（小数点以下の第一位を四捨五入）となる。
【００７５】
図１１は、乱流評価指標を用いて蓄積疲労を計算した場合を示す。実測データが１２方位であるのに対し、シミュレーションが１６方位で行われる。
【００７６】
この例では、乱流評価指標が０．２を超える方位として１２方位における方位４（東）を、０．２を超えず０．１を超える方位として１２方位における方位１２を、０．１を超えない方位として１２方位における方位１（北）を選択し、直線近似により疲労荷重式を計算している。方位４、１２、１に対応する疲労荷重式は、それぞれ、実測データより得て、傾きが０．１５２７、０．０９４５、０．０７８２であり、切片が−０．４０８９、−０．２３０４、−０．１４５８である。図１１（ａ）は、これらの疲労荷重式と、設計ＤＥＬ線を示す。横軸は風速（ｍ／ｓ）であり、縦軸は疲労（ＤＥＬ）である。
【００７７】
図１１（ｂ）は、各方位の乱流評価指標を示す。０．２を超える１６方位における方位５は、１２方位における方位４に対して得られた疲労荷重式を用いる。０．２を超えず０．１を超える１６方位における方位２、４、６、７、８、９，１０、１２、１３、１４、１５、１６は、１２方位における方位１２に対して得られた疲労荷重式を用いる。０．１を超えない１６方位における方位１、３、１１は、１２方位における方位１に対して得られた疲労荷重式を用いる。図１１（ｃ）は、実測データにおいて、各方位の風向で、風速に応じた時間を積算することにより得られた、実測ＤＥＬと設計ＤＥＬを示す。１年間の実測ＤＥＬ及び設計ＤＥＬの蓄積疲労は、それぞれ、13,184及び20,425である。20年間の蓄積疲労は、それぞれ、263,684及び408,496である。ストレス（蓄積疲労）割合を、実測ＤＥＬ／設計ＤＥＬとする。寿命（年）は、２０（年、耐用年数）／ストレス割合であり、２０／０．６５＝３１（小数点以下の第一位を四捨五入）となる。
【００７８】
図１０及び図１１を比較して、乱流評価指標を用いて寿命予測した結果、寿命は３１年であった。１２方位の実測疲労（ＤＥＬ）線を用いて寿命予測した結果（寿命２９年）とおおむね一致した。
【００７９】
以上より、風速と風速標準偏差が、風車ストレス（蓄積疲労）に密接に関連していることが裏付けられた。さらに、風況シミュレーション、風況データ及び運転データの分析を行うことにより、風車ストレス（蓄積疲労）が大きい風向を特定することができる。さらに、実測風速と風況シミュレーション結果で得られた乱流評価指標に対応した疲労荷重式により、ブレード寿命予測が可能となる。
【００８０】
この実施例を踏まえ、風力発電装置に対する乱流の影響を分析する手法の一例を説明する。風力発電装置では、通常、風況データを取得している。そのため、過去の風速、風向、風速標準偏差を得ることができる。また、過去の運転データとして、ＡＮＮ、発電出力、加速度などを得ることができる。これらを利用して、風況・運転データを分析することができる。本実施例によれば、これに加えて、風況シミュレーションを実施し、複数の方位ごとに乱流評価指標を得ることができるとともに、過去の風況データを利用して、乱流評価指標に基づき疲労荷重式を選択して、風車の寿命を予測することができる。そして、風況・運転データ分析及び風車の寿命により乱流の影響を判定し、現状の保守を継続すべきか、新たな保守・運転を提案するかを判断することができる。
【符号の説明】
【００８１】
１，２１  乱流評価装置、３  地形データ記憶部、５  境界入力条件記憶部、７  基準値記憶部、９  乱流評価指標記憶部、１１  シミュレーション部、１２  制御部、１３  標準偏差演算部、１５  乱数評価指標演算部、１７  判定部、１９  表示部、２３  実測値記憶部、２５  設計ＤＥＬ記憶部、２７  疲労荷重式演算部、２９  蓄積疲労演算部、３１  寿命評価部
【図１】
2018091328.tif 000003
【図２】
2018091328.tif 000004
【図３】
2018091328.tif 000005
【図４】
2018091328.tif 000006
【図５】
2018091328.tif 000007
【図６】
2018091328.tif 000008
【図７】
2018091328.tif 000009
【図８】
2018091328.tif 000010
【図９】
2018091328.tif 000011
【図１０】
2018091328.tif 000012
【図１１】
2018091328.tif 000013