Finite element modeling of convective pore-fluid flow in fluid-saturated porous rocks within upper crust: An overview

Abstract

Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.

Keywords

convective flow / steady-state approach / transient-state approach / numerical modeling / upper crust / porous rock

References

[1] HobbsB E, ZhangY, OrdA, ZhaoC. Application of coupled deformation, fluid flow, thermal and chemical modelling to predictive mineral exploration [J]. Journal of Geochemical Exploration, 2000, 69: 505-509

[2] GowP, UptonP, ZhaoC, HillK. Copper-gold mineralization in the New Guinea: Numerical modeling of collision, fluid flow and intrusion-related hydrothermal systems [J]. Australian Journal of Earth Sciences, 2002, 49(4): 753-771

[3] OrdA, HobbsB E, ZhangY, BroadbentG C, BrownM, WillettsG, Sorjonen-WardP, WalsheJ, ZhaoC. Geodynamic modelling of the Century deposit, Mt Isa Province, Queensland [J]. Australian Journal of Earth Sciences, 2002, 49(6): 1011-1039

[4] SchaubsP, ZhaoC. Numerical modelling of gold-deposit formation in the Bendigo-Ballarat zone, Victoria [J]. Australian Journal of Earth Sciences, 2002, 49(6): 1077-1096

[5] Sorjonen-WardP, ZhangY, ZhaoC. Numerical modelling of orogenic processes and mineralization in the south eastern part of the Yilgarn Craton, Western Australia [J]. Australian Journal of Earth Sciences, 2002, 49(6): 935-964

[6] Harcouet-MenouV, Guillou-FrottierL, BonnevilleA, AdlerP M, MourzenkoV. Hydrothermal convection in and around mineralized fault zones: Insights from two- and three-dimensional numerical modeling applied to the Ashanti belt, Ghana [J]. Geofluids, 2009, 9(2): 116-137

[7] JuM, DaiT, YangJ. Finite element modeling of pore-fluid flow in the Dachang ore district, Guangxi, China: Implications for hydrothermal mineralization [J]. Geoscience Frontiers, 2011, 2(3): 463-474

[8] NieldD A, BejanAConvection in porous media [M], 1992, New York, Springer-Verlag: 356

[9] ZhaoC, HobbsB E, OrdAConvective and advective heat transfer in geological systems [M], 2008, Berlin, Springer: 255

[10] LapwoodE R. Convection of a fluid in a porous medium [J]. Proceedings of the Cambridge Philosophical Society, 1948, 44(4): 508-521

[11] HortonC W, RogersF T. Convection currents in a porous medium [J]. Journal of Applied Physics, 1945, 16(6): 367-370

[12] TurcotteD L, SchubertGGeodynamics: Applications of continuum physics to geological problems [M], 1982, New York, John Wiley & Sons

[13] GasserR D, KazimiM S. Onset of convection in a porous medium with internal heat generation [J]. ASME Journal of Heat Transfer, 1976, 98(1): 49-54

[14] HorneR N, CaltagironeJ P. On the evaluation of thermal disturbances during natural convection in a porous medium [J]. Journal of Fluid Mechanics, 1980, 100(2): 385-395

[15] BauH H, TorranceK E. Low Rayleigh number thermal convection in a vertical cylinder filled with porous materials and heated from below [J]. ASME Journal of Heat Transfer, 1982, 104(1): 166-172

[16] KavianyM. Thermal convective instabilities in a porous medium [J]. ASME Journal of Heat Transfer, 1984, 106(1): 137-142

[17] CaltagironeJ P, BoriesS. Solutions and stability criteria of natural convective flow in an inclined porous layer [J]. Journal of Fluid Mechanics, 1985, 155: 267-287

[18] LebonG, ClootA. A thermodynamical modeling of fluid flows through porous media: Application to natural convection [J]. International Journal of Heat and Mass Transfer, 1986, 29(3): 381-390

[19] PillatsisG, TaslimM E, NarusawaU. Thermal instability of a fluid-saturated porous medium bounded by thin fluid layers [J]. ASME Journal of Heat Transfer, 1987, 109(3): 677-682

[20] BjorlykkeK, MoA, PalmE. Modelling of thermal convection in sedimentary basins and its relevance to diagenetic reactions [J]. Marine and Petroleum Geology, 1988, 5(4): 338-351

[21] AlavyoonF. On natural convection in vertical porous enclosures due to prescribed fluxes of heat and mass at the vertical boundaries [J]. International Journal of Heat and Mass Transfer, 1993, 36(10): 2479-2498

[22] ChevalierS, BernardD, JolyN. Natural convection in a porous layer bounded by impervious domains: From numerical approaches to experimental realization [J]. International Journal of Heat and Mass Transfer, 1999, 42(4): 581-597

[23] TournierC, GenthonP, RabinowiczM. The onset of natural convection in vertical fault planes: Consequences for the thermal regime in crystalline basements and for heat recovery experiments [J]. Geophysical Journal International, 2000, 140(3): 500-508

[24] LinG, HobbsB E, OrdA, MuhlhausH B. Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity [J]. Communications in Numerical Methods in Engineering, 2003, 19(10): 787-799

[25] PHILLIPS O M. Flow and reactions in permeable rocks [M]. Cambridge: Cambridge University Press, 1991: 286.

[26] RaffenspergerJ P, GarvenG. The formation of unconformity-type uranium ore deposits: Coupled hydrochemical modeling [J]. American Journal of Science, 1995, 295(6): 639-696

