Continued work on Fervo model; Produces 3.35 MW (Paper says 2-3.5 MW). Average Production Temp is 174.2C (paper says 169C). Average Pumping Power is 0.00 MW (Paper says .5 to 1 MW).

Author: malcolm-dsider

src/geophires_x/WellBores.py

@@ -1196,11 +1196,13 @@ def Calculate(self, model: Model) -> None:
         # making all thr values available to your methods. but you had better have set all the parameters!
 
         # calculate the reservoir pressure as a function of time
-        self.production_reservoir_pressure.value = get_hydrostatic_pressure_kPa(model.reserv.Trock.value, model.reserv.Tsurf.value,
+        if self.usebuiltinhydrostaticpressurecorrelation:
+            self.production_reservoir_pressure.value = get_hydrostatic_pressure_kPa(model.reserv.Trock.value, model.reserv.Tsurf.value,
                                                                                 model.reserv.depth.quantity().to('m').magnitude,
                                                                                 model.reserv.averagegradient.value,
-                                                                                model.reserv.lithostatic_pressure()) if self.usebuiltinhydrostaticpressurecorrelation else self.Phydrostatic.quantity().to(
-            self.production_reservoir_pressure.CurrentUnits).magnitude
+                                                                                model.reserv.hydrostatic_pressure())
+        else:
+            self.production_reservoir_pressure.value = self.Phydrostatic.quantity().to(self.production_reservoir_pressure.CurrentUnits).magnitude
 
         self.production_reservoir_pressure.value = ReservoirPressurePredictor(model.surfaceplant.plant_lifetime.value,
                                                                               model.economics.timestepsperyear.value,

tests/examples/Fervo_Norbeck_Latimer_2024.txt

@@ -1,22 +1,38 @@
 --- Case Study: Fervo Energy Whitepaper
---- Norbeck, J.H. and Latimer, T.M. 2024.
+--- Norbeck J.H. and Latimer T.M. 2024.
 ---   "Commercial-Scale Demonstration of a First-of-a-Kind Enhanced Geothermal System."
 ---   https://doi.org/10.31223/X52X0B
+--- Drilling costs from:
+---- Kareem El-Sadi; Brittany Gierke; Elliot Howard; Christian Gradl; 2024
+---- "Review Of Drilling Performance In A Horizontal EGS Development"
+---- https://pangea.stanford.edu/ERE/db/GeoConf/papers/SGW/2024/Elsadi.pdf?t=1706897516
+----
+---- This case study is a commercial-scale demonstration of a first-of-a-kind enhanced geothermal system.
+---- The reservoir is a 7700-foot deep; 2800 kg/m3 density; 1000 J/kg/K heat capacity; 2.7 W/m/K thermal conductivity reservoir.
+---- The reservoir has 87 fractures; each 750 feet high; 1000 feet wide; and 30 feet apart.
+---- The reservoir has a 72.86 degree temperature gradient.
+---- The second paper figure 4 shows a drilling cost that has plateaued and about $400 per foot. GEOPHIRES cannot
+---- handle a drilling cost estimate on a per foot basis but drilling cost correlation of VERTICAL_LARGE_INT1 (10)
+---- has a cost of $397.12 per foot at 7700 feet so we will use that.
+---- NOTES:
+-------- The paper does not discuss stimulation costs which could be substantially more than the default values used here.
+-------- The paper does not discuss makeup water - how much what cost etc.
+-------- The paper does not discuss redrilling or rentry to refrac etc.
 
 ***Reservoir Parameters***
 Reservoir Model, 1
 Reservoir Volume Option,1,					---Should be 1 2 3 or 4. See manual for details.
-Reservoir Density, 2700
-Reservoir Depth, 2347 meter
+Reservoir Density, 2800, about the average for phyllite, quartzite, diorite, and granodiorite, per paper
+Reservoir Depth, 7700 feet
 Reservoir Heat Capacity,1000,					---[J/kg/K]
 Reservoir Thermal Conductivity, 2.7
-Number of Fractures, 100
+Number of Fractures, 87, between 75 and 100 per paper
 Fracture Shape, 4
-Fracture Height, 100
-Fracture Width, 991 meter
-Fracture Separation, 10 meter
+Fracture Height, 750 feet, per figure 14 of paper
+Fracture Width, 1000 feet, per figure 14 of paper
+Fracture Separation, 30 feet, per paper
 Number of Segments,1,						---[-]
-Gradient 1, 68
+Gradient 1, 72.86
 Drawdown Parameter,.00002,					---[kg/s/m2]
 Maximum Temperature, 500
 
@@ -27,13 +43,14 @@ Well Geometry Configuration, 4,       ----L
 Has Nonvertical Section, True
 Multilaterals Cased, True
 Number of Multilateral Sections, 2,  ---Two parallel horizontal sections
-Total Nonvertical Length, 990 meter
-Well Drilling Cost Correlation,3,				--- [-] Use built-in well drilling cost correlation #3 = vertical large diameter
-Horizontal Well Drilling Cost Correlation,4,				--- [-] Use built-in well drilling cost correlation #4 = deviated large diameter
-Production Flow Rate per Well, 63
-Production Well Diameter, 7
-Injection Well Diameter, 7
-Injection Temperature, 50
+Well Separation, 365 feet, per paper
+Total Nonvertical Length, 3250 feet
+Well Drilling Cost Correlation,10,			See above -this gets us about $400/foot
+Horizontal Well Drilling Cost Correlation,10,			See above -this gets us about $400/foot at 7700 feet
+Production Flow Rate per Well, 47.32, this is 750 gpm, range from paper was 650 gpm to 850 gpm (41-54 kg/s)
+Production Well Diameter, 7 inch
+Injection Well Diameter, 7 inch
+Injection Temperature, 100 degF, The paper says is is between 75 and 125 F
 Injection Wellbore Temperature Gain,3,				---[deg.C]
 Injectivity Index, 5,						---[kg/s/bar]
 Productivity Index, 5
@@ -43,13 +60,14 @@ Water Loss Fraction,0.02,						---[-]
 Maximum Drawdown,1,						---[-] no redrilling considered
 
 ***Surfaceplant Parameters***
-Ambient Temperature, 10.27,						---[deg.C]
-Surface Temperature, 10.27
+Ambient Temperature, 8 degC, per the paper, but the mean annual temperature in Milford is 10.27		---[deg.C]
+Surface Temperature, 8 degC, per the paper, but the mean annual temperature in Milford is 10.27		---[deg.C]
 End-Use Option,1,						---electricity as the main product
-Plant Lifetime, 30
-Power Plant Type,4,						---Double-Flash
+Plant Lifetime, 23, not referenced in paper but performance drops off substantially after 23 years due to heat mined
+Power Plant Type,2,   per the paper, it will be a supercritical ORC
 End-Use Efficiency Factor,.9,					---[-]
 Circulation Pump Efficiency,.80,				---[-]
+Plant Outlet Pressure, 1450 psi,  paper says 1000 psi to 2000 psi (1500 will exceed max value, so we use 1450)
 
 *** Economic Parameters***
 Economic Model,3,						---BICYCLE Levelized Cost Model