이질적인 다공성 매체에서 변형탄성 및 투과성 변화를 위한 강화 갈킨 방법

초록

본 논문은 이질적인 다공성 매체의 유체 흐름과 고체 변형을 동시에 고려하는 Biot 방정식을 풀기 위해 강화 갈킨(Enriched Galerkin, EG) 유한요소법을 제안한다. EG는 연속 갈킨(CG)과 불연속 갈킨(DG)의 장점을 결합해, 압력의 국부 질량 보존과 계면에서의 진동 억제를 달성하면서도 자유도는 DG보다 2~3배 적게 요구한다. 저자는 EG의 블록 구조 구현, 투과성 변화 모델링, 그리고 다양한 1D·2D·3D 테스트와 실제 북해 유전 사례를 통해 CG·DG와의 성능을 비교·검증한다.

상세 분석

이 논문은 다공성 매체에서 발생하는 변형탄성 현상과 투과성 변화를 동시에 모델링하는 데 있어 기존 연속 갈킨(CG)과 불연속 갈킨(DG) 방법이 갖는 한계를 명확히 짚고 있다. CG는 자유도가 적고 구현이 간단하지만, 고투과성 대비 저투과성 구역 사이의 계면에서 압력 진동과 질량 보존 실패가 빈번히 발생한다. 반면 DG는 내부 페널티와 면적 기반 플럭스 계산을 통해 질량 보존을 보장하지만, 각 요소마다 상수 자유도가 추가돼 자유도가 크게 늘어나고 계산 비용이 급증한다.

강화 갈킨(EG)은 이러한 문제를 해결하기 위해 CG의 연속 다항식 공간에 각 요소 중심에 상수 함수를 추가하는 방식으로 설계되었다. 즉, 압력 공간을 P_EG^k = P_CG^k + P_DG^0 로 정의함으로써, 요소 내부에서는 연속성을 유지하면서도 면 내부에서는 불연속적인 상수 자유도를 제공한다. 이 구조는 DG와 동일한 내부 페널티 형태(대칭 내부 페널티)를 사용하지만, 상수 자유도만 추가되므로 전체 자유도는 DG의 약 1/2~1/3 수준에 머문다.

수학적으로는 Biot 방정식의 선형 모멘텀 방정식은 CG 기반 변위 공간 U_CG^k에 대해 전통적인 약식 형태로 이산화하고, 질량 보존 방정식은 EG 기반 압력 공간 P_EG^k에 대해 대칭 내부 페널티와 평균·점프 연산자를 활용한다. 특히, 투과성 텐서 κ는 변형에 따라 k_m = k_m0·(1+ε_v φ)/(1+ε_v) 형태로 업데이트되며, 이는 변형에 따른 투과성 변화 모델을 자연스럽게 통합한다. 시간 전진은 완전 암시적 1차 스키마(Backward Euler)로 처리해 안정성을 확보한다.

구현 측면에서는 FEniCS와 multiphenics 프레임워크를 활용해 블록 구조를 정의하고, 각 물리량(변위·압력)을 별도 블록으로 구성한다. 이를 통해 기존 DG 코드에 최소한의 수정만으로 EG 솔버를 구축할 수 있다. 또한, 오픈소스 형태로 제공함으로써 재현성과 확장성을 높였다.

성능 평가에서는 Terzaghi 1차 압축 해석, 2차원·3차원 이질 매체(구조적 및 무작위) 시뮬레이션, 그리고 실제 북해 유전 모델을 대상으로 CG·DG·EG를 비교하였다. 결과는 다음과 같다. 첫째, EG는 압력과 체적 변형률에서 계면 진동을 거의 보이지 않으며, DG와 동일한 질량 보존 특성을 유지한다. 둘째, 플럭스 근사에서는 CG가 크게 차이 나는 반면, EG와 DG는 거의 일치한다. 셋째, 자유도 면에서 EG는 2D에서 DG 대비 약 2배, 3D에서는 약 3배 적은 자유도를 요구한다. 넷째, 격자 재세분화 실험에서 EG는 매우 거친 메쉬에서도 수렴성을 유지해, 대규모 시뮬레이션에 유리함을 확인했다.

이러한 분석을 통해 EG가 고투과성·저투과성 이질 매체의 복합 변형탄성 문제에 최적의 균형을 제공한다는 결론을 내렸다. 특히, 자유도 절감과 질량 보존, 계면 진동 억제라는 세 가지 핵심 요구사항을 동시에 만족시키는 점이 큰 강점으로 부각된다.