Doing a Monte Carlo Simulation for Fervo. Also, extending the maximum…

… possible range for Plant Outlet Pressure
malcolm-dsider · May 15, 2024 · 0086c8f · 0086c8f
1 parent 4025de9
commit 0086c8f
Show file tree

Hide file tree

Showing 3 changed files with 22 additions and 5 deletions.
diff --git a/src/geophires_x/SurfacePlant.py b/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}')

diff --git a/tests/examples/Fervo_Norbeck_Latimer_2024.txt b/tests/examples/Fervo_Norbeck_Latimer_2024.txt
@@ -25,7 +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, 170 degC, per the paper
+Maximum Temperature, 300 degC, per the paper
 Plant Outlet Pressure, 1450 psi, per the paper
 
 *** Wellbore Parameters***

diff --git a/tests/examples/MC_Fervo_Norbeck_Latimer_2024.txt b/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