after-tax LCOE

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -471,7 +471,8 @@ def CalculateLCOELCOHLCOC(econ, model: Model) -> tuple:
                 model.surfaceplant.annual_heating_demand.value * discount_vector) * 1E2  # cents/kWh
             LCOH = LCOH * 2.931  # $/Million Btu
     elif econ.econmodel.value == EconomicModel.SAM_SINGLE_OWNER_PPA:
-        LCOE = calculate_sam_economics(model)['LCOE']['value']
+        # FIXME TODO designate nominal (as opposed to real) in client result
+        LCOE = calculate_sam_economics(model)['LCOE (nominal)']['value']
     else:
         # must be BICYCLE
         # average return on investment (tax and inflation adjusted)

src/geophires_x/EconomicsSam.py

@@ -56,7 +56,7 @@ def calculate_sam_economics(model: Model) -> dict[str, dict[str, Any]]:
     cash_flow = _calculate_cash_flow(model, single_owner)
 
     data = [
-        ('LCOE', single_owner.Outputs.lcoe_real, 'cents/kWh'),
+        ('LCOE (nominal)', single_owner.Outputs.lcoe_nom, 'cents/kWh'),
         ('IRR', single_owner.Outputs.project_return_aftertax_irr, '%'),
         ('NPV', single_owner.Outputs.project_return_aftertax_npv * 1e-6, 'MUSD'),
         ('CAPEX', single_owner.Outputs.adjusted_installed_cost * 1e-6, 'MUSD'),
@@ -261,6 +261,24 @@ def single_value_row(row_name: str, single_value: float) -> list[Any]:
     data_row('After-tax cumulative IRR (%)', _soo.cf_project_return_aftertax_irr)
     data_row('After-tax cumulative NPV ($)', _soo.cf_project_return_aftertax_npv)
 
+    blank_row()
+
+    category_row('AFTER-TAX LCOE AND PPA PRICE')
+    data_row('Annual costs ($)', _soo.cf_annual_costs)
+    data_row('PPA revenue ($)', _soo.cf_energy_value)  # TODO config-ify repeated
+    data_row('Electricity to grid (kWh)', _soo.cf_energy_sales)  # TODO config-ify repeated
+
+    blank_row()
+    single_value_row('Present value of annual costs ($)', _soo.npv_annual_costs)
+    single_value_row('Present value of annual energy nominal ($)', _soo.npv_energy_nom)  # TODO config-ify repeated
+    single_value_row('LCOE Levelized cost of energy nominal (cents/kWh)', _soo.lcoe_nom)
+
+    blank_row()
+
+    single_value_row('Present value of PPA revenue ($)', _soo.npv_ppa_revenue)
+    single_value_row('Present value of annual energy nominal ($)', _soo.npv_energy_nom)  # TODO config-ify repeated
+    single_value_row('LPPA Levelized PPA price nominal (cents/kWh)', _soo.lppa_nom)
+
     return profile
 
 

tests/geophires_x_tests/generic-egs-case-2.txt

@@ -43,7 +43,7 @@ Surface Plant Capital Cost Adjustment Factor, 0.71, -- $1900/kW (peak production
 Exploration Capital Cost, 30
 Well Drilling and Completion Capital Cost, 3.84, -- Assume cost continues to decrease from $4.8M/well (https://houston.innovationmap.com/fervo-energy-drilling-utah-project-2667300142.html)
 Reservoir Stimulation Capital Cost, 234, -- 78 wells @ $3M/well
-Economic Model, 3
+Economic Model, 5, -- SAM Single Owner PPA
 Starting Electricity Sale Price, 0.15
 Ending Electricity Sale Price, 1.00
 Electricity Escalation Rate Per Year, 0.004053223

tests/geophires_x_tests/test_economics_sam.py

@@ -39,7 +39,7 @@ def _npv(r: GeophiresXResult) -> float:
 
         base_result = self._get_result({})
         base_lcoe = _lcoe(base_result)
-        self.assertGreater(base_lcoe, 6)
+        self.assertGreater(base_lcoe, 7)
 
         npvs = [_npv(self._get_result({'Starting Electricity Sale Price': x / 100.0})) for x in range(1, 20, 4)]
         for i in range(len(npvs) - 1):
@@ -71,13 +71,20 @@ def test_cash_flow(self):
         def get_row(name: str) -> list[float]:
             return next(r for r in cash_flow if r[0] == name)[1:]
 
+        def get_single_value(name: str) -> list[float]:
+            return get_row(name)[0]
+
         self.assertListEqual(get_row('PPA revenue ($)'), get_row('Total revenue ($)'))
 
-        tic = get_row('Total installed cost ($)')[0]
+        tic = get_single_value('Total installed cost ($)')
         self.assertLess(tic, 0)
-        self.assertAlmostEqual(get_row('Cash flow from investing activities ($)')[0], tic, places=2)
+        self.assertAlmostEqual(get_single_value('Cash flow from investing activities ($)'), tic, places=2)
+
+        self.assertAlmostEqual(get_single_value('Cash flow from financing activities ($)'), -1.0 * tic, places=2)
 
-        self.assertAlmostEqual(get_row('Cash flow from financing activities ($)')[0], -1.0 * tic, places=2)
+        self.assertAlmostEqual(
+            m.economics.LCOE.value, get_single_value('LCOE Levelized cost of energy nominal (cents/kWh)'), places=2
+        )
 
     def test_only_electricity_end_use_supported(self):
         with self.assertRaises(RuntimeError):