› 風速プロファイル IEC 61400 シミュレーター 戻る
風力発電・大気境界層

風速プロファイル IEC 61400 シミュレーター

大気境界層内の風速プロファイル（べき乗則・対数則）と IEC 61400 風車クラスを組み合わせ、ハブ高風速・ロータ上下の鉛直シア・等価風速・発電出力をリアルタイム計算します。地表面粗度や大気安定度を変えて、立地に応じた風車の最適設計を検討できます。

パラメータ設定
IEC 風車クラス（地形分類）
基準風速 V_ref と年平均 V_avg を自動設定
表面粗度 z₀
対数則で使用。地表の凹凸スケール
参照高 z_ref
m
参照風速 V_ref
m/s
ハブ高 z_hub
m
ロータ径 D
m
べき指数 α
中立 0.143、安定 0.30、不安定 0.10 が目安
大気安定度
α の推奨値を自動表示
計算結果
—
べき乗則ハブ風速 (m/s)
—
対数則ハブ風速 (m/s)
—
ロータ上下風速差 (m/s)
—
等価風速 V_eq (m/s)
—
ハブ高発電出力 (kW)
—
等価発電損失 (%)
—
風速プロファイルとロータ — 模式図

右側の曲線が風速プロファイル V(z)。タワー上のロータ円とハブ位置、上下風速ベクトルが表示されます。色は鉛直シアの強さを表します。

風速プロファイル V(z) — べき乗則 vs 対数則
α 値別ハブ高風速比較
理論・主要公式

$$V(z) = V_{\mathrm{ref}}\left(\frac{z}{z_{\mathrm{ref}}}\right)^{\alpha} = V_{\mathrm{ref}}\,\frac{\ln(z/z_{0})}{\ln(z_{\mathrm{ref}}/z_{0})}$$

べき乗則と対数則による高度 z での風速。α：べき指数（中立 0.143）、z₀：表面粗度長。

$$V_{\mathrm{eq}} = \left(\frac{V_{\mathrm{top}}^{3}+V_{\mathrm{hub}}^{3}+V_{\mathrm{bot}}^{3}}{3}\right)^{1/3}$$

等価風速 V_eq（REWS）。ロータ上端・ハブ・下端の 3 点を 3 乗平均（cube-weighted）で平均化したエネルギー等価の風速。

$$P = \tfrac{1}{2}\,\rho\,A\,V^{3}\,C_p, \qquad A=\tfrac{\pi}{4}D^{2}$$

発電出力 P（kW）。ρ：空気密度 1.225 kg/m³、A：ロータ掃過面積、Cp：パワー係数 0.45（典型値）。

風速プロファイル べき乗則・対数則 — IEC 61400 風力タービン

🙋
風って高いところほど強くなるって聞きますけど、どれくらい違うものなんですか？
🎓
地表近くは地面の摩擦で風が遅くなる「大気境界層」になっていて、上にいくほど自由大気の速度に近づいていくんだ。デフォルトの設定（草地 z₀=0.03、参照高 10m で 8 m/s）でハブ高 100m に外挿すると、べき乗則で約 11.1 m/s、対数則で約 11.2 m/s。地表 10m の値と比べて 40% も増える。風車の出力は風速の 3 乗で効くから、これだけ違うとざっと 2.7 倍のエネルギーが取れる計算になるんだよ。
🙋
だからタワーをどんどん高くしてるんですね。じゃあ最近の Haliade-X みたいに 220m まで上げれば、もっと取れるんですか？
🎓
そう。でも代わりに「鉛直シア」という別の問題が出てくる。ロータ径が 150m あるとブレードの上端は z=175m、下端は z=25m、その風速差はデフォルト設定でも約 2.89 m/s もある。これを 1 回転 ⼀周ごとに翼が受けるから、ブレード根元や主軸に周期的な疲労荷重が乗る。シア/V_hub が 0.2 を超えると IEC 規格上も疲労評価が厳しくなるよ。本ツールはこの比を自動チェックして verdict 表示する。
🙋
出力の見積もりも、ハブ高の風速で計算してはダメなんですか？
🎓
それが「等価風速（REWS）」の話だね。風車が実際に受けるパワーはロータ面平均の V³ で決まる。デフォルト条件だと V_top=12.0、V_hub=11.1、V_bot=9.1 で、3 乗平均すると V_eq=10.86 m/s。ハブ高でそのまま計算すると 6691 kW、等価で計算すると 6262 kW、差が 6.4% も出る。これを無視すると年間収益を過大評価することになるんだ。最近の風資源評価では LiDAR でプロファイルを直接測って REWS を使うのが標準だよ。
🙋
大気安定度って、α にどれくらい影響するんですか？
🎓
かなり大きい。安定大気（夜、地表が冷えて空気が混ざらない）だと α は 0.3〜0.5 まで跳ね上がる。逆に晴天の昼間で対流が起きてる不安定大気だと 0.10 くらい。同じ立地でも夜と昼で全然違うプロファイルになるから、年間の発電量予測には Monin–Obukhov 長を使った安定度補正をかけるのが本格的なやり方だ。ツールでは安定度セレクタを切り替えると、α の推奨値が verdict に表示されるよ。
🙋
沿岸（IEC I）と内陸（II/III）でクラスを分けてるのは、こういう違いを取り入れるためなんですね。
🎓
そう、IEC 61400-1 は基準風速 V_ref（50年再現期間の 10 分平均最大風速）で I=50、II=42.5、III=37.5、IV=30 m/s と区切っていて、年平均は V_avg = 0.2·V_ref。沿岸のような風強・低粗度サイトには I、低風速の内陸丘陵には III/IV を使う。逆に Class III の機体を沿岸に持ってくると、定格を超える強風で疲労寿命を消費しすぎる。立地→クラス選定→ハブ高・ロータ径という順番で詰めるのが王道だよ。

