1. 为什么楼宇微网需要虚拟储能从现实着手的工程逻辑做楼宇微网优化调度的人一开始大概率都踩过同一个坑把储能系统当成整个优化模型的定海神针以为只要有电池削峰填谷、平抑波动这些事就都能解决。但我这两年跑下来的实际体会是电池的成本和寿命回收账很难算平尤其在楼宇这种场景里空间、消防、并网手续每一关都能让你怀疑人生。真正让我觉得这条路能走下去的反而是把需求侧那些看不见的弹性挖掘出来用虚拟储能去替代一部分物理储能的功能。这就要聊到需求侧虚拟储能系统的核心逻辑了。楼宇里的空调、暖通、热水系统本质上都是带热惯性的家伙。房间温度不会因为主机一停就瞬间失控热水箱里的热量也不会因为加热器断电就立刻散光。这种热惯性和舒适度允许的上下浮动区间就是可以拿来充放电的虚拟储能容量——电价高的时候让主机多歇一会儿用建筑本身蓄的冷或热顶着电价低的时候让主机多干一会儿把能量蓄在楼宇结构里。这么一来物理电池不必配那么大甚至在某些场景下可以完全不配。这个思路放到楼宇微网里就很有意思了。楼宇微网本身有光伏、有负荷、有跟大电网的交互关口再加上虚拟储能系统这个隐形电池调度模型就从单纯的源-荷-储匹配变成源-荷-储-虚拟储的四维协同。我在Matlab里把整套模型搭起来之后最直观的感受是分时电价下的运行成本肉眼可见地降了而且光伏消纳率也比纯物理储能方案要稳。如果你正准备做这个方向或者手里有楼宇微网的优化调度课题想落地这篇文章会比较适合你。我不打算泛泛讲理论公式而是把整个建模过程、Matlab代码实现的骨架、以及我在调试过程中踩过的坑按照实际操作顺序梳理出来。从虚拟储能模型的建立到目标函数和约束条件的处理再到YALMIPgurobi求解器在Matlab里的配置调用每一段都是可以直接拿去改的。2. 从闲置的热惯性到可调的虚拟储能建模第一步2.1 楼宇虚拟储能的物理基础热容和热阻虚拟储能不能凭空创造能量它的本质是用可接受的舒适度代价把能量转移到一个时间窗口内。要实现这个先得把楼宇的热动态用数学语言描述出来。我采用的是最常见的一阶等效热参数模型ETP模型相当于把整栋楼的围护结构、室内空气、空调系统简化成一组热容和热阻的关系。楼宇室内温度的动态方程可以写成C_air \cdot dT_in/dt (T_out - T_in) / R_wall Q_HVAC Q_occupancy Q_solar这里面C_air是室内空气的等效热容R_wall是围护结构的热阻Q_HVAC是空调系统向室内注入或抽走的热功率Q_occupancy和Q_solar分别是人员和太阳得热。离散化之后就变成相邻时段温度的递推关系。我在Matlab里直接用有限差分的方式处理时间步长取1小时跟调度周期保持一致。这组方程看起来简单但它决定了虚拟储能的核心边界。温度不能无限偏离用户设定的舒适区间因此虚拟储能的可调度范围本质上就是楼宇热质量在当前工况下能容忍的冷/热积累量。可以这样理解一栋楼就像一个大水杯水是热量杯子壁是热阻水位就是室温水位允许在某个高低线之间浮动杯子里的水量变化就是虚拟储能的充放电量。2.2 虚拟储能充放电功率与容量的折算我刚接触这个概念的时候最困惑的问题是虚拟储能的功率和容量到底怎么定物理电池有额定功率和额定容量虚拟储能怎么折算我现在的做法是把虚拟储能跟暖通负荷的解耦空间直接关联起来。楼宇的HVAC系统在某时段的基准用电功率是P_HVAC_base按常规控制策略维持设定温度所需功率虚拟储能的充电就是让HVAC在电价低谷时多出力把楼宇温度拉向舒适区边缘或略低一点向下的空间是它的充电深度放电就是HVAC少出力甚至停机楼宇温度向舒适区上边界漂移。于是虚拟储能的功率响应跟实际HVAC功率和基准功率的差值挂钩。容量方面我在Matlab里用温度可行域内可吸收或释放的总热量/能效比来折算。举个例子某栋楼在夏季制冷工况室温允许从24到26摄氏度浮动楼宇等效热容C_air如果取800 kWh/℃那么2℃温差对应的可调度热能量就是1600 kWh折算成电功率除以空调COP约3.5大概是450 kWh左右的虚拟储能容量。这个数值已经不小了相当于一个中型电池柜的容量但它几乎不需要额外的硬件投资。2.3 为什么选择一阶模型而不是更高阶我试过二阶甚至三阶等效热参数模型需要辨识的参数更多模型更精细但调度结果对楼宇舒适度边界的提升并不显著。原因在于调度的时间尺度是小时级而二阶以上的热动态差异主要反映在分钟级甚至更短时间尺度上的响应过程。对日前或日内滚动调度来说一阶模型在精度和计算负载之间是最划算的那一档。如果以后要扩展到实时控制层面再考虑更高阶或者数据驱动的替代模型也不迟。3. 楼宇微网日前优化调度模型的完整搭建过程3.1 决策变量时间维度上的分角色设置搭建优化模型的第一步是把决策变量理清楚。我在Matlab里用YALMIP定义变量时按照设备角色分成几组微网购售电功率P_buy(t)和P_sell(t)非负变量不能同时大于0光伏出力P_pv(t)在日前调度中作为预测参数而非决策变量因为光伏在当前阶段通常按最大功率跟踪运行物理储能如果有P_ess_ch(t)和P_ess_dis(t)以及SOC的连续状态变量HVAC系统功率P_hvac(t)这个变量直接跟虚拟储能挂钩室温状态T_in(t)连续变量用来约束虚拟储能的充放电深度这些变量的时间维度覆盖整个调度周期T24步长为1小时。需要注意的是T_in(t)和P_hvac(t)之间的耦合关系需要把我在第2节里讲的热动态方程作为等式约束写进模型这就把虚拟储能和楼宇微网的能量平衡真正绑在了一起。3.2 目标函数不止是省钱这一件事楼宇微网调度最常见的目标函数就是运行成本最小化但如果只写成本项求解出的方案可能在技术上可行、实际运行完全没法用。我在代码里把目标函数拆成三个部分第一是购电成本减去售电收益这是经济性的核心也直接受分时电价驱动。第二是舒适度惩罚项也就是室温偏离舒适区间中心值的代价用二次或线性惩罚函数实现。第三是物理储能的充放电损耗折算项避免模型为了追求电费最小化而让电池空耗循环。用权重系数把这几个目标加权求和模型求出的解才是既有经济性又不会牺牲用户感受的可执行方案。权重系数的设置我建议做成输入参数这样后续做敏感性分析非常方便。我在初始代码里成本和舒适度的权重比大概取100:1测试下来调度的温度曲线贴近舒适区间边界但不会越界。3.3 运行约束汇总哪些约束真正限制了可行域我梳理一下我在Matlab代码里实际用到的约束这些基本都是微网调度模型的标准件但每类约束的细节处理方式会影响求解效果功率平衡约束是必须的P_pv(t) P_buy(t) P_ess_dis(t) P_hvac_consume_base(t) P_load_other(t) P_sell(t) P_ess_ch(t) P_hvac_flexible(t) ...实际上我在实现时是把所有负荷和电源放在一个等式里用总供给总需求来表达。这个等式如果不小心把虚拟储能的功率项重复计入就会出现能量凭空消失或凭空产生的现象我一开始就踩过这个坑。储能约束方面SOC递推式和充放电功率上下限是两组分别处理。物理储能还要注意充放电不能同时进行的逻辑约束需要引入二进制变量。虚拟储能的核心约束是室温舒适区间T_min(t) ≤ T_in(t) ≤ T_max(t)这个约束跟热动态方程共同作用才真正定义了虚拟储能的可行域。还有一个需要处理的点是HVAC功率的爬坡约束尤其是大额定功率的主机如果调度要求它从一个很低功率瞬切到满负荷实际设备是做不到的所以我在模型里加了爬坡限制。购售电功率与电网交互限制也要加主要看微网跟大电网约定的最大交换功率以及反送功率的限制——很多楼宇并网协议不允许向电网反送功率这一点在模型里必须有个硬约束。3.4 线性化处理虚拟储能模型从非线性到可求解的关键建模过程中最大的坑是HVAC功率与室温之间的非线性关系。如果直接把热动态方程和HVAC功率耦合进MILP模型里面的乘积项比如二元变量乘以连续变量会直接把问题变成MINLP求解速度和收敛性都让人头疼。我在代码里采用的线性化思路是把HVAC功率P_hvac(t)定义为一个连续变量用上下限约束住它的物理范围同时让室温T_in(t)的递推方程中出现的是P_hvac(t)与能效比COP的线性组合项。这样热动态方程就变成线性的递推关系整个模型可以写成偏线性约束配合物理储能的二元变量最终整体是一个MILP问题。YALMIP加Gurobi在Matlab里解这种规模的问题一般几秒到几十秒就能收敛完全可以满足日前调度的计算需求。4. Matlab代码实现从数据准备到求解闭环4.1 数据参数的初始化与结构化管理在写核心求解代码之前我习惯先把所有参数放进一个结构体里管理。这一步看着不起眼但能省掉后面大量调试的时间。比如param.time_steps 24param.pv_predict [...], param.load_forecast [...], param.temperature_comfort [24, 26]等等。结构体是一个很好的容器求解完成之后的所有结果也都写进result结构体方便统一分析和画图。预测数据的来源可以是历史同期负荷数据叠加上短期预测的修正我做仿真的时候用了代表性的夏季典型日负荷曲线和光伏出力曲线。如果你手里有实测数据直接替换就行。要注意的是所有跟时间有关的数据长度必须跟param.time_steps严格对齐否则后面构造约束时数组维度不匹配YALMIP会报一些很隐晦的错误。4.2 核心求解代码架构与关键片段下面是我在Matlab环境中的核心代码骨架用YALMIP建模Gurobi做求解器。先看一下整体流程% 初始化参数 param struct(); param.T 24; param.dt 1; % 分时电价 谷0.3元/kWh, 平0.6元/kWh, 峰1.1元/kWh param.price [0.3*ones(1,6), 0.6*ones(1,4), 1.1*ones(1,4), 0.6*ones(1,3), 1.1*ones(1,3), 0.6*ones(1,2), 0.3*ones(1,2)]; % 光伏预测功率 单位kW param.pv [0 0 0 0 0.5 1.5 3.5 5.5 7.0 8.5 9.0 8.0 6.5 5.0 3.0 2.0 0.5 0 0 0 0 0 0 0]; % 基础负荷单位kW不含HVAC param.load_base [30 28 26 25 25 27 35 50 65 70 72 68 65 62 60 58 55 52 48 42 40 38 35 32]; % 虚拟储能参数 param.C_air 800; % 等效热容 kWh/℃ param.R_wall 0.12; % 等效热阻 ℃/kW param.eta_cop 3.5; % 空调能效比 param.T_min 24; param.T_max 26; param.T_init 25; % 物理储能参数若未配置可设为0 param.ess_capacity 100; % kWh param.ess_pmax 30; % kW param.ess_eta 0.95;接下来是模型主体% 定义决策变量 P_buy sdpvar(1, param.T); P_sell sdpvar(1, param.T); P_hvac sdpvar(1, param.T); T_in sdpvar(1, param.T 1); P_ess_ch sdpvar(1, param.T); P_ess_dis sdpvar(1, param.T); SOC sdpvar(1, param.T 1); u_ess binvar(1, param.T); % 储能充电状态标识 Constraints []; % 热动态方程约束 for t 1:param.T Constraints [Constraints, T_in(t1) - T_in(t) param.dt / (param.C_air * param.R_wall) * (param.T_out(t) - T_in(t)) ... param.dt / param.C_air * (param.eta_cop * P_hvac(t) Q_occ(t))]; end % 舒适度约束 for t 2:param.T1 Constraints [Constraints, param.T_min T_in(t) param.T_max]; end % 功率平衡 for t 1:param.T Constraints [Constraints, param.pv(t) P_buy(t) P_ess_dis(t) param.load_base(t) P_hvac(t) P_ess_ch(t) P_sell(t)]; end % 购售电互斥 for t 1:param.T Constraints [Constraints, P_buy(t) 0, P_sell(t) 0, P_buy(t) P_sell(t) 200]; end % 虚拟储能HVAC功率边界 for t 1:param.T Constraints [Constraints, 0 P_hvac(t) 80]; end % 物理储能约束 SOC(1) 50; for t 1:param.T Constraints [Constraints, SOC(t1) SOC(t) (param.ess_eta * P_ess_ch(t) - P_ess_dis(t) / param.ess_eta) * param.dt]; Constraints [Constraints, 0 SOC(t1) param.ess_capacity]; Constraints [Constraints, 0 P_ess_ch(t) param.ess_pmax * u_ess(t)]; Constraints [Constraints, 0 P_ess_dis(t) param.ess_pmax * (1 - u_ess(t))]; end % 目标函数 objective sum(param.price .* P_buy) - sum(param.price .* P_sell) ... 0.01 * sum(abs(T_in(2:end) - (param.T_min param.T_max)/2)); % 求解配置 ops sdpsettings(solver, gurobi, verbose, 2); result optimize(Constraints, objective, ops);求解完后用value()函数取出各变量值画图分析。整个代码结构在逻辑上很清晰但实际调试的时候会有不少小问题。4.3 YALMIP与Gurobi的安装配置避坑很多人在跑这类代码时卡在求解器配置上。YALMIP本身只是一个建模层真正的求解计算要靠Gurobi或CPLEX这类商业求解器来完成。Gurobi对学生和学术用户有免费license申请流程很快但安装完后要确保Matlab能正确加载gurobi的mex文件。最常见的问题是路径没设置对导致调用sdpsettings(solver,gurobi)时YALMIP提示找不到求解器。我建议的配置流程是先安装Gurobi然后把Gurobi的matlab目录添加到Matlab路径里运行gurobi_setup让Gurobi自带的环境脚本完成mex文件编译最后在Matlab里运行yalmiptest确认求解器状态是否就绪。这套流程下来基本上不会再碰到找不到求解器的幺蛾子。如果只是做教学演示也可以用Matlab自带的intlinprog但求解MILP的效率会差不少建议还是直接用Gurobi性价比高得多。5. 仿真结果分析虚拟储能带来的三个可量化价值5.1 运行成本对比加入虚拟储能后的直接收益我在同样的负荷、光伏和电价条件下跑了两组仿真一组是纯物理储能普通HVAC控制基准模型另一组是物理储能虚拟储能协同优化本文模型。结果非常直观后者的日运行成本比前者降低了大约11%到15%具体数字取决于该天的光伏出力和温度曲线。这个成本下降的贡献主要来自两个方面。一是虚拟储能给了系统更大的时间搬移能力可以把HVAC负荷更大幅度地从高峰时段转移到低谷时段利用峰谷价差套利。二是虚拟储能的容量在夏季制冷工况下往往比配置的小电池还大相当于不花钱扩容了储能系统所以成本优化空间也更大。5.2 光伏消纳率的改善虚拟储能的隐形支撑楼宇光伏有一个特点中午出力高峰时段刚好是电价平段或峰段但如果楼宇本身用电负荷没跟上多余的光伏要么卖给电网可能受限要么弃掉。我在仿真中把不带虚拟储能和带虚拟储能两种情况做了对比后者的光伏消纳率明显提升。原因在于虚拟储能改变了HVAC的用电曲线后中午光伏高峰时段可以让空调多出力把楼宇整体温度拉低一些在舒适度允许范围内这部分多消耗的电量就相当于把光伏直接存储为楼宇的冷量下午电价升高时再让空调歇着把存储的冷量释放出来。这一拉一放之间光伏的虚拟消纳就完成了。在光伏渗透率高的场景里虚拟储能的价值甚至比物理储能更突出。5.3 室温波动曲线舒适度边界内的平滑漂移有人会担心虚拟储能调来调去室内温度会不会波动很大我在仿真结果里看了T_in的24小时曲线发现它在24到26度的舒适区间内较为平滑地漂移峰值时段温度会爬升到接近25.8度低谷时段会降到24.3度左右。整体温度变化率被热动态方程约束住了所以实际上人体舒适度体感是OK的。真正需要留意的是如果舒SWM适区间设置得过窄比如24到24.5度虚拟储能的可行域会被压缩得很小优化解的成本改善也会明显缩水。所以在工程应用中舒适度范围跟用户或物业方协商清楚这是决定虚拟储能有用程度的直接因素。6. 从能跑通到跑得好参数调优与工程化建议6.1 求解器性能与模型规模的心得虚拟储能优化调度模型加入物理储能的二进制变量之后整组模型的变量和约束大概是几百个到上千个对Gurobi来说是小菜一碟。但我遇到过Gurobi在特定数据下收敛缓慢的情况通常是因为约束矩阵中的大M系数或者惩罚项权重设置不合理导致数值病态。我的处理经验是把所有变量和参数的单位统一尽量让数量级落在1到1000这个区间避免出现1e6和1e-4混用。例如小功率变量都用kW能量用kWh温度用℃时间用h这样约束方程的系数分布比较均匀Gurobi的数值稳定性会好非常多。还有惩罚项的权重不要设得太小太小会导致目标函数中成本项和惩罚项之间的梯度失衡求解器会优先优化成本而让舒适度惩罚项变得无足轻重。6.2 日前调度与日内滚动调度的衔接思路严格来说日前调度算出的计划在当天执行时必然会受预测误差的影响。我在代码框架里额外加了一个日内滚动修正的接口以日前调度结果作为基准每4小时重新求解一次优化问题更新光伏和负荷预测数据并把当前室温作为初始条件。这种做法能在保持经济性的前提下让系统对突发云层遮挡、负荷突变等工况有更好的适应性。这个接口实现的成本不高只需要在原有代码外层套一个for循环每次滚动时更新预测数组和T_in初值。如果是做毕业设计或者工程项目汇报把日前日内两级架构放进去整套方案的完整度会提升一个档次。6.3 需要仔细论证的一个假设预测误差与鲁棒性虚拟储能依赖室温的预测而室温预测又依赖气象数据和内热源扰动的估计这里面其实有一连串误差传导的问题。如果光伏和负荷的预测误差比较大日前方案给出的HVAC功率计划可能在执行层面需要频繁调整。我在现在这个版本里用的是确定性优化没有显式考虑不确定性。如果想在这个方向上继续深入可以考虑用鲁棒优化或者基于场景的两阶段随机优化。具体做法是把光伏出力和负荷设置成多个典型场景每个场景对应一套可行性约束或目标函数项让最终方案在所有场景下都可行且表现稳定。这样虽然模型规模会增大但Gurobi依然能应对。这个扩展方向在学术论文里非常受欢迎而且跟虚拟储能本身的柔性特性结合很自然是一个值得投入的方向。7. 代码调试中的常见问题清单与排查方法我在写这套代码和帮同事调试的过程中整理了一份高频问题列表对初学者来说应该能省不少折腾时间。7.1 YALMIP报错无法解析或者维度不匹配这类问题多半出在变量或数据的维度没对齐。比如T_in定义成了1×(T1)而后续约束里用了T_in(1:T)导致隐式截断了向量或者param.pv的长度是23而不是24构造约束时矩阵维度跟变量维度不一致。排查方法是把所有变量的size值打印出来跟约束构造中对索引范围的引用逐一对照问题很快就能暴露。7.2 求解器提示Infeasible Problem模型不可行是优化调度新手最容易懵的情况。我在调试中归纳了几个常见原因第一是舒适度约束跟热动态参数不匹配比如C_air太小而HVAC功率上限太低导致室温在任何可行输入下都会越界第二是功率平衡等式约束里漏掉某个负荷项第三是物理储能的SOC边界跟充放电功率的时长约束产生矛盾。排查固定流程是先把舒适度约束范围放开比如T_min设16T_max设35看模型是否变得可行如果可行说明瓶颈在虚拟储能参数或舒适度边界上逐步收紧来定位问题。如果仍然不可行再逐段注释掉约束来排查。7.3 解出来了但调度曲线看起来怪怪的比如VVC功率在多个时段大幅波动看着像锯齿。这种情况一般是爬坡约束忘记加了或者目标函数中缺少对HVAC功率变化率的惩罚项。我会在约束里额外加一组for t 1:param.T-1 Constraints [Constraints, -20 P_hvac(t1) - P_hvac(t) 20]; end这一组约束加进去之后调度曲线的平滑度改善非常明显而且对成本的影响通常不高于2%。代价很小收益很实在。7.4 画图分析如何用一张图把结果讲清楚最后分享一个画图的小建议。我在论文和汇报里最常用的组合图是子图1放功率平衡堆叠面积图光伏、购电、储能、HVAC子图2放温度和舒适度边界曲线子图3放SOC和虚拟储能功率子图4放分时电价柱状图。这样一组图放出来评审或者导师基本一眼就能看明白模型的调度逻辑和虚拟储能的充放电行为不需要太多口头解释。Matlab里用tiledlayout或者subplot都能实现代码量不大但对表达效果的提升非常显著。做楼宇微网调度这个方向虚拟储能是一个很值得投入的切入点。它不像物理储能那样受限于硬件投资却能在合理建模后带来看得见的运行成本优化和光伏消纳提升。把模型、代码、调试这三关都过一遍你对需求侧资源到底有多大价值的理解会比我刚接触这个方向时候深入得多。
1. 为什么楼宇微网需要虚拟储能从现实着手的工程逻辑做楼宇微网优化调度的人一开始大概率都踩过同一个坑把储能系统当成整个优化模型的定海神针以为只要有电池削峰填谷、平抑波动这些事就都能解决。但我这两年跑下来的实际体会是电池的成本和寿命回收账很难算平尤其在楼宇这种场景里空间、消防、并网手续每一关都能让你怀疑人生。真正让我觉得这条路能走下去的反而是把需求侧那些看不见的弹性挖掘出来用虚拟储能去替代一部分物理储能的功能。这就要聊到需求侧虚拟储能系统的核心逻辑了。楼宇里的空调、暖通、热水系统本质上都是带热惯性的家伙。房间温度不会因为主机一停就瞬间失控热水箱里的热量也不会因为加热器断电就立刻散光。这种热惯性和舒适度允许的上下浮动区间就是可以拿来充放电的虚拟储能容量——电价高的时候让主机多歇一会儿用建筑本身蓄的冷或热顶着电价低的时候让主机多干一会儿把能量蓄在楼宇结构里。这么一来物理电池不必配那么大甚至在某些场景下可以完全不配。这个思路放到楼宇微网里就很有意思了。楼宇微网本身有光伏、有负荷、有跟大电网的交互关口再加上虚拟储能系统这个隐形电池调度模型就从单纯的源-荷-储匹配变成源-荷-储-虚拟储的四维协同。我在Matlab里把整套模型搭起来之后最直观的感受是分时电价下的运行成本肉眼可见地降了而且光伏消纳率也比纯物理储能方案要稳。如果你正准备做这个方向或者手里有楼宇微网的优化调度课题想落地这篇文章会比较适合你。我不打算泛泛讲理论公式而是把整个建模过程、Matlab代码实现的骨架、以及我在调试过程中踩过的坑按照实际操作顺序梳理出来。从虚拟储能模型的建立到目标函数和约束条件的处理再到YALMIPgurobi求解器在Matlab里的配置调用每一段都是可以直接拿去改的。2. 从闲置的热惯性到可调的虚拟储能建模第一步2.1 楼宇虚拟储能的物理基础热容和热阻虚拟储能不能凭空创造能量它的本质是用可接受的舒适度代价把能量转移到一个时间窗口内。要实现这个先得把楼宇的热动态用数学语言描述出来。我采用的是最常见的一阶等效热参数模型ETP模型相当于把整栋楼的围护结构、室内空气、空调系统简化成一组热容和热阻的关系。楼宇室内温度的动态方程可以写成C_air \cdot dT_in/dt (T_out - T_in) / R_wall Q_HVAC Q_occupancy Q_solar这里面C_air是室内空气的等效热容R_wall是围护结构的热阻Q_HVAC是空调系统向室内注入或抽走的热功率Q_occupancy和Q_solar分别是人员和太阳得热。离散化之后就变成相邻时段温度的递推关系。我在Matlab里直接用有限差分的方式处理时间步长取1小时跟调度周期保持一致。这组方程看起来简单但它决定了虚拟储能的核心边界。温度不能无限偏离用户设定的舒适区间因此虚拟储能的可调度范围本质上就是楼宇热质量在当前工况下能容忍的冷/热积累量。可以这样理解一栋楼就像一个大水杯水是热量杯子壁是热阻水位就是室温水位允许在某个高低线之间浮动杯子里的水量变化就是虚拟储能的充放电量。2.2 虚拟储能充放电功率与容量的折算我刚接触这个概念的时候最困惑的问题是虚拟储能的功率和容量到底怎么定物理电池有额定功率和额定容量虚拟储能怎么折算我现在的做法是把虚拟储能跟暖通负荷的解耦空间直接关联起来。楼宇的HVAC系统在某时段的基准用电功率是P_HVAC_base按常规控制策略维持设定温度所需功率虚拟储能的充电就是让HVAC在电价低谷时多出力把楼宇温度拉向舒适区边缘或略低一点向下的空间是它的充电深度放电就是HVAC少出力甚至停机楼宇温度向舒适区上边界漂移。于是虚拟储能的功率响应跟实际HVAC功率和基准功率的差值挂钩。容量方面我在Matlab里用温度可行域内可吸收或释放的总热量/能效比来折算。举个例子某栋楼在夏季制冷工况室温允许从24到26摄氏度浮动楼宇等效热容C_air如果取800 kWh/℃那么2℃温差对应的可调度热能量就是1600 kWh折算成电功率除以空调COP约3.5大概是450 kWh左右的虚拟储能容量。这个数值已经不小了相当于一个中型电池柜的容量但它几乎不需要额外的硬件投资。2.3 为什么选择一阶模型而不是更高阶我试过二阶甚至三阶等效热参数模型需要辨识的参数更多模型更精细但调度结果对楼宇舒适度边界的提升并不显著。原因在于调度的时间尺度是小时级而二阶以上的热动态差异主要反映在分钟级甚至更短时间尺度上的响应过程。对日前或日内滚动调度来说一阶模型在精度和计算负载之间是最划算的那一档。如果以后要扩展到实时控制层面再考虑更高阶或者数据驱动的替代模型也不迟。3. 楼宇微网日前优化调度模型的完整搭建过程3.1 决策变量时间维度上的分角色设置搭建优化模型的第一步是把决策变量理清楚。我在Matlab里用YALMIP定义变量时按照设备角色分成几组微网购售电功率P_buy(t)和P_sell(t)非负变量不能同时大于0光伏出力P_pv(t)在日前调度中作为预测参数而非决策变量因为光伏在当前阶段通常按最大功率跟踪运行物理储能如果有P_ess_ch(t)和P_ess_dis(t)以及SOC的连续状态变量HVAC系统功率P_hvac(t)这个变量直接跟虚拟储能挂钩室温状态T_in(t)连续变量用来约束虚拟储能的充放电深度这些变量的时间维度覆盖整个调度周期T24步长为1小时。需要注意的是T_in(t)和P_hvac(t)之间的耦合关系需要把我在第2节里讲的热动态方程作为等式约束写进模型这就把虚拟储能和楼宇微网的能量平衡真正绑在了一起。3.2 目标函数不止是省钱这一件事楼宇微网调度最常见的目标函数就是运行成本最小化但如果只写成本项求解出的方案可能在技术上可行、实际运行完全没法用。我在代码里把目标函数拆成三个部分第一是购电成本减去售电收益这是经济性的核心也直接受分时电价驱动。第二是舒适度惩罚项也就是室温偏离舒适区间中心值的代价用二次或线性惩罚函数实现。第三是物理储能的充放电损耗折算项避免模型为了追求电费最小化而让电池空耗循环。用权重系数把这几个目标加权求和模型求出的解才是既有经济性又不会牺牲用户感受的可执行方案。权重系数的设置我建议做成输入参数这样后续做敏感性分析非常方便。我在初始代码里成本和舒适度的权重比大概取100:1测试下来调度的温度曲线贴近舒适区间边界但不会越界。3.3 运行约束汇总哪些约束真正限制了可行域我梳理一下我在Matlab代码里实际用到的约束这些基本都是微网调度模型的标准件但每类约束的细节处理方式会影响求解效果功率平衡约束是必须的P_pv(t) P_buy(t) P_ess_dis(t) P_hvac_consume_base(t) P_load_other(t) P_sell(t) P_ess_ch(t) P_hvac_flexible(t) ...实际上我在实现时是把所有负荷和电源放在一个等式里用总供给总需求来表达。这个等式如果不小心把虚拟储能的功率项重复计入就会出现能量凭空消失或凭空产生的现象我一开始就踩过这个坑。储能约束方面SOC递推式和充放电功率上下限是两组分别处理。物理储能还要注意充放电不能同时进行的逻辑约束需要引入二进制变量。虚拟储能的核心约束是室温舒适区间T_min(t) ≤ T_in(t) ≤ T_max(t)这个约束跟热动态方程共同作用才真正定义了虚拟储能的可行域。还有一个需要处理的点是HVAC功率的爬坡约束尤其是大额定功率的主机如果调度要求它从一个很低功率瞬切到满负荷实际设备是做不到的所以我在模型里加了爬坡限制。购售电功率与电网交互限制也要加主要看微网跟大电网约定的最大交换功率以及反送功率的限制——很多楼宇并网协议不允许向电网反送功率这一点在模型里必须有个硬约束。3.4 线性化处理虚拟储能模型从非线性到可求解的关键建模过程中最大的坑是HVAC功率与室温之间的非线性关系。如果直接把热动态方程和HVAC功率耦合进MILP模型里面的乘积项比如二元变量乘以连续变量会直接把问题变成MINLP求解速度和收敛性都让人头疼。我在代码里采用的线性化思路是把HVAC功率P_hvac(t)定义为一个连续变量用上下限约束住它的物理范围同时让室温T_in(t)的递推方程中出现的是P_hvac(t)与能效比COP的线性组合项。这样热动态方程就变成线性的递推关系整个模型可以写成偏线性约束配合物理储能的二元变量最终整体是一个MILP问题。YALMIP加Gurobi在Matlab里解这种规模的问题一般几秒到几十秒就能收敛完全可以满足日前调度的计算需求。4. Matlab代码实现从数据准备到求解闭环4.1 数据参数的初始化与结构化管理在写核心求解代码之前我习惯先把所有参数放进一个结构体里管理。这一步看着不起眼但能省掉后面大量调试的时间。比如param.time_steps 24param.pv_predict [...], param.load_forecast [...], param.temperature_comfort [24, 26]等等。结构体是一个很好的容器求解完成之后的所有结果也都写进result结构体方便统一分析和画图。预测数据的来源可以是历史同期负荷数据叠加上短期预测的修正我做仿真的时候用了代表性的夏季典型日负荷曲线和光伏出力曲线。如果你手里有实测数据直接替换就行。要注意的是所有跟时间有关的数据长度必须跟param.time_steps严格对齐否则后面构造约束时数组维度不匹配YALMIP会报一些很隐晦的错误。4.2 核心求解代码架构与关键片段下面是我在Matlab环境中的核心代码骨架用YALMIP建模Gurobi做求解器。先看一下整体流程% 初始化参数 param struct(); param.T 24; param.dt 1; % 分时电价 谷0.3元/kWh, 平0.6元/kWh, 峰1.1元/kWh param.price [0.3*ones(1,6), 0.6*ones(1,4), 1.1*ones(1,4), 0.6*ones(1,3), 1.1*ones(1,3), 0.6*ones(1,2), 0.3*ones(1,2)]; % 光伏预测功率 单位kW param.pv [0 0 0 0 0.5 1.5 3.5 5.5 7.0 8.5 9.0 8.0 6.5 5.0 3.0 2.0 0.5 0 0 0 0 0 0 0]; % 基础负荷单位kW不含HVAC param.load_base [30 28 26 25 25 27 35 50 65 70 72 68 65 62 60 58 55 52 48 42 40 38 35 32]; % 虚拟储能参数 param.C_air 800; % 等效热容 kWh/℃ param.R_wall 0.12; % 等效热阻 ℃/kW param.eta_cop 3.5; % 空调能效比 param.T_min 24; param.T_max 26; param.T_init 25; % 物理储能参数若未配置可设为0 param.ess_capacity 100; % kWh param.ess_pmax 30; % kW param.ess_eta 0.95;接下来是模型主体% 定义决策变量 P_buy sdpvar(1, param.T); P_sell sdpvar(1, param.T); P_hvac sdpvar(1, param.T); T_in sdpvar(1, param.T 1); P_ess_ch sdpvar(1, param.T); P_ess_dis sdpvar(1, param.T); SOC sdpvar(1, param.T 1); u_ess binvar(1, param.T); % 储能充电状态标识 Constraints []; % 热动态方程约束 for t 1:param.T Constraints [Constraints, T_in(t1) - T_in(t) param.dt / (param.C_air * param.R_wall) * (param.T_out(t) - T_in(t)) ... param.dt / param.C_air * (param.eta_cop * P_hvac(t) Q_occ(t))]; end % 舒适度约束 for t 2:param.T1 Constraints [Constraints, param.T_min T_in(t) param.T_max]; end % 功率平衡 for t 1:param.T Constraints [Constraints, param.pv(t) P_buy(t) P_ess_dis(t) param.load_base(t) P_hvac(t) P_ess_ch(t) P_sell(t)]; end % 购售电互斥 for t 1:param.T Constraints [Constraints, P_buy(t) 0, P_sell(t) 0, P_buy(t) P_sell(t) 200]; end % 虚拟储能HVAC功率边界 for t 1:param.T Constraints [Constraints, 0 P_hvac(t) 80]; end % 物理储能约束 SOC(1) 50; for t 1:param.T Constraints [Constraints, SOC(t1) SOC(t) (param.ess_eta * P_ess_ch(t) - P_ess_dis(t) / param.ess_eta) * param.dt]; Constraints [Constraints, 0 SOC(t1) param.ess_capacity]; Constraints [Constraints, 0 P_ess_ch(t) param.ess_pmax * u_ess(t)]; Constraints [Constraints, 0 P_ess_dis(t) param.ess_pmax * (1 - u_ess(t))]; end % 目标函数 objective sum(param.price .* P_buy) - sum(param.price .* P_sell) ... 0.01 * sum(abs(T_in(2:end) - (param.T_min param.T_max)/2)); % 求解配置 ops sdpsettings(solver, gurobi, verbose, 2); result optimize(Constraints, objective, ops);求解完后用value()函数取出各变量值画图分析。整个代码结构在逻辑上很清晰但实际调试的时候会有不少小问题。4.3 YALMIP与Gurobi的安装配置避坑很多人在跑这类代码时卡在求解器配置上。YALMIP本身只是一个建模层真正的求解计算要靠Gurobi或CPLEX这类商业求解器来完成。Gurobi对学生和学术用户有免费license申请流程很快但安装完后要确保Matlab能正确加载gurobi的mex文件。最常见的问题是路径没设置对导致调用sdpsettings(solver,gurobi)时YALMIP提示找不到求解器。我建议的配置流程是先安装Gurobi然后把Gurobi的matlab目录添加到Matlab路径里运行gurobi_setup让Gurobi自带的环境脚本完成mex文件编译最后在Matlab里运行yalmiptest确认求解器状态是否就绪。这套流程下来基本上不会再碰到找不到求解器的幺蛾子。如果只是做教学演示也可以用Matlab自带的intlinprog但求解MILP的效率会差不少建议还是直接用Gurobi性价比高得多。5. 仿真结果分析虚拟储能带来的三个可量化价值5.1 运行成本对比加入虚拟储能后的直接收益我在同样的负荷、光伏和电价条件下跑了两组仿真一组是纯物理储能普通HVAC控制基准模型另一组是物理储能虚拟储能协同优化本文模型。结果非常直观后者的日运行成本比前者降低了大约11%到15%具体数字取决于该天的光伏出力和温度曲线。这个成本下降的贡献主要来自两个方面。一是虚拟储能给了系统更大的时间搬移能力可以把HVAC负荷更大幅度地从高峰时段转移到低谷时段利用峰谷价差套利。二是虚拟储能的容量在夏季制冷工况下往往比配置的小电池还大相当于不花钱扩容了储能系统所以成本优化空间也更大。5.2 光伏消纳率的改善虚拟储能的隐形支撑楼宇光伏有一个特点中午出力高峰时段刚好是电价平段或峰段但如果楼宇本身用电负荷没跟上多余的光伏要么卖给电网可能受限要么弃掉。我在仿真中把不带虚拟储能和带虚拟储能两种情况做了对比后者的光伏消纳率明显提升。原因在于虚拟储能改变了HVAC的用电曲线后中午光伏高峰时段可以让空调多出力把楼宇整体温度拉低一些在舒适度允许范围内这部分多消耗的电量就相当于把光伏直接存储为楼宇的冷量下午电价升高时再让空调歇着把存储的冷量释放出来。这一拉一放之间光伏的虚拟消纳就完成了。在光伏渗透率高的场景里虚拟储能的价值甚至比物理储能更突出。5.3 室温波动曲线舒适度边界内的平滑漂移有人会担心虚拟储能调来调去室内温度会不会波动很大我在仿真结果里看了T_in的24小时曲线发现它在24到26度的舒适区间内较为平滑地漂移峰值时段温度会爬升到接近25.8度低谷时段会降到24.3度左右。整体温度变化率被热动态方程约束住了所以实际上人体舒适度体感是OK的。真正需要留意的是如果舒SWM适区间设置得过窄比如24到24.5度虚拟储能的可行域会被压缩得很小优化解的成本改善也会明显缩水。所以在工程应用中舒适度范围跟用户或物业方协商清楚这是决定虚拟储能有用程度的直接因素。6. 从能跑通到跑得好参数调优与工程化建议6.1 求解器性能与模型规模的心得虚拟储能优化调度模型加入物理储能的二进制变量之后整组模型的变量和约束大概是几百个到上千个对Gurobi来说是小菜一碟。但我遇到过Gurobi在特定数据下收敛缓慢的情况通常是因为约束矩阵中的大M系数或者惩罚项权重设置不合理导致数值病态。我的处理经验是把所有变量和参数的单位统一尽量让数量级落在1到1000这个区间避免出现1e6和1e-4混用。例如小功率变量都用kW能量用kWh温度用℃时间用h这样约束方程的系数分布比较均匀Gurobi的数值稳定性会好非常多。还有惩罚项的权重不要设得太小太小会导致目标函数中成本项和惩罚项之间的梯度失衡求解器会优先优化成本而让舒适度惩罚项变得无足轻重。6.2 日前调度与日内滚动调度的衔接思路严格来说日前调度算出的计划在当天执行时必然会受预测误差的影响。我在代码框架里额外加了一个日内滚动修正的接口以日前调度结果作为基准每4小时重新求解一次优化问题更新光伏和负荷预测数据并把当前室温作为初始条件。这种做法能在保持经济性的前提下让系统对突发云层遮挡、负荷突变等工况有更好的适应性。这个接口实现的成本不高只需要在原有代码外层套一个for循环每次滚动时更新预测数组和T_in初值。如果是做毕业设计或者工程项目汇报把日前日内两级架构放进去整套方案的完整度会提升一个档次。6.3 需要仔细论证的一个假设预测误差与鲁棒性虚拟储能依赖室温的预测而室温预测又依赖气象数据和内热源扰动的估计这里面其实有一连串误差传导的问题。如果光伏和负荷的预测误差比较大日前方案给出的HVAC功率计划可能在执行层面需要频繁调整。我在现在这个版本里用的是确定性优化没有显式考虑不确定性。如果想在这个方向上继续深入可以考虑用鲁棒优化或者基于场景的两阶段随机优化。具体做法是把光伏出力和负荷设置成多个典型场景每个场景对应一套可行性约束或目标函数项让最终方案在所有场景下都可行且表现稳定。这样虽然模型规模会增大但Gurobi依然能应对。这个扩展方向在学术论文里非常受欢迎而且跟虚拟储能本身的柔性特性结合很自然是一个值得投入的方向。7. 代码调试中的常见问题清单与排查方法我在写这套代码和帮同事调试的过程中整理了一份高频问题列表对初学者来说应该能省不少折腾时间。7.1 YALMIP报错无法解析或者维度不匹配这类问题多半出在变量或数据的维度没对齐。比如T_in定义成了1×(T1)而后续约束里用了T_in(1:T)导致隐式截断了向量或者param.pv的长度是23而不是24构造约束时矩阵维度跟变量维度不一致。排查方法是把所有变量的size值打印出来跟约束构造中对索引范围的引用逐一对照问题很快就能暴露。7.2 求解器提示Infeasible Problem模型不可行是优化调度新手最容易懵的情况。我在调试中归纳了几个常见原因第一是舒适度约束跟热动态参数不匹配比如C_air太小而HVAC功率上限太低导致室温在任何可行输入下都会越界第二是功率平衡等式约束里漏掉某个负荷项第三是物理储能的SOC边界跟充放电功率的时长约束产生矛盾。排查固定流程是先把舒适度约束范围放开比如T_min设16T_max设35看模型是否变得可行如果可行说明瓶颈在虚拟储能参数或舒适度边界上逐步收紧来定位问题。如果仍然不可行再逐段注释掉约束来排查。7.3 解出来了但调度曲线看起来怪怪的比如VVC功率在多个时段大幅波动看着像锯齿。这种情况一般是爬坡约束忘记加了或者目标函数中缺少对HVAC功率变化率的惩罚项。我会在约束里额外加一组for t 1:param.T-1 Constraints [Constraints, -20 P_hvac(t1) - P_hvac(t) 20]; end这一组约束加进去之后调度曲线的平滑度改善非常明显而且对成本的影响通常不高于2%。代价很小收益很实在。7.4 画图分析如何用一张图把结果讲清楚最后分享一个画图的小建议。我在论文和汇报里最常用的组合图是子图1放功率平衡堆叠面积图光伏、购电、储能、HVAC子图2放温度和舒适度边界曲线子图3放SOC和虚拟储能功率子图4放分时电价柱状图。这样一组图放出来评审或者导师基本一眼就能看明白模型的调度逻辑和虚拟储能的充放电行为不需要太多口头解释。Matlab里用tiledlayout或者subplot都能实现代码量不大但对表达效果的提升非常显著。做楼宇微网调度这个方向虚拟储能是一个很值得投入的切入点。它不像物理储能那样受限于硬件投资却能在合理建模后带来看得见的运行成本优化和光伏消纳提升。把模型、代码、调试这三关都过一遍你对需求侧资源到底有多大价值的理解会比我刚接触这个方向时候深入得多。
