فهرست مطالب عنوان صفحه فصل اول: مقدمه1-1- پیشگفتار....... 21-2- اهمیت موضوع...... 31-3- اهداف تحقیق........ 41-4- نوآوری پایان نامه........ 5 1-5- ساختار پایان نامه....... 5 فصل دوم: مبانی نظری تحقیق و مدلهاي شبيهسازي کمی و کیفی جریان آبهاي زيرزميني2-1- تعریف مدل آب زیرزمینی.......72-2- انواع مدلها........ 82-2-1- مدلهای فیزیکی....... 92-2-2- مدلهای آنالوگ........92-2-3- مدلهای ریاضی......... 102-2-3-1- مدلهای تجربی....... 112-2-3-2- مدلهای احتمالاتی........ 122-2-3-3- مدلهای علت ومعلولی...... 122-3-نرم افزارهای مدلسازی جریان و پخش و انتقال آلودگی در آبهای زیرزمینی..... 15عنوان صفحه 2-3-1- MODFLOW................................................................................................... 152-3-2- Visual MODFLOW................................................................................... 162-3-3- MIKE SHE...................................................................................................... 172-3-4- FEFLOW..........................................................................................................172-3-5- CTRAN/W......................................................................................................182-3-6- MT3D................................................................................................................182-3-7- MODPATH.................................................................................................... 192-3-8- Groundwater Modeling System; GMS.............................................. 202-3-8-1- Map Module................................................................................... 222-3-8-2- TIN Module.................................................................................... 232-3-8-3- 2D Scatter Point Module.......................................................... 242-3-8-4- Solid Module................................................................................. 252-3-8-5- 3D Scatter Point Module ........................................................ 252-3-8-6- Modflow........................................................................................... 262-4- مبانی فیزیکی و معادلات حاكم بر جريان آب زيرزميني.............................................. 272-4-1- قانون دارسي................................................................................................................. 272-4-2- معادله پيوستگي.......................................................................................................... 292-4-3- حل معادلات حاكم بر جريان در آبهاي زيرزميني............................................332-4-3-1- روش عناصر محدود.................................................................................... 342-4-3-2- روش تفاضلات محدود................................................................................ 342-5- آلودگي آبهاي زيرزميني................................................................................................... 412-5-1- منابع آلودگي آبهاي زيرزميني.............................................................................. 41 عنوان صفحه 2-5-2- حرکت آلودگی در آبهای زیرزمینی......................................................................... 432-5-2-1- حرکت همراه با جریان................................................................................ 442-5-2-2- انتشار هیدرودینامیکی................................................................................ 442-5-2-3- تأخیر...............................................................................................................462-5-2-4- واکنش شیمیایی........................................................................................... 47 فصل سوم: مروری بر تحقیقات انجام شده3-1- مروري بر مطالعات گذشته.................................................................................................. 49 فصل چهارم: روش تحقیق4-1- محدوده مطالعاتی.................................................................................................................. 614-1-1- هواشناسی...................................................................................................................... 624-1-2- ایستگاههای هواشناسی............................................................................................... 634-1-3- دما.................................................................................................................................... 664-1-4- باران................................................................................................................................. 684-1-5- تبخیر............................................................................................................................... 694-2- زمینشناسی منطقه...............................................................................................................714-3- مراحل ساخت و آماده سازي مدل منطقه مورد مطالعه................................................. 744-3-1- مدل مورد استفاده در اين مطالعه............................................................................. 744-3-2- تهيه مدل مفهومي........................................................................................................764-3-3- شبکه بندی مدل..........................................................................................................774-3-4- اعمال توپوگرافي سطح و كف آبخوان به محدوده مدل........................................ 784-3-5- اعمال شرايط مرزي به محدوده مدل....................................................................... 80عنوان صفحه 4-3-6- اعمال سطح آب مشاهدهاي........................................................................................ 824-3-7- هدايت هيدروليكي و آبدهي ويژه آبخوان............................................................... 844-3-8- برآورد ميزان تغذيه سطحي........................................................................................ 874-3-9- ايجاد لايه اطلاعاتي مربوط به پيزومترها................................................................. 874-4- اجرا و واسنجي مدل.............................................................................................................. 874-5- آناليز حساسيت.......................................................................................................................924-6- مدل کيفي............................................................................................................................... 924-6-1- تعیین ناحیه گیرش چاههای با اهمیت.................................................................... 934-6-2- اثر احداث یک لندفیل فرضی.................................................................................... 94 فصل پنجم: نتایج وبحث5-1- نتايج واسنجي مدل در شرايط پايدار................................................................................. 985-1-1- مقادیر محاسباتی جهت پارامتر هدایت هیدرولیکی............................................. 985-1-2- مقادیر محاسباتی جهت پارامتر تغذیه سطحی ..................................................... 995-1-3- توزیع سطح ایستابی و جهت جریان...................................................................... 1005-1-4- مقادير محاسباتي و مشاهداتي سطح آب زيرزميني...........................................1015-2- نتايج واسنجي مدل در شرايط ناپايدار........................................................................... 1035-2-1- مقادیر محاسباتی جهت پارامتر هدایت هیدرولیکی......................................... 1035-2-2- مقادیر محاسباتی جهت پارامتر آبدهی ویژه...................................................... 1045-2-3- مقادير محاسباتي و مشاهداتي سطح آب زيرزميني......................................... 1055-2-4- مقادیر خطا در شرایط ناپایدار............................................................................... 1085-2-5- بیلان........................................................................................................................... 1085-3- آنالیز حساسیت.................................................................................................................... 110عنوان صفحه 5-3-1- حساسیت مدل نسبت به تغییرات پارامتر هدایت هیدرولیکی....................... 1105-3-2- حساسیت مدل نسبت به پارامتر آبدهی ویژه..................................................... 1105-3-3- حساسیت مدل نسبت به پارامتر تغذیه سطحی................................................ 1115-4- صحت سنجي....................................................................................................................... 1125-5- نتایج مدلسازی کیفی....................................................................................................... 1155-5-1- تعیین ناحیه گیرش چاههای حائز اهمیت.......................................................... 1155-5-1-1- نتایج بررسی چند چاه.............................................................................. 1165-5-2- نتایج شبیهسازی لندفیل فرضی............................................................................ 1195-5-2-1- مکانیابی حرکت ذرات شیرابه در آب زیرزمینی............................... 1195-5-2-2- روند پخش و انتقال آلودگی ................................................................... 120 فصل ششم: نتیجه گیری و پیشنهادات6-2- نتیجه گیری.......................................................................................................................... 1396-3- پیشنهادات............................................................................................................................ 140- فهرست منابع............................................................................................................................. 141- چکیده به زبان انگلیسی فهرست جدولها عنوان و شماره صفحه جدول 4-1: موقعیت ایستگاههای هواشناسی منطقه مورد مطالعه.................................................63جدول 4-2: معادلات گرادیانهای ماهانه حرارتی محدوده..............................................................67جدول 4-3: متوسط دمای ماهانه محدوده به تفکیک دشت و ارتفاعات محدوده ایج (درجه سانتیگراد).......................................................................................................................................67جدول 4-4: متوسط بارش ماهانه و سالانه در ایستگاه معرف، دشت و ارتفاعات محدوده ایج (میلیمتر).................................................................................................................................................... 69جدول 4-5: متوسط تبخیر ماهانه محدوده به تفکیک دشت و ارتفاعات محدوده مطالعاتی ایج (میلیمتر).....................................................................................................................................................71جدول 4-6: مقدار هدايت هيدروليکي تشکيلات مختلف................................................................ 85جدول 4-7: مقدار آبدهي ويژه تشکيلات مختلف...............................................................................86جدول 5-1: مقادیر بیلان جریان آب زیرزمینی منطقه مورد مطالعه در طول دوره واسنجی شرایط ناپایدار (سالهای آبی 1390-1388).................................................................................. 109جدول 5-2: طول کمینه، متوسط و بیشینه محدودهی گیرش مهمترین چاه آب شرب شهر ایج .....................................................................................................................................................................117جدول 5-3: طول کمینه، بیشینه، متوسط و مساحت محدودهی گیرش چاه شماره یک.......118جدول 5-4: طول کمینه، بیشینه، متوسط و مساحت محدودهی گیرش چاه شماره دو.........118 فهرست شکلها عنوان صفحه شکل 2-1- تقسیم بندی مدلهای آب زیرزمینی .............................................................................. 8شکل 2-2- نحوه تغییرات هد در ستون ماسه در آزمایش دارسی ............................................... 28شکل 2-3- جریان ورودی و خروجی از المان حجمی ................................................................... 30شکل 2-4- هیدروگراف برای سلول i, j, k ....................................................................................... 39شکل 4-1- موقعیت دشت ایج فارس ................................................................................................ 62شکل 4-2- موقعیت چاههای بهره برداری و مشاهداتی منطقه مورد مطالعه .............................77شکل 4-3- شبکه بندی منطقه مورد مطالعه .................................................................................. 78شکل 4-4- نقشه DEM توپوگرافی سطح زمین منطقه مورد مطالعه ....................................... 79شکل 4-5- نقشه DEM رقوم ارتفاعی سنگ بستر منطقه مورد مطالعه .................................. 80شکل 4-6- هیدروگراف واحد دشت ایج طی سالهای 1375-1391 ...................................... 83شکل 4-7- نقشه سطح آب (متر) مهرماه 1388 دشت مورد مطالعه ....................................... 83شکل 4-8- میله رنگی نمایانگر خطای واسنجی ............................................................................ 90شکل 4-9- محل دفن زباله فرضی برای منطقه مورد مطالعه ................................................... 94شکل 5-1- نقشه زونبندي و مقادير به دست آمده براي پارامتر هدايت هيدروليکي (متر بر روز) طي کاليبراسيون مدل در شرايط پايدار (مهرماه 1388)..................................................... 98شکل 5-2- نقشه زون بندي و مقادير به دست آمده براي پارامتر تغذیه سطحی (متر بر ماه) طي کاليبراسيون مدل در شرايط پايدار (مهرماه 1388) ............................................................... 99شکل 5-3- نقشه توزيع سطح ايستابي به دست آمده طي کاليبراسيون مدل در شرايط پايدار (مهرماه 1388) ..................................................................................................................................... 100شکل 5-4- جهت حرکت جریان آب زیرزمینی در منطقه مورد مطالعه ................................. 101شکل 5-5- مقادير محاسباتي و مشاهداتي بار هيدروليکي در چاههاي مشاهداتي مختلف محدوده مطالعاتي در دوره واسنجي شرايط پايدار (مهرماه 1388) ........................................... 102شکل 5-6- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در دوره واسنجي شرايط پايدار (مهرماه 1388) ............................................................................................................ 102شکل 5-7- مقادير مشاهداتي بار هيدروليکي در مقابل مقادير باقيمانده (تفاوت مقادير محاسباتي و مشاهداتي) در دوره واسنجي شرايط پايدار (مهرماه 1388) ............................... 103شکل 5-8- نقشه زونبندي و مقادير به دست آمده براي پارامتر هدايت هيدروليکي (متر بر روز) طي کاليبراسيون مدل در شرايط ناپايدار (سالهای آبي 1390-1388) ...................... 104شکل 5-9- نقشه زون بندي و مقادير به دست آمده براي پارامتر آبدهي ويژه طي کاليبراسيون مدل در شرايط ناپايدار (سالهای آبي 1390-1388) ............................................................... 105شکل 5-10- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره یک در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) ................. 106شکل 5-11- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره دو در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) ................... 106شکل 5-12- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره سه در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) ................. 107شکل 5-13- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره چهار در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) .............. 107شکل 5-14- حساسیت مدل نسبت به تغییرات پارامتر هدایت هیدرولیکی ........................ 110شکل 5-15- حساسیت مدل نسبت به تغییرات پارامتر آبدهی ویژه .................................... 111شکل 5-16- حساسیت مدل نسبت به تغییرات پارامتر تغذیه سطحی ................................. 112شکل 5-17- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره یک در طول دوره صحت سنجی (سال آبي 1391-1390) ........................................ 113شکل 5-18- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره دو در طول دوره صحت سنجی (سال آبي 1391-1390) ........................................... 113شکل 5-19- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره سه در طول دوره صحت سنجی (سال آبي 1391-1390) .......................................... 114شکل 5-20- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره چهار در طول دوره صحت سنجی (سال آبي 1391-1390) ........................................ 114شکل 5-21- محدوده گیرش چاههای با اهمیت منطقه مورد مطالعه پس از گذشت بیست سال ....................................................................................................................................................................116شکل 5-22- ناحیه گیرش چاه آب شرب شهر ایج پس از گذشت ده سال ............................ 117شکل 5-23- ناحیه گیرش دو چاه نمونه در دشت ایج ............................................................... 118شکل 5-24- مکانیابی حرکت ذرات شیرابه در آب زیرزمینی ناشی از وجود لندفیل ......... 119شکل 5-25- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از یک سال، مکانيسم انتقال جرم: پخش، انتشار) .................................................................................................................................................... 121شکل 5-26- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم:پخش، انتشار) .................................................................................................................................................... 122شکل 5-27- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم:پخش، انتشار) ................................................................................................................................................... 122شکل 5-28- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم:پخش، انتشار) ............................................................................................................................ 123شکل 5-29- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از یک سال، مکانيسم انتقال جرمپخش، انتشار و تأخیر) ..................................................................................................................................... 123شکل 05-30- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ..................................................................................................................................... 124شکل 5-31- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ..................................................................................................................................... 124شکل 5-32- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) .............................................................................................................. 125شکل 5-33- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از یک سال، مکانيسم انتقال جرم: پخش و انتشار) ..................................................................................................................................................... 125شکل 5-34- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم: پخش و انتشار) ..................................................................................................................................................... 126شکل 5-35- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم: پخش و انتشار) ................................................................................................................................................. 126شکل 5-36- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم: پخش و انتشار) ........................................................................................................................... 127شکل 5-37- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از یک سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) .................................................................................................................................... 127شکل 5-38- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ...................................................................................................................................... 128شکل 5-39- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ..................................................................................................................................... 128شکل 5-40- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) .............................................................................................................. 129شکل 5-41- توزیع ابر آلودگی پس از گذشت بیست و هشت سال در شرایط نرخ نفوذ زیاد و مکانیسم انتقال پخش و انتشار به همراه موقعیت چاه شماره شش ........................................... 130شکل شماره 5-42- موقعیت سلولهای دلخواه مورد نظر جهت بررسی جزئی تر نتایج تغییر غلظت شیرابه در طول زمان ............................................................................................................... 131شکل 5-43- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ بالا، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 132شکل 5-44- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ بالا، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 133شکل 5-45- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ کم، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 134شکل 5-46- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ کم، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 134شکل 5-47- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ بالا، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 135شکل 5-48- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ بالا، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 136شکل 5-49- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ کم، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ............................................................................ 136شکل 5-50- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ کم، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 137 آب مایه حیات و فراوانترین ماده مرکب برروی سطح کره زمین و بستر اولیه حیات به شکلی که امروزه میشناسیم، میباشد. بیش از ۷۰٪ سطح کره زمین را آب پوشانده است (نزدیک به ۳۶۰ میلیون از ۵۱۰ میلیون کیلومتر مربع)؛ با وجود این حجم عظیم آب، تنها 2% از آبهای کره زمین شیرین و قابل شرب است و باقی آن به علت محلول بودن انواع نمکها خصوصاً نمک طعام غیر قابل استفاده است. از همین 2% آب شیرین، بیش از ۹۰% به صورت منجمد در دو قطب زمین و دور از دسترس بشر واقع شده است (Davie, 2002). به علاوه، منابع آب شیرین به طور یکنواخت در سطح زمین پراکنده نشدهاند. در حال حاضر، 60% کل منابع آب شیرین در 9 کشور جهان وجود دارد؛ در مقابل حدود 80 کشور با کمبود آب مواجهاند که برخی از آنها تقریباً به هیچ منبع آب شیرین قابل توجهی دسترسی ندارند (http://fa.wikipedia.org/wiki/%D8%A2%D8%A8#cite_note-autogenerated4-9). طبق آمار برنامه عمران سازمان ملل متحد در سال 2006:بنابراین، نحوه و چگونگی استفاده از منابع آب موجود در اکوسیستم طبیعی در سالهای اخیر موضوع بحث و بررسی علوم مختلف بوده است. امروزه مدیریت منابع آب[1] با در نظر گرفتن آنالیز همزمان سیستمهای آبی، اعم از سطحی و زیرزمینی و بررسی اثرات متقابل این منابع بر یکدیگر که از آن با نام مدیریت تلفیقی[2] منابع یاد میشود، به عنوان روش مناسب کارآمد مدیریتی در زمینه شناخت و برنامه ریزی هیدروسیستمها مطرح میشود. 1-2- اهمیت موضوع براساس آمار آژانس اطلاعات مرکزی (CIA[3]) در سال 2011 میلادی، کشور ایران با مساحت 1648195 کیلومتر مربع در رتبهی هجدهم در لیست وسیعترین کشورهای دنیا قرار دارد. براساس این آمار 3/0% خشکیهای کرهی زمین در محدودهی مرزهای ایران است. این در حالی است که تنها 03/0% آبهای موجود به ایران تعلق دارد. اهمیت این موضوع هنگامی بیشتر میشود که بدانیم 09/1% جمعیت جهان در ایران قرار دارد. از طرفی در آمار منتشر شده بین سالهای 1930 تا 1960 میلادی، کشور ما با 230 میلیمتر میانگین بارندگی سالانه در رتبه هشتاد و چهارم در میان کل کشورها است. دردناکتر اینکه میانگین بارندگی سالانه در میان کل کشورها 963 میلیمتر است(CIA, 2011).متأسفانه در سالهای اخیر برداشتهای بیش از ظرفیت از منابع آبهای زیرزمینی باعث ایجاد حالت بحرانی در اکثر دشتهای کشور ایران شدهاست. از طرف دیگر همزمانی آن با دورههای خشکسالی که دورههای تناوب آن بتدریج کوتاه شده است، مشکلاتی را برای توسعه بخشهای کشاورزی، اقتصادی و حتی تامین آب شرب به وجود آورده است. بدین ترتیب سطح آبهای زیرزمینی در منطقه روز به روز افت کرده و خشکیدن قنوات و چشمهها، کاهش آبدهی چاهها، فرونشست زمین، توسعه بیابان و در نهایت از دست رفتن سرمایه گذاریهای انجام شده در بخش کشاورزی، صنعت و مهاجرت را به همراه خواهد داشت (اکبرپور و همکاران، 1389).با نگاهی به مطالب ذکر شده به نظر میرسد اگر شرایط به همین روند ادامه یابد و مدیریت صحیح در زمان مناسب اعمال نشود، سرانجام به جایی خواهد رسید که آب مناسبی برای استخراج وجود نخواهد داشت و کشور با بحران شدید منابع آبی در قسمتهای مختلف روبرو خواهد شد.همچنین افزایش برداشت از آب زیرزمینی، ورود انواع پسابهای صنعتی، آب برگشتی کشاورزی آلوده به کودهای فسفاته و نیتراته، نشت از مخازن و خطوط فراوردههای نفتی و ... سبب شده است که کیفیت این منابع رو به نامناسب شدن پیش رود. اگر روند آلودگی، پخش و انتشار آن در آبهای زیرزمینی که سهم عمدهای در مصارف شرب و کشاورزی دارد مدیریت و کنترل نشود مخاطرات جبران ناپذیری برای سلامت بشر و محیط زیست خواهد داشت.
شبيهسازي كمي و كيفي جريان آبهای زيرزمينی دشت ايج با استفاده از نرمافزار GMS
فهرست مطالب عنوان صفحه فصل اول: مقدمه1-1- پیشگفتار....... 21-2- اهمیت موضوع...... 31-3- اهداف تحقیق........ 41-4- نوآوری پایان نامه........ 5 1-5- ساختار پایان نامه....... 5 فصل دوم: مبانی نظری تحقیق و مدلهاي شبيهسازي کمی و کیفی جریان آبهاي زيرزميني2-1- تعریف مدل آب زیرزمینی.......72-2- انواع مدلها........ 82-2-1- مدلهای فیزیکی....... 92-2-2- مدلهای آنالوگ........92-2-3- مدلهای ریاضی......... 102-2-3-1- مدلهای تجربی....... 112-2-3-2- مدلهای احتمالاتی........ 122-2-3-3- مدلهای علت ومعلولی...... 122-3-نرم افزارهای مدلسازی جریان و پخش و انتقال آلودگی در آبهای زیرزمینی..... 15عنوان صفحه 2-3-1- MODFLOW................................................................................................... 152-3-2- Visual MODFLOW................................................................................... 162-3-3- MIKE SHE...................................................................................................... 172-3-4- FEFLOW..........................................................................................................172-3-5- CTRAN/W......................................................................................................182-3-6- MT3D................................................................................................................182-3-7- MODPATH.................................................................................................... 192-3-8- Groundwater Modeling System; GMS.............................................. 202-3-8-1- Map Module................................................................................... 222-3-8-2- TIN Module.................................................................................... 232-3-8-3- 2D Scatter Point Module.......................................................... 242-3-8-4- Solid Module................................................................................. 252-3-8-5- 3D Scatter Point Module ........................................................ 252-3-8-6- Modflow........................................................................................... 262-4- مبانی فیزیکی و معادلات حاكم بر جريان آب زيرزميني.............................................. 272-4-1- قانون دارسي................................................................................................................. 272-4-2- معادله پيوستگي.......................................................................................................... 292-4-3- حل معادلات حاكم بر جريان در آبهاي زيرزميني............................................332-4-3-1- روش عناصر محدود.................................................................................... 342-4-3-2- روش تفاضلات محدود................................................................................ 342-5- آلودگي آبهاي زيرزميني................................................................................................... 412-5-1- منابع آلودگي آبهاي زيرزميني.............................................................................. 41 عنوان صفحه 2-5-2- حرکت آلودگی در آبهای زیرزمینی......................................................................... 432-5-2-1- حرکت همراه با جریان................................................................................ 442-5-2-2- انتشار هیدرودینامیکی................................................................................ 442-5-2-3- تأخیر...............................................................................................................462-5-2-4- واکنش شیمیایی........................................................................................... 47 فصل سوم: مروری بر تحقیقات انجام شده3-1- مروري بر مطالعات گذشته.................................................................................................. 49 فصل چهارم: روش تحقیق4-1- محدوده مطالعاتی.................................................................................................................. 614-1-1- هواشناسی...................................................................................................................... 624-1-2- ایستگاههای هواشناسی............................................................................................... 634-1-3- دما.................................................................................................................................... 664-1-4- باران................................................................................................................................. 684-1-5- تبخیر............................................................................................................................... 694-2- زمینشناسی منطقه...............................................................................................................714-3- مراحل ساخت و آماده سازي مدل منطقه مورد مطالعه................................................. 744-3-1- مدل مورد استفاده در اين مطالعه............................................................................. 744-3-2- تهيه مدل مفهومي........................................................................................................764-3-3- شبکه بندی مدل..........................................................................................................774-3-4- اعمال توپوگرافي سطح و كف آبخوان به محدوده مدل........................................ 784-3-5- اعمال شرايط مرزي به محدوده مدل....................................................................... 80عنوان صفحه 4-3-6- اعمال سطح آب مشاهدهاي........................................................................................ 824-3-7- هدايت هيدروليكي و آبدهي ويژه آبخوان............................................................... 844-3-8- برآورد ميزان تغذيه سطحي........................................................................................ 874-3-9- ايجاد لايه اطلاعاتي مربوط به پيزومترها................................................................. 874-4- اجرا و واسنجي مدل.............................................................................................................. 874-5- آناليز حساسيت.......................................................................................................................924-6- مدل کيفي............................................................................................................................... 924-6-1- تعیین ناحیه گیرش چاههای با اهمیت.................................................................... 934-6-2- اثر احداث یک لندفیل فرضی.................................................................................... 94 فصل پنجم: نتایج وبحث5-1- نتايج واسنجي مدل در شرايط پايدار................................................................................. 985-1-1- مقادیر محاسباتی جهت پارامتر هدایت هیدرولیکی............................................. 985-1-2- مقادیر محاسباتی جهت پارامتر تغذیه سطحی ..................................................... 995-1-3- توزیع سطح ایستابی و جهت جریان...................................................................... 1005-1-4- مقادير محاسباتي و مشاهداتي سطح آب زيرزميني...........................................1015-2- نتايج واسنجي مدل در شرايط ناپايدار........................................................................... 1035-2-1- مقادیر محاسباتی جهت پارامتر هدایت هیدرولیکی......................................... 1035-2-2- مقادیر محاسباتی جهت پارامتر آبدهی ویژه...................................................... 1045-2-3- مقادير محاسباتي و مشاهداتي سطح آب زيرزميني......................................... 1055-2-4- مقادیر خطا در شرایط ناپایدار............................................................................... 1085-2-5- بیلان........................................................................................................................... 1085-3- آنالیز حساسیت.................................................................................................................... 110عنوان صفحه 5-3-1- حساسیت مدل نسبت به تغییرات پارامتر هدایت هیدرولیکی....................... 1105-3-2- حساسیت مدل نسبت به پارامتر آبدهی ویژه..................................................... 1105-3-3- حساسیت مدل نسبت به پارامتر تغذیه سطحی................................................ 1115-4- صحت سنجي....................................................................................................................... 1125-5- نتایج مدلسازی کیفی....................................................................................................... 1155-5-1- تعیین ناحیه گیرش چاههای حائز اهمیت.......................................................... 1155-5-1-1- نتایج بررسی چند چاه.............................................................................. 1165-5-2- نتایج شبیهسازی لندفیل فرضی............................................................................ 1195-5-2-1- مکانیابی حرکت ذرات شیرابه در آب زیرزمینی............................... 1195-5-2-2- روند پخش و انتقال آلودگی ................................................................... 120 فصل ششم: نتیجه گیری و پیشنهادات6-2- نتیجه گیری.......................................................................................................................... 1396-3- پیشنهادات............................................................................................................................ 140- فهرست منابع............................................................................................................................. 141- چکیده به زبان انگلیسی فهرست جدولها عنوان و شماره صفحه جدول 4-1: موقعیت ایستگاههای هواشناسی منطقه مورد مطالعه.................................................63جدول 4-2: معادلات گرادیانهای ماهانه حرارتی محدوده..............................................................67جدول 4-3: متوسط دمای ماهانه محدوده به تفکیک دشت و ارتفاعات محدوده ایج (درجه سانتیگراد).......................................................................................................................................67جدول 4-4: متوسط بارش ماهانه و سالانه در ایستگاه معرف، دشت و ارتفاعات محدوده ایج (میلیمتر).................................................................................................................................................... 69جدول 4-5: متوسط تبخیر ماهانه محدوده به تفکیک دشت و ارتفاعات محدوده مطالعاتی ایج (میلیمتر).....................................................................................................................................................71جدول 4-6: مقدار هدايت هيدروليکي تشکيلات مختلف................................................................ 85جدول 4-7: مقدار آبدهي ويژه تشکيلات مختلف...............................................................................86جدول 5-1: مقادیر بیلان جریان آب زیرزمینی منطقه مورد مطالعه در طول دوره واسنجی شرایط ناپایدار (سالهای آبی 1390-1388).................................................................................. 109جدول 5-2: طول کمینه، متوسط و بیشینه محدودهی گیرش مهمترین چاه آب شرب شهر ایج .....................................................................................................................................................................117جدول 5-3: طول کمینه، بیشینه، متوسط و مساحت محدودهی گیرش چاه شماره یک.......118جدول 5-4: طول کمینه، بیشینه، متوسط و مساحت محدودهی گیرش چاه شماره دو.........118 فهرست شکلها عنوان صفحه شکل 2-1- تقسیم بندی مدلهای آب زیرزمینی .............................................................................. 8شکل 2-2- نحوه تغییرات هد در ستون ماسه در آزمایش دارسی ............................................... 28شکل 2-3- جریان ورودی و خروجی از المان حجمی ................................................................... 30شکل 2-4- هیدروگراف برای سلول i, j, k ....................................................................................... 39شکل 4-1- موقعیت دشت ایج فارس ................................................................................................ 62شکل 4-2- موقعیت چاههای بهره برداری و مشاهداتی منطقه مورد مطالعه .............................77شکل 4-3- شبکه بندی منطقه مورد مطالعه .................................................................................. 78شکل 4-4- نقشه DEM توپوگرافی سطح زمین منطقه مورد مطالعه ....................................... 79شکل 4-5- نقشه DEM رقوم ارتفاعی سنگ بستر منطقه مورد مطالعه .................................. 80شکل 4-6- هیدروگراف واحد دشت ایج طی سالهای 1375-1391 ...................................... 83شکل 4-7- نقشه سطح آب (متر) مهرماه 1388 دشت مورد مطالعه ....................................... 83شکل 4-8- میله رنگی نمایانگر خطای واسنجی ............................................................................ 90شکل 4-9- محل دفن زباله فرضی برای منطقه مورد مطالعه ................................................... 94شکل 5-1- نقشه زونبندي و مقادير به دست آمده براي پارامتر هدايت هيدروليکي (متر بر روز) طي کاليبراسيون مدل در شرايط پايدار (مهرماه 1388)..................................................... 98شکل 5-2- نقشه زون بندي و مقادير به دست آمده براي پارامتر تغذیه سطحی (متر بر ماه) طي کاليبراسيون مدل در شرايط پايدار (مهرماه 1388) ............................................................... 99شکل 5-3- نقشه توزيع سطح ايستابي به دست آمده طي کاليبراسيون مدل در شرايط پايدار (مهرماه 1388) ..................................................................................................................................... 100شکل 5-4- جهت حرکت جریان آب زیرزمینی در منطقه مورد مطالعه ................................. 101شکل 5-5- مقادير محاسباتي و مشاهداتي بار هيدروليکي در چاههاي مشاهداتي مختلف محدوده مطالعاتي در دوره واسنجي شرايط پايدار (مهرماه 1388) ........................................... 102شکل 5-6- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در دوره واسنجي شرايط پايدار (مهرماه 1388) ............................................................................................................ 102شکل 5-7- مقادير مشاهداتي بار هيدروليکي در مقابل مقادير باقيمانده (تفاوت مقادير محاسباتي و مشاهداتي) در دوره واسنجي شرايط پايدار (مهرماه 1388) ............................... 103شکل 5-8- نقشه زونبندي و مقادير به دست آمده براي پارامتر هدايت هيدروليکي (متر بر روز) طي کاليبراسيون مدل در شرايط ناپايدار (سالهای آبي 1390-1388) ...................... 104شکل 5-9- نقشه زون بندي و مقادير به دست آمده براي پارامتر آبدهي ويژه طي کاليبراسيون مدل در شرايط ناپايدار (سالهای آبي 1390-1388) ............................................................... 105شکل 5-10- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره یک در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) ................. 106شکل 5-11- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره دو در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) ................... 106شکل 5-12- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره سه در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) ................. 107شکل 5-13- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره چهار در طول دوره واسنجي شرايط ناپايدار (سالهای آبي 1390-1388) .............. 107شکل 5-14- حساسیت مدل نسبت به تغییرات پارامتر هدایت هیدرولیکی ........................ 110شکل 5-15- حساسیت مدل نسبت به تغییرات پارامتر آبدهی ویژه .................................... 111شکل 5-16- حساسیت مدل نسبت به تغییرات پارامتر تغذیه سطحی ................................. 112شکل 5-17- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره یک در طول دوره صحت سنجی (سال آبي 1391-1390) ........................................ 113شکل 5-18- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره دو در طول دوره صحت سنجی (سال آبي 1391-1390) ........................................... 113شکل 5-19- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره سه در طول دوره صحت سنجی (سال آبي 1391-1390) .......................................... 114شکل 5-20- مقادير محاسباتي در مقابل مقادير مشاهداتي بار هيدروليکي در چاه مشاهدهاي شماره چهار در طول دوره صحت سنجی (سال آبي 1391-1390) ........................................ 114شکل 5-21- محدوده گیرش چاههای با اهمیت منطقه مورد مطالعه پس از گذشت بیست سال ....................................................................................................................................................................116شکل 5-22- ناحیه گیرش چاه آب شرب شهر ایج پس از گذشت ده سال ............................ 117شکل 5-23- ناحیه گیرش دو چاه نمونه در دشت ایج ............................................................... 118شکل 5-24- مکانیابی حرکت ذرات شیرابه در آب زیرزمینی ناشی از وجود لندفیل ......... 119شکل 5-25- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از یک سال، مکانيسم انتقال جرم: پخش، انتشار) .................................................................................................................................................... 121شکل 5-26- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم:پخش، انتشار) .................................................................................................................................................... 122شکل 5-27- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم:پخش، انتشار) ................................................................................................................................................... 122شکل 5-28- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم:پخش، انتشار) ............................................................................................................................ 123شکل 5-29- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از یک سال، مکانيسم انتقال جرمپخش، انتشار و تأخیر) ..................................................................................................................................... 123شکل 05-30- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ..................................................................................................................................... 124شکل 5-31- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ..................................................................................................................................... 124شکل 5-32- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ زیاد (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) .............................................................................................................. 125شکل 5-33- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از یک سال، مکانيسم انتقال جرم: پخش و انتشار) ..................................................................................................................................................... 125شکل 5-34- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم: پخش و انتشار) ..................................................................................................................................................... 126شکل 5-35- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم: پخش و انتشار) ................................................................................................................................................. 126شکل 5-36- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم: پخش و انتشار) ........................................................................................................................... 127شکل 5-37- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از یک سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) .................................................................................................................................... 127شکل 5-38- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پنج سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ...................................................................................................................................... 128شکل 5-39- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از پانزده سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) ..................................................................................................................................... 128شکل 5-40- چگونگي توزيع ابر آلودگي ناشي از شيرابه آلاينده محل دفن زباله فرضي در آب زيرزميني در شرایط نرخ نفوذ کم (مدت زمان: پس از بیست و هشت سال، مکانيسم انتقال جرم: پخش، انتشار و تأخیر) .............................................................................................................. 129شکل 5-41- توزیع ابر آلودگی پس از گذشت بیست و هشت سال در شرایط نرخ نفوذ زیاد و مکانیسم انتقال پخش و انتشار به همراه موقعیت چاه شماره شش ........................................... 130شکل شماره 5-42- موقعیت سلولهای دلخواه مورد نظر جهت بررسی جزئی تر نتایج تغییر غلظت شیرابه در طول زمان ............................................................................................................... 131شکل 5-43- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ بالا، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 132شکل 5-44- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ بالا، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 133شکل 5-45- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ کم، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 134شکل 5-46- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ کم، غلظت اولیه آلودگی 20000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 134شکل 5-47- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ بالا، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 135شکل 5-48- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ بالا، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 136شکل 5-49- نمودار غلظت-زمان سلول شماره یک در شرایط نفوذ کم، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ............................................................................ 136شکل 5-50- نمودار غلظت-زمان سلول شماره دو در شرایط نفوذ کم، غلظت اولیه آلودگی 30000ppm و دو نوع مکانیسم انتقال آلودگی ........................................................................... 137 آب مایه حیات و فراوانترین ماده مرکب برروی سطح کره زمین و بستر اولیه حیات به شکلی که امروزه میشناسیم، میباشد. بیش از ۷۰٪ سطح کره زمین را آب پوشانده است (نزدیک به ۳۶۰ میلیون از ۵۱۰ میلیون کیلومتر مربع)؛ با وجود این حجم عظیم آب، تنها 2% از آبهای کره زمین شیرین و قابل شرب است و باقی آن به علت محلول بودن انواع نمکها خصوصاً نمک طعام غیر قابل استفاده است. از همین 2% آب شیرین، بیش از ۹۰% به صورت منجمد در دو قطب زمین و دور از دسترس بشر واقع شده است (Davie, 2002). به علاوه، منابع آب شیرین به طور یکنواخت در سطح زمین پراکنده نشدهاند. در حال حاضر، 60% کل منابع آب شیرین در 9 کشور جهان وجود دارد؛ در مقابل حدود 80 کشور با کمبود آب مواجهاند که برخی از آنها تقریباً به هیچ منبع آب شیرین قابل توجهی دسترسی ندارند (http://fa.wikipedia.org/wiki/%D8%A2%D8%A8#cite_note-autogenerated4-9). طبق آمار برنامه عمران سازمان ملل متحد در سال 2006:بنابراین، نحوه و چگونگی استفاده از منابع آب موجود در اکوسیستم طبیعی در سالهای اخیر موضوع بحث و بررسی علوم مختلف بوده است. امروزه مدیریت منابع آب[1] با در نظر گرفتن آنالیز همزمان سیستمهای آبی، اعم از سطحی و زیرزمینی و بررسی اثرات متقابل این منابع بر یکدیگر که از آن با نام مدیریت تلفیقی[2] منابع یاد میشود، به عنوان روش مناسب کارآمد مدیریتی در زمینه شناخت و برنامه ریزی هیدروسیستمها مطرح میشود. 1-2- اهمیت موضوع براساس آمار آژانس اطلاعات مرکزی (CIA[3]) در سال 2011 میلادی، کشور ایران با مساحت 1648195 کیلومتر مربع در رتبهی هجدهم در لیست وسیعترین کشورهای دنیا قرار دارد. براساس این آمار 3/0% خشکیهای کرهی زمین در محدودهی مرزهای ایران است. این در حالی است که تنها 03/0% آبهای موجود به ایران تعلق دارد. اهمیت این موضوع هنگامی بیشتر میشود که بدانیم 09/1% جمعیت جهان در ایران قرار دارد. از طرفی در آمار منتشر شده بین سالهای 1930 تا 1960 میلادی، کشور ما با 230 میلیمتر میانگین بارندگی سالانه در رتبه هشتاد و چهارم در میان کل کشورها است. دردناکتر اینکه میانگین بارندگی سالانه در میان کل کشورها 963 میلیمتر است(CIA, 2011).متأسفانه در سالهای اخیر برداشتهای بیش از ظرفیت از منابع آبهای زیرزمینی باعث ایجاد حالت بحرانی در اکثر دشتهای کشور ایران شدهاست. از طرف دیگر همزمانی آن با دورههای خشکسالی که دورههای تناوب آن بتدریج کوتاه شده است، مشکلاتی را برای توسعه بخشهای کشاورزی، اقتصادی و حتی تامین آب شرب به وجود آورده است. بدین ترتیب سطح آبهای زیرزمینی در منطقه روز به روز افت کرده و خشکیدن قنوات و چشمهها، کاهش آبدهی چاهها، فرونشست زمین، توسعه بیابان و در نهایت از دست رفتن سرمایه گذاریهای انجام شده در بخش کشاورزی، صنعت و مهاجرت را به همراه خواهد داشت (اکبرپور و همکاران، 1389).با نگاهی به مطالب ذکر شده به نظر میرسد اگر شرایط به همین روند ادامه یابد و مدیریت صحیح در زمان مناسب اعمال نشود، سرانجام به جایی خواهد رسید که آب مناسبی برای استخراج وجود نخواهد داشت و کشور با بحران شدید منابع آبی در قسمتهای مختلف روبرو خواهد شد.همچنین افزایش برداشت از آب زیرزمینی، ورود انواع پسابهای صنعتی، آب برگشتی کشاورزی آلوده به کودهای فسفاته و نیتراته، نشت از مخازن و خطوط فراوردههای نفتی و ... سبب شده است که کیفیت این منابع رو به نامناسب شدن پیش رود. اگر روند آلودگی، پخش و انتشار آن در آبهای زیرزمینی که سهم عمدهای در مصارف شرب و کشاورزی دارد مدیریت و کنترل نشود مخاطرات جبران ناپذیری برای سلامت بشر و محیط زیست خواهد داشت.