考虑公平分配与风险管理的社区微电网多阶段优化策略——基于VCG机制与GlueVaR方法

《Journal of Modern Power Systems and Clean Energy》：Multi-Stage Optimization of Community Microgrid Considering Fair Allocation and Risk Management

　　

摘要：本文提出一种考虑公平分配与风险管理的社区微电网（Community Microgrid, CMG）多阶段优化调度策略。该策略将碳配额（Carbon Quota）引入传统社区共享消费（Community Shared Consumption, CSC）框架，使配有光伏（Photovoltaic, PV）单元、储能系统（Energy Storage System, ESS）及燃气机组（Gas Unit/Turbine）的用户可在日前市场进行电能与碳配额的内部交易。为应对PV出力不确定性带来的运行风险，研究人员采用Glue Value-at-Risk（GlueVaR）——一种融合不同置信水平下VaR与条件VaR（Conditional Value-at-Risk, CVaR）的组合风险度量方法——构建计及风险偏好（Risk Preference）的日前阶段优化模型。实时阶段通过二阶锥规划（Second-Order Cone Programming, SOCP）松弛的交流最优潮流（Alternating Current Optimal Power Flow, AC-OPF）进行电压调节与偏差校正。月度结算阶段引入Vickrey–Clarke–Groves（VCG）机制确定内部能—碳分配系数，抑制虚假申报并确保激励相容（Incentive Compatibility）与分配公平。模型采用交替方向乘子法（Alternating Direction Method of Multipliers, ADMM）分布式求解。基于IEEE欧洲低压（Low-Voltage, LV）测试馈线55节点算例验证表明：所提策略较传统仅电能调度方案显著提升社区总收益约8%，促进可再生能源消纳，VCG机制使Jain公平指数达0.7156（购电）与0.6431（售电），GlueVaR提供灵活的风险—成本权衡手段，实时阶段可有效维持节点电压于0.95~1.05 p.u.安全范围内。

论文解读：考虑公平分配与风险管理的社区微电网多阶段优化策略——基于VCG机制与GlueVaR方法

一、研究背景与意义

随分布式电源（Distributed Generation, DG）与产消者（Prosumer）大量接入社区微电网（Community Microgrid, CMG），传统集中式调度难以兼顾隐私保护与个体理性，且单纯电能共享未充分利用碳减排激励机制。现有社区共享消费（Community Shared Consumption, CSC）框架忽略：（1）光伏（Photovoltaic, PV）出力不确定性引发的运行风险；（2）碳配额（Carbon Quota / Allowance）可内部交易带来的额外收益与减排激励；（3）成员虚假申报电量造成的分配不公与外部性。为此，研究人员在CSC基础上耦合碳配额交易与风险敏感优化，引入GlueVaR（Glue Value-at-Risk）处理多置信水平风险谱，并在结算阶段采用Vickrey–Clarke–Groves（VCG）机制实现真实申报诱导与公平分配，最终形成日前—实时—月度三阶段优化架构。该成果发表于Journal of Modern Power Systems and Clean Energy。

二、主要关键技术方法

研究人员以IEEE欧洲低压（Low-Voltage, LV）测试馈线55节点含PV/ESS/燃气机组（Gas Turbine / Gas Unit）的典型社区微网格为算例，PV出力取自ECMWF与NASA历史再分析数据经蒙特卡洛模拟与K-means聚类生成5种场景。日前阶段建立含GlueVaR的风险—成本加权目标函数（式(41)），燃气机组按式(9)~(12)、PV按式(13)~(15)、ESS按SOC（State of Charge）动力学式(16)~(21)建模，节点功率平衡按式(22)~(24)，碳配额初始分配按式(25)，碳排放按式(26)，净买入/卖出碳量按式(27)~(28)，买卖成本与收益按式(29)~(32)，GlueVaR由式(35)~(37)定义并以辅助变量线性化（式(38)~(40)）。实时阶段以最小化日前—实时偏差为目标（式(47)~(48)），通过SOCP松弛DistFlow支路模型（式(49)）及节点电压约束（式(50)）由PMO（Platform Management Operator）协调求解（式(51)）。月度结算按VCG机制计算各成员边际贡献（式(55)~(57)），导出VCG分配系数（式(58)），结算成本按式(59)~(61)，公平性与盈余分别按Jain指数（式(62)）与社区结余（式(63)）评估。整体采用交替方向乘子法（Alternating Direction Method of Multipliers, ADMM）分布式迭代（式(66)~(68)），Gurobi求解。

三、研究结果

Section II. 社区微网成员与网络模型（Local Community Member Model）

