• Tunnel and Underground Space
• /
• v.31 no.5
• /
• pp.333-359
• /
• 2021

Gases such as hydrogen and radon can be generated around the canister in high-level radioactive waste disposal systems due to several reasons including the corrosion of metal materials. When the gas generation rate exceeds the gas diffusion rate in the low-permeability bentonite buffer, the gas phase will form and accumulate in the engineered barrier system. If the gas pressure exceeds the gas entry pressure, gas can migrate into the bentonite buffer, resulting in pathway dilation flow and advective flow. Because a sudden occurrence of dilation flow can cause radionuclide leakage out of the engineered barrier of the radioactive waste disposal system, it is necessary to understand the gas migration behavior in the bentonite buffer to quantitatively evaluate the long-term safety of the engineered barrier. Experimental research investigating the characteristics of gas migration in saturated bentonite and research developing numerical models capable of simulating such behaviors are being actively conducted worldwide. In this technical note, previous gas injection experiments and the numerical models proposed to verify such behaviors are introduced, and the future challenges necessary for the investigation of gas migration are summarized.

• Tunnel and Underground Space
• /
• v.31 no.5
• /
• pp.309-332
• /
• 2021

Natural barrier systems surrounding the geological repository for the high-level radioactive waste should guarantee the hydraulic performance for preventing or delaying the leakage of radionuclide. In the case of the behavior of a crystalline rock, the hydraulic performance tends to be decided by the existence of discontinuities, so the coupled hydro-mechanical(HM) processes on the discontinuities should be characterized. The discontinuum modelling can describe the complicated behavior of discontinuities including creation, propagation, deformation and slip, so it is appropriate to model the behavior of a crystalline rock. This paper investigated the coupled HM processes in discontinuum modelling such as UDEC, 3DEC, PFC, DDA, FRACOD and TOUGH-UDEC. Block-based discontinuum methods tend to describe the HM processes based on the fluid flow through the discontinuities, and some methods are combined with another numerical tool specialized in hydraulic analysis. Particle-based discontinuum modelling describes the overall HM processes based on the fluid flow among the particles. The discontinuum methods that are currently available have limitations: exclusive simulations for two-dimension, low hydraulic simulation efficiency, fracture-dominated fluid flow and simplified hydraulic analysis, so it could be improper to the modelling the geological repository. Based on the concepts of various discontinuum modelling compiled in this paper, the advanced numerical tools for describing the accurate coupled HM processes of the deep geological repository should be developed.

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

• Tunnel and Underground Space
• /
• v.31 no.3
• /
• pp.167-183
• /
• 2021

It is essential to understand the complex thermo-hydro-mechanical-chemical (THMC) coupled behavior in the engineered barrier system and natural barrier system to secure the high-level radioactive waste repository's long-term safety. The heat from the high-level radioactive waste induces thermal pressurization and vaporization of groundwater in the repository system. Groundwater inflow affects the saturation variation in the engineered barrier system, and the saturation change influences the heat transfer and multi-phase flow characteristics in the buffer. Due to the complexity of the coupled behavior, a numerical simulation is a valuable tool to predict and evaluate the THMC interaction effect on the disposal system and safety assessment. To enhance the knowledge of THMC coupled interaction and validate modeling techniques in geological systems. DECOVALEX, an international cooperation project, was initiated in 1992, and KAERI has participated in the projects since 2008 in Korea. In this study, we introduced the main contents of all tasks in the DECOVALEX-2023, the current DECOVALEX phase, to the rock mechanics and geotechnical researchers in Korea.

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

• Tunnel and Underground Space
• /
• v.31 no.1
• /
• pp.52-65
• /
• 2021

In this study, authors analyzed the vapor cloud explosion induced by propane leak at the PEMIX Terminal, which is the propane storage facility outside of Mexico City. TNT equivalence mass for the leaked 4750 kg propane was estimated to be 9398 kg. Blast parameters such as peak overpressure, positive phase duration, and impact at 40-400 (m) away from the center of the explosion were calculated by applying TNT Equivalency Method and Multi-Energy Method. The probability of damage due to lung damage, eardrum rupture, head impact, and whole-body displacement impact by applying the probit function obtained using blast parameters was evaluated. The peak overpressure obtained using Multi-Energy Method was found to be greater than the peak overpressure obtained by applying the TNT Equivalency Method at all distances considered, but it was evaluated that there was no significant difference from the points above 200 m. The peak overpressure obtained by Multi-Energy Method was computed to assess the extent of damage to the structure, and it was shown that structures within 100 m of the explosion center would collapse completely, and that the glasses of the structures 400 m away would be almost broken. The probability of death due to lung damage was shown to vary depending on a human body's position located in the propagating direction of shock wave, and if there is a reflecting surface in the immediate surroundings of a human body, the probability of death was estimated to be the greatest. The impact of shock wave on lung damage, eardrum rupture, head impact, and whole-body displacement impact was evaluated and found to affect whole-body impact < lung damage < eardrum rupture

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.509-518
• /
• 2020

This study constructed a new simulation platform, including mesh generation process, numerical simulation, and post-processing for results analysis based on exploration data to perform real-scale numerical analysis considering the actual geological structure efficiently. To build the simulation platform, we applied for open-source programs. The source code is open to be available for code modification according to the researcher's needs and compatibility with various numerical simulation programs. First, a three-dimensional model(3D) is acquired based on the exploration data obtained using a drone. Then, the domain's mesh density was adjusted to an interpretable level using Blender, the free and open-source 3D creation suite. The next step is to create a 3D numerical model by creating a tetrahedral volume mesh inside the domain using Gmsh, a finite element mesh generation program. To use the mesh information obtained through Gmsh in a numerical simulation program, a converting process to conform to the program's mesh creation protocol is required. We applied a Python code for the procedure. After we completed the stability analysis, we have created various visualization of the study using ParaView, another open-source visualization and data analysis program. We successfully performed a preliminary stability analysis on the full-scale Dokdo model based on drone-acquired data to confirm the usefulness of the proposed platform. The proposed simulation platform in this study can be of various analysis processes in future research.

• Tunnel and Underground Space
• /
• v.31 no.1
• /
• pp.41-51
• /
• 2021

The hydraulic fracturing developed to improve permeability of tight reservoir is one of key stimulation technologies for developing unconventional resources such as shale gas and deep geothermal energy. The experimental study was conducted to improve disadvantage of hydraulic fracturing which has simple fracture pattern and poor fracturing efficiency. The fracturing experiments was conducted for tight rock using various fracturing fluids, water, N2, and CO2 and the created fracture pattern and fracturing efficiency was analyzed depending on fracturing fluids. The borehole pressure increased rapidly and then made fractures for hydraulic fracturing with constant injection rate, however, gas fracturing shows slowly increased pressure and less fracture pressure. The 3D tomography technic was used to generate images of induced fracture using hydraulic and gas fracturing. The stimulated reservoir volume (SRV) was estimated increment of 5.71% (water), 12.72% (N2), and 43.82% (CO2) respectively compared to initial pore volume. In addition, permeability measurement was carried out before and after fracturing experiments and the enhanced permeability by gas fracturing showed higher than hydraulic fracturing. The fracture conductivity was measured by increasing confining stress to consider newly creating fracture and closing induced fracture right after fracturing. When the confining stress was increased from 2MPa to 10MPa, the initial permeability was decreased by 89% (N2) and 50% (CO2) respectively. This study shows that the gas fracturing makes more permeability enhancement and less reduction of induced fracture conductivity than hydraulic fracturing.

• Journal of the Korean Geosynthetics Society
• /
• v.20 no.1
• /
• pp.1-8
• /
• 2021

In urban areas, structures are installed deep underground in the lower part of the structure to utilize space. Therefore, a retaining wall is used to prevent earth pressure from the ground when constructing a structure. Due to the development of construction technology, retaining wall applied to excavation work are used to prevent danger such as falling rocks and landslides in temporary facilities when construction or retaining walls are installed. In general, the application of a retaining wall to a temporary facility during the embankment construction is the case of expanding an existing roads or railways. Therefore, it is necessary to study the retaining wall applied to the embankment construction such as the double-track site of the high-speed railway. In this study, two types of common one row H-pile retaining wall and two types of IER retaining wall were analyzed, and the stability of the retaining wall applied to the construction of double-track of the high-speed railway was analyzed. The earth retaining wall is a construction method that combines forced pile applied to the stabilization of the slope with the wall of the earth retaining wall. As a result of the analysis, the IER retaining wall had maximum lateral displacement of 19.0% compared to the type with H-plie installed only in the front while dynamic load was applied. In addition, the slower the speed of high-speed railway, the more displacement occurred, and the results show that more caution is needed when designing the ground in low-speed sections.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.