Infeasibility analyzer : HydroPower constraint (#2227) [ANT-1825]

Add constraint that can possibily lead to infeasibilities. HydroPower is the constraint on the sum of generated hydro energy over the week for a given area. The provided RHS in some cases makes the problem infeasible. --------- Co-authored-by: Florian OMNES <[email protected]>
AntaresSimulatorTeam · Jul 8, 2024 · 0742047 · 0742047
1 parent e9d7c45
commit 0742047
Show file tree

Hide file tree

Showing 4 changed files with 14 additions and 1 deletion.
diff --git a/src/solver/infeasible-problem-analysis/constraint.cpp b/src/solver/infeasible-problem-analysis/constraint.cpp
@@ -26,6 +26,7 @@
 #include <sstream>
 
 #include <boost/algorithm/string/regex.hpp>
+#include <boost/regex.hpp>
 
 namespace
 {
@@ -143,6 +144,10 @@ ConstraintType Constraint::getType() const
     {
         return ConstraintType::hydro_reservoir_level;
     }
+    if (nameComponents_.at(0) == "HydroPower")
+    {
+        return ConstraintType::hydro_production_weekly;
+    }
     if (nameComponents_.at(0) == "Level")
     {
         return ConstraintType::short_term_storage_level;
@@ -194,6 +199,8 @@ std::string Constraint::prettyPrint() const
     case ConstraintType::hydro_reservoir_level:
         return "Hydro reservoir constraint at area '" + getAreaName() + "' at hour "
                + getTimeStepInYear();
+    case ConstraintType::hydro_production_weekly:
+        return "Hydro weekly production at area '" + getAreaName() + "'";
     case ConstraintType::short_term_storage_level:
         return "Short-term-storage reservoir constraint at area '" + getAreaName() + "' in STS '"
                + getSTSName() + "' at hour " + getTimeStepInYear();

diff --git a/...m-analysis/include/antares/solver/infeasible-problem-analysis/constraint-slack-analysis.h b/...m-analysis/include/antares/solver/infeasible-problem-analysis/constraint-slack-analysis.h
@@ -57,7 +57,8 @@ class ConstraintSlackAnalysis: public UnfeasibilityAnalysis
 
     std::vector<const operations_research::MPVariable*> slackVariables_;
     const std::string constraint_name_pattern = "^AreaHydroLevel::|::hourly::|::daily::|::weekly::|"
-                                                "^FictiveLoads::|^Level::";
+                                                "^FictiveLoads::|^Level::|"
+                                                "^HydroPower::";
 };
 
 } // namespace Antares::Optimization
diff --git a/...feasible-problem-analysis/include/antares/solver/infeasible-problem-analysis/constraint.h b/...feasible-problem-analysis/include/antares/solver/infeasible-problem-analysis/constraint.h
@@ -32,6 +32,7 @@ enum class ConstraintType
     binding_constraint_weekly,
     fictitious_load,
     hydro_reservoir_level,
+    hydro_production_weekly,
     short_term_storage_level,
     none
 };

diff --git a/src/solver/infeasible-problem-analysis/report.cpp b/src/solver/infeasible-problem-analysis/report.cpp
@@ -92,6 +92,10 @@ void InfeasibleProblemReport::logSuspiciousConstraints()
         Antares::logs.error() << "* Hydro reservoir impossible to manage with cumulative options "
                                  "\"hard bounds without heuristic\"";
     }
+    if (nbConstraintsByType_[ConstraintType::hydro_production_weekly] > 0)
+    {
+        Antares::logs.error() << "* impossible to generate exactly the weekly hydro target";
+    }
     if (nbConstraintsByType_[ConstraintType::fictitious_load] > 0)
     {
         Antares::logs.error() << "* Last resort shedding status,";