حل مشکلات

سلام دوستان.
بالاخره پس از یک وقفه طولانی، اگرچه از برنامه کمی عقب موندم ولی باید عرض کنم که خیلی از مشکلات حل شد. که گزارش اون به شرح زیره.

1. در ابتدا با کمک نمودار نسبت به هر یک از پیزومترها و کل مدل ساخته شده حساسیت سنجی انجام گرفته شد. و با کمک آنها طرح اولیه محدوده هدایت هیدرولیکی تعریف گردید.
   
2. تعیین دقیق محدوده بندی Hc انجام گرفت و مدل در حالت Steady تا حد واریانس 0.65 کالیبره شد.
3. از آنجائیکه رواداریهای مجاز برای کنترل رقوم سطح آب برای پیزومترهای موجود در حد 5% توسط دکتر رضایی اعلام شده بود در مورد پیزومتر شماره 2 کمی مشکل وجود داشت که آخر سر با تخفیف حدود 0.001 مورد قبول واقع شد.
4. با کمک حساسیت سنجی مدل برای ضریب ذخیره، در ابتدا به کمک Hc و فرمول مربوطه مقادیر اولیه Sy تخمین زده شد.
5. بعد از تخمین اولیه Sy، مطابق Hc نمودارهای حساسیت سنجی برای هر یک از پیزومترها و کل مدل رسم گردید.
6. تعیین دقیق محدوده های ضریب ذخیره با آزمون و خطا امتحان گردید که در نهایت مدل با واریانس 0.061 مدل شد و بعلت پایین بودن واریانس و محدوده های رواداری مدل مورد قبول واقع شد.
7. زمانهای محدوده طرح با توجه به وفور دیتا های موجود و داده سازی از شهریور 83 لغایت مرداد 83 به مدت 365 روز جهت دوره آموزش مدل و کالیبراسیون در نظر گرفته شد.
8. دیتاهای مربوط به Recharge و Evapotranspiration با کمک نرم افزار ساخته شده محاسبه گردید و در زمانهای مربوطه وارد مدل گردید.
9. تنها مشکل فعلی آن است که مدل زمانها را نمی شناسد و بعد از Run فقط اولین Period رو نمایش می دهد. که البته قبلاً با این مساله برخورد کرده بودم و گمان می کنم که به راحتی قابل حل است.

در نهایت فازهای بعدی جهت انجام در روزهای آینده به شرح زیر خواهند بود.

1. اجرای مدل و کالیبراسیون دقیقتر در صورت نیاز
2. ایجاد مدل Key جهت حفاظت از اطلاعات
3. احداث سد در محل مورد نظر
4. فعال کردن Polyline مربوط به زهکش سد و ایجاد ارتباط با قنات امیر کهریزی
5. اجرای مدل و کنترل Water budget
6. نوشتن فصل 6 و 7 باقیمانده.
7. اصلاح و تحریر نهایی پایان نامه و ....

دوستان مثل اینکه دیگه آروم آروم قصعه این پایان نامه هم داره به سر می رسه. حالا خدا کریمه ببینیم چی میشه.

با تشکر واحترام
واحددوست

وضعیت مدل

سلام.
این عکس رو اینجا قرار میدم تا اگه دکتر Chiang بخواهد به وبلاگ سر بزنه ببینه.



همونطور که دیده میشه خط آبی رودخانه منه، خطهای زرد زهکشها DRAIN هستن و خطوط پیزومتریک دیده میشن اون نقاطی هم که آبی کم رنگه در حقیقت اون قسمت تز زهکش هستش که بعد از اجرای سد ساخته خواهد شد والان غیر فعاله.

با تشکر
واحددوست.
Baba;.V.doost

فقط مشکل

سلام دوستان عزیز.

