3. 微电网元件建模

微电网EMS的核心任务是在秒级到分钟级的时间尺度上完成功率平衡、经济调度与安全约束管理。EMS要做出正确决策,前提是它“认识”每一个可控与不可控元件——这种认识就是元件模型。元件模型不是仿真软件里的黑箱,而是EMS内部用于预测、优化和校验的数学表达。同一个物理设备,在稳态调度和动态稳定分析中需要不同粒度的模型:稳态模型回答“这个功率点能不能长期维持”,动态模型回答“从当前状态切换到目标状态需要多久、会不会越限”。

一、建模的基本约定与统一接口

在微电网EMS中,所有元件模型最终都要落到三个接口量上:有功功率P、无功功率Q、以及反映能量状态的变量(如SOC、转速、温度)。本文统一约定:功率正方向为元件向交流母线注入功率,负荷吸收功率为负;时间步长在稳态调度中取15 min或1 min,在动态分析中取10 ms~1 s。建模时先确定“状态变量—代数变量—控制变量”三类量的划分,再写出功率平衡方程和状态演化方程。

稳态模型通常写成代数方程 P=f(u)、Q=g(u),其中u是环境或指令变量;动态模型写成微分-代数方程组 dx/dt=h(x,u)、P=k(x,u)。EMS中常用的简化策略是:对慢过程(电池SOC、储热)保留动态,对快过程(逆变器电流环)用一阶惯性或纯增益近似。

二、光伏与风机:不可控源的稳态与动态

光伏的稳态模型以最大功率点跟踪(MPPT)后的出力为核心。给定辐照度G(W/m²)和电池温度T(℃),单块组件的直流功率可近似为:

P_dc = P_stc × (G/1000) × [1 + γ×(T - 25)]
γ ≈ -0.0035 /℃(晶硅组件温度系数)

再经逆变器效率η_inv(典型0.96~0.98)得到交流侧P_ac。EMS中常把光伏聚合为“辐照度—出力”分段线性曲线,便于线性规划求解。动态方面,光伏本身无机械惯性,但逆变器有功指令响应可用一阶惯性描述:τ_pv·dP/dt + P = P_ref,τ_pv约0.1~0.5 s。当EMS下发限功率指令时,实际出力按该惯性逼近目标值,这是判断AGC指令能否被跟踪的依据。

风机稳态模型以风速v为输入,经功率曲线得到P_wt:

v < v_cut_in(3 m/s):P = 0
v_cut_in ≤ v < v_rated(12 m/s):P = P_rated × (v³ - v_cut_in³)/(v_rated³ - v_cut_in³)
v_rated ≤ v ≤ v_cut_out(25 m/s):P = P_rated
v > v_cut_out:P = 0

动态上,风机转子惯量较大,功率响应时间常数约1~5 s,且存在塔影、风剪切引起的功率波动。EMS做超短期预测时,通常用风速的滑动平均加波动分量来近似。

微电网元件模型:稳态约束与动态响应分层 输入量 辐照度 G 风速 v 负荷需求 P_L 电价 / 调度指令 SOC 初值 燃料指令 母线电压/频率 元件模型(稳态 + 动态) 光伏 PV 稳态:P=f(G,T) 动态:τ_pv≈0.1~0.5s 风机 WT 稳态:分段功率曲线 动态:τ_wt≈1~5s 储能电池 BESS 稳态:SOC 约束 动态:一阶 RC / 内阻 超级电容 SC 稳态:E=½CU² 动态:τ≈0.01~0.1s 柴油发电机 / 燃气轮机 / 燃料电池 稳态:P_min ≤ P ≤ P_max,爬坡约束 动态:柴油 τ≈1~3s;燃气 τ≈5~20s;燃料电池 τ≈0.5~2s 负荷 Load 稳态:ZIP 模型 P=P₀(a·V²+b·V+c);动态:电动机再启动 EMS 使用 功率平衡方程 经济调度 SOC 轨迹预测 爬坡可行性校验 频率支撑评估 备用容量计算 切换逻辑判断 时间尺度: SC 10ms 电池 100ms 柴发 1s 风机 5s 光伏 0.5s

三、储能电池与超级电容:能量型与功率型的分工

储能电池的稳态模型核心是SOC演化与功率约束。以额定容量E_rated(kWh)、充放电效率η_ch/η_dis表示,离散时间步Δt内:

SOC(k+1) = SOC(k) + (η_ch × P_ch × Δt) / E_rated      (充电)
SOC(k+1) = SOC(k) - (P_dis × Δt) / (η_dis × E_rated)  (放电)
约束:SOC_min ≤ SOC ≤ SOC_max,|P| ≤ P_max

