• 지질공학
• /
• 제26권3호
• /
• pp.325-330
• /
• 2016

현재 경주 중저준위방사성폐기물 처분장(이하 '경주방폐장')에서는 2단계 표층처분시설이 계획중에 있으며, 포화대에 위치한 1단계 처분시설과는 달리 불포화대 상부에 위치하게 된다. 단열을 포함하는 불포화대의 특성상 지하수 및 용질의 대부분이 단열을 통해 이동할 것으로 예상된다. 따라서 불포화 암반 매질에 대한 정밀한 해석을 위하여 단열망 연속체와 암반 매질 연속체를 구분하여 해석하는 TOUGH2 전산코드의 meshmaker 모듈의 MINC 기법을 활용하였다. TOUGH2 MINC 기법의 기존 국내 연구 사례가 미미하여 본 연구에서는 MINC를 이용한 mesh 구성 방법에 대한 절차를 개발하였으며, 단열 연속체와 암반매질 연속체의 k-field를 생성하였다. 이와 같이 생성된 도메인은 향후 이중 연속체를 기반으로 경주방폐장의 지하수 유동 및 오염물질 이동 등에 활용될 뿐만 아니라 단열이 발달한 암반에서 단열-암반매질 연결성을 고려한 단열망 유동 특성을 분석하는데 참고가 될 것으로 기대한다.

• 터널과지하공간
• /
• 제34권2호
• /
• pp.127-142
• /
• 2024

수치해석 기법의 강건성은 다양한 모델링 조건에서 계산 성능이 유지되는 것을 의미하며, 새로운 해석기법 또는 수치코드는 벤치마크 테스트를 통해 강건성이 평가될 필요가 있다. TOUGH-FLAC 모델링 기법은 국내외적으로 이산화탄소 지중저장, 사용후핵연료 지층처분, 지열 개발 등 다양한 분야에 적용되었으며, 실험 계측자료, 다른 수치코드들과의 결과 비교를 통해 모델링 유효성이 분석되었다. 본 연구에서는 해석해를 갖는 열-수리-역학적 복합거동 문제를 토대로 TOUGH-FLAC 기법의 벤치마크 테스트를 수행하였다. 적용된 해석해는 완전히 포화된 지반에 점열원 작용 시 주변매질의 온도, 간극수압, 역학적 거동과 관계되며, 해석해와 수치모사 결과를 비교하여 TOUGH-FLAC 기법의 강건성이 평가되었다. 또한, 열-수리-역학 해석의 연계항, 유체 상변화, 시간증분이 복합거동 계산에 미치는 영향을 조사하였다.

• 터널과지하공간
• /
• 제30권4호
• /
• pp.335-358
• /
• 2020

본 연구에서는 TOUGH-FLAC 연동해석기법을 통해 단층의 수리역학적 거동을 평가할 수 있는 수치해석 모델을 제안하였다. 이는 국제공동연구 DECOVALEX-2019 Task B의 일환으로 수행되었으며, 불투수성 암반 내 유체 주입으로 인한 단층의 재활성을 예측하고 주변 암반의 수리역학적 안정성을 평가할 수 있는 해석모델을 개발하는 데에 그 목적이 있다. 본 연구에서는 TOUGH2 수리유동모델과 FLAC3D의 역학적 인터페이스 모델의 연동을 통해 단층의 역학적 거동을 보다 합리적으로 구현할 수 있는 해석기법을 제안하고, 벤치마크 해석과 스위스 Mont Terri 지하연구시설 현장시험에 적용하여 그 타당성과 유효성을 검증하였다. 개발된 해석 모델은 유체의 주입으로 인한 단층 내 압력 분포의 발달, 역학적 변형에 따른 수리간극의 변화, 변위와 응력 등 단층의 수리역학적 거동을 적절히 재현할 수 있는 것으로 나타났다.

• 터널과지하공간
• /
• 제30권1호
• /
• pp.39-62
• /
• 2020

