Mathematical Modeling of Salt Precipitation and Multi-Phase Flow in High Enthalpy Fractured Geothermal Systems

본 논문은 고엔탈피 균열형 지열 저류층에서의 비등온 다상 유동 및 암염 침전을 시뮬레이션하기 위해, 견고한 주요 변수 정식화와 이산 균열-매질 접근법을 활용하여 투과도 손상과 운영상의 과제를 정확하게 예측하는, PorePy 프레임워크 내에 구현된 새로운 오픈 소스 조성 유동 모델을 제시한다.

핵심

고온의 지하 저류층을 뜨거운 소금물을 가득 채운 거대한 천연 압력솥이라고 상상해 보십시오. 이것이 지열 에너지의 원천입니다. 이 에너지를 얻기 위해 엔지니어들은 우물을 파고 차가운 물을 주입하여 뜨거운 물을 다시 밀어 올립니다. 하지만 이 과정은 물속에 포함된 '소금' 때문에 까다롭습니다.

소금(암염)을 뜨거운 차에 들어있는 설탕이라고 생각해 보세요. 차를 식히거나 물을 증발시키면, 설탕은 더 이상 녹아있지 못하고 다시 고체 결정으로 변하기 시작합니다. 지열 우물에서도 이와 같은 현상이 발생하는데, 주입정 근처에서 물이 식거나 생산정 근처에서 물이 끓어 증발할 때 그렇습니다. 그 결과, 고체 소금 결정이 형성되어 물이 흐르는 암석의 미세한 구멍과 균열(틈새)을 막아버립니다. 이는 마치 설탕 결정이 고속도로를 막아버려 발생하는 교통 체증과 같습니다.

이 논문은 복잡하게 갈라진 암석 속에서 이러한 "설탕 정체"가 정확히 어디서, 어떻게 발생하는지 예측하도록 설계된 새로운 컴퓨터 시뮬레이션 도구를 소개합니다.

다음은 이 도구가 어떻게 작동하는지, 그리고 무엇을 발견했는지에 대한 내용을 쉬운 비유를 사용하여 정리한 것입니다.

1. 지도: 균열을 명확하게 보기

지하 암석의 전통적인 지도는 모든 것을 매끄럽게 처리하여 암석을 하나의 단단한 스펀지처럼 취급하곤 합니다. 하지만 실제로는 물이 콘크리트 자체를 통과하기보다는 갈라진 보도블록 사이를 흐르는 것처럼, 균열 네트워크를 통해 흐릅니다.

• 혁신 요소: 이 새로운 모델은 "이산 균열-매질(Discrete Fracture-Matrix)" 접근 방식을 사용합니다. 이는 균열을 배경에 흐릿하게 섞는 대신, 지도 위에 뚜렷하고 가는 선으로 직접 그리는 것과 같습니다. 이를 통해 컴퓨터는 균열이 어떻게 연결되는지(혹은 연결되지 않는지)를 정확히 파악하고, 소금이 주변 암석이 아닌 특정 균열을 어떻게 막을 수 있는지 볼 수 있습니다.

2. 엔진: 물리학을 위한 "만능 리모컨"

끓는 물, 증기, 그리고 고체 소금을 동시에 시뮬레이션하는 것은 컴퓨터에게 매우 어려운 일입니다. 보통 컴퓨터는 계속해서 "모드"를 전환해야 합니다(예: "좋아, 이제 액체야; 이제 기체야; 이제 고체야"). 이 과정에서 계산이 충돌하거나 멈출 수 있습니다.

• 혁신 요소: 저자들은 "통합된" 시스템을 만들었습니다. 이것은 배터리를 바꾸거나 모드를 변경할 필요 없이 모든 기기를 작동시키는 만능 리모컨과 같습니다. 이 모델은 물이 액체든, 증기든, 혹은 고체 소금으로 변하든 상관없이 동일하게 유지되는 세 가지 고정된 "다이얼"(압력, 열에너지, 소금 양)을 사용합니다. 덕분에 시뮬레이션은 상태 변화 사이의 혼란스러운 전환을 훨씬 더 부드럽고 안정적으로 처리할 수 있습니다.

