• 터널과지하공간
• /
• 제29권3호
• /
• pp.197-213
• /
• 2019

본 연구에서는 TOUGH-FLAC 연동해석기법을 이용하여 Mont Terri 지하연구시설에서 수행된 단층 내 물 주입시험을 수치적으로 모델링하고, 단층의 재활성과 수리역학적 거동 특성을 살펴보았다. TOUGH2 해석에서는 단층을 Darcy의 법칙과 삼승법칙(Cubic law)을 따르는 연속체 요소로 모델링하였으며, FLAC3D 해석에서는 미끄러짐과 개폐가 허용되는 불연속 인터페이스 요소를 통해 모사하였다. 현장에서 획득한 단층의 균열개방압력(fracture opening pressure), 주입율, 모니터링 압력, 변위 곡선 등을 바탕으로, 단층의 탄성적 변형과 파괴에 의한 수직팽창 특성을 반영할 수 있는 수리간극모델과 수리역학 커플링 관계를 해석모델에 반영하였다. 한편, 현지응력 조건, 단층의 강도 및 변형 특성에 따른 파라미터 해석을 실시하여 각 입력변수가 해석 결과에 미치는 영향을 분석하였으며, 이를 통해 현장시험 결과를 가장 잘 재현할 수 있는 파라미터 조합을 선정하였다. 해석 결과, 균열개방압력에서 단층의 주입율과 모니터링 압력이 크게 증가하는 현상을 합리적으로 재현할 수 있었다. 하지만, 동일한 입력 변수 조건에서 단층의 전단변위와 파괴영역의 범위는 현장시험 결과에 비해 과대평가되는 결과를 보였다. 이는 해석모델에서는 고압의 주입조건에서 단층의 지속적인 전단파괴가 유도되는 반면, 현장에서는 수리간극의 변화가 전단 미끄러짐보다는 인장력에 의한 단층면의 개방(tensile opening)에 크게 의존하는 것으로 추정되기 때문이다.

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

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

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

본 수치해석연구에서는 국제공동연구프로젝트 DECOVALEX2019의 Task B의 일환으로 PFC3D를 기반으로한 수리역학연계모델을 개발하여 스위스 Mont Terri 지하연구시설에서 수행된 단층의 유체주입으로 인한 슬립시험을 모사하였다. 이를통해, 개발한 PFC3D 수리역학연계모델이 가진 한계점과 향후 보완할 점을 검토하고자 하였다. PFC3D를 기반으로한 3차원 입자결합모델 내 공극-유동통로모델을 생성하였으며 이를 사용하여 Mont Terri Step 2 단층내 유체주입실험을 모사하였다. 모델링결과 단층대를 따라 주입유체의 유동에 의한 단층대의 변형을 확인하였지만, 관측정에서의 시간에 따른 수압변화는 현장측정치와 부분적으로 일치하는 경향을 확인하였다. 현장측정 관측수압은 초기 유체주입 압력증가에 거의 변화를 보이지 않고 주입수압이 최대치에 도달할때쯤 급격한 증가를 보이는반면, 모델링에서는 주입압력이 증가함에 따라 관측수압도 부드럽게 증가하는 경향을 보였다. 이러한 부분적으로 일치하는 결과의 원인으로는 Mont Terri 현장의 단층을 모사하는 방법에 기인하는 것으로 판단하다. PFC3D에서는 단층을 손상대와 코어균열의 조합으로 모사하였고 단층대의 두께가 약 2 m로 주입유체가 단층대를 통해 유동하도록 모사하였기에 현장에서의 주입유체의 단층내 유동보다 그 유동범위가 크게 모사되었다고 판단한다. 또한, 현장단층에서와 같이 단층내부에 존재하는 충진물질로 인해 단층내 수리유동이 제한되어 국부적으로 과잉공급수압이 형성될 수 있는 기재를 모사하지 못한 점 또한 모델링 결과와 현장측정결과가 부분적으로 일치하는 원인일 수 있다. 단층변형의 경우는 모델링결과와 현장측정결과 유사한 수준으로 일치하는 결과를 확인하였다. 수치모델을 변형하여 단층대의 두께를 감소시키고 단층내 충진 물질의 비균질적인분포를 모사할 수 있는 방법론에 대한 후속 연구를 통해 PFC3D 수리역학연계모델의 유체주입으로 인한 단층활성화 연구로의 적용성을 향상시키는 것을 제안하고 한다.

• 터널과지하공간
• /
• 제32권6호
• /
• pp.518-529
• /
• 2022

DECOVALEX-2023 Task C에서는 6개국 9개 참여 기관들이 스위스 Mont Terri 지하처분연구시설에 서 수행된 FE 실험을 대상으로 열-수리-역학 복합거동 모사를 위한 해석코드 개발을 수행하고 있다. 현재 현장시험 결과와 비교 분석을 위한 Step 1이 진행되고 있으며, 본 연구진은 OGS-FLAC 해석 시뮬레이터를 활용하여 일련의 해석을 진행하였다. 해석 결과 히터 가열에 따른 온도 상승이 잘 구현되었고, 상 변화에 따른 완충재 내 포화도 변화를 관측할 수 있었다. 반면 완충재 흡입력의 과대평가로 완충재 내 상대습도, 온도 변화 및 Opalinus 점토암 내 압력 변화가 현장 결과와 다소간 차이를 나타내는 것을 확인하였다. 이를 통해 완충재 흡입력이 처분시스템 해석 시 유동 해석 결과에 지배적인 영향을 미침을 확인할수 있었으며, 향후 지보재 및 초기 수압 모사 개선을 통해 향상된 결과를 도출하고자 한다. 또한, Opalinus 점토암의 열, 수리, 역학적 이방성이 잘 구현되었으며 해석 결과를 통해 OGS-FLAC 시뮬레이터의 처분시스템 해석 적용성을 확인하였다.

• 터널과지하공간
• /
• 제31권6호
• /
• pp.610-622
• /
• 2021

DECOVALEX 프로젝트는 고준위방사성폐기물 처분 시스템에서 발생하는 복잡한 열-수리-역학-화학적(THMC) 복합거동에 대해 수치해석을 통해 보다 더 깊이 이해하기 위해 수행되고 있는 대표적 국제공동 연구이다. 현재 DECOVALEX-2023이 7개 task를 기반으로 진행 중이며, 이 중 Task C는 Mont-Terri 지하연구시설에서 수행된 실규모 정치(FE) 시험을 대상으로 처분시스템 내 THM 복합거동을 모사하는 것을 목표로 하고 있다. 본 연구에서는 자체개발 해석 코드인 OGS-FLAC을 활용하여 수치해석 연구를 수행하였다. 수치모델에서는 FE 시험과 같이 일정 출력의 히터를 수평으로 위치시켰으며, 주어진 계측지점에서 압력 분포, 온도 변화, 역학적 변형을 계측하였다. 완충재 내부로 유입되는 유체 흐름은 완충재의 흡입력으로 인해 발생하였으며, 주변 영역에서는 열 팽창 및 열 압력이 지배적으로 작용함을 확인하였다. 해석 결과는 향후 타 참여 그룹 및 실험 결과와 비교 검증을 수행할 계획이다.