• 지구물리와물리탐사
• /
• 제26권4호
• /
• pp.268-274
• /
• 2023

유한 단층 미끌림 역산에는 지진 변위 측지 자료와 그린 함수 행렬(Green's function matrix)을 주로 사용한다. 그린 함수 행렬은 일반적으로 오카다 모형(Okada, 1985)을 기반으로 한다. 그러나 최근 물리 기반 지진 모델링을 활용하여 그린 함수 행렬을 제작하고 유한 단층 미끌림 역산을 수행하는 연구가 활발하다. 물리 기반 지진 모델링은 다양한 물성(탄성, 점탄성, 탄소성 등)을 고려하여 현실적인 환경에서 지진을 모사할 수 있다는 장점이 있다. 물리 기반 유한요소 소프트웨어 PyLith는 단층을 구성하는 절점을 두 개로 나누어 지진을 모사할 수 있으므로 지진 모사 모델링에 적합하다. 하지만 PyLith는 격자망 생성 기능을 자체 제공하지 않아, 모형 내부에 수십~수백 개의 소단층과 관측점을 설정해야 하는 유한 단층 미끌림 역산 수행에는 어려움이 있다. 본 연구에서는 소단층과 관측점을 포함한 수치 모형을 제작하고, 지진 모사 모델링을 수행하여 그린 함수 행렬을 제작하는 일련의 과정을 연계하여 유한 단층 미끌림 역산의 편리성을 높이기 위해 CPInterface (COMSOL-PyLith Interface)를 개발하였다. CPInterface는 COMSOL의 격자 생성 능력과 PyLith의 지진 모사 능력을 결합하여 그린 함수 행렬을 자동으로 생성할 수 있다. CPInterface는 간단한 변수들로 모형 및 단층 정보를 조절할 수 있고, 지하 탄성 이상체와 GPS 관측점을 자유롭게 배치할 수 있다. 또한, 그린 함수 행렬을 생성하는 복잡한 과정을 간소화하여 더욱 편리하게 유한 단층 미끌림 역산을 할 수 있게 한다.

• 지질공학
• /
• 제32권4호
• /
• pp.643-659
• /
• 2022

최근 지중저장기술(예, 온실가스 심지층 처분, 인공지열저류층 발전 등)이 활발히 수행됨에 따라, 유체 주입과 저장부지 안정성 사이의 역학적 관계에 관한 정량적 이해의 중요성이 인지되고 있다. 지중 유체 주입은 공극압 및 지중응력 교란과 지층의 역학적 불안정성을 야기할 수 있어, 유체 주입에 대한 다공탄성 수치 모형 구축이 요구된다. 본 연구에서는 순차적인 COMSOL-PyLith-COMSOL 유체 주입-유발지진 다공탄성 수치 모사를 수행한다. 유한요소 상용 소프트웨어인 COMSOL을 이용해 단층에 가해지는 쿨롱 파괴 응력(CFS) 변화를 시간에 따라 추적하였고, CFS 변화량이 임계값(예, 0.1 MPa)을 초과할 경우, 모형의 정보(기하구조, 물성 등)를 유한요소 오픈소스 코드인 PyLith로 이동시키는 알고리즘을 구축했다. PyLith는 단층의 미끄러짐을 모사하고, 미끌림에 의한 변위장을 획득한다. 이후 변위장을 COMSOL로 이동시켜 지진에 의한 응력 및 표면 변위를 계산한다. 수치 모사 결과, 주입 기간 중엔 주입정 근거리에서 큰 변화(공극압, CFS 변화 등)를 보였고, 주입 종류 후에는 잔류 응력이 원거리 영역으로 확산하는 양상이 나타났다. 이는 주입 종료 후 지속적인 모니터링의 필요성을 제안한다. 또한, 단층과 주입층 물성(예, 투수계수, Biot-Willis 계수)에 따른 CFS 변화량 비교는 주입정 위치 선정 시 주입층 및 주변 지층에 대한 물성 파악이 중요함을 의미한다. 단층 미끄러짐 양에 따른 표면 변위 및 이암층에 가해지는 편차응력은 다양한 단층 미끌림 시나리오 설정의 필요성을 지시한다.