2.1. 대상유역 하천개황

남강유역은 남강댐을 경계로 남강 상류유역과 하류유역 으로 구분되며 남강댐 하류유역은 남강댐을 기점으로 낙동 강 합류 전까지의 남강 본류 77.6 km에 해당하며 유역 인 구는 343,000명 정도이고 진주시와 진주상평공단 등 중·소 도시 및 공업지역이 입지해 있어 오염물질 발생이 상당히 높다(Fig. 1). 낙동강에 대한 남강 하류유역의 BOD, TP 오 염부하기여율은 각각 44.2%, 23.0%로서 매우 높은 편이다 (^{NRERC, 2009}; ^{Yu et al., 2012}).

현재 남강 하류유역의 환경기초처리시설은 하수종말처리 장이 6개소(총 시설용량 167,000 m³/day), 소규모 마을하수처 리시설이 32개소(총 시설용량 1,792 m³/day)가 있으며, 그 외 에도 폐수종말처리시설 5개소 등이 운영 중이다. Table 1에 서 보는 바와 같이 2007년~2011년까지 최근 5년간 하수종 말처리장 방류수 농도를 조사한 결과, 진주하수처리장 BOD, TP 방류수 농도가 2007년 평균 6.9 mg/L, 0.931 mg/L에서 2011년 각각 평균 3.1 mg/L, 0.380 mg/L로 대부분의 처리장 방류수 농도는 목표수질 달성을 위해 지속적으로 개선되었 음을 알 수 있다.

2.2. 대상유역 분석자료

본 연구에서는 남강 하류유역에서 조사된 2008년~2012년 까지 총 5년간의 모니터링 자료를 분석에 사용하였다. 먼저 수질 자료는 국립환경과학원 낙동강물환경연구소 측정자료 를 이용하였으며, 분석항목은 BOD, COD_Mn, TOC, TN, TP, Chl-a, Water temperature, pH, DO, Conductivity, Ammonia nitrogen, Nitrate nitrogen, Phosphorate, outflow 등 총 14개 항목을 사용하였고, 유량자료는 국토부(환경부) 정암수위국 (남강D)과 거룡강수위국(남강E) 지점의 각 연도별 유출분석 을 통해 산출된 일유량자료를 분석에 사용하였다.

강우량자료는 남강유역에 구축된 기상청 티센망의 티센 계수값을 이용한 티센평균 일강우량을 분석에 사용하였고, 남강댐 방류량 자료는 국가수자원관리종합정보 인터넷홈페 이지(www.wamis.go.kr)에서 수집된 남강댐 본류 일방류량 자료를 분석하여 사용하였다.

2.3. 수질인자간 통계 분석

상관성 분석은 수질인자간 연관성이 있는 수질항목을 파 악하기 위하여, 인자간의 상관관계를 분석하여 상호 연관성 높은 인자와 항목을 추출하는 방법이다. 상관계수의 크기가 0에 가까우면 두 변수사이에 상관성이 거의 없으며, 크기 가 1에 가까우면 두 변수 사이의 상관성이 매우 높다고 해 석한다. 다음으로 상관계수의 부호가 양(Positive)이면 두 변수가 같은 쪽으로 변화하고 음(Negative)이면 두 변수는 반대로 변화한다. 요인분석은 주성분분석법과 최소 고유값 기준(고유값 1이상)으로 요인의 수를 결정하였으며, 각 요 인의 특성을 알고자 할 때 유용한 베리맥스(Varimax) 회전 방식을 적용하였다.

2.4. 모델의 정확도 평가 방법

모델의 정확도를 평가하기 위해서 Table 2에 제시된 통 계값을 사용하였다. 결정계수(Coefficient of determination; R²)는 예측값과 실측값의 선형 회귀관계의 적합도를 결정 하는 계수로써 0과 1사이의 값을 가지며 1에 가까울수록 적합도 및 모형의 정확도가 높다. 평균제곱근오차(Root mean square error; RMSE)는 예측결과의 평균오차를 나타내는 지표로써 그 값이 0에 가까울수록 즉 값이 작을수록 모형의 정확도가 높고 모형의 성능도 우수하다고 평가할 수 있다 (^{Choi. et al., 2008}).

3.1. 남강댐 현황 및 이수기능