امروز رفتم خدمت جناب آقای دکتر رضایی و کمی در مورد پایان نامه حرف زدیم قرار شد برای Hc ضریب حساسیت در بیارم که بهترین باشه ولی متاسفانه مدل اصلاً درست جواب نمیده، البته باید بگم که variance در حد زیر عدد 10 هستش ولی چون هنوز برای تمام مدل یک Hc در نظر گرفتم این مساله بزرگی نیست. مشکل اینجاست که اصلا PMPATH برای مدل من Error میده ومنو حسابی کلافه کرده برای همین هم یک Email به دکتر Chiang زدم، امیدوارم که جوابم رو بده. اصلاً کلاً نمی دونم چیکار باید بکنم البته حدود هدایت هیدرولیکی من در حدود 5.5 m/day هستش ولی از قسمت 2D-Visulization که نگاه می کنم هدهای سطح آب درست بنظر نمی رسند و فکر کنم باید عدد رو در حد 0.055 در نظر بگیرم که با داده ها اصلاً جور نیست.
به همین دلیل مطمئن میشم که یکجا لنگه!.
از طرفی دیگه وقت زیادی باقی نمونده شاید بهتر باشه با خانم مهندس حسنی هم یک مشورتی بکنم.

با تشکر
واحددوست
Babak.Vaheddoost

Horizantal Hydraulic Conductivity

سلام.
فكر كنم منشا اصلي مشكل رو پيدا كردم. توي داده هايي كه من گرفتم با استفاده از لوگ چاه مقدار هدايت هيدروليكي رو برابر 0.0005 در نظر گرفته كه واحدش m/sec هستش ، اگه بخواهيم اونو تبديل كنيم به m/day مقدارش ميشه 43.2كه منطقي نيست چون مثل اينه كه از يك واحد مورد بررسي بستر ما روزانه 42 متر آب توليد ميشه!!!.
تازه براي روش دانه بندي وضع خرابتره، چون براي GW‌مقدار 7.1E-2 و براي GM مقدار 7.8E-4 رو در نظر گرفته كه واحدشون همون m/sec هستش و اگه بخواهيم ببريم به واحد m/day ميشنه 6134 و 67.4 كه واحدشون m/day هستش و مطمئناً‌ اين ضرايب مشكل دارند. در مورد GW كه شكي وجود نداره و لي در مورد GM و لوگ چاه احتمالاً واحدشون cm/sec هستش.
از طرفي برنامه هم اينو نشون ميده وبراي اين ارقام Error ميده، توي مثالهاش هم رنج براي چند نمونه خاك مدل ساززي براي اشل m/day در حدود يك رقمي يا بعد از اعشاره.
امروز عصر مي برم ، 5 فصلي رو كه نوشتم دكتر رضايي ببينه و در اين مورد هم صحبت مي كنم باهاش.

با تشكر
واحددوست
Babak.V.Dst

p2

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

با تشکر
واحددوست
Babak.V.Dst

یک قدم به جلو

سلام،
بالاخره در شرایط Steady مدل بدون حضور Recharge و Evapotranspiration برای شهریور 1383به واریانس 4.39 رسید. حالا دیگر نمی توان با تغییر پارامترهای رودخانه ویا زهکش مدل نسب به کوچکتر کردن واریانس اقدام کرد چراکه هر کدام از آنها مقدار قابل توجهی از هدهای مدل را جابجا می کنند و اصولاً جهت Fine Tunning باید از تغییر داده های با دامنه تاثیر زیاد اجتنای کرد و رفت سراغ پارامتر هدایت هیدرولیکی و با تعریف محدوده های مختلف نسبت به کالیبراسیون دقیقتر اقدام کرد.
یکی دیگر از مشکلات موجود در سر راه مدل آن است که قنات دوم طرح که در پایین دست سد زیر زمینی قرار دارد بعد از اجرای سد با کمک یک زهکش سرتاسری تغذیه می شود که این خود شرایط زهکش را تغییر می دهد، لذا باید کل مسیر را از ابتدا تعریف کرد اما سلولهای زهکش بالا دست سد را در Period های قبلی غیر فعال نمود.

با تشکر
واحددوست
VDst

بازهم مشکل

