A synthesis of approaches for modelling coupled thermal–hydraulic–mechanical–chemical processes in a single novaculite fracture experiment

Abstract The geological formation immediately surrounding a nuclear waste disposal facility has the potential to undergo a complex set of physical and chemical processes starting from construction and continuing many years after closure. The DECOVALEX project (DEvelopment of COupled models and their VALidation against EXperiments) was established and maintained by a variety of waste management organisations, regulators and research organisations to help improve capabilities in experimental interpretation, numerical modelling and blind prediction of complex coupled systems. In the present round of DECOVALEX (D-2015), one component of Task C1 has considered the detailed experimental work of Yasuhara et al. (Earth Planet Sci Lett 244:186–200, 2006), wherein a single artificial fracture in novaculite (micro- or crypto-crystalline quartz) is subject to variable fluid flows, mechanical confining pressure and different applied temperatures. This paper presents a synthesis of the completed work of six separate research teams. A range of approaches are presented including 2D and 3D high-resolution coupled thermo–hydro–mechanical–chemical models. The results of the work show that while good, physically plausible representations of the experiment can be obtained using a range of approaches, there is considerable uncertainty in the relative importance of the various processes, and that the parameterisation of these processes can be closely linked to the interpretation of the fracture surface topography at different spatial scales.

[1]  Xia-Ting Feng,et al.  Coupled THM processes in EDZ of crystalline rocks using an elasto-plastic cellular automaton , 2009 .

[2]  B. Berkowitz Characterizing flow and transport in fractured geological media: A review , 2002 .

[3]  F. Schmitt,et al.  Scaling invariance of crack surfaces , 1995 .

[4]  H. Eyring The Activated Complex in Chemical Reactions , 1935 .

[5]  D. Elsworth,et al.  Evolution of permeability in a natural fracture: Significant role of pressure solution , 2004 .

[6]  S. Norris,et al.  Modelling reactions between alkaline fluids and fractured rock: The Maqarin natural analogue , 2016 .

[7]  Gang Xu,et al.  FEM Analyses for T-H-M-M Coupling Processes in Dual-Porosity Rock Mass under Stress Corrosion and Pressure Solution , 2012, J. Appl. Math..

[8]  I. Neretnieks Stress‐mediated closing of fractures: Impact of matrix diffusion , 2014 .

[9]  C. C. Pain,et al.  A high‐order TVD transport method for hybrid meshes on complex geological geometry , 2005 .

[10]  B. Atkinson Fracture Mechanics of Rock , 1987 .

[11]  H. Zheng,et al.  An approach for simulating the THMC process in single novaculite fracture using EPCA , 2016, Environmental Earth Sciences.

[12]  Joshua Taron,et al.  Numerical simulation of thermal-hydrologic-mechanical-chemical processes in deformable, fractured porous media , 2009 .

[13]  T. A. Jones,et al.  Surface measurement and fractal characterization of naturally fractured rocks , 1992 .

[14]  Patricia M. Dove,et al.  Reply to Comment on “Kinetics of quartz dissolution in electrolyte solutions using a hydrothermal mixed flow reactor” , 1990 .

[15]  Special I Ssue,et al.  OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical (THM/C) processes in porous media , 2012 .

[16]  Gudmundur S. Bodvarsson,et al.  Lubrication theory analysis of the permeability of rough-walled fractures , 1991 .

[17]  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 .

[18]  S. Norris,et al.  Evaluating the importance of different coupled thermal, hydraulic, mechanical and chemical process simulations during fluid flow experiments in fractured novaculite and fractured granite , 2016, Environmental Earth Sciences.

[19]  Olaf Kolditz,et al.  Development of approaches for modelling coupled thermal–hydraulic–mechanical–chemical processes in single granite fracture experiments , 2016, Environmental Earth Sciences.

[20]  Phillip M. Halleck,et al.  Permeability reduction of a natural fracture under net dissolution by hydrothermal fluids , 2003 .

[21]  Patricia M. Dove,et al.  Geochemical controls on the kinetics of quartz fracture at subcritical tensile stresses , 1995 .

[22]  Stephen R. Brown,et al.  Closure of random elastic surfaces in contact , 1985 .

[23]  Johannes M. Miocic,et al.  Experimental investigation and hybrid numerical analytical hydraulic mechanical simulation of supercritical CO2 flowing through a natural fracture in caprock. , 2016 .

[24]  Xiating Feng,et al.  Modelling fluid flow through a single fracture with different contacts using cellular automata , 2011 .

[25]  J. Tester,et al.  Correlating quartz dissolution kinetics in pure water from 25 to 625°C , 1994 .

[26]  Ingvi Gunnarsson,et al.  Amorphous silica solubility and the thermodynamic properties of H4SiO°4 in the range of 0° to 350°C at Psat , 2000 .

[27]  A. Paluszny,et al.  Hydraulic sealing due to pressure solution contact zone growth in siliciclastic rock fractures , 2015 .

[28]  Eduardo F. D'Azevedo,et al.  HBGC123D: a high-performance computer model of coupled hydrogeological and biogeochemical processes , 2001 .

[29]  Hui Zhou,et al.  Simulation of the rock microfracturing process under uniaxial compression using an elasto-plastic cellular automaton , 2006 .

[30]  Shinichiro Nakashima,et al.  Temporal alteration of fracture permeability in granite under hydrothermal conditions and its interpretation by coupled chemo-mechanical model , 2011 .

[31]  Wenqing Wang,et al.  OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical/chemical (THM/C) processes in porous media , 2012, Environmental Earth Sciences.

[32]  E. Oelkers,et al.  SUPCRT92: a software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000 ° C , 1992 .

[33]  Christopher McDermott,et al.  Application of hybrid numerical and analytical solutions for the simulation of coupled thermal, hydraulic, mechanical and chemical processes during fluid flow through a fractured rock , 2015, Environmental Earth Sciences.

[34]  W. M. Heston,et al.  The Solubility of Amorphous Silica in Water , 1954 .

[35]  Yousif K. Kharaka,et al.  A Compilation of Rate Parameters of Water-Mineral Interaction Kinetics for Application to Geochemical Modeling , 2004 .

[36]  D. Elsworth,et al.  Compaction of a Rock Fracture Moderated by Competing Roles of Stress Corrosion and Pressure Solution , 2008 .

[37]  T. Cao,et al.  MODELING THE COUPLED HYDROGEOCHEMICAL AND MECHANICAL BEHAVIOR OF A SINGLE ARTIFICIAL FRACTURE IN NOVACULITE ROCK , 2015 .

[38]  J. B. Walsh,et al.  EFFECT OF PORE PRESSURE AND CONFINING PRESSURE ON FRACTURE PERMEABILITY , 1981 .

[39]  Yasuhiro Mitani,et al.  Evolution of fracture permeability through fluid–rock reaction under hydrothermal conditions , 2006 .

[40]  S. P. Neuman,et al.  Trends, prospects and challenges in quantifying flow and transport through fractured rocks , 2005 .

[41]  D. Elsworth,et al.  A numerical model simulating reactive transport and evolution of fracture permeability , 2006 .

[42]  Xiating Feng,et al.  Numerical study on coupled thermo-mechanical processes in Äspö Pillar Stability Experiment , 2013 .

[43]  J. Donald Rimstidt,et al.  QUARTZ SOLUBILITY AT LOW TEMPERATURES , 1997 .