남강댐의 주요 시설물은 높이 43 m, 길이 1,126 m의 콘 크리트 표면차수벽 석괴댐인 본댐과 문비 3문을 갖춘 여수 로, 시설용량 14,000 kW의 발전소, 그리고 홍수시 사천만으 로 홍수를 방류하기 위한 높이 31 m, 길이 258 m의 콘크리 트 중력식 제수문 구조물 등으로 구성된다. 남강댐의 일반 제원으로 유역면적은 2,285.0 km², 저수면적 28.2 km², 유역 연평균유입량 64.4 m³/s, 계획홍수위 46.0 EL.m, 상시만수위 41.0 EL.m, 저수위 32 EL.m이다.

남강댐을 취수원으로 광역상수도 사업을 통하여 용수를 공급받는 곳은 통영시, 사천시, 고성군, 진주시, 거제시, 고 성군, 하동군, 남해군 등이다. Table 3과 같이 기본계획공급 량에 따르면 남강댐은 총 573.3 백만톤 중에서 생·공업용수 224.4 백만톤, 농업용수 226.8 백만톤, 하천유지용수 122.1 백만톤을 공급하는 것으로 계획되어 있다.

물공급량(Water yield)은 특정 근원으로부터 획득할 수 있는 물의 양을 말한다. 대개 침투나 증발이 없다고 가정 하면 하천으로부터 최대 공급량은 연평균 유출량과 같게 된다. 이에 따라 남강댐에서 공급가능한 최대 물공급량은 기본계획상으로는 64.4 m³/s이고, Table 4의 최근 5년간 유 입량 분석결과로는 73.5 m³/s (25.2 ~ 106.2 m³/s)까지 공급할 수 있는 것으로 분석되었다.

3.2. 수문 및 수리 특성분석

Table 4에 요약한 남강하류 유역 강우량을 살펴보면 2012년 기준으로 6월 하순부터 8월 중순까지의 여름철 장 마기와 9월에 강우가 집중되어 연강우량(1,788.1 mm)의 약 70.7% (1,264.3 mm)가 이 기간 동안 내렸다. 이에 따라 남 강댐 방류량도 증가하였고, 그 결과 하류의 유량도 이때 가장 많았다. 한편 2011년 12월과 2012년 1월 ~ 2월까지의 강우량은 매우 적어 가뭄현상을 보였으며 2009년과 2010 년 봄철에는 강우 빈도가 산발적으로 많았으나 그 양은 많 지 않아 남강댐 방류량과 유량은 적은 양이 방류되고 실측 되었다.

남강댐은 남강본류와 사천만방수로를 통하여 방류하고 있다. 남강댐은 농업용수 226.8 백만톤(7.2 m³/s) 및 하천유 지용수 122.1 백만톤(3.9m³/s)을 연간 348.9 백만톤(11.1 m³/year) 을 공급하는 것으로 계획되어 있다. Table 5와 같이 남강댐 이 본류로 방류하는 물의 양은 발전방류량과 여수로방류량 의 합이며, 일반적으로 발전방류는 농업용수와 하천유지용 수를 공급하며, 여수로방류는 홍수피해저감을 위하여 방류 한다. Table 5의 남강댐 본류 연간 방류량을 살펴보면, 발 전방류량의 평균은 915.0 백만톤(29.0 m³/s)으로 농업용수와 하천유지용수 연간 계획공급량 348.9 백만톤(11.1 m³/s)을 훨씬 넘는 것으로 나타났다. 하지만 Table 4의 월별 방류량 을 살펴보면 위 계획공급량 이하의 월이 1월~ 4월과 10월 ~ 12월에 집중되어 있는 것을 알 수 있다. 여수로 방류는 주로 비가 많이 내리는 우기(6월~ 9월)에 이루어지기 때문 에 그 외 월의 방류량은 모두 발전방류량에 해당하는 것으 로 남강하류 수질개선을 위해서는 추가적인 댐 방류량 증 가가 필요할 것으로 판단된다.

하천의 수질은 물리적 특성, 즉 수리·수문학적 특성에 의해 많은 영향을 받는다. 오염총량제 2단계 기준유량은 1995년 ~ 2004년까지 통합모형화 기법으로 구축한 SWAT 모형을 이용하여 자연유량 및 댐 방류로 인한 효과와 물수지 요소 에 의한 영향 등을 고려하여 10년 평균 저수량(Q₂₇₅)을 기 준으로 낙동강 수계 41개 오염총량관리 단위유역별로 기준 유량이 산정되었으며, 275번째 해당하는 유량 즉 365일 중 275일 이상을 유지하는 최소유량으로 초과확률(Exceedance Probability) 75%에 해당하는 유량을 의미한다. 저수량을 기준 유량으로 선정한 것은 저수량보다 유량이 풍부한 경우에는 목 표수질을 만족할 수 있다는 의미로 볼 수 있다(^{GNDI, 2010}).