سلام،
دوستان از تیتر تعجب نکنید، آخه الان ساعت 1:10 am هست و بنده در نهایت یاس و ناامیدی دارم از نو برنامه رو اونهم تو شرایط Steady از اول Run می کنم.
بله دیگه اینجوریاست، پایان نامه اینش سخته دیگه !!!!!!!
به هر حال مساله اینه که من بعد از استخراج نفوذ و تبخیر وتعرق بعد از کالیبراسیون شهریور ماه 1383 رفتم دنبال مهرماه ولی اصلاً جوابها بی ربط بودن واسه همین با آقای دکتر رضایی و خانم مهندس حسنی مشورت کردم (بعد از کلی سوزاندن فسفر) آخر سر برگشتیم به اول و باید کارهای زیر رو انجام بدهیم (به نقل از خودم، دکتر رضایی و خانم مهندس حسنی):

1. قبل از هرچیز باید برم (که الان رفتم) توی شرایط Steady تا جای ممکن مدل رو کالیبره کنم، البته بدون Recharge و Evapotraspiration تا پارامترهای ثابت کالیبره بشن.
   
2. بعد باید Recharge وEvapotranspiration رو وارد کنم و باز کالیبره کنم واز حالت Steady برم روی حالت Transient.
3. بعد باید چند Period تعریف کنم (متاسفانه نه تو سال 83 و نه 84 یکسری ماه کامل از یک سال رو نداریم.) تا ماههای مورد مطالعه رو وارد کنم و بعد باز کالیبره کنم، احتمالاً بعد اون باز با سال 84 هم کالیبراسیون انجام بدم.
4. آخرسر برم رو پش بینی برای سالهای 85 و 86

البته گفتنش راحته، بنده فعلاً که تو شماره 1 گیر هستم. راستی معلوم نیست که در مراحل پیشرفت چه مشکلاتی پیش بیاد، خدا رفتگان این خانم مهندس رو بیامرزه واقعاً چند جا کمک درست و حسابی برام کرده. به هر حال مسولیت اصلی به عهده خودمه، متاسفانه بعلت کمی وقت (آموزش می گه تا آخر بهمن ماه، یهنی 2 ماه تمام) می ترسم نتونم تموم کنم. البته تدابیری هم اندیشیدم از جمله اینکه وقتی برنامه پیشرفت نداره شروع می کنم فصلها رو می نویسم، حالا کی میخواد بره مقاله هم کار کنه!!!! انشا الله که بتونم برسم به هر حال خواستم درد دلی کرده باشم. هرکی موضوعی شبیه این کار می کنه می دونه که بنده چی میگم. این مدل واقعاً چیز حساسیه به چند دهم یا صدم یکدفعه می بینی 2 تا از پیزومترها خشک شدن و....
راستی روی نرم افزار محاسبه تبخیر، راواناب و نفوذم چند تا تغییر گرافیکی دادم ولی متاسفانه نتونستم تو صفحه Download وبلاگ بزارمش واسه همین بجاش اینجا برچسب لینک رو هم میزنم و ازینجا می زارم واسه DOWNLOAD.

با تشکر
بابک واحددوست
Babak.Vaheddoost

همچنان مشکل تبخیر و تعرق

سلام،
متاسفانه هنوز مشکل داده های ورودی برای Recharge و Evapotranspiriation، وجود دارد. چونکه داده هاییکه از نرم افزار خانم مهندس حسنی برای Recharge بدست آمده مساوی صفر شد و بنظر جناب آقای دکتر رضایی این احتمال بسیار کم می باشد. با این وجود در مذاکرات تلفنی با خانم مهندس حسنی از احتمال وجود خشکسالی سراسری در حین سالهای 82 الی 87 صحبت به میان آمد که احتمالی قریب به یقین است. در نتیجه دست بکار شدم و خودم نرم افزار محاسبه تبخیر و تعرق و نفوذ نوشتم که از روش Thorenthwaite عمل می کند.
بنظر می رسد که جوابهای بدست آمده منطقی هستند. با این وجو قرار است که نتایج بدست آمده از دو نرم افزار تحت Excel رو به دکتر رضایی نشان بدهم.

