بررسی حساسیت پارامترهای موثر بر روندیابی هیدروگراف سیل با روش موج پخشی دیفیوژن با مدل هیدرولوژیکی توزیعی WetSpaدر حوزه آبخیز زیارت گرگان

نوع مقاله : مقالات پژوهشی

نویسندگان

• نرگس جاویدان
• عبدالرضا بهره مند

دانشگاه علوم کشاورزی و منابع طبیعی ساری

https://doi.org/10.22067/jsw.v30i3.40498

چکیده

روندیابی سیل در رودخانه روشی است که به‌وسیله آن و با در دست داشتن هیدروگراف ورودی در نقطه‌ای از بالادست رودخانه، می‌توان هیدروگراف خروجی را در نقطه‌ای از پایین دست محاسبه کرد. در روندیابی رواناب در حوزه های آبخیز بارش مازاد با تکنیک های روندیابی به هیدروگراف جریان در پایین دست تبدیل می شود. در این تحقیق با استفاده از مدل هیدرولوژیکی – توزیعی WetSpa، روندیابی جریان سیل در حوزه آبخیز زیارت صورت گرفت. جهت اجرای مدل از داده های هیدرومتئورولوژی ساعتی از سال 1386 تا 1390 شامل داده‌های بارش، تبخیر و تعرق و دما به عنوان داده ورودی مدل استفاده شد. نتایج شبیه سازی انطباق خوبی بین هیدروگراف شبیه‌سازی شده و هیدروگراف مشاهداتی نشان می-دهد. در این مدل روندیابی جریان سطحی و جریان آبراهه با استفاده از روش معادله خطی موج پخشی صورت می گیرد. در مطالعه حاضر، بررسی حساسیت پارامترهایی از قبیل: شعاع هیدرولیکی، ضریب زبری کانال، آستانه شیب حداقل و آستانه سطح برای ترسیم شبکه آبراهه انجام شد. بررسی حساسیت پارامترها نشان داد که تاثیر تغییر فراوانی سیل و ضریب زبری نسبت به آستانه شیب حداقل و آستانه سطح بر روی هیدروگراف خروجی و هیدروگراف واحد حوزه بیشتر است.

کلیدواژه‌ها

• آستانه شیب حداقل
• داده‌های هیدرومتئورولوژی
• روندیابی سیل
• هیدروگراف شبیه‌سازی شده
• هیدروگراف واحد

20.1001.1.20084757.1395.30.3.11.5

عنوان مقاله [English]

Sensitivity Test of Parameters Influencing Flood Hydrograph Routing with a Diffusion-Wave Distributed using Distributed Hydrological Model, Wet Spa, in Ziarat Watershed

نویسندگان [English]

• narges javidan
• Abdolreza Bahremand

Sary University of Agricultural Sciences and Natural Resources

چکیده [English]

