1. 서 론
2. 연구대상지역
3. SWAT 모형을 이용한 신속회귀수량 산정방법
3.1 SWAT에서의 관개량 적용
3.2 SWAT에서 pothole을 이용한 담수 논 모의
3.3 자료수집 및 분석방법
4. 회귀수량 산정결과
4.1 SWAT의 검·보정
4.2 신속회귀수량 산정 결과
5. 결 론
1. 서 론
2019년 6월 물관리기본법이 시행됨으로서 통합물관리를 위한 법적 기반이 마련되었다. 그러나, 물관리일원화는 수량 확보, 수질 보전, 재해방지 등의 통합물관리 정책이 포함되어 있음에도 불구하고 농업용수 부문이 포함되지 않았다. 농업용수부문은 관행수리권으로부터의 허가수리권 조정 및 체계적 관리, 물 이용 시설관리 비용부담 구조개선, 농업용수 효율화를 통한 적정 용수 공급 등 다양한 과제에 직면하고 있다(KREI, 2020). 특히, 부족한 수자원을 효율적으로 관리하기 위한 측면에서, 농업용수는 우리나라 전체 물 이용량 중에서 약 40% (2016년 기준)에 해당하는 수량(MCT, 2016)으로 농업용수 효율 증대에 따른 다량의 가용수량 확보가 가능하기 때문에 통합물관리에서 농업용수부문을 배제할 수 없다. 이러한 농업용수 통합물관리 부문중에서 합리적인 농업용수의 수요량, 공급량 및 회귀량을 정량화 하기위한 노력은 하천유량관리 측면에서 가장 큰 의미가 있다고 할 수 있다(KRC, 2020).
농업용수는 4월말부터 9월 중순까지 논벼 재배기간에 집중적으로 사용되지만 전량 작물이 소비하는 것은 아니다. 용수원에서 취수하여 관개지구로 도수하는 동안 발생하는 수로손실, 논에서의 관리손실과 지하로 침투하는 침투손실 등과 같이 논에서 증발산으로 소비되지 않는 관개용수는 다소 시차는 있지만 결국은 상당부분 하천으로 회귀하게 된다. 이와 같이 관개용수 중에서 증발산에 의해 소비되지 않고 하천으로 다시 유입하는 용수를 회귀용수 혹은 회귀수량(Return flow)이라고 한다(Kim et al., 2010). 관개용수 회귀수량은 시설물에서 직접 공급되는 하천유지용수와 함께 하천유량의 증대와 갈수기 하천유지유량을 확보할 수 있는 물관리 방안이다. 또한, 수계 하천실무업무에서는 하천유지유량이 부족한 시기에 배수간선으로부터 하천으로 배수되는 회귀수량을 정량화하고 이를 이용하여 하천유량을 증대시키는 통합 유역 물관리 계획 수립을 목표로 하고 있다. 이 뿐만 아니라, 하천수허가업무상에 하천 가용수량을 산정하기 위해서 사용되는 회귀율은 획일적인 35%값을 이용하여 계산되고 있다. 따라서, 회귀수량의 정확한 추정은 유역의 용수공급계획, 하천유지유량 증대, 하천수허가 가용수량 증대, 하천유황의 예측, 관개용수의 사용량 결정, 하천수질관리 및 중소 농업유역의 수문모델링 등에 대단히 중요한 인자로 대표되지만 이에 대한 자료 및 연구가 미비하고 활발하지 못한 관계로 회귀수량의 영향을 고려하지 않거나 현실적이지 않은 회귀율을 사용하여 유역 물 수지 분석을 하는 경우가 대부분이다(Kim et al., 2022).
농업용수 회귀수량 산정에 대한 조사는 1990년 이전까지 미비한 실정이었으며, 최근에서야 회귀수량 산정에 대한 체계적인 연구가 일부 진행되었다. 국내에서 수행된 연구의 경우, Choi and Choi (2002)는 강원도 춘천시 신동, 천천 유역에 대해 용수로 유입부의 관개량과 하천수위 계측자료를 이용하여 논물수지 모델링을 통해 회귀율을 산정한 바 있으며, Noh et al. (2010)은 DAWST 모형을 이용하여 금강유역 내 5개 용수구역에 대한 농업용수 회귀율 35%로 산정한 바 있다. Kim et al. (2011)은 TANK 모형을 이용하여 하이저수지에 대한 회귀율을 추정하였고, Kim (2014)는 HOMWRS 모형과 저수지 유입부, 방류량, 관개량 및 배수량 계측자료를 이용하여 충북 진천군 내 무수저수지에 대한 신속회귀와 지연회귀율을 각각 28.4%, 24.3%로 도출하였다. Park and Kim (2016)은 BROOK90-K 모델을 이용하여 덕유산 내 구랑천에 대한 회귀율을 산정하였고, Kim et al. (2022)는 SWMM 모델을 이용하여 강원도 원주시 흥업저수지에 대한 신속회귀수량을 산정하였다.
국내에서 진행된 회귀수량 산정에 대한 대부분의 연구는 필지단위 또는 소규모 논 및 논 관개지구에서 물수지방법에 의해 이루어졌으며, 침투량을 모두 강하침투량으로 가정하여 연구를 진행하였다. 또한 일부 표본지역을 대상으로 한 단기 모니터링 결과는 전국단위의 유역별 대표 회귀율 산정에는 자료 등의 적용에 한계가 존재해왔다. 농업용수는 타 목적별 용수에 비해 사용량이 연중 고르지 못해 관개기에 주로 소요되는 용수이다. 농업용수의 필요수량은 기후, 토양, 관개방식, 용배수로, 양수, 지하수위등의 다양한 인자들의 지배를 받으며, 시기적, 지역적으로 변동성이 매우 커 정확하게 그 양을 측정하는 것은 매우 어렵기 때문에(Kang, 2016) 생활용수 회귀수량에 비해 연구사례가 많지 않다. 현재까지 조사된 회귀율을 살펴보면, 2000년 이전, 2000년 이후, 그리고 2010년대에 들어 변화하고 있다. 2000년대 이전의 회귀율을 살펴보면 32%-35%를 나타내었으며, 2000-2010년대는 41%-70%로 증가하였다. 2010년대에 들어서는 회귀율이 최대 84%까지 증가하는 추세이다(Kim et al., 2022).
반면 국외에서는 모델링을 이용한 회귀수량 산정 기법에 대한 연구가 꾸준히 진행되어왔다. Yalcin (2019) SWAT을 이용하여 터키 남동부 지역에 대한 밭관개 회귀수량을 모델링한 바 있으며, Wu et al. (2019)는 회귀수량의 재사용 규모 효과를 분석하기 위해 SWAT을 이용하여 중국 내 대규모 관개지역에 대한 관개량 operation이 포함된 수문모델링을 실시하였다. Dash et al. (2020)은 아시아의 논 물수지를 효과적으로 모의하기 위해 SWAT의 pothole을 이용한 담수논 모델링을 제안하였으며, Song et al. (2020)은 개념적 및 물리적 기반의 매개변수화 접근방식이 회귀수량의 불확실성에 미치는 영향을 분석하였으며, 제약적인 조건에서 정확도가 향상된 회귀수량 산정 기법을 제안하였다.
궁극적으로, 합리적이며 경제적인 수자원 이용 및 관리를 위해서는 수계별로 정확한 농업용수 회귀율을 실험을 통하여 측정하고 이를 물이용 계획에 반영해야 한다. 이를 위해서 유역단위에 대한 조사 분석과 모델링 기반의 유역단위 물수지 분석을 연계한 연구방법이 제시되어야 한다. 따라서, 본 연구에서는 유역 수문모델인 SWAT을 이용하여 소유역 내 담수 논에 대한 물 공급과 수문기작을 모의할 수 있는 모델링 기법을 개발하여 유역단위의 회귀수량을 정량화 해보고자 한다.
2. 연구대상지역
신속회귀수량 산정을 위한 연구대상지역은 유역을 대표할 수 있는 대표적인 농업지역을 선정하고, 가급적 단일 영농작업이 이뤄지는 중소규모의 유역으로 선정한다. 특히 논, 밭, 과수원, 주거지역 등 복합적인 토지이용이 이루어지는 유역에서는 용수의 이용패턴이 다르고 특히 강우 시 강우의 효과를 측정할 수 없기 때문에 정확한 회귀율을 산정할 수 없다. 또한 외부 유역으로부터 지표수 유입이나 지하수 용출이 없는 지구를 선정한다(Chung and Park, 2004). 이에 근거하여 본 연구의 대상지역은 영산강 유역 내 만봉천 표준유역(97.34 km2)을 연구 대상지역으로 선정하였다. 만봉천 표준유역은 유역 내 타 하천 유입이 없고, 나주호로부터 농업용수를 공급받고 있으며(Sung et al., 2021), 유역의 하류는 영산강 본류와 합류한다(Fig. 1).
3. SWAT 모형을 이용한 신속회귀수량 산정방법
3.1 SWAT에서의 관개량 적용
SWAT (Soil and Water Assessment Tool)은 미국 농무성 농업연구소(USDA Agricultural Research Service, ARS)에서 개발된 일 단위의 모의가 가능한 유역단위의 준 분포형 모형이다. 이 모형은 물리적 기반의 장기 강우-유출 모형으로서, 대규모의 복잡한 유역에서 장기간에 걸친 다양한 종류의 토양과 토지이용 및 토지관리 상태에 따른 유출의 영향을 예측하기 위해 개발되었으며 물수지 방정식에 근거를 두고 강수, 증발산, 지표유출, 기저유출, 지하수, 토양수분 등에 대한 모의를 각 수문반응단위(Hydrologic Response Units, HRU) 별로 계산할 수 있다(Arnold et al., 1996).
SWAT에서 관개는 사용자(Manual) 또는 자동(Auto)으로 계획된 관개일정에 따라서 적용되며, 관개모의는 .mgt의 Irrigation operation을 이용하여 모의·조작 할 수 있다. 관개 일정과 관개량을 구체적으로 적용하기 위해서는 관개용수의 수원공을 지정해야 하며, SWAT에서 적용 가능한 수원공의 형태는 하천취수(reach), 충적층(shallow aquifer), 심층(deep aquifer), 저수지(reservoir), 유역 외 공급(outside source)의 5개로 구분되어 있다. 하천취수(reach) 적용 시에는 최소 하천 유지유량과 최대 관개용수량을 조절하여 관개량을 계산하며, 충적층과 심층에서의 취수 방식은 물수지방정식을 통해 산정된 토양수분 및 지하수량에서 물을 취수되며, 저수지의 경우 SWAT 내 입력된 제원 및 방류량을 통해 산정된 담수량에서 물공급을 모의할 수 있다.
물공급 일정은 작물의 생육 시점(Plant / begin, growing season)과 수확 시기(Harvest and kill season)에 따라 조절이 가능하며, 작물 생육 모의를 위해서는 작물의 종류, LAI (Leaf area index) 초기값, 수확비율, 작물의 Curve Number (CN)에 따라 필요수량을 산정한다. 또한 SWAT은 관개 일정에 따라 수원공에서 공급 가능한 물의 양(관개량)을 산정하고, 관개량 모의를 통해 HRU에 적용된 물은 물수지방정식을 통해 지표, 중간, 기저유출을 계산한다. 만약 물공급량이 포장용수량을 초과하면 하천으로 회귀 된다.
3.2 SWAT에서 pothole을 이용한 담수 논 모의
SWAT에서 pothole로 지정된 HRU의 경우 관개용수는 토양 표층에 담수된 형태로 모의되며, 지정된 담수높이를 초과할 경우 하천으로 배수된다. SWAT에서 관개용수가 입력된 HRU pothole 담수상태와 수문성분 해석을 위한 기본구조는 물수지방정식을 기본으로 하여 선담수 후 손실량 산정의 방식으로 구성되어 있다. 먼저 HRU에 pothole 설정을 하여 pothole의 제원을 입력한 후, 해당 HRU에 대해 .mgt 모듈을 이용하여 물공급량이 입력된 후 HRU에 대한 수문성분 해석이 이뤄진다. 수문성분 해석은 HRU (output.hru) 및 pothole (.pot, .wtr) 단위로 산출된다. SWAT 모형 내 pothole 에서의 회귀수량은 다음 Eq. (1)을 이용하여 계산된다.
여기서, IRflow는 가능 관개회귀량(mm/day), Irrday는 적용한 관개량(mm/day), Pday는 강수량(mm/day), qday는 하천에 유입된 유출량(m3/day), lday는 하천에 유입된 측방흐름(m3/day), gwday는 하천에 유입된 기저 유출량(m3/day), 은 관개 효율이다.
Eq. (1)을 토대로 Pothole에서 하천으로 유출되는 양은 다음 Eq. (2)을 이용하여 계산된다.
여기서, qdqyout은 포트홀에서 하천으로 유입되는 총 유출량(mm/day)이다.
SWAT Pothole 모듈에서는 유입되는 강수량과 물 공급량 및 전날 담수심은 현재의 담수심을 결정하게 되는데 SWAT 소스코드 상에서는 전날 담수량을 고려하여 현재 부족한 담수심 만큼만 물 공급을 하고, 나머지는 하천으로 방류되도록 산정되도록 설계되어 있다. 이때 하천으로 방류되는 유량은 SPILLO라는 output file로 출력된다. 수로를 통해 공급되는 농업용수는 수혜구역에 무조건 공급이 되는 것이 아니라, 농민들이 필요한 수량 만큼 수로에서 취수하여 사용하는 형태로 물공급이 이뤄지고 있다. 또한 저수지의 말단부 까지 물 공급을 원활하게 하기위해, 수혜구역에서 필요한 수량 이외에 배분관리 용수량을 추가로 더 공급하는 것이 일반적이다. 이러한 경우 취수되지 않은 물 공급량은 하천으로 바로 배수되게 되며, SWAT에서는 이를 SPILLO 라는 변수로 결과값을 산정한다. 따라서 본 연구에서는 .mgt를 통해 입력된 IRR_AMT를 공급량으로 설정하고, Output.hru에서 산정되는 논 지역 내 증발산량(ET)과 삼투량(Percolation)을 논 사용량으로 가정하였다. 또한 Output.pot에서 출력되는 SPILLO를 공급량 중 취수되지 않고 바로 방류되어 하천으로 회귀되는 양으로 가정하여 다음과 같은 수식으로 회귀수량을 산정하였다. 또한 SPILLO의 단위 환산을 위해 수표면적 산정결과를 활용하였다.
여기서, Returnflow는 신속회귀수량(mm/day), IRRAMT는 물 공급량(mm/day), ET는 증발산량(mm/day), Percolation은 삼투량(mm/day), SPILLO는 공급량 중 취수되지 않고 방류되는 유량(mm/day)를 나타낸다.
IRRAMT의 단위는 mm/day로 본 연구에서는 Kim et al. (2007)이 적용한 IRRAMT 환산방식과 동일하게 물공급량을 수혜구역(HRU 면적)으로 나누어 입력하였다(Eq. (4)).
Fig. 2는 Pothole이 적용된 HRU에서의 수문성분 및 회귀수량 산정을 나타낸 것으로, 앞서 언급한 바와 같이 HRU에 입력된 pothole 비율, 최대, 최소 담수 높이(mm), pothole의 증발계수, Pothole의 투수계수를 변수로 입력한다. 이를 토대로 해당 HRU에 대한 .mgt에 작물 생육시작일, LAI, 용수공급원, 공급량, 관개효율 등의 관개정보가 입력된 상태에서 SWAT 모의를 실시한다. 해당 HRU에서는 기상인자와 .mgt를 통해 입력된 물공급량에 대한 물수지 계산이 이뤄지며, 담수된 pothole에서의 수면증발량과 논에서의 증발산량과 pothole에서 천층으로의 삼투량을 합하여 농지에서의 사용량으로 계산하고, 관개용수로 사용되지 않아 하천으로 배수되는 유량(배분관리 용수량 등)을 합하여 회귀수량을 산정한다.
HRU의 담수면적은 증발산량과 담수량 변화에 큰 영향을 주는 인자로, 기존의 SWAT에서는 pothole의 담수면적은 HRU의 평균경사를 이용하여 원뿔형으로 가정하여 일별로 산정하고 있다(Eq. (5)).
여기서, Surface area: 담수면적(ha), V: 담수량(m3), slp: HRU의 평균경사(m/m)이다.
그러나 원뿔형태는 우리나라 농경지의 대부분이 계단식의 형태를 띄고 있는 것과는 다소 차이가 있으므로 이를 직육면체로 가정하여 실제 HRU의 면적에 대한 관개지와 비관개지의 면적비를 이용하도록 다음 Eq. (6)과 같이 소스코드를 변경하였다.
여기서 areahru는 HRU의 면적(ha)이다.
3.3 자료수집 및 분석방법
3.3.1 GIS, 기상 입력자료 및 수문 검보정 자료
SWAT의 GIS 입력자료로 DEM (Digital Elevation Model)과 토양도는 국가수자원관리종합정보시스템(WAter Management Information System, WAMIS,)의 30 m×30 m DEM (Fig. 3(a)), 토지이용도는 환경부의 2020년 세분류 자료를 이용하였으며, 이를 9개 항목(활엽수림, 침엽수림, 혼효림, 주거지, 논, 밭, 초지, 나지, 수역)으로 분류하여 모형에 입력하였다(Fig. 3(b)). 만봉천 유역의 토지이용은 논이 27.1%, 밭이 20.5%를 차지하는 대표적인 농업유역이다. 토양도는 1:25,000 정밀토양도를 사용하였다(Fig. 3(c)). SWAT 모형의 검보정을 위한 기상자료는 4개(광주, 장흥, 강진군, 목포) 강수량(mm), 최고 및 최저기온(℃), 풍속(m/sec), 상대습도(%), 일사량(MJ/m2) 기상관측소의 14년(2005-2017)동안의 일자료를 사용하였다.
3.3.2 논 필지 HRU 생성을 위한 토양도 수정
SWAT은 DEM, 토지피복도, 토양도를 중첩하여 HRU를 생성하게 된다. SWAT에서의 관개모듈은 HRU 단위로 수행되며, 논으로 분류된 HRU를 기준으로 관개용수 수원결정을 위한 판단됨에 따라 HRU의 설정은 매우 중요하다. SWAT 모의 시 적용된 정밀토양도는 논 필지별 HRU 설정이 매우 어렵다. 이에 Kim et al. (2005, 2007)은 담수 논 모의를 위해 논에 해당하는 토양통을 다양한 별개의 HRU보다 대표 논 토양통으로 설정하여 모의하는 것이 관개지구를 벗어나는 지역에 대한 과잉 물 공급 모의를 방지할 수 있음을 분석한 바 있어, 본 연구에서도 동일한 방법으로 적용하였다. 논 필지 설정을 위해 기설수혜논 지역 중 저수지 수혜구역도를 기준으로 HRU을 생성하고자 하였다. 기설수혜논 지역의 구분은 농어촌용수 이용합리화계획에서 제시한 기설수혜논 지역을 기준으로 하였으며, 토양DB는 정밀토양도 자료를 이용하여 토양통을 조사하였다. ArcGIS를 이용하여 기설수혜 논 지역과 토양통을 중첩하였으며, 기설 수혜 논 지역 내 토양통은 수혜구역 내 가중평균을 통해 대표 토양통을 선정하여 토양도를 수정하였다.
3.3.3 물공급량 입력자료
SWAT에는 용수공급 시작·종료일, 공급일수, 시설물 효율 등의 정보가 입력됨에 따라, 유역 내 시설물 운영패턴 분석이 필요하다. 이를 위해 만봉천 유역 내 저수지 수리시설물의 운영자료를 농어촌공사로부터 제공받아 분석을 수행하였다. 유역 내 농업용 저수지는 총 20개소로 조사되었으나, 이중 일별 물공급량 모니터링자료를 제공하는 저수지는 6개소(곡치(GH), 입석(IS), 토동1(TD1), 지초(JC), 식산(SS), 연보(YB))로 조사되었다. 또한, 만봉천 유역 외에 위치한 나주호로부터 농업용수를 공급받음에 따라 나주호의 일 공급량 자료를 합계에 포함시켰다. Sung et al. (2021)은 나주호의 간선수로별 용수공급량을 분석하였으며, 운곡(UG)간선 수로를 이용하여 만봉천 유역에 나주호 용수량의 11%의 용수를 공급하고 있음을 확인하였다. 따라서 본 연구에서도 나주호의 용수공급량 중 11%에 대한 용수량을 저수지를 통한 용수공급량에 포함하였다.
2012년부터 2017년까지 분석기간 동안 최소 9,096.9 천톤에서 최대 20,043.2 천톤까지 농업용수를 공급한 것으로 분석되었다(Table 1). 다만 2015년 이전 입석(IS)과 연보(YB) 저수지의 둑높임 사업으로 인해 물공급량이 크게 차이가 발생하여, 본 연구에서는 2015~2017년 자료만 분석에 활용하였다.
Table 1.
Irrigation amount from agricultural reservoir (Unit: 1,000 m3)
4. 회귀수량 산정결과
4.1 SWAT의 검·보정
본 연구에서는 회귀수량 산정에 앞서 모형 보정을 위해, 대상지구 내 나주 유량관측소 지점의 유량자료를 이용하여 모형의 검증 및 보정을 수행하였다. 매개변수는 만봉천 유역에 대한 물공급량을 모두 입력한 후 조정되었으며, 검증과 보정에 사용한 매개변수는 민감도 분석결과(Sensitivity analysis)와 선행연구(Kim et al., 2020)를 참고하여 매개변수를 선정하였다. 선정된 매개변수 중 ESCO는 증발 보상계수로 증발산량을 조정하는 매개변수이다. LAT_TIME은 중간유출 흐름 시간이며, GW_DELAY와 GWQMN은 지하수 관련 매개변수로 각각 지하수 흐름시간, 기저유출이 발생하는 충적층 내 수위를 조정하는 매개변수이다.
또한, 농업용수 회귀수량 산정을 위해 작물의 생육 시점, 수확 시기, 그리고 물 공급량을 입력하였으며, 입력 변수는 Table 2와 같다. RDA (2018)에 따르면, 우리나라 남부지방의 이앙시기는 평균 5월 25일에 시작하여 수확은 9월 25일 전후로 이뤄진다. 따라서 본 과업에서는 생육시작일을 5월 25일로 설정하였으며, 수확은 9월 25일로 설정하였다. 작물 생육시작 시 작물의 CN 값은 SWAT DB에서 사용하는 73을 적용하였고, 초기 엽면적 지수는 선행연구를 참고하여 1로 입력하였다(Ha et al., 2008; Park et al., 2012). LAI는 5월 25일부터 9월 25일까지 1부터 10 사이의 값으로 누적되도록 설정되었다. 작물의 수확 시 CN값은 생육 시작 시 CN값과 동일하게 73을 적용하였다.
Table 2.
Calibration parameter for irrigation operation
| Parameter | Definition | Range / Reference | Adjusted | |
| Evapotranspiration | ESCO | Soil evaporation compensation coefficient | 0~1 | 0.2 |
| Lateral flow | LAT_TIME | Lateral flow travel time (days) | 0~180 | 5 |
| Groundwater flow | GW_DELAY | Delay time for aquifer recharge (days) | 0~500 | 40 |
| GWQMN |
Threshold water level in shallow aquifer for base flow (mm) | 0~5,000 | 500 | |
|
Planting / beginning of growing season | Planting date | Planting start date | RDA (2018) | 5/25 |
| CNOP | Crop SCS CN | SWAT DB | 73 | |
| LAI initial | Initial leaf area index | Park et al. (2012), Ha et al. (2008) | 1 | |
| Irrigation operation | Start date | Irrigation start date | Observed data | 4/1 |
| End date | Irrigation end date | Observed data | 9/29 | |
| IRR_AMT | Depth of irrigation water applied on HRU | Observed data | Table 1 | |
| IRR_EFM | Irrigation efficiency (fraction) |
Conveyance loss in canal (concrete canal 5~7% ) | 0.93 | |
| IRR_SQ | Surface runoff ratio (fraction) | Sakaguchi et al. (2014) | 1 | |
| Harvest and kill operation | Harvest date | Harvest date | NICS | 9/25 |
물공급 시작일과 종료일은 모니터링 자료와 동일하게 각각 4월 1일과 9월 29일을 입력하였으며, 송수손실율(IRR_ EFM)은 콘크리트 수로손실율을 고려하여 0.93으로 입력하였다. 공급용수 중 지표유출의 비율인 IRR_SQ는 선행연구를 토대로 1로 입력하였다. Sakaguchi et al. (2014)은 SWAT 모형에서 Pothole을 이용하여 담수 논 모의 시, pothole은 배수로 없이 바로 하천으로 유입된 것으로 모의 되기 때문에 공급용수는 용수 사용 후 모두 지표유출로 반영됨에 따라 1로 설정하였다. 물공급량 입력 변수인 IRR_AMT는 Table 1과 같이, 수집한 2015년부터 2017년까지의 일별 물공급량 실적자료를 분석하여 일단위로 입력하였다.
본 연구에서는 담수 논 모의를 위해 SWAT의 Pothole 모듈을 이용하였다. 변수 입력을 위해 선행연구의 결과를 조사하였으며, 이를 토대로 매개변수를 설정하였다(Table 3). Pothole의 제원 관련 매개변수 3개(POT_FR, POT_VOLX, POT_VOL)와 pothole의 수문반응 관련 매개변수 2개(EVPOT, POT_K)를 선정하였다. HRU에 대한 Pothole 면적비를 나타내는 매개변수인 POT_FR은 HRU 생성 시 농어촌용수 이용합리화계획의 기설수혜논 지역을 고려하였기 때문에 입력변수의 값을 1로 설정하였다. 또한 담수논의 최대담수심을 의미하는 POT_VOLX 매개변수는 Kim et al. (2007)이 적용한 40 mm를 동일하게 적용하였으며, 관개용수 공급 시작일의 초기 pothole 담수심높이(POT_VOL)는 선행연구와 동일하게 0으로 적용하였다.
Table 3.
Pothole input parameters in HRU
| Parameter | Definition | Reference | Adjusted | ||
| Kim et al. (2007) | Sakaguchi et al. (2014) | Wu et al. (2019) | |||
| POT_FR |
Fraction of HRU area that drains into the pothole | 0.75 | 1 | 0.95 | 1 |
| POT_VOLX |
Maximum volume of water stored in the pothole (mm) | 40 | 30 | 60 | 40 |
| POT_VOL |
Initial volume of water stored in the pothole (mm) | 0 | 0 | 0 | 0 |
| EVPOT | Pothole evaporation coefficient | - | - | - | 2 |
| POT_K |
Hydraulic conductivity of soil surface of pothole (mm/hr) | - | - | - | 1 |
Pothole에서의 증발계수(EVPOT)와 토양수리전도도(POT_K)는 각각 0~2, 0.02~10 mm/hr의 범위 내에서 입력 가능하며, 두 매개변수는 선행연구사례가 존재하지 않음에 따라 매개변수 민감도 분석을 수행하였다. 매개변수 민감도 분석 결과, EVPOT는 0에서 2로 증가할 수록 증발산량이 1.7% 감소하는 경향을 나타내었고, POT_K는 0.02에서 10으로 증가할수록 삼투량이 최대 23.3% 감소하는 것으로 나타났다. 특히 0.02에서 1로 변화시 가장 민감하게 변화하는 것으로 분석되었다. 따라서 본 연구에서는 민감도 분석결과를 바탕으로 EVPOT와 POT_K를 각각 2와 1로 입력하였다.
SWAT의 검보정은 물공급량 입력자료의 기간을 고려하여 보정기간은 2년(2015-2016), 검증기간은 1년(2017)으로 설정하였다. 모형의 적용성 평가는 coefficient of determination (R2), Nash and Sutcliffe (1970) model efficiency (NSE), root-mean square error (RMSE), percent bias (PBIAS)를 사용했다. R2은 1에 가까울수록 관측 값과 모의 값이 완벽하게 일치함을 나타내고, NSE는 -∞에서 1까지의 범위를 통해 모형의 효율을 평가하는데 0보다 클수록 모형의 적용성이 높다고 설명할 수 있다(Mkhwanazi et al., 2012). RMSE는 실측치와 모의치의 오차를 나타내는 것으로 0에 가까울수록 오차가 작음을 의미한다.
Table 4와 Fig. 4는 SWAT 나주(NJ) 유량관측소에 대한 유량 검·보정결과를 나타낸 것으로 전체 기간 동안의 R2, NSE, RMSE는 각각 0.76, 0.69 0.50 mm/day로 분석되었다. 또한 관개기간의 R2, NSE, RMSE는 각각 0.73, 0.70, 0.64 mm/day로 분석되었다. NSE의 경우 전체기간과 관개기간의 효율은 큰 차이는 없었으나, R2의 정확도는 0.03 감소하였으며, RMSE는 0.14 mm/day 증가하였다. Gassman et al. (2014)은 SWAT 검보정 시 NSE ≥ 0.5, R2 ≥ 0.6 이상일 경우 통계적으로 유의하다는 결과를 제시한 바 있다. 선행연구를 토대로 본 연구에서의 검보정 결과는 통계적으로 유의한 것으로 분석되었다.
Table 4.
Calibration and validation results for streamflow (NJ)
4.2 신속회귀수량 산정 결과
본 연구에서는 2015~2017년까지의 모의 결과를 토대로 물공급 기간(5/1~9/30)에 대한 회귀수량과 공급량 대비 신속회귀율을 산정하였다. 결과분석에 사용된 변수는 재료 및 방법에서 제시하였던 저수지 물공급량을 공급량(Qi)으로 설정하고, 논 지역 내 증발산량(ET)과 삼투량(Percolation)을 논 사용량(Qu)으로 가정하였으며 하천으로 배수되는 유출량을 취수되지 않고 바로 방류되어 하천으로 신속회귀되는 양(Qo)으로 가정하였다.
Table 5는 월별 신속회귀수량 분석결과를 나타낸 것으로, 여기서 Return flow 1은 저수지 물공급량(Qi)에서 논 사용량(Qu)을 뺀 회귀량이며, Return flow2는 취수되지 않고 바로 하천으로 방류되는 신속회귀량(Qo) 모의 결과이다. 4월의 물공급량은 전체적으로 낮아 회귀율을 산정하지 않았고, 본격적으로 이앙기가 시작되는 5월부터 신속회귀수량을 산정하였다. 분석기간(2015~2017) 동안의 산정된 평균 신속회귀량은 494.9 mm, 회귀율은 53.4%로 분석되었으며, 2016년이 60.2%로 회귀율이 가장 높았다. 신속회귀율이 가장 낮은 연도는 2017년으로 회귀율이 42.8% 였다. 월별로 회귀율을 분석한 결과, 8월의 회귀율이 68.0~83.9%로 가장 높았으며 5월의 회귀율이 36.1~46.7%로 가장 낮게 나타났다.
그러나 본 연구의 대상 지역인 만봉천 유역에 대한 신속회귀수량 조사가 이뤄지지 않음에 따라 영산강 유역 내 다른 표본조사 지점에 대한 회귀수량조사 자료와의 비교를 통해 도출된 회귀율의 적정성을 평가하였다. 1997년 수자원 계획의 최적화 연구에서 조사된 영산강 유역 저수지와 양수장 및 보에 대한 신속회귀율은 각각 51.7%, 23%, 50%로 나타났으며(Chung et al., 1997), 2006년 영산강·섬진강 유역조사 중 영산강 유역 내 학야지구에 대한 회귀율은 22.5%로 조사된 바 있으나 이는 저수지에 대한 신속회귀율이 아닌 학야지구 내 필지에 관한 회귀율 산정 결과이다(K-water, 2006). 본 연구에서 산정된 회귀율은 신속회귀량으로, 1997년 조사된 영산강 유역 저수지에 대한 회귀율인 51.7%와 유사하게 나타난 것으로 분석되었다.
Table 5.
Results of water supply and quick return flow from irrigation canal during 2015-2017 (unit: mm)
5. 결 론
국내에서 농업용수 회귀수량에 대한 연구 및 조사는 주로 필지단위에 대하여 이뤄져 왔으며 용수공급현황과 물관리 방식에 따라 다양하게 산정 및 추정되어왔다(Song et al., 2015; Kim et al., 2021). 유역 중심의 통합물관리의 목적을 달성하기 위해서는 수계별로 농업용수의 합리적인 사용에 대한 중요성이 강조되고 있으며 특히, 정확한 농업용수 회귀율의 산정은 하천수 허가업무와 장래 물이용 계획 등 다양한 계획에 반영되어야 하는 매우 중요한 요소이다. 이에 본 연구에서는 유역 모델링을 이용한 유역별 회귀율 산정을 위한 방법을 모색하기 위해 물리적기반의 유역수문 모형을 이용한 유역단위 회귀수량 산정 방법론을 제안하고 이를 대표적인 농업유역에 적용하여 회귀율을 산정하였다.
SWAT에서 논 관개는 Irrigation operation을 이용하여 모의·조작 하였으며, 물공급 일정은 작물의 생육 시점(Plant/begin, growing season)과 수확 시기(Harvest and kill season)에 따라 필요수량을 산정하고 실제 저수지로부터 공급된 물공급량을 적용하였다. 우리나라 담수논 형태의 모의를 위해 토양통을 논 수혜구역 단위로 수정하여 하나의 HRU로 만든 후, pothole로 지정하였다. 기존의 SWAT에서는 pothole의 담수면적은 HRU의 평균경사를 이용하여 원뿔형으로 가정하여 일별로 산정됨에 따라 본 연구에서는 직육면체의 형태로 계산되도록 소스코드를 일부 수정하였다.
회귀수량 산정에 앞서 모형 보정을 위해, 대상지구 내 나주 유량관측소 지점의 유량자료를 이용하여 모형의 검증 및 보정을 수행하였으며, 보정결과 R2, NSE, RMSE는 각각 0.76, 0.69, 0.50 mm/day로 분석되어 통계적으로 유의한 수준으로 나타났으며, 2015년부터 2017년까지의 회귀수량 산정 결과 회귀율은 평균 53.4%로 분석되었다. 산정된 회귀율은 1997년 조사된 영산강 유역 저수지에 대한 회귀율인 51.7%와 유사하게 나타난 것으로 분석되었다.
실제 하천수허가업무에 있어서 물수지 분석을 통한 가용수량 산정시 사용되는 농업용수 회귀율 값은 35%로 획일적으로 사용하고 있다(MOLTM, 2011). 본 연구결과를 활용하여 회귀율 값을 지역별로 사용할수 있다면 현재 가용수량이 없는 지역에서의 물 사용 가능성을 증대 시킬수 있으며 실제 물관리 업무에 효율성을 개선 시킬 것으로 기대된다. 다만, 본 연구결과는 SWAT을 이용한 신속회귀량 산정 가능성을 검토하였고, 회귀수량 산정에 있어 가장 필수적인 요소에 대해서만 모델이 개선되었다. 즉, 신속회귀수량에 대한 직접유출량이 기여할 수 있는 수량을 정량화한 연구로 지하수로 유출되는 지연회귀에 대한 추가 모델링 방법이 정립되어야 실무에서 최종적으로 활용 될 수 있다. 따라서, 향후 연구로 지연회귀를 포함한 총 회귀수량에 대한 방법 및 조사를 통한 관측 회귀수량과의 검증이 반드시 이루어져야 하며 용수의 재이용 관련 자료, 취입보, 양수장 등의 운영자료가 보완된다면, 보다 신뢰성 있는 회귀수량이 산정될 것으로 기대된다. 또한 이를 활용할 수 있는 농업용수 최적 운영방안에 대해서 제시되어야 할 것으로 판단된다.