با تشکر
واحددوست
Babak.V.dst

تصمیم گیری در مورد روش محاسبات تبخیر و تعرق و نفوذ

سلام،
بالاخره پس از کلی کلنجار رفتن با مسال و افراد و با کمکهای بی دریغ جناب آقای دکتر رضایی و همچنین خانم مهندس حسنی به یکسری نتایج اساسی دست پیدا کردم و فکر می کنم دیگه کار تمومه و دیگه می تونم پایان نامه رو نهایی کنم. به این دلیل تو این قسمت در مورد داده های ایستگاه هام صحبت می کنم:
در این طرح که واقع در شمال استان آذربایجان غربی و تقریباً میان دو محدوده "33.09'49'44 East و "27.68'7'39 North هست، بنده (همانطور که قبلاً هم گفته بودم) رفتم سراغ اداره هواشناسی استان و از اونجا تو نستم داده های ایستگاههای چالدران و ماکو و قره ضیا الدین را بدست بیارم و بعداً آوردم و دیدم که متاسفانه در روش پولیگونهای تیسن منطقه بنده داخل محدوده تاثیر ایستگاهها قرار نمی گیره لذا پس از کلی دردسر رفتم سراغ آب منطقه ای استان آذربایجان غربی و از اونجا داده های 3 ایستگاه دیگه به نامهای جمالکندی، پلدشت و بازرگان رو بدست آوردم.
بعد از اینکه 6 ایستگاه هواشناسی را پیش هم گذاشتم دیدم که تمامی داده ها در آنها یکسان نیست مثلاً آمارهای موجود برای 3 ایستگاه اول شامل موارد زیر بود:

1. مختصات جغرافیایی ایستگاه، ارتفاع ایستگاه از سطح دریا
2. بارش متوسط ماهیانه
3. درجه حرارت Max ماهانه
4. درجه حرارت Min ماهانه
5. درجه حرارت متوسط ماهانه
6. درجه حرارت مکزیمم مطلق ماهانه
7. درجه حرارت مینیمم مطلق ماهانه
8. رطوبت نسبی Max ماهانه
9. رطوبت نسبی Min ماهانه
10. متوسط رطوبت نسبی ماهانه
11. ساعات متوسط آفتابی در هر ماه
12. میزان تبخیر پتانسیل ماهانه
13. سمت و سرعت باد روزانه
14. تعداد روزهای یخبندان

که آمار های موجود در فایل ارائه شده به بنده در محدوده سالهای 82 تا 89 برای ایستگاه چالدران و سالهای 64 تا 88 برای ایستگاه ماکو و سالهای 47 تا 89 برای ایستگاه قره ضیاالدین می بود که بخوبی سالهای مورد بررسی طرح را یعنی 83 تا 86 را پوشش می داد.
داده هاییکه از طریق سازمان آب منطقه ای گرفته بودم شامل موارد زیر بود:

1. مختصات جغرافیایی ایستگاه، ارتفاع ایستگاه از سطح دریا
2. درجه حرارت Max ماهانه
3. درجه حرارت Min ماهانه
4. درجه حرارت معدل ماهانه
5. درجه حرارت حداکثر مطلق ماهانه
6. درجه حرارت حداقل مطلق ماهانه
7. بارش متوسط ماهانه
8. میزان تبخیر از تشتک ماهانه
9. نم نسبی 6:30 ساعته
10. نم نسبی 12:30 ساعته
11. نم نسبی 18:30 ساعته
12. مجموع سرعت باد در ماه به Km

