Parameterize Reservoir Stimulation Indirect Capital Cost Percentage

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -629,6 +629,20 @@ def __init__(self, model: Model):
             Valid=True,
             ToolTipText="Multiplier for reservoir stimulation capital cost correlation"
         )
+        self.stimulation_indirect_capital_cost = \
+          self.ParameterDict[self.stimulation_indirect_capital_cost.Name] = floatParameter(
+            'Reservoir Stimulation Indirect Capital Cost Percentage',
+            DefaultValue=5,
+            Min=0,
+            Max=100,
+            UnitType=Units.PERCENT,
+            PreferredUnits=PercentUnit.PERCENT,
+            CurrentUnits=PercentUnit.PERCENT,
+            ToolTipText=f'The indirect capital cost for reservoir stimulation, '
+                        f'calculated as a percentage of the direct cost. '
+                        f'(Not applied if {self.ccstimfixed.Name} is provided.)'
+        )
+
         self.ccexplfixed = self.ParameterDict[self.ccexplfixed.Name] = floatParameter(
             "Exploration Capital Cost",
             DefaultValue=-1.0,
@@ -1638,8 +1652,11 @@ def __init__(self, model: Model):
             CurrentUnits=EnergyCostUnit.DOLLARSPERMMBTU
         )
 
-        # TODO https://github.com/NREL/GEOPHIRES-X/issues/383?title=Parameterize+indirect+cost+factor
-        stimulation_contingency_and_indirect_costs_tooltip = 'plus 15% contingency plus 5% indirect costs'
+        stimulation_contingency_and_indirect_costs_tooltip = (
+            f'plus 15% contingency ' # TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+            f'plus {self.stimulation_indirect_capital_cost.quantity().to(convertible_unit("%")).magnitude}% '
+            f'indirect costs'
+        )
 
         # noinspection SpellCheckingInspection
         self.Cstim = self.OutputParameterDict[self.Cstim.Name] = OutputParameter(
@@ -2382,14 +2399,17 @@ def Calculate(self, model: Model) -> None:
             stim_cost_per_production_well = self.stimulation_cost_per_production_well.quantity().to(
                 self.Cstim.CurrentUnits).magnitude
 
-            # 1.15 for 15% contingency and 1.05 for 5% indirect costs
-            # TODO https://github.com/NREL/GEOPHIRES-X/issues/383?title=Parameterize+indirect+cost+factor
-            self.Cstim.value = ((
-                                 stim_cost_per_injection_well * model.wellbores.ninj.value
-                                 + stim_cost_per_production_well * model.wellbores.nprod.value
-                                )
-                                * self.ccstimadjfactor.value
-                                * 1.05 * 1.15)
+            stimulation_indirect_cost_fraction = (self.stimulation_indirect_capital_cost.quantity()
+                                                  .to('dimensionless').magnitude)
+            self.Cstim.value = (
+                (
+                    stim_cost_per_injection_well * model.wellbores.ninj.value
+                    + stim_cost_per_production_well * model.wellbores.nprod.value
+                )
+                * self.ccstimadjfactor.value
+                * (1 + stimulation_indirect_cost_fraction)
+                * 1.15  # 15% contingency TODO https://github.com/NREL/GEOPHIRES-X/issues/383
+            )
 
         # field gathering system costs (M$)
         if self.ccgathfixed.Valid:

src/geophires_x_schema_generator/geophires-request.json

@@ -1422,6 +1422,15 @@
       "minimum": 0,
       "maximum": 10
     },
+    "Reservoir Stimulation Indirect Capital Cost Percentage": {
+      "description": "The indirect capital cost for reservoir stimulation, calculated as a percentage of the direct cost. (Not applied if Reservoir Stimulation Capital Cost is provided.)",
+      "type": "number",
+      "units": "%",
+      "category": "Economics",
+      "default": 5,
+      "minimum": 0,
+      "maximum": 100
+    },
     "Exploration Capital Cost": {
       "description": "Total exploration capital cost",
       "type": "number",

tests/test_geophires_x.py

@@ -975,3 +975,31 @@ def _get_result(prod_well_stim_MUSD: Optional[int] = None) -> GeophiresXResult:
             result_prod_stim.result['CAPITAL COSTS (M$)']['Stimulation costs']['value'],
             places=1,
         )
+
+    def test_stimulation_indirect_cost(self):
+        def _get_result(indirect_cost_percent: Optional[int] = None) -> float:
+            p = {}
+            if indirect_cost_percent is not None:
+                p['Reservoir Stimulation Indirect Capital Cost Percentage'] = indirect_cost_percent
+
+            return (
+                GeophiresXClient()
+                .get_geophires_result(
+                    ImmutableGeophiresInputParameters(
+                        from_file_path=self._get_test_file_path('geophires_x_tests/generic-egs-case.txt'),
+                        params=p,
+                    )
+                )
+                .result['CAPITAL COSTS (M$)']['Stimulation costs']['value']
+            )
+
+        result_default_indirect_cost: GeophiresXResult = _get_result()
+
+        higher_indirect = 12
+        result_higher_indirect_cost: GeophiresXResult = _get_result(higher_indirect)
+
+        self.assertAlmostEqual(
+            result_default_indirect_cost / 1.05,
+            result_higher_indirect_cost / (1 + (higher_indirect / 100)),
+            places=1,
+        )