Subject Areas : geographical and environmental
1 - Assistant Professor of Islamic Azad University, Larestan branch
Keywords:
Abstract :
منابع
1 - آبیل ، ا . طاووسی ، ت . خسروی ، م ( 1398 ) به بررسی و تحلیل مناطق بالقوه در معرض مخاطرات سیلاب شهری زاهدان ، جغرافیا و توسعه ، شماره 54 ، ص 106 – 91 .
2 - امیرانی ، محمد هادی ( 1371 ) ، اثرات و پیامدهای بلایای طبیعی بر منابع طبیعی ، ماهنامه جهاد ، شماره 230 و 231 .
3-تاجبخش فخرآبادی ، سید محمد، چزگی ،جواد ،1402؛ اولویتبندی سیلخیزی زیرحوزههای آبخیز شمال دشت بیرجند با استفاده از عوامل مورفومتری و مدل ویکور،مدل سازی و مدیریت آب و خاک دوره 3، شماره 3، مهر
4 - شریفی پیچون ، م . امیدوار ، ک . متذکر ، ک ( 1398 ) استفاده از روش تحلیل خوشه ای و رگرسیون چند متغیره در ارزیابی پتانسیل سنجی سیلاب با تاکید بر پارامترهای هیدروژئومورفولوژیکی (مورد مطالعه: حوضه آبخیز رودخانه مارون) ، مخاطرات محیط طبیعی ، شماره 21 ، ص 92 – 75 .
5 - حسام ، ر . ضرابی ، ا . تقوایی ، م ( 1398 ) پتانسیل سنجی خطر سیلاب شهری با رویکرد توسعه شهری ایمن شهر گنبد کاووس، مدیریت مخاطرات طبیعی ، دوره 6 ، شماره 1 ، ص 32 – 17 .
6 - حاتمی نژاد ، ح . آتش افروز ، ن . آروین ، م ( 1396 ) پهنهبندی خطر سیل با استفاده از تحلیل چندمعیاره و GIS مطالعه موردی: شهرستان ایذه ، دانش پیشگیری و مدیریت بحران ، دوره 7 ، شماره 2 ، ص 57 – 44 .
7 - دهقانی ، حسن ( 1394 ) پهنه بندی خطر سیلاب در شهر شیراز با استفاده از GIS ، پایان نامه کارشناسی ارشد رشته سنجش از دور ، دانشگاه آزاد لارستان .
8 - کریمی فیروزجانی ، م . نیسانی سامانی ، ن ( 1397 ) پهنه بندی ریسک پذیری سیل با استفاده از روس تصمیم گیری چند معیاره مبتنی بر استدلال شهودی دمپستر – شافر درحوزه آبخیز نکا ، پژوهش های جغرافیای طبیعی ، دوره 50 ، شماره 1 ، ص 194 – 77 .
9 - مجرد کاهانی ، م . خواهنده کارنما ا . زهره نژاداکبری راوری ، ز ( 1393 ) مدیریت بحران سیل شهری، اولین کنگره تخصصی مدیریت شهری و شوراهای شهر، ساری، مرکز همایشهای توسعه ایران .
10.Amoon P. H. (2015). “From hydraulic modelling to urban flood risk”, Procedia Engineering, 115, pp: 37-44.
Chinese Soil & Water Conservation28(1),11-20.
11.Morellin, S.; Segoni, S.; Manzo, G.; Ermini, L.; & Catani, F, (2018), “Urban planning, flood risk and public policy: the case of the Arno River, Firenze”, Italy, Applied Geography, 34, pp: 205-218.
12.Parker, A.; Mignot, E.; & Bazin, P. H. (2015). “From hydraulic modelling to urban flood risk”, Procedia Engineering, 115, pp: 37-44.
13.Plat,E.G (2010). Flood risk and flood management , Journal of hydrology 267 (1-2), PP. 2-11.
14.Stefan , D.S.,Lutza,M.A. (2011),Urban Flood Hazard Zoning in Tucumán Province, Argentina,Using GIS and Multicriteria Decision Analysis, Engineering Geology,111(1–4),90–98
15.Telford, G. (2012). A GIS-based spatial multi-criteria approach for flood risk assessment in the Dongting Lake Region, Hunan, Central China, Water resources management, 25(13): 3465-3484
مجله علوم جغرافيايي، دانشگاه آزاد اسلامي واحد مشهد، دوره 20، شماره 48، پاییز 1403، صص 120-99
پتانسیلسنجی خطر سیلاب شهری با رویکرد توسعه شهری ایمن با استفاده از مدل Knn ( مطالعه موردی : شهر لار)
دکتر محمد ابراهیم عفیفی: دانشیار گروه جغرافیا، واحد لارستان، دانشگاه آزاد اسلامی، لارستان،ایران(نویسنده مسئول)
Email: afifi.ebrahim6353@gmail.com
دريافت: 18/1/1402 پذيرش: 30/7/1402
چکیده
سیلابها، یکی از پدیده ها و مخاطرات طبیعی بسیار مهم و مخرب در سطح جهان محسوب میشوند. سیلاب ها نه تنها در نواحی برون شهری ،بلکه در نواحی شهری نیز خسارات فراوانی به بار می آورند . تغییرات انسانی ایجاد شده در الگوی طبیعی سیستمهای زهکشی در نواحی شهری وحاشیه شهرها عامل اصلی بروز سیلاب وخسارات وارده است. روند رو به رشد رخداد سیل در سال های اخیر حاکی از آن است که اکثر مناطق کشور نیز در معرض تهاجم سیلاب قرار دارند . هدف این نوشتار ، بررسی وضعیت سیل خیزی در شهر لار استان فارس است.برای این منظور پس از رقومی سازی داده های مورد نیاز در سیستم اطلاعات جغرافیایی ، لایه های نقشه ای مختلف شامل بارندگی ، زمین شناسی ، جاده ، رودخانه ، خاکشناسی و ... تولید شد.سپس ارزیابی و پهنه بندی سیل خیزی با استفاده از مدل knn به عمل آمد. نقشه نهایی خطر با طبقه بندی متنوع ترسیم شد. که بر اساس نتایج تحقیق قسمت های شمالی و مرکزی حریم رودخانه لار در معرض پتانسیل خطر بالا قرار دارد و قسمت های جنوبی در معرض با پتانسیل خطر پایین قرار داد .
کلمات کلیدی : سیلاب ، پتانسیل سنجی ،توسعه، خطر،مدل Knn
مقدمه
برخی این موضوع را مطرح می کنند که جهان نسبت به آن چه سابقاً و پیش از این بوده مخاطره انگیز تر شده است. اما به هر حال دوران و زمان ما خطرناک تر از گذشته نیست،بلکه تعادل بین مخاطرات و خطرات تغییر کرده است .ما در جهانی زندگی می کنیم که مخاطراتی که خود ایجاد شده است به همان اندازه یا بیش از آن چه که از دنیای خارج به ما تحمیل می شود، است.مخاطرات محیطی عبارتند از رویدادها یا شرایط فیزیکی که توانایی ایجاد مرگ و میر ها ،آسیب ها ،خسارات مالی ،خسارات زیر بنایی ،ضرر و زیان های مربوط به کشاورزی ،خسارت به محیط ،وقفه های تجاری و دیگر ضررها و خسارت ها را دارند (امیرانی ، 1371 ). بلایای طبیعی رویدادهای ناخواسته و اغلب ناگهانی هستند که باعث بوجود آمدن خسارات جانی ،مادی ،اقتصادی و یا محیط می شوند و فراتر از قدرت مقابله کشور یا جامعه تحت تاثیر آن ها هستند.در طی چند دهه گذشته آسیب ها و اثرات اقتصادی بلایای طبیعی و انسان ساخت با سرعت زیادی رو به افزایش بوده و همچنان ادامه دارد .دامنه این بلایا و ارتباط جهانی ما تا نقطه ای رشد کرده که در آن بلایا(طبیعی و انسان ساخت ) ممکن است پیامدهای منطقه ای ،ملی و حتی پیامدهای جهانی داشته باشد. نتایج تحقیقات نشان می دهد که جوامع جهانی نسبت به قبل در برابر بلایای طبیعی و انسان ساخت حساس تر شده اند . از بین مخاطرات طبیعی شاید سیلاب ها گستره و فراوانی زیادی در ایران داشته و همیشه ایران را با مشکلات فراوانی مواجه می سازند(عفیفی،1401) .سیلاب ها در ايران معمولاً،در حين وقوع توجه مسئولين را به خود جلب کرده است، و با قطع سیلاب ها به بوته فراموشي سپرده شده است. شايد همين موضوع يکي از دلايل عدم برنامه ريزي موثر در مسائلي نظير پيش بيني1،پايش و مقابله و تخفيف آثار زیانیار اين پديده شده است ( مجرد کاهانی، 1393 ) .
سیل یکی از پدیدههای پیچیده و مخرب طبیعی است که هر ساله خسارات فراوانی را به دنبال دارد. به منظور برنامهریزی جهت بهرهبرداری بهینه و کنترل پدیدههای طبیعی سیل که از مسائل موجود در کره زمین و از جمله کشور ما نیز می باشد، مطالعات و تحقیقات زیادی مورد نیاز است. در کشور ایران، در طی چند سال اخیر 1260 سیل با میانگین 30 سیل در سال رخ داده است، که این آمار حاکی از رشد وقوع سیل در سالهای اخیر به میزان 40% می باشد. آمارها حاکی از آن است که سیلاب چه از نظر تلفات جانی و چه از نظر تلفات مالی مقام اول را در میان حوادث دیگر داراست. این امر کشور ایران را به لحاظ سیل خیزی در رتبه هفتمین کشور دنیا قرار داده است. سیل علل مختلفی می تواند داشته باشد( عفیفی،1401)
- دهقانی ( 1394 ) در پایان نامه ارشد خود به پهنه بندی خطر سیلاب شهر شیراز با استفاده از GISپرداختند که نتایج حاصل از مدل SCS نشان می دهد که در اراضی بالا دست حوضه که دارای شیب بالا، اراضی لخت، خاکهای نفوذ ناپذیر پراکندگی بیشتری دارد ظریب نگهداشت خاک پایین بوده و در نتیجه میزان رواناب تولیدی نیز بیشتر است. با افزایش رواناب در سطح حوضه بالا دست حداکثر دبی نیز افزایش یافته و با این افزایش پهنه های سیلابی با وسعت بیشتری زیر اب فرو می روند .
- حاتمی نژاد و همکاران ( 1396 ) به پهنهبندی خطر سیل با استفاده از تحلیل چندمعیاره و GIS مطالعه موردی: شهرستان ایذه پرداختند که نقشه نهایی پهنهبندی سیلاب نشان داد که حدود 43 درصد از شهرستان در پهنه خطر خیلی زیاد ، 16 / 14 درصد در پهنه خطر زیاد ، 12 / 46 درصد در پهنه خطر متوسط، 10 / 91 درصد در پهنه خطر کم و 8 / 3 درصد خیلی کم از لحاظ سیلگیری قرار دارد. بهطورکلی حدود60 درصد از پهنههای شهرستان درخطر سیلگیری خیلی زیاد و زیاد است . طبق نقشه ارزیابی خسارت، بیشترین خسارت مربوط به مناطق شهر ایذه، سوسن شرقی و غربی و روستاهای دنباله رود جنوبی است. میتوان گفت مناطقی که در موقعیت کوهستانی واقه شدهاند، درخطر سیلگیری قرار دارند. زیرا تمام این مناطق در بین کوهها و یا محل تلاقی درهها و تنگهها واقع شدهاند. بهطورکلی میتوان گفت نقشههای ریز پهنهبندی خطر سیلاب، به مدیریت غیر سازهای سیل کمک میکند و به مدیران برنامهریزی این امکان را میدهد تا بخشهای امن، از نظر سیلخیزی را برای توسعه انتخاب کنند.کریمی فیروزجانی ( 1397 ) به بررسی و پهنهبندی ریسکپذیریسیل با استفاده ازروس تصمیم گیری چند معیاره مبتنی بر استدلال شهودی دمپستر – شافر درحوزه آبخیز نکا پرداختند که بر اساس نتایج بدست آمده از بررسی طبقات مختلف ریسک پذیری زیر حوزه پایاب نکا و سراب نکا به ترتیب بیشترین و کمترین مساحت مربوط به طبقه ریسک زیاد و خیلی زیاد را به خود اختصاص داده اند . حسام و همکاران ( 1398 ) به پتانسیل سنجی خطر سیلاب شهری با رویکرد توسعه شهری ایمن شهر گنبد کاووس پرداختند که نتایج حاصل از نقشه حاصل از خطر پذیری شهرستان گنبد کاووس که با استفاده از 9 پارامتر تعیین شد بیانگر قرار گیری این شهر در پهنه با پتانسیل خطر زیاد است که نتایج حاصل و اتفاقات سیل اخیر گنبد کاووس ناشی از ریسک بیرونی سیلاب است . پلات(2010 ) در تحقیقات خود با نام خطر سیل و مدیریت آن ، پهنه بندی خطر را به عنوان شیوه کاملاً مشخصی برای ساماندهی و مدیریت خطرهای ناشی از عوامی طبیعی ، زیست محیطی یا انسانی که از بین آنها سیل نیز بارزتر است ، تعریف کرد . استفان(2011) در تحقیقات خود به بررسی سیلهای به وقوع پیوسته در سالهای 1966و 1977 در حوزۀ آبخیز رودخانه سنگ زرد واقع در ایالت مونتانای آمریکا پرداخت . وی همچنین برای تعیین پهنههای سیل رخ داده در این 2 سال مشخصات 25 مقطع را در طول یک مسیر 18 کیلومتری برداشت و پهنههای سیل با دورههای بازگشت مختلف را تعیین کرد . تل فورد(2012) در تحقیقی با عنوان تجزیه و تحلیل سیلابها در کشور آمریکا به بررسی عوامل به وجود آورنده سیلابها از دیدگاه ژئومورفولوژیک پرداخته است .وی در روش خود با استفاده از سیستم اطلاعات جغرافیایی (GIS) و با مدلAHPبه مقوله فوق پرداخته و در پایان به این نتیجه رسید که عوامل ژئومورفولوژیک نقش بسزایی در پیدایش سیلابها در منطقه مورد نظر دارند.آمون( 2013 ) با بررسی حوزه های مختلف در کالیفرنیا از لحاظ سیلابی بودن با استفاده از GIS وRS به این نتیجه رسید که عوامل بسیاری در پیدایش سیلابها نقش دارند که میتوان سهم بیشتری را به عوامل ژئومورفولوزیک داد.وی با مدلی جامعه به بررسی مقوله فوق پرداخته است و به نتایجی مطلوب دیگری نیز دست یافته است .پارکر( 2014 ) در پژوهشی پس از بررسی و مطالعه سیلابهای کشور فرانسه به اثرات اقتصادی و اجتماعی آنها پرداخته است .وی با مطالعه بیش از صد سیلاب به این پژوهش پرداخته است .در روش خود با استفاده از مدلهای اقتصادی اجتماعی مطالعه سیلابها به نتایجی دست یافته که یکی از آنها مداخله انسان در طبیعت و تغییر اشکال ژئومورفولوژیکی میباشد و این تغییر باعث تبدیل شده و به جان بشر امروزی افتاده است . موریلس ( 2018 ) در تحقیقی پتانسیل خیزی رودخانه آرنو در ایتالیا را با استفاده از سنجش از دور و GIS مورد بررسی قرار دادند که نتایج تحقیق حاکی از آن است که مناطق شهری توسعه یافته تر بیشتر در معرض ریسک قرار دارند . پديده سیل و سیلاب طي دهه های گذشته در منطقه مورد مطالعه ( لارستان ) رخ داده است و آثار و پيامدهاي نامطلوب خود را نشان داده است . وضعيت بارندگي های شدید طی سالیان گذشته و محدوديت ناشی از آن و شرايط اقليمي منطقه مورد مطالعه اين واقعيت را نشان میدهد که بايد براي مقابله پديده سیلاب برنامه داشت و در زمان وقوع با آن آثار مقابله نمود.از اين رو شناخت وضعيت منطقه به لحاظ سیلاب و به خصوص پي بردن به رفتار خاص دورههایسیلابی آن به لحاظ وقوع، شدت و گسترش مکاني هر یک از سیلاب ها بدون شک ابزار بسيار سودمندي و کارآمدي در اختيار برنامه ریزان شهری خواهد داد که از اين طريق به اهدافي که در خصوص زیانهای ناشي از اين حالت منظور است ، نائل گرديد.برای شناسایی و پتانسیل سیلخیزی حوضههای آبریز سیستم اطلاعات جغرافیایی ابزار توانمندی هستند و بوسیله آن می توان تغییرات محیطی مانند،آلودگی آب، فرسایش و... پایش نمود . لذا در این تحقیق به بررسی پتانسیل سنجی خطر سیلاب شهری با رویکرد توسعه شهری ایمن در شهر لار پرداخته می شود که این تحقیق در صدد پاسخگویی به سوال آیا کاهش پوشش گیاهی و بارندگی از دلایل اصلی افزایش پتانسیل سیل خیزی شهر لار می باشد ؟ است .
موقعیت جغرافیایی منطقه مورد مطالعه
شهرستان لار در فاصله۳۳۰کیلومتری جنوب شرقی شهر شیراز و ۱۹۰ کیلومتری شمال غرب بندرعباس قرار دارد. اين شهرستان در جنوب استان فارس، در حدفاصل شهرستانهاي داراب ، جهرم ، فيروزآباد و استان بوشهر و هرمزگان واقع شده و مساحت آن 33/667 كيلومتر مربع است ، از نظر وسعت بزرگترين شهرستان استان فارس به شمار ميرود و ارتفاع آن از سطح دريا 915 متر ميباشد . اين شهرستان ، برابر سرشماري سال 1355 داراي 183921 نفر و مطابق سرشماري سال 1390 داراي 226879 نفر جمعيت است . از نظر تقسيمات كشوري شامل شش بخش : مركزي ، جويم ، خنج ، اوز ، گراش و بيرم ميباشد . درجه حرارت و بارندگي حكايت از خشكي حاكم بر اين شهر دارد ، بطوري كه مجموع تبخير ساليانه به بيش از سه هزار ميليمتر ميرسد كه حدود 16 برابر بارش متوسط ساليانه (حدود 188 ميليمتر) ميباشد . بنابراين در طراحي و ساخت شهرهاي مناطق گرم و خشك مانند لار ، ميبايست از تكنيكها و امكاناتي استفاده نمود كه نه تنها اثر اين عوامل نامساعد جوي را تشديد نكند ، بلكه از آثار مضّر آن نيز بكاهد ، شهرستان لارستان از '5،ْ27 تا '30، ْ28 عرض شمالي و '17،ْ50 تا '44 ، ْ55 شرقي وسعت دارد و در جنوب شرقي فارس قرار گرفته است . لارستان از شمال به زرين شهر و جهرم و فيروزآباد و از شرق و جنوب شرقي به بندرعباس و بندرلنگه واز جنوب غربي و غرب به لامرد محدود است . لارستان با وسعت 20964كيلومتر مربع وسيع ترين شهرستان استان فارس است و حدود 17 درصد از مساحت استان را شامل مي شود ( اسدیانی و همکاران ،1387) . مرکز شهرستان لارستان شهر لار است که در جنوب استان فارس و مرکز شهرستان لارستان که در فاصله ۳۳۰ کیلومتری واقع شده است.
شکل 1 تصویر ماهواره ای لارستان ( ترسیم کننده : نگارنده:1402 ) |
---|
– عوامل موثر بر سیلاب های شهر لار
بر اساس بررسي هاي انجام شده ، مهمترين عوامل مؤثر بر سيلابهاي حوضة آبريز لارستان عبارتند از :
الف . عوامل فيزيكي :
علت اصلي و قوع سيل ريزش بيش از اندازه نزولات جوي مي باشد هر چند براي به وقوع پيوستن سيل ، باران زياد در مدت كم و شدت زياد و تداوم آن مورد نياز مي باشد ، اما چگونگي ريزش به محيط هم بستگي دارد . به همين دليل بررسي شيب ، تراكم زهكشي ، زمان تمركز ، ارتفاع حوضه ، مساحت و محيط حوضه ، اشكال و فرآيندهاي ژئومورفيك ، خصوصيات زمين ساختي ، خصوصيات خاك و پوشش گياهي هركدام به نوبة خود ميتواند تأثيرگذاريرابرحجم سيلاب و خصوصيات سيل داشته باشند .
ب . جنس زمين :
در منطقة لارستان ، تأثير جنس زمين در واحدهاي هيدرولوژيك يكسان نمي باشد و هر كدام داراي خصوصيات متفاوتي هستند . وجود درزوشكاف در زمين هاي آهكي يا كارستيك كه خود سبب نفوذ آب به درون زمين مي شود ، از شدت سيلاب ها مي كاهد و آن را تعديل ميكند . و درزمين هاي مارني و رسي ، شرايط عكس آن حاكم است .
ج . پوشش گياهي :
وجود زمين هاي باير و يا با پوشش گياهي تنك در مقابل ريزش هاي جوي بازتاب ي متفات با زمين هايي كه داراي پوشش گياهي متراكم هستند ، دارند . از طرفي با افزايش تراكم پوشش گياهي از شدت روانابهاي سطحي كه نهايتا"منجر به ايجاد سيل ميشوند كاسته مي گردد . بيشترين وسعت حوضه وربند و زيرحوضه هاي آن را مراتع تشكيل مي دهند كه از پوشش گياهي آنچناني برخوردار نيست به همين دليل بيشتر آب جاري مي گردد و منجر به تشكيل سيلاب مي شود .
روش تحقیق
سیل یکی از مهم ترین مخاطراتی است که در مخروط افکنه ها ایجاد می شود و تحلیل آن با پیچیدگی های زیادی همراه است . از طرفی دیگر تعداد زیادی از شهرها و روستاهای ایران بر روی مخروط افکنه ها واقع شده اند و بالقوه در خطر وقوع سیلاب قرار دارند . در این تحقیق خطر سیل در شهر لار در استان فارس به روش بررسی فرضیه ها ، تحلیل محتوا ، اسنادی و تاریخی صورت گرفته است . روش تحقیق توصیفی – تحلیلی مبتنی بر مطالعات میدانی و کتابخانه ای و آنالیز نقشه ای بوده و تحقیق از نوع تحقیقات کاربردی است . اطلاعات لازم از طریق تکنیک های دور سنجی و GIS برای پهنه بندی خطر سیل فراهم گردید . روش تحقیق با توجه به اهداف و ماهیت کاری متفاوت میباشد . در اين تحقیق ابتدا به جمع آوري آمار و اطلاعات و مطالعات انجام شده در سطح منطقه ( شامل : پارامترهای اقلیمی ، ژئومورفولوژی ، زمين شناسي ، هیدروژئولوژی ، انسانی و محیطی ) . و سپس بازديد از منطقه مورد مطالعه جهت شناخت ويژگیهای طبيعی ، محیطی ،اجتماعی وا قتصادی و تطبيق نقشه هاي موجود با منطقه. مدل مورد استفاده در تحقیق الگوریتم knn می باشد . ابتدا نقشه ها با استفاده از سنجش از دور و GIS ترسیم گردیدند و سپس اطلاعات با استفاده از نرم افزار orange مورد تجزیه و تحلیل قرار گرفتند .
نقشه 1 موقعیت و مسیل های رخ داده شهر لار ( ترسیم کننده : نگارنده:1402 )
نتایج حاصل از بکار بردن GIS
الف . محدوده ارتفاعی DEM
نقشه ارتفاعی یا DEM(Digital Elevation Model) در مطالعات و پژوهش های مهم که می تواند بازتاب جهانی داشته باشد استفادهمیشود . نقشه ارتفاعی یا DEM محصول وزارت صنایع، تجارت کشور ژاپن(METI) و ایلات متحده امریکا سازمانNASA (national Aeronautics and Space Administration) شروع به تولید داده های ارتقاعی با ارتفاع 30 متر تولید کرد و ما در این پژوهش مهم هم برای محدوده مورد نظر از نقشه ارتفاعی 30 متر استفاده شده است.
نقشه 2 ارتفاع منطقه مورد مطالعه ( ترسیم کننده : نگارنده ) نقشه 3 شیب منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
ب . شیب منطقه
شیب به عنوان یکی از عوامل موثر در تعیین شناسایی نقاط خطر سیل مطرح می باشد. شیب بسیار بالای منطقه علاوه بر فرسایش ، مانع از نفوذ آب به زمین می شود و شیب کمتر فرصت نفوذ آب را فراهم آورده و به ذخیره مناسب آب کمک می کند و بر اساس بررسی ها و مطالعات شیب بین 5 الی 15 از آنجا که درصد بهترین شیب می باشد. مقدار یا ویژگییک معیار تصمیم گیری مانند مقدار شیب زمین در نقاط مختلف یکGISلایه تغییر میکند به نقاطی که دارای مقدار یا ویژگییکسان هستند ارزش (ضریب تاثیر) مساوی اختصاص مییابد.
برای تهیه نقشه شیب منطقه ، نقشه ارتفاعی منطقه تهیه گردید سپس با استفاده از نرم افزار در قسمت Spatial Analyst tool سپس گزینه Surface گزینه Slope و مقدار را بر حسب درصد بیان شده است . سپس در قسمت طبقه بندی محدوده مورد نظر را به 5 قسمت تقسیم انجام شده است . سپس مرزهای بین واحد های مجاور و با مساحت زیر 300 متر از نظر شیب حذف و نقشه شیب تهیه شده است . در نقشه شیب منطقه مطالعاتی نشان دهنده آن است که شیب کمتر از 10 درصد بیشترین عرصه را در سطح منطقه دارد و شیب های بالاتر از 10 درصد با مجموع مساحت 1131 هکتار از کل منطقه مورد مطالعه و بیشترین پتانسیل خطرات سیلابی شهری را در منطقه مورد مطالعه را در بر می گیرد . شیب زمین کمتر از 10% باشد این خصوصیت بیشتر در دره ها یا رودخانه ای باریک و یکدست می شود در شرایط استثنایی از شیب 10- 16% نیز استفاده شده است . جهت بررسی شیب منطقه ، ابتدا مدل ارتفاعی رقومی منطقه تهیه شده و سپس با استفاده از آن نقشه شیب منطقه تهیه شده ، که در این مطالعه شیب های بالاتر از 10% استفاده گردیده است.
جدول 1 مشخصات شیب منطقه مورد مطالعه نگارنده1402 )
ردیف | طبقات | کد | مساحت (هکتار) | درصد سطح نسبت به کل منطقه |
---|---|---|---|---|
1 | 0-87/3 | 1 | 964 | 46/0 |
2 | 87/3–25/12 | 2 | 401 | 19/0 |
3 | 25/12- 09/27 | 3 | 130 | 06/0 |
4 | 09/27- 83/44 | 4 | 295 | 14/0 |
5 | 83/44–25/82 | 5 | 305 | 15/0 |
پ . جهت شیب منطقه
جهت شیب نشان دهنده تاثیر متفاوت میزان دریافت نور خورشید ، بادهای گرم و خشک و میزان بارش در جهات مختلف است . رده بندی جهت شیب با توجه به حضور عوامل مختلف در جهات مختلف شیب دامنه و تفاوت در روند گسترش دامنه صورت می گیرد . جهت شیب منطقه مورد مطالعه با توجه به مساحت های بدست آمده بیشترین جهت غربی بوده و واضح است که به دلیل قرار گیری منطقه در نیم کره شمالی ، دامنه های جنوبی نسبت به دامنه های شمالی آن تابش خورشید ، تابش بیشتری از نور خورشید را در طول روز دریافت می کنند و آب و هوای خشک تری دارند و در برابر آن دامنه های مرطوب هستند .
نقشه 4 نقشه جهتشیب منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
جدول 2 مشخصات جهت شیب منطقه مورد مطالعه
ردیف | جهت شیب | مساحت | درصد سطح نسبت به کل منطقه | زوایا |
---|---|---|---|---|
1 | Flat | 5 | 0 | 1- < |
2 | North | 323 | 15/0 | 5/22-0 |
3 | NorthEast | 295 | 14/0 | 5/67-5/22 |
4 | East | 165 | 08/0 | 5/112-5/67 |
5 | SouthtEast | 108 | 05/0 | 5/157-5/112 |
6 | South | 167 | 08/0 | 5/202-5/157 |
7 | SouthWest | 153 | 07/0 | 5/247-5/202 |
8 | West | 215 | 10/0 | 5/292-5/247 |
9 | NorthWest | 363 | 17/0 | 5/337-5/292 |
10 | North | 301 | 14/0 | 360-5/337 |
ج . نقشه خاک شناسی
خاکیک پدیده طبیعی است که عموما در همه جا و در هر نوع شرایط آب و هوایی می تواند وجود داشته باشد یکی از شرایط مهم پیدایش حیات و تمدن مادی بر روی کره زمین است . از آنجا که منطقه مورد مطالعه در ناحیه خشک و کم رطوبت قرار گرفته است . کمبود رطوبت به همراه شدت بارندگی های منطقه که اغلب بصورت رگباری می باشد دست به دست هم داده و باعث می شود که خطر سیل در مناطقی رخ دهد و حتی ممکن است که از خاک از مناطق با ارتفاع بالا به مناطق مرتفع به پایین دست بیاورد و عدم پوشش گیاهی مناسب در مسیر شرایط را برایسیل و سیلاب و حتی فرسایش فراهم کند.
نقشه 5 خاکمنطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
جدول 3 مشخصات خاک منطقه مورد مطالعه
ردیف | کد | مساحت (هکتار) | مساحت کل (درصد) |
---|---|---|---|
1 | XBL | 1249 | 60/0 |
2 | XRO/ent | 21 | 01/0 |
3 | Inc | 825 | 39/0 |
چ . زمین شناسی محدوده مورد مطالعه
بررسی خصوصیات زمین شناسی و ژئومورفولوژی منطقه مورد مطالعه از جمله نیازهای اولیه وضروری برای مطالعات زمین شناسی و بعنوان پایه و اساس سایر مطالعات دیگر می باشد . منطقه مورد مطالعه معمولا از سازندهای متنوع زمین شناسی تشکیل شده که هر کدام از آنها تاثیر خاص خود را بر سایر خصوصیات دیگر مانند لغزش و رانش ، نفوذپذیری، ژئوموفولوژی منطقه و عوامل موثر بوجود آورنده اشکال مختلف زمین را دارا می باشد. همچنین با توجه به اینکه هر چه جنس لایه نفوذناپذیرتر باشد در سیل و سیلاب زمین جلوگیری می کند، چرا که از نفوذپذیری کمتری برخوردارند لذا قابلیت بالاتری جهت تشکیل مخزن نفوذناپذیر دارند بررسی زمین شناسی بیانگر این است که با توجه به نوع سنگ های این منطقه و میزان نفوذپذیری آنها کم است.در کل مطالعاتی که در گذشته صورت گرفته و همچنین مطالعه حاضر امتیاز دهی به انواع واحد های زمین شناسی و تعیین مناطق مناسب به دلیل گسسته بودن معیار زمین شناسی بر اساس نظر کارشناسی بوده است.
نقشه 6 زمین شناسی منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
جدول 4 مشخصات زمین شناسی منطقه مورد مطالعه
ردیف |
| کد | مساحت (هکتار) | مساحت کل (درصد) | امتیاز |
---|---|---|---|---|---|
1 | میشان | Mm | 12 | 0 | 1 |
2 | مارن | Mr | 229 | 11 | 2 |
3 | ابرفت | OE | 222 | 11 | 3 |
4 | ابرفت کواترنری | Qf | 1624 | 78 | 4 |
نگارنده1402 )
خ . نقشه کاربری اراضی محدوده مورد مطالعه
کاربری اراضی به عنوان یک عامل مهم در پتانسیل سنجی سیلاب در نظر گرفته است. این عامل ویژگی های زمین را تحت تاثیر قرار داده و سبب تغییر رفتار می شود . واحد اراضی به تقسیماتی در درون تیپهای اراضی اطلاق می شود که بر اساس تفاوت در پوشش فیزیکی سطح زمین در یک فیزیوگرافی می توان آنها را تشخیص داد و از هم تفکیک کرد . لذا هر یک از تیپهای فیزیوگرافی دارای شیب زیاد نظیر کوهستان و تپه ، براساس رخساره های ژئومورفولوژیکی بیرون زدگی سنگی ، دامنه پوشیده از نهشته های منفصل ، به واحدهای اراضی کوچتر تقسیم می شوند. در نقشه کاربری اراضی انواع کاربریهای موجود شامل مرتع نیمه متراکم، مرتع متراکم، مرتع کم زراعت آبی و باغات، جنگل دست کاشت، مناطق مسکونی، بیرون زدگی سنگی ، تراکم، زراعت دیم سطوح آبی مشخص گردیدند. با توجه به کاربریهای مختلف و مناسب بودن آنها نشان داده شده است . با اضافه نمودن ابزارArcSDMدر محیطGIS و با استفاده از Spatial Dara Modeller Tool قسمتFuzzy Logic دستورCategorical &Reclass نقشه شاهد وزندار کاربری اراضی تولید گردید که وزنهای کاربری آن بر اساس اهمیت آنها اولویت بندیشدهاند،با توجه به نقشه طبقه بندی کاربری اراضی و وزنهای مربوط به آنها به اراضی با بیشترین نفوذ و کمترین محدودیت بیشترین امتیاز، و به اراضی که کمترین نفوذپذیری و بیشترین محدودیت را دارا می باشند کمترین امتیاز به آنها اختصاص خواهد داشت .بررسی اثر اینعامل در پژوهش حاضر در مناطق مختلف در مطالعات انجام شده گزارش شده است که وزن های اختصاص یافته ، در تمام مطالعات انجام شد و نیز مطالعه حاضر بر اساس سلیقه و نظر کارشناسی بوده است.
نقشه 7 کاربری زمین منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
جدول 6 کاربری زمین منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
ردیف | نوع کاربری | مساحت (هکتار) | امتیاز | مساحت کل (درصد) |
---|---|---|---|---|
1 | Urban | 1080 | 1 | 51/0 |
2 | Agri | 1 | 2 | 01/0 |
3 | mix(agri_follow) | 675 | 3 | 32/0 |
4 | Poorrange | 339 | 4 | 16/0 |
ه . مسیر رودخانه منطقه مورد مطالعه
خصوصیات هندسی یا مورفومتری مسیر رودخانه ها به مجموع عوامل فیزیکی گفته می شود که مقادیر آنها برای هر برای هر حوضه نسبتا ثابت بوده و نشان دهنده وضع ظاهری است . این پارامتر از این جهت حائز اهمیت اند . به منظور تعیین محدوده مورد مطالعه و برآورد پارامترهای موفومتریک با دقت بالا از توابع تحلیل GIS استفاده شد . ابتدا بر اساس نقشه های توپوگرافی مربوط به سازمان نقشه بردارری کشور مرز و محدوده اولیه حوضه تعیین گردید و پس از مرزبندی روی نقشه های توپوگرافی با عکس های هوایی مطابقت گردیده و محدوده مورد مطالعه بطور دقیق تعیین شد. مجموع طول آبراهه در این منطقه محدوده مورد مطالعه 30995 متر است .
نقشه 8 نقشه مسیر رودخانه منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
ع . مسیر جاده محدوده مورد مطالعه
از پای شیبهای طبیعی ، وضعیت تعادلی تنش را در شیبهای مشرف بر جاده بهم می زند و در نتیجه شرایط مناسب برای بروز حرکات دامنه ای در این مناطق را فراهم می آورد. از آنجاییکه راههای در منطقه اخیرا احداث گردیده اند ، تعداد حرکات دامنه ای بوقوع پیوسته در حاشیه این مسیر ها کم می باشد و پراکندگی در این منطقه چندان از مسیرهای پیروی نمی کند لیکن در حاشیه مسیر های قدیمی قابل مشاهده اند . هر چند عدم تطابق بین سیل و سیلاب بوقوع پیوسته لیکن بروز ناپایداری های آتی در حاشیه مسیرهای قابل انتظار است، لذا در پهنه بندی خطر سیلاب شهری در این منطقه مورد مطالعه تاثیر عامل ناپایدار کننده را بصورت نیمی از طول راه بصورت از طول راه در هر واحد از شبکه های انتخاب در نظر گرفته شده است . با توجه به نحوه تاثیر عامل ناپایدار کننده را به صورت نیمی از طول راه در هر واحد از شبکه های انتخاب در نظر گرفته شده است . با توجه به نحوه تاثیر عوامل ناپایدار کننده راه و بدین منظور مجموع جاده ها و نیمی از طول راه برای واحدهای شبکه محاسبه گردید . مجموع طول جاده در این منطقه محدوده مورد مطالعه 360881 متر است . سپس فاصله اقلیدسی جاده حساب شده است .(عفیفی،1401)
نقشه 9 فاصله اقلیدسی مسیر رودخانه منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
غ . بارندگی
یکی از فاکتور های مهم در بحت پتانسیل سنجی سیلاب ، میزان بارندگی در منطقه است که با توجه به ایستگاه های هواشناسی در منطقه و یا روش های درون یابی می توان میزان بارندگی در منطقه می توان داشت که با توجه به میزان بارندگی می توان خطوط بارندگی را نمایش داد .
نقشه 10 نقشه خطوط بارندگی منطقه مورد مطالعه ( ترسیم کننده : نگارنده1402 )
الگوریتم KNN یکی از متداول ترین الگوریتم ها طبقه بندی است . این الگوریتم مبتنی بر نمونه است و بر اساس K همسایه نزدیک طبقه بندی را انجام می دهد . الگوریتم KNN به عنوان الگوریتم تنبل شناخته می شود ، زیرا مبتنی بر تقریب محلی است و همه محاسبات تا انجام طبقه بندی معوق می ماند. این روش بر اساس شباهت داده ها طبقه بندی را انجام می دهد . در واقع برای هر دادهی آزمایشی جدید ، فواصل K همسایه نزدیک را محاسبه و بر چسبی مشابه بر چسب غالب این k همسابه برای نقطه مورد را تعیین می کند. طبقه بندی کننده K نزدیکترین همسابه ، یکی از الگوریتم های طبقه بندی شناخته شده و ساده می باشد. این اولین بار توسط فیکس و هادجس به عنوان یک الگوریتم ناپارامتری معرفی شد که هیچ فرضی بر توزیع داده های ورودی ایجاد نمی کند ؛ بنابراین به طور گسترده در کاربردهای مختلف استفاده می شود.در طبقه بنـدي كننـده KNN يـك نمونـه ناشـناخته براسـا س شـباهت بـين نمونـه هـاي ، شناخته شده آموزش ديده يا برچسب دار برمبناي محاسـبه فاصـله بـين نمونـه هـاي ناشـناس بـا نمونه هاي برچسب دار، شناخته می شود. سپس K نزدیکترین نمونه هـا بـه عنـوان پايـه بـراي طبقه بندي انتخاب می شود و نمونه های نامشخص به کلاسی اختصاص می یابد که بیشترین نمونه را در میان نزدیکترین نمونه ها دارد. به همین منظور الگوریتم طبقه بندی کننده KNN بستگی دارد به 1) تعداد همسایه عدد صحیح و تغییر مقدار پارامتر K که ممکن است نتایج طبقه بندی را تغییر دهد. 2)مجموعه داده های برچسب دار بنابر این اضافه کردن یا حذف هر گونه نمونه به نمونه های آموزشی، بر تصمیم نهایی طبقه بندی کننده KNN تاثیر می گذارد و 3) معیار فاصله در KNN از فاصله اقلیدسی معمولا به عنوان معیار فاصله برای اندازه کیری فاصله بین دو نمونه استفاده می شود. طبقه بندی کننده KNN به صورت تحلیلی قابل ردیابی است و به سادکی پیاده سازی می شود اما یکی از مشکلات اصلی الگوریتم KNN این است که به همه نمونه های آموزشی نیاز دارد که در زمان اجرا در حافظه باشند، به همین دلیل طبقه بندی مبتنی بر حافظه نامیده می شود. در روشK نزدیک ترین همسایه (K-Nearest Neighbors : KNN) دسته یا کلاس یک نمونه (شئء یا مورد) جدید مشابه دسته K نزدیک ترین نمونه از مجموعه داده های آموزشی به نمونه مورد نظر است. در واقع می توان گفت پارامتر K اشاره به تعداد نمونه هایی دارد که در همسایگی نمونه مورد نظر از جهت نزدیکی بررسی می شوند. معمولا معیار اقلیدوسی برای سنجش نزدیکی مورد استفاده قرار می گیرد.(آبیل و همکاران،1398) الگوریتم روش K نزدیک ترین همسایه به صورت زیر قابل بیان است.
گام1) داده آموزشی (نمونه های آموزشی) شامل نمونه هایی که دارای بر چسب دسته یا کلاس هستند را بگیرد و K را مشخص کنید.
گام2) فاصله نمونه جدید که قرار است بر چسب دسته یا کلاس آن پیش بینی شود را با همه نمونه های آموزشی حساب کنید.
گام3) نمونه های آموزشی را بر اساس فاصله مرتب نمایید و K همسایه نزدیک را انتخاب کنید.
گام4) دسته ای که اکثریت را در K نمونهی همسایه دارد، به عنوان تخمینی برای دسته نمونه جدید انتخاب کنید.
جدول 7 ماتریس عوامل موثر در پتانسیل سنجی خطر سیلاب نگارنده1402 )
ارتفاع)متر) | شیب | جهت شیب | خاکشناسی | زمین شناسی | کاربری اراضی | مسیر رودخانه | مسیر جاده | بارندگی(میلیمتر) | |
---|---|---|---|---|---|---|---|---|---|
1 | 821 | 3 | Northeast | XBL | QF | Urban | 192 | 223 | 188 |
2 | 815 | 11 | North | XBL | QF | Urban | 155 | 372 | 187 |
3 | 816 | 32/1 | North | XBL | QF | Urban | 155 | 156 | 188 |
4 | 810 | 3 | South | XBL | QF | MIX | 174 | 80 | 196 |
5 | 851 | 9 | Southeast | Inc | QF | Urban | 160 | 266 | 243 |
6 | 820 | 2 | South | XBL | MR | Urban | 135 | 123 | 266 |
7 | 815 | 2 | East | XBL | QF | Urban | 124 | 266 | 186 |
8 | 815 | 2 | SouthWest | XBL | QF | MIX | 99 | 106 | 190 |
9 | 811 | 3 | Northeast | XBl | QF | Urban | 62 | 376 | 185 |
10 | 810 | 4 | Southeast | XBL | QF | Poorrange | 78 | 0 | 212 |
11 | 819 | 2 | Northeast | XBL | QF | Urban | 116 | 0 | 187 |
12 | 817 | 4 | Southeast | XBL | QF | Urban | 174 | 266 | 187 |
همان گونه که مشخص است، روش KNN منجر به یک مدل و سپس استفاده از آن برای پیش بینی نمی شود و در هر بار پیش بینی باید به کل داده آموزشی رجوع کرد. این امر می تواند برای داده های با حجم بالا دردسر ساز گردد. از طرفی چنانچه عدد K زوج انتخاب شود، احتمال داشتن تعداد رای مساوی برای دو برچسب دسته وجود دارد. برای جلوگیر از بوجود آمدن این مساله معمولا توصیه می شود K عددی فرد انتخاب شود. اما چنانچه زوج هم باشد، در صورتی که تعداد رای برای دو برچسب دسته مساوی باشد، بر چسب نمونه جدید را می توان با برچسب نمونه هایی که نزدیک ترین نمونه با نمونه جدید در آن حضور دارد برار گرفت. روش K نزدیک ترین همسایه یک گروه شامل K رکورد از مجموعه رکورد های آموزشی که نزدیک ترین رکوردها یا برچسب مربوط به آن ها در مورد دسته رکورد آزمایشی مزبور تصمیم گیری می نماید. به عبارت ساده تر این روش رده ای را انتخاب می کند که در همسایگی انتخاب شده بیشترین تعداد رکورد منتسب به آن دسته باشند. بنابراین رده ای که از همه ردهها بیشتر در بین K نزدیک ترین همسایه مشاهده شود، به عنوان رده رکورد جدید در نظر گرفته می شود. K نزدیکترین همسایگی یک الگوریتم سادهی طبقه بندی است. که همه موارد موجود را نگهداری و مورد جدید را بر اساس محاسبه شباهت طبقه بندی می کند.
(رابطه1)
2
Kنزدیک ترین همسایه یا همان KNN در خانواده الگوریتم هایی قرار میگیرد که ما نیاز داریم تا داده هایی را با ویژگی های خاص که با X مشخص میکنیم و با برچسب مربوطه که با Y مشخص میکنیم، را داشته باشیم تا بتوانیم در قبال نمونه جدید که برچسب آن را نمیدانیم تصمیم گیری کنیم.
KNNClassifier یک الگوریتم یادگیری مبتنی بر نمونه و ناپارامتری است که در تنظیمات دسته بندی یا همان Classification، با توجه به مقدار مشخص شده برای K، به محاسبه فاصله نقطه ای که میخواهیم برچسب آن را مشخص کنیم با نزدیک ترین نقاط میپردازد و با توجه به تعداد رای حداکثری این نقاط همسایه، در رابطه با برچسب نقطه مورد نظر تصمیم گیری میکنیم. برای محاسبه این فاصله میتوان از روش های مختلفی استفاده کرد که یکی از مطرح ترین این روش ها، فاصله اقلیدسی است که از طریق رابطه زیر میتوان آن را محاسبه کرد.
(رابطه2)
حال با توجه به اینکه مجموعه نزدیک ترین همسایه های ما، چه برچسبی دارند، با توجه به فرمول زیر به محاسبه مقدار احتمال رخداد یک برچسب پرداخته و نهایتا با توجه به اینکه احتمال وقوع کدام برچسب، بیشتر از بقیه است، به تعیین برچسب نقطه مورد نظر میپردازیم. در فرمول زیر، (I(x یک تابع شاخص است که در صورتی که x گزاره درستی باشد، مقدار آن برابر با یک و در غیر اینصورت برابر با صفر می باشد.
(رابطه3)
تعیین مقدار K در الگوریتم K نزدیک ترین همسایگی مهم ترین فاکتور تاثیرگذار بر روی نتیجه دسته بندی و برچسب گذاری است که با تغییر آن، ممکن است برچسب نمونه نیز تغییر کند. همانطور که در شکل زیر میبینید، در صورتی که مقدار K برابر با یک باشد، برچسب X به صورت "آبی" و وقتی که مقدار K برابر سه باشد، برچسب X به صورت "قرمز" پیشیبنی میشود. در پاسخ به این سوال که چطور میتوانیم بهترین مقدار را برای K در نظر بگیریم تا پیشبینی دقیق تری انجام دهیم، باید گفت تعیین مقدار K کاملا بستگی به داده هایی دارد که به بررسی آن میپردازیم و تحلیلگر باید با بهکارگیری مقادیر مختلف برای K ، بهترین مقدار از نظر دقت را مشخص نماید.
البته همانطور که در قسمت (b) شکل بالا مشهود است، بکارگیری اعداد زوج برای K مناسب به نظر نمیرسد. چون در صورتی که تعداد داده های هر کلاس با هم برابر باشند، برچسب گذاری بر روی نمونه مورد نظر سخت تر میشود. در برخی موارد نیز ممکن است برخی نقاط همسایه، فاصله بسیار کمتری، نسب به باقی همسایه ها داشته باشند که در این موارد نیز توصیه میشود به همسایه های نزدیک تر، وزن بیشتری دهیم تا از این طریق بتوانیم KNN را با دقت بالاتری پیاده سازی کنیم.
جدول 8 وزن هر کدام از عوامل موثر در پتانسیل سنجی خطر سیلاب نگارنده1402 )
ردیف | نام نقشه | وزن نقشه |
---|---|---|
1 | نقشه بارندگی | 11 |
2 | نقشه جهت شیب | 15 |
3 | نقشه رودخانه ها | 15 |
4 | نقشه کاربری اراضی | 10 |
5 | نقشه زمین شناسی | 12 |
6 | نقشه خاکشناسی | 10 |
7 | نقشه شیب | 15 |
8 | مسیر جاده | 12 |
جمع کل |
| 100 |
تلفیق نقشه ها
هدف از این نوشتار، بررسی و مقایسه و ترکیب و تلفیق نقشه های مختلف در پتانسیل سنجی سیلاب با استفاده از ترکیبات منطقی بوده و برای ترکیب نقشه به هر نقشه ای درصد کاربردی داده شده است که شامل موارد زیر می باشد . به هر یک از نقشه های اطلاعاتی بر اساس اهمیت آن در پتانسیل سنجی سیلاب در مطالعات ارزش خاصی داده می شود . برای پتانسیل سنجی خطر ، عوامل موثر اولویت بندی می شوند که در آن بارندگی و شیب مهم ترین عوامل در نظر گرفته شده است .
جدول 9 مقادیر شاخص های محاسبه شده برای عوامل موثر در پتانسیل سنجی خطر سیلاب نگارنده1402 )
شاخص های محاسبه شده | شیب | جهت شیب | بارندگی | فاصله از جاده | فاصله از رود | کاربری اراضی | خاکشناسی | زمین شناسی |
---|---|---|---|---|---|---|---|---|
| 908/0 | 873/0 | 968/0 | 951/0 | 931/0 | 911/0 | 897/0 | 942/0 |
| 564/0 | 985/0 | 836/0 | 973/0 | 896/0 | 933/0 | 894/0 | 951/0 |