با همکاری انجمن آبخیزداری ایران

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

نویسندگان

1 گروه مهندسی عمران، واحد خرم آباد، دانشگاه آزاد اسلامی، خرم آباد، ایران

2 گروه مهندسی عمران، واحد استهبان، دانشگاه آزاد اسلامی، استهبان، ایران

چکیده

مقدمه
تخمین سیلاب حوزه‌های آبخیز بدون آمار برای طراحی سازه‌های هیدرولیکی بسیار حائز اهمیت است. جریان کل سیلاب در حوزه‌های آبخیز از دو بخش جریان سطحی و زیرسطحی تشکیل شده است. در حوضه‌های با نفوذپذیری بالا جریان زیرسطحی سهم مهمی در رواناب کل حوضه دارد و در این موضوع تحقیقات کمی انجام شده است. مکانیسم رواناب در حوضه‌ها به دو شکل مکانیسم هورتونی که خاک از بالا اشباع می‌شود و مکانیسم دانی بلاک که خاک از زیر اشباع می‌شود و در این تحقیق از مکانیسم دانی بلاک استفاده شده است.
 
مواد و روش‌ها
یکی از روش‌های پیش‌بینی هیدروگراف جریان سطحی و زیرسطحی حوضه‌ها، استفاده از روش هیدروگراف واحد لحظه‌ای ژئومورفولوژیک بر اساس اطلاعات ژئومورفولوژیک حوضه‌های فاقد آمار است. اطلاعات ژئومورفولوژیکی حوضه‌ها با استفاده از نرم‌افزار ArcGIS و الحاقیه‌های هیدرولوژیکی محاسبه می‌شوند. روش GIUH، قابلیت جداسازی جریان سطحی و زیرسطحی را از کل جریان دارد. مدل GIUH در تحقیقات گذشته بیشتر برای تخمین جریان سطحی استفاده ‌شده است، لذا معادلات این مدل گسترش یافته تا بتوان هیدروگراف جریان زیر سطحی کل حوضه را تخمین زد. در این پژوهش از مدل GIUH برای تخمین جریان سطحی و زیرسطحی دو حوضه کسیلیان در ایران و حوضه گاجاس در کشور هند استفاده شد.
 
نتایج و بحث
مدل هیدروگراف واحد لحظه‌ای ژئومورفولوژیکی (GIUH) ارائه‌ شده در این پژوهش و برای تحلیل جریان سطحی و زیرسطحی دو حوضه کسیلیان (ایران) و گاجاس (هند) به‌کار گرفته شد. نتایج حاصل از شبیه‌سازی برای چهار رویداد بارندگی-رواناب در هر دو حوضه نشان داد که مدل توانایی مناسبی در تخمین هیدروگراف کل و مؤلفه‌های آن دارد. مقایسه دبی‌های پیک برآوردی و مشاهده‌ای در حوضه کسیلیان نشان داد که مقدار دبی پیک شبیه‌سازی‌شده در تاریخ ۱۰ می ۱۹۹۲ برابر با 10.1 مترمکعب بر ثانیه بوده، درحالی‌که مقدار مشاهده‌ای آن 11.8 مترمکعب بر ثانیه است. این اختلاف در دیگر رویدادها نیز در محدوده سه تا ۱۶ درصد قرار داشت که نشان‌دهنده دقت قابل‌قبول مدل در تخمین رواناب است. در حوضه گاجاس، مقدار حداکثر دبی کل جریان بین ۴۴ تا ۱۱۰ مترمکعب بر ثانیه متغیر بود و جریان زیرسطحی حدود پنج تا شش درصد از کل جریان را تشکیل می‌داد. در حوضه کسیلیان، مقدار دبی کل بین 1.6 تا ۱۲ مترمکعب بر ثانیه و مقدار دبی پیک جریان زیرسطحی بین ۳۵ تا ۶۰ لیتر بر ثانیه برآورد شد. بررسی ارتباط بین میزان بارندگی و دبی پیک جریان زیرسطحی نشان داد که کاهش بارندگی منجر به کاهش دبی پیک جریان زیرسطحی شده است. تحلیل حساسیت مدل نشان داد که ضریب هدایت هیدرولیکی خاک یکی از مؤثرترین پارامترها بر روی شبیه‌سازی جریان زیرسطحی است. در خاک‌های با ضریب هدایت هیدرولیکی بالا، جریان زیرسطحی سهم بیشتری از کل جریان را تشکیل داده و زمان تأخیر هیدروگراف افزایش می‌یابد. برای نمونه، در صورتی که ضریب هدایت هیدرولیکی خاک حوضه کسیلیان از 0.0025 به 0.0009 متر بر ثانیه کاهش یابد، مقدار دبی پیک جریان زیرسطحی از 0.35 به 1.3 مترمکعب بر ثانیه افزایش خواهد یافت. همچنین، بررسی تأثیر ضریب مانینگ صفحات نشان داد که کاهش مقدار این ضریب از 0.2 به 0.4، موجب کاهش ۳۱ درصدی مقدار دبی پیک سیلاب می‌شود. این یافته‌ها بر اهمیت ویژگی‌های هیدرولوژیکی و ژئومورفولوژیکی حوضه در برآورد دقیق رواناب و طراحی سازه‌های کنترل سیلاب تأکید دارد.
 
