Record THP limit re-specification in SimulatorUpdate

The new thp_respec_wells set carries the WELL_THP_UPDATE information
through ACTIONX processing, where report-step events do not reach the
simulator.

Author: GitPaean

opm/input/eclipse/Schedule/Action/SimulatorUpdate.hpp

@@ -44,6 +44,7 @@ struct SimulatorUpdate
         simulatorUpdate.well_structure_changed = true;
         simulatorUpdate.affected_wells = {"test"};
         simulatorUpdate.welpi_wells.insert("I-45");
+        simulatorUpdate.thp_respec_wells.insert("P-3");
         simulatorUpdate.new_frac_wconns.assign({
                 std::pair { "I-45"s, std::vector { std::size_t{11}, std::size_t{22}, std::size_t{33} } },
                 std::pair { "RA-MAN"s, std::vector { std::size_t{1}, std::size_t{7}, std::size_t{29} } },
@@ -57,6 +58,7 @@ struct SimulatorUpdate
     {
         serializer(affected_wells);
         serializer(welpi_wells);
+        serializer(thp_respec_wells);
         serializer(new_frac_wconns);
         serializer(tran_update);
         serializer(well_structure_changed);
@@ -70,6 +72,15 @@ struct SimulatorUpdate
     /// its internal notion of the connection transmissibility factors.
     std::unordered_set<std::string> welpi_wells{};
 
+    /// Wells whose THP limit and/or VFP table has been (re)specified by
+    /// WCONPROD or WCONHIST for the named wells, by WELTARG with THP or VFP
+    /// control, or by WTMULT with THP control. Unlike affected_wells, this
+    /// set also includes wells for which the entered values are unchanged,
+    /// since re-specifying either input cancels a THP limit imposed
+    /// dynamically by the simulator (e.g. by network balancing). See also
+    /// the WELL_THP_UPDATE schedule event.
+    std::unordered_set<std::string> thp_respec_wells{};
+
     /// New well connections created as a result of a geomechanical
     /// fracturing process.
     ///
@@ -116,6 +127,9 @@ struct SimulatorUpdate
         this->welpi_wells.insert(otherSimUpdate.welpi_wells.begin(),
                                  otherSimUpdate.welpi_wells.end());
 
+        this->thp_respec_wells.insert(otherSimUpdate.thp_respec_wells.begin(),
+                                      otherSimUpdate.thp_respec_wells.end());
+
         this->new_frac_wconns.insert(this->new_frac_wconns.end(),
                                      otherSimUpdate.new_frac_wconns.begin(),
                                      otherSimUpdate.new_frac_wconns.end());
@@ -128,6 +142,7 @@ struct SimulatorUpdate
         this->well_structure_changed = false;
         this->affected_wells.clear();
         this->welpi_wells.clear();
+        this->thp_respec_wells.clear();
         this->new_frac_wconns.clear();
     }
 
@@ -137,6 +152,7 @@ struct SimulatorUpdate
             && (this->well_structure_changed == that.well_structure_changed)
             && (this->affected_wells == that.affected_wells)
             && (this->welpi_wells == that.welpi_wells)
+            && (this->thp_respec_wells == that.thp_respec_wells)
             && (this->new_frac_wconns == that.new_frac_wconns)
             ;
     }

opm/input/eclipse/Schedule/HandlerContext.cpp

@@ -48,6 +48,13 @@ void HandlerContext::affected_well(const std::string& well_name)
     }
 }
 
+void HandlerContext::thp_respec_well(const std::string& well_name)
+{
+    if (sim_update != nullptr) {
+        sim_update->thp_respec_wells.insert(well_name);
+    }
+}
+
 void HandlerContext::welpi_well(const std::string& well_name)
 {
     if (sim_update != nullptr) {

opm/input/eclipse/Schedule/HandlerContext.hpp

@@ -107,6 +107,10 @@ class HandlerContext
     //! \brief Mark that a well is affected by WELPI.
     void welpi_well(const std::string& well_name);
 
+    //! \brief Mark that the THP limit and/or VFP table of a well has been
+    //! (re)specified, also when the entered values are unchanged.
+    void thp_respec_well(const std::string& well_name);
+
     //! \brief Mark that transmissibilities must be recalculated.
     void record_tran_change();
 

opm/input/eclipse/Schedule/Well/WellKeywordHandlers.cpp

@@ -192,6 +192,7 @@ void handleWCONHIST(HandlerContext& handlerContext)
             // whether or not any value changed.
             handlerContext.state().events().addEvent( ScheduleEvents::WELL_THP_UPDATE );
             handlerContext.state().wellgroup_events().addEvent( well2.name(), ScheduleEvents::WELL_THP_UPDATE);
+            handlerContext.thp_respec_well(well2.name());
 
             if (update_well) {
                 handlerContext.state().events().addEvent( ScheduleEvents::PRODUCTION_UPDATE );
@@ -446,6 +447,7 @@ void handleWCONPROD(HandlerContext& handlerContext)
             // by every WCONPROD, whether or not any value changed.
             handlerContext.state().events().addEvent( ScheduleEvents::WELL_THP_UPDATE );
             handlerContext.state().wellgroup_events().addEvent( well2.name(), ScheduleEvents::WELL_THP_UPDATE);
+            handlerContext.thp_respec_well(well2.name());
 
             if (update_well) {
                 handlerContext.state().events().addEvent( ScheduleEvents::PRODUCTION_UPDATE );
@@ -851,6 +853,7 @@ void handleWELTARG(HandlerContext& handlerContext)
             {
                 handlerContext.state().events().addEvent( ScheduleEvents::WELL_THP_UPDATE );
                 handlerContext.state().wellgroup_events().addEvent( well_name, ScheduleEvents::WELL_THP_UPDATE);
+                handlerContext.thp_respec_well(well_name);
             }
 
             if (update) {

opm/input/eclipse/Schedule/Well/WellPropertiesKeywordHandlers.cpp

@@ -642,6 +642,7 @@ void handleWTMULT(HandlerContext& handlerContext)
                     handlerContext.state().wellgroup_events()
                         .addEvent(well_name, ScheduleEvents::WELL_THP_UPDATE);
 
+                    handlerContext.thp_respec_well(well_name);
                 }
 
                 if (update_well) {