最近在做一个海上风电直流并网的预研项目,把一台5MW永磁直驱发电机在Simulink里完整搭建起来,直流母线电压定在1200V。整套模型从风力机气动特性、永磁同步发电机本体参数、机侧变流器控制,到直流并网侧的控制全部串了起来,验证功率传输、母线电压稳定性和低电压穿越策略都非常方便。这篇文章把我实际搭模型的过程、参数计算思路、控制策略选型逻辑和踩过的坑全部写出来,适合正在做风电变流器仿真、或者想快速掌握永磁直驱发电机直流并网控制逻辑的工程师和研究生参考。
1. 整体设计与思路拆解
1.1 机型选型:为什么用永磁直驱而不是双馈
风电领域最主流的两条技术路线,一个是双馈异步发电机(DFIG)加齿轮箱,另一个是永磁同步发电机(PMSG)直驱。做直流并网方案时,我毫不犹豫选了永磁直驱,理由很直接:DFIG的转子绕线需要滑环和碳刷,机组在海上高湿高盐环境下维护成本非常高;而永磁直驱没有齿轮箱,电机主轴直接连接风力机,传动链简化了一大截,故障率低很多。更关键的是,永磁直驱必须配全功率变流器,这意味着发电机和电网之间完全解耦,在直流并网场景下天然适配。
全功率变流器这个词听起来复杂,但落地到模型里就很好理解:发电机出来的电能,全部经过机侧变流器整流成直流,再通过并网侧变流器转换成所需电压等级的电能。发电机频率变化、转速波动,都不会直接传导到直流侧,控制上可以分两边独立处理。用生活化的类比,DFIG像自行车链条传动,齿轮箱把脚蹬转速调整到合适范围;永磁直驱更像轮毂电机电动车,电机直接驱动轮子,省掉中间传动机构,结构简单了,但驱动控制电路要全功率承受电能。
在Simulink里模拟永磁直驱还有一个额外好处:PMSG的转子磁场恒定,可以用转子磁场定向控制简化成直流电机特性,控制算法直观得多,不像异步电机还要考虑励磁电流和转子磁场估算。
1.2 直流并网方案:1200V母线电压是怎么定下来的
传统风电场通常是交流汇集,每台机组升压到35kV再集中升压到110kV或更高。直流并网则完全不同,机组发出的交流电就地整流成直流,经过直流汇集线路送到海上换流站或者直接在微电网里用。直流方案在海上风电里特别有优势,因为交流海缆存在大量容性无功电流,传输距离一长损耗急剧上升,而直流电缆没有这个问题,还可以通过直流母线天然形成多端互联网络。
1200V这个电压等级不是拍脑袋定的,核心是变流器功率器件和电缆绝缘的综合权衡。5MW功率等级下,如果用690V交流并网,母线电压通常在1000V左右,电流接近4000A,母线铜排和电缆截面会非常恐怖,开关器件并联数量也多。把直流母线提升到1200V,按照容量不变、功率等于电压乘电流来算,电流能降到3500A以下,损耗随电流平方下降,这是个很可观的收益。再往上提到1500V甚至2000V,当然还能进一步降低电流,但1200V的系统采用1700V耐压的IGBT模块时,电压裕量在1.4倍左右,工作应力合理,过压保护也有充足空间,属于工程上比较成熟的平衡点。
从绝缘系统角度看,1200V直流母线的绝缘设计比低压690V系统要高一个等级,但还不需要引入中压绝缘设计的复杂度。如果做到3000V以上,电缆接头、断路器和变压器绝缘方案全都要换一套体系,系统成本反而上去了。所以5MW这个级别,1200V直流母线是性价比最优的选择。
1.3 Simulink模型整体架构
整套Simulink模型采用模块化分层设计,顶层模型包含四个子系统:
- 风力机与传动链模型:输入风速,输出机械转矩。
- PMSG电机模型:电磁关系,输出定子电流和电磁转矩。
- 机侧变流器及控制系统:实现最大功率跟踪和发电机转速控制。
- 直流并网系统:包含直流母线电容、并网变换器和并网控制策略。
这样划分的好处是调试时可以单独走通每个子系统。我习惯先开环测试电机模型,确认空载反电动势和额定转矩正确无误,再闭环控制,最后再整机联调。如果模型一上来就全部堆在一起,出了问题很难定位是电机参数错了、控制器发散还是并网侧振荡,排查成本非常高。
2. 核心参数设计与数学模型
2.1 5MW风力机关键参数计算
仿真模型要让人信服,第一步是风机参数基本靠谱。5MW风力机可以用以下参数作基准:
- 叶轮直径:126m,扫风面积约12469平方米;
- 额定风速:11.4m/s;
- 切入风速:3m/s,切出风速:25m/s;
- 额定转速:约12rpm,对应1.26rad/s;
- 齿轮箱:无,直驱连接。
风力机气动模型的核心是贝兹理论。捕获的机械功率为:
P_m = 0.5 × ρ × A × V³ × Cp(λ, β)
其中ρ取1.225kg/m³,A是扫风面积,V是风速,Cp是风能利用系数,λ是叶尖速比,β是桨距角。仿真里Cp可以用经验公式拟合,典型的最大值在0.48左右,对应最佳叶尖速比约8.5。
桨距角在额定风速以下保持最优值不变,额定风速以上需要通过变桨系统限制捕获功率。做额定风速以下仿真时,可以把桨距角固定,只做最大功率跟踪控制;做高风速工况模拟时,才考虑加入变桨控制。
额定转速对应的发电机机械角速度,5MW直驱机组大约在1.0到1.5rad/s。永磁电机极对数设计得比较多,比如取100极对,那么电气频率约为:
f = n × p / 60
转速按12.6rpm算,极对数100,频率就是21Hz。频率不高,但转矩非常大,因为功率等于转矩乘以转速,转速低时转矩必然要上去。5MW除以1.32rad/s,机械额定转矩接近3790kN·m,这个数据对后面电机参数选择和变流器设计很关键。
2.2 PMSG在dq坐标系下的数学模型
永磁同步电机在转子磁场定向的dq坐标系下,定子电压方程可以写成:
u_d = R_s × i_d + L_d × di_d/dt - ω_e × L_q × i_q
u_q = R_s × i_q + L_q × di_q/dt + ω_e × (L_d × i_d + ψ_f)
其中u_d、u_q是定子d、q轴电压,i_d、i_q是定子d、q轴电流,L_d、L_q是d、q轴电感,R_s是定子电阻,ψ_f是永磁磁链,ω_e是电气角速度。
电磁转矩方程为:
T_e = 1.5 × p × [ψ_f × i_q + (L_d - L_q) × i_d × i_q]
对于表贴式永磁电机,L_d等于L_q,转矩简化为:
T_e = 1.5 × p × ψ_f × i_q
这说明id控制到零时,转矩完全由iq决定,控制和直流电机一模一样。对内置式永磁电机,L_d小于L_q,可以注入负的id来利用磁阻转矩,但5MW直驱大电机通常采用表贴式或半内置结构,id=0策略足够用。
用之前参数估算:3790kN·m额定转矩,极对数100,那么ψ_f × i_q约为25.27kN·m/A。如果额定电流的q轴分量取1000A,永磁磁链大约0.25Wb左右。实际项目中电机厂家会给出精确磁链值,仿真这里用这些倒推的值初始化模型,特性已经非常接近真实机组。
2.3 直流母线电容选取
直流母线电容是仿真中必须认真对待的元件。它的作用有两方面:一是稳定母线电压,缓冲机侧和并网侧之间的瞬时功率不平衡;二是为高频开关纹波提供低阻抗通路。
电容容量的工程估算,通常按母线电压纹波要求来。假设在功率阶跃过程中允许母线电压波动5%,即1200V母线的60V,一周期内能量波动导致的电压变化满足:
C ≥ P × T / (U_bus × ΔU)
取输出功率5MW,开关周期0.5ms(2kHz开关频率),算出来是:
C ≥ 5e6 × 0.5e-3 / (1200 × 60) ≈ 34.7mF
实际会留一些裕量,取40到50mF比较合理。有些仿真模型为了加快仿真速度,把电容值设得远小于物理值,仿真确实能跑通,但电压纹波明显偏大,稳压特性失真。我建议还是按物理参数设置,这样后续做故障穿越和负载扰动时结果才有参考价值。
电容等效串联电阻ESR也要加上,否则母线电压会出现高频尖峰。Simulink里电解电容模型可以直接填等效串联电阻,一般取0.1mΩ量级。
3. 控制策略设计与原理
3.1 机侧变流器:最大功率跟踪控制
额定风速以下,风机要尽可能捕获最大风能,这就是最大功率跟踪(MPPT)控制的含义。从2.1节公式能看出,在桨距角最优时,Cp是叶尖速比λ的函数,而叶尖速比等于叶尖线速度除以风速:
λ = ω_r × R / V
要保持λ在最优值,转速必须跟随风速变化。也就是说,风速变了,发电机转速要跟着变,这需要通过控制电磁转矩来实现。转子运动方程:
J × dω_r/dt = T_m - T_e
其中J是转动惯量,T_m是风力机机械转矩,T_e是电磁制动转矩。当电磁转矩等于机械转矩时,转速稳定。
工程上最常用的MPPT方法是功率曲线法。在最优叶尖速比下,最优功率和转速有确定关系:
P_opt = k_opt × ω³
其中k_opt可以由风力机参数计算得到:k_opt = 0.5 × ρ × A × R⁵ × Cp_max / λ_opt³。在Simulink里,用这个公式生成转速-功率曲线,转速反馈作为输入,查表得到功率参考值,再除以转速得到转矩参考值,作为电流环的指令。
传递函数的角度看,机侧控制采用双闭环结构:外环是转速环,输出q轴电流参考值;内环是电流环,输出dq轴电压指令。电流环带宽一般设计在100到200Hz,转速环带宽在1到5Hz,内外环带宽差一个数量级,保证两个环不会互相干扰。我在仿真里把转速环PI参数整定到带宽2Hz左右,电流环带宽150Hz,实际调节效果比较理想。
3.2 直流母线电压控制策略
直流并网系统的母线电压稳定是整个系统的核心控制目标。母线电压的波动本质反映的是功率不平衡,母线电容上的功率守恒关系为:
C × U_bus × dU_bus/dt = P_gen - P_grid
其中P_gen是机侧注入母线的功率,P_grid是并网侧输出的功率。要让母线电压稳定,并网侧必须快速平衡机侧功率变化。常规做法是并网侧采用电压外环加电流内环的双闭环控制:电压外环检测母线电压误差,输出电流指令,电流内环跟踪电流指令并输出调制电压。
在设计电压环时,要注意带宽必须低于电流环响应速度。否则母线电压有一点波动,电流指令立即大幅跳变,电流环来不及跟随,整个系统就容易振荡。我习惯将电压环带宽设在10到20Hz,比电流环低5到10倍,这样稳压效果和动态响应兼顾得比较好。
3.3 并网侧功率传输控制
直流并网的输出侧取决于接入的对象。如果接入交流电网,需要通过逆变器将直流变换为交流,并采用电网电压定向矢量控制,实现有功和无功的解耦控制。如果接入的是直流微网或者另一个直流母线,那就由直流变换器(如DC/DC)控制输出电流或功率。
从标题来看,这个项目关注的是“1200V直流并网”,最典型的情景是:
- 经过DC/DC升压后接入更高电压等级的直流汇集母线;
- 或者通过DC/AC逆变器接入中压交流电网。
如果输出侧是单向直流变换器,电流环指令由功率外环或者电压外环给出。如果输出侧是双向AC/DC变换器,网侧还需要锁相环(PLL)同步电网电压相位。这里特别提醒一句,PLL在弱电网下容易出相位抖动问题,仿真里可以把电网阻抗加进去观察一下,别只看理想电压源的情况。
4. Simulink模型搭建实操
4.1 顶层模型搭建
打开Simulink,新建一个空白模型,先把顶层框架搭出来。我建议顺序是:
- 左侧放风速输入和风力机子系统;
- 风力机输出机械转矩连接到PMSG电机子系统;
- 电机定子三相输出接机侧变流器;
- 机侧变流器直流侧接直流母线电容;
- 母线再接并网侧变换器和电网模型。
模型里我习惯用Simscape Electrical组件库里的永磁同步电机和三相桥式变流器,因为电气物理模型比纯数学方程更接近真实系统。如果不用Simscape,也可以用SimPowerSystem里的电机模块。前者适合做精细电气仿真,后者接口更简洁、仿真速度更快,看个人需求选。
4.2 风力机气动模型
在Simulink里可以用MATLAB Function模块写一个风力机气动模型函数,输入是风速V、转速ω_r和桨距角β,输出是机械转矩T_m。核心代码如下:
function Tm = wind_turbine(V, wr, beta) % 风力机参数 rho = 1.225; % 空气密度 kg/m^3 R = 63; % 叶轮半径 m (直径126m) A = pi * R^2; % 扫风面积 m^2 % 叶尖速比计算 lambda = wr * R / V; % Cp经验公式(典型参数) % 系数根据实际机组数据拟合 c1 = 0.5176; c2 = 116; c3 = 0.4; c4 = 5; c5 = 21; c6 = 0.0068; lambda_i = 1 / (1/(lambda + 0.08*beta) - 0.035/(beta^3 + 1)); Cp = c1 * (c2/lambda_i - c3*beta - c4) * exp(-c5/lambda_i) + c6*lambda; % 捕获机械功率 Pm = 0.5 * rho * A * V^3 * Cp; % 机械转矩 Tm = Pm / wr; if wr < 0.01 Tm = 0; % 防止除零 end end这段代码写完以后,有必要先用脚本扫一遍风速-转速-输出功率关系,确认在额定风速附近能输出5MW。如果偏差较大,检查Cp函数峰值系数是否正确,特别是c1到c6这组系数,它们直接决定了拟合出来的Cp曲线高度和位置。
在机械连接上,风力机输出转矩接到电机的机械端口。为了模拟传动轴的柔性,我在中间加了一段等效弹簧阻尼模型,用Simscape里的Rotational Spring Damper模块实现。直驱机组虽然省掉齿轮箱,但叶片和轮毂的扭振模态仍然存在,不加柔性传动模型的话,动态响应会显得过于“刚硬”,做低电压穿越仿真时差异更明显。
4.3 PMSG电机及变流器参数设置
电机参数在Simscape Electrical的PMSG模块中设置:
- 定子电阻Rs:5mΩ;
- d轴电感Ld:0.3mH,q轴电感Lq:0.3mH;
- 永磁磁链ψf:0.25Wb;
- 极对数:100;
- 额定机械转矩:约3790kN·m。
变流器采用三相两电平IGBT桥,IGBT通态压降设置为2V左右,导通电阻设1mΩ。直流侧电容按2.3节计算的45mF设置。
需要注意,Simscape Electrical里的电机模型默认是国际单位制,转矩单位是N·m,功率单位是W,转速单位是rad/s。如果直接把工程单位(转速用rpm,功率用MW)填进去,你会看到仿真结果莫名其妙差了好几倍。我在第一次搭建模型时踩过这个坑,后来养成习惯,所有参数先统一单位换算,再填入模块对话框。
4.4 控制器搭建与PI参数初调
电流环控制器用Simulink标准PID模块即可。在dq旋转坐标系下,PI控制器的输出还要加入交叉耦合补偿项和反电动势前馈,这一步很多人会忽略,但不加前馈的结果就是系统动态响应慢、静态误差大。
q轴电压指令为:
u_q_ref = PI(i_q_ref - i_q) + ω_e × L_d × i_d + ω_e × ψ_f
d轴电压指令为:
u_d_ref = PI(i_d_ref - i_d) - ω_e × L_q × i_q
这样处理之后,d、q轴电流环之间实现了解耦。控制器输出的u_d_ref和u_q_ref经过反Park变换得到三相调制波,输入PWM发生器。PWM频率设2kHz,直流母线电压1200V,载波幅值按正弦调制规则设置。
PI参数初调我用工程近似法:电流环的开环截止频率设为150Hz,按最优阻尼比0.707计算,比例系数Kp约等于(L × ω_c),积分系数Ki约等于(R × ω_c)。L取0.3mH,R取5mΩ,算出来Kp在0.28左右,Ki在4.7左右。把这两个数填进去跑起来,再根据仿真波形微调即可。
5. 常见问题与排查技巧实录
5.1 仿真速度极慢,怎么优化
5MW直驱模型的仿真速度是很多人第一次跑时最崩溃的问题。PWM频率2kHz,最小步长必须低到5到10微秒才能正确捕捉开关事件,仿真10秒工况,Simulink可能要跑好几个小时。
我的优化手段有这么几条:一是把IGBT开关器件改成平均值模型或者理想开关模型,Simscape里可以选择“Average model”模式,它忽略单个开关动作,只保留调制波的传递特性,仿真步长可以放大到微秒秒级;二是把控制环和电气环分开,用不同采样时间执行,电流环用10微秒步长,转速环和电压环用1毫秒步长,这样大幅减少控制器计算节点数;三是如果只需要验证控制策略,可以直接用平均值模型加简化的母线电压方程替代全功率级电气模型,仿真速度能提升好几个数量级。
5.2 母线电压出现低频振荡,怎么定位
母线电压低频振荡通常表现为电压在几百毫秒周期内上下波动,幅度逐步增大。遇到这种情况,先别急着调参数,用Scope分别记录发电机功率、并网功率和母线电压波形,看功率不平衡的相位关系。
最常见的原因是电压外环带宽和机侧MPPT环带宽太接近,两个环相互作用。比如电压环带宽20Hz,转速环带宽5Hz,如果转速环响应过快,机侧功率快速跳变,电压环又急于修正,就会形成闭环振荡。我的解决方法是把转速环带宽降到1Hz以下,电压环带宽保持15到20Hz,两个环拉开差距,振荡通常就消失了。
另一种可能性是直流母线电容等效串联电阻设置过大,导致控制系统看到的母线电压含有大量纹波,反馈量失真。检查一下ESR,控制在0.1mΩ以内。
5.3 电磁转矩的模型发散问题
如果仿真一启动,电磁转矩就出现大幅高频振荡,甚至数值达到额定值的几十倍,大概率是电流环的PI参数过冲或者解耦前馈算错了。先冻结转速,给定阶跃iq指令,看电流响应曲线。电流响应如果出现严重的超调振荡,把Kp调小一倍,Ki保持不变。如果电流跟踪正常但转矩仍然振荡,检查反Park变换的角度对不对,角度错了一个编码器脉冲,转矩就会周期性波动。
另外一种情况是初始工况设置不匹配。并网前发电机转速为零,磁链初始值也是零,如果直接启动满风速仿真,相当于给电机一个巨大的阶跃冲击。建议先把风速降到切入风速附近跑通初始段,再用信号发生器逐步增大风速,让系统平滑切入最大功率跟踪。
5.4 仿真结果与理论值偏差大
如果稳态功率和理论值偏差超过5%,建议从以下几个方面检查:
- 电机参数填入是否有单位错误,特别是磁链和电感;
- Cp函数是否在额定工况点输出最大值,扫一遍工作区间;
- 直流母线电压是否稳定在1200V,如果电压偏低,变流器调制深度可能进入了非线性区;
- PWM调制方式是否正确,过调制会导致输出电压畸变。
我做一个简单的自检清单,跑任何新模型前都过一遍:
| 检查项目 | 预期结果 | 常见错误 |
|---|---|---|
| 空载反电动势 | 频率对应转速,幅值合理 | 磁链参数错误、极对数错误 |
| 额定风速机械功率 | 约5MW | 扫风面积或Cp公式错误 |
| 母线电压稳态波动 | 小于2% | 电容偏小、电压环PI不合适 |
| 发电机电磁转矩 | 与风力机转矩平衡 | 电流环方向反了、dq轴映射错误 |
上手实操的小技巧分享
最后分享一个我自己的经验:拿到这类模型,先别急着看控制策略怎么写的,第一件事是把模型从头到尾按信号流路径点一遍,搞清楚每个模块的输入输出关系。然后在稳态工况下先把机侧回路单独验证,确认发电机转速和转矩匹配、MPPT输出正确,再把并网侧加进来调母线电压。分段验证方式看着慢,但实际问题排查起来特别快,比整体模型跑崩了再慢慢找原因高效很多。
这套5MW永磁直驱直流并网模型搭完之后,后续可以扩的方向还挺多的:加入低电压穿越控制逻辑、故障工况模拟、多机并联协调控制,都是比较成熟的研究方向。如果后面有时间,我再单独写一篇低电压穿越和故障仿真的实现细节。