2.pdf1.06MB

1. 서 론

 중·저준위 방사성폐기물과는 달리 고준위 방사성폐기물은 인간 및 생태계에 매우 위협적인 양의 방사능을 갖는다. 이와 같은 고준위 폐기물을 생태계와 격리시키기 위한 심지층처분장에는 다중방벽의 개념을 도입하여 폐기물로부터 나오는 유독물질의 이동을 늦추거나 차단한다. 특히 천연 방벽인 모암은 낮은 수리전도도와 높은 수착 능력으로 핵종의 이동을 지연시키는 역할을 한다. 하지만 암반에는 수많은 균열이 있는데 이는 지하수의 이동통로의 역할을 함과 동시에 방사성 핵종을 생태계로 이어주는 통로의 역할을 하게 된다. 균열을 따라 이동하는 방사성 핵종들은 암반 매질, 균열 충진 물질, 다양한 종류의 콜로이드, 유기물, 미생물 등과 상호작용을 함으로써 그 이동성이 빨라지거나 늦춰질 수 있다. 특히 콜로이드의 구성종류 중 벤토나이트는 많은 핵종들에 대해 뛰어난 수착능을 가지고 있음을 고려할 때 점토입자들은 핵종 이동에 영향을 미칠 수 있다[1]. 콜로이드는 지하수가 흐르는 어느 곳에나 존재하고 지하수와 폐기물 처분장의 완충재와의 반응으로 생성될 수 있으며 이들은 방사성 핵종의 이동을 촉진시킨다[2]. 실례로 미국 Nevada 주의 실험 장소에서 1.3 km 떨어진 곳에서 플루토늄이 발견된 바 있다[3]. 이는 지하수에서 용질 이동시 핵종의 수착반응을 고려한다면 일어나기 어려운 일이다. 하지만 콜로이드는 지하수 흐름보다 빠르게 이동하고 또 어떤 핵종은 콜로이드에 흡착되어 이동할 수 있기 때문에 가능한 일이다. 즉, 위에서 발견된 플루토늄은 콜로이드가 플루토늄 이동을 촉진시켰기 때문에 먼 거리까지 도달한 것이다. 이와 같이 콜로이드는 핵종의 이동을 촉진시키기 때문에 고준위 폐기물 처분의 안전성을 확보하기 위하여 이와 관련된 실험을 진행하여 핵종 이동에 대한 정보를 수집할 필요가 있다.

 고준위 폐기물은 처분장의 폐쇄 후에도 장기간 안전성이 보장되어야 하는데 실험만으로는 장기간 안전성 보장에 어려움이 따르기 때문에 이를 보완해 줄 수치모델이 동반되어야 한다. 이미 해외에서는 스웨덴의 Äspö와 스위스의 GTS 실험장 같이 핵종 이동 실험과 수치모델 개발이 수행되고 있는 실정이다. Äspö의 연구 프로젝트 중 삼차원 추적자 이동실험은 현장에서 핵종이동 실험을 수행하는 팀, 실험을 해석하고 예측하기위한 이동모델을 개발하는 팀으로 분반해서 연구하고 있고 이미 1, 2단계 연구가 진행되어 2007년에 최종보고서가 발간되었다. 스위스 GTS에서도 콜로이드 형성 및 이동 실험을 통하여 핵종이동에 대한 성능평가에 사용될 컴퓨터 코드의 검증 및 개선을 위해 현재 2단계(2008~2013)연구를 진행 중에 있다. 우리나라에서도 지난 2006년 11월 한국원자력연구원 부지 내에 고준위폐기물 처분을 위한 소규모 지하연구시설인 KURT(KAERI Underground Research Tunnel)이 건설 되어 활발한 연구가 진행 중에 있다[4]. 하지만 KURT는 약 100 m 수준으로 심도가 낮고 방사성 물질을 이용할 수 없다는 단점이 있다. 또 실제 처분장 수준의 연구시설을 국내에 건설하기에는 부지와 예산확보, 장기간의 건설기간 등의 문제점이 있기 때문에 국제공동연구 참여와 수치모델 개발의 필요성이 있다.

 따라서 본 연구에서는 단일 균열 암반에서 콜로이드가 핵종이동에 미치는 영향을 규명하기 위하여 수치모델을 제시하고 가상의 주어진 조건에서 모델의 실행 예를 기술하였다.

