Skip to content

Commit 29a03c0

Browse files
Merge pull request #12 from softwareengineerprogrammer/vapor-pressure-calc-with-pressure-2
Calculate vapor_pressure_water_kPa with hydrostatic pressure
2 parents fc37843 + 935abd5 commit 29a03c0

11 files changed

+1114
-1118
lines changed

src/geophires_x/WellBores.py

Lines changed: 6 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -364,10 +364,7 @@ def ProdPressureDropAndPumpingPowerUsingIndexes(
364364
# Minimum production pump inlet pressure and minimum wellhead pressure
365365
Pminimum_kPa = vapor_pressure_water_kPa(
366366
Trock_degC,
367-
# TODO pass pressure https://github.com/NREL/GEOPHIRES-X/issues/118; either
368-
# pressure=model.reserv.lithostatic_pressure(),
369-
# or
370-
# pressure=quantity(Phydrostaticcalc_kPa, 'kPa'),
367+
pressure=quantity(Phydrostaticcalc_kPa, 'kPa'),
371368
) + Pexcess_kPa
372369

373370
if usebuiltinppwellheadcorrelation:
@@ -499,10 +496,7 @@ def InjPressureDropAndPumpingPowerUsingIndexes(
499496
# Minimum production pump inlet pressure and minimum wellhead pressure
500497
Pminimum_kPa = vapor_pressure_water_kPa(
501498
Trock_degC,
502-
# TODO pass pressure https://github.com/NREL/GEOPHIRES-X/issues/118; either
503-
# pressure=model.reserv.lithostatic_pressure(),
504-
# or
505-
# pressure=quantity(Phydrostaticcalc_kPa, 'kPa'),
499+
pressure=quantity(Phydrostaticcalc_kPa, 'kPa'),
506500
) + Pexcess_kPa
507501

508502
if usebuiltinppwellheadcorrelation:
@@ -1035,7 +1029,8 @@ def Calculate(self, model: Model) -> None:
10351029
model, self.Tinj.value, self.prodwellflowrate.value, self.injwelldiam.value, self.impedancemodelused.value,
10361030
model.reserv.depth.value, self.nprod.value, self.ninj.value, model.reserv.waterloss.value)
10371031

1038-
if self.impedancemodelused.value: # assumed everything stays liquid throughout, based on TARB in Geophires v1.2
1032+
if self.impedancemodelused.value:
1033+
# assumed everything stays liquid throughout, based on TARB in Geophires v1.2
10391034
self.DPOverall.value, self.PumpingPower.value, self.DPProdWell.value, self.DPReserv.value, self.DPBouyancy.value = \
10401035
ProdPressureDropsAndPumpingPowerUsingImpedenceModel(f3, vprod, self.rhowaterinj, self.rhowaterprod,
10411036
model.reserv.rhowater.value,
@@ -1052,7 +1047,8 @@ def Calculate(self, model: Model) -> None:
10521047
model.surfaceplant.pump_efficiency.value,
10531048
self.DPOverall.value)
10541049

1055-
else: # PI and II are used
1050+
else:
1051+
# PI and II are used
10561052
self.PumpingPower.value, self.PumpingPowerProd.value, self.DPProdWell.value, self.Pprodwellhead.value = \
10571053
ProdPressureDropAndPumpingPowerUsingIndexes(model,
10581054
self.productionwellpumping.value,

tests/example1_addons.csv

Lines changed: 399 additions & 399 deletions
Large diffs are not rendered by default.

tests/examples/example1.out

Lines changed: 85 additions & 85 deletions
Large diffs are not rendered by default.

tests/examples/example10_HP.out

Lines changed: 50 additions & 50 deletions
Original file line numberDiff line numberDiff line change
@@ -4,17 +4,17 @@
44

55
Simulation Metadata
66
----------------------
7-
GEOPHIRES Version: 3.4.1
7+
GEOPHIRES Version: 3.4.3
88
GEOPHIRES Build Date: 2022-06-30
9-
Simulation Date: 2024-02-14
10-
Simulation Time: 15:13
11-
Calculation Time: 0.152 sec
9+
Simulation Date: 2024-02-18
10+
Simulation Time: 09:03
11+
Calculation Time: 0.143 sec
1212

1313
***SUMMARY OF RESULTS***
1414

1515
End-Use Option: Direct-Use Heat
1616
Average Direct-Use Heat Production: 16.64 MW
17-
Direct-Use heat breakeven price: 14.49 USD/MMBTU
17+
Direct-Use heat breakeven price: 18.82 USD/MMBTU
1818
Number of production wells: 2
1919
Number of injection wells: 2
2020
Flowrate per production well: 70.0 kg/sec
@@ -29,10 +29,10 @@ Simulation Metadata
2929
Accrued financing during construction: 5.00
3030
Project lifetime: 30 yr
3131
Capacity factor: 90.0 %
32-
Project NPV: 5.39 MUSD
33-
Project IRR: 0.08 %
34-
Project VIR=PI=PIR: 1.17
35-
Project MOIC: 0.98
32+
Project NPV: 4.04 MUSD
33+
Project IRR: 0.07 %
34+
Project VIR=PI=PIR: 1.13
35+
Project MOIC: 0.92
3636

3737
***ENGINEERING PARAMETERS***
3838

@@ -67,8 +67,8 @@ Simulation Metadata
6767
Fracture separation: 80.00 meter
6868
Reservoir volume: 176000000 m**3
6969
Reservoir hydrostatic pressure: 20849.48 kPa
70-
Plant outlet pressure: 416.74 kPa
71-
Production wellhead pressure: 485.69 kPa
70+
Plant outlet pressure: 21125.23 kPa
71+
Production wellhead pressure: 21194.18 kPa
7272
Productivity Index: 10.00 kg/sec/bar
7373
Injectivity Index: 5.00 kg/sec/bar
7474
Reservoir density: 3000.00 kg/m**3
@@ -84,8 +84,8 @@ Simulation Metadata
8484
Initial Production Temperature: 104.5 degC
8585
Average Reservoir Heat Extraction: 11.89 MW
8686
Wellbore Heat Transmission Model = Constant Temperature Drop: 5.0 degC
87-
Average Injection Well Pump Pressure Drop: 1661.4 kPa
88-
Average Production Well Pump Pressure Drop: 830.3 kPa
87+
Average Injection Well Pump Pressure Drop: -19047.1 kPa
88+
Average Production Well Pump Pressure Drop: 21538.8 kPa
8989

9090

9191
***CAPITAL COSTS (M$)***
@@ -95,20 +95,20 @@ Simulation Metadata
9595
Stimulation costs: 3.02 MUSD
9696
Surface power plant costs: 7.70 MUSD
9797
of which Heat Pump Cost: 3.74 MUSD
98-
Field gathering system costs: 2.55 MUSD
99-
Total surface equipment costs: 10.25 MUSD
98+
Field gathering system costs: 3.73 MUSD
99+
Total surface equipment costs: 11.43 MUSD
100100
Exploration costs: 3.85 MUSD
101-
Total capital costs: 31.03 MUSD
101+
Total capital costs: 32.21 MUSD
102102

103103

104104
***OPERATING AND MAINTENANCE COSTS (M$/yr)***
105105

106-
Wellfield maintenance costs: 0.23 MUSD/yr
106+
Wellfield maintenance costs: 0.24 MUSD/yr
107107
Power plant maintenance costs: 0.31 MUSD/yr
108108
Water costs: 0.07 MUSD/yr
109-
Average Reservoir Pumping Cost: 0.25 MUSD/yr
109+
Average Reservoir Pumping Cost: 2.12 MUSD/yr
110110
Average Heat Pump Electricity Cost: 3.63 MUSD/yr
111-
Total operating and maintenance costs: 4.49 MUSD/yr
111+
Total operating and maintenance costs: 6.37 MUSD/yr
112112

113113

114114
***SURFACE EQUIPMENT SIMULATION RESULTS***
@@ -119,44 +119,44 @@ Simulation Metadata
119119
Initial Net Heat Production: 17.24 MW
120120
Average Annual Heat Production: 130.78 GWh
121121
Average Annual Heat Pump Electricity Use: 51.90 GWh/year
122-
Average Pumping Power: 0.45 MW
122+
Average Pumping Power: 3.85 MW
123123

124124
************************************************************
125125
* HEATING, COOLING AND/OR ELECTRICITY PRODUCTION PROFILE *
126126
************************************************************
127127
YEAR THERMAL GEOFLUID PUMP NET HEAT PUMP
128128
DRAWDOWN TEMPERATURE POWER HEAT ELECTRICITY USE
129129
(deg C) (MWe) (MWt) (MWe)
130-
0 1.0000 104.50 0.4443 17.2387 6.8408
131-
1 1.0000 104.50 0.4443 17.2387 6.8408
132-
2 1.0000 104.50 0.4443 17.2387 6.8408
133-
3 1.0000 104.50 0.4443 17.2387 6.8408
134-
4 1.0000 104.50 0.4443 17.2387 6.8407
135-
5 1.0000 104.50 0.4443 17.2380 6.8405
136-
6 1.0000 104.50 0.4443 17.2352 6.8394
137-
7 0.9999 104.49 0.4444 17.2283 6.8366
138-
8 0.9997 104.47 0.4444 17.2150 6.8314
139-
9 0.9995 104.44 0.4445 17.1937 6.8229
140-
10 0.9991 104.41 0.4446 17.1632 6.8108
141-
11 0.9986 104.36 0.4447 17.1231 6.7949
142-
12 0.9980 104.29 0.4448 17.0733 6.7751
143-
13 0.9973 104.22 0.4450 17.0143 6.7517
144-
14 0.9965 104.14 0.4452 16.9466 6.7249
145-
15 0.9956 104.04 0.4454 16.8711 6.6949
146-
16 0.9946 103.94 0.4457 16.7886 6.6621
147-
17 0.9936 103.83 0.4459 16.6999 6.6269
148-
18 0.9924 103.71 0.4462 16.6057 6.5896
149-
19 0.9913 103.59 0.4465 16.5070 6.5504
150-
20 0.9900 103.46 0.4468 16.4044 6.5097
151-
21 0.9888 103.33 0.4471 16.2984 6.4676
152-
22 0.9875 103.19 0.4474 16.1898 6.4245
153-
23 0.9862 103.05 0.4478 16.0791 6.3806
154-
24 0.9848 102.91 0.4481 15.9666 6.3360
155-
25 0.9835 102.77 0.4484 15.8529 6.2908
156-
26 0.9821 102.63 0.4488 15.7383 6.2453
157-
27 0.9807 102.48 0.4491 15.6230 6.1996
158-
28 0.9793 102.34 0.4495 15.5075 6.1538
159-
29 0.9780 102.20 0.4498 15.3918 6.1079
130+
0 1.0000 104.50 3.8462 17.2387 6.8408
131+
1 1.0000 104.50 3.8462 17.2387 6.8408
132+
2 1.0000 104.50 3.8462 17.2387 6.8408
133+
3 1.0000 104.50 3.8462 17.2387 6.8408
134+
4 1.0000 104.50 3.8462 17.2387 6.8407
135+
5 1.0000 104.50 3.8462 17.2380 6.8405
136+
6 1.0000 104.50 3.8462 17.2352 6.8394
137+
7 0.9999 104.49 3.8462 17.2283 6.8366
138+
8 0.9997 104.47 3.8462 17.2150 6.8314
139+
9 0.9995 104.44 3.8462 17.1937 6.8229
140+
10 0.9991 104.41 3.8462 17.1632 6.8108
141+
11 0.9986 104.36 3.8462 17.1231 6.7949
142+
12 0.9980 104.29 3.8462 17.0733 6.7751
143+
13 0.9973 104.22 3.8462 17.0143 6.7517
144+
14 0.9965 104.14 3.8462 16.9466 6.7249
145+
15 0.9956 104.04 3.8461 16.8711 6.6949
146+
16 0.9946 103.94 3.8461 16.7886 6.6621
147+
17 0.9936 103.83 3.8461 16.6999 6.6269
148+
18 0.9924 103.71 3.8461 16.6057 6.5896
149+
19 0.9913 103.59 3.8461 16.5070 6.5504
150+
20 0.9900 103.46 3.8460 16.4044 6.5097
151+
21 0.9888 103.33 3.8460 16.2984 6.4676
152+
22 0.9875 103.19 3.8460 16.1898 6.4245
153+
23 0.9862 103.05 3.8460 16.0791 6.3806
154+
24 0.9848 102.91 3.8459 15.9666 6.3360
155+
25 0.9835 102.77 3.8459 15.8529 6.2908
156+
26 0.9821 102.63 3.8459 15.7383 6.2453
157+
27 0.9807 102.48 3.8458 15.6230 6.1996
158+
28 0.9793 102.34 3.8458 15.5075 6.1538
159+
29 0.9780 102.20 3.8458 15.3918 6.1079
160160

161161

162162
*******************************************************************

0 commit comments

Comments
 (0)