تحلیل نا پایداری بستر کانال با استفاده از تنش برشی بحرانی و دبی واحد بحرانی

الموضوعات :

میلاد رستمی ¹ ( دانشجو )
محمد مهدی حسین زاده ² ( دانشیار دانشگاه شهید بهشتی )

تاريخ الإرسال : 14 الثلاثاء , ذو الحجة, 1438 تاريخ التأكيد : 25 الأربعاء , ربيع الأول, 1439 تاريخ الإصدار : 01 الأربعاء , ذو الحجة, 1438

الکلمات المفتاحية: آستانه حرکت ذرات رسوبی, تنش برشی بحرانی, دبی واحد بحرانی, رود خانه گلالی,

ملخص المقالة :

از مسائل مهم و کلیدی در بررسی فرسایش و پایداری رودخانه‌ها مسئله آغاز حرکت ذرات رسوبی می‌باشد. جریانی که سبب شروع حرکت ذره خواهد شد، جریان بحرانی است. آستانه حرکت اولیه ذرات رسوبی با چند روش مختلف بررسی می‌شود که مهمترین آنها میانگین سرعت بحرانی جریان، روش تنش برشی بحرانی و روش دبی واحد جریان. برای بررسی پایداری رسوبات بستر رودخانه گلالی قروه با توجه به شرایط رودخانه از روش تنش برشی بحرانی و دبی واحد بحرانی استفاده کردیم. با هدف بررسی آستانه حرکت ذرات در بستر این رودخانه شرایط دو مقطع عرضی از این رودخانه بررسی شد. در ادامه تنش برشی بحرانی و دبی واحد بحرانی برای چهار اندازه ذره D16، D50، D84 و D95 در هر دو مقطع محاسبه شد و مشخص گردید که در مقطع عرضی شماره 1و2 در دبی بانکفول مقطع تمامی ذرات با اندازه‌های D16، D50، D84 و D95 شروع به حرکت می‌کنند و بنابراین در دبی بانکفول شرایط رسوبی بستر رودخانه تغییر کلی خواهد کرد. در مقطع شماره 3 ذرات با قطر D95 و بیشتر از آن در بستر رودخانه پایدار خواهند ماند و سایر ذرات کوچکتر از این اندازه شروع به حرکت می‌کنند. وجود شیب بیشتر بستر و عرض کمتر کانال فعال در مقطع عرضی شماره 1و2 سبب شده که تمامی ذرات مورد بررسی ناپایدار باشند. در نهایت نتایج به دست آمده از هر دو روش همخوانی بالایی داشتند و این خود تأییدی بر استفاده از این روش‌ها برای مطالعات پایدار سازی و تولید رسوب رودخانه‌ها می‌باشد.

نص كامل:



تحلیل ناپایداری بستر کانال با استفاده از تنش برشی بحرانی و دبی واحد بحرانی

( مطالعه موردی رودخانه گلالی قروه)

چکیده:

از مسائل مهم و کلیدی در بررسی فرسایش و پایداری رودخانه‌ها مسئله آغاز حرکت ذرات رسوبی می‌باشد. تولید رسوب رودخانه  سبب کدورت آب، مشکلات مواد مغذی و آلودگی‌های آب، مدفون نمودن تأسیسات انحراف آب و همچنین پر شدن دریاچه‌های ذخیره آب می‌گردد. جریانی که سبب شروع حرکت ذره خواهد شد، جریان بحرانی نامیده می‌شود. آستانه حرکت اولیه ذرات رسوبی با چند روش مختلف بررسی می‌شود که مهمترین آنها عبارتند از میانگین سرعت بحرانی جریان، روش تنش برشی بحرانی و روش دبی واحد جریان. برای بررسی پایداری رسوبات بستر رودخانه گلالی قروه با توجه به شرایط رودخانه از روش تنش برشی بحرانی و دبی واحد بحرانی استفاده شد. با هدف بررسی آستانه حرکت ذرات در بستر این رودخانه شرایط سه مقطع عرضی از این رودخانه بررسی شد. در ادامه تنش برشی بحرانی و دبی واحد بحرانی برای چهار اندازه ذره D16، D50، D84 و D95 در هر دو مقطع محاسبه شد و مشخص گردید که در مقطع عرضی شماره 1و2 در دبی بانکفول مقطع تمامی ذرات با اندازه‌های D16، D50، D84 و D95 شروع به حرکت می‌کنند و بنابراین در دبی بانکفول شرایط رسوبی بستر رودخانه تغییر کلی خواهد کرد. در مقطع شماره 3 ذرات با قطر D95 و بیشتر از آن در بستر رودخانه پایدار خواهند ماند و سایر ذرات کوچکتر از این اندازه شروع به حرکت می‌کنند. وجود شیب بیشتر بستر و عرض کمتر کانال فعال در مقطع عرضی شماره 1و2 سبب شده که تمامی ذرات مورد بررسی ناپایدار باشند. در نهایت نتایج به دست آمده از هر دو روش همخوانی بالایی داشتند و این خود تأییدی بر استفاده از این روش‌ها برای مطالعات پایدار سازی و تولید رسوب رودخانه‌ها می‌باشد.

