Understanding the long-term geochemical evolution of engineered barrier system is crucial for conducting safety assessment in high-level radioactive waste disposal repository. One critical scenario to consider is the intrusion of seawater into the engineered barrier system, which may occur due to global sea level rise. Seawater is characterized by its high ionic strength and abundant dissolved cations, including Na, K, and Mg. When seawater infiltrates an engineered barrier, such dissolved cations displace interlayer cations within the montmorillonite and affect to precipitation/ dissolution of accessory minerals in bentonite buffer. These geochemical reactions change the porewater chemistry of bentonite buffer and influence the reactive transport of radionuclides when it leaked from the canister. In this study, the adaptive process-based total system performance assessment framework (APro), developed by the Korea Atomic Energy Research Institute, was utilized to simulate the geochemical evolution of engineered barrier system resulting from seawater intrusion. Here, the APro simulated the geochemical evolution in bentonite porewater and mineral composition by considering various geochemical reactions such as mineral precipitation/dissolution, temperature, redox processes, cation exchange, and surface complexation mechanisms. The simulation results showed that the seawater intrusion led to the dissolution of gypsum and partial precipitation of calcite, dolomite, and siderite within the engineered barrier system. Additionally, the composition of interlayer cation in montmorillonite was changed, with an increase in Na, K, and Mg and a decrease in Ca, because the concentrations of Na, K, and Mg in seawater were 2-10 times higher than those in the initial bentonite porewater. Further studies will evaluate the geochemical sorption and transport of leaked uranium-238 and iodine-129 by applying TDB-based sorption model.
During and after the construction of LILW disposal facilities, the decrease of groundwater head potential has been monitored. In addition, an increase of the electrical conductivity (EC) has been observed in several monitoring wells installed along the coastal coastline. Monitoring activity for groundwater head potential and hydrogeochemical properties is important to reduce the uncertainty in the evaluation of groundwater flow characteristics. However, the data observed in the monitoring wells are spatial point data, so there is a limit to the dimension. Several researchers evaluated groundwater head potential changes and seawater intrusion (SWI) potential for disposal sites using groundwater flow modeling. In case of groundwater flow modeling results for SWI, there is a spatial limit in directly comparing the EC observed in the monitoring wells with the modeling results. In a recent study, it was confirmed that the response of the long-range ground penetraiing radar (GPR) system was severely attenuated in the presence of saline groundwater. In order to reduce the spatial constraint of the groundwater monitoring wells for SWI, the characteristics of SWI within the disposal facility site by using the the results of a recent study of the long-range GPR system were investigated and evaluated in this study.
우리나라 서, 남해안에 지하수의 해수침투 현상을 관측하기 위해 설치된 관측망 자료를 이용하여, 해당지역별 대수층을 통한 해수침투의 유형을 분석하였다. 분석에 이용된 관측정은 25개 해안지역 소유역의 총 45개 암반관정으로, 지하수위, 수온, 전기전도도를 대상으로 기본통계분석, 상관성 분석 및 변동유형 분류를 수행하였다. 분석 결과 지하수위의 경우 강우 영향을 받거나 양수에 의한 수위강하가 나타나는 관측정에서 큰 폭으로 변동하고 있다. 지하수온은 대부분의 경우 변동폭이 0.2˚C 이내로 안정적인 특징을 보여주고 있으며, 전체적으로 평균 온도가 15˚C 이상으로 나타났다. 전기전도도의 경우 상대적으로 변동폭이 크고 불규칙한 특징을 가지고 있으며, 평균값이 2,000 μS/cm 이내의 관측정이 28개, 10,000 μS/cm를 초과하는 관측정은 9개소로 나타났다. 교차상관도 분석에 의하면, 지하수위는 강우의 영향을 받아 변동하는 형태가 많았지만 수온과 전기전도도는 상대적으로 강우의 영향이 없는 것으로 나타났다. 조석의 경우에는 일부 관측정에서 강우에 의한 영향보다 높은 교차상관도가 나타났다. 본 연구에서는 장기관측자료를 변동형태에 따라 정상형, 추세형, 주기형, 충동형, 계단형, 경사형 등 6가지로 분류하였다. 지하수위의 경우 강우나 양수의 영향을 받는 충동형이 가장 많은 73.3%이며, 조석의 영향은 13.4%에 해당되었다. 지하수온의 경우 변동폭이 일정한 정상형태가 51.2로 가장 높게 나타났으며, 강우나 양수에 의한 충동형이 26.7%에 해당되었다. 전기전도도의 경우는 지하수위나 지하수온과 달리 추세형, 계단형, 경사형 등의 변동형태가 우세한 것으로 나타났다. 전체적으로 동일한 해안지역에 설치된 관측정별로 관측자료의 변동 특성이 상이한 경우가 나타남에 따라, 향후 각 소유역별로 관측정 자료의 정기적인 변동 경향성을 분석하는 것이 필요하다.
본 연구에서는 SEAWAT을 이용하여 여수지역의 해수침투 피해 면적을 파악하고, 기후변화 시나리오 적용에 따른 미래의 해수침투 피해 예상 면적을 산출하였으며 해수침투 피해 저감 대책의 피해면적 저감 능력을 분석하였다. 2015년 기준 여수지역의 해수침투 피해 면적은 14.90 km2로 나타났고, 기후변화 시나리오를 적용하였을 때 2099년 여수지역의 예상 해수침투 피해 면적은 RCP 4.5의 경우 19.19 km2이며, RCP 8.5의 경우 20.43 km2로 나타났다. 이에 대한 저감대책으로 인공함양을 고려하였을 때, 총 300 m3/d, 100 m3/d, 50 m3/d의 함양 시나리오를 설정하였을때 RCP 4.5의 경우 해수침투 면적은 평균 7.03%, RCP 8.5의 경우 8.32% 감소하였다. 물리적 차수벽 대책의 경우는 차수벽의 두께를 0.8 m, 1.3 m, 1.8 m 로 설정하였을 때 RCP 4.5의 경우 해수침투 면적은 평균 9.80%, RCP 8.5의 경우 10.30% 감소하였다. 본 연구는 연안지역의 해수에 의한 지하수 오염과 그에 수반한 2차적인 피해를 막기 위한 대비책을 결정하기 위한 정량적인 근거로서 사용될 수 있을 것으로 기대된다.
기후변화에 따른 해수면 상승으로 인하여 해안지역의 지하수계에 해수침투가 가중된다. 지하수의 염분농도가 증가하면 지하수면 상부의 불포화 토양에서도 염분 농도가 증가할 수 있으며, 이는 농경지에서 작물피해를 일으킬 수 있다. 해수면이 상승함에 따라 내륙의 지하수위도 함께 상승한다. 이는 불포화 토양층의 두께를 감소시켜 해안 저지대의 경작에 피해를 끼칠 수 있다. 본 연구에서 지하수 해수침투는 3차원 모델, 토양 염류화 평가는 연직 1차원 모델을 합성 적용하여 해안 농경지에 대한 해수면 상승 피해를 평가하는 방법을 개발하였다. 3차원 해수침투 모델에서 지하수면의 수위 와 농도분포를 계산하고 최상부 절점 중에서 염분 농도가 기준 값 이상인 절점에서 지하수면과 지표면 사이의 토양층에서 연직 1차원 모델링으로 토양층의 염분 농도와 불포화대 두께를 계산하였다. 농경지의 토양 염류화는 작물의 뿌리 심도에서 보통 작물의 생육한계 염분농도를 기준으로 판 단하였다. 개발된 모델링 방법을 가상의 간척농경지에 적용하였다. 해수면 상승자료로 IPCC의 RCP 4.5와 8.5 시나리오를 사용하였다. 평가 결과 는 2050년과 2100년에 대하여 제시하였다. 연구결과 대상지역에서 기후변화 시나리오 RCP 8.5에서 2100년에는 지하수 염류화 피해 면적은 간 척지 육지면적 대비 7.8%, 염류화 토양 면적은 6.0%, 불포화층의 두께가 뿌리심도보다 적은 지역의 면적은 1.6% 증가하는 것으로 분석되었다.
본 연구에서는 해수-담수-해안대수층의 비선형 상호작용을 직접 해석할 있는 PBM(Porous Body Model) 기반의 3차원 N-S Solver인 LES-WASS-3D ver 2.0을 적용하며, 해안대수층의 해수침투모의를 수행하였다. 이와 같은 N-S Solver를 적용한 해안대수층의 해수침투모의는 국내 최초 수행되는 것일 뿐만 아니라, 국외적으로도 찾아보기 어려운 새로운 수치해석방법이라고 할 수 있다. 먼저 적용하는 수치모델을 검증하기 위하여 해안대수층의 해수-담수 경계면에 관한 수리모형실험결과와 비교·검토 하여 수치모델의 타당성 및 유효성을 확인하였다. 그리고 해수위 및 지하수위 변화를 고려한 해안대수층 내의 해수침투모의를 수행하여 해수위-지하수위 차와 해수위의 비(△h/h)의 증가에 따른 해안대수층 내의 유동장 그리고 해수-담수 경계면 분포 특성에 관하여 논의하였다. 또한 기존의 비확산 수치모델에서 도출할 수 없었던, △h/h에 따른 해안대수층 내의 연직 염분농도로 부터 해수침투 특성을 파악하였으며, 최종적으로 지표화 할 수 있는 △h/h가 해안대수층 내의 해수침투거리에 미치는 영향에 대하여 분석하였다. 이 결과로부터 △h/h가 작을수록 해안대수층 내의 해수침투가 약해지는 메커니즘을 이해할 수 있었다.