2.1 Flow-R 모형의 활용

실제 대규모 토석류 흐름을 모의하기 위해서는 침식, 퇴적에 대한 상당량의 질량 변화를 고려해야 하기 때문에 정확한 흐름을 모의는 어렵다(Iverson and Denlinger, 2001). 대규모 댐의 경우에는 댐이 파괴 됨에 따라 발생하는 인명 및 경제적 피해가 상당하기 때문에 물리적모형을 활용한 댐 파괴에 대한 영향력 분석에 필요한 정보를 획득하는 것이 필요하다. 그러나 산지 곳곳에 분포하고 있는 소규모 농업용 저수지의 경우에는 그 규모와 개수로 인하여 댐 위치, 저류용량, 제방고 등 기본적인 정보만이 있는 상태이기 때문에 제방 파괴로 인한 영향력을 분석하는데 필요한 자료의 확보에 한계가 있다.

Flow-R은 Matlab을 기반으로 개발된 경험적 공간분포 모형으로 토석류가 발생할 수 있는 발생원을 식별하고, 식별된 토석류 발생원으로부터 마찰법칙과 흐름방향 알고리즘을 통해 토석류의 확산을 산정할 수 있다. 경험식을 통해 발생원 및 영향범위를 산정하기 때문에 실제 발생한 토석류의 흐름이라고 설명하기에는 한계가 있으나, 수치지형지도만을 사용하여 토석류의 발생원 및 영향범위를 산정할 수 있는 장점을 가지고 있다. 따라서 제한적인 자료를 보유하고 있는 소규모 저수지의 제방 파괴에 대한 하류의 영향을 분석하기 위해서 Flow-R을 활용할 수 있다.

Flow-R을 활용하여 토석류 범위를 산정하기 위해 필수적으로 입력되어야 하는 자료는 Raster로 구성된 유역 범위, DEM이다. 대상으로 하는 유역에 대하여 토석류의 흐름에 가장 큰 영향을 미치는 지형정보인 DEM을 활용하여 흐름방향 및 범위가 결정된다. DEM은 국토지리정보원을 통해 30 m by 30 m 해상도의 자료를 획득할 수 있으며, 필요에 따라 항공영상 등을 활용하여 보다 고해상도의 DEM자료를 획득하여 활용할 수 있다. 토석류 발생원의 경우 DEM을 보함하여, 토지이용, 경사, 곡률 등 다양한 지형정보를 활용하여 식별해 낼 수 있지만, 발생원 자체를 모의를 위해 필요한 입력자료로 활용할 수 있다. 본 연구에서도 DEM을 활용하여 발생원을 식별하기보다는 대상으로 하는 저수지의 제방의 위치를 발생원 입력자료로 활용하였다.

3.1 대상지역 선정

전라북도 완주군 하천 상류의 산악지역에 위치한 백여저수지, 탐저수지, 운지저수지, 오성저수지를 대상으로 붕괴 영향을 분석하였다(Fig. 1). 연구대상 4개 저수지는 모두 총 저수량이 30만 m³ 이하인 농업용 저수지로 완주군이 관리주체이다. 4개 저수지 모두 법적 EAP 수립 대상이 아니나, 저수지 하류부 인근에 마을이 위치하여 재해발생 시 인명 및 재산 피해의 위험이 높은 지역이다(Fig. 2). 백여저수지는 전북 완주군 구이면에 위치하고 있으며, 소하천인 백여천으로 방류하여 호동천, 삼천을 지나 구이저수지로 유입되며 직하류에 원백여마을이 있다. 탑저수지는 전라북도 완주군 봉동읍에 위치하고 있으며, 지방하천인 석탑천으로 방류하여 만경강으로 합류되며, 저수지 직하류에는 민가와 농경지로 이루어져 있다. 운지저수지는 전라북도 완주군 용진읍에 위치하고 있으며, 소하천인 봉계천으로 방류하여 양지천을 지나 만경강으로 합류하며, 직하류에 봉계마을, 지동마을 등 비교적 큰 마을이 위치하고 있다. 오성저수지는 전라북도 완주군 소양면에 위치하여 지방하천인 오도천으로 방류하며 저수지 직하류에 오성한옥마을이 위치해 있다. Table 1은 모의 대상 4개 저수지의 제원을 나타낸다.

Fig. 1

Location of Research Area

Fig. 2

Location Map of the Research Reservoirs

Table 1

Research Reservoirs Facility Specifications

3.2 수위별 저수지용량 산정

