🔬 condensed matter

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

본 논문은 높은 점도비 조건에서 다공성 매질 내 비혼합성 점성 핑거링(viscous fingering)의 발생과 진화를 정확하게 시뮬레이션하기 위한 물리적 메커니즘과 모델링 요구 사항을 조사하며, 실험적인 핑거 스케일과 포화 패턴을 일치시키기 위해서는 작은 초기 불안정 파장, 후방 안정성을 방해하기 위한 소규모 채널링 효과의 포함, 그리고 우회된 오일을 포착하기 위한 약한 유성(oil-wet) 모세관 압력 함수의 포함이 필요함을 입증한다.

원저자: Paulo L. K. Caetano Chang, Kundan Kumar, Arne Skauge, Kenneth S. Sorbie

게시일 2026-08-24

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Paulo L. K. Caetano Chang, Kundan Kumar, Arne Skauge, Kenneth S. Sorbie

원본 논문은 CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

돌로 만든 스펀지를 상상해 보십시오. 그 안은 걸쭉하고 끈적한 기름으로 가득 차 있습니다. 이제 이 기름을 밀어내기 위해 물을 밀어 넣는다고 상상해 보십시오. 완벽하고 균일한 세상이라면, 물은 매끄럽고 고른 벽처럼 움직이며 빗자루처럼 기름을 앞으로 쓸어 넘길 것입니다. 하지만 자연은 결코 완벽하지 않습니다. 실제 지하의 석유 저장소 세계에서는 유체들의 점도가 서로 다르고, 암석 자체도 결코 완벽하게 균일하지 않습니다. 물처럼 묽은 유체가 기름처럼 훨씬 더 걸쭉한 유체를 밀어내려 할 때, 두 유체의 경계면은 불안정해집니다. 물은 매끄러운 벽이 되는 대신, 들쭉날쭉한 손가락 모양의 경로를 만들며 돌파하고, 그 뒤에는 손을 쓸 수 없는 거대한 기름 주머니들을 남겨둡니다. '점성 핑거링(viscous fingering)'이라고 알려진 이 현상은 석유를 추출하거나 지하수를 관리하려는 엔지니어들에게 큰 골칫거리입니다. 또한, 복잡하고 숨겨진 다공성 암석의 미로 속에서 유체가 어떻게 이동하는지 예측하려는 과학자들에게는 하나의 난제이기도 합니다.

연구팀은 이 퍼즐의 특정 조각을 해결하기 위해 나섰습니다. 그들은 이 '손가락(fingers)'들이 정확히 어떻게 형성되고, 성장하고, 합쳐지는지, 그리고 왜 표준 컴퓨터 모델들이 실험에서 관찰되는 무질서한 현실을 예측하는 데 자주 실패하는지를 이해하고자 했습니다. 그들의 초점은 물을 벤트하이머 사암(Bentheimer sandstone) 판에 주입하여 물보다 2,000배나 더 걸쭉한 기름을 밀어내는 특정 실험에 맞춰져 있었습니다. 이 극단적인 시나리오에서, 이론적으로 물은 암석을 통해 혼란스러운 분사 형태의 손가락들을 만들며 빠르게 통과해야 하며, 이 과정에서 대부분의 기름을 고립된 상태로 남겨두어야 합니다. 그러나 과학자들이 이 과정을 컴퓨터 시뮬레이션으로 실행했을 때, 결과는 너무나 깔끔했습니다. 시뮬레이션 속의 손가락들은 너무 넓게 자랐고, 손가락 바로 뒤에는 실험에서는 전혀 존재하지 않는 기이하고 안정적인 미접촉 구역(stable zone)이 나타나곤 했습니다. 연구자들은 디지털 모델이 실제 암석처럼 작동하게 만들기 위해 어떤 물리적 요소들이 빠져 있는지 알아내야 했습니다.

