跳到正文
1836 · 开源社区

1836 · 开源社区

AI时代的人本基建

登录/注册
← 返回项目《初焰 • 纳宙 第一赛季 - 风象议题 09 水瓶座》

严寒地区 AIDC 微电网钠/锂储能配置研究 · 进展报告 - 20261017

参与共建发布

robinbg · 发布于 2026/10/17 · 更新于 2026/10/17

RO

robinbg2026/10/17

版本：2026-10-17 · 覆盖里程碑 M1–M3（全部验收） 规格基线：SPEC.md v2.9（唯一权威规格，经十一轮复核与两轮阶段零独立审查） 开源仓库：https://github.com/robinbg/aidc-microgrid（本报告引用的全部代码、 配置、验证脚本与产物均已入库）

⚠️ 全部结果为确定性仿真结果，非真实设备测量或现场运行数据。 气象与负荷序列由固定种子的确定性生成器生成；电化学与经济参数中仍有 PLACEHOLDER 项待公开数据标定（见 §6）。本报告不做任何收益承诺： 所有数字仅表示"在该仿真情景下"的结果。

0. 摘要

本项目研究布置在户外严寒环境的 AI 数据中心（AIDC）微电网中，钠离子电池与 磷酸铁锂电池的容量配置问题：内层是小时级调度 MILP（五种策略 × 正常/故障双模式）， 外层是四目标（LCOE、单位 Token 碳排、韧性缺口、SLA 违约率）容量寻优。

截至目前完成三个里程碑并通过 268 项自动化验收测试（0 失败、0 跳过）：

里程碑内容状态关键数字
M1数据与参数层、纯函数✅ 验收11 份 YAML 冻结基线；单位换算只在加载层发生一次
M2内层 MILP 双模式（正常/故障）✅ 验收190 项测试；手算逐项对照通过
M2.5韧性评估（六类故障场景）+ 求解提速✅ 验收244 次真实故障求解 + 84 个显式跳过；168h/2200 任务实例 34.7 s，认证全局 gap 0.2884%
M3单配置闭环（四目标全部可用）✅ 验收268 项测试；高保真 52 窗滚动年度 + 82 场景：A0 f3=0.0694、A3 f3=0.0536（全部为删失下界口径）；试点标定 250 个可行周值产出 S-balance 归一化范围

M3 闭环结果一览（R-NW、S-cost、种子 20260816，全部经独立手算公式复核）：

目标（最小化口径）A0 低保真A3 低保真A0 高保真A3 高保真
f1 LCOE（元/MWh）1349.981476.191282.531417.90
f2 碳排（kgCO₂/MTok）0.028730.028730.0096240.009508
f3 韧性缺口NaN（规格要求）NaN（规格要求）0.06940.0536
f4 SLA 违约率0008.7×10⁻⁶
储能年化成本（元/年）0398.8 万0398.8 万
求解5 窗全可行5 窗全可行52/52 窗可行，21.4 min52/52 窗可行，50.7 min

解读注意（详见 §4.7）：f1 的 A3 > A0 是占位造价下的量化关系而非选型结论； f3 全部由右删失下界构成，当前不可用于配置排序；f2 的高低保真差异来自 年化方法（代表周加权 vs 年度直接求和），两口径不可混比。

最重要的三条发现（详见 §4）：

  1. 低温折减的正确建模位置是"可用能量带"，不是 SOC 动态。把温度折减系数 γ 写进 SOC 更新式会产生 γ² 的凭空往返损耗（γ=0.5 时充 1 MWh 只能放出 0.25 MWh）； 而"收缩 SOC 上界"的写法会让锂电在 −25℃ 下上界（≈0.48）低于初始 SOC（0.50）， MILP 开局即不可行——这正是北方冬季的典型工况。最终规格：SOC 窗口固定， γ 只作用于以 soc_ref 为中心的对称可用带与能量上限。
  2. 当前未标定数据下，储能定容对孤岛续航没有可测影响——所有孤岛观测全部 右删失（连无储能臂 A0 也零失供），韧性目标 f3 是平坦的。这是一条如实 报告的负结果：平坦目标比有噪声目标更危险（看起来干净），在延长故障时长 或完成标定之前，f3 不得用于配置排序。
  3. 大规模调度 MILP 的"认证提速"可以不牺牲任何精度：用确定性规则排程做 可行上界 + 保留全部物理二元变量的辅助模型做有效下界，把原问题 84.4% 的 gap 夹逼到 0.29%——不缩短窗口、不截断任务、不放宽整数性。