2. 연구 방법

2.1 개념 모델

본 연구에서는, 단일 균열암반에서 핵종과 콜로이드의 복합반응을 고려한 이동현상을 재연하기 위해서 여러 가지 현상들을 모델 개발에서 고려하였고, 이는<Fig. 1>에서 도식화되었고 <Fig. 1>에서 사용된 축약어는 <Table 1>에 정리하였다. 이러한 현상들은 크게 균열 공간, 균열 표면, 암반 공극, 그리고 암반 매질에서 일어나는 현상들로 분류 할 수 있다. 고준위 폐기물 처분장으로부터 유출된 핵종은 완충재로부터 유출되는 콜로이드와 함께 암반의 균열을 따라 이동하게 된다. 핵종은 균열을 따라 단독 혹은 콜로이드에 흡착되어 이동할 수 있고 균열 표면에 흡착 혹은 균열 표면에 흡착된 콜로이드에 흡착될 수 있다. 핵종은 암반 공극으로의 확산이 될 수 있으나 콜로이드는 일반적인 암반 매질의 공극의 크기보다 그 크기 때문에 공극으로의 확산을 고려하지 않는다[5]. 방사성 핵종은 시간에 따라서 그 방사능이 줄어들기 때문에 이도 같이 고려해 주었다.

2.2 수학 모델

2.2.1 이동 지배방정식

 단일 균열 암반에서의 이동 지배방정식은 크게 콜로이드의 이동 지배방정식과 핵종(용질) 이동 지배방정식, 두 가지로 나누었다. 이들 지배 방정식은 단일 균열에서 1차원 이류-분산 방정식(advection - dispersion equation)에 콜로이드와 핵종의 동역학적(kinetic) 반응들이 더하여진다.

Fig. 1. Conceptual model for colloid-radionuclide interaction in fractured rock.

Table 1. Explanations of Abbreviations for the reaction terms used in <Fig. 1>

2.2.1.1 콜로이드 이동 지배방정식

1차원 이류-분산 방정식에 콜로이드의 동역학적 반응들을 더한 콜로이드의 이동 지배방정식은 다음과 같다. 

여기서, F_c는 암반 균열에서 콜로이드의 농도 (g/m³), 는 암반 균열에서 지하수 유속 (m/s), D_c는 암반 균열에서 콜로이드의 수리분산 계수 (m²/s), t 는 시간 (s), 그리고 z 는 암반 균열 방향으로의 거리 (m) 이다. 여기서V_c는 콜로이드의 유속으로 지하수보다 평균적으로 약 1.1배 빠르다고 가정하였다. 암반 균열에서 콜로이드의 수리분산 계수(D_c)는 유속에 의한 이산(advection)으로 다음과 같이 정의된다.

 여기서, α_c 는 암반 균열에서 콜로이드의 분산계수(dispersivity) (m) 이다. 또한, 암반 균열 표면에 흡착된 콜로이드의 농도 변화는 이류 및 분산에 영향을 받지 않는다고 가정하면 다음 식과 같다.

 여기서, S_c는 암반 균열 표면에 흡착된 콜로이드의 농도(g/m³), 그리고 는 암반 균열 간극의 거리 (m) 이다.

식(1)과 (2)에서 사용된 동역학적 반응들의 속도를 나타낸 Q (g/m³/s) 항들의 정의는 다음과 같다.

 는 암반 균열에 존재하는 콜로이드가 암반 균열 표면으로 흡착하는 속도, 는 암반 균열 표면에 흡착된 콜로이드가 암반 균열로 탈착하는 속도이다. 이와 같은 동역학적 반응들은 <Fig. 1>에 형상화 되어있다.