3. 속도의 비결: "컨닝 페이퍼"

고온에서 소금물이 겪는 정확한 물리학을 계산하려면 컴퓨터가 복잡한 수학 문제를 반복해서 풀어야 하며, 이는 매우 느립니다.

• 혁신 요소: 연구팀은 미리 계산된 "컨닝 페이퍼"(룩업 테이블)를 만들었습니다. 시뮬레이션을 시작하기 전, 다양한 조건 하에서 소금이 어떻게 행동하는지에 대한 가능한 모든 결과를 미리 계산하여 저장해 두었습니다. 시뮬레이션 중에 컴퓨터는 매번 어려운 수학 문제를 푸는 대신, 이 시트에서 답을 찾아보기만 하면 됩니다. 이 방식은 정확도를 유지하면서도 시뮬레이션 속도를 훨씬 빠르게 만듭니다.

4. 막힘 효과: "공극 크기의 축소"

소금 결정이 형성되면 공간을 차지하게 됩니다.

• 혁신 요소: 모델은 소금이 쌓임에 따라 "파이프"(공극률 및 투과율)를 자동으로 축소시킵니다. 모델은 "만약 소금이 구멍의 10%를 채운다면, 파이프는 상당히 좁아진다"라는 규칙(코죄니-카르만 법칙)을 사용합니다. 이를 통해 소금 정체가 심해짐에 따라 흐름이 어떻게 느려지거나 완전히 멈추는지 예측할 수 있습니다.

시뮬레이션 결과

연구팀은 이 도구를 두 가지 주요 시나리오에서 테스트했습니다.

시나리오 A: 끊어진 고속도로 (연결되지 않은 균열)

• 설정: 균열들이 서로 연결되지 않아 물이 균열 사이의 단단한 암석을 뚫고 지나가야 하는 저류층을 가정합니다.
• 결과: 차가운 물을 주입하자, 생산정 근처의 뜨거운 물이 급격히 끓어올랐습니다. 이로 인해 암석 근처에서 소금이 결정화되어 암석을 막았습니다.
• 반전: 만약 물을 더 빠르게 주입하면, 막힘 현상이 훨씬 심해졌고 에너지 출력도 크게 떨어졌습니다. 모델은 "교통 체증"이 단순히 균열 내부뿐만 아니라 생산정 근처의 암석에서도 발생한다는 것을 보여주었습니다.

시나리오 B: 연결된 고속도로 (연결된 균열)

• 설정: 균열이 주입정에서 생산정까지 이어지는 연속적인 고속도로를 형성하고 있는 저류층을 가정합니다.
• 결과: 차가운 물이 균열을 통해 빠르게 이동했습니다. 물이 매우 빠르게 움직이고 차가운 상태를 유지했기 때문에, 오히려 생산정 근처의 소금을 녹여버리는 결과가 나타났습니다!
• 반전: 소금 침전은 우물을 막는 대신, 차가운 물의 경계 지점 바로 옆에서 발생했습니다. 이는 연결된 균열 네트워크가 우물을 막히는 것으로부터 오히려 보호할 수도 있다는 점을 시사합니다.

결론

이 논문은 엔지니어들이 지열 우물 내의 열, 압력, 소금 사이의 복잡한 상호작용을 이해할 수 있도록 돕는 새로운 오픈 소스 소프트웨어 도구를 소개합니다. 균열이 어떻게 연결되는지와 소금이 어떻게 균열을 막는지 정확하게 매핑함으로써, 이 도구는 다음을 예측하는 데 도움을 줄 수 있습니다:

1. 우물이 소금에 의해 어디서 막힐 수 있는지.
2. "파이프"가 막히기 전까지 얼마나 많은 에너지를 안전하게 추출할 수 있는지.
3. 지하 균열의 구조가 생산 과정에 도움이 될지 혹은 해가 될지.

