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做超短期预测时,通常用风速的滑动平均加波动分量来近似。
三、储能电池与超级电容:能量型与功率型的分工
储能电池的稳态模型核心是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、γ、η_inv | 0.1~0.5 s | 日前/日内出力预测 |
| 风机 | v_cut_in、v_rated、v_cut_out | 1~5 s | 超短期功率预测 |
| 储能电池 | E_rated、SOC_min/max、η | 0.1~1 s | 能量搬移、调峰 |
| 超级电容 | C、U_max、ESR | 0.01~0.1 s | 高频功率支撑 |
| 柴油发电机 | P_min、爬坡率、油耗系数 | 1~3 s | 旋转备用、孤岛主源 |
| 燃气轮机 | 启停时间、最小运行时间 | 5~20 s | 热电联供、长期调度 |
| 燃料电池 | 电效率、最小运行时间 | 0.5~2 s | 清洁基荷、备用 |
| 负荷 | ZIP系数、再启动倍数 | 秒级~分钟级 | 平衡约束、切换校验 |
实际工程中,元件模型参数应来自设备铭牌、出厂试验报告或现场辨识,而不是随意假设。EMS投运前,建议用历史运行数据对模型做一轮回归校验:把实际辐照度、风速、负荷代入模型,比较预测功率与实测功率的偏差,若光伏模型日均偏差超过5%、风机超过10%,就需要重新标定功率曲线或效率系数。只有模型能复现历史运行,EMS的调度决策才具备可信的物理基础。