• 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.26 no.3
• /
• pp.166-180
• /
• 2016

The present study numerically simulated the CO2 injection into the saline aquifer of CO2CRC Otway pilot project and the resulting hydrological-mechanical coupled process in the storage site by TOUGH-FLAC simulator. A three-dimensional numerical model was generated using the stochastic geological model which was established based on well log and core data. It was estimated that the CO2 injection of 30,000t over a period of 200 days increased the pressure near the injection point by 0.5 MPa at the most. The pressure increased rapidly and tended to approach a certain value at an early stage of the injection. The hydrological and mechanical behavior observed from the CO2 flow, effective stress change and stress-strength ratio revealed that the CO2 injection into the saline aquifer under the given condition would not have significant effects on the mechanical safety of the storage site and the hydrological state around the adjacent fault.

• 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.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.

• 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.26 no.1
• /
• pp.12-23
• /
• 2016

For a successful accomplishment of Carbon Capture Sequestration (CCS) projects, appropriate injection conditions should be designed and optimized for site specific geological conditions. In this study, we evaluated the effect of injection conditions such as injection temperature and injection rate on the geomechanical stability of CCS system in terms of TOUGH-FLAC simulator, which is one of the well-known T-H-M coupled analysis methods. The stability of the storage system was assessed by a shear slip potential of the pre-existing fractures both in a reservoir and caprock, expressed by mobilized friction angle and Mohr stress circle. We demonstrated that no tensile fracturing was induced even in the cold CO2 injection, where the injected CO2 temperature is much lower than that of the reservoir and tensile thermal stress is generated, but shear slip of the fractures in the reservoir may occur. We also conducted a scenario analysis by varying injected CO2 volume per unit time, and found out that it was when the injection rate was decreasing in a step-wise that showed the least potential of a shear slip.

• Tunnel and Underground Space
• /
• v.33 no.4
• /
• pp.265-280
• /
• 2023

Water leakage in tunnels is a defect that can affect tunnel stability and the ground movement by changing the stress and pore water pressure of the surrounding ground. Long-term or large-scale water leaks may lead to damage of tunnel structure and the surrounding environment, such as tunnel lining instability and ground surface settlement. The present study numerically investigated the effects of water leakage on the structural stability of a tunnel and the ground behavior. The tunnel was assumed to be under undrained conditions for preventing the inflow of the surrounding water and leaks occurred in the concrete lining after completion of the tunnel construction. A coupled hydro-mechanical analysis using a TOUGH-FLAC simulator developed in Python was conducted for assessing the leakage induced-behavior of the tunnel structure and ground under different conditions of the amount and location of water leak. Additionally, the effect of hydro-mechanical coupling terms on the results of coupled response was investigated and discussed.

• 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.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.