1. 研究的问题

AIDC 的负荷特征（高功率密度、Token 产量驱动的 SLA、部分算力任务可迁移）与 严寒地区电力条件（风电反调峰、冬季极端低温、电网碳因子偏高）叠加后，储能 配置需要同时回答四个互相牵制的问题：

  • 经济性：LCOE（含储能全生命周期成本摊销，分母为年服务电量——保证无储能 基线 A0 也有定义、四臂可比）；
  • 低碳性：每百万 Token 净碳排（购电 + 锅炉燃气 − 工业供热替代，制造碳按 项目寿命摊销；碳资产不入账）；
  • 韧性：故障场景下的期望未供电量 EENS 与孤岛续航 T_island 合成的缺口指标 f3 （参考值为全局配置常数 50 MWh / 12 h，不随画像漂移）；
  • 服务质量：SLA 达标率（分母只含可规划的 S2/S3 任务，S1 在线推理已并入 基础负荷，不重复计算）。

低温机制是设计前提：储能户外布置时，钠电与锂电各自的 γ(T) 曲线决定可用容量带 的收窄幅度，进而影响夜间风电消纳、购电量与孤岛续航。本项目只陈述机制， 不预设"钠电更优"的结论方向——若实测显示机制未导致预期差异，同样如实报告。

2. 方法概览

外层（M3–M5）  x = [E_Na, P_Na, E_Li, P_Li]，4 维连续
  ├─ 低保真：5 个 168h 代表周（含全年最冷连续周 W-WINX），权重 0.18/0.16/0.22/0.22/0.22
  │           年化 = Σ 周结果 × 权重 × 365/7；f3 = NaN（代表周承载不了连续故障轨迹）
  ├─ 高保真：8760h 滚动时域（52 窗，SOC 逐窗传递，首末窗固定 soc_init）
  │           + 年度序列上的六类故障场景 → f3（§5.4 唯一公式）
  └─ 对照：离散网格扫描（同预算 hypervolume 对比，如实报告谁更好）

内层（M2，已完成）  小时级 MILP，HiGHS，单线程 + 固定种子
  ├─ 正常模式：需求全满足、任务全完成、SOC 首末相等、市政供热完全补足
  ├─ 故障模式：失供变量显式建模、可迁移任务暂停、四级字典序降级策略
  │             （第一级 gap=0 求到证明最优，否则标 feasible_not_proven，不得冒充）
  └─ 五策略：S-cost / S-carbon / S-heat / S-balance / S-base（各自独立标量目标，
              不是同一组权重的不同取值；S-base 为规则驱动，不调用优化器）

验证体系（与实现同权重）：268 项 pytest（手算对照、四角参数验算、故障链 一致性、退化行为、确定性 1e-9/1e-12）；11 份 YAML 与关键模块 SHA-256 冻结； 每个里程碑产出独立验证脚本与 report.json，全部可复现。

3. 已建成的系统能力

  • 物理模型：功率平衡（PUE 随环境温度变化，非常数）、储能双化学体系 （温度相关效率/倍率/可用带、充放互斥、面积硬约束：钠 45 vs 锂 25 m²/MWh）、 余热回收（κ=0 只计 IT 侧排热，热泵双侧上限，产热上界系数 COP/(COP−1)·η_rec 由配置现算，禁止截断常数）、市政/工业供热分配的精确 min 编码（单二进制 y(t)， 常数差额上界 M）、燃气锅炉备用热源。
  • 算力-电力耦合：S1（在线推理，并入 IT_base）/ S2（硬截止）/ S3（容忍窗口 12–24 h）三级任务；柔性占比按电量口径（分母含 S1）由目标比例反解任务规模， 带可达性预检（四画像 u_req = 0.372/0.468/0.273/0.333，均 ≤ 0.60）。
  • 韧性评估：六类场景（电网失电 2/6/12h、长时失电 24h、风光低分位骤降、 N-1 储能失效、两类组合故障），每场景概率加权；T_island 右删失显式标记； 故障初始 SOC 取正常调度真实衔接值；N-1 按"关键负荷 + 市政供热完整满足"判定。
  • 求解管线：正常模式 168h/2200 任务实例 34.7 s 返回带认证区间的解 （0.2884% 是"距真实最优不超过 0.2884%"的证书，不是"已证明最优"的宣称）； 故障模式四级字典序，第一级强制 gap=0。

