• Tunnel and Underground Space
• /
• v.34 no.2
• /
• pp.127-142
• /
• 2024

The robustness of a numerical method means that its computational performance is maintained under various modeling conditions. New numerical methods or codes need to be assessed for robustness through benchmark testing. The TOUGH-FLAC modeling approach has been applied to various fields such as subsurface carbon dioxide storage, geological disposal of spent nuclear fuel, and geothermal development both domestically and internationally, and the modeling validity has been examined by comparing the results with experimental measurements and other numerical codes. In the present study, a benchmark test of the TOUGH-FLAC approach was performed based on a coupled thermal-hydro-mechanical behavior problem with an analytical solution. The analytical solution is related to the temperature, pore water pressure, and mechanical behavior of a fully saturated porous medium that is subjected to a point heat source. The robustness of the TOUGH-FLAC approach was evaluated by comparing the analytical solution with the results of numerical simulation. Additionally, the effects of thermal-hydro-mechanical coupling terms, fluid phase change, and timestep on the computation of coupled behavior were investigated.

• 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.30 no.1
• /
• pp.39-62
• /
• 2020

In this study, Barcelona Basic Model (BBM) was implemented into TOUGH2-MP/FLAC3D for the numerical analysis of coupled thermo-hydro-mechanical (THM) behavior of unsaturated soils and the prediction of long-term behaviors. Similar to the methodology described in a previous study for the implementation of BBM into TOUGH-FLAC, the User Defined Model (UDM) of FLAC based on the Modified Cam Clay Model (MCCM) and the FISH function of FLAC3D were used to extend the existing MCCM module in FLAC3D for the implementation of BBM into TOUGH2-MP/FLAC3D. In the developed BBM module in TOUGH2-MP/FLAC3D, the plastic strains due to change in suction increase (SI) in addition to mean effective stress are calculated. In addition to loading-collapse (LC) yield surface, suction increase (SI) yield surface is changed by hardening rules in the developed BBM module. Several numerical simulations were conducted to verify and validate the implementation of BBM: using an example presented in the FLAC3D manual for the standard MCCM, simulation results using COMSOL, and experimental data presented in SKB Reports. In addition, the developed BBM analysis module was validated by simultaneously performing a series of modeling tests that were performed for the validation of the Quick tools developed for the purpose of effectively deriving BBM parameters, and by comparing the Quick tools and Code_Bright results reported in a previous study.

• 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.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.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.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.21 no.5
• /
• pp.394-405
• /
• 2011

In this paper, we applied coupled non-isothermal, multiphase fluid flow and geomechanical numerical modeling using TOUGH-FLAC coupled analysis to study the complex thermodynamic and geomechanical performance of underground lined rock caverns (LRC) for compressed air energy storage (CAES). Mechanical stress in concrete linings as well as pressure and temperature within a storage cavern were examined during initial and long-term operation of the storage cavern for CAES. Our geomechanical analysis showed that effective stresses could decrease due to air penetration pressure, and tangential tensile stress could develop in the linings as a result of the air pressure exerted on the inner surface of the lining, which would result in tensile fracturing. According to the simulation in which the tensile tangential stresses resulted in radial cracks, increment of linings' permeability and air leakage though the linings, tensile fracturing occurred at the top and at the side wall of the cavern, and the permeability could increase to $5.0{\times}10^{-13}m^2$ from initially prescribed $10{\times}10^{-20}m^2$. However, this air leakage was minor (about 0.02% of the daily air injection rate) and did not significantly impact the overall storage pressure that was kept constant thanks to sufficiently air tight surrounding rocks, which supports the validity of the concrete-lined underground caverns for CAES.

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