3. 风速与风轮特性建模

在转矩与功率控制中，风速不是“一个数”，而是同时决定可用风功率上限、气动转矩大小和控制器工作点的外部输入。风轮也不是固定效率的能量转换器，它的转换能力随风速、转速和桨距角变化。因此，本节把风速模型、风能利用系数 \(C_p\)、叶尖速比 \(\lambda\) 和桨距角 \(\beta\) 放在同一条气动链路上理解：风速与转速共同决定 \(\lambda\)，\(\lambda\) 与 \(\beta\) 共同决定 \(C_p\)，\(C_p\) 再与风速立方共同决定可捕获功率，最终通过转矩—转速关系影响机组运行。

一、风速模型：从自然风到控制器可用输入

自然风可分解为四部分：平均风速决定当前风资源水平；阵风表现为短时幅值抬升；渐变风表现为缓慢爬升或下降；湍流表现为围绕平均值的随机脉动。工程建模常用“平均风 + 湍流分量”的形式：

v(t) = v_mean + v_turb(t) + v_gust(t) + v_ramp(t)

其中 \(v_{mean}\) 可由威布尔分布或现场测风数据统计得到，湍流分量常通过冯·卡门或凯马尔谱生成。对控制而言，关键不是复现每一秒风速，而是明确控制器看到的是轮毂高度风速估计值，它经过滤波、限幅和延迟后参与转矩或桨距控制。若风速模型只给平均风速，仿真会低估转矩波动和变桨动作频次；若湍流强度设置过大，又会掩盖控制策略本身的差异。因此，风速模型必须与后续气动模型、传动链模型和控制器采样周期匹配。

二、叶尖速比与风能利用系数：风轮特性的核心关系

叶尖速比定义为叶尖线速度与风速之比：

λ = ωR / v

其中 \(ω\) 为风轮角速度，\(R\) 为风轮半径，\(v\) 为轮毂高度风速。\(\lambda\) 的物理含义是“叶尖相对于风跑得有多快”。\(\lambda\) 过小，叶片像一块挡风板，气流分离严重；\(\lambda\) 过大，叶片相对来流角度过小，升力下降且阻力增加。因此存在一个最佳叶尖速比 \(\lambda_{opt}\)，使风能利用系数 \(C_p\) 达到最大值 \(C_{p,max}\)。

风能利用系数定义为风轮捕获功率与来流风功率之比：

P_mech = 0.5 · ρ · A · v³ · C_p(λ, β)
A = πR²

其中 \(ρ\) 为空气密度，\(A\) 为风轮扫风面积。贝兹极限说明理想风轮 \(C_{p,max}\) 不超过 0.593，实际三叶片机组通常在 0.45～0.50 之间。\(C_p\) 不是常数，而是 \(\lambda\) 和 \(\beta\) 的二元函数。常见工程拟合形式为：

C_p(λ, β) = c1·(c2/λ_i - c3·β - c4)·exp(-c5/λ_i) + c6·λ
1/λ_i = 1/(λ + 0.08β) - 0.035/(β³ + 1)

这组系数因机型而异，只用于说明 \(C_p\) 随 \(\lambda\) 和 \(\beta\) 连续变化。实际建模时，应优先使用厂家提供的 \(C_p\) 曲线或查表数据，而不是随意套用通用系数。

三、桨距角对转矩与功率的调节作用

桨距角 \(\beta\) 是叶片绕变桨轴旋转的角度。\(\beta=0°\) 附近通常对应最大升力与最大 \(C_p\)；\(\beta\) 增大时，叶片迎风角度改变，升力下降、阻力上升，\(C_p\) 整体下移。这正是变桨控制能够限制功率的物理基础：当风速超过额定值，控制器不再追求最大 \(C_p\)，而是主动增大 \(\beta\)，使气动转矩与发电机反力矩重新平衡，把转速和功率稳定在额定值附近。

从控制角度看，\(\beta\) 对转矩的影响具有明显非线性：小桨距角变化对 \(C_p\) 影响较小，大桨距角附近变化更敏感。因此变桨控制器通常采用增益调度，在不同桨距角区间使用不同 PI 参数。若把 \(\beta\) 当作线性执行器处理，在高风速段容易出现超调或振荡。

四、把模型串起来：从风速到转矩的完整计算

假设某 2 MW 机组：\(R=40\) m，\(ρ=1.225\) kg/m³，\(A=πR²≈5027\) m²。取风速 \(v=10\) m/s，风轮转速对应 \(\lambda=8\)，\(\beta=0°\)，查 \(C_p\) 曲线得 \(C_p≈0.47\)。则：

P_mech = 0.5 × 1.225 × 5027 × 10³ × 0.47
       ≈ 0.5 × 1.225 × 5027 × 1000 × 0.47
       ≈ 1.447 × 10⁶ W ≈ 1.45 MW

若风速升到 12 m/s，仍保持 \(\lambda=8\)、\(\beta=0°\)，则 \(v³\) 从 1000 变为 1728，功率约升至 2.50 MW，超过额定值。此时控制器必须增大 \(\beta\)，把 \(C_p\) 从 0.47 降到约 0.38，使功率回落到 2 MW 附近。这个例子说明：风速决定“有多少能量”，\(C_p\) 决定“能拿多少”，\(\beta\) 决定“要不要少拿”。

五、建模时的关键检查点

• 风速参考高度：轮毂高度与测风高度不一致时，需用风切变指数换算，否则 \(v\) 系统性偏差会直接放大到 \(v³\)。
• \(C_p\) 数据来源：优先使用厂家曲线或查表，拟合公式只作插值补充，并检查 \(\lambda\) 和 \(\beta\) 的定义范围。
• λ 的转速口径：确认 \(ω\) 是风轮侧还是发电机侧，传动比是否已折算，避免 \(\lambda\) 整体偏移。
• β 的零位定义：明确 \(0°\) 对应最大升力位置还是顺桨位置，不同定义下 \(C_p\) 曲线方向相反。
• 空气密度修正：高原或高温环境 \(ρ\) 下降，相同风速下转矩和功率同步降低，控制参数需相应调整。

综上，风速模型提供输入，\(\lambda\) 描述风轮与风速的匹配程度，\(C_p\) 描述气动转换效率，\(\beta\) 提供主动调节手段，四者共同决定风轮气动转矩和输出功率。后续转矩控制与变桨控制策略，都是在这条气动链路基础上设计工作点和动态补偿。