Discrete element modelling of fracture transmissivity evolution around a deep geological repository for high-level radioactive waste

Abstract

원자력발전에서 발생하는 고준위 방사성 폐기물은 계속해서 방사선과 열을 발생하기 때문에, 충분한 시간 동안 생태계로부터 격리되어야 한다. 고준위 방사성 폐기물을 영구적으로 처분하는 방법으로 공학적방벽과 천연방벽을 통해 격리시키는 심지층처분방식이 제시되었다. 심지층처분 시스템의 성능을 평가하기 위하여 공학적방벽과 천연방벽의 장기적 거동을 조사하는 다양한 실험 및 수치해석이 진행되고 있다. 특히 천연방벽의 투수량계수는 핵종의 누출을 평가할 수 있는 수단으로 다양한 열, 수리. 역학적 영향 요인에 대한 분석이 필요하다. 본 논문에서는 고준위 방사성 폐기물 심층처분장의 건설 및 운영 중에 천연방벽, 특히 결정질 균열암반의 투수량계수 변화에 대하여 검토하였다. 결정질 균열암반의 실제 거동을 모사하기 위하여 스웨덴 애스푀 지하연구시설에서 취득한 암석 및 균열 물성을 활용하여 삼차원 개별요소법 기반 균열암반모델을 구축하였다. 균열투수량계수에 영향을 주는 요인인 굴착손상영역, 열, 빙하, 지진에 의한 응력 변화를 수치해석 모델에 반영하였다. 굴착손상영역의 영향은 터널 자유면에 의한 응력재분배가 야기한 균열의 수직변형 및 전단 미끄러짐을 통하여 투수량계수의 변화를 분석하였으며, 이렇게 구축된 터널 모델에 각각 열원, 빙하에 의한 경계조건변화, 지진에 의한 동적 하중을 추가하여 각 영향요인에 의한 투수량계수 변화를 분석하였다. 굴착손상영역으로 인한 영향은 크게는 1,000배가량의 투수량계수 증진으로 나타났고, 수직 열림/닫힘과 전단 미끄러짐에 의한 영향이 균열의 방향에 따라서 복합적으로 작용한 것으로 나타났다. 열, 빙하, 지진에 의한 영향 분석에서는 모두 하중의 재하와 제하가 반복되는 비슷한 응력변화 양상을 수반하였으며 영향의 분석이 끝난 후에는 모두 초기 응력조건으로 회귀하는 과정이 포함되어 있었다. 열, 빙하, 지진에 의한 영향 분석 모두 일부 균열에서 비가역적인 투수량계수의 증진이 확인되었다. 시간에 따른 투수량계수 및 응력 상태를 분석한 결과 하중의 재하 시에는 균열의 수직 닫힘이 지배적으로 나타났지만, 제하 시에는 닫혔던 균열이 가역적으로 열림과 동시에 수직 닫힘에 의해 가려졌던 비가역적 전단팽창의 영향이 드러났다. 이러한 가역적 수직변형과 비가역적 전단팽창의 영향은 균열의 기하학적 특성에 따라 선별적으로 나타났다. 투수량계수의 변화를 터널 주변 절리의 기하학적 특성에 대하여 보다 일반적으로 분석하기 위하여, 균일절리군을 터널 주변에 추가하여 각각 굴착손상영역, 열, 빙하, 지진에 의한 영향을 정리하였다. 투수량계수의 증진은 주로 최대주응력과 평행인 절리에서 크게 나타났으며, 이는 최대주응력의 영향이 작아 수직방향으로 작용하는 압축응력이 상대적으로 낮았기 때문이었다. 또한 터녈면으로부터 약간 기울어진 방향의 절리에서 투수량계수의 증진이 크게 발생했으며, 이는 역학적으로 자유로운 터널면 때문에 추가적인 응력이 절리면 상에서 전단응력으로 집중되어 비가역적인 전단팽창이 크게 발생한 것으로 확인되었다. 본 논문의 결과는 고준위 방사성 폐기물 심층처분장의 계획 및 건설 시 굴착손상영역, 열, 빙하, 지진 시나리오를 고려한 성능 평가의 기초 자료로 활용될 수 있으며, 투수량계수와 같은 천연방벽의 수리학적 인자가 시간에 따라서 변화할 수 있다는 것을 제시하였다.

High-level radioactive waste produced from nuclear power plant continuously emits radioactivity and heat after it has been disposed of, so it should be isolated from the biosphere for a sufficient amount of time. To dispose of high-level radioactive waste, deep geological repository systems consisting of engineered barriers and natural barriers are suggested. As a performance assessment, numerous experiments and numerical simulations are performed to characterize the long-term behaviors of engineered and natural barriers. Especially, the transmissivity of natural barriers is considered an important parameter related to radionuclide leakage. Therefore, the relation between transmissivity and a variety of thermal, hydraulic, and mechanical factors are raised, and quantitative studies to determine the effects on transmissivity are suggested for performance assessment. This thesis aims to quantify the transmissivity of natural barriers, especially crystalline fractured rock, under a variety of factors possibly happening during the construction and operation of a geological repository. Three-dimensional discrete element models were constructed based on the rock and fracture characteristics extracted from the Äspö Hard Rock Laboratory in Sweden to describe the realistic behaviors of the crystalline fractured rock. The effects of excavation damage zone, thermal loading, glaciation, and earthquake were applied to numerical models as the factors disturbing the fracture transmissivity. The excavation damage zone was considered the effect of stress re-distribution inducing the fracture normal deformation and shear slip. The heat sources, glacial boundary conditions, and dynamic load were separately applied to the tunnel model after the stress re-distribution to determine the transmissivity changes. The transmissivity changes induced by the excavation damage zone appeared as three-order increases due to the combined effects of normal opening/closure and shear dilation, depending on the fracture orientation. Thermal, glacial, and earthquake scenarios included the loading and unloading cycles and recovered the initial stress conditions after each scenario ended. Irreversible transmissivity increases were found for several fractures under thermal, glacial, and earthquake effects. According to the transient analysis of transmissivity and stress path, the normal closure dominated the transmissivity during loading cycles, while the irreversible effects of shear dilation were revealed after the dissipation of normal loads. These reversible normal deformations and irreversible shear dilations appeared depending on the geometrical characteristics of fractures. To define the relation between the geometry of fractures and transmissivity changes, the uniformly distributed joints were applied on the tunnel models to analyze the effects of excavation damage zone, thermal loading, glaciation, and earthquake. The transmissivity increases were larger on joints parallel to the direction of maximum horizontal stress due to the absence of effects from the largest principal stress inducing the normal stress on fractures. At the vicinity of the tunnel, the joints that are slightly inclined from the tunnel surface accompanied a large amount of permanent transmissivity increase, because the additional stresses were converted to the shear stress on fractures due to the mechanically free surface. The results extracted in this research can be the basis data for performance assessments of geological repositories under excavation damage zone, thermal loading, glaciation, and earthquake scenarios which can happen during the lifetime of repositories. This thesis implies that the hydraulic parameters of natural barriers should be considered as dynamic variables that change as a repository operates.