Add Reservoir Stimulation Capital Cost per Production Well parameter (default value = $0)

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -585,18 +585,38 @@ def __init__(self, model: Model):
             Valid=False,
             ToolTipText="Total reservoir stimulation capital cost"
         )
+
+        max_stimulation_cost_per_well_MUSD = 100
         self.stimulation_cost_per_injection_well = \
           self.ParameterDict[self.stimulation_cost_per_injection_well.Name] = floatParameter(
             'Reservoir Stimulation Capital Cost per Injection Well',
             DefaultValue=1.25,
             Min=0,
-            Max=100,
+            Max=max_stimulation_cost_per_well_MUSD,
             UnitType=Units.CURRENCY,
             PreferredUnits=CurrencyUnit.MDOLLARS,
             CurrentUnits=CurrencyUnit.MDOLLARS,
             Provided=False,
             ToolTipText='Reservoir stimulation capital cost per injection well'
         )
+
+        stimulation_cost_per_production_well_default_value_MUSD = 0
+        stimulation_cost_per_production_well_default_value_note = \
+            '. By default, only the injection wells are assumed to be stimulated unless this parameter is provided.' \
+                if stimulation_cost_per_production_well_default_value_MUSD == 0 else ''
+        self.stimulation_cost_per_production_well = \
+          self.ParameterDict[self.stimulation_cost_per_production_well.Name] = floatParameter(
+            'Reservoir Stimulation Capital Cost per Production Well',
+            DefaultValue=stimulation_cost_per_production_well_default_value_MUSD,
+            Min=0,
+            Max=max_stimulation_cost_per_well_MUSD,
+            UnitType=Units.CURRENCY,
+            PreferredUnits=CurrencyUnit.MDOLLARS,
+            CurrentUnits=CurrencyUnit.MDOLLARS,
+            ToolTipText=f'Reservoir stimulation capital cost per production well'
+                        f'{stimulation_cost_per_production_well_default_value_note}'
+        )
+
         self.ccstimadjfactor = self.ParameterDict[self.ccstimadjfactor.Name] = floatParameter(
             "Reservoir Stimulation Capital Cost Adjustment Factor",
             DefaultValue=1.0,
@@ -607,7 +627,7 @@ def __init__(self, model: Model):
             CurrentUnits=PercentUnit.TENTH,
             Provided=False,
             Valid=True,
-            ToolTipText="Multiplier for built-in reservoir stimulation capital cost correlation"
+            ToolTipText="Multiplier for reservoir stimulation capital cost correlation"
         )
         self.ccexplfixed = self.ParameterDict[self.ccexplfixed.Name] = floatParameter(
             "Exploration Capital Cost",
@@ -1629,8 +1649,9 @@ def __init__(self, model: Model):
             CurrentUnits=CurrencyUnit.MDOLLARS,
             ToolTipText=f'Default correlation: ${self.stimulation_cost_per_injection_well.value}M '
                         f'per injection well {stimulation_contingency_and_indirect_costs_tooltip}. '
-                        f'Provide {self.stimulation_cost_per_injection_well.Name} to set the correlation '
-                        f'cost per injection well. '
+                        f'Provide {self.stimulation_cost_per_injection_well.Name} and '
+                        f'{self.stimulation_cost_per_production_well.Name} to set the correlation '
+                        f'costs per well. '
                         f'Provide {self.ccstimadjfactor.Name} to multiply the correlation-calculated cost. '
                         f'Provide {self.ccstimfixed.Name} to override the correlation and set your own '
                         f'total stimulation cost.'
@@ -2358,11 +2379,15 @@ def Calculate(self, model: Model) -> None:
         else:
             stim_cost_per_injection_well = self.stimulation_cost_per_injection_well.quantity().to(
                 self.Cstim.CurrentUnits).magnitude
+            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
+            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)
 

src/geophires_x_schema_generator/geophires-request.json