4. 关键发现与 insight

4.1 低温机制：γ 的正确位置（两次真实建模错误的教训）

  • 错误一（v1.2）：γ 进入 SOC 动态 → 充放各乘一次 γ，产生 γ² 的凭空往返 损耗。物理上 SOC 是守恒状态量，不存在"温度套利"空间。
  • 错误二（规格初版）："收缩 SOC 上界"→ 锂电 −25℃ 上界 ≈0.48 < 初始 0.50， 北方冬季典型工况开局不可行。
  • 最终形式：SOC 窗口固定；γ 收窄以 soc_ref 为中心的对称可用带（含上下 截断），能量上限带正部截断 [·]⁺ 与效率换算（e_ch ≤ (E/η_ch)[·]⁺、 e_dis ≤ η_dis·E·[·]⁺），功率上限 min(P_nom, c_max(T)·E_nom) 不乘 γ。
  • 物理含义：低温是"装得下但放不出来"（极化增大、可放电深度收窄）， 不是"电池变小"。守门测试专门验证：η=1 时充放往返无损、可用带随 γ 对称收缩、 低温充电裕度约束真实生效。

4.2 负结果：当前数据不能支撑"储能续航优劣"的任何结论

M2.5 对 4 个配置臂 × 82 个场景做了 244 次真实故障求解：全部有效孤岛观测 右删失（censoring_rate = 1.0），延长到 168 h 长时失电诊断仍然如此。原因已从 输入数据核实：该代表周最低风光出力 4.494 MW，而故障期最大基础 IT+制冷需求仅 0.851 MW，燃气锅炉未切断继续补热——电网+风光本身已足够，储能根本没有机会 成为瓶颈。由此：

  • f3 在当前参数下是平坦目标：看起来干净，实则不可分辨配置优劣；
  • 启用 f3 排序前必须先延长故障时长或把风光/负荷标定到储能真正成为瓶颈的区间；
  • 这条负结果与正结果同等重要——它防止了"用平坦指标排出假名次"这一最隐蔽的错误。

4.3 认证提速：84.4% gap → 0.2884%，且没有放松任何东西

168h/2200 任务的 S-cost 实例曾在 60 s 时限下只收敛到 84.4% gap（目标值被 SLA 惩罚项 5:1 主导）。解法不是放宽 gap 或缩短窗口，而是两侧夹逼：

  • 上界：确定性规则排程（S-base 前向规则）产生的完整整数任务排程，经原模型 逐约束独立验证后以 c·x ≤ c·x₀ 割注入——全局最优不会被切掉；
  • 下界：保留全部物理二元变量与物理约束、只把逐任务时序替换为必要总量界的 辅助模型，只用后端对偶界认证，绝不用松弛解的目标值冒充下界；
  • 两侧合并得到包含真实最优的区间，34.7 s 完成，状态如实标 feasible。

该方法论可迁移到任何"整数排程 + 物理调度"耦合的大规模 MILP。

4.4 基线策略的内建不对称（比较时必须携带的注脚）

S-base 的"盈余只充到 soc_ref"规则（BS3）导致其弃电随储能容量放大 （A0 = 0，A1 = 2010.8 MWh/168h）——这是策略约束的后果，不是物理限制。 与其余四策略对比时，基线携带此劣势是内建的；审计字段 s_base_curtailed_from_soc_cap_mwh 使其始终可见而非被掩盖。

4.5 规格工程：独立复核的收益率

阶段零的"只读不写、逐条复述、列出矛盾"流程两轮各发现 9 条与 20 条问题， 全部成立，其中 8 条是规格真实错误（YAML 序列/映射混用的语法错误、 碳公式漏 Δt、SLA 周惩罚被逐时累计 168 次、A0 的 LCOE 分母为零、S3 软截止 与任务窗口矛盾等）。经验："复核过多轮"不等于没有错；结论对但理由错 比结论错更危险（能通过"看起来合理"的审查）；每条质疑都必须附反例或验算。