动态模型常用一阶RC等效电路:开路电压U_oc(SOC)串联内阻R₀和RC并联支路,端电压U = U_oc - I·R₀ - U_rc。EMS中若只做分钟级调度,可用“SOC—可用功率”查表代替完整RC模型;若要做一次调频评估,则必须保留R₀和极化环节,否则会高估电池的瞬时支撑能力。

超级电容的稳态能量为E=½C(U_max²-U_min²),其突出特点是功率密度高、循环寿命长,但能量密度低。动态上它可用理想电容串联等效串联电阻ESR描述,时间常数通常在10~100 ms量级。在微电网中,超级电容常与电池组成混合储能:高频功率波动由超级电容承担,低频能量搬移由电池承担。EMS的功率分配可用一阶低通滤波:

P_bat_ref = LPF(P_total, T_f),T_f 取 10~60 s
P_sc_ref  = P_total - P_bat_ref

四、柴油发电机、燃气轮机与燃料电池

柴油发电机稳态模型由燃料消耗曲线和电气约束组成。典型油耗曲线为F = a·P_rated + b·P_out(L/h),其中a为空载油耗系数,b为边际油耗系数。EMS做经济调度时,常把柴油机表示为分段线性成本函数,并施加最小出力P_min(通常30%额定)和爬坡率限制(如每分钟20%额定)。动态上,柴油机有调速器和励磁系统,有功响应时间常数约1~3 s,适合承担旋转备用。

燃气轮机响应更慢,从启动到满功率可能需要数分钟到数十分钟,且存在最小启停时间约束。其稳态模型与柴油机类似,但爬坡率更低(约每分钟5%~10%额定),EMS在日前调度中必须把启停成本纳入目标函数。燃料电池(如质子交换膜PEMFC)效率高、响应较快,电效率约40%~50%,热电联供时综合效率可达80%以上。其动态可用一阶惯性加纯延迟近似:τ_fc约0.5~2 s,延迟约0.1~0.5 s。EMS需注意燃料电池不宜频繁启停,通常设定最小运行时间。

五、负荷模型:ZIP与动态再启动

负荷不是简单的常数P。稳态负荷常用ZIP模型描述电压依赖性:

P = P₀ × [a × (V/V₀)² + b × (V/V₀) + c]
a + b + c = 1

其中a为恒阻抗分量,b为恒电流分量,c为恒功率分量。居民负荷中恒功率分量较大,工业电动机负荷中恒阻抗分量较大。动态负荷方面,空调、冰箱等温控负荷有热惯性,电动机类负荷在电压跌落恢复后会出现再启动冲击,冲击功率可达额定值的3~5倍。EMS在孤岛切换或大功率机组退出时,必须预留这部分再启动容量,否则会触发二次低频减载。

六、建模成果的校验与使用

模型建好后,至少要做两类校验。第一类是稳态校验:给定一组典型工况(如正午光伏满发、晚间无风),检查功率平衡方程是否闭合、SOC是否越限、机组是否在爬坡能力内。第二类是动态校验:施加一个阶跃扰动(如负荷突增10%),观察各元件响应时间是否与设定时间常数一致,母线频率最低点是否在允许范围内。下表给出一组教学示例参数,用于说明不同元件在EMS中的典型建模粒度。

元件稳态关键参数动态时间常数EMS典型用途
光伏P_stc、γ、η_inv0.1~0.5 s日前/日内出力预测
风机v_cut_in、v_rated、v_cut_out1~5 s超短期功率预测
储能电池E_rated、SOC_min/max、η0.1~1 s能量搬移、调峰
超级电容C、U_max、ESR0.01~0.1 s高频功率支撑
柴油发电机P_min、爬坡率、油耗系数1~3 s旋转备用、孤岛主源
燃气轮机启停时间、最小运行时间5~20 s热电联供、长期调度
燃料电池电效率、最小运行时间0.5~2 s清洁基荷、备用
负荷ZIP系数、再启动倍数秒级~分钟级平衡约束、切换校验
建模粒度要与EMS的时间尺度匹配:做15 min经济调度时,光伏和风机可用稳态功率曲线,电池用SOC递推即可;做一次调频或孤岛切换分析时,必须引入各元件的动态时间常数和功率限幅,否则优化结果会给出物理上无法执行的指令。

实际工程中,元件模型参数应来自设备铭牌、出厂试验报告或现场辨识,而不是随意假设。EMS投运前,建议用历史运行数据对模型做一轮回归校验:把实际辐照度、风速、负荷代入模型,比较预测功率与实测功率的偏差,若光伏模型日均偏差超过5%、风机超过10%,就需要重新标定功率曲线或效率系数。只有模型能复现历史运行,EMS的调度决策才具备可信的物理基础。