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.

In recent years, the importance of the thermo-hydraulic-mechanical-chemical coupled processes is increasing in the performance assessment (PA) of the high-level radioactive waste repository. In the case of mechanical behavior, it is very important because it can affect fluid flow and radionuclide transport by changing the porosity and permeability of the medium. In particular, Excavation Damaged Zone (EDZ) should be considered essential in PA because the migration of radionuclide is affected by the enhanced hydraulic transmissivity and altered geomechanical behavior of EDZ. Furthermore, due to various thermo-hydraulic behaviors such as decay heat generated from radioactive waste, pore water pressure increase, and swelling pressure of bentonite buffer material, mechanical evolution is occurred which may change the size and physical properties of EDZ. Therefore, to solve this problem, analysis of coupled thermal-hydraulic-mechanical (THM) processes with the effect of long-term evolution of EDZ due to the mechanical behavior should be accompanied. In this study, numerical model for the long-term evolution due to mechanical behavior considering EDZ using the Adaptive Process-based total system performance analysis framework for a geological disposal system (APro) proposed by the Korea Atomic Energy Research Institute (KAERI). In the case of EDZ, the concept of Mazars’ damage evolution model was applied to simulate the behavior using the continuum model, and the change in hydraulic properties according to the degree of damage was considered. To investigate the importance of mechanical behavior in PA, the results were compared by performing numerical analysis according to the presence or absence of mechanical analysis. Finally, numerical analysis considering the mechanical evolution of EDZ was conducted using the model developed in this study to investigate the effect of EDZ.