飞轮储能这个方向,这几年在电网调频、UPS、轨道交通能量回收这些场景里越来越热。做仿真模型的时候,最常见也最核心的一套方案就是用永磁同步电机当能量转换核心,通过控制电机来实现飞轮的充放电。我见过不少人卡在这个控制策略的搭建上,尤其是充电外环转速、内环电流这套双闭环逻辑,看起来简单,真要把参数调稳、把充放电切换做顺,还是有不少细节坑。这篇把我自己的建模和调试经验完整拆开讲讲。
1. 项目核心思路与整体方案拆解
1.1 为什么飞轮储能系统选择永磁同步电机
飞轮储能系统的本质,是用一个高速旋转的飞轮转子把电能转化成机械能存起来,需要的时候再把机械能转化回电能。这套系统里,电机既是电动机也是发电机,充电的时候电动运行把飞轮加速,放电的时候发电运行把飞轮的动能变成电能送出去。整个能量转换链路里,电机和它的控制器决定了系统的效率上限、动态响应速度和稳定边界。
选永磁同步电机而不是异步电机,核心原因有三点:第一是效率,永磁电机的转子没有铜耗和铝耗,励磁损耗几乎为零,高效率区间宽,这对储能系统来说很关键,飞轮储能本来拼的就是循环效率;第二是功率密度,同样的体积和重量,永磁电机能输出的功率和转矩更大,飞轮系统为了减小风阻损耗,通常会把电机和飞轮封装在真空腔体里,体积约束很紧;第三是控制精度,永磁同步电机的转子磁场位置可以通过编码器精确获得,矢量控制下电流环和转速环都能做到很高的控制带宽,这对频繁充放电切换的动态性能是刚需。
1.2 飞轮储能充放电控制的双闭环结构选型
这套模型采用的控制架构是很标准的转速-电流双闭环结构。充电过程外环控制转速,内环控制电流,这是大部分飞轮储能充电控制的基本骨架。外环转速环根据目标转速和实际转速的偏差,输出一个转矩电流指令(q轴电流参考值),内环电流环跟踪这个指令,通过电压源逆变器输出电压矢量,驱动电机加速。放电的时候,控制目标从转速维持变成功率输出或母线电压维持,控制逻辑需要切换。
整体系统框图如下:
转速指令 → [转速环PI] → iq_ref → [电流环PI] → dq电压指令 → SVPWM → 逆变器 → PMSM ↑ ↓ 转速反馈 ← 编码器 电流采样 → iq/id反馈这套结构的好处是动静分离,转速环负责宏观的能量管理,电流环负责微观的转矩快速响应。电流环带宽通常整定到几百赫兹,转速环带宽在几十赫兹的量级,两者差一个数量级以上,才能保证不互相耦合振荡。
放电模式的切换逻辑也在这套模型里做了处理,当系统收到放电指令后,外环从转速环切换为直流母线电压环,控制目标从飞轮转速变为母线电压恒定,内环仍然控制电流。这个时候电机会自然进入发电工况,转速下降,动能转化为电能输出。
1.3 这套控制方案解决了什么问题
飞轮储能最怕的就是充放电切换瞬间的大电流冲击和转矩突变,直接表现就是母线电压跌落、电流尖峰、甚至机械轴系扭振。双闭环控制通过两个层面的协调把这些问题控制在合理范围:
转速环让转速按设定好的斜坡或S型曲线平滑升降,避免了直接给全转矩导致飞轮转子承受过大的剪切应力。飞轮的储能量和转速的平方成正比,转速控制得稳,能量管理就稳。电流环则把电磁转矩限制在设定范围内,防止大电流冲击损坏功率器件和永磁体退磁。
另外,这套模型里面把电机模型、逆变器模型、控制算法全部放进了Simulink环境,可以方便地验证不同工况下的控制性能。我在搭这个模型之前,也参考过不少现成的飞轮储能仿真例子,但很多模型要么控制策略太简单,没有完整的充放电切换逻辑;要么电机模型用理想模型替换,控制效果失真。这套模型好就好在完整保留了永磁同步电机的电磁动态和逆变器的开关特性,控制算法也按实际工程的做法来设计,仿真结果有说服力。
2. 核心控制参数设计与关键细节
2.1 永磁同步电机在Simulink中的建模方式选择
Simulink里做永磁同步电机的建模有三条路:用电气专业的永磁同步电机模块、用Simscape Electrical里的物理模型、或者自己搭状态方程。
我在这套模型里用的是电气专业(Simscape Electrical,老版本叫SimPowerSystems或者Specialized Power Systems)的永磁同步电机模块,配合三相两电平逆变器和直流母线构成主电路。这个选择的原因很简单:电机模块内部已经考虑了凸极效应、永磁体磁场谐波、阻尼绕组等实际因素,比自己搭状态方程准确得多。而且逆变器的IGBT开关动作可以通过PWM信号驱动的MOSFET/IGBT桥臂来建模,比用理想的受控电压源替代要真实得多,能观察到电枢电流的谐波成分和开关纹波。
用基于物理的建模有个问题是仿真速度慢,尤其是开关频率高的时候。我试验过Continuous和Discrete两种模式,连续模式下仿真精度高但速度较慢,离散模式下需要设定合理的采样时间,推荐用1e-6秒的步长和离散求解器,速度提升明显,而且便于后续部署到实时仿真平台。
2.2 转速环和电流环PI参数的计算过程
双闭环控制的核心就是两个PI调节器。我直接给出我这里用的参数整定过程,整个过程全部可以套公式算,不需要靠感觉调。
电流环采用PI控制,其控制对象是永磁同步电机的定子绕组,简化为一阶惯性环节。电流环的目标是快、准、稳,因此采用典型I型系统整定。以d轴电流环为例,给定电流环PI参数为Kp_i和Ki_i,电流环控制对象传递函数近似为:
G_i(s) = (Kp_i * s + Ki_i) / s * (1 / (Ls + R))
其中L是d轴电感(表贴式PMSM中Ld=Lq),R是定子电阻。按典型I型系统整定,取Kp_i / Ki_i = L / R,即零极点对消后,电流环开环增益为Kp_i/L,要求闭环带宽为ω_ci。最终参数计算公式:
Kp_i = ω_ci * L,Ki_i = ω_ci * R
我这里电机的参数是:L = 0.000835H(8.35mH),R = 0.0325Ω,等效到标幺值后按p.u.系统整定更方便。实际在模型里我用的是标幺值系统,基值电流取峰值电流,基值电压取额定相电压峰值。电流环带宽ω_ci取500Hz,即3140rad/s:
Kp_i = 3140 * 0.000835 = 2.62(标幺约等于0.52) Ki_i = 3140 * 0.0325 = 102
转速环的整定思路类似,但需要把电机运动方程考虑进去,控制对象变成电流环闭环加上机械转动惯量。转速环带宽ω_cn一般取电流环带宽的1/5到1/10,我这里取50Hz,即314rad/s。转速环PI参数计算公式:
Kp_n = J * ω_cn / (1.5 * p * Ψf) Ki_n = Kp_n * ω_cn / 5
其中J是飞轮加电机的总转动惯量,我这里飞轮转子和电机的合成转动惯量约为0.052kg·m²,p是极对数(这里取4),Ψf是永磁体磁链(约0.192Wb)。代入:
Kp_n = 0.052 * 314 / (1.5 * 4 * 0.192) = 7.09 Ki_n = 7.09 * 314 / 5 = 445.7
实际上在Simulink里做标幺化之后,PI参数要除以对应的基值折算系数,但比例关系不变。初版模型里就按这个计算值先填进去,出来的波形基本稳,然后微调。
提示:电流环带宽不能取太高,受限于逆变器开关频率,一般不超过开关频率的1/10。如果开关频率是10kHz,电流环带宽取500Hz-1kHz是合理的。超过这个范围,数字控制延迟会带来相位裕度不足,系统就会振荡。
2.3 SVPWM调制与坐标变换的工程细节
模型里电流环输出的dq轴电压指令,需要经过反Park变换和SVPWM调制才能生成逆变器的开关信号。SVPWM的实现逻辑我直接用了Simulink里的SVPWM模块,但在参数设置上有几个点值得注意:
直流母线电压的利用率:SVPWM的线性调制范围比SPWM高15%,最大输出相电压基波幅值可以达到Udc/√3。这意味着同样的直流母线电压,SVPWM能让电机跑得更高转速,这对飞轮储能来说很有意义,因为飞轮转速通常很高。
矢量切换的零矢量分配:SVPWM采用七段式还是五段式调制,对开关损耗和电流谐波影响很大。七段式每个开关周期内零矢量对称分布,电流谐波抑制效果更好,但开关损耗大;五段式通过减少开关次数降低损耗,但谐波略大。我的模型里用的是七段式,追求控制精度优先。
过调制:当电机运行在弱磁区域时,电压指令超过线性调制范围,如果不过处理直接限幅,输出波形会畸变。模型里做了过调制处理,把电压矢量限幅到SVPWM的最大线性范围(即Udc/√3对应的圆内),同时保留一个弱磁控制器来调整d轴电流,确保高速区也能稳定工作。
坐标变换环节看起来简单,但转子位置角对齐关系到整个控制能否建立。我习惯把电机的零位对齐到A相绕组轴线,然后通过编码器信号(模型里用增量编码器模块模拟)获得转子位置,再换算成电角度。有一点容易被忽略:如果电机的极对数是4,机械角度和电角度的换算系数是4,在反馈环节漏乘这个系数,电流环直接发散,这个坑我踩过。
2.4 飞轮转动惯量与转速环动态特性的匹配
飞轮储能系统的机械时间常数通常很大,因为飞轮的转动惯量不小,转速响应相对电流会很慢。这个特性既是好事也是需要处理的问题。
好处在于,转速环对被控对象参数变化不敏感,转动惯量的测量误差对转速环稳定性的影响不大,PI参数有较大的裕量。坏处在于,如果按转速环带宽50Hz去整定,会出现转速响应过快、电流限幅导致加速转矩不够的情况。
所以我在模型里对转速环做了两个处理。第一,转速指令发生器采用斜坡/梯形波,限制转速变化率,让转速环的给定变化平缓。这个限速率的做法在工程上等同于给飞轮系统一个最大加速度限制,保护机械轴系。第二,转速环PI输出端串联了一个限幅器,把iq指令幅值限制在逆变器和电机的最大允许电流范围内。限幅值按电机峰值转矩电流设置,一般取额定电流的1.2到1.5倍。
这两个处理合在一起,相当于在原有双闭环外面又加了一层动态特性和机械约束协调层。很多初学者只盯着PI参数调,忽略了限幅和给定整形,结果仿真出来转速环不稳定,电流波形全是限幅饱和的方波,这个现象我后面在常见问题里详细讲。
3. 仿真模型搭建与充放电实现
3.1 主电路与电机参数的设置
模型的主电路结构是:直流电压源(模拟母线或是前级变换器输出)→ 三相两电平逆变器 → 永磁同步电机 → 飞轮负载(转动惯量+摩擦阻尼)。
直流母线电压我这里设的是700V,这是一个比较典型的中小型飞轮储能系统母线电压等级。逆变器的IGBT导通压降设为1.8V,开关频率10kHz,死区时间2μs。这些参数都按照实际工程器件水平来设,不能设成理想开关,否则仿真结果偏乐观,后续做硬件在环会有落差。
永磁同步电机参数按一台额定功率约55kW、额定转速12000rpm的高速电机来设置,具体参数如下:
| 参数 | 数值 | 说明 |
|---|---|---|
| 定子电阻R | 0.0325Ω | 常温下的绕组电阻 |
| d轴电感Ld | 0.000835H | 表贴式结构,Ld=Lq |
| q轴电感Lq | 0.000835H | 与Ld相等 |
| 永磁体磁链 | 0.192Wb | 以转子位置0度时A相反电动势峰值为准 |
| 极对数 | 4 | 高速电机常用4-8极 |
| 转动惯量 | 0.052kg·m² | 电机转子+飞轮折算后的合成惯量 |
这里我需要提醒一下,不同文献里永磁同步电机的参数基准可能不一样,有的用线电压有效值,有的用相电压峰值。从电机厂家拿到的参数表,需要先统一化到Simulink模块使用的坐标系中。最简单的验证办法:空载反电动势仿真结果和理论值对比,设一个500rpm的转速,测A相反电动势波形幅值,和0.19242π*500/60对比一下,误差在3%以内说明参数输入正确。
3.2 充电过程的转速外环电流内环控制实现
充电过程的控制目标是把飞轮从静止加速到额定转速,并在达到目标转速后维持转速不变,此时系统处于待机状态,飞轮以额定转速空转,仅克服风阻摩擦损耗维持转速。
转速外环的输入是转速指令与实际转速的偏差,输出是q轴电流指令。具体到Simulink里的实现:
转速指令发生器:用阶跃+斜坡限制模块,设定加速度限制为500rpm/s,从0加速到12000rpm需要24秒。这个加速度限制是根据电机最大转矩电流和飞轮转动惯量计算出来的:最大加速转矩大约140Nm,折算到转速变化率就是约600rpm/s,留10%裕量后取500rpm/s。
转速环PI:含积分限幅和输出限幅,iq指令上限设为200A(对应约280Nm的峰值转矩),下限-200A。
电流内环:d轴电流指令恒为0(表贴式PMSM最大转矩电流比控制退化为id=0控制),q轴电流跟踪转速环输出。两个电流环各自使用PI控制器,输出为dq轴电压。
解耦补偿:在电流环输出处加入交叉解耦项,补偿dq轴之间的耦合电动势。解耦项的计算方法:Vd_comp = -ω_e * Lq * iq,Vq_comp = ω_e * (Ld * id + Ψf)。这一项不加的话,高速时dq轴电流控制会互相干扰,动态响应会明显变差。
整套逻辑跑下来,充电过程三相电流波形呈逐步增大的正弦波包络,转速线性上升,转矩电流维持在设定值附近,直流母线功率从零爬升到稳定值再下降,整体过程清晰稳定。
3.3 放电过程的控制切换与能量回馈实现
放电过程是整个模型里最容易出问题的地方。控制目标从转速恒定切换到母线电压恒定,本质上是从调速模式切到电压控制模式,两者需要平滑衔接。
我的实现方案是:把转速环的PI控制器旁路,电机进入发电模式,外环改为直流母线电压环。电压环的输入是母线电压给定与实际母线电压的偏差,输出是q轴电流指令(发电工况下iq为负值),内环电流环仍然不变,跟踪这个负的电流指令,电机转矩反向,转速下降,动能转化为电能回馈到直流母线。
切换的关键细节是切换瞬间的电流指令缓冲。当发电控制使能时,电压环PI的输出初始值不能从零开始,而应该从当前iq值开始,否则控制器积分器充电过程会造成数十毫秒的电流断档,母线电压会瞬间跌落。我在模型里用了一个带外部初值设定的积分器,切换信号到来时,把积分初值设置为切换前的iq值,这样就实现了无扰切换。
放电工况的功率大小由负载决定。模型中用了一个可变功率负载来模拟功率吸收需求,母线电压闭环会自动增大放电电流来维持电压。如果飞轮转速下降到最低允许转速(模型设定为6000rpm),触发欠速保护,停止放电,飞轮系统进入待机状态。这个最低转速与最大放电深度相关,放电深度越深,能量释放越多,但转速过低会导致母线电压无法维持,需要提前切出。
3.4 转速-功率协调与S型速度曲线规划
基础的转速斜坡指令可以做,但充放电过程更柔顺的做法是用S型速度曲线。原因在于飞轮系统的机械轴承和轴系在加速度突变时会有冲击,长期运行会造成疲劳损伤。
我在模型的速度指令发生器里做了一个简单的S型曲线规划:加速过程中先以线性增加的加速度启动,中间恒加速度,末端加速度线性减小到零。对应到转速曲线上就是一段S型过渡。这个曲线可以通过Simulink里的信号生成模块实现,或者直接用Matlab函数写一段数值积分。
S型曲线的好处不仅是机械冲击小,对电网侧的能量冲击也小。充电开始时如果直接给满转矩,直流母线相当于短时功率冲击,对前端整流器和电网都不友好。而S型曲线让功率缓慢爬升,对系统友好得多。
我做了一个简化处理:设定加速时间为24秒,但把前3秒作为加速度递增阶段,后3秒作为加速度递减阶段。仿真波形上可以看到转速指令在起始和结束阶段呈现出平滑的弧度,而不是折线拐角。这个细节虽然不影响控制稳定性,但对后续的并网功率质量很有意义。
4. 常见问题与排查技巧
4.1 电流环发散或振荡怎么办
最典型的现象是仿真开始后iq或者id电流振荡发散,直接飞掉。排查顺序是这样的:
先看坐标变换方向,尤其是转子位置角的正负号。把d轴方向和转子磁极方向对齐,如果符号反了,电流反馈就变成正反馈,系统必发散。验证方法:给一个常数iq指令,看电机是否向预期方向加速,同时看id是否保持在零附近。
再看电流环PI参数,如果计算出来的Kp_i和Ki_i是按照连续域算的,而Simulink里用的是离散PI模块,需要先做离散化转换。直接用连续域的增益放到离散控制器里,采样频率不够高时系统会振荡。
最后看解耦项,有没有加对。解耦项的符号错了会引入额外的正反馈。常见错误是交叉项漏乘转速电角速度,导致低速时正常、高速时振荡。
4.2 转速环响应慢且超调
转速环的PI参数按带宽50Hz计算,但实际效果发现转速跟踪慢,且到达目标转速时超调量达到3%-5%。这个问题多半是电流环限幅提前饱和了。
加速过程中,转速误差大,转速环PI输出想要一个很大的iq指令,但限幅器把它截在200A,这时候转速环相当于开环控制,等转速接近目标时误差变小,PI才从饱和状态退出来。退出过程中积分器还在持续累积,所以会产生超调。
解决办法有两个:一是给转速环PI的积分器做抗积分饱和设计,当输出达到限幅值时冻结积分器;二是在转速接近目标时提前切换控制策略,比如在目标转速的95%开始减小加速度。我在模型里用的是积分限幅+条件冻结,切换后超调量压到1%以内。
4.3 充电到放电切换瞬间母线电压跌落
这个问题在不少做飞轮储能仿真的朋友那里反复出现。现象是母线电压在切换后的前50ms内跌落了20V以上,然后才慢慢恢复。原因是切换时刻电压环PI的初始输出和当前工况不匹配。
调试方法:把放电开始的时刻记录下来,观察切换前瞬间的iq值、母线电压和直流负载功率。然后把电压环PI的积分初始值设置成切换前iq对应的值(或者切换前的实际占空比对应的等效值),切换就能做到基本无扰。另一种做法是切换后短时间内保持切换前的电流指令,用软启动的方式让电压环慢慢接管,也能改善。
4.4 SVPWM电压指令饱和
高速运行时,反电动势接近直流母线电压能输出的极限,电流环输出的dq电压指令会超出SVPWM的线性调制范围,导致实际输出电压不够,电流跟踪失败。这个现象通常发生在额定转速以上。
处理办法是给电流环输出电压加限幅,保证电压矢量幅值在最大线性范围内。同时需要考虑是否需要弱磁。如果飞轮的工作转速不超过额定转速,不加弱磁也能跑;如果要超速运行到120%额定转速,就必须引入弱磁控制。我的模型里加了一个简单的弱磁控制器:当电压指令超过限幅值时,减少d轴电流指令(负值),通过d轴去磁效应来降低反电动势。
弱磁控制器在调试时要格外小心,负的d轴电流很容易把电流限幅值占满,反而挤占q轴电流空间。一般设一个优先级,在总电流限幅的约束下,q轴电流优先,弱磁电流在剩余空间里调整。
4.5 仿真速度过慢的优化经验
高速PMSM仿真,尤其是带开关器件和PWM调制的时候,仿真速度很容易让人崩溃。一次完整的充放电仿真要模拟几十秒的物理时间,如果步长太小,几个小时都跑不完。
我优化之后的配置是:功率电路和解算器用离散模式,采样时间设为1μs(对应10kHz开关频率的PWM分辨率和足够的电流环带宽),机械部分(转速、位置)用较慢的采样率模块单独处理。Simulink里可以设置不同的采样时间,这样做之后仿真速度提升了3-5倍,波形精度损失很小。
另外就是避免在循环里面加连续积分器,用离散积分器替代,数值稳定性也更好。如果模型确实太大,还可以考虑把充放电过程拆分成两个独立的模型来跑:一个是充电过程模型,一个是放电过程模型,中间用初始状态继承来衔接。这样每个模型的计算量小,调试更方便,最后再把两个模型合并验证。
4.6 快速排查表
我把做这个仿真过程中最容易踩的坑整理成一个速查表,按优先级排序排查:
| 现象 | 可能原因 | 排查手段 |
|---|---|---|
| 仿真初瞬即发散 | 坐标变换符号错误、电流反馈极性反 | 检查Park变换角度符号,断开转速环只测电流环 |
| 电流振荡 | 电流环PI增益过高、解耦项符号错误 | 降低Kp_i,关闭解耦项验证 |
| 转速超调大 | 转速环积分饱和、限幅过低 | 给PI加抗饱和,提高限幅到1.5倍额定电流 |
| 母线电压跌落 | 模式切换无初始值匹配 | 切换前记录iq,预置电压环PI初值 |
| 高速时转矩不足 | 电压饱和、未加弱磁 | 检查SVPWM线性范围,启用弱磁控制 |
| 仿真极慢 | 连续求解器步长过小 | 改用离散模式,合理设置采样时间 |
5. 充电阶段深度解析与功率转矩特性
5.1 充电过程的三段式控制策略
把充电过程进一步拆开看,完整的充电控制其实是三个阶段,每个阶段对控制器有不同的要求:
第一阶段是最大转矩加速阶段。转速从静止到某一预设定转速的过程中,转速环一直处于饱和状态,iq指令稳定在限幅值,电机输出恒定最大转矩,转速线性上升。这个阶段电流环的跟踪性能和直流母线的供电能力是主要约束。如果直流母线电压不够高,大电流会造成母线电压跌落,反而限制了最大转矩输出。
第二阶段是恒功率加速阶段。随着转速升高,反电动势增大,电压余量逐渐减少,电机从恒转矩区过渡到恒功率区。这时候q轴电流会自然回落,因为同样的电压极限下,高转速能输出的最大转矩变小了。我在模型里观察到的现象是:转速达到约8000rpm后,iq指令从200A开始回落,转矩功率乘积基本保持恒定。
第三阶段是转速维持阶段。达到目标转速后,转速环从饱和退出,iq减小到维持摩擦力矩和风阻损耗的数值。对于高速飞轮系统,这个维持电流并不为零,因为高速下的风阻损耗虽然真空环境可以降低,但不完全为零,轴承摩擦损耗同样存在。
充电策略里还有一个放电深度管理的问题需要考虑,也就是不能让飞轮一直充到过高的转速,要给最高转速留一个保护余量,通常控制在额定转速的105%以内。
5.2 充电功率特性与直流侧需求
充电过程的功率特性对前端变换器选型有直接影响。最大转矩加速阶段,机械功率Pm = T * ω,转速上升的同时功率上升,直流侧输入功率从零平滑爬升到最大值。这个特性对前级AC/DC变换器比较友好,因为不存在阶跃式的功率冲击。
我算过一组数据验证功率匹配关系:最大转矩140Nm,在8000rpm时机械功率Pm = 140 * (2π * 8000/60) ≈ 117kW。考虑电机和逆变器损耗后,直流侧需要输入的功率大约在120kW左右。这就是为什么这台55kW额定功率的电机,在峰值转矩状态下短时输入功率会明显超过额定值,这是飞轮快速充能的特性需求,也是电机和功率器件需要短时过载能力的原因。
模型里把直流母线的电容和电压源内阻也做了建模,目的是验证大电流充电时母线电压的波动幅度。实际仿真下来,700V母线在峰值充电电流条件下电压波动控制在3%以内,如果有超出,需要增加母线电容或前级响应速度。
5.3 充电结束的转速保持策略
飞轮储能系统充满电之后,不是立刻放电,而是要维持在额定转速待命。这个阶段有个参数值得关注:转速保持的精度和待机损耗。
转速保持仍然用转速环,但这时候转速指令保持不变,干扰来自风阻和轴承摩擦的缓慢变化。转速环的积分作用可以把稳态误差消除到接近零,但代价是持续的电流输出。我做了一个估算:待机时所需的维持转矩大概是额定转矩的1%-2%,对应的iq维持电流约2-4A,直流侧输入功率约1-2kW。这个待机损耗决定了飞轮储能的"自放电"时间常数。
为了降低待机损耗,工程上会在真空腔体内维持高真空度,把风阻降到忽略不计。仿真模型里我给了一个风速相关的阻尼系数,并且在模型里留了接口可以在不同真空等级下切换。实际调试中发现,如果真空度足够好(气压低于1Pa),风阻转矩可以忽略,待机时只需要补偿轴承损耗,这样储存的能量可以保持较长时间,这也是飞轮储能系统适合做短时高频次储能的根本原因。
6. 放电阶段控制精讲与母线支撑逻辑
6.1 母线电压闭环的实现逻辑
放电时飞轮储能系统的角色从电压接受方变成电压支撑方。在独立直流母线系统中,飞轮系统通过控制逆变器输出功率来维持母线电压稳定。
母线电压环的传递链是:Udc_ref与Udc_feedback的误差经过PI控制器输出iq_ref,电流环跟踪iq,电机发电产生负转矩,转速下降,机械能转化为电能输入直流母线。这里有一个因果关系值得注意:母线电压环控制的本质是功率平衡,并不是直接控制电压。当外部负载功率增大时,母线电压会先下降一个小的偏差,电压环感知到偏差后增大电流输出,补充功率缺口,电压恢复。所以这个闭环的物理延时取决于母线电容的大小和电压环带宽。
放电时电流环的dq轴控制和充电时是一样的,只是q轴电流参考值是负的。为保护电机和逆变器,放电电流同样需要限幅。因为发电工况下,如果突然甩负荷,机械功率无法瞬间降为零,会产生过电压。所以模型中给放电工况做了直流母线过压保护,当母线电压超过740V时立即封锁PWM,并切到能耗制动模式,防止飞轮系统失控。
6.2 放电深度与转速下限保护
放电过程飞轮转速从12000rpm下降,存储动能E = 0.5 * J * (ω²_max - ω²_min)。这个公式说明放电深度对释放能量有平方关系影响,转速从12000rpm放到6000rpm可以释放75%的动能,而放到9000rpm只能释放43.75%的动能。
但转速下限不能为零,原因有三条:
第一,永磁同步电机反电动势与转速成正比,转速过低时反电动势太小,无法维持母线电压,发电能力丧失。第二,转子永磁体在高速旋转下如果突然降到很低转速,轴承润滑油膜和动压轴承的支承可能失效,机械安全性存疑。第三,飞轮系统的控制策略通常在深度放电末端需要切换回电机模式来维持飞轮基础转速,避免停转后需要重新从零加速。
模型里我把转速下限设为6000rpm,对应最大放电深度50%。这个值对应放电结束后的母线电压维持能力足够(反电动势还有额定值的一半),而且再次充电的加速时间相对较短。实际工程中,这个下限要根据具体系统和工况来设定,有的飞轮储能系统为了追求更大的放电能量,会把下限压到额定转速的40%。
6.3 放电工况的功率平滑策略
模拟真实负载工况时,负载往往是波动的,母线电压闭环需要应对这种波动。为了验证模型的动态响应能力,我设置了一个可编程负载序列,模拟功率在20kW到100kW之间阶跃变化。
实测母线电压的恢复时间大约在10-20ms之间,超调量小于5%。这个动态响应指标对于飞轮储能系统来说是可以接受的。如果把母线电压环带宽再调高,理论上可以获得更快的恢复速度,但过高的带宽会和电流环带宽接近,引发控制环路之间的耦合振荡。
在功率平滑策略上,我建议使用一阶低通滤波对负载功率参考值做平滑处理,限制功率变化率。这样做的好处是可以减小储能单元的冲击电流,延长功率器件寿命,代价是动态响应慢一些。对于电网调频这类需要快速响应的应用场景,平滑时间常数要取得很小,或者干脆不用滤波,直接让电压环快调。
7. 系统性能评估与待机损耗分析
7.1 充放电效率的仿真评估方法
储能系统的效率是最核心的性能指标。仿真中评估效率的方法是记录整个充放电过程中的能量值变化,然后对比放电输出能量和充电输入能量的比值。
我在模型中加了若干个积分器来记录能量。充电过程中直流侧的输入能量通过积分Udc * Idc计算,放电过程中负载消耗的能量通过积分负载功率得到。仿真结束后用能量比值就能算出系统效率。
以我这组参数为例,充放电深度50%的循环效率仿真结果大约是88%-90%。这个效率的量级比较合理,铜损、铁损、开关损耗、摩擦损耗都算进去了。如果把中间待机时间拉长,待机损耗会让循环效率明显下降,这也解释了为什么飞轮储能适合高频次短时储能而不是长时间储能。
7.2 待机损耗的来源与参数敏感性
待机损耗主要有三块:风阻损耗、轴承摩擦损耗、电机铁损和控制系统的自耗电。仿真模型中,风阻损耗用转速相关的阻尼系数模拟,轴承摩擦损耗用恒定摩擦转矩模拟,铁损通过电机模块内部参数计算。
我特意做了参数敏感性分析,看待机损耗对哪些参数最敏感。结论是风阻损耗对转速最敏感,转速降10%,风阻损耗下降约30%。所以如果应用场景允许,降低待机转速比维持满转速待机要划算得多,哪怕需要快速升速时再加速回来。这也是有些飞轮系统采用"低转速待机、高速充能"策略的原因。
7.3 模型的扩展方向
这套模型搭建完成之后,可以往几个方向做扩展。
第一个方向是接入上层调度。把充放电控制封装成一个可接受功率指令的储能单元模型,然后放到微电网系统仿真里,验证飞轮储能在削峰填谷、调频辅助服务等场景中的作用。
第二个方向是并联扩容。实际工程中单台飞轮机组功率不够时,会把多台飞轮并联接入同一个直流母线。并联系统的控制难点在于多台飞轮的荷电状态不一致,需要做功率分配策略。基于这套模型,可以复制出多台机组,在功率分配算法上做研究。
第三个方向是做硬件在环。用实时仿真机运行电机和主电路模型,控制器跑在实际的控制板上,验证实际控制器代码的正确性。硬件在环的重点就是模型离散化程度和I/O配置,基于这套模型做离散化改造是可行的。
7.4 模型验证与数据校准的实操心得
仿真模型做得再精细,终究是简化模型,要保证仿真结果对工程有参考价值,必须做验证。我的做法是分三步验证:
第一步是空载验证。把飞轮负载去掉,只保留电机空载运行,对比仿真反电动势波形和电机参数表中给出的数值是否一致。这一步如果对不上,先查坐标变换和参数换算。
第二步是阶跃响应验证。给一个小的转速阶跃,记录iq响应波形和转速响应波形,检查超调量、上升时间是否和理论计算一致。这一步验证PI参数整定是否正确。
第三步是能量闭环验证。做完整的充放电循环,检查初始转速和最终转速是否一致,母线电压是否维持稳定,能量积分结果是否符合效率预期。这一步验证系统级的能量管理是否正确。
三次验证都通过之后,仿真模型才算是可信度较高的工具,后面调控制策略、做工况分析才能有依据。盲目相信仿真结果的后果,我在做硬件在环时遇到过,控制器参数在仿真里表现很好,上到实际系统就振荡,原因就是仿真模型忽略了数字控制的离散延迟。所以做仿真的时候,要刻意把离散采样、PWM死区这些非理想因素加进去。
最后再说一个容易被忽略的细节:飞轮储能系统的仿真最好能加入低转速下的数值稳定处理。有些求解器在转速接近零时,电机反电动势接近零,状态变量可能出现数值奇异。解决的办法是给机械方程式加一个很小的阻尼项,或者在仿真初始化时给出一个初始转速值,比如100rpm,避免从完全静止开始仿真时的数值抖动。这个小技巧能让仿真的启动过程顺畅很多,尤其对后续要扩展到硬件在环的朋友来说,这个处理几乎是必须的。