Table 6과 같이 지난 5년간(2008년 ~ 2012년)의 실측유량 자료를 바탕으로 유황분석을 실시한 결과 남강D(남강2A) 지점의 저수량은 Q₂₇₅ = 7.3m³/s이며, 남강E(남강4-1) 지점 의 저수량은 Q₂₇₅ = 13.8 m³/s으로 나타났다. 오염원으로부터 배출되는 오염물질의 양이 동일하다고 가정한다면, 유량이 저수량(Q₂₇₅) 이하일 때 측정한 수질은 높게 나타날 가능성 이 있다. 실제로 상기 조건에서 저수량 이하에서 측정된 수질을 선별하여 조사한 결과 남강D 지점의 평균 수질농도 는 BOD = 3.9 mg/L, 남강E 지점의 평균 수질농도는 BOD = 4.2 mg/L로 목표수질 남강D(BOD = 2.5 mg/L)와 남강E(BOD = 3.1mg/L)를 초과하는 것으로 나타났다. 또한 저수기에 해당하는 평균유속을 선별하여 조사한 결과 남강댐 하류 평균유속은 0.19 m/s 정도의 저유속으로 나타났다. 또한 본 연구에서는 남강댐 하류로의 하상경사 및 유속분포 그리고 유하시간을 산정하고자 HEC-RAS 모형을 이용하였으며, 오염총량제 2단계 기본계획에 적용된 단면자료를 활용하였다.

경계조건(Boundary Condition)은 앞서 수집·분석된 남강 댐 방류량 자료 및 수위국 지점들에서 분석된 일평균유량 을 통해 산출된 유황분석 결과값을 상류단 경계값으로 사 용하였으며, 하류단 경계조건은 Nomal Depths(=0.000271) 을 적용하였다. 각 단위유역별 산출된 최종 유황분석결과는 Table 6과 같다.

HEC-RAS 모형을 이용하여 남강댐 하류의 유속을 모의한 결과 저수기(Q₂₇₅)기준 남강D(남강2A) 부근은 0.10 m/s, 남강 E(남강4-1) 부근은 0.16 m/s로 실측한 자료와 매우 근사한 값을 나타내었다. 산정된 유하시간은 5.6 day로 2009년 안 동조정지댐으로부터 하구까지 물이 도달하는데 대략 12 day이 소요되는 것으로 산정된 것을 감안할 때 상당히 긴 시간으로 나타났다(Fig. 2 참조)(^{Choi et al., 2008}).

Fig. 2. Hydraulic residence time simulation in the Nam River watershed.

Fig. 3에서 보듯 남강댐 하류로 내려 갈수록 하상 경사가 매우 완만하며 다수의 수중보와 교량 등 인공 구조물과 만 곡흐름의 특징으로 저유속으로 인한 강의 유하시간이 증가 하는 것으로 분석되었다. 하천에서의 유하시간 증가는 곧 체류시간의 증가를 의미한다. 따라서 남강댐의 방류량이 적 은 늦가을부터 다음해 봄까지의 갈수기에 남강 하류 수역 의 유속도 매우 느려져 이런 조건은 남강 하류 수역에서 식물플랑크톤의 증식을 촉진 할 것으로 예상된다.

Fig. 3. Location of dam and example vertical grid of a selected in the Nam River watershed.

3.3. 시·공간적 변동특성

남강댐 직하류인 남강C(남강)의 수질은 2008년~ 2012년 까지 거의 변화가 없었고, Fig. 4에 나타난 남강댐 하류 수 질측정망 지점에서의 최근 5년간 연평균 분포를 보면 남강 하류로 갈수록 BOD, TP, TN, Chl-a의 농도가 높아지는 것 을 알 수가 있다.

Fig. 4. Comparison of water quality mean values of each sampling stations during January, 2008 - December, 2012.

특히, Chl-a의 경우 단위유역 남강D(남강2A)지점을 지나 면서 급격히 증가하는 경향을 보이는 것으로 나타났다. 이 러한 현상은 지형적으로 경사가 완만해지는 중·하류부에서 강의 흐름(유속)이 느려지면서 체류시간이 길어져 미처리된 유기오염물질과 질소, 인 등 과량의 영양염류가 집적되면서 수체 내 부유성 식물플랑크톤이 대량 증식하는 현상을 초 래한 것으로 추정된다.