Introduction: Flood routing is a procedure to calculate flood stage and water depth along a river or to estimate flood hydrograph at river downstream or at reservoir outlets using the upstream hydrography . In river basins, excess rainfall is routed to the basin outlet using flow routing techniques to generate flow hydrograph.
A GIS-based distributed hydrological model, Wet Spa, has been under development suitable for flood prediction and watershed management on a catchment scale. The model predicts outflow hydrographs at the basin outlet or at any converging point in the watershed, and it does so in a user-specified time step. The model is physically based, spatially distributed and time-continuous, and simulates hydrological processes of precipitation, snowmelt, interception, depression, surface runoff, infiltration, evapotranspiration, percolation, interflow, groundwater flow, etc. continuously both in time and space, for which the water and energy balance are maintained on each raster cell. Surface runoff is produced using a modified coefficient method based on the cellular characteristics of slope, land use, and soil type, and allowed to vary with soil moisture, rainfall intensity and storm duration. Interflow is computed based on the Darcy’s law and the kinematic approximation as a function of the effective hydraulic conductivity and the hydraulic gradient, while groundwater flow is estimated with a linear reservoir method on a small subcatchment scale as a function of groundwater storage and a recession coefficient. Special emphasis is given to the overland flow and channel flow routing using the method of linear diffusive wave approximation, which is capable to predict flow discharge at any converging point downstream by a unit response function. The model accounts for spatially distributed hydrological and geophysical characteristics of the catchment.
Determination of the river flow hydrograph is a main target in hydrology. Hydrological modeling deals with calculation of watershed hydrograph using hydro-meteorological information and terrain data, and processes of transforming rainfall into a flood hydrograph and the translation of hydrographs throughout a watershed. Flow routing subjects hydrography transformation and translation throughout a river basin. The Wet Spa model used in this study is a simple grid-based distributed runoff and water balance simulation model that runs on an hourly time step. It predicts hourly overland flow occurring at any point in a watershed, hydrography at the outlet, and provides spatially distributed hydrologic characteristics in the basin, in which all hydrologic processes are simulated within a GIS framework (Bahremand, 2007). The Wet Spa model was originally developed by Wang et al. (1997) and adapted for flood prediction by De Smedt et al. (2000) and Liu et al. (2003).
Materials and Methods: The outlet is accomplished using the first passage time response function based on the mean and variance of the flow time distribution, which is derived from the advection–dispersion transport equation. The flow velocity is location dependent and calculated in each cell by the Manning equation based on the local slope, roughness coefficient and hydraulic radius. The hydraulic radius is determined according to the geophysical properties of the catchment and the flood frequency. The total direct runoff at the basin outlet is obtained by superimposing all contributions from every grid cell.
The routing of overland flow and channel flow is implemented by the method of the diffusive wave approximation. This method has been used in some recent GIS-based flood models (Fortin et al., 2001; Olivera and Maidment, 1999). Liu et al 2003 has presented the flow routing method of the WetSpa model in detail. A two-parameter response function, based on the average flow time and the standard deviation of the flow time, is proposed in this study. The flow time and its variance are determined by the local slope, surface roughness and the hydraulic radius for each grid cell. The flow path response function at the outlet of the catchment or any other downstream convergence point is calculated by convoluting the responses of all cells located within the drainage area in the form of the probability density function (PDF) of the first passage time distribution. This routing response serves as an instantaneous unit hydrograph and the total discharge is obtained by a convolution integral of the flow response from all generated spatially distributed runoff.
Starting from the continuity equation and the St.Venant momentum equation, assuming the one-dimensional unsteady flow, and neglecting the inertial terms and the lateral inflow to the flow element, the flow process can be modeled by the diffusive wave equation (Cunge et al., 1980)
In this study, flood routing is done as the main part of flow simulation of the distributed hydrological Wet Spa model in the Ziarat watershed. In order to execute the model, hourly hydrometeorological data for a period of four years (2007-2010) including rainfall, evapotranspiration, temperature, and discharge are used as inputs. Additionally, three main maps of the digital elevation model, soil map (texture), and land use are also applied and converted to digital formats. The result of the simulation shows a good agreement between the simulated hydrography and the observed one. The routing of overland flow and channel flow is implemented by the method of the diffusive wave approximation.
Results and Disscasion: The Wet Spa model has been applied in several studies, e.g. the Barebeek catchment in Belgium, the Alzette river basin in Luxembourg, the Hornad watershed in Slovakia, In this study, flood routing is done as the main part of flow simulation of the distributed hydrological WetSpa model in the Ziarat watershed. In order to execute the model, hourly hydrometeorological data for a period of four years (2007-2010) including rainfall, evapotranspiration, temperature, and discharge are used as inputs. Additionally, three main maps of the digital elevation model, soil map (texture), and landuse are also applied and converted to digital formats. The result of the simulation shows a good agreement between the simulated hydrography and the observed one. The routing of overland flow and channel flow is implemented by the method of the diffusive wave approximation. A sensitivity test shows that the parameter of flood frequency and the channel roughness coefficient have a large influence on the outflow hydrography and the calculated watershed unit hydrograph, while the threshold of minimum slope and the threshold of drainage area in delineating channel networks have a marginal effect.

