Honor 'Capital Cost for Power Plant for Electricity Generation' parameter for all electricity end-use, not just AGS. Affects Fervo_Project_Cape-4 and example_SAM-single-owner-PPA, which have this parameter in their inputs (but was not having an effect on results, which was not noticed until today)

Author: softwareengineerprogrammer

src/geophires_x/AGSEconomics.py

@@ -74,19 +74,7 @@ def __init__(self, model: Model):
             Required=False,
             ErrMessage="assume default Capital cost for surface plant for direct-use system (100 $/kWth)"
         )
-        self._default_Power_plant_cost_USD_per_kWe = 3000
-        self.Power_plant_cost_per_kWe = self.ParameterDict[self.Power_plant_cost_per_kWe.Name] = floatParameter(
-            "Capital Cost for Power Plant for Electricity Generation",
-            DefaultValue=self._default_Power_plant_cost_USD_per_kWe,
-            Min=0.0,
-            Max=10000.0,
-            UnitType=Units.ENERGYCOST,
-            PreferredUnits=EnergyCostUnit.DOLLARSPERKW,
-            CurrentUnits=EnergyCostUnit.DOLLARSPERKW,
-            Required=True,
-            ErrMessage=f'assume default Power plant capital cost per kWe '
-                       f'({self._default_Power_plant_cost_USD_per_kWe} USD/kWe)'
-        )
+        self.Power_plant_cost_per_kWe.Required = True
 
         # results are stored here and in the parent ProducedTemperature array
 

src/geophires_x/Economics.py

@@ -700,6 +700,7 @@ def __init__(self, model: Model):
             Valid=True,
             ToolTipText="Multiplier for built-in wellfield O&M cost correlation"
         )
+
         self.ccplantfixed = self.ParameterDict[self.ccplantfixed.Name] = floatParameter(
             "Surface Plant Capital Cost",
             DefaultValue=-1.0,
@@ -724,6 +725,19 @@ def __init__(self, model: Model):
             Valid=True,
             ToolTipText="Multiplier for built-in surface plant capital cost correlation"
         )
+        self._default_Power_plant_cost_USD_per_kWe = 3000
+        self.Power_plant_cost_per_kWe = self.ParameterDict[self.Power_plant_cost_per_kWe.Name] = floatParameter(
+            "Capital Cost for Power Plant for Electricity Generation",
+            DefaultValue=self._default_Power_plant_cost_USD_per_kWe,
+            Min=0.0,
+            Max=10000.0,
+            UnitType=Units.ENERGYCOST,
+            PreferredUnits=EnergyCostUnit.DOLLARSPERKW,
+            CurrentUnits=EnergyCostUnit.DOLLARSPERKW,
+            ErrMessage=f'assume default Power plant capital cost per kWe '
+                       f'({self._default_Power_plant_cost_USD_per_kWe} USD/kWe)'
+        )
+
         self.ccgathfixed = self.ParameterDict[self.ccgathfixed.Name] = floatParameter(
             "Field Gathering System Capital Cost",
             DefaultValue=-1.0,
@@ -2826,13 +2840,19 @@ def calculate_plant_costs(self, model:Model) -> None:
                 self.CAPEX_cost_electricity_plant = self.Cplant.value * self.CAPEX_heat_electricity_plant_ratio.value
                 self.CAPEX_cost_heat_plant = self.Cplant.value * (1.0 - self.CAPEX_heat_electricity_plant_ratio.value)
             else:
-                # 1.02 to convert cost from 2012 to 2016
+                if self.Power_plant_cost_per_kWe.Provided:
+                    nameplate_capacity_kW = np.max(model.surfaceplant.ElectricityProduced.quantity().to('kW'))
+                    direct_plant_cost_MUSD = (nameplate_capacity_kW.magnitude *
+                                              model.economics.Power_plant_cost_per_kWe
+                                              .quantity().to('MUSD / kW').magnitude)
+                else:
+                    # 1.02 to convert cost from 2012 to 2016
+                    # factor 1.10 to convert from 2016 to 2022
+                    direct_plant_cost_MUSD = self.ccplantadjfactor.value * self.Cplantcorrelation * 1.02 * 1.10
 
                 # factor 1.15 for 15% contingency and 1.12 for 12% indirect costs.
                 # TODO https://github.com/NREL/GEOPHIRES-X/issues/383?title=Parameterize+indirect+cost+factor
-
-                # factor 1.10 to convert from 2016 to 2022
-                self.Cplant.value = 1.12 * 1.15 * self.ccplantadjfactor.value * self.Cplantcorrelation * 1.02 * 1.10
+                self.Cplant.value = 1.12 * 1.15 * direct_plant_cost_MUSD
                 self.CAPEX_cost_electricity_plant = self.Cplant.value
 
         # add direct-use plant cost of co-gen system to Cplant (only of no total Cplant was provided)

src/geophires_x_schema_generator/geophires-request.json

@@ -37,8 +37,7 @@
     "Dead-state Pressure",
     "Economic Model",
     "Time steps per year",
-    "Operation & Maintenance Cost of Surface Plant",
-    "Capital Cost for Power Plant for Electricity Generation"
+    "Operation & Maintenance Cost of Surface Plant"
   ],
   "properties": {
     "Reservoir Model": {
@@ -1486,6 +1485,15 @@
       "minimum": 0,
       "maximum": 10
     },
+    "Capital Cost for Power Plant for Electricity Generation": {
+      "description": "",
+      "type": "number",
+      "units": "USD/kW",
+      "category": "Economics",
+      "default": 3000,
+      "minimum": 0.0,
+      "maximum": 10000.0
+    },
     "Field Gathering System Capital Cost": {
       "description": "Total field gathering system capital cost",
       "type": "number",
@@ -2230,15 +2238,6 @@
       "minimum": 0.0,
       "maximum": 10000.0
     },
-    "Capital Cost for Power Plant for Electricity Generation": {
-      "description": "",
-      "type": "number",
-      "units": "USD/kW",
-      "category": "Economics",
-      "default": 3000.0,
-      "minimum": 0.0,
-      "maximum": 10000.0
-    },
     "AddOn Nickname": {
       "description": "If using multiple add-ons: either (1) specify this value as an array or (2) use multiple parameters suffixed with a number e.g. 'AddOn Nickname 1', 'AddOn Nickname 2', etc.",
       "type": "array",