남강댐 하류 단위유역 남강D(남강2A)지점의 BOD농도는 2009년 평균 3.6 mg/L에서 2012년 평균 1.9 mg/L로 남강E (남강4-1)지점의 BOD농도는 2009년 평균 4.3 mg/L에서 2012년 평균 2.4 mg/L로 수질이 개선되어진 것으로 조사되 었다. 하지만 Fig. 5에 도시한 바와 같이 월평균 BOD값을 살펴보면 두 지점 모두 대체로 1월~ 5월까지 높은 BOD 값 을 보여 봄철 유기물 오염도가 높은 것으로 나타났다. 또 한 강우가 집중되는 6월 하순에서 8월 중순 풍수기와 9월 에서 11월 평수기를 제외하고 12월부터 저·갈수기인 봄철 5월 까지 지속적으로 수질이 악화되어 목표수질 남강D(BOD = 2.5 mg/L)와 남강E(BOD = 3.1 mg/L)를 초과하는 것으로 나 타났다. 일반적으로 강우량이 많을 때 하천은 세척효과로 정화되지만 선행무강우일수가 길면 하천은 대부분 하·폐수 로 이루어지면서 오염물질이 하천에 집적되어 수질오염은 심화될 수밖에 없다.

Fig. 5. Monthly variation of BOD in Nam River D and Nam River E stations during January, 2008 - December, 2012.

Fig. 4에서 BOD와 Chl-a의 경우 상류에서 하류로의 수질 분포 특성을 보았을 때 남강댐 하류 수역은 농도가 비교적 낮은 구간 (a)상류(남강, 남강1, 남강2), 농도가 증가하는 구 간 (b)중류(남강2A, 남강3), 높은 농도가 유지되는 구간 (c) 하류(남강4, 남강4-1) 구간으로 구분되어진다. 이를 바탕으 로 수질조사 기간 동안 각 구간별 이화학적 수질값의 최대, 최소값 및 평균값을 Box Plot 형태의 통계값으로 표시하여 나타내었다(Fig. 6 참조). 유기물 오염지표인 BOD, COD_Mn, TOC 분포 특성을 살펴보면 BOD의 중앙값은 1.0 ~ 2.6 mg/L, 제1사분위수(25%)는 0.8 ~ 1.7 mg/L, 제3사분위수(75%)는 1.3 ~ 4.2mg/L의 농도 범위를 보이고 있으며, 중 하류로 내려 가면서 농도가 증가하고, 농도 변동 폭도 커지는 것으로 나타났다. COD_Mn, TOC의 경우 상류에서 하류까지의 농도 변화는 BOD와 비슷한 경향을 보이며 농도분포가 중앙값에 밀집된 분포를 보이고 있다. TP의 제1사분위수는 0.01 ~ 0.072 mg/L(7.2배), 제3사분위수는 0.03 ~ 0.135 mg/L(4.5배), 중앙 값은 0.02 ~ 0.1 mg/L(5배) 증가한 농도를 보이고 있다. TN 은 사분위범위가 중앙값을 기준으로 저 농도로 분포되어 있는 밀집형태를 보이고 있다. 식물플랑크톤 현존량 지표인 Chl-a의 최소 10 %percentile 값과 최대 90 %percentile 값은 상류에서 하류로 각각 2.40 ~ 12.87 μg/L, 5.44 ~ 76.83 μg/L, 6.10 ~ 103.6 μg/L로 다량의 영양염 유입과 완만한 하상구배의 느린 흐름으로 Chl-a 농도의 증가 범위가 연중 2.5 ~ 185.8 μg/L로 매우 큰 변동 폭을 보이고 있으며, 중앙값을 기준 으로 흩어짐 정도가 매우 크게 나타났다.

Fig. 6. Box-plots of water quality parameter in each group during January, 2008 - December, 2012.

3.4. 상관성 분석

남강댐 하류 수역에서의 수질관측자료의 통계분석을 위 해 SPSS 12.0 통계 프로그램을 사용하여 항목 간 상관관 계를 분석하였다. 통계분석을 실시하기 전에 주요 수질항목 이 정규분포를 따르는지를 검토하였으며, 정규분포를 따르 지 않는 항목에 대해서는 적절한 자료의 변환을 검토하였 다. 하지만 원 자료의 관측치들이 정규분포를 따르는 경우 에는 보다 더 의미가 있겠으나, 이러한 정규성의 가정이 반드시 필요한 것은 아니다(^{Kim et al., 2007}). 그러므로 본 연구에서는 자료의 변환 없이 원 자료를 그대로 표준화 시 킨 후 분석하였다.