저자들은 자신들의 도구를 기존의 산업 표준 모델과 비교 검증하였으며, 모델이 완벽하게 일치함을 확인했습니다. 이는 이 도구가 고온의 짠물과 갈라진 암석이 존재하는 환경을 시뮬레이션하는 데 있어 신뢰할 수 있는 방법임을 입증합니다.

기술 요약

기술 요약: 다상 파쇄 고엔탈피 시스템 시뮬레이션을 위한 통합 조성 유동 모델

문제 정의
고엔탈피 지열 저류층을 시뮬레이션하는 것은 비등온, 다상, 다성분 유동의 복잡한 결합과 강력한 비선형 열역학, 그리고 파쇄 네트워크의 지배적인 영향으로 인해 상당한 어려움을 수반한다. 이러한 시스템의 핵심적인 운영 문제는 광물 스케일링, 특히 할라이트(NaCl) 침전이다. 생산 중 발생하는 압력 강하(drawdown)는 비등(boiling)을 유도하여 잔류 염수를 농축시키고 용해도를 초과하게 만들며, 냉수 주입은 국부적인 용해도를 낮춘다. 두 메커니즘 모두 공극 공간과 균열 내에 염을 침전시켜 공극률과 투과율을 심각하게 저해하고, 우물 생산성을 감소시키며, 잠재적으로 우물 폐쇄를 유발할 수 있다. 기존의 연속체 모델은 국부적인 침전 및 균열 규모의 투과율 진화를 포착하는 데 실패하는 경우가 많으며, 변수 전환(variable-switching) 공식에 의존하는 전통적인 조성 시뮬레이터는 상 경계 근처에서 수치적 안정성이 떨어지는 문제가 있다.

방법론
저자들은 오픈 소스 PorePy 프레임워크 내에 구현된 새로운 통합 조성 유동 모델을 제시한다. 이 방법론은 세 가지 핵심 구성 요소를 통합한다:

1. 지속 변수 공식화(Persistent-Variable Formulation): 상의 존재 여부에 따라 주요 변수를 전환하는 전통적인 방식과 달리, 이 모델은 모든 상 구성(액체, 증기, 고체 할라이트)에 대해 고정된 주요 변수 집합인 압력($p$), 비엔탈피($h$), 그리고 전체 염 질량 분율($z$)을 사용한다. 이 접근 방식은 수동적인 전환 없이도 상 전이를 자연스럽게 처리하여 수치적 견고성과 안정성을 향상시킨다.
2. 이산 균열-매질(DFM) 접근법: 물리적 도메인은 균열과 그 교차점이 2D 다공성 매질 내에 임베딩된 저차원 매니폴드(1D 및 0D)로 취급되는 혼합 차원 기하 구조를 사용하여 표현된다. 이러한 명시적 표현은 균열의 연결성과 불균질성을 보존하여, 침전에 의한 투과율 감소와 같은 국부적 과정을 해상할 수 있게 한다.
3. 열역학적 폐쇄 및 선형화: 모델은 마그마 조건을 포함하는 H2O–NaCl 열역학에 대한 Driesner 및 Christoph (2007)의 종합적인 상관관계를 채택한다. 계산 비용을 줄이기 위해 모델은 연산자 기반 선형화(OBL) 전략을 사용한다. 열역학적 특성과 그 미분값은 주요 변수 공간의 구조화된 그리드 상에서 사전 계산되며, 시뮬레이션 중에 다선형 보간(multilinear interpolation)을 통해 호출된다. 이는 비용이 많이 드는 실시간 상 분리(flash) 계산을 제거한다.
4. 동적 투과율 및 간극(Aperture) 진화: 모델은 할라이트 포화도에 따라 매질 공극률과 투과율을 동적으로 업데이트하는 Kozeny-Carman 유형의 관계식을 포함한다. 균열의 경우, 간극은 할라이트 포화도에 따른 멱법칙(power-law) 관계에 따라 진화하며, 이는 다시 3차 법칙(cubic law)을 통해 균열 투과율을 업데이트한다.