کلیدواژه‌ها [English]

• Flood routing
• Hydrometeorological data
• Hydrogrph simulation
• Threshold of minimum slope
• Unit Hydrograph

مراجع

1-Azinmehr M., Bahremand A., and Kabir A. 2013. Analysis of flood flows using wave approximation (diffusion) in the spatial distribution of hydrologic model WetSpa. the International River Engineering Conference 9 shahid Chamran University,22-24 jan 2013, Ahwaz, 1-9.

2-Bahremand A., Corluy J., Liu Y.B., and De Smedt F. 2006. Stream flow simulation by WetSpa model in Hornad river basin, Slovakia. In: J. van Beek, M. Taal (eds.), floods, from Defence to Management, Toylor and Francis group, London, 415-422.

3-Bahremand A. 2007. Simulation the effects of reforestation on floods using spatially distributed hydrologic modelling and GIS, Phd Thesis, Brussel Belgium: Vrij Universiteit, 1-150.

4-Bahremand A., and De Smedt F. 2010. Predictive analysis and simulation uncertainty of a distributed hydrological model, Water Resouce Management, 24: 2869-2880.

5-Bahremand A., Yaghubi F., and Kabir A. 2010. Introduce WetSpa Model Water and Energy Transfer between Soil, Plant and Atmosphere, the International Watershed Management Science and Engineering Conference 6 Tarbiat Modarres University, 1-8.

6-Bayat B. 2009. Simulated river flow and analysis Use change effects on river flow using WetSpa models in GIS environment, A thesis submitted in partial fulfillment of the requirements for the degree of M.Sc Remote sensing and GIS, Shahid Beheshti university, School of Geosciences,1-142.

7-Chow V.T. (Ed). 1964. Handbook of Applied Hydrology, McGraw- Hill Book Company, New York, 7-25.

8-De Smedt F., Liu Y.B., and Gebremeskel S. 2000. Hydrological modeling on a catchment scale using GIS and remote sensed land use information, in: C.A. Brebbia (ed.), Risk Analysis II, WTI press, Boston, 295-304.

9-Eaguson P.S. 1970. Dynamic Hydrology, McGraw-Hill, New York, 1- 364.

10- Henderson F.M .1966. Open Channel Flow, McMillan, New York, 1- 522.

11-Kabir A., Mahdavi M., Bahremand A., and Noora N. 2010. Application of the distributed hydrological WetSpa model for runoff simulation in Gorganrood river basin, Iran, the International Watershed Management Science and Engineering Conference 6 Tarbiat Modarres University,1-8.

12-Liu Y.B., Gebremeskel S., De Smedt F., Hoffmann L. and Pfister L. 2003. A diffusive transport approach for flow routing in GIS-based flood modelling, Journal of Hydrology, 91-106.

13-Miller W.A., and Cunge J.A. 1975. Simplified equation of unsteady flow, eds, K. Mamood & V. Yevjevich, Unsteady flow in open channels, Water Resources Publications, Fort Collins, CO.

14-Molnar P., and Ramirez J.A. 1998. Energy dissipation theories and optimal channel characteristics of river networks, Water Resour, Research, 34:1809–1818.

15-Shafii M., and De Smedt F. 2009. Multi-objective calibration of a distributed hydrological model (WetSpa) using a Genetic Algoritm. Hydrol, Earth Syst, 13: 2137-2149.

16- Olivera F., and Maidment D. 1999. Geographical information system (GIS)-based spatially distributed model for runoff routing, Water Resources Research, 35: 1155–1164.

17-Wang Z., Batelaan O., and De Smedt F. 1997. A distributed model for water and energy transfer between soil, plants and atmosphere (WetSpa), Physics and Chemistry of the Earth, 21: 189–193.

18-Yaghubi F. 2010. Simulated river flow using WetSpa models ,Case study Chehe Chai watershed in Golestan Province, Journal of Soil and Water Conservation, Gorgan university of agricultural sciences and natural resources, 3:185-207.