연구팀은 X선 영상을 사용하여 손가락들이 어떻게 처음 나타나는지 관찰하며 실제 실험의 초기 순간들을 면밀히 살피는 것부터 시작했습니다. 그들은 이 초기 손가락들 사이의 간격을 측정했고, 그것들이 상당히 좁다는 것을 발견했습니다. 하지만 표준 시뮬레이션을 실행했을 때, 손가락들은 훨씬 더 넓었습니다. 세심한 분석을 통해 그들은 직관에 반하는 규칙을 발견했습니다. 시뮬레이션이 올바르고 좁은 손가락을 만들어내려면, 컴퓨터 모델이 최종적인 손가락 크기보다 훨씬 더 작은 규모에서 불안정성을 생성하도록 설정되어야 한다는 것이었습니다. 이는 마치 모델이 아주 작은, 불안정한 물결들로 이루어진 혼란스러운 덩어리에서 시작해야 하는 것과 같습니다. 시뮬레이션이 진행됨에 따라, 이 작은 물결들은 자연스럽게 합쳐지고 서로를 보호하며, 실험에서 보이는 더 크고 안정적인 손가락들로 결합됩니다. 만약 모델이 이미 최종 손가락 크기만한 물결들로부터 시작했다면, 그 결과는 현실보다 너무 매끄럽고 느린 시뮬레이션이 되었을 것입니다.

하지만 손가락의 크기를 맞추는 것은 절반의 성공에 불과했습니다. 실제 실험에서는 손가락 뒤로 평온하고 안정적인 구역이 나타나지 않았으며, 손가형들은 주입 지점까지 길게 뻗어 있었습니다. 반면, 시뮬레이션은 계속해서 이 평온한 구역을 만들어냈습니다. 연구자들은 실험 속의 암석이 육안으로는 균일해 보일지라도, 구조적으로 미세하고 숨겨진 변동성을 가지고 있음이 틀림없다는 사실을 깨달았습니다. 그들은 유체가 통과할 수 있는 능력, 즉 투과율(permeability)의 미세한 변화를 도입했습니다. 이러한 미세한 변동성은 방해 요소로 작용하여, 평온한 구역을 깨뜨리고 손가락이 뒤틀리고 갈라지도록 강제하여 실제 암석에서처럼 작동하게 했습니다. 그러나 그들은 주의를 기울여야 했습니다. 만약 변동성이 너무 강하면 물이 점박이 패턴의 포화 상태를 만들며 갇히게 되는데, 이는 실험에서 관찰되지 않은 현상이었습니다. 핵심은 평온한 구역을 깨뜨릴 만큼은 충분히 강력하면서도, 점박이 형태의 엉망인 상태를 만들지는 않는 '골디락스(Goldilocks)' 수준의 미세한 거칠기를 찾는 것이었습니다.

마지막 퍼즐 조각은 작업이 끝난 후 암석에 남아 있는 갇힌 기름이었습니다. 실험에서는 상당한 양의 기름이 우회되어, 오랜 시간이 지난 후에도 큰 주머니 형태로 뒤처져 있었습니다. 연구자들은 이것이 암석의 습윤성(wettability)과 숨겨진 변동성 사이의 미묘한 상호작용 때문에 발생한다는 것을 발견했습니다. 실험 속의 암석은 약간 '오일 웨트(oil-wet)', 즉 기름이 암석 표면에 달라붙는 성질을 띠고 있었습니다. 이를 미세한 투과율 변동성과 결합하자, 물이 저투과 구역으로부터 밀려나면서 그곳에 기름이 갇히는 상황이 만들어졌습니다. 이러한 특정 유형의 암석 거동을 모델에 포함함으로써, 연구자들은 마침 finally 실험에서 관찰된 우회된 기름의 큰 주머니들을 재현할 수 있었습니다.

결과는 이러한 복잡한 유체 이동을 모델링하기 위한 새로운 가이드라인이 되었습니다. 연구자들은 지표 아래에서 기름과 물이 어떻게 상호작용하는지를 정확하게 예측하려면, 암석이나 유체의 평균값만을 사용해서는 안 된다는 것을 보여주었습니다. 대신, 가장 작은 불안정성이 가장 큰 손가락으로 성장한다는 점, 미세하게 숨겨진 변동성이 흐름이 너무 질서 정연해지는 것을 방지한다는 점, 그리고 유체가 암석 표면에 달라붙는 방식이 겉보기에 균일한 암석 속에서도 기름을 가둘 수 있다는 점을 반드시 고려해야 합니다. 이러한 요소들을 결합함으로써, 팀은 실제 실험과 거의 완벽하게 일치하는 시뮬레이션을 만들어냈으며, 손가락의 혼돈스러운 아름다움, 평온한 후방 구역의 부재, 그리고 뒤처진 끈질긴 기름 주머니들을 모두 포착해 냈습니다. 이 연구는 지각 속에서 유체가 어떻게 이동하는지에 대한 더 명확한 경로를 제시하며, 자원을 관리하고 지하 저장소의 거동을 예측하는 데 더 나은 도구를 제공합니다.

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

문제 정의
다공성 매질에서의 점성 핑거링 시뮬레이션은 석유 공학 및 수문 지질학의 치환 과정을 이해하는 데 매우 중요하다. 특히 혼합 시스템(miscible systems)에 대해서는 광범위한 이론적 및 수치적 연구가 존재하지만, 불혼합(immiscible) 점성 핑거링을 다르시(Darcy) 척도에서 모델링하는 것은 여전히 도전적인 과제이다. Salmo 등의 선행 연구는 무작위 상관 관계를 가진 투과도 장(random correlated permeability fields)과 특정 상대 투과도 함수를 사용하여 파과(breakthrough) 전까지의 핑거 진화를 성공적으로 모델링하였다. 그러나 이러한 모델들은 1회 이상의 공극 부피(pore volume)가 주입된 후의 후기 단계 실험에서 관찰되는 상당한 "바이패스된 오일(bypassed oil)" 영역을 포착하는 데 실패하였으며, 핑거 뒤쪽에 안정적이고 교란되지 않은 "희박 구역(rarefaction zone)"을 생성하는 경향이 있었다. 이는 Skauge 등이 Bentheimer 사암 슬래브를 대상으로 수행한 특정 고점도비 실험( $\mu_o/\mu_w = 2000$ )에서는 나타나지 않는 특징이다. 본 연구가 다루는 핵심 문제는, 비현실적인 불균질성 크기에 의존하지 않으면서 초기 단계의 핑거 스케일과 후기 단계의 포화 패턴(바이패스된 오일 포함)을 모두 재현할 수 있는 물리 기반의 수치 모델을 개발하는 것이다.

방법론
저자들은 상용 블랙 오일 시뮬레이터인 IMEX를 사용하여 E2000 실험(30×30×2.55 cm Bentheimer 슬래브 내에서 2000 mPa·s 오일을 물이 치환하는 실험)에 초점을 맞춘 수치 시뮬레이션 프레임워크를 개발하였다. 방법론은 세 가지 주요 구성 요소의 통합으로 이루어진다:

선형 안정성 분석 (LSA): 저자들은 불안정성의 발생을 특성화하기 위해 LSA를 활용하였다. 그들은 실험 X-선 이미지를 분석하여 핑거의 지배적인 파수(dominant wavenumber, $\hat{\nu} \approx 0.6 \text{ cm}^{-1}$ )를 결정하였다. 결정적으로, 저자들은 비선형 영역(차폐 및 병합 이후)에서 실험적 핑거 스케일과 일치시키기 위해서는 LSA에 의해 예측된 가장 불안정한 파수( $\nu_{max}$ )가 관찰된 지배적 파수보다 훨씬 커야 함(약 10배)을 입증하였다.
유동 함수 선택: 두 가지 세트의 상대 투과도(RP) 및 모세관 압력( $P_c$ ) 함수를 테스트하였다. 선택된 함수(RP2 및 $P_c2$ )는 불안정성을 극대화(높은 $\nu_{max}$ )하면서도 생산 데이터(수분율 및 회수율)와 일관성을 유지하도록 선택되었다. $P_c$ 함수는 미세하게 유층 습윤(oil-wet)인 것으로 모델링되었으며, 모세관 불균질성 효과를 강화하기 위해 미분값과 크기 비율( $|J'_c|/|J_c|$ )이 낮은 Leverett J-함수를 사용하였다.
다중 스케일 투과도 모델링: 실험적 관찰을 재현하기 위해 저자들은 다중 스케일 투과도 장을 구축하였다:
- 소규모 불균질성: 핑거 형성을 촉진하고 핑거 뒤쪽의 안정적인 희박 구역을 파괴하기 위해 짧은 상관 길이를 가진 무작위 가우시안 시뮬레이션을 사용하였다(이는 균질한 시뮬레이션에서 흔히 나타나는 안정적인 꼬리 현상을 방지한다).
- 대규모 불균질성: 관찰된 포화 패턴과 고/저 투과도 구역을 정렬하기 위해 X-선 스캔 데이터에 조건화된 더 큰 규모의 상관 관계 장을 중첩시켜 바이패스된 오일 형성을 유도하였다.
- 모세관 불균질성: 모델은 포화도와 투과도 모두에 의존하는 $P_c$ 함수를 포함하여, 모세관력이 저투과도 영역에 젖음상(유층/오일)을 가두는 것을 허용하였다.

주요 기여 및 결과
본 연구는 E2000 실험의 포화도 장과 생산 데이터를 일치시키는 수치 모델을 성공적으로 구축하였다. 주요 결과는 다음과 같다:

파수 격차 (Wavenumber Disparity): 불안정성 발생 시의 가장 불안정한 파장은 원하는 최종 핑거 폭보다 몇 배 더 작아야 한다. 이를 통해 차폐 및 병합이라는 비선형 과정이 파수를 실험적으로 관찰된 스케일로 감소시킬 수 있다.
희박 구역의 파괴: 실험에서 나타나는 안정적인 후방 구역의 부재는 소규모 투과도 불균질성에 기인한다. 이러한 불균질성이 없다면 수치 시뮬레이션은 일반적으로 핑거 뒤에 안정적인 영역을 보여주지만, 소규모 변동성을 도입하면 이 구역을 교란하고 핑거 끝 분열(tip splitting)을 촉진한다.
오일 바이패스 메커니즘: 모델은 비교적 균질한 매질 내에서의 바이패스된 오일이 약한 유층 습윤 모세관 압력 함수와 중간 규모의 대규모 투과도 불균질성의 결합으로 설명될 수 있음을 보여준다. 모세관 불균질성 항은 장벽 역할을 하여, 침입하는 수상이 치환된 유상(wetting phase, 오일)이 갇혀 있는 저투과도 영역으로 들어가는 것을 방지한다.
정량적 일치: 최종 모델은 돌파 시간, 수분율, 오일 회수율 및 다양한 단계(바이패스된 오일이 존재하는 후기 단계 포함)에서의 공간적 포화도 분포(히스토그램을 통해 검증됨)를 재현하였다.
습윤성에 대한 민감도: 모세관 불균질성이 없거나 수층 습윤(water-wet) $P_c$ 함수를 사용한 시뮬레이션은 바이패스된 오일을 재현하는 데 실패하였으며, 이는 습윤성과 투과도 변화 사이의 특정한 상호작용이 이 현상에 필수적임을 확인해 준다.

의의 및 주장
본 논문은 불혼합 점성 핑거링 역학을 포착하기 위한 "물리 기반 제약 조건"과 모델링 가이드라인을 제공한다고 주장한다. 저자들은 이러한 복잡한 시스템에 대한 유일한 예측 모델을 찾는 것이 여전히 어렵지만, 본 연구의 접근 방식이 높은 이동도 대비(mobility contrast)를 가진 다상 치환을 일치시키는 강력한 템플릿을 제공한다고 주장한다.

본 연구의 의의는 극단적인 불균질성을 끌어들이지 않고도 바이패스된 오일을 실험적 관찰과 일치시킬 수 있는 연속체 규모 시뮬레이션을 구축했다는 점에 있다. 저자들은 필드 규모에서는 무시될 수 있는 모세관 불균질성 효과가 거시적 규모에서 젖음상을 가두는 데 중요한 역할을 하며, 결과적으로 전체 치환 효율에 영향을 미친다고 주장한다. 적당한 투과도 대비를 가진 약한 유층 습윤 시스템이 복잡한 바이패스 패턴을 재현할 수 있음을 입증함으로써, 본 연구는 불안정한 치환을 지배하는 물리적 메커니즘에 대한 정교한 이해를 제공한다.

유사한 논문