Pyomo电池储能优化：特定时段禁充约束的异常问题排查

核心问题分析 #

变量全为0或None的本质是模型不可行，大概率是禁充约束的构建逻辑错误，或者该约束与原有核心约束（如SOC首尾相等、功率平衡）产生了冲突。以下是具体排查和修正方案：

1. 修正禁充约束的构建逻辑（避免索引匹配错误） #

别用numpy生成索引后直接匹配，直接在Pyomo约束规则中对每个时间点做判断，从根源避免索引类型不匹配、时段范围错误的问题：

def no_charge_rule(model, t):
    # 假设t为datetime类型，直接提取小时判断
    hour = t.hour
    # 8:00-15:00禁止充电，注意边界：若15:00允许充电则用8<=hour<15，若15:00也禁则用8<=hour<=15
    return model.Charge[t] == 0 if 8 <= hour < 15 else Constraint.Skip

# 给所有时间点绑定约束规则，符合条件的自动生成禁充约束
model.NoCharge = Constraint(model.t, rule=no_charge_rule)

如果你的时间集t是整数索引（如0代表0点，1代表1点...），可以这样计算禁充索引：

total_days = len(model.t) // 24
no_charge_indices = []
for day in range(total_days):
    # 对应每日8点到14点的索引（步长1小时）
    no_charge_indices.extend(range(day*24 +8, day*24 +15))

# 仅对目标时段添加禁充约束
model.NoCharge = Constraint(
    [t for t in model.t if t in no_charge_indices],
    rule=lambda model, t: model.Charge[t] == 0
)

2. 验证模型可行性 #

如果修正约束后仍无有效解，说明禁充规则与原有核心约束冲突：

• 临时移除SOC首尾相等约束，重新求解，看是否能得到非零结果。如果可以，说明原约束要求电池必须在禁充时段充电才能满足首尾SOC一致，需要调整SOC约束的松弛度或目标函数权重。
• 打印所有禁充约束，确认仅对目标时段生效：

  for idx, con in model.NoCharge.items():
      print(f"时段{idx}约束：{con}")
  

3. 排查求解器输出 #

运行求解器时开启日志，查看是否有明确的不可行提示：

solver = SolverFactory('glpk')  # 或你使用的求解器
result = solver.solve(model, tee=True)
print(result)

如果日志显示infeasible，说明需要调整约束逻辑（比如放宽SOC上下限、允许部分时段的功率约束松弛）。

内容的提问来源于stack exchange，提问作者MiPre