做风电并网仿真的人大概率都遇到过这么一种尴尬风机模型跑得挺欢逆变器控制也调得挺顺但一到并网那一下系统就跟你翻脸——电压闪变、功率抖动、频率掉坑甚至直接数值发散Simulink整块报红。原因很简单传统跟网型PQ控制本质上就是一个受控电流源电网给什么频率电压它就跟着跑电网一旦偏软它就撑不住。我最近完整做了一个风储VSG项目也就是基于虚拟同步发电机VSG的风储并网系统Simulink仿真从控制原理、参数整定到模型搭建一条龙走下来踩了不少坑也把很多关键细节真正想明白了。这篇就来把整个项目和盘托出给正在做新能源并网、微电网仿真或者VSG相关毕设和课题的朋友一份可以直接抄作业的实操参考。1. 方案选型为什么非得给风储系统配一个“虚拟发电机”1.1 风电并网的核心痛点双馈风机或者直驱永磁风机本身都是电力电子变换器接入电网转子运动系统和电网之间隔着直流母线机械惯量几乎没有传递到电网上。传统控制用锁相环盯住电网相位然后输出指定的有功无功电网频率波动时它基本不参与支撑。比例高起来之后整个电网系统惯量下降得厉害频率变化率RoCoF变大一次调频能力变弱这是老生常谈的问题但却是VSG方案存在的根本理由。我刚开始做这个课题的时候想的是能不能直接在风电机组侧改控制把变流器的PQ控制换成VSG控制。后来发现风机的功率本身受风速影响波动很大如果强行让风机承担虚拟惯量和调频任务功率指令波动和虚拟惯量所需功率叠加在一起直流母线电压根本稳不住。这是一个很容易被新手忽略的核心点VSG要求功率源能够在短时间内提供或吸收额外功率而风机的最大功率跟踪MPPT控制响应速度远不够快。1.2 VSG和传统并网逆变器控制的本质区别先看一个最直观的对比表对比维度传统PQ控制VSG控制外特性受控电流源受控电压源带内阻频率支撑不参与锁相环跟随具备虚拟惯性和一次调频电压支撑无功指令控制无功-电压下垂自动调压并网切换需要预同步处理同样需要但支撑能力更强参数整定PI参数为主惯量J、阻尼D、下垂系数等VSG的本质就是把同步发电机的转子运动方程搬进逆变器控制里。控制回路里人为构造一个“虚拟转子”用有功功率偏差去积分出虚拟角频率再用频率偏差积分出虚拟功角整个电压源的相位由此生成。这样电网频率跌落的时候VSG会因为它“有惯量”而自然地减速并释放存储的能量表现出来就是输出有功功率增加对电网形成支撑。这里有个很重要的地方要说透VSG的“惯量”不是凭空产生的物理飞轮它本质上是控制算法制造的对外动态特性功率的来源最终还是直流侧。所以风储储系统里“储能”不是辅助项而是必不可少的能量缓冲池。我做的项目里采用锂电池储能、双向DC-DC加逆变器的架构风机功率经过整流后和储能双向变换器一起挂在直流母线上再统一经过VSG控制的逆变器并网。这样VSG对外表现出来的惯量、阻尼、调频能力都由储能单元来兜底。1.3 储能系统在VSG中的能量平衡逻辑这就解释了为什么做风储VSG而不是单纯做“风机VSG”。虚拟同步机的惯量响应功率计算公式可以粗略写成ΔP_J ≈ J·ω_n·(dω/dt)频率变化率越大需要注入或吸收的功率就越大。一台2MW级别的风机如果设置惯性时间常数H为4秒在0.5 Hz/s的频率变化率下瞬时功率差额会非常可观风机MPPT不可能提供这种瞬态功率。但储能DC-DC的响应速度在毫秒级完全可以通过双向Buck-Boost快速拉高或降低直流母线功率。我实际搭仿真的时候储能模型用的是简化一阶惯性环节加荷电状态SOC计算没有用复杂的电化学模型。原因是在VSG并网系统级仿真里重点考察的是功率响应和电能质量电池内部细节对结果影响不大反而会增加仿真时间和数值刚性问题。当然如果课题侧重储能优化那就需要换更细的电池模型这是另一个方向。仿真初期用简化模型把控制链路跑通后期再逐步细化是我一直推荐的做法。2. VSG数学模型与参数整定先把公式吃透再动手搭模型2.1 转子运动方程和虚拟惯量参数设计VSG控制的核心公式转子运动方程J·(dω/dt) T_m - T_e - D·(ω - ω_n)功率形式下可以写为J·ω_n·(dω/dt) P_m - P_e - D·ω_n·(ω - ω_n)其中P_m是机械功率输入对应值这里是有功功率参考值PrefP_e是实际电磁功率D是阻尼系数ω_n是额定角频率。控制框图里有功功率参考Pref与实际输出功率Pe的差值经过惯性环节生成角频率ω再积分得到相位角θ。J的物理意义我习惯于用飞轮来理解转动惯量大相当于飞轮重转速变化慢对电网频率波动的抵抗能力强但动态响应慢功率超调大J小响应快但支撑能力弱。合适的J取值原则是让VSG的振荡频率落在合理范围内。VSG的虚拟功角特性近似二阶振荡环节自然振荡角频率可以估算为ω_n_vsg ≈ √(K·P_max / (J·ω_n))其中K是功角灵敏度系数在额定工作点附近近似为电网电压和逆变器输出电压乘积与同步电抗的比值。工程上更直接的做法是根据期望的振荡周期T和阻尼比ξ反推参数。我项目中设J对应的惯性时间常数H为3~6秒切换测试阻尼比取0.7左右对应的D为D 2·ξ·√(J·ω_n·K)这个推导过程说来简单实际调参的时候如果不理解纯靠盲试参数可能一天都试不出稳定波形。我在项目里用MATLAB脚本算了一组初值再做小范围扫描调整效率高很多。2.2 有功-频率下垂与一次调频VSG除了虚拟惯量还必须有稳态下的下垂特性也就是有功-频率下垂控制模拟发电机组的调速器。稳态时频率偏差和功率偏差满足ω ω_n - K·(P_e - P_ref)这里K就是下垂系数单位rad/(s·W)工程上常用%表示比如2%下垂意味着满载功率变化对应2%的频率偏移。并网接口处Simulink里通常把这个关系叠加上去替代转子方程中的阻尼项。实际上阻尼D和下垂K在稳态特性上有重叠初学时容易混淆。我自己的理解阻尼项影响的是动态过程的振荡衰减速度下垂项决定最终的频率-功率静态关系。两者可以统一成一个等效阻尼加死区环节但为了物理意义清晰我在模型里分开实现。仿真结果显示负荷突增2%时VSG频率短暂下降然后恢复到一个比额定值低约0.03 Hz的稳态值这与下垂设计值基本吻合。2.3 无功-电压下垂与励磁调节无功控制回路相对简单模拟同步发电机的励磁调节器表达式E E_0 K_q·(Q_ref - Q_e)E是虚拟励磁电动势Q是无功功率。这里K_q是电压下垂系数取值通常让无功功率从0变化到额定值时电压偏差在2%以内。逆变器输出电压的幅值就由E决定再通过电压电流双环控制最终生成调制波。无功回路有一个非常容易被忽视的坑初始无功参考值必须和实际负载无功匹配否则并网瞬间会有巨大的无功冲击。做预同步的时候理想情况是把VSG输出电压幅值、相位、频率都调到和电网一致然后再合闸。但即使相位幅值都对齐无功功率参考值不对还是会冲击。我调试的时候反复遇到并网瞬间无功震荡好半天才稳定后来发现是给定的Q_ref不是系统稳态无功把Q_ref改成负载消耗的无功值后冲击基本消失。2.4 滤波器和线路阻抗对VSG性能的影响VSG作为电压源型控制输出阻抗特性直接影响功率传输能力和稳定性。逆变器输出端一般加L或者LCL滤波器。LCL在滤除高频开关纹波方面效果好但在VSG控制下容易和虚拟同步电抗产生谐振交互。我踩过这个坑后来把LCL改成了单纯的L滤波器配合虚拟阻抗问题就清爽很多。虚拟阻抗的引入是VSG控制的一个常用技巧在控制回路里人为叠加一个阻抗压降E_ref E - Z_v·I_o这样做的效果是改变逆变器等效输出阻抗的比例削弱电网阻抗变化对控制性能的影响同时还能抑制环流。我设置的虚拟感抗为0.2 Ω标幺值。从波形上看加入虚拟阻抗后并网电流的谐波含量和振荡幅度都有明显改善。3. Simulink仿真模型搭建从整体架构到各模块实现细节3.1 系统顶层架构设计整个仿真系统在Simulink里按模块化搭我分成五个子系统风电机组模型、储能系统模型、VSG控制器、主电路和并网断路器、电网模型。这样做的好处是每个子系统可以单独调试长时间仿真时也能快速定位问题。主电路的直流侧我用了理想直流电压源和受控电流源的简化方式来表示储能加风机的合成外特性后期再替换成完整的Buck-Boost加电池模型。并网侧电网用三相电压源串联阻抗模拟。断路器用三相Breaker模块实现预同步完成后闭合。3.2 风电模型风速、气动、机械传动、发电机整流风机部分我用了经典的两质量块模型风速模型可以考虑基本风、阵风、斜坡风和噪声风组合。气动模型按标准气动公式P_w 0.5·ρ·π·R²·C_p(λ, β)·v³其中叶尖速比λ ω_r·R/v。C_p曲线我直接用了MATLAB里常用的多项式拟合系数Bezier那种查表也可以。发电机部分为了简化把永磁同步机整流环节等效为受控电流源接直流母线重点关注功率响应而非电机内部电磁细节。这里给一个我做仿真时经常提的注意点风机的MPPT控制会让风机输出功率随风速变化而波动这种波动对VSG来说是一种干扰输入。如果不希望论文里解释频率动态时被风机波动干扰可以在特定工况下把风速设成恒定值只保留阵风扰动去测试系统鲁棒性。项目里我做了两种工况恒风速下测试VSG惯量响应变风速下测试整体协调控制。3.3 储能双向变换器模型储能侧我搭了一个双向Buck-Boost的平均模型控制目标是把直流母线电压稳定在额定值。控制器是个双环外环电压环内环电流环。参数整定上带宽分离原则很重要电压环带宽设为电流环的1/5到1/10避免环路交互振荡。实际上这里有一个逻辑闭环储能稳压VSG控制功率两者分工明确。VSG输出功率变化时直流母线电压会产生波动储能变换器立刻响应去补功率。这个响应速度决定了VSG虚拟惯量能否真正实现。如果储能变换器带宽不够VSG的虚拟惯量响应就会被直流母线电压变化拖垮。我模拟的时候储能DC-DC的开关频率设在10 kHz电流环PI带宽2 kHz电压环带宽400 Hz整体响应速度完全够用。3.4 VSG核心控制回路实现这是整个项目最核心的部分。VSG控制器在Simulink里我分四步实现三相电压电流采样的Park变换旋转坐标系角度来自VSG自身生成的相位瞬时功率计算得到Pe和Qe转子运动方程求解生成虚拟角频率和相位角无功-电压下垂生成电压幅值参考再送入电压电流双环。用代码生成的函数形式容易入门后面我为了加速仿真部分环节改用S-Function和Stateflow配合。实际上用纯Simulink框图搭也完全没问题关键是积分器要设置好初值。相位生成环节特别要注意初始化的顺序。预同步前VSG的相位应该由锁相环追踪电网相位获得而不是从零开始积分。如果直接把初始相位设成0然后等它追赶电网相位动态过程又长又难看。我的做法是预同步阶段把VSG角频率参考绑定到PLL输出频率初始角度设置成电网相位初值合闸前再切到自洽的转子方程模式。3.5 预同步与并网控制策略并网瞬间的电流冲击是衡量控制好坏的重要指标。VSG的预同步逻辑主要有三个维度幅值对齐、频率对齐、相位对齐。幅值对齐通过无功-电压下垂的E0初值设定频率对齐就是把虚拟角频率的初值设成电网频率相位对齐则用一个微调积分器持续产生一个小的频率补偿量把相位差逼到零。我用的方法是在预同步PLL旁边加一个PI调节器输入是VSG输出电压和电网电压的相位差输出是一个频率微调叠加量作用到虚拟转子方程上。当相位差小于0.01 rad时给并网断路器发合闸信号。实测下来并网瞬间电流冲击被控制在额定电流的5%以内这个结果对常规工程来说已经很理想。3.6 并网之后控制模式的平滑切换合闸成功后预同步的PI调节器要退出否则它会持续干扰调频过程。这里存在一个模式切换的问题处理不好会引入切换扰动。最简单可靠的做法是用一个条件切换或多路复用器合闸前用预同步调节器输出合闸后用标准转子方程。切换时把预同步PI的输出清零或者保持最后一个输出值避免跳变。我实际遇到过一次切换后频率振荡的情况原因是合闸瞬间两个频率源电网和虚拟转子之间存在细微偏差切换时相当于给转子方程加了一个阶跃扰动。解决办法是把切换做成斜坡过渡让预同步调节器在0.2秒内逐渐退出而不是瞬间退出。这个细节在论文里讲出来是很加分的点。4. 关键参数整定与典型工况仿真分析4.1 VSG参数初值计算实例以一台额定功率2 MW、额定电压690 V的VSG逆变器为例。电网额定频率50 Hz直流母线电压1100 V开关频率5 kHz。设虚拟惯量惯性时间常数H 4 s则J 2H·S_n/ω_n² 2×4×2e6/(2π×50)² ≈ 8.1 kg·m²。阻尼系数D按期望阻尼比0.7和功角特性斜率计算。假设同步电抗等效为0.3 p.u.额定功角30度功角灵敏度K E·U·cosδ/X ≈ 0.126 p.u.则D 2×0.7×(J·ω_n·K)^0.5 ≈ 2×0.7×(8.1×314×0.126×2e6/314)^0.5。简化算下来D大约是几千牛米秒每弧度。实际MATLAB脚本里这些参数我都是写成可调变量便于扫描分析。下垂系数Kp按2%调差率设定对应的功率频率下垂为0.02×314 rad/s除以额定功率2 MW大概是3.14e-6 rad/(s·W)。这个值不算大但在仿真里如果不加单位换算很容易出错我这里特别提醒参数单位一定要统一到SI单位制。4.2 惯量和阻尼对动态性能的影响我做了一个参数扫描J从2 kg·m²扫到20 kg·m²记录负载阶跃下频率恢复时间和最大频率偏差。结果很符合经典的二阶系统规律J越大频率变化率越小最大频率偏差越小但恢复到稳态的振荡时间越长阻尼比太小会有明显的超调和振荡环阻尼比太大响应就变得迟缓。这里面有一个设计取舍和同步发电机调速系统里的思维是一样的。我做仿真的时候设定的目标是在10%负载阶跃下频率最大偏差不超过0.1 Hz恢复时间小于2秒。最终调出来的参数是H 5 s阻尼比约0.9。如果你期望的频率指标更严格虚拟惯量就还要加大但相应地储能系统的功率容量也要跟上否则直流母线电压会出现明显波动。4.3 典型场景仿真结果风速突变和负载阶跃风机风速从8 m/s阶跃到10 m/s的情况风机输出功率上升直流母线首先出现压升储能变换器吸收多余功率VSG输出功率缓慢上升频率有小幅震荡但不失控。这个场景验证了风储协调的核心逻辑储能做缓冲VSG做支撑。负载阶跃场景则更直接并网点新增20%的负荷系统频率跌落VSG依靠虚拟惯量快速提供额外功率频率变化率明显受限随后通过下垂环节建立新的稳态频率点。同时无功功率回路响应电压跌落维持并网点电压在允许范围。对比PQ控制的仿真结果后发现同样负荷阶跃下PQ控制频率跌落深度大、恢复慢而VSG把最大频率偏差缩小了40%以上这个对比能直接体现方案的工程价值。5. 常见问题与排查技巧实录5.1 Simulink模型频繁发散或数值爆炸Simulink报错一般集中在两类一是微分方程数值发散二是开关动作引发的代数环问题。VSG控制的强非线性在并网瞬间最容易触发数值问题。处理手段我的经验排序是先检查积分器步长仿真步长固定设置为1e-5秒甚至更小不要用可变步长因为高频开关和断续控制会导致步长剧烈波动其次检查是否形成了代数环可以尝试在反馈通道上增加单位延迟或者用Memory模块打断最后再排查参数极端值比如阻尼系数为0或者J过小导致的数值奇异。5.2 Bus Selector选不到信号的问题热搜词里有人搜“Simulink Bus Selector没有可选信号”这个我在项目里也遇到过。原因基本是仿真模型没更新或者信号类型没被Simulink识别。解决办法是先CtrlD更新模型再选如果是从Matlab Function输出的Bus数据需要确认总线的数据类型定义、输出端口的信号名是否已经正确设置。更推荐的做法是尽量用Goto/From标签或者StartSelector替代老式Bus Selector在大型模型里维护起来省心得多。5.3 并网瞬间电流冲击过大如果预同步相位偏差做不好合闸时冲击电流轻松超过2倍额定电流。排查顺序是先检查相位差是否真的在合闸时刻归零用Scope直接观察VSG电压和电网电压的波形重合情况再检查电压幅值是否一致虚拟励磁电压E0设置是否匹配最后检查无功参考值是否合理。经验表明相位差控制在0.01 rad以内才能得到理想的合闸效果。5.4 求解器稳定性和仿真速度的矛盾VSG控制系统里既有快速开关动态kHz级又有慢速机械动态Hz级这种多时间尺度系统用单一固定步长仿真会很慢。三分钟的动态可能要多仿真好几个小时。我在模型调试后期用了一个技巧把电力电子开关部分改成平均值模型电网侧用Thevenin等效这样求解器可以用更大的步长仿真速度能提升好几倍。平均值模型在商业模式验证层面精度完全够用只有需要看开关纹波的时候才换回详细模型。5.5 常被忽略的参数单位和标幺化问题VSG控制里国际单位制和标幺值混用是新手最容易出问题的地方。如果在一些模块用了标幺值另一些用了实际值转子运动方程输出的单位不一致整个控制就会偏移甚至失真。我的建议是整个控制系统统一使用标幺值功率基值取VSG额定功率电压基值取额定电压频率基值取额定角频率这样转子方程的参数和计算会干净很多下垂系数也更符合工程习惯。6. 一点项目复盘与扩展建议整个项目做下来我最深的感受是VSG不是一个现成落地的控制模块而是一套需要结合系统能量关系才能正确发挥作用的控制哲学。风储缺一不可惯性参数不是越大越好储能也不是一个可有可无的缓冲池。仿真模型里每一处参数的取舍背后都要回答“物理上能不能实现、经济上要不要这代价、控制上会不会失控”这三个问题。如果后续想在这个基础上扩展可以考虑的方向有VSG与电网电压暂降时LVRT控制的融合、多台VSG并联时虚拟阻抗对功率均分的影响、以及把Simulink模型通过C代码生成部署到DSP控制器做硬件在环测试。我目前正在做多机并联方向等有阶段成果了再来更新。
风储VSG并网Simulink仿真：从控制原理到工程实现
做风电并网仿真的人大概率都遇到过这么一种尴尬风机模型跑得挺欢逆变器控制也调得挺顺但一到并网那一下系统就跟你翻脸——电压闪变、功率抖动、频率掉坑甚至直接数值发散Simulink整块报红。原因很简单传统跟网型PQ控制本质上就是一个受控电流源电网给什么频率电压它就跟着跑电网一旦偏软它就撑不住。我最近完整做了一个风储VSG项目也就是基于虚拟同步发电机VSG的风储并网系统Simulink仿真从控制原理、参数整定到模型搭建一条龙走下来踩了不少坑也把很多关键细节真正想明白了。这篇就来把整个项目和盘托出给正在做新能源并网、微电网仿真或者VSG相关毕设和课题的朋友一份可以直接抄作业的实操参考。1. 方案选型为什么非得给风储系统配一个“虚拟发电机”1.1 风电并网的核心痛点双馈风机或者直驱永磁风机本身都是电力电子变换器接入电网转子运动系统和电网之间隔着直流母线机械惯量几乎没有传递到电网上。传统控制用锁相环盯住电网相位然后输出指定的有功无功电网频率波动时它基本不参与支撑。比例高起来之后整个电网系统惯量下降得厉害频率变化率RoCoF变大一次调频能力变弱这是老生常谈的问题但却是VSG方案存在的根本理由。我刚开始做这个课题的时候想的是能不能直接在风电机组侧改控制把变流器的PQ控制换成VSG控制。后来发现风机的功率本身受风速影响波动很大如果强行让风机承担虚拟惯量和调频任务功率指令波动和虚拟惯量所需功率叠加在一起直流母线电压根本稳不住。这是一个很容易被新手忽略的核心点VSG要求功率源能够在短时间内提供或吸收额外功率而风机的最大功率跟踪MPPT控制响应速度远不够快。1.2 VSG和传统并网逆变器控制的本质区别先看一个最直观的对比表对比维度传统PQ控制VSG控制外特性受控电流源受控电压源带内阻频率支撑不参与锁相环跟随具备虚拟惯性和一次调频电压支撑无功指令控制无功-电压下垂自动调压并网切换需要预同步处理同样需要但支撑能力更强参数整定PI参数为主惯量J、阻尼D、下垂系数等VSG的本质就是把同步发电机的转子运动方程搬进逆变器控制里。控制回路里人为构造一个“虚拟转子”用有功功率偏差去积分出虚拟角频率再用频率偏差积分出虚拟功角整个电压源的相位由此生成。这样电网频率跌落的时候VSG会因为它“有惯量”而自然地减速并释放存储的能量表现出来就是输出有功功率增加对电网形成支撑。这里有个很重要的地方要说透VSG的“惯量”不是凭空产生的物理飞轮它本质上是控制算法制造的对外动态特性功率的来源最终还是直流侧。所以风储储系统里“储能”不是辅助项而是必不可少的能量缓冲池。我做的项目里采用锂电池储能、双向DC-DC加逆变器的架构风机功率经过整流后和储能双向变换器一起挂在直流母线上再统一经过VSG控制的逆变器并网。这样VSG对外表现出来的惯量、阻尼、调频能力都由储能单元来兜底。1.3 储能系统在VSG中的能量平衡逻辑这就解释了为什么做风储VSG而不是单纯做“风机VSG”。虚拟同步机的惯量响应功率计算公式可以粗略写成ΔP_J ≈ J·ω_n·(dω/dt)频率变化率越大需要注入或吸收的功率就越大。一台2MW级别的风机如果设置惯性时间常数H为4秒在0.5 Hz/s的频率变化率下瞬时功率差额会非常可观风机MPPT不可能提供这种瞬态功率。但储能DC-DC的响应速度在毫秒级完全可以通过双向Buck-Boost快速拉高或降低直流母线功率。我实际搭仿真的时候储能模型用的是简化一阶惯性环节加荷电状态SOC计算没有用复杂的电化学模型。原因是在VSG并网系统级仿真里重点考察的是功率响应和电能质量电池内部细节对结果影响不大反而会增加仿真时间和数值刚性问题。当然如果课题侧重储能优化那就需要换更细的电池模型这是另一个方向。仿真初期用简化模型把控制链路跑通后期再逐步细化是我一直推荐的做法。2. VSG数学模型与参数整定先把公式吃透再动手搭模型2.1 转子运动方程和虚拟惯量参数设计VSG控制的核心公式转子运动方程J·(dω/dt) T_m - T_e - D·(ω - ω_n)功率形式下可以写为J·ω_n·(dω/dt) P_m - P_e - D·ω_n·(ω - ω_n)其中P_m是机械功率输入对应值这里是有功功率参考值PrefP_e是实际电磁功率D是阻尼系数ω_n是额定角频率。控制框图里有功功率参考Pref与实际输出功率Pe的差值经过惯性环节生成角频率ω再积分得到相位角θ。J的物理意义我习惯于用飞轮来理解转动惯量大相当于飞轮重转速变化慢对电网频率波动的抵抗能力强但动态响应慢功率超调大J小响应快但支撑能力弱。合适的J取值原则是让VSG的振荡频率落在合理范围内。VSG的虚拟功角特性近似二阶振荡环节自然振荡角频率可以估算为ω_n_vsg ≈ √(K·P_max / (J·ω_n))其中K是功角灵敏度系数在额定工作点附近近似为电网电压和逆变器输出电压乘积与同步电抗的比值。工程上更直接的做法是根据期望的振荡周期T和阻尼比ξ反推参数。我项目中设J对应的惯性时间常数H为3~6秒切换测试阻尼比取0.7左右对应的D为D 2·ξ·√(J·ω_n·K)这个推导过程说来简单实际调参的时候如果不理解纯靠盲试参数可能一天都试不出稳定波形。我在项目里用MATLAB脚本算了一组初值再做小范围扫描调整效率高很多。2.2 有功-频率下垂与一次调频VSG除了虚拟惯量还必须有稳态下的下垂特性也就是有功-频率下垂控制模拟发电机组的调速器。稳态时频率偏差和功率偏差满足ω ω_n - K·(P_e - P_ref)这里K就是下垂系数单位rad/(s·W)工程上常用%表示比如2%下垂意味着满载功率变化对应2%的频率偏移。并网接口处Simulink里通常把这个关系叠加上去替代转子方程中的阻尼项。实际上阻尼D和下垂K在稳态特性上有重叠初学时容易混淆。我自己的理解阻尼项影响的是动态过程的振荡衰减速度下垂项决定最终的频率-功率静态关系。两者可以统一成一个等效阻尼加死区环节但为了物理意义清晰我在模型里分开实现。仿真结果显示负荷突增2%时VSG频率短暂下降然后恢复到一个比额定值低约0.03 Hz的稳态值这与下垂设计值基本吻合。2.3 无功-电压下垂与励磁调节无功控制回路相对简单模拟同步发电机的励磁调节器表达式E E_0 K_q·(Q_ref - Q_e)E是虚拟励磁电动势Q是无功功率。这里K_q是电压下垂系数取值通常让无功功率从0变化到额定值时电压偏差在2%以内。逆变器输出电压的幅值就由E决定再通过电压电流双环控制最终生成调制波。无功回路有一个非常容易被忽视的坑初始无功参考值必须和实际负载无功匹配否则并网瞬间会有巨大的无功冲击。做预同步的时候理想情况是把VSG输出电压幅值、相位、频率都调到和电网一致然后再合闸。但即使相位幅值都对齐无功功率参考值不对还是会冲击。我调试的时候反复遇到并网瞬间无功震荡好半天才稳定后来发现是给定的Q_ref不是系统稳态无功把Q_ref改成负载消耗的无功值后冲击基本消失。2.4 滤波器和线路阻抗对VSG性能的影响VSG作为电压源型控制输出阻抗特性直接影响功率传输能力和稳定性。逆变器输出端一般加L或者LCL滤波器。LCL在滤除高频开关纹波方面效果好但在VSG控制下容易和虚拟同步电抗产生谐振交互。我踩过这个坑后来把LCL改成了单纯的L滤波器配合虚拟阻抗问题就清爽很多。虚拟阻抗的引入是VSG控制的一个常用技巧在控制回路里人为叠加一个阻抗压降E_ref E - Z_v·I_o这样做的效果是改变逆变器等效输出阻抗的比例削弱电网阻抗变化对控制性能的影响同时还能抑制环流。我设置的虚拟感抗为0.2 Ω标幺值。从波形上看加入虚拟阻抗后并网电流的谐波含量和振荡幅度都有明显改善。3. Simulink仿真模型搭建从整体架构到各模块实现细节3.1 系统顶层架构设计整个仿真系统在Simulink里按模块化搭我分成五个子系统风电机组模型、储能系统模型、VSG控制器、主电路和并网断路器、电网模型。这样做的好处是每个子系统可以单独调试长时间仿真时也能快速定位问题。主电路的直流侧我用了理想直流电压源和受控电流源的简化方式来表示储能加风机的合成外特性后期再替换成完整的Buck-Boost加电池模型。并网侧电网用三相电压源串联阻抗模拟。断路器用三相Breaker模块实现预同步完成后闭合。3.2 风电模型风速、气动、机械传动、发电机整流风机部分我用了经典的两质量块模型风速模型可以考虑基本风、阵风、斜坡风和噪声风组合。气动模型按标准气动公式P_w 0.5·ρ·π·R²·C_p(λ, β)·v³其中叶尖速比λ ω_r·R/v。C_p曲线我直接用了MATLAB里常用的多项式拟合系数Bezier那种查表也可以。发电机部分为了简化把永磁同步机整流环节等效为受控电流源接直流母线重点关注功率响应而非电机内部电磁细节。这里给一个我做仿真时经常提的注意点风机的MPPT控制会让风机输出功率随风速变化而波动这种波动对VSG来说是一种干扰输入。如果不希望论文里解释频率动态时被风机波动干扰可以在特定工况下把风速设成恒定值只保留阵风扰动去测试系统鲁棒性。项目里我做了两种工况恒风速下测试VSG惯量响应变风速下测试整体协调控制。3.3 储能双向变换器模型储能侧我搭了一个双向Buck-Boost的平均模型控制目标是把直流母线电压稳定在额定值。控制器是个双环外环电压环内环电流环。参数整定上带宽分离原则很重要电压环带宽设为电流环的1/5到1/10避免环路交互振荡。实际上这里有一个逻辑闭环储能稳压VSG控制功率两者分工明确。VSG输出功率变化时直流母线电压会产生波动储能变换器立刻响应去补功率。这个响应速度决定了VSG虚拟惯量能否真正实现。如果储能变换器带宽不够VSG的虚拟惯量响应就会被直流母线电压变化拖垮。我模拟的时候储能DC-DC的开关频率设在10 kHz电流环PI带宽2 kHz电压环带宽400 Hz整体响应速度完全够用。3.4 VSG核心控制回路实现这是整个项目最核心的部分。VSG控制器在Simulink里我分四步实现三相电压电流采样的Park变换旋转坐标系角度来自VSG自身生成的相位瞬时功率计算得到Pe和Qe转子运动方程求解生成虚拟角频率和相位角无功-电压下垂生成电压幅值参考再送入电压电流双环。用代码生成的函数形式容易入门后面我为了加速仿真部分环节改用S-Function和Stateflow配合。实际上用纯Simulink框图搭也完全没问题关键是积分器要设置好初值。相位生成环节特别要注意初始化的顺序。预同步前VSG的相位应该由锁相环追踪电网相位获得而不是从零开始积分。如果直接把初始相位设成0然后等它追赶电网相位动态过程又长又难看。我的做法是预同步阶段把VSG角频率参考绑定到PLL输出频率初始角度设置成电网相位初值合闸前再切到自洽的转子方程模式。3.5 预同步与并网控制策略并网瞬间的电流冲击是衡量控制好坏的重要指标。VSG的预同步逻辑主要有三个维度幅值对齐、频率对齐、相位对齐。幅值对齐通过无功-电压下垂的E0初值设定频率对齐就是把虚拟角频率的初值设成电网频率相位对齐则用一个微调积分器持续产生一个小的频率补偿量把相位差逼到零。我用的方法是在预同步PLL旁边加一个PI调节器输入是VSG输出电压和电网电压的相位差输出是一个频率微调叠加量作用到虚拟转子方程上。当相位差小于0.01 rad时给并网断路器发合闸信号。实测下来并网瞬间电流冲击被控制在额定电流的5%以内这个结果对常规工程来说已经很理想。3.6 并网之后控制模式的平滑切换合闸成功后预同步的PI调节器要退出否则它会持续干扰调频过程。这里存在一个模式切换的问题处理不好会引入切换扰动。最简单可靠的做法是用一个条件切换或多路复用器合闸前用预同步调节器输出合闸后用标准转子方程。切换时把预同步PI的输出清零或者保持最后一个输出值避免跳变。我实际遇到过一次切换后频率振荡的情况原因是合闸瞬间两个频率源电网和虚拟转子之间存在细微偏差切换时相当于给转子方程加了一个阶跃扰动。解决办法是把切换做成斜坡过渡让预同步调节器在0.2秒内逐渐退出而不是瞬间退出。这个细节在论文里讲出来是很加分的点。4. 关键参数整定与典型工况仿真分析4.1 VSG参数初值计算实例以一台额定功率2 MW、额定电压690 V的VSG逆变器为例。电网额定频率50 Hz直流母线电压1100 V开关频率5 kHz。设虚拟惯量惯性时间常数H 4 s则J 2H·S_n/ω_n² 2×4×2e6/(2π×50)² ≈ 8.1 kg·m²。阻尼系数D按期望阻尼比0.7和功角特性斜率计算。假设同步电抗等效为0.3 p.u.额定功角30度功角灵敏度K E·U·cosδ/X ≈ 0.126 p.u.则D 2×0.7×(J·ω_n·K)^0.5 ≈ 2×0.7×(8.1×314×0.126×2e6/314)^0.5。简化算下来D大约是几千牛米秒每弧度。实际MATLAB脚本里这些参数我都是写成可调变量便于扫描分析。下垂系数Kp按2%调差率设定对应的功率频率下垂为0.02×314 rad/s除以额定功率2 MW大概是3.14e-6 rad/(s·W)。这个值不算大但在仿真里如果不加单位换算很容易出错我这里特别提醒参数单位一定要统一到SI单位制。4.2 惯量和阻尼对动态性能的影响我做了一个参数扫描J从2 kg·m²扫到20 kg·m²记录负载阶跃下频率恢复时间和最大频率偏差。结果很符合经典的二阶系统规律J越大频率变化率越小最大频率偏差越小但恢复到稳态的振荡时间越长阻尼比太小会有明显的超调和振荡环阻尼比太大响应就变得迟缓。这里面有一个设计取舍和同步发电机调速系统里的思维是一样的。我做仿真的时候设定的目标是在10%负载阶跃下频率最大偏差不超过0.1 Hz恢复时间小于2秒。最终调出来的参数是H 5 s阻尼比约0.9。如果你期望的频率指标更严格虚拟惯量就还要加大但相应地储能系统的功率容量也要跟上否则直流母线电压会出现明显波动。4.3 典型场景仿真结果风速突变和负载阶跃风机风速从8 m/s阶跃到10 m/s的情况风机输出功率上升直流母线首先出现压升储能变换器吸收多余功率VSG输出功率缓慢上升频率有小幅震荡但不失控。这个场景验证了风储协调的核心逻辑储能做缓冲VSG做支撑。负载阶跃场景则更直接并网点新增20%的负荷系统频率跌落VSG依靠虚拟惯量快速提供额外功率频率变化率明显受限随后通过下垂环节建立新的稳态频率点。同时无功功率回路响应电压跌落维持并网点电压在允许范围。对比PQ控制的仿真结果后发现同样负荷阶跃下PQ控制频率跌落深度大、恢复慢而VSG把最大频率偏差缩小了40%以上这个对比能直接体现方案的工程价值。5. 常见问题与排查技巧实录5.1 Simulink模型频繁发散或数值爆炸Simulink报错一般集中在两类一是微分方程数值发散二是开关动作引发的代数环问题。VSG控制的强非线性在并网瞬间最容易触发数值问题。处理手段我的经验排序是先检查积分器步长仿真步长固定设置为1e-5秒甚至更小不要用可变步长因为高频开关和断续控制会导致步长剧烈波动其次检查是否形成了代数环可以尝试在反馈通道上增加单位延迟或者用Memory模块打断最后再排查参数极端值比如阻尼系数为0或者J过小导致的数值奇异。5.2 Bus Selector选不到信号的问题热搜词里有人搜“Simulink Bus Selector没有可选信号”这个我在项目里也遇到过。原因基本是仿真模型没更新或者信号类型没被Simulink识别。解决办法是先CtrlD更新模型再选如果是从Matlab Function输出的Bus数据需要确认总线的数据类型定义、输出端口的信号名是否已经正确设置。更推荐的做法是尽量用Goto/From标签或者StartSelector替代老式Bus Selector在大型模型里维护起来省心得多。5.3 并网瞬间电流冲击过大如果预同步相位偏差做不好合闸时冲击电流轻松超过2倍额定电流。排查顺序是先检查相位差是否真的在合闸时刻归零用Scope直接观察VSG电压和电网电压的波形重合情况再检查电压幅值是否一致虚拟励磁电压E0设置是否匹配最后检查无功参考值是否合理。经验表明相位差控制在0.01 rad以内才能得到理想的合闸效果。5.4 求解器稳定性和仿真速度的矛盾VSG控制系统里既有快速开关动态kHz级又有慢速机械动态Hz级这种多时间尺度系统用单一固定步长仿真会很慢。三分钟的动态可能要多仿真好几个小时。我在模型调试后期用了一个技巧把电力电子开关部分改成平均值模型电网侧用Thevenin等效这样求解器可以用更大的步长仿真速度能提升好几倍。平均值模型在商业模式验证层面精度完全够用只有需要看开关纹波的时候才换回详细模型。5.5 常被忽略的参数单位和标幺化问题VSG控制里国际单位制和标幺值混用是新手最容易出问题的地方。如果在一些模块用了标幺值另一些用了实际值转子运动方程输出的单位不一致整个控制就会偏移甚至失真。我的建议是整个控制系统统一使用标幺值功率基值取VSG额定功率电压基值取额定电压频率基值取额定角频率这样转子方程的参数和计算会干净很多下垂系数也更符合工程习惯。6. 一点项目复盘与扩展建议整个项目做下来我最深的感受是VSG不是一个现成落地的控制模块而是一套需要结合系统能量关系才能正确发挥作用的控制哲学。风储缺一不可惯性参数不是越大越好储能也不是一个可有可无的缓冲池。仿真模型里每一处参数的取舍背后都要回答“物理上能不能实现、经济上要不要这代价、控制上会不会失控”这三个问题。如果后续想在这个基础上扩展可以考虑的方向有VSG与电网电压暂降时LVRT控制的融合、多台VSG并联时虚拟阻抗对功率均分的影响、以及把Simulink模型通过C代码生成部署到DSP控制器做硬件在环测试。我目前正在做多机并联方向等有阶段成果了再来更新。
