最近在梳理孤岛微电网的二次控制方案手头这个基于事件触发机制的Simulink仿真模型值得好好聊一聊。孤岛微电网在脱离主网后电压幅值和频率很容易偏移传统的一次控制只能靠下垂特性硬扛想要做到无差调节就必须上二次控制。但二次控制如果按周期广播数据通信开销和计算压力都不小。我这次选择事件触发机制让控制器只在系统偏差超过触发阈值时才交换邻居信息既保证了恢复精度又把通信次数降了一个数量级。模型里同时兼顾了二次电压恢复和频率恢复二者通过一套协同框架耦合在一起在Simulink里完整跑通了。这篇文章会把主电路拓扑、事件触发逻辑、协同控制器设计、模型搭建步骤和踩坑记录都放出来给准备做微电网二次控制仿真的朋友一个参考。1. 项目背景与核心问题拆解1.1 为什么孤岛微电网必须做二次控制孤岛微电网相当于一个独立的小型电力系统由分布式电源、储能、负载构成。并网状态下电网电压和频率由大电网支撑本地逆变器只要跟住电网就行。一旦静态开关断开进入孤岛系统就失去“靠山”所有电压频率都由本地DER分布式能源共同维持。最常用的维持方式是下垂控制也就是模拟同步发电机的一次调频调压特性。有功功率和频率挂钩无功功率和电压幅值挂钩。下垂控制不需要通信动态响应快但它的本质是有差调节负载一变频率和电压就偏离额定值。比如一条典型的P-F下垂曲线空载频率50.2Hz满载频率49.8Hz中间是线性下倾的。在50%负载下频率会稳定在50.0Hz附近似乎还行但负载变化剧烈时偏差可能超出允许范围。GB/T 33593-2017对微电网电能质量有明确要求频率偏差不能超过±0.5Hz电压偏差不能超过±7%。下垂控制很难一直满足这个指标所以必须有二次控制来“扶正”。二次控制的任务很直白检测到电压幅值或频率偏离额定值后通过调整参考量让下垂曲线整体平移从而把系统运行点拉回到额定工作点。传统的二次控制是集中式的由一个中央控制器采集所有DER的信息计算修正量后下发。集中式的优点是结构简单但中心节点一旦故障整个二次控制就瘫痪了而且通信时延和带宽压力都集中在中心扩展性差。分布式二次控制则让每个DER只跟邻居交换信息通过一致性算法逐步收敛到同一个修正量这是目前微电网二次控制的主流研究方向。1.2 事件触发机制解决了什么痛点分布式二次控制的前提就是通信。如果每个控制器都按照固定周期把自己的状态广播出去系统能稳定但通信资源被白白浪费。我在做仿真时统计过在一个4节点微电网模型里固定周期通信模式下每个仿真秒大概要交换2400个数据包实际工程中通信链路带宽有限而且无线通信的功耗、丢包、时延都会随数据量上升。事件触发机制的核心思想是“按需通信”。每个控制器本地维护一个触发条件只有当测量值和上次触发时刻的存储值之差超过设定阈值时才通过网络发送当前数据。可以类比成家里装了智能摄像头画面没变化时不推消息只有检测到运动物体才推送警报。这样通信次数大幅减少但系统状态依然能被感知到。对于微电网二次控制事件触发机制能在保证稳定性的前提下把通信资源占用降一个量级。把这个机制放到二次电压与频率协同控制里还有一个额外的挑战电压和频率是两个不同的物理量变化速率不同频率差异主要受有功功率影响电压幅值受无功功率影响。如果对二者使用同一个触发阈值和同一个触发条件容易出现“频率还没超差但电压已经需要更新”或者反过来。所以我的模型里把事件触发条件做成复合形式既考虑电压偏差也考虑频率偏差只要其中一个超过阈值就触发通信。这也是仿真里最值得关注的设计点。2. 系统总体架构与仿真模型设计2.1 微电网主电路拓扑与参数选型模型主电路采用典型孤岛微电网结构4台三相逆变器并联通过LC滤波器接入公共交流母线再统一向负载供电。额定线电压380V/50Hz直流侧电压用理想直流源800V代替这样可以忽略直流动态对控制的影响把注意力集中在二次控制策略上。逆变器开关频率设为10kHzLC滤波器参数按常规设计滤波电感L1.8mH滤波电容C50μF阻尼电阻Rd2Ω。这个参数组合下截止频率约在530Hz左右能有效滤除开关谐波同时保证工频信号顺利通过。4台DER的容量可以不完全一致我仿真时设置了两种配置等比容量和不等比容量。等比容量下4台都是50kVA不等比容量下分别是50kVA、40kVA、30kVA、30kVA。为什么不等比容量更能说明问题因为在真实微网里光伏、储能、柴油机容量天生不同二次控制不仅要恢复电压频率还要保证有功和无功功率按容量比例均分。如果容量不等二次控制的一致性算法设计不好就容易出现“大容量带得多、小容量带得少”的失衡情况。负载端我放了三个并联负荷额定有功30kW加无功10kvar的平衡负载加一组可在1.2秒投入的冲击负载20kW再加一组非线性负载来模拟谐波干扰。冲击负载用来测试频率和电压的动态恢复能力非线性负载用来考察控制策略对波形畸变的容忍度。负载全部用三相RL或带整流桥模块Simulink里直接用Three-Phase Series RLC Branch和Three-Phase Rectifier就能搭出来。2.2 控制层级划分一次控制与二次控制的边界在下垂控制框架下每个DER的控制器分为内环电流环、电压环和一次功率环。功率环计算输出有功和无功经过低通滤波器后生成电压幅值参考和频率参考频率参考ω_i ω_n - m_i * P_i δω_i电压幅值参考V_i V_n - n_i * Q_i δV_i这里的δω_i和δV_i就是二次控制修正量它们是分布式控制器计算出来叠加到Droop Controller的给定参考上的。内环一般采用电压电流双闭环电压外环用PI控制器电流内环也用PI控制器输出SPWM脉冲驱动逆变器。这些内环模型每个DER内部自己闭环不依赖通信属于一次控制范畴。二次控制只在更高层起作用。它采集各个DER的本地频率/电压信息通过一致性算法估算出全局平均偏差然后输出修正量叠加给一次控制。也就是说二次控制和一次控制之间的信号边界就是那两个修正量。我在Simulink里用Goto/From标签跨子系统传递信号避免一层层连线把模型搞乱。2.3 协同控制的数学模型与坐标变换要实现电压和频率协同不能各自为政。频率与有功、电压与无功之间有天然的数学联系所以协同控制的模型通常写成动态平均一致性形式。用λ表示某个待恢复的量每个DER的本地值记为λ_i通过邻居间交换λ_j按照以下动态更新λ_hat_i(t) λ_i(t) ∫₀ᵗ Σ_{j∈N_i} a_ij (λ_hat_j(s) - λ_hat_i(s)) ds这其实就是一致性算法里的平均一致项。所有节点的λ_hat_i最终收敛到全局平均值。电压和频率协同的巧妙之处在于控制器修正量可以写成δV_i k_pv * (V_ref - λ_v_hat_i) k_iv * ∫(V_ref - λ_v_hat_i)dtδω_i k_pf * (ω_ref - λ_f_hat_i) k_if * ∫(ω_ref - λ_f_hat_i)dt这里的λ_v_hat和λ_f_hat分别是通过通信估计出的电压和频率平均值。如果只做电压控制不管频率或者反过来都会让另一个量继续漂移。协同的核心是用同一个通信网络同时交换电压和频率信息在统一的触发条件下并行计算两个修正量。在Simulink里我把状态量封装成二维向量[V_i; ω_i]来传递整个一致性计算的实现就简化为对向量的加权求和。3. 事件触发机制原理与Simulink实现3.1 事件触发条件设计与阈值选择事件触发机制的关键是“什么时候更新通信”。我采用的触发条件为e_i(t) σ_i * (∥状态∥的基准量)其中e_i(t)是当前本地状态与最近一次触发时缓存的状态之间的差值。触发后发送端更新缓存值接收端更新邻居状态。阈值σ_i直接控制通信频率和精度的折中σ越小通信越频繁控制效果越接近连续通信σ越大通信越稀疏但性能会下降。在原仿真里我对电压和频率分别设置阈值。电压阈值为0.5V额定380V的0.13%频率阈值为0.01Hz额定50Hz的0.02%。你可能觉得电压阈值有点小其实电压本身波动比频率大且电压测量值传送的是经过park变换后的dq轴分量相当于进入了直流域用0.5V的阈值已经能筛掉大部分纹波干扰。频率信号我直接用PLL输出的数字频率比较干净所以阈值设得更小。为了验证阈值的影响我做了三组对比测试阈值放大5倍通信次数降到原来的30%但频率恢复时间从0.8秒拉长到1.5秒阈值缩小5倍通信次数增加约3倍恢复时间缩短到0.4秒曲线更平滑。实际工程中阈值应根据允许的动态偏差、通信带宽和传感器噪声来定没有绝对最优值只能做折中。3.2 触发信号生成模块的搭建在Simulink里实现事件触发核心是构造一个单位脉冲信号在触发时刻输出1其他时刻输出0。这个脉冲信号控制一个逻辑开关决定是否将当前测量值写入零阶保持器。我的做法是用Memory模块存储上一次发送的值注意Memory模块的初值需要设置为初始测量值否则第一个仿真步就会误触发。当前测量值减去Memory模块输出得到差值e。用Abs模块取绝对值再用Compare To Constant模块和阈值比较输出触发脉冲。触发脉冲接到Switch模块的控制端当触发脉冲为1时开关输出当前测量值否则输出前一个值并通过环回到Memory更新。这其实是一个数据流“自刷新”结构不需要S-Function也能实现。如果想更清晰也可以用Triggered Subsystem创建一个子系统触发端口类型改为rising edge内部直接用Input端口接收测量值输出到零阶保持器。我实测过两种方案功能等效但在大模型里Memory的代数环问题更少。3.3 通信链路与数据封装处理模型里4台DER之间的通信拓扑采用环形结构DER1和DER2通信DER2和DER3通信DER3和DER4通信DER4和DER1通信。环形拓扑是最简单的连通图能保证一致性算法收敛同时通信链路数量最少。在搭建时我把每条通信链路单独建模为一个“发送-接收”对用Goto/From标签模拟点对点传输。当然这只是一个理想通信模型没有考虑时延和丢包。如果想加通信延时可以在Goto/From链路里插入一个Variable Transport Delay模块延时值设为0.02-0.1秒之间。实测下来事件触发机制对0.05秒以内的时延不太敏感超过0.1秒会出现振荡。数据封装上我把每个DER发送给邻居的信息打包成一个1×4的向量[本机电压估计值, 本机频率估计值, 本机有功功率, 本机无功功率]。这个向量用Bus Creator打包后通过Goto标签传出接收端用Bus Selector解包。这里要特别注意Bus Selector的信号名称必须和Bus Creator里的信号名完全一致拼写错误或者层级不匹配都会报错“No signal available”。热词里提到“simulink bus selector 没有可选信号”就是这个原因。我在第6章还会专门说这个坑。4. 二次电压与频率协同控制器设计4.1 分布式二次控制器的拓扑结构协同控制器在模型里被封装为一个独立子系统输入是本地测量值和邻居接收值输出是δV和δω两个修正量。控制器内部结构分成三层通信层、一致性算法层、PI调节层。通信层负责处理事件触发后的数据交互一致性算法层通过加权积分计算全局平均估计值。具体到Simulink实现一致性算法层的核心是一个积分器输入是邻居估计值与本地估计值之差的加权和输出是本地估计值。这个积分器和测量值相加就是最终的平均估计值。权重矩阵A的元素a_ij按照通信拓扑设置连通且非负的权重矩阵能让算法稳定收敛。我的权重设置是相邻节点权重0.2自环权重0.6没有直接通信的节点为0。这个权重满足行随机矩阵的条件可以加快收敛速度。PI调节层就比较直白了两个PI控制器分别处理电压和频率。PI输出限幅要根据下垂控制的安全范围来设定。我做的是δV限制在±20Vδω限制在±2.5rad/s。限幅值太小会导致二次控制“使不上劲”太大则可能在启动阶段造成电压过冲。4.2 电压恢复控制环的实现电压恢复的目标是让公共母线电压平均值回到额定值。但每个DER测量的电压是本地端电压需要滤波处理。我在模型里用一阶低通滤波器带宽150Hz对电压幅值进行滤波避免触发逻辑被高频纹波反复触发。电压环的PI控制器参数比例系数kpv取0.05积分系数kiv取8。为什么比例系数这么小因为电压偏差经过一致性估计后本身已经被放大到平均值比例项太大容易在触发瞬间产生较大的修正量突变进而引起无功功率环振荡。积分系数则可以适当大一些保证稳态误差为零。实际调试时可以先断开事件触发通信把系统强制设为连续通信模式整定好PI参数后再切换到事件触发模式。这个顺序很重要很多新手上来就用事件触发调参数结果波形一塌糊涂还以为是通信机制的问题其实单纯是PI参数没调好。4.3 频率恢复与有功功率均分环节频率恢复比电压恢复更微妙因为频率是全局量所有DER通过公共母线天然连在一起无法单独改变某台变压器的频率。如果某台DER的有功出力不足它的频率就会偏低但整个系统频率在任意时刻几乎相等所以需要各台DER的有功功率按容量比例分配。有功功率均分通常是在一次控制的下垂系数m_i中体现容量大的DER下垂系数小容量小的下垂系数大。二次控制则是在一致性估计的基础上给出有功修正量。更高级的做法是直接控制有功功率参考值使P_i/S_i趋于一致。我在模型里采用的方法是在频率PI调节器的输出基础上叠加一个功率均分误差项误差系数kp_p0.2。这个误差项能动态调整下垂曲线加快均分收敛。调试时需要注意频率PI的输出限幅大小会直接影响均分效果。如果限幅太紧即使频率恢复得很完美功率均分也可能卡在某个误差上。我试验了限幅±1.0Hz与±0.5Hz两种情况前者均分偏差在5%以内后者均分偏差超过12%。所以在这类协同控制里不能只看频率波形要同时观察有功均分指标。4.4 协同控制器的参数整定方法整定协同PI参数我的通用方法是先按经典频域法估算初值再结合仿真微调。以频率环为例开环截止频率设定为2Hz相位裕度不小于45度。对于积分型的一致性控制器可以近似根轨迹法估算。实际仿真中我给出一个更接地气的整定顺序将事件触发阈值设置为极小值让模型近似连续通信。固定电压环比例系数为0.01频率环比例系数为0.2积分系数先设1。观察系统在1秒投入冲击负载后的响应如果频率振荡持续超过1.5秒降低比例系数如果稳态误差不为零增大积分系数。电压环重复同样的流程重点关注无功功率超调量目标超调不超过5%。参数稳定后再把事件触发阈值调回设定值根据通信次数和性能指标做最终折中。这套方法虽然不如严格的理论整定那么严密但在工程仿真里很实用。我最终调整出一组能兼顾动态性和通信频率的参数电压环kp0.05, ki8频率环kp0.3, ki5功率均分项kp_p0.2事件触发电压阈值0.5V频率阈值0.01Hz。5. 完整Simulink仿真模型搭建流程5.1 模型分层搭建与子系统封装我的Simulink模型根目录下分成三个大块主电路、控制层、通信与触发层。主电路直接放在根层方便看到母线电压和电流。控制层按照4台DER封装成四个子系统子系统里包含下垂控制、电压电流双闭环、PWM发生器和状态测量。通信与触发层单独建立一个子系统把所有事件触发模块、一致性算法模块都放进去输入来自四个DER的状态信号输出修正量再送回各个DER子系统。有人可能问为什么控制层不把分布式控制器也集成到DER子系统里我的经验是分布式控制器每个节点都要用到邻居信息如果放在DER内部就需要跨子系统的信号连线一旦拓扑变化就得重画很多线。单独做成一个“集中通信分布式计算”的模块可以清晰地区分本地控制和网络控制。当然从严格意义上说真正的分布式通信不应该通过Goto/From全局标签来实现但仿真阶段这样做调试效率最高。5.2 事件触发模块与控制器的信号连接在连接事件触发模块和控制器时最需要注意的是防止代数环。事件触发模块内部用Memory存储上一次发送值如果控制器计算得到的修正量又被用来影响测量值而测量值又反馈到事件触发模块就可能形成代数环。Simulink会报“Algebraic loop detected”错误。我的解决办法是在测量值进入事件触发模块之前先串一个Unit Delay单位延迟模块打破瞬时反馈回路。这个延迟步长很小比如1e-4秒对控制效果基本没有影响却能让仿真稳定解算。另外在控制器输出叠加回一次控制参考时也建议加一个Rate Limiter防止触发瞬间输出突变导致PWM饱和。信号连接顺序我建议按以下步骤明确各DER的测量信号来源电压幅值、频率、有功、无功。将测量信号接入事件触发子系统的输入端口。事件触发子系统输出打包后的向量到一致性算法模块。一致性算法模块返回估计的电压平均值和频率平均值。计算偏差、经过PI调节、限幅后输出δV和δω。将δV和δω通过Goto标签发送到DER子系统内部的加法器叠加到下桂参考上。5.3 仿真参数配置与解算器选择这个模型涉及PWM开关平均开关频率10kHz所以仿真步长必须足够小。我使用的是可变步长ode23tb最大步长设置为1e-4秒相对容差1e-4。用固定步长的离散解算器也能跑但速度慢很多而且大模型容易发散。ode23tb对带开关器件和间断信号的控制系统比ode45更稳定这也是微电网仿真圈里比较公认的经验。仿真时长设置为4秒。0-0.5秒为启动阶段0.5秒时投入第一组负载1.2秒时投入冲击负载2.5秒时切除冲击负载4秒结束。这个时间轴能覆盖动态响应、稳态恢复和负载突变多个工况。为了减少数据存储压力可以只在需要观察的模块上添加To Workspace记录点比如母线电压、各DER频率、事件触发次数计数器。模型初始化脚本我单独放在一个.m文件里定义基础参数物理量、控制参数、事件触发阈值、通信权重矩阵。脚本先运行再打开模型这样Simulink可以直接读取工作区变量。我习惯用结构体打包参数比如sys.L1.8e-3, sys.C50e-6这样模型里直接写sys.L就能引用改动参数也很方便。热词里提到“matlab simulink 输入变量是结构体的形式”这就是典型应用场景。6. 仿真结果分析与问题排查6.1 关键观察波形与性能指标跑完仿真后我重点观察四组波形公共母线三相电压、4台DER的频率曲线、有功功率分配曲线、事件触发次数记录。三相电压波形主要看幅值和THD。在事件触发机制下二次控制修正量是离散跳变的电压幅值会有小台阶但正常情况下不影响THD。我的模型稳态THD在2.8%左右加入非线性负载后升到4.1%还在可接受范围。频率曲线在负载突增时会下跌约0.4Hz二次控制大约在0.8秒内将其恢复到50Hz超调小于0.05Hz。有功功率曲线在2秒左右达到均分均分误差在6%以内。事件触发次数方面4秒仿真时长内4个DER总共触发约320次相比固定周期通信的9600次减少了96.7%。这个数据很能说明事件触发机制的价值。6.2 常见问题与排查技巧我在搭建模型时遇到了几个典型的坑都记录下来。第一个坑是Bus Selector没有可选信号。出现这个问题的原因通常是Bus Creator和Bus Selector之间的Goto/From标签在跨子系统时没有正确关联或者信号名中包含空格/中文导致名称匹配失败。解决办法是在Bus Creator的输出端用Bus Object.bus类型显式定义信号结构并在Goto/From属性里设置Signal Copy而不是Global就能稳定传递。第二个坑是Memory模块初值错误。如果Memory模块初值设为零而测量值初始不为零那么第一步就会触发事件导致系统启动时就有一堆无意义的通信。解决方法是初值设定为额定值比如380/50或者用Initial Condition模块预处理。第三个坑是Men希望事件触发逻辑中“面”的问题用Compare To Constant和Switch组合时如果触发脉冲为1的时间只有一个仿真步Switch输出可能会因为步长太短而无法更新。解决方法是触发后加一个MinMax或Hold模块延长逻辑电平保持一个控制周期或者把Switch的输出接回零阶保持器ZOH而不是普通Memory。第四个坑是ode23tb求解器在事件触发瞬间容易报“Newton iteration failed”。这是因为事件触发引起的信号跳变很锐利导致积分步长收缩极度频繁。我通过把最大步长从1e-4改到2e-5虽然能解决但仿真速度大幅下降。更好的办法是在跳变信号输出端加一个小惯性环节时间常数5e-5秒既不影响控制效果又能让求解器稳定通过。6.3 事件触发参数对系统性能的影响这个主题值得单独拿出来说。因为很多人做完基础模型后会想问“阈值到底调多大才合适”。我在第3.1节提过三组测试这里进一步解释。触发阈值本质上是给控制器“感知误差”设置了一个死区。死区小控制器反应敏捷但通信频繁可能会放大测量噪声。我在电压测量链路上加了高斯噪声SNR 40dB来模拟真实传感器当电压阈值设为0.5V时噪声引起的误触发率只有2%当阈值降到0.1V时误触发率飙升到35%通信次数直接翻倍而恢复速度提升不到10%。这说明阈值不能一味减小还要匹配传感器噪声统计特性。触发阈值也会影响一致性算法的收敛精度。事件触发机制相当于在通信中引入了“事件量化误差”这个误差的方差和阈值平方成正比。如果阈值过大分布式估计值会停留在某个稳态偏差上二次控制无法完全消除静态误差。我的测试结果阈值设为1V时稳态电压误差达到2V勉强算合格阈值设为0.5V时稳态误差0.3V效果满意。频率阈值0.02Hz时稳态频率误差0.002Hz差别不大。所以频率阈值可以适当放宽电压阈值要保守一些这正好对应我选择0.5V/0.01Hz组合的原因。7. 实测下来的一些心得做这个模型最大的收获不是把控制算法跑通而是明白了“通信资源和控制性能”之间的权衡在仿真里如何量化。事件触发机制并不是简单地把固定周期通信变稀疏而是要同时设计触发条件、阈值、通信拓扑和控制器结构。在Simulink里这四者相互耦合任何一环没弄好波形都会给你颜色看。单独提醒一个细节事件触发计数器的统计。如果想要在仿真报告里展示通信次数建议用Triggered Subsystem套一个计数器模块触发边沿设为rising每次事件触发计数器加1。通过To Workspace输出到MATLAB工作区后用numel(find(diff(cnt)0))就能得到总触发次数。这个统计比手动数示波器上的脉冲靠谱得多。如果后续要做更深入的扩展可以尝试把事件触发机制和模型预测控制结合让控制器在触发瞬间计算出一条最优的修正轨迹而不是仅仅输出一个PI增量。另外通信时延和丢包加入后事件触发机制会面临额外的稳定性挑战这也是目前学术界很关注的方向。我计划下一步就在这个基础上把异步触发和时延鲁棒性做进模型里。到时候有心得再继续分享。