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.

고준위방사성폐기물의 처분터널 및 처분공 간격을 결정하고 처분시스템의 성능을 평가하기 위해서는 열-수리-역학적인 복합 거동 변화에 대한 이해가 반드시 필요하고 이를 반영하여 해석해야만 한다. 하지만 한국형 기준 처분시스템에서의 처분 터널 및 처분공 간격을 결정하기 위해 수행된 기존의 연구들은 이러한 복합거동 특성을 반영하지 않고 열 해석 결과만을 근거로 처분시스템을 설계하였다. 따라서 본 연구에서는 열-수리-역학적인 복합거동 특성을 반영하여 한국형 기준 처분시스템의 성능을 TOUGH2-MP/FLAC3D를 이용하여 평가하였다. 고준위방사성폐기물이 처분된 이후 방사성 붕괴열에 의해 처분 시스템의 온도는 급격히 증가하다가 붕괴열의 감소로 온도는 서서히 감소하였으며, 해석 기간 1,000년 동안 벤토나이트 완충재의 최고 온도는 설계 기준인 100℃ 이하로 유지되는 것으로 나타났다. 처분용기와 벤토나이트 완충재의 계면에서의 최고 온도는 약 3.21년이 지난 시점에 용기의 중간 지점에서 약 96.2℃로 나타났으며, 암반에서의 최고 온도는 폐쇄 후 약 17년 이 지난 시점에서 약 68.2℃로 계산되었다. 처분용기 부근 벤토나이트 완충재는 처분 초기에 온도 변화에 따른 건조현상이 발생하여 포화도가 감소하지만, 시간이 지남에 따라 주변 암반으로부터의 지하수 유입에 의해 포화도가 증가하는 것으로 계 산되었다. 이후, 벤토나이트 완충재 및 뒷채움재 모두 약 266년 이후 완전히 포화되는 것으로 계산되었다. 처분시스템에서의 온도 변화에 따른 열응력 그리고 벤토나이트 완충재 및 뒷채움재의 팽윤압으로 인한 응력 변화가 처분장 주변 암반에 미치는 영향을 평가하고자 수치해석에서 계산된 응력을 스폴링 강도(spalling strength)와 Mohr-coulomb 파괴 기준식과 비교 하였다. 계산 결과 일축압축강도와 스폴링 강도에 도달하지 않는 것으로 나타나 처분시스템이 스폴링에 의한 파괴는 나타나지 않을 것으로 판단되며, Mohr-coulomb 파괴 기준 역시 충족하는 것으로 나타났다. 본 연구에서 사용된 수치해석 코드와 방법론은 다양한 조건에서의 한국형 기준 처분시스템에 대한 성능평가뿐만 아니라, 복층 처분시스템에 대한 설계와 성능평가에 활용될 수 있을 것으로 판단된다.

고준위방사성폐기물 처분시스템에서는 방사성 핵종의 붕괴열과 암반으로부터의 지하수 유입으로 열응력 및 팽윤압의 발생으로 열-수리-역학적 복합거동(coupled thermo-hydro-mechanical behavior)이 예상되기 때문에 한국원자력연구원은 처분시스템 및 근계암반에서의 열-수리-역학적인 복합거동 특성을 평가하기 위해서 지하처분연구시설(KAERI Underground Research Tunnel, KURT)에서 2016년부터 현장시험(In-situ Demonstration of Engineered Barrier System, In-DEBS)을 수행 중에 있다. 본 연구에서는 In-DEBS 현장시험 데이터 분석하고 벤토나이트 완충재와 화강암반에서의 열-수리-역학적 복합거동 특성을 평가하기 위해 TOUGH2-MP/FLAC3D을 이용하여 수치해석을 수행하였다. 또한 벤토나이트 블록과 KURT 화강암의 열-수리-역학적 복합거동 특성을 평가하기 위해 사용된 각각의 열, 수리, 그리고 역학적 모델의 적합성을 평가하고 자 현장시험에서 계측된 온도, 상대습도, 그리고 변위의 결과와 수치해석으로 계산된 결과를 비교하였다. 온도와 상대습도의 계산 결과를 현장 데이터와 비교·분석한 결과, 전체적으로 유사한 경향을 보일 뿐만 아니라 시간에 따라 변화하는 정량적인 값 역시 유사하게 나타났다. 역학적 해석 결과를 살펴보면, 계산된 변위의 전반적인 경향은 유사하지만 해석 결과가 계측 값에 비해 상대적으로 작게 나타났다. 축대칭 모델을 이용하여 In-DEBS 현장시험에서 관측된 열-수리-역학적 복합거동 특성을 전반적으로 평가할 수 있었지만, 벤토나이트 블록 및 KURT 암반에서의 열-수리-역학적 복합거동을 면밀히 살펴보기 위해서는 추후 터널의 형상과 주변 KURT 터널의 영향을 반영한 3차원 해석이 필요할 것으로 판단된다. 본 연구에서 사용된 입력 물성과 열-수리-역학적 모델은 추후 In-DEBS 장기 거동 및 처분시스템에서의 열-수리-역학적 복합거동 특성을 평가하고 예측하는데 활용될 수 있을 것으로 기대된다.