Pyomo可再生能源最优容量选择优化模型问题求助

问题描述 #

• 开发Pyomo优化模型，输入客户用电负荷与多个可再生能源项目发电数据，目标是求解最低成本的可再生能源项目组合，以最小化电力消耗
• 编写第一版代码后遇到两个关键问题：
  1. d_spill_var规则执行失败，触发PyomoException：无法转换非恒定Pyomo表达式，推测根源是在依赖表达式中使用了max(expr, 0)
  2. 运行时无法生成合法LP文件，原因是目标函数包含非二次非线性项

核心解决思路 #

针对两个问题的针对性处理方向：

1. 线性化max(expr, 0)表达式：Pyomo的线性求解器不支持直接解析非线性的max函数，需通过引入辅助变量将其转化为线性约束。例如定义辅助变量spill，添加约束spill >= expr和spill >= 0，用该变量替代原max(expr, 0)的位置
2. 修正目标函数非线性问题：检查目标函数中的非线性项，将非线性逻辑拆分出来定义为独立的线性约束，确保目标函数仅保留线性或二次项，避免求解器无法生成LP文件

最终可行方案 #

通过以下操作彻底解决问题，更新后的代码可正常运行：

• 对max约束进行线性化处理，用辅助变量替代原非线性表达式
• 将依赖表达式重定义为显式模型约束，不再嵌入变量规则中
• 移除模型中的冗余变量，简化模型结构
• 将最后两个约束修改为求和形式，确保所有约束均为线性格式

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