做了几年电力电子仿真要说储能单元里最让我印象深刻的**超导磁能储存系统SMES**绝对算一个。名字听着高冷实际在Simulink里搭起来反而不复杂一个线圈、一套功率变换器、一组控制算法剩下全是数学处理。和锂电池、飞轮储能相比SMES的响应速度可以达到毫秒级功率密度极高非常适合电网暂态支撑、电能质量治理这种“瞬间要出力”的场景。这也是为什么我建议做电力系统方向的同学或者搞新能源并网的朋友把SMES建模当成一次很好的练手项目。本文我会从原理、数学模型、控制策略到Simulink实现把整个流程完整走一遍。你不需要提前掌握超导物理细节只需要有基本的电路和电力电子知识配合MATLAB/Simulink基础操作就能按步骤把模型搭出来。我会把我在调试过程中踩过的坑一并写出来尤其是那些一跑就发散、一调就震荡的经典问题。1. 先看懂SMES是什么再谈建模1.1 用一个例子理解SMES的储能逻辑SMES的核心物理基础很简单电流流过线圈会建立磁场磁场本身携带能量。能量公式是 E ½LI² L是线圈电感I是流过线圈的电流。传统线圈因为有电阻电流会逐渐衰减能量存不住但超导线圈的直流电阻几乎为零只要维持低温环境线圈电流就可以长期保持能量也就一直“锁”在磁场里。我用一个生活化类比来解释超导线圈就像一个“超级飞轮”只不过飞轮储存的是旋转动能SMES储存的是磁场能量。飞轮靠惯性维持转动SMES靠零电阻维持电流。需要释放能量时让线圈电流减小需要吸收能量时让线圈电流增大。整个储能过程不涉及化学反应所以不存在电池老化问题充放电循环寿命极长响应速度还比化学电池快好几个数量级。1.2 建模之前先确定系统的硬件结构SMES系统从电网侧到线圈侧典型结构包含四个部分电网、功率变换系统PCS、超导线圈、低温与保护系统。在仿真里低温系统不用建模我们用等效模型替代即可重点放在电气主回路和控制环节上。功率变换系统通常采用两级结构。第一级是三相电压源型PWM变换器VSC负责连接交流电网第二级是双向DC-DC斩波器负责连接超导线圈和直流母线。为什么要两级因为VSC需要维持直流母线电压稳定而线圈电流又需要独立控制两级结构可以把这两件事解耦控制逻辑清晰工程实现也成熟。我给仿真模型设定的参数如下这些参数并不代表某个真实机组但结构上按实际SMES系统映射研究结论可以推广。部件参数数值交流电网线电压/频率380V / 50Hz直流母线额定电压800V直流母线电容C10mF超导线圈电感L5H超导线圈额定电流50A仿真比例缩小开关频率VSC/斩波器10kHz / 5kHz仿真采样步长离散控制器10μs线圈额定电流50A看起来不大配合5H的电感存储能量约6.25kJ。用于机理验证足够你可以按比例放大到实际机组级别。2. 数学模型把线圈和变换器翻译成微分方程2.1 超导线圈的等效模型在Simulink里超导线圈用一个大电感串联一个极小电阻来建模。超导状态下电阻理论上是零但仿真中如果不加任何电阻直流路径可能会出现数值奇异所以我习惯在参数里保留一个微欧级别的电阻比如1e-4Ω既不影响物理规律又能避免数值问题。线圈的电气方程是L * (dI_coil / dt) R * I_coil V_coil其中V_coil是斩波器施加在线圈两端的电压。当电网需要吸收能量时V_coil为正线圈电流上升需要释放能量时V_coil为负线圈电流下降。线圈吸收的瞬时功率为P_coil V_coil * I_coil LI_coil(dI_coil/dt)功率流动方向完全由V_coil极性决定。还有一点容易忽略线圈电流不能突变这是电感的天然属性。所以任何要求SMES“瞬时放出几千焦能量”的控制逻辑落实到电流指令上都必须有限斜率或限幅处理否则实际系统里线圈根本跟不上你的指令。这一点在后面控制器设计时我会再强调。2.2 三相PWM整流器的dq模型VSC是SMES连接交流系统的桥梁。三相坐标系下的VSC方程很直观但不好控制因为三相量都在正弦变化。工程上最常用的做法是把三相静止坐标系转换到dq同步旋转坐标系让基波分量变成直流量PI控制器才有发挥空间。Park变换后的VSC交流侧方程如下Ls * (disd / dt) -Rs * isd ω * Ls * isq usd - vsd Ls * (disq / dt) -Rs * isq - ω * Ls * isd usq - vsq其中usd、usq是电网电压的dq分量vsd、vsq是VSC交流侧电压的dq分量isd、isq是交流侧电流ω是电网角频率Ls和Rs是进线电感和等效电阻。这不是我随手写的它本质上是KVL在dq坐标系下的投影只要做完坐标变换剩下的就是代数操作。dq坐标系下还有一个好处有功功率主要由d轴电流决定无功功率主要由q轴电流决定。分析公式如下P 1.5 * (usd * isd usq * isq) Q 1.5 * (usq * isd - usd * isq)如果锁相环让q轴电压usq0那P就正比于usdisdQ就反比于usdisq。这就实现了有功和无功的解耦控制d轴电流管有功控制q轴电流管无功互不干扰。2.3 直流母线功率平衡关系直流母线电容上的电压变化本质是流入和流出母线的功率不平衡造成的。母线电容方程为C * (dVdc / dt) Idc_vsc - Idc_chopIdc_vsc是VSC注入母线的直流电流Idc_chop是斩波器从母线抽取的电流。VSC注入功率多母线电压就上升斩波器抽取功率多母线电压就下降。整个系统的能量流比很多人想象得复杂一点电网先通过VSC把交流电转成直流再通过斩波器把直流能量“泵”进线圈磁场里释放时倒过来。如果两级控制配合不好母线电压就会剧烈波动你说不清是电网问题还是线圈问题。这一段的建模是整个仿真的地基地基打不牢后面控制调得再精细也白搭。尤其是dq变换的符号约定不同教材正负方向定义不一样一旦弄错方向反馈极性就会反系统一闭环就崩。我在第三节会详细讲控制结构怎么和这个模型配合。3. 控制系统设计3.1 VSC双闭环控制结构VSC采用经典的双闭环控制外环是直流母线电压环内环是交流电流环。外环计算“母线电压离目标还差多少”输出d轴电流参考值内环快速跟踪这个参考值实际输出PWM占空比。我设计的控制规律如下isd_ref Kp_v * (Vdc_ref - Vdc) Ki_v * ∫(Vdc_ref - Vdc)dt vsd usd ω*Ls*isq - (Kp_i*(isd_ref - isd) Ki_i*∫(isd_ref - isd)dt) vsq usq - ω*Ls*isd - (Kp_i*(isq_ref - isq) Ki_i*∫(isq_ref - isq)dt)注意vsd、vsq的表达式里有 ωLsisq 和 -ωLsisd 两项这是d轴和q轴之间的耦合补偿项。如果没有这两项d轴电流变化会影响q轴q轴变化会影响d轴调试时表现为两个电流环互相牵制怎么调都调不出理想阶跃响应。加了补偿项两个电流环近似解耦调试难度瞬间下降。无功功率一般设定为0即isq_ref0让系统单位功率因数运行。如果你想仿真无功补偿能力直接给isq_ref一个阶跃就行后面的功率响应曲线会非常直观。3.2 斩波器控制超导线圈电流斩波器的任务是控制线圈电流跟踪参考值Iref。这个参考值怎么定通常由上层能量管理策略给出检测到电网需要支撑时Iref下降线圈放能电网有余量时Iref上升线圈吸能。斩波器电流环的原理是比较Iref和实际线圈电流误差经过PI调节器输出占空比D。D的正负决定了线圈两端电压的极性从而控制电流上升下降。我用的式子V_coil_ref Kp_chopper * (Iref - I_coil) Ki_chopper * ∫(Iref - I_coil)dt D V_coil_ref / Vdc 0.5为什么会加0.5偏移这是为了把双向斩波器的控制信号映射到PWM比较器的工作区间里。这个映射关系取决于你用的斩波器拓扑如果用的是全桥电路D的解释会略有不同但总体思路一致让平均输出电压等于V_coil_ref。3.3 PI参数初步整定思路很多人拿到PI就习惯性试凑我不反对试凑但效率太低。对电流环我习惯用带宽设计法。电流环控制对象近似为一阶惯性环节闭环带宽选在开关频率的1/10到1/5之间。假设开关频率10kHz电流环带宽取1kHz进线电感Ls取5mH则Kp_i Ls * 2π * f_bandwidth ≈ 0.005 * 6280 ≈ 31.4 Ki_i Kp_i * Rs / Ls ≈ 31.4 * 0.1 / 0.005 ≈ 628电压环带宽再比电流环低一个数量级取100Hz左右避免内外环之间产生频带重叠和交互振荡。斩波器电流环同理控制对象是大电感动态慢带宽可以取低一些200-300Hz即可PI参数根据L_coil和等效串联电阻计算。这里面有个我反复强调的原则内环必须比外环快5到10倍。内环慢、外环快就会出现“外环发出指令、内环还没跟上”的失控状态母线电压过冲大甚至直接发散。参数计算完之后再在Simulink里微调5%-15%基本就能稳定运行。4. Simulink仿真搭建实操4.1 顶层模型框架与模块选型打开Simulink新建模型我习惯先在脑子里把顶层分成几个子系统电网与变压器、VSC主电路、直流母线与斩波器、超导线圈、三相PLL与坐标变换、VSC控制器、斩波器控制器、测量与示波器。每个子系统封装成一个Subsystem内部线路不会干扰外部逻辑排查问题时可以逐级双击查看。主电路部分我建议用Simscape Electrical以前叫SimPowerSystems里的Universal Bridge模块来建模VSC开关管选IGBT/Diode。直流母线电容直接用Series RLC Branch参数设为电容。斩波器可以用两个IGBT构成半桥或者更简单一点用Controlled Voltage Source配合电流源来表示后者属于平均值模型避开高频开关细节仿真速度大幅提升。超导线圈用Series RLC Branch电感5H、电阻1e-4Ω即可。如果你对电力电子开关过程不感兴趣只想研究上层控制策略我强烈推荐第一版用平均值模型。等系统逻辑调通之后再换成全开关模型验证波形细节这样做能省掉很多时间。4.2 坐标变换与控制器模块坐标变换需要三相电压和电流作为输入。方法是用PLL模块锁定电网相位θ然后根据Park变换矩阵把abc转成dq。Simscape Electrical自带PLL模块直接用三相锁相环的作用是给坐标变换提供准确的相位基准相位不准dq分量算出来就是歪的控制自然不对。VSC控制器我用四个PI控制器模块加一堆数学运算搭的。之所以不推荐直接使用Simulink自带的“PID Controller”模块是因为在离散控制系统中我需要精确控制积分器的抗饱和逻辑。自定义PI模块内部如下err ref - feedback P_out Kp * err I_out Ki * err * Ts I_out saturate(I_out, Limit_min, Limit_max) out P_out I_out积分限幅特别重要。线圈电流参考值如果因为某种原因超出超导线圈的临界电流实际系统会发生失超仿真模型里则表现为电流飞升。所以我在斩波器电流环的输出和参考值两处都加了Saturation限幅模块限幅值取超导临界电流的80%作为保护裕量。PWM生成部分用PWM Generator模块载波频率设置成10kHz载波幅值范围映射到[-1,1]直接把调制波送入即可。调制波就是电压环和电流环算出来的vsd、vsq经过反Park变换后的三相调制信号。4.3 仿真参数与求解器配置这是最容易踩坑的地方。我首次搭建时直接在默认设置下运行结果系统发散得一塌糊涂。后来总结的经验是主电路开关器件必须用离散化求解器不能选默认的变步长ode45。我的配置如下配置项设置值说明求解器类型离散固定步长匹配数字控制器采样逻辑固定步长1e-6 s1μs1us可以覆盖10kHz开关过程仿真时长2s覆盖稳态暂态恢复全过程离散采样时间10μs控制器控制器和主电路解耦采样功率模块积分方式默认Tustin可避免部分数值振荡固定步长为什么是1μs因为开关频率10kHz一个开关周期100μs为了让PWM边沿解析足够精细步长至少要小于开关周期的百分之一1μs正好。你要是把步长放宽到10μs有些开关组合波形的细节会丢失仿真结果看着还行但换成实际硬件就会对不上。线圈初始电流也要设置。我建议把初始电流设为额定值50A这样系统启动时已经有储能可以直接观察SMES的放电支撑性能。从零电流启动也不是不行但那样前面几百毫秒都在给线圈充电你观察到的动态过程就不纯粹了。5. 仿真结果分析与验证5.1 典型工况电网电压跌落支撑设定工况系统稳态运行到t0.5s时电网电压突降20%持续0.3s后恢复。仿真目的就是看SMES能不能在电压跌落期间维持直流母线电压并向系统提供有功支撑。我在模型中用三相可编程电压源实现电压跌落PCC点接一个受控有功负荷。仿真结果里最核心的几条曲线是直流母线电压Vdc、线圈电流I_coil、VSC交流侧有功功率P、电网电压幅值Us。实测下来电压跌落瞬间VSC的功率平衡被打破母线电压有下降趋势。此时斩波器控制器检测到母线电压跌落立刻降低线圈电流指令线圈释放储能有功功率在几个毫秒内上升把母线电压“顶”回接近额定值。电压恢复后斩波器反向工作线圈电流回升到50A系统回到初始状态。整个过程没有切机没有低频振荡这就是SMES“毫秒级功率支撑”能力的直观体现。5.2 验证模型正确性的几个指标仿真跑通不等于仿真正确。我通常用三个指标来审计模型第一能量守恒。把VSC吸收的电能、线圈磁场能变化、线路损耗积分分别算出来校验能量守恒关系成立。用MATLAB的trapz函数对功率曲线积分最终差值不超过2%就算合格。第二功率响应时间。从功率指令阶跃到实际功率达到90%目标值这个时间应该在十几毫秒以内。如果响应明显迟滞先检查电流环带宽是否太低再检查线圈电流限幅是否触发。第三电流纹波。线圈电流在稳态时应该基本平滑纹波幅度控制在1%以内。如果纹波偏大说明斩波器开关频率太低或电感值太小。当然这跟实际SMES的大电感参数有关我在仿真里故意把电感设成5H纹波天然很小这符合超导线圈“大电感、慢电流变化”的特点。另外一个容易被忽略的验证点无功支撑能力。我做过一组测试保持有功不变把isq_ref从0阶跃到20A观察输出无功是否快速跟踪、d轴电流有没有串扰。由于我在控制器里加了dq解耦补偿q轴阶跃对d轴的扰动非常小双闭环的独立性验证良好。如果你的模型里d轴电流跟着q轴一起跳那十有八九是耦合补偿项出了问题。6. 常见问题排查与调试心得6.1 仿真一跑就发散问题出在哪我刚开始做SMES仿真时最常见的发散原因有三个。第一个是步长太大开关过程解析不足数值振荡越来越烈把步长1e-5改到1e-6一下就稳定了。第二个是初始条件不匹配线圈初始电流和控制器输出占空比不对应启动瞬间产生巨大冲击表现为Vdc和I_coil剧烈跳变。解决办法是让Iref的初始值和线圈初始电流一致同时让VSC双闭环的积分器初始值设为0。第三个是代数环。如果你从VSC测量模块直接取电流再反馈回同一模块的调制波输入端Simulink会提示代数环求解容易卡死。解决办法是给控制器输出加一个Memory或Unit Delay模块打断代数环或者用离散PID模块自带的“公式化”选项避免直接代数依赖。6.2 PI参数为什么总是调不稳调PI最忌讳一次改两三个参数根本不知道是谁在影响稳定性。我推荐的顺序是先断开外环给内环一个固定参考值只调内环PI调到阶跃无超调、无稳态误差然后接回外环在低增益下启动逐步增加外环Kp观察母线电压动态最后微调积分系数消除稳态误差。还有一点积分限幅和输出限幅必须匹配。如果积分器积分到特别大的值而输出被限幅卡住下一次误差反向时积分还要很久才能回退这就是“积分饱和”现象表现为大阶跃后系统恢复极慢。解决办法就是我前面说的给积分器本身加限幅并且适当让输出限幅略大于积分限幅防止两者互相憋住。6.3 开关模型选择的两难如果是纯粹研究控制策略平均值模型速度更快、更好调但如果要观察PWM谐波和电流纹波全开关模型才是真实反映。我建议分两阶段调策略用平均值模型出波形用全开关模型。全开关模型仿真速度确实慢2秒仿真时长在普通电脑上可能要跑十几分钟这时候合理设置数据记录范围也能省不少内存比如只在暂态窗口内记录高速波形其他时间只记录低速采样数据。6.4 我自己的一套排查顺序当仿真结果异常时我有一套固定的排查顺序这里分享给你1. 先看模型能不能“跑”而不是“跑得快”——发散先查步长和代数环 2. 看稳态值对不对——母线电压是否落在800V附近线圈电流是否稳态不飘 3. 看扰动响应——给定一个小阶跃看响应过程有没有异常振荡 4. 逐级断开——先断斩波器只调VSC再断VSC只调斩波器 5. 对比文献或简化模型——如果你能找到公开的SMES论文波形哪怕参数不同形状趋势也应该一致实践证明按照这个顺序排查绝大多数问题都能在半小时内定位。最忌讳的是看着屏幕上乱七八糟的波形瞎猜参数那不叫调试叫碰运气。最后再分享一个小技巧在搭建过程中每完成一个子系统就单独测试一个子系统。我先单独搭VSC接一个纯电阻负载看它能不能把直流母线稳在800V再单独搭斩波器看电流能不能跟踪方波指令最后才把两个子系统连起来联调。这样每一步出了问题你心里都清楚“问题只在本环节内”而不是满屏找错误。SMES这个课题本身并不神秘它的建模思路完全可以平移到其他储能设备上你就当是用一个高功率密度的储能系统练一遍从物理到数学再到仿真的完整链路这个经验绝对不亏。
超导磁能储存系统SMES的Simulink建模与仿真实现
做了几年电力电子仿真要说储能单元里最让我印象深刻的**超导磁能储存系统SMES**绝对算一个。名字听着高冷实际在Simulink里搭起来反而不复杂一个线圈、一套功率变换器、一组控制算法剩下全是数学处理。和锂电池、飞轮储能相比SMES的响应速度可以达到毫秒级功率密度极高非常适合电网暂态支撑、电能质量治理这种“瞬间要出力”的场景。这也是为什么我建议做电力系统方向的同学或者搞新能源并网的朋友把SMES建模当成一次很好的练手项目。本文我会从原理、数学模型、控制策略到Simulink实现把整个流程完整走一遍。你不需要提前掌握超导物理细节只需要有基本的电路和电力电子知识配合MATLAB/Simulink基础操作就能按步骤把模型搭出来。我会把我在调试过程中踩过的坑一并写出来尤其是那些一跑就发散、一调就震荡的经典问题。1. 先看懂SMES是什么再谈建模1.1 用一个例子理解SMES的储能逻辑SMES的核心物理基础很简单电流流过线圈会建立磁场磁场本身携带能量。能量公式是 E ½LI² L是线圈电感I是流过线圈的电流。传统线圈因为有电阻电流会逐渐衰减能量存不住但超导线圈的直流电阻几乎为零只要维持低温环境线圈电流就可以长期保持能量也就一直“锁”在磁场里。我用一个生活化类比来解释超导线圈就像一个“超级飞轮”只不过飞轮储存的是旋转动能SMES储存的是磁场能量。飞轮靠惯性维持转动SMES靠零电阻维持电流。需要释放能量时让线圈电流减小需要吸收能量时让线圈电流增大。整个储能过程不涉及化学反应所以不存在电池老化问题充放电循环寿命极长响应速度还比化学电池快好几个数量级。1.2 建模之前先确定系统的硬件结构SMES系统从电网侧到线圈侧典型结构包含四个部分电网、功率变换系统PCS、超导线圈、低温与保护系统。在仿真里低温系统不用建模我们用等效模型替代即可重点放在电气主回路和控制环节上。功率变换系统通常采用两级结构。第一级是三相电压源型PWM变换器VSC负责连接交流电网第二级是双向DC-DC斩波器负责连接超导线圈和直流母线。为什么要两级因为VSC需要维持直流母线电压稳定而线圈电流又需要独立控制两级结构可以把这两件事解耦控制逻辑清晰工程实现也成熟。我给仿真模型设定的参数如下这些参数并不代表某个真实机组但结构上按实际SMES系统映射研究结论可以推广。部件参数数值交流电网线电压/频率380V / 50Hz直流母线额定电压800V直流母线电容C10mF超导线圈电感L5H超导线圈额定电流50A仿真比例缩小开关频率VSC/斩波器10kHz / 5kHz仿真采样步长离散控制器10μs线圈额定电流50A看起来不大配合5H的电感存储能量约6.25kJ。用于机理验证足够你可以按比例放大到实际机组级别。2. 数学模型把线圈和变换器翻译成微分方程2.1 超导线圈的等效模型在Simulink里超导线圈用一个大电感串联一个极小电阻来建模。超导状态下电阻理论上是零但仿真中如果不加任何电阻直流路径可能会出现数值奇异所以我习惯在参数里保留一个微欧级别的电阻比如1e-4Ω既不影响物理规律又能避免数值问题。线圈的电气方程是L * (dI_coil / dt) R * I_coil V_coil其中V_coil是斩波器施加在线圈两端的电压。当电网需要吸收能量时V_coil为正线圈电流上升需要释放能量时V_coil为负线圈电流下降。线圈吸收的瞬时功率为P_coil V_coil * I_coil LI_coil(dI_coil/dt)功率流动方向完全由V_coil极性决定。还有一点容易忽略线圈电流不能突变这是电感的天然属性。所以任何要求SMES“瞬时放出几千焦能量”的控制逻辑落实到电流指令上都必须有限斜率或限幅处理否则实际系统里线圈根本跟不上你的指令。这一点在后面控制器设计时我会再强调。2.2 三相PWM整流器的dq模型VSC是SMES连接交流系统的桥梁。三相坐标系下的VSC方程很直观但不好控制因为三相量都在正弦变化。工程上最常用的做法是把三相静止坐标系转换到dq同步旋转坐标系让基波分量变成直流量PI控制器才有发挥空间。Park变换后的VSC交流侧方程如下Ls * (disd / dt) -Rs * isd ω * Ls * isq usd - vsd Ls * (disq / dt) -Rs * isq - ω * Ls * isd usq - vsq其中usd、usq是电网电压的dq分量vsd、vsq是VSC交流侧电压的dq分量isd、isq是交流侧电流ω是电网角频率Ls和Rs是进线电感和等效电阻。这不是我随手写的它本质上是KVL在dq坐标系下的投影只要做完坐标变换剩下的就是代数操作。dq坐标系下还有一个好处有功功率主要由d轴电流决定无功功率主要由q轴电流决定。分析公式如下P 1.5 * (usd * isd usq * isq) Q 1.5 * (usq * isd - usd * isq)如果锁相环让q轴电压usq0那P就正比于usdisdQ就反比于usdisq。这就实现了有功和无功的解耦控制d轴电流管有功控制q轴电流管无功互不干扰。2.3 直流母线功率平衡关系直流母线电容上的电压变化本质是流入和流出母线的功率不平衡造成的。母线电容方程为C * (dVdc / dt) Idc_vsc - Idc_chopIdc_vsc是VSC注入母线的直流电流Idc_chop是斩波器从母线抽取的电流。VSC注入功率多母线电压就上升斩波器抽取功率多母线电压就下降。整个系统的能量流比很多人想象得复杂一点电网先通过VSC把交流电转成直流再通过斩波器把直流能量“泵”进线圈磁场里释放时倒过来。如果两级控制配合不好母线电压就会剧烈波动你说不清是电网问题还是线圈问题。这一段的建模是整个仿真的地基地基打不牢后面控制调得再精细也白搭。尤其是dq变换的符号约定不同教材正负方向定义不一样一旦弄错方向反馈极性就会反系统一闭环就崩。我在第三节会详细讲控制结构怎么和这个模型配合。3. 控制系统设计3.1 VSC双闭环控制结构VSC采用经典的双闭环控制外环是直流母线电压环内环是交流电流环。外环计算“母线电压离目标还差多少”输出d轴电流参考值内环快速跟踪这个参考值实际输出PWM占空比。我设计的控制规律如下isd_ref Kp_v * (Vdc_ref - Vdc) Ki_v * ∫(Vdc_ref - Vdc)dt vsd usd ω*Ls*isq - (Kp_i*(isd_ref - isd) Ki_i*∫(isd_ref - isd)dt) vsq usq - ω*Ls*isd - (Kp_i*(isq_ref - isq) Ki_i*∫(isq_ref - isq)dt)注意vsd、vsq的表达式里有 ωLsisq 和 -ωLsisd 两项这是d轴和q轴之间的耦合补偿项。如果没有这两项d轴电流变化会影响q轴q轴变化会影响d轴调试时表现为两个电流环互相牵制怎么调都调不出理想阶跃响应。加了补偿项两个电流环近似解耦调试难度瞬间下降。无功功率一般设定为0即isq_ref0让系统单位功率因数运行。如果你想仿真无功补偿能力直接给isq_ref一个阶跃就行后面的功率响应曲线会非常直观。3.2 斩波器控制超导线圈电流斩波器的任务是控制线圈电流跟踪参考值Iref。这个参考值怎么定通常由上层能量管理策略给出检测到电网需要支撑时Iref下降线圈放能电网有余量时Iref上升线圈吸能。斩波器电流环的原理是比较Iref和实际线圈电流误差经过PI调节器输出占空比D。D的正负决定了线圈两端电压的极性从而控制电流上升下降。我用的式子V_coil_ref Kp_chopper * (Iref - I_coil) Ki_chopper * ∫(Iref - I_coil)dt D V_coil_ref / Vdc 0.5为什么会加0.5偏移这是为了把双向斩波器的控制信号映射到PWM比较器的工作区间里。这个映射关系取决于你用的斩波器拓扑如果用的是全桥电路D的解释会略有不同但总体思路一致让平均输出电压等于V_coil_ref。3.3 PI参数初步整定思路很多人拿到PI就习惯性试凑我不反对试凑但效率太低。对电流环我习惯用带宽设计法。电流环控制对象近似为一阶惯性环节闭环带宽选在开关频率的1/10到1/5之间。假设开关频率10kHz电流环带宽取1kHz进线电感Ls取5mH则Kp_i Ls * 2π * f_bandwidth ≈ 0.005 * 6280 ≈ 31.4 Ki_i Kp_i * Rs / Ls ≈ 31.4 * 0.1 / 0.005 ≈ 628电压环带宽再比电流环低一个数量级取100Hz左右避免内外环之间产生频带重叠和交互振荡。斩波器电流环同理控制对象是大电感动态慢带宽可以取低一些200-300Hz即可PI参数根据L_coil和等效串联电阻计算。这里面有个我反复强调的原则内环必须比外环快5到10倍。内环慢、外环快就会出现“外环发出指令、内环还没跟上”的失控状态母线电压过冲大甚至直接发散。参数计算完之后再在Simulink里微调5%-15%基本就能稳定运行。4. Simulink仿真搭建实操4.1 顶层模型框架与模块选型打开Simulink新建模型我习惯先在脑子里把顶层分成几个子系统电网与变压器、VSC主电路、直流母线与斩波器、超导线圈、三相PLL与坐标变换、VSC控制器、斩波器控制器、测量与示波器。每个子系统封装成一个Subsystem内部线路不会干扰外部逻辑排查问题时可以逐级双击查看。主电路部分我建议用Simscape Electrical以前叫SimPowerSystems里的Universal Bridge模块来建模VSC开关管选IGBT/Diode。直流母线电容直接用Series RLC Branch参数设为电容。斩波器可以用两个IGBT构成半桥或者更简单一点用Controlled Voltage Source配合电流源来表示后者属于平均值模型避开高频开关细节仿真速度大幅提升。超导线圈用Series RLC Branch电感5H、电阻1e-4Ω即可。如果你对电力电子开关过程不感兴趣只想研究上层控制策略我强烈推荐第一版用平均值模型。等系统逻辑调通之后再换成全开关模型验证波形细节这样做能省掉很多时间。4.2 坐标变换与控制器模块坐标变换需要三相电压和电流作为输入。方法是用PLL模块锁定电网相位θ然后根据Park变换矩阵把abc转成dq。Simscape Electrical自带PLL模块直接用三相锁相环的作用是给坐标变换提供准确的相位基准相位不准dq分量算出来就是歪的控制自然不对。VSC控制器我用四个PI控制器模块加一堆数学运算搭的。之所以不推荐直接使用Simulink自带的“PID Controller”模块是因为在离散控制系统中我需要精确控制积分器的抗饱和逻辑。自定义PI模块内部如下err ref - feedback P_out Kp * err I_out Ki * err * Ts I_out saturate(I_out, Limit_min, Limit_max) out P_out I_out积分限幅特别重要。线圈电流参考值如果因为某种原因超出超导线圈的临界电流实际系统会发生失超仿真模型里则表现为电流飞升。所以我在斩波器电流环的输出和参考值两处都加了Saturation限幅模块限幅值取超导临界电流的80%作为保护裕量。PWM生成部分用PWM Generator模块载波频率设置成10kHz载波幅值范围映射到[-1,1]直接把调制波送入即可。调制波就是电压环和电流环算出来的vsd、vsq经过反Park变换后的三相调制信号。4.3 仿真参数与求解器配置这是最容易踩坑的地方。我首次搭建时直接在默认设置下运行结果系统发散得一塌糊涂。后来总结的经验是主电路开关器件必须用离散化求解器不能选默认的变步长ode45。我的配置如下配置项设置值说明求解器类型离散固定步长匹配数字控制器采样逻辑固定步长1e-6 s1μs1us可以覆盖10kHz开关过程仿真时长2s覆盖稳态暂态恢复全过程离散采样时间10μs控制器控制器和主电路解耦采样功率模块积分方式默认Tustin可避免部分数值振荡固定步长为什么是1μs因为开关频率10kHz一个开关周期100μs为了让PWM边沿解析足够精细步长至少要小于开关周期的百分之一1μs正好。你要是把步长放宽到10μs有些开关组合波形的细节会丢失仿真结果看着还行但换成实际硬件就会对不上。线圈初始电流也要设置。我建议把初始电流设为额定值50A这样系统启动时已经有储能可以直接观察SMES的放电支撑性能。从零电流启动也不是不行但那样前面几百毫秒都在给线圈充电你观察到的动态过程就不纯粹了。5. 仿真结果分析与验证5.1 典型工况电网电压跌落支撑设定工况系统稳态运行到t0.5s时电网电压突降20%持续0.3s后恢复。仿真目的就是看SMES能不能在电压跌落期间维持直流母线电压并向系统提供有功支撑。我在模型中用三相可编程电压源实现电压跌落PCC点接一个受控有功负荷。仿真结果里最核心的几条曲线是直流母线电压Vdc、线圈电流I_coil、VSC交流侧有功功率P、电网电压幅值Us。实测下来电压跌落瞬间VSC的功率平衡被打破母线电压有下降趋势。此时斩波器控制器检测到母线电压跌落立刻降低线圈电流指令线圈释放储能有功功率在几个毫秒内上升把母线电压“顶”回接近额定值。电压恢复后斩波器反向工作线圈电流回升到50A系统回到初始状态。整个过程没有切机没有低频振荡这就是SMES“毫秒级功率支撑”能力的直观体现。5.2 验证模型正确性的几个指标仿真跑通不等于仿真正确。我通常用三个指标来审计模型第一能量守恒。把VSC吸收的电能、线圈磁场能变化、线路损耗积分分别算出来校验能量守恒关系成立。用MATLAB的trapz函数对功率曲线积分最终差值不超过2%就算合格。第二功率响应时间。从功率指令阶跃到实际功率达到90%目标值这个时间应该在十几毫秒以内。如果响应明显迟滞先检查电流环带宽是否太低再检查线圈电流限幅是否触发。第三电流纹波。线圈电流在稳态时应该基本平滑纹波幅度控制在1%以内。如果纹波偏大说明斩波器开关频率太低或电感值太小。当然这跟实际SMES的大电感参数有关我在仿真里故意把电感设成5H纹波天然很小这符合超导线圈“大电感、慢电流变化”的特点。另外一个容易被忽略的验证点无功支撑能力。我做过一组测试保持有功不变把isq_ref从0阶跃到20A观察输出无功是否快速跟踪、d轴电流有没有串扰。由于我在控制器里加了dq解耦补偿q轴阶跃对d轴的扰动非常小双闭环的独立性验证良好。如果你的模型里d轴电流跟着q轴一起跳那十有八九是耦合补偿项出了问题。6. 常见问题排查与调试心得6.1 仿真一跑就发散问题出在哪我刚开始做SMES仿真时最常见的发散原因有三个。第一个是步长太大开关过程解析不足数值振荡越来越烈把步长1e-5改到1e-6一下就稳定了。第二个是初始条件不匹配线圈初始电流和控制器输出占空比不对应启动瞬间产生巨大冲击表现为Vdc和I_coil剧烈跳变。解决办法是让Iref的初始值和线圈初始电流一致同时让VSC双闭环的积分器初始值设为0。第三个是代数环。如果你从VSC测量模块直接取电流再反馈回同一模块的调制波输入端Simulink会提示代数环求解容易卡死。解决办法是给控制器输出加一个Memory或Unit Delay模块打断代数环或者用离散PID模块自带的“公式化”选项避免直接代数依赖。6.2 PI参数为什么总是调不稳调PI最忌讳一次改两三个参数根本不知道是谁在影响稳定性。我推荐的顺序是先断开外环给内环一个固定参考值只调内环PI调到阶跃无超调、无稳态误差然后接回外环在低增益下启动逐步增加外环Kp观察母线电压动态最后微调积分系数消除稳态误差。还有一点积分限幅和输出限幅必须匹配。如果积分器积分到特别大的值而输出被限幅卡住下一次误差反向时积分还要很久才能回退这就是“积分饱和”现象表现为大阶跃后系统恢复极慢。解决办法就是我前面说的给积分器本身加限幅并且适当让输出限幅略大于积分限幅防止两者互相憋住。6.3 开关模型选择的两难如果是纯粹研究控制策略平均值模型速度更快、更好调但如果要观察PWM谐波和电流纹波全开关模型才是真实反映。我建议分两阶段调策略用平均值模型出波形用全开关模型。全开关模型仿真速度确实慢2秒仿真时长在普通电脑上可能要跑十几分钟这时候合理设置数据记录范围也能省不少内存比如只在暂态窗口内记录高速波形其他时间只记录低速采样数据。6.4 我自己的一套排查顺序当仿真结果异常时我有一套固定的排查顺序这里分享给你1. 先看模型能不能“跑”而不是“跑得快”——发散先查步长和代数环 2. 看稳态值对不对——母线电压是否落在800V附近线圈电流是否稳态不飘 3. 看扰动响应——给定一个小阶跃看响应过程有没有异常振荡 4. 逐级断开——先断斩波器只调VSC再断VSC只调斩波器 5. 对比文献或简化模型——如果你能找到公开的SMES论文波形哪怕参数不同形状趋势也应该一致实践证明按照这个顺序排查绝大多数问题都能在半小时内定位。最忌讳的是看着屏幕上乱七八糟的波形瞎猜参数那不叫调试叫碰运气。最后再分享一个小技巧在搭建过程中每完成一个子系统就单独测试一个子系统。我先单独搭VSC接一个纯电阻负载看它能不能把直流母线稳在800V再单独搭斩波器看电流能不能跟踪方波指令最后才把两个子系统连起来联调。这样每一步出了问题你心里都清楚“问题只在本环节内”而不是满屏找错误。SMES这个课题本身并不神秘它的建模思路完全可以平移到其他储能设备上你就当是用一个高功率密度的储能系统练一遍从物理到数学再到仿真的完整链路这个经验绝对不亏。
