ARTICLE DETAIL

资讯详情

深耕网站建设与运营推广的一线实战洞察。

5MW直驱永磁风机并网仿真:从气动到电网的建模与调参实践

5MW直驱永磁风机并网仿真:从气动到电网的建模与调参实践 把一台5MW的直驱永磁海上风机整机并网模型跑通最熬人的其实不是控制算法本身而是把气动、机械、电机、电力电子、电网这五个层面的模型拼在一起还能稳定仿真。这套模型我在Matlab/Simulink上搭了将近三个月直流母线电压按1200V设计机侧整流加网侧逆变全功率变流最终跑出了从风速阶跃到电网电压跌落的完整工况。如果你也在做风电并网仿真或者刚接触Simulink里的电力电子建模这篇笔记应该能帮你少走几个弯路。下面我把整个建模思路、控制策略、调参过程和平时的坑全部摊开讲。1. 为什么选这套拓扑直驱PMSG与1200V母线的前因后果1.1 直驱路线和双馈路线的取舍海上风电的机组选型最核心的矛盾是可靠性和可维护性。双馈机组用齿轮箱加转子侧部分功率变流器变流器容量只需要整机功率的三成左右成本优势明显但齿轮箱在海上高盐雾、高湿环境里是故障高发点一旦坏了更换成本高到让人头疼。直驱永磁同步发电机PMSG把齿轮箱整个去掉传动链缩短机械故障点大幅减少而且永磁励磁不需要转子励磁电流低风速段效率更高全功率变流器还能从结构上把电网侧的故障与发电机侧隔离开来。对5MW这个档位来说直驱PMSG已经是非常成熟的路线既有大量工程数据可以参考算力上又不像10MW以上机组那么吃紧非常适合在普通工作站上做整机仿真验证。1.2 1200V直流母线怎么来的以及它的边界条件1200V这个母线电压不是拍脑袋定的。背靠背全功率变流器要正常向690V交流系统馈电直流母线电压必须高于交流线电压的峰值。690V线电压峰值大约是975V留出10%到20%的调制裕量直流母线取1100到1200V是比较合理的区间。所以我建模时直接把直流母线参考值设为1200V对应的交流侧就按690V设计。这里必须提醒一句5MW在1200V直流母线下的电流接近4.2kA单个IGBT模块在这个电流等级下已经非常吃紧。实际工程往往会走中压方案或者多个功率单元并联仿真模型里出于控制策略验证的目的按单机全功率变流器抽象是没问题的电压等级只影响标幺基准控制结构完全一致。用标幺值系统建模还有一个附带好处把电机的伏安单位和母线电容等参数归一化之后数值刚度会明显下降仿真更容易稳定。1.3 Simulink顶层架构五个子系统的分工我的模型分成五个子系统风速与气动模块、机械轴系、PMSG电机本体、背靠背变流器机侧整流加直流母线加网侧逆变、LCL滤波与三相电网。整体架构建议先用平均值模型把控制逻辑跑稳再用开关模型验证波形细节。平均值模型里的开关桥臂用受控电压源和受控电流源代替没有PWM纹波仿真速度快、不挑求解器控制器调参阶段非常好用。开关模型用Simscape Electrical或传统电力电子库都能搭步长小、仿真慢但能看到THD、并网冲击电流这些平均值模型看不到的东西。我踩过的坑是如果一上来就上开关模型PI参数没调好时波形全是毛刺根本分不清是控制器问题还是数值问题。分层推进是排坑最快的方法。2. 风轮与PMSG模型三层物理的递进2.1 风速模型从恒风到湍流的场景选择风速模型直接决定你要验证什么。恒风速只能看稳态工作点验证MPPT动态性能必须给阶跃风速而想体现海上风况的随机性就得用湍流风谱。我仿真里常用组合风基本风8m/s叠加一个持续6s、幅值3m/s的阵风再叠加一段渐变风用来验证转速环的跟踪能力和变桨逻辑。做并网动态测试时我反而不太用强湍流因为风轮惯量很大湍流中的高频分量传到电功率侧时已经被机械系统的大惯性滤掉了主要影响在低频段。海上平均风速高、湍流强度通常低于10%这也解释了为什么海上机型普遍把叶轮做大、功率密度做低——风资源更平稳就能在低风速段多捞能量。2.2 Cp(λ,β)气动模型与最大功率跟踪气动模型的核心是功率系数函数Cp(λ,β)。我在Matlab里用的是一组经典经验公式直接写成函数可以到处复用function Cp cp_5MW(lambda, beta) % 经典风力机Cp经验公式lambda为叶尖速比beta为桨距角(rad) c [0.5176 116 0.4 5 21 0.0068]; if lambda 0 Cp 0; return; end lambda_i 1 / (1/(lambda 0.08*beta) - 0.035/(beta^3 1)); Cp c(1)*(c(2)/lambda_i - c(3)*beta - c(4))*exp(-c(5)/lambda_i) c(6)*lambda; end叶尖速比定义为λ ωR/v其中ω是风轮转速R是叶轮半径v是风速。额定风速以下保持桨距角为0只要让转速跟随风速变化维持最优叶尖速比就能让Cp停在最大值附近。最大功率跟踪MPT的给定转速有两种做法一是按λ_opt直接实时计算ω_ref二是预存功率-转速查找表。我习惯用前者参数可调性强代码就是一行ω_ref λ_opt * v / R。如果用查找表记得给切入风以下加最小转速保护否则低速区转速给定会趋近于0电流环跟着算出来的转矩给定会非常离谱。2.3 轴系直驱也要用两质量块模型直驱没有齿轮箱很多人就顺手把轴系当成单质量块。实际上风轮和发电机转子之间依然存在柔性连接扭振模态频率一般在1到10Hz之间电网扰动或者风速剧变时这个模态会被激发。单质量块模型看不到这个过程等于把机械动态中最有价值的一段信息丢失了。两质量块模型用轴刚度K_shaft和阻尼D_shaft把风轮惯量和发电机惯量连起来。风轮惯量比发电机转子惯量大好几个数量级两者之间的速度差会产生轴系扭振。验证轴系参数是否合理有个小技巧在轴系上故意加一个脉冲扰动转矩观察发电机转速曲线有没有出现明显的低频振荡如果振幅不衰减说明阻尼参数设小了需要加大D_shaft。仿真里把这条曲线录下来对后面解释电网故障下的功率振荡非常有帮助。2.4 PMSG的dq方程与id0控制的适用边界永磁同步电机在dq旋转坐标系下的电压方程是ud R·id Ld·d(id)/dt - ωe·Lq·iquq R·iq Lq·d(iq)/dt ωe·(Ld·id ψf)电磁转矩Te 1.5·p·[ψf·iq (Ld - Lq)·id·iq]。表贴式PMSG的Ld等于Lq磁阻转矩为零用id0控制最省事每安培电流的转矩产出最大。内置式PMSG因为有磁阻转矩要做最大转矩电流比MTPA控制逻辑会更复杂一些。我的模型用的是表贴式参数定子电阻0.002Ωd轴和q轴电感都取0.35mH极对数84额定转速14.5rpm对应的电频率约20Hz。这里最容易翻车的点是电角速度和机械转速之间的极对数换算一不留神就会让转矩方向反掉或者弱磁曲线彻底乱套。建议把PMSG模型单独封装成子系统内部统一用电角度计算对外接口用机械角度这样上层控制策略不用时刻惦记极对数。3. 机侧变流器转速环、电流环与变桨配合3.1 级联控制结构转速外环加电流内环机侧整流器的任务有两个把PMSG发出的有功功率稳定送到直流母线以及通过控制电磁转矩改变转速让机组运行在最大功率点。控制上采用转速外环加电流内环的级联结构。转速环输入是来自MPT模块的ω_ref输出是q轴电流参考iq_ref电流环在dq轴分别做PI控制d轴电流参考id_ref设为0。dq轴电流之间存在耦合项需要加前馈解耦。解耦不是不加也能跑而是加了之后电流环带宽可以更接近物理极限动态响应更好。机侧电流环的PI参数可以按带宽法起步电流环带宽取开关频率的1/10到1/5Kpc Ls·ω_bwKic Rs·ω_bw。按2kHz开关频率、500Hz电流环带宽算下来Kpc大约是1.1Kic大约是6.3。这个组合先跑起来再根据波形微调基本不会出大问题。3.2 抗积分饱和风电并网调参的头号避坑点转速环的输出就是iq_ref它必须限幅一般取1.1倍额定电流。大工况下转速环非常容易饱和比如风速从8m/s阶跃到12m/s转速误差突然拉大PI积分项疯狂累加一旦深入饱和区等误差反过来时积分器还要慢慢释放于是转速超调严重恢复过程拖得很长极端情况下直接发散。处理办法是限幅逻辑和积分项联动当实际输出不等于限幅前输出时冻结积分或者用反计算法把超出限幅的量按倍数反馈回积分器输入端。这个逻辑在Simulink里用几个基本运算模块就能搭但在风电并网仿真里它就是你调半天都稳不下来的头号嫌疑犯。我看到很多初学者的模型PI参数看着没问题波形就是稳不住一问限幅做了但没有抗积分饱和。3.3 额定风速以上的变桨配合额定风速以下桨距角保持0转速环自动跟踪MPT曲线。超过额定风速后转速给定锁定在额定值不再上升桨距角环开始接管功率限制。桨距角环同样用PI控制转速误差经过PI输出β_ref再经过一个一阶惯性环节模拟变桨执行机构时间常数一般取0.02到0.2秒。取太小变桨动作生硬功率波形毛糙取太大高风速下转速容易超调。我给变桨执行器加了一阶惯性加限速限速5°/s限幅0到30°。30°以上气动功率已经减得差不多但安全位置仍然要给到90°。仿真里变桨和转矩的配合时序值得反复试先增大桨距角降低气动转矩再让转速环跟进电流给定。这个先后顺序对功率波动幅度影响很大顺序反了母线电压会出现明显的凹陷。4. 网侧变流器与1200V直流母线功率出去的前提是母线稳住4.1 母线电容取值与预充电环节直流母线是背靠背系统的枢纽机侧往母线注入功率网侧从母线抽出功率瞬时不平衡由母线电容吸收。电容越大电压纹波越小、抗扰动越强但体积成本都上去了。工程经验公式是C ≈ P / (2ω_grid·ΔV·Vdc)按5MW、1200V、允许5%纹波和50Hz电网计算理论值会到0.11F级别。实际电压环自己能抑制一部分波动折中取30mF到50mF起步先看动态波形再调。我模型里初始设30mF电网电压跌落工况下观察母线过压幅度不够再往上加。预充电这个点我这次特别想提醒。如果模型里母线电容初始电压设成0然后直接合闸启动你会看到母线电流冲到几十倍额定值控制器疯狂饱和误以为是参数问题其实只是初始条件不对。物理系统里通常会先通过预充电电阻把母线电容充到接近1200V再启动机侧。仿真里两招都能用一是直接给电容的初始电压条件写1200V这样最省事二是搭一个预充电电阻加旁路开关的序列更接近物理过程但要额外处理开关时序。新手建议先选第一种把控制逻辑调稳了再回来补细节。4.2 PLL与并网同步合闸冲击的根源网侧逆变器要并网第一件事是锁住电网相位。同步旋转坐标系锁相环SRF-PLL的原理是把q轴电压控到0让d轴和电网电压同相。PLL的PI带宽一般取10到30Hz带宽太高会把电网谐波引进来导致相位抖动。并网断路器合闸的条件是逆变器输出电压的幅值、频率、相位都要接近电网。仿真里最容易犯的错是PLL初始相角设0但电网电压源本身已经有相角合闸瞬间相位差大冲击电流非常吓人。正确做法是让PLL先运行一段时间完成锁定或者直接把逆变器侧的电压源初始相角也设成匹配值。我在调试时习惯记录PLL输出的q轴电压如果合闸前q轴电压一直不归零绝不闭合并网开关。4.3 电压外环加电流内环有功无功怎么协调网侧变流器控制的核心是稳定Vdc。直流电压环把母线电压误差经过PI变成d轴有功电流参考q轴电流参考决定无功功率井网友好型通常设0。电流内环和机侧类似但因为控制对象是带LCL滤波器的逆变器电流环带宽宜取低一些250到350Hz比较合适。控制器输出经过SVPWM送到IGBT桥臂。联调时最容易忽略的问题是坐标系约定。机侧、网侧、PLL各自都有旋转坐标系d/q轴方向必须统一。我吃过这个亏两套控制的iq符号相反结果机侧往母线送功率网侧又从母线抽回来母线电压呈现出周期性的充放电锯齿差点被误导去调母线电容。先花半小时把坐标方向理清比调试三天都有效。4.4 LCL滤波器参数设计与有源阻尼LCL比单L滤波器高频衰减能力强但它本身是欠阻尼二阶系统在谐振频率处有尖峰。设计法则是让谐振频率f_res落在0.3到0.5倍开关频率之间既保证高频衰减又避开控制带宽。我用的参数是L1取0.15pu约0.45mHL2取0.10pu约0.30mHCf取0.03pu约280uF谐振频率约700Hz在2kHz开关频率下大约0.35倍位置比较合适。谐振尖峰必须处理。最省事的办法是给滤波电容串联无源阻尼电阻缺点是持续损耗更推荐电容电流反馈有源阻尼——在电流环输出里减去kc倍的电容电流等效增加阻尼却不产生实打实的线路损耗。Simulink里就是多一路测量信号乘一个系数效果立竿见影。4.5 电网电压跌落时的行为LVRT本质是功率平衡电网跌落后逆变器输出电流能力受限机侧还在持续输入功率不平衡功率会灌入直流母线导致电压过冲。所以低电压穿越LVRT逻辑的本质就是功率平衡管理检测到电压跌落立刻限制有功电流同时机侧限功率或者投入Chopper卸荷电路。简化模型里我用一个受控电阻模拟ChopperVdc超过阈值就投入把多余能量泄放。这次只做了对称跌落的验证工况是0.1s跌落到70%持续0.6s再恢复。非对称跌落还有负序抑制、正负序分离控制这是下一步要补的但先把对称跌落下的功率管理跑稳母线电压不超1300V就已经覆盖了这个直流母线等级下最核心的安全约束。5. 调试现场那些让模型从发散到稳定的关键操作5.1 模型第一次跑就发散先查什么仿真第一次运行就数值爆炸90%的情况和PI参数无关而是初始条件矛盾。检查顺序我固定是四步母线电容初始电压有没有设到1200VPMSG初始转速有没有和风速匹配PLL有没有输出有效相位断路器初始状态是开还是合。我踩得最深的一个坑是PMSG初始机械转速设成0MPT模块却给了一个接近额定的转速参考。转速环误差瞬间拉到最大积分饱和还没做电流内环也跟着饱和模型在第一个仿真步长里就崩了。解决方法是把电机初速预设到额定转速附近再让转速环自己收敛比从0硬拉要稳得多。记住这个教训之后我每次新建模型都会先把所有初始条件列一张表检查一遍。5.2 PI参数起步法一套可以直接抄的基准值给一套可以复现的试凑参数标幺制下机侧电流环Kp1.0、Ki6.3转速环Kp0.2、Ki0.5网侧电流环Kp0.5、Ki5电压环Kp2.0、Ki20。这套组合在多数5MW级模型里能跑稳再按带宽公式微调就行。调参时固定一个环、只看另一个环的动态是定位振荡源头的基本功。把转速给定改成阶跃先断开电流环饱和因素只观察电流环跟踪效果再恢复观察转速环响应。更关键的是机侧和网侧分开测试机侧接电阻负载网侧接理想直流源分别验证系统最后再联调。两边同时振荡的时候你根本不知道是谁在惹事分开测是唯一高效的路。5.3 求解器与步长决定仿真结果可信度开关模型建议用离散定步长求解器步长不超过开关周期的1/20到1/30。2kHz开关频率对应0.5ms周期步长取10到15微秒够用取1微秒能看更干净的波形但仿真速度太慢。平均值模型可以用ode45或ode23tb变步长即可。还有一个容易被忽视的细节SVPWM调制波和PWM发生器如果不在同一个离散任务里采样时刻对不上会额外引入数字谐波。建议把PWM生成固定在某个离散周期比如1e-5s其他连续部分用零阶保持器对接录出来的波形才可信。5.4 典型故障现象的快速对照现象常见原因处理方向母线电压持续缓慢上升网侧有功送出不足或电压环增益太小加大电压环Kp检查网侧电流限幅并网瞬间电流尖峰PLL相位未同步或断路器相角不匹配查看PLL q轴电压等待锁定后再合闸转速低频振荡轴系阻尼太小或转速环带宽过高增大D_shaft降低转速环带宽电流波形含强烈谐振尖峰LCL没有阻尼或谐振频率接近控制带宽加有源阻尼检查f_res位置波形高频毛刺明显求解器步长过大或控制采样率过低缩小步长统一离散任务这张表是我调试时最常翻的一张基本能覆盖九成以上的发散去向。5.5 波形评估指标什么叫并网质量合格定性判断之外需要一组定量指标。电流THD用Simulink的FFT分析工具看LCL参数合理时THD应低于5%母线电压纹波不超过±5%转速跟踪误差在稳态下不超过0.05pu并网合闸冲击电流峰值不超过额定值的30%。这些指标不是随意拍的它们直接对应并网电能质量的主流工程要求。仿真阶段把这几项录下来后面做实验台或者出报告时都有依据。6. 一组典型工况的波形解读6.1 风速阶跃下的MPPT过程风速从8m/s阶跃到11m/s时转速环应该在3到5秒内把转速平滑升到新的平衡点叶尖速比维持在最优值附近。关键看两个量转速是否有超调电磁转矩是否有瞬时冲击。风轮惯量非常大转速本就不该快速响应阶跃控制目标是平滑地逼近最优叶尖速比而不是像伺服系统那样快速精确跟踪。如果转速超调超过5%我会先降低转速环带宽再检查轴系阻尼。这个工况同时能验证机侧电流环是否饱和如果电流给定瞬时冲到限幅值说明带宽或者给定斜率还要再放慢一点。6.2 并网合闸瞬间的冲击控制PLL锁定、逆变器输出电压与电网同频同相后合闸理想情况冲击电流应该接近0实际允许一小段衰减振荡峰值不超过额定电流的30%。我这次合闸冲击峰值大约15%额定波形在半周波内衰减到稳态这就算合格。合闸成功后母线电压基本不变只有微小跌落这是变流器已经接管功率馈送的直观证据。如果合闸后母线电压出现明显凹陷说明网侧电压环还没来得及接管或者合闸时PLL相位仍有点偏差需要回头查同步逻辑。合闸冲击这个指标是所有并网研究里最容易拿出来对比的参数建议每次改动模型后都重录一次这个工况。6.3 电网电压对称跌落30%的响应0.1s时电网电压跌落到70%持续0.6s恢复。波形上的典型顺序是电网电流瞬间限幅逆变器输出有功功率下降母线电压短暂上升后被Chopper或机侧限功率抑制在1250V以内恢复瞬间电流会有一段无功支撑过程。这段工况验证的不是稳态性能而是功率管理逻辑能不能在毫秒级把机侧和网侧功率重新平衡。如果母线电压过冲超过1300V要么Chopper投入阈值太高或泄放功率不够要么机侧限功率响应太慢。两个方向都要查。做完这个工况基本可以认为这套1200V背靠背系统的并网逻辑闭环了。我个人的体会是Simulink里跑风电模型最怕的不是模型复杂而是把物理假设忘在模型外。5MW直驱、1200V母线这套方案真正有价值的不在那些公式而在你如何用仿真逼问出工程矛盾直流母线电压多高合适、电流大到什么程度不可接受、并网冲击控制在什么水平才算合格。这套模型跑通之后我后续又往里面加了湍流风速功率谱密度分析和全风速段功率曲线验证把控制策略的前期验证基本都搬到了仿真环境里。如果你也在搭类似的模型建议先搭平均值模型把控制逻辑跑稳再切换开关模型看细节遇到问题先查初始条件再查PI限幅最后才去怀疑算法本身。这条路走通之后你会觉得那些发散、振荡、毛刺其实都是模型在跟你说话关键是你有没有耐心听它说完。
返回列表