남강하류 유역에서 조사된 인자별 상관분석 결과를 Table 7 에 (a)상류, (b)중류, (c)하류 순으로 나타내었다. 상류 구간 에서 WT(r = 0.421), BOD(r = 0.340)가 Chl-a 농도와 작지만 통계적으로 유의한 양의 상관관계(p < 0.01)를 나타냈다. 중 류 구간과 하류구간에서는 pH(r = 0.547), BOD(r = 0.718), COD_Mn(r = 0.766), TOC(r = 0.606) 등이 Chl-a 농도와 높은 양의 상관관계(p < 0.01)를 가짐으로써 조류성장에 따른 자 생 유기오염물질과의 관련성을 추정해 볼 수 있다. 또한 댐(본류)방류량과 BOD(r = 0.388)와는 작지만 통계적으로 유의한 음의 상관관계(p = 0.01)로 나타나 방류량에 의한 오 염원 제어 효과가 있는 것으로 판단된다. pH와 Chl-a 농도 가 높은 상관성을 보인 것은 조류의 성장과 관련한 pH의 변화로 설명할 수 있다. 조류의 경우 광합성이 활발하면 수중의 용존산소를 증가시키는 한편 조류의 증식으로 인한 CO₂의 소모에 따라 탄산염 시스템(Carbonate system)의 역 반응으로 H⁺이온이 감소하면서 pH를 상승시키는 결과로 나타난다. 즉, 대기 중 CO₂가 수중으로 녹아들면, 탄산 (H₂CO₃)이 형성되며, 이는 다시 수소이온(H⁺)과 중탄산염 (HCO₃^-)으로 분해되는데, 이 중 중탄산염은 수중식물(조류) 의 광합성 시 CO₂흡수의 주된 형태이므로 식물플랑크톤이 증식하면서 광합성 양이 증가할수록 중탄산염 흡수가 증가 되어 수소(H⁺)이온의 양은 계속 감소하게 된다. 이는 결국 수체의 pH 증가를 유발하게 된다. 또한 BOD, COD_Mn, TOC 등이 Chl-a와 높은 인자 부하량을 가짐으로써 조류성 장에 따른 자생 유기오염물질의 발생을 나타내고 있다. 그 리고 전기전도도와 NH₃-N(r = 0.519), NO₃-N(r = 0.503)가 양의 상관관계(p < 0.01)로 나타나 수중에서 유기물과 함께 전해질로 작용함을 알 수 있다. 한편 NO₃-N와 PO₄-P는 음 의 관계를 보여 식물플랑크톤 증식과정에서 우선적으로 소 비되는 이들 무기영양염류는 식물플랑크톤의 증식 과정에 서 소모되어 음의 상관을 나타낸 것으로 추정된다. 또한 남강댐 하류 구간에서 영양염류 농도와 Chl-a가 뚜렷한 상 관성을 보이지 않고 있다. 이는 하류 수역에 영양염류 농 도가 식물플랑크톤이 증식하기에 이미 충분한 상태이기 때 문인 것으로 판단된다. 조류발생의 필요조건인 인 농도는 미국 EPA 기준으로 0.02 mg/L, OECD 기준으로 0.035 mg/L로(^{Kim, 2012}) 앞서 구간에서 조사된 인 농도는 (a)상 류(0.028 ~ 0.107 mg/L), (b)중류(0.073 ~ 0.152 mg/L), (c)하류 (0.071 ~ 0.128 mg/L)로 이미 충분한 상태이다. 한편, 상류구 간에서 하류구간까지 TN/TP의 질량비(TN:TP mass ratio)를 보면 각 지점에서 조사기간 평균 54.7, 37.4, 38.1 (범위: 6.2 ~ 757.3)의 비를 보였고, 특히 상대적으로 TP 농도가 낮 았던 하천구간에서 더 높은 비를 보여 전 지점에서 모두 인 이 제한영양염인 것으로 조사되었다(^{Forsberg and Ryding, 1980}). 일반적으로 TN과 TP의 질량비가 16 이하일 경우 총질소가 제한요인으로 작용하고, 29 이상일 경우 TN에 비 하여 총인의 부족으로 1차 생산력에 영향을 미칠 수 있는 잠재성을 갖고 있는 것으로 판단되며 세 구간 모두 TN과 TP의 질량비가 평균 29 이상으로 인이 제한(P - limitation) 요인으로 나타났다(^{Park et al., 2013}). 또한, 조사기간 평균 값으로 볼 때 각 지점에서 TP농도는 0.06 mg/L, 0.11 mg/L, 0.1 mg/L로 부영양 기준을 초과하였다.

