다중상 반응 광상과 이방성 다공성 매체에서의 반응성 용질 전이 분석

초록

본 연구는 Deng et al. (2013)의 다중상 반응 광상 모델과 라그랑지안 확률론을 결합해 3차원 이방성 다공성 매체에서의 반응성 용질 확산을 정량화한다. 수리전도도(K)와 평형흡착계수(Kd) 사이의 상관관계가 용질 전이에 큰 영향을 미치며, 이방성 비율은 상대적으로 미미한 영향을 보인다. 평균, 분산, 적분 스케일 등 K와 Kd의 통계적 특성을 변화시킨 민감도 분석을 통해 주요 제어 요인을 도출하였다.

상세 분석

이 논문은 물리적·화학적 이질성이 반응성 수송에 미치는 복합 효과를 정량적으로 해석하기 위해 라그랑지안 입자 추적 기반의 확률론적 프레임워크를 채택하였다. 핵심은 Deng et al. (2013)이 제시한 ‘다중상 반응 광상(multimodal reactive mineral facies)’ 개념으로, 서로 다른 광상마다 고유의 K와 Kd 값을 갖는 계층적 구조를 가정한다. 이러한 가정 하에 K와 Kd의 공간적 변동성을 각각 평균(μK, μKd), 변동성(σ²K, σ²Kd), 그리고 적분 스케일(LK, LKd)으로 파라미터화하고, 두 변수 간 상관계수 ρKH를 도입해 공동 변동성을 기술한다. 라그랑지안 접근법은 입자 경로를 확률적 흐름으로 모델링함으로써, 전단 확산(전단계수)과 반응성 확산(반응계수) 사이의 상호작용을 명시적으로 포함한다.

수치 실험에서는 실제 현장 데이터를 기반으로 3차원 이방성 매체를 설정하고, 시간에 따른 반응성 용질의 변동계수(DR(t))를 계산하였다. 결과는 ρKH가 양(양의 상관)일 때 DR이 크게 감소하고, 음(음의 상관)일 때는 DR이 증가함을 보여준다. 이는 고전적인 비반응성 입자 확산과 달리, K와 Kd가 동시에 높은 영역에서는 용질이 빠르게 이동하면서 동시에 강한 흡착에 의해 지연되는 복합 효과가 발생하기 때문이다.

또한, 이방성 비율(α = Kx/Kz)의 변화가 DR에 미치는 영향은 통계적으로 유의미하지 않음이 확인되었다. 이는 입자 이동 경로가 장축 방향으로 지배적이지만, 반응성 확산은 주로 K와 Kd의 공동 변동성에 의해 좌우되기 때문이다.

민감도 분석에서는 μK, σ²K, LK, μKd, σ²Kd, LKd 각각을 독립적으로 변동시켰을 때 DR에 미치는 영향을 정량화하였다. K의 평균이 증가하면 전반적인 흐름 속도가 빨라져 DR이 감소하고, K의 변동성이 커질수록 전단 확산이 확대되어 DR이 증가한다. Kd의 평균이 상승하면 흡착 강도가 강화돼 DR이 크게 감소하지만, Kd의 변동성이 커지면 흡착 강도의 공간적 불균일성이 증가해 DR이 다시 상승한다. 적분 스케일은 장거리 상관성을 반영하는데, LK와 LKd가 모두 커질수록 큰 규모의 연속된 고/저 전도·흡착 구역이 형성돼 DR의 시간적 진화가 느려진다. 이러한 결과는 현장 설계와 관리에서 K와 Kd의 통계적 특성을 정확히 파악하는 것이 반응성 오염물 전이 예측에 핵심임을 시사한다.