这段时间一直在折腾一个5MW永磁直驱风电的直流并网仿真模型核心目标就一句话把风轮捕获的机械功率通过永磁同步发电机和全功率变流器稳定地送进一条1200V直流母线。模型涉及的环节不少——风力机空气动力学、PMSG的dq模型、机侧PWM整流器、MPPT、双闭环PI控制、SVPWM调制最后还得保证直流侧功率流平稳。这套东西放Simulink里搭起来不算难但真正跑出稳定不振荡的波形还是有不少门道。这篇文章把我搭模型时的一些设计取舍、参数计算和踩坑记录整理出来希望能给正在做风电直流并网仿真的同行省点时间。1. 模型整体设计与控制思路1.1 为什么是永磁直驱加直流并网永磁直驱在风电领域已经很主流了——没有齿轮箱传动链短可靠性高低风速段效率也好。尤其是直驱结构的发电机转速很低极对数做得多同步电感相对较大配合全功率变流器控制上很灵活。相比双馈异步机组PMSG天然适合做“全功率变频隔离”机侧和网侧用直流母线解耦两边各自控制就行。这次模型选1200V直流并网而不是传统的690V交流并网主要考虑的是一种直流汇集场景。现在海上风电、陆上大基地经常会遇到直流汇集母线、中压直流微网这类拓扑风机整流之后直接挂在直流母线上省掉网侧逆变器那一大套。直流并网的好处很明显机组之间不需要严格的交流同步控制上只剩机侧变流器这一级效率更高系统也简化了。但代价是直流母线电压必须稳得住否则功率怎么送出去都成问题。我在模型里把直流电网处理成1200V恒压源这个思路在验证阶段是合理的先把发电侧和变流器的行为搞对再考虑直流电缆阻抗、母线电压动态这些细节。1.2 模型的四个子系统整个模型我习惯分四块搭每个子系统职责单一排错的时候不至于满地图找问题风机子系统输入是风速信号和发电机转速输出是机械转矩。核心是那个风能利用系数Cp(λ,β)的计算。发电机子系统用Simulink里的永磁同步电机模块机械端接收风机转矩电端口接变流器。机侧变流器子系统三相两电平PWM整流器把发电机输出的交流电变成直流电挂到1200V母线上。控制子系统这是最核心的包含MPPT转速计算、转速外环、电流内环、坐标变换、SVPWM调制。数据流向也很直观风速进来风力机算功率带动PMSG旋转机侧变流器控制发电机的电磁转矩让风机始终工作在最佳叶尖速比附近直流侧输出功率自然就等于风机捕获的功率减掉损耗。1.3 控制思路是怎么分层的控制上用的是经典的双闭环矢量控制分层很清晰。外环是转速环目标值由MPPT给出也就是根据实时风速算出最优转速外环的输出是电磁转矩参考值再折算成功率电流指令。内环是dq轴电流环d轴电流控制为零q轴电流控制电磁转矩。电流环的输出是dq轴电压指令经过Park反变换和SVPWM生成六路PWM波去驱动IGBT。这里要特别说明一下直流母线电压在这个模型里是由外部恒压源“箝位”的所以机侧只负责功率传输不需要管电压。如果后面要做成真正独立的直流并网接口就需要加一级DC-DC变换器来专门管母线电压这个我放到文章后面讲扩展方向的时候再说。系统层级上电压控制是靠并网侧或直流源提供的电流和转速控制由机侧完成这个分工在风电直流并网里是主流做法。2. 关键参数解算与准备2.1 风机和发电机的额定参数参数是模型的骨架不能随便写。我参考了常见5WM级直驱机型的量级做了一组可以和模型直接对应的参数。叶轮半径按63米算扫风面积约12470平方米额定风速取11.5米/秒左右。风机侧的关键参数我列在下面参数数值说明额定功率5 MW仿真验证目标叶轮半径63 m典型5MW级切入风速3 m/s低于这个转速风机不转额定风速11.5 m/s对应输出接近额定功率最佳叶尖速比8.5用于MPPT最大风能利用系数约0.48桨距角为0时的峰值发电机参数是按1200V直流母线反推的。三相PWM整流器本身有升压能力交流侧线电压有效值取690VSVPWM调制比约0.82整流后直流电压可以稳在1200V。具体参数如下参数数值说明额定功率5 MW与风机匹配极对数60低速直驱特征定子电阻8 mΩ极低大电流低电压机组定子电感Ld/Lq0.7 mH表贴式PMSGLd≈Lq转子磁链约5.8 Wb和大功率直驱电机匹配额定转速约15 rpm对应最优叶尖速比转速反推一下λ_opt8.5v11.5m/sR63m机械角速度ωrλ_opt·v/R≈1.55rad/s也就是约14.8rpm和15rpm对上了。电气频率等于极对数乘以机械角速度60×1.55≈93rad/s约14.8Hz属于典型的直驱低频工况。2.2 直流母线电容怎么选直流母线电容的取值经常被随手一填但实际它对直流电压纹波和动态响应影响很大。粗算公式可以用能量法在负载阶跃或者功率波动时电容要能维持母线电压不掉太多。工程上模拟用简单方式估算按额定功率5MW、直流电压1200V允许电压波动1%即12V开关频率2500Hz来算电容大概需要几个毫法的量级。我在模型里取的是6mF仿真结果直流侧电压波动在允许范围内动态响应也还能接受。启动前记得给电容设置初始电压1200V否则仿真一开始会出现很大的充电冲击电流直接把波形弄花。这个小细节经常有人忽略但非常影响仿真体验。2.3 MPPT参考转速的推导MPPT有很多种做法爬山搜索、功率信号反馈、最优叶尖速比法。这个模型风速是已知输入所以最优叶尖速比法最直接不需要做扰动搜索也不用查功率曲线直接算出参考转速。公式很简单ω_ref λ_opt · v / R。风机的叶尖线速度等于λ_opt乘以风速除以半径就是机械角速度。比如风速8米/秒时ω_ref8.5×8/63≈1.079rad/s也就是10.3rpm。风速10米/秒时参考转速变成约12.9rpm。风速变了参考转速跟着变转速环再去追这个参考值风机就能一直停在最大功率点附近。但要注意一点这个控制策略只适合额定风速以下。超过额定风速以后叶尖速比要保持最优的话转速会超过额定值而且功率会超发这时候必须切到桨距角限功率。我的模型把工作范围限定在额定风速以下桨距角固定为0这是简化处理符合大部分变流器仿真验证的需求。2.4 坐标变换最容易栽的坑坐标变换看着简单实际中出BUG的概率很高。用Park变换把三相电流转到dq轴时必须用转子的电角度而且转子的初始角度要跟模型里永磁同步电机模块的设置一致。Simulink的PMSG模块出来的转子角度通常默认是d轴和A相轴线的夹角。如果你用D轴对齐转子磁链的约定那标准Park变换用cos/sin函数就行但有些文档里用的是带初始偏移的变换直接套标准公式就会导致d轴q轴解耦失败电流内环完全没法工作。我的处理办法是直接用模块自带的转子角度输出接进Park变换不要自己额外加偏移如果电流波形不对就检查旋转方向——仿真里PMSG转速如果是负的角度方向也要跟着反过来。坐标变换的代码建议用函数块封装方便统一改符号不要散落在各个Fcn块里排查起来太痛苦。3. Simulink建模实操3.1 风力机转矩模块风力机模块可以用Simscape自带的风机模型也可以自己写。我更喜欢自己写一个MATLAB Function块因为参数可控公式透明。核心是这三步先算叶尖速比再算风能利用系数最后算机械转矩。具体代码长这样function Tm wind_turbine(v, wm) R 63; % 叶轮半径 rho 1.225; % 空气密度 lambda_opt 8.5; beta 0; % 桨距角固定为0 if v 0.1 v 0.1; % 防除零 end lambda wm * R / v; % 标准Cp曲线公式 lambda_i 1 / (1/(lambda 0.08*beta) - 0.035/(beta^3 1)); Cp 0.5176 * (116*lambda_i - 0.4*beta - 5) * exp(-21*lambda_i) 0.0068*lambda; % 机械功率和转矩 Pm 0.5 * rho * pi * R^2 * v^3 * Cp; Tm Pm / wm; end这段代码里有个细节lambda_i的计算当叶尖速比很小时公式会出现奇异点加了小风速保护依然可能出现lambda0时除零的问题所以建议在lambda小于0.1时直接给Cp赋值0。仿真跑的时候如果看到转矩跳成NaN基本就是这里的问题。实际运行时风机还有一个转速下限风速太低时风轮不应转起来这可以在控制器里加判断如果风速小于切入风速直接让电磁转矩参考值为零发电机不发电。3.2 机侧变流器双闭环控制机侧控制是模型的核心。我搭控制子系统时严格分了三层逻辑。第一层是MPPT层输入风速输出转速参考值这个上面已经说了。第二层是转速外环。转速误差进PI控制器输出电磁转矩参考值。注意转矩参考值要限幅不能让转矩超过发电机的极限。我设的上限是3.5e6牛米略大于额定转矩。转矩参考值再换算成功率电流指令公式是iq_ref 2 * Tref / (3 * p * ψf)三相合成转矩公式里的系数要看电流是用峰值还是有效值我这里用的是峰值坐标系下的标准写法。把额定转矩3.18e6N·m代进去算出额定q轴电流约在5500A量级和前面估算的定子相电流峰值基本一致。第三层是电流内环。d轴电流指令为零q轴电流指令来自转速环输出。电流误差经过PI调节再叠加前馈解耦项vd -ωe * Lq * iq PI_d vq ωe * Ld * id ωe * ψf PI_q这里前馈补偿很重要。如果不加解耦d轴和q轴之间会互相影响电流响应会变得很慢甚至出现振荡。我试过把前馈去掉只靠PI硬扛结果是转速一变化电流就开始发散。3.3 调制与仿真步长PWM调制可以直接用Simulink的PWM Generator模块选二电平模式载波频率取2500Hz。很多教程喜欢自己搭三角波比较器其实没必要模块里还能直接设置死区时间比自己写方便得多。SVPWM比SPWM多了约15%的电压利用率这个模型直流1200V下SPWM最多输出约734V交流线电压SVPWM可以到约844V。690V的线电压在SVPWM下调制比约0.82留了裕量这样仿真中不容易出现调制饱和。如果发现调制波满头就要检查直流电压是不是不够或者交流侧参数是不是偏离了690V太多。仿真步长方面我用的求解器是ode23tb最大步长设20e-6秒。开关频率2500Hz载波周期400微秒如果步长太大PWM边沿会失真电流纹波和实际不符控制器采样到的电流也会有台阶。有些教程让你把步长设为1e-6那仿真会慢到让人崩溃其实20微秒在工程验证层面足够前提是求解器别用ode45碰到电力电子这种刚性系统会非常痛苦。3.4 PI参数整定的工程方法PI参数是新手最容易卡住的地方我给一套能直接抄的初值和方法。电流内环被控对象近似为一个电阻加电感的惯性环节G(s)1/(LsR)。带宽法整定选择电流环带宽ωc约为1000rad/s约160Hz则Kp ωc · L 1000 × 0.0007 0.7 Ki ωc · R 1000 × 0.008 8这只是纯被控对象算出来的理想值实际仿真中因为存在采样延迟和PWM滞后我直接把Kp调到了2Ki调到了40响应速度更快稍微有一些超调但可控。电流环的调试标准是dq电流指令阶跃时实际电流能在1~2毫秒内跟上。转速外环因为风力机的转动惯量太大风力机本身是慢系统转速环带宽只取电流环的1/10甚至更低。我在自然单位制下用了Kp10000Ki1000效果是转速响应没有明显振荡。要知道这个数值换一台机组就不一定能用因为转动惯量和磁链参数都变了。更稳妥的办法是先设Kp为0只加Ki慢慢往上加积分看到转速开始缓慢追参考值后再加Kp消除超调。这个顺序能少走很多弯路。4. 仿真结果怎么看4.1 启动和风速阶跃场景我设计了一个三段式仿真来验证模型性能0到5秒风速固定在8米/秒10秒时风速阶跃到10米/秒18秒时风速降到6米/秒。整个仿真时长25秒。8米/秒风速下理论上风轮捕获功率是P 0.5 × ρ × A × v³ × Cp 0.5 × 1.225 × 12469 × 512 × 0.48 ≈ 1.87MW模型跑稳后直流侧功率应该在1.8~1.9MW之间考虑到PMSG和变流器的损耗会有几kW的出入。10米/秒时理论功率约3.66MW直流输出功率相应提升。6米/秒时理论功率约0.79MW。这个对照非常重要——模型输出必须跟理论值保持在合理误差内否则说明MPPT或者参数设置有系统性问题。仿真启动时我做了两个处理直流电容初始电压设1200V避免启动冲击发电机初始转速设成风速8米/秒对应的最优转速附近避免经历一个特别长的启动爬坡过程。如果你不设初始转速仿真前几秒转速要从零爬上去MPPT参考值又一直比实际转速大转速外环会长时间饱和波形看起来会非常难看。4.2 从波形判断模型是否正常看波形不是光看功率多大要按下面几条逐项检查转速波形应该平滑地追着MPPT参考值走8米/秒时参考转速约10.3rpm阶跃到10米/秒后目标变成约12.9rpm转速环响应大概在2~3秒内稳定。如果转速出现持续振荡多半是转速环PI增益偏大或者风速信号没滤波。d轴电流应该在0附近小幅波动幅值不超过几十安培q轴电流应该随着功率变化有明显改变10米/秒时q轴电流应该比8米/秒时大接近一倍。如果d轴电流明显偏置检查坐标变换角度是否准了。相电流波形应该是光滑的正弦叠加高频纹波畸变太大就看开关频率设置和步长。直流侧电流是判断功率流的关键。1200V母线1.87MW对应约1560A直流电流3.66MW对应约3050A。直流电流跟着风速走稳定后纹波应该在几十安以内。如果直流电流振荡回头看直流电容和电流环参数。有功传输方向是正的——从风机往直流电网送电电流方向反了说明模型里电源极性接错。5. 调试实录与避坑清单5.1 四个高频问题仿真跑着跑着发散是常态我遇到过的问题基本可以归成下面四类。第一类是启动就报错提示初始化失败。这是PMSG的转速初始值或者变流器状态没设对。解决方法是给电机模块设置正确的初始机械转速给直流电容设初始电压再不行就把仿真最大步长调小一点先让模型跑通再调参数。第二类是电流环振荡。现象是dq轴电流大幅波动相电流畸形。九十的概率是PI增益有问题。我调试这个模型时把电流环Kp从理论值0.7调到2Ki从8调到40马上稳了。如果增大增益反而振荡那还要看电感参数和实际模型是否一致——SimuSCAPE里的永磁同步电机模块定子电感是基于相电压定义的和dq模型里用的电感有时候差个系数要对着文档换算清楚。第三类是转速追不上MPPT参考值。风速阶跃后转速响应特别慢像拉不动的样子。这个基本是转速环PI增益太小或者转矩限幅卡住了。把转速环Ki往上拉一到两个数量级看到转速缓慢追上去再加Kp。第四类是仿真速度慢到无法忍受。这个要从两个角度解决一是用离散化PWM载波周期设为开关频率的倒数仿真器会自动用固定步长处理开关边沿二是减少波形记录点示波器存成千上万个点对内存消耗很大。还有一个很有效的技巧把风力机和变流器的损耗忽略掉模型更简单跑得更快。我把这些问题整理成了速查表现象可能原因处理办法启动初始化失败电机初始转速未设置设置机械转速初始值接近MPPT参考值启动电流冲击过大直流电容初值未设电容初始电压设1200Vdq电流振荡电流环PI过大或电感参数换算错误降低Kp核对电感定义转速响应过慢转速环KI偏小增大Ki再适当增加Kp直流电流纹波过大电容偏小或开关频率偏低加大电容到6~10mF提高载频仿真极慢求解器选择不当用ode23tb最大步长20微秒5.2 调参顺序和我的习惯调参这个事顺序对了比什么都重要。我的习惯是从内到外调先调电流环再调转速环最后才看功率外特性。电流环调试的时候先给固定的iq参考值不接转速环直接看电流响应曲线能不能跟上。这里可以用Signal Builder或者一个阶跃信号代替转速环输出等电流环调好了再闭合转速环。转速环闭合以后先给一个风速阶跃观察转速的超调和振荡最后再上完整的MPPT逻辑。这个流程看起来慢实际上最省时间。我见过很多人从外环开始调结果电流内环根本没稳定外环再怎么调都是白费。还有一个小习惯是把风速信号在进入MPPT和风机模块之前做一个一阶惯性滤波时间常数取1秒左右。真实风速里有很多湍流成分虽然模型里是阶跃变化但加一个滤波能模拟出风场的平滑特性也避免MPPT参考转速剧烈跳动。我在模型里用的就是1秒的一阶滤波器效果很稳。如果你想把模型往深了做还有几个现成的扩展方向一是加一级Boost升压变换器让发电机侧变流器只管MPPT由Boost专门控制1200V母线电压这样就能脱离恒压源假设二是把直流母线换成带电阻电感的直流电缆模型做直流侧稳定性分析三是把单机扩展成多机并联研究直流汇集风场里的功率分配和电压下垂控制。这几个方向我后续都会继续折腾到时候再把心得整理出来分享。
