WIP - display redrilling costs in O&M

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -1775,6 +1775,14 @@ def __init__(self, model: Model):
             UnitType=Units.CURRENCYFREQUENCY,
             PreferredUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR,
             CurrentUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR
+            # TODO TooltipText to document how this is calculated
+        )
+        self.redrilling_annual_cost = self.OutputParameterDict[self.redrilling_annual_cost.Name] = OutputParameter(
+            Name="Redrilling costs",
+            UnitType=Units.CURRENCYFREQUENCY,
+            PreferredUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR,
+            CurrentUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR,
+            # FIXME WIP TooltipText
         )
         self.Cplant = self.OutputParameterDict[self.Cplant.Name] = OutputParameter(
             Name="Surface Plant cost",
@@ -1809,7 +1817,7 @@ def __init__(self, model: Model):
             UnitType=Units.CURRENCYFREQUENCY,
             PreferredUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR,
             CurrentUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR,
-            ToolTipText='Assumes $3.5/1,000 gallons of water'
+            ToolTipText='Assumes $3.5/1,000 gallons of water'  # TODO parameterize
         )
         self.CCap = self.OutputParameterDict[self.CCap.Name] = OutputParameter(
             Name="Total Capital Cost",
@@ -2521,8 +2529,10 @@ def Calculate(self, model: Model) -> None:
 
         if model.wellbores.redrill.value > 0:
             # account for well redrilling
-            self.Coam.value = self.Coam.value + \
-                              (self.Cwell.value + self.Cstim.value) * model.wellbores.redrill.value / model.surfaceplant.plant_lifetime.value
+            redrilling_costs: PlainQuantity = self.calculate_redrilling_costs(model)
+            self.redrilling_annual_cost.value = redrilling_costs.to(self.redrilling_annual_cost.CurrentUnits).magnitude
+            self.Coam.value += redrilling_costs.to(self.Coam.CurrentUnits).magnitude
+
 
         # Add in the AnnualLicenseEtc and TaxRelief
         self.Coam.value = self.Coam.value + self.AnnualLicenseEtc.value - self.TaxRelief.value
@@ -2751,6 +2761,11 @@ def calculate_stimulation_costs(self, model: Model) -> PlainQuantity:
 
         return quantity(stimulation_costs, self.Cstim.CurrentUnits)
 
+    def calculate_redrilling_costs(self, model) -> PlainQuantity:
+        return ((self.Cwell.quantity() + self.Cstim.quantity())
+                * model.wellbores.redrill.quantity()
+                / model.surfaceplant.plant_lifetime.quantity())
+
     def calculate_field_gathering_costs(self, model: Model) -> None:
         if self.ccgathfixed.Valid:
             self.Cgath.value = self.ccgathfixed.value
@@ -3164,5 +3179,7 @@ def _calculate_derived_outputs(self, model: Model) -> None:
                 (model.wellbores.nprod.value + model.wellbores.ninj.value)
             )
 
+
+
     def __str__(self):
         return "Economics"

src/geophires_x/Outputs.py

@@ -475,10 +475,6 @@ def PrintOutputs(self, model: Model):
                         f.write(f'         District Heating System Cost:                  {model.economics.dhdistrictcost.value:10.2f} {model.economics.dhdistrictcost.CurrentUnits.value}\n')
                     f.write(f'         Total surface equipment costs:                 {(model.economics.Cplant.value+model.economics.Cgath.value):10.2f} ' + model.economics.Cplant.CurrentUnits.value + NL)
                     f.write(f'         {model.economics.Cexpl.display_name}:                             {model.economics.Cexpl.value:10.2f} {model.economics.Cexpl.CurrentUnits.value}\n')
-                if model.economics.totalcapcost.Valid and model.wellbores.redrill.value > 0:
-                    f.write(f'         Drilling and completion costs (for redrilling):{model.economics.Cwell.value:10.2f} ' + model.economics.Cwell.CurrentUnits.value + NL)
-                    f.write(f'      Drilling and completion costs per redrilled well: {(model.economics.Cwell.value/(model.wellbores.nprod.value+model.wellbores.ninj.value)):10.2f} ' + model.economics.Cwell.CurrentUnits.value + NL)
-                    f.write(f'         Stimulation costs (for redrilling):            {econ.Cstim.value:10.2f} {econ.Cstim.CurrentUnits.value}\n')
                 if model.economics.RITCValue.value:
                     f.write(f'         {model.economics.RITCValue.display_name}:                         {-1*model.economics.RITCValue.value:10.2f} {model.economics.RITCValue.CurrentUnits.value}\n')
 
@@ -506,6 +502,10 @@ def PrintOutputs(self, model: Model):
                         f.write(f'         Annual District Heating O&M Cost:              {model.economics.dhdistrictoandmcost.value:10.2f} {model.economics.dhdistrictoandmcost.CurrentUnits.value}\n')
                         f.write(f'         Average Annual Peaking Fuel Cost:              {model.economics.averageannualngcost.value:10.2f} {model.economics.averageannualngcost.CurrentUnits.value}\n')
 
+                    if model.wellbores.redrill.value > 0:
+                        redrill_label = Outputs._field_label(econ.redrilling_annual_cost.display_name, 47)
+                        f.write(f'         {redrill_label}{econ.redrilling_annual_cost.value:10.2f} {econ.redrilling_annual_cost.CurrentUnits.value}\n')
+
                     f.write(f'      {econ.Coam.display_name}:            {(econ.Coam.value + econ.averageannualpumpingcosts.value + econ.averageannualheatpumpelectricitycost.value):10.2f} {econ.Coam.CurrentUnits.value}\n')
                 else:
                     f.write(f'      {econ.Coam.display_name}:            {econ.Coam.value:10.2f} {econ.Coam.CurrentUnits.value}\n')

src/geophires_x_client/geophires_x_result.py

@@ -290,6 +290,7 @@ class GeophiresXResult:
                 # SUTRA
                 'Average annual auxiliary fuel cost',
                 'Average annual pumping cost',
+                'Redrilling costs',
                 'Total average annual O&M costs',
             ],
             'SURFACE EQUIPMENT SIMULATION RESULTS': [

tests/examples/Fervo_Project_Cape-4.out

@@ -4,10 +4,10 @@
 
 Simulation Metadata
 ----------------------
- GEOPHIRES Version: 3.9.28
- Simulation Date: 2025-07-02
- Simulation Time:  12:19
- Calculation Time:      1.773 sec
+ GEOPHIRES Version: 3.9.36
+ Simulation Date: 2025-07-25
+ Simulation Time:  10:59
+ Calculation Time:      1.734 sec
 
                            ***SUMMARY OF RESULTS***
 
@@ -113,6 +113,7 @@ Simulation Metadata
          Wellfield maintenance costs:                         6.20 MUSD/yr
          Power plant maintenance costs:                      25.43 MUSD/yr
          Water costs:                                        24.86 MUSD/yr
+         Redrilling costs:                                   70.46 MUSD/yr
       Total operating and maintenance costs:                126.95 MUSD/yr