In KAERI, a site descriptive model for stress field estimation had already been constructed by using integrated field data within KURT site scale. A sub-divided rock block domain containing major fracture zones has spatial rock mass and fault properties. The properties were decided based on the rock classification results of several borehole investigations. Modeling for maximum and minimum horizontal stress field estimation was performed and compared with the in-situ data. As a result, a depth-dependent stress ratio was adopted to obtain numerical results closer to actual in-situ data. Although the results were suitable at a relatively low depth (~500 m), there is still some deviation trend at a deep depth. This study aims to improve these modeling results by incorporating not only depth-dependent stress ratio but also changes in rock mass properties along the depth. The deep borehole of DB2 in the KURT site indicated fracture distribution corresponding to the property changes. Natural fractures are typically randomly oriented, and the fracture frequency decreases with increasing depth. The increase in P-wave velocity log data accompanies these features. A discrete fracture network (DFN) model can be used to simulate fractured rock explicitly, but DFN modeling is not feasible for site scale analysis because of its numerical efficiency. Therefore, as a preliminary model in this study, the effect of fracture distribution was considered by substituting the influence for the depth-dependent property. The properties were estimated from the fracture frequency and P-wave velocity log data. The influence of elastic modulus and density on the site stress field was dominant, with decreasing the deviation trend between modeling and in-situ data at a deep depth. Considering that the depth of the repository construction is within about 500 m, it may not be necessary to consider the change of rock properties with depth. However, it was determined that the rock property effect might need to be considered when the loading conditions change due to subsidence in the long-term evolution scenario. Continuously, this site descriptive modeling will be interdependently conducted with a representative DFN block model for deriving equivalent properties in fractured rock.

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