水力机组组合优化模型不可行：Gurobi/IPOPT求解差异及可行化咨询

我正在复现某学术论文的水力机组组合优化问题，已核对所有问题数据与约束，确认无误。目标是使用Gurobi求解器，但该求解器判定模型不可行无法求解；改用IPOPT求解器后，虽同样提示模型不可行，却得到了正确的目标函数值。我已通过Infeasible Irreducible System (IIS)工具定位到导致模型不可行的约束，相关结果如下：

Model unknown_copy
\ LP format - for model browsing. Use MPS format to capture full model detail.
Minimize
 
Subject To
 constraints(1474): - 470 u(5_2_17) + w(5_2_17) <= 0
 constraints(1522): - 470 u(5_3_17) + w(5_3_17) <= 0
 constraints(1570): - 470 u(5_4_17) + w(5_4_17) <= 0
 constraints(2290): g(5_1_17) - 215 u(5_1_17) <= 0
 constraints(2561): g(5_1_17) + g(5_2_17) + g(5_3_17) + g(5_4_17) = 830
 constraints(9042): g(5_2_17) + 37.5685717831504 u(5_2_17)
   - 0.461258179047618 w(5_2_17) <= 30.60699513446293
 constraints(9138): g(5_3_17) + 37.5685717831504 u(5_3_17)
   - 0.461258179047618 w(5_3_17) <= 30.60699513446293
 constraints(9140): g(5_3_17) + 260.913290772527 u(5_3_17)
   - 0.197897437142857 w(5_3_17) <= 366.4855962350004
 constraints(9234): g(5_4_17) + 37.5685717831504 u(5_4_17)
   - 0.461258179047618 w(5_4_17) <= 30.60699513446293
 constraints(9236): g(5_4_17) + 260.913290772527 u(5_4_17)
   - 0.197897437142857 w(5_4_17) <= 366.4855962350004
Bounds
 g(5_1_17) free
 g(5_2_17) free
 g(5_3_17) free
 g(5_4_17) free
 w(5_2_17) free
 w(5_3_17) free
 w(5_4_17) free
Binaries
 u(5_1_17) u(5_2_17) u(5_3_17) u(5_4_17)
End here

可行化处理方案 #

1. 约束松弛策略 #

针对IIS中的约束按优先级松弛：

• 等式约束优先松弛：对约束(2561)的功率平衡等式，引入非负松弛变量s_low和s_high，将原等式改为830 - s_low ≤ g(5_1_17) + g(5_2_17) + g(5_3_17) + g(5_4_17) ≤ 830 + s_high，同时在目标函数中加入惩罚项λ*(s_low + s_high)，λ取远大于目标函数系数的值（如1e6），保证仅在必要时松弛，且解贴近原问题要求。
• 出力类约束松弛：对约束(2290)、(9042)等出力关联约束，引入单侧非负松弛变量，例如将g(5_1_17) - 215u(5_1_17) ≤ 0改为g(5_1_17) - 215u(5_1_17) ≤ s1，目标函数添加惩罚项μ*s1，μ根据出力偏差的经济成本设定。
• 启停关联约束松弛：对约束(1474)等w与u的关联约束，可引入小范围松弛，如-470u(5_2_17) + w(5_2_17) ≤ s2，s2≥0，惩罚项系数参考机组启停成本设置。

2. Gurobi内置工具应用 #

• 调用Model.computeIIS()再次确认最小不可行约束集，确保没有遗漏关键约束。
• 适当调整FeasibilityTol参数（如从默认1e-6调整至1e-4），放宽可行性容忍度，此方法仅作临时调试，需注意解的精度损失。
• 使用FeasRelax()函数自动生成松弛模型，可指定仅松弛等式约束或特定类型约束，工具会自动添加合理惩罚项。

3. 变量边界优化 #

检查g、w等自由变量的物理意义，若实际场景中存在隐含边界（如机组出力g不能为负），补充g(...) ≥ 0这类边界约束，避免变量取值极端导致约束冲突。

4. 反向验证IPOPT结果 #

将IPOPT给出的目标值对应的变量取值，代入IIS约束计算每个约束的违反量，优先松弛违反量最大的约束，针对性调整松弛系数，确保Gurobi求解结果贴近IPOPT的可行解。

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