فهرست مطالب
مجله پژوهش آب ایران
پیاپی 22 (پاییز 1395)
- تاریخ انتشار: 1395/09/27
- تعداد عناوین: 20
-
- مقالات پژوهشی
-
صفحه 1در این پژوهش اثر زبری بستر روی توزیع سرعت در بدنه جریان غلیظ در دو بخش جت و دیواره بررسی شده است. بدین منظور آزمایش هایی به صورت جریان غلیظ نمکی با 3 دبی ورودی 7 /0، 1 و 3 /1 لیتر بر ثانیه در محدوده شیب 0 تا 2 /2 درصد انجام شد. برای زبر کردن بستر از ذرات شن طبیعی با اندازه های 4، 8، 12 و 15 میلی متر (که در کف کانال چسبانده شد) استفاده شد. همچنین برای بررسی اثر تغییرات غلظت جریان ورودی بر مشخصات جریان غلیظ، از سه غلظت 10، 16 و 20 گرم بر لیتر استفاده شد. توزیع سرعت در بدنه و بررسی روند تغییرات آن در زبری و شیب های مختلف با یک سرعت سنج پروفایلر آلتراسونیک اندازه گیری شد. نتایج نشان می دهد که ارتفاع سرعت حداکثر در پروفیل سرعت بدنه با افزایش زبری کف افزایش می یابد. نسبت ارتفاع سرعت حداکثر به ضخامت بدنه جریان غلیظ زمانی که زبری نسبی (نسبت ارتفاع زبری به ضخامت بدنه جریان) ، به سمت صفر میل می کند، برابر مقدار ثابت 3 /0 خواهد بود و با افزایش زبری نسبی ، ارتفاع نسبی محل سرعت حداکثر ، به سمت عددی ثابت 65 /0 همگرا می شود. نتایج نشان می دهد که زبری عامل مهمی در توزیع سرعت در بدنه جریان غلیظ است. بررسی تشابه توزیع سرعت در ناحیه جت و دیواره، به ترتیب با توزیع های نیمه گوسی و توانی، نشان می دهد که ضرایب توزیع سرعت در بدنه با افزایش زبری بستر تغییر می کند.کلیدواژگان: زبری بستر، توزیع سرعت، جریان غلیظ، ناحیه جت، دیواره
-
صفحه 11تولید پتانسیل محصول زراعی از پارامترهایی است که با محاسبه آن می توان خلاء عملکرد را به دست آورد و بر اساس عوامل تاثیرگذار در تولید، به مدیریت منابع آب و خاک پرداخت. در این پژوهش با هدف محاسبه تولید پتانسیل گندم آبی در 30 مرکز استانی کشور، ابتدا با استفاده از تحلیل عاملی و خوشه ایو با در نظر گرفتن 9 پارامتر اقلیمی، مراکز استانی در چهار گروه طبقه بندی شدند. سپس با روش ناحیه بندی بر اساس بوم شناسی زراعی، تولید پتانسیل محصول و خلاء عملکرد تعیین شد. بر اساس نتایج به دست آمده، روش مذکور در مناطق گرم خشک و کم بارش کشور نتوانست تولید پتانسیل محصول را برآورد کند؛ اما در مناطق نیمه خشک تا مرطوب کشور، خلاء عملکرد مشاهده شد. همچنین نتایج نشان داد که مدل از بین عوامل اقلیمی، بیشترین حساسیت را نسبت به میانگین درجه حرارت دارد. این امر سبب شده تا کمترین خلاء عملکرد در گروه های اقلیمی مرطوب کم بارش و خشک با 15- و 5 /27- درصد دیده شود. همچنین وجود اقلیم معتدله در گروه اقلیمی پربارش کشور سبب رشد بیشتر گونه ها و طولانی تر شدن فصل رشد شده که این امر سبب اختلاف بین عملکرد متوسط و عملکرد پتانسیل شده است و خلاء عملکرد در استان های ساحلی شمالی کشور را به 8 /55- درصد رسانده است.کلیدواژگان: تغییرات اقلیمی، تحلیل عاملی و خوشه ای، داده های هواشناسی، راندمان محصول گندم
-
صفحه 21آب پاک یکی از عوامل مهم توسعه هر منطقه است. با توجه به قرارگیری ایران در منطقه گرم و خشک و کمبود منابع آب، حفاظت و تامین کیفیت آب لازم برای مصارف مختلف اهمیتی دو چندان دارد. به طور معمول ارزیابی کیفی آب های سطحی پرهزینه و زمان بر بوده و انتخاب روشی که در آن با حداقل پارامترهای هیدروشیمیایی بتوان پیش بینی به نسبت دقیقی از کیفیت آب داشت، ترجیح داده میشود. یکی از مهم ترین پارامترهای کیفی آب در زمینه فعالیت های کشاورزی، نسبت جذبی سدیم (SAR) است که تخمین و ارزیابی دقیق مقدار آن بسیار ضروری است. در این بررسی، امکان سنجی تخمین شاخص کیفی SAR در رودخانه آجی چای در منطقه آذربایجان شرقی با استفاده از پارامترهای هیدروشیمیایی مختلف با مدل درختی قوانین M5 و ماشین بردار پشتیبان بررسی شد. برای بررسی دقت مدل های M5 و ماشین بردار پشتیبان از چهار آماره ضریب همبستگی (R) ، نش- ساتکلیف (NSC) ، جذر میانگین مربعات خطا (RMSE) و میانگین خطای مطلق مقادیر (MAE) استفاده شد. مقادیر این آماره ها برای روش ماشین بردار پشتیبان (98 /0R=، 97 /0N-SC=، (mg/l) 22 /6RMSE= و (mg/l) 06 /6MAE=) و برای مدل M5(98 /0R=، 96 /0N-SC=، (mg/l) 33 /7RMSE= و (mg/l) 9 /3MAE=) محاسبه شد. نتایج مقایسه نشان داد هر دو روش عملکرد خوبی در تخمین میزان SAR داشته اند، اما مدل درختی قوانین M5 در محدوده داده های مورد استفاده روابط خطی ساده و کاربردی تر ارائه می کند.
چکیده (انگلیسی):
Clean water is one of the important factors in any region''s development. Since Iran is located in an arid and semi-arid area with scarce water resources، preservation of water required for various uses and maintenance of its quality takes redoubles this importance. Evaluation of surface water is normally a costly and time-consuming process. Therefore، a method is preferred which has the minimum number of hydrochemical parameters and can yield a relatively accurate prediction of water quality. One of the most significant qualitative parameters of water for agricultural uses is the sodium absorption ratio (SAR)، the factor which should be estimated and evaluated accurately. This research employed various hydrochemical parameters، a model tree using the M5-Rules، and a Support Vector Machine to study the feasibility of estimating the qualitative index SAR in the Ajichai River located in East Azerbaijan Province. The four statistics of correlation coefficient (R)، Nash-Sutcliffe coefficient (NSC)، Root Mean Square Error (RMSE)، and Mean Absolute Error (MAE) were used to determine the accuracy of both M5 model and the Support Vector Machine.
The study region was the Ajichai River on the northern hillsides of the Sahand Mountain. Hydrochemical data from the Hydrometric Station in Vanyar was used to evaluate and predict the SAR in the river. The Vanyar Station has the longitude of 46 ̊24 east، the latitude of 38 7́ north، and the altitude of 1460 meters. Effects of Total Dissolved Solids (TDS)، Electrical Conductivity (EC)، PH، chlorine (Cl-)، sulfate (SO42+)، calcium (Ca2+)، magnesium ( Mg2+) and sodium ( Na+) parameters on SAR were determined in SAR estimation. The model tree M5-Rules is a new data mining method. The main goal of this model is derived from regression trees. The difference is that this model has regression functions in its leaves instead of constant values and classification tags. The major advantage of the model tree M5-Rules over regression trees is that the model tree M5-Rules is much smaller than regression trees. Furthermore، regression functions normally do not include many parameters. A decision tree usually consists of four parts of root، branches، nodes، and leaves. Each node corresponds to a certain characteristic، and the branches represent values of the intervals. These intervals consider known values for each of the characters. The branching operation takes place with one of the predictor variables. The branching intervals are selected in a way that the sum of squared deviations from the mean of the data in each node is minimized. The branching criterion indicates the amount of the error in the related node، and the model calculates the minimum expected error as a result of each characteristic testing in the related node. The model error is generally assessed by measuring the predicted unobserved target values accuracy. In this research، the WEKA software which is developed at Waikato University in New Zealand was used to model the M5 method. Modeling was performed with this software using the option of M5-Rules which present simple and linear rules. Support Vector Machines are data mining algorithms similar to the model tree M5 and the artificial neural network. There are two groups of Support Vector Machines: Support Vector Classification (SVC) and Support Vector Regression. Furthermore، Support Vector Machines are based on the concept of decision planes that define decision boundaries، i.e. a decision plane separates data with different tags from each other. The goal in a linearization algorithm with the help of a Support Vector Machine، the assumptions of the input value of xi، and the output value of yi is to find a function with the minimum deviation (ε) from the yis (ε is the amount of deviation). In this research، the Statistica software is used for modeling the SAR values employing Support Vector Regression.
In the modeling of the SAR values by using the tree model M5-Rules، the best answer was obtained when 66 percent of the data was allocated to training and the rest to testing. To model the SAR values using the Support Vector Machine، various functions were tested as kernel functions، and it was found that the RBF function exhibited the best performance in the modeling of the SAR values. Among the 10 scenarios studied in this research، the best one was selected. The four statistics of correlation coefficient (R)، Nash-Sutcliffe coefficient (NSC)، Root Mean Square Error (RMSE)، and Mean Absolute Error (MAE) were used to determine the accuracy of both M5 model and the Support Vector Machine. The obtained values of these calculations were: R =0.98، N-SC=0.97، RMSE=6.22 (mg/l)، MAE=6.06 (mg/l) for the Support Vector Machine method; and R=0.98، N-SC=0.96، RMSE=7.33 (mg/l)، and MAE=3.9 (mg/l) for the M5 model. Results of the comparison indicated that both methods studies in this work، i.e. Support Vector Regression and the M5 model، were highly capable of predicting the SAR values in the Ajichai River، using the available data. However، the M5 model is recommended to be used due to the fact that the formulas employed in this method are so simple and linear.کلیدواژگان: کیفیت آب، مدل درختی قوانین M5، نسبت جذبی سدیم، _ ماشین بردار پشتیبان -
صفحه 31با توجه به اهمیت انتقال سیالات حیاتی از جمله آب در دنیای امروز، تلاش برای یافتن راهی که بتواند با کمترین هزینه سبب کاهش مصرف انرژی در هنگام انتقال سیال شود، مهم است. از جمله مهم ترین و ساده ترین راه ها، استفاده از مواد کاهنده اصطکاک در جریان است. هدف از این پژوهش، افزودن مقادیر کم از سه پلیمر کارآاکریل30، کارآاکریل 58 و سدیم کربوکسی متیل سلولز در خطوط لوله سیستم آبیاری بارانی برای بررسی تغییرات میزان تلفات هد ناشی از اصطکاک درخطوط لوله ها و تاثیر آن بر شعاع پاشش آبپاش ها است، در این آزمایش کارآاکریل30 و 58 در غلظت های 500،100 و 1000 و سدیم کربوکسی متیل سلولز در غلظت های 50، 100 و 150 میلی گرم بر لیتر به آب داخل مخزن اضافه و بعد از حل شدن، با پمپ وارد خط لوله و از طریق آبپاش ها پاشش انجام شده است. پلیمرهای سدیم کربوکسی متیل سلولز و کارآاکریل30 عملکرد بهتری نسبت به کارآاکریل 58 نشان دادند. مقدار افت در دبی 423 /0 لیتر بر ثانیه برای آب خالص، کارآاکریل30 در غلظت 1000 میلی گرم بر لیتر، کارآاکریل 58 در غلظت 1000 میلی گرم بر لیتر و سدیم کربوکسی متیل سلولز در غلظت 150 میلی گرم بر لیتر به ترتیب 04 /20، 8 /12، 09 /16 و 39 /12 متر در لوله ای به طول 40 متر و به قطر 20 میلی متر اندازه گیری شد.کلیدواژگان: کاهش دراگ، شعاع پاشش آبپاش، ماده کاهنده دراگ، _ خطوط لوله، پلیمر
-
صفحه 39اگرچه تاکنون پژوهش های زیادی در رابطه با پرش در سطوح صاف شیب دار انجام شده اما بررسی ای در سطوح شیب دار با بستر زبر انجام نشده است. هدف اصلی این پژوهش ارائه رابطه ای کلی برای محاسبه نسبت عمق های مزدوج پرش در سطوح شیب دار با بستر زبر بود. بدین منظور انواع مختلف پرش در سطوح شیب دار در دو حالت بستر صاف و زبر در محدوده وسیعی از زبری نسبی (032 /0-546 /0)، شیب های مختلف تندآب (5 /14، 5 /20 و 5 /27 درجه) و عدد فرود (3 /1 تا 4 /7) به صورت آزمایشگاهی بررسی شد. با استفاده از نظریه باکینگهام تحلیل ابعادی، پارامتر های بی بعد موثر بر نسبت عمق های مزدوج پرش در سطوح شیب دار به دست آمد و با استفاده از تئوری خود- تشابهی ناقص، رابطه ای کلی استخراج شد. سپس، ضریب های معادله کلی استخراج شد و دقت آن با استفاده از داده های آزمایشگاهی این بررسی و منابع منتشر شده قبلی ارزیابی شد. نتایج بررسی ها نشان داد که معادله دقت بالایی داشته و درصد متوسط خطای آن 5 /7 درصد است. همچنین نتایج نشان داد که افزایش زبری نسبی سبب کاهش نسبت عمق های مزدوج و عدد فرود اولیه جریان می شود به گونه ای که در زبری نسبی 54 /0 کاهش 50 درصدی عدد فرود و نسبت عمق های مزدوج را در پی داشت.کلیدواژگان: پرش هیدرولیکی، خود تشابهی ناقص، حوضچه آرامش، زبری، تندآب، استهلاک انرژی
-
صفحه 49مطالعات مربوط به شکست سد به ویژه موج سیلاب ناشی از آن، برای مدیریت ایمنی پایین دست سدها ضروری است. دانستن سرعت امواج ناشی از شکست سد و اثرات زبری بستر در تغییرات کمی و کیفی گسترش سیلاب در پایاب سدها در هیدرولیک جریان های با سطح آزاد مهم است. پدیده شکست سد با نصب دریچه در مسیر کانال آزمایشگاهی و باز کردن ناگهانی دریچه قابل بررسی است. در این پژوهش با به کارگیری نرم افزار Flow-3D به مدل سازی سه بعدی بررسی مشخصات هیدرولیکی سیلاب ناشی از شکست سد و تاثیر پذیری موج شکل گرفته به دلیل وجود زبری های مصنوعی، نحوه آرایش آن ها، تاثیر عمق پایین دست، و مقایسه آن ها با نتایج آزمایشگاهی پرداخته شده است. نتایج نشان داد که وجود زبری در کف سبب افزایش مقاومت در برابر جریان شده و از سرعت موج پیش رونده می کاهد که این کاهش سرعت در آرایش زیگزاگ زبری ها بیشتر از آرایش منظم است که دلیل آن تفاوت در فاصله عرضی بین زبری ها در دو نوع آرایش است که نحوه پخش جریان را متفاوت می کند. به طور کلی وجود زبری با آرایش زیگزاگ در کنار افزایش عمق پایین دست بیشترین تاثیر را در استهلاک سرعت موج پیش رونده دارد.کلیدواژگان: مدل سازی عددی، شکست سد، سیلاب، مهندسی رودخانه
-
صفحه 57شاخص های مستخرج از منحنی تداوم جریان، همواره یکی از موضوعات مهم در هیدرولوژی بوده و اطلاع از میزان آن، نقش به سزایی در برنامه ریزی و مدیریت بهینه منابع آب دارد. در این پژوهش، با استفاده از داده های جریان روزانه رودخانه، شاخص های تداوم جریان در هجده ایستگاه آب سنجی حوضه خزری، تعیین شد. در محیط GIS پارامترهای فیزیوگرافی زیرحوضه ها استخراج و عوامل مختلف اقلیمی، زمین شناسی و کاربری اراضی نیز محاسبه شد. روش آزمون تحلیل عاملی از بین بیست ویک پارامتر موثر، شش عامل به عنوان مهم ترین عوامل مستقل انتخاب شدند. سپس روابط رگرسیون تک و چند متغیره خطی، توانی، نمایی در سطح معنی داری 1 درصد تعیین شدند. در روش رگرسیون های چندمتغیره از بین عامل های زمین شناسی به طور مشخص سازندهای سخت کربناته و آبرفتی، از بین عوامل مورفومتری، طول رودخانه اصلی و شیب متوسط حوضه و از جنبه کاربری اراضی نیز عامل سطح کاربری کشاورزی، به عنوان بهترین برآورد کننده ی شاخص تداوم جریان، تشخیص داده شدند. همچنین با توجه به مدل های تک متغیره توانی و نمایی، پارامتر بارش متوسط حوضه به عنوان موثرترین عامل اقلیمی در تخمین شاخص های تداوم جریان رودخانه ها تشخیص داده شد. نتایج اعتبار سنجی روابط رگرسیونی نشان داد که توزیع خطا ها در تمامی مدل های استخراج شده، نرمال است. با توجه به روابط منطقه ای به دست آمده، مشارکت عامل سازندهای سخت کربناته و نیز همراهی آن ها با نهشته های کواترنری، نشان از نقش موثر آن ها در تامین تداوم جریان رودخانه های منطقه دارد.کلیدواژگان: شاخص های تداوم جریان، مدل منطقه ای، آبرفت، کربنات ها، سنگ شناسی
-
صفحه 69مشخص بودن رابطه هیدرولیکی آب های سطحی و زیرزمینی در مدیریت صحیح منابع آب نقش بسیار مهمی دارد و باید در برنامه ریزی ها درنظر گرفته شود. یکی از بهترین روش ها برای بررسی ارتباط آب های زیرزمینی و سطحی، مدل سازی آبخوان ها است. در بیشتر پژوهش های انجام شده به ویژه در داخل ایران، نقش رودخانه ها و آب های سطحی در مدل سازی آب زیرزمینی نادیده گرفته می شود. از این رو در بیشتر مناطق رابطه بین آب های سطحی و زیرزمینی همچنان نامشخص باقیمانده است. در این پژوهش به بررسی رابطه بین آب سطحی و زیرزمینی در دشت دزفول- اندیمشک پرداخته شده است. بدین منظور، مدل آب زیرزمینی منطقه در طول شش سال 91-1385 به کمک نرم افزارGMS تهیه شد. نتایج نشان دهنده تفاوت رابطه هیدرولیکی آبخوان و رودخانه ها در نقاط مختلف دشت است. بیشترین تبادل هیدرولیکی با آبخوان در رودخانه دز در بازه بین سد تنظیمی دز تا محل ورود رودخانه بالارود دیده می شود. به طور کلی رودخانه ها در منطقه بیشتر تغذیه کننده آبخوان هستند. تنظیمی بودن جریان رودخانه ها تاثیر زیادی بر ارتباط هیدرولیکی دارد. به طوری که رودخانه های تنظیمی در فصل خشک، به مقدار بیشتر و در فصل تر به میزان کمتری تغذیه کننده آبخوان هستند و نسبت به رودخانه های طبیعی، وضعیت متفاوتی رخ می دهد.کلیدواژگان: نرم افزار GMS، _ رابطه هیدرولیکی، مدل آب زیرزمینی، رودخانه دز، جریان تنظیمی رودخانه
-
صفحه 79در این بررسی اثر تنش آبی دوره ای و شوری ناشی از آب دریای خزر بر عملکرد پنبه به صورت آزمایش فاکتوریل در قالب طرح بلوک کامل تصادفی مشتمل بر 20 تیمار در سه تکرار در بهار 1391 در مزرعه تحقیقات کشاورزی دانشگاه علوم کشاورزی و منابع طبیعی گرگان انجام شد. تنش آبی دوره ای در 4 سطح شامل I1 تنش اولیه، I2 تنش توسعه مرحله اول، I3 تنش توسعه مرحله دوم، I4 تنش مرحله میانی رشد پنبه و تنش شوری در 5 سطح S1، S2، S3، S4 و S5 به ترتیب اختلاط 0، 25، 50، 75 و 100 درصد آب دریای خزر با آب معمولی با شوری 7 /0، 1 /4،1 /11، 1 /18 و 4 /25 اعمال شد. نتایج نشان داد تنش آبی دوره ای و شوری ناشی از آب دریای خزر بر ارتفاع بوته، وزن غوزه، عملکرد وش در سطح یک درصد و بر تعداد غوزه در سطح پنج درصد معنی دار است. به دلیل این که اختلاف عملکرد وش در اختلاط 25% آب دریا و تنش دوره اولیه مرحله توسعه (I2S2)، با تیمار شاهد نسبت به دیگر تیمارها کمتربود به عنوان یک گزینه مناسب برای استان های شمالی توصیه می شود.کلیدواژگان: تنش خشکی، دریای خزر، آب شور، وش
-
صفحه 87پرش هیدرولیکی یکی از مهم ترین پدیده ها در جریان متغیر سریع است که در نتیجه تغییر حالت جریان از فوق بحرانی به زیربحرانی در جریان های با سطح آزاد ایجاد می شود و با پدید آمدن آن عمق جریان در مسیر به نسبت کوتاهی به میزان قابل توجهی افزایش می یابد. در این پژوهش خصوصیات پرش هیدرولیکی تشکیل یافته در تبدیل های همگرا و واگرا با دیواره های مستقیم و انحنادار به ازای پنج عدد فرود مختلف در محدوده 1 /9-8 /5 بررسی و مقایسه شد. بدین منظور مقادیر عمق و سرعت نقطه ای در نقاط مختلف پرش هیدرولیکی در تبدیل ها اندازه گیری شد. میزان استهلاک انرژی پرش در تبدیل های واگرا با دیواره مستقیم و انحنادار به ترتیب 95 /24 و 23 /21 درصد و مقدار متناظر آن در تبدیل های همگرا به ترتیب 93 /15 و 49 /12 درصد حاصل شد. همچنین نتایج نشان داد که با افزایش فاصله طولی از ابتدای پرش مقادیر حداکثر تنش های رینولدز در تبدیل همگرا افزایش و در تبدیل واگرا کاهش می یابد.کلیدواژگان: پرش هیدرولیکی، تنش های رینولدز، تبدیل همگرا، تبدیل واگرا، استهلاک انرژی
-
صفحه 97پایش کیفیت منابع آب برای توسعه کشاورزی و صنعتی ضروری است. هدف این پژوهش ارزیابی تغییرات مکانی و زمانی غلظت عناصر سنگین در آب و رسوبات بستر رودخانه زاینده رود است. برای رسیدن به این هدف 10 ایستگاه نمونه برداری بین سد زاینده رود و زرین شهر انتخاب شده است. شش نمونه مرکب آب و رسوب از تیرماه 1390 تا اردیبهشت 1391 هر دو ماه یکبار برداشت شدند. نمونه های رسوب پس از تیمار با اسید عصاره گیری و غلظت عناصر سنگین در عصاره و غلظت عناصر محلول در آب رودخانه به کمک دستگاه جذب اتمی اندازه گیری شدند. نتایج نشان می دهند که حداکثر غلظت های مس، کادمیم، روی و سرب به ترتیب در رسوبات ایستگاه های روستای چلوان، زرین شهر، بند ذخیره ای و پل کله وجود دارند در حالی که حداکثر میزان کادمیم، مس و سرب رسوبات در اردیبهشت ماه و برای روی در آذرماه مشاهده شد. از نظر زمانی حداکثر غلظت مس محلول در آب در دی ماه و برای کادمیم، روی و مس در آذرماه مشاهده شدند. به نظر می رسد کارخانه ذوب آهن و مزارع کشاورزی منابع اصلی عناصر سنگین در حاشیه رودخانه هستند. این پژوهش نشان داد که غلظت عناصر سنگین بررسی شده محلول در آب رودخانه کمتر از آستانه مجاز برای آشامیدن بودند.کلیدواژگان: کیفیت آب، رسوب، عناصر سنگین، رودخانه زاینده رود، آلود
-
صفحه 107نفوذ آب در خاک، از حساس ترین پارامترهای هیدرولیکی موثر بر آبیاری سطحی است که به دلیل تاثیر پذیری از عوامل متعدد خاک و شیوه آبیاری، پیچیدگی زیادی دارد. بررسی اثر ارتفاع آب آبیاری بر نفوذ تجمعی و لحاظ کردن این پارامتر در مدل های نفوذ، منجر به کارآیی بهتر این مدل ها در شرایط کاربری مختلف می شود. هدف در پژوهش حاضر، بررسی تاثیر ارتفاع های مختلف آب روی سطح خاک در فرآیند نفوذ، در برخی کلاس های بافتی خاک بخشی از دشت فومنات در استان گیلان در نظر گرفته شد. بدین منظور 17 نمونه از خاک منطقه ی مورد مطالعه برداشت شد. آزمایش های نفوذ برای سه تیمار ارتفاع آب 3، 5 و 7 سانتی متر در سه تکرار به صورت فاکتوریل در قالب طرح تصادفی و در مدل فیزیکی که برای این منظور ساخته شده بود، انجام شد. برای ارزیابی نتایج، تجزیه واریانس و مقایسه میانگین ها به روش توکی در سطح 5 درصد صورت گرفت. مقایسه تیمار ارتفاع آب در 17 نمونه خاک، نشان داد که این تیمار اثر قابل توجهی را در مقدار نفوذ تجمعی چهار نمونه شامل بافت های Silty Clay Loam و Silty Clay داشت. بررسی سه تیمار زمان نفوذ تجمعی 30، 360 و 600 دقیقه نشان داد که در 14 نمونه از خاک ها بین نفوذ تجمعی در زمان 30 دقیقه با زمان های 360 و 600 اختلاف معنی داری (05 /0p<) وجود دارد که این اختلاف در بافت خاک Clay Loam بیشتر بود.کلیدواژگان: مدت زمان نفوذ، ارتفاع آب آبیاری، آبیاری سطحی، خاک رسی، نفوذ تجمعی
-
صفحه 117در این پژوهش برای تعیین الگوی توزیع زمانی بارش در حوزه آبخیز سنگانه استان خراسان رضوی، ابتدا با بررسی کلیه رگبارهای مشاهداتی در طول دوره ی آماری 2013- 2005، رگبارهای منفرد در تداوم های 5 /0>، 1-5 /0، 3-1، 6-3، 9-6 و 9< ساعته استخراج شده و برای ارایه الگوهای توزیع زمانی رگبارها از روش محاسباتی پیل گریم استفاده شده است. نتایج حاصله از این پژوهش نشان می دهد که بیشتر رگبارهای رخ داده در منطقه سنگانه از نوع چارک اولی بوده و با افزایش زمان رگبار، به تعداد رگبارهای چارک دوم تا چهارم افزوده می شود. علاوه بر این، 61 درصد از رگبارهای فرساینده تداومی کمتر از سه ساعت داشته و اکثرا از نوع چارک دومی و سومی هستند. به منظور ارزیابی آماری داده های توزیع با زمان مشخص نیز از آزمون کای اسکوئر استفاده شده است. نتایج به دست آمده نشان داد که در تداوم های کمتر از 5 /0، 1-5 /0 و بیشتر از 9 ساعت، رابطه معنی داری بین رتبه و چارک وجود داشته و الگوها قابلیت بالایی برای استفاده دارند؛ اما در سایر تداوم ها، الگوها قابل اعتماد نیستند. علاوه بر این، مقایسه الگوهای منطقه ای به دست آمده با الگوهای تیپ SCS نیز نشان داد که تنها الگوی توزیع زمانی بارش های چارک دومی با تیپ IA ضریب همبستگی بالایی دارند؛ بنابراین با توجه به این که بیشتر رگبارهای منطقه سنگانه، از نوع چارک اولی و با تداوم 3-1 ساعته هستندکلیدواژگان: آزمون کای اسکوئر، روش پیل گریم، الگوی بارش، حوزه آبخیز سنگانه، چارک، SCS
-
صفحه 127امروزه، آلودگی محیط متخلخل اشباع و غیراشباع، از جمله جدی ترین مشکلات زیست محیطی در سراسر جهان است. از مهم ترین این آلاینده ها، آلاینده های نفتی و ترکیبات آن است که در اثر نشت؛ تراوش هنگام ذخیره سازی یا انتقال، با نفوذ به خاک، منابع آب زیرزمینی را آلوده می کند. ترکیبات بنزن، تولوئن، اتیل بنزن و زایلن (BTEX) ، از معمولترین منابع آلاینده های نفتی در محیط متخلخل زیرزمین هستند. درک دقیق شیوه حرکت و انتشار پلوم آلودگی در محیط متخلخل، نقش مهمی در موفقیت و انتخاب تصمیم های صحیح برای پالایش آلودگی محل های آلوده دارد. در این پژوهش، در محیط آزمایشگاه یک تانک دوبعدی به منظور شبیه سازی فیزیکی حرکت BTEX طراحی شد و توزیع درجه اشباع پلوم آلاینده در محیط متخلخل با استفاده از روش پردازش تصاویر تعیین گردید. مقایسه نتایج حاصل از روش پردازش تصاویر، حاکی از مطابقت مقدار درجه اشباع با تصاویر گرفته شده از پلوم رنگی BTEX بود. همچنین بررسی مقدار درجه اشباع پلوم در زمان های مختلف از شروع تزریق، نشان داد که در زمان های اولیه به علت زیاد بودن سهم نیروی گرانشی در گرادیان کل، حرکت پلوم به صورت عمودی است. پس از رسیدن پلوم BTEX به ناحیه موئینگی، به دلیل کاهش کشش سطحی بین هوا و آب در این ناحیه، پلوم BTEX در بالای سطح ایستابی متوقف می شود و در ناحیه موئینگی به صورت افزایشی تجمع می یابد.کلیدواژگان: محیط متخلخل، غیراشباع، BTEX، پردازش تصاویر
-
صفحه 135از فاکتورهای مهم در برآورد بیلان آبی یک حوضه و برنامه ریزی و مدیریت منابع آب آن، تعیین پتانسیل تغذیه آب زیرزمینی (Groundwater Recharge Potential، GRP) است. در این زمینه، روش های تهیه نقشه های موضوعی و تلفیق آن ها در محیط GIS به منظور تهیه یک نقشه پهنه بندی شده GRP در کلاس های مختلف، می تواند در صرفه جویی زمان و منابع مالی موثر باشد. در این مطالعه، نقشه های موضوعی پارامترهای هفت گانه تراکم خطواره ها، تراکم طولی شبکه آبراهه ها، لیتولوژی، زمین های کارستی، کاربری اراضی- پوشش زمین، شیب زمین و بارش برای حوضه آبریز رودخانه کر تهیه شد. سپس با در نظر گرفتن بر هم کنش این پارامترها در تغذیه حوضه، نقشه پهنه بندی پتانسیل تغذیه در پنج کلاس تغذیه خیلی بالا تا خیلی پایین با استفاده از وزن نهایی پارامترها تهیه گردید. در نهایت با مقایسه ای که بین برآورد تغذیه با استفاده از روش مورد استفاده در این مقاله با مقدار تغذیه برآوردی با روش بیلان هیدروکلیماتولوژی (تورنت وایت) صورت گرفت، مشخص شد حدود 14 درصد اختلاف بین مقدار تخمینی با هر یک از این روش ها وجود دارد.کلیدواژگان: GRP، روش ارجاع دادن، حوضه آبریز رودخانه کر، بیلان هیدروکلیماتولوژی
-
صفحه 147استفاده از روش های اقتصادی برای جلوگیری یا کاهش خطر تخریب پل ها در زمان سیلابی شدن رودخانه ها، امری ضروری است. از جمله این روش ها می توان به استفاده از طوقه، شکاف، سنگ چین، کابل پیچیده شده دور پایه و یا ترکیب آن ها اشاره کرد. در این پژوهش تاثیر استفاده همزمان از کابل و طوقه به عنوان روشی جدید برای کاهش آب شستگی در کنار گروه پایه های سه تایی با فاصله بین پایه دو برابر قطر پایه (2=S/b) و سه برابر آن (3=S/b) ، بررسی شد. آزمایش ها در شرایط آب شستگی آب زلال برای رسوبات یکنواخت انجام شد. از طوقه با قطر دو برابر قطر پایه و کابل با قطر 15 /0 قطر پایه و زاویه پیچش 15 درجه استفاده شد. نتایج نشان داد که استفاده هم زمان از کابل و طوقه کارایی بسیار بهتری نسبت به کابل و طوقه تنها دارد به طوری که سبب کاهش 20، 45 و 53 درصدی عمق نهایی آب شستگی به ترتیب در کنار پایه های اول، دوم و سوم می شود. همچنین مشاهده شد که با افزایش فاصله بین پایه ها زمان شروع آب شستگی با تاخیر بیشتری همراه است به طوری که 234 دقیقه تاخیر برای پایه سوم با 3=S/b ثبت شد.کلیدواژگان: کنترل آب شستگی موضعی، کابل، عمق آب شستگی موضعی، گروه پایه، طوقه
- یادداشت های فنی
-
صفحه 157تاکنون برای محاسبه هیدروگراف واحدلحظه ای ، روابط مختلفی ارائه شده است. یکی از روش های محاسبه و پیش بینی سیلاب خروجی از حوضه، استفاده از مفهوم مخازن متوالی است. گامای معکوس یکی از توزیع های احتمالاتی است که به دو پارامتر و وابسته است و برای تخمین پارامترهای آن روش های مختلفی وجود دارد. دقت تخمین پارامترهای و در مدل گامای معکوس، در بهبود نتایج شبیه سازی جریان تاثیر زیادی دارد. در این مقاله برای برآورد پارامترهای مدل گامای معکوس از روش های گشتاورگیری، روسو و بهینه سازی استفاده شده است. نتایج نشان می دهد که استفاده از روش بهینه سازی، درمحاسبه پارامترهای و بیشترین هماهنگی را با هیدروگراف های مشاهداتی دارد. روش های روسو و گشتاورگیری در مرتبه های بعدی قرار داردند. هر سه روش مذکور، حجم سیلاب خروجی را با خطای کمتر از 11 درصد، محاسبه کرده اند. دبی اوج و زمان دبی اوج به ترتیب با روش های بهینه سازی، روسو و گشتاورگیری با دقت بالایی برآورد شده است. ولی روش گشتاورگیری در 5 رویداد بارندگی، در محاسبه دبی اوج با مقدار مشاهداتی بیش از 50 درصد اختلاف دارد. بنابراین در حوضه آبخیز جنگ، روش بهینه سازی به عنوان روش مناسب در تعیین مقدار و در مدل گامای معکوس معرفی می شود.کلیدواژگان: روش گشتاورگیری، مدل روسو، _ حوضه، روش بهینه سازی، آبریز جنگ، مدل گامای معکوس، برآورد پارامترهای و
-
صفحه 163هدف از این مطالعه ارزیابی مدل SRM در برآورد رواناب روزانه اندازه گیری شده در ایستگاه آب سنجی گلینک می باشد. در این پژوهش ورودی های مدل شامل متغیرهای هیدرولوژی، هواشناسی و وضعیت فیزیکی حوضه آبخیز محاسبه شد و 72 گذر از تصاویر ماهواره MODIS، برای تعیین سطح پوشش برف استفاده شد و بعد از اعمال ورودی ها به مدل، شبیه سازی رواناب انجام گرفت. مقایسه آب نمودها از طریق ضریب کارآیی نش- ساتکلیف و درصد اختلاف حجم دبی انجام شد. مقادیر ضریب کارآیی و درصد اختلاف حجم دبی در ایستگاه گلینک برای دوره واسنجی 95 /0 و 07 /2 درصد و برای دوره صحت سنجی 86 /0 و 81 /4 درصد محاسبه گردید که نشان می دهد مدل با کارآیی بالا و حداقل خطا رواناب حاصل از ذوب برف را برآورد نموده است.کلیدواژگان: رواناب، ذوب برف، تصویر MODIS، مدل سازی
-
صفحه 169در این پژوهش، پتانسیل های بیوفیزیکال محصولات زراعی در محدوده اراضی شبکه آبیاری و زهکشی دشت قزوین ارزیابی شد. پس از ترسیم واحدهای اراضی، شبیه سازی رشد و تولید محصولات زراعی مورد بررسی با مدل های گیاهی بسته نرم افزاری DSSAT انجام شد. در این شبیه سازی برای هر محصول، دو رژیم آبیاری (معرف آبیاری سطحی و تحت فشار) به همراه 2 تا 3 تاریخ کشت مختلف در نظر گرفته شد. تحلیل نتایج نشان از تفاوت عملکرد محصولات در تیمارهای مختلف تاریخ کشت و رژیم های آبیاری بود. در بین محصولات عمده منطقه، مقدار ضریب تغییرات مکانی گندم و جو در سطح متوسط قرار دارد، در حالیکه محصول ذرت دانه ای دارای مقدار کوچک تر ضریب تغییرات است. در ارتباط با محصولاتی مانند گوجه فرنگی، پنبه و سویا، اختلاف زیادی بین عملکرد محصول در واحدهای مختلف اراضی وجود دارد و این مساله باید در تعیین الگوی کشت این مناطق مورد توجه قرار گیرد. شاخص ضریب تغییرات زمانی محصولات نیز می تواند به عنوان شاخصی از ریسک موجود در ارتباط با کاشت این محصولات مورد توجه قرار گیرد. نتایج این بررسی، اطلاعات ورودی لازم برای یک مدل پشتیبان برنامه ریزی مکانی را فراهم می کند.کلیدواژگان: شبیه سازی، ریسک، واحدهای اراضی، _ ایران، برنامه ریزی
-
صفحه 175افزایش کارایی ایستگاه های پمپاژ آب کشاورزی، به دلیل سهم بالا در مصرف انرژی و هزینه، مهم است. تاکنون در ارزیابی شبکه های آبیاری و زهکشی، عملکرد ایستگاه های پمپاژ نادیده گرفته شده است، در حالی که هر اختلالی در ایستگاه پمپاژ، سبب اختلال در آبگیری و در نتیجه عملکرد شبکه خواهد شد. از آنجا که، روش ارزیابی مشخصی برای این تاسیسات وجود ندارد، در این پژوهش روش ارزیابی کمی تحلیل پوششی داده ها برای آنها معرفی و پیشنهاد شده است. در این راستا، پس از معرفی روش و ارائه مدل های آن، کارایی زمانی ایستگاه پمپاژ اصلی سبیلی در سال های 1380 تا 1389، با نرم افزار DEAP2 محاسبه شده است. در این روش، با معرفی شاخص های ارزیابی عملکرد، کارایی متوسط ایستگاه سبیلی، 69 درصد به دست آمد. همچنین بالاترین کارایی برابر 100 درصد در سال های 1380، 1381 و 1389 و کمترین کارایی برابر 21 درصد در سال 1388، محاسبه شد. عمده این ناکارایی، به دلیل هزینه بالای برق مصرفی است، که با توجه به آن در برنامه ریزی و اقدامات کاهش مصرف انرژی، می توان به حداکثر پتانسیل کارایی در بهره برداری سیستم پمپاژ رسید.کلیدواژگان: شبکه آبیاری و زهکشی، کارایی، ارزیابی عملکرد، تحلیل پوششی داده ها، ایستگاه پمپاژ اصلی سبیلی
-
Page 1Gravity currents form when a heavier fluid propagates into a lighter one in a predominantly horizontal direction. They are frequently encountered both in the environment and engineering applications. Gravity currents can be driven by density differences of the fluids involved, or by differential particle loading. In natural or human-made aquatic settings such as lakes, oceans or reservoirs, there is a myriad of possible contributors to these density differences including temperature differences, salinity contrasts, suspended material, both organic and inorganic, as well as combinations of these mechanisms. In water resources management, to prevent sedimentation and drain sediment from dam reservoirs, in most cases, is tried to remove fine sediments by using of hydrodynamic forces. One of these methods is removing sediment by the gravity current dynamics. The gravity currents are the most important effective events on reservoir's sedimentary processes (transmission, distribution and deposition of particles). Therefore, identifying the factors affecting on this type of flow is very important. Both bed slope and surface roughness are the most important parameters on gravity current dynamics and have a considerable effect on velocity distribution of the body.
This study investigates the effect of bed roughness on the gravity currents characteristics and obtaining relations between effective parameter on gravity currents and finding their influences in body velocity distribution on rough and slope beds. In this study to determine the bed roughness effect on the velocity distribution, a series of gravity currents experiments was performed on beds with the rough natural element sizes of 4, 8, 12 and 15 mm and slopes from 0 to 2.2%. The experiments of saline gravity currents were carried out in three inflow discharges with amounts of 0.7, 1 and 1.3 (lit/s) and with three concentrations of 10, 16 and 20 (gr/lit). All Experiments was done on a Plexiglas flume with the length of 10 meters and width of 35 cm and in hydraulic models laboratory in Shahid Chamran university of Ahvaz (Iran). In all experiment, the velocity profile in the body of gravity current was acquired by an Ultrasonic Doppler Velocity meter device.
The results of this study can be summarized as an investigation on variations of velocity profile characteristics and analysis of velocity distribution in tow (wall and jet) regions of gravity current's body. Determination of general equations for velocity distribution in gravity current's body on rough beds showed that by increasing bed roughness the height of maximum velocity rose up to a higher position and the coefficient of velocity distribution equation changed consequently. Based on results of this study, the relationship between relative height of maximum velocity location in the velocity profile ( ), versus relative roughness ( ), was rational type; in which, hm is the height of maximum velocity location, ks is the bed roughness and h is the thickness of gravity current. The relative height of maximum velocity location was equal to 0.3 when was near to zero and by increasing , the value of converged to the constant of 0.65. The investigation on acquired velocity profile showed that the bed roughness was one of the most important factors on velocity distribution. The similarity of the velocity distribution, in both jet and wall flow regions, by Semi-Gaussian and power distribution equations showed that the coefficient of velocity distribution was varied by increasing the bed roughness.
For the velocity profile in jet region, the coefficients of Semi-Gaussian distribution (α and β) were varied from 1.2 to 2 and 2.2 to 1.8, respectively, by increasing the bed roughness. Also, for the velocity profile in wall region, the coefficient of Power distribution (n) was changed from 0.2 to 0.6.Gravity currents form when a heavier fluid propagates into a lighter one in a predominantly horizontal direction. They are frequently encountered both in the environment and engineering applications. Gravity currents can be driven by density differences of the fluids involved, or by differential particle loading. In natural or human-made aquatic settings such as lakes, oceans or reservoirs, there is a myriad of possible contributors to these density differences including temperature differences, salinity contrasts, suspended material, both organic and inorganic, as well as combinations of these mechanisms. In water resources management, to prevent sedimentation and drain sediment from dam reservoirs, in most cases, is tried to remove fine sediments by using of hydrodynamic forces. One of these methods is removing sediment by the gravity current dynamics. The gravity currents are the most important effective events on reservoir's sedimentary processes (transmission, distribution and deposition of particles). Therefore, identifying the factors affecting on this type of flow is very important. Both bed slope and surface roughness are the most important parameters on gravity current dynamics and have a considerable effect on velocity distribution of the body.
This study investigates the effect of bed roughness on the gravity currents characteristics and obtaining relations between effective parameter on gravity currents and finding their influences in body velocity distribution on rough and slope beds. In this study to determine the bed roughness effect on the velocity distribution, a series of gravity currents experiments was performed on beds with the rough natural element sizes of 4, 8, 12 and 15 mm and slopes from 0 to 2.2%. The experiments of saline gravity currents were carried out in three inflow discharges with amounts of 0.7, 1 and 1.3 (lit/s) and with three concentrations of 10, 16 and 20 (gr/lit). All Experiments was done on a Plexiglas flume with the length of 10 meters and width of 35 cm and in hydraulic models laboratory in Shahid Chamran university of Ahvaz (Iran). In all experiment, the velocity profile in the body of gravity current was acquired by an Ultrasonic Doppler Velocity meter device.Keywords: Gravity Current, Jet, Wall Region, Velocity Distribution., Bed Roughness -
Page 11By calculating the potential yield of crops, the yield gap can be determined. Also, the management of water and soil resources can be done based on effective factors in production. Attending to limitations in agricultural programs, for increasing the wheat yield, applying a predictive model of yield potential rate could be effective in the yield management based on various parameters including climate factors. A suitable model for this goal is regionalism method based on cropped ecology, in which it is assumed the climatic requirements are prepared and water, nutrients, salinity, pests and diseases do not effect on growth and yield potential. In this study, with the goal of estimating the yield potential in 30 province centers in Iran, by using factor and cluster analysis and considering nine climate factors, the stations were classified. Then, climate data was collected for a period of 30 years (up to 2010) from synoptic stations of province's centers. Also, the average values of actual yield for winter wheat were collected from the site of Centre for Information and Communication Technology, Ministry of Agriculture Jahad, in each province. Then the potential yield and gap yield (kilograms per hectare) were determined, using the regionalism method based on cropped ecology.
According to the variance percentage and special values of effective factors, KMO value was equal to 0.598; forwhich, four climatic parameters including mean temperature, relative humidity, wind speed and precipitation were selected in classifying process. Therefore, nine provinces were categorized in dry climate including Birjand, Kerman, Yazd, Zahedan, Qom, Isfahan, Shiraz, Tehran and Semnan. In semi-arid and semi-humid climate there were 15 provinces including Mashhad, Qazvin, Bojnoord, Sanandaj, Kermanshah, Tabriz, Arak, Shahrekord, Zanjan, Hamedan, Orumieh, Ilam, Khorramabad, Yasooj and Ardebil. Three provinces were in humid with low precipitation climate including Busheher, Bandarabbas and Ahvaz. Also, three provinces were in humid-rainy climate including Rasht, Sari and Gorgan. Then, the estimated values of parameters in regionalism method based on cropped ecology and the actual and potential of wheat yield were determined in four regional groups. Also, the percentage of yield gap which caused by the difference between the actual and potential yield due to the lack of climatic requirements was determined in each group.
As a result, regionalism method based on cropped ecology in arid -low precipitation regions of the country could not estimate potential yield of wheat. In the region with semi-arid and semi-humid climate, the lowest and highest amount of yield gap was 30.9 and 53.5 percent in Ardebil and Ilam stations, respectively. Also, the average percentage of yield gap in this group was 43.2 percent. In regions with arid climate, the minimum and maximum of yield gap percentage were in Semnan and Zahedan stations with 10.2 and 52.1 percent, respectivelyand the average percentage of yield gap in this group was 27.5 percent. In humid-rainy region, including the stations of the northern coasts of Iran, the percentage of yield gap was very high. In this group, the maximum yield gap was in the Rasht station with 68.9 percent and the minimum was in Gorgan station with 41.2 percent. Also, the average percentage of yield gap in this group was 55.8 percent. Among the four regional groups, in terms of required climatic parameters, the maximum of yield gap occurred in humid-rainy regions and the lowest was in low precipitation humid climate. Furthermore, among 30 provinces which were studied, the Rasht station from humid-rainy regions had the highest yield gap and the Semnan station from dry regions had the lowest. The results of factor analysis showed that temperature was important for classifying the provinces by climate. Also, the results showed that model had the greatest sensitivity to the average temperature, comparing to the other climatic factors. This caused the lowest yield gap occured in the regions with low precipitation and dry climates with 15 and 27.5 percent, respectively. Also, the moderate climate in rainy climate condition caused more growing and a longer season of growth in the most species. Therefore, the difference between average yield and potential yield has been occcured which caused the average of yield gap in the northern coastal provinces of the country to be 55.8 percent. The conclusion of results represented that the regionalism method based on cropped ecology is one of the appropriate tools for predicting wheat yield and determining yield gap in different climatic regions.Keywords: Meteorological data, Climate changes, Factor, cluster analysis, Wheat yield efficiency -
Page 31On the one hand, the limited energy resources and on the other hand increasing energy demand are the factors that enhance the importance of energy management in micro and macro programs, in developed and developing countries. The head loss in different parts of pipelines causes an increase in the amount of energy consumed by the pumps. Therefore, reducing it is critical. One of the most important and easiest ways to reduce energy consumption during fluid supply is using friction-reducing materials in the flow.
In this study, a setup was conducted to find out the effect of adding a small dose of three different polymers (Karaakryl30, Karaakryl58 and Sodium carboxymethyl cellulose) on changing friction losses and its effect on the spray radius of sprinklers in a pressurized irrigation system; in which three Hunter type sprinklers were used. To investigate the effects of these three polymers, the irrigation system consisted of a galvanized tank, a valve to control the flow rate at the beginning of the pipeline, polyethylene pipes with a diameter of 20 mm and a length of 40 meters, three Hunter sprinklers. And two amount of pressure were used. In the first stage, pure water was used to measure the flow rate and the head loss. In the first stage, the experiments were related to the control treatments in which pure water was used to measure the flow rate and energy loss. First, the tank was filled with water. Then, the initial and final sprinkler's output were closed. The sprinkler's nozzle at the end of the path was opened to measure the discharges. And then the control valve fully opened. By switching the pump, water from the sprinkler's output flowed out. After this stage, the discharge valve fully opened and the head loss was measured. Then, by adjusting the control valve (reducing the flow), the amount of flow and pressure were measured again. Adjusting the control valve was performed in 5 different conditions with pure water. The polymer Karaakryl30 was tested with the concentration of 100 mg per liter which was once dissolved in water and then stirred after half an hour again. After that, the pump started immediately. In the same way of the pure water measuring, the pressure and flow rate were measured. Then the Pump was turned off and all three sprinklers were installed, adjusted and fasten to fulcrums in the given situations. After pump working, the radius of spraying was measured. In the next stage, without any changes to the system compared to its previous state, the tank was filled with water again and the pump was turned on. After pumping all the water in the tank, the spraying radius was measured and compared with a spraying radius of Karaakryl30 solution (100 mg per liter). All the above steps for each of the three materials were tested at three concentrations. By fitting a curve for discharge and head loss data, the relationship between them was obtained, for pure water and the three polymers in three concentrations. Using these relationships, the head loss values were obtained in the same discharges. Considering that friction-reducing material reduces head loss in the pipeline, it increases the pressure in sprinkler nozzles and as a result increases the spraying radius. In this aspect, the effect of adding the polymers on spraying radius was studied, in three concentrations.
The results showed that Karaakryl30 and Sodium carboxymethyl cellulose had the better effect on reducing the head loss and increasing the spraying radius, comparing to the Karaakryl58. The amount of head loss in the discharge of 0.423 lit/s for pure water, Karaakryl30, Karaakryl58 with the concentration of 1000 ppm and Sodium carboxymethyl cellulose with the concentration of 150 ppm was 20.04, 12.8, 16.09 and 12.39 meter, respectively, at the pipe length of 40m and diameter of 20mm. Karaakryl30 with 1000 ppm concentration led to maximum reduction of friction (36.1%) in pipelines and maximum increase in sprinkler spray radius (11.6%), comparing to pure water. Using polymer Karaakryl58 at the same level of concentration increased the same parameters 19.6 and 5.9 percent, respectively. Sodium carboxymethyl cellulose (150 ppm) increased the same measured parameters to 38.1 and 15.6 percent, respectively.
Thus, Karaakryl30 and Sodium carboxymethyl cellulose could be recommended to use in sprinkler systems in short-time and emergency situations such as high evapotranspiration for which it is not economical to consider the pipes with larger diameters in design.Keywords: Friction loss, Sprinkler spray radiation., Polymer, Pipeline, Drag Reducing Agent -
Page 39Hydraulic jumps in sloping channels were first classified by Kindsvater (1944). An A-jump is the jump that starts at the foot of the chute. For the B-jump, the toe of the jump forms on a positive slope and the roller ending on the downstream stilling basin. The C-jump begins on a positive slope and the roller length ending at the foot of the chute and finally the D-jump that its roller length entirely occurs on the chute. The most common type of jump in practice is the B-jump which from the point of hydraulic calculation, is more complicated than others (Hager, 1988). For a jump on sloping beds, application of the one-dimensional (1-D) momentum equation is not easy to compute the sequent depth ratio because some additional information is needed to estimate both the weight component of the jump and the bottom pressure acting on the sloping channel portion (Carollo et al, 2011). Hence, for B-jump sequent depth ratio, previous researchers developed several empirical equations to compute the sequent flow depth of the jump. Although, many researches have been done about jump on smooth sloping beds, so far there has not been a study on roughened sloping beds. The primary aim of the present study was to evaluate a general equation to calculate the sequent depth ratio of different types of jump on the roughened sloping bed. For this purpose, different types of jump on both the smooth and the roughened sloping bed were investigated experimentally in a wide range of relative roughness (0.032 0.54), different slopes of chute (14.5, 20.5 and 27.5 degree) and Froude number (1.3 7.4). Using -Buckingham theory of dimensional analysis, the effective parameters for determination of sequent depth ratio on smooth and rough bed were obtained and by applying of the incomplete self - similarity theory, a general equation was derived. Then, the coefficients of the general equation were derived and its accuracy was determined, using experimental data from present study and previous published resources.
One of the main reasons of using roughened beds is to dissipate some of the kinetic energy of the flowing water over the chute. In order to examine the effect of the roughness size, the Froude number at the start of the jump and the sequent depth were calculated. In this study, The variation between the sequent depth ratio and upstream Froude number for different types of jumps and for relative roughness sizes (ks/yc) was figured. A parameter E is introduced by Hager (1944) and defines as the location where the jump begins. The results showed for a given range of E, by increase of the relative roughness, the upstream Froude number and sequent depth ratio decreased. In order to develop a functional relationship for computing the B-jump sequent depth ratio on sloping smooth beds, the function m should be determined from experimental data. Therefore, the procedure of Carollo et al (2011) was followed and the values of m(α,E) were calculated for different slopes. In this function, for a given chute angle α, the function m decreases as E increases, and for a fixed E value, m decreases as the angle α increases.
The form of relation between the function m and the variable E is exponential. The empirical relationship ( ) was fitted to each experimental series corresponding to a given α value. In which a and b are coefficients that have to be experimentally determined. For smooth beds, the coefficients a and b depend only on the chute slope (α). In order to illustrate the effect of roughness height on the function m (or the sequent depth ratio), for each slope, the m values were plotted versus E values for various roughness heights. For example, the results of chute with bed angle of 14.5 illustrated that the coefficients a and b depended not only on the slope of the chute, but also on the roughness height. Also, the measurements, for different roughness sizes, indicated that m decreased as E increased and, for a fixed E value, m slightly decreased when the roughness size increased. A comparison between the rough and hydraulically smooth bed conditions shows that boundary roughness reduced m and the sequent depth ratio. This result confirmed the findings of previous investigations. The coefficients a and b was obtained from the investigations of Carollo et al (2011) in the analysis of the data in this study. Using experimental data, two equations obtained for a and b coefficients. Thus, a general equation for sequent depth ratio (Y) obtained. A comparison between experimental Y values and those calculated by main Equation was illustrasted in this study. The lines representing the 15% error range are also illustrated.
The results showed that the general equation has high accuracy and the average percentage error is about 7.5%. Also, the result showed that increasing relative roughness decreased the sequent depth ratio and the Froude number of the upstream flow and when relative roughness is about 0.54 the upstream Froude number and sequent depth ratio decreased about 50%.Keywords: Energy dissipation, Hydraulic jump, Self-similarity., Chute, Roughness, Stilling basin -
Page 49Being aware of such phenomena like dam break may seriously prevent human and environmental disasters. Dam's breakage may lead to devastating waves in downstream valleys, which cause casualties, severe damages to infrastructure and human loss of life. Therefore, leaves huge economic burdens to the governments. Studies related to dam break, especially the flood wave, are necessary for safety management in the downstream of dams. Wave velocity and bed roughness are the two main parameters affecting the qualitative and quantitative changes of flood spreading in downstream. Bed roughness can produce more resistance against flood flow and thus, reduces the wave velocity and increases the depth of flow. It can change the form of wave profile, as well.
The flood wave caused by dam break phenomenon can be simulated experimentally and numerically. In the laboratory, a dam break can be produced by an abrupt opening of sluice gate, which is installed in a flume. In this research study, the dam break has been numerically modeled by using the three- dimensional Flow-3D model. The measured water elevations and velocities of flood wave were provided from another research study, which was carried out in the laboratory of Water Sciences Engineering Faculty, Shahid Chamran University of Ahvaz. Flow-3D model is powerful software in CFD, which has been developed by Flow Science Inc. and is generally used for open channels flow modeling. In this software, five turbulence models are applied and each of them is able to produce a calculation mechanism to quantify the effect of turbulence fluctuations on flow parameters such as the average velocity and the depth. In this study, the geometry model was produced in AutoCAD software and the appropriate grids in upstream and downstream of the sluice gate were developed by the Flow-3D model. The positions of the small triangular pieces used to generate rough bed in downstream were defined as the fixed points in the model. Defining the boundary conditions, the wall conditions were used for the bed and symmetry condition assumed for walls. The initial conditions were defined as constant water elevation in upstream and downstream of the sluice gate.
In this laboratory study, 4 values for downstream water depths (0, 3, 6, 9cm), 3 values for upstream head (24, 28 and 32cm) and 3 types of bed roughness (soft, regular and zigzag patterns of small triangular pieces) were the variables and the water elevation and velocity profiles were provided using the image processing method and a suitable software. In this software, the user was able to specify the wave profile at any time of wave propagation and any place. In this numerical simulation, the upstream head was just 24 cm, so for all 12 water elevations, the profiles were simulated and were compared with the corresponding measured profiles. All five turbulence models were used, to complete the numerical simulation. Comparison of the simulated and measured water elevation profiles showed that the Prandtl and one equation turbulence models were not able to simulate the flow conditions, accurately. The RNG turbulence model showed the best performance among all the proposed turbulence models. Therefore, after verifying the RNG turbulence model, it was possible to print out the water elevation and velocity profiles at any time and for any type of bed roughness and downstream water depths.
Analyzing the simulated results showed that the bed roughness was an important parameter in progressive wave velocity and acted as a resistant factor. It was found that, for the zigzag pattern of bed roughness, the longitudinal decreasing slope of the wave velocity was higher than the regular pattern. And this could be due to the form of ordering the small triangular pieces. The turbulence in the flow, near the bed, was greater for the zigzag pattern in comparison to the regular pattern, due to the interface of pieces and water which caused more resistance, lower velocities and higher depths in longitudinal profile, in the same hydraulic boundaries. Downstream depth was another parameter which affected the form of longitudinal velocity and the water elevation profiles. The average wave velocity decreased with increasing downstream water depth. It was obvious that increasing the downstream water depth amplify the effect of bed roughness in decreasing the velocity. The maximum effectiveness was simulated for the zigzag pattern with 9cm downstream water depth. The other effect of downstream water depth was in the form of the wave break, immediately after the dam break.
According to the results of model predictions and measured values, it was found that the effect of downstream water elevation on reducing wave velocity was much more than the bed roughness. For example, the average of reduction in velocity, soft bed to rough bed, was calculated %24. However, this value for 0 and 9cm downstream depth was %71.Keywords: Numerical modeling, Dam break, River engineering, Flood -
Page 57The flow duration curve (FDC) is one of the most important and widely used signatures of catchment runoff response. It has been used in numerous hydrological applications as a part of planning activities related to agriculture, industry, urban water supply, navigation, flood control, reservoir and sedimentation studies, in-stream flow assessment, water quality management, and impacts of land use changes. The FDC is a graph of river discharge plotted against exceedance frequency and is normally derived from the complete time series of recorded river flows. Representation of the relationship between the frequency and magnitude of stream flows, making it a compact signature of a catchments functioning; it can be used to diagnose the rainfall-runoff responses in gauged catchments at a holistic functional level, as well as to regionalize them to ungauged catchments. For these reasons, in the past few decades, considerable effort has been expended towards detailed studies of FDCs, especially in the context of predictions in ungauged basins (Smakhtin, 1998, 2001, Castellarina, 2004, Niadas, 2005, Santhi, 2008. Li, 2010). However, most recent studies on the FDCs have been empirically based, which generally fall into two categories: graphical and statistical. The graphical approach focuses on exploring the controls of catchment climatic and physiographic characteristics on the shape of the FDC (Smakhtin et al), while the statistical approach employs statistical distributions to fit the FDC and then relates the parameters of the distribution to the catchments physical characteristics (Castellarin et al., 2004; Li et al., 2010). Several studies have been conducted on the impact of geometric and hydrologic parameters on the FDC. (Booker and Snelder, 2012. Brown, 2013). In unregulated streams, the shape of FDC is controlled chiefly by the geology of the basin and research on the effects of lithology on FDC in a certain area, leading to obtaining information for use in water resources management (Bloomfield, 2009).
Although both graphical and statistical methods demonstrated that different climatic and landscape characteristics impact the shape of the FDCs in different regions of the world (Li et al., 2010), it has been difficult to generalize the results from one region to another.
The aim of this research was to determine the factors affecting the flow duration curve with emphasis on the role of lithology in the estimation of flow duration indices.
The study area comprises 18 basins with areas from 69.1 to 3941 km2 which are located in the North part of Iran and is underlain mostly by rocks belonging to the hard formations. Carbonate formations in this area are appropriate distribution and cover about 28% of the studied area. Mean annual rainfall varies from 427 to 1131 mm/ year. Median slopes vary from 2.85 to 15.9%.
In this research, using data from the daily flow of rivers, flow duration curve indices, in eighteen hydrometric stations of the Caspian Basin were determined. Physiographic parameters extracted from the GIS environment and the various factors of climate, geology and land use were also calculated. Lithological units using digital geological map, with a scale of 1: 250,000 and based on expert opinions were divided into 4 classes and the area covered by each unit was calculated. Land use map based on satellite images with the division of forest, pasture, agriculture, surface water resources, and rock outcrops were prepared. Factor analysis and cluster analysis were conducted and the regression equations were derived. The accuracy of models using the coefficient of determination and standard errors were investigated. Finally, the most appropriate models in 1% and 5%, significant of levels extracted.
Factor analysis of 21 variables showed, six factors were chosen as the most important independent factors. In the multivariate regressions between agents of geologic formations, particularly carbonate hard formation and alluvial and between morphometric parameters, slope, and main river length and from land use classes, agricultural land use was determined as the best predictor of duration curve indices. Also, according to the power and exponential models, annual precipitation as the most important factors between climatologically factors was diagnosed. The validation results showed that the distributions of errors in all models were normalized and Durbin-Watson index was varied between 1.5 to 2.5.According to regional relations obtained, the contribution of hard carbonate formations along with Quaternary deposits, indicating their role in ensuring the continuity of the region's rivers.Keywords: Flowduration Curve indices, Carbonates, Regional Model, Lithology, Alluvium -
Page 69Groundwater is one of the main sources to supply drinking water, especially in arid and semi-arid areas. The interaction of surface water and groundwater is one of the important issues in groundwater management and modeling, and it is essential to understand. This research aims at investigating river and groundwater interactions in the Dezful-Andimesh plain, which is about 2073 Km2 as the largest agricultural plain in the northern part of Khuzestan province, southwest of Iran. This plain contains five smaller plains: western Dez, eastern Dez, Sabili, Deymche and Lor. The Dez, Karkhe, Balaroud, Kohnak and Shavoor rivers are located in the study area. Thus this area is very important for supplying water needs and requires an efficient water resource management, a correct planning and reconsideration of water policy. Only 4% of the region is classified as the urban area and the remaining are agricultural plains. All the plains have irrigation networks except the Lor plain. Although, the irrigation networks supply a great part of the crops demand, in some regions groundwater use is also significant. Wheat, sugarcane and maize (32%, 21% and 16%) are the most common crops in this region. The cultivation period of wheat is usually from November to May, so the maximum recharge of groundwater by irrigation networks occurs in this period. The study area has a semi-arid climate with the mean monthly temperature of 36.5°C in July and 11.8°C in January. The precipitation varies from 296 mm in eastern Dez to 394 mm in Lor. The Dezful aquifer with an average thickness about 100 m is an unconfined aquifer system. There are over 2700 wells which pump about 500 million m3 per year, in this area. The general direction of groundwater flow in the area is from the north to the south. Recharging the aquifer is via direct rainfall infiltration and return water from irrigation networks. Hydraulic conductivity range is from 14 m/day for clayey sediments in Sabili to 49 m/day for sandy deposits, particularly in Lor. In order to model the groundwater, a MODFLOW programming code was implemented in a GMS environment. MODFLOW is a numerical model that is one of the best methods for assessing the quantity and quality of groundwater. The numerical models are difficult and time-consuming. However, in the recent decade's research using simulation models have been developed due to the improvement of high-speed computers. Information related to steady-state flow along river information was used in the model. Water level and bed river data were measured in river-cross sections and hydrometric stations in different months. The model was configured to a steady-state condition in October (2009) and was calibrated using 62 groundwater wells. Also, the model was calibrated in an unsteady-state condition for five years (October 2006 to September 2011). Then, the model was validated for the period of 2011-2012.
According to the results, the regulated rivers such as Karkhe and Dez were the major source of recharging aquifers; however, there was a decrease in their contribution to aquifers in wet seasons. There was an increase in groundwater level during wet seasons due to the groundwater recharge from rainfall and irrigation systems. However, there were no significant changes in the water level of the regulated rivers, and the minimum water level was observed in November. Therefore, an increase in groundwater level could not be related to rivers when the river levels were constant, during the wet season. The Balaroud non-regulated rivers generally drained the aquifers during dry seasons, but not during wet seasons (February), when recharge occurred from these rivers to aquifers. The rivers such as Shavoor and Kohnak and their tributaries drained aquifers in every month. There were not increasing or decreasing patterns in groundwater discharge when there was a change in both rivers and groundwater levels with different patterns in different months of a year.
According to the simulation results over six years, the most river- groundwater interaction was belong to Dez and Karkhe rivers with the highest amount of recharge to aquifers. The Shavoor and Kohnak rivers and their tributaries generally drained groundwater. However, there was no specific pattern of river-groundwater interaction for the Balaroud seasonal rivers. There was more groundwater recharge through regulated rivers in dry seasons compared to wet seasons which was due to the high amount of rainfall (along with the amount of water from irrigation systems) in wet seasons.Keywords: Dez River, Hydraulically connected, Groundwater model, GMS, Regulation of river flows. -
Page 79Water is undoubtedly one of the factors of productivity and reducing water consuming reduces the yield. But the decrease in yield due to water stress is inevitable. Because of limited water resources in Iran, deficit irrigation and water stress on crops are inevitable. Also in some regions, saline water is being used for irrigation. The saline lands of Iran are about 44.5 million hectares with various salinity and alkalinity, which are 27% of the area of the country and they are mainly located in the central plateau, southern coasts of Iran and the Caspian Sea coasts. Saline water has been successfully used for irrigating the crops such as cotton, wheat, barley and sugar beet. Deficit Irrigation (DI) is a method in irrigation management and an engineering technique which is used to optimize the yield in water scarcity and increasing water efficiency. In this method, some irrigation events which have less effect on productivity could be omitted or the volume of water in each irrigation event could be reduced. The salinity of water resources and soils in Iran is one of the problems barricades the effective application of these resources in agriculture. However, it has been tried to use unconventional water resources such as Caspian Sea to ease water shortage, especially in droughts in Northern provinces. Various reports have been presented on the effect of salinity on cotton. Some research expressed that irrigating by 6,000 mg of salt per liter had acceptable yield without any bad effect on the quality of cotton fiber and seeds. In this study, the simultaneous effects of periodic water stress and salinity stress due to using the Caspian Sea water on cotton yield was investigated, which was carried out in the agricultural research farm of Gorgan University of Agricultural Sciences and Natural Resources. This research is based on a factorial experiment in completely block design (CRD) including 20 treatments and 3 replicates in the farm of Gorgan University of Agriculture on spring 2012. Three cotton seed of Golestan variety were planted in lysimeters with the diameter of 35 cm and an average height of 70 cm. The distance between Lysimeters was one meter. To measure the cotton water requirement, pan (class A) was used. Also the effective rain during the cultivation was measured with a rain gauge which was installed near the research farm. Periodic water stress in four stages (Initial, early stage of development, the second stage of development, midseason) was applied and the salinity stress was performed at 5 levels (0% 25%, 50%, 75% and 100% mixture of Caspian sea water with ordinary water) .The Salinity was 0.7, 4.1, 11.1, 18.11, 25.1 (dS/m), respectively.
Results of this study showed that the effect of Periodic water stress and salinity from Caspian Sea water on Plant's height, boll weight, seed cotton yield was significant at the 1% level. But its effect on boll number was significant at 5% level. Also, the interaction of water stress and salinity on boll weight, seed cotton yield and plant height was significant at 1% level but the interaction of water stress and salinity at the boll number was not significant. Comparison of mean of the salinity stress in soil was done at each level of water stresses in order to reach a better understanding of the interactions. Depths of the measurements were 0-10 cm and 10-20 cm from the soil surface. By increasing salinity of irrigating water and applying water stress, soil salinity was increased. The highest soil salinity among the different level was occurred in S5 level. Because the soil salinity was increased in all treatments which were irrigated by saline water, considering leaching water in water events is necessary.
Periodic water stress also increased the salinity in soil. This increase was higher in the midseason period due to the length of this period.
By the results, the cotton yield was generally decreased by increasing salinity of irrigating water. The interaction of periodic water stress and salinity had significant effects on the cotton yield. But at the level of 25% mix of Caspian Sea water with ordinary water, in the early stage of crop development, the difference of results were lower than other treatments. So this planning could be an appropriate choice for Northern provinces of Iran in droughts and limitation of water resources.Keywords: Caspian Sea, Salt water, Water stress., Cotton -
Page 87Hydraulic jump is one of the most important phenomena in rapid varied flow which occurs in free surface flow by transforming from supercritical into sub critical flow and causes the depth to be increased in a small distance of flow. Changes in the geometrical parameters of the bed including roughness, the slope of the floor, as well as shape change and their effects on the characteristics of hydraulic jump have always drawn the attention of the researchers.
In the present study, the characteristics of the hydraulic jump formed in contractions and expansions with straight and curved walls were investigated for five Froude numbers of 5.8-9.1. In order to create supercritical flow in the contractions, in all experiments, the opening value of the sliding gate was 2 cm and the width of the upper and lower canals were 80 and 40 cm, respectively, (the convergence ratio equals 0.5). Also, in the experiments concerned with expansions, the opening of the sliding gate was 3 cm and the width of the upper and lower canals equaled to 40 and 80 cm (the divergence ration equals 0.5). An audio flume meter model UFM610P was used to measure the discharge flow of the flume. The accuracy of the flow meter was ±0.02 liters per second and its sensors were installed on the hydration pipe of the flume. Also, a bathometer with the accuracy of ±0.1 was utilized to measure the height of the waves surface. In addition, a horizon surface, two-dimensional electromagnetic speedometer with the accuracy of ±0.5 was employed to measure the wave velocity of the different parts of the wave.
The velocity profiles showed that, in all sections, the velocity increased as it got away from the floor and after reaching to its maximum, its value decreased. In fact, the occurrence of hydraulic jump causes intense blending of air and water and formation of eddy currents near the water surface. As a result, the energy of the flow decreases and its velocity value reduces, too. Besides, it was observed that by getting distance from the origin of the jump, the boundary layer extends while the maximum values reduce. On the other hand, negative velocity values were observed in the velocity profiles of contractions near the water surface. This phenomenon is a sign of the overcome of the return flows to the main flow of the canal (near the water surface). Indeed, the entrance of the hydraulic jump to the tightening section caused an increase in the intensity of water and air blending and occurrence of the return flow phenomenon. This situation was observed more with the increase of the Froude number. Furthermore, the velocity value in contraction was higher than its counterpart in expansion for a fixed Froude number in the vertical direction. Also, the results showed that the secondary depth and jump length in contractions were more than expansions for the fixed Froude number while the formed hydraulic jump in contractions had less energy dissipation than expansions. In fact, widening of section intensifies flow energy loss in expansions. However, in contractions, the gradual decrease of the width of the section prevents complete flow energy dissipation. For this reason, energy dissipation or flow from supercritical to subcritical takes longer time and consequently the length of hydraulic jump and the depth of secondary jump increase. On the other hand, the results indicated that the amount of tension faced a sudden drop and reached to a negative value after getting to its maximum value. Considering that Reynolds tensions have a direct relationship with fluctuations of point velocity, the sudden drop of the tension value can be due to the intense blending of air and water in jump surface and after that the reduction of range of fluctuations of point velocity. Therefore, a higher negative value was observed in contractions compared to expansions because of more intense water and air blending. It needs to be noted that the contraction with straight wall had higher negative value. Besides, in both types of the conversions, the highest negative value of the tension was observed in the middle section. Also, in contractions, by getting distance from the origin of the jump, the maximum values of Reynolds tension increase. In fact, the progress of the jump in tightening section increases the fluctuations of point velocity and through this increases the intensity of the tension. However, in expansions, the opposite case happened. On the other hand, by comparing the vertical distributions of velocity and tension in both conversions, it can be found that the site of the occurrence of maximum velocity in different sections of the jump relatively matches the site of the occurrence of the maximum intensity. The values of energy dissipation in the expansions with straight and curved walls were averagely obtained 24.95% and 21.23%, relatively. For the contractions, they were averagely obtained 15.93% and 12.49%. Also, the result showed Reynolds tensions increased in the contraction and decreased in the expansion by increasing the longitudinal distance from the origin of the jump.Keywords: Contraction, Dissipation energy, Hydraulic jump, Reynolds stresses, Expansion -
Page 97Different human activities depend on water availability so most of the civilizations have been emerged and developed in the proximities of rivers. However, such developments lead to polluting water bodies (Alizade, 2004). For sustainable agricultural and industrial developments monitoring the quality of water resources is inevitable. Heavy metals negatively impact water quality as they form complex and poisonous compounds when entering water bodies (Ostan, 2005). Heavy metals exchange between sediments and water and the equilibrium exists between these phases in rivers. Such equilibrium changes by any variation in media reaction, oxidation and reduction condition, the amount of organic matter, microorganism's activities, concentration of soluble elements and cation exchange capacity of suspended particles. The sources of heavy metals in water resources are anthropogenic (from human activities) or naturally released from mineral weathering (Khosheqbal et al. 2012). Hoba et al. (1983) reported that severe erosion of stones in the Vender river watershed and releasing a lot of Ca2 led to formation of carbonates of calcium and some heavy metals and also Fe and Mn hydroxides in sediments of the river. Zhao et al. (2012) analyzed heavy metal concentration of fish and crab's tissues and found a significant and positive correlation between the concentration in legs of crabs, skins and heads of fish and the concentrations in water. Sediments are considered pollutants as they absorb and transport heavy metals. Akhondi et al. (1989) studied heavy metal concentration of the Qomroud River and reported that the concentrations of Cd, Zn and Pb in water are below the permissible threshold for drinking. They expressed that sources of heavy metals are petroleum refineries, metal mills, electrical and textile industries. The Zayanderoud River is one of the main rivers of interior plateau of Iran also a water resource for drinking, industries and croplands in Chaharmahal & Bakhtiari and Esfahan provinces. The river origins from the Zagros Mountains and passes through cultivated lands, orchards, fish farms and proximity of steel mill industries and pours into Gavkhoni wetland in the southeast of Esfahan City. Due to the importance of the river, the objectives of this study were monitoring spatial and temporal variations of heavy metal concentrations in water and sediments of the Zayanderoud River.
Mean annual rainfall of the Zayanderoud watershed is 450mm and the surface area of the watershed is 27100 km2 (Johari, 2010). To fulfill the objective 10 points were selected between the Zayanderoud Dam in the west and Zarrinshahr in the east based on overlook land uses and the accessibility of the river for sediment and water samplings. Water and sediment samples were collected bimonthly according to the Geological Survey of US (1978) in 2012. Sediment samples were collected from the surficial (0-10 cm) of river bed, air dried, grinded and passed through 2 mm sieve. Total heavy metals Pb, Zn, Cu and Cd of sediments were extracted by acid (Sposito et al., 1982) and the extracted heavy metals and the concentration of water samples were measured using Atomic Absorption Spectrometer. Landsat ETM data of 2012 were used and the land use map of the study area was delineated.
According to the land use map, the overlooked areas to the river were rangelands, stone outcrops, settlements, rice cultivated areas, orchards and industrial zones. By the temporal and spatial variations of mean Pb concentration in water and sediments, the highest mean Pb concentration belonged to point 10 which was significantly different from of all points except point 9. The results also indicated that the maximum mean Pb concentration of sediments was 20.63 ppm and belongs to point 9. Wang et al. (2012) studied the water quality of the Lancang River and expressed that the exhausted smoke of cars was the source of Pb. The main sources of Pb in the study area could be steel mill industries as point 9 locates in the vicinity of these industries. A significant and positive correlation existed between Pb concentration of water and of sediments in all points. Also, the Pb concentration of water and sediments increased in rainy months which showed possible transportation of Pb by runoff from polluted areas to the river. Significant intercorrelation coefficients among Pb, Zn and Cu indicated possible common sources of these elements which are different from Cd resources. Spatial and temporal variations of Zn mean concentration in water were significant. The maximum mean concentration of Zn in water was 6.2 ppb which was below the permitted threshold for drinking according to the Institute of Standards and Industrial Researches of Iran (2010). Maximum concentration of Zn in sediments was 34.36 ppm. Zn used as a plant micronutrient and a component of insecticides and herbicides in croplands and orchards in the area. On the other hand, maximum Zn mean concentration observed in point 9 (5.6 ppb) and in sediments of point 2. The land use map showed that the overlooked land use of point 2 were orchards while of point 9 were steel mill industries. Therefore, the source of Zn could be fertilizers and insecticides which applied in orchards. Likewise, Wang et al. (2012) reported that Zn application as fertilizer and insecticide in a part of China led to river pollution. Namazi et al. (2014) analyzed the dust in the proximities of steel mill industries and reported that they contain heavy metals consist of Zn. Mean Cd concentration in water increased from the west to the east of studied points. The highest and lowest Cd mean concentration sediments were in points 10th and 1st, respectively. Maximum mean Cu concentration was observed in the sediments of the 6th point.
Regarding the permissible thresholds of studied heavy metals for drinking water (Institute of Standards and Industrial Researches of Iran, 1388), the concentrations of the studied heavy metals were below the threshold and the water in all points was safe for drinking. The results also indicated that industrial areas, croplands and orchards were the main non-point resources of heavy metal for the river and runoff was one of the main ways of transferring pollutants to the river.Keywords: Heavy metal, Zayanderoud River., Sediment, Water quality, Pollution -
Page 107Infiltration is one the most sensitive hydraulic parameters that has more effect on the efficiency of water, especially in surface irrigation design and management. The fact that infiltration is affected by several factors of soil and irrigation method makes it a complex parameter. Evaluating the impacts of irrigation water head on cumulative infiltration and using this parameter in infiltration models will lead to better performance of the models which are widely used in different conditions. So far, many studies have been conducted, in the field of infiltration. Forman et al (2006) investigated the effect of water height on parameters of Kostiakov soil infiltration equation in the soils with Sandy loam, loam and Silt showed that the parameters k and a had the most variability in sandy loam soil texture. However, the parameter k, in the Silty soil texture, showed little variation due to water height changes. Valyntzas et al (2009) studied the effect of water height as an upstream boundary condition on the parameters of Kostiakov equation in Sandy Clay Loam soil texture and showed that by modifying the infiltration coefficient based on water height, the simulation results of infiltration by HYDRUS improved. The purpose of this study was to investigate the effect of different water heights on the infiltration process, in some soil texture in Foumanat plain in Guilan province. In this way, 17 soil samples that had more than 15% Clay and less than 30% Sand were gathered. Infiltration tests were carried out with three replications in a completely randomized design, in a physical model. The factors included soil texture with 17 levels: the height of water with three levels (3, 5 and 7 centimeters) and three levels of the time. The results also showed that height of water on the soil surface had more impact on Loam soil texture than Clay soil texture. Infiltration testing time was 600 minutes and cumulative infiltration values were calculated in 30, 360 and 600 minutes. The physical model includes a water tank, manometer for reading the water level, reservoir water for fixing water level on the surface of soil and polyethylene cylinder with an internal diameter of 15.23 cm and a height of 35 cm, as a soil box. To maintain the soil water level, a float was used. Analysis of variance hypotheses for normal distribution was done. After passing normal distribution test, Analysis of variance and the comparison of averages were done by Tukeys method at 5% level using SAS software. Results showed that treatment of water height had a significant effect on cumulative infiltration of the four soil samples in Silty Clay and Loam Silty Clay soil textures. Although, the results of statistical tests showed no significant effect between water height treatments (3, 5 and 7 cm) in 13 soil samples, the difference between the amount of infiltrated water in the samples should be considered due to the critical situation of water resources and necessity of increasing in irrigation efficiency. Comparison of the difference percentage of cumulative infiltration in different soil textures showed that, apart from the 3 cm water height treatment, in the two other water height treatments, Clay Loam, Silty Clay Loam, Silty loam and Silty Clay soil textures had the highest percentage of difference between the two time periods, respectively. On the others in estimating of cumulative infiltration, Clay Loam texture has more sensitivity to the cumulative time of infiltration. The difference of cumulative infiltration between treatments 3 and 5 cm was more than the other two treatments (5 and 7 cm), despite the equal difference of 2 cm in the water height. The reason for this result could be explained by reducing the influence of gravity due to increasing the height of water on the soil surface. Investigating the three infiltration time treatments (30, 360 and 600 minutes) showed that in 14 soil samples, there were significant differences between cumulative infiltration in 30 min with infiltration time in 360 and 600 minutes (pKeywords: Cumulative infiltration, Irrigation water depth, Clay soil, Infiltration duration, Surface irrigation
-
Page 117Comprehensive management of water resources nowadays demands accurate and complete meteorological and hydrological data. Simulation of runoff hydrograph for hydraulic structures design often requires rainfall characteristics such as the amount, duration, and temporal rainfall distribution. The design storm is one of the main factors in flood estimation in watersheds. Considering the temporal variations in heavy storms, precise recognition of the rainfall temporal pattern is very important. Temporal rainfall distribution is one of the significant features of rainfall affecting the quantity and quality of runoff and erosion, which has unfortunately attracted less attention. Temporal rainfall distribution means the amount of rainfall volume per temporal periods of precipitation. In fact, the temporal rainfall distribution pattern indicates the changes in rainfall intensity. The effect of rainfall characteristics, especially temporal rainfall distribution and its intensity, on simulation of runoff events is more understandable. A huge part of errors in the simulation of rainfall-runoff models is due to the uncertainties regarding regions rainfall distribution patterns. Therefore, provision of rainfall distribution pattern for each region by using rainfall data is very important since the pattern varies in regions with different climates. In the studies in Iran, with regard to the variety of mechanisms causing rainfall events (altitude, proximity to the sea, and general circulation systems), it is impossible to consider a predominant distribution pattern for all parts of Iran. In this regard, presented 6 and 24-hour rainfall patterns of American Soil Conservation Society have been used. The SCS patterns which are broadly used in Iran and other countries in order to simulate runoff and designate hydraulic structures.
In this research temporal rainfall distribution has been evaluated in the Sanganeh research watershed. The Sanganeh watershed is located in the north-east of Khorasan Razavi province (6015' E, 3641' N), with the area of 50 ha. According to the collected data from the climatologic station and Domartan coefficient of 10.2, the climate of the watershed is classified as semi-arid. The area receives less than 180 mm annual precipitation. In this study, for determination of temporal rainfall distribution Pilgrim method has been used. In this method, in every time duration based on this fact that in which quarter the maximum rain has happened, the first, second, third and fourth rainstorms were defined and temporal rainfall distribution pattern was evaluated by use of the data obtained from rain gauge station. In this study, by investigating all observed storms in the Sanganeh research station in the period of 2005-2013, all individual storms have been derived in the different durations of 9 hr. Then, the Pilgrim method was used based on storms separation in four quarters. By averaging ranks in each quartile, the index rank the final pattern for each time-duration was specified. To statistically evaluate the obtained patterns with durations of interest, Chi square test was used after establishing a contingency table and estimating the rainfall percentage for each quartile and ranking. By averaging ranks in each quartile, index rank of the final pattern of each time-duration was determined.
The obtained results showed that most of the storms which occurred in Sanganeh region were in the first quadrate. So the assumption that the highest amount of precipitation occurs in the middle of time duration is doubtable. Also, with increasing of the storms duration, the number of the storms decreased in the first quadrate and increased in the second to fourth quadrates. Moreover, 61% of erosive storms had the duration of fewer than 3 hours and most of them were in second and third quadrate. Chi square test was used for statistical evaluation of the obtained patterns with distinct durations. The results showed that in the duration of 9 hours, there was a significant relationship between the rank and the quadrate. These patterns had high ability to use while in the others durations, the patterns were not reliable. In addition, comparison of the obtained logical pattern with SCS type patterns showed that only rainfall temporal pattern of the second quadrate rainfalls had a high correlation coefficient with IA type. Thus, as the most of the storms in the Sanganeh watershed were in the first quadrate and had the duration of 1 to 3 hours, using SCS type pattern in this region should be precautionary.Keywords: Quadrate, Precipitation Pattern, Pilgrim Methods, SCS, Chi Square Test, Sanganeh Watershed -
Page 127The groundwater contamination by hydrocarbon products has been a common problem in industrial areas in recent decades. Leakage from storage tanks or during the transportation seeps through the soil and pollutes the groundwater resources. Materials such as benzene, toluene, ethyl-benzene and xylene (BTEX) are common compounds of gasoline and can be found frequently in subsurface layers. In order to understand the behavior and to analyze the distribution of this kind of contamination in subsurface layers, a laboratory experimental were carried out as part of this research. An experimental sand tank was constructed in University of Tehran Laboratory to simulate the movement of BTEX in unsaturated and saturated porous media. The dimensions of the tank were 70, 50, and 8 cm for length, height, and thickness, respectively. The walls of the tank were made of Plexiglas to make the movement of the plume visible and possible for photography during the experiment. Since BTEX is colorless, SUDAN III, a red soluble substance, was added to BTEX. Image processing technique was used to estimate the concentration of BTEX in porous media in different time steps. In this experiment, 200 milliliters of BTEX with a flow rate of 15 milliliters per minute injected into the soil layer 2 cm below the surface. A groundwater flow in the saturated zone was applied with the rate of 10.4 cm per day. The upgradient depth of water was 12 cm and that for the downgradient was kept at 11.8 cm.
The profiles showed that the maximum degree of saturation happened 50 seconds after the start of the experiment with 0.92. The thickness of BTEX plume in the unsaturated zone reached 4 cm after 650 seconds. The degree of saturation and BTEX velocity were decreased with depth and increase in water saturation. After stopping the BTEX injection, the degree of saturation decreased in upper layer until it reached the residual saturation. The degree of saturation was increased above the water table and since it could not penetrate the saturated zone, it spread laterally.
The image processing technique using MATLAB software was used to estimate the BTEX degree of saturation at different location of the plume and for different time steps. The concept behind image processing algorithm is that the intensity of the pixel in a photograph is directly related to the degree of saturation of the plume. High degree of saturation means that the presence of the BTEX is higher in that location and the red color has higher intensity. The correlation between the pixel color intensity and BTEX saturation degree was obtained by calibrating the pixel intensity verses the measured saturation. The calibration process showed that there is a linear relationship between saturation degree and pixel intensity. The saturations extracted by image processing technique were used to produce the iso-lines of BTEX saturation for different time steps using Surfer 10.0 software. These plumes were compared to photographs for each time step to additionally confirm the accuracy of the technique.
The results of the image processing technique and comparison with the laboratory results showed that the proposed technique can accurately estimate the plume distribution. In addition, the results showed that in the initial stages of the BTEX injection, the gravitational force is the dominate force and lateral movement is negligible compare to vertical movement. Once the plume reaches the capillary fringe, BTEX may migrate laterally as a continuous, free-phase along the upper boundary of the water saturated zone due to gravity and capillary forces. BTEX saturation and accumulated volume increase in the capillary fringe and they depend on the depth of the capillary zone, BTEX viscosity, hydraulic head of the injection source, and volume of the injected BTEX. BTEX depressed the capillary fringe by approximately 4 cm. Ultimately, BTEX plume stopped above the water table.
The importance of groundwater, as the source of agricultural, industrial, and drinking water, makes it imperative to protect this precious resource. Understanding the behavior of the petroleum hydrocarbons in saturated and unsaturated zone can help to better prevent the pollution or remediate the contaminated groundwater. The successful application of the image processing technique in BTEX saturation estimate proved that this method can reduce the cost of laboratory experiment and can be used in evaluation of the numerical models, as wellKeywords: Unsaturation, BTEX, Image analysis, Porous media -
Page 135The world's growing population has caused irrecoverable variations on natural resources. Groundwater as one of these resources has an important role in the health and economy of countries. Nowadays, groundwater is the first source of drinking water in arid and semi-arid regions such as Iran. In order to manage the consumption of groundwater resources, knowledge of groundwater recharge in spatial and temporal scales before detailed exploration such as pumping test and geoelectrical studies is required. Hence, preparing and analyzing the factors affecting groundwater recharge in the GIS (Geographic information system) can be very helpful. The existence of a groundwater recharge potential (GRP) map is vital in regions where groundwater is the most important source of water supply. The aim of this study is to provide a GRP map for Kore river basin in the Fars province considering the interaction between the factors affecting groundwater recharge. Quantitative estimate of recharge due to precipitation using remote sensing (RS) and GIS is another aim of this study. Kore river basin is located in the south of Iran and north of Fars province, between geographic coordinates of 560000 to 740000 and 3250000 to 3450000 (zone 39) in the Mercator scale (UTM). This basin has 5 sub-basins including Khosro-Shirin, Bakan, Dozkord-Kamfirooz, Asopas and Marvdasht-kherameh with an area of 8284 km2. The annual rainfall of the study area varies from 300 mm in the southeast (plains) to more than 800 mm in the northwest (mountainous region). Geologically, this area includes two different carbonated and non-carbonated parts. Carbonated formations with fractures or karst systems cover 40 percent of the study area. These formations are Asmari, Sarvak, Ilam, Tarbur, Daryan in the mountainous regions while non-carbonated formations cover other parts of the study area (lowland regions).
Satellite data (Landsat images) of the study area were interpreted in the Erdas Imaging, Envi and Arc GIS software. In this study, thematic maps of 7 factors including lineament, and longitudinal density drainage, lithology, karst terrains, land use-land cover, land slope and rainfall for Kor River basin were prepared. Then, each of the factors was classified in 9 classes and their maps were designed. These classes are in the numerical range for lineaments, drainage, lithology, slope and rainfall while they are in the descriptive manner for karst domain and land use. The interactions between the factors were considered using a schematic pattern and the internal rate of each factor was extracted. Then this rate was multiplied by the specific weight of each factor. In this way, the effect of each factor on GRP was obtained in terms of percentage which is known as refereeing method. Then, by assigning an internal weight and considering the interaction of these parameters with each other in regard to GRP a weighted rating was calculated for every parameter in GIS software. By applying these weighted rates and using the final weight of the parameters a classified map of GRP in five classes from very high to very low grades was prepared. Taking into account more factors and their interaction is the advantage of this study compared to previous studies.
Based on the obtained recharge rates using RS and GIS, the areas of recharge zones and annual rainfall, the volume of recharged water was estimated. This estimate was equal to 1090.59 MCM which is about 20 percent of annual rainfall in the Kore river basin. In order to verify this estimated recharge, hydroclimatological balance was applied to recalculate the recharge volume. Considering the lack of hydrometery data in the Kore river basin, the SCS method was used to calculate the runoff volume regarding different land use. By these methods, the percentage of recharge from rainfall was estimated to be 18 percent. Accordingly, recharge volume was equal to 954.13 MCM. Finally, by comparing recharge estimation in this article with that of hydroclimatological balance method, it was revealed that there is about 14% difference in the estimated values of each of these methods. It should be mentioned that in this study, the effect of unsaturated zone properties on the GRP has not been considered. This methodology can be applied in similar basins in arid and semiarid regions in order to estimate the percentage of infiltration resulted from monthly excess water in the hydroclimatological balance method.Keywords: Kor River basin, Hydroclimatological balance, GIS, Groundwater recharge potential -
Page 147This study is presented to investigate the effectiveness of simultaneous use of cable and collar in reduction of local scour depth in groups of two and three piers. The experiments were carried out in a 20m long, 0.60m wide, and 0.60m deep, Plexiglas-sided flume. The working section, in which piers were located, was 8m long with a recess on the bed 0.2m deep and was situated 8m downstream from the entrance of the flume in order to generate fully developed flows. The recess was filled with uniform sediment with the mean particle sizes of 0.71mm and geometric standard deviation of particles was σg=1.135. Cylindrical Perspex pipes with diameter of 60mm were used as pier models. The tests were conducted on a group of three piers in line. For each group, the spacing (S/b where S is the distance between the piers sections and b is the pier diameter) of the piers was varied. The collar with the diameter of 3 times wider than the pier diameter was developed at the bed level around group of three piers. The cable with the diameter equal to 0.15 of the pier diameter with a torsion angle equal to 15° and combination of them around each pier was studied.
In experiment of group piers with collar, two different spacing of S/b=2 and 3 were tested for a group of three piers in line. The experiment with spacing S/b=3 the space between two piers were not completely covered with collar. The scour hole was developed at back of the third pier. Also, two grooves were gradually developed at the down-stream rim of the collar and gradually developed and extended upstream and eventually reached the upstream edge. At this moment, the flow was intensified through the groove, reducing the side slope of the groove and with sediment removal from the grooves, the scour hole extended upstream of the pier and below the collar. By penetration of the flow beneath the front edge of the collar the rate of sediment removal increased. Maximum depth of scour hole is equal in front of the first pier with and without collar. In the first pier, both spaces had the same delay time and maximum scour depth because the space between piers had no effect on efficiency of collar. The results of experiment with cable showed that reduction of scour depth in the second and third pier is more than the first pier, due to sheltering effect of the first pier on the second and third piers. Also, in this case, with increasing the space between the piers the reduction of scour depth was decreased slightly. The simultaneous use of cable and collar showed that reduction of scour depth in upstream of group of three piers with S/b=2 is about 20%, 45% and 53% respectively. Since the strength of wake vortex reduces, reduction of scour depth in upstream of the pier protected with cable and collar is more than the reduction of scour depth of the pier protected with collar and pier protected with cable. In this case, comparison of delay time reveals that using collar caused the delay in scouring mechanism and when the collar and cable combined the delay time is increased. In the first pier, cable and collar cannot cause salient reduction in maximum scour depth due to the reinforcement effects counteract the effect of collar and cable. In the second and the third piers when the collar and cable combined, the delay time and equilibrium time are increased. In these piers, the maximum scour depth is decreased as well. The efficiency of cable and collar to reduce the scour depth in the third pier is greater than the second pier due to increasing the sheltering effect.
The results indicate that the scouring is postponed and reduced in protected pier with cable and collar more than protected pier with collar; therefore cable increases the efficiency of collar. The space between the piers has a considerable effect in reducing the scour depth, so that the scour depth was decreased and the delay time was increased, with increasing the space between the piers. The most reduction of scour depth was observed for group of three piers with S/b=2, in which the scour depth of the first, second and third piers reduced to about 20%, 45% and 53% respectively. In situation that a pier protected with cable and collar with S/b=3, the scouring postponed was maximum. So, use the cable with the collar can greatly reduce the risk of bridge destruction.Keywords: Control of local scour, Collar, Cable, Group piers, Local scour depth -
Page 157Several formulas have been proposed for the calculation of the Instantaneous Unit Hydrograph (IUH). The concept of series reservoirs is one of the methods to calculate and predict the flow Hydrograph in the catchment. Inverse gamma model is one of probability distributions, which depends on two parameters: n and k. And several formulas have existed for the estimation of the Inverse gammas instantaneous unit hydrograph. The accuracy of estimating the parameters n and k in Inverse gamma model directly affect simulation results. This study for the evaluating parameters of the Inverse gamma model used different methods such as moment, optimization and Rosso. The results showed that using optimization method in estimating parameters n, k had the most compatibility with the observed hydrographs. In all these methods, output runoff volume has been calculated with the error less than 11%. The optimization, Rosso and moment methods were successful in computing the time peak and the peak discharge. But in 5 rainfall event, the moment method, in computing the peak discharge has more than 50% of the difference with the observed value. Thus regionalization of Inverse gammas model parameters (n and k), by the mentioned methods, is advised for development and application of this hydrological model, in similar catchments.Several formulas have been proposed for the calculation of the Instantaneous Unit Hydrograph (IUH). The concept of series reservoirs is one of the methods to calculate and predict the flow Hydrograph in the catchment. Inverse gamma model is one of probability distributions, which depends on two parameters: n and k. And several formulas have existed for the estimation of the Inverse gammas instantaneous unit hydrograph. The accuracy of estimating the parameters n and k in Inverse gamma model directly affect simulation results. This study for the evaluating parameters of the Inverse gamma model used different methods such as moment, optimization and Rosso. The results showed that using optimization method in estimating parameters n, k had the most compatibility with the observed hydrographs. In all these methods, output runoff volume has been calculated with the error less than 11%. The optimization, Rosso and moment methods were successful in computing the time peak and the peak discharge. But in 5 rainfall event, the moment method, in computing the peak discharge has more than 50% of the difference with the observed value. Thus regionalization of Inverse gammas model parameters (n and k), by the mentioned methods, is advised for development and application of this hydrological model, in similar catchmentsSeveral formulas have been proposed for the calculation of the Instantaneous Unit Hydrograph (IUH). The concept of series reservoirs is one of the methods to calculate and predict the flow Hydrograph in the catchment. Inverse gamma model is one of probability distributions, which depends on two parameters: n and k. And several formulas have existed for the estimation of the Inverse gammas instantaneous unit hydrograph. The accuracy of estimating the parameters n and k in Inverse gamma model directly affect simulation results. This study for the evaluating parameters of the Inverse gamma model used different methods such as moment, optimization and Rosso. The results showed that using optimization method in estimating parameters n, k had the most compatibility with the observed hydrographs. In all these methods, output runoff volume has been calculated with the error less than 11%. The optimization, Rosso and moment methods were successful in computing the time peak and the peak discharge. But in 5 rainfall event, the moment method, in computing the peak discharge has more than 50% of the difference with the observed value. Thus regionalization of Inverse gammas model parameters (n and k), by the mentioned methods, is advised for development and application of this hydrological model, in similar catchments..Several formulas have been proposed for the calculation of the Instantaneous Unit Hydrograph (IUH). The concept of series reservoirs is one of the methods to calculate and predict the flow Hydrograph in the catchment. Inverse gamma model is one of probability distributions, which depends on two parameters: n and k. And several formulas have existed for the estimation of the Inverse gammas instantaneous unit hydrograph. The accuracy of estimating the parameters n and k in Inverse gamma model directly affect simulation results. This study for the evaluating parameters of the Inverse gamma model used different methods such as moment, optimization and Rosso. The results showed that using optimization method in estimating parameters n, k had the most compatibility with the observed hydrographs. In all these methods, output runoff volume has been calculated with the error less than 11%. The optimization, Rosso and moment methods were successful in computing the time peak and the peak discharge. But in 5 rainfall event, the moment method, in computing the peak discharge has more than 50% of the difference with the observed value. Thus regionalization of Inverse gammas model parameters (n and k), by the mentioned methods, is advised for development and application of this hydrological model, in similar catchments.Keywords: Inverse gamma model, Optimization method, Rosso model., Moment method, Jonog watershed, N, K parameters estimation
-
Page 163The purpose of this study was to evaluate the Snowmelt Runoff Model (SRM) to estimate of daily runoff at Galinak hydrometric station in Taleghan watershed, Iran. In this study, the input variables of the model such as hydrology and meteorology variables and physical condition of the catchment area were calculated and a total of 72 pass from satellite images of MODIS was used to determine the snow cover area. Comparison of measured and calculated hydrographs was carried out using the efficiency coefficient of Nash-Sutcliffe (R2) and difference of volume percent (Dv). The R2 and Dv values were found equal to 0.95 and 2.07 percent in the calibration period and equal to 0.86 and 4.81 percent in the validation period for Galinak hydrometric station. The results show that the model proved to be an excellent alternative, performing high-quality daily snowmelt runoff prediction.
The purpose of this study was to evaluate the Snowmelt Runoff Model (SRM) to estimate of daily runoff at Galinak hydrometric station in Taleghan watershed, Iran. In this study, the input variables of the model such as hydrology and meteorology variables and physical condition of the catchment area were calculated and a total of 72 pass from satellite images of MODIS was used to determine the snow cover area. Comparison of measured and calculated hydrographs was carried out using the efficiency coefficient of Nash-Sutcliffe (R2) and difference of volume percent (Dv). The R2 and Dv values were found equal to 0.95 and 2.07 percent in the calibration period and equal to 0.86 and 4.81 percent in the validation period for Galinak hydrometric station. The results show that the model proved to be an excellent alternative, performing high-quality daily snowmelt runoff prediction.
The purpose of this study was to evaluate the Snowmelt Runoff Model (SRM) to estimate of daily runoff at Galinak hydrometric station in Taleghan watershed, Iran. In this study, the input variables of the model such as hydrology and meteorology variables and physical condition of the catchment area were calculated and a total of 72 pass from satellite images of MODIS was used to determine the snow cover area. Comparison of measured and calculated hydrographs was carried out using the efficiency coefficient of Nash-Sutcliffe (R2) and difference of volume percent (Dv). The R2 and Dv values were found equal to 0.95 and 2.07 percent in the calibration period and equal to 0.86 and 4.81 percent in the validation period for Galinak hydrometric station. The results show that the model proved to be an excellent alternative, performing high-quality daily snowmelt runoff prediction.
The purpose of this study was to evaluate the Snowmelt Runoff Model (SRM) to estimate of daily runoff at Galinak hydrometric station in Taleghan watershed, Iran. In this study, the input variables of the model such as hydrology and meteorology variables and physical condition of the catchment area were calculated and a total of 72 pass from satellite images of MODIS was used to determine the snow cover area. Comparison of measured and calculated hydrographs was carried out using the efficiency coefficient of Nash-Sutcliffe (R2) and difference of volume percent (Dv). The R2 and Dv values were found equal to 0.95 and 2.07 percent in the calibration period and equal to 0.86 and 4.81 percent in the validation period for Galinak hydrometric station. The results show that the model proved to be an excellent alternative, performing high-quality daily snowmelt runoff prediction.
The purpose of this study was to evaluate the Snowmelt Runoff Model (SRM) to estimate of daily runoff at Galinak hydrometric station in Taleghan watershed, Iran. In this study, the input variables of the model such as hydrology and meteorology variables and physical condition of the catchment area were calculated and a total of 72 pass from satellite images of MODIS was used to determine the snow cover area. Comparison of measured and calculated hydrographs was carried out using the efficiency coefficient of Nash-Sutcliffe (R2) and difference of volume percent (Dv). The R2 and Dv values were found equal to 0.95 and 2.07 percent in the calibration period and equal to 0.86 and 4.81 percent in the validation period for Galinak hydrometric station. The results show that the model proved to be an excellent alternative, performing high-quality daily snowmelt runoff prediction.Keywords: Snowmelt, SRM model, Snow cover area, Runoff, MODIS image -
Page 169In this research, the biophysical potentials of agricultural crops were analyzed in Qazvin irrigation and drainage network using DSSAT crop models. These models were calibrated and validated using field data or proper input coefficients were selected for the crops without such data. Land Units (LUs) were extracted by overlaying soil mapping units (SMUs) and weather grids using GIS. Weather grids were mapped using Thiessen method assuming each unit as a homogeneous area. Crop growth simulation was carried out using the DSSAT crop models. For each crop, two irrigation regimes (surface and pressurized) and 2-3 planting dates were considered. The results revealed the difference of crop production for each planting date and irrigation regime. Among the major crops in the study area, the values of spatial coefficient of variation were mediocre for Wheat and Barley, but lower for Corn. In relation to Tomato, Soybean and Cotton, significant differences were found between LUs and this issue must be taken in consideration when determining cropping pattern. The index of time coefficient of variation can also be used as a measure of risk in relation to the studied crops. The results of this research can be used as inputs of a spatial planning support system (PSS).In this research, the biophysical potentials of agricultural crops were analyzed in Qazvin irrigation and drainage network using DSSAT crop models. These models were calibrated and validated using field data or proper input coefficients were selected for the crops without such data. Land Units (LUs) were extracted by overlaying soil mapping units (SMUs) and weather grids using GIS. Weather grids were mapped using Thiessen method assuming each unit as a homogeneous area. Crop growth simulation was carried out using the DSSAT crop models. For each crop, two irrigation regimes (surface and pressurized) and 2-3 planting dates were considered. The results revealed the difference of crop production for each planting date and irrigation regime. Among the major crops in the study area, the values of spatial coefficient of variation were mediocre for Wheat and Barley, but lower for Corn. In relation to Tomato, Soybean and Cotton, significant differences were found between LUs and this issue must be taken in consideration when determining cropping pattern. The index of time coefficient of variation can also be used as a measure of risk in relation to the studied crops. The results of this research can be used as inputs of a spatial planning support system (PSS).In this research, the biophysical potentials of agricultural crops were analyzed in Qazvin irrigation and drainage network using DSSAT crop models. These models were calibrated and validated using field data or proper input coefficients were selected for the crops without such data. Land Units (LUs) were extracted by overlaying soil mapping units (SMUs) and weather grids using GIS. Weather grids were mapped using Thiessen method assuming each unit as a homogeneous area. Crop growth simulation was carried out using the DSSAT crop models. For each crop, two irrigation regimes (surface and pressurized) and 2-3 planting dates were considered. The results revealed the difference of crop production for each planting date and irrigation regime. Among the major crops in the study area, the values of spatial coefficient of variation were mediocre for Wheat and Barley, but lower for Corn. In relation to Tomato, Soybean and Cotton, significant differences were found between LUs and this issue must be taken in consideration when determining cropping pattern. The index of time coefficient of variation can also be used as a measure of risk in relation to the studied crops. The results of this research can be used as inputs of a spatial planning support system (PSS).In this research, the biophysical potentials of agricultural crops were analyzed in Qazvin irrigation and drainage network using DSSAT crop models. These models were calibrated and validated using field data or proper input coefficients were selected for the crops without such data. Land Units (LUs) were extracted by overlaying soil mapping units (SMUs) and weather grids using GIS. Weather grids were mapped using Thiessen method assuming each unit as a homogeneous area. Crop growth simulation was carried out using the DSSAT crop models. For each crop, two irrigation regimes (surface and pressurized) and 2-3 planting dates were considered. The results revealed the difference of crop production for each planting date and irrigation regime. Among the major crops in the study area, the values of spatial coefficient of variation were mediocre for Wheat and Barley, but lower for Corn. In relation to Tomato, Soybean and Cotton, significant differences were found between LUs and this issue must be taken in consideration when determining cropping pattern. The index of time coefficient of variation can also be used as a measure of risk in relation to the studied crops. The results of this research can be used as inputs of a spatial planning support system (PSS).Keywords: risk, Iran, Planning, Land units, Simulation
-
Page 175Increasing the efficiency of the agricultural water pumping stations is important because of the high percent of cost and energy consumption. The performance of pumping stations in the evaluation of irrigation and drainage networks has been ignored up to now, whereas any impairment in pumping station will cause the destruction on priming and thereupon disruption of network performance. Since, there is not specific assessment for pump's performance, therefore, in this study, a quantitative evaluation method named data envelopment analysis, was introduced and applied. In this regard, after the introduction of this technique and its various models, the performance of Sabili main station located in Khuzestan province was calculated during 2001 to 2010 with DEAP2 software. In this method, the average performance of Sabili main station was calculated 69 percent. This station had the highest efficiency in 2001, 2002 and 2010 (100 percent) and it had the lowest efficiency in 2009 (21 percent). Most of the inefficiency was due to the high cost of electricity consumption. Hence, maximum potential efficiency could be accessible with energy consumption reduction activities.
Increasing the efficiency of the agricultural water pumping stations is important because of the high percent of cost and energy consumption. The performance of pumping stations in the evaluation of irrigation and drainage networks has been ignored up to now, whereas any impairment in pumping station will cause the destruction on priming and thereupon disruption of network performance. Since, there is not specific assessment for pump's performance, therefore, in this study, a quantitative evaluation method named data envelopment analysis, was introduced and applied. In this regard, after the introduction of this technique and its various models, the performance of Sabili main station located in Khuzestan province was calculated during 2001 to 2010 with DEAP2 software. In this method, the average performance of Sabili main station was calculated 69 percent. This station had the highest efficiency in 2001, 2002 and 2010 (100 percent) and it had the lowest efficiency in 2009 (21 percent). Most of the inefficiency was due to the high cost of electricity consumption. Hence, maximum potential efficiency could be accessible with energy consumption reduction activities.Increasing the efficiency of the agricultural water pumping stations is important because of the high percent of cost and energy consumption. The performance of pumping stations in the evaluation of irrigation and drainage networks has been ignored up to now, whereas any impairment in pumping station will cause the destruction on priming and thereupon disruption of network performance. Since, there is not specific assessment for pump's performance, therefore, in this study, a quantitative evaluation method named data envelopment analysis, was introduced and applied. In this regard, after the introduction of this technique and its various models, the performance of Sabili main station located in Khuzestan province was calculated during 2001 to 2010 with DEAP2 software. In this method, the average performance of Sabili main station was calculated 69 percent. This station had the highest efficiency in 2001, 2002 and 2010 (100 percent) and it had the lowest efficiency in 2009 (21 percent). Most of the inefficiency was due to the high cost of electricity consumption. Hence, maximum potential efficiency could be accessible with energy consumption reduction activities.
Increasing the efficiency of the agricultural water pumping stations is important because of the high percent of cost and energy consumption. The performance of pumping stations in the evaluation of irrigation and drainage networks has been ignored up to now, whereas any impairment in pumping station will cause the destruction on priming and thereupon disruption of network performance. Since, there is not specific assessment for pump's performance, therefore, in this study, a quantitative evaluation method named data envelopment analysis, was introduced and applied. In this regard, after the introduction of this technique and its various models, the performance of Sabili main station located in Khuzestan province was calculated during 2001 to 2010 with DEAP2 software. In this method, the average performance of Sabili main station was calculated 69 percent. This station had the highest efficiency in 2001, 2002 and 2010 (100 percent) and it had the lowest efficiency in 2009 (21 percent). Most of the inefficiency was due to the high cost of electricity consumption. Hence, maximum potential efficiency could be accessible with energy consumption reduction activities.
Increasing the efficiency of the agricultural water pumping stations is important because of the high percent of cost and energy consumption. The performance of pumping stations in the evaluation of irrigation and drainage networks has been ignored up to now, whereas any impairment in pumping station will cause the destruction on priming and thereupon disruption of network performance. Since, there is not specific assessment for pump's performance, therefore, in this study, a quantitative evaluation method named data envelopment analysis, was introduced and applied. In this regard, after the introduction of this technique and its various models, the performance of Sabili main station located in Khuzestan province was calculated during 2001 to 2010 with DEAP2 software. In this method, the average performance of Sabili main station was calculated 69 percent. This station had the highest efficiency in 2001, 2002 and 2010 (100 percent) and it had the lowest efficiency in 2009 (21 percent). Most of the inefficiency was due to the high cost of electricity consumption. Hence, maximum potential efficiency could be accessible with energy consumption reduction activities.Keywords: Efficiency, Performance evaluation, Sabili main pumping station., Data Envelopment Analysis method, Irrigation, drainage network
©2000-2024 magiran.com All Rights Reserved.