• 터널과지하공간
• /
• 제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 수리역학연계모델의 유체주입으로 인한 단층활성화 연구로의 적용성을 향상시키는 것을 제안하고 한다.

• 터널과지하공간
• /
• 제28권5호
• /
• pp.400-425
• /
• 2018

본 논문에서는 국제공동연구인 DECOVALEX-2019 프로젝트 Task B의 연구결과와 현황을 소개하였다. Task B의 주제는 'Fault slip modelling'으로 유체의 주입으로 인해 발생하는 단층의 재활성(미끄러짐, 전단파괴)과 수리역학적 거동을 예측할 수 있는 해석기법을 개발하는 데에 그 목적이 있다. 1단계 연구는 참가팀들이 연구주제에 대해 숙지하고, 벤치마크 모델을 대상으로 단층의 투수특성과 역학적 거동의 상호작용을 모사할 수 있는 해석코드를 개발할 수 있도록 하는 준비 단계의 연구이다. 본 연구에서는 TOUGH-FLAC 연동해석 기법을 사용하여 물 주입으로 인한 단층의 수리역학적 연계거동을 모사하였다. TOUGH2 해석에서는 단층을 Darcy의 법칙과 삼승법칙을 따르는 연속체 요소로 모델링하였으며, FLAC3D 해석에서는 미끄러짐과 개폐가 허용되는 불연속 인터페이스 요소를 통해 모사하였다. 두 가지 수리간극모델에 대하여 수리역학적 커플링 관계식을 수치화하였으며, 연속체 요소(수리모델)와 인터페이스 요소(역학모델)의 거동을 연계할 수 있는 해석기법을 제시하였다. 또한, 단층의 역학적 변형(간극의 변화)으로 인한 수리물성 변화와 기하학적 변화(해석 메쉬의 변형)를 수리해석에 반영할 수 있는 해석기법을 개발하였다. 다양한 압력의 물을 단계적으로 주입하고 이로 인해 유도되는 단층의 탄성거동 및 전단파괴(미끄러짐)에 대해 살펴보았으며, 수리간극의 변화 양상과 원인, 압력 분포와 주입율의 관계 등을 면밀히 검토하였다. 해석 결과, 본 연구에서 개발한 해석기법이 물 주입으로 인한 단층의 미끄러짐 거동을 합리적인 수준에서 재현할 수 있는 것으로 판단할 수 있었다. 본 연구의 해석모델은 Task B에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 연구의 현장시험에 적용하여 타당성을 검증할 예정이다.