1. 서 론
탁수는 강우에 의해 많은 양의 토사가 일시에 하천에 유입됨에 따라 발생하게 된다. 이는 저수지에서 수개월 체류하면서 하류 하천으로 방류된다. 하절기에 강우가 집중되어있는 우리나라의 강우 특성과 최근 기후변화에 따라 이상강우 및 기상이변으로 인한 집중강우의 영향이 크게 나타난다(Chung and Oh, 2006). 특히 과거 2000년대 ‘루사’, ‘매미’, ‘에위니아’부터 2020년대까지 ‘마이삭’, ‘하이선’ 등 태풍으로 인해 소양강, 임하댐에서 탁수 및 장기화 문제가 발생하였다. 이러한 탁수 발생은 댐의 탁수 장기화 문제로 이어져 댐 하류의 농·공업, 수생태, 등 사회·환경적 다양한 문제를 야기시킨다(Kim and Jung, 2007; Chung et al., 2008). 이러한 탁수 발생에 따른 대응 방안 마련을 위해 호 내 탁수 예측 정확도 개선을 위한 연구가 활발히 진행됨에(Kim et al., 2006) 따라 상류 유입지점의 탁수 데이터를 필요로 한다. 또한 탁수 발생 시 댐 상류 탁수 발생과 분포에 따라 댐의 효과적인 운영 및 대응(Chung et al., 2007)이 중요함에 따라 댐 상류의 지속적인 탁수 데이터뿐만 아니라 공간적인 탁수 분포 등을 파악하는 것 또한 중요하다. 특히 소양강댐의 경우 인북천 및 내린천의 합류 이후 하류 지점에 위치함에 따라 두 하천의 합류 이후 지점의 탁수 발생 현황 및 거동의 파악이 중요하다. 현재 부유사 모니터링 방법은 부유사 채집기, 채수 등을 통한 재래식 방식이 주로 이용되고 있으며 최근 레이저부유사측정기(LISST-200X)를 활용한 센서 기반 측정을 활용한 연구가 증가하고 있다(Kwak et al., 2017; Kim et al., 2012). 특히 LISST의 활용을 위한 검ㆍ보정 기법(Mikkelsen and Pejrup, 2001), 입경크기를 통한 SSC 산정, 채수를 통한 실험실 분석 기반 LISST의 적용성 검토 및 SSC (Suspended Sediment Concentratio), PSD (Particle Size Distribution)(Ehrbar et al., 2017) 등의 연구가 활발히 진행되고 있다. 하지만 이러한 점, 선 단위의 측정은 시·공간적으로 탁수의 현황 및 분포를 파악하기에는 한계점이 있다. 따라서, 시간적 해상도를 개선하기 위해 탁도-SS 관계식을 활용하여 실시간 모니터링이 가능한 탁도를 SS로 변환하는 방식이 활발히 사용되고 있으며(Chung et al., 2011), 관계식 개선을 위한 연구 또한 활발히 진행되고 있다(Landers and Sturm, 2013).
탁수 분포의 공간적 해상도 개선은 주로 분광 영상을 활용하여 공간 분포를 취득하는 방식으로 많은 연구가 진행되고 있다(Kwon et al., 2022; 2023). 원격 탐사를 위해 초분광 센서 및 다중분광센서는 위성이나 항공기 및 드론 등에 탑재되어 운용하여 토지 피복, 지질 상태, 작물의 생장 및 질병의 유무 및 토양의 수분 함유도 등의 연구에서 시작되어 최근 하천환경에서 하천의 수리량, 수질, 유사량에도 활용되기 시작하였다(Gwon et al., 2020). 이에 초분광영상을 활용한 원격탐사를 통해 하천 수질 모니터링 방법론을 개발하는 연구(Kisevic et al., 2016), 수체 반사율의 현장 측정을 기반으로 수질의 초분광 원격감지 알고리즘 개발(Fan, 2014), 드론 기반 초분광영상 활용 하천 수심 계측 기술 개발(You, 2018) 등 하천을 대상으로 다양한 연구가 활발히 이루어지고 있다. 하지만 탁도와 SS의 경우, 다른 물리량들에 비해 연구 사례가 부족한 실정이다.
따라서 본 연구에서는 SS 항목을 측정 가능한 LISST-200X와 탁도를 측정하는 YSI를 활용하여 연속적인 데이터를 취득하여 시간적 해상도를 높일 수 있다는 장점을 활용하여 관계식을 개선하고자 체 거름을 이용한 탁수 시료 제작을 통해 탁도-SS 관계식을 산정하였다. 접촉식 센서로 실측된 탁도, SS와 더불어 탁도-SS 관계식으로부터 산정된 탁도를 UAV 및 초분광센서를 활용한 면 단위의 탁수 모니터링을 비교 분석하여 탁도-SS 관계식의 적용성 평가와 댐 상류 유입지점에 대한 탁수의 공간적 분포를 파악하고자 하였다.
2. 연구 방법
2.1 실내 실험 기반 탁도-SS 관계식 산정
탁도-SS 관계에서 부유사의 입경 크기가 탁도에 큰 영향을 끼침에 따라 시중 황토의 체거름 분석을 통해 실제 하천 적용성을 검토하고자 하였다. 탁도-SS 관계식 산정을 위해 체 거름 분석에 이용한 체는 200, 100, 60번 체를 통해 75㎛ 이하, 75㎛ 이상 150㎛이하, 150㎛ 이상의 3 분류로 나누어 분율(Clay : Silt : Sand), 농도 별 시료를 제작하였다(Table 1). 시료에 대하여 LISST-200X를 통해 입경크기 별 농도와 YSI-EXO를 이용하여 탁도를 측정하였다(Fig. 1(a)). LISST-200X 측정 시 안정성을 높이기 위하여 소형 챔버를 활용하였다(Fig. 1(b)). 취득한 데이터의 경우 사분위 범위(Interquartile range, IQR)를 이용하여 IQR의 1.5배 차를 이상치로 판단하여 제거 후 탁도, SS에 대한 상관성 및 관계식을 산정하였다.
Table 1.
Sample production by fraction
| Case | Proportion (Clay : Silt : Sand) |
| 1 | 0.03 : 0.69 : 0.28 |
| 2 | 0.04 : 0.82 : 0.14 |
| 3 | 0.06 : 0.76 : 0.18 |
| 4 | 0.08 : 0.79 : 0.13 |
| 5 | 0.1 : 0.72 : 0.18 |
| 6 | 0.13 : 0.71 : 0.16 |
2.2 실내 실험 기반 탁도-SS 관계식 현장 계측 적용
취득한 관계식을 적용하고자 한 대상지역은 소양호 상류 사구미교 지점으로 소양강댐 상류 약 59km에 위치 해있으며 인북천과 내린천 합류 이후 지점인 사구미교 직하류 지점에서 측정하였다(Fig. 2). 소양강댐은 북한강 수계에 위치하여 서울을 포함한 수도권 지역의 생활, 공업, 농업 용수 공급과 홍수조절을 담당하는 인공호이자 다목적댐으로 소양강 댐 유역 면적은 2,703 km2으로 총 저수용량은 2,900 × 106 m3, 유효 저수용량은 1,900 × 106 m3, 저수면적은 70 km²이다.
실험실 기반 탁도-SS 관계식 적용을 위해 소양강 상류 사구미교 지점에서 YSI-EXO, LISST-200X, 초분광 센서를 통해 탁도, SS, 초분광영상을 취득하고 LISST-200X 측정 결과를 통해 입경 크기별 비율을 분석하여 탁도-SS 관계식 결과를 이용하여 탁도를 산정하였다. 탁도-SS 관계식으로부터 산정된 탁도, 실측 탁도와 부유사농도를 기반으로 초분광영상을 활용하여 각각 맵핑을 진행하였다. 실측 탁도와 탁도-SS 관계식으로부터 산정된 탁도를 비교 분석하여 탁도-SS 관계식에 대한 적용성을 검토하고 부유사, 탁도를 포함한 탁수 모니터링을 진행하였다(Fig. 3).
초분광 센서 기반 탁수 모니터링을 위해 소양강 상류 맵핑 대상 지점을 초분광 센서를 탑재한 UAV를 활용하여 촬영을 진행하였다. 탁도-SS 관계식 적용 및 탁도의 공간분포를 제시하고자 LISST-200X와 YSI-EXO를 통해 일정 지점에 대하여 연속 데이터를 취득하였으며 UAV의 호버링을 통해 실측지점에 대하여 초분광영상을 연속적으로 취득하였다(Fig. 4). 활용한 장비는 4 nm 간격으로 400~1000 nm의 측정이 가능한 초분광 센서(Corning사의 Micro HSITM 410 Shark)(Corning, 2017), 센서 탑재를 위한 UAV (DJI 사 Matrice 600), 촬영 각, 흔들림을 보정 해주기 위한 짐벌(DJI 사의 Ronin-MX)을 이용하여 측정을 진행하였다(Fig. 5(a)). 또한 취득 영상의 분광 특성을 비교하기 위해 LISST-200X, YSI-EXO를 사구미교 교량에서 밧줄을 연결하여 교량법으로 측정하였다(Fig. 5(b)). 측정 시기는 2023년 8월 31일로 사전 누적 강수량은 2023년 8월 28일~2023년 8월 30일까지 총 43.3mm로 강수량이 많지는 않았으나 육안상 탁수가 확인되었으며 수위는 3.57m로 누적 강수량으로 인해 상승기로 나타났다.
측정방법은 초분광 센서의 경우 라인스캐닝 기법으로 촬영함에 따라 바람, 진동, 촬영각 등 흔들림에 취약한 센서의 움직임을 보정 하기 위해 짐벌을 활용하였다. 또한 취득 영상의 방사 보정을 위해 일정한 반사율을 가지는 4종류(6%, 12%, 48%, 64%)의 반사천의 분광정보를 함께 측정하여 진행하였다. 취득된 초분광영상 분석을 위해 전처리 과정이 필요하다. 첫 번째로 초분광 센서가 라인스캐닝 기반으로 영상을 축적하여 면 단위의 촬영이 됨에 따라 6개의 영상을 정합한 RGB 이미지 기반으로 영상정합을 진행하였다. 두 번째로는 영상의 대기흡수 및 산란 효과 제거를 위한 방사보정을 수행하였으며 식은 다음과 같다(Kim et al., 2023)(Eq. (1)).
이때, R: 에 대한 반사율, : 파장, : 반사율로 보정하기 위한 기울기, : 절편값이다.
초분광영상 측정된 탁도, 부유사농도, 탁도-SS 관계식으로부터 산정된 탁도의 농도 추정을 위해서 호버링 기법을 통한 분광 정보를 수집하여 접촉식 센서 실측 데이터 및 탁도-SS 관계식으로부터 산정된 데이터와 연계하였다. 이때, 초분광영상의 경우 1초에 120개의 데이터를 수집하며 접촉식 센서의 경우 1초당 1개의 데이터를 수집함에 따라 시간 동기화를 위해 초분광영상 1초 동안 측정된 120개의 초분광영상 데이터 평균을 통한 다운 스케일링을 진행하였다. 최적밴드탐색기법(Optimal Band Ratio Analysis, OBRA)을 적용하여 부유사 농도와 최대 상관계수를 나타낸 두 개의 파장을 선택하여 분광 정보와 맵핑 항목(실측 탁도, 부유사 농도, 탁도-SS 관계식으로부터 산정된 탁도) 값 간의 회귀분석을 통해 오차 및 상관계수를 산정하였다. 선정된 두 파장의 광량을 서로 나눈 밴드비(Eq. (2))를 산정할 수 있으며 이를 이용하여 부유사에 대한 공간적인 농도 추정을 진행하고자 하였다(Fig. 6).
3. 연구 결과
실험실 분석을 통한 Table 1의 탁수 샘플 케이스별 탁도-부유사 농도 그래프를 통해 상관성 및 탁도-SS 회귀식을 산정하였다(Fig. 7)(Table 2). 실제 하천의 LISST-200X를 통해 취득된 SS 데이터의 분율(Clay : Silt : Sand)이 0.03 : 0.62 : 0.35 나타남에에 따라 Case1번의 관계식을 적용하여 탁도를 산정하였다. 실측 탁도와 산정된 탁도의 결정계수는 0.887(Fig. 8)로 유의미하게 나타났으며 탁도-SS 관계식에 있어서 입경 크기에 따른 분율(Clay : Silt : Sand)의 영향이 매우 큰 것으로 판단되었다.
Table 2.
laboratory analysis results
초분광영상 정합 과정에서 이미지 겹침으로 인한 줄이 생성된다는 한계점이 있었으며 전체 영상 전처리과정 결과는 다음과 같다(Fig. 9). 접촉식 센서를 통해 취득한 부유사 농도는 최대 221.50 ㎕, 최소 114.15 ㎕, 탁도는 최대 29.5 FNU, 최소 25.8FNU로 나타났으며 탁도-SS 관계식을 통해 산정된 탁도는 최대 32.9 FNU, 최소 25.2 FNU로 나타났다. 접촉식 센서 계측 결과 및 탁도-SS 관계식을 통해 산정된 값과 초분광영상 촬영 결과를 이용한 OBRA 수행 결과는 Fig. 10과 같으며, 분석 결과로부터 최대 상관계수를 나타낸 두 파장을 선택하였다(Fig. 10). SSC에 대한 두 파장에서의 상관성은 0.854(Fig. 10(a)), 실측 탁도에 대한 두 파장의 상관성은 0.911(Fig. 10(b)), 탁도-SS 관계식으로부터 산정된 탁도에 대한 두 파장 상관성은 0.854(Fig. 10(c))로 유의미한 결과를 나타냈으며 이를 통해 농도 산정 회귀식을 산정하였다(Table 3)(Eq. (3)).
여기서, 는 구하고자 하는 항목의 농도이며은 반사율로 보정하기 위한 기울기 계수, 는 A의 반사율, 은 절편값이다.
Table 3.
Turbid water measurement OBRA results
| SSC | 837.88 nm | |
| 770.17 nm | ||
| -1563.67 ㎕ | ||
| 13.27 ㎕ | ||
| 0.854 | ||
|
Turbidity (Obs) | 646.70 nm | |
| 571.03 nm | ||
| -20.62 FNU | ||
| 23.35 FNU | ||
| 0.911 | ||
|
Turbidity (Estimated) | 837.88 nm | |
| 770.17 nm | ||
| 21.35 FNU | ||
| 26.49 FNU | ||
| 0.854 |
정합된 초분광영상에 대하여 OBRA 결과 Eq. (3)를 적용하여 부유사와 탁도의 공간분포에 대한 맵핑을 진행하였다(Figs. 10(a) and 10(b)). 공간분포 결과 검증을 위해 센서 측정지점에 해당하는 초분광영상의 PIxel 농도와 실측값을 비교하였다. 동 시간대 실측 가능한 장비가 1대로 검증 데이터의 한계가 있었다. 부유사 실측 농도 132.2 ㎕, 맵핑 농도 127.3 ㎕로 오차율은 4%로 나타났으며 전체 픽셀의 평균 농도는 108.52 ㎕로 나타났다. 탁도의 경우 실측 값 25.8 FNU, 맵핑 값 24.6 FNU로 오차율 5% 정도로 나타났으며 탁도-SS 관계식 기반 탁도 맵핑 결과 부유사 데이터에 대한 탁도-SS 관계식 적용으로 인해 비슷한 데이터 경향을 나타냈다. 이에 따라 부유사 OBRA 결과와 동일한 파장대를 선정함으로써 부유사 맵핑 결과와 유사한 결과로 나타났다. 관계식을 활용하여 산정된 탁도의 맵핑 시 실측 값 25.8 FNU, 맵핑 값 23.5 FNU로 9%의 오차로 나타났으며 실측 탁도의 경우 전체 픽셀 평균 농도 23.8 FNU로 전반적으로 부유사 및 산정된 탁도 맵핑 결과에 비해 낮은 농도 분포로 나타났다. 탁도-SS 관계식 검토를 위해 실측 탁도와 관계식 활용하여 산정한 탁도 맵핑 결과에 대해 비교하였다(Figs. 11(b) and 11(c))(Table 4). 실측 탁도를 활용한 맵핑과 동일 지점 10개의 픽셀에 대해서 추출하여 비교하였을 때 최소 2%에서 우안 측 1번의 경우 최대 52%까지 오차가 발생하였으나 우안을 제외하고는 유사한 것으로 판단되었다. 전체적으로 하천 중앙지점 좌안 측이 높게 나타난 이유는 강우 이후 양안에 비하여 중앙지점의 유속이 높으며 좌측방향 만곡부 직상류 지점에 위치함으로서 좌안 측이 높게 나타난 것으로 판단된다.
4. 결 론
본 연구에서는 소양강댐 탁수 문제에 있어 상류 유역의 탁수에 대한 시·공간적인 데이터의 해상도를 취득하는 것이 매우 중요함에 따라 실험실 분석 기반 입경크기 별 비율에 따른 탁도-SS 관계식을 산정하였다. 산정된 탁도-SS 관계식에 대한 검증 및 적용성 검토를 위해 소양강댐의 상류 유입지점인 사구미교 지점에서 초분광 센서와 LISST, YSI를 통해 탁도, SS 데이터를 취득하였다. LISST 센서를 이용하여 취득한 SS 데이터와 유사한 입경 크기 별 분율(Clay : Silt : Sand)을 가진 케이스의 탁도-SS 관계식 적용을 통해 탁도의 공간분포를 산정하였다. 접촉식 센서를 이용하여 계측한 SS 및 탁도와 탁도-SS 관계식으로부터 산정된 탁도 총 3 항목에 대해 각각의 맵핑을 진행하여 탁수의 공간분포 파악를 파악하고 비교하여 탁도-SS 관계식의 적용성을 검토하고자 하였다. 실제 사구미교 직 하류에서 사전 누적 강수량 43.3mm, 약 12시간 이후 수위 상승기 시점에서 진행하였다. 접촉식 센서를 이용하여 점 단위의 연속 데이터를 취득하였으며 UAV에 탑재한 초분광 센서를 이용하여 센서 측정지점에 대한 영상을 취득하였다. 접촉식 센서를 이용하여 취득한 탁수 데이터와 초분광영상의 분광 정보를 이용하여 OBRA 분석을 수행하였으며 초분광영상의 분광 정보를 OBRA 분석을 통해서 최대 상관계수를 나타낸 두 개의 파장을 선정하였다. 이때, 두 파장과 실측한 부유사농도 및 탁도, 탁도-SS 관계식 기반으로 산정된 탁도의 상관계수는 각각 0.854, 0.911, 0.854로 높은 상관관계를 나타냈다. 선정된 두 파장을 이용한 밴드비를 통해 부유사 농도 회귀식을 산정하였으며, 초분광영상을 활용한 부유사 농도의 공간분포 맵핑을 위해 활용하였다. 측정 장비의 수량 및 안전성 문제로 인해 탁수 모니터링 결과 검증 데이터가 부족하다는 한계점이 있었으나 실측 데이터 비교하였을 시, 오차율 4%로 나타났다. 또한 탁도 실측 값 25.8 FNU, 맵핑 값 24.6 FNU로 오차율 5%로 나타났으며 탁도-SS 관계식 기반 탁도 맵핑 값은 23.5 FNU로 9%의 오차율을 보였다. 탁도-SS 관계식을 기반으로 산정된 탁도 맵핑 결과의 경우 SSC에 대한 회귀식을 통한 변환으로 인해 부유사 맵핑 결과와 공간적인 분포 특성이 유사하게 나타났다. 실측 탁도 기반 분석 결과 및 탁도-SS 관계식 기반 탁도 분석 결과 우안에서 다른 지점에 비해 오차가 크게 나타났으나 우안을 제외한 탁수 현황에 대한 결과는 유사하게 나타나 활용 가능할 것으로 판단된다. 해당 결과의 경우 사전강우량이 매우 많지 않았으며 강우 종료 후 약 12시간이 경과함에 따라 수체가 안정화 된 시점에서 측정됨에 따라 실험실 분석결과를 적용하였을 때 잘 일치한 것으로 판단된다. 초분광영상을 활용한 탁수 모니터링 기법의 경우 표면 부근 데이터만 취득 가능함에 따라 표층 탁수 현황에 대해서만 제시 가능하다는 한계점이 있으나 접촉식 센서를 함께 활용하여 개선 가능할 것으로 판단된다. 추후 측정 조건에 따른 추가 실험 및 실측 농도, 입경분포에 대한 데이터 축적과 영상 전처리 및 보정 방법 개선 시 하천 전반적인 탁수 발생 현황에 대한 공간적 분포를 재현 하는데 있어 정확도가 상승할 것으로 판단된다.
