Computational modeling of convective seepage flow in fluid-saturated heterogeneous rocks: Steady-state approach

Abstract Convective seepage flow plays a key role not only in naturally forming mineral deposits and oil reservoirs, but also in the carbon-dioxide sequestration in fluid-saturated rocks. This paper presents a steady-state numerical solver, which is based on the finite element method (FEM) and progressive asymptotic approach procedure (PAAP), to solve steady-state convective seepage flow problems in fluid-saturated heterogeneous rocks. Although this kind of flow problem is commonly solved by using transient-state numerical solvers, the use of the proposed steady-state numerical solver, which is the main characteristic of this study, can have certain advantages. After the proposed steady-state numerical solver is verified by a benchmark testing problem, for which analytical solutions are available for comparison, it is further applied to simulate steady-state convective seepage flow in the Australia North West Shelf basin, which is composed of naturally-formed heterogeneous porous rocks. Through the simulated numerical results, it can be found that: (1) compared with the analytical solutions of the benchmark testing problem, the use of the steady-state numerical solver can produce correct and accurate numerical simulation results for dealing with convective seepage flow problems in fluid-saturated heterogeneous rocks. (2) The main advantage of using the steady-state numerical solver is to avoid the need of considering initial conditions of the problem, which are commonly difficult to be determined because the convective seepage flow within the upper crust is a past event in geology. (3) In the Australia North West Shelf basin, the convective seepage flow and temperature distribution patterns depend strongly on both the permeability contrast of the faults and the basin geometry (including the layer structures and fault locations) in the computational model.

[1]  Alison Ord,et al.  Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity , 2003 .

[2]  R. D. Gasser,et al.  Onset of Convection in a Porous Medium With Internal Heat Generation , 1976 .

[3]  Alison Ord,et al.  Investigating dynamic mechanisms of geological phenomena using methodology of computational geosciences: An example of equal-distant mineralization in a fault , 2008 .

[4]  Pierre Genthon,et al.  The onset of natural convection in vertical fault planes:consequences for the thermal regime in crystalline basementsand for heat recovery experiments , 2000 .

[5]  Stephan K. Matthäi,et al.  Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks , 2007 .

[6]  Mohammad E. Taslim,et al.  Thermal Instability of a Fluid-Saturated Porous Medium Bounded by Thin Fluid Layers , 1987 .

[7]  Jeff P. Raffensperger,et al.  The formation of unconformity-type uranium ore deposits; 2, Coupled hydrochemical modeling , 1995 .

[8]  P. Adler,et al.  Hydrothermal convection in and around mineralized fault zones: insights from two‐ and three‐dimensional numerical modeling applied to the Ashanti belt, Ghana , 2009 .

[9]  Roland N. Horne,et al.  On the evolution of thermal disturbances during natural convection in a porous medium , 1980, Journal of Fluid Mechanics.

[10]  Jianwen Yang,et al.  Three-dimensional numerical modeling of salinity variations in driving basin-scale ore-forming fluid flow: Example from Mount Isa Basin, northern Australia , 2010 .

[11]  Thomas Graf,et al.  Combined inter- and intra-fracture free convection in fracture networks embedded in a low-permeability matrix , 2015 .

[12]  A. Bejan,et al.  Convection in Porous Media , 1992 .

[13]  G. Lebon,et al.  A thermodynamical modelling of fluid flows through porous media: application to natural convection , 1986 .

[14]  Chongbin Zhao,et al.  Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below , 1997 .

[15]  E. Palm,et al.  Modelling of thermal convection in sedimentary basins and its relevance to diagenetic reactions , 1988 .

[16]  D. Bernard,et al.  Natural convection in a porous layer bounded by impervious domains: from numerical approaches to experimental realization , 1999 .

[17]  Alison Ord,et al.  Convective and Advective Heat Transfer in Geological Systems , 2008 .

[18]  Alison Ord,et al.  Numerical modelling of fluids mixing, heat transfer and non‐equilibrium redox chemical reactions in fluid‐saturated porous rocks , 2006 .

[19]  Stuart W. Bull,et al.  Basin-Scale Numerical Modeling to Test the Role of Buoyancy-Driven Fluid Flow and Heat Transfer in the Formation of Stratiform Zn-Pb-Ag Deposits in the Northern Mount Isa Basin , 2006 .

[20]  Massoud Kaviany,et al.  Thermal Convective Instabilities in a Porous Medium , 1984 .

[21]  Alison Ord,et al.  Fundamentals of Computational Geoscience: Numerical Methods and Algorithms , 2009 .

[22]  Farid Alavyoon,et al.  On natural convection in vertical porous enclosures due to prescribed fluxes of heat and mass at the vertical boundaries , 1993 .

[23]  A. A. Pek,et al.  Onset of Fault-Bounded Free Thermal Convection in a Fluid-Saturated Horizontal Permeable Porous Layer , 2015, Transport in Porous Media.

[24]  J. Z. Zhu,et al.  The finite element method , 1977 .

[25]  C. W. Horton,et al.  Convection Currents in a Porous Medium , 1945 .

[26]  Jianwen Yang.  Full 3-D numerical simulation of hydrothermal fluid flow in faulted sedimentary basins: Example of the McArthur Basin, Northern Australia , 2006 .

[27]  Chongbin Zhao,et al.  Finite element modeling of pore-fluid flow in the Dachang ore district, Guangxi, China: Implications for hydrothermal mineralization , 2011 .

[28]  Jean-Paul Caltagirone,et al.  Solutions and stability criteria of natural convective flow in an inclined porous layer , 1985, Journal of Fluid Mechanics.

[29]  Kenneth E. Torrance,et al.  Low Rayleigh Number Thermal Convection in a Vertical Cylinder Filled With Porous Materials and Heated From Below , 1982 .

[30]  A. A. Pek,et al.  Linked thermal convection of the basement and basinal fluids in formation of the unconformity‐related uranium deposits in the Athabasca Basin, Saskatchewan, Canada , 2016 .

[31]  Wolfgang Kinzelbach,et al.  Simulation of reactive processes related to biodegradation in aquifers: 1. Structure of the three-dimensional reactive transport model , 1998 .

[32]  Alison Ord,et al.  Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones. , 2004 .

[33]  O. M. Phillips,et al.  Flow and Reactions in Permeable Rocks , 1991 .