이중공극 매질 습윤상 상대 투과성 예측 임계경로 분석의 우수성

초록

본 연구는 이중공극 구조를 가진 토양에서 물 상대 투과성(krw)을 예측하기 위해 임계경로 분석(CPA)을 적용하고, 기존의 연속·병렬 관 모델(vG‑M)과 비교하였다. 두 종류의 경작 방식(전통 경작 CT와 무경작 NT)에서 수집된 13개 시료의 모세관 압력 곡선과 실험적 krw 데이터를 이용해 두 모델을 적용했으며, CPA가 전반적으로 더 정확한 krw 추정치를 제공함을 확인하였다. 특히 구조와 텍스처 영역을 구분하는 교차점의 정확한 설정과 모세관 압력 곡선의 정밀한 측정이 CPA의 성능에 결정적임을 강조한다.

상세 분석

이 논문은 이중공극 매질, 즉 거시적 구조(domain 1)와 미세 텍스처(domain 2)로 구분되는 토양에서 물 흐름을 기술하는 상대 투과성 모델링에 초점을 맞춘다. 전통적인 연속·병렬 관 모델(vG‑M)은 각 공극 집단을 독립적인 관으로 가정하고, 전체 투과성을 두 집단의 기여를 선형적으로 합산하는 방식이다. 그러나 이러한 접근은 공극 간 연결성 및 임계 흐름 경로를 충분히 반영하지 못한다는 한계가 있다.

저자들은 통계물리학에서 차용한 임계경로 분석(CPA)을 도입한다. CPA는 전체 매질을 무작위 네트워크로 간주하고, 물이 흐를 수 있는 최소 저항 경로, 즉 임계 경로를 찾아내어 전체 투과성을 추정한다. 이때 중요한 변수는 모세관 압력 곡선(Pc–S)이며, 이를 통해 각 공극 크기별 포화도와 압력 관계를 정량화한다. 특히, 구조(domain 1)와 텍스처(domain 2)를 구분하는 교차점(crossover point)을 정확히 정의해야 두 영역의 임계 경로가 올바르게 연결된다.

실험 데이터는 두 가지 경작 방식, 전통 경작(CT)과 무경작(NT)에서 각각 13개의 토양 시료를 선정해 500개의 포화도 단계에서 측정한 모세관 압력 곡선과 물 상대 투과성 데이터를 사용하였다. CPA와 vG‑M을 동일한 입력 데이터에 적용한 결과, 평균 절대 오차와 결정계수(R²) 측면에서 CPA가 현저히 우수하였다. 특히, 고포화 구간에서 vG‑M은 과도하게 높은 krw 값을 예측하는 반면, CPA는 실험값에 근접한 곡선을 제공한다. 이는 고포화 구간에서 거시적 구조가 지배적인 흐름 통로가 되며, 임계 경로가 이를 정확히 포착하기 때문이다.

또한, 교차점 위치를 약간씩 변동시켰을 때 CPA의 예측 정확도가 급격히 변하는 것을 확인하였다. 이는 구조와 텍스처 영역의 경계가 물 흐름 메커니즘을 전환하는 ‘스위치’ 역할을 함을 의미한다. 따라서 실험적 Pc–S 곡선의 고해상도 측정과 교차점 추정이 CPA 적용의 전제조건으로 강조된다.

이 연구는 이중공극 매질의 물 흐름을 모델링할 때, 단순한 병렬·연속 관 가정보다 통계적 네트워크 접근이 더 현실적이며, 특히 토양 관리(경작 방식)와 같은 실제 현장 적용에 있어 정확한 투과성 예측이 가능함을 보여준다. 향후 연구에서는 비균질성, 비정상적인 입자 배열, 그리고 다상 흐름(예: 오일·수)까지 확장할 수 있는 CPA 기반 프레임워크 개발이 기대된다.