که محدوده داده های این ایستگاهها در سالهای 78 تا 86 برای حوزه ایستگاه تازه کند بازرگان و سالهای 45 تا 86 برای ایستگاه پلدشت و سالهای 81 تا 86 برای ایستگاه جمالکندی بود.
همانطور که مشاهده می شود این 6 ایستگاه در صورتیکه در نظر داشته باشیم که از داده های هر 6 تا از آنها استفاده کنیم دچار مشکلات زیادی خواهیم شد، چراکه علاوه بر اینکه از هر فرمولی برای تمامی آنها نمی توان استفاده کرد و از طرفی دیگر بعلت اینکه منطقه مورد نظر در میان ترکیبی از ایستگاههای اداره هواشناسی و سازمان آب منطقه ای قرار داده شده است باید از درون یابی و یا رسم منحنی های همباران، همدما و... استفاده کرد که کاری بود بس زمانگیر.
جهت اینکه هر یک از مطالب از حوصله خواننده خارج نباشد ادامه آنرا در مطلب بعدی می نویسم.

با تشکر
واحددوست
Babak.V.Dst

تبخیر و بارش

سلام.
امروز با شماره ای که آقای دکتر داده بودند تماس گرفتم و با جناب آقای مهندس مهدیزاده (دانشجوی کارشناسی ارشد آبیاری) در مورد تبخیر و تعرق صحبت کردیم. ایشون پایان نامشون در مورد تبخیر و تعرق بود ایشون می فرمودند که:

1. تبخیر عبارتست از از دست دادن آب از سطح محلهایی مانند رودخانه، دریاچه و خاک.
2. تعرق عبارتست از از دست دادن آب از سطح بدنه گیاهان