نتیجه‌گیری
در این پژوهش، مدل هیدروگراف واحد لحظه‌ای ژئومورفولوژیکی (GIUH) برای تخمین رواناب سطحی و زیرسطحی در دو حوضه کسیلیان (ایران) و گاجاس (هند) ارزیابی شد. نتایج نشان داد که مدل دقت مناسبی در برآورد دبی پیک دارد، به‌طوری که خطای تخمین در کسیلیان بین سه تا ۱۶ درصد و در گاجاس بین 1.6 تا ۱۲ درصد بود. همچنین، جریان زیرسطحی در حوضه‌های با نفوذپذیری بالا سهم قابل‌توجهی دارد و در حوضه گاجاس به‌طور میانگین پنج تا شش درصد از کل جریان را تشکیل داد. تحلیل حساسیت نشان داد که افزایش ضریب هدایت هیدرولیکی، پیک جریان زیرسطحی را افزایش می‌دهد، درحالی‌که کاهش ضریب مانینگ صفحات موجب افزایش پیک سیلاب شد. این یافته‌ها اهمیت ویژگی‌های ژئومورفولوژیکی و هیدرولوژیکی را در مدل‌سازی رواناب تأیید می‌کند. درنهایت، مدل GIUH می‌تواند ابزار مناسبی برای مدیریت سیلاب و ارزیابی پاسخ هیدرولوژیکی حوزه‌های آبخیز باشد.

کلیدواژه‌ها

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

Estimating the surface and subsurface flow of catchments using geomorphological instantaneous unit hydrograph method

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

  • Hamidreza Babaali 1
  • Mostafa Akbari Kheirabadi 2
  • touraj sabzevari 2

1 Department of civil Engineering, Khorramabad Branch, Islamic Azad University, Khoramabad, Iran

2 Department of civil Engineering, Estahban Branch, Islamic Azad University, Estahban, Iran

چکیده [English]

Introduction
Flood estimation in ungauged watersheds is crucial for designing hydraulic structures. The total flood flow in watersheds consists of surface and subsurface flows. In highly permeable watersheds, subsurface flow significantly contributes to total runoff, yet limited studies have addressed this aspect. Runoff mechanisms in watersheds generally follow two models: the Hortonian mechanism, where the soil saturates from the top, and the Dunne mechanism, where saturation occurs from below. This study adopts the Dunne mechanism for runoff generation.
 
