• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.6
• /
• pp.119-130
• /
• 2010

The universal model for the parameter estimation of the Storage Function Model(SFM) was developed through the applications of the distributed model for various hypothetical watersheds and runoff conditions. The existing parameter estimation equations are based on observations and these equations which are derived from the restricted conditions are not sensitive to the variation of physical characteristics of a watershed. This study developed the universal model for the parameter estimation through the runoff simulations of 35,000 times. As the simulation results, we have known that the lag time is related to the longest stream channel characteristics and the storage coefficient is related to the watershed characteristics.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.143-143
• /
• 2015

기존의 빈도해석에서는 자료의 정상성을 가정하며, 이에 따라 적정모형 선정 시에 $x^2$ 검정이나 PPCC(Probability Plot Correlation Coefficient)검정과 같은 적합도 검정방법을 사용한다. 하지만 자료에서 경향성이 나타나거나 평균, 분산, 매개변수 등이 시간에 따라 변하는 등의 비정상성 현상들이 관측됨에 따라 비정상성 빈도해석에 관한 연구들이 활발히 진행되고 있다. 비정상성 빈도해석에서는 시간항과 같은 공변량이 포함된 매개변수를 가지는 비정상성 모형을 적용하게 되는데, 시간에 따라 매개변수가 계속 변하므로 매개변수에 따라 검정통계량이 고정되어 있는 기존의 적합도 검정방법의 적용이 어렵다. 따라서 비정상성 빈도해석의 적정 모형 선정에 적용할 수 있는 방법으로 최우도 함수에 기반한 모형 평가 방법인 AIC와 BIC가 추천되고 있으며 자료길이가 충분하지 않은 경우에는 AIC 대신하여 AICc의 사용이 추천되고 있다. 본 연구에서는 극치사상을 나타내는데 적합한 분포형인 GEV분포형의 위치, 규모 매개변수를 시간항으로 나타낸 다양한 비정상성 GEV모형에 대하여 Monte-Carlo 모의실험을 통해 AICc와 BIC의 적용성을 검토하였으며, 비정상성이 관측되는 실측 자료에 적용해보았다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.288-288
• /
• 2016

유역단위 오염부하량 산정에는 SWAT, HSPF 등의 물리적 매개변수 기반 분포형 모형이 주로 사용되고 있으나, 공간분포형 입력자료로 인한 많은 매개변수는 모의 과정을 복잡하게 하며, 보정 과정에 있어 많은 시간과 노력을 요구하는 단점이 있다. 이로 인해 실무에서는 원단위법이나 유량-부하량 관계식과 같은 통계적 분석에 의한 회귀식이 주로 사용되고 있다. 그 중 LOADEST는 회귀식 기반 프로그램으로, 다양한 연구자들에 의해 연구되고 있으나, 수질 모형과의 모의능력을 비교하는 연구는 부족하다. 본 연구에서는 청미천 상류유역을 대상으로 유역특성에 따른 LOADEST 기반 회귀식의 매개변수를 추정하여 오염부하량을 모의하고, SWAT 모형에 의한 오염부하량 모의결과와 비교 평가하고자 한다. 모형의 구동 및 회귀식 매개변수 추정에 필요한 입력 자료는 용인시 백암면 일대에서 2013년부터 2015년까지 모니터링한 수질, 유량 및 기상자료와 지형자료 (토지이용도, 토양도, 수치표고자료)를 이용하여 구축하였다. LOADEST 기반 회귀식의 매개 변수 추정은 김계웅 (2015)이 개발한 방법을 사용하였으며, 유역면적, 토지이용비율 등은 지형자료를 이용하여 산정하였다. SWAT 모형의 보정은 2013년부터 2014년까지의 자료를 이용하였으며, 2015년 자료를 이용하여 검정하였다. 본 연구의 결과는 비점오염원 모델에 대한 이해를 넓히고, 오염부하량 모의를 위한 모형 선정에 있어 도움이 될 수 있을 것으로 기대한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.249-253
• /
• 2008

