초록
본 수치해석연구에서는 국제공동연구프로젝트 DECOVALEX2019의 Task B의 일환으로 PFC3D를 기반으로한 수리역학연계모델을 개발하여 스위스 Mont Terri 지하연구시설에서 수행된 단층의 유체주입으로 인한 슬립시험을 모사하였다. 이를통해, 개발한 PFC3D 수리역학연계모델이 가진 한계점과 향후 보완할 점을 검토하고자 하였다. PFC3D를 기반으로한 3차원 입자결합모델 내 공극-유동통로모델을 생성하였으며 이를 사용하여 Mont Terri Step 2 단층내 유체주입실험을 모사하였다. 모델링결과 단층대를 따라 주입유체의 유동에 의한 단층대의 변형을 확인하였지만, 관측정에서의 시간에 따른 수압변화는 현장측정치와 부분적으로 일치하는 경향을 확인하였다. 현장측정 관측수압은 초기 유체주입 압력증가에 거의 변화를 보이지 않고 주입수압이 최대치에 도달할때쯤 급격한 증가를 보이는반면, 모델링에서는 주입압력이 증가함에 따라 관측수압도 부드럽게 증가하는 경향을 보였다. 이러한 부분적으로 일치하는 결과의 원인으로는 Mont Terri 현장의 단층을 모사하는 방법에 기인하는 것으로 판단하다. PFC3D에서는 단층을 손상대와 코어균열의 조합으로 모사하였고 단층대의 두께가 약 2 m로 주입유체가 단층대를 통해 유동하도록 모사하였기에 현장에서의 주입유체의 단층내 유동보다 그 유동범위가 크게 모사되었다고 판단한다. 또한, 현장단층에서와 같이 단층내부에 존재하는 충진물질로 인해 단층내 수리유동이 제한되어 국부적으로 과잉공급수압이 형성될 수 있는 기재를 모사하지 못한 점 또한 모델링 결과와 현장측정결과가 부분적으로 일치하는 원인일 수 있다. 단층변형의 경우는 모델링결과와 현장측정결과 유사한 수준으로 일치하는 결과를 확인하였다. 수치모델을 변형하여 단층대의 두께를 감소시키고 단층내 충진 물질의 비균질적인분포를 모사할 수 있는 방법론에 대한 후속 연구를 통해 PFC3D 수리역학연계모델의 유체주입으로 인한 단층활성화 연구로의 적용성을 향상시키는 것을 제안하고 한다.
This study presents an application of hydro-mechanical coupled Particle Flow Code 3D (PFC3D) to simulation of fluid injection induced fault slip experiment conducted in Mont Terri Switzerland as a part of a task in an international research project DECOVALEX-2019. We also aimed as identifying the current limitations of the modelling method and issues for further development. A fluid flow algorithm was developed and implemented in a 3D pore-pipe network model in a 3D bonded particle assembly using PFC3D v5, and was applied to Mont Terri Step 2 minor fault activation experiment. The simulated results showed that the injected fluid migrates through the permeable fault zone and induces fault deformation, demonstrating a full hydro-mechanical coupled behavior. The simulated results were, however, partially matching with the field measurement. The simulated pressure build-up at the monitoring location showed linear and progressive increase, whereas the field measurement showed an abrupt increase associated with the fault slip We conclude that such difference between the modelling and the field test is due to the structure of the fault in the model which was represented as a combination of damage zone and core fractures. The modelled fault is likely larger in size than the real fault in Mont Terri site. Therefore, the modelled fault allows several path ways of fluid flow from the injection location to the pressure monitoring location, leading to smooth pressure build-up at the monitoring location while the injection pressure increases, and an early start of pressure decay even before the injection pressure reaches the maximum. We also conclude that the clay filling in the real fault could have acted as a fluid barrier which may have resulted in formation of fluid over-pressurization locally in the fault. Unlike the pressure result, the simulated fault deformations were matching with the field measurements. A better way of modelling a heterogeneous clay-filled fault structure with a narrow zone should be studied further to improve the applicability of the modelling method to fluid injection induced fault activation.