شنبه، اسفند ۲۱، ۱۳۸۹
سلام
متاسفانه باید عرض کنم که بعلت فیلتر شدید بنده مدتی هست که نمیتونم بیام و مطلب بنویسم و کلاً Blogger بسته بود.
به هر حال امروز بالاخره به نحوی موفق شدم که وارد بشم و از خبرهای جدید مطلعتون کنم.
بله، خوشبختانه بنده در مورخ 12/12/1389 مصادف با 03/03/2011 موفق شدم که از پایان نامه ام دفاع کنم.
استاد راهنما: جناب آقای دکتر رضایی
استاد مشاور: جناب آقای دکتر یوسف زاده
داوران: جنای آقای پروفسور بهروش و جناب آقای دکتر قادری
مدیر گروه و ناظر جلسه: جناب آقای دکتر ارمغانی
متاسفانه چون مقاله نداشتم فقط تونستم نمره کامل تو قسمت خود پایان نامه رو بدست بیارم، 18
به هر حال ناراضی هم نیستم. اخیراً هم دارم یک سری ویرایش نهایی انجام میدهم که دیگه تکثیر کنم و برم تو خط صحافی.
بعداً حتماً البته اگر این فیلتراسیون هم اجازه بده، مطلب می نویسم مخصوصاً در مورد نتایج نهایی، عکسها ، نمودارهاو...
با تشکر
واحددوست
دوشنبه، بهمن ۱۱، ۱۳۸۹
فصل بندي جديد
همانطور كه قبلاً گفته بودم بعد از تايپ قسمتي از پايان نامه، اونو بردم دادم كه آقاي دكتر رضايي بررسي كنند.
ايشون نظرشون اينه كه تغييراتي به شرح زير بديم و فصل بندي رو اينطوري كنيم:
- فصل اول: مقدمه
- فصل دوم: تئوري سدهاي زيرزميني و پيشينه تحقيق
- فصل سوم: تعريف ومعرفي Modflow
- فصل چهارم:معرفي طرح و ابعاد مساله
- فصل پنجم:آناليز و داده ها
- فصل ششم:نتيجه گيري وپيشنهادات
با تشكر
واحددوست
Babak.V.Dst
اتمام مدل
خوشبختانه پس از كلي كلنجار و مشكل برنامه Run شد و مدل را بستم.
نتايج در حد قابل قبولي بود كه با توجه به آناليزهاي انجام شده براي سال پيش بيني به اين نتيجه رسيدم كه عليرغم اينكه در سال پيش بيني تبخير و تعرق بيشتر بوده و آورده Recharge و رودخانه كمتر شده و قنات نيز همچنان در حال تخليه مخزن هستند. باز ميزان بيلان ذخيره منطقه 4 برابر شده كه از حدود 13 متر مكعب در سال به 60 مترمكعب در سال رسيده كه عليرغم كم بودن آن بايد به اين نكته توجه كرد كه اين مخزن مدام در حال تخليه شدن و برآوردن نيازهاي ما مي باشد.
در مورد واريانسها،پيزومترهاي بالا دست من به اعداد خوبي (زير 0.5) رسيدم ولي پيزومتر پايين دست با واريانس حدود 7 كمي دردسر ساز شد.چون در گزارش جهاد نيز بدان اشارهه شده بايد گفت كه احتمالاً در طول مسير (بعد از سد) و يا از بدنه سد نشت وجود دارد. چراكه ارتفاع سطح آب محاسبه شده در پيزومتر بسيار پايينتر از سطح آب برداشت شده مي باشد.
در ادمه مفصلاً در مورد نتايج صحب خواهم كرد.
با تشكر
واحددوست
پنجشنبه، دی ۳۰، ۱۳۸۹
حل مشکلات
بالاخره پس از یک وقفه طولانی، اگرچه از برنامه کمی عقب موندم ولی باید عرض کنم که خیلی از مشکلات حل شد. که گزارش اون به شرح زیره.
- در ابتدا با کمک نمودار نسبت به هر یک از پیزومترها و کل مدل ساخته شده حساسیت سنجی انجام گرفته شد. و با کمک آنها طرح اولیه محدوده هدایت هیدرولیکی تعریف گردید.
- تعیین دقیق محدوده بندی Hc انجام گرفت و مدل در حالت Steady تا حد واریانس 0.65 کالیبره شد.
- از آنجائیکه رواداریهای مجاز برای کنترل رقوم سطح آب برای پیزومترهای موجود در حد 5% توسط دکتر رضایی اعلام شده بود در مورد پیزومتر شماره 2 کمی مشکل وجود داشت که آخر سر با تخفیف حدود 0.001 مورد قبول واقع شد.
- با کمک حساسیت سنجی مدل برای ضریب ذخیره، در ابتدا به کمک Hc و فرمول مربوطه مقادیر اولیه Sy تخمین زده شد.
- بعد از تخمین اولیه Sy، مطابق Hc نمودارهای حساسیت سنجی برای هر یک از پیزومترها و کل مدل رسم گردید.
- تعیین دقیق محدوده های ضریب ذخیره با آزمون و خطا امتحان گردید که در نهایت مدل با واریانس 0.061 مدل شد و بعلت پایین بودن واریانس و محدوده های رواداری مدل مورد قبول واقع شد.
- زمانهای محدوده طرح با توجه به وفور دیتا های موجود و داده سازی از شهریور 83 لغایت مرداد 83 به مدت 365 روز جهت دوره آموزش مدل و کالیبراسیون در نظر گرفته شد.
- دیتاهای مربوط به Recharge و Evapotranspiration با کمک نرم افزار ساخته شده محاسبه گردید و در زمانهای مربوطه وارد مدل گردید.
- تنها مشکل فعلی آن است که مدل زمانها را نمی شناسد و بعد از Run فقط اولین Period رو نمایش می دهد. که البته قبلاً با این مساله برخورد کرده بودم و گمان می کنم که به راحتی قابل حل است.
- اجرای مدل و کالیبراسیون دقیقتر در صورت نیاز
- ایجاد مدل Key جهت حفاظت از اطلاعات
- احداث سد در محل مورد نظر
- فعال کردن Polyline مربوط به زهکش سد و ایجاد ارتباط با قنات امیر کهریزی
- اجرای مدل و کنترل Water budget
- نوشتن فصل 6 و 7 باقیمانده.
- اصلاح و تحریر نهایی پایان نامه و ....
با تشکر واحترام
واحددوست
جمعه، دی ۲۴، ۱۳۸۹
فقط مشکل
امروز رفتم خدمت جناب آقای دکتر رضایی و کمی در مورد پایان نامه حرف زدیم قرار شد برای Hc ضریب حساسیت در بیارم که بهترین باشه ولی متاسفانه مدل اصلاً درست جواب نمیده، البته باید بگم که variance در حد زیر عدد 10 هستش ولی چون هنوز برای تمام مدل یک Hc در نظر گرفتم این مساله بزرگی نیست. مشکل اینجاست که اصلا PMPATH برای مدل من Error میده ومنو حسابی کلافه کرده برای همین هم یک Email به دکتر Chiang زدم، امیدوارم که جوابم رو بده. اصلاً کلاً نمی دونم چیکار باید بکنم البته حدود هدایت هیدرولیکی من در حدود 5.5 m/day هستش ولی از قسمت 2D-Visulization که نگاه می کنم هدهای سطح آب درست بنظر نمی رسند و فکر کنم باید عدد رو در حد 0.055 در نظر بگیرم که با داده ها اصلاً جور نیست.
به همین دلیل مطمئن میشم که یکجا لنگه!.
از طرفی دیگه وقت زیادی باقی نمونده شاید بهتر باشه با خانم مهندس حسنی هم یک مشورتی بکنم.
با تشکر
واحددوست
Babak.Vaheddoost
سهشنبه، دی ۰۷، ۱۳۸۹
آخرین وضعیت فصول
با توجه به آخرین وضعیت موجود باید عرض کنم که بنده از امروز یعنی 1389/10/06 نهایتاً تا تاریخ 1389/11/30 زمان دارم تا دفاع از پایان نامه رو انجام بدهم ولی ظاهراً تا اواسط اسفند ماه نیز می توان از پایان نامه دفاع نمود لذا با فرض دفاع در اواخر بهمن ماه 54 روز وقت دارم تا دفاع کنم لذا زمانبندی ای که در نظر گرفته ام بصورت زیر است.
- 14 روز جهت اتمام مدل
- با توجه به وضعیت فصول گذشته که قبلاً نوشته شده، نوشتن فصول 6و7 نهایتاً 10 روز
- ویرایش فصول 5 روز
- ارئه نسخه به اساتید جهت ویرایش وتائید، 7 روز
- تهیه وساخت فایل دفاع PowerPoint ، 5 روز
- کارهای متفرقه 7 روز
- فصل اول: سدهای زیر زمینی
- فصل دوم: پیشینه تحقیق
- فصل سوم: معرفی نرم افزار Modflow
- فصل چهارم: معرفی محل
- فصل پنجم: مواد وروش تحقیق
- فصل ششم: توسعه مدل
- فصل هفتم: نتیجه گیری و پیشنهادات
خدا کمک کند!!
با تشکر
بابک واحددوست
Babak.V.Doost
جمعه، دی ۰۳، ۱۳۸۹
بازهم مشکل
دوستان از تیتر تعجب نکنید، آخه الان ساعت 1:10 am هست و بنده در نهایت یاس و ناامیدی دارم از نو برنامه رو اونهم تو شرایط Steady از اول Run می کنم.
بله دیگه اینجوریاست، پایان نامه اینش سخته دیگه !!!!!!!
به هر حال مساله اینه که من بعد از استخراج نفوذ و تبخیر وتعرق بعد از کالیبراسیون شهریور ماه 1383 رفتم دنبال مهرماه ولی اصلاً جوابها بی ربط بودن واسه همین با آقای دکتر رضایی و خانم مهندس حسنی مشورت کردم (بعد از کلی سوزاندن فسفر) آخر سر برگشتیم به اول و باید کارهای زیر رو انجام بدهیم (به نقل از خودم، دکتر رضایی و خانم مهندس حسنی):
- قبل از هرچیز باید برم (که الان رفتم) توی شرایط Steady تا جای ممکن مدل رو کالیبره کنم، البته بدون Recharge و Evapotraspiration تا پارامترهای ثابت کالیبره بشن.
- بعد باید Recharge وEvapotranspiration رو وارد کنم و باز کالیبره کنم واز حالت Steady برم روی حالت Transient.
- بعد باید چند Period تعریف کنم (متاسفانه نه تو سال 83 و نه 84 یکسری ماه کامل از یک سال رو نداریم.) تا ماههای مورد مطالعه رو وارد کنم و بعد باز کالیبره کنم، احتمالاً بعد اون باز با سال 84 هم کالیبراسیون انجام بدم.
- آخرسر برم رو پش بینی برای سالهای 85 و 86
راستی روی نرم افزار محاسبه تبخیر، راواناب و نفوذم چند تا تغییر گرافیکی دادم ولی متاسفانه نتونستم تو صفحه Download وبلاگ بزارمش واسه همین بجاش اینجا برچسب لینک رو هم میزنم و ازینجا می زارم واسه DOWNLOAD.
با تشکر
بابک واحددوست
Babak.Vaheddoost
دوشنبه، آذر ۲۹، ۱۳۸۹
نرم افزار محاسبه تبخیرو تعرق و نفوذ
بالاخره خودم آستینها رو زدم بالا و تحت Excel یک نرم افزار نوشتم که از روش Thornthwaite تبخیر و تعرق رو محاسبه کرده و با کمک روش SCS رواناب رو بدست می آره، بعد از اون هم نفوذ رو حساب میکنه. متاسفانه عکسشو نتونستم رو هاست بزنم پس مستقیم قرار می دهم رو وبلاگ در ضمن می تونین برین از قسمت دانلود وبلاگ نرم افزار رو که گزاشتمش تو rapidshare دانلود کنین. اگه هم زمانش گذشته بود به من یک Email بزنین براتون بفرستم.
با تشکر
بابک واحددوست
Babak.V.Dst
سهشنبه، آذر ۲۳، ۱۳۸۹
تصمیم گیری در مورد روش محاسبات تبخیر و تعرق و نفوذ
بالاخره پس از کلی کلنجار رفتن با مسال و افراد و با کمکهای بی دریغ جناب آقای دکتر رضایی و همچنین خانم مهندس حسنی به یکسری نتایج اساسی دست پیدا کردم و فکر می کنم دیگه کار تمومه و دیگه می تونم پایان نامه رو نهایی کنم. به این دلیل تو این قسمت در مورد داده های ایستگاه هام صحبت می کنم:
در این طرح که واقع در شمال استان آذربایجان غربی و تقریباً میان دو محدوده "33.09'49'44 East و "27.68'7'39 North هست، بنده (همانطور که قبلاً هم گفته بودم) رفتم سراغ اداره هواشناسی استان و از اونجا تو نستم داده های ایستگاههای چالدران و ماکو و قره ضیا الدین را بدست بیارم و بعداً آوردم و دیدم که متاسفانه در روش پولیگونهای تیسن منطقه بنده داخل محدوده تاثیر ایستگاهها قرار نمی گیره لذا پس از کلی دردسر رفتم سراغ آب منطقه ای استان آذربایجان غربی و از اونجا داده های 3 ایستگاه دیگه به نامهای جمالکندی، پلدشت و بازرگان رو بدست آوردم.
بعد از اینکه 6 ایستگاه هواشناسی را پیش هم گذاشتم دیدم که تمامی داده ها در آنها یکسان نیست مثلاً آمارهای موجود برای 3 ایستگاه اول شامل موارد زیر بود:
- مختصات جغرافیایی ایستگاه، ارتفاع ایستگاه از سطح دریا
- بارش متوسط ماهیانه
- درجه حرارت Max ماهانه
- درجه حرارت Min ماهانه
- درجه حرارت متوسط ماهانه
- درجه حرارت مکزیمم مطلق ماهانه
- درجه حرارت مینیمم مطلق ماهانه
- رطوبت نسبی Max ماهانه
- رطوبت نسبی Min ماهانه
- متوسط رطوبت نسبی ماهانه
- ساعات متوسط آفتابی در هر ماه
- میزان تبخیر پتانسیل ماهانه
- سمت و سرعت باد روزانه
- تعداد روزهای یخبندان
داده هاییکه از طریق سازمان آب منطقه ای گرفته بودم شامل موارد زیر بود:
- مختصات جغرافیایی ایستگاه، ارتفاع ایستگاه از سطح دریا
- درجه حرارت Max ماهانه
- درجه حرارت Min ماهانه
- درجه حرارت معدل ماهانه
- درجه حرارت حداکثر مطلق ماهانه
- درجه حرارت حداقل مطلق ماهانه
- بارش متوسط ماهانه
- میزان تبخیر از تشتک ماهانه
- نم نسبی 6:30 ساعته
- نم نسبی 12:30 ساعته
- نم نسبی 18:30 ساعته
- مجموع سرعت باد در ماه به Km
همانطور که مشاهده می شود این 6 ایستگاه در صورتیکه در نظر داشته باشیم که از داده های هر 6 تا از آنها استفاده کنیم دچار مشکلات زیادی خواهیم شد، چراکه علاوه بر اینکه از هر فرمولی برای تمامی آنها نمی توان استفاده کرد و از طرفی دیگر بعلت اینکه منطقه مورد نظر در میان ترکیبی از ایستگاههای اداره هواشناسی و سازمان آب منطقه ای قرار داده شده است باید از درون یابی و یا رسم منحنی های همباران، همدما و... استفاده کرد که کاری بود بس زمانگیر.
جهت اینکه هر یک از مطالب از حوصله خواننده خارج نباشد ادامه آنرا در مطلب بعدی می نویسم.
با تشکر
واحددوست
Babak.V.Dst
شنبه، آذر ۰۶، ۱۳۸۹
گزارش پيشرفت
خوشبختانه بالاخره اولین مدل کالیبره شده بنده با داده های نیمه فرضی تقریباً به اتمام رسید. (چه برای قبل و چه برای بعد از احداث)
حالا وقتشه که با پیدا کردن مقادیر دقیق Recharge و Evaporation شروع کنم به fine tuning و کالیبراسیون دقیق ولی حقیقت این است که در کالیبراسیون طرحهایی مانند پایان نامه بنده بعلت کوچک بودن ابعاد و نیز تغییرات در حد سانتیمتر و میلیمتر و با توجه به در نظر گرفتن ناهمواری و نا هماهنگی های موضعی و احتمال تاثیر پدیده موئینگی و... دقت زیاد میسر نخواهد بود و مشاهد واریانسهایی در حد دهم و صدم دستاوردی بزرگ به شمار می آید. چراکه در ابتدای امرعده زیادی معتقد بودند که این طرح اصلاً جواب نخواهد داد. از طرفی دیگر جواب ندادن داده ها علل منطقی دیگری هم دارد و آن اینکه در طرح اصلی نیز احتمال وجود شکاف و ترک در بستر یا دیواره مطرح شده چراکه جوابها نا متقارن می بوده اند.
با تشکر
بابک واحددوست
Babak.Vaheddoost
پنجشنبه، آذر ۰۴، ۱۳۸۹
موفقیت در کالیبراسیون شماره 1
با سلام.
بله دوستان، همونطور که از تیتر مشخص میشه بنده برای اولین مدل کالیبراسیون شده سد زیرزمینی صوفی مربوط به اطلاعات شهریور 81 با داده های رودخانه و زهکش (قنات) با واریانس 0.16 موفق به این کامیابی اولیه شدم اما در قسمت زهکش (قنات) و در در بخش Water budget که دبی قنات بنده در آن تاریخ حدود 13 lit/sec بود داده های مدل نزدیک 14.8 lit/sec را نشان میدهد که بنظر خودم موفقیت اولیه خوبی به حساب میاد. نکات اساسی در تغییر حالت از واریانسهای بالا به این مقدار موارد زیرهستند.
زمان را از حالت Steady به transient تبدیل کردم. با توجه به تغییر زمان به غیر ماندگار، ضریب ذخیره ای رو که توی طرح بدست آمده را در مدل دخیل نمودم. به کمک آقای دکتر رضایی مقدار Bulk Density را هم قرار دادیم چون بنده در اوایل درک دقیقی از Bulk Density نداشتم. با نظر اجمالی به مدل پایان نامه و همفکری با جناب خانم مهندس حسنی پکیج های رودخانه و زه کشم رو هم از سیستم cell by sell input method تبدیل کرد به سیستم poly line input method. چون تو این حالت کنترل بعضی از موارد مانند عرض رودخانه یا زهکش راحتتر بود چون سیستم گرید بندی بنده ثابت بود و آب باریکه جاری در کف مخزن بنده به هیچ وجه با عرض 5 متر مطعلق به مش بندی بنده همخوانی نداشت. محل سد را با کد 1- در نظر گرفتم که البته این مربوط به قسمت کالیبراسیون نمی شود.خوشبختانه نتایج جالبی هم گرفته شد؛ با کمک PMPATH مشاهده کردم که همانطوریکه در طرح هم به آن اشاره شده بود. قنات دیرسک که واقع در دل مخزن ما می بود بعلت عدم بررسی های اولیه و بالاتر بودن Elevation آن از روی تاج سد بیشتر دارد زه پس میدهد تا اینکه مخزن ما را زهکشی کند از طرفی دیگر رودخانه در بالادست آب را به سفره داده و در پایین دست آنرا زهکشی می نماید (بعلت شرایط توپوگرافی). و چون در ابتدا من کل مسیر اولیه رودخانه را در نظر گرفته بودم ، رودخانه از کنار پیزومتر شماره 3 بنده هم رد میشد درحالیکه در عمل همانطور که در گزارش بازدید از محل هم مشخص بود آب باریکه جاری را قبل از رسیدن به محل سد در بالادست سوار کانالی کرده و از محدوده طرح خارج می نمایند که این می تواند دلیل منطقی ای برای عدم تطبیق پیزومتر 3 که در پایین دست سد قرار دارد باشد چراکه حدود 90% از اینکه واریانس زیاد است مربوط به پیزومتر3 می باشد و در آن نقطه داده های زیادی در دست نیست و از طرفی دیگر بعد از مشاهده شرایط محیطی احنتمال آنرا می دهم که در محل تراوش آب یا جریان آب به طرف آن در قسمتهای پایین دست برقرار باشد یا که خود حوزه آبریز سد ما آب هدایت شده ارتفاعات خود را به سمت این پیزومتر بصورت موضعی هدایت کند. و با این منوال مشاهده گردید که در کف مخزن ما جریانهای عمومی بصورت سرگردان در جهتای مختلفی سیر کرده و گاهی در محلی (علیرغم شرایط توپو گرافی) انباشته می شوند و در ادامه پیشروی زمانی آن محل را رد می کنند.
حال در فاز بعدی جهت Fine tunning باید مقادیر Recharge و Evaporatian رو قرار بدم تا دقیق کالیبره بشه و بعد از اون باید برای ماههای سالهای 83 و 84 به تعداد 24 مدل کالیبراسیونی انجام بدهم (البته ضروری نیست و برای افزایش دقت و کنترل صحت مدل انجام می شود). و بعد برم سر سالهای 85 و 86 که نتایج احداث سد را مقایسه کنم. بعد از اون هم در نظر دارم چند تا مقاله کار کنم که احتمالاً اولیش برای یک کنفرانس داخلی باشه.
با تشکر و احترام.
بابک واحددوست
Babak.V.Dst
یکشنبه، آبان ۳۰، ۱۳۸۹
گزارش پيشرفت
الان نشستم دارم رو مدل کار می کنم. چند تا تصمیم گرفتم که دارم اعمال می کنم.
- دارم river و drain رو به کمک polyline کار می کنم. چون اینجوری خیلی منطقی تره تازه کنترل بیشتری روی داده ها بخصوصو تو قنات دارم چون موقعیت چاهها رو دقیقاً می دونم
- برای ساده کردن تخمین Initial Heads از surfer استفاده کردم.
- می خواهم تو قسمت Head Observation از منوی PEST قبل از اجرا کردنش برای چاهک مشاهداتی P2 بعلت نزدیکتر بودن به محل سد و قرار گرفتن در دل مخزن مقدار Weight بیشتری رو در نظر بگیرم.
با تشکر از توجهتون
بابک واحددوست
Babak.V.Dst
تغییرات
امروز سر ظهر حوالی ساعت 2 رفتم به دانشکده کشاورزی و در قسمت گروه آب، اونجا با جناب آقای دکتر رضایی و خانم مهندس حسنی جلسه ای داشتیم در مورد مدل بندی و پس از مشاهده مدل خانم مهدس کمی هم در مورد مدل بنده حرف زده شد که نتایج آن بدین شرح است.
- تعیین Kv برای کف رودخانه
- تعیین سد زیر زمینی با هد ثابت 1-
- کالیبره کردن مدل، یک بار با steady state ویک بار با transient
- کنترل Head و واریانسها از طریق PEST
با تشکر
واحددوست
Babak.V.Doost
چهارشنبه، آبان ۲۶، ۱۳۸۹
وضعیت مدل
وقتی مطلب قبلی رو می نوشتم در حقیقت چند تا مدل ساخته بودم و Run کرده بودم منتهی نمی خواستم همه چیز قاطی بشه، واسه همین الان آخرین وضعیت رو قرار می دهم تا ببینین.
- زه کش رو دوباره اضافه کرد. چون چاه دقیقاً اون چیزی رو که می خواستم جواب نداد. (خیلی هم جالب شد).
- PEST هم Run شد ولی دقیقاً نمی دونم کجاشو باید کنترل کنم، چون Discrepancy برار 0.01% هست!!!!!!
- سد رو با Horizatal flow barriers مدل کردم.
این حالت قبل از PMPATH هستش. همونطور که می بینی سیستم زه کش، رودخانه و سلولهای 1- مشخص هستند. چند تا المان هم تو جاهای مختلف در نظر گرفتم تا Particel tracking کنمشون
تو عکس زیر هم مشاهده می کنین، همه چیز گویا هستش. تازه همونجور که گفتم می بینی که زه کش عمل نمی کنه و فقط داره از خودش تراوش می ده که تو عمل هم مهندس جعفری اینا با این مشکل مواجه شده بودند. رودخانه هم تو بالا دست داره تراوش می ده ولی تو پایین دست داره آب رو زهکشی میکنه. خیلی جالبه نه؟ اگه تو مقطع کناری هم ببینین میشه دید که آب میاد از روی تاج سد سرریز میشه!!!!!!
با تشکر / واحددوست/ B.V.Dst
'گزارش تغییرات و پیشرفت
از آخرین دفعه ای که مطلب ننوشتم کلی اتفاق خوب و بد افتاده. مثلاً وقت گذاشتم و رفتم دانشکده کشاورزی دانشگاه ارومیه و چند تا پایان نامه در مورد Modflow خوندم. بعد آقای دکتر رضایی بنده رو معرفی کردند برم با آقای مهنس خالق زاده که پایان نامه ایشون در مورد Modflow بود صحبت بکنم. متاسفانه بعلت اینکه زمان زیادی از فارغ التحصیلص ایشون گذشته بود نتونستند کمکی بکنند. (البته لازم به ذکر هستش که ایشون بسیار مودب و با نزاکت بنده را تحویل گرفتند). بعد از اون زنگ زدم با آقای دکتر دادمهر هم کمی صلاح و مشورت نمودم. و بعد از این همه ورجه وورجه فعلاً تصمیم دارم کارهای زیر رو انجام بدم:
- برای تست ورودی و خروجی رودخانه رو کلاً هد ثابت بدم (1-) تا گرادیان هیدرولیکی برقرار بشه (دکتر دادمهر صلاح نمی دونستند).
- رودخانه کاملاً محصور در داخل مدل میشه
- زه کش رو حذف کنم چراکه در عمل نیز بعلت بالاتر بودن سطح آن از سطح تاج زیرزمینی، زهکش بازده خود را از دست داده است و بعد به جای همون یک چاه با تخلیه معادل دبی زه کش عمل کنم.
- تست کنم ببینم که خود سد را با چی مدل کنم؟ با سلولهای 0 یا با Horizantal flow barriers؟
- به پیشنهاد آقای دکتر دادمهر ابعاد سلولهامو بزرگتر انتخاب کنم؟؟؟
به هر حال کلاً تو استرس هستم چون حد اکثر تا اواسط بهمن باید دفاع کنم ولی هنوز اولین مدل هم درست و حساب run نشده.
برام دعا کنید.
با تشکر
واحددوست.
جمعه، آبان ۱۴، ۱۳۸۹
فصول پایان نامه
الان داشتم مقاله ای رو که گفته بودم میخوندم همون
Methods and Guide lines for Effective Model Calibration
خیلی خفنه کلاً منطق ریاضی Calibration هستش، هر کی علاقه مند به مدل ریاضی مطلب هستش توصیه بنده اینه که این مقاله رو بخونه تقریباً 90 صفحه هستش.
بعد از خوندن چند صفحه از این مقاله سنگین، اومدم تا کم کم در مورد فصول پایان نامه تصمیم بگیرم، فعلاً چیزی که تو ذهنم هستش به قراره زیره:
- فصل اول : سدهای زیر زمینی ( تعریف، انواع، کاربرد، مزایا، معایب ، تاریخچه و ...).
- فصل دوم: پیشینه تحقیق و کارهای انجام شده و همچنین ذکر پروژه ها.
- فصل سوم: لزوم تحقیق و پیش بینی های اولیه.
- فصل چهارم: معرفی طرح، پارامترهای اولیه، بررسی های اولیه، محل و موقعیت.
- فصل پنجم: روش تحقیق، معرفی Modflow و PMWIN و گزینه های بررسی.
- فصل ششم: محاسبات، چالشهای تحقیق، چارتها، دیاگرامها و روشهای ابتکاری.
- فصل هفتم: نتایج، بررسی ها و مقایسه مدل با برداشتها.
- فصل هشتم: پیشنهادات و ژرف نگری.
با تشکر
بابک واحددوست
Babak.V.Dst
شنبه، آبان ۰۸، ۱۳۸۹
تغییرات مدل
ببخشید یادم رفت که بگم آخرین تغییرات فیزیکی مدل به قرار شکل زیره:
منطقش اینه که شیب هیدرولیکی رو باید یکجورایی بین دو سر رودخانه ایجاد کرد، سلولهای خالی هر طرف هم مربوط به محل رودخانه هستش چون اگر کد 1- بگیره دیگه نمیشه سلول رودخانه باشه چون نرم افزار خود رودخانه رو بعنوان Time variant میشناسه ولی 1- یعنی هد ثابت. البته به غیر اینها یکسری تغییراتی رو هم به TOL, BOL دادم که بعداً سعی میکنم راجع به اونها چیزایی بنویسم.
با تشکر
واحددوست
Babak.V.Dst
چهارشنبه، مهر ۲۱، ۱۳۸۹
گزارش پيشرفت كار
بالاخره همچي رو وارد كردم و اولين بار Runكردم، با كلي Error مواجه شدم. خيلي دلخور شدم چون نمي دونستم مشكل چيه و رفتم با آقاي دكتر رضايي تماس گرفتم. بعد از اون آقاي دكتر شماره خانم مهندس حسني رو به بنده داد كه بتونم از ايشون كمك بخوام. واسه همين زنگ زدم و از خانم مهندس حسني كمي كمك خواستم. كه علي الحق كمكم كرد و چند تا پيشنهاد اوليه داد. با اينكه ايشون از جزئيات پايان نامه بنده خبر نداشت، چيزهاي خوبي براي شروع گيرم اومد.
واسه همين چند تا كار كردم:
- بجاي modflow از مدل PEST استفاده مي كنم، همونطور كه دكتر كهربايي هم گفته بود
- كلاً زه كش و Evaporationرو حذف كردم تا مرحله به مرحله وارد كنم
- تنها از داده هاي River package و چاههاي مشاهداتي استفاده كردم
- در ضمن مدل رو هم تغيير دادم (براساس سلولهاي 1-)ه بعداً دلايلش رو عرض ميكنم.
واحددوست
Babak.V.Doost
چهارشنبه، مهر ۱۴، ۱۳۸۹
گزارش پيشرفت
بالاخره بعد اينكه منوي Gridبه اتمام رسيد. رفتم سر وقت منوي Parameters اونجا مجبور شدم كه باز ضخامت لايه ها رو كاليبره كنم، آخه يكم اختلاف تو رقوم سطح آب با رقوم كف لايه ها بود كه برطرف شد. بعد از اون زمان رو براي شهريور ماه 1383 معادل 30 روز در نظر گرفتم البته به خاطر 30 روز ماه بطور كلي البته ممكنه تغييرش بدم بكنمش 31 روزه. بعد از اون چاهكهاي مشاهداتي رو مشخص كردم،
حالا بايد 3 تا كار اساسي بكنم.
- مشخص كردن Horizatal Hydraulic Conductivity
- Effective Porosity
- ببينم كه آيا لازمه Bulk Densityرو هم بدم يا نه (احتمالاً لازم نيست)
با تشكر
واحددوست
Babak.V.Doost
جمعه، شهریور ۲۶، ۱۳۸۹
مش بندی محیط
بعد از جلسه هماهنگی با آقای دکتر شروع به ساخت وتکمیل جدول مشخصات سلولها کرده، یک محیط با ابعاد 250 متر در 250 متر مشخصاً در محدوده طرح انتخاب شد. سلولها با ابعاد 5 متر در 5 متر تقسیم بندی شدند ( 50 ستون، 50 ردیف).
با توجه به عمق لایه در سلول ( مقطع ژئوالکتریک)و تپوگرافی هر سلول، شروع به مشخص کردن عکس العمل هر سلول،TOL و BOL آن نموده شد تا عملاً از مرحله منوی Grid گذر کنیم.
در مورد مشخصات و نوع لایه مشخصاً نوع Unconfined انتخاب شد، چون تغیر ارتفاع سفره در نظر ما قرار ندارد و تغییرات فشار در اینجا کنترل نمی شود.
با تشکر - واحددوست