Please wait a minute...
Frontiers in Energy

ISSN 2095-1701

ISSN 2095-1698(Online)

CN 11-6017/TK

Postal Subscription Code 80-972

2018 Impact Factor: 1.701

Front. Energy    2019, Vol. 13 Issue (3) : 439-449    https://doi.org/10.1007/s11708-019-0612-4
RESEARCH ARTICLE
A numerical study of non-Darcy flow in EGS heat reservoirs during heat extraction
Wenjiong CAO1, Wenbo HUANG1, Guoling WEI2, Yunlong JIN2, Fangming JIANG1()
1. Laboratory of Advanced Energy Systems, Guangdong Key Laboratory of New and Renewable Energy Research and Development, CAS Key Laboratory of Renewable Energy, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, China
2. Guangdong Hydrogeology Battalion, Guangzhou 510510, China
 Download: PDF(1289 KB)   HTML
 Export: BibTeX | EndNote | Reference Manager | ProCite | RefWorks
Abstract

Underground non-Darcy fluid flow has been observed and investigated for decades in the petroleum industry. It is deduced by analogy that the fluid flow in enhanced geothermal system (EGS) heat reservoirs may also be in the non-Darcy regime under some conditions. In this paper, a transient 3D model was presented, taking into consideration the non-Darcy fluid flow in EGS heat reservoirs, to simulate the EGS long-term heat extraction process. Then, the non-Darcy flow behavior in water- and supercritical CO2 (SCCO2)-based EGSs was simulated and discussed. It is found that non-Darcy effects decrease the mass flow rate of the fluid injected and reduce the heat extraction rate of EGS as a flow resistance in addition to the Darcy resistance which is imposed to the seepage flow in EGS heat reservoirs. Compared with the water-EGS, the SCCO2-EGS are more prone to experiencing much stronger non-Darcy flow due to the much larger mobility of the SCCO2. The non-Darcy flow in SCCO2- EGSs may thus greatly reduce their heat extraction performance. Further, a criterion was analyzed and proposed to judge the onset of the non-Darcy flow in EGS heat reservoirs. The fluid flow rate and the initial thermal state of the reservoir were taken and the characteristic Forchheimer number of an EGS was calculated. If the calculated Forchheimer number is larger than 0.2, the fluid flow in EGS heat reservoirs experiences non-negligible non-Darcy flow characteristic.

Keywords enhanced geothermal system      non-Darcy flow      heat extraction      Reynolds number      Forchheimer number     
Corresponding Author(s): Fangming JIANG   
Online First Date: 25 February 2019    Issue Date: 04 September 2019
 Cite this article:   
Wenjiong CAO,Wenbo HUANG,Guoling WEI, et al. A numerical study of non-Darcy flow in EGS heat reservoirs during heat extraction[J]. Front. Energy, 2019, 13(3): 439-449.
 URL:  
