同步电机与构网型变流器频率稳定性：Simulink虚拟同步机仿真研究

1. 同步电机与构网型变流器频率稳定性研究的“同频共振”起点这两年做新能源并网仿真尤其是在Matlab/Simulink里做微电网或储能PCS控制一个绕不开的话题就是构网型变流器。它之所以火最根本的原因是传统火力发电和水力发电里的同步电机用一种近乎“物理必然”的方式维持着电网频率。但现在新能源占比上来了风机、光伏走的是电力电子变流器接口它们以前绝大多数是“跟网型”的靠锁相环去追外部电网的频率和相位本身不提供惯量和频率支撑。电网一旦出现大扰动频率跌落速度会比以前快得多这就是频率稳定性问题的来源。把同步电机和构网型变流器摆在一起研究频率稳定性本质上是在新能源变流器上“复刻”同步电机那种天然的惯量响应和阻尼特性所以也叫虚拟同步机技术。这类研究在Simulink里做非常适合的一类目标是电气工程专业研究生、系统工程师或者是在做储能变流器、新能源场站调频控制的一线开发人员。通过Simulink中的电气模型和控制算法你可以把同步电机的摇摆方程、下垂控制、虚拟惯量这些抽象概念变成一条条还能复用的模型和可测的波形曲线。这类课题看起来并不新但真正做到位的人不多。我见过不少初学者把同步电机和构网型变流器当成两个孤立的模型来做最后只是拿个示波器看频率没跌落就不管了这离真正的“频率稳定性研究”还很远。本文我想从一个相对完整的研究视角聊聊怎么在Simulink里把同步电机和构网型变流器放到同一个平台下做频率稳定性研究包括控制原理、模型搭建顺序、参数初始化脚本以及我在实测仿真中踩过的不少坑。2. 把“转子惯量”写进控制器构网型变流器的虚拟同步机映射与参数映射逻辑2.1 同步电机为什么会天然稳定同步电机的频率稳定能力核心来自机械旋转部分。转子有质量转动惯量用 ( J ) 表示当电磁功率和机械功率不平衡时转子的运动方程决定了转速变化不会突变而是按惯性缓慢变化。用摇摆方程表示[ J\frac{d\omega}{dt} T_m - T_e - D(\omega - \omega_0) ]其中 ( T_m ) 是机械转矩( T_e ) 是电磁转矩( D ) 是阻尼转矩系数( \omega ) 是电角速度。将这个方程乘以转速 ( \omega )就可以转化为功率形式的方程。工程上更习惯用惯性时间常数 ( H ) 来表示惯量的大小 ( H ) 表示在额定转矩作用下转子从静止加速到额定转速所需的时间。这带来了两个关键特性惯性响应频率变化率RoCoF被限制在一定范围内不会瞬间掉下去一次调频转子转速下降的同时调速器动作改变原动机输入功率使频率在新的平衡点恢复。构网型变流器要“仿”的主要就是这两点。2.2 构网型变流器的控制结构构网型变流器听起来很玄乎但控制结构上其实可以清晰拆成有功-频率下垂、无功-电压下垂、虚拟阻抗、电压/电流双闭环这几个环节。如果去掉锁相环通过功率计算和下垂方程生成相位和幅值就构成了典型的构网型控制如Droop控制、VSG控制。有功-频率下垂关系[ f f_0 - k_p(P - P_0) ]无功-电压下垂关系[ U U_0 - k_q(Q - Q_0) ]这与同步电机的一次调频特性异曲同工。如果要更贴近同步电机的动态行为就在有功-频率下垂环前面再加上惯性-阻尼环节这就成了虚拟同步机VSG控制。此时控制器里的 ( H ) 和 ( D )替代了同步电机转子里的 ( J ) 和阻尼绕组的作用。经过这样一番映射你在Simulink里搭建两个模型作对比研究时就比较清楚该调哪些参数为什么调。同步电机模型改惯性时间常数直接改机器参数构网型变流器模型则改控制器的 ( H ) 和 ( D )两者放在同一个频率扰动场景下对比波形上的惯性响应、频率跌峰、稳态偏差就变成了可量化、可对比的指标。2.3 “构网”和“跟网”的本质差异把构网型变流器与跟网型变流器放在一起对比也有利于理解频率稳定性的本质。跟网型控制必须有电网电压作为参考无法自主建立电压和频率在弱电网和孤岛工况下较难生存构网型控制本质上是构建一个可控电压源输出电压的幅值和频率由自身给定。构网型由于能够瞬时建立电压源特性在电网频率跌落时可以主动提供有功支撑这是它参与频率稳定研究最值得关注的原因。注意构网型并不等于一定带惯量。有些构网型控制只是单纯的静态下垂没有惯性环节还有一些构网型控制为了应对限流问题引入了虚拟阻抗和功率同步环。做研究时不要笼统地把“构网型”等同于“高惯量”这一点会在仿真的动态响应上体现得很明显。3. Simulink仿真环境与主电路建模从容量、电压到滤波参数的实现顺序3.1 仿真模型的架构选择在Simulink里做这类研究一个实用且不过度复杂的模型架构是储能或直流电压源模拟机侧。三相两电平电压源型逆变器。逆变器出口经L或LCL滤波器接入交流母线。交流母线接负载并通过接触器与三相电网电压源连接。这种“变流器负荷电网”的架构足以覆盖多数频率稳定性研究场景。比直接构建一个大电网系统来得清晰也可以很容易地改变负荷大小、电网频率扰动时刻、电网阻抗数值来生成不同评估工况。同步电机模型则单独并列一套也接在同一交流母线上。在Simulink库中同步电机用“Synchronous Machine”标准模块需要设置为发电机模式或电动机模式根据研究对象设定原动机输入。一般情况研究频率稳定性时可以设置同步电机作为主力电源构成一个离网微电网场景也可以把它当作调频机组并联在电网上与构网型变流器共同承担负荷波动。模型整体不需要过于庞杂重点是两侧都能调节关键参数并且能输出有功功率、无功功率、母线频率、电压幅值这四类信号用于分析频率稳定特性。3.2 仿真模型中的核心参数设置模型参数设置决定了仿真是否稳定可靠。下面是我常用的一组基础参数实际应用时可以在此基础上扩展额定容量( S_{rated} 10 ) kVA额定线电压( V_{LL} 380 ) V三相交流额定频率( f_0 50 ) Hz直流母线电压( V_{dc} 800 ) V滤波电感( L_f 3 ) mH滤波电阻( R_f 0.1 ) Ω滤波电容( C_f 20 ) μF载波频率( f_{sw} 10 ) kHz控制周期( T_s 100 ) μs注意Sentences里的直流母线电压要高于交流侧线电压幅值否则整流或逆变环节会产生耦合失真。通常 ( V_{dc} \sqrt{2} \cdot V_{LL} )所以这里取800V是够用的。滤波电感的取值会影响电流纹波和动态响应太小电流纹波大太大电压降明显实际调参时要兼顾。3.3 构建构网型变流器主电路主电路搭建顺序在Simulink中放置“Universal Bridge”选择器件类型为IGBT/Diode。为了准确模拟死区和谐波要合理设置开关管导通电阻、缓冲电路参数。逆变器输出端串联三相电感和电阻。如果研究弱电网特性可以在电网侧加入可调电感模拟线路阻抗。在公共连接点PCC接入电压电流测量模块。不少初学者为了省事会把功率模块和控制模块用近似的一阶惯性环节直接简化这在研究“控制策略”时可能勉强可行但一旦涉足“频率稳定性”主电路的动态特性特别是滤波电感和电容的交互会直接影响频率跌落和恢复过程所以从完整主电路起步是比较合理的。3.4 构建同步电机模型侧同步电机侧相对成熟直接使用Simulink基础库的同步电机模块设置好转子类型隐极或凸极、阻尼绕组、定子电阻、漏感等参数其原动机可采用调速器模型如燃气轮机或水轮机调速器。如果是简化的对比研究用给定额定机械功率输入也是常见的做法。有两个参数对频率稳定性影响较大惯性时间常数 ( H )同步电机模拟时模型库用转动惯量 ( J (kg·m^2) ) 形式输入需要根据极对数和额定转速折算。阻尼绕组等效 ( D )阻尼绕组的作用是抑制转子振荡即便同步电机自身物理阻尼很小在动态仿真中也可以通过附加阻尼绕组参数来模拟实际阻尼效果。4. 频率稳定性测试场景的搭建与Matlab初始化脚本从惯量常数到负荷突变事件4.1 Matlab初始化脚本的功能分区在仿真前最好将关键参数都写在Matlab初始化脚本中这样修改参数后更新基座即可不需要频繁在Simulink框图里翻找对参数扫描和敏感性分析也友好很多。下面给出一个可直接使用的初始化脚本框架%% 构网型变流器与同步电机频率稳定性研究 % 基本电气参数 clear; clc; close all; %% 主电路参数 f0 50; % 额定频率Hz V_LL 380; % 线电压有效值V V_dc 800; % 直流母线电压V P_nom 10e3; % 变流器额定有功功率W Q_nom 0; % 变流器额定无功功率Var % LCL滤波参数 L_f 3e-3; % 滤波电感H R_f 0.1; % 滤波电阻Ohm C_f 20e-6; % 滤波电容F % 电网阻抗用于弱电网模拟 L_g 1e-3; % 电网等效电感H R_g 0.15; % 电网等效电阻Ohm %% 构网型变流器控制参数 Dp 1e-3; % 有功-频率下垂系数 Dq 1e-3; % 无功-电压下垂系数 H_vsg 4; % 虚拟惯量时间常数s D_vsg 20; % 虚拟阻尼系数 % VSG控制方程用到的相关变量 Ki_p 0.1; % 有功功率环积分增益 Kp_i 8; % 电流内环比例增益 Ki_i 200; % 电流内环积分增益 Kp_u 0.5; % 电压外环比例增益 Ki_u 50; % 电压外环积分增益 %% 同步电机参数 S_sm 10e3; % 同步电机容量VA V_sm 380; % 同步电机额定电压V H_sm 3; % 同步电机惯性时间常数s D_sm 5; % 同步电机等效阻尼系数 %% 仿真控制 T_stop 10; % 仿真时长s T_load_change 5.0; % 负荷突变时刻s P_load_step 4e3; % 负荷突变功率W这段脚本能做的事情是统一管理所有关键参数。我建议你在脚本末尾再添加一个assert()用于检查直流母线电压是否满足逆变条件以及控制周期是否小于载波周期这能有效避免参数设置式错误导致的仿真异常。4.2 在Simulink中构造频率扰动事件频率稳定性研究最重要的支撑是扰动事件注入。常见事件类型有负荷阶跃在原有负载基础上突增有功负载比如在( t 5s )时投入4kW电炉负载观察频率跌落幅度和恢复时间。电网频率跌落模拟大电网频率波动如从50Hz跌到49.5Hz观察构网型变流器与同步电机的功率输出响应。孤岛切换断开并网开关使微网转入孤岛运行这是最严苛的频率稳定性测试工况。在Simulink中可通过Step模块设置时间变量来控制断路器闭合也可以直接通过可变电阻模块模拟负荷投入。这里我习惯用Breaker模块相比可控电压源切换更接近真实电气行为。4.3 测量指标与分析逻辑仿真不应只停留在波形好看。频率稳定性研究需要明确指标频率跌落最低点Nadir频率跌落深度单位Hz越低越危险。频率变化率的最大值Max RoCoF( \frac{df}{dt} )表示频率变化的剧烈程度。恢复稳态时间和振荡次数频率重新回到允许误差带的时间以及是否出现振荡。在Simulink中频率一般由母线三相电压经过锁相环测量获得或者使用PLL模块测出相位后求导。要注意简单的过零检测算频率容易在谐波较大时抖动一般用二阶广义积分器或者同步参考坐标系锁相环更稳定。5. 仿真中的隐性壁垒频率波动指标、求解器设置与参数整定的避坑手册5.1 求解器设置对频率稳定性仿真结果的影响仿真频率稳定性属于动态过程仿真步长配合敏感度高。有些初学者一上来就用ode45结果模型运行速度慢且波形噪声大。我的实践是主电路为连续系统控制为离散系统时用ode23tbstiff/TR-BDF2能较好兼顾步长收敛性。如果控制系统完全离散化且采用定步长控制此时主电路的开关频率和控制周期不同易产生数值发散。建议把控制系统的采样周期设为载波周期的整数倍比如载波从10kHz控制周期设为100μs即整数倍采样。另外要打开“Zero-Crossing Detection”或改为“Non-adaptive”零交叉检测防止断路器开断瞬间产生大量计算步长堆积直接仿真到时间耗尽都推不完一个暂态过程。这类问题在正弦电压源切换场景中特别常见。5.2 虚拟惯量 ( H ) 与阻尼 ( D ) 的参数整定陷阱虚拟惯量和阻尼越往大调整频率跌落深度会不会改善实测效果并非如此简单。在一定范围内增加惯量 ( H )可以将RoCoF压得比较低让频率变化变慢。但是 ( H ) 过大会导致功率环响应变慢频率跌落到最低点反而可能更深因为支撑功率释放慢。阻尼 ( D ) 过小时系统容易振荡有功功率响应会出现正弦衰减振荡( D ) 过大时频率恢复能力变差存在稳态频差。所以调参场景静态来看应该是一个“先调 ( H ) 限制RoCoF、再调 ( D ) 抑制振荡、最后配合下垂系数消除稳态偏差”的过程。我在实际仿真中一般先令 ( D 0 ) 观测功率是否产生等幅振荡再逐步增加 ( D ) 观察振荡衰减时间直到功率响应满足衰减比要求。5.3 锁相环与大扰动环境下的可靠性问题构网型变流器的频率观测有时候仍然依赖锁相环或频率测量环节。在电网电压畸变严重的场景下锁相环的动态带宽和滤波能力会影响频率反馈的准确性进而影响频率稳定性控制的性能。建议在Simulink中设置两套频率观测途径一套用标准同步参考坐标系锁相环用于对比一套用自适应带通滤波器加过零检测观测扰动工况下的差异。如果只想要构网型变流器在弱电网下挂网运行可以不用锁相环完全依靠功率下垂控制的方向但频率指标反馈链路建议保留便于后续做自动发电控制或二次调频。5.4 滤波电容与弱电网阻抗的谐振问题LCL滤波器和电网阻抗之间可能形成谐振在频率跌落瞬间容易出现高频振荡表现为电压电流波形出现明显纹波。对频率稳定性研究而言这种高频纹波会污染测量值导致锁相环输出频率抖动控制量抖动又反过来放大纹波。解决谐振的基本方案有两个方向增加无源阻尼在滤波电容支路串联电阻典型值20Ω至50Ω带来一定功耗但效果好直观可靠。采用有源阻尼在电流环中增加电容电流前馈或虚拟阻抗Simulink中实现也不复杂。从工程角度宁可为频率稳定性研究多一层有源阻尼因为纯无源阻尼会降低系统效率且寄生参数较多时计算容易偏离。6. 构网型变流器与同步电机联合仿真的工程启示与延伸空间6.1 从仿真结果看两者协同调频的分工当我按上述模型将同步电机和构网型变流器并联在同一母线后负载突变下的频率跌落实验表现出了一些有意义的分工特征。同步电机频率响应曲线表现为惯量响应和阻尼振荡动态响应相对慢但稳态带载能力强频率恢复依赖调速器。构网型变流器的虚拟同步机控制则表现为下垂外特性功率响应较快能在频率跌落初期迅速释放功率在频率恢复阶段可能会略微增加稳态功率输出。两者并联时整个系统频率跌落深度明显小于单一电源带同等负载改变。这种互补性在新能源接入比例高的新型电力系统中表现为储能变流器作为一种快速频率响应资源与常规火电机组的慢速调频形成时间尺度上的配合。6.2 拓展为参数扫描与敏感性分析频率稳定性研究中往往需要回答一些问题某一台机组的惯量低至多少会导致系统失稳负荷突变幅值与频率跌落深度是否有线性关系在Simulink中通过Monte Carlo仿真或参数扫描可以更系统地探索。使用Matlab脚本驱动Simulink仿真时可以通过batchsim或者循环修改模型参数、保存波形、提取指标的方式批量完成。比如在for循环里修改H_vsg从1到8然后运行模型获取MaxRoCoF和Nadir值并绘制关系曲线就能非常直观地看出惯量配置规律H_list 1:0.5:8; Nadir_list zeros(size(H_list)); for i 1:length(H_list) H_vsg H_list(i); % 更新工作区变量 sim(gfc_sm_freq_stability.slx); % 提取频率信号计算Nadir freq_signal freq_log.signals.values; Nadir_list(i) min(freq_signal); end figure; plot(H_list, Nadir_list, bo-, LineWidth, 1.5); xlabel(虚拟惯量 H (s)); ylabel(频率跌落最低点 (Hz)); title(虚拟惯量对频率稳定性的影响); grid on;这个脚本可以帮助快速判断构网型变流器的惯量参数对系统频率稳定性的影响规律很适合毕设和工程前期选型参考。6.3 当前构网型变流器频率稳定性研究的一些真实挑战虽然有了以上基础模型但这种仿真研究仍有两个较大的简化假设一级控制器设计完成后没有考虑二次调频AGC与三次调频经济调度的时间尺度交互虚拟同步机虽然有惯性模拟但物理上不存储动能其能量来自直流母线或储能电池因此构网型变流器提供惯量时背后必须有充足的储能支撑。在仿真中如果直流母线幅值波动过大通常意味着储能已经处于出力极限此时的频率支撑能力并不可持续。对严肃的研究项目建议把直流母线电容的尺寸和储能单元的功率约束引入模型而不是一直把直流电压当作恒定值。这样会更多涉及到双向DC/DC和电池模型模型复杂度上来了但仿真结论在工程上更有意义。6.4 我所经历过的一条现实路径在亲手做这个课题的过程中我越发意识到好用的仿真模型一定不是把参考模块一股脑拖进来完事关键是建模前想清楚下面的链路研究问题定义在哪个时间尺度频率稳定性主要指惯量响应、一次调频还是二次调频构网型变流器和同步电机各自扮演什么角色参数变化对结论敏感度如何如果这些在模型搭建前有了清晰答案后续在Simulink里的工作就会顺畅很多。同步电机侧明明白白告诉你什么叫做“物理惯性”构网型变流器则告诉你通过控制算法“虚拟惯性”也可以做到同样的频率支撑效果。真正难的部分是让这两种惯性来源在同一个系统里协调配合既不冲突又不失稳。从更广阔的视角来看这种研究也会帮助做储能系统设计的人更清楚变流器在系统频率风险面前的责任边界。毕竟在不久的将来面对极低系统惯量的运行场景构网型变流器很可能不再是一个可选项而是一个必选基础功能。另外如果后续做工程应用建议多在硬件在环测试平台上验证一下仿真模型的控制参数。毕竟Simulink仿真中对理想电压源、理想断路器等抽象环境的依赖在实际样机中都会变成死区时间、采样延时和通信延迟这些细节同样会影响频率稳定性的真实表现。