최근 국내에서 사용되고 있는 격자분포형 강우-유출 모형 $Vflo^{TM}$을 실무에 효율적으로 적용하기 위하여 매개변수값의 결정과 그 민감도를 분석하고, 중랑천 유역의 토지이용 및 지형특성을 고려한 적정 크기의 격자에 대한 연구를 수행하였다. 연구대상유역으로는 토지이용현황이 상이한 중랑천 상류(자연유역), 우이천(복합유역) 및 청계천(도시유역) 유역을 선택하였다. $2003{\sim}2006$년까지 5개 호우사상의 실측자료로부터 매개변수값의 최적화를 실시하였으며 수치모의된 결과로 부터 조도계수, 투수계수, 초기함수비 등의 매개변수에 의한 계산치의 민감도를 분석하였다. 검토결과 지표면 조도계수가 가장 예민했으며, Green-Ampt식에 적용되는 투수계수의 영향이 큰 것으로 분석되었다. 또한, 수문학적 지형특성 반영을 위한 적정 격자크기는 자연유역(중랑천 상류)과 복합유역(우이천)에서는 100과 200 m에서 관측 유출수문곡선과 대체로 일치하나, 격자크기가 300 m에서는 지형특성을 반영하지 못해 유출수문곡선이 왜곡되는 경향을 나타내어 불투수율 80% 이하인 자연 및 복합유역에서는 DEM 격자크기를 $200{\sim}300\;m$이하로 적용하는 것이 적정하다고 판단된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.904-904
• /
• 2012

강우-유출모형중 단기 홍수 사상 모형은 주로 집중형으로 하천유량 모의에 중점을 두어 왔고, 장기유출모형은 분포형으로 기저유출 모의에 중점을 두어 왔다. 단기 홍수사상 모형은 유역 전반에 걸친 수자원의 분포특성을 파악하기 힘들다는 단점이 있고, 장기유출모형은 단기 홍수사상을 분석하기 힘들다는 단점이 있다. 이와 같은 단점은 두 모형을 함께 이용하여 보완할 수 있다. 따라서 본 연구에는 단기 홍수사상모형인 HEC-HMS모형과 장기유출모형인 SWAT모형을 연동하기 위하여 두 모형에 공통으로 적용되는 매개변수 CN값을 적절히 구해 이용하여야 한다. 두 모형에서 CN값을 얻는 방법은 첫 번째, HEC-HMS모형의 홍수기 호우사상별 최적평균 CN값을 구하는 것이고, 두 번째, SWAT모형 내에서 적용되는 일별 CN값을 구하는 것이다. 그리고 두 모형의 매개변수를 이용하여 보정 3년(2002년~2004년), 검증 2년(2005년~2006년), 적용 3년(2007년~2009년)을 모델링하여 SWAT-HMS모형의 유출분석 신뢰성을 판단하였다. 보정 3년 결과로 선행 5일강우량 기준으로 150mm 이하이면 모의가능하고, 150mm 이상이면 과소 모의되는 SWAT모형의 단점을 알 수 있었다. 검증을 실시한 2년은 SWAT-HMS모형의 중요 매개변수인 CN값을 HEC-HMS모형의 CN값과 SWAT모형의 CN값을 적용한 결과를 비교하여 상관계수 0.90과 0.95를 얻었으며, SWAT-HMS모형에 적용되는 매개변수 CN값은 SWAT모형의 CN값을 사용할 경우 더욱 신뢰할 수 있음을 알았다. 적용 3년(2007년~2009년)은 SWAT모형의 CN값을 SWAT-HMS모형에 적용하여 유출분석을 하였으며, 관측치와 모의치의 상관계수 0.95라는 신뢰할 수 있는 결과를 얻었다. 그리고 관측치와 모의치의 첨두유량은 평균적으로 5% 이내의 차이로 정확한 모의결과를 얻을 수 있었다. 따라서 본 연구는 격자분포형 유출모형인 SWAT모형과 집중형 유출모형인 HEC-HMS모형을 두 모형에 공통으로 적용되는 매개변수 CN값을 이용하여 두 모형이 가지는 단점을 해소하여 보다 정확한 유출분석 결과를 얻을 수 있는 방법을 제시 하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.366-366
• /
• 2017

