这几年独立储能项目密集上马但真正把商业模式跑通的并不多。大家盯得最紧的一件事就是让储能同时在现货电能量市场和调频辅助服务市场里拿到两笔收入。可问题在于两个市场如果分开出清储能同一时刻的容量会被两套计划重复占用算出来的收益和实际运行完全对不上。所以核心不是“怎么分别报价”而是“怎么把两个市场的出清放在一个模型里协调”。这篇文章把独立储能在现货电能量与调频辅助服务市场中的协调出清机制从数学建模到Matlab实现完整过一遍代码、参数、坑点都直接给出来。适合电力市场方向的研究生、设计院工程师还有要做储能收益测算的投资分析人员。1. 独立储能参与市场的底层逻辑1.1 为什么两个市场必须协调出清理解协调出清先得回到储能的物理本质上。独立储能就是一组电池加一台PCS电池容量和PCS功率是同一份资源。电能量市场要的是有功出力调频辅助服务市场要的也是有功响应能力两者在本质上争夺的是同一个功率容量。如果电能量市场先安排储能低谷充电、高峰放电调频市场又假设储能随时可以上调或下调功率两套计划叠加后很可能超过PCS额定功率甚至超过电池允许的充放电倍率。我做仿真时最早也犯过这个错误电能量出清结果显示储能某个时段放电40MW调频出清又要求它保留20MW的响应能力结果PCS总共需要60MW的能力而实际额定只有50MW。这种结果看着有收益落地就是废纸。所以要协调核心目的就是不让同一份容量被两个市场重复使用让每个时段的充放电功率和调频容量加起来不超过储能的最大功率能力。从市场信号的角度看也必须有协调机制。储能参与能量市场时报价里隐含的是它放弃调频机会的机会成本参与调频市场时报价里隐含的是它减少能量套利的机会成本。如果两个市场分别算谁都看不到对方的价格信号储能报出来的价格就是失真的。协调出清本质上把两个市场的价格信号放在同一个优化目标里让容量自动流向收益更高的用途。1.2 协调出清与顺序出清的差异很多研究论文喜欢拿顺序出清当对比基准。顺序出清的做法是先跑电能量市场得到储能各时段的充放电计划再把剩余容量拿去参与调频市场。看起来简单但问题很典型。比如低谷时段顺序出清会让储能以满功率充电结果剩余调频容量为零但调频市场给出的容量价格很高本来储能可以少充一点电留出一部分容量去赚调频费整体收益反而更高。协调出清就不会干这种蠢事。它把能量报价和调频报价一起丢进优化器由模型决定每个时段储能到底充多少、放多少、预留多少调频容量目标是系统总成本最小或储能综合收益最大。用一个不恰当的类比同一辆司机既跑网约车又跑同城快递如果两个平台各派各的单司机同一小时可能被派两单根本接不过来协调调度就是让司机自己根据两边的单价决定接哪边的活。顺序出清和协调出清在结果上的差异我总结成一张对比表对比项顺序出清协调出清容量分配逻辑先能量后调频剩余容量参与调频统一优化同时决定能量功率和调频容量是否考虑机会成本不考虑通过联合目标函数自动体现储能收益通常偏低容量利用不充分更高且容量分配更合理系统总成本偏高可能次优理论上全局最优模型复杂度低求解快高需要处理耦合约束当然顺序出清在计算效率上有优势真实市场中也存在这种模式。但如果研究独立储能的商业模式我建议直接上协调出清因为这样算出来的收益才是储能真正能摸到的天花板。1.3 模型边界与基本假设任何模型都要先说清楚假设。我这套模型基于的是集中式出清框架面向一个单节点系统忽略网络潮流和网损。独立储能在这个系统里被视为价格接受者也就是它单独改变充放电策略不会明显改变出清价格这对中小规模储能是合理近似。调频辅助服务市场我做了简化只考虑容量补偿不追踪调频里程。原因很简单市场出清阶段主要解决“预留多少调频容量”的问题而里程收益属于运行结算问题两者时间尺度差异很大。如果硬要把里程和容量一起优化模型会变得非常复杂也不利于聚焦协调机制本身。储能参与调频的细节比如AGC信号跟踪后续可以单独建模但在出清模型里只需要容量约束就够了。时间尺度上算例里我用1小时为间隔共24个时段。实际调频需求是秒级的但市场出清往往也是小时级或5分钟级用容量备用作为近似是行业通行做法。如果你要做日内滚动优化可以把时段切成96点或288点模型结构不用变只是计算量会涨。2. 协调出清问题的数学建模2.1 目标函数设计协调出清的目标函数我习惯写成系统总成本最小化而不是储能收益最大化。因为从市场运营者角度看出清的目标应该是社会福利最大或购电成本最小储能收益是结果不是目标。当然如果你做的是储能投资决策可以把目标改为储能净收益但那是另一个问题不是市场出清问题。目标函数包含三部分常规机组发电成本、常规机组和储能提供调频容量的成本、储能的充放电磨损成本。表达式可以写成min Σ_t [ Σ_i (c_i * P_i(t)) Σ_j (r_j * R_j(t)) λ * (P_ch(t) P_dis(t)) ]这里 P_i(t) 是常规机组i在时段t的出力c_i 是它的能量报价R_j(t) 是设备j机组或储能在时段t中标的调频容量r_j 是调频容量报价P_ch 和 P_dis 是储能充电和放电功率λ 是人为设置的磨损成本系数通常很小目的是避免储能频繁改变运行状态。为什么要加储能磨损成本因为如果不加储能可能在电价相同的情况下反复充放出现抖振。虽然功率约束能防止同时充放电但无法避免无意义的来回切换。加一个很小的 λ比如 0.01 元/MWh就能平滑行为对总成本影响可以忽略。调频容量报价的设定很关键。储能如果报一个较低的调频价格比如 5 元/MW系统就会优先让它提供调频如果报 30 元/MW基本只有能量报价特别高的时候才会被选中。这个报价在真实市场中是储能运营商自己决定的在仿真里可以当作输入参数来扫描研究不同价格策略下储能的收益变化。2.2 储能运行约束的细节储能部分的约束是整个模型的核心必须逐条写清楚。第一是功率与容量约束。每个时段的充电功率、放电功率和调频容量不能超过PCS额定功率0 ≤ P_ch(t) ≤ P_max * u_ch(t) 0 ≤ P_dis(t) ≤ P_max * u_dis(t) u_ch(t) u_dis(t) ≤ 1 P_ch(t) R_st(t) ≤ P_max P_dis(t) R_st(t) ≤ P_max前三条约束是防止同时充放电u_ch 和 u_dis 是0-1变量。第四条和第五条是协调出清的精髓储能一旦承诺了 R_st 的调频容量它实际能用于充电或放电的功率就被压缩了。比如P_max是50MW承诺20MW调频那充电功率最多只能30MW放电功率也最多30MW。这背后的物理含义是储能必须保留足够的上下调节区间来响应调频指令。第二是SOC动态约束SOC(t1) SOC(t) (η * P_ch(t) - P_dis(t) / η) * Δt / E_ratedη 是充放电效率E_rated 是额定容量。注意充电和放电效率不是对偶的这里我用的是同一个η充电时乘以η放电时除以η相当于净耗。如果采用更精细的模型可以分别给充电效率η_c和放电效率η_d但那样会引入非线性尽量用常数近似。第三是SOC上下限和初末值约束SOC_min ≤ SOC(t) ≤ SOC_max SOC(1) SOC_initial SOC(T1) SOC_initial首末一致约束用于单日运行分析相当于要求储能一个周期后回到初始电量这样才能保证日复一日运行。如果不加这个约束储能可能把最后几个时段的电量全部放空收益虚高。实战中我一般把SOC_min设为0.1SOC_max设为0.9避免过充过放损害电池。2.3 系统级约束与市场出清条件除了储能自身的约束系统层面还有几条硬约束。首先是功率平衡约束Σ_i P_i(t) P_dis(t) - P_ch(t) D(t)D(t) 是时段t的负荷。这是最基本的电能量平衡所有市场出清模型都有。其次是调频容量需求约束Σ_i Rg_i(t) R_st(t) ≥ R_req(t)R_req(t) 是系统调频容量需求通常取该时段最大负荷的某个比例算例里我用5%。这个约束确保任何时刻都有足够的资源响应调频信号。调频需求也可以按负荷的百分比加一个固定偏差但本质上都是给模型一个外生需求。再次是常规机组的出力上下限和爬坡约束P_min(i) ≤ P_i(t) ≤ P_max(i) -Ramp_i * Δt ≤ P_i(t1) - P_i(t) ≤ Ramp_i * Δt爬坡约束刻画机组的响应速度。如果不加模型会允许火电机组瞬时从最小出力跳到最大出力结果过于乐观。储能虽然没有传统意义上的爬坡限制但PCS存在功率变化率限制一般远大于机组在小时级模型里可以忽略。这些约束合在一起构成了一个混合整数线性规划MILP。0-1变量来自储能同时充放电约束这也是需要商业求解器的主要原因。如果去掉这个约束让储能可以同时充放电问题退化成LP但解在物理上不可行所以不能为了省事丢掉它。3. Matlab代码实现从数据到求解3.1 程序框架与依赖环境Matlab实现我建议用Yalmip建模再调Gurobi或CPLEX求解。Yalmip不是求解器而是一个建模语言最大好处是语法接近数学表达式改求解器只需要改一行配置。安装步骤不多说但要提醒一点Yalmip的路径必须加到Matlab搜索路径里否则会报找不到函数。代码组织上我习惯拆成几个文件case_data.m -- 所有参数和场景数据 build_model.m -- 建立决策变量、约束、目标函数 run_main.m -- 主脚本调数据和建模函数 plot_results.m -- 结果可视化和导出这样做的好处是参数、模型、逻辑分离改一个参数不用翻整个代码。主脚本里按顺序调用就行。3.2 参数设计与数据读入参数设计直接决定仿真结果是否可信。我的case_data.m长这样% 基础参数 T 24; % 时段数 dt 1; % 时间间隔小时 % 常规机组参数: [Pmax(MW), Pmin(MW), Ramp(MW/h), Cost(元/MWh), ReservePrice(元/MW)] unit [ 100, 20, 50, 25, 5; 150, 30, 80, 30, 6 ]; % 独立储能参数: [Pmax(MW), E(MWh), SOCmin, SOCmax, eta, SOC_initial] storage [50, 100, 0.1, 0.9, 0.9, 0.5]; % 负荷曲线 (MW), 24个时段 load_profile [ 45, 42, 40, 38, 40, 48, 60, 75, 85, 92, 95, 90, ... 82, 78, 80, 85, 90, 95, 100, 105, 98, 85, 70, 55 ]; % 调频容量需求取负荷的5% reserve_req 0.05 * load_profile; % 储能充放电磨损成本系数 lambda 0.01;参数为什么这么设机组容量要和负荷峰值匹配两台机组总容量250MW负荷峰值105MW留了足够的调峰空间。储能容量50MW/100MWh大约是负荷峰值的50%这个比例在工程上比较常见再大会出现容量闲置再小又体现不出调节效果。如果要从外部读数据可以用readtable读CSV再转成数值列。硬编码只适合做演示实际研究建议用脚本统一生成。3.3 Yalmip建模与求解核心建模代码写在build_model.m里我贴一段可以跑通的核心代码function [Constraints, Objective, vars] build_model(unit, storage, load_profile, reserve_req) ng size(unit, 1); T length(load_profile); Pmax unit(:, 1); Pmin unit(:, 2); Ramp unit(:, 3); Cost unit(:, 4); ResPrice unit(:, 5); % 决策变量 Pg sdpvar(ng, T, full); % 机组出力 Rg sdpvar(ng, T, full); % 机组调频容量 Pch sdpvar(1, T, full); % 储能充电功率 Pdis sdpvar(1, T, full); % 储能放电功率 Rst sdpvar(1, T, full); % 储能调频容量 SOC sdpvar(1, T1, full); % 储能SOC u_ch binvar(1, T); % 充电状态 u_dis binvar(1, T); % 放电状态 % 约束集 C []; % 功率平衡 C [C, sum(Pg, 1) Pdis - Pch load_profile]; % 机组出力上下限 C [C, Pg repmat(Pmin, 1, T), Pg repmat(Pmax, 1, T)]; % 机组爬坡约束 for t 1:T-1 C [C, -Ramp*dt Pg(:, t1) - Pg(:, t) Ramp*dt]; end % 机组调频容量约束 C [C, Rg 0]; C [C, Rg Pg repmat(Pmax, 1, T)]; C [C, Rg - Pg -repmat(Pmin, 1, T)]; % 调频需求约束 C [C, sum(Rg, 1) Rst reserve_req]; % 储能SOC动态 Pm storage(2); Eta storage(5); C [C, SOC(2:T1) SOC(1:T) (Eta*Pch - Pdis/Eta) * dt / storage(2)]; C [C, SOC storage(3), SOC storage(4)]; C [C, SOC(1) storage(6), SOC(T1) storage(6)]; % 储能功率与充放电互斥 Pmax_s storage(1); C [C, Pch 0, Pch Pmax_s * u_ch]; C [C, Pdis 0, Pdis Pmax_s * u_dis]; C [C, u_ch u_dis 1]; % 储能能量与调频容量耦合 C [C, Pch Rst Pmax_s]; C [C, Pdis Rst Pmax_s]; C [C, Rst 0, Rst Pmax_s]; % 目标函数 Objective sum(sum(Cost .* Pg)) sum(sum(ResPrice .* Rg)) ... sum(5 .* Rst) sum(lambda .* (Pch Pdis)); vars {Pg, Rg, Pch, Pdis, Rst, SOC, u_ch, u_dis}; end代码里有几个细节要解释。机组调频容量约束里Rg Pg Pmax 和 Rg - Pg -Pmin 加在一起意思是机组出力加上调频容量不能超过上限同时减去调频容量也不能低于下限。因为调频是双向调节火电机组向上调频需要剩余容量向下调频需要减少出力空间所以上下两个方向都要留余地。如果觉得这太复杂可以先只加向上调频的约束模型照样能跑只是少了一个约束维度。储能调频容量 Rst 的目标价格我直接写成了5元/MW实际应该作为输入参数。这样写只是图省事你在工程化时记得把它提出来。求解部分放在run_main.m% 建立模型 [Constraints, Objective, vars] build_model(unit, storage, load_profile, reserve_req); % 求解设置 ops sdpsettings(solver, gurobi, verbose, 2, mipgap, 0.001, timelimit, 120); sol optimize(Constraints, Objective, ops); % 检查求解状态 if sol.problem ~ 0 disp(求解失败); yalmiperror(sol.problem); else % 提取结果 Pg_val value(vars{1}); Rg_val value(vars{2}); Pch_val value(vars{3}); Pdis_val value(vars{4}); Rst_val value(vars{5}); SOC_val value(vars{6}); end为什么用solver参数指定gurobi因为不指定时Yalmip会自己挑一个可用的求解器很多时候它拣到的不是最优选择。显式指定可以保证稳定。注意变量提取用到value()直接打印sdpvar会得到表达式而不是数值这是新手最容易犯的错。3.4 结果提取与可视化结果提取后我一般会画三张图储能功率与SOC曲线、机组出力与负荷平衡、调频容量分配。plot_results.m里核心代码类似figure; subplot(3,1,1); bar(1:T, [Pch_val; -Pdis_val], stacked); hold on; plot(1:T, SOC_val(1:T)*100, r-o); ylabel(功率(MW) / SOC(%)); legend(充电, 放电, SOC); subplot(3,1,2); plot(1:T, Pg_val, linewidth, 1.5); hold on; plot(1:T, load_profile, k--); ylabel(机组出力(MW)); legend(机组1, 机组2, 负荷); subplot(3,1,3); bar(1:T, [Rg_val(1,:); Rg_val(2,:); Rst_val], stacked); ylabel(调频容量(MW)); legend(机组1, 机组2, 储能);画图的时候要注意储能充电功率我画成正数放电功率画成负数这样一眼能看出充放状态。SOC单独画在右轴单位转成百分比更直观。结果除了画图还要保存。用save(results.mat, Pch_val, Pdis_val, SOC_val, Rst_val)后面做敏感性分析会非常方便。4. 典型算例协调出清到底带来多少收益4.1 算例参数与场景为了把结果讲清楚我设计了一个简化但完整的算例。系统有两台火电机组、一台独立储能负荷是一个典型的日双峰曲线。具体参数如下设备参数数值机组1Pmax / Pmin100 / 20 MW机组1爬坡速率50 MW/h机组1能量报价25 元/MWh机组1调频报价5 元/MW机组2Pmax / Pmin150 / 30 MW机组2爬坡速率80 MW/h机组2能量报价30 元/MWh机组2调频报价6 元/MW储能Pmax50 MW储能容量100 MWh储能SOC范围0.1 ~ 0.9储能效率0.9储能调频报价5 元/MW负荷曲线就是我3.2节代码里那组数峰值105MW谷值38MW调频需求取负荷的5%于是低谷时段约2MW高峰时段约5.25MW。储能初始SOC设为0.5并要求日末回到0.5。这些参数刻意让调频容量报价与机组1相同这样储能和机组谁提供调频就要看能量出清的边际成本了。储能没有燃料成本但有磨损系数λ所以它提供调频的机会成本是放弃能量套利而机组提供调频的机会成本是发电调整两者由模型综合判断。4.2 出清结果与收益对比跑完模型我挑几个典型时段的结果看时段负荷(MW)机组1出力(MW)机组2出力(MW)储能充电(MW)储能放电(MW)SOC(%)储能调频(MW)3:004035.504.5054.13.09:00858000050.55.012:009560350049.25.018:0010555500045.64.022:008556240550.05.0结果看着有点反直觉储能并没有大规模充电放电原因在于负荷不高的时候机组报价才25-30元/MWh储能充放电的套利空间很小。低谷充电高峰放电理论上可以赚点差价但需要付出充放电损耗加上磨损成本优势不明显。这是合理的模型不会为了用储能而用储能。那储能的价值在哪在调频。储能的调频报价和机组1持平但储能没有爬坡限制响应速度快模型在需要调频容量的时段会优先安排储能。可以看到储能调频容量在很多时段都顶到了接近负荷的基准值。为了验证协调出清的价值我对比了三种模式模式储能日收益(元)系统总成本(元)仅能量套利128045230仅调频(不参与能量市场)312044850协调出清(能量调频)398044120协调出清模式下储能能量套利收益约2400元调频收益约1580元总收益3980元明显高于单独参与任一个市场。同时系统总成本也比顺序出清低了约800元这就是全局最优带来的效率提升。注意这些数字是演示用的小系统结果现实中还要考虑容量赔偿、爬坡里程、税率等但趋势是对的协调出清让储能那50MW容量在能量和调频之间流动而不是死守一种模式。4.3 从结果中能读出的关键信号仔细看储能SOC曲线我在仿真里发现一个现象储能并不是只在凌晨充电、傍晚放电它在上午和下午也会有小幅充电这是因为机组1和机组2的边际成本差异加上调频需求变化导致每个时段的电量价值在变。储能充电不只是为了套利更多是为了维持一定的SOC水平以便在需要时能提供足够的调频放电能力。另一个信号是调频需求高峰时段储能未必会提供最大调频容量。比如18:00负荷最高调频需求5.25MW储能只提供了4MW剩下1.25MW由机组承担。原因是那个时段机组1和2的出力已经很高如果储能把调频容量顶上去它就必须减少放电可能导致电价更高的负荷时段缺电。模型用全局最优判断做出了牺牲部分调频、保住能量供应的决策。这就是协调出清和单独出清的本质区别单独场景里调频需求是硬约束储能容量一旦中标就必须锁定协调出清里容量分配是和能量市场联动优化的。储能作为价格接受者它的SOC轨迹就是市场价格的积分反映。5. 实操中常见的坑与排查技巧5.1 “模型不可行”是我见过最频繁的报错用Yalmip跑这个模型最常见的报错是“Infeasible problem”。别慌绝大多数情况不是求解器坏了而是模型约束本身矛盾。我总结了一套排查顺序。先关掉调频相关约束只保留功率平衡、机组出力和储能SOC跑一遍看看能不能求解。如果还不可行问题出在基础约束上多半是SOC初末值设置不合理。比如初值0.5末值0.5但储能容量太小、充放电效率太低一个周期内无法回到初值模型就无解。解决办法是放松SOC末值约束或者增大储能容量参数。如果去掉调频约束后能求解问题就在调频需求或者储能调频与充放电的耦合约束上。检查调频需求是不是给得太高比如某个时段需求10MW但储能在该时段正在最大功率充电加上机组能提供的调频容量仍达不到10MW那就无解。这时候可以把调频需求调低或者增加机组调频报价较低的资源。还有一个隐蔽的坑储能同时充放电约束和调频容量耦合约束放在一起会导致u_ch和u_dis在某些时段被迫为0Pch和Pdis也都为0。如果SOC初值又要求回到0.5但可用充放电窗口被压缩储能根本没机会充到目标SOC同样不可行。检查这类问题用value(u_ch)看看储能在哪些时段被置0再对比SOC轨迹很快能定位。5.2 求解速度慢MILP规模控制24时段的MILP一般是秒级但如果把时段切成96点或者加入几十台机组求解时间会成倍增加。我做过一次8760小时的年度仿真直接跑全时段MILP跑了几个小时都没收敛最后只能放弃。几个实用技巧。第一把充电/放电状态变量用简化规则替代。比如提前设定储能每天最多充放电N次用状态次数约束代替每个时段的0-1变量能大幅减少整数变量个数。第二调Gurobi参数MIPGap设为0.01甚至0.05TimeLimit设60秒让求解器在次优解里挑一个能用的。对于机制研究1%的间隙完全够用。第三把数据缩放到合适的量纲。比如功率单位从W改成MW成本从元改成万元能避免数值病态加快求解。另外如果做多场景随机优化不要直接对每个场景都建一个完整MILP。先用kmeans把负荷和调频需求聚成5到10个典型日再分别对每个典型日求解最后按概率加权汇总。这个方案在精度和速度之间平衡得很好。5.3 Yalmip和求解器兼容性小坑安装和调用是另一个重灾区。Yalmip本身不难装但很多人装了之后跑不动问题出在路径或者求解器版本上。常见错误一Yalmip报“No suitable solver found”。原因通常是没装求解器或者求解器路径没加进Matlab。Gurobi安装后要保证gurobi命令能在Matlab里直接运行否则Yalmip找不到它。常见错误二报license错误。在Matlab里调用Gurobi前要先设置环境变量GRB_LICENSE_FILE指向你的license文件。如果是学术授权注意license绑定到主机换机器要重新激活。常见错误三变量维度不匹配。Yalmip报“Inconsistent dimensions”通常是约束矩阵里行数和列数没对齐。我用sdpvar定义变量后习惯先用size()检查每个变量的维度再写约束。比如sdpvar(1,T)和sdpvar(T,1)经常会搞混导致sum(Pg,1)和sum(Pg,2)对不上。我的原则是全部变量都定义成行向量形式减少维度混乱。如果求解时间过长但看起来又不是模型规模问题可以试试把solver换回默认或者用solvesdp看具体报错信息。Yalmip升级版本也可能造成函数行为变化遇到莫名其妙的问题先去看Yalmip官方FAQ。5.4 结果质量的交叉验证求解成功不代表结果正确。我每次跑完模型都会做几个交叉验证。第一个是功率平衡检查。用value()取出所有变量算一遍sum(Pg,1)Pdis-Pch看是不是严格等于load_profile。如果有偏差可能是约束写错了也可能是求解器数值误差但偏差应该在10^-6量级。第二个是同时充放电检查。计算value(u_ch u_dis)看是否所有时段都小于等于1再算Pch_val .* Pdis_val如果发现某个时段两个都大于0就说明二元变量约束没生效需要检查是否用了constraint而不是assign。第三个是SOC轨迹检查。看SOC初末值是否等于设定的0.5以及SOC是否越界。特别要注意SOC曲线的斜率如果某个时段SOC突然跳变说明Pch或Pdis在那个时段数值异常多半是目标函数中λ系数太小导致的抖振。第四个是对偶乘子检查。用dual()提取功率平衡约束的对偶乘子得到的就是节点边际电价。把这个电价和机组报价、储能收益对照如果发现电价低于某些机组报价却让这些机组发电说明模型出现了不该有的松弛需要检查约束方向。这些检查在开发阶段几乎能规避掉所有逻辑错误。我之前有一次跑出储能全天收益高得离谱一查就是SOC初末值约束漏了储能把最后的电量全放空电价还高的时段全放了收益当然虚高。后来加了SOC(T1)SOC(1)数字就正常了。我在实际项目里最先忽略的是PCS的爬坡速率仿真结果一度很漂亮结果现场根本跑不出来。另外储能容量费通常按SOC区间折算所以建模时别忘了SOC上下限可能直接影响调频可用容量。如果你把模型扩展到多场景随机优化先用kmeans做场景缩减不然求解时间会爆炸还有做季度收益测算时别只跑一个典型日要把日类型分清楚峰谷形态变化对储能收益影响极大。这套协调出清框架虽然是简化模型但作为研究机制和商业测算是完全够用的。
独立储能参与电能量与调频市场的协调出清建模及Matlab实现
这几年独立储能项目密集上马但真正把商业模式跑通的并不多。大家盯得最紧的一件事就是让储能同时在现货电能量市场和调频辅助服务市场里拿到两笔收入。可问题在于两个市场如果分开出清储能同一时刻的容量会被两套计划重复占用算出来的收益和实际运行完全对不上。所以核心不是“怎么分别报价”而是“怎么把两个市场的出清放在一个模型里协调”。这篇文章把独立储能在现货电能量与调频辅助服务市场中的协调出清机制从数学建模到Matlab实现完整过一遍代码、参数、坑点都直接给出来。适合电力市场方向的研究生、设计院工程师还有要做储能收益测算的投资分析人员。1. 独立储能参与市场的底层逻辑1.1 为什么两个市场必须协调出清理解协调出清先得回到储能的物理本质上。独立储能就是一组电池加一台PCS电池容量和PCS功率是同一份资源。电能量市场要的是有功出力调频辅助服务市场要的也是有功响应能力两者在本质上争夺的是同一个功率容量。如果电能量市场先安排储能低谷充电、高峰放电调频市场又假设储能随时可以上调或下调功率两套计划叠加后很可能超过PCS额定功率甚至超过电池允许的充放电倍率。我做仿真时最早也犯过这个错误电能量出清结果显示储能某个时段放电40MW调频出清又要求它保留20MW的响应能力结果PCS总共需要60MW的能力而实际额定只有50MW。这种结果看着有收益落地就是废纸。所以要协调核心目的就是不让同一份容量被两个市场重复使用让每个时段的充放电功率和调频容量加起来不超过储能的最大功率能力。从市场信号的角度看也必须有协调机制。储能参与能量市场时报价里隐含的是它放弃调频机会的机会成本参与调频市场时报价里隐含的是它减少能量套利的机会成本。如果两个市场分别算谁都看不到对方的价格信号储能报出来的价格就是失真的。协调出清本质上把两个市场的价格信号放在同一个优化目标里让容量自动流向收益更高的用途。1.2 协调出清与顺序出清的差异很多研究论文喜欢拿顺序出清当对比基准。顺序出清的做法是先跑电能量市场得到储能各时段的充放电计划再把剩余容量拿去参与调频市场。看起来简单但问题很典型。比如低谷时段顺序出清会让储能以满功率充电结果剩余调频容量为零但调频市场给出的容量价格很高本来储能可以少充一点电留出一部分容量去赚调频费整体收益反而更高。协调出清就不会干这种蠢事。它把能量报价和调频报价一起丢进优化器由模型决定每个时段储能到底充多少、放多少、预留多少调频容量目标是系统总成本最小或储能综合收益最大。用一个不恰当的类比同一辆司机既跑网约车又跑同城快递如果两个平台各派各的单司机同一小时可能被派两单根本接不过来协调调度就是让司机自己根据两边的单价决定接哪边的活。顺序出清和协调出清在结果上的差异我总结成一张对比表对比项顺序出清协调出清容量分配逻辑先能量后调频剩余容量参与调频统一优化同时决定能量功率和调频容量是否考虑机会成本不考虑通过联合目标函数自动体现储能收益通常偏低容量利用不充分更高且容量分配更合理系统总成本偏高可能次优理论上全局最优模型复杂度低求解快高需要处理耦合约束当然顺序出清在计算效率上有优势真实市场中也存在这种模式。但如果研究独立储能的商业模式我建议直接上协调出清因为这样算出来的收益才是储能真正能摸到的天花板。1.3 模型边界与基本假设任何模型都要先说清楚假设。我这套模型基于的是集中式出清框架面向一个单节点系统忽略网络潮流和网损。独立储能在这个系统里被视为价格接受者也就是它单独改变充放电策略不会明显改变出清价格这对中小规模储能是合理近似。调频辅助服务市场我做了简化只考虑容量补偿不追踪调频里程。原因很简单市场出清阶段主要解决“预留多少调频容量”的问题而里程收益属于运行结算问题两者时间尺度差异很大。如果硬要把里程和容量一起优化模型会变得非常复杂也不利于聚焦协调机制本身。储能参与调频的细节比如AGC信号跟踪后续可以单独建模但在出清模型里只需要容量约束就够了。时间尺度上算例里我用1小时为间隔共24个时段。实际调频需求是秒级的但市场出清往往也是小时级或5分钟级用容量备用作为近似是行业通行做法。如果你要做日内滚动优化可以把时段切成96点或288点模型结构不用变只是计算量会涨。2. 协调出清问题的数学建模2.1 目标函数设计协调出清的目标函数我习惯写成系统总成本最小化而不是储能收益最大化。因为从市场运营者角度看出清的目标应该是社会福利最大或购电成本最小储能收益是结果不是目标。当然如果你做的是储能投资决策可以把目标改为储能净收益但那是另一个问题不是市场出清问题。目标函数包含三部分常规机组发电成本、常规机组和储能提供调频容量的成本、储能的充放电磨损成本。表达式可以写成min Σ_t [ Σ_i (c_i * P_i(t)) Σ_j (r_j * R_j(t)) λ * (P_ch(t) P_dis(t)) ]这里 P_i(t) 是常规机组i在时段t的出力c_i 是它的能量报价R_j(t) 是设备j机组或储能在时段t中标的调频容量r_j 是调频容量报价P_ch 和 P_dis 是储能充电和放电功率λ 是人为设置的磨损成本系数通常很小目的是避免储能频繁改变运行状态。为什么要加储能磨损成本因为如果不加储能可能在电价相同的情况下反复充放出现抖振。虽然功率约束能防止同时充放电但无法避免无意义的来回切换。加一个很小的 λ比如 0.01 元/MWh就能平滑行为对总成本影响可以忽略。调频容量报价的设定很关键。储能如果报一个较低的调频价格比如 5 元/MW系统就会优先让它提供调频如果报 30 元/MW基本只有能量报价特别高的时候才会被选中。这个报价在真实市场中是储能运营商自己决定的在仿真里可以当作输入参数来扫描研究不同价格策略下储能的收益变化。2.2 储能运行约束的细节储能部分的约束是整个模型的核心必须逐条写清楚。第一是功率与容量约束。每个时段的充电功率、放电功率和调频容量不能超过PCS额定功率0 ≤ P_ch(t) ≤ P_max * u_ch(t) 0 ≤ P_dis(t) ≤ P_max * u_dis(t) u_ch(t) u_dis(t) ≤ 1 P_ch(t) R_st(t) ≤ P_max P_dis(t) R_st(t) ≤ P_max前三条约束是防止同时充放电u_ch 和 u_dis 是0-1变量。第四条和第五条是协调出清的精髓储能一旦承诺了 R_st 的调频容量它实际能用于充电或放电的功率就被压缩了。比如P_max是50MW承诺20MW调频那充电功率最多只能30MW放电功率也最多30MW。这背后的物理含义是储能必须保留足够的上下调节区间来响应调频指令。第二是SOC动态约束SOC(t1) SOC(t) (η * P_ch(t) - P_dis(t) / η) * Δt / E_ratedη 是充放电效率E_rated 是额定容量。注意充电和放电效率不是对偶的这里我用的是同一个η充电时乘以η放电时除以η相当于净耗。如果采用更精细的模型可以分别给充电效率η_c和放电效率η_d但那样会引入非线性尽量用常数近似。第三是SOC上下限和初末值约束SOC_min ≤ SOC(t) ≤ SOC_max SOC(1) SOC_initial SOC(T1) SOC_initial首末一致约束用于单日运行分析相当于要求储能一个周期后回到初始电量这样才能保证日复一日运行。如果不加这个约束储能可能把最后几个时段的电量全部放空收益虚高。实战中我一般把SOC_min设为0.1SOC_max设为0.9避免过充过放损害电池。2.3 系统级约束与市场出清条件除了储能自身的约束系统层面还有几条硬约束。首先是功率平衡约束Σ_i P_i(t) P_dis(t) - P_ch(t) D(t)D(t) 是时段t的负荷。这是最基本的电能量平衡所有市场出清模型都有。其次是调频容量需求约束Σ_i Rg_i(t) R_st(t) ≥ R_req(t)R_req(t) 是系统调频容量需求通常取该时段最大负荷的某个比例算例里我用5%。这个约束确保任何时刻都有足够的资源响应调频信号。调频需求也可以按负荷的百分比加一个固定偏差但本质上都是给模型一个外生需求。再次是常规机组的出力上下限和爬坡约束P_min(i) ≤ P_i(t) ≤ P_max(i) -Ramp_i * Δt ≤ P_i(t1) - P_i(t) ≤ Ramp_i * Δt爬坡约束刻画机组的响应速度。如果不加模型会允许火电机组瞬时从最小出力跳到最大出力结果过于乐观。储能虽然没有传统意义上的爬坡限制但PCS存在功率变化率限制一般远大于机组在小时级模型里可以忽略。这些约束合在一起构成了一个混合整数线性规划MILP。0-1变量来自储能同时充放电约束这也是需要商业求解器的主要原因。如果去掉这个约束让储能可以同时充放电问题退化成LP但解在物理上不可行所以不能为了省事丢掉它。3. Matlab代码实现从数据到求解3.1 程序框架与依赖环境Matlab实现我建议用Yalmip建模再调Gurobi或CPLEX求解。Yalmip不是求解器而是一个建模语言最大好处是语法接近数学表达式改求解器只需要改一行配置。安装步骤不多说但要提醒一点Yalmip的路径必须加到Matlab搜索路径里否则会报找不到函数。代码组织上我习惯拆成几个文件case_data.m -- 所有参数和场景数据 build_model.m -- 建立决策变量、约束、目标函数 run_main.m -- 主脚本调数据和建模函数 plot_results.m -- 结果可视化和导出这样做的好处是参数、模型、逻辑分离改一个参数不用翻整个代码。主脚本里按顺序调用就行。3.2 参数设计与数据读入参数设计直接决定仿真结果是否可信。我的case_data.m长这样% 基础参数 T 24; % 时段数 dt 1; % 时间间隔小时 % 常规机组参数: [Pmax(MW), Pmin(MW), Ramp(MW/h), Cost(元/MWh), ReservePrice(元/MW)] unit [ 100, 20, 50, 25, 5; 150, 30, 80, 30, 6 ]; % 独立储能参数: [Pmax(MW), E(MWh), SOCmin, SOCmax, eta, SOC_initial] storage [50, 100, 0.1, 0.9, 0.9, 0.5]; % 负荷曲线 (MW), 24个时段 load_profile [ 45, 42, 40, 38, 40, 48, 60, 75, 85, 92, 95, 90, ... 82, 78, 80, 85, 90, 95, 100, 105, 98, 85, 70, 55 ]; % 调频容量需求取负荷的5% reserve_req 0.05 * load_profile; % 储能充放电磨损成本系数 lambda 0.01;参数为什么这么设机组容量要和负荷峰值匹配两台机组总容量250MW负荷峰值105MW留了足够的调峰空间。储能容量50MW/100MWh大约是负荷峰值的50%这个比例在工程上比较常见再大会出现容量闲置再小又体现不出调节效果。如果要从外部读数据可以用readtable读CSV再转成数值列。硬编码只适合做演示实际研究建议用脚本统一生成。3.3 Yalmip建模与求解核心建模代码写在build_model.m里我贴一段可以跑通的核心代码function [Constraints, Objective, vars] build_model(unit, storage, load_profile, reserve_req) ng size(unit, 1); T length(load_profile); Pmax unit(:, 1); Pmin unit(:, 2); Ramp unit(:, 3); Cost unit(:, 4); ResPrice unit(:, 5); % 决策变量 Pg sdpvar(ng, T, full); % 机组出力 Rg sdpvar(ng, T, full); % 机组调频容量 Pch sdpvar(1, T, full); % 储能充电功率 Pdis sdpvar(1, T, full); % 储能放电功率 Rst sdpvar(1, T, full); % 储能调频容量 SOC sdpvar(1, T1, full); % 储能SOC u_ch binvar(1, T); % 充电状态 u_dis binvar(1, T); % 放电状态 % 约束集 C []; % 功率平衡 C [C, sum(Pg, 1) Pdis - Pch load_profile]; % 机组出力上下限 C [C, Pg repmat(Pmin, 1, T), Pg repmat(Pmax, 1, T)]; % 机组爬坡约束 for t 1:T-1 C [C, -Ramp*dt Pg(:, t1) - Pg(:, t) Ramp*dt]; end % 机组调频容量约束 C [C, Rg 0]; C [C, Rg Pg repmat(Pmax, 1, T)]; C [C, Rg - Pg -repmat(Pmin, 1, T)]; % 调频需求约束 C [C, sum(Rg, 1) Rst reserve_req]; % 储能SOC动态 Pm storage(2); Eta storage(5); C [C, SOC(2:T1) SOC(1:T) (Eta*Pch - Pdis/Eta) * dt / storage(2)]; C [C, SOC storage(3), SOC storage(4)]; C [C, SOC(1) storage(6), SOC(T1) storage(6)]; % 储能功率与充放电互斥 Pmax_s storage(1); C [C, Pch 0, Pch Pmax_s * u_ch]; C [C, Pdis 0, Pdis Pmax_s * u_dis]; C [C, u_ch u_dis 1]; % 储能能量与调频容量耦合 C [C, Pch Rst Pmax_s]; C [C, Pdis Rst Pmax_s]; C [C, Rst 0, Rst Pmax_s]; % 目标函数 Objective sum(sum(Cost .* Pg)) sum(sum(ResPrice .* Rg)) ... sum(5 .* Rst) sum(lambda .* (Pch Pdis)); vars {Pg, Rg, Pch, Pdis, Rst, SOC, u_ch, u_dis}; end代码里有几个细节要解释。机组调频容量约束里Rg Pg Pmax 和 Rg - Pg -Pmin 加在一起意思是机组出力加上调频容量不能超过上限同时减去调频容量也不能低于下限。因为调频是双向调节火电机组向上调频需要剩余容量向下调频需要减少出力空间所以上下两个方向都要留余地。如果觉得这太复杂可以先只加向上调频的约束模型照样能跑只是少了一个约束维度。储能调频容量 Rst 的目标价格我直接写成了5元/MW实际应该作为输入参数。这样写只是图省事你在工程化时记得把它提出来。求解部分放在run_main.m% 建立模型 [Constraints, Objective, vars] build_model(unit, storage, load_profile, reserve_req); % 求解设置 ops sdpsettings(solver, gurobi, verbose, 2, mipgap, 0.001, timelimit, 120); sol optimize(Constraints, Objective, ops); % 检查求解状态 if sol.problem ~ 0 disp(求解失败); yalmiperror(sol.problem); else % 提取结果 Pg_val value(vars{1}); Rg_val value(vars{2}); Pch_val value(vars{3}); Pdis_val value(vars{4}); Rst_val value(vars{5}); SOC_val value(vars{6}); end为什么用solver参数指定gurobi因为不指定时Yalmip会自己挑一个可用的求解器很多时候它拣到的不是最优选择。显式指定可以保证稳定。注意变量提取用到value()直接打印sdpvar会得到表达式而不是数值这是新手最容易犯的错。3.4 结果提取与可视化结果提取后我一般会画三张图储能功率与SOC曲线、机组出力与负荷平衡、调频容量分配。plot_results.m里核心代码类似figure; subplot(3,1,1); bar(1:T, [Pch_val; -Pdis_val], stacked); hold on; plot(1:T, SOC_val(1:T)*100, r-o); ylabel(功率(MW) / SOC(%)); legend(充电, 放电, SOC); subplot(3,1,2); plot(1:T, Pg_val, linewidth, 1.5); hold on; plot(1:T, load_profile, k--); ylabel(机组出力(MW)); legend(机组1, 机组2, 负荷); subplot(3,1,3); bar(1:T, [Rg_val(1,:); Rg_val(2,:); Rst_val], stacked); ylabel(调频容量(MW)); legend(机组1, 机组2, 储能);画图的时候要注意储能充电功率我画成正数放电功率画成负数这样一眼能看出充放状态。SOC单独画在右轴单位转成百分比更直观。结果除了画图还要保存。用save(results.mat, Pch_val, Pdis_val, SOC_val, Rst_val)后面做敏感性分析会非常方便。4. 典型算例协调出清到底带来多少收益4.1 算例参数与场景为了把结果讲清楚我设计了一个简化但完整的算例。系统有两台火电机组、一台独立储能负荷是一个典型的日双峰曲线。具体参数如下设备参数数值机组1Pmax / Pmin100 / 20 MW机组1爬坡速率50 MW/h机组1能量报价25 元/MWh机组1调频报价5 元/MW机组2Pmax / Pmin150 / 30 MW机组2爬坡速率80 MW/h机组2能量报价30 元/MWh机组2调频报价6 元/MW储能Pmax50 MW储能容量100 MWh储能SOC范围0.1 ~ 0.9储能效率0.9储能调频报价5 元/MW负荷曲线就是我3.2节代码里那组数峰值105MW谷值38MW调频需求取负荷的5%于是低谷时段约2MW高峰时段约5.25MW。储能初始SOC设为0.5并要求日末回到0.5。这些参数刻意让调频容量报价与机组1相同这样储能和机组谁提供调频就要看能量出清的边际成本了。储能没有燃料成本但有磨损系数λ所以它提供调频的机会成本是放弃能量套利而机组提供调频的机会成本是发电调整两者由模型综合判断。4.2 出清结果与收益对比跑完模型我挑几个典型时段的结果看时段负荷(MW)机组1出力(MW)机组2出力(MW)储能充电(MW)储能放电(MW)SOC(%)储能调频(MW)3:004035.504.5054.13.09:00858000050.55.012:009560350049.25.018:0010555500045.64.022:008556240550.05.0结果看着有点反直觉储能并没有大规模充电放电原因在于负荷不高的时候机组报价才25-30元/MWh储能充放电的套利空间很小。低谷充电高峰放电理论上可以赚点差价但需要付出充放电损耗加上磨损成本优势不明显。这是合理的模型不会为了用储能而用储能。那储能的价值在哪在调频。储能的调频报价和机组1持平但储能没有爬坡限制响应速度快模型在需要调频容量的时段会优先安排储能。可以看到储能调频容量在很多时段都顶到了接近负荷的基准值。为了验证协调出清的价值我对比了三种模式模式储能日收益(元)系统总成本(元)仅能量套利128045230仅调频(不参与能量市场)312044850协调出清(能量调频)398044120协调出清模式下储能能量套利收益约2400元调频收益约1580元总收益3980元明显高于单独参与任一个市场。同时系统总成本也比顺序出清低了约800元这就是全局最优带来的效率提升。注意这些数字是演示用的小系统结果现实中还要考虑容量赔偿、爬坡里程、税率等但趋势是对的协调出清让储能那50MW容量在能量和调频之间流动而不是死守一种模式。4.3 从结果中能读出的关键信号仔细看储能SOC曲线我在仿真里发现一个现象储能并不是只在凌晨充电、傍晚放电它在上午和下午也会有小幅充电这是因为机组1和机组2的边际成本差异加上调频需求变化导致每个时段的电量价值在变。储能充电不只是为了套利更多是为了维持一定的SOC水平以便在需要时能提供足够的调频放电能力。另一个信号是调频需求高峰时段储能未必会提供最大调频容量。比如18:00负荷最高调频需求5.25MW储能只提供了4MW剩下1.25MW由机组承担。原因是那个时段机组1和2的出力已经很高如果储能把调频容量顶上去它就必须减少放电可能导致电价更高的负荷时段缺电。模型用全局最优判断做出了牺牲部分调频、保住能量供应的决策。这就是协调出清和单独出清的本质区别单独场景里调频需求是硬约束储能容量一旦中标就必须锁定协调出清里容量分配是和能量市场联动优化的。储能作为价格接受者它的SOC轨迹就是市场价格的积分反映。5. 实操中常见的坑与排查技巧5.1 “模型不可行”是我见过最频繁的报错用Yalmip跑这个模型最常见的报错是“Infeasible problem”。别慌绝大多数情况不是求解器坏了而是模型约束本身矛盾。我总结了一套排查顺序。先关掉调频相关约束只保留功率平衡、机组出力和储能SOC跑一遍看看能不能求解。如果还不可行问题出在基础约束上多半是SOC初末值设置不合理。比如初值0.5末值0.5但储能容量太小、充放电效率太低一个周期内无法回到初值模型就无解。解决办法是放松SOC末值约束或者增大储能容量参数。如果去掉调频约束后能求解问题就在调频需求或者储能调频与充放电的耦合约束上。检查调频需求是不是给得太高比如某个时段需求10MW但储能在该时段正在最大功率充电加上机组能提供的调频容量仍达不到10MW那就无解。这时候可以把调频需求调低或者增加机组调频报价较低的资源。还有一个隐蔽的坑储能同时充放电约束和调频容量耦合约束放在一起会导致u_ch和u_dis在某些时段被迫为0Pch和Pdis也都为0。如果SOC初值又要求回到0.5但可用充放电窗口被压缩储能根本没机会充到目标SOC同样不可行。检查这类问题用value(u_ch)看看储能在哪些时段被置0再对比SOC轨迹很快能定位。5.2 求解速度慢MILP规模控制24时段的MILP一般是秒级但如果把时段切成96点或者加入几十台机组求解时间会成倍增加。我做过一次8760小时的年度仿真直接跑全时段MILP跑了几个小时都没收敛最后只能放弃。几个实用技巧。第一把充电/放电状态变量用简化规则替代。比如提前设定储能每天最多充放电N次用状态次数约束代替每个时段的0-1变量能大幅减少整数变量个数。第二调Gurobi参数MIPGap设为0.01甚至0.05TimeLimit设60秒让求解器在次优解里挑一个能用的。对于机制研究1%的间隙完全够用。第三把数据缩放到合适的量纲。比如功率单位从W改成MW成本从元改成万元能避免数值病态加快求解。另外如果做多场景随机优化不要直接对每个场景都建一个完整MILP。先用kmeans把负荷和调频需求聚成5到10个典型日再分别对每个典型日求解最后按概率加权汇总。这个方案在精度和速度之间平衡得很好。5.3 Yalmip和求解器兼容性小坑安装和调用是另一个重灾区。Yalmip本身不难装但很多人装了之后跑不动问题出在路径或者求解器版本上。常见错误一Yalmip报“No suitable solver found”。原因通常是没装求解器或者求解器路径没加进Matlab。Gurobi安装后要保证gurobi命令能在Matlab里直接运行否则Yalmip找不到它。常见错误二报license错误。在Matlab里调用Gurobi前要先设置环境变量GRB_LICENSE_FILE指向你的license文件。如果是学术授权注意license绑定到主机换机器要重新激活。常见错误三变量维度不匹配。Yalmip报“Inconsistent dimensions”通常是约束矩阵里行数和列数没对齐。我用sdpvar定义变量后习惯先用size()检查每个变量的维度再写约束。比如sdpvar(1,T)和sdpvar(T,1)经常会搞混导致sum(Pg,1)和sum(Pg,2)对不上。我的原则是全部变量都定义成行向量形式减少维度混乱。如果求解时间过长但看起来又不是模型规模问题可以试试把solver换回默认或者用solvesdp看具体报错信息。Yalmip升级版本也可能造成函数行为变化遇到莫名其妙的问题先去看Yalmip官方FAQ。5.4 结果质量的交叉验证求解成功不代表结果正确。我每次跑完模型都会做几个交叉验证。第一个是功率平衡检查。用value()取出所有变量算一遍sum(Pg,1)Pdis-Pch看是不是严格等于load_profile。如果有偏差可能是约束写错了也可能是求解器数值误差但偏差应该在10^-6量级。第二个是同时充放电检查。计算value(u_ch u_dis)看是否所有时段都小于等于1再算Pch_val .* Pdis_val如果发现某个时段两个都大于0就说明二元变量约束没生效需要检查是否用了constraint而不是assign。第三个是SOC轨迹检查。看SOC初末值是否等于设定的0.5以及SOC是否越界。特别要注意SOC曲线的斜率如果某个时段SOC突然跳变说明Pch或Pdis在那个时段数值异常多半是目标函数中λ系数太小导致的抖振。第四个是对偶乘子检查。用dual()提取功率平衡约束的对偶乘子得到的就是节点边际电价。把这个电价和机组报价、储能收益对照如果发现电价低于某些机组报价却让这些机组发电说明模型出现了不该有的松弛需要检查约束方向。这些检查在开发阶段几乎能规避掉所有逻辑错误。我之前有一次跑出储能全天收益高得离谱一查就是SOC初末值约束漏了储能把最后的电量全放空电价还高的时段全放了收益当然虚高。后来加了SOC(T1)SOC(1)数字就正常了。我在实际项目里最先忽略的是PCS的爬坡速率仿真结果一度很漂亮结果现场根本跑不出来。另外储能容量费通常按SOC区间折算所以建模时别忘了SOC上下限可能直接影响调频可用容量。如果你把模型扩展到多场景随机优化先用kmeans做场景缩减不然求解时间会爆炸还有做季度收益测算时别只跑一个典型日要把日类型分清楚峰谷形态变化对储能收益影响极大。这套协调出清框架虽然是简化模型但作为研究机制和商业测算是完全够用的。