본 연구에서는 불포화토의 열-수리적 변화로 야기되는 역학적 거동 특성을 분석하고 장기 거동을 예측하기 위해 Barcelona Basic Model(BBM) 해석 모듈을 TOUGH2-MP/FLAC3D에서 구동할 수 있도록 개발하였다. 기본적으로는 TOUGH-FLAC에 BBM 해석 모듈을 개발하기 위해 사용된 선행연구의 방법과 마찬가지로 FLAC에서 제공하고 있는 Modified Cam Clay Model(MCCM)을 기반으로 User Defined Model(UDM)과 FLAC3D에서 제공하는 FISH function을 이용하였다. 본 연구에서 개발한 BBM 해석 모듈에서는 평균유효응력뿐만 아니라 흡입력의 변화에 따른 소성변형률을 모두 고려하였으며, 평균유효응력 및 흡입력 증가에 따른 항복면의 변화를 모두 반영할 수 있도록 하였다. 개발된 BBM 해석 모듈을 FLAC3D 매뉴얼에 기술되어 있는 MCCM 예제, BBM을 처음 제안한 선행연구에 언급되어 있는 예제들, 그리고 스웨덴 SKB 보고서에 언급되어 있는 실험실 시험 결과를 이용하여 검증하였다. 뿐만 아니라, 효과적인 BBM의 파라미터 도출을 목적으로 개발된 Quick tools의 검증을 위해 수행된 일련의 모델링 시험들을 동일하게 수행하고 선행연구에서 보고된 Quick tools와 Code_Bright의 결과와 비교함으로써 개발된 BBM 해석 모듈을 검증하였다.

• 터널과지하공간
• /
• 제25권6호
• /
• pp.556-567
• /
• 2015

고준위방사성폐기물 처분장의 장기 안전성 및 안정성 평가를 위해서는 심지층 처분장 환경에서의 열-수리-역학적 복합 거동에 대한 평가가 필요하다. 신뢰도 높은 THM 해석 기법의 개발을 효과적으로 하기 위해 DECOVALEX 국제공동연구가 수행되고 있다. DECOVALEX-2015 Task B2에서는 일본 JAEA에서 계획중인 현장 THM 실험에 대하여 세계 각국의 연구진이 모델링을 실시하였다. 본 연구에서는 TOUGH2-FLAC을 연계한 THM 해석 기법을 개발하고 이를 이용하여 암반, 완충재, 히터로 구성되는 1차원 THM 모델링을 수행하였으며 해석 결과를 외국의 연구팀 결과와 비교하였다.

• 지질공학
• /
• 제26권4호
• /
• pp.461-471
• /
• 2016

수주지열정의 지하수 유동 및 지열 이동 현상을 모사할 수 있는 수치 모델링 기법을 제시하기 위하여 열-수리학적 모델인 TOUGH2를 기반으로 한 모듈을 개발하였고 두 유형의 수주지열정에서의 냉방 및 난방 운영을 수치 모의하였다. 수치 모델링 결과들은 수주지열정의 냉난방 운영이 주변 지층 및 지하수와의 열교환 및 유체 혼합에 크게 영향을 받으며, 특히 PVC 또는 PE 파이프 내부와 수주지열정 내부에서 서로 반대방향으로 유동하는 유입수와 유출수 간의 열교환이 중요한 요소임을 보여준다. 또한 본 연구에서 개발 및 제안된 수치 모델링 기법이 수주지열정의 내부 구조 및 유체 유동과 열 이동 현상을 합리적으로 모사할 수 있음을 보여준다. 이러한 수치 모델링 기법은 수주지열정의 설계, 시공 및 운영 단계에서 열교환 성능을 정량적으로 분석할 때 유용하게 활용될 수 있을 것으로 기대된다.

• Geomechanics and Engineering
• /
• 제16권4호
• /
• pp.363-373
• /
• 2018

The evaluation of Thermo-Hydro-Mechanical (THM) coupling behavior is important for the development of underground space for various purposes. For a high-level radioactive waste repository excavated in a deep underground rock mass, the accurate prediction of the complex THM behavior is essential for the long-term safety and stability assessment. In order to develop reliable THM analysis techniques effectively, an international cooperation project, Development of Coupled models and their Validation against Experiments (DECOVALEX), was carried out. In DECOVALEX-2015 Task B2, the in situ THM experiment that was conducted at Horonobe Underground Research Laboratory(URL) by Japan Atomic Energy Agency (JAEA), was modeled by the research teams from the participating countries. In this study, a THM coupling technique that combined TOUGH2 and FLAC3D was developed and applied to the THM analysis for the in situ experiment, in which rock, buffer, backfill, sand, and heater were installed. With the assistance of an artificial neural network, the boundary conditions for the experiment could be adequately implemented in the modeling. The thermal, hydraulic, and mechanical results from the modeling were compared with the measurements from the in situ THM experiment. The predicted buffer temperature from the THM modelling was about $10^{\circ}C$ higher than measurement near by the overpack. At the other locations far from the overpack, modelling predicted slightly lower temperature than measurement. Even though the magnitude of pressure from the modeling was different from the measurements, the general trends of the variation with time were found to be similar.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제22권3호
• /
• pp.42-49
• /
• 2017

