1. 서 론
2. 연구방법
2.1 연구대상지역
2.2 지상 LiDAR 및 드론
2.3 FLO-2D 모형이론
3. 결과 및 분석
3.1 유역정보 분석
3.2 지상 LiDAR 및 드론 데이터 분석
3.3 FLO-2D모형 분석
4. 결 론
1. 서 론
기후변화로 인하여 태풍과 집중호우가 반복적으로 발생하고 이로 인해 산사태나 토석류와 같은 산지재해 발생 확률이 증가함에 따라 산악지역에 위치한 도로나 주택가에 심각한 위험을 초래하며 많은 재산과 인명피해를 발생시킨다. 그러므로 무엇보다 산지재해의 피해를 줄이기 위해서는 재해 발생 예상 지역이나 피해 정도 및 규모에 대한 예측이 중요하며 재해의 피해를 최소화하고 효율적으로 산지를 관리하기 위해서는 산사태가 발생한 지점이나 발생이 예상되는 지점의 환경적인 원인을 규명하는 것이 중요하며, 이는 산사태 발생에 대한 사전 예방적 접근에 가장 핵심적인 요소이다(Sidle et al., 2006; Lee et al., 2015). 또한 LiDAR와 GPS장비를 도입하여 현지조사를 정밀하게 수행하고 정확도를 높이는 분석이 필요하며 무인항공사진측량 방법으로 소규모 지형에 대한 정사영상과 수치지도를 제작하고 그 정확도를 검증함으로써 향후 수치지도 갱신에 적용가능성을 실험(Lim et al., 2015)하였으며 드론을 활용한 산림공간정보 조사 및 분석 방법은 기존 방법보다 인력과 시간을 절감할 수 있으며, 정량적인 분석이 가능하기 때문에 국가산림계획 및 정책 수립에 크게 기여 할 것이라고 하였다(Park and Jung, 2018). O'Brien (1993)은 이차원으로 토석류 거동을 해석한 연구를 수행하여 FLO-2D라는 홍수 및 토석류 해석 상용 프로그램을 개발하였으며 토석류의 피해를 저감하기 위해서는 토석류가 어떻게 이동하고 퇴적되어 하류부에 피해를 주는지를 예측할 수 있는 기술이 필요하며 토석류에 대한 유동 매개변수가 토석류의 유동에 미치는 영향 즉 민감도 분석 결과를 제시하였다(Kim et al., 2013).
본 연구에서는 산악지역에서 발생한 토석류를 분석하기 위해 2개의 유역을 선정하고 공간자료를 구축을 위해 두 가지 방법 사용하여 공간 데이터를 취득하였으며 첫 번째 유역의 경우 지상 LiDAR를 활용하여 토석류 발생 구간을 스캔하고 지형자료를 구축하였으며 두 번째 유역은 드론을 활용하여 유역의 퇴적부를 중심으로 촬영하고 데이터 처리를 통해 정사영상과 포인트 클라우드 형태의 DSM을 생성하였다. 그리고 토석류 발생이 하류부에 미치는 영향을 분석하기 위해 2차원 상용 모델인 FLO-2D를 사용하여 토석류의 흐름 영역을 시뮬레이션하고 지상 LiDAR 및 드론 측정데이터와 퇴적부를 비교하였다.
2. 연구방법
2.1 연구대상지역
2006년 7월 태풍 에위니아로 인해 강원도 인제군 북면 한계리 지역은 많은 토석류 피해가 발생하였으며 설악산 국립공원을 동서 방향으로 횡단하는 국도 44호선을 따라 토석류가 집중되었으며 연구대상 유역은 Fig. 1과 같다.
연구대상 지역은 설악산 국립공원내에 장수 5교 상류부와 하류부에 위치하고 있는 두 개의 유역으로 사진과 같이 지금까지 토석류 발생 이력이 그대로 존재하고 있어 재해 발생 시점과 시간의 흐름에 따른 지형의 변화 분석이 가능하다. 연구지역의 경우 경사가 심하고 높은 고도로 인하여 GPS 및 지형자료 획득의 한계점을 가지고 있어 지상 LiDAR 및 드론을 활용한 2가지 형태로 유역의 지형자료를 취득하였다. Fig. 2는 연구 흐름도로 먼저 1:5000 수치 지도를 활용하여 DEM을 생성하고 수문분석을 통해 두 개의 유역정보를 추출하였으며 GIS 공간정보를 구축하였다. 이후, 지상 LiDAR를 활용하여 Site 1의 토석류 발생 구간을 따라 스캔하고 지형데이터를 구축하였다. Site 2의 경우는 드론 촬영을 통하여 유역의 지형데이터를 구축하였다. 그리고 FLO-2D 모형을 이용하여 각각 유역의 토석류 흐름을 모의하고 지상 LiDAR 및 드론 촬영데이터와 FLO-2D에서 모의 된 확산면적을 비교 분석하였다.
2.2 지상 LiDAR 및 드론
지상 LiDAR의 측정 원리는 레이저를 이용한 지형·지물의 측정은 측정 대상에 레이저를 쏘아 표면에서 반사되어 되돌아오는 레이저광을 광센서(light-detecting sensor)를 이용하여 감지, 분석하는 방법으로 장치는 송신부(transmitter), 수신부(receiver), 처리부(processor)로 구성된다. Fig. 3(a)는 현장측량에 사용된 장비로 RIEGL사의 LMS-Z210ii로 반사율 80%기준에서 최대 측정거리 650 m, 수직 0°~80°, 수평 0°~360° 범위로 15 mm의 오차수준으로 포인트를 획득할 수 있으며 카메라와 GPS를 장착하면 측정 포인트의 RGB값 및 절대좌표를 획득할 수 있다. GPS는 HUACE사의 X90 2주파 RTK (Real Time Kinematic) 장비로 국토지리원에서 구축한 VRS (Virtual Reference System) RTK시스템을 이용하면 수평방향 ±10 mm+1 ppm, 수직방향 ±20 mm+1 ppm의 정확도를 가진다.
Fig. 3(b)의 드론은 DJI사의 Phantom3 Professional 모델로 4K의 영상을 촬영하고 배터리를 포함한 무게는 1280 g으로 최대 16 m/s의 속도로 약 23분 동안 촬영이 가능하다. 1200만 화소의 카메라를 장착하였으며 F/2.8렌즈는 화각이 94°로 넓은 화각은 항공 촬영 시 발생할 수 있는 불필요한 왜곡을 제거한다.
2.3 FLO-2D 모형이론
FLO-2D모델(O'Brien et al., 1993)은 미국 콜로라도 대학에서 개발되었으며 그리드 기반의 물리적 모델로써 지표면과 하도내의 홍수-수문곡선과 강우-유출을 추적할 수 있다. 이 모델은 운동학, 확산, 운동파 방정식을 이용하여 운동량방정식으로 접근한다.
FLO-2D는 격자와 시간에 관계없이 유동심, 속도, 압력이 예상가능하고, 그리드 기반의 모델로써 GIS 및 기타 응용 프로그램들과 연동 가능하다. 하지만 하상 침식에 의해 유발된 토석류 체적의 증가는 FLO-2D 모델에서 고려할 수 없기 때문에 토석류의 퇴적이나 확산영역에서의 토석류 모의에만 사용되어야 하는 단점이 있다. FLO-2D는 2차원 수치해석 프로그램으로 8개의 흐름 방향으로 정의되며 각 방향의 연속방정식과 2개의 운동량방정식으로 구성되어있으며 지배방정식은 다음과 같다(FLO-2D, 2009).
| $$\frac {\partial h} {\partial t} + \frac {\partial (uh)} {\partial x} + \frac {\partial (vh)} {\partial y} =i$$ | (1) |
| $$S_{fx}=S_{ax}-\frac{\partial h}{\partial x}-\frac{\partial u}{g\partial t}-u\frac{\partial u}{g\partial x}-v\frac{\partial u}{g\partial y}$$ | (2) |
| $$S_{fy}=S_{ay}-\frac{\partial h}{\partial y}-\frac{\partial v}{g\partial t}-u\frac{\partial v}{g\partial x}-v\frac{\partial v}{g\partial y}$$ | (3) |
여기서, h= 토석류의 유동심, = 유동심을 통해 산정된 x, y축 이동속도, = 강우강도, = 흐름 방향과 흐름직각 방향의 마찰 경사, = 바닥 경사이며, 일반 유체와 달리 고밀도인 토석류 전단응력은 다음 식 (4)와 같이 나타낼 수 있다.
| $$\tau = \tau _{y} + \tau _{v} + \tau _{t} + \tau _{d} = \tau _{y} + \eta ( \frac {\partial u} {\partial y} )+C( \frac {\partial u} {\partial y} ) ^{2}$$ | (4) |
여기서, 는 전단응력, 는 항복전단응력으로 Mohr-Coulomb 전단저항(shear resistance)을 나타내며, 는 점성전단응력, 는 난류전단응력, 는 분산(dispersive)전단응력, 는 동점성(dynamic viscosity)계수, 그리고 는 내부(internal)전단계수를 나타낸다. 총 전단응력을 수심 적분하여 경사형식(gradient form)으로 표현하면 다음 식 (5)과 같다.
| $$S _{f} =S _{y} +S _{v} +S _{td} = \frac {\tau _{y}} {\gamma _{m} h} + \frac {K \eta u} {8 \gamma _{m} n ^{2}} + \frac {n ^{2} u ^{2}} {n ^{4/3}}$$ | (5) |
여기서, 는 전체 마찰 경사, 는 항복경사, 는 점성경사, 는 난류-분산경사, 은 유사혼합물(sediment mixture)의 비중량, K는 저항매개변수, 은 등가 Manning 계수값이다.
3. 결과 및 분석
3.1 유역정보 분석
지형특성 분석을 위해 공간해상도 5 m× 5m 수치지도를 활용하여 지형분석을 수행하였으며 Fig. 4의 (a) 연구대상 유역으로 토석류 발생 이후의 항공사진이며 화살표를 따라 토석류가 발생하였다. (b)는 대상유역의 등고선이며 (c)는 DEM으로 해발 고도는 대략 510~1485 m로 나타났으며 (d) 유역 및 집수면적을 분석한 자료이다. 각 유역의 상세 정보는 Table 1과 같으며 공간 데이터는 FLO-2D 모의 시 유역경계 구현 등의 기초자료로 입력된다.
Table 1. Watershed Information
| Area (km2) | Maximum elevation (m) | Minimum elevation (m) | Average elevation (m) | Average slope (°) | Deposition area (km2) | |
| Site 1 | 1.23 | 1280 | 510 | 858 | 31 | 0.064 |
| Site 2 | 1.98 | 1485 | 528 | 979 | 39 | 0.021 |
3.2 지상 LiDAR 및 드론 데이터 분석
지상 LiDAR지형 자료는 토석류 발생 구간을 따라 올라가며 스캔하였으며 기존 1:5000 수치지도의 지형자료와 비교 분석하였다. Fig. 5의 (a)는 토석류 발생지점을 화살표를 따라 스캔한 항공사진 자료이며 (b)는 토석류 발생 이전의 지형데이터로 수치지도를 활용하여 DEM을 생성한 자료이고 (c)는 수치지도와 지상 LiDAR 측량데이터를 활용하여 생성한 정밀 DEM이다. (b)와 (c) 지형데이터를 비교하였을 경우 토석류는 계곡을 따라 약 500 m 정도 이동하였고 도로를 지나 한계천으로 유입되었으며 계곡 중심부에는 많은 침식이 계곡주변부에는 많은 퇴적이 발생하였다.
드론 촬영을 수행하기 이전에 RTK-GPS측량을 통해 24곳의 지상기준점(GCP, Ground Control Point)을 선정하고 측량을 진행하였으며 4개의 기지점(CP, Check Point)에 대하여 비교하여 오차율을 확인하였다. 오차율은 평균 제곱근 오차(RMSE, Root Mean Square Deviation)를 이용하여 계산하였다. 계산 결과는 Table 2와 같이 평균 오차 RMSE의 값이 X축 0.021 m, Y축 0.031 m, Z축 0.028 m로 나타났다.
Table 2. Result of accuracy by GCP survey
드론 촬영 영상은 PIX4D Mapper 소프트웨어를 활용하여 0.064 km2의 면적에 1401장의 촬영된 사진을 중첩하고 Fig. 6과 같이 정사영상과 0.79 cm의 샘플링 거리를 갖는 DSM을 생성하였다.
3.3 FLO-2D모형 분석
FLO-2D 모형 분석을 위해 토석류 모의 시 확률강우량은 토석류 발생 당일의 실측강우량을 반영하여 토석류 흐름과 퇴적면적 등을 산정하고 각 유역의 지상 LiDAR 및 드론 측량데이터와 이동 경로 및 퇴적구간을 비교하였다. Fig. 7의 (a), (b)의 경우 각 유역을 대상으로 FLO-2D를 모의한 결과로 이동 경로와 확산면적을 비교한 것이며 (c)는 (a), (b)유역을 같이 모의한 것으로 site 2의 유역에서 발생한 토석류가 하천으로 유입되어 site 1의 토석류와 합쳐지는 결과를 보였다. 각 유역의 퇴적면적의 경우는 토석류가 더 이상 이동하지 않는 마지막 유동심을 이용하여 계산하였으며 site 1의 경우 실제 토석류 퇴적면적은 51,336 m2이고 FLO-2D모의 결과는 59,543 m2로 나타났다. site 2의 경우는 드론 측량의 토석류 퇴적면적은 64,000 m2이고 모의 결과 퇴적량은 66,996 m2로 분석 되었다. 두 유역 모두 실제 토석류 퇴적량보다는 약간 높게 분석되었지만 비슷한 양상을 보이고 있다. 그리고 (d)는 토석류의 흐름 속도를 나타낸 것으로 하천으로 이동 시에는 약 1~7 m/s의 속도를 보이며 site 1과 합쳐지는 부분에는 약 5~20 m/s의 속도를 보여 주고 있다. 토석류가 한계천으로 유입되면서 유속이 빨라지고 위험도가 높아지는 양상으로 보여진다.
4. 결 론
본 연구에서는 실제 산악지역에서 발생한 토석류를 모의하기 위해 설악산 국립공원에 속해있는 한계천 일대를 대상 유역으로 선정하고 수치지도와 지상LiDAR 드론 등을 활용한 지형분석과 FLO-2D 모형을 적용하였으며 그 결과는 다음과 같다.
대상 유역 분석 및 공간데이터 생성을 위해 1:5000 수치지도를 이용하여 DEM을 생성하고 수문분석을 통하여 유역정보를 추출하였다. 지상 LiDAR와 드론을 활용하여 고해상도의 지형 정보를 구축하고 토석류 퇴적면적을 산정하였으며 site 1의 경우 51,336 m2으로 site 2의 경우 64,000 m2으로 분석되었다.
FLO-2D 모형을 적용하여 토석류 발생지역의 토석류 흐름과 퇴적면적을 산정하였으며 실제 측정 구간과 비교하였다. FLO-2D 모의 결과 토석류 퇴적면적은 site 1의 경우 59,543 m2, site 2의 경우 66,996 m2으로 측정구간보다 조금 높게 분석되었으나 흐름 형태나 퇴적양상은 실제 발생현장과 비슷하게 분석되었다. 이는 FLO-2D 모형을 활용함으로써 토석류 위험지역 선정 시 유용하게 활용 가능할 것으로 판단된다.