2.2.1.2 핵종 이동 지배방정식

 콜로이드는 입자가 크기 때문에 암반 균열부에서만 이동을 하지만 지하수에 용존하는 핵종은 암반 공극으로 확산을 한다. 게다가 핵종의 방사능 붕괴까지 고려해야하므로 핵종의 이동 지배방정식은 콜로이드의 이동 지배방정식보다 약간 더 복잡하다. 암반 균열에서의 핵종 농도에 대한 방정식은 다음과 같다.

 여기서, F_r  는 암반 균열에서 핵종의 농도 (g/m³), 그리고 D_r 는 암반 균열에서 핵종의 수리분산 계수 (m²/s) 이다. 핵종의 수리분산 계수(D_r )는 콜로이드의 수리분산 계수와 비슷하지만 분자의 확산 항이 더해지면서 다음과 같이 정의된다 [6].

 여기서, σ_r는 암반 균열에서 핵종의 분산계수 (m), 그리고 D_r^*는 물과 같은 자유공간에서 핵종의 분자확산 계수 (m²/s) 이다.

 암반 균열 표면에 흡착된 핵종의 농도 변화가 이류와 분산에 영향을 받지 않는다고 가정하면 다음 식과 같다.

 여기서, S_r는 단위면적당 암반 균열 표면에 흡착된 핵종의 농도 (g/m²) 이다.

 암반 균열에서 이동하는 콜로이드에 흡착된 핵종은 콜로이드의 이동 특성에 따라 이동하게 되므로 식으로 정리하면 식(5)와 같다.

여기서, σ_Fcr는 암반 균열에 존재하는 콜로이드 표면에 흡착된 핵종의 질량비이다.

 다음으로 암반 균열 표면에 흡착된 콜로이드에 흡착되는 핵종에 대한 식은 이류와 분산에 영향을 받지 않는다고 가정하면 다음과 같다.

 여기서, σ_Scr 는 암반 균열 표면에 흡착된 콜로이드에 흡착된 핵종의 질량비이다.

암반 공극에서의 핵종 농도 변화는 화강암의 공극의 크기가 수 nm 이하임을 고려해서 지하수의 흐름은 고려하지 않고 오직 확산에 의해서만 핵종이 이동한다고 가정하였다. 핵종의 암반 공극으로의 확산은 식(7)과 같이 균열방향의 직각 방향으로 1차원적으로 나타내었다.
 

 여기서 P_r 는 암반 공극에서 핵종의 농도 (g/m³), D_Mr는 암반 공극에서 핵종의 수리확산 계수 (m²/s), 그리고 는 암반 균열과 수직방향으로의 거리 (m) 이다. 암반 공극에서 핵종의 수리확산 계수(D_Mr)는 분자확산만을 고려하여 다음과 같이 정의된다[6].

여기서, T는 암반 공극의 굴곡비(totuosity)로 경로의 직선거리와 실제 유체가 이동한 거리의 비의 제곱이다.

 마지막으로, 암반 매질 표면에서의 핵종 농도 변화는 다음과 같다.

여기서, M_r는 암반 매질 표면에 흡착된 핵종의 농도 (g/m³), 그리고 θ는 암반의 공극률이다.

 핵종의 이동 지배방정식에서 사용된 동역학적 반응들은 Q로 나타내었고, 이는 <Fig. 1>에 형상화 되어있다.

2.2.2 동역학적 반응

 콜로이드와 핵종의 이동 지배방정식에서의 동역학적 반응들은 모두 1차 반응속도(1st-order kinetics)로 가정하였다.

2.2.2.1 암반 균열에서 암반 공극으로의 핵종 확산

 암반 균열에서 암반 공극으로의 핵종 이동은 농도차로 인한 확산(diffusion)에 의해서 이루어진다. 따라서 이러한 이동을 Fick‘s first law에 의해서 나타내었다.

2.2.2.2 콜로이드의 흡착과 탈착

암반 균열에서 균열 표면으로 콜로이드의 흡착과 탈착의 속도는 암반 균열에 존재하는 콜로이드의 농도에 1차원적으로 비례한다고 가정하였다. 식(10)과 식(11)은 각각 콜로이드의 흡착과 탈착 속도에 관한 식이다.
 

여기서, 는 콜로이드가 암반 균열에서 암반 균열 표면으로 흡착하는 속도 계수 (1/s), 그리고 는 콜로이드가 암반 균열 표면에서 암반 균열로 탈착하는 속도 계수 (1/s)이다.

콜로이드의 흡착과 탈착은 콜로이드 표면에 흡착된 핵종의 이동에 영향을 미치게 된다. 이러한 흡착(식(12))과 탈착(식(13))에 대하여 아래와 같이 나타내었다.
 

