• Geomechanics and Engineering
• /
• v.16 no.4
• /
• pp.363-373
• /
• 2018

The evaluation of Thermo-Hydro-Mechanical (THM) coupling behavior is important for the development of underground space for various purposes. For a high-level radioactive waste repository excavated in a deep underground rock mass, the accurate prediction of the complex THM behavior is essential for the long-term safety and stability assessment. In order to develop reliable THM analysis techniques effectively, an international cooperation project, Development of Coupled models and their Validation against Experiments (DECOVALEX), was carried out. In DECOVALEX-2015 Task B2, the in situ THM experiment that was conducted at Horonobe Underground Research Laboratory(URL) by Japan Atomic Energy Agency (JAEA), was modeled by the research teams from the participating countries. In this study, a THM coupling technique that combined TOUGH2 and FLAC3D was developed and applied to the THM analysis for the in situ experiment, in which rock, buffer, backfill, sand, and heater were installed. With the assistance of an artificial neural network, the boundary conditions for the experiment could be adequately implemented in the modeling. The thermal, hydraulic, and mechanical results from the modeling were compared with the measurements from the in situ THM experiment. The predicted buffer temperature from the THM modelling was about $10^{\circ}C$ higher than measurement near by the overpack. At the other locations far from the overpack, modelling predicted slightly lower temperature than measurement. Even though the magnitude of pressure from the modeling was different from the measurements, the general trends of the variation with time were found to be similar.

• 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.26 no.2
• /
• pp.59-67
• /
• 2016

In this paper, we perform a numerical analysis to evaluate the potential of fault reactivation during gas production from hydrate bearing sediments and the moment magnitude of induced seismicity. For the numerical analysis, sequential coupling of TOUGH+Hydrate and FLAC3D was used and the change in effective stress and consequent geomechanical deformation including fault reactivation was simulated by assuming that Mohr-Coulomb shear resistance criterion is valid. From the test production simulation of 30 days, we showed that pore pressure reduction as well as effective stress change hardly induces the fault reactivation in the vicinity of a production well. We also investigated the influence of stress state conditions to a fault reactivation, and showed that normal fault stress regime, where vertical stress is relatively greater than horizontal, may have the largest potential for the reactivation. We tested one simulation that earthquake can be induced during gas production and calculated the moment magnitude of the seismicity. Our calculation presented that all the magnitudes from the calculation were negative values, which indicates that induced earthquakes can be grouped into micro-seismic and as small as hardly perceived by human beings. However, it should be noted that the current simulation was carried out using the highly simplified geometric model and assumptions such that the further simulations for a scheduled test production and commercial scale production considering complex geometric conditions may produce different results.

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

• Geomechanics and Engineering
• /
• v.13 no.1
• /
• pp.1-23
• /
• 2017

A coupled liquid-gas-solid three-phase model, linking two numerical codes (TOUGH2/EOS3 and $FLAC^{3D}$), was firstly established and validated by simulating an in-situ air flow test in Essen. Then the coupled model was employed to investigate responses of multiphase flow and soil skeleton deformation to compressed air or freshwater injection using the same simulation conditions in an aquifer of Tianjin, China. The simulation results show that with injecting pressurized fluids, the vertical effective stress in some area decreases owing to the pore pressure increasing, an expansion of soil skeleton appears, and land uplift occurs due to support actions from lower deformed soils. After fluids injection stops, soil deformation decreases overall due to injecting fluids dissipating. With the same applied pressure, changes in multiphase flow and geo-mechanical deformation caused by compressed air injection are relatively greater than those by freshwater injection. Furthermore, the expansion of soil skeleton induced by compressed air injection transfers upward and laterally continuously with time, while during and after freshwater injection, this expansion reaches rapidly a quasi-steady state. These differences induced by two fluids injection are mainly because air could spread upward and laterally easily for its lower density and phase state transition appears for compressed air injection.