• Journal of the Korean Geotechnical Society
• /
• v.32 no.1
• /
• pp.19-33
• /
• 2016

Stability analysis based on the limit equilibrium method combined with the result of infiltration analysis is commonly used to evaluate the effect of rainfall infiltration on the slope stability. Soil is a three-phase mixture composed of solid particle, water and air. Therefore, a fully coupled mixture theories of stress-deformation behavior and the flow of water and air should be used to accurately analyze the process of rainfall infiltration through soil slope. The purpose of this study is to study the effect of interaction of air and water flow on the mechanical stability of slope. In this study, stability analyses based on the coupled hydro-mechanical model of three-phases were conducted for slope of weathered granite soil widespread in Korea. During the process of hydro-mechanical analysis strength reduction technique was applied to evaluate the effect of rainfall infiltration on the slope stability. The results showed an increase of air pressure during infiltration because rain water continuously displaced the air in the unsaturated zone. Such water-air interaction in the pore space of soil affects the stress-deformation behavior of slope. Therefore, the results from the three-phase model showed different behavior from the solid-water model that ignores the transport effect of air in the pores.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.20 no.4
• /
• pp.429-453
• /
• 2022

Numerical modeling and scenario composition are needed to characterize the geological environment of the disposal site and analyze the long-term evolution of natural barriers. In this study, processes and features of the hydro-mechanical behavior of natural barriers were categorized and represented using the interrelation matrix proposed by SKB and Posiva. A hydro-mechanical coupled model was evaluated for analyzing stress field changes and fracture zone re-activation. The processes corresponding to long-term evolution and the hydro-mechanical mechanisms that may accompany critical processes were identified. Consequently, practical numerical methods could be considered for these geological engineering issues. A case study using a numerical method for the stability analysis of an underground disposal system was performed. Critical stress distribution regime problems were analyzed numerically by considering the strata's movement. Another case focused on the equivalent continuum domain composition under the upscaling process in fractured rocks. Numerical methods and case studies were reviewed, confirming that an appropriate and optimized modeling technique is essential for studying the stress state and geological history of the Korean Peninsula. Considering the environments of potential disposal sites in Korea, selecting the optimal application method that effectively simulates fractured rocks should be prioritized.

• Structural Engineering and Mechanics
• /
• v.71 no.6
• /
• pp.699-712
• /
• 2019

The reinforced concrete lining of hydraulic pressure tunnels tends to crack under high inner water pressure (IWP), which results in the inner water exosmosis along cracks and involves typical hydro-mechanical interaction. This study aims at the development, validation and application of an indirect-coupled method to simulate the lining cracking process. Based on the concrete damage plasticity (CDP) model, the utility routine GETVRM and the user subroutine USDFLD in the finite element code ABAQUS is employed to calculate and adjust the secondary hydraulic conductivity according to the material damage and the plastic volume strain. The friction-contact method (FCM) is introduced to track the lining-rock interface behavior. Compared with the traditional node-shared method (NSM) model, the FCM model is more feasible to simulate the lining cracking process. The number of cracks and the reinforcement stress can be significantly reduced, which matches well with the observed results in engineering practices. Moreover, the damage evolution of reinforced concrete lining can be effectively slowed down. This numerical method provides an insight into the cracking process of reinforced concrete lining in hydraulic pressure tunnels.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

• Structural Engineering and Mechanics
• /
• v.86 no.1
• /
• pp.139-156
• /
• 2023

The reinforced concrete lining in the hydraulic pressure tunnel tends to crack during the water-filling process. The lining will be detached from the surrounding rock due to the inner water exosmosis along concrete cracks. From the previous research achievements, the cohesive element is widely adopted to simulate the concrete crack but rarely adopted to simulate the lining-rock interface. In this study, the zero-thickness cohesive element with hydro-mechanical coupling property is not only employed to simulate the traditional concrete crack, but also innovatively introduced to simulate the lining-rock interface. Combined with the indirect-coupled method, the hydro-mechanical coupling algorithm of the reinforced concrete lining in hydraulic pressure tunnels is proposed and implemented in the finite element code ABAQUS. The calculated results reveal the cracking mechanism of the reinforced concrete lining, and match well with the observed engineering phenomenon.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.21 no.1
• /
• pp.23-41
• /
• 2023

