다중 스케일 침투 퍼콜레이션을 통한 암석 균열망 연구

초록

본 논문은 이질적인 균열 개구면을 가진 암석 균열망에서 유체 침투 과정을 다중 스케일로 모델링한다. 미세 스케일에서는 균열 내부 흐름이 개구면 변동에 의해 채널화되는 현상을 확인하고, 이를 벡터화된 인공 채널망(ACN)으로 축소한다. ACN을 이용해 다양한 균열 크기와 교차부의 효율적 개구면을 정량화하고, 균열 경사, 흐름 경로 길이, 교차 길이 등의 매개변수를 스케일링 관계로 정립한다. 최종적으로 이러한 관계를 적용해 포화된 이산 균열망에 대한 비혼합 유체 침투 시뮬레이션을 수행한다.

상세 분석

이 연구는 암석 균열망의 복잡한 물리적 특성을 정량화하기 위해 두 단계의 스케일링 전략을 도입한다. 첫 번째 단계는 개별 균열 내부의 미세 흐름을 직접 해석하는 것으로, 균열 표면의 개구면 분포가 랜덤하게 변동함에 따라 유체는 고전적인 포어 흐름이 아니라 제한된 고속 채널을 따라 이동한다는 점을 실험적·수치적으로 입증한다. 이러한 채널화 현상은 전통적인 등방성 평균 개구면 모델이 과소평가하는 전도성을 보정하는 데 핵심적인 역할을 한다. 저자들은 채널 네트워크를 그래프 형태로 추출하고, 각 엣지를 흐름 저항과 연결된 효과적 개구면(효율적 aperture)으로 매핑한다. 이때 채널의 폭, 길이, 곡률 등 기하학적 파라미터가 전도성에 미치는 비선형 영향을 다항식 형태의 스케일링 법칙으로 정리한다.

두 번째 단계에서는 이러한 미세 채널 정보를 인공 채널망(Artificial Channel Network, ACN)이라는 벡터화된 모델로 축소한다. ACN은 실제 균열의 통계적 특성을 보존하면서도 계산 비용을 크게 낮춘다. 다양한 크기의 균열에 대해 다수의 ACN 시뮬레이션을 수행함으로써, 균열 크기(L), 입·출구 교차 길이(L_int), 그리고 균열 경사(θ)와 같은 매크로 파라미터가 효율적 개구면에 미치는 영향을 정량화한다. 특히, 입구와 출구 영역에서의 개구면 비대칭성은 침투 전선의 전파 속도와 포화도 분포에 결정적인 영향을 미친다. 저자들은 이러한 관계를 경험적 상수와 지수 형태의 식으로 정리하고, 이를 통해 임의의 균열망에 대한 전도성 매개변수를 빠르게 추정할 수 있는 알고리즘을 제시한다.

마지막으로, 정립된 스케일링 관계와 ACN 기반 전도성 모델을 이용해 비혼합 유체(예: 물과 비점착성 가스)의 침투 퍼콜레이션을 이산 균열망에 적용한다. 침투 과정은 전통적인 침투 퍼콜레이션(Invasion Percolation, IP) 알고리즘에 비해, 각 균열의 효율적 개구면을 가중치로 사용함으로써 실제 물리적 압력 구배와 유동 저항을 보다 정확히 반영한다. 시뮬레이션 결과는 기존 연구에서 보고된 침투 전선의 비선형 성장, 클러스터 형성, 그리고 임계 침투 압력의 분포와 일치하면서도, 채널화 효과가 포함된 경우 침투 속도가 현저히 증가하고, 특정 경사와 교차 길이를 가진 균열이 침투 경로를 주도한다는 새로운 통찰을 제공한다. 전체적으로 이 논문은 미세-거시 스케일을 연결하는 체계적인 방법론을 제시함으로써, 암석 공학, 지하수 흐름, 그리고 탄산가스 저장 등에서 요구되는 정밀한 다공성 매체 모델링에 중요한 기여를 한다.