[27] ZhaoC, MuhlhausH B, HobbsB E. Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1997, 21(12): 863-881

[28] OrdA, PengS, MuhlhausH B, LiuL. Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones [J]. Tectonophysics, 2004, 387(1–4): 47-64

[29] HornbyP, OrdA, PengS. Numerical modelling of fluids mixing, heat transfer and non-equilibrium redox chemical reactions in fluid-saturated porous rocks [J]. International Journal for Numerical Methods in Engineering, 2006, 66(7): 1061-1078

[30] SchaferD, SchaferW, KinzelbachW. Simulation of reactive processes related to biodegradation in aquifers: 1. Structure of the three-dimensional reactive transport model [J]. Journal of Contaminant Hydrology, 1998, 31(12): 167-186

[31] YangJ W. Full 3D numerical simulation of hydrothermal fluid flow in faulted sedimentary basins: Example of the Mcarthur basin, northern Australia [J]. Journal of Geochemical Exploration, 2006, 89(1–3): 440-444

[32] YangJ W, FengZ, LuoX, ChenY. Three-dimensional numerical modeling of salinity variations in driving basin-scale ore-forming fluid flow: Example from Mount Isa Basin, northern Australia [J]. Journal of Geochemical Exploration, 2010, 106(1–3): 236-243

[33] YangJ W, LargeR, BullS, ScottD. Basin-scale numerical modelling to test the role of buoyancy driven fluid flow and heat transport in the formation of stratiform Zn-Pb-Ag deposits in the northern Mt Isa basin [J]. Economic Geology, 2006, 101(6): 1275-1292

[34] PalusznyA, MatthaiS K, HohmeyerM. Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks [J]. Geofluids, 2007, 7(2): 186-208

[35] Alt-EppingP, ZhaoC. Reactive mass transport modeling of a three-dimensional vertical fault zone with a finger-like convective flow regime [J]. Journal of Geochemical Exploration, 2010, 106(1–3): 8-23

[36] ZhaoC, HobbsB E, OrdA. Investigating dynamic mechanisms of geological phenomena using methodology of computational geosciences: An example of equal-distant mineralization in a fault [J]. Science in China Series D: Earth Sciences, 2008, 51(7): 947-954

[37] ZhaoC, HobbsB E, OrdAFundamentals of computational geoscience: Numerical methods and algorithms [M], 2009, Berlin, Springer: 285

[38] DierschH J GFEFLOW Reference manual [M], 2002, Berlin, Wasy GmbH: 225

[39] ZhaoC, ReidL B, Regenauer-LiebK. Some fundamental issues in computational hydrodynamics of mineralization: A review [J]. Journal of Geochemical Exploration, 2012, 112: 21-34

[40] ZienkiewiczO CThe finite element method [M], 1977, London, McGraw-Hill: 586

[41] KuhnM, DobertbF, GessnerK. Numerical investigation of the effect of heterogeneous permeability distributions on free convection in the hydrothermal system at Mount Isa, Australia [J]. Earth and Planetary Science Letters, 2006, 244(34): 655-671

[42] ZhaoC, HobbsB E, OrdA. Modeling of mountain topography effects on hydrothermal Pb-Zn mineralization patterns: Generic model approach [J]. Journal of Geochemical Exploration, 2018, 190: 400-410

[43] SchaubsP, HobbsB E. Acquisition of spatially-distributed geochemical data in geoinformatics: Computational simulation approach [J]. Journal of Geochemical Exploration, 2016, 164: 18-27

[44] ZhaoC, SchaubsP, HobbsB E. Computational simulation of seepage instability problems in fluid-saturated porous rocks: Potential dynamic mechanisms for controlling mineralization patterns [J]. Ore Geology Reviews, 2016, 79: 180-188

[45] MalkovskyV I, PekA A. Onset of fault-bounded free thermal convection in a fluid-saturated horizontal permeable porous layer [J]. Transport in Porous Media, 2015, 110(1): 25-39

[46] VujevicK, GrafT. Combined inter- and intra-fracture free convection in fracture networks embedded in a low-permeability matrix [J]. Advances in Water Resources, 2015, 84: 52-63

[47] PekA A, MalkovskyV I. Linked thermal convection of the basement and basinal fluids in formation of the unconformity-related uranium deposits in the Athabasca Basin, Saskatchewan, Canada [J]. Geofluids, 2016, 16(5): 925-940

[48] ZhaoC. Advances in numerical algorithms and methods in computational geosciences with modeling characteristics of multiple physical and chemical processes [J]. Science China Technological Sciences, 2015, 58(5): 783-795

[49] ZhaoC, HobbsB, OrdA. A new alternative approach for investigating acidization dissolution front propagation in fluid-saturated rocks [J]. Science China Technological Sciences, 2017, 60(8): 1197-1210

[50] ZhaoC, SchaubsP, HobbsB. Effects of porosity heterogeneity on chemical dissolution-front instability in fluid-saturated rocks [J]. Journal of Central South University, 2017, 24(3): 720-725

[51] ZhaoC, HobbsB, OrdA. A unified theory for sharp dissolution front propagation in chemical dissolution of fluid-saturated porous rocks [J]. Science China Technological Sciences, 2019, 62(1): 163-174

[52] ZhaoC, HobbsB, OrdA. Effects of different numerical algorithms on simulation of chemical dissolution-front instability in fluid-saturated porous rocks [J]. Journal of Central South University, 2018, 25(8): 1966-1975