It is crucial to understand the hydro-mechanical behavior of rock mass to assess the performance of natural barriers. As rock fractures serve as both mechanically weak planes and prominent pathways for hydraulic flow, they significantly influence the hydro-mechanical behavior of the rock mass. Hence, understanding the characteristics of rock fractures is necessary to analyze the long-term behavior of natural barriers. In particular, fracture apertures are crucial parameters directly associated with groundwater flow and consequently hold significant importance in determining the hydro-mechanical behavior of natural barriers. Fracture apertures are defined as mechanical and hydraulic apertures, and various studies have been conducted to measure and analyze them. However, direct measurement of mechanical aperture according to changes in normal stress is known to be a challenging task. For this reason, there has been a scarcity of direct comparative findings between mechanical and hydraulic apertures under various normal stress conditions. This study aims to analyze the characteristics of the mechanical and hydraulic apertures according to changes in normal stress based on experimental results. A digital analysis technique using a pressure film image was applied to analyze the mechanical aperture characteristics of the fracture. This technique can be applied by performing a pressure film compression test and a normal stiffness test on a fracture specimen, and has the advantage of being able to derive mechanical apertures under various normal stress conditions. The hydraulic aperture characteristics of the fracture were analyzed based on Cubic law after measuring the flow rate by performing a constant pressure injection test under triaxial compression conditions. By applying various confining pressures, it was possible to examine the hydraulic apertures according to changes in normal stress conditions. Through the experimental results, the relationship between the mechanical and hydraulic apertures of the fracture was summarized under various normal stress conditions. In addition, the experimental results were used to examine the applicability of various empirical equations for mechanical and hydraulic apertures proposed in previous studies. The characteristics of the fracture aperture resulting from this study are significant because they are required in the hydro-mechanical model of natural barriers. Future studies will entail further experiments, with the objective of establishing novel relationships based on the accumulation of experimental data.

A low- and intermediate-level radioactive waste repository contains different types of radionuclides and organic complexing agents. Their chemical interaction in the repository can result in the formation of radionuclide-ligand complexes, leading to their high transport behaviors in the engineered and natural rock barriers. Furthermore, the release of radionuclides from the repository can pose a significant risk to both human health and the environment. This study explores the impact of different experimental conditions on the transport behaviors of 99Tc, 137Cs, and 238U through three types of barrier samples: concrete, sedimentary rock, and granite. To assess the transport behavior of the samples, the geochemical characteristics were determined using X-ray diffraction (XRD), X-ray fluorescence (XRF), Fouriertransform infrared spectroscopy (FTIR), scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDS), and Brunauer-Emmett-Teller (BET) analysis. The adsorption distribution coefficient (Kd) was used as an indicator of transport behavior, and it was determined in batch systems under different conditions such as solution pH (ranging from 7 to 13), temperature (ranging from 10 to 40°C), and with the presence of organic complexing agents such as ethylenediaminetetraacetic acid (EDTA), nitrilotriacetic acid (NTA), and isosaccharinic acid (ISA). A support vector machine (SVM) was used to develop a prediction model for the Kd values. It was found that, regardless of the experimental parameters, 99Tc may migrate easily due to its anionic property. Conversely, 137Cs showed low transport behaviors under all tested conditions. The transport behaviors of 238U were impacted by the order of EDTA > NTA> ISA, particularly with the concrete sample. The SVM models can also be used to predict the Kd values of the radionuclides in the event of an accidental release from the repository.

The distribution characteristics of rock fractures determine the hydro-mechanical behavior of natural barriers. Rock fractures are defined by various parameters, which are analyzed as the probability distribution from observation results by surveying the exposed rock surface or borehole. The size is known to have the most uncertainty among the fracture parameters because it cannot be directly measured. Therefore, various estimation methods have been proposed for fracture size distribution using the fracture traces observable on the rock surface. However, most methods are based on a planar survey area, limiting their applicability to the underground research laboratory (URL) excavated in the form of tunnels. This study aims to review a method that can be applied to estimate the size distribution of fractures in deep rock masses at the URL site. The estimation method using the joint center volume (JCV) has recently been extended to be applicable regardless of the geometry of the survey area, which means that it can be applied to the URL site with complex structures. To apply the JCV-based estimation method to non-planar survey areas, JCV calculation using Monte Carlo simulation and estimation of fracture size distribution using the maximum likelihood method are required. In this study, we applied the JCV-based estimation method to a tunnel-shaped survey area to examine its applicability to the URL site. The error rates were analyzed when there were fracture sets with various orientations, size distributions, and maximum fracture sizes in the rock mass, and it was found to be less than 10% in all cases. This result indicates that the JCV-based estimation method can be used to estimate the fracture size distribution of the surrounding rock mass if accompanied by a reliable survey of fracture traces on the tunnel surface inside the URL site. Also, since there are no restrictions on the geometry of the survey area, we can continuously update the estimation results during the URL excavation process to increase reliability. The fracture size distribution is essential for constructing the discrete fracture network (DFN) model of the rock mass units at the URL site. In the future, the uncertainty for the fracture size in the DFN model is expected to be reduced by applying the JCV-based estimation method.

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.