하천의 유량관측 자료는 지표수자원의 확보와 수공구조물의 설계를 위해 가장 기초적인 수문자료로써 정밀하고 지속적인 관측을 요구한다. 최근 유량 관측법은 접촉식 유속측정 방법의 단점을 보완한 전자파 표면유속계나 영상분석기법을 적용한 표면영상유속계(SIV)가 활용되고 있다. 이들 관측장비는 표면유속 관측법에 의해 유량을 측정하므로 보다 정밀한 유량자료를 확보하기 위해서는 측정 영역의 표면유속과 단면의 평균유속에 대한 해석이 필요하다. 이 연구에서는 제주도 남부지역에 위치한 강정천을 대상으로 2011년 7월부터 2015년 6월 까지 월 1~2회 현장관측한 ADCP 자료를 활용하여 Chiu(1987)가 제안한 2차원 유속분포식의 매개변수를 추정하여 정밀한 유량을 산정하였다. 또한 표면영상유속계(SIV)로 산정된 표면유속을 Chiu-2 차원 유속분포식에서 평균유속으로 환산하여 기존의 표면유속을 일률적으로 적용한 수심평균유속 환산계수인 0.85의 적용 값과 비교 분석하였다. 대상하천의 현장에서 72회 관측된 ADCP 자료를 활용하여 각각의 최대유속과 평균유속을 분석하고 엔트로피 계수(M)를 산정한 결과와 유속의 공간적 분포를 모델링하기 위해 제시되는 $_{urf}$를 산정하였으며, 산정된 계수 값을 이용하여 표면유속을 계산한 결과와 ADCP의 관측된 표면유속의 $^2$는 0.874로 나타났다. 이러한 결과는 Chiu-2차원 유속분포식을 연구대상하천에 적용하는 과정에서 추정되는 매개변수의 평균값 사용에 대한 타당성을 보여준다. 이 후 추정된 하천매개변수를 하천현장에 적용성 확인을 위해 강정천의 동일 관측지점에서 표면영상유속계(SIV)를 사용한 표면유속과 유량을 산정함과 동시에 ADCP에 의한 유속 및 유량과 비교 분석하였다. 표면영상유속계(SIV)로 분석된 유속 벡터를 Chiu-2차원 유속분포식에 적용하여 산정된 유량과 기존의 수심평균유속환산계수 0.85를 적용한 유량은 각각 $0.7171m^3/s$과 0$0.5758m^3/s$였다. ADCP 평균유량 $0.6664m^3/s$과의 오차율은 각각 7.63%, 13.64%로 나타나 Chiu-2차원 유속분포식을 적용한 유량이 수심평균유속환산계수 0.85를 적용한 유량에 비해 작은 오차율을 보였다.

• Journal of Korea Water Resources Association
• /
• v.43 no.12
• /
• pp.1011-1027
• /
• 2010

In hydrologic modeling, prediction uncertainty generally stems from various uncertainty sources associated with model structure, data, and parameters, etc. This study aims to assess the parameter uncertainty effect on hydrologic prediction results. For this objective, a distributed rainfall-sediment yield-runoff model, which consists of rainfall-runoff module for simulation of surface and subsurface flows and sediment yield module based on unit stream power theory, was applied to the mesoscale mountainous area (Cheoncheon catchment; 289.9 $km^2$). For parameter uncertainty evaluation, the model was calibrated by a multi-objective optimization algorithm (MOSCEM) with two different objective functions (RMSE and HMLE) and Pareto optimal solutions of each case were then estimated. In Case I, the rainfall-runoff module was calibrated to investigate the effect of parameter uncertainty on hydrograph reproduction whereas in Case II, sediment yield module was calibrated to show the propagation of parameter uncertainty into sedigraph estimation. Additionally, in Case III, all parameters of both modules were simultaneously calibrated in order to take account of prediction uncertainty in rainfall-sediment yield-runoff modeling. The results showed that hydrograph prediction uncertainty of Case I was observed over the low-flow periods while the sedigraph of high-flow periods was sensitive to uncertainty of the sediment yield module parameters in Case II. In Case III, prediction uncertainty ranges of both hydrograph and sedigraph were larger than the other cases. Furthermore, prediction uncertainty in terms of spatial distribution of erosion and deposition drastically varied with the applied model parameters for all cases.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.33 no.3
• /
• pp.122-130
• /
• 2021