https://academic.hep.com.cn/fie/EN/10.1007/s11708-019-0612-4
https://academic.hep.com.cn/fie/EN/Y2019/V13/I3/439
Fig.1  Model geometry and numerical mesh system
B correlation B value/cm-1
β=4.1 ×109 K1.5 [20] 1.3 × 108
β=1.82 ×108 K1.25ε0.75 [21] 3.2 × 108
β=4.8 ×1010K1.176 [22] 3.2 × 109
Tab.1  Correlations for the non-Darcy coefficient b
Group # Case # Working fluid b/cm-1
1 Case 1 H2O Darcy
Case 2 H2O 1.3 × 108
Case 3 SCCO2 Darcy
Case 4 SCCO2 1.3 × 108
2 Case 5 H2O 3.2 × 108
Case 6 H2O 3.2 × 109
Case 7 SCCO2 3.2 × 108
Case 8 SCCO2 3.2 × 109
Tab.2  Simulated cases
Fig.2  Fluid pressure distribution in the heat reservoir (The plots from left to right correspond to time instants of 1, 5, and 10 years, respectively)
Fig.3  Seepage speed distribution in the heat reservoir (The plots from left to right correspond to time instants: 1, 5, and 10 years, respectively)
Fig.4  Fluid temperature fields in the heat reservoir (The plots from left to right correspond to time instants of 1, 5, and 10 years, respectively)
Fig.5  Mass flow rate as a function of EGS operation time in water- and SCCO2-based EGS
Fig.6  Non-Darcy effect on the production temperature curve for water- and SCCO2-based EGSs
Fig.7  Non-Darcy effects on the real-time heat extraction rate versus EGS operation time for water- and SCCO2-based EGSs
Fig.8  Relative velocity reduction in the heat reservoir due to non-Darcy flow after 5 years of EGS operation
Case # Heat transfer fluid β/cm1 Foch γ Re Fo
Point A Point B Point C Point A Point B Point C Point A Point B Point C
2 H2O 1.3 × 108 0.03 2.0% 1.1% 2.1% 0.14 5.8 × 104 0.13 0.03 1.3 × 104 0.03
5 H2O 3.2 × 108 0.07 2.5% 1.7% 2.6% 0.08 3.3 × 104 0.07
6 H2O 3.2 × 109 0.74 13.1% 10.1% 14.9% 0.78 3.2 × 103 0.74
4 SCCO2 1.3 × 108 0.35 11.0% 7.2% 11.1% 2.25 5.4 × 103 1.51 0.51 1.2 × 103 0.35
7 SCCO2 3.2 × 108 0.85 21.1% 13.6% 21.8% 1.27 3.1 × 103 0.85
8 SCCO2 3.2 × 109 8.54 53.9% 49.7% 59.1% 12.70 3.1 × 102 8.54
Tab.3  Relative velocity reduction, Reynolds number, and Forchheimer number at the three monitoring points after 5 years of EGS operation
A Area/m2
cp Specific heat capacity/(J·kg -1·K-1)
Fo Forchheimer number
Foch Characteristic Forchheimer number
g Acceleration of gravity/(m·s -2)
ha Volumetric heat exchange coefficient/(W·m-3·K-1)
k Thermal conductivity/(W·m-1·K-1)
K Permeability/m2
p Pressure/Pa
Q Mass flow rate/(kg·s-1)
Re Reynolds number
t Time/s
SND Non-Darcy term
T Temperature/K
Tg Ground temperature/K
u Velocity vector/(m·s-1)
x, y, z Cartesian coordinates
Greek symbols
ρ Density/(kg·m-3)
ε Porosity
μ Dynamic viscosity/(Pa·s)
β Non-Darcy coefficient
γ Relative velocity reduction
Subscripts/superscripts
ch Characteristic
eff Effective
f Fluid
H2O Water
inj Injection
ini Initial
out Production, outlet
ref Reference
s Solid or rock
SCCO2 Supercritical carbon dioxide
  
