Add Contingency Percentage parameter to Economics. (TODO to incorporate/consolidate SBTEconomics)

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -1044,6 +1044,21 @@ def __init__(self, model: Model):
             ErrMessage="assume default inflation rate during construction (0)",
             ToolTipText='For SAM Economic Models, this value is treated as an indirect EPC capital cost percentage.'
         )
+
+        self.contingency_percentage = self.ParameterDict[self.contingency_percentage.Name] = floatParameter(
+            'Contingency Percentage',
+            DefaultValue=15.,
+            Min=0.,
+            Max=100.,
+            UnitType=Units.PERCENT,
+            PreferredUnits=PercentUnit.PERCENT,
+            CurrentUnits=PercentUnit.PERCENT,
+            ToolTipText='The contingency percentage applied to the direct capital costs for stimulation, '
+                        'field gathering system, exploration, and surface plant. '
+                        '(Note: well drilling and completion costs do not have contingency applied and are not '
+                        'affected by this parameter.)'
+        )
+
         self.wellcorrelation = self.ParameterDict[self.wellcorrelation.Name] = intParameter(
             "Well Drilling Cost Correlation",
             DefaultValue=WellDrillingCostCorrelation.VERTICAL_LARGE_INT1.int_value,
@@ -2366,8 +2381,8 @@ def Calculate(self, model: Model) -> None:
             if self.ccexplfixed.Valid:
                 self.Cexpl.value = self.ccexplfixed.value
             else:
-                self.Cexpl.value = 1.15 * self.ccexpladjfactor.value * self._indirect_cost_factor * (
-                    1. + self.cost_one_production_well.value * 0.6)  # 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+                self.Cexpl.value = self._contingency_factor * self.ccexpladjfactor.value * self._indirect_cost_factor * (
+                    1. + self.cost_one_production_well.value * 0.6)
 
             # Surface Piping Length Costs (M$) #assumed $750k/km
             self.Cpiping.value = 750 / 1000 * model.surfaceplant.piping_length.value
@@ -2635,6 +2650,10 @@ def _wellfield_indirect_cost_factor(self) -> float:
     def _stimulation_indirect_cost_factor(self) -> float:
         return 1 + self.stimulation_indirect_capital_cost_percentage.quantity().to('dimensionless').magnitude
 
+    @property
+    def _contingency_factor(self) -> float:
+        return 1 + self.contingency_percentage.quantity().to('dimensionless').magnitude
+
     def calculate_wellfield_costs(self, model: Model) -> None:
         if self.per_production_well_cost.Valid:
             self.cost_one_production_well.value = self.per_production_well_cost.value
@@ -2727,7 +2746,7 @@ def calculate_stimulation_costs(self, model: Model) -> PlainQuantity:
                 )
                 * self.ccstimadjfactor.value
                 * self._stimulation_indirect_cost_factor
-                * 1.15  # 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+                * self._contingency_factor
             )
 
         return quantity(stimulation_costs, self.Cstim.CurrentUnits)
@@ -2766,8 +2785,8 @@ def calculate_field_gathering_costs(self, model: Model) -> None:
                             1750 * injpumphpcorrected ** 0.7) * 3 * injpumphpcorrected ** (-0.11)
                 self.Cpumps = Cpumpsinj + Cpumpsprod
 
-            # Based on GETEM 2016: 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
-            self.Cgath.value = 1.15 * self.ccgathadjfactor.value * self._indirect_cost_factor * (
+            # Based on GETEM 2016
+            self.Cgath.value = self._contingency_factor * self.ccgathadjfactor.value * self._indirect_cost_factor * (
                     (model.wellbores.nprod.value + model.wellbores.ninj.value) * 750 * 500. + self.Cpumps) / 1E6
 
     def calculate_plant_costs(self, model: Model) -> None:
@@ -2777,21 +2796,32 @@ def calculate_plant_costs(self, model: Model) -> None:
             if self.ccplantfixed.Valid:
                 self.Cplant.value = self.ccplantfixed.value
             else:
-                # 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
-                self.Cplant.value = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
-                    model.surfaceplant.HeatExtracted.value) * 1000.
+                self.Cplant.value = (self._indirect_cost_factor
+                                     * self._contingency_factor
+                                     * self.ccplantadjfactor.value
+                                     * 250E-6
+                                     * np.max(model.surfaceplant.HeatExtracted.value)
+                                     * 1000.)
 
         # absorption chiller
         elif model.surfaceplant.enduse_option.value == EndUseOptions.HEAT and model.surfaceplant.plant_type.value == PlantType.ABSORPTION_CHILLER:  # absorption chiller
             if self.ccplantfixed.Valid:
                 self.Cplant.value = self.ccplantfixed.value
             else:
                 # this is for the direct-use part all the way up to the absorption chiller
-                self.Cplant.value = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
-                    model.surfaceplant.HeatExtracted.value) * 1000.  # 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+                self.Cplant.value = (self._indirect_cost_factor
+                                     * self._contingency_factor
+                                     * self.ccplantadjfactor.value
+                                     * 250E-6
+                                     * np.max(model.surfaceplant.HeatExtracted.value)
+                                     * 1000.)
                 if self.chillercapex.value == -1:  # no value provided by user, use built-in correlation ($2500/ton)
-                    self.chillercapex.value = self._indirect_cost_factor * 1.15 * np.max(
-                        model.surfaceplant.cooling_produced.value) * 1000 / 3.517 * 2500 / 1e6  # $2,500/ton of cooling. 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+                    self.chillercapex.value = (
+                        self._indirect_cost_factor
+                        * self._contingency_factor
+                        * np.max(model.surfaceplant.cooling_produced.value)
+                        * 1000 / 3.517 * 2500 / 1e6 # $2,500/ton of cooling.
+                    )
 
                 # now add chiller cost to surface plant cost
                 self.Cplant.value += self.chillercapex.value
@@ -2802,11 +2832,11 @@ def calculate_plant_costs(self, model: Model) -> None:
                 self.Cplant.value = self.ccplantfixed.value
             else:
                 # this is for the direct-use part all the way up to the heat pump
-                self.Cplant.value = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
-                    model.surfaceplant.HeatExtracted.value) * 1000.  # 1.15 for 15% contingency
+                self.Cplant.value = self._indirect_cost_factor * self._contingency_factor * self.ccplantadjfactor.value * 250E-6 * np.max(
+                    model.surfaceplant.HeatExtracted.value) * 1000.
                 if self.heatpumpcapex.value == -1:  # no value provided by user, use built-in correlation ($150/kWth)
-                    self.heatpumpcapex.value = self._indirect_cost_factor * 1.15 * np.max(
-                        model.surfaceplant.HeatProduced.value) * 1000 * 150 / 1e6  # $150/kW. 1.15 for 15% contingency  TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+                    self.heatpumpcapex.value = self._indirect_cost_factor * self._contingency_factor * np.max(
+                        model.surfaceplant.HeatProduced.value) * 1000 * 150 / 1e6  # $150/kW - TODO parameterize
 
                 # now add heat pump cost to surface plant cost
                 self.Cplant.value += self.heatpumpcapex.value
@@ -2816,8 +2846,7 @@ def calculate_plant_costs(self, model: Model) -> None:
             if self.ccplantfixed.Valid:
                 self.Cplant.value = self.ccplantfixed.value
             else:
-                # 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
-                self.Cplant.value = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
+                self.Cplant.value = self._indirect_cost_factor * self._contingency_factor * self.ccplantadjfactor.value * 250E-6 * np.max(
                     model.surfaceplant.HeatExtracted.value) * 1000.
 
                 # add 65$/KW for peaking boiler
@@ -2986,23 +3015,28 @@ def calculate_plant_costs(self, model: Model) -> None:
                     # factor 1.10 to convert from 2016 to 2022
                     direct_plant_cost_MUSD = self.ccplantadjfactor.value * self.Cplantcorrelation * 1.02 * 1.10
 