@@ -1404,8 +1404,17 @@
       "minimum": 0,
       "maximum": 100
     },
+    "Reservoir Stimulation Capital Cost per Production Well": {
+      "description": "Reservoir stimulation capital cost per production well. By default, only the injection wells are assumed to be stimulated unless this parameter is provided.",
+      "type": "number",
+      "units": "MUSD",
+      "category": "Economics",
+      "default": 0,
+      "minimum": 0,
+      "maximum": 100
+    },
     "Reservoir Stimulation Capital Cost Adjustment Factor": {
-      "description": "Multiplier for built-in reservoir stimulation capital cost correlation",
+      "description": "Multiplier for reservoir stimulation capital cost correlation",
       "type": "number",
       "units": "",
       "category": "Economics",

src/geophires_x_schema_generator/geophires-result.json

@@ -351,7 +351,7 @@
         "Drilling and completion costs per redrilled well": {},
         "Stimulation costs": {
           "type": "number",
-          "description": "Default correlation: $1.25M per injection well plus 15% contingency plus 5% indirect costs. Provide Reservoir Stimulation Capital Cost per Injection Well to set the correlation cost per injection well. Provide Reservoir Stimulation Capital Cost Adjustment Factor to multiply the correlation-calculated cost. Provide Reservoir Stimulation Capital Cost to override the correlation and set your own total stimulation cost.",
+          "description": "Default correlation: $1.25M per injection well plus 15% contingency plus 5% indirect costs. Provide Reservoir Stimulation Capital Cost per Injection Well and Reservoir Stimulation Capital Cost per Production Well to set the correlation costs per well. Provide Reservoir Stimulation Capital Cost Adjustment Factor to multiply the correlation-calculated cost. Provide Reservoir Stimulation Capital Cost to override the correlation and set your own total stimulation cost.",
           "units": "MUSD"
         },
         "Stimulation costs (for redrilling)": {},

tests/test_geophires_x.py

@@ -12,6 +12,7 @@
 from geophires_x_client import _get_logger
 from geophires_x_client.geophires_input_parameters import EndUseOption
 from geophires_x_client.geophires_input_parameters import GeophiresInputParameters
+from geophires_x_client.geophires_input_parameters import ImmutableGeophiresInputParameters
 from geophires_x_tests.test_options_list import WellDrillingCostCorrelationTestCase
 from tests.base_test_case import BaseTestCase
 
@@ -941,3 +942,36 @@ def test_sbt_coaxial_raises_error(self):
             )
             client.get_geophires_result(params)
         self.assertIn('SBT with coaxial configuration is not implemented', str(e.exception))
+
+    def test_production_well_stimulation_cost(self):
+        def _get_result(prod_well_stim_MUSD: Optional[int] = None) -> GeophiresXResult:
+            p = {}
+            if prod_well_stim_MUSD is not None:
+                p['Reservoir Stimulation Capital Cost per Production Well'] = prod_well_stim_MUSD
+
+            return GeophiresXClient().get_geophires_result(
+                ImmutableGeophiresInputParameters(
+                    from_file_path=self._get_test_file_path('geophires_x_tests/generic-egs-case.txt'),
+                    params=p,
+                )
+            )
+
+        result_no_prod_stim: GeophiresXResult = _get_result()
+
+        result_prod_stim: GeophiresXResult = _get_result(1.25)
+
+        # TODO https://github.com/NREL/GEOPHIRES-X/issues/383?title=Parameterize+indirect+cost+factor
+        indirect_and_contingency = 1.05 * 1.15
+
+        self.assertAlmostEqual(
+            (
+                2
+                * (
+                    result_no_prod_stim.result['CAPITAL COSTS (M$)']['Stimulation costs']['value']
+                    / (indirect_and_contingency)
+                )
+            )
+            * indirect_and_contingency,
+            result_prod_stim.result['CAPITAL COSTS (M$)']['Stimulation costs']['value'],
+            places=1,
+        )