1 J W Tester, B J Anderson, A S Batchelor , D D Blackwell , R Di Pippo , E M Drake, J Garnish, B Livesay , M C Moore , K Nichols , S Petty , M N Toksoz , R W Veatch. The Future of Geothermal Energy: Impact of Enhanced Geothermal Systems (EGS) on the United States in the 21st Century. DOE Contract DE-AC07-05ID14517 Final Report, Massachusetts Institute of Technology, 2006
2 Y Yang, H Yeh. Modeling heat extraction from hot dry rock in a multi-well system. Applied Thermal Engineering, 2009, 29(8–9): 1676–1681
https://doi.org/10.1016/j.applthermaleng.2008.07.020
3 P Olasolo, M C Juárez, M P Morales, S D´Amico, I A Liarte. Enhanced geothermal systems (EGS): a review. Renewable & Sustainable Energy Reviews, 2016, 56: 133–144
https://doi.org/10.1016/j.rser.2015.11.031
4 P Forchheimer. Water moving through soil. Zeitschrift des Vereines Deutscher Ingenieuer, 1901, 45: 1782–1788
5 R Gelet, B Loret, N Khalili. A thermo-hydro-mechanical coupled model in local thermal non-equilibrium for fractured HDR reservoir with double porosity. Journal of Geophysical Research, 2012, 117: B07205
https://doi.org/10.1029/2012JB009161
6 S Held, A Genter, T Kohl, T Kölbel, J Sausse, M Schoenball. Economic evaluation of geothermal reservoir performance through modeling the complexity of the operating EGS in Soultz-sous-Forêts. Geothermics, 2014, 51: 270–280
https://doi.org/10.1016/j.geothermics.2014.01.016
7 D W Brown. A hot dry rock geothermal energy concept utilizing supercritical CO2 instead of water. In: The 25th Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, CA, USA, 2000
8 K Pruess. On production behavior of enhanced geothermal systems with CO2 as working fluid. Applied Thermal Engineering, 2008, 49: 1446–1454
9 F M Jiang, L Luo, J L Chen. A novel three-dimensional transient model for subsurface heat exchange in enhanced geothermal systems. International Communications in Heat and Mass Transfer, 2013, 41: 57–62
https://doi.org/10.1016/j.icheatmasstransfer.2012.11.003
10 T Kohl, K F Evans, R J Hopkirk, R Jung, L Rybach. Observation and simulation of non-Darcian flow transients in fractured rock. Water Resources Research, 1997, 33(3): 407–418
https://doi.org/10.1029/96WR03495
11 J Zhang, H L Xing. Numerical modeling of non-Darcy flow in near-well region of a geothermal reservoir. Geothermics, 2012, 42: 78–86
https://doi.org/10.1016/j.geothermics.2011.11.002
12 M Haghshenas Fard, K Hooman, H T Chua. Numerical simulation of a supercritical CO2 geothermosiphon. International Communications in Heat and Mass Transfer, 2010, 37(10): 1447–1451
https://doi.org/10.1016/j.icheatmasstransfer.2010.08.016
13 C S Xu, P Dowd, Q Li. Carbon sequestration potential of the Habanero reservoir when carbon dioxide is used as the heat exchange fluid. Journal of Rock Mechanics and Geotechnical Engineering, 2016, 8: 50–59
14 A B Borgia, C M Oldenburg, R Zhang, L Pan, T M Daley, S Finsterle, T S Ramakrishnan. Simulations of CO2 injection into fractures and faults for improving their geophysical characterization at EGS sites. Geothermics, 2017, 69: 189–201
https://doi.org/10.1016/j.geothermics.2017.05.002
15 C L Wang, W L Cheng, Y L Nian, L Yang, B B Han, M H Liu. Simulation of heat extraction from CO2-based enhanced geothermal systems considering CO2 sequestration. Energy, 2018, 142: 157–167
https://doi.org/10.1016/j.energy.2017.09.139
16 W J Cao, W B Huang, F M Jiang. Numerical study on variable thermophysical properties of heat transfer fluid affecting EGS heat extraction. International Journal of Heat and Mass Transfer, 2016, 92: 1205–1217
https://doi.org/10.1016/j.ijheatmasstransfer.2015.09.081
17 J L Chen, F M Jiang. A numerical study of EGS heat extraction process based on a thermal non-equilibrium model for heat transfer in subsurface porous heat reservoir. Heat and Mass Transfer, 2016, 52(2): 255–267
https://doi.org/10.1007/s00231-015-1554-y
18 J L Chen, F M Jiang. Designing multi-well layout for enhanced geothermal system to better exploit hot dry rock geothermal energy. Renewable Energy, 2015, 74: 37–48
https://doi.org/10.1016/j.renene.2014.07.056
19 F M Jiang, J L Chen, W B Huang, L Luo. A three-dimensional transient model for EGS subsurface thermo-hydraulic process. Energy, 2014, 72: 300–310
https://doi.org/10.1016/j.energy.2014.05.038
20 T Friedel, H D Voigt. Investigation of non-Darcy flow in tight-gas reservoirs with fractured wells. Journal of Petroleum Science Engineering, 2006, 54(3-4): 112–128
https://doi.org/10.1016/j.petrol.2006.07.002
21 J D Janicek, D L Katz. Applications of unsteady state gas flow calculations. In: Research Conference on Flow of Natural Gas Reservoirs, University of Michigan, 1955
22 H Pascal, R G Quillian, J Kingston. Analysis of vertical fracture length and non-Darcy flow coefficient using variable rate tests. SPE 9348, 1980
23 Ф И Koтяxoв. The Physical Foundation of Reservoir. Beijing: Petroleum Industry Press, 1958
24 C Li. The research of natural gas small leakage & dispersion rule. Dissertation for the Doctoral Degree. Beijing: China University of Petroleum, 2011
25 H Saboorian-Jooybari, P Pourafshary. Significance of non-Darcy flow effect in fractured tight reservoirs. Journal of Natural Gas Science and Engineering, 2015, 24: 132–143
https://doi.org/10.1016/j.jngse.2015.03.003
26 Z Wen, G Huang, H Zhan. An analytical solution for non-Darcian flow in a confined aquifer using the power law function. Advances in Water Resources, 2008, 31(1): 44–55
https://doi.org/10.1016/j.advwatres.2007.06.002
27 Z Y Zhang, J Nemcik. Fluid flow regimes and nonlinear flow characteristics in deformable rock fractures. Journal of Hydrology (Amsterdam), 2013, 477: 139–151
https://doi.org/10.1016/j.jhydrol.2012.11.024
28 H C Brinkman. A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles. Applied Scientific Research, 1949, 1: 27–34
[1] Jian GUO, Wenjiong CAO, Yiwei WANG, Fangming JIANG. A novel flow-resistor network model for characterizing enhanced geothermal system heat reservoir[J]. Front. Energy, 2019, 13(1): 99-106.
[2] Yufei MAO, Liejin GUO, Bofeng BAI, Ximin ZHANG. Convective heat transfer in helical coils for constant-property and variable-property flows with high Reynolds numbers[J]. Front Energ Power Eng Chin, 2010, 4(4): 546-552.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed