Influence of Fluid Rheology on Fluid Flow in a Natural Fracture Network

본 논문은 비뉴턴 유체의 점탄성 특성이 자연 균열 네트워크 내 유동 경로, 점도 변화 및 압력 강하에 미치는 복잡한 영향을 수치 시뮬레이션을 통해 규명하여, 균열 유동 모델링 시 유체 레올로지를 반드시 고려해야 함을 강조합니다.

핵심

이 논문은 지하의 복잡한 갈라진 틈 (단층) 사이를 흐르는 액체의 움직임을 연구한 것입니다. 보통 우리는 물처럼 점성이 일정하게 유지되는 액체만 생각하지만, 실제 지하에서는 고분자 용액, 원유, 혹은 점성이 매우 높은 진흙 같은 액체가 흐릅니다. 이 액체들은 '뉴턴 유체'가 아니라 **'비뉴턴 유체'**라고 불리는데, 흐르는 속도에 따라 점성이 변하는 특이한 성질을 가지고 있습니다.

이 연구는 이 복잡한 액체들이 지하의 거미줄처럼 얽힌 갈라진 틈 (단층 네트워크) 을 어떻게 통과하는지 컴퓨터 시뮬레이션으로 분석했습니다.

이 내용을 일상적인 언어와 비유로 쉽게 설명해 드릴게요.

1. 핵심 비유: "꿀과 물의 차이"와 "교통 체증"

생각해 보세요. 물은 흐르는 속도가 빨라지든 느려지든 항상 똑같이 미끄러집니다. 하지만 꿀이나 치약은 다릅니다.

• 꿀 (비뉴턴 유체): 젓가락으로 천천히 저으면 끈적하게 붙어있다가 (점성 증가), 세게 빠르게 저으면 뚝뚝 떨어집니다 (점성 감소).
• 치약: 튜브를 살짝 짜면 나오지 않지만, 세게 짜면 쏙 나옵니다.

이 연구는 지하의 갈라진 틈을 복잡한 도시의 도로망으로, 액체를 차량으로 비유합니다.

2. 연구의 주요 발견 (세 가지 상황)

이 연구는 액체가 흐르는 속도에 따라 세 가지 다른 현상이 일어난다는 것을 발견했습니다.

① 느릴 때: "고체처럼 굳어버리는 교통 체증" (항복 응력)

액체가 아주 천천히 흐를 때, 특히 치약처럼 짜야 나오는 액체 (항복 응력이 있는 유체) 는 문제가 됩니다.

• 현상: 액체가 흐르려면 일정 힘 (압력) 을 받아야 하는데, 힘이 부족하면 액체가 고체처럼 굳어 버립니다.
• 비유: 도로에 차가 너무 느리게 움직이다 보니, 차들이 서로 밀착되어 아예 움직이지 않는 고체 덩어리가 되어버린 상황입니다.
• 결과: 지하의 갈라진 틈 중 약 65% 가량이 아예 막혀버려 액체가 통할 수 있는 길이 매우 좁아집니다. 마치 도로의 3 분의 2 가 공사 중이라 차가 한 줄로만 지나가야 하는 것과 같습니다.

② 빠를 때: "미끄러워지는 액체와 소용돌이" (전단 박리)

액체가 빠르게 흐르기 시작하면 이야기가 달라집니다.

• 현상: 빠르게 흐르는 액체 (특히 전단 박리 성질이 있는 것) 는 속도가 빨라질수록 점성이 급격히 낮아져 훨씬 미끄러워집니다.
• 비유: 도로 위를 질주하는 스포츠카처럼, 액체가 미끄러워져서 구석구석까지 빠르게 퍼져나갑니다.
• 결과: 오히려 느릴 때보다 더 많은 갈라진 틈으로 액체가 분산되어 흐르게 됩니다. 하지만 속도가 너무 빠르면 교차로 (갈라진 틈이 만나는 곳) 에서 **소용돌이 (와류)**가 생기며 에너지 손실이 발생합니다.

③ 복잡한 네트워크의 특징: "다양한 속도의 혼재"

지하의 갈라진 틈은 크기와 모양이 제각각입니다.

• 현상: 같은 액체라도 좁고 빠른 길에서는 미끄러워지고, 넓고 느린 길에서는 끈적해집니다.
• 비유: 한 도시 안에서 어떤 차는 고속도로를 200km/h 로 달리고, 어떤 차는 골목길에서 10km/h 로 기어가는 혼란스러운 교통 상황이 동시에 발생합니다.
• 결과: 액체의 속도가 한 가지가 아니라 여러 가지로 나뉘어 분포하게 됩니다.

3. 왜 이 연구가 중요한가요?

기존의 연구들은 대부분 물처럼 점성이 변하지 않는 액체만 가정하고 지하를 분석했습니다. 하지만 실제 석유 시추나 지하수 관리에서는 **고분자 용액 (폴리머)**이나 점성 원유를 사용합니다.

• 기존의 오해: "물이 흐르는 것처럼 계산하면 돼."
• 이 연구의 경고: "아닙니다! 액체가 고체처럼 굳어 길을 막거나, 반대로 너무 미끄러져 예상치 못한 곳으로 퍼질 수 있습니다."

4. 결론: 무엇을 배울 수 있을까요?

이 연구는 지하 자원을 개발할 때 (예: 석유를 더 많이 뽑아내거나, 지열 에너지를 얻거나, 오염된 물을 정화할 때) **액체의 '성격' (점성 변화)**을 무시하면 큰 실수를 할 수 있음을 보여줍니다.

• 천천히 주입할 때: 액체가 굳어서 길이 막힐 수 있으니, 압력을 충분히 줘서 굳지 않게 해야 합니다.
• 빠르게 주입할 때: 액체가 미끄러져서 예상치 못한 경로로 퍼지거나, 교차로에서 에너지가 많이 소모될 수 있으니 이를 고려해야 합니다.

요약하자면, **"지하의 복잡한 길에서 액체가 어떻게 움직이는지 알기 위해서는, 액체가 흐르는 속도에 따라 점성이 변하는 '성질'을 반드시 고려해야 한다"**는 것이 이 논문의 핵심 메시지입니다.

기술 요약

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 지하 유체 (석유, 가스, 지열 등) 의 이동은 암석 내 균열 네트워크를 통해 발생하며, 실제 지하 유체는 점탄성, 전단박화 (shear-thinning), 항복응력 (yield stress) 등을 갖는 비뉴턴 유체인 경우가 많습니다. 특히 중유, 폴리머 주입 (EOR), 수압파쇄 등에서 이러한 유체가 널리 사용됩니다.
• 문제점:
  • 기존 균열 흐름 연구의 대부분은 유체를 뉴턴 유체로 가정하거나, 단순한 단일 균열 (single-fracture) 또는 교차점 (intersection) 모델에 국한되어 있었습니다.
  • 복잡한 자연 균열 네트워크 (Discrete Fracture Network, DFN) 에서 유체의 관성 (inertia) 과 비뉴턴 특성이 상호작용하는 메커니즘은 충분히 연구되지 않았습니다.
  • 기존 모델은 균열 내부의 복잡한 유동 구조 (순환, 와류, 국부적 압력 강하) 와 비선형적인 흐름 거동을 정확히 포착하지 못했습니다.
• 연구 목적: 실제 자연 균열 네트워크 구조에서 비뉴턴 유체 (항복응력과 전단박화를 동시에 갖는) 의 흐름을 시뮬레이션하여, 유체 레올로지가 균열 네트워크 내 유동 분포, 속도 분포, 압력 강하에 미치는 영향을 규명하는 것.

2. 방법론 (Methodology)

• 수치 모델:
  • 지배 방정식: 비압축성 Navier-Stokes 방정식을 기반으로 한 RANS (Reynolds-Averaged Navier-Stokes) 접근법을 사용하여 관성 효과가 지배적인 흐름을 모사했습니다.
  • 유체 모델: 항복응력과 전단박화 특성을 모두 포함하는 Herschel-Bulkley-Papanastasiou (HBP) 모델을 사용했습니다. 이는 항복응력 ($\tau_0$) 이 있는 유체가 전단응력 임계값을 넘을 때만 흐르기 시작하고, 흐를 때는 전단율에 따라 점도가 변하는 특성을 정교하게 표현합니다.
  • 격자 및 해석 도구: Ansys Fluent 2023R1 을 사용하였으며, 복잡한 균열 기하구조를 삼각형 유한요소로 이산화했습니다.
• 모델 설정:
  • 단일 교차점 모델: 모델 검증 및 개별 효과 (항복응력 vs 전단박화) 분석을 위해 X 자형 단일 교차점 모델을 사용했습니다.
  • 자연 균열 네트워크 (DFN) 모델: 노르웨이 Devonian Hornelen 분지의 실제 노두 (outcrop) 데이터를 기반으로 한 8m x 8m 규모의 복잡한 균열 네트워크를 사용했습니다. 이 네트워크는 431 개의 교차점 (Y, X, T 자형 등) 과 불규칙한 개구부 (aperture) 를 포함합니다.
• 유체 조건:
  • 실제 EOR 에 사용되는 잔탄검 (Xanthan gum) 수용액의 레올로지 파라미터 (농도별) 를 적용했습니다.
  • 다양한 주입 속도 ($u_0$) 를 시나리오로 설정하여 저속 (점성 지배) 에서 고속 (관성 지배) 까지 흐름 체제를 분석했습니다.

3. 주요 기여 및 발견 (Key Contributions & Results)

A. 유동 패턴 및 점도 분포의 복잡성

