احتراماً چند لينك بدرد بخور قرار ميدم اگه علاقه مند بوديد برين ببينين:
Finite difference
Finite element method
History of calculus
Water
Groundwater
Computer simulation
GIS
Subsurface flow
واحددوست
مطالعه بیلان آبی مخازن سدهای زیرزمینی با استفاده از نرم افزار MODFLOW (مطالعه موردی: سد زیرزمینی مسیل روستای صوفی-ماکو) Study of water balance in subsurface dam's cachment, using MODFOLW software
با سلام.
بله دوستان، همونطور که از تیتر مشخص میشه بنده برای اولین مدل کالیبراسیون شده سد زیرزمینی صوفی مربوط به اطلاعات شهریور 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
این حالت قبل از PMPATH هستش. همونطور که می بینی سیستم زه کش، رودخانه و سلولهای 1- مشخص هستند. چند تا المان هم تو جاهای مختلف در نظر گرفتم تا Particel tracking کنمشون
تو عکس زیر هم مشاهده می کنین، همه چیز گویا هستش. تازه همونجور که گفتم می بینی که زه کش عمل نمی کنه و فقط داره از خودش تراوش می ده که تو عمل هم مهندس جعفری اینا با این مشکل مواجه شده بودند. رودخانه هم تو بالا دست داره تراوش می ده ولی تو پایین دست داره آب رو زهکشی میکنه. خیلی جالبه نه؟ اگه تو مقطع کناری هم ببینین میشه دید که آب میاد از روی تاج سد سرریز میشه!!!!!!
با تشکر / واحددوست/ B.V.Dst
به هر حال کلاً تو استرس هستم چون حد اکثر تا اواسط بهمن باید دفاع کنم ولی هنوز اولین مدل هم درست و حساب run نشده.
برام دعا کنید.
با تشکر
واحددوست.
ρs . g.Z
و فشاری که از بالا وارد می شود معادل فشار آی شیرین و عبارت است از:
ρf . g.(Z+h)
بنا بر این خواهیم داشت
(ρs . g.Z = ρf . g.(Z+h
که در آن ρs : دانسیته آب شور، ρf: دانسیته آب شیرین، h: فاصله سطح ایستابی تا تراز آب سطح دریا (بار پیزومتری نسبت به سطح دریا) و Z: عمق جبهه آب شور نسبت به سطح دریا در هر نقطه از فرمول زیر می توان Z را محاسبه کرد.
Z= [(ρf) / (ρs- ρf)].h
چون دانسیته آب شور 1/025 و دانسیته آب شیرین 1 گرم در سانتیمتر مکعب می باشد لذا:
Z = [(1.0) / (1.025-1.0)]*h
Z=40h
بدین ترتیب مشاهده می شود که اگر سطح ایستایی در نواحی ساحلی به اندازه h متر بالاتر از تراز سطح آب دریا قرار گیرد تا عمق 40h در زیر آب شیرین وجود خواهد داشت و هر زمان که این آب به پایین تر از سطح دریا برسد آب شور به سمت ساحل حرکت خواهد کرد."
برای ادامه مطلب رجوع شود به صفحه 443 کتاب.
و دوستان این نشان دهنده ترسناکترین و خطرناکترین و مطرحترین مساله در موضوع نفوذ آب شور است، چراکه این بدان معنی است که اگر فقط 1متر وتنها 1 متر از سطح آب لایه آبدار در اثر براشت بی رویه کاسته شود، آب شور 40 متر نسبت به مرز قبلی خود بیشتر به سمت داخل نفوذ خواهد کرد.
با تشکر
واحددوست
Babak.V.Dst