Materials and methods
One method for predicting surface and subsurface hydrographs in ungauged watersheds is the Geomorphologic Instantaneous Unit Hydrograph (GIUH), which utilizes geomorphological data. These data were derived using ArcGIS and hydrological extensions. The GIUH model can separate surface and subsurface flow components. While previous studies have primarily used GIUH for surface flow estimation, the model equations were expanded in this study to estimate the total watershed subsurface hydrograph. The GIUH model was applied to estimate surface and subsurface runoff in two watersheds: Kasillian in Iran and Gagas in India.
 
Results and discussion
The proposed GIUH model was used to analyze surface and subsurface flows in the Kasillian (Iran) and Gagas (India) watersheds. Simulation results for four rainfall-runoff events in each watershed demonstrated that the model effectively estimated total hydrographs and their components. A comparison of estimated and observed peak discharges in the Kasillian watershed showed that the simulated peak discharge on May 10, 1992, was 10.1 m³/s, whereas the observed value was 11.8 m³/s. The error margin across different events ranged from 3% to 16%, indicating an acceptable model accuracy in runoff estimation. In the Gagas watershed, total peak discharge varied between 44 and 110 m³/s, with subsurface flow contributing approximately 5%–6% of the total flow. In the Kasillian watershed, total discharge ranged from 1.6 to 12 m³/s, while peak subsurface discharge was estimated between 35 and 60 L/s. The relationship between rainfall and subsurface peak discharge revealed that lower rainfall led to reduced subsurface peak discharge. Sensitivity analysis showed that hydraulic conductivity was one of the most influential parameters in subsurface flow simulation. In soils with high hydraulic conductivity, subsurface flow accounted for a larger portion of total flow, and the hydrograph lag time increased. For example, reducing the hydraulic conductivity in the Kasillian watershed from 0.0025 to 0.0009 m/s increased subsurface peak discharge from 0.35 to 1.3 m³/s. Additionally, reducing the Manning’s roughness coefficient from 0.2 to 0.4 resulted in a 31% decrease in flood peak discharge. These findings highlight the importance of hydrological and geomorphological characteristics in accurate runoff estimation and flood control structure design.
 
Conclusion
This study evaluated the GIUH model for estimating surface and subsurface runoff in the Kasillian (Iran) and Gagas (India) watersheds. Results showed that the model provided accurate peak discharge estimates, with estimation errors ranging from 3% to 16% in Kasillian and 1.6% to 12% in Gagas. The subsurface flow played a significant role in highly permeable watersheds, contributing 5%–6% of total runoff in the Gagas watershed. Sensitivity analysis revealed that increasing hydraulic conductivity led to higher subsurface peak discharge, whereas reducing the Manning’s coefficient increased flood peaks. These findings confirm the importance of geomorphological and hydrological characteristics in runoff modeling. Ultimately, the GIUH model can serve as a useful tool for flood management and watershed hydrological response assessment.
 

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

  • Gagas
  • Kasillian watershed
  • Peak discharge
  • Subsurface runoff
  • Surface runoff