A wind-speed estimation at the arbitrary elevations is key component for the design of the offshore wind energy structures and the computation of the wind-wave generation. However, the wind-speed estimation of the target elevation has been carried out by using the typical functions and their typical parameters, e.g., power and logarithmic functions because the available wind speed data is limited to the specific elevation, such as 2~3m, 10 m, and so on. In this study, the parameters of the vertical profile functions are estimated with optimal and analyzed the parameter ranges using the HeMOSU-1 platform wind data monitored at the eight different locations. The results show that the mean value of the exponent of the power function is 0.1, which is significantly lower than the typically recommended value, 0.14. The values of the exponent, the friction velocity, and the roughness parameters are in the ranges 0.0~0.3, 0~10 (m/s), and 0.0~1.0 (m), respectively. The parameter ranges differ from the typical ranges because the atmospheric stability condition is assumed as the neutral condition. To improve the estimation accuracy, the atmospheric condition should be considered, and a more general (non-linear) vertical profile functions should be introduced to fit the diverse profile patterns and parameters.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.180-180
• /
• 2018

가능최대강수량(PMP)은 대규모 수공구조물의 설계 시 기준이 되는 강수량으로, 최근 대규모 거대재난에 대비한 대피계획수립에 PMP를 활용하려는 움직임이 있으며 PMP에 대한 국내 연구가 활발히 수행되고 있다. PMP를 추정하기 위해 Hershfield의 통계적 방법에 대한 간단한 대안이 제안되었다. PMP는 물리적인 강우량 상한계로, 확률론적 개념과는 모순적이다. 또한, Hershfield의 PMP는 연 최대 시계열 평균의 선형함수로 주어지는 모양 매개변수를 가지는 GEV 분포의 약 60,000년 빈도임이 밝혀졌다. 따라서 본 연구에서는 Hershfield의 방법을 확률론적으로 해석하는 것이 바람직할 것으로 판단하였고, 기상청 ASOS 및 AWS 자료를 이용하여 우리나라 각 지점자료 중 10년 이상의 자료를 사용하여 Hershfield 방법을 적용하여 PMP를 산정하였다. 각 지점의 빈도계수를 구하여 우리나라 자료에 적합한 확률분포의 형태를 적용하였고, 분포형의 매개변수 값을 추정하였다. 또한, Hershfield의 빈도계수와, 우리나라 자료에 해당하는 빈도계수가 몇 년 빈도로 계산되는지 각각 확인해 보았다. ASOS 및 AWS 자료를 이용하여 연 최대 강수량 시계열 평균과 모양 매개변수의 관계 공식 또한 구성하였다. 본 연구의 방법을 검증하기 위하여 우리나라에서 제일 오래된 자료(57년)인 서울지점 자료를 이용하여 경험적인 분포함수와 본 연구에서 제안하고 있는 방법을 비롯한 다양한 방법을 통하여 구한 분포함수를 비교하여 도시하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.1137-1141
• /
• 2004

본 연구에서는 Probability Weighted Moments의 새로운 선형조합기법인 LQ-moments를 이용하여 GEV 분포의 매개변수를 추정하고 L, LH, LQ-moments를 사용하여 뉴욕주의 Donnattsburg에 위치한 Independence River의 홍수량을 빈도 해석하였다. LH, LQ-moments가 제시된 근본적인 이유는 L-moments가 극치값에 내해 지나치게 민감한 단점을 보완하기 위해서인데, 이번 연구의 결과에 의하면 오히려 LH, LQ-moments가 극치값에 대해 민감하게 반응하여 부정확한 결과가 도출되었다. 그러므로 항상 LH, LQ-moments가 L-moments의 대안이 될 수 있는 것은 아님을 알게 되었다. 그리고 수학적 유도에서 L, LH, LQ-moments는 좀더 쉽고 간편한 메개변수 추정을 위해 Probability Weighted Moments의 선형조합을 통해 고안되었다는 공통점을 가지고 있지만, 이 점을 제외한 나머지 부분의 수식 유도에서는 서로 많은 차이가 있어서 지역적인 특성과 확률분포형의 특성을 고려하여 L, LH, LQ-moments 중에서 선별 사용해야 할 것이다.