تجزیه و تحلیل پدیده فرونشست زمین در دشت نجف آباد با شبیه سازی عددی آب زیرزمینی

1- بررسی جهت عمومی جریان آب زیرزمینی منطقه بر اساس منحنی تراز سطح آب زیرزمینی نشان می دهد که آبخوان مورد مطالعه با آبخوان های اطراف از طرف شرق، غرب و جنوب در ارتباط است و فقط در ناحیه ای از شمال آبخوان است که ارتباط آبخوان با آبخوان های اطراف قطع شده است. 2- روند تقریباً یکسان هیدروگراف آبخوان و حجم آب رودخانه زاینده رود و منحنی تراز سطح آب زیرزمینی نشان می دهد که رودخانه زاینده رود نقش تغذیه کننده بر آبخوان دارد. 3- بررسی هیدروگراف ماهانه، سالانه، میانگین متحرک سه ماهه و پنج ماهه پیزومترهای منطقه و تفکیک پیزومترهای هم رفتار، نشان می دهد که منطقه به 4 زون تقسیم می شود. 4- با واسنجی مدل در حالت پایدار و ناپایدار، ضرایب هیدرودینامیک، شامل هدایت هیدرولیکی(در حالت پایدار) و آبدهی ویژه (در حالت ناپایدار) محاسبه گردید. پهنه بندی آبخوان بر اساس ضرایب هیدرودینامیک به دست آمده از واسنجی مدل، با پهنه بندی بر اساس هیدروگراف پیزومترها ، انطباق قابل قبولی دارد. 5- با تغییر پارامترها در مرحله واسنجی مدل جریان در حالت پایدار و ناپایداراین نتیجه اتخاذ گردید که مدل نسبت به تغییرات هدایت هیدرولیکی و سپس به پمپاژ از آبخوان حساسیت دارد. 6- با صحت سنجی مدل برای 6 سال آینده، بیلان دشت برای این بازه زمانی به صورت زیر به دست آمد: مولفه های بیلان حجم (میلیون متر مکعب) مربوط به سال 82-81 حجم (میلیون متر مکعب) مربوط به سال 83-82 حجم (میلیون متر مکعب) مربوط به سال 84-83 حجم (میلیون متر مکعب) مربوط به سال 84-85 حجم (میلیون متر مکعب) مربوط به سال 85-86 حجم (میلیون متر مکعب) مربوط به سال 86-87 برداشت از آب زیرزمینی توسط پمپاژ 3053- 2766- 2849- 2771- 799- 1110- جریان های خروجی از مرزها 503- 445- 917- 885- 474- 903- جریان های ورودی زیرزمینی از مرزها 1084 2283 2589 2081 588 877 نفوذ سطحی و آب برگشتی 5/9 5/7 5/3 22 0 0 تغذیه از طریق رودخانه 511 480 450 525 399 235 تغییرات حجم مخزن 1951- 442- 772- 1026- 285- 901- 7- نتایج بالا بیان می کند که بیش ترین تغذیه آبخوان از طریق آبخوان های مجاور و سپس رودخانه انجام شده و بیش ترین تخلیه نیز توسط چاه های بهره برداری بوده است. 8- با استفاده از تصاویر ماهواره ای فرونشستی حداکثر برابر 41 سانتی متر را در یک دوره 6 ساله و 30 سانتی متر را در یک دوره 2 ساله بر آورد شد. 9- با استفاده از تصاویر ماهواره ای سازمان فضایی اروپا، ضریب ذخیره ارتجاعی و ضریب ذخیره غیرارتجاعی واسنجی گردید، پهنه بندی آبخوان براساس این ضرایب نیز با پهنه بندی آبخوان براساس ضرایب هیدرودینامیک و پیزومترهای هم رفتار مطابقت دارد. 10- مقدار میانگین فرونشست برآورد شده توسط مدل برای دو سال واسنجی 17 سانتی متر بوده است. این مقدار فرونشست با توجه به کسری حجم آب در حدود 1186 میلیون متر مکعب در این دو سال قابل توجیه است. 11- با اجرای مدل برای 2 سال پیش بینی شده با شرایط مدل واسنجی، به طور متوسط فرونشستی معادل 25 سانتی متر برای منطقه به دست آمد. 12- هم چنین با اجرای مدل برای چنین شرایطی، فقط با 50 از پمپاژ کنونی به طور متوسط فرونشستی معادل 04/9 سانتی متر خواهیم داشت.