کلید واژه:  آستانه حرکت ذرات رسوبی، تنش برشی بحرانی، دبی واحد بحرانی، رود خانه گلالی

 

 

 

 

 

 

مقدمه

 

رودخانه فراهم کننده آب و انرژی برای طبیعت و انسان است و می توان گفت تآمین آب مهمترین نقش اقتصادی رودخانه است. ناپایداری مجرا و سواحل رودخانه نه تنها موجب خسارات به اراضی مستعد کشاورزی و تآسیسات مجاور رودخانه شده بلکه رسوبات حاصل از فرسایش در برخی از موارد درصد قابل توجهی از مجموع کل رسوبات انتقالی توسط جریان رودخانه را شامل می شوند[12].کاهش ظرفیت انتقال رودخانه ناشی از عوامل مختلفی است که هر یک به گونه ای با تآثیر گذاری در رفتار جریان و مشخصه های هندسی تشدید حالت سیلابی را سبب می گردد. شکل ظاهری و رفتار رودخانه متآثر از توان جریان (دینامیک جریان) و مقاومت رسوبات در برابر فرسایش است .[9] درک درست از انتقال رسوب و حساسیت کانال به فرسایش یا رسوب در مدیریت فرایندهای رودخانه ای ضروری است[10]. ناپایداری مجرا و سواحل رودخانه نه تنها موجب خسارات به اراضی مستعد کشاورزی و تآسیسات مجاور رودخانه شده بلکه رسوبات حاصل از فرسایش در برخی از موارد درصد قابل توجهی از مجموع کل رسوبات انتقالی توسط جریان رودخانه را شامل می شوند[5]. نا پایداری مجرا و سواحل رودخانه ها در توسعه پهنه سیلابی و مدیریت منابع آب اهمیت دارد[12]. مواد بستر و کناره رودخانه نقش مؤثری در تغییرات مورفولوژی رودخانه ایفا می کند. رودخانه های با مواد بستری ریزدانه نسبت به تغییرات دبی و بار رسوبی، حساسیت بیشتری از خود نشان می دهند[4]. در حالت کلی معمولآ رودخانه ای که بار رسوبات آن رس و سیلت باشد، نسبت به آبراهه ای که دارای بار رسوبی ماسه و ریگ است، فعال و ناپایدار تر بوده و کانال اصلی معمولآ باریک تر و عمیق تر است. با توجه به این موارد بافت رسوبات موجود در یک رودخانه نقش به سزایی در بررسی میزان فرسایش و میزان پایداری یک رودخانه می تواند داشته باشد[14]. درک حرکت رسوب و حساسیت کانال نسبت به فرسایش یا رسوبگذاری، ضرورتی مهم برای مدیریت رودخانه است. برای مثال غلبه فرآیند رسوبگذاری در یک بازه از رودخانه می‌تواند انتقال رسوب را کاهش دهد و درنتیجه افزایش قدرت رود و غلبه فرسایش در بازه پایین دست شود، و این فرسایش سبب ناپایداری و خسارت به زیرساخت‌ها می‌گردد. طرح‌های مدیریتی اخیر جهت ارزیابی فرآیندهای ژئومورفولوژیکی با هدف طراحی مناسب و پیشنهاد راه‌حل‌های مؤثر ارائه می‌شود[13,15]. راه حل های متعددی جهت تحلیل ناپایداری کانال رودخانه وجود دارد از جمله آنها روش تنش برشی بحرانی و دبی واحد بحرانی می باشد که در این پژوهش از این دو مدل جهت برآورد میزان ناپایداری در رود خانه گلالی قروه استفاده شده است. تاکنون تحقیقات زیادی توسط محققین داخلی و خارجی در زمینه ناپایداری و فرسایش کانال انجام شده است از جمله: کوریت (2014) با روش پیامدهای منبع غیر نقطه­ای رسوب(BANCS) برای پیش‌بینی فرسایش و پایداری کرانه رودخانه در استونی کلوو کریک¹ در کت اسکیلز² پرداخته است.

---

[1]  - Stony Clove Creek

[2]  - Catskills

این مطالعه نشان داد که نتایج بدست آمده از طریق روش BANCS ممکن است منجر به بهبود مدیریت آینده حوضه استونی کلوو شود. بیزی و لرنر(2015) به بررسی فرآیند فرسایش و رسوب با استفاده از شاخص قدرت رود به عنوان شاخص حساسیت کانال پرداخته اند. که در پایان منطقه مورد مطالعه را بر اساس تنش برشی بحرانی و دبی واحد بحرانی به چهار کلاس فرسایشی طبقه بندی کرده اند. قمشی و وجدانی(1384) به بررسی تنش برشی بحرانی فرسایش رسوبات چسبنده و نقش آن در کانال های رو باز پرداخته اند که باتوجه به نتیجه به دست آمده برای بازه مورد مطالعه این انتظار می رود که دبی طراحی کانال، رسوباتی راکه در بستر کانال ته نشین می شود را فرسایش می دهد. حداد چی و همکاران (1390) در پژوهشی معادلات تجربی برآورد دبی بار بستر در رودخانه کوهستانی با بستر شنی را در رودخانه چهل چای در استان گلستان مورد ارزیابی قرار داده اند. نتایج نشان می دهد که معادله ایکرز و وایت بهترین نتایج را با 3/46 درصد نرخ حمل پیش بینی شده در رنج نسبت ناجوری بین 5/0 تا 2 ارائه می دهد. معادلات مییر- پیتر و مولر و ون راین نیز نتایج خوبی را به ترتیب با 43درصد و 36 درصد نرخ حمل پیش بینی شده ارائه می دهد. فولادفر و همکاران (1391) در پژوهشی به بررسی آزمایشگاهی و مدل سازی تنش برشی بحرانی برای فرسایش سطحی بستر چسبنده پرداخته اند که نتایج نشان داد تنش برشی دارای حساسیت زیادی نسبت به متغییر بعد فرکتالی می باشد و نشان داده شده که هر چه بعد فرکتالی کمتر باشد، تغییرات تنش برشی بحرانی برای فرسایش نسبت به تغییرات غلظت حجمی کمتر خواهد بود. رضایی مقدم و همکاران(1391) در پژوهشی به تحلیل وضعیت پایداری مجرای رودخانه ی قزل اوزن با استفاده از روش های تنش برشی، شاخص مقاومت نسبی بستر و مطالعات صحرایی پرداخته اند. که این بررسی نشان داده است روش های تنش برشی و شاخص مقاومت نسبی بستر، تنها در مجرا های آبرفتی اهمیت داشته و به نتایج آنها می توان استناد کرد. اصغری سراسکانرود و زینالی(1392) به بررسی مقاومت ذرات رسوبی در برابر فرسایش رودخانه سر اسکند چای هشترود با استفاده از روش های ریاضی پرداخته اند که نتایج پژوهش نشان داده که در اکثر نقاط رودخانه میزان تنش برشی بحرانی بیشتر از تنش برشی کل می باشد که نشان دهنده این است که ذرات رسوبی موجود در بستر رودخانه نسبت به فرسایش در بستر رودخانه مقاوم می باشند بنا بر این رودخانه حالت پایداری دارد و میزان تغییرات رودخانه در آینده بسیار کم و به صورت پایدار خواهد بود. بروجنی و همکاران (1392) به بررسی آزمایشگاهی تنش برشی آستانه نهشته شدن رسوبات رسی و سیلتی در غلظت های مختلف پرداخته اند، که نتایج نشان داده که تنش برشی در آستانه نهشته شدن رسوبات. معلق ریز دانه مورد آزمایش برابر 01/1نیوتن بر متر مربع و تنش برشی آستانه برای نهشته شدن کامل این رسوبات 038/0 نیوتن بر متر مربع به دست آمد و در نهایت محاسباتی برای محاسبه تنش برشی آستانه؛ غلظت تعادلی و درصد نهشته شدن ارایه شد. مددی و همکاران(1393) در پژوهشی به آشکار سازی تغییرات بستر رودخانه وبررسی مورفولوژی رودخانه ای با توجه به تآثیر ساختار های زمین شناسی در رودخانه زرینه رود پرداخته اند. که میزان فرسایش در کناره و بستر رودخانه در هر مقطع را به دست آورده اند. 

 

 

مواد و روش

منطقه مورد مطالعه

حوضه آبخیز گلالی در فاصله 43 کیلومتری شهرستان قروه، در جنوب شرقی استان کردستان و در موقعیت جغرافیایی 47 درجه و 59 دقیقه تا 47 درجه و 48 دقیقه طول جغرافیایی ، 34درجه و 56 دقیقه تا 35 درجه و 2 دقیقه عرض جغرافیایی واقع شده است (شکل1). مساحت حوضه 27/75 کیلو متر مربع، مرتفع ترین قله آن 3124 متر و کمترین ارتفاع آن 1827 متر است. شش منطقه مسکونی به نام های گلالی، شیروانه، تکیه علیا، تکیه سفلی، درمه و قورمه دره به اضافه یک معدن آهن و یک کارخانه ماسه شویی و آسفالت در منطقه مورد مطالعه قرار دارد. دمای متوسط سالانه منطقه 3/16 درجه سانتیگراد، متوسط حداکثر دمای ماهانه در مرداد ماه برابر 8/28درجه سانتیگراد و متوسط حداقل دمای ماهانه 6/11 درجه سانتیگراد است. مقدار میانگین بارش سالانه منطقه 4/418 میلیمتر در سال می باشد. اقلیم منطقه مورد مطالعه در روش آمبرژه، اقلیم نیمه خشک سرد است. بازه مورد مطالعه در این مقاله به طول 3 کیلومتر و از پل کنار روستای شیروانه تا پل بعدی که در روستای گلالی قرار دارد.

 

شکل1: حوضه گلالی و بازه مورد مطالعه

روش مطالعه:

جهت بررسی تحرک ذرات بستر رودخانه قروه تعداد 3 مقطع عرضی در بازه مورد نظر از پل روستای شیروانه تا پل روستای گلالی انتخاب شد (شکل1) و اقدام به تهیه مقاطع عرضی کانال با استفاده از متر پارچه ای و شاخص نقشه برداری گردید و پارامتر هایی نظیر عرض دبی لبالبی، عرض دبی فعال، عمق متوسط دبی لبالبی، شیب بستر، اندازه ذرات بستر و سرعت جریان و در پایان دبی اندازه گیری شد. بر پایه داده های اندازه گیری شده از مقاطع، پارامترهای مورد نیاز در اندازه گیری تنش برشی بحرانی و دبی واحد به شرح ذیل محاسبه گردید. برای بررسی اندازه رسوباتی که بیشترین فراوانی رادارند از روش شمارش پبل (ریگ) که به‌وسیله ولمن در سال 1952 ارائه شده است، استفاده شد. در این روش 100 ذره به صورت اتفاقی انتخاب شده و محور متوسط(d) آن‌ها اندازه‌گیری شد. سپس داده‌های به دست‌ آمده از اندازه ذرات به صورت نمودار تجمعی ترسیم شد و 50D  (اندازه ذراتی که 50 درصد نمونه مساوی یا کوچکتر از آن

 

 

 

 

 

هستند) محاسبه می‌شود (حسین زاده و همکاران، 1394). همچنین اندازه ذرات D95(اندازه ذراتی که 95 درصد نمونه مساوی یا کوچکتر از آن هستند)،  D84(اندازه ذراتی که 84 درصد نمونه مساوی یا کوچکتر از آن هستند) و D16 (اندازه ذراتی که 16 درصد نمونه مساوی یا کوچکتر از آن هستند) رسوبات در هر سه مقطع عرضی تعیین شد. بعد از بررسی اولیه مشخصات استخراج شده از مقاطع عرضی، مشخص شده که

دربازه مورد مطالعه شیب بستر بیش از یک درصد است، غوطه‌وری نسبی کمتر از 10، و بستر کانال دارای

ریفل و واریزه‌های چوبی می‌باشد، بنابراین یکی از روش‌های مناسب برای آنالیز تحرک و پایداری ذرات بستر

روش دبی واحد بحرانی می‌باشد. علاوه براین روش، با استفاده از معادله اصلاح شده تنش برشی بحرانی، تنش برشی بحرانی برای ذرات D16، D50، D84، D95 نیز محاسبه گردید. علت انتخاب این ذرات برای آنالیز این است که D84، D95 ذرات کلیدی می‌باشند که زمانی که شروع به حرکت می‌کنند نشان‌دهنده حرکت بیشتر اجزاء بستر بوده و در واقع شروع تغییر در ساختار بستر را نشان می‌دهند. D50  میانه ذرات بوده و متوسط تنش برشی را نشان می‌دهد و همچنین اندازه این ذرات برای اکثر محاسبات لازم است. اندازه D16 در بسترهای نسبتاً پایدار اهمیت دارد.

معادله تنش برشی بحرانی اصلاح شده(Komar 1987, 1997; Komar and Carling 1991) براساس رابطه بین اندازه ذره مورد آنالیز و D50 رسوبات بستر به صورت زیر بیان می‌شود:

رابطه شماره 1: τci= τD50(ɤs-ɤ) Di0.3D500.7

