همانگونه که آلودگی نفتی در دریاها و اقیانوس ها بارها رخ داده است، احتمال وقوع این حوادث در رودخانه ها و دریاچه ها نیز وجود دارد. به عنوان مثال طی سال های 1979 تا 1986، 367 مورد آلودگی نفتی در رودخانه ها در سطح جهان گزارش شده است که در 65% موارد حجم این آلودگی ها کمتر ازm3 100 می باشد. وقوع آلودگی نفتی علاوه بر ضایع کردن آب شرب و کشاورزی و به خطر انداختن سلامت انسان ها، موجب وارد آمدن خسارت به محیط زیست و از بین رفتن زیستگاه جانوری و گیاهی منطقه می شود. لذا تعیین نحوه انتشار آلودگی از نظر مسائل زیست محیطی و همچنین در مدیریت بحران حائز اهمیت می باشد هدف اصلی این پایان نامه تحلیل نحوه مدلسازی انتقال و انتشار ماده آلاینده پیچیده (نفت) در محیط رودخانه با تمرکز بر حادثه روی داده در فروردین 87 رودخانه زاینده رود بعنوان مطالعه موردی است. در این پایان نامه از مدل wasp با توجه به قابلیت های مناسب آن نظیر توانایی شبیه سازی مواد آلی سمی استفاده گردیده است. حادثه نفتی رودخانه زاینده رود یکی از سوانحی است که طی آن در اثر برخورد بولدوزر راهسازی با لوله انتقال نفت از منطقه مسجد سلیمان به پالایشگاه نفت اصفهان حدود 4 میلیون لیتر نفت خام از لوله خارج و عمده آن به رودخانه زاینده رود وارد شده است. همچنین در این پایان به بررسی روش های مختلف جداسازی نفت از محیط رودخانه پرداخته شده است. براساس نتایج حاصله از شبیه سازی حادثه مشخص شد که دو پارامتر دبی جریان و دمای هوا در میزان غلظت آلودگی بیشترین تأثیر را دارا می باشند و دیگر آنکه برش های سنگین تر نفت اثرات نامطلوب و طولانی تری از خود در محیط رودخانه برجای خواهند گذاشت. بنابراین در این پایان نامه به منظور حفظ سلامت و پاکیزگی آب و حفاظت از محیط زیست، تلاش جدی برای افزایش ایمنی و جلوگیری از وقوع حوادث مشابه (به عنوان پیشگیری) توصیه شده است. اما در صورت وقوع حادثه باید نسبت به انجام فعالیت های به موقع و اصولی نظیر اعلام به موقع به تصفیه خانه و قطع جریان آب توسط سد جهت انجام عملیات پاکسازی اقدام شود که در غیر اینصورت پیامدهای نامطلوب حادثه افزایش خواهد داشت. براساس نتایج حاصل از سناریوهای تعریف شده نیز مشخص شد که در صورت وقوع حادثه در فصول تابستان و پائیز با بهترین و بحرانی ترین حالت روبرو خواهیم بود.

در این پایان نامه، از مدلqual2e به منظور شبیه سازی کیفیت آب رودخانه استفاده شد. متغیرهای کیفی انتخاب شده برای شبیه سازی کیفی، بر اساس آمار در دسترس عبارتند از do, bod, tds,nh3, no2, no3,t . دوره های زمانی اردیبهشت ماه 1380 و بهمن ماه 1379 به ترتیب برای واسنجی و صحت سنجی نتایج استفاده گردید تا ضرایب مورد نیاز مدل بدست آید. پس از شبیه سازی کیفی، با استفاده از زیربرنامه qual2e-uncas (که در مدل qual2e به منظور تحلیل عدم قطعیت تعبیه شده است) سهم متغیرهای ورودی در عدم قطعیت نتایج شبیه سازی مشخص گردید. این کار با استفاده از روش تحلیل خطای مرتبه اول (f.o.e.a) انجام شده است. در کنار استفاده از روش f.o.e.a، شبیه سازی مونت کارلو نیز انجام شد و با مقایسه انحراف معیار بدست آمده از این روش با روش قبلی، میزان خطای فرض خطی بودن در روش f.o.e.a مشخص گردید.

