DC微电网这几年在实验室和工程里都特别热但真正动手在Simulink里搭一套完整的直流微电网模型时很多人第一反应是这玩意全搭在一起跑得起来吗。风机、光伏、蓄电池、并网变流器每一块单独拿出来都有现成demo可一旦把它们挂到同一条直流母线上就出现各种奇怪问题母线电压振荡、MPPT不收敛、蓄电池功率方向乱跳、仿真步长小到跑不动。我这次把整套模型从零搭了一遍风机用MPPT、光伏用MPPT、蓄电池做母线稳压、再加一路并网接口把控制逻辑和参数整定都打通了。这篇就记录搭建过程中的关键设计思路和踩坑点给要做直流微电网仿真的朋友一条能直接参考的路径。阅读建议这篇偏实操适合已经会Simulink基本操作、但刚开始接触微电网或分布式发电控制的人。如果你只是需要交一份能跑的仿真模型里面每个模块的搭建顺序、参数初值、调试方法都可以直接照搬。1. 直流微电网的接线图母线架构与各单元接口选型1.1 母线电压等级怎么定直流微电网说白了就是一条公共直流母线所有电源和负荷都往里送功率、从里面取功率。第一步不是拉模块而是定母线电压。选多少伏直接决定后面每台变换器的拓扑和控制方案。实验室小规模系统常用48V模拟真实现场用400V或750V。我这次选的是400V原因很实际市面上600V等级的IGBT/MOSFET模块好找400V母线对器件的耐压要求适中而且和常见交流380V系统的整流输出电平能自然接轨后面做并网变流器时电压匹配不用额外加变压器。定了母线电压各单元接口方式基本就清楚了单元接口拓扑方向风机交流输出不控整流 Boost升压单向源→母线光伏直流输出Boost升压单向源→母线蓄电池双向Buck-Boost双向可充可放并网变流器三相PWM整流/逆变双向母线↔交流电网如果你的风机选的是永磁直驱型输出是频率变化的交流电那前面得加一个不控整流桥变成直流再进Boost。如果风机模型你已经简化成机械功率输入发电模块输出直流那可以直接省掉不控整流这一级但牺牲的是风速动态时的真实响应。我后面详说。1.2 公共接口设计别让子系统之间互相不知道对方是谁Simulink建模最忌讳把所有模块铺在一个画布里。每个单元都应该封装成独立子系统对外只暴露必要的物理接口和信号接口。我的习惯是给每个子系统定义统一的接口协议电气接口母线的电压、电流用物理信号线连接控制信号MPPT输出的占空比、使能信号、运行模式测量输出往里输出功率、电压、电流、SOC等观测信号。这样每个子系统可以单独调试单独跑通后再挂到公共母线上。我见过很多人省掉这一步直接在总模型里堆模块最后一个问题查三天查不出来。1.3 公共母线电容的作用母线电容是整个微电网的压舱石。瞬态功率不平衡时母线电压波动主要靠电容吸收或释放能量来缓冲。电容选多大有一个粗略估算公式C P × Δt / (ΔU × U_bus)比如系统额定功率5kW允许瞬时功率缺额持续0.1秒母线电压400V允许波动5%20V那电容量级大约在C 5000 × 0.1 / (20 × 400) ≈ 62.5mF实际仿真中这个值会偏大但一开始可以按这个数量级给等调试稳定以后再降下来。电容太小母线电压会高频抖动电容太大电压建立变慢而且并网启动瞬间的冲击电流也大。2. 风机MPPT的Simulink实现从气动特性到Boost占空比2.1 风轮机建模功率系数的决定性作用风机MPPT听起来是控制问题但控制的对象是风轮机发电机的机电特性所以模型要先有风轮机。风轮机输出的机械功率表达式是P_mech 0.5 × ρ × π × R² × V_w³ × Cp(λ, β)其中ρ是空气密度R是风轮半径V_w是风速Cp是风能利用系数。Cp不是常数它和叶尖速比λ、桨距角β有关。最常见的简化模型是λ ω_r × R / V_wCp(λ, β) 0.5176 × (116/λ_i - 0.4β - 5) × exp(-21/λ_i) 0.0068λ这里λ_i是一个修正后的中间变量1/λ_i 1/(λ 0.08β) - 0.035/(β³ 1)这几个公式看着复杂其实就是在复现一条以λ为横轴、Cp为纵轴的单峰曲线。峰值对应最优叶尖速比λ_opt在这个点运行捕获的风能最大。Simulink里直接在Fcn模块里写公式即可也可以用Matlab Function块封装一段function。2.2 MPPT算法选择TSR法和PO法对比风机MPPT主流有两种思路一是最优叶尖速比法TSR。它需要实时测量风速V_w和风轮转速ω_r算当前λ然后和λ_opt比较误差经过PI调节器输出转速或功率参考值。这个方法响应快但问题是要测风速——实际风机机舱上的风速仪测的是受扰流影响的局部风速并非真实来流风速另外风速仪故障率高。仿真里不存在这个问题所以TSR法特别适合验证模型逻辑。二是扰动观察法PO。它不测风速而是不断对电压或占空比施加小扰动观察功率变化方向功率增加就继续同方向扰动功率减小就反向扰动。好处是不需要传感器坏处是稳态时功率点会在最大值附近来回振荡风速突变时还可能朝错误方向爬。我在模型里用的是TSR法做基础验证因为逻辑清晰、收敛快适合先把整体系统跑通后续要模拟真实工况再改成PO。TSR法的风机MPPT控制框图风速 V_w → 计算 λ → 与 λ_opt 比较得到误差 → PI 控制器 → 输出功率参考值 P_refP_ref 与母线电压测得的 V_bus 相除得到电流参考 I_refI_ref 与 Boost 输入电流比较 → 第二个 PI → 输出占空比 d或者更常见的做法直接用 MPPT 输出一个功率参考通过功率/电压关系折算成电流环指令再进电流内环。Simulink里用两三个PI控制器就能搭完。2.3 系统的惯性折算风电转动惯量不能省很多初学者搭风机MPPT时直接在风轮机输出端挂一个受控电流源节省了大量建模时间但代价是忽略了机械惯性。风机本身的转动惯量是很大的缓冲风速突变时转速不会突变功率输出是缓慢变化的。如果你省掉惯性风速阶跃时功率会瞬间跳变这会给母线电压带来巨大的冲击仿真出来的控制效果差其实不是MPPT的问题而是你丢掉了物理惯性。我在风轮机模型里加了传动链的一阶惯性环节用时间常数T_mech J / D近似模拟J是实际风机折算后的转动惯量D是阻尼系数。这样风速从8m/s阶跃到10m/s时功率是斜坡式上升的和现场实测更像。2.4 风机侧Boost的占空比限制风机MPPT输出的占空比一定要限幅否则启动瞬间或有扰动时占空比冲到100%会导致电感电流失控。我的限幅范围是0.1~0.9而且加了rate limiter限制占空比变化率不超过每秒0.2左右。这个变化率限制就是给MPPT的功率变化套一个软约束避免大幅度占空比跳变。实际仿真中如果你不加占空比变化率限制MPPT在风速剧烈变化时会出现明显的功率过冲甚至触发代数环收敛问题。这不是控制器参数不好调而是执行机构变换器的物理响应速度被忽略了。3. 光伏阵列与Boost升压MPPT控制在工程模型上的落地3.1 光伏组件的工程等效模型光伏阵列的I-V特性曲线是强非线性的Simulink里有Simscape Electrical的PV模块可以直接用但我更推荐自己用公式建模。原因有二一是自建模型可以精确控制参数适配你拿到手的实际组件数据二是运算速度快不会因为Simscape的细节模型让整个系统仿真变慢。光伏电池单体的输出方程I I_ph - I_0 × [exp(q(V IR_s)/(nkT)) - 1] - (V IR_s)/R_sh这个方程是隐式方程工程上一般用简化双参数模型或四参数模型通过开路电压V_oc、短路电流I_sc、最大功率点电压V_mp和电流I_mp这四个厂家标称参数拟合。在Simulink里用Matlab Function块实现输入光照S和温度T输出I-V曲线上对应的工作点。光照变化对Isc影响大温度变化对Voc影响大这是调MPPT时必须记住的物理规律。光照从1000W/m²降到600W/m²最大功率点电压变化不大但短路电流大幅下降温度从25°C升到45°C开路电压明显下降最大功率点的电压区间会整体左移。3.2 扰动观察法的状态机实现我在光伏侧用的是电导增量法Incremental Conductance它比PO的稳态振荡小误判概率低。核心判据是dP/dV 0 → 在最大功率点dP/dV 0 → 工作在左侧需要增大电压dP/dV 0 → 工作在右侧需要减小电压实际离散化实现时判断dI/dV和-I/V的关系即可。Simulink里实现时注意采样保持整个MPPT以固定频率执行比如每0.01秒执行一次电压电流测量值更新后才允许占空比改变。3.3 PWM调制和Boost电路参数Boost电路的电感和电容参数直接影响纹波和动态响应。我按连续导通模式设计L V_out × (1-D) / (2 × f_sw × I_L_min)其中f_sw是开关频率。比如母线400V光伏侧200VD0.5开关频率20kHz最小电流5A那么L要大于L 400 × 0.5 / (2 × 20000 × 5) 1mH实际取3~5倍裕量我用的是5mH。输出电容按电压纹波约束取C P_out / (f_sw × ΔV × V_bus)开关频率我选在20kHz这个频率在仿真步长和真实器件特性之间取了个平衡。低于10kHz纹波大MPPT会因采样噪声而误判高于50kHz仿真步长必须压到很小模型整体跑不动。3.4 光伏MPPT的启动行为光伏系统启动时光照刚施加输出电压从零开始建立。这时候如果MPPT立刻工作占空比会乱跳。我的做法是加一个启动判断只有当光伏输出电压大于一个阈值比如50V后MPPT才开始运行。这模拟的是现实中变换器的软启动过程避免启动瞬间冲出一个巨大的浪涌电流。4. 蓄电池模型与双向Buck-Boost储能端如何稳定母线4.1 电池等效电路选型蓄电池是这个系统里唯一能吸收和释放能量的单元。模型太简单SOC算不准模型太复杂仿真速度拉胯。我采用的是带内阻和一阶RC网络的戴维南等效模型V_bat V_oc(SOC) - I × R_0 - V_1其中R_0是欧姆内阻V_1是极化电容C_1两端的电压受极化电阻R_1和电流决定。V_oc作为SOC的函数通过查表给定常见锂离子电池的V_oc-SOC曲线大致是SOC在10%~90%区间平缓、两端陡峭。SOC计算采用安时积分法SOC SOC_0 - (1/Q_nom) × ∫I_bat dtQ_nom是电池额定容量单位Ah。放电时I为正SOC下降。要注意仿真开始时SOC初值要设对比如0.5或者0.8如果从0开始充控制逻辑很容易触碰充电下限。4.2 双向Buck-Boost的四种工作状态双向变换器本质上就是Buck和Boost两种拓扑的合并。蓄电池端电压低比如300V母线端400V那么放电蓄电池→母线Buck-Boost工作在Boost模式能量升压流向母线充电母线→蓄电池变换器工作在Buck模式降压给电池充电。实现上有两种做法一种是两个开关管互补导通用同一个占空比控制电流方向自动决定模式另一种是两个独立开关分别控制。工程上第一种更常用控制结构简单。但为了避免上下桥臂直通必须设置死区时间比如2微秒。在Simulink的离散仿真里死区逻辑可以用On/Off Delay模块实现也可以用Stateflow搭。4.3 蓄电池的控制策略恒压稳定母线、限流保护我把蓄电池作为整个微电网的平衡节点——谁多了谁少了蓄电池负责补差。控制结构是电压电流双闭环外环母线电压参考值400V和实测母线电压比较经过PI输出电流参考内环蓄电池电流跟随电流参考经过PI输出占空比。电压环PI参数整定时要注意母线电容的存在让电压环被控对象近似一阶惯性积分I增益过大会导致母线电压低频振荡。我现在的参数是电压环Kp2Ki40电流环Kp1.5Ki100采样周期都是0.1ms跑起来母线电压波动在±2V以内。限流保护必须加。蓄电池最大充放电电流一般为0.5C比如100Ah的电池最大电流50A。我在电流参考后面加了一个Satulation模块上下限设为±50A。同时加了一个SOC上下限判断逻辑SOC低于10%时只允许充电高于95%时只允许放电。这个逻辑用Matlab Function块写很简单但一定要和主控制器的使能信号联动否则保护动作后控制器还在输出占空比会顶死在饱和限。4.4 为什么蓄电池侧容易出功率倒灌调试中一个典型问题是光伏和风机出力大于负荷需求母线电压被抬高蓄电池电压外环会给出负的电流参考充电这没问题。但如果你并网变流器还在向电网输送功率两者抢着消化盈余母线电压就可能出现拍频振荡。解决办法是给并网变流器和蓄电池设置优先级正常并网运行时优先通过并网变流器送回电网蓄电池只处理电网来回波动的小功率分量。这个逻辑需要在协调层做后面专门讲。5. 并网运行的核心母线电压调节与多单元协调策略5.1 并网变流器在直流微电网里的定位直流微电网并网指的是通过一个AC/DC变流器把直流母线和交流大电网连接起来。这个变流器可以工作在整流状态从电网吸收功率补充直流侧也可以工作在逆变状态把直流侧盈余功率送回电网。控制上常用电压定向矢量控制VOC直流侧电压外环交流电流内环外环直流母线电压参考V_bus_ref与实测V_bus比较经过PI得到有功电流参考i_d_ref内环d轴电流和q轴电流PI调节输出d/q轴电压指令最后经过Park反变换和SVPWM产生开关信号。这套框图如果你做过三相整流器就很熟悉。在微电网里它就是直流母线和交流母线之间的功率通道。5.2 并网模式下的母线电压谁说了算这是整个系统协调的核心问题。直流母线电压的稳定必须由某一个单元来保证否则各单元都按自己的MPPT跑母线电压没人管系统一定会崩。我把控制模式分成两种离网模式并网变流器不工作蓄电池作为主控单元定电压控制母线风机和光伏都做功率源并网模式并网变流器作为主控单元稳压母线蓄电池退居辅助只做短时功率缓冲。这两种模式的切换通过一个逻辑判断实现检测到交流电网侧电压、频率正常则并网变流器使能同时蓄电池从定电压模式切换为定功率/下垂模式。这个切换如果处理不好母线电压会瞬间跳变。我在切换时加了一个比例积分前馈补偿切换前记录当前母线电压和蓄电池电流切换后把电流参考初始化为切换前的值避免控制器饱和导致冲击。5.3 下垂控制在多单元协调中的作用如果未来要扩展多台蓄电池或多台并网变流器并联运行定电压控制就不合适了——多个主控单元会因为参考电压差异互相环流。这时候要用下垂控制Droop ControlV_bus_ref V_0 - k_d × I_out下垂系数k_d决定了功率分配比例。比如两台变换器k_d相同则输出电流自动均分。这个思想的本质是让每个单元的电压参考值随输出功率增大而略微下降利用母线电压作为通信总线协调各单元功率。我在模型里把蓄电池预设了一个下垂系数初始值k_d0.5V/A即输出20A时电压参考下降10V。下垂控制的好处是不需要单元间通信即插即用坏处是母线电压静态偏差变大所以后面通常还要加二次调压。仿真时可以先只跑下垂观察母线电压的静差是否符合预期再加入二次调节。5.4 并网启动的预充电问题直流母线电容直接接到400V上启动瞬间的冲击电流非常大仿真中会让模型进入数值奇异。我的做法是加预充电电阻先通过限流电阻给母线电容充电母线电压升到95%以后再让主接触器吸合切除预充电回路。Simulink里实现很简单用Breaker或可变电阻时序交给Stateflow状态机控制。这一步跳过的人后面多半会遇到电压明明建立了但总报警过流的怪问题。6. 参数整定、求解器选择与调试避坑让模型从能跑到跑稳6.1 求解器选型连续域和离散域不要混战Simulink微电网模型最烦人的问题就是仿真速度慢。变流器工作在20kHzPWM周期50微秒如果你用固定步长求解器步长至少0.5~1微秒才能看清开关细节仿真10秒就要跑几千万步基本跑不动。我的实用方案是场合推荐设置说明整体控制策略验证变步长 ode23tb误差1e-3开关管用平均值模型速度快开关级波形分析定步长1e-6sode3只跑几毫秒观察细节含电力电子开关的离散模型Powergui离散化采样1μs适合硬件在环前的验证第一阶段的策略验证里我把所有DC-DC变流器都用平均值模型替代不再建模开关管和PWM用一个受控电压源/电流源来等效变换器的平均行为。比如Boost变换器输入侧用受控电流源模拟电感电流输出侧用受控电压源模拟经过升压后向母线输送的电压。这样模型步长可以放到50微秒10秒仿真十几分钟就能跑完速度提升了两个数量级。6.2 代数环微电网模型最常见的仿真杀手代数环的典型表现是仿真刚开始就报Simulink cannot solve algebraic loop或者卡住不动。微电网里代数环最常见的位置是电压测量 → 控制计算 → 占空比 → 功率源 → 又影响电压测量整个回路里没有延迟环节。破解方法有三种在反馈测量信号后面加一个Memory模块或者Unit Delay打破代数环。用变步长求解器时在配置参数里勾选Allow algebraic loop solvers。控制链路上用离散采样器Zero-Order Hold让测量值以固定频率更新。我推荐第三种因为这和真实数字控制器的行为一致——DSP采样是离散的控制也是离散的。把控制环改成离散域以后代数环基本消失了而且模型的物理意义更贴近真实系统。代价是需要整定离散PI参数和连续域PI参数略有差异。6.3 PI参数初值的实用整定法微电网控制里PI参数整定我很少用理论公式硬推更常用的是基于物理对应关系的试凑以蓄电池电压外环为例思路是这样先只保留内环电流环给定一个固定的I_ref调k_p让电流跟随不用管超调先记录响应速度然后单独用比例环节调电压环逐渐增加k_p直到母线电压出现等幅振荡此时k_p约为临界增益k_cr按照Ziegler-Nichols经验取k_p≈0.45k_crT_i≈0.83T_crT_cr是振荡周期然后微调。实测下来这个流程比随机试参数快得多。要注意的是离散控制器的采样时间越小稳定裕度越大但抗噪声能力越差。我一般控制在0.1ms采样也就是5个PWM周期更新一次控制量这样实现成本低稳定性也够。6.4 数据观测与故障复现把Scope换成Logger调试微电网时Scope窗口只能看个大概波形要分析功率流向、SOC变化、多组信号对比时非常不方便。我后期全部改成Simulink的数据记录器Simulation Data Inspector每个子系统把关键变量导出为结构化数据母线电压、光伏功率、风机功率、蓄电池电流、SOC、并网电流。一次仿真跑完可以在SDI里同时拖拽6条曲线用游标精确读取数据方便发现振荡频率和相位关系。如果想记录被试控制策略在某个故障下的行为可以设置多个故障触发时刻比如t2s: 光照从1000W/m²阶跃到900W/m²t4s: 风速从9m/s阶跃到8m/st6s: 负荷突增2kWt8s: 并网开关断开切换离网模式把这几个工况放在同一个仿真脚本里跑观察母线电压的动态恢复时间和波动范围基本能完整评估控制策略的性能。6.5 我自己调试时踩过的一个最深的坑最后必须分享一个我这次搭建中踩得最深的坑蓄电池SOC积分漂移导致系统越跑越假。原因很简单安时积分法计算SOC时电流测量中的直流偏置会被不断累积10秒仿真看起来没事跑到100秒时SOC已经漂了3%。在长工况仿真里这个漂移会让蓄电池的充放电限制逻辑提前触发然后整个系统的功率分配就乱了。解决办法有两个方向一是定期用开路电压法校正SOC但运行中不能开路二是限制安时积分器的输出范围同时给电流测量信号加一个高速滤波滤掉直流偏置。工程上更稳妥的是用扩展卡尔曼滤波估计SOC但仿真阶段用限幅滤波就够了没必要上来就上EKF。另外一个容易忽略的是蓄电池初始SOC和开路电压必须匹配。如果你设SOC0.5但V_oc表格里0.5对应的电压是310V而仿真初始化时电池电压给的是280V那么启动瞬间会有一个很大的过渡电流冲击。检查和初始化状态一致性是仿真跑不跑得稳的隐藏关键点。数据记录和参数整定做扎实以后这套模型基本就稳定了。最后再分享一个个人体会微电网仿真里系统能跑和系统跑得对之间隔着一整条调试路不要一看到波形稳定就认为控制策略正确一定要在工况切换瞬间仔细看电压恢复时间和功率分配是否合理。后续如果你想扩展这个模型可以在现有基础上加负荷预测、加故障穿越逻辑、甚至把某个单元换成硬件在环整个架构都是开放的不会推倒重来。