よくある質問

実務ではべき乗則 V(z)=V_ref·(z/z_ref)^α が広く使われます。指数 α 一つで上空への外挿ができ、IEC 61400 でも中立で α=0.143（陸上）、α=0.20（簡易設計）などの値が規定されています。一方、対数則 V(z)=V_ref·ln(z/z0)/ln(z_ref/z0) は表面粗度長 z0 が明確なときに物理的妥当性が高く、空港や港湾の風観測でよく用いられます。本ツールは両方を同時に計算し、ハブ高での差を比較できます。
ロータ上下で風速差 ΔV があると、1 回転中に翼が受ける空力荷重が周期的に変動し、ヨー軸・主軸・ブレード根元に疲労荷重を与えます。IEC 61400-1 では正規シアモデル（NWP, EWS）で α=0.20 を疲労評価の基準にしています。本ツールでは ΔV/V_hub > 0.20 を NG（疲労過大）、>0.15 を警告として表示します。直径 150m・ハブ高 100m の機体だと α=0.143 でも上下差は約 3 m/s 程度に達し、疲労寿命の主要因の一つとなります。
Rotor-equivalent wind speed（REWS）は、ロータ面内の風速分布から取得される総運動エネルギーを「ロータ面平均の風速」に換算した値で、V_eq = ((V_top³+V_hub³+V_bot³)/3)^(1/3) のように 3 乗平均（cube weighted）で求めます。風車出力は V³ に比例するため、シアが強い場合はハブ高風速で見積もると数 % 過大評価になります。本ツールではハブ高出力 P_hub と等価出力 P_eq の差を「等価発電損失（%）」として表示します。
中立大気（風が強く混合が活発）で α≈0.14、安定大気（夜間・地表が冷える）で α≈0.30〜0.50、不安定大気（晴天日中・対流活発）で α≈0.10 程度になります。風車の出力評価で年平均を用いるなら α≈0.14〜0.20、夜間ピーク負荷の疲労評価では安定大気の高 α を採用します。本ツールの「大気安定度」セレクタは α の推奨値を表示するので、現場観測がないときの初期値として活用してください。

実世界での応用

洋上・沿岸ウィンドファーム：北海・台湾海峡などの洋上サイトは表面粗度 z₀≈0.0002（海面）と極小で α≈0.10、ハブ高 100〜150m でハブ高風速が 10〜12 m/s に達します。Vestas V236 や Siemens Gamesa SG 14-222 のような大型機が高効率で稼働でき、年平均稼働率（capacity factor）が陸上の倍近い 45〜55% に達します。一方、波浪・塩害・洗掘などの固有の課題があり、設計には IEC 61400-3（洋上）が併用されます。

陸上低風速サイトの最適化：内陸の低風速地帯（年平均 5〜6 m/s、IEC III/IV）では、より長いブレード（D=160〜170m）と高いタワー（140〜160m）で掃過面積と高高度風を稼ぐのが定石です。GE 6.1-158、Vestas V162 などはまさにこのセグメント向けに開発されました。本ツールでクラス III、ハブ高 140m、D=160m と入れると、低風速サイトの設計感が掴めます。

