• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.335-358
• /
• 2020

We present a numerical model to simulate coupled hydro-mechanical behavior of fault using TOUGH-FLAC simulator. This study aims to develop a numerical method to estimate fluid injection-induced fault reactivation in low permeability rock and to access the relevant hydro-mechanical stability in rock as part of DECOVALEX-2019 Task B. A coupled fluid flow and mechanical interface model to explicitly represent a fault was suggested and validated from the applications to benchmark simulations and the field experiment at Mont Terri underground laboratory in Switzerland. The pressure build-up, hydraulic aperture evolution, displacement, and stress responses matched those obtained at the site, which indicates the capability of the model to appropriately capture the hydro-mechanical processes in rock fault.

• Tunnel and Underground Space
• /
• v.25 no.2
• /
• pp.155-167
• /
• 2015

The thermal-hydrological-mechanical (T-H-M) behavior of rock mass surrounding a high-temperature cavern thermal energy storage (CTES) operated for a period of 30 years has been investigated by TOUGH2-FLAC3D simulator. As a fundamental study for the development of prediction and control technologies for the environmental change and rock mass behavior associated with CTES, the key concerns were focused on the hydrological-thermal multiphase flow and the consequential mechanical behavior of the surrounding rock mass, where the insulator performance was not taken into account. In the present study, we considered a large-scale cylindrical cavern at shallow depth storing thermal energy of $350^{\circ}C$. The numerical results showed that the dominant heat transfer mechanism was the conduction in rock mass, and the mechanical behavior of rock mass was influenced by thermal factor (heat) more than hydrological factor (pressure). The effective stress redistribution, displacement and surface uplift caused by heating of rock and boiling of ground-water were discussed, and the potential of shear failure was quantitatively examined. Thermal expansion of rock mass led to the ground-surface uplift on the order of a few centimeters and the development of tensile stress above the storage cavern, increasing the potential of shear failure.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.2
• /
• pp.183-202
• /
• 2019

For design and performance assessment of a high-level radioactive waste (HLW) disposal system, it is necessary to understand the characteristics of coupled thermo-hydro-mechanical (THM) behavior. However, in previous studies for the Korean Reference HLW Disposal System (KRS), thermal analysis was performed to determine the spacing of disposal tunnels and interval of disposition holes without consideration of the coupled THM behavior. Therefore, in this study, TOUGH2-MP/FLAC3D is used to conduct THM modeling for performance assessment of the Korean Reference HLW Disposal System (KRS). The peak temperature remains below the temperature limit of $100^{\circ}C$ for the whole period. A rapid rise of temperature caused by decay heat occurs in the early years, and then temperature begins to decrease as decay heat from the waste decreases. The peak temperature at the bentonite buffer is around $96.2^{\circ}C$ after about 3 years, and peak temperature at the rockmass is $68.2^{\circ}C$ after about 17 years. Saturation of the bentonite block near the canister decreases in the early stage, because water evaporation occurs owing to temperature increase. Then, saturation of the bentonite buffer and backfill increases because of water intake from the rockmass, and bentonite buffer and backfill are fully saturated after about 266 years. The stress is calculated to investigate the effect of thermal stress and swelling pressure on the mechanical behavior of the rockmass. The calculated stress is compared to a spalling criterion and the Mohr-Coulumb criterion for investigation of potential failure. The stress at the rockmass remains below the spalling strength and Mohr-Coulumb criterion for the whole period. The methodology of using the TOUGH2-MP/FLAC3D simulator can be applied to predict the long-term behavior of the KRS under various conditions; these methods will be useful for the design and performance assessment of alternative concepts such as multi-layer and multi-canister concepts for geological spent fuel repositories.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Tunnel and Underground Space
• /
• v.33 no.5
• /
• pp.399-413
• /
• 2023

The present study conducted a numerical modeling of CO₂ injection at the Janggi Basin using the TOUGH-FLAC simulator, and examined the hydro-mechanical stability of the aquifer and the fault. Based on the site investigations and a 3D geological model of the target area, we simulated the injection of 32,850 tons of CO₂ over a 3-year period. The analysis of CO₂ plume with different values of the aquifer permeability revealed that assuming a permeability of 10^-14 m² the CO₂ plume exhibited a radial flow and reached the fault after 2 years and 9 months. Conversely, a higher permeability of 10^-13 m² resulted in predominant westward flow along the reservoir, with negligible impact on the fault. The pressure changes around the injection well remained below 0.6 MPa over the period, and the influence on the hydro-mechanical stability of the reservoir and fault was found to be insignificant.

• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.446-461
• /
• 2020

Numerical simulations of CIEMAT column test in Spain are performed to investigate the coupled thermo-hydro-mechanical (THM) behavior of MX80 bentonite pellets using TOUGH2-FLAC3D. The heater power and injection pressure of water in the numerical simulations are identical to those in the laboratory test. To investigate the applicability of the thermo-hydraulic (TH) model used in TOUGH2 code to prediction of the coupled TH behavior, the simulation results are compared with the observations of temperature and relative humidity with time. The tendencies of the coupled behavior observed in the test are well represented by the numerical models and the simulator in terms of temperature and relative humidity evolutions. Moreover, the performance of the models for the reproduction and prediction of the coupled TH behavior is globally satisfactory compared with the observations. However, the calculated stress change is relatively small and slow due to the limitations of the simple elastic and swelling pressure model used in numerical simulations. It seems that the two models are insufficient to realistically reproduce the complex coupled THM behavior in the bentonite pellets.

• Tunnel and Underground Space
• /
• v.25 no.2
• /
• pp.168-185
• /
• 2015

The thermal-hydrological-mechanical (T-H-M) behavior of rock mass surrounding a large-scale high-temperature cavern thermal energy storage (CTES) at a shallow depth has been investigated, and the effects of hydrological conditions such as water table and rock permeability on the behavior have been examined. The liquid saturation of ground water around a storage cavern may have a small impact on the overall heat transfer and mechanical behavior of surrounding rock mass for a relatively low rock permeability of $10^{-17}m^2$. In terms of the distributions of temperature, stress and displacement of the surrounding rock mass, the results expected from the simulation with the cavern below the water table were almost identical to that obtained from the simulation with the cavern in the unsaturated zone. The heat transfer in the rock mass with reasonable permeability ${\leq}10^{-15}m^2$ was dominated by the conduction. In the simulation with rock permeability of $10^{-12}m^2$, however, the convective heat transfer by ground-water was dominant, accompanying the upward heat flow to near-ground surface. The temperature and pressure around a storage cavern showed different distributions according to the rock permeability, as a result of the complex coupled processes such as the heat transfer by multi-phase flow and the evaporation of ground-water.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.221-242
• /
• 2021

A numerical study of the performance assesment of coupled thermo-hydro-mechanical (THM) processes in improved Korean reference disposal system (KRS⁺) for high-level radioactive waste is conducted using TOUGH2-MP/FLAC3D simulator. Decay heat from high-level radioactive waste increases the temperature of the repository, and it decreases as decay heat is reduced. The maximum temperature of the repository is below a maximum temperature criterion of 100℃. Saturation of bentonite buffer adjacent to the canister is initially reduced due to pore water evaporation induced by temperature increase. Bentonite buffer is saturated 250 years after the disposal of high-level radioactive waste by inflow of groundwater from the surrounding rock mass. Initial saturation of rock mass decreases as groundwater in rock mass is moved to bentnonite buffer by suction, but rock mass is saturated after inflow of groundwater from the far-field area. Stress changes at rock mass are compared to the Mohr-Coulomb failure criterion and the spalling strength in order to investigate the potential rock failure by thermal stress and swelling pressure. Additional simulations are conducted with the reduced spacing of deposition holes. The maximum temperature of bentonite buffer exceeds 100℃ as deposition hole spacing is smaller than 5.5 m. However, temperature of about 56.1% volume of bentonite buffer is below 90℃. The methodology of numerical modeling used in this study can be applied to the performance assessment of coupled THM processes for high-level radioactive waste repositories with various input parameters and geological conditions such as site-specific stress models and geothermal gradients.