ارزیابی حساسیت فرونشست زمین در آبخوان دشت سلماس با استفاده از منطق فازی

مقدمهدشت سلماس، یکی از بحرانی‌‌ترین دشت‌‌های کشور از نظر فرونشست است. به‌طورکلی، عوامل مختلفی سبب فرونشست زمین می‌شوند، ولی در بسیاری مناطق برداشت بی‌رویه آب از سفره‌های زیرزمینی باعث فرونشست زمین می‌شود. افزایش روزافزون بهره‌برداری از آب‌های زیرزمینی به‌ویژه در حوضه‌هایی که با نهشته‌های آبرفتی، دریایی کم‌عمق یا دریاچه‌ای تحکیم نیافته انباشته شده‌اند، به نشست و یا فروریزش زمین منجر می‌شود. با برداشت بیش از حد آب زیرزمینی، سطح ایستابی لایه آبدار پایین رفته و فشار هیدرواستاتیک کاهش می‌یابد که درنتیجه امکان نشست زمین را به‌صورت تدریجی فراهم می‌نماید. فرونشست ماندگار دشت‌ها عمدتاً در اثر همین عامل یعنی برداشت بیش از حد مجاز آب و فشرده‌شدن لایه‌های رسی و سیلتی میان آبخوان صورت می‌گیرد. در این صورت، حتی اگر سطح ایستابی مجدداً بالاتر آید، امکان بازگشت زمین به تراز اولیه خود وجود ندارد. مواد و روش‌‌هادر این پژوهش، حساسیت فرونشست دشت سلماس با بهره‌‌گیری از لایه‌‌های عوامل موثر در فرونشست با نرم‌افزار arcgis و منطق فازی بررسی شد. در مرحله اول، آمار و اطلاعات برخی عوامل که موجب ایجاد فرونشست می‌‌شوند، شامل افت سطح آب زیرزمینی، دبی بهره‌‌برداری از چاه‏ها، ضریب ذخیره آبخوان، ضریب انتقال‌‌پذیری، بارش، نقشه dem، بافت خاک و عمق سنگ بستر، جمع‌‌آوری و نقشه رستری هر کدام از این عوامل در سطح آبخوان تهیه شد. در مرحله بعد، فازی‌‌سازی لایه‌‌ها با استفاده از توابع عضویت فازی، با توجه به تاثیر کاهش یا افزایش هر کدام از این عوامل در فرونشست زمین انجام گرفت. سپس، نقشه‌‌ها با عملگرهای فازی (gamma، or، and، sum، product) با هم تلفیق شده و نقشه واحد حساسیت فرونشست آبخوان به‌دست آمد. در نهایت، برای انتخاب بهترین عملگر ترکیب، نتایج را با داده‌‌های مشاهده‌‌ای صحرایی مقایسه شد و با شاخص عملکرد منحنی roc کنترل و صحت‌‌سنجی صورت گرفت. نتایج و بحثنتایج نشان داد که عملگر or با auc= 0.693 کمترین انطباق را با فرونشست رخ‌داده در منطقه دارد. عملگرهای گاما با سطح زیر منحنی بالای 70 درصد، بیشترین همپوشانی یا انطباق را با فرونشست مشاهده‌‌ای در دشت داشتند و در این پژوهش عملگر گامای 0.9 با  0.805=auc به‌عنوان بهترین عملگر فازی انتخاب شد. نتایج نشان می‌‌دهد که قسمت شرق آبخوان از نظر فرونشست بحرانی است. از کل مساحت دشت سلماس، 25 درصد معادل 93 کیلومترمربع دارای فرونشست با حساسیت خیلی بالا است. نتیجه‌‌گیریبا توجه به نتایج به‌دست آمده می‌توان بیان کرد که علی‌رغم اینکه مقدار auc اپراتور فازی sum بالاتر است؛ ولی در نقشه فازی گاما با مقدار 0.9 بیشترین انطباق را با واقعیت زمینی دارد، هرچند مقدار auc کمتری دارد. ذکر این نکته ضروری است، عملگرهای حداقل‌سازی and و product با ایجاد منطقه با حساسیت کم و عملگرهای حداکثرسازی or و sum با حداکثرسازی منطقه، نمی‌‌توانند کارایی مطلوبی در تهیه نقشه حساسیت فرونشست داشته باشند و در اینجا فقط برای نمایش‌‌دادن ضعف کارایی اپراتورهای فازی در حداکثر یا حداقل‌سازی حساسیت فرونشست زمین، استفاده شده است.

assessment of land subsidence susceptibility in the salmas plain aquifer northwestern lake urmia, utilizing fuzzy logic

introductionsalmas plain represents one of the most critical areas in the country experiencing subsidence. in general, various factors cause land subsidence, but in many areas, the excessive extraction of ground water from aquifers causes land subsidence. the increasing use of ground water, especially in the sites that are accumulated with alluvial deposits, shallow sea or unconsolidated lake, leads to subsidence or collapse of the land. with the excessive extraction of ground water, the water level of the aquifer decreases and the hydrostatic pressure decreases, which makes it possible for the land to subside gradually. subsidence in plains mostly occurs due to this factor, namely excessive groundwater extraction and compaction of clay and silt layers between aquifers. in this case, even if the water table level rises again, the land cannot return to its original level. materials and methodsin this study, the susceptibility of land subsidence in salmas plain was investigated using layers of influential factors in subsidence with arcgis software and fuzzy logic. in the first stage, statistical information on some factors causing subsidence, including groundwater level decline, well extraction rate, aquifer storage coefficient, transmissivity coefficient, precipitation, dem map, soil texture, and bedrock depth, was collected and raster maps of each of these factors at the aquifer level were prepared. in the next stage, fuzzy layering was performed using fuzzy membership functions based on the impact of decreasing or increasing each of these factors on land subsidence. subsequently, the maps were combined using fuzzy operators (gamma or, and, sum, product) to obtain a unified map of aquifer subsidence susceptibility. finally, to select the best combination of operators, the results were compared and evaluated with field observation data and the roc curve performance index. results and discussionthe results showed that the or operator had the lowest conformity with observed subsidence in the area with an auc of 0.693. gamma operators with an auc above 70% had the highest overlap or conformity with observed subsidence in the plain. in this study, the gamma 0.9 operator was selected as the best fuzzy operator with an auc of 0.805. the results indicate that the eastern part of the aquifer is critical in terms of subsidence. approximately 25% of the total area of salmas plain, equivalent to 93 square kilometers, has subsidence with very high susceptibility. conclusionbased on the results obtained, it can be said that although the auc value of the fuzzy operator sum is higher, the gamma operator with a value of 0.9 has the highest conformity with the ground reality on the fuzzy map, even though it has a lower auc value. it is essential to mention that the minimum operator and and product create a region with low susceptibility, while the maximum operator or and sum maximize the susceptible area. they cannot achieve satisfactory performance in preparing a subsidence susceptibility map. here, they have only been used to demonstrate the inefficiency of fuzzy operators in maximizing or minimizing subsidence susceptibility.