اندکی درباره مدل MODFLOW:

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

از جمله ویژگی های دیگر MODFLOW می توان به ساختار ماژولی آن اشاره نمود. بدین ترتیب که می توان برای فرآیندهای خاص هیدرولوژیکی ماژول خاصی را فعال یا غیرفعال نمود. همچنین مدول های جدیدی برای مسایل مربوط به جریان (مثل اندرکنش جریان رودخانه با آبخوان) یا روش های عددی جدید در حال توسعه می باشند. این دلایل سبب شده MODFLOW از سوی بسیاری از سازمان ها به عنوان نرم افزار برتر مورد استفاده قرار گیرد.

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

- فرم هندسی محدوده آبخوان

- نوع تشکیلات زمین شناسی آبخوان از نظر همگنی و ناهمگنی

- نحوه بررسی مساله (به صورت یک، دو و یا سه بعدی)

- تعیین رژیم جریان به صورت ورقه ای و یا متلاطم

- تعیین نوع منبع آلودگی (نقطه ای، خطی، توزیعی)

- نوع آلاینده ها و نحوه توزیع آن

- بررسی تغییر مکانی و زمانی متغیر حالت (غلظت) در آبخوان

- تغذیه و تخلیه مواد محلول در سامانه کیفی آبخوان

- شرایط مرزی و ارتباطی که از نظر کیفیت، آبخوان، با خارج از محدوده خود دارد.

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

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

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

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

میان یابی داده ها و توزیع آنها بر روی گره های مدل

انواع داده های مورد نیاز در مدل سازی آب های زیرزمینی در بخش های گذشته این نوشتار مورد اشاره قرار گرفته است. داده های صحرایی پس از جمع آوری و اصلاح، باید به عنوان مقادیر پارامتری یا متغیر با روش های مناسب برای گره ها یا سلول ها/المان ها تعریف شوند. چگونگی توزیع داده های گوناگون در گستره مدل بسیار مهم می باشد.

تطبیق پارامترها با نوع مدل

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

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

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

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

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

روش الکترومغناطیس

روش قطبش القایی

آژانس های آب و شرکت های خصوصی آب با چالش توسعه و حفظ منابع آب پایدار، اغلب در آب های زیرزمینی به عنوان محور وظایف خود روبرو می باشند. از سال 1987، سیستم Westbay برای مدیران آب، فضای انعطاف پذیری را فراهم کرده است که به توصیف آب های زیرزمینی در چهار بعد پرداخته، و منجر به بهینه سازی فعالیت های مدیریت و توسعه شیوه های مدیریت پایدار گشته است.

نرم افزار مدیریت منابع آب زیرزمینی

• مدیریت حوضه آب های زیرزمینی
• عملیات MAR
• نفوذ آب دریا، نظارت، پیشگیری و کاهش
• بهینه سازی مدل های عددی

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

نرم افزار تصویر سازی، تجزیه و تحلیل و مدیریت داده های محیطی

نرم افزار Hydro GeoAnalyst یک فضای یگانه برای مدیریت سیستمی داده های آب های زیرزمینی و محیطی است؛ که سبب سازماندهی داده های محیطی شده و همچنین به شما برای افزایش سرعت ایجاد نتایج تحلیلی-جامع که قابل اعتماد و دارای مفاهیم ساده ای جهت درک می باشد، کمک می نماید. با Hydro GeoAnalyst، متخصصان آب های زیرزمینی، علوم زمین و سازمان های دولتی تمامی ابزارهای مورد نیازی که به طور موثری برای مدیریت حجم زیادی از داده ها، و تصمیم گیری آگاهانه در مورد محیط زیست و منابع آب کارامد است، در اختیار خواهند داشت.

واسنجی مدل در شرایط ماندگار و غیرماندگار

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

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

شرایط اولیه

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

ارتباط شرط مرزی و روش محاسبات

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

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

انواع اصلی شرایط مرزی

شرایط مرزی اصلی که در واقع شرایط مرزی خارجی مدل های آب زیرزمینی هستند، به سه دسته زیر تقسیم می شوند:

نوع اول - شرط مرزی بار هیدرولیکی معلوم:

این نوع شرط مرزی در نهایی وجود دارد که تراز آب و یا بار هیدرولیکی مشخصی به صورت ثابت یا متغیر با زمان وجود داشته باشد. مثل نواحی ساحلی، دریاچه ها یا مخازن و رودخانه ها که دارای بار هیدرولیکی معلوم هستند.

در بحث مدل آبهای زیرزمینی این نوع مرزها مرز با بار ثابت و یا مرز وابسته به بار هیدرولیکی خوانده می شوند. در بحث ریاضی، این نوع مرزها به نام های Dirichlet یا Stable یا Essential خوانده می شوند.

شرایط مرزی و شرایط اولیه

در این قسمت با معرفی شرایط مرزی و شرایط اولیه، مسایل مهم مربوط به تعیین این شرایط در هنگام شبیه سازی آب های زیرزمینی تشریح می شود.

شرایط مرزی

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

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

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

ویژگی های نرم افزار های مورد نیاز برای مدل سازی آب زیرزمینی

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

از نظر میزان استفاده از برنامه های رایانه ای، دو نوع کد وجود دارد: کدهای عمومی با قابلیت دسترسی برای همگان مانند برنامه 22MT3 و کدهای اختصاصی که معمولا برای مصارف ویژه و به سفارش مراکز تحقیقاتی خاصی توسعه یافته و در دسترس همه نمی باشد.

مدل مفهومی

معمولا پس از تعیین اهداف مدل سازی دومین گام در فرآیند مدل سازی و در حقیقت مهمترین مرحله در مطالعات مدل تهیه یک مدل مفهومی می باشد. مدل مفهومی یک تصویر ساده شده از دنیای واقعی است. این مدل خلاص های از ویژگی های سامانه هیدروژئولوژیکی است که ضمن برخورداری از دقت لازم معمولا به صورت توصیفی و گرافیکی (نیمرخ زمین شناسی هیدروژئولوژی). مدل مفهومی، پایه مدل ریاضی است و خود بر پایه - و یا نمودار بلوکی) ارائه دهنده ویژگی های اصلی سامانه می باشد اطلاعات اولیه داده های صحرایی و تعبیر و تفسیرهای هیدروژئولوژیکی استوار می باشد.

پروتکل مدل سازی

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