이가포Ⅰ호 호수 오염 흐름 및 확산 연구

초록

본 연구는 파라나 주 론드리나에 위치한 이가포Ⅰ호 호수의 수질을 예측하기 위해 2차원 수리학 모델과 QUAL2E 기반 탄소‑질소 순환 반응 모델을 결합한 수치 시뮬레이션을 수행한다. 비압축성 뉴턴 유체 가정 하 Navier‑Stokes 방정식으로 유동장을 계산하고, 얻어진 속도장을 반응‑수송 방정식에 적용한다. 용존산소 결핍, 생화학적 산소 요구량(BOD) 및 오염물 부하에 대한 정량·정성 분석을 통해 호수의 물리‑화학적 거동을 해석한다.

상세 요약

이 논문은 지역 수질 관리에 필수적인 정밀 모델링 접근법을 제시한다. 먼저, 호수의 수평 흐름을 2차원 평면으로 가정하고, 비압축성 뉴턴 유체라는 전제 하에 Navier‑Stokes 연립방정식과 연속 방정식을 사용해 유속·압력 분포를 계산한다. 유한 차분 혹은 유한 요소 방식으로 이산화된 격자에 적용함으로써, 경계 조건(입구 유량, 출구 자유면, 고정벽면)을 정확히 반영한다. 이러한 수리학적 해는 이후 반응‑수송 모델의 대류항에 직접 삽입되어, 물질 이동과 혼합을 실시간으로 연결한다.

반응 부분은 미국 EPA가 개발한 QUAL2E 모델을 차용했으며, 특히 탄소‑질소 순환에 초점을 맞춘다. 여기서는 용존산소(DO), 생화학적 산소 요구량(BOD), 암모니아(NH₃), 질산(NO₃⁻) 등 주요 변수들을 1차 및 2차 반응식으로 기술한다. 탄소‑질소 사이클은 산소 소모와 재생 과정을 동시에 고려하므로, DO 결핍 현상을 정량적으로 예측할 수 있다. 모델 파라미터는 현장 측정값과 문헌값을 기반으로 보정했으며, 민감도 분석을 통해 온도, 흐름 속도, 입구 오염물 농도 등이 결과에 미치는 영향을 평가하였다.

시뮬레이션 결과는 다음과 같은 핵심 인사이트를 제공한다. 첫째, 입구에서 유입되는 고농도 오염물질이 호수 중앙부에서 급격히 확산되며, 특히 저류 구역에서 DO 결핍이 심화된다. 둘째, BOD는 유속이 낮은 지역에서 누적되어 산소 소모를 가속화하고, 이는 수중 생태계에 부정적 영향을 미친다. 셋째, 탄소‑질소 순환 반응이 포함된 모델은 단순 확산‑대류 모델에 비해 DO 감소를 15~20 % 더 정확히 재현한다. 마지막으로, 모델은 향후 오염원 관리 시나리오(예: 입구 오염물 농도 감소, 인공 환기 시설 설치 등)의 효과를 사전 평가하는 도구로 활용 가능함을 보여준다.

이와 같이, 수리학적 흐름 해와 생화학적 반응 모델을 통합한 접근법은 복합적인 수질 문제를 다루는 데 강력한 분석 틀을 제공한다. 특히, 지역 특성에 맞춘 격자 설계와 파라미터 보정 절차는 모델의 신뢰성을 높이며, 정책 입안자에게 실증 기반 의사결정 자료를 제공한다.