남강댐 하류유역 총 세 구간에서 인(P)이 제한영향인자인 것으로 조사되었고, 식물플랑크톤 증식에 영향을 미치는 인 자인 것으로 판단이 되나, 강우 및 상류 댐의 방류량 증가 로 인한 유속의 증가(체류시간 감소) 또한 중요한 영향을 미치는 것으로 추정된다. 조류의 발생원인은 영양염류의 과 다 이외에도 수리학적 인자인 유속, 체류시간, 강우량과 이 화학적 인자인 수온, 태양광량, 영양염류, 그리고 생물학적 인자들이 복합적으로 작용하여 발생한다. 따라서 저·갈수기 시에 남강댐 방류량이 하천유지유량인 발전방류량이 줄어 들어 유속이 느려지면, 하류로는 상류에서 흘러내려온 하·폐수가 집적되어 수질이 악화되는 현상을 초래할 것으로 판단된다. 앞서 언급한 수질 특성분석에서 월별 BOD농도를 분석한 결과 봄철 유기물 오염도가 높은 것으로 나타났는데 이는 상기시기와 일치하는 현상으로 볼 수 있다.

한편 월별 방류량에 따른 BOD와 Chl-a의 상관관계 분석 을 실시한 결과 1월(r = 0.840)(r = 0.783)과 2월(r = 0.833)(r = 0.757), 5월(r = 0.751)(r = 0.658)이 가장 높은 음의 상관관계(p < 0.01)로 나타났으며, 3월(r = 0.586) (r = 0.430)과 4월(r = 0.509)(r = 0.472)이 그 다음 순으로 나타났다(Table 8 참조). 강우가 집중되고 유량이 많아 유속이 빠른 7월과 8월은 음 의 관계이지만 유의성은 작은 것으로 분석되었으며, 강우 후 수체가 안정화된 9월과 10월은 유의성이 없는 것으로 나타 났다. 대체로 하천에 유량이 작고 유기물 오염도가 높은 시 기에 상관성이 매우 높은 것으로 나타났다. 위 결과로 볼 때 남강댐 방류량 증가에 따라 오염농도가 감소하는 것은 분명 한 것으로 판단된다. 따라서 남강댐 중·하류 구간의 유기물 오염관리를 위해서는 봄철 수중에서 식물플랑크톤의 증식과 더불어 오염부하량의 제어가 매우 중요하며 이를 위해서는 남강댐의 수문조작을 통한 적정한 유량에 의한 유속을 유지 해 주는 것이 효과적일 것으로 판단된다. 미국 Patuxent River 하구의 연구에서 겨울철 조류 발생이 질소(N) 제한조 건과 관련이 있으며, 유속증가에 따라 조류 발생이 사라짐을 보고하였다(^{Sohn et al., 2013}).

3.5. 주성분 및 요인 분석

분류된 세 구간은 서로 다른 수질특성을 나타내므로 각 구간의 수질특성을 파악하기 위하여 주성분 분석 및 요인 분석을 실시하여 그 결과를 Table 9와 10에 나타내었다 (^{Bernard et al., 2004}; ^{Kim et al., 2007}; ^{Shin et al., 1997}; ^{Yoo, 2002}). 주성분 수를 결정하기 위해서 고유치(eignvalue) 가 1.0 이상인 값을 갖는 주성분 축만을 고려하여 상류 4 개, 중·하류 3개의 요인을 추출하였다. 주성분 전체 수질 변동에 대하여 (a)상류 78.531%, (b)중류 78.394%, (c)하류 75.669%를 설명하고 있다.

상류 구간에서 제 1요인은 31.161%, 제 2요인은 52.722%, 제 3인은 68.591%, 제 4요인은 78.531%를 기여하고 있고, 제 1요인은 NH₃-N, Cond, TN로 분류되었고, 제 2요인은 수온(WT), 제 3요인은 BOD로 분류되었다. 중류 구간에서 제 1요인은 27.000%, 제 2요인은 53.461%, 제 3인은 78.394% 를 기여하고 있고, 제 1요인은 NH₃-N, 제 2요인은 Chl-a, BOD, COD_Mn, pH, 제 3요인은 TP로 분류되었다. 하류 구 간에서 제 1요인은 31.835%, 제 2요인은 58.616%, 제 3인 은 75.669%를 기여하고 있고, 제 1요인은 NO₃-N, TN, DO, 제 2요인은 Chl-a, COD_Mn, BOD, TOC 제 3요인은 PO₄-P, TP로 분류되었다.