Amini, A., Hosseini, M., Mohammad khan, S., Moghimi, E., Yamani, A.,  2021. Evaluation of methods estimating overland/streamflow travel-times for predicting flood hydrograph based on PGIUH model (case study: Amameh Watershed, Tehran Province). Quantit. Geomorphol. Res (in Persian).
Aryal, S.K., O’Loughlin, E.M., Mein, R.G., 2005. A similarity approach to determine response times to steady-state saturation in landscapes. Adv. Water Resour. 28, 99-115.
Chang, C.H., Lee, K.T., 2008. Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic-index threshold and geomorphology-based runoff model. Hydrol. Process. 22, 802-812.
Derakhshan. 2010. Study of flood potential of Kasilian watershed using geographic information system. J. Applied Res. Geograph. Sci. (in Persian), 16(13), 51-63.‎
Gupta, V.K., Waymire, E., Wang, C.T., 1980. A representation of an instantaneous unit hydrograph from geomorphology. Water Resour. Res. 16)5(, 863-870.
Henderson, F.M., Wooding, R.A., 1964. Overland flow and ground water flow from a steady rainfall of finite duration, J. Geophys. Res., 69(7), 1531-1540.
Jabbari, I., Asadi, M., Hesadi, H., 2023. Evaluation and comparison of GIUH model hydrograph simulation results in different rainfall seasons. Geograph. Environ. Hazards 12(3), 241-254 (in Persian).
Keshtkaran, P., Sabzevari, T., Karami Moghadam, M., 2018. Estimation of runoff in ungauged catchments using the Nash non-dimensional unit hydrograph (case study: Ajay and Kasilian Catchments). Water Resource. Eng. 11(36), 1-10 (in Persian).
Kumar, A., Kumar, D., 2008. Predicting direct runoff from hilly watershed using geomorphology and stream-order law ratios: case study. Hydrol. Eng. 13(7), 570-576.
Lee, K.T., Chang, C.H., 2005. Incorporating subsurface-flow mechanism into geomorphology-based IUH modeling. J. Hydrol. 311)1-4(, 91-105.
Lee, K.T., Yen, B.C., 1997. Geomorphology and kinematic-wave based hydrograph derivation. J. Hydrol. Eng. 123(1),73-80.
Mahmudi, F., Yamani, M., Bahrami, S., 2007. Evaluation of Geomorphological Instantaneous Unit Hydrograph Model (GIUH) in Kheez Kangir Watershed (Ivan Gharb). Geograph. Res. 60(39), 1-14 (in Persian).
Najafi, M., Behbahani, S.M.R., Abdollahi, J., Hosseini, S.M., 2009. Comparative study of geomorphologic artificial intelligent model and GIUH For direct runoff computations. Iran-Water Resour. Res. 5(2), 10-9 (in Persian).
Rodríguez‐Iturbe, I., González‐Sanabria, M., Bras, R.L., 1982. A geomorphoclimatic theory of the instantaneous unit hydrograph. Water Resour. Res. 18(4), 877-886.
Rodriguez-Iturbe, I., Valdes, J.B., 1979. The geomorphologic structure of hydrologic response. Water Resour. Res. 15(6), 1409-1420.
Sabzevari, T., Fattahi, M.H., Mohammadpour, R., Noroozpour, Sh., 2013. Prediction of surface and subsurface flow in catchments using the GIUH, under publication. J. Flood Risk Manage. 6(2), 135-145.
Sabzevari, T., Rezaeian, M., 2011. Flood estimation of  ungauged catchments by GIUH model.  Research Project of Islamic Azad University, Estehban Branch and Fars Regional Water Company (in Persian).
Sabzevari, T., Talebi, A., Ardakanian, R., Shamsai A., 2010. A steady-state saturation model to determine the subsurface travel time (STT) in complex hillslopes, Hydrol. Earth Syst. Sci. 14, 891-900.
Sallary jazi, M., Adib, A., Mamudian Shushtari, M., Akhunali, M.A., 2009. Evaluation of GIUH-NASH and GCIUH-CLARK models in Kesilian Basin. 8th International Congress of Civil Engineering (in Persian).
Shokoohi, A., Azizian, A., Jemaat, R., Singh, V., 2017. Sensitivity analysis of KW-GIUH rainfall runoff model with respect to infiltration methods and roughness coefficients. Manage. Engineer. Watershed J. 9(3), 162-175 (in Persian).
Tarahi, M., Sabzevari, T., Fattahi, M.H., Derikvand, T., 2022. Estimating runoff in ungauged catchments by Nash-GIUH model using image processing and fractal analysis. Stochas. Environ. Res. Risk Assess. 36(1), 51-66.
Troch, P.A., van Loon, A.H., Hilberts, A.G.J., 2002. Analytical solutions to a hillslope storage kinematic wave equation for subsurface flow. Adv. Water Resour. 25(6), 637-649.