انجمن هیدرولیک ایران نشریه علمی هیدرولیک 2345-4237 16 1 2021 04 21 Uncertainty analysis of the HEC-RAS results in hydraulic simulation of Karoon River flow by Monte-Carlo approach تحلیل عدم قطعیت نتایج مدل HEC-RAS در شبیه‌سازی پارامترهای هیدرولیکی جریان رودخانه کارون با رویکرد مونت کارلو 1 22 129403 10.30482/jhyd.2021.253266.1483 FA حسین ریاحی مدوار 1- استادیار گروه علوم و مهندسی آب، دانشکده کشاورزی، دانشگاه ولی‌عصر (عج) رفسنجان،ایران 0000-0002-5902-4985 بهمن فکوری دانشجوی دکتری سازه‌های آبی، دانشکده کشاورزی، دانشگاه تربیت مدرس، تهران Journal Article 2020 10 18 Introduction <br /> The roughness parameter in hydraulic modelling of natural river and channel flows is not measurable easily and accurate point determination of roughness coefficient, its spatial and temporal variations includes several uncertainties that acts as the main source of error and uncertainties in hydraulic modeling. These drawbacks restricts the applicability of hydraulic modelling in river engineering projects, flood control and management, re-habitation and river restoration. Due to the uncertainties in rating curve, stage, top water width, stream power, shear stress, Froude number and velocity. Because of these drawbacks, in the current study uncertainty analysis by Monte-Carlo Simulation (MCS) combined with HEC-RAS model.<br /> <br /> Methodology <br /> In the current study uncertainty analysis, by Monte-Carlo Simulation (MCS) combined with HEC-RAS model is used to study a 105 km reach of Karoon River from Mollasani to Farsiat as shown in Fig. 1. The model is calibrated and verified using two year daily data of river flow and stage levels in Ahvaz station at the middle of the river reach. A computational control module is developed and combined with computational core of HEC-RAS to perform MCS automatically and the flowchart of modeling strategy and uncertainty analysis is presented in Fig.2. The MCS approach is coupled with computational core of HEC-RAS model by developing a subprogram that create and modifies the input files of HEC-RAS, run it automatically based on random samples of n Manning, and extracting the results of HEC-RAS model in each execution for further analysis in an automatic procedure. By using probability distribution of Manning roughness, 3000 simulations performed and graphical and quantitative indices used to evaluate the uncertainties of model results. In order to refine proper MCS from non-proper ones, the NSE>0.75 index is used to objectively sample n Manning from uncertainty analysis. The uncertainty analysis of proper MCS evaluated by 5 and 95% uncertainty bounds. The uncertainty analysis of model results are evaluated based on the six parameters of water surface elevation, top width of water, flow velocity, Froude number, stream power and shear stress in 3000 runs of peak flow and mean flow discharges respectively and quantified by two indices of 95PPU and d-factor.<br /> Results and Discussion <br /> The calibration and verification results of the HEC-RAS model in Figs 3-4 shows that in the calibration data set the R2 and RMSE of model in discharge are 0.94 and 21 (m3/s); and 093 and 0.6 (m) for water stage respectively. These values in the verification stage were 0.94 and 25.2(m3/s) for discharge; and 0.91 and 0.1(m) for water stage respectively. The results in 105 km length of Karoon River reveals high level of uncertainties with d-factor greater than 1 up to 11 in peak discharge of 3000 and mean daily discharge of 457 m3/s. These results revealed that using conditional evaluations based on NSE>0.75 reduced the uncertainty of d-factor in results of rating curve, stage, top water width, stream power, shear stress, Froude number and velocity. The d-factor of water stage reduced from 2 to 0.07 in peak discharge, and from 0.96 to 0.02 in average flows. These uncertainty reductions in top width of water were 2.5 to 0.19 in peak discharge and 1.3 to 0.078 in average flows of Karoon river. The highest uncertainty of HEC-RAS model results observed in water velocity and Froude number with di-factor 10.85 and 7.44 in peak discharge respectively. This trend of uncertainty reduction observed for water velocity, Froude number, stream power and shear stress along the river, as provided in Tables 1-3. The spectral responses of hydraulic parameters in model result that presented in Figs 5-10, indicate that although the HEC-RAS model produced high uncertainty values, especially in the complex domain of Karoon river, but these uncertainties dos not deviates the hydraulic patters of river flow in the study reach. The peak and maximum values and the zones of high vales of parameters, show high level of uncertainty than the small or moderate hydraulic situations. These indicates the inherent uncertainty in model results that causes high extents of spectral responses for model simulations. The provided findings necessitates the accurate determination of roughness coefficient according to its spatial and dynamic variations along the river reach. <br /> Conclusion <br /> The uncertainty results revealed high level of latent uncertainties in HEC-RAS model results and probabilistic analysis of models results is required for river re-habitation and management practices of large rivers such as Karoon River to provide certain and reliable results. The presented methodology and framer in the current study that uses automatic control and automation of HEC-RAS runs, strengths the modeling capability of one dimensional river flows for probabilistic analysis and automatic calibration of this mode. پارامتر ضریب زبری در مدل‌سازی هیدرولیکی جریان رودخانه‌ها به سادگی قابل اندازه‌گیری نیست و تعیین آن همواره با عدم قطعیت و خطا در نتایج همراه است. به همین منظور در این مقاله از رویکرد شبیه‌سازی مونت-کارلو برای تحلیل عدم قطعیت نتایج مدل هیدرولیکی HEC-RAS در بازه ملاثانی تا فارسیات رودخانه کارون به طول 105 کیلومتر استفاده شده است. با توسعه یک ماژول محاسباتی کنترلی و ترکیب آن با هسته محاسباتی HEC-RAS اجرای خودکار فرآیند مونت-کارلو فراهم شد. شبیه‌سازی3۰۰۰ نمونه مونت کارلو براساس توزیع احتمال ضریب مانینگ انجام شد و تحلیل گرافیکی و کمی‌سازی نتایج عدم قطعیت روی پارامترهای هیدرولیکی خروجی مدل صورت گرفت. نتایج نشان‌دهنده عدم قطعیت زیاد با پهنای باند اطمینان بزرگتر از 1 تا 11، در دبی حداکثر سیلاب 3000 و دبی متوسط روزانه 457 مترمکعب بر ثانیه است. برای پالایش شبیه‌سازی‌های مونت-کارلوی کارآمد و غیرکارآمد از معیار شاخص NSE>0.75 در تحلیل هدفمند عدم قطعیت زبری بر مبنای کران‌های عدم قطعیت 5 و 95 درصد استفاده شد. در این حالت پهنای باند عدم قطعیت (d-factor ) هر شش پارامتر تراز سطح آب، عرض سطح آب، عدد فرود، سرعت جریان، تنش برشی و توان جریان، بعنوان متغیرهای پاسخ، کمتر از 1 شد که نشان دهنده کارآمدی رویکرد به‌گزینی ضریب زبری است. نتایج تحلیل عدم قطعیت نشان داد عدم قطعیت‌های پنهان در نتایج مدل HEC-RAS بالا است و در صورت تحلیل احتمالاتی روی نتایج می‌توان در مطالعات بهسازی، لایروبی و احیای رودخانه‌هایی همچون کارون به نتایج اطمینان و اعتماد بالاتری بخشید و پهنه‌های سطوح سیل‌گیری احتمالاتی را استخراج نمود. https://jhyd.iha.ir/article_129403_e94e45271ce3d4aaacff9955f1551c8c.pdf