که در آن τci تنش برشی بحرانی برای ذره مورد نظر(N/m2)، τD50 پارامتر بی بعد شیلدز برای ذرات D50(از جدول شیلدز به دست می‌آید)، ɤs وزن مخصوص رسوب، ɤ وزن مخصوص آب، Di اندازه ذره مورد آنالیز(m(، D50 اندازه ذراتی که 50 درصد نمونه مساوی یا کوچکتر از آن هستند(m) می‌باشد.

پس از محاسبه تنش برشی بحرانی برای چهار اندازه مختلف، تنش برشی بحرانی در دبی لبالبی(τbf) با استفاده از رابطه شماره 2 محاسبه شد. در صورتی که تنش برشی بحرانی برای ذره مورد نظر بیشتر از تنش برشی در دبی لبالبی باشد، ذره مورد نظر پایدار خواهد بود و در صورتی که کوچکتر باشد ذره در دبی لبالبی حرکت خواهد نمود.

رابطه شماره 2: τbf = gRS

τbf  تنش برشی بحرانی در عرض دبی بانکفول(N/m2)، R عرض کانال در حالت دبی بانکفول(m)، S شیب بستر رودخانه(m/m) و g وزن مخصوص آب است.

برای کانال‌های با شیب بیشتر از 1 درصد (01/0 متر بر متر)، جایی که نسبت عمق جریان به اندازه ذرات کمتر است(غوطه‌وری نسبی: R/D50 کمتر از 10)، عمق آب کاملاً متغیر است، زیرا سنگ‌های بزرگ یا واریزه‌های چوبی نزدیک سطح، عمق آب را کاملاً تحت تأثیر قرار می‌دهند. (Bathurst et al,1987) برای چنین شرایطی با استفاده از داده‌های فلوم معادله‌ای برای پیش‌بینی دبی واحد بحرانی موردنیاز برای آغاز حرکت ذرات D50 ارائه نمودند (رابطه 3):

رابطه شماره 3: qcD50 =

 

که در آن qcD50 دبی واحد بحرانی برای حمل ذرات D50(m2/s)، D50 معادل قطر میانه و s شیب

بستر(m/m) می‌باشند.

 (Bathurst et al1987) برای محاسبه دبی واحد بحرانی برای آغاز حرکت ذرات در اندازه‌های مختلف رابطه شماره 4 را ارائه

نمودند:

رابطه شماره 4: qcD50(Di / D50)b qci =        

که در آن qci دبی واحد بحرانی برای حمل ذرات در اندازه مورد نظر(m2/s)، qcD50 دبی واحد بحرانی برای حمل ذرات D50(m2/s)، Di اندازه ذره مورد نظر و D50 اندازه ذراتی که 50 درصد نمونه مساوی یا کوچکتر از آن هستند می‌باشد.

توان b دامنه‌ای از اندازه ذراتی است که ترکیب بستر کانال را تشکیل می‌دهند. با وارد کردن توان b در محاسبات تأثیر اندازه ذرات کوچکتر مخفی در پناه ذرات بزرگتر نیز در محاسبات وارد می‌شود. توان b از رابطه 5 محاسبه می‌شود:

رابطه شماره5: b = 1.5(D84/D16)

که در آن  D84 اندازه84 درصد اندازه ذرات و D16 اندازه 16 درصد ذرات است.

برای بازه مورد مطالعه پس از اینکه مشخص شد نسبت R/ D50 کمتر از 10 می‌باشد(جدول شماره1) برای محاسبه دبی واحد بحرانی با این روش مراحل زیرا اجرا شده است:

1-       ابتدا با استفاده از رابطه شماره 1 دبی واحد بحرانی با استفاده از قطر ذرات D50 و شیب بستر محاسبه می شود.

2-       با استفاده از رابطه شماره 5 مقدار توان b را محاسبه نموده و سپس با استفاده از رابطه شماره 3 دبی واحد بحرانی برای اندازه ذرات D16، D85 و D95 محاسبه می گردد.

3-       در ادامه دبی واحد درون کانال در هر مقطع با استفاده از رابطه شماره 6 محاسبه می‌ شود. و در نهایت برای تعیین اینکه اندازه ذرات داده شده در دبی واحد درون کانال حرکت خواهد کرد یا خیر؟ دبی واحد بحرانی ذره مورد نظر با دبی واحد درون کانال فعال(q) مقایسه می‌شود. در صورتی که دبی واحد بحرانی ذره(qci) کمتر از دبی واحد درون کانال فعال(q) باشد، ذره کشیده خواهد شد. اگر دبی واحد بحرانی ذره (qci) بیش از دبی واحد(q) اعمال شده بر روی ذره باشد، ذره حرکت نخواهد کرد.

رابطه شماره 6:                                 q = Q/w

در این رابطه q دبی واحد درون کانال(دبی در واحد عرض)(m2/s)، Q دبی حداکثر لحظه‌ای رودخانه(m2/s) و W عرض کانال فعال(m) در مقطع عرضی مورد نظر می‌باشد.

پس از تعیین نسبت R/ D50، (غوطه وری نسبی) مشخص شد مناسب‌ترین روش برای آنالیز تحرک و پایداری ذرات بستر رودخانه گلالی قروه در بازه مورد نظر روش دبی واحد بحرانی می‌باشد. علاوه بر این روش با استفاده از معادله تنش برشی بحرانی، میزان تنش برشی بحرانی برای ذرات در چهار اندازه مختلف نیز محاسبه گردید. بدین منظور متغیرهای مورد نیاز برای تعیین دبی واحد بحرانی و تنش برشی بحرانی در سه مقطع عرضی در این بازه برداشت شد (جدول شماره1)  با توجه به جنس رسوبات بستر در بازه مورد مطالعه که کوارتز است وزن مخصوص رسوب معادل 65/2 و چگالی آب نیز با توجه به دمای آب رودخانه مورد نظر که در هنگام نمونه برداری اندازه گیری شده است 1 در نظر گرفته شده است.

   

بحث و نتایج

 

 

    

 

جدول1: مشخصات مقاطع عرضی رودخانه گلالی قروه

 

شماره مقطع	شیب بستر(S)	عمق جریان(R)	عرض کانال فعال(W)	D16(mm)	D50(mm)	D86(mm)	D95(mm)	R/ D50
1	023/0	5/0	5/7	26	64	128	180	62/6
2	024/0	49/0	4/9	32	90	180	256	16/5
3	017/0	4/0	9/5	26	90	180	256	27/4

 

با استفاده از رابطه شماره 1 تنش برشی بحرانی برای سه اندازه ذره کلیدی محاسبه شد(جدول شماره2). مقدار تنش برشی ذرات با قطر D50 که برای محاسبات اولیه لازم بود از جدول پارامتر شیلدز محاسبه شد. همچنین تنش برشی کانال در حالت دبی لبالبی نیز با استفاده از رابطه شماره 2 محاسبه شده است.

جدول2: تنش برشی محاسبه شده برای اندازه‌های مختلف ذرات در مقاطع عرضی رودخانه گلالی قروه

شماره مقطع	τD16(N/m2)	τD84(N/m2)	τD95(N/m2)	τbf(N/m2)
1	79/0	26/1	87/1	99/1
2	22/1	95/1	18/2	27/2
3	1/1	95/1	18/2	33/1

 

 

 

همانطور که از جدول شماره 2 قابل ملاحظه است در مقطع عرضی شماره 1 تنش برشی بحرانی ایجاد شده توسط دبی لبالبی از تنش برشی لازم برای شروع حرکت ذرات با قطر D16، D84 و D95 بیشتر است. بنابراین تمامی ذرات دارای این اندازه‌ها یا کوچکتر از آن در دبی لبالبی این مقطع حرکت خواهند کرد. در مقطع عرضی شماره 2 تنش برشی بحرانی ایجاد شده توسط دبی لبالبی رودخانه بیشتر از تنش برشی لازم برای شروع حرکت ذرات با قطر D16،D95 و D84 می‌باشد، بنابراین در دبی لبالبی ذرات ذکر شده حرکت خواهند کرد. در مقطع عرضی شماره 3 تنش برشی ایجاد شده توسط دبی لبالبی این مقطع عرضی از تنش برشی لازم برای شروع حرکت ذرات با اندازه, D84 D95 کمتر است، بنابراین در دبی لبالبی ذرات با اندازه

D95 و D84 بزرگتر از آن پایدار خواهند بود. در ادامه کار با استفاده از روابط شماره 3 تا 6 مقادیر qcd50، qcD16، qcD84، qcD95 و q محاسبه شد(جدول شماره 3).

 

 

 

جدول3: دبی واحد بحرانی محاسبه شده برای اندازه‌های مختلف ذرات در مقاطع عرضی رودخانه گلالی قروه

شماره مقطع	Q(m2/s)	qcD16(m2/s)	qcd50(m2/s)	qcD84(m2/s)	qcD95(m2/s)
1	2/1	42/0	51/0	59/0	94/0
2	15/1	69/0	83/0	94/0	1
3	14/1	02/1	22/1	35/1	42/1

 

 

با توجه به نتایج به‌دست‌آمده مشخص شد که در مقطع عرضی شماره 1 دبی واحد بحرانی لبالبی در عرض رودخانه در مقایسه با دبی واحد بحرانی لازم برای حرکت تمامی ذرات مورد بررسی (اندازه‌های D16، D50، D84 و D95) بیشتر است. بنابراین ذراتی که اندازه آنها مساوی یا کمتر از این اندازه‌ها باشد در بستر رودخانه پایدار نخواهند ماند و در دبی واحد بحرانی لبالبی رودخانه حرکت می‌کنند. در مقطع شماره 2 دبی واحد بحرانی لبالبی در عرض رودخانه در مقایسه با دبی واحد بحرانی لازم برای شروع حرکت ذرات با اندازه‌های D16، D50 ، D84 و D95بیشتر است. بنابراین ذراتی که اندازه آنها مساوی یا کمتر از این اندازه‌ها باشد در بستر رودخانه پایدار نخواهند ماند و در دبی واحد بحرانی لبالبی رودخانه حرکت می‌کنند. در مقطع شماره 3 دبی واحد بحرانی لبالبی در عرض رودخانه در مقایسه با دبی واحد بحرانی لازم برای شروع حرکت ذرات با اندازه‌ی D16، بیشتر است. در این مقطع عرضی دبی بحرانی لازم برای شروع حرکت ذرات با اندازه D95 D50 و D84  و بزرگتر از آن بیشتر از دبی واحد بحرانی در عرض مقطع می‌باشد و در شرایط دبی لبالبی این ذرات پایدار خواهند ماند و حرکت نمی‌کنند.

همان‌گونه که قابل مشاهده است در مقطع عرضی شماره 1و2 ذرات نسبت به مقطع عرضی شماره 3 ناپایدارترند و در شرایط دبی لبالبی تمامی ذرات مورد بررسی شروع به حرکت خواهند کرد، بنابراین در دبی لبالبی شرایط کلی بستر رودخانه در این مقطع تغییر خواهد کرد. با توجه به شرایط مقطع شماره 1و2 انتظار می‌رود که به‌علت شیب بیشتر بستر و عرض کمتر کانال فعال در این مقطع ذرات در اندازه‌های کوچکتر نسبت به مقطع عرضی شماره 3 شروع به حرکت کنند. در مقطع شماره 3 نیز در دبی لبالبی ذرات  کوچکتر حرکت خواهند کرد و ذرات بزرگتر از D50 ثابت خواهند ماند.

 

بحث و نتیجه‌گیری:

درک حرکت رسوب و حساسیت کانال نسبت به فرسایش یا رسوبگذاری، ضرورتی مهم برای مدیریت رودخانه است. روش‌های متفاوتی برای تعیین آستانه حرکت ذرات رسوبی بستر رودخانه توسط محققین ارائه شده است از جمله این روش‌ها روش تنش برشی بحرانی و روش دبی واحد بحرانی است که در این مقاله به این دو روش پرداخته شده است. متاسفانه، تلاش برای محاسبه یا اندازه‌گیری مقادیر تنش برشی در رودخانه‌ها، بصورت پیچیده‌ای توسط زبری بستر کانال و نوسانات سرعت و آشفتگی جریان تحت تأثیر قرار می‌گیرد(whole,2000). با این حال هنوز هم یکی از دقیق‌ترین روش‌های برآورد آستانه حرکت ذرات رسوبی بستر رودخانه است. جهت بررسی تحرک ذرات بستر رودخانه گلالی قروه مقطع عرضی در بازه مورد نظر انتخاب شد و پارامترهای مورد نیاز برای برآورد تنش برشی بحرانی و دبی واحد بحرانی در هر سه مقطع استخراج شد. در این مقاطع تنش برشی بحرانی و دبی واحد بحرانی برای شروع حرکت چهار اندازه کلیدی از ذرات محاسبه شد و با تنش برشی ایجاد شده توسط دبی لبالبی و دبی بحرانی در شرایط دبی لبالبی در واحد عرض کانال مقایسه گردید. نتایج نشان داد که در مقطع عرضی شماره 1و 2 در دبی لبالبی رودخانه تمامی ذرات مورد بررسی ناپایدار می‌باشند و حرکت خواهند کرد. درحالی‌که در مقطع عرضی شماره 3 ذرات با اندازه D95 D50 و D84پایدار هستند و سایر ذرات در دبی لبالبی در واحد عرض رودخانه حرکت خواهند کرد. این امر به شیب بیشتر بستر و عرض کمتر کانال رودخانه در مقطع شماره 1و2برمی‌گردد. این پژوهش نشان داد که استفاده از هر دو روش تنش برشی بحرانی و دبی واحد بحرانی نتایج مشابهی در ارتباط با آستانه حرکت ذرات به دست می‌دهد و این خود تأییدی بر استفاده از این روش‌ها برای مطالعات پایداری رودخانه است. استفاده از روش تنش برشی بحرانی و دبی واحد بحرانی برای بررسی آغاز حرکت ذرات رسوبی کانال رودخانه‌ها روش مناسب، مقرون ‌به‌صرفه و ساده‌ای است که به‌راحتی می‌توان از نتایج آن برای برنامه‌ریزی‌های پایدارسازی کانال و پیشگیری از فرسایش کناری استفاده نمود.

 

 

منابع

1-حدادچی، آرمان؛ امید، محمد حسین؛ دهقانی، امیر احمد. 1390. ارزیابی معادلات تجربی برآورد دبی بار بستر در رودخانه کوهستانی با بستر شنی (مطالعه موردی: رودخانه چهل چای در استان گلستان)، مجله پژوهش های آب و خاک جلد 18، ش 3 صص149-177

1خواستار بروجنی، میلاد؛ صمدی بروجنی، حسین؛ اسماعیلی، کاظم. 1392. بررسی آزمایشگاهی تنش برشی آستانه نهشته شدن رسوبات رسی – سیلتی در غلظت های مختلف. مجله پژوهشی آب ایران سال هفتم. ش13. صص131-138

2-حسین زاده، محمدمهدی؛ اسماعیلی، رضا. (1394). ژئومورفولوژی رودخانه‌ای. مفاهیم، اشکال و فرآیندها. انتشارات دانشگاه شهید بهشتی. تهران.

 4-رضایی مقدم، محمد حسین؛ ثروتی، محمدرضا؛ اصغری سراسکانرود؛ صیاد.1391. تحلیل وضعیت پایداری مجرای رودخانه قزل اوزن با استفاده از روش های تنش برشی، شاخص مقاومت نسبی بستر و مطالعات صحرایی، مجله پژوهش های کمی، ش 1، صص 33-46

5-سر اسکانرود، صیاد؛ زینالی، بتول. 1392،بررسی مقاومت ذرات رسوبی در برابر فرسایش رودخانه سر اسکند چای هشترود با استفاده از روش های ریاضی،فصل نامه علمی پژوهشی پژهش های فرسایش محیطی هرمزگان. سال سوم. ش12،صص42-30

6-فولادفر، حسام؛ شفاعی بجستان، فتحی مقدم، منوچهر. 1391. بررسی آزمایشگاهی و مدل سازی تنش برشی بحرانی برای فرسایش سطحی بستر چسبنده، فصلنامه علمی پژوهشی مهندسی آبیاری و آب، سال دوم ش8، صص 78-94

7-قمشی، مهدی؛ وجدانی، ناصر.1384. تنش برشی بحرانی فرسایش رسوبات چسبنده و نقش آن در طراحی کانال های رو باز، مجله بین المللی علوم مهندسی، ش3. صص 83-93. دانشگاه علم و صنعت ایران

8-مددی، عقیل؛ بهشتی جاوید، ابراهیم؛ فتحی، محمد حسین.1393. آشکار سازی تغییرات بستر رودخانه و بررسی مورفولوژِی رودخانه ای با توجه به تآثیر ساختار های زمین شناسی (مطالعه موردی: زرینه رود)، مجله هیدروژئومورفولوژی، ش2. صص 25-42

9-نیری، ه، 1389. تحلیل دینامیک وشکل مجرا در حوضه آبریز مهاباد، رساله دکتری جغرافیای طبیعی، گرایش جغرافیای طبیعی، گرایش ژئومورفولوژی، دانشکده جغرافیا-دانشگاه تبریز.

 

10-Bizzi.,s and Lerner.,D.N. 2015, The use of stream power as an indicator of channel sensitivity to erosion and deposition processes, River research and application River Res.applic.31: 16-27

11- Coryat., M, 2014, Analysis of the Bank Assessment for Non-point Source consequences of Sediment (BANCS) Approach for the Prediction of Streambank Stability and Erosion along

stony Clove Creek in the Catskills, Master Theses, Syracuse University.

12-Mighui, Y., Hongyan, W.,Yanjie, L., and Chunyan, H., 2010, Study on Stability of Noncohesive River Bank, International Journal of Sediment Research, Vol. 25, No. 4,pp.391-398

13- Rinaldi M; Simoncini C; Piegay H.(2009). Scientific design strategy for promoting sustainable sediment management: the case of the magra river(central northern Italy) river research and application 25. pp. 607-625

14-Tokaldany E.A., Darby S. E., and P. Tosswell.,2007. Coupling bank stability and deformation models to predict equilibrium bad topography in river bends. Journal of Hydraulic Engineering, Vol. 133, No10,pp 1167-1170

15-Wallerstein N.P; Soar P.J; Thorne C.R.(2006). River energy auditing scheme(REAS) for catchment flood management planning, International conference on fluvial hydraulics. Pp. 1923-1932. Lisbon portugal

16-Wohl, E. (2000). Mountain Rivers. American Geophysical Union. Washington, D.C.

 

 

 

 

Channel bed unsteady analysis using critical shear stress and critical unit discharge (Case Study: Rivers Galali, Kurdistan)

 

 

Abstract

 

One of the key issues of erosion and stability of the river is the beginning the movement of sediment particles. River sediment yield caused water turbidity, nutrients and pollution problems of water, buried in the water diversion facilities as well as filling the water storage lakes. Flow that will start moving particle, called the critical flow. The initial move threshold sediment particles has been investigated in several different methods, most important of which are the average critical velocity, critical shear stress and critical unit discharge. To investigate the stability of the Galali river bed sediments in the city of Qorveh, according to river conditions used the critical shear stress and critical unit discharge methods. To determine initial move threshold sediment particles in the bed of the river were two cross-section of the river, two cross-section conditions were assessed. After that, the critical shear stress and critical unit discharge for four particle size D16, D50, D84 and D95 in each cross-section was calculated, and showed that in cross-section (1,2) in the bankfull discharge conditions all particles with sizes D84 and D95 begin to moving, and so in the bankfull discharge, sediment conditions will generally change. In cross-section (3) particles with a diameter D95, D16, D50 and the more it will be remain stable in the river bed, and other smaller particles of this size are starting to move. Due to further slope of the bed and width less active channel in the cross-section (1,2) has caused that all particles are unstable. Finally, the results from both methods were highly concordant, and this confirms that the use of these methods for stabilizes studies and sediment yield of river are appropriate.

Key words: Move threshold sediment particles, Critical shear stress, critical unit discharge, Galali river.

---