风速剖面 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 点采用三次方加权平均,得到能量等效风速。

$$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。每次转动叶片都经历这个差,导致叶根和主轴产生周期性疲劳荷载。IEC 规范中风切变/轮毂风速超过 0.2 时疲劳评估就变严格了。本工具会自动检查这个比值并显示结果。
🙋
发电量估算能直接用轮毂高风速吗?
🎓
这就是「等效风速(REWS)」的作用。风机实际获得的功率由转子面平均 V³ 决定。默认条件下 V_top=12.0、V_hub=11.1、V_bot=9.1,三次方加权平均得 V_eq=10.86 m/s。用轮毂高算出 6691 kW,用等效风速算出 6262 kW,差了 6.4%。忽视这点会高估年收益。现代风资源评估标准是用激光雷达直接测转子风速分布,然后用 REWS 计算。
🙋
大气稳定度对 α 影响有多大?
🎓
影响很大。稳定大气(夜间、地表冷却混合弱)α 可达 0.3~0.5。反之不稳定大气(晴朗白天、对流活跃)α≈0.10。同一立地夜间和白天剖面完全不同,年均发电需要用 Monin-Obukhov 长度做稳定度修正才准确。工具的「大气稳定度」选择器会显示 α 推荐值,无现场观测时可作为初始参考。
🙋
沿海(IEC I)和内陆(II/III)分等级,就是为了考虑这些差异吧。
🎓
没错。IEC 61400-1 用 50 年重现期 10 分钟平均最大风速定义参考风速:I=50、II=42.5、III=37.5、IV=30 m/s,年平均约 0.2 倍参考值。沿海风强粗度低用 I,低风速内陆丘陵用 III/IV。反过来把 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 时,叶片在每次旋转中受到周期性变化的气动荷载,对偏航轴、主轴、叶根产生疲劳荷载。IEC 61400-1 在标准风切变模型(NWP, EWS)中使用 α=0.20 作为疲劳评估基准。本工具对 ΔV/V_hub > 0.20 显示 NG(疲劳过大),>0.15 显示警告。直径 150m、轮毂高 100m 的风机即使在 α=0.143 下,上下风速差也约 3 m/s,是疲劳寿命的主要因素之一。
转子等效风速(REWS)是将转子面内风速分布获得的总动能转换为「转子面平均风速」的值,计算公式为 V_eq = ((V_top³+V_hub³+V_bot³)/3)^(1/3),采用三次方加权平均。由于风机输出与 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 等中尺度气象模型或 OpenFOAM ABL、WindSim、Meteodyn WT 等微尺度 CFD 算三维风场。本工具定位是快速一维剖面把握,最终评估应结合 CFD 和激光雷达观测校准。

激光雷达遥感: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 附录 E 规定)或三维 CFD 耦合分析。

第三个陷阱:「从地形图片决定对数律 z₀ 就行」太草率。实际 z₀ 随植被高度、密度、季节(落叶前后)、农作物生育期变化。春天裸地(z₀≈0.01)和夏季麦田(z₀≈0.05)轮毂高 100m 预测风速差 5% 以上。年均应用需季节加权 z₀,月度运营用月别 z₀。

使用指南

  1. 输入参考高(z_ref)和参考风速(V_ref)。例如地面 10m 处实测 12 m/s
  2. 输入风力涡轮机的轮毂高(例 100m)和转子直径(例 120m),设置幂指数 α(根据地表粗度 0.14~0.40)
  3. 点击「计算」,自动演算轮毂高垂直风切变修正风速、上下相对风速差、等效风速 V_eq,显示预期输出

具体计算例

IEC 61400-1 II 级(年均参考风速 V_ref=10 m/s at z_ref=10m)、轮毂高 90m、转子直径 110m、地表粗度 B 类(α=0.19),幂律则轮毂风速为 V_hub=10×(90/10)^0.19≈15.2 m/s。转子上端(z=145m)和下端(z=35m)风速差约 3.8 m/s,等效风速 V_eq≈15.6 m/s。定格 3.6MW 机组在此风切剖面下发电输出约 2.95MW,等效损失约 18%。

实务注意

  1. 地表粗度分类必须现场调查。都市(α=0.35)与开放地(α=0.11)同风速下风切变差异大,下部叶片应力变化超 25%
  2. IEC 61400 「等效风速」是转子整体运动能量平均,忽视与轮毂单点风速差异会产生出力预测偏差
  3. 台风、寒流等突发风切变事件,对数律不适用,推荐直接输入实测值
  4. 转子面风速分布非均(如龙卷风),本工具不适用,需数值流体力学详细分析