Coupled thermo-hydraulic-mechanical (THM) processes are essential for the long-term performance of deep geological disposal of high-level radioactive waste. In this study, a numerical sensitivity analysis was performed to analyze the effect of rock properties on THM responses after the execution of the heater test at the Kamaishi mine in Japan. The TOUGH-FLAC simulator was applied for the numerical simulation assuming a continuum model for coupled THM analysis. The rock properties included in the sensitivity study were the Young's modulus, permeability, thermal conductivity, and thermal expansion coefficients of crystalline rock, rock salt, and clay. The responses, i.e., temperature, water content, displacement, and stress, were measured at monitoring points in the buffer and near-field rock mass during the simulations. The thermal conductivity had an overarching impact on THM responses. The influence of Young's modulus was evident in the mechanical behavior, whereas that of permeability was noticed through the change in the temperature and water content. The difference in the THM responses of the three rock type models implies the importance of the appropriate characterization of rock mass properties with regard to the performance assessment of the deep geological disposal of high-level radioactive waste.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.670-691
• /
• 2018

This study presents the research results of the BMT(Benchmark Model Test) simulations of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to predict fault reactivation and the coupled hydro-mechanical behavior of fault. BMT scenario simulations of Task B were conducted to improve each numerical model of participating group by demonstrating the feasibility of reproducing the fault behavior induced by water injection. The BMT simulations consist of seven different conditions depending on injection pressure, fault properties and the hydro-mechanical coupling relations. TOUGH-FLAC simulator was used to reproduce the coupled hydro-mechanical process of fault slip. A coupling module to update the changes in hydrological properties and geometric features of the numerical mesh in the present study. We made modifications to the numerical model developed in Task B Step 1 to consider the changes in compressibility, Permeability and geometric features with hydraulic aperture of fault due to mechanical deformation. The effects of the storativity and transmissivity of the fault on the hydro-mechanical behavior such as the pressure distribution, injection rate, displacement and stress of the fault were examined, and the results of the previous step 1 simulation were updated using the modified numerical model. The simulation results indicate 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 interaction and collaboration with other research teams of DECOVALEX-2019 Task B and validated using the field experiment data in a further study.

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

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.4 no.2
• /
• pp.117-131
• /
• 2006

A coupled T-H-M(Thermo-Hydro-Mechanical) analysis was carried out for KENTEX (KAERI Engineering-scale T-H-M Experiment for Engineered Barrier System), which is a facility for validating the coupled T-H-M behavior in the engineered barrier system of the Korean reference HLW(high-level waste) disposal system. The changes of temperature, water saturation, and stress were estimated based on the coupled T-H-M analysis, and the influence of the types of mechanical constitutive material laws was investigated by using elastic model, poroelastic model, and poroelastic-plastic model. The analysis was done using ABAQUS, which is a commercial finite element code for general purposes. From the analysis, it was observed that the temperature in the bentonite increased sharply for a couple of days after heating the heater and then slowly increased to a constant value. The temperatures at all locations were nearly at a steady state after about 37.5 days. In the steady state, the temperature was maintained at $90^{\circ}C$ at the interface between the heater and the bentonite and at about $70^{\circ}C$ at the interface between the bentonite and the confining cylinder. The variation of the water saturation with time in bentonite was almost same independent of the material laws used in the coupled T-H-M processes. By comparing the saturation change of T-H-M and that of H-M(Hydro-Mechanical) processes using elastic and poroelastic material mod31 respectively, it was found that the degree of saturation near the heater from T-H-M calculation was higher than that from the coupled H-M calculation mainly because of the thermal flux, which seemed to speed up the saturation. The stresses in three cases with different material laws were increased with time. By comparing the stress change in H-M calculation using poroelasetic and poroelasetic-plastic model, it was possible to conclude that the influence of saturation on the stress change is higher than the influence of temperature. It is, therefore, recommended to use a material law, which can model the elastic-plastic behavior of buffer, since the coupled T-H-M processes in buffer is affected by the variation of void ratio, thermal expansion, as well as swelling pressure.

• Geomechanics and Engineering
• /
• v.15 no.3
• /
• pp.823-831
• /
• 2018

The most important issue during shield TBM tunneling in soft ground formations is to appropriately control ground surface settlement. Among various operational conditions in shield TBM tunneling, the face pressure and backfill pressure should be the most important and immediate measure to restrain surface settlement during excavation. In this paper, a 3-D hydro-mechanical coupled FE model is developed to numerically simulate the entire process of shield TBM tunneling, which is verified by comparing with real field measurements of ground surface settlement. The effect of permeability and stiffness of ground formations on tunneling-induced surface settlement was discussed in the parametric study. An increase in the face pressure and backfill pressure does not always lead to a decrease in surface settlement, but there are the critical face pressure and backfill pressure. In addition, considering the relatively low permeability of ground formations, the surface settlement consists of two parts, i.e., immediate settlement and consolidation settlement, which shows a distinct settlement behavior to each other.