Low- and intermediate-level radioactive wastes have been disposed of in the first-phase deep underground silo disposal at Gyeongju in South Korea. These radioactive wastes contain harmful radionuclides such as Uranium-238 (238U), which can pose long-term and deleterious effects on humans and the natural environment. Ethylenediaminetetraacetic acid and isosaccharinic acid, which can be formed via cellulosic waste degradation under high alkaline conditions might considerably enhance the transport behavior of 238U with the intrusion of rainwater and groundwater. In this study, the engineered barriers (concrete and grout) and natural barriers (sedimentary rock and granite) were used to investigate the 238U transport behavior in artificial cementitious porewater of State I (pH 13.3) and State II (pH 12.5) based on groundwater or rainwater. The surface properties and geochemical compositions of barrier samples were characterized using XRD, XRF, SEM-EDX, and BET. The transport behaviors of 238U in various solution conditions were observed by sorption distribution coefficient (Kd) at a range of initial chelating agents concentration (10-5-10-2 M). The sorption behavior of 238U was retarded more in the engineered rock barriers than in the natural rock barriers. The mobility enhancement of 238U was more significant in State I than in State II. In comparison with the absence of chelating agents, negligible changes in the Kd values of 238U were observed at less than initial chelating agent concentrations of 10-4 M. However, the Kd values of 238U were significantly reduced at initial chelating agent concentrations higher than 10-3 M. Therefore, these experimental findings show that the transport behavior of 238U into the geo- and bio-sphere could be accelerated by the presence of chelating agents and the type of cement degradation states.

Engineered barriers (concrete and grout) in Low- and Intermediate-Level Waste (L/ILW) disposal facilities tend to degrade by groundwater or rainfall water over a long period of time. During the degradation process, radionuclides stored in the disposal facility might be released into the pore water, which can pass through the natural rock barriers (granite and sedimentary rock) and may reach the near-field and far-field. In this transportation, radionuclide might be sorbed onto the engineered and natural rock barriers. In addition, the organic complexing agent such as ethylenediaminetetraacetic acid (EDTA) and α-isosaccharinic acid (ISA), is also present in pore water, which may affect the sorption and mobility of radionuclide. In this study, the sorption and mobility of 90Sr under different conditions such as two pHs (7 and 13), different initial concentrations of organic complexing agents (from 10-5 M to 10-2 M), and solutions (groundwater, pore water, and rainfall water) were investigated in a batch system. The groundwater was collected at the L/ILW disposal facility located at Gyeongju in South Korea. The pore water and rainfall water were artificially made in the laboratory. The concrete, grout, granite, and sedimentary rock samples were collected from the same study sites from where the groundwater was collected. The rock samples were crushed to 53-150 micrometers and were characterized by XRD, XRF, SEM-EDS, BET, and zeta potential analyzer. 90Sr concentration was determined using liquid scintillation counting. The sorption of 90Sr was described by distribution coefficients (Kd) and sorption reduction factor (SRF). In the case of EDTA, the Kd values of 90Sr remained constant from 10-5 M to 10-3 M and tended to decrease at 10-2 M, while in case of ISA the Kd values decreased steadily as the concentration of ISA was increased from 10-5 M to 10-3 M; However, a sudden reduction in the Kd values were observed above 10-2 M. In comparison to EDTA, ISA gave a higher SRF of 90Sr. Therefore, from the above results, it can be concluded that the presence of ISA has a greater effect on the sorption and mobility of radionuclide in the solutions than EDTA, and the radionuclide may reach near- and far-field of the L/ILW disposal facility.

Radionuclides stored in a radioactive waste repository over a long period of time might be leached through the barriers such as engineered rock (cement) and natural rock (granite). Organic complexing agents such as ethylenediaminetetraacetic acid (EDTA) and isosaccharinic acid (ISA) may also influence the mobility of radionuclides. In this study, a continuous fixed-column reactor packed with engineered and natural rocks was designed to investigate the effect of organic complexing agents on cesium mobility through cement and granite under anaerobic conditions. The influent flow rate of the mixed solution (organic complexing agent and cesium) at the column bottom was 0.1 mL/min, while that of groundwater was 0.2 mL/min, which was introduced between cement and granite layers in the middle of the column. The hydraulic properties such as diffusion coefficient and retardation factor were derived by a bromide tracer test. The effects of different operating parameters, such as initial cesium concentrations, initial EDTA or ISA concentrations, and bed size, on the cesium adsorption were investigated. The Thomas, Adams-Bohart, and Yoon-Nelson models were applied to the experimental data to predict the breakthrough curves using non-linear regression. These results suggest that organic complexing agents such as EDTA and ISA significantly influence the mobility of cesium in the barriers, indicating that the presence of complexing agents enhances the migration of cesium to the geosphere.

Radionuclides can be leached into groundwater or soil over a long period of time due to unexpected situations even after being permanently disposed of in a repository. Therefore, it is necessary to investigate the mobility of radionuclides for the safety assessment of radioactive waste disposal. In this study, the effects of organic complexing agents such as ethylenediaminetetraacetic acid (EDTA) and isosaccharinic acid (ISA) on the sorption behavior of 239Pu and 99Tc over cementitious (concrete and grout) and natural rock samples (granite and sedimentary rock) were investigated in batch sorption experiments. For characterization of rock samples, XRD, XRF, FT-IR, FE-SEM, BET, and Zeta-potential analyses were performed. For the evaluation of mobility, the distribution coefficient (Kd) was selected and compared. The adsorption experiment was carried out at two pHs (7 and 13), a temperature of 20°C, and a range of organic complexing agents concentrations (10-7~10-2 M and 10- 5~10-2 M for 239Pu and 99Tc, respectively). The radionuclides concentrations in adsorption samples were analyzed using ICP-MS. The Kd values for 239Pu in all rock samples reduced significantly due to the presence of EDTA, even at low concentrations such as 10-5 M. In the case of ISA, the limiting noeffect concentration was much higher than that of EDTA. On the other hand, 99Tc showed relatively lower Kd values than 239Pu, and the sorption behavior of 99Tc was almost unaffected by the organic complexing agents for all rock samples. Therefore, it is possible to assume that the increased mobility of radionuclides, especially, 239Pu, in groundwater caused by the lowering of sorption at even low concentrations of organic complexing agents may result in the transport of radionuclides to the nearand far-field location of the repository.