본 연구에서는 Muskingum 하도추적모형을 수문학적으로 재해석하여 지체효과만을 고려하는 선형하천모형과 저류효과만을 고려하는 선형저수지모형의 선형결합으로 나타내었다. 유도된 모형은 일종의 순간단위도의 형태가 되며, 그 매개 변수는 Muskingum 모형의 매개변수와 동일하다. 즉, 추적시간간격 ${\Delta}t$ 또는지체시간 $T_c$ 후에 최초의 유출이 발생하게 되고, 총 유입량 중 x 만큼은 선형하천모형에 의해 저류효과 없이 빠져나가고 나머지(1-x) 만큼은 선형저수지모형에 의해 저류상수 $K_c$로 대변되는 저류효과를 나타내며 빠져나가는 형태이다. Muskingum 하도추적 모형과 그에 대응하는 순간단위도를 가상하도에 적용해 본 결과, 두 모형이 근본적으로 하도추적결과가 동일함을 확인하였다. 이러한 결과는 대청댐 방류량에 대한 금남 및 공주지점까지의 하도추적결과에서도 확인할 수 있었다.
홍수 예측을 위한 분포형 수문모형의 유출해석에서 하도추적은 수리학적 하도 추적과 수문학적 하도 추적 방법이 있다. 수리학적 하도 추적은 운동파 방정식, 확산파 방정식 등을 이용하여 수리현상을 시간과 공간으로 편미분하여 홍수량 예측을 한다. 수리적 하도 추적은 시간적, 공간적 안정조건(stability condition)을 만족해야된다. 면적이 큰 유역에서 적용할 때에는 계산에 소요되는 시간이 크다. 그러므로 국지호우로 인한 돌방홍수 예 경보를 위해서는 준실시간 또는 실시간 홍수 감시 및 예측이 필요하므로 계산에 소요되는 시간이 큰 수리학적 하도추적을 이용한 홍수 예측은 한계를 가진다. 본 연구에서는 유역면적이 큰 유역의 준실시간 홍수 감시 및 예측을 위하여 수문학적 하도추적 기법은 하천차수별 저류상수를 적용한 multi-Muskingum방법을 개발하여 모의하였다. multi-Muskingum 적용한 결과 모의시간이 상당히 단축되었으며 자료동화 기법을 통하여 모형의 정확도를 개선하였다.
수문학적 하도추적법의 하나인 Muskingum 모형은 미 육군공병단(U.S. Army Corps of Engineers)에 의해서 미국 Ohio 주의 Muskingum 유역에 홍수조절계획으로 처음 사용되었으며 모형의 구조 및 입력자료의 단순성에 비하여 비교적 우수한 결과를 모의할 수 있는 것으로 알려져 있다. 1938년 McCarthy에 의해서 개발되었고 구간내 총저류량은 prism 저류와 wedge 저류로 구분하여 prism 저류는 유출량에 wedge 저류는 유입량과 유출량의 차에 직접 비례한다는 가정하에 추적식을 개발하였다. 이후 지속적인 연구가 이뤄져 1985년 O'Donnel은 측방유입량(lateral inflow)을 상류단의 유입량에 비례하는 형태로 3-매개변수 muskingum 모형을 제안하여 추적계수의 결정을 선형대수(linear algebra)에서 동차(homogeneous)연립방정식 해를 구하는 Cramer 법칙인 matrix 기법을 적용하였다. 본 연구에서는 홍수사상으로부터 측방유입량이 고려되고 추적계수 결정에 있어서 직접 계산이 가능한 O'Donnel(1985)이 제안한 3-매개변수 muskingum 모형을 적용하였다. 추적계수들의 결정은 직접 matrix 기법을 적용하였고 적용대상은 낙동강 유역의 낙동 지점을 상류단으로 구미 지점을 하류단으로 선정하였다. 홍수사상은 낙동강 유량측정 조사사업 2005년${\sim}$2007년 보고서에 수록된 수문자료를 선정하여 관측치와 계산치를 비교하였고 홍수사상에 적용하여 수문곡선을 추정하였으며, 각각의 매개변수가 추적구간에 어떠한 영향을 미치는지 변수간의 관계를 분석하였다. 또한, 관측치와 계산치의 적합도 검증은 평균제곱근오차(root mean squar error; RMSE)와 모형 효율성 계수(model efficiency; ME)를 산정하여 분석하였으며, 하도 구간내 저류량은 대상구간에 대한 유입량과 유출량의 가중합에 비례한다는 선형모형을 적용하였다.
본 연구에서는 선형시스템 가정에 근거하여 하도구간에 대한 Muskingum 하도추적모형의 매개변수 결정방법을 제안하였다. 제안된 모형은 충주댐 유역에 적용되어 검토되었다. 추가적으로 영춘-충주댐 유역에 대해 총 7개의 호우사상을 대상으로 유출해석을 실시하고 그 결과를 검토하였다. 그 결과를 정리하면 다음과 같다. 먼저, 유역분할에 의해 발생하는 하도의 집중시간 및 저류상수는 상류 분할유역과 상류 분할유역을 포함한 하류 유역의 집중시간 및 저류상수의 차로써 표현가능하다. 이와 같은 방법으로 산정된 하도구간에서의 저류상수는 Muskingum 하도추적모형의 저류상수와 동일하며, 가중인자 역시 집중시간과 저류상수와의 비를 이용하여 간단히 산정할 수 있다. 둘째, Russel 계수와 Muskingum 모형의 가중인자는 서로 반비례 관계에 있으며 일반적으로 적용되고 있는 Russel 계수의 범위에 해당하는 가중인자의 범위는 0.4166-0.625이다. 마지막으로, 영춘-충주댐 구간을 대상으로 한 적용에서는 관측자료의 불확실성과 같은 한계에도 불구하고 제안된 방법의 유효성을 확인할 수 있었다.
하도 홍수추적과 관련하여 하천에서의 시 공간적 홍수파를 해석하는데 수리학적 방법과 수문학적 방법이 일반적으로 많이 이용되어 왔다. 수문학적 홍수추적 방법은 수리학적 방법에 비해 수행하기에는 비교적 간단하면서도 합리적인 정확성을 지닌다. 수문학적 홍수추적 방법 중 광범위하게 적용되어지고 있는 Muskingum 모형의 중요 변수인 저류상수는 유하시간과 매우 유사한 값을 가진다. 이러한 점에 착안하여 본 연구에서는 저류상수를 산정하기 위해 HEC-RAS를 이용한 유하시간을 산정하고, 하도거리, 하도경사, 유량 자료를 이용하여 유하시간에 대한 비선형 회귀곡선식을 개발하였다. 비선형 회귀곡선에 의해서 산정된 저류상수를 Muskingum 모형에 대입하여 구한 유출량은 HEC-RAS 1차원 부정류 모의를 적용하여 산정된 유출량과 비교하였다. 이와 함께 본 연구에서는 가중인자에 대한 영향 및 상하류 사이의 구간 분할에 대해서 검토하였는데, 그 결과 가중인자 값이 클수록 첨두홍수량이 올라가는 것으로 나타났으며, 구간 분할을 많이 할수록 RMSE가 감소하는 것으로 나타났다.
하도의 횡단 및 종단 지형자료와 조도계수를 이용하여 자연하천에 대한 Muskingum-Cunge 모형의 매개변수들을 추정하는 방법을 제안하였다. 우선 각 단면에서의 다양한 수위에 대하여 통수 단면 및 동수반경을 계산한 후, Manning 공식을 이용하여 유량을 산정한다. 이러한 과정은 하도에서의 모든 단면에 대하여 반복되며, 최종적으로 통수단면과 유량을 통한 회귀 분석에 의하여 매개변수들을 추정한다. 앞서 설명한 Muskingum-Cunge 모형의 매개변수 추정과정을 남한강 구간에 적용하였다. 추정된 매개변수들을 사용한 Muskingum-Cunge 모형의 계산결과가 HEC-1의 Muskingum-Cunge 모형에 비하여 동역학적 모형의 계산결과와 잘 일치하는 것으로 나타났다.
측방유입수가 고려되는 3변수 Muskingum하도추적모형을 낙동강 수계중 왜관에서 적포교구간의 12개 홍수사상에 대하여 적용하였고, 기존 방법인 2변수 Muskingum방법의 저류상수 K와 가중계수 x에 추가된 α는 측방유입수를 고려해주는 변수이다. 3변수모형의 추적방정식은 유한차분방정식으로 표현되며, 추적상수 결정은 Matrix Inversion에 의하여 직접 계산 가능하며, 이로부터 각 홍수사상의 K x 값을 구할 수 있다. 본 연구를 실유역에 적용하여 실측치와 비교하여본 결과 비교적 잘 맞음을 알 수 있어쓰ㅇ며 K와 x값은 하도특성인자로서 홍수규모와 관계되고 측방유입인자 α는 강우 특성에 의하여 지배되는 변수로 판단되었으며 α값은 측방유입량이 클수록 값이 커지는 성향으로 나타났다.
측방유입수가 고려되는 3변수 Muskingum하도추적모형을 낙동강수계중 왜관에서 적포교구간의 12개 홍수사상에 대하여 적용하였고, 기존방법인 2변수 Muskingum방법의 저류상수 K와 가중계수 x에 추가된 $\alpha$는 측방유입수를 고려해주는 변수이다. 3변수모형의 추적방정식은 유한차분 방정식으로 표현되며, 추적상수 결정은 Matrix Inversion에 의하여 직접 계산가능하며, 이로부터 각홍수사상의 K x $\alpha$값을 구할수 있다. 본 연구를 실유역에 적용하여 실측치와 비교하여본 결과 비교적 잘 맞음을 알 수 있었으며, K와 x값은 하도특성인자로서 홍수규모와 관계되고 측방유입인자 $\alpha$는 항특성에 의하여 지배되는 변수로 측방유입량이 클수록 값이 커지는 성향으로 나타났다.
실제 유역에서 지류유입량(tributary inflow)의 형태로든 사면류의 형태로든 측방유입은 반드시 존재한다. 측방유입이 하도의 지배적인 흐름이 되는 경우 이는 유출수문곡선의 종거값과 형태를 변화시키는데 중요한 역할을 하게 된다. 따라서 측방유입의 형태를 저류상수 및 집중시간과 같은 수문학적 특성으로 적절하게 표현할 수 있다면, 전체 유역 내에서 측방유입의 지체효과 및 저류효과를 파악하는데 크게 기여하게 된다. 측방유입과 관련된 선행연구들을 살펴보면, Saint-Venant 방정식을 근간으로 하는 연구가 주를 이루고, Muskingum 하도추적모형 또는 Muskingum-Cunge 하도추적 모형을 확장한 연구가 나머지 부분을 차지한다(Hayami, 1951; Dooge et al. 1982). 지금까지 수행된 대다수의 연구들은 수치해석적으로 측방유입의 유출량을 모의하거나 혹인 관측값이 존재하는 경우 역으로 측방유입의 특성을 유추한 것들로 다소 복잡하고, 관측값이 존재하지 않는 경우에는 적용이 어려운 문제점을 가지고 있다. 그러나 유역의 물리적인 특성과 주하도의 특성만을 이용하여 측방유입의 특성을 대략적으로 유추하는 것이 가능하다면, 전체유역과 각 소유역의 관계는 전체유역의 물리적인 특성과 주하도의 수문학적 특성만으로 충분히 파악할 수 있게 된다. 다시 말해 유역분할 시 각 소유역 사이의 관계를 고려하여 전체유역의 유출량을 파악하는 것이 가능해진다. 이에 본 연구에서는 Muskingum 하도추적모형을 재해석한 순간단위도를 이용하여 측방유입의 수문학적 해석을 시도하였다. 대상유역으로는 격자형태의 사각형과 삼각형 유역을 임의로 가정하였으며, 각각의 유역에서의 순간단위도를 선형하천모형과 선형저수지모형의 합으로 유도하였다. 이때 저류상수는 하도길이와 비례한다는 가정을 바탕으로 사각형과 삼각형 유역에서의 저류상수 및 집중시간을 유도하였다. 특히 유역 출구에서 최원점에 위치한 격자에서 유출이 발생시간을 집중시간으로 가정하였으며, 이 시점에서의 종거값과 기울기를 이용하여 저류상수를 유도하였다. 그 결과, 선형하천과 선형저수지모형 각각은 집중시간과 저류상수로 특징지어짐을 알 수 있었으며, 결정된 측방유입의 저류상수 및 집중시간이 적절한 것을 확인하였다. 이러한 결과는 향후 대유역에서 유역분할의 효과뿐만 아니라 홍수량 할당문제를 입증하는데 크게 기여할 것으로 기대된다.
1차원 Saint-Venant 방정식을 이용한 동력학적 부정류 계산모형은 하도의 지형과 흐름에 대한 물리적 특성을 잘 반영하고 있지만 하폭이 좁고 경사가 큰 하천에서는 해의 안정성 문제를 갖고 있기 때문에 홍수 예 경보를 목적으로 하는 하도추적 방법으로 항상 적합하다 할 수 없다. 본 연구에서는 이러한 동력학적 부정류 계산모형의 대안으로서 확산파 모형의 수치해법으로 볼 수 있는 Muskingmum-Cunge 모형을 남한강에 적용하여 분포형 매개변수를 추정하는 방법을 제시하였다. 먼저 남한강의 상류인 충주 조정지댐부터 여주 수위관측소까지 가용한 각각의 횡단면에 대해 최심하상고에서 제방고까지 다양한 수위를 가정하고 수면폭과, 통수단면적 그리고 동수반경을 계산한 후 Manning 식을 이용하여 유량을 산정한다. 이때 각 단면별로 유량에 대한 통수단면적과 수면폭에 대하여 회귀분석 방법을 이용해 단면별 매개변수를 추정한다. 추정된 매개변수를 Muskingum-Cunge 모형에 적용하여 상류로부터 하류까지 각 단면별로 하도추적을 수행한다. 계산된 결과는 HEC-HMS Musking-Cunge 모형과 첨두유량의 크기 및 무차원 RMS 등을 비교하였을 때 동력학적 모형의 계산결과에 잘 일치함을 알 수 있다.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.