지하 석회암 광산의 풍화·퇴화 연쇄 파괴 수치모델링

초록

본 논문은 지하 석회암 광산에서 물 응축·침투·CO₂ 확산·칼사이트 용해가 연쇄적으로 진행되어 암석의 포화도·공극률이 증가하고 강도가 저하되는 과정을 정량화한다. 수리·화학·기계적 현상을 각각의 지배 방정식으로 기술하고, 이들을 결합한 다중물리 모델을 제시한다. 간단한 3차원 수치 실험을 통해 모델 구현 가능성을 검증하고, 풍화가 진행될수록 압축·인장 강도가 급격히 감소함을 보여준다.

상세 분석

이 연구는 지하 석회암 광산 붕괴 메커니즘을 ‘응축‑침투‑CO₂ 확산‑용해‑강도 저하’라는 5단계 연쇄 반응으로 정의하고, 각 단계에 대한 물리·화학·역학 방정식을 체계적으로 전개한다. 첫 번째 단계인 물 전이에서는 토양 물리학에서 차용한 수두·흡입압(식 1·2)과 다르시 법칙(식 4)을 이용해 비포화 매질 내 물 흐름을 기술하고, 연속 방정식(식 5)으로 질량 보존을 확보한다. 두 번째 단계인 확산‑대류 전이는 용질(주로 CO₂·Ca²⁺)의 이동을 대류 항(v·C)과 확산 항(D·∇C)으로 구분한 전송 방정식(식 6·7)으로 모델링한다. 여기서 반응 속도 R은 칼사이트 용해에 기인한다.

세 번째 단계인 용해 메커니즘은 CO₂가 물에 용해되어 탄산을 형성하고, 이 탄산이 칼사이트와 반응해 Ca²⁺와 HCO₃⁻를 방출하는 화학 평형식(식 8·9·10·11)으로 표현한다. 온도·압력 의존성을 반영한 평형 상수 K₁, K₂, K_C, K_H를 도입함으로써 실제 지하 환경의 변동을 반영한다. 용해 속도는 경험적 반응식(식 12)으로, 포화도와 반응 차수 n_cal, 반응 계수 k_cal에 따라 비선형적으로 변한다.

네 번째 단계인 기계적 열화는 두 가지 경로—포화도 상승에 의한 강도 감소와 공극률 증가에 의한 구조 약화—를 구분한다. 실험 데이터(그림 2·3·4)를 기반으로 압축 강도 R_c와 인장 강도 R_t가 포화도 S_r 및 공극률 n에 대해 선형·지수적으로 감소함을 확인하고, 이를 식 13·14·15에 통합한다. 특히, 전체 강도 저하의 80 %를 포화도 증가, 20 %를 용해에 의한 공극률 증가로 가정함으로써 모델 단순화를 시도한다.

수치 적용 사례에서는 20 m 깊이의 지하 공동을 9 m × 9 m × 3 m 영역으로 제한하고, 초기 조건으로 Ca²⁺ 농도 1.607 × 10⁻⁴ mol/L, CO₂ 농도 거의 0, 포화도 1인 루프 표면을 설정한다. 경계 조건은 루프에서의 완전 포화와 CO₂ 용해, 하부에서는 수직 흐름 차단을 적용한다. 이러한 설정을 통해 시간에 따른 포화도·공극률·강도 변화를 시뮬레이션하고, 결과적으로 몇십 년에서 수백 년에 걸쳐 ‘폰티스’라 불리는 원형 함몰이 발생함을 예측한다.

본 모델의 강점은 수리·화학·역학을 하나의 연속체로 통합해 풍화 진행 속도와 파괴 시점을 정량적으로 추정할 수 있다는 점이다. 그러나 몇 가지 한계도 존재한다. 첫째, 투과성(K)와 확산계수(D)를 고정값으로 가정했는데, 실제로는 포화도·공극률 변화에 따라 가변적이다. 둘째, 강도 저하를 포화도와 용해 효과로 80 %/20 % 비율로 단순 분배했으나, 실험적 검증이 부족하다. 셋째, 열전달·기체 대류에 의한 온도 구배가 모델에 포함되지 않아, 실제 지하 온도 변동이 미치는 영향을 과소평가할 가능성이 있다. 향후 연구에서는 비선형 투과성 모델, 온도·기체 흐름 연계, 현장 관측 데이터를 통한 파라미터 보정 등을 통해 모델 정밀도를 높일 필요가 있다.