-                # factor 1.15 for 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
-                self.Cplant.value = self._indirect_cost_factor * 1.15 * direct_plant_cost_MUSD
+                self.Cplant.value = self._indirect_cost_factor * self._contingency_factor * direct_plant_cost_MUSD
                 self.CAPEX_cost_electricity_plant = self.Cplant.value
 
         # add direct-use plant cost of co-gen system to Cplant (only of no total Cplant was provided)
-        if not self.ccplantfixed.Valid:  # 1.15 below for contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+        if not self.ccplantfixed.Valid:
             if model.surfaceplant.enduse_option.value in [EndUseOptions.COGENERATION_TOPPING_EXTRA_ELECTRICITY,
                                                           EndUseOptions.COGENERATION_TOPPING_EXTRA_HEAT]:  # enduse_option = 3: cogen topping cycle
-                self.CAPEX_cost_heat_plant = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
-                    model.surfaceplant.HeatProduced.value / model.surfaceplant.enduse_efficiency_factor.value) * 1000.
+                self.CAPEX_cost_heat_plant = (
+                    self._indirect_cost_factor
+                    * self._contingency_factor
+                    * self.ccplantadjfactor.value
+                    * 250E-6
+                    * np.max(model.surfaceplant.HeatProduced.value / model.surfaceplant.enduse_efficiency_factor.value)
+                    * 1000.
+                )
             elif model.surfaceplant.enduse_option.value in [EndUseOptions.COGENERATION_BOTTOMING_EXTRA_HEAT,
                                                             EndUseOptions.COGENERATION_BOTTOMING_EXTRA_ELECTRICITY]:  # enduse_option = 4: cogen bottoming cycle
-                self.CAPEX_cost_heat_plant = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
+                self.CAPEX_cost_heat_plant = self._indirect_cost_factor * self._contingency_factor * self.ccplantadjfactor.value * 250E-6 * np.max(
                     model.surfaceplant.HeatProduced.value / model.surfaceplant.enduse_efficiency_factor.value) * 1000.
             elif model.surfaceplant.enduse_option.value in [EndUseOptions.COGENERATION_PARALLEL_EXTRA_ELECTRICITY,
                                                             EndUseOptions.COGENERATION_PARALLEL_EXTRA_HEAT]:  # cogen parallel cycle
-                self.CAPEX_cost_heat_plant = self._indirect_cost_factor * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
+                self.CAPEX_cost_heat_plant = self._indirect_cost_factor * self._contingency_factor * self.ccplantadjfactor.value * 250E-6 * np.max(
                     model.surfaceplant.HeatProduced.value / model.surfaceplant.enduse_efficiency_factor.value) * 1000.
 
             self.Cplant.value = self.Cplant.value + self.CAPEX_cost_heat_plant

src/geophires_x_schema_generator/geophires-request.json

@@ -1719,6 +1719,15 @@
       "minimum": 0.0,
       "maximum": 1.0
     },
+    "Contingency Percentage": {
+      "description": "The contingency percentage applied to the direct capital costs for stimulation, field gathering system, exploration, and surface plant. (Note: well drilling and completion costs do not have contingency applied and are not affected by this parameter.)",
+      "type": "number",
+      "units": "%",
+      "category": "Economics",
+      "default": 15.0,
+      "minimum": 0.0,
+      "maximum": 100.0
+    },
     "Well Drilling Cost Correlation": {
       "description": "Select the built-in well drilling and completion cost correlation: 1: vertical small diameter, baseline; 2: deviated small diameter, baseline; 3: vertical large diameter, baseline; 4: deviated large diameter, baseline; 5: Simple (per-meter cost); 6: vertical small diameter, intermediate1; 7: vertical small diameter, intermediate2; 8: deviated small diameter, intermediate1; 9: deviated small diameter, intermediate2; 10: vertical large diameter, intermediate1; 11: vertical large diameter, intermediate2; 12: deviated large diameter, intermediate1; 13: deviated large diameter, intermediate2; 14: vertical open-hole, small diameter, ideal; 15: deviated liner, small diameter, ideal; 16: vertical open-hole, large diameter, ideal; 17: deviated liner, large diameter, ideal. Baseline correlations (1-4) are from NREL's 2025 cost curve update. Intermediate and ideal correlations (6-17) are from GeoVision.",
       "type": "integer",