투과성 이질성이 반응성 대류 용해에 미치는 영향

초록

본 연구는 이질적인 투과성 구조가 이중분자 반응 A + B → C에 의해 유도되는 부력‑구동 대류 용해에 미치는 영향을 수치적으로 조사한다. 균일, 수평·수직 층상, 로그정규 분포 투과성 매체를 비교하고, 반응 전·후의 지문 전진, 혼합 길이, 용해 플럭스 등 핵심 변수들을 분석한다. 수직 층상 및 로그정규 매체는 전반적으로 혼합·용해 효율을 높이며, 수평 층상은 흐름을 억제해 가장 낮은 용해 플럭스를 보인다.

상세 분석

이 논문은 다중 스케일 이질성을 포함한 다공성 매체에서 부력‑구동 대류 용해와 동시에 일어나는 A + B → C 반응을 고차원 수치 모델(rhoDarcyFoam 기반)로 구현하였다. 기본 방정식은 질량 보존(∇·q = 0)과 세 물질의 대류‑확산‑반응식(식 2‑4)이며, 투과성 k(x,z)와 밀도 ρ = R_A c_A + (R_B − R_C) c_B + β R_C 로 정의된 라일리 수(R_i)들이 비선형 밀도 구배를 형성한다. 세 가지 반응 시나리오(R1, R2, R3)는 각각 비단조·비단조·반응에 의해 비안정성이 유도되는 경우를 대표한다.

이질성은 (1) 수평 층상(λ_z = 50, λ_x = ∞), (2) 수직 층상(λ_x = 100, λ_z = ∞), (3) 다중 가우시안 로그정규 매체(λ_x, λ_z 다양) 로 구현했으며, 로그정규 변동성 σ²_log k = 1,2,3을 적용해 통계적 강건성을 확보하였다. 각 경우를 10번씩 난수 초기조건으로 반복해 평균값을 분석하였다.

결과적으로, 수직 층상과 로그정규 매체는 고투과성 채널이 대류 손가락을 가늘고 길게 성장시키며, 혼합 길이와 용해 플럭스를 모두 증가시켰다. 특히 λ_x/λ_z < 1(수직 지배)일 때는 손가락이 수직으로 빠르게 전파돼 파동길이가 짧아지고, 반응 전선이 급격히 전진한다. 반대로 λ_x/λ_z > 1(수평 지배)에서는 손가락이 옆으로 퍼지며 병합이 빨라지고, 전선 전진 속도가 감소한다.

수평 층상은 저투과성 층이 장벽 역할을 해 손가락이 두껍고 수직 이동이 억제되며, 전반적인 용해 플럭스가 가장 낮았다. 반응에 따라 밀도 구배가 비단조(예: R1)일 경우, 반응이 흐름을 안정화시켜 전선이 상부에 머무는 현상이 관찰되었다. 반면 R2와 같이 반응 생성물 C가 가장 높은 라일리 수(R_C = 2)를 갖는 경우, 밀도 구배가 단조 증가해 대류가 급격히 촉진되고, 혼합 길이와 용해 플럭스가 최대가 된다. R3는 초기에는 안정하지만 반응이 진행되면서 비단조 밀도 구배가 형성돼 중간 정도의 대류 강도를 보였다.

또한, 로그정규 매체에서 σ²_log k가 커질수록 수직 층상에서는 혼합 길이가 늘어나지만, 수평 층상에서는 오히려 감소한다는 상반된 효과가 나타났다. 이는 고투과성·저투과성 영역의 공간 배치가 대류 손가락의 성장 방향과 병합 메커니즘을 결정하기 때문이다.

전반적으로, 투과성 이질성은 반응‑대류 상호작용을 복합적으로 조절한다. 짧은 수평 상관길이(λ_x < λ_z)는 수직 채널화를 촉진해 혼합·용해 효율을 높이며, 반대로 긴 수평 상관길이는 수평 확산을 강화해 전반적인 전진 속도를 저하시킨다. 이러한 결과는 CO₂‑암수 저장과 같은 실제 지하공정에서 이질성을 고려한 설계·예측의 중요성을 강조한다.