2.2.2.3 핵종의 흡착과 탈착

 암반 균열에서 표면으로 핵종의 흡착과 탈착의 속도는 암반 균열에 존재하는 핵종의 농도와 암반 표면에 존재하는 핵종의 농도에 각각 1차원적으로 비례한다고 가정하였다. 식(14)는 핵종의 흡착, 식(15)는 핵종의 탈착 속도를 나타낸다.

 여기서, 는 핵종이 암반 균열에서 암반 균열 표면으로 흡착하는 속도 계수 (1/s), 그리고 는 핵종이 암반 균열 표면에서 암반 균열로 탈착하는 속도 계수 (1/s) 이다.

 암반 균열에서 핵종이 콜로이드로 흡착하는 정도는 핵종과 콜로이드의 농도 모두를 고려해 주어야 한다. 따라서 암반 균열에서 핵종이 콜로이드에 흡착하는 속도를 암반 균열에 존재하는 콜로이드에 흡착하는 속도(식(16))와 암반 균열 표면에 흡착된 콜로이드에 흡착하는 속도(식(17))에 대해서 아래와 같이 정의하였다.

 여기서, 는 핵종이 암반 균열에서 암반 균열에 존재하는 콜로이드 표면으로 흡착하는 속도 계수 (1/s), 는 핵종이 암반 균열에서 암반 균열 표면에 흡착된 콜로이드 표면으로 흡착하는 속도 계수 (1/s) 이다.

 암반 균열에 존재하는 콜로이드로부터 핵종의 탈착 속도는 흡착된 핵종의 농도에 1차원 적으로 비례한다고 가정하여 암반 균열에 존재하는 콜로이드(식(18))와 균열 표면에 흡착된 콜로이드(식(19))에 대해서 각각 다음과 같이 정의하였다.

 여기서, 는 핵종이 암반 균열에 존재하는 콜로이드 표면에서 암반 균열로 탈착하는 속도 계수 (1/s), 그리고 는 핵종이 암반 균열 표면에 흡착된 콜로이드 표면에서 암반 균열로 탈착하는 속도 계수 (1/s) 이다.

 암반 공극에서 핵종 흡착과 탈착 속도는 암반 공극에 존재하는 핵종의 농도와 암반 매질에 흡착된 핵종의 농도에 1차원적으로 비례한다고 가정하였다. 흡착(식(20))과 탈착 (식(21)) 속도는 아래와 같다.

 여기서, 는 핵종이 암반 공극에서 암반 매질 표면으로 흡착하는 속도 계수 (1/s), 그리고 는 핵종이 암반 매질 표면에서 암반 공극으로 탈착하는 속도 계수 (1/s)이다.

2.2.2.4 핵종의 감쇠

 방사성 핵종은 반감기를 가지므로, 핵종의 감쇠를 다음과 같이 암반 균열과 암반 공극, 그리고 콜로이드 표면에서 각각 1차 반응으로 가정하였다.

여기서, λ_r 는 핵종의 감쇠 계수 (1/s) 로, 핵종의 반감기(t_r^1/2  )를 이용해 다음과 같이 계산된다.
 

2.3 수치해석

 앞에서 설명된 이류-분산(advection-dispersion) 방정식과 동역학적(kinetic) 반응식들의 합은 비선형 방정식이기 때문에 이론적 해석해를 구하는 것이 매우 어렵거나 불가능하다. 따라서 적절한 수치해석 방법을 사용하여 수치해를 구할 수 밖에 없다. 본 연구에서는 이류-분산 방정식과 동역학적 반응식들을 별개로 계산하는 Operator Splitting Method를 사용하여 해당 방정식의 수치해를 구하였다.

2.3.1 Operator Splitting Method(OSM)

