• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.609-619
• /
• 2018

This study was carried out at a local limestone mine to analyze the ventilation efficiency of the shaft equipped with a main fan. The results show that its ventilation efficiency is clearly verified for the natural as well as the mechanical ventilation. The airflow rate of $11.7m^3/s$ was induced by the natural ventilation force and the maximum quantity is almost same as the airflow rate estimated by monitoring the average temperatures in the upcast and downcast air columns. Meanwhile, the airflow rate exhausted by the main fan through the shaft was $20.3{\sim}24.8m^3/s$; variation of the quantity was caused by the upward shift of the mine ventilation characteristic curve due to the frequent movement of the equipment. This indicates efforts are required to reduce the ventilation resistance and raise the quantity supplied by the main fan. The turbulent diffusion coefficients along the 1912 m long airway from the portal to the shaft bottom was estimated to be $15m^2/s$ and $18m^2/s$. Since these higher coefficients imply that contaminants will be dispersed at a faster velocity than the airflow, prompt exhaust method should be planned for the effective air quality control. The ventilation shaft and main fan are definitely what local limestone mines inevitably need for better working environment and sustainable development.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.635-650
• /
• 2018

Many people have been recognized that the Korean Peninsula is no longer safe area from the earthquake by the recent earthquakes occurred in the country. The earthquakes that occurred at Pohang and Gyeongju appeared differently from them considered in the seismic design and researches on the seismic design method have been also conducted by many researchers. Studies on seismic loads are mainly focused on existing superstructures, and research involving them has been actively carried out in reality. However, paper regarding structural stability of reinforcement from seismic load such as soil-nails, rock-bolts, ground anchors which were constructed to ensure stability of serviced structure have been published rarely. In this study, ground anchor been effected by static load and seismic load which is settled in the weathered rock is analyzed. Results for static load are obtained from field test and seismic load is from numerical analysis. In this study, the behavioral characteristics of the ground anchor were analyzed by numerical analysis in case of seismic loading based on the result of the in-situ tensile test of the ground anchor settled weathered rock. As a result, settlement of concrete block due to application of tension force for ground anchor occurred as well as following loss of axial force for ground anchor. Also, as bond length and period of seismic load are longer, increasement of displacement is greater.

• Tunnel and Underground Space
• /
• v.29 no.4
• /
• pp.262-279
• /
• 2019

In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the environment can be exposed to radioactivity, it is very important to clarify gas migration in terms of long-term integrity of the engineered barrier system. In particular, it is necessary to identify the hydro-mechanical mechanism for the dilation flow, which is a very important gas flow phenomenon only in medium containing large amounts of clay materials such as bentonite buffer, and to develop and validate new numerical approach for the quantitative evaluation of the gas migration phenomenon. Therefore, in this study, we developed a two-phase flow model considering the mechanical damage model in order to simulate the gas migration in the engineered barrier system, and validated with 1D gas flow modelling through saturated bentonite under constant volume boundary conditions. As a result of numerical analysis, the rapid increase in pore water pressure, stress, and gas outflow could be simulated when the dilation flow was occurred.

• Tunnel and Underground Space
• /
• v.29 no.4
• /
• pp.230-242
• /
• 2019

While various automatic rock fracture survey methods have been researched, the evaluation of the accuracy of these methods raises issues due to the absence of a metric which fully expresses the similarity between automatic and manual fracture maps. Therefore, this paper proposes a geometry similarity metric which is especially designed to determine the overall similarity of fracture maps and to evaluate the accuracy of rock fracture survey methods by a single number. The proposed metric, Scanline Intersection Similarity (SIS), is derived by conducting a large number of scanline surveys upon two fracture maps using Python code. By comparing the frequency of intersections over a large number of scanlines, SIS is able to express the overall similarity between two fracture maps. The proposed metric was compared with Intersection Over Union (IoU) which is a widely used evaluation metric in computer vision. Results showed that IoU is inappropriate for evaluating the geometry similarity of fracture maps because it is overly sensitive to minor geometry differences of thin elongated objects. The proposed metric, on the other hand, reflected macro-geometry differences rather than micro-geometry differences, showing good agreement with human perception. The metric was further applied to evaluate the accuracy of a deep learning-based automatic fracture surveying method which resulted as 0.674 (SIS). However, the proposed metric is currently limited to 2D fracture maps and requires comparison with rock joint parameters such as RQD.

• 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