这两年做微电网仿真有一个问题越来越绕不开系统里电力电子设备占比高了以后频率抗扰动能力肉眼可见地变差。以前电网频率出问题靠同步电机转子的物理惯量就能扛住第一波冲击现在很多电源都经变流器并网顶不上这个位置负荷一动频率就掉得又快又深。于是构网型变流器被推到了台前——说白了就是让电力电子装置学着同步电机的外特性把惯量和阻尼“装”进控制算法里。这篇文章是我自己完整做了一遍“同步电机构网型变流器频率稳定性”研究的记录重点放在Simulink模型怎么搭、VSG参数怎么算、频率曲线怎么看以及踩过的那些仿真坑。无论你是刚接触构网型控制的初学者还是已经搭过模型但被振荡、发散问题折磨过的人这份内容应该都能帮你省下不少试错时间。1. 频率稳定性为什么成了硬骨头构网型变流器为什么被寄予厚望1.1 同步电机的转子是系统频率的“物理压舱石”要理解构网型变流器在干什么先得把传统同步电机在频率稳定里的角色说透。同步电机的转子是一个实实在在转动的质量块它携带的动能可以近似写成 E_k (1/2) Jω²。当系统里出现负荷突增发电机的电磁功率瞬间变高但原动机输入机械功率来不及跟上转子就会把一部分动能释放出来表现为转速下降。这一过程不需要任何控制器参与完全是物理特性所以叫“惯性响应”。转速下降的速率与转子转动惯量直接相关惯量越大频率变化率越小给一次调频留出的时间就越充裕。频率稳定分析里经常用到一个“惯性常数”H定义是 H E_k / S_n也就是转子在额定转速下的动能除以发电机额定容量单位是秒。一台典型的大型汽轮发电机H通常在4到8秒之间小型柴油发电机可能在1到2秒。这个数字的意义在于当系统遭遇功率扰动时频率变化率大致可以按 RoCoF ≈ ΔP·f₀ / (2H·S_total) 来估算。也就是说如果总惯量足够大同样的功率缺额造成的频率下跌速度就慢保护装置和调度员都更容易处理。等到惯性响应把频率的下跌速度拖慢之后接下来登场的是调速器一次调频。同步电机的调速器检测到转速偏差调整原动机进汽量或进油量逐步把功率补充上来最终让频率稳定在一个新的平衡点。这两层配合——惯性响应决定“跌多快”一次调频决定“跌多深”——是传统电力系统频率稳定的底层逻辑。1.2 构网型变流器如何把惯量“装”进控制算法变流器并网的新能源没有转子自然没有物理惯性响应。常见的跟网型GF变流器靠锁相环盯住电网电压相位然后按指令注入功率本质上是一个“受控电流源”。电网存在稳定的电压和频率它才能工作一旦系统频率快速变化它不但不帮忙反而可能因为锁相环动态跟不上而加剧问题。这就是高比例新能源系统频率稳定性变差的直接原因。构网型变流器GFM的思路换了方向不再去追电网的相位而是自己建立一个电压和频率参考本质上成为一个“受控电压源”。最主流的实现方式是虚拟同步发电机VSG控制也叫VSG算法。它在控制环路里写入了同步电机的转子运动方程J·(dω/dt) P_ref/ω_nom - P_e - D·(ω - ω_nom)其中J是自定义的虚拟转动惯量D是阻尼系数P_ref是有功功率参考P_e是实测输出功率。每个控制周期按这个方程更新一次虚拟角频率再对ω积分得到相位从而决定输出电压的相位和频率。这样一来变流器虽然没有真实的旋转部件但它的外特性——频率随功率变化的动态过程——和同步电机非常接近。有一点必须心里有数虚拟惯量是“用控制算法模拟出来的”但它背后需要真实的能量支撑。如果构网型变流器的直流侧是恒压源那频率响应可以做得和同步电机几乎一样如果直流侧接的是电池或光伏惯量响应本质上是在调用直流侧储能此时必须考虑SOC和功率限幅否则虚拟惯量只是“账面值”而不是实际能力。这一点在做仿真时很容易被忽略但在工程化落地时恰恰最关键。1.3 为什么我选择Matlab/Simulink来做这个课题电网级的频率稳定性分析可以用PSCAD、DIgSILENT PowerFactory也可以用Matlab/Simulink。我做这个课题选择Simulink主要因为三个原因。第一是电气和控制信号的耦合非常方便。这个课题里同步电机模型、构网型变流器控制、负荷阶跃和频率测量都在同一个环境里既有电磁暂态层面的电压电流波形也有控制层面的功率、频率信号Simulink的信号线模型天然适合这种“电气控制”混合仿真。第二是Matlab脚本可以做批量扫参。研究频率稳定性必然要看不同惯量、不同阻尼、不同下垂系数下的响应差异。在Simulink里把关键参数设成工作区变量再用Matlab脚本循环改参数、自动运行仿真、自动提取指标效率比手动一个个跑高太多了。第三是模型复杂度可控。用Simulink里的Simscape Electrical模块库搭同步电机和受控电压源模型既不牺牲物理准确性又比PSCAD更贴近控制算法的调试逻辑。后面我会具体说明频率稳定性分析其实不需要特别精细的开关级模型用平均化受控电压源模型就够了。2. 先从参数说起系统拓扑、电机参数与VSG控制关键值的计算2.1 仿真系统拓扑设计这次频率稳定性研究核心要回答的问题是同步电机和构网型变流器在同一个系统里运行时频率响应会怎么样所以拓扑必须同时包含这两种电源并且要让两者共同承担负荷变化。我采用的是一个小型孤岛微电网结构一套同步发电机、一台构网型变流器、一条连接线路、一组可阶跃投切的负荷。系统额定电压我取400V线电压有效值额定频率50Hz总装机容量10kVA。这个量级对应实验室微电网平台仿真速度快参数也容易理解。对于频率稳定研究量级本身不影响物理规律关键看标幺化的惯量、阻尼和下垂系数取值。构网型变流器的直流侧在仿真里用理想直流电压源代替容量也给到10kVA。这样做是为了聚焦控制策略对频率动态的影响暂时不考虑储能能量约束。如果后面想研究构网型变流器在储能SOC限制下还能不能提供惯量再把直流源换成电池模型即可。负荷用恒功率负荷模型初始有功负荷8kW在t5s时额外投入2kW使总有功负荷阶跃到10kW对应20%的功率突变。这个幅值足够让频率动态清晰可见又不会让系统失稳。2.2 同步电机模型参数同步电机在Simulink里可以用Simscape Electrical模块库中的“Synchronous Machine”模块我选的是SI Fundamental参数类型所有参数用国际单位。这样虽然填参数时麻烦一点但不会出现标幺值基准不一致的问题。表同步电机主要参数10kVA / 400V / 50Hz参数数值说明额定视在功率10 kVA系统基准容量额定线电压400 V50Hz交流侧定子电阻 Rs0.05 ΩSI单位直接填写定子漏感 Ll1.2 mH约0.1标幺d轴同步电感 Ld31.8 mH对应Xd≈2.5 puq轴同步电感 Lq25.5 mH对应Xq≈2.0 pud轴暂态电感 Ld5.1 mH对应Xd≈0.4 pud轴次暂态电感 Ld1.9 mH对应Xd≈0.15 puq轴次暂态电感 Lq1.9 mH对应Xq≈0.15 pu惯性常数 H4 s转子动能/额定容量极对数2四极机转子的转动惯量J由H反算J 2HS_n / ω_m,nom²。这里的ω_m,nom是机械角速度。四极机额定机械转速是1500r/min对应机械角速度约157.1 rad/s。于是J 2×4×10000 / (157.1²) ≈ 3.24 kg·m²。如果模块里要求的是电气侧折算惯量需要按极对数关系换算但SI Fundamental模块的转动惯量是按机械轴填写的直接用这个数就行。同步电机的原动机我配了一个简单的调速器下垂系数R取5%。也就是频率偏差1%时机械功率调整对应额定功率的20%这是常见的调速器设定。有了调速器才能看到最终频率稳定在哪否则在负荷阶跃后频率会一直降下去没有平衡点。2.3 VSG虚拟惯量与阻尼参数整定全过程VSG参数是这次仿真的灵魂。核心计算公式是J 2HS_n / ω_n²先看惯性常数H怎么选。为了和同步电机形成对照我让VSG的虚拟惯量H也取一个可以比较的量级。本例中H_base取4s和同步电机一致。于是J 2×4×10000 / (314.16²) 80000 / 98696 ≈ 0.811 kg·m²要注意的是这里ω是电角频率不是机械角频率。因为VSG控制算法直接在电气坐标系里写转子运动方程所以虚拟转动惯量J的单位是kg·m²折算到电气侧这和同步电机模块里填机械侧J的数值不是一回事。很多第一次做VSG的人会把这两个数搞混结果就是虚拟惯量参数对不上频率响应做出来跟预想完全不符。阻尼系数D的选择稍微麻烦一点。D的物理意义是频率偏差1rad/s需要额外提供多少功率来维持平衡单位可以理解为W·s/rad。一个自然的基准是额定功率除以额定角频率D_nom S_n / ω_n 10000 / 314.16 ≈ 31.8 W·s/rad实际仿真中D取D_nom的1到2倍比较合适。我最后整定D50也就是约1.6倍基准阻尼。如果D取得太小频率扰动后会出现明显的超调和振荡取得太大虽然频率平稳但变流器对功率分配的响应太“软”类似于同步电机阻尼绕组过强的效果反而会让功率共享不均匀。有功-频率下垂系数Kf在VSG里可以有两种实现方式。一种是把下垂直接融合在转子运动方程里即P_ref本身就是有功功率参考不需要额外的频率反馈另一种是显式加上P_ref P_0 Kf·(f_nom - f_meas)。我这套模型用第一种结构更简洁也是工业上常见的做法。无功-电压下垂则用标准Q-V下垂E E_nom Kq·(Q_ref - Q_e)。Kq取5%即无功功率从零变化到额定值时电压幅值调整5%。频率稳定性研究里电压外环的影响相对次要但也要保证电压在合理范围内否则电磁功率计算会偏离预期。表VSG控制参数最终取值参数符号数值虚拟转动惯量J0.811 kg·m²惯性常数H4 s阻尼系数D50 W·s/rad无功下垂系数Kq0.05 pu功率测量滤波时间常数Tf0.02 s额定电压E_nom230 V相电压峰值2.4 频率稳定性评价指标与阈值做了仿真总得有个客观标准来比较不同参数下的结果。频率稳定性分析通常看三个指标最大频率变化率RoCoF、频率最低点值nadir、稳态频率偏差。RoCoF是扰动后频率变化最快时刻的变化率它主要反映系统惯量水平。孤岛微电网里很多保护装置都有RoCoF闭锁或动作逻辑一般要求不超过0.5 Hz/s甚至更严格。按照前面提到的估算公式20%负荷阶跃、H4s时RoCoF大约在1.25Hz/s左右单靠同步电机是压不下来的。这正好能直观展示为什么需要构网型变流器补惯量。频率最低点nadir是扰动后频率降到的谷值它和惯量、阻尼、调速器动作速度都有关系。实际系统里低频减载装置动作频率通常设定在49Hz左右如果nadir低于这个值就会触发切负荷所以nadir是非常硬的安全约束。稳态频率偏差由系统总的有功-频率调节能力决定。同步电机的调速器和构网型变流器的下垂特性共同把频率拉回一个高于扰动前但低于额定值的水平偏差大小取决于总的等效下垂系数。这三个指标在Matlab脚本里可以自动提取找回频率波形后用gradient函数算一阶导数找最大值即RoCoF同段波形找最小值即nadir取末尾稳态数据的平均值算稳态偏差。后面批量扫参的脚本就是这么写的。3. Simulink搭建全流程与Matlab自动化脚本3.1 主电路与同步电机模块的搭建模型我命名为 gfm_freq.slx里面分几个大块同步电机电源区块、构网型变流器电源区块、负荷区块、测量与记录区块。先放同步电机。从Simscape Electrical Specialized Power Systems库拖出“Synchronous Machine SI Fundamental”双击把前面表格里的参数填进去。机械输入选“Shaft torque”在端口接入一个阶跃信号作为原动机输出配合一个调速器子系统。调速器内部就是简单的比例控制测到转速偏差乘以R系数叠加到额定机械功率上再加一个一阶低通模拟原动机时间常数柴油机一般取0.5s左右。电机出口接一个断路器模块和电流电压测量模块。励磁方面我为了简化直接用恒定励磁电压不配AVR。原因是这次研究聚焦有功-频率动态恒定励磁可以让分析更清晰真正要做工程再把AVR加进来看电压耦合的影响。构网型变流器那一路我用了受控电压源而不是完整的IGBT两电平变流器。这是频率稳定性研究里很实用的简化手段受控电压源的输出幅值和相位直接由VSG控制算法决定跳过了直流电容、PWM调制和电流内环这些高频动态仿真步长可以放大速度能快好几倍。如果后续要验证构网型变流器在并网瞬间的冲击电流或者低电压穿越特性再替换成开关级模型不迟。3.2 VSG控制环路的搭建细节VSG子系统是整个模型最核心的部分。输入是参考有功P_ref、实测功率P_e、参考无功Q_ref、实测功率Q_e输出是三相受控电压源的调制参考。每个控制环路单独说明。功率测量环节用“Three-Phase V-I Measurement”模块测构网型变流器出口电压和电流再用“Instantaneous Active Reactive Power”模块算瞬时功率。这个模块输出的P和Q带有测量噪声尤其是用受控电压源模型时电路动态本身很快功率信号容易毛糙所以要加一个一阶低通滤波器时间常数Tf0.02s。这一步不能省略否则VSG的转子运动方程会被功率纹波干扰容易产生数值振荡。有功控制环核心方程是 J·(dω/dt) (P_ref - P_e)/ω_nom - D·(ω - ω_nom)。具体搭建时先做减法得到一个功率偏差信号除ω_nom换算成虚拟转矩再减掉阻尼项D·Δω得到的净转矩除以J就是角加速度。对角加速度积分得到ω再积分一次得到相位θ。积分器初始值设为额定值即ω_0314.16rad/sθ_00这样变流器起始就是50Hz。无功控制环根据Q_e和Q_ref的偏差做下垂输出电压幅值E。公式是E E_nom Kq·(Q_ref - Q_e)。这里E是相电压峰值E_nom取230×√2≈325V。E和θ合成三相参考电压可以写一个Matlab Function或者直接用三角函数模块组合。这里有一个调试经验值得分享VSG的ω积分器和θ积分器初始值必须设置正确且匹配否则模型一启动就产生巨大的频率偏移仿真前期波形完全不可用。而且θ积分器的初始值要按“电压相位为0时刻开始”来设如果接的线路带有阻抗还需要考虑初始相角差否则启动瞬间会有功率冲击。3.3 频率测量与结果导出研究频率稳定性首先要有一套可靠、统一的频率测量。同步电机侧我用PLL模块测频率构网型变流器侧由于相位θ是VSG自己生成的可以直接用dθ/dt换算。两者的频率都要经过滤波后再记录否则PLL的动态会造成虚假的暂态波动。PLL模块在Simulink里可以直接用“PLL”库模块内部是Park变换加PI调节器。PI参数按带宽30Hz左右整定。这里有个矛盾带宽太低测出来的频率响应滞后严重RoCoF数值会被低估带宽太高电网电压畸变和谐波会导致频率抖动nadir附近可能出现假尖峰。30Hz是我在10kVA微电网模型里试出来比较合适的折中。数据记录用“To Workspace”模块变量名freq和time输出格式选Structure With Time。同步电机频率和VSG频率分别记录方便后续脚本处理。仿真时间设30秒负荷在5秒时阶跃前5秒是稳态段用来确认模型稳定中间5到20秒是完整的暂态段后面10秒足够让频率落到新稳态。3.4 用Matlab脚本批量扫参手动每次改参数再运行效率太低。我把所有关键参数写在Matlab脚本开头Simulink模型里的参数一律引用工作区变量。这样批量扫参时只要改工作区变量再循环调sim函数就行了。扫参脚本的核心结构如下% 系统参数初始化 Sn 10000; % 额定容量 VA Vn 400; % 线电压有效值 V fn 50; % 额定频率 Hz wn 2*pi*fn; % 额定角频率 rad/s P_step 2000; % 负荷阶跃 W % VSG控制参数 H 4; % 惯性常数 s D 50; % 阻尼系数 J 2*H*Sn/wn^2; % 虚拟转动惯量 kg*m^2 % 扫参测试 H_set [2, 4, 6, 8, 10]; D_set 50; results table(); for k 1:length(H_set) H H_set(k); J 2*H*Sn/wn^2; % 运行仿真 simOut sim(gfm_freq, StopTime, 30); % 提取频率波形 f simOut.freq.Data; t simOut.freq.Time; % 分析暂态段 idx (t 5.1) (t 20); % 跳过阶跃瞬间的毛刺 rocof max(abs(gradient(f(idx), t(idx)))); nadir min(f(idx)); f_ss mean(f(t 23 t 30)); results [results; table(H, rocof, nadir, f_ss)]; end disp(results);这里有一个很实际的经验数据里的频率波形在前几个周期内往往有数值毛刺直接求gradient会把RoCoF算得离谱大。所以我的脚本把阶跃后0.1秒内的数据剔掉再算RoCoF这样得到的是“真实物理响应的最大变化率”而不是数值积分的人造尖峰。如果你发现算出的RoCoF有几百先检查是不是漏了这个滤除步骤。4. 频率响应结果分析惯量、阻尼和协同运行的规律4.1 基准工况下的频率跌落过程先看基准工况H4sD50系统是同步电机加构网型变流器并列运行负荷在5秒阶跃增加20%。频率响应曲线的典型过程分三个阶段。第一阶段是阶跃后最初约0.5秒频率开始下降但下降速率受系统总惯量约束。同步电机贡献物理惯量构网型变流器贡献虚拟惯量两者叠加后RoCoF明显比单台同步电机时低。第二阶段是下降持续到最低点时间大约在阶跃后1.5到2秒。RoCoF指标更多由惯量决定而nadir还取决于阻尼和调速器开始提升功率的速度。第三阶段是频率从谷底回摆逐渐稳定到新的稳态值这个值由同步电机调速器下垂和VSG有功均衡共同决定。我在模型里同时记录了同步电机电磁功率和VSG输出功率可以看到负荷突增后瞬间同步电机的电磁功率立刻上升把转子动能释放出来对付需求缺口VSG因为虚拟惯量的存在输出功率也在极短时间内快速拉升。两个电源的出力变化几乎是并行的这正是构网型变流器被看作“跟同步机同台竞技”的原因——它能在第一时间参与功率平衡而不是像跟网型变流器那样傻等着锁相环更新。4.2 不同虚拟惯量H下的扫参结果把H从2s扫到10s其余参数不动表格结果如下。表虚拟惯量H对频率稳定指标的影响H (s)RoCoF (Hz/s)Nadir (Hz)稳态频率 (Hz)21.3848.7149.5240.9248.8949.5160.7148.9849.5180.5849.0449.51100.4949.0749.51从数据能读出两个关键规律。第一惯量H对稳态频率几乎没有影响稳态频率完全由下垂特性和调速器决定所以H不需要堆很大来改善稳态控制精度。第二H对RoCoF的影响非线性递减从2s到4s改善明显从8s到10s收益开始变弱这说明只靠增大虚拟惯量不是万能药而且过大的虚拟惯量意味着直流侧需要支撑更长时间的能量输出实际系统未必承受得住。还有一个细节H增大后nadir确实抬高了但回摆过程变慢系统恢复到稳态的时间拉长。这跟同步电机大惯量机组的物理特性完全一致——惯量是“缓冲”但不是“恢复力”恢复速度要靠阻尼和控制作用。如果H取得太大可能出现频率缓慢爬回稳态的“拖尾效应”对调度运行来说也不是好事。4.3 同步电机与构网型变流器协同运行的风险点两个电源都具备频率调节能力协同运行并不是简单地把两者参数加起来就完事。我在仿真里试过把同步电机的H设为4s、VSG的H也设为4s结果系统确实总惯量翻倍但把波形放大看同步电机转子和VSG虚拟转子之间出现了一个低频功率振荡模态频率在0.5到1Hz之间。这个模态的物理来源很容易理解两台“发电机”通过线路阻抗连接一个转子是真铁的一个转子是算法算出来的两者之间有功率交换形成的等效“弹性连接”就好像两个质量块中间夹了一根弹簧。如果线路阻抗大、连接弱这个振荡越明显如果VSG的阻尼D不够振荡衰减很慢极端情况下会与同步电机调速器产生不良交互。处理这个问题的套路是在VSG控制里加强阻尼项或者适当减小虚拟惯量让两者的动态时间尺度拉开。我的最终整定结果是同步电机H4s、VSG H3s、VSG D50这样既能提供可观的惯量支撑又不会和同步电机的转子振荡强耦合。实际的参数整定当然要根据具体容量和阻抗来微调但“两个惯量源不要靠得太近”这个原则是通用的。5. 仿真路上的坑常见问题与调试实录5.1 代数环和数值发散第一次搭VSG模型十有八九会遇到“代数环”警告。原因在于功率测量模块的输出直接参与了转子运动方程的计算而功率本身又由电压电流决定电压又由θ生成θ又依赖功率——这形成了一个无延迟的闭合循环Simulink必须迭代求解严重时直接发散。解决办法有两个按优先级来。最推荐的做法是在功率测量后加一阶低通滤波时间常数取0.01到0.02秒滤波本身让功率信号延迟了一点点代数环被物理上“打破”。另一个做法是在功率测量模块和VSG控制之间加一个Memory模块等于引入一个步长的计算延迟。Memory模块在仿真精度要求不高的场景能用但如果你的模型里有许多高频动态额外的延迟会改变控制相位反而制造振荡不如滤波干净。数值解算器也很关键。带PLL和降压器模型呈现一定的刚性我用ode23tb变步长。相对误差设1e-4最大步长限制在1e-3秒。不要为了追求速度把最大步长设太大否则频率波形上会出现阶梯状伪影直接影响RoCoF计算的准确性。5.2 PLL测量造成的人为误差构网型变流器的频率可以直接从θ微分获得但同步电机侧的频率必须经过PLL测量。如果PLL带宽和VSG的功率外环带宽比较接近PLL的滞后特性会让测量的频率信号“看起来”更平稳从而让你误判系统的真实稳定性。我在调试过程中做过一个对比实验把PLL带宽从30Hz降到5Hz测得的RoCoF数值下降了将近三成nadir数值也偏高。这完全是测量手段造成的假象不是真实物理响应。所以做频率稳定性分析时要么把PLL带宽固定且合理20到40Hz要么直接使用模型内部的转角信号换算频率。混用两种来源的数据时切记先确认测量延迟一致否则后续参数整定都会被带偏。5.3 功率外环振荡怎么压下去VSG有功环在特定参数组合下会出现等幅或缓慢衰减的振荡频率在2到5Hz。排查发现最常见的诱因是虚拟惯量J取得过小导致角频率变化率放得过大阻尼项D相对不足。按控制系统的视角VSG有功环本质是一个二阶振荡环节自然振荡频率大约正比于√(S_n/(J·ω_n))阻尼比正比于D/J。我遇到的几次振荡都是因为J减半后没有同步调整D阻尼比掉到了0.3以下结果频率曲线出现好几个周期的过冲回调。经验做法是把阻尼比目标设为0.5到0.7根据当前J反算D。还有一个简单实用的调试技巧给频率反馈信号加一个带宽稍低于振荡频率的低通滤波器可以从反馈通道上直接把振荡模态压掉大半但要注意不能影响正常的惯性响应带宽。5.4 单位与标幺值使用错误同步电机模块里填的是国际单位线路阻抗也是国际单位但VSG控制里如果用标幺值就必须小心基准。我这套模型选择全部统一的国际单位制J用kg·m²D用W·s/rad功率用W电压用V频率用rad/s。这样虽然数值有时很小有时很大但不会出基准错误。如果你习惯用标幺值那所有和“额定容量”相关的计算像J 2H/ω_n标幺惯性常数和D的标幺值必须放在同一个基准容量下。两个电源的容量不同标幺值参数必须折回公共基准否则功率共享比例完全乱套。我见过太多次因为“总觉得标幺值都是1上下”而把参数填错的情况所以这次直接放弃标幺全部有名值调试时用scope看波形也直观得多。5.5 结果处理建议频率波形从scope里看可能觉得没问题但提交论文或报告时通常还需要平滑处理和定量计算。我建议所有频率指标的计算脚本化不要手动在scope上读数。把To Workspace记录的数据导入Matlab脚本用gradient算RoCoF、用find峰谷算nadir、用均值算稳态偏差结果统一存表格全程可回放。这样调整参数后再跑一遍所有指标自动更新对比实验的效率完全不是手动操作能比的。一点实际操作后的体会这套“同步电机构网型变流器”频率稳定性仿真做下来我最大的体会是构网型技术不是简单地把同步电机的转子运动方程抄进控制环里就完了。物理世界里同步电机的惯量是一种被动存在的能量储备响应天然无延迟而虚拟惯量本质是控制系统对直流侧能量的一种主动调度它受测量延迟、滤波时间常数、PWM载波周期和储能容量的诸多约束。仿真里把H设成10s很简单但真要落地你得回答“能量从哪来”这个基本问题。所以做研究时我建议把模型分两层先用理想恒压源的构网型变流器把控制算法跑通理清楚惯量和阻尼的整定规律再逐步替入电池模型、限幅环节、电流内环检验虚拟惯量在真实约束下还能兑现多少。一步一步对照RoCoF、nadir、稳态频率三个指标你会对构网型变流器的能力边界有非常具体的认识。希望这份记录能帮你少踩几个坑。