Operator Splitting Method는 콜로이드와 핵종의 동역학적 반응식을 포함한 이류-분산 방정식의 해석해를 구하는데 있어서 계산 효율을 상당히 향상시켜 준다. 이 방법은 이류-분산 방정식과 동역학적 반응식을 별개의 계산 방법으로 각각계산하는 방법이다.

 일반적으로 m개의 용질에 대한 반응-이동 방정식은 다음과 같이 주어진다.

 여기서, C_e  는 용질 e의 농도, L(C_e, z )는 용질 e의 이류-분산항으로 다른 용질의 영향을 받지 않고, f_e (C₁, C₂,...C_e,...,C_m)는 용질 e의 반응항으로 다른 용질의 농도에 영향을 받는다.

 이에 대한 Operator Splitting Method는 계산 순서와 반복 계산 여부에 따라 표준 순차 무반복(Sequential non-iterative; SNI) 방식, Strang의 분리(Strang-splitting) 순차 무반복(Sequential noniterative; SNI) 방식, 표준 순차 반복(Sequential iterative; SI) 방식, 그리고 대칭 순차 반복(Symmetric sequential iterative) 방식으로 나누어진다. 본 연구에서는 모델의 안정성에 문제없이 정확한 결과를 도출해 낼 수 있는 Strang의 분리 SNI방식을 채택하여 이동-반응 모델의 수치해를 구하였다[7].

 Strang의 분리 SNI 방식은 단위 시간 간격을 반으로 나누어서 반응식들을 제외한 용질 이류-분산 방정식을 각 용질 e에 대해서 다음과 같이 먼저 계산한다.

 여기서, C_e^{n + 1/2, *}  는 시간 (n + 1/2)·Δt에서 용질 이류-분산 방정식에 의해 계산된 용질 e의 농도(임시값)이다. 식(29)에서 계산된 C_e^{n + 1/2, *}를 이용하여 각 용질 e에 대한 반응식을 순차적으로 계산한다.

 여기서,C_e^{n + 1 , * *}는 시간 (n + 1)·Δt 에서 용질 반응에 의해 계산된 용질 e의 농도(임시값)이다. 식(30)에서 계산된 C_e^{n + 1 , * *} 를 이용하여 나머지 시간 간격 반에 대한 용질 이류-분산 방정식을 각 용질 e에 대해서 다음과 같이 다시 계산한다.

2.3.2 콜로이드 및 용질 이류-분산

2.3.2.1 암반 균열에서 콜로이드의 이류-분산

 암반 균열에서 콜로이드의 이류-분산 방정식에 대하여 본 연구에서는 경계조건으로 유입구에서는 콜로이드가 고정농도로 주입된다고 가정하였고 유출구에서는 0 구배 농도 조건을 사용하였다. 초기조건으로는 최초 암반 균열에 콜로이드가 없다고 설정하였다. 초기조건과 경계조건을 식으로 표현하면 아래와 같다.

2.3.2.2 암반 균열에서 용질의 이류-분산

 암반 균열에서 용질 이류-분산 방정식에 대하여 본 연구에서 사용한 경계조건으로, 유입구에서는 고정농도 조건이, 그리고 유출구에서는 0 구배농도 조건이 설정되었다. 초기 조건으로는, 최초 암반 균열에 용질이 존재하지 않는다고 가정하였다. 이들을 수식으로 표현하면 다음과 같다.

2.3.2.3 암반 균열에서 콜로이드 표면에 흡착된 용질의 이류-분산

 암반 균열에서 콜로이드 표면에 흡착된 용질의 이류-분산방정식에 대하여 본 모델에서 사용한 경계조건으로, 유입구에서는 고정농도 조건이, 그리고 유출구에서는 0구배농도 조건이 설정되었다. 초기조건으로는, 최초 암반 균열에 콜로이드 표면에 흡착된 용질이 존재하지 않는다고 가정하였다. 이들을 수식으로 표현하면 다음과 같다.

2.3.2.4 암반 공극에서 용질의 확산

 암반 공극에서 용질 확산 방정식에 대하여 본 연구에서 사용한 경계조건으로, 유입구에서는 고정농도 조건으로 각 시간별 암반 균열에서의 용질 농도가, 그리고 유출구에서는 0 구배농도 조건이 설정되었다. 초기조건으로는, 최초 암반 공극에 용질이 존재하지 않는다고 가정하였다. 이들을 수식으로 표현하면 다음과 같다.

3. 결과 및 논의

