본 연구에서는 낙동강 유역을 대상으로 유역-하천 흐름 연계 모형을 구축하였다. 본 연구는 기후변화에 따른 하천 흐름 및 수질 영향 분석을 위한 전단계의 연구이다. 낙동강 유역을 대상으로 유역 수문 모형인 SWAT 모형을 구축하고 낙동강 본류 EFDC 모형에 유량 경계자료를 제공하여 본류 흐름을 예측하였다. 유역 모의를 통해 계산된 유출량을 본류 13개, 지류 30개 지점에서 2004년부터 2009년까지의 환경부8일 자료를 이용하여 검보정하였다. 모형의 적용성을 평가하기 위하여 % difference, NSE, $R^2$를 산정하였으며, % difference는 몇몇 지점을 제외하고는 대부분 15% 이내인 것으로 나타났고, NSE와$R^2$ 역시 대부분 적정 수준인 것으로 모의 되었다. 한편, SWAT 모형으로부터 본류의 상류 및 지류 경계의 유출량 값을 제공받아 EFDC모형에 입력하였고 낙동강 본류에서의 2007년도 흐름을 예측하였다. 계산된 결과를 낙동강 본류의 각 수위 관측소에서의 측정 데이터와 비교한 결과, 본 모형이 저수기 및 홍수기시 모두 잘 예측하는 것으로 나타났다. 또한 만곡부 및 하중도 부근에서의 유속 증가와 감소 현상 등을 잘 모의하는 것으로 확인되었다. 한편, 상류 경계로 부터 1월1일 유입된 염료의 이동현상을 모의한 결과 하구언까지 도달하는데 약 65일이 소요되는 것으로 나타났다.
본 연구에서는 낙동강 유역을 대상으로 유역-하천 흐름 연계 모형을 구축하였다. 본 연구는 기후변화에 따른 하천 흐름 및 수질 영향 분석을 위한 전단계의 연구이다. 낙동강 유역을 대상으로 유역 수문 모형인 SWAT 모형을 구축하고 낙동강 본류 EFDC 모형에 유량 경계자료를 제공하여 본류 흐름을 예측하였다. 유역 모의를 통해 계산된 유출량을 본류 13개, 지류 30개 지점에서 2004년부터 2009년까지의 환경부8일 자료를 이용하여 검보정하였다. 모형의 적용성을 평가하기 위하여 % difference, NSE, $R^2$를 산정하였으며, % difference는 몇몇 지점을 제외하고는 대부분 15% 이내인 것으로 나타났고, NSE와$R^2$ 역시 대부분 적정 수준인 것으로 모의 되었다. 한편, SWAT 모형으로부터 본류의 상류 및 지류 경계의 유출량 값을 제공받아 EFDC모형에 입력하였고 낙동강 본류에서의 2007년도 흐름을 예측하였다. 계산된 결과를 낙동강 본류의 각 수위 관측소에서의 측정 데이터와 비교한 결과, 본 모형이 저수기 및 홍수기시 모두 잘 예측하는 것으로 나타났다. 또한 만곡부 및 하중도 부근에서의 유속 증가와 감소 현상 등을 잘 모의하는 것으로 확인되었다. 한편, 상류 경계로 부터 1월1일 유입된 염료의 이동현상을 모의한 결과 하구언까지 도달하는데 약 65일이 소요되는 것으로 나타났다.
In this study, the watershed and water body models were linked for the simulation of the Nakding river flow. This is a pre-step study for the estimation of the effect of the flow and water quality on the climate change. For models of watershed and river flow, the SWAT and EFDC were used, respectivel...
In this study, the watershed and water body models were linked for the simulation of the Nakding river flow. This is a pre-step study for the estimation of the effect of the flow and water quality on the climate change. For models of watershed and river flow, the SWAT and EFDC were used, respectively. The runoff discharge at each boundary points for the simulation of the river flow was provided from the drainage basin model. The calculated runoff discharge by the SWAT model was compared with the measured data of the Ministry of Environment at 13 locations along the Nakdong river and 30 locations along the tributary streams. The computed water discharge was shown to be similar with the measured data. For the model calibration and verification, % difference, NSE, and $R^2$ were computed. The computed % difference was within 15% except of a few points. The NSE and $R^2$ were also within a fair level. The Nakdong river flow of 2007 was simulated by using the EFDC model. The comparison with the measured data showed that the model reflected the actual values of low and high flow well. Also, it was confirmed that the acceleration and deceleration in the curved areas were appropriately simulated. The movement of dye injected at the upstream boundary was simulated. The result showed that the arrival time up to the estuary dam was computed to be about 65 days.
In this study, the watershed and water body models were linked for the simulation of the Nakding river flow. This is a pre-step study for the estimation of the effect of the flow and water quality on the climate change. For models of watershed and river flow, the SWAT and EFDC were used, respectively. The runoff discharge at each boundary points for the simulation of the river flow was provided from the drainage basin model. The calculated runoff discharge by the SWAT model was compared with the measured data of the Ministry of Environment at 13 locations along the Nakdong river and 30 locations along the tributary streams. The computed water discharge was shown to be similar with the measured data. For the model calibration and verification, % difference, NSE, and $R^2$ were computed. The computed % difference was within 15% except of a few points. The NSE and $R^2$ were also within a fair level. The Nakdong river flow of 2007 was simulated by using the EFDC model. The comparison with the measured data showed that the model reflected the actual values of low and high flow well. Also, it was confirmed that the acceleration and deceleration in the curved areas were appropriately simulated. The movement of dye injected at the upstream boundary was simulated. The result showed that the arrival time up to the estuary dam was computed to be about 65 days.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
그러나 각 년도별 최대 유량을 본 모형이 과대 산정하는 것을 볼 수 있는데, 이는 실측값이 8일 간격 유량이 라는 점과 홍수기 보다는 평수기 중심으로 유량을 관측하였기 때문인 것으로 판단된다. 따라서 총량 8일 간격 관측 자료의 홍수기 유량한계를 보완하기 위하여 국토해양부의 주요 본류 및 지류지점의 일간격 유량 관측자료를 이용하여 홍수기 유량에 대한 보정을 수행하였다. Fig.
본 연구는 낙동강 유역을 대상으로 유역 유출 모형과 하천에서의 흐름모형 연계를 통한 본류 흐름 해석을 수행하였다. 사용된 유역과 하천 흐름 모형으로는 각각 SWAT과 EFDC 모형을 적용하였다.
본 연구는 대규모 유역에서의 향후 기후변화에 의한 하천 흐름과 수질 변화를 예측하기 위한 전단계의 연구이다. 본 연구에서는 유출-흐름 연계 모형 구축을 목적으로, 낙동강 유역을 대상으로 모형 구축과 모형 검·보정을 수행하였다.
본 연구에서는 유출-흐름 연계 모형 구축을 목적으로, 낙동강 유역을 대상으로 모형 구축과 모형 검·보정을 수행하였다.
제안 방법
SWAT 모형의 모의기간은 총 8년 (2002∼2009)이며, 모의기간의 안정화를 위해 2002∼2003년을 추가 모의하였다.
먼저 SWAT 모형을 구축하고 2004년-2009년까지 낙동강 본류 및 지류 43개 지점에서의 모형 검·보정을 수행하였다. 또한 낙동강 본류에 대해 EFDC 모형을 구축하고, 상류 경계와 각 12개 주요 지류 경계에서의 유출량은 SWAT 모형의 모의결과를 이용하여 2007년도 낙동강 본류의 각 수위표에서의 측정된 수위값과 비교하였다. 한편, 상류경계 로부터 유입된 염료의 이동을 모의하여 하구언까지의 도달시간을 산정하였다.
건설교통부 (2000)를 참고하여 약 35%가 하천으로 회귀하는 것으로 하였다. 또한 생공용수량과 농업용수량을 합산하여 산정한 각각의 소유역별 취수량자료를 유역모형에 입력하였다.
그림에서 지류 위치는 앞의 유역 모델링과 동일하다. 또한 유역 모의를 통한 2007년 일평균 유량 결과를 12개의 각 지류 및 상류 경계에서의 유량 경계조건으로 사용하였다. Fig.
한편, 본 연구의 대상유역인 낙동강 유역에 대한 연구사례로는 신현석 등 (2007)이 유역내 댐방류 및 인위 적인 물이용 요소를 고려한 물수지 시나리오에 따른 하천 유량에 미치는 영향을 정량적으로 분석하였다. 또한, 낙동강 유역에서의 기후변화에 대한 영향 분석으로는 김문성등 (2009)이 캐나다에서 개발된 CGCM 3.1 T63의 A2 기후변화 시나리오를 대상으로 한 미래기후정보를 MSPG에 의해 지역 상세화시킨 기후자료로 변환한 후 수정 TANK 강우-유출모형에 적용시켜 기후변화가 낙동강 본류 유량에 미치는 영향에 대하여 분석하였다. 손태석 등 (2010)은 낙동강 유역을 대상으로 기후변화로 인한 하천 유량의 영향을 분석한 결과, 하천 및 호소의 수온상승과 유출량 증대가 예상되며 이에 따른 수리 구조물의 치수능력 증대 방안에 대한 필요성을 제시한 바 있다.
1은 유역모형 구축을 위한 소유역분할, 수치고도모델 (DEM), 토지이용도, 토양 분포도를 보여준다. 먼저 DEM은 국가수자원관리 종합정보 (resolution 30 m)에서 제공받았으며, SWAT 모형의 소유역 자동분할 기능을 사용하여 하천도 및 대상유역의 경계를 추출하였다. 오염총량관리 단위유역과 댐 방류지점 등을 기준으로 하여 총 67개의 소유역을 추출하였다.
먼저 SWAT 모형을 구축하고 2004년-2009년까지 낙동강 본류 및 지류 43개 지점에서의 모형 검·보정을 수행하였다.
사용된 유역과 하천 흐름 모형으로는 각각 SWAT과 EFDC 모형을 적용하였다. 먼저 SWAT 모형을 이용하여 2004~2009년도의 유출해석을 수행하고 낙동강 본류와 지류 43개 지점에서의 측정값과 비교하였다. 한편, 구축된 유출모형으로부터 낙동강 본류 경계와 지류 경계에서의 계산된 유출량을 EFDC 모형의 경계유량으로 입력하여 계산하고, 실측 수위와 비교하였다.
매개변수 보정 방법에는 크게 시행착오법과 최적화기법에 의한 자동보정법이 있다. 본 연구에서 사용한 자동보정법의 경우 한 개의 보정지점 당 계산시간이 상당히 길기 때문에 시행착오를 최소화하기 위하여 국내외 연구자들이 보고한 매개변수를 참고하여 유량모의에 적절한 매개변수를 선정하여 허용범위 내에서 변화시키면서 단순시행착오법에 의해 보정하였다 (Table 1). 보정은 본류구간 13개 지점 (낙본 A∼M)과 지류구간 30개 지점 총 43개 지점에 대하여 각각 수행되었으며, 각각의 위치는 Fig.
유역경사도 범주는 일반적으로 이용되는 0∼7%, 7∼15%, 15∼30%, 30∼60%, 60%이상 등 5개로 구분하였으며, 수문반응단위 (HRU: Hydrologic Response Unit) 생성을 위해 소유역별 토양특성과 토지이용특성을 중첩하고 토지이용, 토양면적, 유역경사도 비율을 각각 10%로 설정하여 총 1,256개의 HRUs를 생성하였다.
한편, 유역 모형과 하천 및 호소 모형을 결합하여 흐름과 수질을 예측한 연구로는 정제호 (2010), 박재충 등 (2010), 정광욱 (2007) 등의 연구가 있다. 이들 연구는 대부분 댐유역에서의 수질변화를 분석하는 것이 목적이었고, 사용된 유역모형으로는 SWAT과 HSPF 모형을, 호소내에서의 수질예측은 EFDC와 WASP 모형을 이용하였다. 이상의 최근 연구동향을 살펴보면, 우선, 낙동강과 같은 대규모 유역에서의 유역 모형 및 하천 모형이 구축된 사례가 많지 않으며, 유역-하천모형의 연계를 통한 해석 역시 많지 않은 것을 볼 수 있다.
또한 낙동강 본류에 대해 EFDC 모형을 구축하고, 상류 경계와 각 12개 주요 지류 경계에서의 유출량은 SWAT 모형의 모의결과를 이용하여 2007년도 낙동강 본류의 각 수위표에서의 측정된 수위값과 비교하였다. 한편, 상류경계 로부터 유입된 염료의 이동을 모의하여 하구언까지의 도달시간을 산정하였다.
대상 데이터
SWAT 모형에서 이용 가능한 방류량 자료는 크게 환경기초시설과 댐으로 구분되며 환경기초시설로부터 방류되는 유량은 2007년 전국 오염원 기초자료 (환경부, 2008a)를 이용하였다. 각각의 소유역에 위치한 환경기초시설의 방류량을 산정하였고, 댐 방류량은 한국수자원공사의 각 댐별 제원과 일방류량 (1994∼2009)을 수집하여 관측 일유량법 (observed daily outflow method)으로 모델에 적용하였다.
5(b)의 각 수위표 위치는 본 연구에서 계산된 수위 결과를 비교하기 위한 지점 수위표이다. 각 수위표 지점에서의 관측 수위는 한국수문조사연보 (국토해양부, 2007)의 자료를 이용하였다. 또한 하류단 경계지점인 하구언에서의 일평균 수위는 한국수문조사연보의 자료를 이용하였으며, 하구언에서의 연평균수위는 0.
강우 데이터 구축을 위해 국토해양부 산하 35개 지점과 기상청 산하 21개 지점 등 총 56개 지점에서의 일강우량 자료를 수집하였다. 또한 강우량을 제외한 온도, 풍속, 일조량, 상대습도 등의 기상자료는 낙동강 유역 내의 합천, 밀양, 구미, 남해, 대구, 마산, 문경, 부산, 산청, 안동, 영덕, 영주, 영천, 울산, 의성, 진주, 추풍, 통영, 포항, 거창 그리고 춘양(봉화)의 총 21개소의 기상관측소에 대한 자료를 수집하여 입력 자료를 구축하였다.
강우 데이터 구축을 위해 국토해양부 산하 35개 지점과 기상청 산하 21개 지점 등 총 56개 지점에서의 일강우량 자료를 수집하였다. 또한 강우량을 제외한 온도, 풍속, 일조량, 상대습도 등의 기상자료는 낙동강 유역 내의 합천, 밀양, 구미, 남해, 대구, 마산, 문경, 부산, 산청, 안동, 영덕, 영주, 영천, 울산, 의성, 진주, 추풍, 통영, 포항, 거창 그리고 춘양(봉화)의 총 21개소의 기상관측소에 대한 자료를 수집하여 입력 자료를 구축하였다.
또한 하류단 경계지점인 하구언에서의 일평균 수위는 한국수문조사연보의 자료를 이용하였으며, 하구언에서의 연평균수위는 0.86 m이고, 연 최고 수위는 7월 22일 3.09 m을 기록하였다.
SWAT 모형의 모의기간은 총 8년 (2002∼2009)이며, 모의기간의 안정화를 위해 2002∼2003년을 추가 모의하였다. 모형의 보정은 2007년부터 2009년까지, 검증은 2004년부터 2006년까지 수행되었으며, 낙동강 유역의 총량 8일 간격 관측자료를 사용하여 연중 일별로 이루어졌다. 매개변수 보정 방법에는 크게 시행착오법과 최적화기법에 의한 자동보정법이 있다.
보정은 본류구간 13개 지점 (낙본 A∼M)과 지류구간 30개 지점 총 43개 지점에 대하여 각각 수행되었으며, 각각의 위치는 Fig. 2와 같다.
먼저 DEM은 국가수자원관리 종합정보 (resolution 30 m)에서 제공받았으며, SWAT 모형의 소유역 자동분할 기능을 사용하여 하천도 및 대상유역의 경계를 추출하였다. 오염총량관리 단위유역과 댐 방류지점 등을 기준으로 하여 총 67개의 소유역을 추출하였다. 토지이용도는 환경부 환경지리정보서비스 (http://egis.
토양특성 등의 정보는 국립식량과학원에서 제공하는 1:25,000 해상도의 정밀토양도를 이용하여 모형을 구축하였으며, 토양의 물리적 특성 및 화학적 특성은 한국토양정보시스템 (http://asis.rda.go.kr)의 자료를 사용하였고, 토양층의 유효수분량, 포화수리전도도 등의 물리적 값들은 한국건설기술연구원에서 제공하는 HyGIS Soil DB에서 추출하여 사용하였다. 낙동강 유역에서는 삼암 (17.
오염총량관리 단위유역과 댐 방류지점 등을 기준으로 하여 총 67개의 소유역을 추출하였다. 토지이용도는 환경부 환경지리정보서비스 (http://egis.me.go.kr/egis)에서 제공하는 2007년 기준으로 작성된 1:25,000 해상도의 중분류 토지 피복도를 사용하였다. 한편, 낙동강유역의 토지이용형태를 살펴보면 산림 69.
하천 단면 데이터는 낙동강수계 하천기본계획 (변경) 보고서 (부산지방국토관리청, 2009a, b)의 하상고 자료를 이용하였다. 계산 구간은 낙동강 하구언 (St.
023인 것으로 나타났다. EFDC 모의를 위한 초기 수심 자료는 2007년 1월 1일 상류단 유량-하류단 수위 조건에 대한 CCHE 모형 (Zhang and Jia, 2009)의 2차원 흐름 모의를 수행하고, 그 결과를 EFDC 모형의 초기 조건으로 이용하였다.
각각의 소유역에 위치한 환경기초시설의 방류량을 산정하였고, 댐 방류량은 한국수자원공사의 각 댐별 제원과 일방류량 (1994∼2009)을 수집하여 관측 일유량법 (observed daily outflow method)으로 모델에 적용하였다.
낙동강 구간에서의 흐름 모델링을 위하여 EFDC (Envi-ronmental fluid Dynamics Code) 모형을 이용하였다. EFDC 모델은 Virginia Institute of Marine Science의 Hamrick에 의해 1990년대 초에 개발되었고, 현재는 USEPA 와 Tetra Tech Inc.
본 연구는 낙동강 유역을 대상으로 유역 유출 모형과 하천에서의 흐름모형 연계를 통한 본류 흐름 해석을 수행하였다. 사용된 유역과 하천 흐름 모형으로는 각각 SWAT과 EFDC 모형을 적용하였다. 먼저 SWAT 모형을 이용하여 2004~2009년도의 유출해석을 수행하고 낙동강 본류와 지류 43개 지점에서의 측정값과 비교하였다.
본 모형은 수직 방향으로 σ좌표계를 사용하고, 수평 방향으로는직선 또는 직교곡선좌표계를 사용하고 있다. 수직 와동점성계수의 산정에는 Mellor and Yamada (1982)의 level 2.5 turbulence closure scheme을 사용하고, 수평와동 점성계수는 Smagorinsky type diffusivity를 사용한다. 또한 Mass conservation scheme을 사용한 마름/젖음 구간 처리가 가능하다.
본 연구에서는 유출-흐름 연계 모형 구축을 목적으로, 낙동강 유역을 대상으로 모형 구축과 모형 검·보정을 수행하였다. 유역 모형으로는 SWAT 모형을 이용하였고, 하천 흐름 모형으로는 EFDC 모형을 적용하였다. 먼저 SWAT 모형을 구축하고 2004년-2009년까지 낙동강 본류 및 지류 43개 지점에서의 모형 검·보정을 수행하였다.
각각의 소유역에 위치한 환경기초시설의 방류량을 산정하였고, 댐 방류량은 한국수자원공사의 각 댐별 제원과 일방류량 (1994∼2009)을 수집하여 관측 일유량법 (observed daily outflow method)으로 모델에 적용하였다. 취수량의 경우 생활용수 및 공업용수, 관개용수에 대한조절유량을 반영하기 위해 2007년도 상수도 통계연보(환경부, 2008b)를 이용하여 소유역별 생공용수량을 산정하였다. 관개용수 즉, 농업취수량의 경우 하천 유량과 강우량 및 수정 Penman식에 의해 산정된 증발산량 등을 이용하여 추산한 자료를 적용하였으며, 농업용수의 회귀율은 특징상 다른 용도의 용수에 비해 유효우량의 비중이 높다는 것과 회귀수량이 반복 이용됨으로써 반복 이용효율이 높으며 물의 소비기구가 복잡하여 기상, 토양, 작물의 종류, 생육조건, 수원공의 종류 등의 관련 인자에 크게 영향을 받기 때문에 조사자와 연구방법에 따라 최소 1.
는 실측 값의 평균이다. 한편, % difference와 NSE에 대해 Table 2와 같이 Donigian (2000)이 제시한 일반적인 모델효율의 범위와 신뢰 구간을 이용하였다. SWAT 모형에 의한 유량모의 결과는 본류와 지류구간에 대하여 Tables 3 and 4와 같이 나타내었다.
성능/효과
한편, 지류 구간 역시 일부 구간을 제외하고는 대부분 실측값의 경향을 적절히 예측하는 것을 확인하였으며, 지류보다는 본류 지점에서 실측값의 경향을 좀 더 잘 반영하는 것으로 나타났다.
8개 비교 지점에서의 결정계수는 0.72~0.94로서 실측과 모델링 예측 값이 비교적 잘 일치함을 확인할 수 있다.
94로 분석되어 두 자료간의 유의성이 있는 것으로 분석되었다. 검증기간 동안 % difference는 용전A (29.00)와 묘산A (109.69)를 제외한 모든 지점에서 25% 이하를 나타내었고, NSE의 경우 미천A (0.50), 위천A (0.55), 금호B (0.50), 묘산A (-1.16), 황강A (0.57), 남강C (0.55)를 제외 하고 0.62~0.97로 분석되었다. 결정계수는 묘산A를 제외 하고 0.
95의 범위로 분석되어 실측값과 모의값 간의 유의성이 높은 것으로 나타났다. 검증기간에 대한 % difference는 2.55~24.00, NSE는 0.60~0.95로 적절하게 모의되었음을 나타내었다. 또한 결정계수는 0.
관개용수 즉, 농업취수량의 경우 하천 유량과 강우량 및 수정 Penman식에 의해 산정된 증발산량 등을 이용하여 추산한 자료를 적용하였으며, 농업용수의 회귀율은 특징상 다른 용도의 용수에 비해 유효우량의 비중이 높다는 것과 회귀수량이 반복 이용됨으로써 반복 이용효율이 높으며 물의 소비기구가 복잡하여 기상, 토양, 작물의 종류, 생육조건, 수원공의 종류 등의 관련 인자에 크게 영향을 받기 때문에 조사자와 연구방법에 따라 최소 1.34%에서 최대 130.2%까지 다양하며 심한 편차를 보이고 있다(추태호, 2004).
먼저, SWAT모형을 이용하여 유출량을 모의한 결과, 본류와 지류 구간에 걸쳐 실측값과 비교적 잘 일치하는 것으로 나타났다. 본류 구간에서의 % difference는 낙본B와 M을 제외하고는 대부분 15% 이내인 것으로 나타났고, NSE과 R2 역시 대부분 적정 수준인 것으로 모의 되었다.
먼저, 전반적인 계산결과를 보면, 전 구간에서의 계산결과와 측정 데이터가 잘 일치함을 볼 수 있다.
먼저, SWAT모형을 이용하여 유출량을 모의한 결과, 본류와 지류 구간에 걸쳐 실측값과 비교적 잘 일치하는 것으로 나타났다. 본류 구간에서의 % difference는 낙본B와 M을 제외하고는 대부분 15% 이내인 것으로 나타났고, NSE과 R2 역시 대부분 적정 수준인 것으로 모의 되었다. 한편, 지류 구간 역시 일부 구간을 제외하고는 대부분 실측값의 경향을 적절히 예측하는 것을 확인하였으며, 지류보다는 본류 지점에서 실측값의 경향을 좀 더 잘 반영하는 것으로 나타났다.
4는 일간격 고유량에 대한 모델링과 실측값을 비교한 것이다. 비교결과 결정계수 R2의 값이 상류와 하류 유역 각각에 대해 0.66과 0.69인 것으로 나타났으며, 홍수기 고유량에 대한 실측과 예측 값이 비교적 잘 일치함을 알 수 있다. 결정계수의 정의에 대해서는 다음 부분에서 언급하였다.
한편, 지류 구간 역시 일부 구간을 제외하고는 대부분 실측값의 경향을 적절히 예측하는 것을 확인하였으며, 지류보다는 본류 지점에서 실측값의 경향을 좀 더 잘 반영하는 것으로 나타났다. 이상의 낙동강 유역과 같이 대규모이고, 각각의 중소규모 유역특성이 다양하고 광범위하다는 것을 고려해 볼 때 본 연구에서 구축된 모형이 허용 가능한 범위 내의 결과를 잘 예측하는 것으로 판단된다.
97을 나타내었다. 이상의 유량모의 결과를 정리하면, 일부 지류구간의 보정지점을 제외하면 본류 및 지류 전 지점에 걸쳐 구축된 SWAT 모형이 실측값의 경향을 적절히 잘 반영하는 것으로 나타났으며, 지류구간보다 본류구간에서 실측값의 일관성이 더 양호하고 실측자료의 경향과 더 유사한 결과를 나타냈다. 이는 저수기에도 본류하천의 유량이 풍부하며 특히 홍수기의 유출량 모의가 잘 반영된 결과로 판단된다.
특히, 하중도 좌안 측에서 통수단면적 감소에 의한 유속 증가 현상 및 그 하류단의 만곡 외측 부분에서의 유속 증가, 그리고 수면폭이 증가하는 구간에서의 유속 감소 현상 등 하천 흐름의 물리적 현상을 잘 예측하고 있음을 확인할 수 있다.
하천 흐름 모형은 2007년도 낙동강 본류 각 지점별 관측 수위를 비교한 결과 저수기와 홍수기를 비교적 잘 예측하는 것으로 나타났으며, 만곡부 및 하중도 부근에서의 유속 감속과 가속 현상 등 물리적 현상을 구축된 모형이 잘 재현하고 있음을 확인하였다.
또한 수산 수위표에서 계산된 저수기 수위가 국부적으로 큰 값을 보이는데, 이는 유입유량이 국부적으로 큰 값을 갖기 때문인 것으로 보인다. 한편, 각 관측소에서의 최대 홍수위를 비교해 보면, 지보, 구미, 화원 관측소에서 최대 홍수위를 다소 과소 산정하는 것을 볼 수 있으나 그 차이는 크지 않으며, 다른 관측소에서는 관측 데이터와 잘 일치하는 것을 확인할 수 있다. 한편, Table 5는 계산과 실측 수위에 대한 결정계수, % difference, NSE의 값을 보여준다.
후속연구
좀 더 명확한 모형 검보정을 위해서는 환경부의 8일 자료 뿐만 아니라 국토해양부 일일자료에 대해서도 검보정이 필요하다. 그러나 본 모형이 오염총량단위유역을 기반으로 구축된 것이고, 모형의 소유역 경계지점과 국토해양부의 일일 관측지점이 서로 일치하지 않는 지점이 많아 검보정 자료로 이용하기에는 한계가 있는 것으로 판단된다.
이상의 최근 연구동향을 살펴보면, 우선, 낙동강과 같은 대규모 유역에서의 유역 모형 및 하천 모형이 구축된 사례가 많지 않으며, 유역-하천모형의 연계를 통한 해석 역시 많지 않은 것을 볼 수 있다. 그러나 향후 대규모 유역에서의 수자원, 토지 이용, 환경 변화 등을 반영할 수 있고, 기후변화에 따른 가뭄과 홍수 등의 예측과 수질 문제 등 하천환경에 미치는 영향을 분석하기 위해서는 대상 유역에 대한 유역-하천의 통합모형 구축이 필요할 것으로 보인다.
한편, 1월 1일 주입된 염료의 이동을 계산한 결과, 상류경계에서 하구언까지의 염료 도달시간은 약 65일이 소요되는 것으로 나타났다. 본 연구는 유역-하천 연계모형을 이용한 기후변화에 따른 흐름과 수질 등 수환경 변화를 분석하기 위한 전단계의 연구로서, 낙동강 유역과 같은 대유역 평가를 위한 모형 구축에 그 의의가 있으며, 향후 본 모형이 다양한 연구에 적용될 수 있을 것으로 기대된다.
그러나 대부분이 치수목적으로 본류 중심의 수문관측이 이루어졌기 때문에 소유역 및 지류에서의 정확한 유출량을 산정하는 것은 어렵다. 이를 위해서는 경험식에 의한 방법과 비유량법, 유역내 유출모형을 이용하는 방법 등을 고려할 수 있을 것이다. 그러나 경험식에 의한 것과 비유량법에 의한 유출량 산정은 단순하고, 계산하기 쉽다는 장점이 있지만, 단순히 유역 크기에 의해 유출량을 산정하는 것으로서 강우의 분포 특성을 적절히 반영하지 못한다.
결정계수의 정의에 대해서는 다음 부분에서 언급하였다. 좀 더 명확한 모형 검보정을 위해서는 환경부의 8일 자료 뿐만 아니라 국토해양부 일일자료에 대해서도 검보정이 필요하다. 그러나 본 모형이 오염총량단위유역을 기반으로 구축된 것이고, 모형의 소유역 경계지점과 국토해양부의 일일 관측지점이 서로 일치하지 않는 지점이 많아 검보정 자료로 이용하기에는 한계가 있는 것으로 판단된다.
질의응답
핵심어
질문
논문에서 추출한 답변
국내에서 유출모형으로 광범위하게 사용되는 것은?
국내에서 유출모형으로 광범위하게 사용되는 것 중의 하나가 SWAT 모형 (Neitsch et al., 2004)이다.
경험식에 의한 방법과 비유량법, 유역내 유출모형을 이용하는 방법의 장점과 단점은?
이를 위해서는 경험식에 의한 방법과 비유량법, 유역내 유출모형을 이용하는 방법 등을 고려할 수 있을 것이다. 그러나 경험식에 의한 것과 비유량법에 의한 유출량 산정은 단순하고, 계산하기 쉽다는 장점이 있지만, 단순히 유역 크기에 의해 유출량을 산정하는 것으로서 강우의 분포 특성을 적절히 반영하지 못한다. 따라서 대규모 유역에서의 각 소유역별 유출량 해석을 위해서는 유역단위 유출모형에 의한 평가가 필요할 것이다.
SWAT 모형이란?
, 2004)이다. SWAT 모형은 미농림부에서 개발한 유역기반의 장기유출모형으로서 토양수분을 중심으로 다양한 토지이용에 대하여 적용되어 왔으며, 국내에서도 유출해석을 위해 광범위하게 적용되고 있다 (김남원 등, 2008; 김대근 등, 2007; 김정곤 등, 2006; 손태석 등, 2010; 정상만 등, 2008; 신현석 등, 2007). 먼저, 김남원 등 (2008)은 충주댐 유역을 대상으로 유출곡선지수를 산정하는 방안을 검토하였으며, 김대근등 (2007)은 SWAT의 적설 및 융설기작을 이용하여 강원도 쌍천 유역을 대상으로 적설 및 융설의 영향을 분석한 결과, 적설 및 융설의 영향을 고려할 경우 하천의 유황이 다소 증가하였고 실측유량과 유사한 경향을 나타내었다.
참고문헌 (23)
건설교통부(2000). 수자원장기종합계획, 건설교통부,한국수자원공사.
김남원, 이진원, 이정우, 이정은(2008). "SWAT을 이용한 충주댐 유역의 유출곡선지수 산정 방안." 한국수자 원학회논문집, 한국수자원학회, 제41권, 제12호, pp. 1231-1244.
정광욱(2007). Development of 3D integratted water quality simulation system by applying BASINS/ WinHSPF and EFDC in Han River basin and Paldang Reservoir. 박사학위논문, 건국대학교.
정상만, 서형덕, 김형수, 한규하(2008). "기후인자의 변화에 따른 대청댐유역의 유출민감도 모의평가." 한국수자원학회논문집, 한국수자원학회, 제41권, 제11호, pp. 1095-1106.
추태호(2004) 낙동강 유역의 월곡저수지 내 논 농업용수 회귀율 조사.분석. 대한토목학회논문집, 대한토목학회, 제24권, 제2b호, pp. 123-129.
환경부(2008a). 2007년 전국오염원 조사.
환경부(2008b). 2007년 상수도통계연보.
ASCE (1993). "ASCE task committee on definition of criteria for evaluation of watershed models, Criteria for Evaluation of Watershed Models." Irrigation Drainage Engineering, Vol. 119, No. 3, pp. 429-442.
Donigian, Jr. A.S. (2000). HSPF Training Workshop Handbook and CD. Lecture #19. Calibration and Verification Issues, Slide #L19-22, EPA Headquarters, Washington Information Center, 10-14 January, 2000, Presented and prepared for U.S. EPA, Office of Water, Office of Science and Technology, Washington, DC.
Mellor, G.L., and Yamada, T. (1982). "Development of a turbulence closure model for geophysical fluid problems." Rev Geophys. Space Phys, Vol. 20, pp. 851-875.
Neitsch, S.L., Arnold, J.G., Kiniry, J.R., Srinivasan, J.R., and Williams, J.R. (2004). Soil and Water Assessment Tool; The theoretical documentation (Version 2005), U.S. Agricultural Research Service.
Zhang, Y., and Jia, Y. (2009). 2D structured mesh generator user's manual, School of Engineering, The University of Mississippi.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.