孤岛微电网的二次控制在仿真里做过的人都有感受下垂控制一投进去负荷波动后频率和电压必然偏离额定值波形是稳了但50Hz和380V就是回不去。这时候就需要二次控制出手。但这几年大家做的二次控制大多是周期采样、周期通信每个控制周期都广播一次数据。通信网络空闲的时候无所谓如果DG数量一多、通信带宽受限或者你在做多微网互联、虚拟电厂这类大场景周期通信的成本和拥塞问题就会变得非常刺眼。事件触发机制就是在这种背景下被反复提到的——它不是按时间周期发送数据而是状态偏差超过某个阈值才发一次目的是大幅降低通信频次同时把电压和频率协同拉回额定值。这篇文章用一个Simulink仿真模型来聊整个实现链路从孤岛微电网的一次下垂控制为什么需要二次修正到事件触发条件怎么设计、触发逻辑在Simulink里怎么落地再到协同控制律的参数怎么调、波形怎么看、常见坑有哪些。适合研究生做微电网方向课题、搞DG并网控制的工程师以及正在做Simulink仿真但被二次控制参数和通信逻辑折磨的同学。下面这些内容和参数都是我实际调模型、跑波形时逐步定下来的可以直接参考也可以在这个基础上改。1. 孤岛微网仿真里二次电压/频率控制为什么值得单独写一篇很多人刚开始搭微电网仿真时第一版模型往往就是DG加下垂控制。确实下垂控制能实现多台DG之间的功率分配也不需要DG之间通信结构简单仿真跑起来稳定又直观。但注意它的代价是电压和频率的稳态偏差。1.1 下垂控制带来的偏差是结构性的不是参数没调好以两台DG并联带一个30kW负荷为例。采用P-f下垂时有功功率增大频率会线性下降采用Q-V下垂时无功功率增大电压幅值会线性下降。这是下垂特性曲线的自然结果。你把下垂系数调小偏差会小一些但功率分配的精度又会变差因为下垂系数本质上就是各DG之间的功率分配比例。频率偏差和功率分配精度是一对矛盾光靠一次控制无法同时满足。所以在孤岛微电网中频率稳态偏差可能达到0.3~0.5Hz电压偏差可能达到1%~3%。对电能质量要求不高的场景还能凑合但如果后面还带着精密负载或者并网转离网切换这种偏差就不能接受。二次控制的作用就是通过平移下垂特性曲线把这些稳态偏差重新修正回额定值对应的控制量是一个补偿项叠加到一次控制的参考值上。1.2 二次控制为什么必须协同而不能每台DG各干各的有同学会问那每台DG本地放一个PI控制器把本地的频率偏差和电压偏差积分回零不就行了吗在理想情况下确实可以但这有两个问题。第一电压幅值在微电网里不是一个全局一致的物理量。因为线路阻抗的存在不同位置的DG测量到的电压幅值本身就是不一样的。如果每台DG都试图把本地电压拉到一个完全相同的额定值而无功功率的流动又会影响电压那么各DG的电压控制器之间会产生明显的耦合振荡。第二频率虽然在全网是基本一致的但频率测量本身有噪声和相位差各自为战会导致稳态后各DG的补偿量不一致进而让有功功率分配结果偏离下垂设计值。所以现在的主流做法是分布式协同控制每台DG只需要和邻居DG交换信息通过一致项把大家的状态拉到同一个目标值同时用本地参考项把该目标值拉回额定值。这就是电压与频率协同控制里协同二字的含义。1.3 仿真验证的典型拓扑与基础参数我用的模型是4台DG组成一个辐射状孤岛微电网每台DG由一个直流电压源、VSC换流器、LC滤波器、耦合电抗和本地负荷构成DG之间通过线路阻抗连接。额定频率50Hz额定线电压380V直流侧电压750V。仿真里比较省时间的做法是把VSC先用平均模型替代也就是不搭IGBT和PWM高频开关直接由控制环输出调制电压。这样可以用较大的固定步长跑先把二次控制逻辑和通信逻辑调通再换回开关模型看谐波和纹波。两种模型的二次控制部分完全一致不需要改动。2. 事件触发机制到底解决什么问题从周期通讯到按需通讯周期通信的缺点是显而易见的。每个仿真步长或者说每个控制周期每台DG都要把本地状态打包发送到通信网络广播给邻居。4台DG的时候看不出来但如果节点数是10台、20台网络报文数量是随连接数增长的。事件触发机制的核心思路是在满足控制性能的前提下降低通信次数。它是一种按需发送的调度策略只有当某台DG的当前状态与它上一次发送给邻居的状态之间的偏差超过某个阈值时才触发一次新的通信。2.1 事件触发条件的具体形式最常用的事件触发条件我取的是这种形式t_{k1} min{ t t_k | | x(t) - x(t_k) | sigma * | x(t_k) | epsilon }其中x(t_k)是最近一次触发时刻发送出去的状态值x(t)是当前时刻的实际状态二者的差称为事件触发误差。sigma是相对阈值epsilon是绝对阈值用来防止状态接近零时因为相对阈值过小而频繁触发。这里用一个相对阈值加一个绝对阈值而不是只用一个绝对阈值是有原因的。比如频率偏差在恢复过程中是逐渐变小的如果只设绝对阈值0.1Hz那么当偏差从0.5Hz减小到接近0.1Hz的过程中触发器会连续触发好几次并没有比周期通信节省多少。而如果加上相对阈值当状态本身已经很小时触发条件自动放宽通信频次自然降下来。2.2 两次触发之间控制量如何保持这里有个关键点事件触发机制下控制器并不是每个步长都收到新鲜数据。两次事件之间邻居DG接收到的状态是之前触发时刻的旧值控制量需要保持这个旧值不变直到下一次触发发生。这在控制理论里叫零阶保持器ZOH特性在Simulink里实现起来就是一个采样保持模块加上一个触发脉冲控制。这个保持特性带来的好处是通信资源被省下来了坏处是事件触发控制下的系统变成一个混合系统连续动力学过程加上离散跳变过程。稳定性分析比连续系统复杂这也是为什么很多论文在事件触发后面都要挂一个李雅普诺夫函数证明。但仿真层面我们更关注实际行为只要触发间隔合理、参数不激进模型一般都能稳定收敛。2.3 通信量的量化指标为了评估事件触发到底省了多少通信我在每个DG的输出路径上加了一个计数器模块用触发脉冲的上升沿累计触发次数。在同样的仿真时长内周期通信的方案总通信次数等于仿真总时长除以控制周期事件触发方案的通信次数就是计数器读出的数字。我实测的一个典型工况仿真时间5秒控制周期10ms周期通信总次数是4台DG各500次共2000次报文。换成事件触发后相对阈值取2%时4台DG的总触发次数大概在300到500次之间通信量下降了70%到80%。而且在这个通信量下频率恢复时间和电压收敛波形与周期通信相比差距很小只有平方级别的微小延缓。3. Simulink整体模型怎么搭四层架构、模块清单与信号走向建Simulink仿真模型时最怕东一块西一块控制逻辑和主电路信号混在一起改一个参数满图红线。我把模型按功能分成四层每一层在子系统内部管理层与层之间只有明确的信号接口。3.1 四层架构总览模型的层级从上到下依次是主电路层、一次控制层、二次协同控制层、事件触发与通信层。主电路层放DG的电源、换流器与滤波电路输出三相电压电流信号同时测量各DG出口的有功功率P、无功功率Q、频率f和电压幅值V。一次控制层做的是本地下垂控制和电压电流双闭环控制它接收二次层的补偿量输出调制电压。二次协同控制层接收本地状态和邻居状态计算频率和电压的二次补偿量。事件触发与通信层负责判断是否需要发数据以及维护每个DG发送给邻居的触发状态值。接口设计上我让主电路层直接输出一个包含所有测量值的总线信号一次控制层和二次控制层分别从这个总线取自己需要的信号。二次控制层的输出是一个二维补偿量数组分别接入对应DG的一次控制参考值上。通信层输出的不是状态本身而是该DG最新触发时发出的状态值这一点要特别注意否则模型里很容易把实时值当成通信值用。3.2 模块清单与关键参数以4台DG为例我用的核心模块和参数如下模块说明关键参数DC电压源模拟直流侧750VVSC换流器PWM或平均模型开关频率按实际选型平均模型无需设置LC滤波器滤波L1.5mHC50uF线路阻抗DG间连线R0.1Ω/kmL1mH/km实际按距离折算本地负荷电阻性/感抗性30kW 5kVar下垂控制器P-f/Q-V下垂kp0.0002 rad/s/Wkq近似按额定功率折算电压电流双环内环电流外环电压带宽约为开关频率1/10二次PI频率恢复/电压恢复频率环kp0.8ki15电压环kp0.5ki8事件触发器MATLAB Functionsigma2%epsilon0.02采样保持模块零阶保持按触发脉冲驱动这里最容易被忽视的是下垂控制器的输出限幅。二次补偿量叠加到参考值上之后如果一次层没有加饱和模块模型在大扰动时可能出现频率或电压超调甚至发散。Simulink里的饱和模块是廉价的保险丝加上它既能避免积分饱和问题也更容易调参。3.3 信号在层间是怎么流动的我画信号流向时遵循一条规则二次控制层只和一次控制层的参考值输入打交道不直接修改主电路参数。一次控制层只向下输出调制电压和PWM信号不向上反馈任何东西。事件触发通信层挂在二次层一侧不参与主电路运算。具体来说每台DG的测量模块计算出本地频率w_i和电压幅值V_i后送入事件触发器。触发器输出经过通信通道仿真里直接用Simulink信号线模拟也可以加一个传输延迟模块模拟网络时延送到邻居侧。二次协同控制器收到本地实时值和邻居触发值后计算出补偿量。最后在DG的下垂方程中加入补偿量w_ref_i w_nom m_p_i * P_i Δw_i V_ref_i V_nom n_q_i * Q_i ΔV_i其中Δw_i和ΔV_i就是二次协同控制的输出。这样设计的好处是模型层次清晰后续如果要升级到一致性控制器或者加上通信拓扑切换只需要改二次层和通信层一次层和主电路层完全不需要动。4. 协同控制律的写法本地偏差矫正加邻居同步项的赛跑逻辑这一节单独拎出来说因为协同控制律是整个模型的控制核心也是论文里公式最多、仿真里最容易出问题的地方。我用的控制律是经典的分布式二次控制结构包含两个部分本地偏差矫正项和邻居一致性同步项。4.1 频率协同控制律频率部分的控制律写成Δw_dot_i k_f * (w_nom - w_i) k_c * Σ_{j∈N_i} (w_i - w_j)这里Δw_dot_i是二次补偿量的变化率对它积分就得到Δw_i。第一部分k_f * (w_nom - w_i)的作用是把本机频率拉向额定值第二部分k_c * Σ(w_i - w_j)是邻居一致性项作用是让所有DG的频率最终达到同一个值。如果只有第一部分那么各DG会因为频率测量噪声和控制器响应差异各自收敛到额定值附近的微小不同点。如果只有第二部分各DG频率会收敛到一致但不一定等于额定值。两部分合在一起同时实现一致性和额定值跟踪这就是协同的内涵。4.2 电压协同控制律的区别电压部分形式上类似ΔV_dot_i k_v * (V_nom - V_i) k_c * Σ_{j∈N_i} (V_i - V_j)但要注意电压幅值本身不是全局一致的物理量。由于线路阻抗和负荷分布DG1的出口电压和DG4的出口电压天然就有差异。所以这里的V_nom - V_i项如果用的偏差过大的话不同DG会产生方向不同的电压调整动作无功功率可能会出现来回涌动。处理办法有两个。一种做法是在电压偏差项里加入一个额定值范围判断偏差在一定范围比如±2%内时弱化额定值项的权重主要依赖一致性项让各DG电压趋于各自合理的值。另一种做法是改用电网平均电压作为参考值但这需要额外的平均观测器模型复杂度上来了。我仿真时用的是第二种思路的简化版在电压环的额定期望项前面加了一个可调的权重系数收敛过程让一致项先起作用等电压接近一致后再让额定项起作用。这个时序通过一个很小的权重系数就实现了模型里就是个增益模块而已。4.3 通信拓扑的影响协同控制律的效果强依赖于通信拓扑的连通性。如果通信网络是链状信息只能相邻传递那么DG4的扰动传到DG1需要经过多跳收敛速度比全连接拓扑慢。我仿真了两种拓扑全连接和环形。全连接收敛明显快但通信次数也更多。环形收敛慢一些但触发次数更少。如果做课题想突出事件触发的价值建议用环形或稀疏拓扑做对比这样通信减少的百分比更明显。如果只是想展示协同控制的效果全连接拓扑更省心因为邻居数量多一致性项的作用更充分参数容错率更高。5. 事件触发器里那段MATLAB Function源码与Zeno防护这是整个模型里代码量最小但最关键的模块。我用一个MATLAB Function实现事件触发判断、状态保持和最短触发间隔限制。Simulink里其他模块也能搭出同样的逻辑比如用比较器和RS触发器但MATLAB Function更直观改阈值也方便。5.1 源码与模块接线每个DG配一个独立的事件触发器子系统。输入有三个当前测量值x_cur、上一次触发后保持的发送值x_hold、使能信号enable。输出有两个x_send经过采样保持的待发送状态值和trigger_pulse一个布尔脉冲上升沿表示触发。function [x_send, trigger_pulse] event_trigger(x_cur, x_hold, enable) sigma 0.02; % 相对阈值 epsilon 0.02; % 绝对阈值 min_interval 0.05; % 最小触发间隔秒 trigger_pulse false; x_send x_hold; if enable % 触发条件判断 if abs(x_cur - x_hold) sigma * abs(x_hold) epsilon trigger_pulse true; x_send x_cur; end else x_send x_cur; end end这只是单纯触发判断但实际模型中必须处理最小触发间隔的问题。在事件触发控制领域有个经典问题叫Zeno现象指的是触发时间点无限逼近导致在有限时间内产生无限多次触发。在Simulink仿真里表现为某个瞬间触发脉冲密集到几乎每个仿真步长都触发计数器疯狂跳变系统退化为周期通信甚至因为步长限制导致仿真卡死。Zeno的一个常见来源是测量噪声。当状态值在阈值边界附近抖动时abs(x_cur - x_hold)反复越过阈值触发器会连续触发。解决方法是加入强制的最小触发间隔也就是在一个触发发生后至少间隔固定时间才允许下一次触发。这就是上面代码里min_interval的作用——不过上面的代码还只是写了参数具体实现需要把它接进循环我实际用的是下面这种带循环逻辑的写法。function [x_send, trigger_pulse] event_trigger(x_cur, x_hold, t) sigma 0.02; epsilon 0.02; min_interval 0.05; trigger_pulse false; x_send x_hold; persistent last_trigger_time; if isempty(last_trigger_time) last_trigger_time -inf; end if t - last_trigger_time min_interval if abs(x_cur - x_hold) sigma * abs(x_hold) epsilon trigger_pulse true; x_send x_cur; last_trigger_time t; end else x_send x_hold; end endt是仿真时间通过把仿真时钟接入MATLAB Function来获取。last_trigger_time用persistent保存防止每次调用都被重置。这样即使测量值在阈值附近抖动相邻两次触发之间也至少间隔0.05秒Zeno被彻底挡在门外。5.2 采样保持回路怎么接x_hold的来源决定了整个回路能不能正常闭环。正确接法是把x_send通过一个Memory模块或Unit Delay模块反馈回触发器输入作为x_hold。这里的逻辑是只有当本次触发发生时x_send才更新为新值更新后作为新的x_hold等待下一次偏差超限。如果不加这个反馈回路x_hold始终等于初始值那么当状态值收敛到额定值附近时abs(x_cur - x_hold)会一直很大触发器会持续误触发通信量反而比周期通信还高。接法可以这样x_send输出后并联两路一路去通信网络另一路经过一个Memory模块回到触发器输入端的x_hold。Memory模块用于离散时刻的采样保持和Unit Delay的区别是它不引入额外的离散采样周期适合放在这种由事件驱动的回路里。5.3 触发次数统计与可视化为了直观看到通信节省效果我会在每个事件触发器的输出端再接一个计数器子系统。用trigger_pulse作为使能信号的上升沿每触发一次计数加1最后用一个Display模块读数值。还有一个更直观的方法把trigger_pulse信号接到Scope里纵轴显示为布尔量就能直接看到触发时刻的稀疏程度。频率恢复初期触发脉冲密集恢复后期脉冲几乎消失这个波形比任何统计数字都更有说服力。我在调整sigma参数的时候就是看着这个脉冲波形来定的如果脉冲在整个仿真时长内一直密集出现说明阈值太小如果脉冲在开始阶段就非常稀少说明阈值太大可能影响收敛性能。6. 调参复盘从振荡、慢收敛到波形正常的三个关键阶段二次协同控制的调参我实际走过了三个阶段每个阶段都有典型的失败波形和对应的解决方法。直接给经验省得你们再踩一遍。6.1 第一阶段二次环和一次环带宽打架频率波形像锯齿我最初把二次控制器的PI参数调得和一次下垂环相当比例项kp2积分项ki50。结果波形出来频率刚开始恢复时振幅先冲到50.3Hz然后拉回接着又冲过来回振荡衰减看起来就像锯齿波。原因在于二次控制环和一次下垂环的带宽没有拉开。一次下垂环的响应时间常数大概在几十毫秒的量级如果二次PI的积分时间也在这个量级就会形成一个双环抢带宽的局面系统在两个环的交互处产生共振。解决方法很直接把二次环的速度明显放慢。频率环kp0.8ki15积分时间大概在1.3秒左右相当于一次环的10倍以上。这样二次环的补偿作用在一个相对慢的时间尺度上进行一次环已经完成了功率分配的快速动态二次环再去缓慢修正稳态偏差两者各管一段波形就平滑了。6.2 第二阶段邻居一致性系数过大DG之间开始互搏刚开始调协同控制的时候为了让各DG频率快速拉齐我把一致性系数k_c设成了比例项的一半也就是0.4。结果4台DG的频率确实很快拉齐了但出现了另一个现象DG1和DG2的频率在做幅度很小的反相振荡DG3和DG4也是振荡频率大概在2Hz左右看不到明显衰减。这是典型的耦合振荡。一致性项的工作原理是如果我和邻居不一样我就往邻居的方向调整但如果比例过强调整量会超过实际偏差需要的量然后邻居又反过来调整形成你拉我扯的往复振荡。在Simulink里表现出来就是相邻DG的状态量成对振荡。解决方法是把一致性系数调低我最终取的是k_c 0.1大约是本地比例项的1/8。这个值下收敛速度比周期通信慢0.2秒左右但完全消除了耦合振荡。调整的时候每次只动20%观察波形至少5秒不要急着下结论。6.3 第三阶段事件触发的阈值权衡省通信和保性能的平衡点事件触发阈值也不是随便设的。sigma设太小比如0.5%触发次数接近周期通信失去事件触发的意义设太大比如10%通信倒是省得极多但二次控制的收敛时间明显变长因为邻居拿到的状态值长期是旧值一致性协调的时效性变差了。我实际测试下来一个工程上合理的范围是2%~5%。在这个范围内通信量可以降低60%~80%同时频率恢复时间与周期通信方案相比只增加5%~15%。如果对性能要求更高可以取1%如果通信带宽极度受限可以放宽到5%并且配合增大epsilon来滤除小信号抖动。还有一个和阈值配套的变量最小触发间隔。如果设得过大比如0.5秒那么在这个间隔内的任何大幅扰动都无法触发通信二次控制会变得迟钝。我建议最小触发间隔设为控制周期的5到10倍即可10ms的控制周期下取0.05秒到0.1秒是比较平衡的区间。6.4 一个容易忽略的细节触发的是状态还是补偿量事件触发通信中应该发送什么量这个选择很影响仿真的收敛效果。我在第一版模型中触发的是二次补偿量Δw_i和ΔV_i本身邻居收到后直接用于一致性计算。后来发现一个问题当负载突然增加时一次下垂先快速响应二次补偿量还在缓慢增加此时邻居之间交换的补偿量偏差很小触发器可能长期不触发但实际频率已经偏离额定值了。更合理的做法是触发本地频率和电压测量值二次控制律在本地计算后输出补偿量。这样触发条件和控制目标直接挂钩频率一旦偏差超限马上触发通信二次补偿量才能及时跟上。这个改动不算大但效果差异非常明显。我在最终模型中通信层发送的始终是w_i和V_i而不是Δw_i和ΔV_i。7. 离线仿真中一定会踩的坑与扩展思路模型搭完后还有一堆仿真和调试层面的坑。这里挑几个最常见的记录下来每一个我都实际碰到过。7.1 代数环问题事件触发回路里用了Memory后一般能避开代数环。但如果在触发器输出到二次控制器再到一次控制器再到主电路测量再到触发器的这条链路中某个模块忘了加Memory或Unit DelaySimulink会在仿真开始时报代数环错误或者自动插入隐式求解器而导致仿真速度骤降。排查方法也很简单仿真结束时在Diagnostics Viewer里看是否有警告提示有代数环或者干脆在信号线上右键查看哪些信号参与了代数环检测。最快的解决方式是在反馈回路任意位置加一个Memory模块。但注意Memory模块会引入一个步长的延迟如果加在控制主回路上可能会让相位裕度下降所以优先加在通信状态保持的那条支路上。7.2 变步长仿真下的事件漏检Simulink的变步长求解器会根据误差自动调整步长步长最大可能达到毫秒甚至秒级。如果事件触发器判断使用的是连续时间信号当步长过大时一个触发点可能恰好落在两次采样之间触发器就会漏掉这次事件。应对方法有两种。一是把事件触发相关的状态更新改成离散模块显式指定采样时间比如10ms这样触发判断和通信频次的统计就是确定的。二是使用固定的较小步长比如1ms或10ms跑完整个仿真代价是仿真时间变长但结果可复现性更好。我做参数扫描时都用固定步长只有最后出结果图时才换回变步长试一下。7.3 数据类型的隐式转换MATLAB Function里输出的trigger_pulse默认是double类型直接接到Enable端口或者计数的上升沿检测模块时可能会因为类型问题报警。在MATLAB Function代码里声明输出变量的类型或者在模块输出后加一个Data Type Conversion模块转成boolean能省掉很多莫名其妙的问题。7.4 从离线仿真走向硬件在环和扩展场景这套模型可以往两个方向扩展。第一个方向是更复杂的微电网场景并网转孤岛、孤岛转并网的切换或者把DG数量扩展到8台并加入非均匀线路参数。事件触发机制在这种更复杂的暂态过程中的表现比简单孤岛工况更有研究价值。第二个方向是硬件在环验证把Simulink模型编译部署到实时仿真平台上和真实的通信网络设备对接。这时候事件触发机制的价值会体现得更充分因为真实通信网络里的带宽限制、延迟和丢包都是离线仿真里无法完全模拟的。另外提一句有些人会把这种微电网仿真和Carsim与Simulink联合仿真相混淆。Carsim联合仿真的是整车动力学模型微电网控制里没有车辆动力学对象两者用的Simulink工具箱和模型架构完全不同。微电网仿真一般用SimPowerSystems也就是现在的Simscape Electrical加上基础的Simulink控制模块。开工具箱之前先看清楚对象别下载一堆用不上的模块库。回到模型本身我实际调通之后最深的体会是事件触发协同控制的难点不在理论公式而在Simulink里如何把触发逻辑、采样保持、通信延时和PI参数串成一个闭环系统。顺序上先调通周期通信的二次控制再替换成事件触发通信这样即使事件触发环节出了问题也能立刻判断是控制参数问题还是触发逻辑问题。这个习惯帮我省了大量排查时间建议你搭模型时也按这个顺序来。