Modelling the onset and evolution of immiscible viscous fingering in porous media

기술 요약: 다공성 매질 내 불침투성 점성 핑거링(Viscous Fingering)의 발생 및 진화 모델링

문제 정의
다공성 매질에서의 불침투성 점성 핑거링 시뮬레이션은 석유 공학 및 수리지질학의 치환 과정을 이해하는 데 매우 중요하다. 상당한 이론적, 수치적 연구가 존재함에도 불구하고, 연속체 규모(continuum-scale)의 시뮬레이션을 실험 데이터와 일치시키는 것은 여전히 어려운 과제로 남아 있다. 구체적으로, 기존 모델들(예: Salmo 등)은 높은 점도비(예: $\mu_o/\mu_w = 2000$)에 대해 초기 단계의 핑거 패턴과 생산 곡선은 성공적으로 재현했으나, 실험에서 관찰된 두 가지 핵심적인 후기 단계 현상을 포착하는 데 실패했다: (1) 핑거 뒤쪽에 안정적인 희박 구역(rarefaction zone)이 존재하지 않는 현상, 그리고 (2) 돌파(breakthrough) 이후 상당한 양의 바이패스된 오일(bypassed oil) 영역이 존재하는 현상이다. 본 논문은 E2000 실험(Bentheimer 사암 슬래브 내에서 물이 2000 mPa·s 오일을 치환하는 실험)에서 관찰된 이러한 특정 특징들을 재현하기 위해 필요한 물리적 메커니즘과 모델링 요구 사항을 다룬다.

방법론
저자들은 상용 블랙오일 시뮬레이터인 IMEX를 사용하여 물리 가이드형 수치 모델링 접근법을 개발하였다. 이 방법론은 세 가지 주요 구성 요소를 통합한다:

1. 선형 안정성 분석 (LSA): 불안정성의 발생을 특성화하는 데 사용된다. 저자들은 초기 단계의 핑거에서 지배적인 파수(wavenumber, $\hat{\nu}$)를 결정하기 위해 E2000 실험의 X-선 이미지를 분석하였다. 이들은 차폐(shielding) 및 병합(merging)에 의한 비선형 효과로 인해 파수가 감소하는 것을 고려하여, 가장 불안정한 파수($\nu_{max}$)가 $\hat{\nu}$보다 유의미하게 크게 되도록 상대 투과도(RP) 및 모세관 압력($P_c$) 함수를 선택하였다.
2. 유동 함수 선택: 두 세트의 RP 함수(RP1, RP2)와 두 개의 $P_c$ 함수(Pc1, Pc2)를 테스트하였다. 선택된 함수(RP2-Pc2)는 실험적 생산 데이터(물 컷 및 오일 회수율)와 일관성을 유지하면서 불안정성을 극대화하도록 선택되었다. 특히 $P_c$ 함수는 미분값과 크기의 비율($|J'_c|/|J_c|$)이 낮은 약한 유성(weakly oil-wet) 시스템을 나타내도록 선택되었으며, 이는 모세관-불균질성 효과를 강화한다.
3. 다중 스케일 투과도 불균질성: 실험적 포화도를 재현하기 위해 저자들은 다중 스케일 투과도 장(field)을 구축하였다. 이 장은 다음을 결합한다:
   • 소규모 불균질성: 핑거 뒤쪽의 안정적인 희박 구역을 파괴하고 팁 스플리팅(tip splitting)을 유도하기 위한 무작위 상관 관계 장(randomly correlated fields).
   • 대규모 불균질성: 핑거 발달을 유도하고 우선 경로를 만들기 위해 실험적 X-선 스캔에서 유도된 "이미지 가이드(image guide)"를 통해 조건화됨.
   • 모세관 불균질성: $P_c$ 함수를 투과도에 의존하게 함($P_c \propto \tau^{-1/2}$), 이는 치환되는 상(phase)이 습윤(wetting)일 때 저투과도 영역에서 침입 상(invading phase)에 대한 장벽으로 작용하는 모세관-불균질성 항을 도입한다.

주요 결과
제안된 모델은 다음과 같은 여러 지표에서 E2000 실험 데이터를 성공적으로 일치시켰다:

• 핑거 스케일: $\nu_{max} \approx 6 \text{ cm}^{-1}$ (실험적 지배 파수인 $\approx 0.6 \text{ cm}^{-1}$의 약 10배)가 되도록 유동 함수를 선택함으로써, 비선형 병합 및 차폐를 고려한 후의 정확한 핑거 폭을 재현하였다.
• 희박 구역의 부재: 소규모 투과도 불균질성을 포함함으로써 균질한 시뮬레이션에서 흔히 나타나는 안정적인 후방 구역을 성공적으로 파괴하였고, 핑거가 입구까지 연장되도록 하여 실험과 일치시켰다.
• 바이패스된 오일: 약한 유성 $P_c$ 함수와 대규모 투과도 불균질성(완만한 대비 $K_{max}/K_{min} \approx 3.5$)의 결합을 통해, 후기 단계(2+ PVI)에서 관찰되는 바이패스된 오일 영역을 성공적으로 재현하였다.
• 생산 데이터: 압력 강하, 유출수 물 컷(effluent water cut), 오일 회수율 곡선을 정확하게 일치시켰다.
• 포화도 분포: 고투과도 영역에서는 높은 포화도를, 저투과도 영역에서는 낮은 포화도(트랩된 오일)를 보이는 포화도 히스토그램의 확장을 포착하였으며, 이는 고정되거나 0인 모세관 압력을 가진 모델에서는 나타나지 않는 특징이다.

의의 및 주장
본 논문은 단순한 히스토리 매칭을 넘어 근본적인 물리적 메커즘을 포착하기 위한 물리 가이드형 제약 조건을 제공한다고 주장한다. 저자들은 다음을 강조한다:

• 가이드로서의 LSA: 선형 안정성 분석은 유동 함수를 선택하기 위한 필수적인 가이드이다. 초기 가장 불안정한 파장은 비선형 진화를 고려할 때 원하는 최종 핑거 폭보다 몇 배 더 작아야 한다.
• 불균질성의 역할: 소규모 불균질성은 희박 구역을 파괴하는 데 필요하며, 대규모 불균질성은 (비록 대비가 완만하더라도) 특히 유성 시스템에서의 모세관-불균질성 효과와 결래되어 바이패스된 오일 영역을 생성하는 데 필요하다.
• 모세한 불균질성 효과: 본 연구는 현장 규모 모델에서 종종 간과되는 모세관 불균질성이 저투과도 영역에서 습윤 상(오일)을 가두는 데 결정적인 역할을 하며, 이는 상대적으로 균질한 매질에서도 치환 효율에 상당한 영향을 미친다는 것을 입증한다.

저자들은 이와 같이 복잡한 시스템에 대한 유일한 예측 모델을 찾는 것은 여전히 어렵지만, 그들의 방법론이 유동 함수, 다중 스케일 불균질성, 그리고 모세관력 사이의 상호작용을 명시적으로 다룸으로써 높은 이동도 대비를 가진 다상 치환을 매칭하는 강력한 템플릿을 제공한다고 결론짓는다.