국내 대부분의 소규모 농업용 저수지 건설된지 오래되어 낙후되고, 관련 자료가 부실한 실정이다. 본 연구의 대상지역인 탑저수지와 운지저수는 1945년에 축조되었으며, 오성저수지도 축조된지 25년이 넘고 있다. 이에 건설 당시에 작성된 수위-저류량 자료를 확인할 수 없어, 붕괴 모의를 위하여 자체적으로 산정해야 한다. 실제 측량을 통해 정확한 수위-저류량 곡선을 산정할 수 있으나, 본 연구는 소규모 저수지에서 간략하게 산정하기 위하여 수치지형도를 사용하였다. 최근에는 지형정보를 활용한 분석 Tool이 점점 발전하고 있으며, ArcGIS, QGIS 등의 프로그램을 활용하여 원하는 지역의 다양한 지형분석이 가능하다. 본 연구에서는 10 m 해상도의 DEM을 이용하여 ArcGIS의 ‘Surface Volume’ 기능을 활용하여 저수지 용량을 산정하였다. Fig. 3은 4개 저수지에서의 수위-저류량 곡선을 나타낸다.

Fig. 3

Stage-Discharge Curve at 4 Reservoirs

4.1 저수지 제방 파괴 시나리오 구축

저수지 제방 파괴정도를 5가지 시나리오로 구분하여 단계별로 하류의 침수피해 영향범위를 분석하였다. Table 2는 5가지 시나리오에 대하여 총 저수지 제방높이 대비 잔여 제방높이를 나타내고 있다. 저수지별로 총 제방높이의 10%가 파괴된 ‘경미한 수준’부터 전체 제방이 완전히 파괴된 ‘완전 파괴수준’까지 시나리오별로 하류의 침수 영향을 분석하였다.

Table 2

Reservoir Levee Failure Scenario

4.2 하류침수 영향 분석

Flow-R 모형을 활용하여 저수지 파괴시 하류침수 분석 모델을 구축하여 저수지 제방 파괴 시나리오에 따른 하류침수 영향을 분석하였다. 저수지 제방 파괴로 인한 홍수파 발생원은 제방에서 홍수 발생이 시작하는 것으로 고려하였다. 이 후, 발생원으로부터 저류된 물이 확산되어가는 과정을 모의하여 영향 범위를 구획하여 침수 범위를 계산하였다. 이때 홍수파의 확산속도는 5~10 m/s로 가정하였다. Table 3은 시나리오별 침수최원점까지의 홍수파 도달시간을 나타낸다. 저수지별 침수피해 발생 시간으로 제방이 완전히 파괴될 경우 백여저수지는 약 1.90~3.80분, 탑저수지는 약 3.29~6.59분, 운지저수지는 약 5.29~10.57분, 오송저수지는 약 2.07~4.14분 사이에 침수가 발생하는 것으로 분석되었다. 산악형 저수지의 특성이 반영되어 제방이 파괴될 경우 깊게 형성된 계곡부를 따라 저류된 물이 매우 빠르게 확산될 것으로 분석되었다. 저수지 파괴가 발생할 경우 10분 이내에 홍수파가 도달해 저수지 하류에 거주하는 주민은 대피를 위한 시간을 확보하기 어려울 것으로 분석되었다.

Table 3

Flood Wave Arrival Time Due to Reservoir Destruction by Scenarios

Figs. 4~7은 제방 파괴 시나리오별 4개 저수지의 침수 영향 범위를 나타내고 있다. 제방의 10%가 파괴되는 시나리오인 Case 01에서는 침수영향 범위가 저수지 직하류의 산지 유역에 국한되나 Case 02부터는 민가를 포함한 하류 전역에 침수를 발생시키는 것을 볼 수 있다. 특히 운지저수지는 제방 파괴시 홍수파가 하류로 가장 광범위하게 전파되었으며 이는 운지저수지의 형상계수가 1.950으로 방사성의 유역 형상을 가지고 있는 것 때문인 것으로 분석된다. 탑저수지는 제방 파괴시 홍수파가 하류 석탑천과 함께 산지의 계곡 양쪽으로 전파되는 것으로 분석되었다. 오성저수지는 제방 파괴시 하류에 침수 범위가 가장 낮았으며 이는 하류의 오도천의 유역평균폭이 1.781 km로 하천폭이 넓어 홍수파의 대부분을 하도에서 분담하기 때문인 것으로 분석되었다. 지형정보만을 이용한 Flow-R 모형의 간단한 분석만으로 저수지 파괴에 따른 침수범위와 침수까지의 시간을 개략적으로 분석할 수 있었다.

Fig. 4

Baekyeo Reservoir Embankment Damage Area of Impact

Fig. 7

Oseong Reservoir Embankment Damage Area of Impact

Fig. 5

Top Reservoir Embankment Damage Area of Impact

Fig. 6

Unji Reservoir Embankment Damage Area of Impact