3.1 모델 실행 영역

 모델 실행 영역은 <Fig. 2>와 같이 평면 단일 암반 균열의 중심축을 기준으로 한쪽 방향만을 고려하여, 균열을 따라 0.2 m (=Z) 그리고 균열과 직각 방향으로 0.001 m (=X) 로 정하였다. 본 모델은 혼합 수직 1차원 모델로, 암반 균열에서는 z-축 방향으로, 암반 매질에서는 x-축 방향으로 계산되었고, y-축 방향으로는 0.052 m의 고정된 너비를 가지는 것으로 가정하였다. 모델 실행을 위한 시간적·공간적 상수들은 <Table 2>에 정리하였다.

3.2 모델 검증

 위에서 기술한 핵종/콜로이드 복합이동의 메커니즘을 포함하는 용질의 이류-분산 방정식의 해석해를 구하는 것은 불가능하거나 매우 어렵다. 따라서, 본 연구에서는 Tang 등[8]이 보고한 단일 암반 균열에서 용질의 이동 및 분산에 대한 해석해와 물질 수지(mass balance) 검사를 통하여 개발된 모델의 검증을 시도하였다.

3.2.1 해석해를 통한 모델 검증

Tang 등[8]의 해석해는 암반 균열에서 이류-분산뿐만 아니라 Fick의 제1법칙(Fick’s first law)을 따르는 암반 매질로의 확산과 용질의 선형 흡착 평형 등을 포함하고 있다. 이러한 용질의 흡착을 저해상수(retardation factor; 암반 균열에서는 R, 암반 매질에서는 R_Mr  로 나타내어 정리된 해석해는 다음과 같다.

Fig. 2. The domain of virtual model.

Table 2. Coefficient values determining temporal and spatial domain of model.

 이고, C는 암반 균열 및 공극에서의 용질 농도로써 본 모델에서는 F_r와 P_r에 해당하겠다. 식(32)에서 무한적분 부분은 MATLAB의‘Mathematics toolbox에’ 포함된 adaptive Gauss-Kronrod quadrature 함수를 이용하여 수치 해석적으로 계산되었다.

 해석해의 조건을 고려하여, 모델 실행시 콜로이드의 이류-분산과 콜로이드의 반응식들은 제외되었고, 오직 용질의 이류-분산과 암반 매질로의 확산만을 고려하였으며 용질의 흡착을 고려하지 않아 저해상수(R 과 R_Mr  ) 는 1로 가정하였다.

아래의 <Fig. 3>에서 모델 실행 결과와 해석해를 유량을 변화시켜가면서 비교하였다. 모델 검증 단계에서는 계산의 편의를 위해 모델의 시·공간적 영역을 결정하는 상수 값들을 변화시켰는데, 암반 균열 간극의 절반(b/2)을 0.0005 m로, 암반 매질방향으로의 거리(X )를 0.003 m ( dx = 0.0005m)로 확장하였고, 시간(T )은 24 h 동안만 확인하였다. 이를 바탕으로 암반 균열에서 유체의 체류시간(retention time)은 유량에 따라서 각각, q = 0.05mL/min 일 때
 

q = 0.01mL/min 일 때 8.7 h, 그리고 q = 0.005mL/min 일 때 17.3 h으로 계산된다. 그 결과, 모델 실행 시간(24h)이 체류시간의 14배에 달하는 조건인 q = 0.05mL/min 일때는 용질이 균열 전반에 포화되었으며(<Fig.3(A)>), 모델 실행 시간이 체류시간의 3배 이하인 경우(q = 0.01, 0.005mL/min )에는 유출부에서 용질의 파과(break-through)가 일어나지 못했다(<Fig. 3(C)>) .q = 0.005mL/min 일 때도 완전한 파과가 일어나지 못했는데, 이는 암반 매질로의 확산이 함께 고려되었기 때문이다. 전반적으로, 모든 유량 조건에서 암반 균열에서의 용질 농도분포 결과는 해석해 결과와 비교하여 피어슨 상관계수의 제곱값(r² )이 0.99 이상으로 모델의 정확성이 입증되었다. 하지만, 암반 매질 공극에서의 용질 농도분포는 모델 결과값과 해석해 값이 다소 상이하게 나타났는데, 이는 암반 매질 방향으로의 모델 영역 이산화(discretization)정도(dx)가 암반 균열과 매질의 경계면에서 용질 농도의 급격한 변화를 나타내기에 부족했기 때문에 나타난 결로 사료된다.

Fig. 3. Model execute results for consequences of different flow rates and comparisons with analytical solution of Tang et al.[8] (A) Spatial distribution of the solute concentration in rock fracture, (B) Spatial distribution of the solute concentration in the direction of outlet of rock fracture toward the matrix rock, (C) Temporal distribution of the solute concentration at outlet of the rock frature.

Table 3. Mass balance check for the model developed in this study.

3.2.2 물질 수지 (mass balance) 검사를 통한 모델 검증

 본 모델을 검증하는 과정에 있어서 해석해는 몇 가지 가정을 통한 간단한 조건에서만 본 모델과의 비교가 가능하므로 추가적인 검증 절차가 필요하다. 따라서 본 논문에서 개발한 모델을 검증하기 위한 추가적인 방법으로 물질 수지 검사를 실행하였다. 물질 수지 검사는 용질의 이동, 용질의 흡착, 용질의 확산, 콜로이드의 이동, 콜로이드의 흡착에 대한 수식 모델 각각의 결과를 알아보기 위해 다음과 같이 다섯 가지 조건에서 실행하였다.

(1) [용질의 이동 방정식]
(2) [(1) + 용질의 흡착과 탈착에 대한 항(term)]
(3) [(2) + 용질의 암반 공극으로의 확산에 대한 항]
(4) [(3) + 콜로이드의 이동 방정식]
(5) [(4) + 콜로이드의 흡착과 탈착에 대한 항(term)] 

위 다섯 가지 조건에 대한 물질 수지 검사 결과를 <Table3>에 정리하였고 이를 통하여 본 모델에 대한 정확성을 다시 입증 하였다.

물질 수지 검사 중 용질의 이동 방정식에 대한 결과와 용질의 이동 방정식과 용질의 암반 공극으로의 확산에 대한 항을 함께 고려한 결과만 ~99 %를 나타내는 것을 확인 할 수 있다. 여기에는 두 가지 원인이 예측 되는데 첫 번째 예측 원인은 컴퓨터 계산으로 인한 절단오차의 발생이다. 두 번째 예측 원인은 해당 모델에서 용질의 질량 유입량에 대한 식에 분산에 관한 항이 없다는 것이다. 따라서 이 식에 따르면 실제 유입량보다 용질 유입량이 적게 계산되므로 물질수지 검사에서 다소 오차가 생길 수 있다. 하지만 용질의 유입량에 관하여 분산 효과는 무시 할 수 있을 정도로 작기 때문에 큰 오차는 나타나지 않았고 전반적으로 모든 조건에서 유입량 대비 회수율이 100%에 근접한 결과를 나타내어 용질 이동 측면에서 개발된 모델의 정확성 및 안정성을 입증할 수 있었다.

4. 결 론

 본 연구에서는 단일 균열 암반에서 콜로이드가 핵종이동에 미치는 영향을 규명하기 위하여 수치모델을 새롭게 자체 개발하였다. 콜로이드의 이동 지배방정식, 핵종 이동 지배 방정식, 그리고 콜로이드와 핵종의 동역학적 반응들을 동시에 고려하였다. 콜로이드와 핵종의 반응-이동 지배 방정식을 풀기 위한 수치해석 방법으로는 Operator Splitting Method 중 Strang의 분리 SNI 방식을 채택하여 컴퓨터 계산 효율 증대를 도모하였고 이는 MATLAB을 이용하여 코드화 되었다.

이렇게 개발된 수치모델은 용질(핵종)의 이동 및 분산만을 고려한 해석해를 통한 검증과정에서 피어슨 상관계수의 제곱값(r² )이 0.99 이상으로 나타나 모델의 정확성이 입증되었다. 또한 본 모델의 추가 검증을 위한 물질 수지 검사를 실행하여 모델의 정확성을 한번 더 입증하였다. 본 연구는 실험실 칼럼 실험을 염두에 둔 가상 조건에서의 모델 실행 결과로써 향후 본 연구에서 개발된 수치모델과 실험을 비교함으로써 수치모델의 정확성 및 안정성을 한층 더 확실히 입증할 수 있을 것으로 보인다.