지배 방정식(질량 및 에너지 균형)은 확산 항에 대한 다점 유속 근사(MPFA)와 대류 항에 대한 업윈딩(upwinding)을 사용하는 셀 중심 유한 체적 기법을 사용하여 이산화된다. 결과적인 비선도 시스템은 백트래킹 Armijo 라인 서치와 적응형 시간 단계(adaptive time stepping)를 갖춘 뉴턴 방법(Newton's method)을 사용하여 해결된다.

주요 결과
모델은 다음과 같은 수치 실험을 통해 검증 및 적용되었다:

• 검증 (1D 벤치마크): 모델은 확립된 폐쇄형 소스 시뮬레이터인 **CSMP++**를 사용하여 1D 염 용해 벤치마크에 대해 검증되었다. 결과는 네 가지 뚜렷한 상 영역(증기+할라이트, 증기+액체, 단상 액체, 액체+할라이트) 및 그 전이 구간에서 강력한 일치를 보였으며, 이를 통해 지속 변수 공식화와 열역학적 구현을 검증하였다.
• 단절된 균열 네트워크 (2D): 고립된 균열이 있는 저류층에서 유동은 매질 지배적이다. 시뮬레이션 결과, 생산정에서의 압력 강하가 비등을 유도하여 주변 매질 및 인접 균열 세그먼트에서 상당한 할라이트 침전을 일으키는 것으로 나타났다. 이 침전은 매질 투과율을 최대 36%까지 감소시켰으며($K/K_0 \approx 0.64$) 유동 경로를 변화시켰다. 주입 속도를 높이면 근처 우물에서의 비등이 심화되어 더 높은 할라이트 포화도($s_{hal} \approx 0.46$)와 더 심각한 투과율 감소($K/K_0 \approx 0.29$)를 초래하여 생산율에 큰 영향을 미쳤다.
• 연결된 균열 네트워크 (2D): 주입정과 생산정을 연결하는 연속적인 균열 체인이 있는 저류층에서, 수송은 균열 지배적이 되었다. 연결된 네트워크는 빠른 압력 통신과 저염도 유체의 생산정으로의 수송을 촉진했다. 결과적으로, 생산정 근처에서 비등이 발생함에도 불구하고, 도착하는 유체는 국부적인 할라이트 포화도를 방지할 만큼 충분히 희석된 상태를 유지했다. 생산정에서의 침전 대신, 냉각 효과로 인해 주입 전방의 매질과 균열 체인을 따라 침전이 재분배되는 현상이 관찰되었다. 이 구성은 단절된 경우에서 나타난 심각한 투과율 저하를 피할 수 있었다.

의의 및 주장
본 논문은 고엔탈피 지열 시스템에서 광물 스케일링 효과를 포함한 복잡한 열 및 질량 수송을 분석하기 위한 견고한 오픈 소스 수치 도구를 제공한다고 주장한다. 본 연구의 주요 의의는 다음과 같다:

• 수치적 안정성: 통합된 지속 변수 공식화는 상 경계 근처에서의 변수 전환과 관련된 수렴 문제를 방지한다.
• 기하학적 충실도: DFM 접근법은 균열의 연결성을 명시적으로 해상하며, 저자들은 이것이 침전의 공간적 분포와 그에 따른 투과율 감소를 제어하는 결정적 요인임을 입증하였다.
• 운영적 통찰: 결과는 균열의 연결성이 운영상의 과제를 결정한다는 점을 강조한다. 즉, 단절된 네트워크는 근처 우물 폐쇄 및 생산성 손실에 취약한 반면, 연결된 네트워크는 생산정에서의 유체 희석을 유지함으로써 이러한 문제를 완화할 수 있다.
• 재현성: 저자들은 도커(Docker) 컨테이너를 통해 완전한 수치 모델과 실행 워크플로우를 제공하여, 지열 커뮤니티 내에서의 재현성과 추가 개발을 용이하게 한다.

저자들은 자신들의 모델이 할라이트 침전 패턴과 그것이 저류층 성능에 미치는 영향을 효과적으로 예측하며, 균열된 지열 저류층에서 유동, 열역학 및 광물 스케일링 사이의 상호작용을 이해하기 위한 가치 있는 도구를 제공한다고 결론짓는다.