Unit test to ensure difference between correlation values and GEOPHIRES-calculated drilling cost per well is due to 5% indirect cost factor

Author: softwareengineerprogrammer

tests/geophires_x_tests/test_options_list.py

@@ -1,3 +1,5 @@
+from typing import ClassVar
+
 from base_test_case import BaseTestCase
 from geophires_x.OptionList import EndUseOptions
 from geophires_x.OptionList import PlantType
@@ -33,23 +35,33 @@ def test_equality(self):
 
 
 class WellDrillingCostCorrelationTestCase(BaseTestCase):
+
     def test_equality(self):
         self.assertFalse(WellDrillingCostCorrelation.VERTICAL_SMALL == WellDrillingCostCorrelation.DEVIATED_SMALL)
         self.assertTrue(WellDrillingCostCorrelation.VERTICAL_SMALL == WellDrillingCostCorrelation.VERTICAL_SMALL)
 
-    def test_baseline_curve_costs(self):
-        # Sanity-check calibration with NREL 2025 Cost Curve Update
-        # https://pangea.stanford.edu/ERE/db/GeoConf/papers/SGW/2025/Akindipe.pdf?t=1740084555
-
-        self.assertAlmostEqual(5.1, WellDrillingCostCorrelation.VERTICAL_SMALL.calculate_cost_MUSD(3500), delta=0.1)
-        self.assertAlmostEqual(13.9, WellDrillingCostCorrelation.VERTICAL_SMALL.calculate_cost_MUSD(6500), delta=0.1)
-        self.assertAlmostEqual(15.9, WellDrillingCostCorrelation.VERTICAL_SMALL.calculate_cost_MUSD(7000), delta=0.1)
-
-        self.assertAlmostEqual(17.2, WellDrillingCostCorrelation.VERTICAL_LARGE.calculate_cost_MUSD(6500), delta=0.1)
-
-        self.assertAlmostEqual(14.9, WellDrillingCostCorrelation.DEVIATED_SMALL.calculate_cost_MUSD(6500), delta=0.1)
-
-        self.assertAlmostEqual(18.3, WellDrillingCostCorrelation.DEVIATED_LARGE.calculate_cost_MUSD(6500), delta=0.1)
+    COST_CORRELATION_TEST_CASES: ClassVar[list[tuple[WellDrillingCostCorrelation, int, float]]] = [
+        (WellDrillingCostCorrelation.VERTICAL_SMALL, 3500, 5.1),
+        (WellDrillingCostCorrelation.VERTICAL_SMALL, 6500, 13.9),
+        (WellDrillingCostCorrelation.VERTICAL_SMALL, 7000, 15.9),
+        (WellDrillingCostCorrelation.VERTICAL_LARGE, 6500, 17.2),
+        (WellDrillingCostCorrelation.DEVIATED_SMALL, 6500, 14.9),
+        (WellDrillingCostCorrelation.DEVIATED_LARGE, 6500, 18.3),
+    ]
+
+    def test_drilling_cost_curve_correlations(self):
+        """
+        Check calibration with NREL 2025 Cost Curve Update.
+        Values in COST_CORRELATION_TEST_CASES derived from graphs on p. 8 of
+        https://pangea.stanford.edu/ERE/db/GeoConf/papers/SGW/2025/Akindipe.pdf?t=1740084555
+        """
+
+        for test_case in WellDrillingCostCorrelationTestCase.COST_CORRELATION_TEST_CASES:
+            correlation: WellDrillingCostCorrelation = test_case[0]
+            depth_m = test_case[1]
+            expected_cost_musd = test_case[2]
+            with self.subTest(msg=str(f'{correlation.name}, {depth_m}m')):
+                self.assertAlmostEqual(expected_cost_musd, correlation.calculate_cost_MUSD(depth_m), delta=0.1)
 
 
 class PlantTypeTestCase(BaseTestCase):

tests/test_geophires_x.py

@@ -5,11 +5,13 @@
 from typing import Optional
 
 from geophires_x.OptionList import PlantType
+from geophires_x.OptionList import WellDrillingCostCorrelation
 from geophires_x_client import GeophiresXClient
 from geophires_x_client import GeophiresXResult
 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_tests.test_options_list import WellDrillingCostCorrelationTestCase
 from tests.base_test_case import BaseTestCase
 
 
@@ -784,3 +786,36 @@ def test_multiple_construction_years(self):
         pre_rev_idx = min(net_rev_idx, net_cashflow_idx)
         self.assertListEqual([0] * pre_rev_idx, rcp[1][:pre_rev_idx])
         self.assertListEqual([1] + [0] * (pre_rev_idx - 1), rcp[2][:pre_rev_idx])
+
+    def test_drilling_cost_curves(self):
+        """
+        Note this is similar to
+        geophires_x_tests.test_options_list.WellDrillingCostCorrelationTestCase.test_drilling_cost_curve_correlations;
+        this test ensures that the indirect cost factor is responsible for the discrepancy between GEOPHIRES-calculated
+        drilling cost per well and the raw value calculated by the curve.
+        """
+
+        indirect_cost_factor = 1.05  # See TODO re:parameterizing at src/geophires_x/Economics.py:652
+
+        for test_case in WellDrillingCostCorrelationTestCase.COST_CORRELATION_TEST_CASES:
+            correlation: WellDrillingCostCorrelation = test_case[0]
+            depth_m = test_case[1]
+            expected_cost_musd = test_case[2]
+            with self.subTest(msg=str(f'{correlation.name}, {depth_m}m')):
+                result = GeophiresXClient().get_geophires_result(
+                    GeophiresInputParameters(
+                        from_file_path=self._get_test_file_path('geophires_x_tests/generic-egs-case.txt'),
+                        params={
+                            'Well Drilling Cost Correlation': correlation.int_value,
+                            'Reservoir Depth': f'{depth_m} meter',
+                            'Number of Production Wells': 1,
+                            'Number of Injection Wells': 1,
+                            'Well Drilling and Completion Capital Cost Adjustment Factor': 1,
+                        },
+                    )
+                )
+
+                cost_per_well_val = result.result['CAPITAL COSTS (M$)']['Drilling and completion costs per well'][
+                    'value'
+                ]
+                self.assertAlmostEqual(indirect_cost_factor * expected_cost_musd, cost_per_well_val, delta=0.1)