TOUGH2 was used to simulate the migration of $CO_2$ injected into a sandy aquifer. A series of numerical simulations was performed to investigate the effects of a low-permeability layer (LPL) embedded in the aquifer on the injection rate and the pressure distribution of $CO_2$. The results show that the size and location of the LPL greatly affected the spread of $CO_2$. The pressure difference between two points in the aquifer, one each below and above the LPL, increased as the size of the LPL increased, showing a critical value at 200 m, above which the size effect was diminished. The location of the LPL with respect to the injection well also affected the migration of $CO_2$. When the injection well was at the center of the LPL, the injection rate of $CO_2$ decreased by 5.0% compared to the case with no LPL. However, when the injection well was at the edge of the LPL, the injection rate was decreased by only 1.6%. The vertical distance between the injection point and the LPL also affected the injection rate. The closer the LPL was to the injection point, the lower the injection rate was, by up to 8.3%. Conclusively, in planning geologic storage of $CO_2$, the optimal location of the injection well should be determined considering the distribution of the LPL in the aquifer.

• 터널과지하공간
• /
• 제33권5호
• /
• pp.399-413
• /
• 2023

본 연구에서는 TOUGH-FLAC 연동해석 기법을 사용하여 장기분지 유망 저장소 부지에서의 CO₂ 주입을 수치적으로 모델링하고, 주변 지층과 단층에 야기되는 수리-역학적 안정성을 검토하였다. 대상부지의 현장조사와 3D 지질모델을 기반으로 3년간 32,850톤의 CO₂ 주입을 가정하여 시뮬레이션을 수행하였다. 저류층의 투수계수에 따른 CO₂ 플럼의 유동 경로를 분석한 결과, 투수계수를 10^-14 m²로 가정한 경우 방사형 유동을 보이며 2년 9개월 경과 시 단층에 도달하는 것으로 나타났다. 투수계수를 10^-13 m²로 가정한 경우에는 저류층을 따라 서쪽 방향의 유동이 우세하게 나타났으며 단층 방향으로의 유동은 거의 발생하지 않았다. 주입공 주변의 압력 변화는 3년 동안 0.6 MPa 이하로, 지층 및 단층의 역학적 안정성에 미치는 영향은 미미한 것으로 분석되었다.

• 지질공학
• /
• 제26권1호
• /
• pp.87-99
• /
• 2016

본 연구에서는 단열을 포함하는 불포화 암반에서의 지하수 유동 흐름을 예측하기 위한 2차원 수치 모델링을 수행하였다. 특히, 불연속 단열망 모델링(Dscrete Fractured Network, DFN)을 통하여 도출된 투수계수의 분포를 나타내는 k-field를 모델링 입력변수로 적용하였다. 투수계수 차이는 불포화대에서 지하수의 이동속도 차이를 야기시키는 중요한 인자로 적용되었다. 연구지역의 지표 토양층으로부터 지하 대수층까지의 불포화대 깊이를 적용하여 지하수 유동 모델링의 초기조건을 실제와 유사하게 설정하였다. 강우의 지표 침투율은 인공구조물과 자연 토양층의 침투율 차이를 적용하여 실제 불포화 암반을 거쳐서 포화대까지 지하수의 이동흐름을 해석하고 시각화하였다. 특히, 오염물질의 이동 시작점이 될 수 있는 인공구조물의 하부에 모니터링 지점을 설정하여 실제 오염물질원이 지하수에 용해되어 불포화대를 이동할 때의 경로를 예측 하였으며 중력방향인 직하부로 이동하는 것을 확인되었다.