서 론
지형의 변화는 지구시스템적으로 자연스러운 현상이다. 이러한 현상은 다양한 형태로 발생하여 인간의 삶에 재난으로써 영향을 주기 때문에 대비 및 대응에 노력하고 있다. 이러한 지형변화는 그 종류가 다양해서, 많은 연구자들이 현상에 따라 구분하고 있다. 대표적으로 Varnes(1978)에 의해 산지재해 의 현상을 구분하는 방법이 많이 인용되고 있으며 이는 지질과 유동 매커니즘을 고려하여 지형 변화의 종류를 구분하고 있다. 다양한 지형변화 중에서도 땅밀림은 완경사의 지반이 산사태나 토석류에 비하여 상대적으로 느린 거동속도로 원형 토체가 큰 파괴없이 이동하거나 침하되는 현상으로 발생원인은 지하수, 파쇄대, 지진 등이 있다. 이러한 땅밀림은 국내뿐만 아니라, 일본, 대만, 필리핀, 인도네시아 등 환태평양지역을 중심으로 플레이트의 침강이나 대륙판이 충돌하는 지역에서 다양하게 발생하고 있으며, 특히 아시아 몬순 지역에서는 강수와 지진으로 인한 크고 작은 땅밀림 재해가 발생하기 쉬운 자연 조건이 갖추어져 있어 재해에 대한 충분한 예방조치와 재해발생이후 신속한 복구가 요구되고 있다.
땅밀림을 포함하여 다양한 지형의 변화에 대한 특성을 분석하고 연구를 수행하기 위해 다양한 접근 방법이 제시되고 있다. 크게 세 가지로 구분하며, GIS를 이용하여 지형 변화가 발생한 위치를 기록하고 데이터베이스를 구축하여 위험성을 평가하는 연구 방법(Carrara et al., 1991; Dai and Lee, 2002; Saha et al., 2010; Zhao et al., 2019), 동역학적 모델을 이용한 산사태나 토석류의 확산 범위, 충격량 등을 예측하는 수치모델링에 관한 연구 방법(Shieh et al., 1996; Rickenmann et al., 2006, Jiang et al., 2011; Mast et al., 2014; Nikooei and Manzari, 2020), 실제 현상을 재현하기 위한 실험 연구 방법(Moriwaki et al., 2004; Francesca et al., 2014; Iverson, 2015)등이 있다. 최근에는 기술이 발전함에 따라 UAV를 활용해 사람이 접근하여 확인하기 힘든 지역에 대한 정보를 획득하기에 유리해졌다. 따라서 사람이 접근하기 힘든 지반활동 지역이나 산지지역에 대해 UAV를 활용한 연구가 활발하게 진행되고 있으며, Kwon(2016)은 한국 제주도의 미지형의 근접 촬영을 통한 단위 지형의 특징을 파악, Choi et al.(2016)은 한국 강릉지역의 사빈지역의 시계열 측정에 따른 지형 변화 탐지, Dewez et al.(2016)은 절벽 지형의 정보의 구축에 따른 위험지 분석, Oh et al.(2017)은 갯벌의 UAV 촬영을 통한 시계열 공간 지형정보 생성 등 다양한 지형의 정보를 구축한 사례가 있다.
본 연구에서는 UAV를 이용하여 땅밀림이 발생하는 사면에 대한 시계열 분석을 수행하였다. 강원도 삼척시 도계읍에 위치한 해당 사면은 과거 선행 연구자에 의해 안전성 조사 및 검토 연구가 수행된 바 있다. Kang et al.(2005)는 사면안정해석을 통해 평상시, 건기시, 우기시의 원호활동을 해석하여 우기시의 안전율이 0.78로 도출되어 붕괴발생의 위험이 높아 옹벽의 설치를 제안하였다. Cho et al.(2011)은 1:5,000 수치지형도를 이용하여 석탄 폐석 적치로 지형 고도가 20~49 m 상승하였음을 분석하였고, 하중으로 인한 하부 사면의 지형변화의 위험성을 감소시키기 위해 사면절취의 필요성을 제안하였다. 이러한 과거의 선행연구의 결과로써 사면의 절취를 통해 안정성을 도모하였으나, 최근 절토사면에 인장균열의 발생과 배수로의 이탈 등이 발생함에 따라 지형의 변화가 지속적으로 발생하는 것으로 보이고 있다. 따라서 본 연구에서는 사면 안정 관리가 이루어지고 있음에도 불구하고 지형의 변화가 이루어지는 모습을 정량적으로 분석하기 위하여 UAV를 이용한 사진측량을 통해 시계열 분석을 수행하였다.
연구지역 및 방법
연구지역
연구지역은 한국 강원도 삼척시 도계읍 지역의 사면을 대상으로 선정하였다(Fig. 1). 연구지역의 사면은 과거에 소규모 붕괴가 발생해 이에 대한 사면안정화가 요구되어 사면경사완화공법을 적용하였다(Cho et al., 2011). 또한, 절토사면의 중심을 기준으로 횡단 및 종단배수로가 설치되었다. 그러나, 사면안정화 공법이 적용된 절토사면에서 지속적인 지형변화로 인해 인장균열, 표토유실, 배수로 이탈 등이 발생되어 사면붕괴의 위험성이 있는 것을 확인하였다. 이러한 원인으로는 연구지역으로부터 상부 약 150 m 지점에 위치한 폐석 적치장으로 인한 하중의 영향일 수도 있다. 토층이 압축이 되고 적치장 경사면의 하부로 지반의 밀림현상이 발생한다고 볼 수 있기 때문이다. 폐석 적치장의 경사각은 약 34°이며, 적치장 경사면 하부부터 정상부까지의 높이는 약 70 m, 수평거리는 약 120 m로 조성되어 있다. 폐석 적치장 아래의 절토사면의 폭은 약 100 m, 경사는 26°, 표면적은 약 7,000 m2이다. 해당 사면의 인장균열과 지반밀림현상으로 인해 사면의 붕괴가 발생하면 연구지역 하부 지역의 주거지역에 큰 피해가 발생 할 수 있다.
연구방법
본 연구에서는 UAV를 활용하여 사진측량을 통해 고해상도 지형자료를 생성하여 절토사면의 땅밀림 현상을 분석하고자 하였다. 이러한 목적을 달성하기 위해 Fig. 2의 연구순서로 진행하였다. 선정된 연구지역의 현장조사를 위하여 UAV를 통해 이미지 자료를 획득할 계획을 준비한다. 이때 UAV를 사용할 수 있는 지역인지 확인해야할 필요가 있다. 군사지역, 항공기 운항지역 등 제한되는 지역일 수 도 있기 때문이다. UAV의 운용의 불안정, 고속이동, 바람의 영향 등으로 인한 GPS 센서의 부정확성이 높아져 측량정보의 오류를 보완하기 위하여 연구지역을 커버하는 GCP를 측정 후 항공사진 측량을 수행한다. 본 연구에서는 시계열 분석을 위하여 2019년 4월과 10월 두차례 걸쳐서 수행하였다. 수집된 자료를 분석하기 위하여 GIS를 이용해 분석하기 위한 사면 이외의 지역을 제거하고 Structure from Model(SfM)기법을 통해 3차원 지형정보를 작성 후, Digital Surface Model(DSM)을 제작하였다. 생성된 지형정보는 GIS 프로그램을 이용하여 두 시점에서의 지형 정보를 중첩하여 고해상도의 지형변화정보를 비교 분석하였다.
사용 장비 및 측량 정보
사면붕괴의 위험성을 가지고 있는 연구지역을 조사하기 위하여 사용된 장비는 크게 두 가지로 구성된다(Fig. 3). Fig. 3a는 GCP를 측정하기 위한 장비세트로 GNSS로부터 위치를 수신하기 위한 RTK 수신기와 컨트롤러 폴대로 구성되어 있다. 30개 정도의 위성신호가 수신되고 이 중 25개 정도의 위성신호를 사용해 RTK측위를 수행하였으며, RTK정확도는 수평으로 8 mm, 수직으로 15 mm로써 측정값의 편차(RMSE)가 최대 15 mm 이내이다. Fig. 3b는 무인항공기와 카메라를 이용하여 지형을 측량하는 장비세트를 보여준다. 무인기, 컨트롤러, 카메라로 구성되어 있으며, 사진촬영을 통해 얻어진 정보를 통하여 공간분석을 수행하였다. UAV의 무게는 1.4 kg, 최대이륙고도는 6,000 m, 최대속도는 72 km/h, 조정범위는 컨트롤러부터 7 km, 비행조정시간 30분이다. UAV에 장착된 카메라는 1인치 CMOS 센서를 사용하고, 픽셀수는 21.Mpx, 렌즈는 화각(Field of view) 84°, 초점거리는 2.8~11 mm, 101.6 × 58.1 × 41.0 mm 크기의 299 g의 무게를 지닌다.
연구지역의 사면에 대하여 2019년 4월 3일과 2019년 10월 13일 두 번의 무인기를 이용한 사진 측량을 수행하였다(Fig. 4). 두 번의 측량 기준이 다른 이유는 4월 3일 측량시 10월 13일보다 높은 풍속을 기록하여 지형정보의 오차를 줄이기 위해 더욱 상세하게 측량을 수행하였다. 두 번의 측량에 대한 기상정보와 무인항공 측량 조건은 Table 1에 나타내었다.
Table 1.
Conditions of photogrammetric flight
무인기를 이용해 얻어진 사진정보에 좌표값을 생성하기 위해서 기준점이 되는 지상기준점(GCP) 측량이 필요하다. 지상기준점은 연구지역에서 직접 측량한 기준점으로 절토사면의 수평위치(X,Y)와 높이(Z)의 기준이 된다. 따라서 기준점의 위치는 모든 사진상에 선명하고 명확하게 인식되어야하기 때문에 절토사면의 외부에 둘러싸인 도로와 절토사면 내부의 배수로 모서리를 기준으로 GCP를 선정하였다. Table 2는 선정된 GCP의 오차율을 나타내며, 4월 3일의 RMSE값은 x, y, z방향으로 각각 0.034 m, 0.015 m, 0.044 m로 나타났으며 평균 RMSE는 0.031로 나타났다. 10월 3일의 RMSE값은 x, y, z방향으로 각각 0.032 m, 0.027 m, 0.032 m로 나타났으며 평균 RMSE는 0.030으로 나타났다. 두 측량 기준점의 오차율은 거의 같은 것으로 나타났으며, 이는 기준점에 의한 두 측량정보의 차이는 없다고 볼 수 있다.
Table 2.
RMSE of GCPs
결과 및 토론
3차원 지형정보 구축
무인기를 활용한 사진영상 취득과 GPS를 이용해 GCP를 획득한 정보를 이용하여 3차원 공간정보를 구축하였다. 공간정보 구축을 위해 사진측량 이미징 기법중의 하나인 Structure from motion(SfM) 기법을 이용하였으며, 2차원 이미지 시퀀스에서 3차원 구조를 추정하기 위한 방법이다.
Fig. 5는 3차원으로 구현된 정보를 나타낸다. 4월 3일과 10월 13일 생성된 지형정보의 해상도는 각각 1.87 cm와 1.98 cm로 나타났다. 이러한 차이는 비행거리, 중첩률, 이미지 수 등 사진 촬영옵션이 다르고 환경적인 변수, 특히 바람의 영향으로 인해 수집된 포인트 클라우드 수가 다르기 때문이다. 4월 3일의 영상이 10월 13일의 영상보다 바람의 영향을 더 많이 받았기 때문에 사진 1장당 2D point의 수는 10월 3일보다 약 11,000개 작았으나, 앞서 설명하였듯이 높은 중첩율과 다수의 사진 촬영으로 인해 10월 13일보다 약 0.1 cm 높은 해상도의 결과를 도출 할 수 있었다(Table 3).
Table 3.
Results of photogrammetric information
지형정보의 시계열 비교 분석 결과
생성된 지형정보를 이용하여 약 6개월동안 변화된 지형을 분석하였다. UAV로 촬영한 전체 면적 중에서 절토사면 부분만 추출하여 4월 3일 DSM과 10월 13일의 DSM을 중첩하여 비교분석하였다(Fig. 6). 두 시점에서의 고도값을 비교하였을 때, Fig. 6a와 같이 절토사면의 왼쪽 하단부에 변화가 높은 것을 확인할 수 있었다. 따라서 지형의 변화가 높은 부분인 절토사면의 하단부를 중심으로 대표적인 수직단면 1개소와 배수로부분이 위치한 수평단면의 2개소의 지형변화를 비교하였다(Fig. 6b).
Fig. 6에서 설명한 대표 수직단면의 변화량은 Fig. 7에 나타내었다. 실선은 4월 3일의 고도값을 의미하고, 점선은 10월 13일의 고도값을 의미한다. 경사면과 수평면이 번갈아 나타나는 모습을 보여주고 있으며, 수평면은 절토사면의 안정을 위한 소단을 의미한다. 소단개수는 총 4개이지만 본 연구에서는, 그 중 아래쪽에 위치한 두 개의 소단에서의 지형변화를 분석하였다. Fig. 7의 소단에 표시된 2번과 3번은 Fig. 6의 2번과 3번의 위치를 의미한다.
도출된 그래프를 통해 4월 3일과 10월 13일의 고도값이 크게 변하고 있다는 것을 정량적으로 보여준다. 다시 말해, 절토사면의 형태를 그대로 유지하면서 지반이 전체적으로 밀리는 모습을 나타낸다. 즉, 지반의 변형이 융기나 침하가 아닌 수평한 방향으로 밀리고 있다. 이러한 변화를 볼 때, 본 연구지역의 지형변화는 땅밀림으로 판단된다. 해당 단면에서의 최대 밀린 거리는 1.39 m로 나타났으며, 6개월간 지형의 변화가 지속되었다는 것을 고려하였을 때, 하루 최대 0.7 cm의 평균적인 변화를 의미하여 매우 빠른 속도의 땅밀림이 진행되는 것으로 판단된다. 이러한 밀림이 발생한 것은 단면의상부에 인장균열이 발생하였을 가능성이 높다. 따라서 수직단면 상부의 인장균열을 확인하고자 약간의 균열이 발생한 지점을 확인하였으나, 방수포로 덮여 있어서 인장균열이 크게 발생한 지점은 확인하기 어려웠다. 방수포로 덮여있는 지점에 인장균열이 발생한 것으로 추측된다.
Fig. 8은 수평단면 2번과 3번에서의 6개월간 지형의 변화량을 의미한다. 그래프에서 x축의 0은 배수로가 시작되는 지점이고 y축은 6개월간 변환된 땅밀림 거리를 의미한다. 실선인 2번 단면은 배수로가 시작되는 지점으로부터 약 14 m 지점에서부터 땅밀림이 시작되어 36 m간 땅밀림이 발생하였으며, 평균적으로 1.17 m의 밀림현상이 나타났다. 겹선인 3번단면은 배수로가 시작되는 지점으로부터 약 15 m 지점에서부터 40 m사이에 땅밀림이 발생하였으며, 평균적으로 1.04 m의 밀림현상이 나타났다.
이러한 땅밀림 현상에 대하여 육안으로도 확인 가능한 부분이 있었다. Fig. 9는 2019년 10월 13일에 촬영한 현장사진으로 인장균열과 배수로가 이탈한 모습을 보여준다. 2019년 7월과 8월의 강우에 의해 지형변화가 급하게 이루어진 것으로 보이며, 추가 인장균열이 발생하기 이전에 방수포로 응급 복구를 하였으나 땅밀림은 지속적으로 진행되고 있는 것으로 판단된다.
결 론
본 연구에서는 사면 붕괴 발생가능성이 있는 지역에 대하여 2019년 4월 3일과 10월 13일 UAV를 활용한 고해상도 영상 분석을 통해 지형의 변화를 비교 분석하고자 하였다. 4월 10일에 촬영한 영상은 10월 13일에 촬영한 영상보다 풍속이 3배정도 높아 외부의 영향으로 인한 지형정보의 정확성이 낮아질 수 있었으나 3배정도 많은 사진정보의 획득을 통해 DSM을 생성하기 위한 포인트의 수가 2배정도 높아 DSM resolution이 1.87 cm로 10월 13일의 DSM resolution인 1.98 cm와 거의 비슷한 것으로 나타났다. 또한 GCP의 오차율은 3.1 cm, 3.0 cm로 약 7,000 m2의 산지사면의 넓이를 고려할 때, 매우 낮은 수치의 오차율로 볼 수 있었다. 생성된 지형정보에 대한 시계열 분석을 통해 약 6개월간의 지형변화는 땅밀림이 발생하고 있는 것으로 판단하였다. 이러한 이유는 사면의 형태를 그대로 유지하며 지형이 변화되고, 사면에 설치된 배수로의 탈락 현상이 발견되었기 때문이다. 향후 추가적인 UAV 사진 측량을 통해 추가적인 땅밀림의 진행 양상을 파악하여 대비책을 마련할 필요가 있다고 판단된다. 하지만 지속적인 관리를 통해 시계열 기반의 변위량, 변위 방향성 등은 쉽게 파악이 가능하지만 지표범위라는 한계성이 있다. 땅밀림은 지표뿐만 아니라 깊은 심도에서 토사가 전체적으로 서서히 이동하기 때문에 심층의 토사 거동의 판단도 필요하다. 이러한 부분을 보완하기 위한 방법으로 계측이나 모니터링의 방법이 있으나 높은 비용이 소모되므로, 사진측량을 통한 심층의 토사 거동을 판단할 수 있는 새로운 기법이 필요하다.
