Abstract:
پیشینه و هدف در سال های اخیر به دلیل تغییرات اقلیمی و خشکسالی و همچنین فقدان نظارت در حفر چاه های کشاورزی بسیاری از سفره های زیرزمینی کشور مورد برداشت غیراصولی قرارگرفته اند که این امر موجب منفی شدن بیلان این آبخوان ها شده است تا جایی که بنا بر گزارش سازمان زمینشناسی بیشتر دشت های کشور دچار افت سطح آب زیرزمینی شده اند. یکی از مخاطراتی که به دلیل عملکرد تدریجی کمتر مورد توجه قرار میگیرد پدیده فرونشست است که در سال های اخیر به علت افزایش بهره برداری از سفره های زیرزمینی مهم در دشت های بارز کشور به صورت مشکلی فراگیر نمایان شده است. در این تحقیق تلاش شده است که امکان وقوع پدیده فرونشست و خطرات احتمالی آن به عنوان یک تهدید در انجام پروژه های انسانی و همچنین سکونتگاه های روستایی در محدوده دشت اردبیل مورد بررسی قرار گیرد. مواد و روش ها منطقه مورد مطالعه دشت اردبیل بین عرض جغرافیای شمالی 38 درجه و 05 دقیقه تا 38 درجه و 03 دقیقه و طول جغرافیای شرقی 48 درجه و 15 دقیقه تا 48 درجه و 35 دقیقه در بخش شمال غربی در مرکز استان اردبیل قرارگرفته است. برای بررسی وضعیت آب زیرزمینی دشت، از داده های 38 مشاهده ای که از سازمان آب منطقه ای اردبیل تهیه شده اند و در سطح دشت قرارگرفته اند استفاده شده است. نخست با استفاده از آمار 30 ساله 65 چاه مشاهده ای و سامانه اطلاعات جغرافیایی، نقشه پتانسیل افت آب برای منطقه تهیه شد. سپس با استفاده از مدل دماتل فازی اقدام به جمع آوری و مدلسازی نظرات کارشناسان شد. این روش ازجمله روش های مفهومی برای ساختار سازی مسئله های تصمیم گیری به شمار میرود. فن دماتل بر پایه نظریه گراف ایجاد شده است و در این روش ما می توانیم معیارها را به دو گروه معیارهای علت و معیارهای معلول تقسیم بندی کنیم تا شناخت بیشتری از روابط بین علت و معلول به دست آوریم و درنهایت بتوانیم یک شبکه ای از روابط متقابل را ایجاد نماییم. درنهایت بعد از ایجاد ماتریس روابط کلی و با توجه به اندازه حد آستانه تعریف شده، ماتریس روابط نهایی را ایجاد می کنیم که در این ماتریس عدد صفر به معنی عدم وجود رابطه و عدد یک به معنی وجود رابطه میان دو معیار هست. با استفاده از ماتریس روابط نهایی اقدام به نظرسنجی از خبرگان در مورد اندازه تأثیر عوامل بر یکدیگر با توجه به وابستگی آنها به هم مینماییم. بعدازآنکه دادهها از سازمانهای مربوط تهیه شد اقدام به تشکیل پایگاه داده برای اطلاعات موجود در محیط نرم افزار ArcCatalog کرده و سپس نقشه های مربوط به این داده ها را در نرمافزار ArcGIS ایجاد شد. بعد از اینکه وزن لایه های مختلف با استفاده از روش تحلیل شبکه ای فازی به دست آمد، وارد نرم افزار ArcGIS شده و وزن های هر زیر معیار را در نقشه ای که برای هر لایه موردنظر ایجاد کرده بودیم ضرب کرده و درنهایت نقشه ها را باهم جمع کردیم تا یک نقشه نهایی به دست آمد. نقشه نهایی نشان دهنده مناطقی از دشت اردبیل است که از نظر خطر فرونشست طبقه بندی شده و در پنج طبقه از لحاظ وضعیت خطر با طیف رنگی نشان داده شد. بعد ناحیه ای که افت شدید آب داشت انتخاب و با نقشه پراکندگی چاه های عمیق مقایسه شد. در مرحله آخر با استفاده از مدل های پیشرفته و ترکیبی فازی و تحلیل شبکه ای در محیط نرم افزار سامانه اطلاعات جغرافیایی هر یک از لایه های حساسیت رسوبات و افت سطح آب تعیین عضویت و با استفاده از همپوشانی خطی فازی نقشه حساسیت منطقه به پدیده فرونشست در پنج کلاس حساسیت بسیار بالا، حساسیت بالا، حساسیت متوسط، حساسیت پایین و حساسیت بسیار پایین تهیه شد. برای تهیه نقشه نهایی وضعیت خطر فرونشست در محدوده دشت اردبیل ابتدا وزن های به دست آمده برای هر زیر فاکتور را در نقشه های فازی شده آن زیر فاکتور ضرب کرده و سپس با استفاده از ابزار RasterCalculator این نقشه های وزندار تجمیع شدند. در نقشه فازی نهایی ارزیابی وضعیت خطر فرونشست دشت اردبیل با طیف های رنگی از آبی که پایین ترین حد را نشان میدهد به سمت قرمز که بیشترین مقدار را نشان می دهد، در تغییر است. رنگ آبی مناطق با خطر پایین و رنگ قرمز مناطق با خطر بالا را از نظر خطر فرونشست در محدوده دشت اردبیل و روستاهای واقع در این محدوده نشان می دهد.نتایج و بحث پس از به دست آوردن نقشه اولیه ارزیابی وضعیت خطر فرونشست و مناطق در معرض خطر فرونشست دشت اردبیل، نقشه نهایی بررسی وضعیت دشت اردبیل از لحاظ خطر فرونشست زمین، با توجه به وزن ها و لایه های به دست آمده تهیه شده است. از آنجایی که تمامی نقشه های لایه های مبناء در پنج طبقه، باز طبقه بندی شده و وزن متناسب با هر طبقه مطابق با وضعیت لایه ها داده شده بود، نقشه نهایی نیز در پنج طبقه، طبقه بندی و وزن دهی شد که بنا بر نظر کارشناسان و اساتید به شرح زیر به دست آمده است؛ 1) مناطق با خطر پایین، 2) مناطق با خطر متوسط، 3) خطر بالا، 4) مناطق در معرض آسیب، 5) مناطق بحرانی. درنهایت با استفاده از نقشه نهایی ارزیابی وضعیت خطر فرونشست دشت اردبیل و نیز رتبه بندی به دست آمده از نظرات کارشناسان ذیربط، نقشه نهایی تحلیل وضعیت خطر فرونشست دشت اردبیل تهیه شد. همچنین نقشه وضعیت چاه های عمیق در محدوده دشت اردبیل و پراکندگی آن در سطح دهستان ها را نشان می دهد، ملاحظه می شود که بیشترین پراکندگی و تمرکز چاه های عمیق در سمت شرق دشت در محدوده دهستانهای ویلکیج مرکزی و فولادلوی شمالی قرارگرفته است که این وضعیت موقعیت پراکندگی چاه های عمیق به خوبی وضعیت متناسب بودن پراکندگی چاه های عمیق با مناطق در معرض خطر فرونشست را نشان می دهد.نتیجه گیری دهستان های ویلکیج مرکزی، فولادلوی شمالی و دهستان شرقی بیشترین سطح آسیب پذیری را ازلحاظ وضعیت خطر فرونشست دارند. وضعیت بحرانی خطر فرونشست در این سه دهستان بیشترین مقدار را دارد. همچنین دهستان کلخوران و آغبلاغ آقاجان خان از وضعیت متوسط به سمت وضعیت در معرض آسیب پیش می روند که نیازمند دقت بیشتری در مدیریت و برنامه ریزی منابع آب این دهستانها هست. همچنین یک رابطه قوی بین نحوه پراکندگی چاه های عمیق در سطح دشت اردبیل و مناطق در معرض خطر بالای فرونشست وجود دارد. با توجه به نتایج به دست آمده، زیر معیار وضعیت آب های زیرزمینی با وزن 0.38 بیشترین تاثیر را در وضعیت خطر فرونشست دشت اردبیل ایجاد کرده است. این وزن با بررسی سایر لایه های مرتبط با وضعیت آب زیرزمینی و لایه پراکندگی جمعیت، نشان دهنده تاثیر زیاد این زیر معیار است.
Background and Objective
In recent years, due to climate change and drought, as well as the lack of supervision in digging agricultural wells, many of the country's aquifers have been harvested improperly, which has led to a negative balance of these aquifers to the extent that, according to the Geological Survey, most of the country's plains have experienced a drop in groundwater levels. Today, the study, planning and planning to reduce the risks of natural hazards is one of the main issues of officials and planners of countries. One of the hazards that are less considered due to gradual performance is the phenomenon of subsidence, which in recent years due to increased use of important aquifers in the plains of the country has become a pervasive problem. In this research, an attempt has been made to investigate the possibility of subsidence and its possible dangers as a threat to human projects as well as rural settlements in the Ardabil plain.tp
Materials and Methods
The study area of Ardabil plain is located between the latitude of thirty-eight degrees and thirty minutes north latitude to thirty-eight degrees and thirty minutes and the longitude of east geography forty-eight degrees and fifteen minutes to forty-eight degrees and thirty-five minutes in the northwestern part of Ardabil province. To investigate the groundwater status of the plain, data from 38 observation wells prepared by the Ardabil Regional Water Organization and located at the plain level have been used. First, using 30-year statistics of 65 observation wells and GIS, a water potential drop map for the region was prepared. Then, using fuzzy Dematel model, experts' opinions were collected and modeled. This method is one of the conceptual methods for structuring decision problems. The Dematel technique is based on graph theory, and in this way we can divide the criteria into two groups of cause and effect criteria to gain a better understanding of the cause-and-effect relationship and finally be able to create a network of interrelationships. Finally, after creating the general relationship matrix and according to the defined threshold size, we create the final relationship matrix in which the number zero means no relationship and the number one means the relationship between two criteria. Using the final relations matrix, we conduct a survey of experts on the extent to which factors affect each other with respect to their interdependence. After the data was obtained from the relevant organizations, a database was created for the information in the ArcCatalog software environment, and then maps related to this data were created in ArcGIS software. After the weights of the different layers were obtained using the fuzzy network analysis method, they entered the ArcGIS software and multiplied the weights of each sub-criterion in the map we created for each layer and finally gathered the maps together to get a final map Came. The final map shows the areas of Ardabil plain that are classified in terms of subsidence risk and in five categories in terms of danger status were shown with color spectrum. Then, the area with severe water loss was selected and compared with the scattering map of deep wells. In the last step, using advanced and fuzzy hybrid models and network analysis in the software environment of GIS, each of the layers of sediment sensitivity and water level drop membership is determined and using fuzzy linear overlap, the area sensitivity map to subsidence in five classes of very sensitivity High, high sensitivity, medium sensitivity, low sensitivity and very low sensitivity were prepared. To prepare the final map of subsidence risk status in Ardabil plain, first the obtained weights for each sub-factor were multiplied in the fuzzy maps of that sub-factor and then these weighted maps were aggregated using the Raster Calculator tool. The final fuzzy map of subsidence risk assessment of Ardabil plain is changing with the color spectrum changes from blue, which represents the lowest, to red, which represents the highest. The blue color indicates low risk areas and the red color indicates high risk areas in terms of subsidence risk in Ardabil plain and villages located in this area.
Results and Discussion
After obtaining the initial plan to assess the subsidence risk status and areas at risk of subsidence in Ardabil plain, the final map to assess the status of Ardabil plain in terms of land subsidence risk has been prepared according to the weights and layers obtained. Since all the base layer maps were reclassified into five layers and the weight corresponding to each layer was given according to the condition of the layers, the final map was classified and weighted into five layers, which according to experts and professors is as follows; 1) Low risk areas, 2) Medium risk areas, 3) High risk, 4) Damaged areas, 5) Critical areas. Finally, using the final map of Ardabil plain subsidence risk assessment, as well as the ranking obtained from the opinions of relevant experts, the final map of Ardabil plain subsidence risk analysis was prepared. Also, the map of deep wells in Ardabil plain and its distribution in rural areas, it can be seen that the highest distribution and concentration of deep wells in the eastern part of the plain is in Wilkij e Markazi and Fooladloo e Shomali villages. This situation shows the scattering position of deep wells showing the proportionality of the scattering of deep wells in areas at risk of subsidence.
Conclusion
Wilkij e Markazi, Fooladloo e Shomali, and Fooladloo e Sharghi have the highest levels of vulnerability in terms of subsidence risk status. The critical situation of landslide risk is the highest in these three villages. Also, Kalkhoran and Aghbalagh Aqajan Khan villages are moving from a moderate to a vulnerable situation, which requires more care in managing and planning the water resources of these villages. Also, there is a strong relationship between the distribution of deep wells in the Ardabil plain and areas at high risk of subsidence. Also, according to the results obtained the groundwater status sub-criterion with a weight of 0.38 has the greatest impact on the subsidence risk of Ardabil plain. This weight shows the high impact of this sub-criterion by examining other layers related to groundwater status and population dispersion layer.