研究人员将成员分为燃气成员（Gas Member）、PV成员（PV Member）与ESS成员（ESS Member）三类。燃气机组视效率η_T与天然气热值Q_CH4将耗气量转为视在功率s_i,t^T（式(9)），分解为有功p_i,t^T（式(10)）与无功q_i,t^T（式(11)），运行成本C_i,t^CH₄按气价π_t^CH₄与耗气量计算（式(12)）。PV出力受上限p_i,t≤p_i,t^─PV、视在功率约束(p_i,t^PV)²+(q_i,t^PV)²≤(s_i^PV)²及无功限值|q_i,t^PV|≤p_i,t^PVtanθ_i^PV（式(13)~(15)）。ESS按SOC递推式(16)，周期始末SOC相等（式(17)），净功率为充放电差（式(18)），SOC及充放电功率受上下限约束（式(19)~(21)）。节点净注入功率由负荷、ESS净功、PV、燃气机组代数和给出（式(22)~(23)），线路容量约束为(p_i,t^lin)²+(q_i,t^lin)²≤(s_i^lin)²（式(24)）。结论：分层成员模型可统一描述含DG与ESS的产消者物理与经济约束。

Section III-A. 阶段Ⅰ：日前管理（Day-Ahead Management）

碳配额初分配按装机w_i,t⁰=D(p_i,t^PV+p_i,t^B+p_i,t^T)（式(25)），燃气机组CO₂排放w_i,t^CO₂=αp_i,t^T（式(26)），净需买/卖碳配额按max(w_i,t^CO₂?w_i,t⁰,0)与max(w_i,t⁰?w_i,t^CO₂,0)区分（式(27)~(28)）。成员目标为最小化[购电成本+购碳成本+燃气成本?售电收入?售碳收入]×时段长（式(33)~(34)，含多场景概率?_k）。GlueVaR_α,β=ω₁CVaR_β+ω₂CVaR_α+ω₃VaR_α（式(35)），CVaR按Rockafellar-Uryasev形式引入辅助变量与超额损失z_i,k^α/β≥f_i,k?VaR（式(36)~(40)）。最终目标为(1?L_i)·期望成本+L_i·GlueVaR（式(41)），L_i∈[0,1]为风险厌恶因子。PMO校核碳平衡Σw^CO₂b=Σw^CO₂s（式(42)），按比例分配系数μ_i,t^b/s=w_i,t^CO₂b/s/Φ_t（式(43)~(44)），全局目标最小化社区总成本（式(45)）。结论：融入碳交易扩大可行域并提高收益，GlueVaR允许成员依自身风险偏好权衡经济性—鲁棒性，ADMM分布式求解保护数据隐私。

Section III-B. 阶段Ⅱ：实时调整与电压调节（Real-Time Adjustment and Voltage Regulation）

实时阶段不含碳交易，ESS SOC被约束于日前值±δ_SOC（式(46)），成员先最小化(p_linDA?p_linRT)²（式(47)~(48)）。若无电压越限直接执行；若越限PMO以最小化Σ(p_linDA?p_linRT)²+ΣI_ij,tr_ij并SOCP松弛I_ij,t≥(p_ij,t²+q_ij,t²)/v_i,t（式(49)），节点电压平方v_i,t∈[v_–,i,t²,v_─,i,t²]，对应0.95~1.05 p.u.（式(50)），整体约束见式(51)。结论：实时阶段在允许SOC微调范围内修正日前计划，SOCP松弛AC-OPF可有效消除电压越限且满足实时性要求。

Section III-C. 阶段Ⅲ：月度结算（Monthly Settlement）

按实时平均功率计算时段内购P_i,t^b、售P_i,t^s（式(52)），内部可匹配电量Ω_t=min(ΣP_i,t^b,ΣP_i,t^s)（式(53)）；碳配额类同得W_i,t^CO₂b/s与Φ_t（式(54)）。VCG效用u_i,t^b=P_i,t^b?λ_i,t^bΩ_t等（式(55)），剔除成员i后其余总效用u_?i,t（式(56)），VCG贡献u_i,t^VCGb=Σu_i,t^b?u_?i,t^b（式(57)），分配系数λ_i,t^VCGb=u_i,t^VCGb/Σu_i,t^b（式(58)）。内部费用C_i,t^cmE/cmCO₂与外网费用C_i,t^gdE/gdCO₂按式(59)~(60)，总费C_i,t^td=C_i,t^cm+C_i,t^gd（式(61)），满足不高于直购零售商成本上界（式(64)），PMO获社区结余ΣC_i,t^cm=C_surplus（式(63)），结算目标