Doing a Monte Carlo Simulation for Fervo. Also, extending the maximum possible range for Plant Outlet Pressure

Author: malcolm-dsider

src/geophires_x/SurfacePlant.py

@@ -342,7 +342,7 @@ def __init__(self, model: Model):
             "Plant Outlet Pressure",
             DefaultValue=100.0,
             Min=0.01,
-            Max=10000.0,
+            Max=15000.0,
             UnitType=Units.PRESSURE,
             PreferredUnits=PressureUnit.KPASCAL,
             CurrentUnits=PressureUnit.KPASCAL,
@@ -576,16 +576,16 @@ def read_parameters(self, model:Model) -> None:
                                 model.wellbores.impedancemodelallowed.value = False
                                 self.setinjectionpressurefixed = True
                     elif ParameterToModify.Name == 'Plant Outlet Pressure':
-                        if ParameterToModify.value < 0 or ParameterToModify.value > 10000:
+                        if ParameterToModify.value < self.plant_outlet_pressure.Min or ParameterToModify.value > self.plant_outlet_pressure.Max:
                                 if self.setinjectionpressurefixed:
                                     ParameterToModify.value = 100
-                                    msg = (f'Provided plant outlet pressure outside of range 0-10000. GEOPHIRES will '
+                                    msg = (f'Provided plant outlet pressure outside of range defined valid range. GEOPHIRES will '
                                            f'assume default plant outlet pressure ({ParameterToModify.value} kPa)')
                                     print(f'Warning: {msg}')
                                     model.logger.warning(msg)
                                 else:
                                     self.usebuiltinoutletplantcorrelation.value = True
-                                    msg = ('Provided plant outlet pressure outside of range 0-10000 kPa. '
+                                    msg = ('Provided plant outlet pressure outside of defined valid range. '
                                            'GEOPHIRES will calculate plant outlet pressure based on production '
                                            'wellhead pressure and surface equipment pressure drop of 10 psi')
                                     print(f'Warning: {msg}')

tests/examples/Fervo_Norbeck_Latimer_2024.out

@@ -6,8 +6,8 @@ Simulation Metadata
 ----------------------
  GEOPHIRES Version: 3.4.29
  Simulation Date: 2024-05-15
- Simulation Time:  10:41
- Calculation Time:      0.337 sec
+ Simulation Time:  11:23
+ Calculation Time:      0.353 sec
 
                            ***SUMMARY OF RESULTS***
 
@@ -52,7 +52,7 @@ Simulation Metadata
 
                          ***RESOURCE CHARACTERISTICS***
 
-      Maximum reservoir temperature:                        400.0 degC
+      Maximum reservoir temperature:                        300.0 degC
       Number of segments:                                     1
       Geothermal gradient:                                    0.0767 degC/m
 

tests/examples/Fervo_Norbeck_Latimer_2024.txt

@@ -25,6 +25,7 @@ Fracture Width, 160, per the paper (distance between wellbores assumes a dipole
 Fracture Separation, 10 meter
 Number of Segments,1
 Gradient 1, 76.74, per the paper
+Maximum Temperature, 300 degC, per the paper
 Plant Outlet Pressure, 1450 psi, per the paper
 
 *** Wellbore Parameters***

tests/examples/MC_Fervo_Norbeck_Latimer_2024.txt

@@ -0,0 +1,17 @@
+INPUT, Reservoir Density, uniform, 2650, 2900
+INPUT, Number of Fractures, uniform, 90, 110
+INPUT, Fracture Height, normal, 100, 30
+INPUT, Fracture Width, normal, 100, 30
+INPUT, Fracture Separation, uniform, 5, 15
+INPUT, Plant Outlet Pressure, uniform, 6895, 13790
+INPUT, Production Flow Rate per Well, triangular, 35, 41, 47
+INPUT, Injection Temperature, uniform, 24, 52
+INPUT, Ambient Temperature, normal, 15, 7
+INPUT, Surface Temperature, normal, 15, 7
+OUTPUT, Average Net Electricity Generation
+OUTPUT, Average Total Electricity Generation
+OUTPUT, Electricity breakeven price
+OUTPUT, Average Production Temperature
+OUTPUT, Average Pumping Power
+ITERATIONS, 5000
+HTML_PATH, D:\TEMP\MC_Fervo_Norbeck_latimer_2024.html