상류 구간에서는 암모니아성 질소, 총질소와 전기전도도 가 제 1요인으로 분류되었고, 수온이 제 2요인으로 유기물 요인인 BOD가 제 3요인으로 분류되어 진주 도시 및 공업 지역 오염특성을 잘 나타내는 것으로 보이며, 중류와 하류 구간은 암모니아성 질소가 제 1요인으로 공통 분류되었고, BOD, COD_Mn, Chl-a가 제 2요인으로 공통 분류되었다. 이 와 같은 결과로 볼 때 남강 중·하류구간의 수질은 조류 증 식에 따른 자생 유기오염물질의 발생이 수질 오염의 요인 인 것을 확인할 수 있었다. 이는 상류 도시지역의 생활하 수가 흘러내려오면서 하류지역에 도달하여 누적됨에 따라 도시특성과 농가의 축산폐수 및 토지이용에 따른 유기오염 물질과 영양염류 유출로 오염도가 극대화됨에 수질관리에 있어서 가장 중요한 변수임을 확인할 수 있었다.

3.6. 방류량 증가에 따른 수질개선 효과

본 연구에서는 2단계(2011년 ~ 2015년) 오염총량관리 기 본계획수립에서 사용된 바 있는 QUALKO2 모형을 이용하 여 남강댐 방류량 변화에 따른 남강댐 하류 지역의 수질개 선 효과를 분석하였다. 먼저 신뢰성 있는 수질모형을 구축 하기 위해서는 대상하천의 유속 및 수심과 같은 수리학적 특징이 모형에 잘 구현되어야 한다. 따라서 Reach별 수리 학적 계수는 낙동강 오염총량관리 2단계 기본계획에 적용 된 단면자료를 이용하여 V = aQ^b, H = αQ^β의 관계식으로부 터 도출된 수리계수 값을 적용하였다. Reach는 수리분석 결과에 근거하여 하천 내 수리특성이 유사한 구간을 일정 하게 구분하여 14개 구간으로 나누었으며, 소구간 Element 길이는 1 km 등간격으로 구분하였다. 남강 본류에 유입되 는 하·폐수처리시설과 취·양수장 및 지류하천은 위치를 확 인하여 모두 점오염원 형태로 모형에 적용하였으며, 그 외 에 비점오염원으로 각 소유역의 유달부하율(유달율)을 고려 하여 수체에 유입하는 것으로 모형을 구성하였다(^{Gyeong-SangNamDo, 2009}; ^{Jung et al., 2013}). 유달율은 2단계 오 염총량관리 기본계획수립에서 산정된 자료를 이용하여 각 소유역의 하단부에서 점오염원으로 유입되는 형태로 모형 에 적용하였다.

모형에서 반응계수의 결정은 실측을 통해 산정하여야 하 나 이는 불가능할 뿐만 아니라 많은 오차가 포함될 가능성 이 크다. 따라서 대부분 수질모형에서 반응계수 추정은 자 연현상과 수학적인 해석간의 차이가 최소가 될 때까지 시 행착오법을 반복하여 최적의 반응계수를 결정하게 된다. 따 라서 모형의 보정 수질은 기준유량(Q₂₇₅)과 가장 유사한 시 기인 2008년 3월 환경부 수질 측정망 자료의 평균값을 이 용하였으며 BOD, TN, TP, Chl-a 항목에 대해 수질모형인 자의 민감도 분석을 통한 보정을 실시하였다. 민감도 분석 은 수질모형인자를 10 ~ 50% 비율로 변화시키면서 모형을 실행 마지막 Element에서 BOD, TN, TP, Chl-a 등 대부분 의 수질항목에서 영향정도가 크게 나타난 변수를 중심으로 실측값을 근거로 보정하였다. 민감도 분석 결과 K₁(BOD 분 해율 계수), K₃(침전에 의한 BOD 제거율), σ₁(조류의 침전 율), μ_max(조류 최대 성장률), ρ(조류의 호흡율), σ₂(용존성 P 의 하상 용출율, 용존인의 저질 용출율) 등의 수질모형인자 에 민감하게 반응하였다(Fig. 7과 Table 11 참조) 이처럼 조류의 성장 및 사멸에 민감하게 반응하는 이유는 QUALKO2 모형이 조류의 성장에 의한 유기물의 내부 재생산을 적절 하게 잘 모의하고 있기 때문으로 판단된다.

