Numerical analysis of FEBEX at Grimsel Test Site in Switzerland

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

스위스 Grimsel Test Site에서 수행된 FEBEX 현장시험에 대한 수치해석적 연구

벤토나이트 완충재에서의 열-수리-역학적 복합거동을 예측하기 위해 TOUGH2-MP/FLAC3D 시뮬레이터를 기반으로 개발된 Barcelona basic 모델(BBM) 해석모듈의 현장 적용성을 검토하고자 국제공동연구 DECOVALEX-2019 Task D에 참여하여 스위스 Grimsel Test Site의 현장시험(full-scale engineered barriers experiment, FEBEX) 모델링을 수행하고 현장시험에서 계측된 히터 파워, 온도, 상대습도, 응력, 포화도, 함수율 그리고 건조밀도를 계산 값과 비교하였다. 수치해석을 이용하여 시간에 따른 히터 파워와 온도 변화는 전반적으로 잘 재현되었지만, 히터 1과 히터 2에서의 파워 차이를 계산할 수는 없었으며 이를 개선하기 위해서는 FEBEX 터널 주변에 분포하는 황반암과 시험장치 및 벤토나이트 블록의 설치 공정을 반영할 필요가 있을 것으로 판단된다. 상대습도 변화와 분포 역시 전반적으로 잘 모사되었으나, 수치해석에서 히터 부근에서의 재포화과정이 상대적으로 빠르게 진행된 것으로 보아 수리모델에 대한 일부 수정이 필요할 것으로 보인다. 현장시험에서는 벤토나이트 완충재와 암반 사이에 틈이 존재하지만 수치해석에서는 완벽하게 접촉하고 있는 것으로 가정하였기 때문에 운영 초기의 응력 변화는 다소 차이를 보였지만, 전반적으로 유사한 경향을 보이는 것으로 나타났다. 해체 이후 측정한 포화도, 함수율, 그리고 건조밀도의 분포 역시 전반적으로 잘 재현되었지만, 건조밀도가 터널 중심과 히터부근에서 조금 크게 계산되어 벤토나이트 블록의 투수계수가 상대적으로 작은 값으로 반영되어 포화도와 함수율이 작게 계산된 것으로 보이며, 이를 개선하기 위해서는 건조밀도에 따른 투수계수 모델에 일부 수정이 필요할 것으로 판단된다. 본 연구의 결과를 토대로 수치모델을 수정하고 추가적인 연구를 수행한다면, 보다 나은 해석 결과와 벤토나이트 완충재에서의 THM 복합거동을 좀 더 현실적으로 예측할 수 있을 것으로 판단된다.