4.6 诚实状态标签体系

  • feasible ≠ optimal；规则驱动的 S-base 其 mip_gap = null 不是最优性宣称；
  • 故障第一级未证明最优 → feasible_not_proven，其 T_island/censoring 记 NaN、 不计入均值，单列 n_unproven；
  • infeasible 原样上报，禁止折算成 0 成本；未知保持 NaN，不静默为零；
  • 删失观测是下界，聚合时按窗口长度计入并注明——由此算出的 f3 是上界估计。

4.7 M3 闭环的情景级结果（占位参数下，全部如实记录）

  • 经济性（f1）：A3 的 LCOE 高于 A0（高保真 1417.9 vs 1282.5 元/MWh）。 分解：储能年化成本 398.8 万元/年（CRF 摊销 + 运维 + 替换），而运营节省 约 155 万元/年（弃电惩罚缓解 + 谷峰套利，A3 运营成本 2150.6 万 vs A0 2305.9 万）。当前占位造价与惩罚系数下投资未回收——这是标定前的 量化关系，不是"储能不划算"的结论；敏感性分析（M6）正是为此设计。
  • 碳（f2）：高保真下 A3 略低于 A0（0.009508 vs 0.009624 kgCO₂/MTok）： 储能吸收本来弃掉的可再生盈余、小幅减少购电。但电网碳因子在生成数据中 为常数（配置未提供波动幅度），时段平移的减碳渠道未激活，差异仅为 消纳渠道；标定含时变碳因子后该目标才有完整区分度。
  • 韧性（f3）：A0 = 0.0694、A3 = 0.0536，两臂 EENS 均为 0，差值全部来自 孤岛项的删失下界均值（10.33 h vs 10.71 h，censoring_rate = 1.0）。 ⚠️ 两个警示必须随数字携带：① 下界口径 → f3 是上界估计且平坦， 不得用于排序（§4.2）；② A0 跳过了 42 个 N-1 族场景（未安装对应储能）， 两臂孤岛项的加权分母不同，跨臂 f3 比较还需共同场景集对照 （M2.5 已为此提供 common_non_n1 聚合口径，M4 沿用）。
  • SLA（f4）：A3 年度 114,733 个任务中 1 个延迟（f4 = 8.7×10⁻⁶）， A0 为零——当前任务密度下该目标接近饱和，区分度同样待标定。
  • 绿电匹配率口径警示：green_matching_rate 的分母含储能充电、 且不追溯储能出力的来源（solver_info 中已声明），A0（0.9997）与 A3 （0.2351）的差异主要是充电负荷进入分母 + 谷电套利购电所致， 该指标不可跨储能规模直接对比。
  • 求解预算如实记录：60 s/窗预算下，A0 高保真 44/52 窗认证 gap ≤1%； A3 为 20/52（中位 1.17%、最大 5.45%）——全部返回经验证的可行整数解， 状态如实为 feasible，不作最优宣称。缓解方向（M4 前决策）：提高单窗 时限、跨窗热启动、或沿用 M2.5 的认证夹逼管线加大下界模型预算。
  • 试点标定（S-balance 归一化前置）：50 组随机配置 × 5 代表周 = 250 个 可行周值（0 个求解不可行；86 次抽样中 36 次被面积约束预筛），产出 cost ∈ [181,652, 648,937] 元、carbon ∈ [−26,450, 37,665] kgCO₂、 heat ∈ [26,659, 105,102] 元 的 min-max 范围。experiments.yaml 未被 修改（哈希冻结是验收不变量），建议回填值存于 outputs/m3_validation/normalization_pilot.json，是否回填由规格方裁定。

5. 潜在应用价值

  1. 寒区储能选址/选型的决策模型：可用带机制 + 面积硬约束（钠电占地约为 锂电 1.8 倍）+ 温度相关倍率，是寒冷地区 BESS"钠 vs 锂"比较中常被忽略的 一组约束；模型把它们全部显式化并经过四角参数验算。
  2. 算力柔性参与电网调节的量化框架：S1/S2/S3 任务分级 + Token 口径 SLA + 电量口径柔性比，把"训练任务可迁移"这一直觉变成可优化的 MILP 变量， 可直接用于 AIDC 参与需求响应/现货市场的容量评估。
  3. 可迁移的认证求解管线：规则排程上界 + 物理下界模型的夹逼方法（§4.3）， 适用于任何含大量整数排程变量的能源调度问题——在不牺牲精度证据的前提下 把求解时间压到外层优化的可承受范围。
  4. 多保真外层优化框架：代表周低保真（f3=NaN、三维搜索）→ 高保真确认 （滚动年度 + 故障场景补齐 f3、重新过滤 Pareto）的流程已打通，BO 与网格 对照将按同一预算如实报告谁更好。
  5. 负结果与删失处理方法学：右删失显式建模、平坦目标检测、"不得预设结论" 的验收条款写法，对任何仿真驱动的选型研究都可复用。

6. 当前限制（如实列出）

  • 参数标定未完成：chemistry_*.yaml 与 generation.yaml 中仍有 PLACEHOLDER （成本、寿命、占地、气象分布）；5 处容量因子与画像参考值偏差超容差 （如 R-NW 光伏模拟 CF 0.104 vs 参考 0.170），按规格不回调参数拟合参考值， 留待公开气象数据标定。
  • f3 当前不可用于排序（§4.2）；启用前需延长故障时长或完成标定。 另：A0 因未安装储能跳过 N-1 族场景，跨臂 f3 的孤岛项加权分母不同， 比较时须用共同场景集口径（M2.5 的 common_non_n1 聚合）。
  • 高保真 A3 在 60 s/窗预算下仅 20/52 窗认证 gap ≤1%（中位 1.17%、 最大 5.45%，全部为经验证的可行整数解）；M4 矩阵前需决策单窗时限或 热启动策略，不得以放宽 gap 口径替代。
  • Na/Li 的 cycle_life.mid 目前都是 6000："先放钠电以保护更长寿命体系"的 降级规则理由暂不成立，标定后须复核。
  • 代表周任务为窗口内重新生成（规格规定），不是年度任务表的映射——跨周任务 迁移不在模型范围内。
  • 高保真年度求解为滚动时域（中间窗口终端自由），存在滚动规划固有的近视性； 每窗可行性单独检查并记录，首个不可行窗即停止并原样上报。

7. 报告口径（按 SPEC §10 红线固定）

口径定义
SLA 达标率1 − 延迟任务数/总任务数，分母只含 S2/S3（S1 并入基础负荷）
关键负荷设施侧单一比例 critical_it_fraction = 0.60 × IT_base；判据 shed_c ≤ 1e-9 MW
降级策略故障模式四级字典序（关键削减 → 非关键削减 → 供热缺口 → 化学体系偏好），T_island 是降级策略下的续航
等效循环Σ max(0, soc(t) − soc(t+1))（SOC 掉幅口径，不除以 2）
LCOE(年化运营成本 + 储能年化成本) / 年服务电量；运营成本含锅炉燃料与弃电惩罚、不含 SLA 惩罚
碳账电网因子 + 锅炉燃气 − 工业供热替代；制造碳按项目寿命摊销计入；碳资产不入账
供热表述余热为低品位热源，经热泵提温后作为市政/工业供热补充；锅炉是主要（峰值）热源
弃热Q_dump（输出侧）与 Q_unused（热源侧）分别报告，不合并

8. 可复现性

python tools/selfcheck.py                                  # 规格一致性自查
python -m pytest tests -q                                  # 268 项验收测试
python -B outputs/m3_validation/validate_m3_low.py         # M3 低保真闭环验证（~6 min）
python -B outputs/m3_validation/validate_m3_high.py A3     # M3 高保真闭环验证（~51 min）
python -B outputs/m3_validation/validate_m3_high.py A0     # M3 高保真 A0 退化臂（~21 min）
python -B outputs/m3_validation/pilot_normalization.py 5   # S-balance 归一化试点（~49 min，5 并行）

产物索引：outputs/m1_validation/（参数层）、outputs/m2_v29_validation/ （内层 MILP，190 项）、outputs/m25_validation/（韧性 + 提速，251 项， report.json 含逐场景删失与认证详情）、outputs/m3_validation/（外层闭环）。 配置冻结：11 份 YAML 的 SHA-256 与 M1 基线一致，任何验证脚本运行前后各核对一次。

本报告由项目实现方撰写；规格方与实现方分离，所有验收测试由规格方预先定义， 实现方不得修改期望值。历史决策链见 CHANGELOG.md（v1.1–v2.9 逐条记录 "哪条结论被哪条取代"）。

登录后可以留言