Fig. 7. Sensitivity analysis of Nam River QUALKO2 according to μmax.

Fig. 8에서 보듯이 BOD와 Chl-a의 경우 하류로 가면서 농도가 증가하는 경향을 보이는데, 이는 남강 중·하류에서 조류의 성장에 의한 유기물질의 내부 생산이 원인인 것으 로 추정된다. 앞서 언급한 Fig. 5와 Fig. 7에 도시한 바와 같은 경향으로 모형 예측에서도 이러한 현상을 잘 모의하 고 있다.

Fig. 8. Calibration results of Nam River QUALKO2.

모형의 검증은 보정된 수질모형의 적정 사용여부를 판단 하는 과정으로, 보정된 모형이 수질예측에 적절한지를 알아 보기 위해 보정에 사용한 매개변수를 그대로 고정하고, 2008년 4월 측정한 유량·수질 자료를 이용하여 모형을 검 증하였다. Fig. 9에 도시한 바와 같이 BOD, TN, TP, Chl-a 모든 항목에서 수질모형 예측값이 실측값을 안정적으로 잘 모의하고 있다.

Fig. 9. Verification results of Nam River QUALKO2.

모형의 검·보정 결과에 대한 정확도를 평가하기 위해 통 계를 이용한 오차 분석을 실시하였다. BOD, TN, TP, Chl-a 항목의 모형 예측 결과의 정확도에 대한 통계분석 결과 모 든 항목에서 양호한 결과를 나타내었다. Chl-a의 RMSE 값 이 다소 높은 이유는 다른 항목과 달리 단위를 μg/L로 설 정하였기 때문이다(Table 12 참조).

최종 보정 및 검정절차를 거친 남강댐 하류유역 QUALKO2 모형을 이용하여 남강댐 방류량 변화에 대한 시나리오를 설정 목표수질 달성여부를 평가하였다. 시나리오 1은 2008 년부터 2012년까지의 남강댐 1월 평균 방류량 인 8.6 m³/sec를 Head water 유량 값으로 설정하였고, 시나리오 2, 3, 4, 5는 시나리오 1의 1.5배, 2배, 2.5배, 3배씩 유량을 증가 시키면서 모의(Simulation)하였다.

Fig. 10에서 보듯이 방류량이 증가할수록 남강댐 방류수 에 의한 희석효과로 남강댐 하류 수역의 BOD농도가 감소 하는 것으로 나타났으며, 남강D 단위유역에서는 17.2 m³/sec 인 경우에 목표수질인 2.5 mg/L을 만족하는 것으로 나타났 고, 남강E 단위유역에서는 21.5 m³/sec인 경우에 목표수질 인 3.1 mg/L을 만족하는 것으로 나타났다.

Fig. 10. The simulated BOD results by water flow level along the Nam River.

Table 13에서 방류량이 8.6 m³/sec에서 25.8 m³/sec로 증가 한 경우 BOD농도가 각각 13%, 13%, 9%, 6% 감소하는 것으로 나타났고, 특히 영양염류 지표인 TP, TN 보다 BOD, Chl-a 같은 유기 오염물이 방류량 증가에 따른 감소 정도가 높은 것으로 나타났다. 한편 시나리오별로 유량을 증가했을 때 감소하는 정도가 점차적으로 줄어들었다. 이는 곧 남강댐 하류 수질을 개선하기위해 남강댐 방류량을 무 한정으로 늘리는 것은 효과적이지 않은 것으로 생각된다. 하지만 남강댐 방류량이 어느 정도 양의 수준에서는 수질 개선 효과에 도움이 되는 것으로 판단된다.

이러한 연구결과로 볼 때 수문조작에 의한 유량 조절은 연중 특정시기에 편중된 강우량으로 유량변동이 매우 큰 우리나라 강우 특성상 한계가 있으며, 남강댐 하류 수역의 영양염류 농도가 식물플랑크톤 증식의 제한인자로 작용하 지 않을 만큼 높다는 점과 식물플랑크톤에 의한 내부 유기 물 생산이 하류 수역의 유기물 오염의 주요 기원으로 나타 난 것을 감안할 때 근본적으로 식물플랑크톤 증식을 제어 하기 위한 영양염류 저감 대책을 장기적으로 마련해야 할 것으로 판단된다.