In many subsurface engineering geoscience applications the impact of thermal, hydraulic, mechanical and chemical (THMC) processes needs to be evaluated. Coupled process models require solution of the partial differential equations describing energy or mass balance. Ignoring the coupling of these processes can lead to a significant oversimplification which may not adequately represent the systems being modelled. Incorporation of coupled processes and associated phenomena inevitably leads to numerical stability issues due to very different scales in terms of spatial distribution, time and parametrical heterogeneity. One approach to simplify the computational demands is to integrate analytical and physical models into standard numerical modelling techniques (in this case finite elements), effectively adding sub-grid scale and sub-time scale information to the model. We present such an approach for the simulation of fluid flow through a fracture validated against experimental data and cross comparison with results of other modelling teams within the DECOVALEX 2015 (development of coupled models and their validation against experiments) project (http:// www.decovalex.org). By replacing the mechanical behaviour and chemical transport processes with physical models, and by utilising the static nature of the temperature changes, only the hydraulic system required numerical solution in a highly coupled problem. Physical models for fracture closure due to pressure solution, fracture opening due to chemical dissolution, the development of channel flow and a change in the reactive transport characteristics with time were implemented and are described here. The main features of the experimental data could be replicated, although lying outside of the parameter range suggested by the literature. Comparison with other teams using different modelling approaches indicated internal consistency.

No 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 cryptocrystalline 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-hydromechanical-chemical models. The results of the work show that while good, physically plausible representations of the experiment can be obtained using a

No abstract

In order to establish sustainable heat loading (heat removal and storage) in abandoned flooded mine workings it is important to understand the geomechanical impact of the cyclical heat loading caused by fluid injection and extraction. This is particularly important where significantly more thermal loading is planned than naturally occurs. A simple calculation shows that the sustainable geothermal heat flux from abandoned coal mines can provide less than a tenth of Scotland's annual domestic heating demand. Any heat removal greater than the natural heat flux will lead to heat mining unless heat storage options are also considered.As a first step, a steady-state, fully saturated, 2D coupled hydromechanical model of a generalized section of pillar-and-stall workings has been created. Mine water rebound was modelled by increasing the hydrostatic pressure sequentially, in line with monitored mine water-level data from Midlothian, Scotland. The modelled uplift to water-level rise ratio of 1.4 mm m −1 is of the same order of magnitude (1 mm m −1 ) as that observed through interferometric synthetic aperture radar (InSAR) data in the coalfield due to mine water rebound. The modelled magnitude of shear stress at the pillar corners, as a result of horizontal and vertical displacement, is shown to increase linearly with water level. Mine heat systems are expected to cause smaller changes in pressure than those modelled but the results provide initial implications on the potential geomechanical impacts of mine water heat schemes which abstract or inject water and heat into pillar-and-stall coal mine workings.

No abstract

In many subsurface engineering geoscience applications the impact of thermal, hydraulic, mechanical and chemical (THMC) processes needs to be evaluated. Coupled process models require solution of the partial differential equations describing energy or mass balance. Ignoring the coupling of these processes can lead to a significant oversimplification which may not adequately represent the systems being modelled. Incorporation of coupled processes and associated phenomena inevitably leads to numerical stability issues due to very different scales in terms of spatial distribution, time and parametrical heterogeneity. One approach to simplify the computational demands is to integrate analytical and physical models into standard numerical modelling techniques (in this case finite elements), effectively adding sub-grid scale and sub-time scale information to the model. We present such an approach for the simulation of fluid flow through a fracture validated against experimental data and cross comparison with results of other modelling teams within the DECOVALEX 2015 (development of coupled models and their validation against experiments) project (http:// www.decovalex.org). By replacing the mechanical behaviour and chemical transport processes with physical models, and by utilising the static nature of the temperature changes, only the hydraulic system required numerical solution in a highly coupled problem. Physical models for fracture closure due to pressure solution, fracture opening due to chemical dissolution, the development of channel flow and a change in the reactive transport characteristics with time were implemented and are described here. The main features of the experimental data could be replicated, although lying outside of the parameter range suggested by the literature. Comparison with other teams using different modelling approaches indicated internal consistency.

A comparative modelling exercise involving several independent teams from the DECOVALEX-2015 project is presented in this paper. The exercise is based on various laboratory experiments that have been carried out in the framework of a French research program called SEALEX and conducted by the IRSN. The program focuses on the long-term performance of swelling clay-based sealing systems that provide an important contribution to the safety of underground nuclear waste disposal facilities. A number of materials are being considered in the sealing systems; the current work focuses on a 70/30 MX80 bentonite-sand mixture compacted at dry densities between 1.67 Mg/m 3 and 1.97 Mg/m 3. The improved understanding of the full set of hydro-mechanical processes affecting the behaviour of an in-situ sealing system requires both experiments ranging from small-scale laboratory tests to full-scale field emplacement studies, and coupled hydro-mechanical models that are able to explain the observations in the experiments. The approach was to build models of increasing complexity starting for the simplest laboratory experiments and building towards the full-scale in situ experiments. Following this approach, two sets of small-scale laboratory experiments have been performed and modelled. The first set of experiments involve characterising the hydro-mechanical behaviour of the bentonite-sand mixture by means of (i) water retention tests under both constant volume and free swell conditions, (ii) infiltration test under constant volume condition, and (iii) swelling and compression tests under suction control conditions. The second, more complex, experiment is a 1/10th scale mock-up of a larger scale in situ experiment. Modelling of the full-scale experiment is described in a companion paper. A number of independent teams have worked towards modelling these experiments using different conceptual models, codes, and input parameters. Their results are compared and discussed. This exercise has enabled an improved modelling of the bentonite-sand mixture behaviour, in particular accounting for the dependence of its retention curve on the dry density. Moreover, it has shown the importance of the technological voids on the short-term behaviour of the sealing system. .