هیدرولیک انتقال رسوبات رودخانه دانشی است که چگونگی حرکت، فرسایش و رسوبگذاری مواد رسوبی در کانالها و رودخانه ها را مورد بحث قرار می دهد. این سه فرآیند فرسایش، انتقال و ته نشینی توسط محققین بوسیله ی مدل های قطعی و استوکستیک مورد تحلیل قرار گرفته است. در مدل های قطعی، معادلات ارائه شده بر اساس مطالعات آزمایشگاهی و مشاهدات صحرایی می باشند اما در مدل های استوکستیک، اساس بر انتخاب متغیر تصادفی و بکار گیری تئوری استوکستیک در مکانیک انتقال رسوب می باشد. متغیر تصادفی متغیری است که دارای عدم قطعیت (تابع چگالی احتمال) بوده و مقدار آن در طول کانال پیوسته در حال تغییر است. استوکستیک فرآیندی است که وابسته به متغیرهای تصادفی است. بنابراین اگر رابطه ای دارای متغیر تصادفی باشد آن رابطه یک رابطه ی استوکستیک است. تحقیقات نشان می دهد که مدل های قطعی خطای بسیار داشته اند. علاوه بر این، دبی رسوب محاسبه شده توسط این روش ها، اختلاف زیادی با هم دارد. یک مدل استوکستیک می تواند با استفاده از تاثیر پارامترهای متفاوت، فرآیند انتقال رسوب را بهتر از مدل های قطعی شبیه سازی کند. بنا به نظر برخی محققین حرکت ذرات رسوب یک پدیده ی قطعی نیست و تنها یک مدل استوکستیک می تواند این رفتار تصادفی ذرات رسوب را توجیه کند. هم چنین در شرایط سیلابی که نیاز مبرم به دستیابی دبی رسوب می باشد، روش های قطعی، انتقال رسوب را به خوبی مدل نمی کنند. این پایان نامه به بررسی روش استوکستیک در مدل کردن انتقال رسوب در شرایط سیلابی و غیر سیلابی می پردازد. برای این منظور از رابطه ی کلی استوکستیک در انتقال بار بستر که توسط سایر محققین به اثبات رسیده، استفاده شده است. در این رابطه سرعت برشی و قطر متوسط ذرات رسوب، متغیر تصادفی می باشند. دو پارامتر ارتفاع پرش و سرعت حرکت ذره، توسط محققین مختلف با روش های گوناگون مورد محاسبه قرار گرفته است. این پایان نامه روش نوینی را برای دستیابی به این دو پارامتر ارائه می کند. برای این منظور متغیرهای تصادفی بیشتری مورد استفاده قرار می گیرد. این متغیرهای تصادفی عبارتند از: شیب کف، شعاع هیدرولیکی، سرعت برشی در کف، سرعت لحظه ای در نزدیکی بستر، عدد رینولدز ذرات، ضریب کششی، نیروی کششی و مساحت تصویر ذره. ابتدا هدف دستیابی به تابع چگالی احتمال قطر متوسط ذرات رسوب، با توجه به دبی های مختلف است. این عمل تنها با نقشه برداری از بازه مورد مطالعه و اندازه گیری سرعت در یکی از این مقاطع میسر است. سرعت حرکت ذره با توجه به احتمال پرش بین دو نقطه و محاسبه ی طول و زمان این پرش قابل محاسبه است. برای دستیابی به احتمال پرش بین دو نقطه، به آنالیز احتمالاتی آستانه حرکت رسوبات در شیب های تند پرداخته می شود. سه پارامتر ارتفاع، طول و زمان یک پرش نیز، در شرایط غیر سیلابی از معادله مسیر حرکت ذره به کمک معادلات استوکستیک قابل تعیین است. این پارامترها در شرایط سیلابی به کمک معادلات استوکستیک و فرایند پواسون بدست می آیند. در واقع این مدل با ترکیب احتمالات و معادلات استوکستیک به محاسبه دبی بار بستر می پردازد. توانایی این مدل را می توان در نظر گرفتن بستر متحرک و متغیرهای تصادفی زیاد عنوان کرد. در این پایان نامه رودخانه امامه به عنوان مطالعه موردی انتخاب شد. دلیل انتخاب رودخانه امامه به عنوان مطالعه موردی، موجود بودن آمار مکفی طی مطالعات انجام شده در این رودخانه در گذشته، برای انجام این پایان نامه بوده است. دبی بار بستر محاسبه شده به کمک روش شبیه سازی مونت کارلو، برای اثبات صحت خود، مورد تجزیه و تحلیل بیشتر قرار گرفت. در ابتدا با سایر روش قطعی و احتمالاتی مورد مقایسه قرار گرفت. دبی بار بستر محاسبه شده با روش انیشتن مطابقت خوبی دارد. سپس با استفاده از آمار موجود، تعداد وقوع یک سیل بخصوص با توجه به هیدروگراف سیل فرضی، مورد محاسبه قرار گرفت. با توجه به متفاوت بودن قطر ذرات و متفاوت بودن جابجایی برای آنها بعد از هر سیل، تعداد وقوع هیدرگراف های سیل فرضی یکسانی برای هر قطر پیش بینی شد. در انتها متغیر تصادفی پرش ذره رسوب، با نتایج سایر محققین مورد مقایسه قرار گرفت. مدل ارائه شده، در تمامی مقایسات مورد قبول قرار گرفت. در پیشنهاد برای ادامه موضوع توسط سایر دانشجویان، توصیه به بررسی دقیق تر فرضیات مورد استفاده در این پایان نامه (مثلا در نظر گرفتن اثر برخورد ذرات به یکدیگر) می باشد.