【风电控制】一张查表背后的生死线：深度拆解转速-转矩限幅曲线
一、前言：代码里最不起眼的1D Lookup Table
在翻阅风机变流器转矩控制代码时，你几乎一定会看到这样一个1D查表模块：
断点: [0, 0.625, 0.8, 1, 1.05, 1.15]
表数据: [0, 0, 1, 1, 1, 0.87]这个模块通常夹在"转矩指令生成"和"电流环"之间，只有一个 clamp 限幅函数的体量，很容易被忽略。
但就是这个不起眼的查表，是保护变流器不被炸毁的最后一道防线。
本文将从DFIG的电压/电流约束出发，用严格的数学推导解释这条曲线每一段的物理原因。
二、核心问题：这条曲线限的是什么？
2.1 误区排雷
❌ "这条曲线是为了保护电机不过载。"
✅ 这条曲线是为了保护机侧变流器（MSC）不因电压/电流饱和而失控。
电机本身的短时过载能力很强（通常是额定值的2~3倍，持续数秒甚至数分钟）。但变流器的IGBT和直流母线电容是"脆弱"的，一旦电压或电流超出能力范围，PWM调制将失效，系统瞬间崩溃。
2.2 两个核心约束
变流器的输出能力受限于两个硬约束：
约束1：电流约束（热极限）
Id2+Iq2≤Imax \sqrt{I_d^2 + I_q^2} \leq I_{max} Id2+Iq2 ≤Imax
其中 ImaxI_{max}Imax 由 IGBT 的额定电流和散热条件决定。
约束2：电压约束（调制极限）
Vd2+Vq2≤Vmax=Vdc3 \sqrt{V_d^2 + V_q^2} \leq V_{max} = \frac{V_{dc}}{\sqrt{3}} Vd2+Vq2 ≤Vmax=3 Vdc
其中 VdcV_{dc}Vdc 是直流母线电压，Vdc3\frac{V_{dc}}{\sqrt{3}}3 Vdc 是SVPWM线性调制区的最大输出电压。
转速越高，反电动势越大，可用于控制电流的"电压裕度"就越小。 这才是转速-转矩曲线存在的根本原因。
三、DFIG稳态电压方程与电压约束
3.1 d-q轴稳态电压方程
在同步旋转坐标系下，DFIG转子侧的稳态电压方程为：
Vrd=RrIrd−sωsLqIrq V_{rd} = R_r I_{rd} - s\omega_s L_q I_{rq} Vrd=RrIrd−sωsLqIrq
Vrq=RrIrq+sωs(LdIrd+ψf) V_{rq} = R_r I_{rq} + s\omega_s (L_d I_{rd} + \psi_f) Vrq=RrIrq+sωs(LdIrd+ψf)
忽略定子电阻（Rs≈0R_s \approx 0Rs≈0），在定子磁链定向（SFO）下简化：
Vrd=RrIrd−sωsσLrIrq V_{rd} = R_r I_{rd} - s\omega_s \sigma L_r I_{rq} Vrd=RrIrd−sωsσLrIrq
Vrq=RrIrq+sωsLmLsψs V_{rq} = R_r I_{rq} + s\omega_s \frac{L_m}{L_s} \psi_s Vrq=RrIrq+sωsLsLmψs
其中：
- σ=1−Lm2LsLr\sigma = 1 - \frac{L_m^2}{L_s L_r}σ=1−LsLrLm2：漏磁系数
- ψs\psi_sψs：定子磁链幅值
- Lm,Ls,LrL_m, L_s, L_rLm,Ls,Lr：互感、定子电感、转子电感
3.2 电压约束的转速依赖性
关键观察：转子电压中与转速相关的项（反电动势项）包含 sωss\omega_ssωs 因子 ：
Vrq,emf=sωsLmLsψs∝∣s∣×ωs V_{rq,emf} = s\omega_s \frac{L_m}{L_s} \psi_s \propto |s| \times \omega_s Vrq,emf=sωsLsLmψs∝∣s∣×ωs
这意味着：
∣s∣↑ ⟹ Vrq,emf↑ ⟹ 电压裕度↓ |s| \uparrow \implies V_{rq,emf} \uparrow \implies \text{电压裕度} \downarrow ∣s∣↑⟹Vrq,emf↑⟹电压裕度↓
当反电动势接近 VmaxV_{max}Vmax 时，PI控制器已经无法提供足够的电压去调节电流，系统失控。
四、曲线分段硬核解析
4.1 第一段：低速死区（ωm∗∈[0,0.625]\omega_m^* \in [0, 0.625]ωm∗∈[0,0.625]）
表数据：Tlim∗=0T_{lim}^* = 0Tlim∗=0
物理原因：
在极低转速（ωm∗<0.625\omega_m^* < 0.625ωm∗<0.625，即 nr<937.5n_r < 937.5nr<937.5 RPM）下：
- 定子电压过低 ：定子感应电动势 es∝ωme_s \propto \omega_mes∝ωm，转速过低时，定子电压远低于电网电压。如果此时并网，将产生巨大的无功涌入电流。
- 转子电压要求异常 ：虽然反电动势低，但定子磁链定向要求定子磁链幅值恒定，这需要极大的转子励磁电流。在 ∣s∣>0.375|s| > 0.375∣s∣>0.375 时：
Vro=∣s∣×Vs,open_circuit>0.375×Vs V_{ro} = |s| \times V_{s,open\_circuit} > 0.375 \times V_s Vro=∣s∣×Vs,open_circuit>0.375×Vs
这已经逼近或超过变流器的电压能力。 - 控制不可观测 ：在极低转差频率下（fr<18.75f_r < 18.75fr<18.75 Hz），定子磁链观测器的精度急剧下降。
设计决策 ：直接将转矩限幅设为 0，禁止出力。但不禁止并网------系统仍可进行空载励磁和并网操作。
4.2 第二段：出力门槛（ωm∗∈[0.625,0.8]\omega_m^* \in [0.625, 0.8]ωm∗∈[0.625,0.8]）
表数据：Tlim∗=0T_{lim}^* = 0Tlim∗=0（但接近0.8时内部已有裕度）
随着转速升高，转子电压要求降低，但此时变流器的大部分容量仍被用于维持定子磁链恒定所需的励磁电流（IrdI_{rd}Ird），留给产生转矩的电流分量（IrqI_{rq}Irq）的裕度仍不足。
4.3 第三段：阶跃------出力解放（ωm∗=0.8\omega_m^* = 0.8ωm∗=0.8）
表数据：Tlim∗T_{lim}^*Tlim∗ 从 0 跃变到 1.0
为什么是 0.8 p.u.？这是一个精确的工程计算结果。
在 ωm∗=0.8\omega_m^* = 0.8ωm∗=0.8 时（即 nr=0.8×1680=1344n_r = 0.8 \times 1680 = 1344nr=0.8×1680=1344 RPM）：
转差率：
s=1500−13441500=0.104 s = \frac{1500 - 1344}{1500} = 0.104 s=15001500−1344=0.104
此时我们可以计算变流器的电压裕度是否足够容纳转矩电流。
假设：
- 直流母线电压 Vdc=650V_{dc} = 650Vdc=650 V
- SVPWM最大输出电压 Vmax=6503=375.3V_{max} = \frac{650}{\sqrt{3}} = 375.3Vmax=3 650=375.3 V
- 额定转子电压 Vr,rated=∣srated∣×Vs,rated≈0.12×380=45.6V_{r,rated} = |s_{rated}| \times V_{s,rated} \approx 0.12 \times 380 = 45.6Vr,rated=∣srated∣×Vs,rated≈0.12×380=45.6 V
在 s=0.104s = 0.104s=0.104 时，转子反电动势：
Vr,emf=0.104×380×LmLs≈35 V (估算) V_{r,emf} = 0.104 \times 380 \times \frac{L_m}{L_s} \approx 35 \text{ V (估算)} Vr,emf=0.104×380×LsLm≈35 V (估算)
电压裕度：
ΔV=Vmax−Vr,emf=375.3−35=340.3 V \Delta V = V_{max} - V_{r,emf} = 375.3 - 35 = 340.3 \text{ V} ΔV=Vmax−Vr,emf=375.3−35=340.3 V
这个裕度远远超过了产生额定转矩电流所需的最大电压（通常不超过 50~80V）。因此，系统在此处"判定"变流器能力充足，允许输出满载转矩。
0.8 p.u. 的工程含义 ：它不是转速的10%余量，而是变流器电压裕度刚好能够容纳满载转矩电流的最低转速点。这是一个由变流器硬件参数精确计算得出的阈值。
4.4 第四段：恒转矩区（ωm∗∈[0.8,1.05]\omega_m^* \in [0.8, 1.05]ωm∗∈[0.8,1.05]）
表数据：Tlim∗=1.0T_{lim}^* = 1.0Tlim∗=1.0
在这个范围内，反电动势虽然持续增长，但变流器仍有足够的电压裕度：
| 转速 (p.u.) | 转速 (RPM) | $ | s | $ | Vr,emfV_{r,emf}Vr,emf (估算V) | 电压裕度 ΔV\Delta VΔV (V) | 转矩限幅 |
|---|---|---|---|---|---|---|---|
| 0.8 | 1344 | 0.104 | 35 | 340 | 1.0 | ||
| 0.9 | 1512 | -0.008 | 2.7 | 373 | 1.0 | ||
| 1.0 | 1680 | -0.12 | 40 | 335 | 1.0 | ||
| 1.05 | 1764 | -0.176 | 59 | 316 | 1.0 |
电压裕度始终远超所需，因此转矩限幅保持 1.0 p.u.。
4.5 第五段：弱磁限幅区（ωm∗∈[1.05,1.15]\omega_m^* \in [1.05, 1.15]ωm∗∈[1.05,1.15]）
表数据：Tlim∗T_{lim}^*Tlim∗ 从 1.0 下降到 0.87
这是变流器电压饱和 导致的必然结果。
在 ωm∗=1.15\omega_m^* = 1.15ωm∗=1.15 时（nr=1932n_r = 1932nr=1932 RPM）：
s=1500−19321500=−0.288 s = \frac{1500 - 1932}{1500} = -0.288 s=15001500−1932=−0.288
转子反电动势急剧增大：
Vr,emf≈0.288×380×LmLs≈97 V (估算) V_{r,emf} \approx 0.288 \times 380 \times \frac{L_m}{L_s} \approx 97 \text{ V (估算)} Vr,emf≈0.288×380×LsLm≈97 V (估算)
加上转矩电流所需的电阻压降和漏抗压降，总电压需求可能已经逼近 VmaxV_{max}Vmax。
弱磁策略启动：
当电压裕度不足时，必须注入负 IdI_dId 电流 来削弱气隙磁链：
Id<0 ⟹ ψgap=ψs+LdId↓ ⟹ Vr,emf↓ I_d < 0 \implies \psi_{gap} = \psi_s + L_d I_d \downarrow \implies V_{r,emf} \downarrow Id<0⟹ψgap=ψs+LdId↓⟹Vr,emf↓
但总电流受到热约束：
Id2+Iq2≤Imax \sqrt{I_d^2 + I_q^2} \leq I_{max} Id2+Iq2 ≤Imax
IdI_dId 占用了电流空间，留给产生转矩的 IqI_qIq 就必须减小：
Iq=Imax2−Id2↓ ⟹ Te∝Iq↓ I_q = \sqrt{I_{max}^2 - I_d^2} \downarrow \implies T_e \propto I_q \downarrow Iq=Imax2−Id2 ↓⟹Te∝Iq↓
计算验证：为什么降到 0.87？
假设弱磁所需的 IdI_dId 占总电流的 ImaxI_{max}Imax 的比例使得可用 IqI_qIq 降至 1−Id2/Imax2=0.87\sqrt{1 - I_d^2/I_{max}^2} = 0.871−Id2/Imax2 =0.87，即：
IdImax=1−0.872=1−0.7569=0.2431≈0.493 \frac{I_d}{I_{max}} = \sqrt{1 - 0.87^2} = \sqrt{1 - 0.7569} = \sqrt{0.2431} \approx 0.493 ImaxId=1−0.872 =1−0.7569 =0.2431 ≈0.493
这意味着约 49.3% 的变流器电流容量被用于弱磁，只有 87% 的容量留给转矩输出。
五、曲线的数学模型
整个转速-转矩限幅曲线可以表示为：
Tlim∗(ωm∗)={0ωm∗<0.81.00.8≤ωm∗≤1.051.0ωm∗ωm∗>1.05(恒功率约束) T_{lim}^*(\omega_m^*) = \begin{cases} 0 & \omega_m^* < 0.8 \\ 1.0 & 0.8 \leq \omega_m^* \leq 1.05 \\ \frac{1.0}{\omega_m^*} & \omega_m^* > 1.05 \quad \text{(恒功率约束)} \end{cases} Tlim∗(ωm∗)=⎩ ⎨ ⎧01.0ωm∗1.0ωm∗<0.80.8≤ωm∗≤1.05ωm∗>1.05(恒功率约束)
验证：在 ωm∗=1.15\omega_m^* = 1.15ωm∗=1.15 时：
Tlim∗=1.01.15=0.8696≈0.87 T_{lim}^* = \frac{1.0}{1.15} = 0.8696 \approx 0.87 Tlim∗=1.151.0=0.8696≈0.87
与查表数据完全吻合！这证明了高转速段的限幅本质是恒功率约束（P=T×ω=constP = T \times \omega = \text{const}P=T×ω=const）。
六、限幅模块在代码中的位置
c
// Step 1: 将机械角速度转换为标幺转速
float omega_pu = RotorMechOmega * K_GAIN; // K_GAIN = 0.005684105
// Step 2: 查表得到当前转速下的最大转矩限幅
float T_max_pu = lookup_table(omega_pu, breakpoints, table_data);
// Step 3: 对转矩指令进行钳位
float T_ref_pu = power_controller_output; // 功率外环输出的原始转矩指令
T_ref_pu = clamp(T_ref_pu, -T_max_pu, T_max_pu);
// Step 4: 将标幺转矩转换为Iq电流指令
float Iq_ref = T_ref_pu * T_rated / (1.5 * pole_pairs * psi_f);七、"能并网"与"能出力"的本质区别
| 功能 | 要求的转速范围 | 对应表数据 |
|---|---|---|
| 空载并网（励磁同步） | ωm∗>0.4\omega_m^* > 0.4ωm∗>0.4 即可 | 与表无关（转矩为0） |
| 并网后出力（发电） | ωm∗≥0.8\omega_m^* \geq 0.8ωm∗≥0.8 | Tlim∗≥1.0T_{lim}^* \geq 1.0Tlim∗≥1.0 |
| 满载出力 | ωm∗∈[0.8,1.05]\omega_m^* \in [0.8, 1.05]ωm∗∈[0.8,1.05] | Tlim∗=1.0T_{lim}^* = 1.0Tlim∗=1.0 |
关键结论 ：在 0.625∼0.80.625 \sim 0.80.625∼0.8 p.u. 区间，转矩限幅为 0，但系统完全可以进行空载并网。并网后处于"同步旋转但零有功输出"状态，等待转速升至 0.8 后才开始发电。
八、结语
这条 [0, 0, 1, 1, 1, 0.87] 的曲线不是拍脑袋定的，它是：
- 低速段（0~0.8） ：变流器电压/电流裕度不足以同时支持励磁和转矩 →\rightarrow→ 禁止出力
- 中速段（0.8~1.05） ：裕度充足 →\rightarrow→ 允许满载转矩
- 高速段（1.05~1.15） ：电压饱和触发弱磁 →\rightarrow→ 恒功率降转矩
在最后一篇中，我们将把这一切串联起来，从增益 KKK 到查表，到限幅，再到最终的 Id/IqI_d/I_qId/Iq 电流指令生成，手撕完整的转矩控制代码逻辑。