WRF / OpenFOAM による微気象シミュレーション：地形が複雑な山岳・丘陵サイトでは、べき乗則の単純な外挿は誤差が大きく、メソスケール気象モデル WRF やマイクロスケール CFD（OpenFOAM ABL、WindSim、Meteodyn WT）で 3 次元の風速場を計算します。本ツールは机上での 1 次元プロファイル把握用と位置づけ、最終評価は CFD と LiDAR 観測で校正するのが実務フローです。

LiDAR リモートセンシング：ZephIR や WindCube などのドップラーライダーは地上から最大 300m までの風速プロファイルを直接測定でき、メットマスト（80〜100m）より高いハブ高の検証に必須のツールとなりました。観測値と本ツールのべき乗則・対数則の差を見れば、その立地の α・z₀ を逆算でき、年間発電量予測の精度を 5〜10% 向上できます。

よくある誤解と注意点

第一の落とし穴は、「α=0.143（1/7 乗則）はどこでも使える普遍値だ」という思い込みです。1/7 乗則は中立大気・滑らかな草地での経験則であり、洋上では 0.10、森林上空では 0.20〜0.30、市街地では 0.30 以上と大きく変動します。立地調査なしに 0.143 を当てはめると、洋上では発電量を過小評価、森林上空では過大評価する致命的なバイアスが出ます。最低 1 年のメットマスト観測を行い、α を立地ごとに同定してください。

第二に、「ハブ高風速で発電量を見積もる」のは大型機では危険です。デフォルト条件（D=150m、α=0.143）ですら REWS 補正で 6.4% の差が出ます。直径 200m を超える次世代機では 10% 以上の差になる場合もあり、年間収益で数億円の見積もり誤差につながります。風資源評価では必ず REWS（IEC 61400-12-1 Annex E に規定）か、3 次元 CFD 結合解析を行ってください。

第三に、「対数則の z₀ は地形写真から決め打ちでよい」という安易な選定。実際の z₀ は植生の高さ・密度・季節（落葉前後）・農作物の生育段階で変動します。例えば春の裸地（z₀≈0.01）と夏の麦畑（z₀≈0.05）では、ハブ高 100m での風速予測に 5% 以上の差が出ます。年間平均を取るなら季節別 z₀ の加重平均を、月別運用予測なら月別 z₀ を使い分けるのが厳密な手法です。

使い方ガイド

  1. 参照高さ（z_ref）と参照風速（V_ref）を入力します。例えば地上10m での実測値 12 m/s を指定します
  2. 風力タービンのハブ高（例：100 m）とロータ径（例：120 m）を入力し、べき乗則指数 α（地表面粗度クラスにより 0.14～0.40）を設定します
  3. 「計算」ボタンで、ハブ高での鉛直シア補正後の風速、上下端の相対風速差、等価風速 V_eq を自動演算し、想定出力を表示します

具体的な計算例

IEC 61400-1 クラス II（平均年間風速 V_ref=10 m/s at z_ref=10 m）の場合、ハブ高 90 m、ロータ径 110 m、地表面粗度カテゴリー B（α=0.19）で計算すると、べき乗則ハブ風速は V_hub=10×(90/10)^0.19≈15.2 m/s となります。ロータ上端（z=145 m）と下端（z=35 m）の風速差は約 3.8 m/s で、等価風速 V_eq≈15.6 m/s に補正されます。定格出力 3.6 MW の機種で、このシアプロファイル下での発電出力は約 2.95 MW となり、等価損失は約 18 % です

実務での注意点

  1. 地表面粗度カテゴリーは現地調査が必須です。都市部（α=0.35）と開放地（α=0.11）で同じ風速でも鉛直シアが大きく異なり、下部翼の応力が 25 % 以上変化します
  2. IEC 61400 の「等価風速」は、ロータ面全体の運動エネルギー平均値であり、ハブ高単点風速との乖離を無視すると出力予測誤差が生じます
  3. 台風・寒冷低気圧など突発的なシアイベント時には、対数則よりも実測値の直接入力が推奨されます
  4. ローター面での風速分布が非均一（竜巻など）の場合、このシミュレーターは適用外であり、数値流体力学解析が必要です