از طرفی دیگر بعلت اینکه ، تعیین دقیق این عدد ممکن نیست لذا با پارامتر جدیدی به نام تبخیر و تعرق مواجه هستیم.
بنا بر عقیده ایشان چونکه داده های بنده بصورت متوسط ماهانه هستند باید از روش تبخیر و تعرق بلانی-کریدل استفاده کرد و باز هم بنا بر نظر ایشان در این فرمول به P (درصد ساعات آفتابی) و T (متوسط درجه حرارت و طول وعرض جغرافیایی احتیاج خواهد بود و از طرفی دیگر بعلت وجود ضرایب تصحیح در مورد استفاده از این فرمول شک داشتند.
در مورد Recharge هم فرمودند که با آقای دکتر رضایی همنظر هستند و استفاده از فرمول C.I.A=Q را پیشنهاد می کنند اما همانطور که مشخص میشود این فرمول دارای پارمتر I بوده و بعلت اینک شدت بارش (بعلت نبود زمان بارش) در دست نیست، لذا دقیقاً نمی توان از آن استفاده کرد. لذا قرار شد امروز برای بنده فایل Excel جهت محاسبه تبخیر و تعرق از روش بلانی - کریدل را بفرستند.
بعد از ایشون با آقای دکتر رضایی تماس گرفتم که متاسفانه موفق به صحبت با ایشون نشدم. پس به ناچار دوباره دست به دامان خانم مهندس حسنی شدم (ایشون در مورد مدل سازی یک دشت وسیع کار می کنند). ایشون می فرمودند که مقدار Recharge رو از روی 700 چاه پمپاژ منطقه و 3% مقدار نفوذ آن بدست آورده اند و از طرفی دیگر روش بلانی - کریدل را جهت استفاده تایید کردند. و قرار شد که ایشون هم فایل Excel در مورد محاسبه نفوذ و تبخیر و تعرق به من بفرستند. از طرفی دیگر بنده سوالی داشتم که کسی جواب نمی داد در مورد سطح تاثیر داده های هواشناسی. یعنی اینکه مثلاً اگر داده ما برای ماه شهریور 10 میلیمتر باشد، این داده بر چه سطحی از منطقه تاثیر گزار است؟ چراکه باید مقدار Recharge و Evaporation را برای هر سلول مجزا یا بصورت منطقه ای وارد کرد تا با ضرب ارتفاع در سطح مقدار حجم (بیلان) بدست آید. که بنظر ایشون نرم افزار ایشون عدد حاصله را از نظر سطح بی بعد خواهد کرد ولذا قابل استفاده بصورت عمومی می باشد.
فعلاً که منتظر هستم ببینم چطور میشه. حتماً نتایج را با شما در میون خواهم گزاشت.

با تشکر
واحددوست
Bababk.VDst

موفقیت در کالیبراسیون شماره 1

با سلام.

بله دوستان، همونطور که از تیتر مشخص میشه بنده برای اولین مدل کالیبراسیون شده سد زیرزمینی صوفی مربوط به اطلاعات شهریور 81 با داده های رودخانه و زهکش (قنات) با واریانس 0.16 موفق به این کامیابی اولیه شدم اما در قسمت زهکش (قنات) و در در بخش Water budget که دبی قنات بنده در آن تاریخ حدود 13 lit/sec بود داده های مدل نزدیک 14.8 lit/sec را نشان میدهد که بنظر خودم موفقیت اولیه خوبی به حساب میاد. نکات اساسی در تغییر حالت از واریانسهای بالا به این مقدار موارد زیرهستند.

1. زمان را از حالت Steady به transient تبدیل کردم.
2. با توجه به تغییر زمان به غیر ماندگار، ضریب ذخیره ای رو که توی طرح بدست آمده را در مدل دخیل نمودم.
3. به کمک آقای دکتر رضایی مقدار Bulk Density را هم قرار دادیم چون بنده در اوایل درک دقیقی از Bulk Density نداشتم.
4. با نظر اجمالی به مدل پایان نامه و همفکری با جناب خانم مهندس حسنی پکیج های رودخانه و زه کشم رو هم از سیستم cell by sell input method تبدیل کرد به سیستم poly line input method. چون تو این حالت کنترل بعضی از موارد مانند عرض رودخانه یا زهکش راحتتر بود چون سیستم گرید بندی بنده ثابت بود و آب باریکه جاری در کف مخزن بنده به هیچ وجه با عرض 5 متر مطعلق به مش بندی بنده همخوانی نداشت.
5. محل سد را با کد 1- در نظر گرفتم که البته این مربوط به قسمت کالیبراسیون نمی شود.

خوشبختانه نتایج جالبی هم گرفته شد؛ با کمک PMPATH مشاهده کردم که همانطوریکه در طرح هم به آن اشاره شده بود. قنات دیرسک که واقع در دل مخزن ما می بود بعلت عدم بررسی های اولیه و بالاتر بودن Elevation آن از روی تاج سد بیشتر دارد زه پس میدهد تا اینکه مخزن ما را زهکشی کند از طرفی دیگر رودخانه در بالادست آب را به سفره داده و در پایین دست آنرا زهکشی می نماید (بعلت شرایط توپوگرافی). و چون در ابتدا من کل مسیر اولیه رودخانه را در نظر گرفته بودم ، رودخانه از کنار پیزومتر شماره 3 بنده هم رد میشد درحالیکه در عمل همانطور که در گزارش بازدید از محل هم مشخص بود آب باریکه جاری را قبل از رسیدن به محل سد در بالادست سوار کانالی کرده و از محدوده طرح خارج می نمایند که این می تواند دلیل منطقی ای برای عدم تطبیق پیزومتر 3 که در پایین دست سد قرار دارد باشد چراکه حدود 90% از اینکه واریانس زیاد است مربوط به پیزومتر3 می باشد و در آن نقطه داده های زیادی در دست نیست و از طرفی دیگر بعد از مشاهده شرایط محیطی احنتمال آنرا می دهم که در محل تراوش آب یا جریان آب به طرف آن در قسمتهای پایین دست برقرار باشد یا که خود حوزه آبریز سد ما آب هدایت شده ارتفاعات خود را به سمت این پیزومتر بصورت موضعی هدایت کند. و با این منوال مشاهده گردید که در کف مخزن ما جریانهای عمومی بصورت سرگردان در جهتای مختلفی سیر کرده و گاهی در محلی (علیرغم شرایط توپو گرافی) انباشته می شوند و در ادامه پیشروی زمانی آن محل را رد می کنند.

حال در فاز بعدی جهت Fine tunning باید مقادیر Recharge و Evaporatian رو قرار بدم تا دقیق کالیبره بشه و بعد از اون باید برای ماههای سالهای 83 و 84 به تعداد 24 مدل کالیبراسیونی انجام بدهم (البته ضروری نیست و برای افزایش دقت و کنترل صحت مدل انجام می شود). و بعد برم سر سالهای 85 و 86 که نتایج احداث سد را مقایسه کنم. بعد از اون هم در نظر دارم چند تا مقاله کار کنم که احتمالاً اولیش برای یک کنفرانس داخلی باشه.

با تشکر و احترام.

بابک واحددوست

Babak.V.Dst

تغییرات

سلام.
امروز سر ظهر حوالی ساعت 2 رفتم به دانشکده کشاورزی و در قسمت گروه آب، اونجا با جناب آقای دکتر رضایی و خانم مهندس حسنی جلسه ای داشتیم در مورد مدل بندی و پس از مشاهده مدل خانم مهدس کمی هم در مورد مدل بنده حرف زده شد که نتایج آن بدین شرح است.

1. تعیین Kv برای کف رودخانه
2. تعیین سد زیر زمینی با هد ثابت 1-
3. کالیبره کردن مدل، یک بار با steady state ویک بار با transient
4. کنترل Head و واریانسها از طریق PEST

در ضمن بقیه مطالب و توضیحات را در ادامه مطلب به خدمتتون عرض خواهم کرد.
با تشکر
واحددوست
Babak.V.Doost

و همچنان مشکل Calibration

سلام دوستان.
میدونم این مطلب خیلی تکراری شده ولی الان بعد از اینکه دوباره تغییراتی دادم و PEST رو Run کردم باز از یک سد گذشتیم و به سد دیگری برخورد کردیم.
ایندفعه اومدم برای تمام سلولهای کد 1 مقدار initial head به مقدار 10 سانتیمتر تعریف کردم، بعد از اون رفتم سراغ pest و تو اونجا گزینه های (مثل شکل زیر) رو انتخاب کردم.

و بعد از اون مدل رو Run کردم ولی بازهم ERROR!!!! ببینید، البته دفعه دوم بهتر شد یعنی کار کرد که این قسمت مال Run record file هستش:

حالا باید برم تو نخش، از هیچی بهتره بالاخره یک قدم به جلو هستش!!!

باتشکر
واحددوست
Babak.V.Dst

یک مشکل اساسی

سلام دوستان،
علیرغم اینکه چند بار مدل را Run کردم هنوز نتونستم اولین کالیبراسیون را اعمال کنم و این اصلاً خوب نیست. فکر میکنم احتیاج به کمک حرفه ای دارم.
متاسفانه نرم افزار قلق خاص خودش را داره و من هم با اینکه هی دارم Toturial رو می خونم و کلی مطلب گیرم اومده، هنوز نتونستم درست و حساب Run کنم. الان هم دارم دنبال کسی می گردم تا یکسری از مهارتها و قلقهاشو یادم بده. اگه پیشنهادی دارین خوشحال میشم که به من ایمیل کنین.
با تشکر
واحددوست.

آماده سازی اولیه محیط

سلام دوستان،
مدتها تو این فکر بودم که فقط قسمت مخزن رو به همراه قسمتهای لازم مدل کنم چون ایجاد شیب طولی و عرضی آبراهه دردسر زیادی درست میکنه، آخه تاثیر خود شیب روی نتیجه فکر کنم زیاد باشه.
ولی متاسفانه یکی از چاههای مشاهداتی من بالاتر از سطح مخزن قرار گرفته که خیلی بده واسه همین تصمیم گرفته بودم که یک حجم 3 بعدی برای تراز بالاترین تا پایین تریت سطحم در نظر بگیرم بعد با پیش بینی لایه صفر در مقاطع عرضی و لایه های زه کش در طول مقطع، شیب منطقه رو مدل کنم تا یکی از نقاط مشاهداتیم کم نشه تا بعداً بتونم منطقه رو کالیبره کنم.
راستش زیاد مطمئن نبودم دارم چیکار میکنم ولی الان داشتم یک PDF رومطالعه میکردم که یک عکس از برش عرضی یک مدل نظرم رو جلب کرد، دیدم می تونم با دستورات Top of layer و Bot of layer لایه های با مقاطع متغیر ایجاد کنم، این به نظر خوب میاد چون همونطور که قبلاً هم گفته بودم یک مشکل دیگه همون لایه نا همگون ماسه داخل بستره که در جناح چپ مخزن در داخل بستر وجود داره که یکم معادلاتم رو به همزده بود فکر کنم الان با انتخاب مش بندی متغیر و استفاده از ارتفاعهای متغیر بتونم هم مساله شیب طولی و عرضی بستر رو حل کنم و هم مساله لایه ها رو البته باید دید که چی میشه و باید بیشتر مطالعه کنم.
متاسفانه یک مدته که عکسها آپلود نمیشن و گرنه میخواستم که عکسهای منطقه و برش رو وهمچنین گزارش تصویریه بازدید 05/05/1389 رو قرار بدم به هر حال هرکی PDF رو گیر آورد به صفحات30 و 32 نگاه کنه. اصلاً بعداً خودم قرار میدمش ببینین.
با تشکر
واحددوست.

دنباله مشكل مدل سازي

دوباره سلام،
راستش در مطلب قبلي بهتر بود كه نقشه كد ها و ترازها رو هم قرار مي دادم. اينجوري بهتر ميشد.
به هر حاب بايد بگم كه تراز روي تاج سد ما 98+ خواهد بود بنا بر اين با توجه به مقاطع ژئو الكتريك طولي در دو مسير A و B كه بعدا نقشه هاشو ميزارم ببينين عملا ما بالاتر از تراز 103+ روي بستر رودخانه تراز ي نخواهيم داشت چون زير اين تراز لايه سنگ به تراز 98+ ميرسه و چون دربالاتر از اين تراز تاج سد زه كش مي شه و عملا براي شرايط محيطي ما لايه مرزي ميشه مگر اينكه بخواهيم نقش كاپيلاريته رو هم در نظر بگيريم. به نظر شخصي خودم لازم نيست ولي بايد با جناب آقاي دكتر مشورت كنم.
حالا اين مساله با احتمال 90% همونجوريه كه فكر ميكنم، مشكل اساسي تر در زمينه اينه كه با حركت در عمق لايه ها هم لايه صفرمون (غير فعال) فضاي كمتري براي سلولهاي فعال باقي ميزارن كه اين مساله بدليل V شكل بودن دره است كه البته خودشم متقارن نيست از طرف ديگر هم جناح چپ كه داراي لايه ماسه لاي دار هستش خودش داراي ضخامت متغير در عمق مقطع هست و اين باعث ميشه نگران اين باشم كه مدل با اين شرايط چطور رفتار ميكنه، شايد بهتر باشه يه حالت متوسط بين شن و ماسه بگيرم. در مورد مش بندي هم فكر مي كنم مربع هاي 10 متر در 10 متر فعلا كافيه، جاهائيكه لازم باشه هم تقسيم بندي موضعي ميكنم.
با تشكر
واحددوست

يك مشكل

سلام دوستان
همونطور كه گفتم الان دارم محيط رو مدل ميكنم. اما يك مشكلي كه وجود داره اينه كه با توجه به مقطع شكل زير كه مقطع ژئوالكتريك در محور سد ما هست چطور مش بندي رو مي تونم خوب انجام بدم تا بهترين جواب رو بدست بيارم. در ضمن يكجورايي بايد اين مقطع را در دوتا مقطع ژئوالكتريك طولي سرشكن بكنم (يكي جناح چپ و يكي راست) چون فقط يك مقطع عرضي داريم.
ببينيد مقطع را: