Merge pull request NREL#414 from softwareengineerprogrammer/fervo_project_cape-4_stimulation-costs-update

softwareengineerprogrammer · web-flow · commit 0370d818a509 · 2025-08-13T07:08:41.000-07:00
`Fervo_Project_Cape-4` Stimulation Costs Update
diff --git a/.bumpversion.cfg b/.bumpversion.cfg
@@ -1,5 +1,5 @@
 [bumpversion]
-current_version = 3.9.48
+current_version = 3.9.49
 commit = True
 tag = True
 
diff --git a/.cookiecutterrc b/.cookiecutterrc
@@ -54,7 +54,7 @@ default_context:
     sphinx_doctest: "no"
     sphinx_theme: "sphinx-py3doc-enhanced-theme"
     test_matrix_separate_coverage: "no"
-    version: 3.9.48
+    version: 3.9.49
     version_manager: "bump2version"
     website: "https://github.com/NREL"
     year_from: "2023"
diff --git a/.gitignore b/.gitignore
@@ -25,6 +25,7 @@ Geothermal_district_heating_system_with_peaking_boilers.png
 
 # TODO may want to add this to source
 regenerate-schemas.sh
+requirements_2025-08-11.txt
 
 # C extensions
 *.so
@@ -112,3 +113,19 @@ _site/
 /src/geophires_x/CLG Simulator/clg_tea_v3_coaxial_water_elec.py
 /src/geophires_x/CLG Simulator/clg_tea_v3_SCO2_elec.py
 /src/geophires_x/CLG Simulator/clgs_v2.py
+
+# TOUGH
+CO2TAB
+Doublet.dat
+Doublet.out
+FOFT_*.csv
+GENER
+GOFT_*.csv
+INCON
+INFILE
+MESH*
+OUTPUT_ELEME.csv
+SAVE
+TABLE
+rcc3.out
+tough_debug_output/
diff --git a/CHANGELOG.rst b/CHANGELOG.rst
@@ -8,6 +8,8 @@ GEOPHIRES v3 (2023-2025)
 3.9
 ^^^
 
+3.9.49: `Fervo_Project_Cape-4 Stimulation Costs Update <https://github.com/NREL/GEOPHIRES-X/pull/414>`__ | `release <https://github.com/NREL/GEOPHIRES-X/releases/tag/v3.9.49>`__
+
 3.9.48: `Well diameter outputs, surface plant cost, non-SAM-EM Inflation costs during construction <https://github.com/NREL/GEOPHIRES-X/pull/412>`__ | `release <https://github.com/NREL/GEOPHIRES-X/releases/tag/v3.9.48>`__
 
 3.9.47: `Add-Ons support for SAM Economic Models <https://nrel.github.io/GEOPHIRES-X/SAM-Economic-Models.html#add-ons>`__ | `release <https://github.com/NREL/GEOPHIRES-X/releases/tag/v3.9.47>`__
diff --git a/README.rst b/README.rst
@@ -58,9 +58,9 @@ Free software: `MIT license <LICENSE>`__
     :alt: Supported implementations
     :target: https://pypi.org/project/geophires-x
 
-.. |commits-since| image:: https://img.shields.io/github/commits-since/softwareengineerprogrammer/GEOPHIRES-X/v3.9.48.svg
+.. |commits-since| image:: https://img.shields.io/github/commits-since/softwareengineerprogrammer/GEOPHIRES-X/v3.9.49.svg
     :alt: Commits since latest release
-    :target: https://github.com/softwareengineerprogrammer/GEOPHIRES-X/compare/v3.9.48...main
+    :target: https://github.com/softwareengineerprogrammer/GEOPHIRES-X/compare/v3.9.49...main
 
 .. |docs| image:: https://readthedocs.org/projects/GEOPHIRES-X/badge/?style=flat
     :target: https://nrel.github.io/GEOPHIRES-X
diff --git a/docs/Fervo_Project_Cape-4.md b/docs/Fervo_Project_Cape-4.md
diff --git a/docs/conf.py b/docs/conf.py
@@ -18,7 +18,7 @@
 year = '2025'
 author = 'NREL'
 copyright = f'{year}, {author}'
-version = release = '3.9.48'
+version = release = '3.9.49'
 
 pygments_style = 'trac'
 templates_path = ['./templates']
diff --git a/setup.py b/setup.py
@@ -13,7 +13,7 @@ def read(*names, **kwargs):
 
 setup(
     name='geophires-x',
-    version='3.9.48',
+    version='3.9.49',
     license='MIT',
     description='GEOPHIRES is a free and open-source geothermal techno-economic simulator.',
     long_description='{}\n{}'.format(
@@ -77,7 +77,9 @@ def read(*names, **kwargs):
         'h5py',
         'scipy',
         'iapws',
-        'coolprop; python_version > "3.8"',
+        # coolprop 7.0.0 fails on macOS + Python > 3.8, see
+        # https://github.com/softwareengineerprogrammer/GEOPHIRES/actions/runs/16888333182/job/47841807973
+        'coolprop == 6.8.0; python_version > "3.8"',
         # coolprop > 6.7.0 fails on Python 3.8, see
         # https://github.com/softwareengineerprogrammer/GEOPHIRES/actions/runs/14599234458/job/40952924048?pr=69#step:5:302
         'coolprop==6.7.0; python_version == "3.8"',
diff --git a/src/geophires_x/__init__.py b/src/geophires_x/__init__.py
@@ -1 +1 @@
-__version__ = '3.9.48'
+__version__ = '3.9.49'
diff --git a/tests/examples/Fervo_Project_Cape-4.out b/tests/examples/Fervo_Project_Cape-4.out
diff --git a/tests/examples/Fervo_Project_Cape-4.txt b/tests/examples/Fervo_Project_Cape-4.txt
@@ -17,10 +17,8 @@ Fraction of Investment in Bonds, .6, -- Based on fraction of CAPEX with $1 billi
 Inflated Bond Interest Rate, .056
 
 Inflation Rate, .023, -- US inflation as of April 2025
-Inflation Rate During Construction, 0.15, -- Models a higher annual inflation rate (4.769%) over the 3-year construction period as a hedge against short-term inflation volatility.
 Construction Years, 1, -- Calibrated to a 2-6 year construction time for a 1 GW plant (Yusifov & Enriquez, 2025)
 
-
 Combined Income Tax Rate, .28
 Investment Tax Credit Rate, 0.3
 Property Tax Rate, 0
@@ -31,7 +29,10 @@ Exploration Capital Cost, 30, -- Estimate significantly higher exploration costs
 Well Drilling Cost Correlation, 3, -- VERTICAL_LARGE (2025 NREL Geothermal Drilling Cost Curve Update)
 Well Drilling and Completion Capital Cost Adjustment Factor, 0.8, -- Adjust correlation-calculated value of $4.72M/well to $3.96M/well per Tim Latimer on 2025-02-12 Volts podcast: less than $4M/well
 
-Reservoir Stimulation Capital Cost Adjustment Factor, 2.66, -- Estimated cost of ~$2M per well. Typical range for Nth-of-kind projects may be $0.5–2M.
+Reservoir Stimulation Capital Cost per Injection Well, 4, -- Baseline stimulation cost of $4.0M, calibrated from high-intensity U.S. shale well analogue (~$39k/frac stage for 102 stages). This is a pre-contingency value and excludes EGS-specific cost premiums such as ceramic proppant and HPHT hardware.
+Reservoir Stimulation Capital Cost per Production Well, 4, -- Baseline stimulation cost of $4.0M, calibrated from high-intensity U.S. shale well analogue (~$39k/frac stage for 102 stages). This is a pre-contingency value and excludes EGS-specific cost premiums such as ceramic proppant and HPHT hardware.
+Reservoir Stimulation Indirect Capital Cost Percentage, 0, -- Baseline stimulation cost includes indirect costs
+
 Field Gathering System Capital Cost Adjustment Factor, 0.54, -- Gathering costs represent 2% of facilities CAPEX per https://www.linkedin.com/pulse/fervo-energy-technology-day-2024-entering-geothermal-decade-matson-n4stc/
 
 # *** SURFACE & SUBSURFACE TECHNICAL PARAMETERS ***
diff --git a/tests/geophires_x_tests/test_fervo_project_cape_4.py b/tests/geophires_x_tests/test_fervo_project_cape_4.py
@@ -72,8 +72,17 @@ def test_case_study_documentation(self):
         inputs_in_markdown = self.parse_markdown_inputs_structured(documentation_file_content)
         results_in_markdown = self.parse_markdown_results_structured(documentation_file_content)
 
-        self.assertEqual(3.96, results_in_markdown['Well Drilling and Completion Cost']['value'])
-        self.assertEqual('MUSD/well', results_in_markdown['Well Drilling and Completion Cost']['unit'])
+        self.assertEqual(3.96, results_in_markdown['Well Drilling and Completion Costs']['value'])
+        self.assertEqual('MUSD/well', results_in_markdown['Well Drilling and Completion Costs']['unit'])
+
+        expected_stim_cost_MUSD_per_well = 4.6
+        self.assertEqual(expected_stim_cost_MUSD_per_well, results_in_markdown['Stimulation Costs']['value'])
+        self.assertEqual('MUSD/well', results_in_markdown['Stimulation Costs']['unit'])
+
+        self.assertEqual(
+            expected_stim_cost_MUSD_per_well, inputs_in_markdown['Reservoir Stimulation Capital Cost per Well']['value']
+        )
+        self.assertEqual('MUSD', inputs_in_markdown['Reservoir Stimulation Capital Cost per Well']['unit'])
 
         class _Q(HasQuantity):
             def __init__(self, vu: dict[str, Any]):
@@ -96,10 +105,16 @@ def __init__(self, vu: dict[str, Any]):
         field_mapping = {
             'LCOE': 'Electricity breakeven price',
             'Project capital costs: Total CAPEX': 'Total CAPEX',
-            'Well Drilling and Completion Cost': 'Drilling and completion costs per well',
+            'Well Drilling and Completion Costs': 'Drilling and completion costs per well',
+            'Well Drilling and Completion Costs total': 'Drilling and completion costs',
+            'Stimulation Costs total': 'Stimulation costs',
         }
 
-        ignore_keys = ['Project capital costs: $/kW', 'Total fracture surface area per production well']
+        ignore_keys = [
+            'Project capital costs: $/kW',  # See https://github.com/NREL/GEOPHIRES-X/issues/391
+            'Total fracture surface area per production well',
+            'Stimulation Costs',  # remapped to 'Stimulation Costs total'
+        ]
 
         example_result = GeophiresXResult(self._get_test_file_path('../examples/Fervo_Project_Cape-4.out'))
         example_result_values = {}
@@ -116,10 +131,10 @@ def __init__(self, vu: dict[str, Any]):
             if ignore_key in results_in_markdown:
                 del results_in_markdown[ignore_key]
 
-        results_in_markdown['Well Drilling and Completion Cost']['unit'] = results_in_markdown[
-            'Well Drilling and Completion Cost'
+        results_in_markdown['Well Drilling and Completion Costs']['unit'] = results_in_markdown[
+            'Well Drilling and Completion Costs'
         ]['unit'].replace('/well', '')
-        self.assertDictAlmostEqual(example_result_values, results_in_markdown, places=3)
+        self.assertDictAlmostEqual(example_result_values, results_in_markdown, percent=0.185)
 
         result_capex_USD_per_kW = (
             _Q(example_result._get_result_field('Total CAPEX')).quantity().to('USD').magnitude
@@ -143,6 +158,12 @@ def __init__(self, vu: dict[str, Any]):
             inputs_in_markdown['Reservoir Volume']['value'],
         )
 
+        expected_stim_cost_total_MUSD = expected_stim_cost_MUSD_per_well * num_doublets * 2
+        self.assertEqual(
+            expected_stim_cost_total_MUSD,
+            example_result.result['CAPITAL COSTS (M$)']['Stimulation costs']['value'],
+        )
+
     def parse_markdown_results_structured(self, markdown_text: str) -> dict:
         """
         Parses result values from markdown into a structured dictionary with values and units.
@@ -192,10 +213,24 @@ def parse_markdown_results_structured(self, markdown_text: str) -> dict:
                 'Project capital costs: Total CAPEX',
                 'Project capital costs: $/kW',
                 'WACC',
-                'Well Drilling and Completion Cost',
+                'Well Drilling and Completion Costs',
+                'Stimulation Costs',
             ]:
                 structured_results[key_] = self._parse_value_unit(value_)
 
+        for result_with_total_key in ['Well Drilling and Completion Costs', 'Stimulation Costs']:
+            entry = structured_results[result_with_total_key]
+
+            unit_str = entry['unit']
+            entry['unit'] = unit_str.split(';')[0]
+
+            # This probably could/should be done with some adaptation of self._parse_value_unit, but one-off parsing
+            # here is fine for now...
+            structured_results[f'{result_with_total_key} total'] = {
+                'value': float(unit_str.split('; ')[1].replace(' total', '').replace('$', '').replace('M', '')),
+                'unit': entry['unit'].split('/')[0],
+            }
+
         return structured_results
 
     def parse_markdown_inputs_structured(self, markdown_text: str) -> dict:

Original file line number	Diff line number	Diff line change
`@@ -1 +1 @@`
`1`		`-__version__ = '3.9.48'`
	`1`	`+__version__ = '3.9.49'`