• 점도 이질성: 네트워크 내 기하학적 복잡성으로 인해 전단율이 수 개 차수 (orders of magnitude) 에 걸쳐 다양하게 분포합니다. 이로 인해 비뉴턴 유체는 네트워크 내에서 극심한 점도 차이를 보입니다.
  • 저전단 영역: 항복응력 효과로 인해 점도가 매우 높음 (물보다 10,000 배 이상).
  • 고전단 영역: 전단박화 효과로 인해 점도가 급격히 감소.

B. 주입 속도에 따른 유동 거동 변화

1. 저유량 조건 (Low Injection Rates):

   • 항복응력의 지배: 유체가 항복되지 않는 '고체와 같은 영역 (rigid zones)'이 네트워크 단면적의 최대 **약 65%**까지 차지합니다.
   • 유동 국소화 (Localization): 이러한 비유동 영역이 사지 (dead-ends) 에 고정되거나, 연결된 가지 균열을 막아 유동 경로를 차단합니다. 결과적으로 유동은 몇 개의 우세한 경로로만 국한되며, 네트워크 전체의 연결성 (connectivity) 이 크게 저하됩니다.
   • 비선형성: Darcy 법칙을 따르지 않는 Pre-Darcy 영역 ($u \propto J^q, q>1$) 에서 흐름이 발생하며, 이는 항복응력에 의한 저항 때문입니다.

2. 고유량 조건 (High Injection Rates):

   • 전단박화와 관성의 지배: 높은 전단율로 인해 항복응력 영역이 사라지고 (dead-ends 제외), 유체가 전 영역에서 항복됩니다.
   • 관성 지배 흐름: 전단박화로 인해 유효 점도가 낮아지면서 관성 효과가 증폭됩니다. 균열 교차점 (intersections) 에서 **순환 와류 (recirculation zones)**가 형성되고, 이는 국부적인 압력 손실을 유발합니다.
   • 분산된 유동: 흥미롭게도, 단일 균열 연구에서는 전단박화가 유동을 특정 경로로 집중시키는 (channeling) 것으로 알려졌으나, 본 연구의 복잡한 네트워크에서는 전단박화가 오히려 유동을 더 넓은 가지로 분산시키는 효과를 보였습니다. 이는 관성 효과가 증가하여 구불구불한 경로에서의 점성 저항을 극복하게 하기 때문입니다.

C. 속도 분포 및 압력 강하

• 다중 모드 속도 분포 (Multimodal Velocity Distribution): 비뉴턴 유체는 뉴턴 유체 (이중 모드) 와 달리 네트워크 내에서 매우 다양한 속도 분포를 보입니다. 이는 항복되지 않은 영역과 고속 유동 영역이 공존하기 때문입니다.
• 압력 강하 - 유량 관계:
  • 저속: 항복응력으로 인한 Pre-Darcy 비선형성.
  • 고속: 관성 효과 (Forchheimer 식) 와 전단박화의 결합으로 인한 강한 비선형성.
  • 뉴턴 유체에 비해 비뉴턴 유체의 압력 강하 - 유량 관계가 훨씬 비선형적입니다.

4. 의의 및 결론 (Significance & Conclusion)

• 연구의 혁신성:
  • 자연 균열 네트워크 전체에서 비뉴턴 유체 (항복응력 + 전단박화) 의 흐름을 시뮬레이션한 최초의 연구 중 하나입니다.
  • 기존 단순화된 모델 (1D 그래프 기반 등) 이 놓치고 있던 균열 내부의 미세 유동 구조와 관성 효과를 정량적으로 규명했습니다.
• 실무적 함의:
  • 석유 회수 (EOR) 및 수압파쇄: 폴리머 주입 시 항복응력으로 인해 유체가 특정 경로로만 흐르거나 (국소화), 네트워크 일부가 완전히 차단될 수 있음을 경고합니다. 이는 회수율 저하로 이어질 수 있으므로, 주입 설계 시 레올로지 특성을 반드시 고려해야 합니다.
  • 모델링의 필요성: 단순한 Darcy 법칙이나 Cubic law 는 비뉴턴 유체의 복잡한 거동을 설명할 수 없으며, Navier-Stokes 기반의 상세한 시뮬레이션이 필수적입니다.
• 한계 및 향후 과제:
  • 현재 연구는 2D 평면 흐름과 정상 상태 (steady-state) 가정을 기반으로 합니다. 실제 3D 구조, 비정상 와류 (LES/DNS 필요), 온도/압력에 따른 유체 특성 변화, 그리고 고전단율에서의 폴리머 분해 (mechanical degradation) 효과는 향후 연구 과제로 남았습니다.

요약하자면, 이 논문은 복잡한 자연 균열 네트워크에서 비뉴턴 유체의 레올로지 (항복응력과 전단박화) 가 유동 경로를 차단하거나 분산시키는 등 유동 거동을 근본적으로 변화시킨다는 것을 증명하며, 지하 자원 개발 및 지하수 흐름 예측 시 유체 특성을 정밀하게 고려할 필요성을 강조합니다.