이중 광물 암석의 용해가 수리특성에 미치는 영향

초록

본 연구는 격자볼츠만법(LBM) 기반의 공극 규모 수치 모델을 이용해, 두 종류의 광물이 혼합된 파쇄 및 다공성 매질에서 용해 반응이 투수성·공극률에 미치는 변화를 조사하였다. 퀸텟 구조 생성 집합(QSGS)으로 광물 분포를 재구성하고, 다양한 Peclet(Pe)와 Damköhler(Da) 수치에서 시뮬레이션을 수행하였다. 비반응성 광물에 의해 형성된 다공층이 용해 전선 뒤에 남아 흐름 저항을 크게 증가시키며, 공극률은 크게 늘어나도 투수성은 미미하게 변하거나 거의 정체된다. 또한 화학·기계·수리적 구멍폭이 서로 크게 차이나는 이질적 용해 현상이 관찰되었으며, 비반응성 광물 층이 존재할 경우 단일 광물에서 흔히 나타나는 웜홀링 현상이 억제된다.

상세 분석

이 논문은 LBM을 활용한 다중 물리·화학 현상 결합 모델을 구축함으로써, 미세공극 수준에서 용해 반응이 매질의 수리특성에 미치는 메커니즘을 정량적으로 해석한다. 먼저, 퀸텟 구조 생성 집합(QSGS) 기법을 통해 두 광물(반응성·비반응성)의 공간적 분포를 통계적으로 재현했으며, 이는 실제 암석의 이질성을 충분히 반영한다는 점에서 의미가 크다. 시뮬레이션에서는 Peclet 수(Pe)를 통해 대류와 확산 비율을, Damköhler 수(Da)를 통해 화학 반응 속도와 물질 이동 속도의 상대적 중요성을 조절하였다. 높은 Pe·Da 조합에서는 반응 전선이 급격히 전진하면서 용해가 집중적으로 진행되지만, 비반응성 광물이 남긴 다공층이 흐름 경로를 차단한다. 이 다공층은 미세공극이지만 연속적인 통로를 형성해 전단 저항을 크게 증가시키며, 결과적으로 투수성(K)은 공극률(ϕ)의 증가에도 불구하고 포화 상태에 가까운 낮은 값을 유지한다. 반면, 낮은 Pe·Da 조건에서는 확산이 지배적이어서 용해가 균일하게 퍼지며, 다공층의 저항 효과가 상대적으로 감소해 투수성 증가가 더 뚜렷하게 나타난다.

또한, 논문은 화학적 구멍폭(aperture), 기계적 구멍폭, 수리적 구멍폭을 각각 정의하고, 이들 사이에 큰 차이가 존재함을 실험적으로 보여준다. 화학적 구멍폭은 반응 전선이 진행한 거리로 측정되지만, 비반응성 광물에 의해 형성된 잔류 다공층은 실제 유체 흐름을 제한하므로 수리적 구멍폭은 훨씬 작다. 이러한 차이는 현장 실험에서 전통적인 투과성 측정만으로는 용해 진행 상황을 정확히 파악하기 어려움을 시사한다.

마지막으로, 복합 다공성 구조에 대한 추가 시뮬레이션 결과는 비반응성 광물 층이 존재할 경우, 단일 광물에서 흔히 관찰되는 웜홀링(채널화) 현상이 억제된다는 점을 강조한다. 웜홀링은 고투과성 경로가 급격히 확장돼 전체 투수성을 급증시키는 현상이지만, 비반응성 광물의 잔류가 이러한 경로 형성을 방해해 용해가 보다 균일하게 진행되도록 만든다. 따라서, 실제 지하수·석유·탄산염 저장 등에서 다중 광물 시스템을 고려한 모델링이 필요함을 강력히 제시한다.