Document Type : پژوهشی

Authors

Abstract

Flood is one of the most common natural worldwide hazards that causes enormous losses of life and property throughout the world. Therefore, the development of flood mapping forecast models is curial in decision making before the flood and for the after flood management. The aim of this study is to determine the flood hazard zones in the khiav Chai basin using the network analysis process. To implement this model in the area under study, various data such as rainfall, land use, slope-morphological characteristics such as convexity (profile curvature) convergence divergence slopes (plan curvature), steep slopes, vegetation index (NDVI), distance from major rivers and drainage network density were considered. The results obtained from the network Analyzer model shows the fact that more than 15 percent of the area affected by the floods with very high potential risk of occurrence are mainly located at the bottom of the basin. These levels are often less than 35% in slope, with low vegetation density profiles, converged areas with concave surfaces, and areas near rivers. Analysis of the final weights derived from the AHP in relation to flood risk shows that the slope of the geological formations operating by 0/99, and 0/822 value, have the highest impact and influence, regarding the high degree of control and influence they have on the amount of runoff and discharge of the surface area. The importance and influence of the SPI with 0/226 and 0/065 STI are of less important compared to other factors.
 

Keywords

مقدمه

سیلاب‌ها در طول تاریخ رایج‌ترین، مرگبارترین و پرهزینه‌ترین خطر در میان مخاطرات طبیعی بوده‌اند (کوسکی[1]، 2008: 61). به­عبارت دیگر سیل یکی از چند مخاطره طبیعی می‌باشد که برآورد دقیق میزان صدمات ناشی از وقوع آن امکان‌پذیر نبوده (روزالیس و همکاران[2]، 2010: 245) و خطر وقوع آن نیز در طی زمان افزایش یافته است (کوسکی، 2008: 61). دلیل این افزایش، عمدتاً گسترش شهرسازی در اطراف رودخانه‌ها و قطع درختان می‌باشد (برونستر[3]، 2003: 545؛ کریستنسن‎ها[4]، 2003: 805) به طوری که در طی سال‌های 2000 تا 2008 هر سال حدود 99 میلیون نفر در جهان تحت تأثیر نتایج حاصل از رخداد سیل قرار گرفته‌اند (شفاپور تهرانی و همکاران، 2013: 69). در ایـران نیـز همانند سایر مناطق سیل‌خیز دنیا در دهه‌های اخیر، شدت وقوع سیلاب‌ها و میـزان خسارت‌های ناشـی از آن به­طور چشمگیری افزایش یافته است (سازمان مدیریت و برنامه‌ریزی کشور، 1380: 6- 7). سطح مناطق سیل‌خیز کشور حدود 91 میلیون هکتار برآورد گردیده است. به عبارتی دیگر حدود 55 درصد از سطح کشور در تولید رواناب مستقیم و سریع نقش داشته که حدود 42 میلیون هکتار آن دارای شدت سیل‌خیزی متوسط تا خیلی زیاد هستند (شعبانلو و همکاران، 1387: 12). در آذربایجان شرقی نیز فراوانی وقوع سیل به­خصوص در ماه‌های اردیبهشت و مرداد زیاد بوده و از 191 نمونه از سیل‌های اتفاق افتاده از سال 1378 تا سال 1389، 14 درصد به شهرستان تبریز تعلق داشته و خسارت‌های زیادی نیز به بار آورده است. میزان خسارات جانی ناشی از وقوع این سیل‌ها در استان شامل: 33 نفر تلفات جانی، 23 نفر مجروح و 661 تلفات دامی نیز شامل: 661 راس دام می‌باشد. انسان با دخالت‌هایی که در چرخه طبیعت مانند تغییر کاربری و پوشش اراضی (رضایی مقدم و همکاران، 1393)، رعایت نکردن اصول شهرسازی و ساخت و ساز بدون توجه به الگوهای مکانی وقوع این پدیده می‌نماید، خسارت‌های حاصل از مخاطره طبیعی مذکور را افزایش می‌دهد. هرچند می‌توان با شناسایی مکان‌های حساس به سیل‌خیزی و مدیریت بهینه و کارآمد در این مناطق از آسیب‌پذیری بیشتر مناطق جلوگیری کرد. در زمینه شناسایی چنین مکان‌هایی مطالعات زیادی صورت گرفته است که در ادامه مطلب به چند نمونه از آن‌ها اشاره می‌گردد. هو و آمیتسو[5] (2011) با تلفیق اشکال ژئومورفولوژی و داده‌های ماهواره‌ای اعم از تصاویر لندست 7 و لایه رقومی ارتفاعی رادار شاتل[6] و با استفاده از دیاگرام‌های سه بعدی در محیط نرم‌افزار GRASS اقدام به پهنه­بندی مناطق حساس به خطر وقوع سیل در استان کنتام ویتنام کردند. زائولی[7] و همکاران (2012) برای تهیه نقشه خطر سیلاب در حوضه رودخانه Beijiang مدل SPA[8] را ارائه دادند و اعلام نمودند که اگر اعتبارسنجی بر روی مدل صورت بگیرد نتیجه مدل بهتر قابل ارائه خواهد بود. شفاپور تهرانی و همکاران (2013) با استفاده از تکنیک­های سنجش از دور و GIS اقدام به شناسایی مناطق حساس به سیل در کلانتان کشور مالزی نمودند. بدین منظور از مدل­های درخت تصمیم‌گیری و رگرسیون لجستیک همچنین معیارهایی چون بارش، کاربری-پوشش اراضی، ارتفاع، شاخص قدرت رودخانه و نوع خاک بهره بردند. کلانتری و همکاران (2014) مناطق سیل‌خیز در غرب سوئد را با استفاده از ویژگی‌های فیزیوگرافی حوضه‌ها که از تصاویر ماهواره‌ای استخراج کرده بودند، مشخص نمودند. امیدوار و همکاران (1389) با به دست آوردن 21 پارامتر ژئومتری، فیزیوگرافی، نفوذپذیری و اقلیمی در 29 زیرحوضه از حوضه آبریز کنجانچم استان ایلام به پهنه­بندی پتانسیل سیل‌خیزی اقدام کردند. در این پژوهش از روش‌های آماری تحلیل عاملی و تحلیل خوشه‌ای و معیارهای شکل حوضه، آبراهه، شیب، زهکشی مورد استفاده قرار گرفتند و بر اساس امتیاز عاملی، منطقه مورد مطالعه به 5 دسته با سیل‌خیزی زیاد تا کم طبقه‌بندی گردید. با توجه به مطالعات صورت گرفته می‌توان ذکر کرد که سنجش از دور و GIS  ابزار با ارزشی برای مدیریت مخاطرات طبیعی بوده و از این ابزارها می‌توان در تولید نقشه‌های خطر وقوع سیل (مناطقی که به­طور بالقوه خطر‌ناک می‌باشند) استفاده نمود. به عبارتی دیگر داده‌ها و تکنیکهای سنجش از دور و GIS به پایش، ارزیابی و تهیه نقشه‌های مناطق سیل‌گیر کمک کرده و در نتیجه داده‌های مطمئنی برای انجام فعالیت‌های مربوطه می‌باشند.

در مطالعه حاضر، حوضه لیقوان چای واقع در محدوده شهرستان تبریز به­عنوان منطقه مورد مطالعه انتخاب شده است. این حوضه در شهرستان تبریز واقع گردیده و رواناب­های حوضه از داخل شهر تبریز جریان می‌یابد. به عبارت دیگر در مطالعه حاضر یکی از سرمنشأ سیل‌های رخ داده در شهرستان تبریز از لحاظ حساسیت به سیل‌خیزی با استفاده از داده‌ها و تکنیک­های سنجش از دور و GIS مورد بررسی قرار گرفته است. با توجه به این­که جریان آب موجود در رودخانه این حوضه از قطب بسیار پرجمعیتی چون شهر تبریز می‌گذرد لذا شناسایی پهنه‌های سیل‌خیز این حوضه و لحاظ کردن اقدامات مدیریتی مناسب برای جلوگیری از وقوع رخداد سیل اجتناب‌ناپذیر می‌باشد.


معرفی منطقه مورد مطالعه

منطقه مورد مطالعه با مختصات جغرافیایی 37 درجه و 43 دقیقه تا 38 درجه و 2 دقیقه و 30 ثانیه عرض شـمالی و 46 درجه و 19 دقیقه و 25 ثانیه تا 46 درجه و 31 دقیقه و 16 ثانیه طول شرقی، از حوضه‌های واگرای دامنه شمالی تـوده کوهستانی سهند به مساحت 192 کیلومتر مربع است (شکل 1). لیقوان چای زهکش اصلی حوضه بوده که با جهت‌گیری جنوبی ـ شـمالی در ابتـدا بـا نـام باغچـا دره‌سی‌ چای، از نقاط و قلل مرتفع کوه سهند، مانند گیروه داغ (3596 متر)، کمال‌داغ (3236متر)، شرشر داغ و ... سرچشمه می‌گیرد. در طول مسیر آبراهه‌های متعددی چون، توله سرچای و بارالی‌چای به آن مـی‌پیونـدد. ایـن رودخانـه بـا گـذر از روستاهای سفیده‌خوان، لیقوان، هربی و دیزج عبدل، به سمت شمال تا شهر باسمنج ادامه می‌یابد و با عبور از شهر تبریـز، در حوالی فرودگاه، به رودخانة آجی‌چای می‌پیوندد (کرمی و همکاران، 1385: 128).

 

شکل (1) موقعیت منطقه مورد مطالعه در استان آذربایجان شرقی

مواد و روش­ها

در تحقیق حاضر از داده، مواد و ابزارهایی به شرح ذیل استفاده گردیده است:

از تصویر ماهواره‌ای لندست 8 (Operational Land Imager) برای استخراج نقشه کاربری و پوشش اراضی مربوط به 10 جولای سال 2013 استفاده گردید. از مدل رقومی ارتفاعی استر با اندازه پیکسل 30 متری برای استخراج حوضه و تصحیح توپوگرافیکی تصویر لندست استفاده گردید. از نقشه‌های زمین‌شناسی  1:100000 برای تولید رقومی لایه سنگ‌شناسی، نقشه‌های توپوگرافی 25000­:1 برای تصحیح هندسی تصویر ماهواره‌ای، استخراج لایه شبکه زهکشی و شیب، از دستگاه سیستم موقعیت جهانی و از سامانه گوگل ارث برای برداشت نقاط تعلیمی و از داده‌های ایستگاه‌های باران­سنجی، تبخیرسنجی و سینوپتیک اشاره شده در جدول (1)، برای تهیه لایه حداکثر بارش 24 ساعته استفاده گردید. لازم به ذکر است که لایه بارش تهیه شده با روش زمین آمار جبری و بر اساس تغییرات بارش برحسب تغییرات ارتفاع به­دست آمد. نرم‌افزارهای تخصصی ENVI4.8، IDRISI17  و ARC GIS10.1 در مطالعه حاضر مورد استفاده قرار گرفتند.

جدول (1) ایستگاه‌های مورد استفاده در تحقیق حاضر

مختصات UTM

سازمان مربوطه

ارتفاع از سطح دریا (متر)

نوع ایستگاه

نام ایستگاه

ردیف

X

Y

626681

4223677

وزارت نیرو

1500

تبخیر سنجی

سد نهند

1

612548

4215839

هواشناسی

1364

سینوپتیک

تبریز

2

611438

4189928

هواشناسی

2110

بارانسنجی

زینجناب

3

639111

4201450

هواشناسی

1950

بارانسنجی

سعید آباد

4

626141

4188470

وزارت نیرو

2200

تبخیر سنجی

لیقوان

5

630394

4197786

وزارت نیرو

1980

بارانسنجی

هروی

6

               

برای طبقه‌بندی تصویر و استخراج لایه کاربری ـ پوشش اراضی از روش نظارت شده ناپارامتریک ماشین بردار پشتیبان[9] استفاده گردید. این روش طبقه‌بندی­کننده به­صورت باینری عمل کرده و دو کلاس را با استفاده از یک فرا صفحه[10] از هم جدا می‌کند. به­طور کلی در این تحقیق از روش نظارت­شده‌ ناپارامتریک پیکسل پایه به­دلیل دقت زیاد آن در طبقه‌بندی (واپنیک[11]، 1999: 139؛ مانتروو همکاران[12]، 2005: 560؛ اندریانی، 1393: 110) استفاده گردید. برای برآورد ارتفاع رواناب از روش SCS و از رابطه (1) استفاده گردید:

رابطه (1)

 

در رابطه فوق Q: رواناب بر حسب میلی‌متر مکعب،P : میزان بارش بر حسب میلی‌متر (در تحقیق حاضر به­جای پارامتر بارش لایه استخراج شده حداکثر بارش 24 ساعته استفاده گردید)، S: عامل مربوط به نگهداشت آب در سطح زمین می‌باشد و با احتساب تلفات اولیه (2/0)، برای سیستم متریک از رابطه (2) محاسبه گردید (علیزاده، 1388: 520).

رابطه (2)

 

در رابطه فوق پارامتر [13]CN عبارت است از شماره منحنی و عددی بی‌بعد می‌باشد. دامنه عددی این پارامتر بین 100 – 0 متغیر می‌باشد، طوری که در عدد 100 تمام باران تبدیل به رواناب شده مثل مخزن آب و در عدد صفر روانابی تشکیل نمی‌شود. به­­عبارت، دیگر چنین چیزی با توجه به رابطه فوق تعریف نشده است. شماره منحنی بر اساس جدول ارائه شده از طرف سازمان حفاظت خاک امریکا که با مدنظر قرار دادن کاربری و پوشش زمین و گروه هیدرولوژی خاک به­دست آمده است، تعیین گردید (سازمان حفاظت خاک امریکا، 1973: 6). به­عبارت دیگر لایه کاربری و پوشش اراضی با گروه هیدرولوژیکی خاک منطقه مورد مطالعه در محیط ARC GIS همپوشانی شدند و مقادیر جدول )2( با توجه به اطلاعات هر ردیف پایگاه داده تشکیل گردید. لازم به ذکر است در تحقیق حاضر لایه گروه‌های هیدرولوژیکی خاک تولید شده توسط «مهندسین مشاور آب اندیشان آذر (1393)» مورد استفاده قرار گرفت.

جدول (2) مقادیر شماره منحنی برای گروه هیدرولوژیکی خاک ـپوشش موجود در منطقه

انواع کاربری و پوشش اراضی

زیر کلاس­های موجود

گروه هیدرولوژیکی خاک

B

C

D

اراضی باغی

پوشش متوسط

60

73

79

اراضی بایر

بدون پوشش و زمین لخت

86

91

94

اراضی دیم

غالباً غلات،با پوشش نسبتاً فقیر

76

84

88

مرتع

مرتع درجه یک

61

74

80

مرتع درجه دو

69

79

84

مرتع درجه سه

79

86

89

مسکونی

مناطق غیر قابل نفوذ و جاده‌ها

84

90

92

بعد از تولید لایه‌های ارتفاع، شیب، سنگ‌شناسی، تراکم زهکشی، ضریب گراویلیوس، ارتفاع رواناب، رتبه‌بندی استراهلر و کاربری و پوشش اراضی در هر یک از 13 زیر حوضه، اقدام به وزن­دهی سلسله مراتبی AHP برای فاکتورهای ذکر شده گردید. سپس زیر معیارهای طبقه‌بندی شده با توجه به اهمیت آن‌ها در تولید سیلاب با نظر کارشناسی فازی­سازی گردیده (طبق جدول 5) و سپس با روش ترکیب خطی وزن­دار مدل‌سازی و پهنه­بندی سیلاب در منطقه مورد مطالعه صورت گرفت. روش‌های ذکر شده از مبانی زیر تبعیت می‌کنند:

روش سلسله­مراتبی در سه گام اصلی: 1- تعریف معیار مورد بررسی برای تصمیم‌گیری 2- مقایسه دو به دوی معیارهای تعریف شده و 3- ارزش‌دهی معیارها با توجه به اهمیت وجودی آن‌ها صورت گرفت. البته لازم به ذکر است که اهمیت و ارزش معیارها بایستی به­صورت کارشناسی انتخاب شود. با وجود این­که ساعتی و ورگاس (1991) شاخصی را برای ناسازگاری تصمیم‌های اتخاذ شده، تعریف کرده‌اند، ولی با این وجود سلایق و دیدگاه‌های کارشناسان نیز متفاوت بوده و این مسئله تأثیر بسزایی در نتیجه تصمیم‌گیری دارد.

در  ماتریس A که معیارها دو به دو مقایسه می‌شوند به­صورت رابطه (3) قابل بیان است:

رابطه (3)

 

در رابطه (3) aijترجیح یک معیار به دیگری می‌باشد که می‌تواند به­صورت  و aij=1،nبرابر تعداد معیارها ‌باشد که در مطالعه حاضر n برابر 8 می‌باشد. wi و wj عناصر نرمال کردن وزن بردار ویژه W می‌باشد به­عبارتی جمع وزن‌های هر مقایسه به‌صورت رابطه 4 قابل بیان است: 

رابطه (4)

 

از مقایسه دو به دوی ماتریس A وزن بردار Wنیز طبق رابطه 5 بیان شد:

رابطه (5)

 

در رابطه (7)،  بزرگترین مقدار ویژه ماتریس A می‌باشد. چنانکه نتیجه ماتریس به‌صورت رابطه 6 باشد به عنوان ماتریس سازگار تلقی می‌شود.

رابطه (6)     

در این روش، نرخ ناسازگاری (CR) کمتر از 1/0 در نظر گرفته می‌شود و در صورت بیشتر شدن از این مقدار بایستی در مقایسه زوجی معیارها تجدیدنظر شود. برای به­دست آوردن نرخ ناسازگاری نیاز به محاسبه شاخص ناسازگاری (CI، در رابطه 7) و شاخص ناسازگاری تصادفی (RI) وجود دارد.

رابطه (7)

 

شاخص ناسازگاری تصادفی از جدول (3) بر اساس تعداد معیارها به دست آمد.

جدول (3) مقدار شاخص ناسازگاری تصادفی برای تعداد مختلف n (معیار) (ساعتی، 1977)

n

2

3

4

5

6

7

8

9

10

RI

0/0

58/0

9/0

12/1

1/1

32/1

41/1

45/1

49/1

نرخ ناسازگاری به‌صورت رابطه (8) بیان می‌شود:

رابطه (8)

 

وزن­دهی به معیارها از طریق روش تحلیل سلسله مراتبی و با استفاده از نظر کارشناسی صورت گرفت. بدین طریق که ابتدا در ماتریس 8 در 8 معیارها بر اساس میزان تأثیر آنها مرتب شده، سپس ما بین 1 تا 9 در مقایسه دو به دو معیارها، امتیاز به آنها تعلق گرفت.

روش ترکیب خطی وزن‌دار: در تحقیق حاضر برای تولید نقشه پهنه بندی مناطق مستعد سیلاب از روش ترکیب خطی وزن‌دار (Weighted Linear Combination) در محیط ARC GISاستفاده گردید. هدفاز ارزیابیچندمعیاری،انتخاببهترین جایگزین (Alternative)برمبنایرتبه‌بندیآنهاازطریقارزیابیچند معیاراصلیاست.در اینارزیابیبرایدستیابیبههدف، معیارهاتعریفوتعیینشدند. معیارهابه‌صورت عامل (Factor)دسته‌بندی گردیدند(مالچوفسکی[14]، 1999 :199). دراینروشتصمیم‌گیری،مقدارهر جایگزین بر اساس رابطه 9 قابل محاسبه است (پنتیوس و همکاران[15]، 2001 :195؛ اوروک[16]، 2003 :2451). 

رابطه (9)

 

در رابطه فوق: =SWمیزان مطلوبیت،  =Wiوزن هر معیار، =Siارزش استاندارد شده هر معیار می‌باشد. لازم به ذکر است که ارزش استاندارد شده هر معیار با روش فازی­سازی و نظر کارشناسی صورت گرفتبا توجه به این­که در این مجموعه‌ها مرز کاملاً مشخصی وجود ندارد، به عبارت دیگر تبدیل آنها بین عضویت و عدم عضویت در یک مجموعه به‌صورت تدریجی است، بنابراین مجموعه فازی به وسیله درجه عضویت فازی آن مشخص شده و دامنه‌ای بین 0 تا 1 به آنها تعلق گرفت که این دامنه نشانگر افزایش پیوسته از عدم عضویت به سمت عضویت کامل می‌باشد.

یافته­های تحقیق

عوامل مؤثر بر رخداد سیلاب

الف) ارتفاع رواناب

هیدروگراف هر رودخانه تحت تأثیر کاربری و پوشش گیاهی آن می‌باشد، به­طوری که اگر حوضه آبریز، فاقد پوشش گیاهی باشد، بارندگی بدون این­که جذب خاک گردد ابتدا چاله‌ها را پر کرده و آنگاه جریان سطحی آغاز می‌گردد. در صورتی که حوضه از پوشش گیاهی خوبی برخودار باشد مقداری از بارندگی توسط شاخه و برگ و ریشه گیاهان جذب می‌گردد. بنابراین CN حوضه که دربرگیرنده وضعیت پوشش گیاهی و نیز گروه‌های هیدرولوژیکی خاک از نقطه نظر سرعت نفوذ آب در آنها می‌باشد مدنظر قرار می‌گیرد. از طرفی میزان نگهداشت خاک نیز به تبعیت از شماره منحنی و در جهت عکس آن عدد مربوطه را به­خود می‌گیرد به عبارتی دیگر در مکان‌هایی که نفوذپذیری دارای مقدار کمی است میزان نگهداشت خاک آن مکان‌ها نیز کاهش می‌یابد که این مسئله با افزایش مقدار شماره منحنی قابل مشاهده است. شکل‌های (2 و 3) نقشه کاربری - پوشش اراضی، گروه هیدرولوژیکی خاک، شماره منحنی و عامل مربوط به نگهداشت خاک منطقه مورد مطالعه را نشان می‌دهد.

 

 

شکل (2) به ترتیب از راست به چپ نقشه کاربری - پوشش اراضی و گروه هیدرولوژیکی خاک

 

شکل (3) به­ترتیب از راست به چپ نقشه شماره منحنی و عامل مربوط به نگهداشت خاک

با توجه به این­که برای محاسبه ارتفاع رواناب نیاز به بارش حداکثر بارش 24 ساعته وجود دارد لذا در تحقیق حاضر بین شش ایستگاه مورد استفاده و ارتفاع آن‌ها رابطه خطی برقرار گردید. شکل (4) میزان ضریب تبیین و معادله مربوطه را نشان می‌دهد.

 

شکل (4) نمودار معادله خطی و ضریب تبیین بین ارتفاع و بارش ایستگاه‌های مورد استفاده

در معادله خطی به­دست آمده x ارتفاع هر پیکسل در لایه رقومی ارتفاعی و y میزان بارش به­دست آمده به ازای ارتفاع هر پیکسل می‌باشد. شکل (5) لایه بارش بر حسب میلی‌متر و ارتفاع رواناب به تبعیت از لایه بارش برحسب میلی‌متر را نشان می‌دهد.

 

شکل (5) به­ترتیب از راست به چپ میزان بارش 24 ساعته و ارتفاع رواناب حاصل از آن

ب) شیب حوضه

شیب حوضه شاید اصلی‌ترین عامل در بروز یا عدم بروز سیلاب در یک حوضه باشد و تأثیر عمده‌ای بر روی زمان تمرکز، مقدار جریان، نفوذ آب و نحوه تغییرات سطح آب زیرزمینی و میزان رطوبت خاک دارد. در حوضه‌هایی که از شیب بالایی برخوردار هستند، هیدروگراف مربوط به سیلاب دارای دامنه کم (مدت پایه هیدروگراف کم) و نقطه اوج بالایی است. یعنی حجم زیاد آب در مدت کوتاهی از یک مقطع عبور می‌کند. شکل (6) شیب منطقه، رتبه آبراهه‌ای منطقه به روش استراهلر[17] و دامنه و تعداد طبقات آن‌ها را نشان می‌دهد. لازم به ذکر است که رتبه‌بندی استراهلر به­صورت یک لایه جداگانه در کنار دیگر لایه‌ها در تصمیم‌گیری مورد استفاده قرار گرفت چرا که وجود رتبه بالاتر در هر یک از زیر حوضه‌ها می‌تواند زیر حوضه‌های تحت اشغال را تحت تأثیر خطر سیلاب قرار دهد.

 

شکل (6) به­ترتیب از راست به چپ نقشه شیب بر حسب درصد و رتبه‌بندی آبراهه‌ها به روش استراهلر در هر یک از زیرحوضه‌ها

ت ـ تراکم شبکه آبراهه‌ها

تراکم شبکه آبراهه‌ها در واحد سطح حوضه، طول کل آبراهه، شکل و ویژگی هندسی آبراهه‌ها و نحوه استقرار آنها در یک حوضه در ایجاد سیلاب نقش مهمی دارند. در حوضه‌ای که نسبت انشعاب آبراهه‌ها زیاد باشد طبیعتاً آب سطح حوضه توسط این کانال‌ها زودتر زهکشی می­شود. در منطقه مورد مطالعه تراکم آبراهه برای کل حوضه 07/1 به­دست آمد ولی برای دخیل دانستن معیار تراکم آبراهه در تک تک زیرحوضه، این عامل در 23 زیرحوضه محاسبه گشت (شکل 7). زیرحوضه‌ها از نظر تراکم زهکشی، دامنه عددی 77/0 الی 67/1 را به خود اختصاص داده‌اند که این امر می­تواند دلیلی بر اختلاف کم زیرحوضه‌ها از لحاظ شیب، لیتولوژی و مساحت باشد. برای دخالت دادن شکل حوضه نیز ضریب گراویلیوس مورد استفاده قرار گرفت. اگر فرم حوضه نزدیک به دایره باشد آبراهه‌ها کمترین مسیر را داشته و در نتیجه در چنین حوضه‌ای سیلاب به وقوع خواهد پیوست. در ضریب گراویلیوس پایین، تخلیه‌ی رواناب به سرعت انجام می‌گیرد و زمان تمرکز کاهش می‌یابد که این خود می‌تواند باعث افزایش خطر وقوع سیل در این مناطق‌ شود. شکل 7 تراکم آبراهه و ضریب گراویلیوس را برای زیرحوضه‌ها نشان می‌دهد. دامنه عددی این ضریب مابین 53/1 الی 95/2 می­باشد که این مقادیر نشان­دهنده دایره کامل و یا کشیدگی کامل زیرحوضه‌ها نمی‌باشد. به­عبارتی دیگر این ضریب نشان دهنده سیل­خیزی منطقه در حد متوسط می­باشد. لازم به ذکر است زیرحوضه‌هایی با طبقه یکسان در هم ادغام گردیده است.

 

شکل (7) به ترتیب از راست به چپ میزان تراکم زهکشی و ضریب گراویلیوس در هر یک از زیرحوضه‌ها

پ ـ سنگ‌شناسی منطقه و فاصله از آبراهه

نوع سنگ و پوشش خاک، ظرفیت نفوذ را تحت تأثیر قرار می‌دهد. خاک یا سنگ قابل نفوذ، شرایط نفوذ آب به داخل زمین را فراهم می‌کند و تخلیه‌ی آن را به داخل آبراهه‌ی اصلی به تأخیر می‌اندازد، از این رو رواناب سطحی کاهش می‌یابد. حوضه‌های با سنگ بستر یا خاک به نسبت غیرقابل نفوذ، حجم بالایی از رواناب سطحی ایجاد می‌کنند (گارد، به نقل از خیری­زاده و همکاران، 1391: 48).لذا در وزن­دهی و فازی­سازی میزان نفوذپذیری واحدها مدنظر قرار گرفته است. از آنجایی که منطقه مورد مطالعه در دامنه شمالی سهند واقع شده است لذا بیشتر واحدها از فعالیت آتشفشانی تأثیر پذیرفته‌اند طوری که  در منطقه مورد مطالعه چینه­شناسی پلیو – کواترنری با واحد PLQC(کنگلومرا با لایه‌هایی از لاهار) با مساحت 211 کیلومتر مربع بیشترین گسترش را در سطح منطقه دارد.به عبارتی 57 درصد از مساحت منطقه از واحد فوق‌الذکر تشکیل شده است. این واحد کنگلومرائی که خاستگاه آبرفتی دارد و در خشکی و یا در مناطق دریاچه‌ای رسوب‌گذاری کرده است دارای قطعاتی با بافت و جنس متفاوت است و گهگاه قطر قطعات آن به بیش از یک متر می‌رسد قطعات آن از سنگ‌های ولکانیکی فعالیت‌های قدیمی سهند تشکیل گردیده و سیمان آن‌ها غالباً سست بوده و از ماسه، رس و خاکسترهای آتشفشانی تشکیل شده است. با این توصیف می‌توان گفت در بیشتر مساحت منطقه نفوذپذیری سنگ‌ها کم بوده و در تولید رواناب با توجه به­وجود پوشش مرتعی ضعیف با درصد زیاد (مرتع درجه 3) بیشترین امتیاز را به­خود اختصاص می‌دهد. شکل (8) سنگ‌شناسی و واحدهای موجود در منطقه مورد مطالعه را نشان می‌دهد.

 

      شکل (8) سنگ‌شناسی منطقه مورد مطالعه                                شکل (9) لایه فاصله از آبراهه‌های منطقه

با توجه به اهمیت بالای آبراهه‌ در وقوع سیلاب، این لایه نیز به­صورت فاصله اقلیدسی در طبقاتی به مدل معرفی گردید. شکل (9) نقشه فاصله از آبراهه‌ها را نشان می‌دهد.

در مطالعه حاضر برای وزن دهی تک‌تک معیارها نسبت به هم به­صورت زوجی از روش سلسله­مراتبی با توجه به نظر کارشناسی استفاده گردید. جدول (4) نتیجه وزن­دهی به معیارهای موثر در بروز خطر سیلاب را نشان می‌دهد.

جدول (4) ماتریس امتیازدهی به معیارهای مؤثر در بروز خطر سیلاب

معیارها

1

2

3

4

5

6

7

وزن‌های مقدار ویژه

ناسازگاری

ارتفاع رواناب (1)

1

 

 

 

 

 

 

34/0

 

شیب (2)

5/0

1

 

 

 

 

 

23/0

 

فاصله از آبراهه (3)

33/0

5/0

1

 

 

 

 

18/0

 

لیتولوژی (4)

33/0

25/0

33/0

1

 

 

 

11/0

 

کاربری-پوشش اراضی(5)

14/0

2/0

25/0

25/0

1

 

 

05/0

 

تراکم آبراهه (6)

12/0

25/0

16/0

2/0

25/0

1

 

03/0

 

ضریب گراویلیوس (7)

11/0

2/0

14/0

16/0

33/0

5/0

1

02/0

 

 

 

 

 

 

 

 

 

 

08/0

جدول (4) و نحوه امتیازدهی نشان‌دهنده تأثیر بیشتر فاکتورهای شیب و سپس میزان ارتفاع رواناب بر رخداد سیلاب می‌باشد. مقدار ناسازگاری 08/ 0 نیز نشان‌دهنده امتیازدهی درست بر معیارها می‌باشد. بعد از امتیازدهی به معیارها، برای زیرمعیارها و یا طبقه‌های هر یک از معیارها (برای وارد شدن به مدل میانگین خطی وزن­دار) با استفاده از روش فازی­سازی و بر اساس میزان اهمیت آنها در رخداد سیلاب، دامنه عددی بین 0 و 1 با تابع تعریف کاربر در نظر گرفته شد. لازم به ذکر است نظر کارشناسی و پرسشنامه‌ای ملاک وزن­دهی در تحقیق حاضر بود. جدول (5) وزن هر یک از زیر معیارها را نشان می‌دهد که بر اساس این وزن‌ها طبقه‌بندی مجدد برای هر یک از معیارها صورت گرفت.


جدول (5) تعریف طبقات هر فاکتور و استانداردسازی فازی آنها

معیار

زیرمعیارها

شماره طبقه

وزن فازی­سازی (1-0)

فاصله از آبراهه (متر)

200

1

1/0

 

400 - 200

2

3/0

 

600 - 400

3

5/0

 

1000 - 600

4

8/0

 

2000 - 1000

5

1

تراکم آبراهه

< 95/0

1

1/0

 

1/1 – 95/0

2

3/0

 

25/1 – 1/1

3

5/0

 

45/1 – 25/1

4

8/0

 

45/1< 

5

1

شیب(درصد)

10-0

1

2/0

 

20-10

2

5/0

 

30-20

3

55/0

 

60 - 30

4

8/0

 

60 < 

5

1

سنگ‌شناسی

آبرفت­های جدید رودخانه‌ای

1

1/0

 

تراس­های جوان، رسوبات دشت و مخروط افکنه­ها

1

1/0

 

خاکستر آتشفشانی، لاهار

2

5/0

 

داسیت-آندزیت

4

8/0

 

داسیتیک آندزیت، کوارتز آندزیت

4

8/0

 

مارن با بین لایه­ای ماسه سنگی

3

6/0

 

کنگلومرا همراه با لایه­هایی از لاهار، توف

3

6/0

ضریب گراویلیوس

 7/1 >

1

1

 

9/1 – 7/1

2

6/0

 

1/2 – 9/1

3

4/0

 

3/2 – 1/2

4

2/0

 

3/2<

5

1/0

کاربری-پوشش اراضی

باغ و زراعت آبی

1

1/0

 

بایر

4

6/0

 

دیم

4

6/0

 

مرتع درجه 1

2

2/0

 

مرتع درجه 2

3

3/0

 

مرتع درجه 3

5

9/0

 

مسکونی

6

9/0

ارتفاع رواناب(میلی‌متر)

5 >

1

2/0

 

10 - 5

2

4/0

 

17 - 10

3

6/0

 

27 - 17

4

8/0

 

27<

5

1

با اعمال مدل تصمیم‌گیری چندمعیاره ترکیب خطی وزن­دار بر روی لایه‌های تولید شده و طبقه‌بندی مجدد شده منطقه مورد مطالعه به پنج طبقه بر اساس شکل (10) طبقه‌بندی گردید. شکل (11) درصد هر یک از پهنه‌های سیلابی در منطقه مورد مطالعه را نشان می‌دهد.

 

شکل (10) شناسایی حساسیت منطقه نسبت به رخداد سیل

 

شکل (11) درصد پهنه‌های سیلابی در منطقه مورد مطالعه

نتیجه‌گیری

شناسایی مناطق حساس به خطر سیل در حوضه‌های آبریز علی‌الخصوص در حوضه‌هایی که رودخانه دائمی داشته و از قطب‌های جمعیتی عبور می‌کنند، با اعمال روش‌های مدیریتی بهینه و کارآمد می‌تواند خسارت‌های وارده از طریق این مخاطره طبیعی را کاهش دهد. در مطالعه حاضر نیز با شناسایی عوامل موثر در رخداد سیلاب از جمله کاربری و پوشش اراضی، ارتفاع، حجم رواناب و ... با بهره‌گیری از تصاویر ماهواره‌ای، داده‌های زمینی و تکنیک­های سنجش از دور و سیستم اطلاعات جغرافیایی، مناطق حساس به سیل‌خیزی مورد پهنه بندی قرار گرفت. بدین منظور ابتدا با استفاده از تصویر ماهواره‌ای لندست 8، نقشه کاربری و پوشش اراضی مربوط به تیرماه سال 2013 استخراج گردید. که نتیجه نشان داد منطقه اکثراً تحت پوشش اراضی مرتعی به خصوص مرتع درجه 3 می‌باشد که این امر نشان‌دهنده سیل‌خیز بودن منطقه و نفوذناپذیری منطقه می‌باشد. از آن جایی که بسترهای گرانیتی نفوذناپذیر می‌باشد در نتیجه سنگ‌های این بستر جریان آب را بدون این که جذب زمین گردد به سیستم رودخانه‌ای منتقل می‌کند، در حالی که بسترهای آهکی و گچی نفوذپذیر بوده و آب را سریعاً جذب می‌کنند. منطقه مورد مطالعه به دلیل قرار گرفتن در دامنه شمالی سهند اکثر دارای لیتولوژی نفوذناپذیر می‌باشد. در ادامه برای تولید لایه مربوط به ارتفاع رواناب نیز از لایه گروه‌های هیدرولوژیکی خاک و کاربری و پوشش اراضی منطقه و ایستگاه‌های درون و نزدیک به منطقه استفاده گردید. در لایه گروه‌های هیدرولوژیکی منطقه نیز بیش از 100 کیلومترمربع جزو گروه‌های C و D می‌باشد که در این گروه­ها نفوذ آب به لایه‌های زیرین کمتر اتفاق افتاده و حجم بیشتری رواناب تولید می‌کند. لایه‌های دیگر نیز به­گونه‌ای سیل‌خیز بودن منطقه را نشان می‌دهند. به هر حال در مطالعه حاضر معیارها و زیر معیارهای آن‌ها با نظر کارشناسی به ترتیب با مدل سلسله مراتبی و فازی وزن­دهی شدند سپس با استفاده از مدل میانگین خطی وزن دار مورد پهنه بندی قرار گرفت. نتایج حاصله نشان‌دهنده سیل‌خیز بودن منطقه در پایین‌دست حوضه و ورودی شهر تبریز می‌باشد به طوری که بیش از 28 درصد از مساحت منطقه جزو مناطق با پتانسیل سیل‌خیزی از نوع زیاد و بسیار زیاد می‌باشد. با توجه به این که هر سال به دلیل وقوع سیل در مناطق ورودی شهر تبریز آسیب‌ها و خسارت‌های جانی و مالی زیادی وارد می‌شود لازم است اقدامات اساسی در بالادست حوضه به ویژه در مناطق با پهنه سیل‌خیزی بالا صورت گیرد این اقدامات را می‌توان در راستای اجرای برنامه‌های مدیریت جامع حوضه آبریز، مدیریت ریسک سیل و آبخیزداری برنامه‌ریزی نمود به گونه‌ای که در دراز مدت علاوه بر بهبود وضعیت پوشش گیاهی در سطح حوضه منجر به کاهش آسیب‌های جانی و مالی وارده گردد. در عین حال با پیگیری و اجرای مستمر قوانین موجود و تهیه و تنظیم قوانین جدید مانع تغییر کاربری اراضی در بستر و حریم رودخانه‌ها و مسیل‌های سیل‌خیز شد.



[1]-  Kusky

[2]- Rozalis et al.,

[3]- Bronstert

[4]- Christensen

[5]- Ho & Umistsu  

[6]- Shuttle Radar Topographic Mission

[7]- Zhaoli et al.,

[8]- Set pair analysis

[9]- Suport Vector Machine

[10]- Hyperplane

[11]-Vapnik

[12]-Mantero et al.,

[13]- Curve Number

[14]- Malczewski

[15]- Portius et al.,

[16]- Oruc

[17]- Sterahler

مقدمه

سیلاب‌ها در طول تاریخ رایج‌ترین، مرگبارترین و پرهزینه‌ترین خطر در میان مخاطرات طبیعی بوده‌اند (کوسکی[1]، 2008: 61). به­عبارت دیگر سیل یکی از چند مخاطره طبیعی می‌باشد که برآورد دقیق میزان صدمات ناشی از وقوع آن امکان‌پذیر نبوده (روزالیس و همکاران[2]، 2010: 245) و خطر وقوع آن نیز در طی زمان افزایش یافته است (کوسکی، 2008: 61). دلیل این افزایش، عمدتاً گسترش شهرسازی در اطراف رودخانه‌ها و قطع درختان می‌باشد (برونستر[3]، 2003: 545؛ کریستنسن‎ها[4]، 2003: 805) به طوری که در طی سال‌های 2000 تا 2008 هر سال حدود 99 میلیون نفر در جهان تحت تأثیر نتایج حاصل از رخداد سیل قرار گرفته‌اند (شفاپور تهرانی و همکاران، 2013: 69). در ایـران نیـز همانند سایر مناطق سیل‌خیز دنیا در دهه‌های اخیر، شدت وقوع سیلاب‌ها و میـزان خسارت‌های ناشـی از آن به­طور چشمگیری افزایش یافته است (سازمان مدیریت و برنامه‌ریزی کشور، 1380: 6- 7). سطح مناطق سیل‌خیز کشور حدود 91 میلیون هکتار برآورد گردیده است. به عبارتی دیگر حدود 55 درصد از سطح کشور در تولید رواناب مستقیم و سریع نقش داشته که حدود 42 میلیون هکتار آن دارای شدت سیل‌خیزی متوسط تا خیلی زیاد هستند (شعبانلو و همکاران، 1387: 12). در آذربایجان شرقی نیز فراوانی وقوع سیل به­خصوص در ماه‌های اردیبهشت و مرداد زیاد بوده و از 191 نمونه از سیل‌های اتفاق افتاده از سال 1378 تا سال 1389، 14 درصد به شهرستان تبریز تعلق داشته و خسارت‌های زیادی نیز به بار آورده است. میزان خسارات جانی ناشی از وقوع این سیل‌ها در استان شامل: 33 نفر تلفات جانی، 23 نفر مجروح و 661 تلفات دامی نیز شامل: 661 راس دام می‌باشد. انسان با دخالت‌هایی که در چرخه طبیعت مانند تغییر کاربری و پوشش اراضی (رضایی مقدم و همکاران، 1393)، رعایت نکردن اصول شهرسازی و ساخت و ساز بدون توجه به الگوهای مکانی وقوع این پدیده می‌نماید، خسارت‌های حاصل از مخاطره طبیعی مذکور را افزایش می‌دهد. هرچند می‌توان با شناسایی مکان‌های حساس به سیل‌خیزی و مدیریت بهینه و کارآمد در این مناطق از آسیب‌پذیری بیشتر مناطق جلوگیری کرد. در زمینه شناسایی چنین مکان‌هایی مطالعات زیادی صورت گرفته است که در ادامه مطلب به چند نمونه از آن‌ها اشاره می‌گردد. هو و آمیتسو[5] (2011) با تلفیق اشکال ژئومورفولوژی و داده‌های ماهواره‌ای اعم از تصاویر لندست 7 و لایه رقومی ارتفاعی رادار شاتل[6] و با استفاده از دیاگرام‌های سه بعدی در محیط نرم‌افزار GRASS اقدام به پهنه­بندی مناطق حساس به خطر وقوع سیل در استان کنتام ویتنام کردند. زائولی[7] و همکاران (2012) برای تهیه نقشه خطر سیلاب در حوضه رودخانه Beijiang مدل SPA[8] را ارائه دادند و اعلام نمودند که اگر اعتبارسنجی بر روی مدل صورت بگیرد نتیجه مدل بهتر قابل ارائه خواهد بود. شفاپور تهرانی و همکاران (2013) با استفاده از تکنیک­های سنجش از دور و GIS اقدام به شناسایی مناطق حساس به سیل در کلانتان کشور مالزی نمودند. بدین منظور از مدل­های درخت تصمیم‌گیری و رگرسیون لجستیک همچنین معیارهایی چون بارش، کاربری-پوشش اراضی، ارتفاع، شاخص قدرت رودخانه و نوع خاک بهره بردند. کلانتری و همکاران (2014) مناطق سیل‌خیز در غرب سوئد را با استفاده از ویژگی‌های فیزیوگرافی حوضه‌ها که از تصاویر ماهواره‌ای استخراج کرده بودند، مشخص نمودند. امیدوار و همکاران (1389) با به دست آوردن 21 پارامتر ژئومتری، فیزیوگرافی، نفوذپذیری و اقلیمی در 29 زیرحوضه از حوضه آبریز کنجانچم استان ایلام به پهنه­بندی پتانسیل سیل‌خیزی اقدام کردند. در این پژوهش از روش‌های آماری تحلیل عاملی و تحلیل خوشه‌ای و معیارهای شکل حوضه، آبراهه، شیب، زهکشی مورد استفاده قرار گرفتند و بر اساس امتیاز عاملی، منطقه مورد مطالعه به 5 دسته با سیل‌خیزی زیاد تا کم طبقه‌بندی گردید. با توجه به مطالعات صورت گرفته می‌توان ذکر کرد که سنجش از دور و GIS  ابزار با ارزشی برای مدیریت مخاطرات طبیعی بوده و از این ابزارها می‌توان در تولید نقشه‌های خطر وقوع سیل (مناطقی که به­طور بالقوه خطر‌ناک می‌باشند) استفاده نمود. به عبارتی دیگر داده‌ها و تکنیکهای سنجش از دور و GIS به پایش، ارزیابی و تهیه نقشه‌های مناطق سیل‌گیر کمک کرده و در نتیجه داده‌های مطمئنی برای انجام فعالیت‌های مربوطه می‌باشند.

در مطالعه حاضر، حوضه لیقوان چای واقع در محدوده شهرستان تبریز به­عنوان منطقه مورد مطالعه انتخاب شده است. این حوضه در شهرستان تبریز واقع گردیده و رواناب­های حوضه از داخل شهر تبریز جریان می‌یابد. به عبارت دیگر در مطالعه حاضر یکی از سرمنشأ سیل‌های رخ داده در شهرستان تبریز از لحاظ حساسیت به سیل‌خیزی با استفاده از داده‌ها و تکنیک­های سنجش از دور و GIS مورد بررسی قرار گرفته است. با توجه به این­که جریان آب موجود در رودخانه این حوضه از قطب بسیار پرجمعیتی چون شهر تبریز می‌گذرد لذا شناسایی پهنه‌های سیل‌خیز این حوضه و لحاظ کردن اقدامات مدیریتی مناسب برای جلوگیری از وقوع رخداد سیل اجتناب‌ناپذیر می‌باشد.


معرفی منطقه مورد مطالعه

منطقه مورد مطالعه با مختصات جغرافیایی 37 درجه و 43 دقیقه تا 38 درجه و 2 دقیقه و 30 ثانیه عرض شـمالی و 46 درجه و 19 دقیقه و 25 ثانیه تا 46 درجه و 31 دقیقه و 16 ثانیه طول شرقی، از حوضه‌های واگرای دامنه شمالی تـوده کوهستانی سهند به مساحت 192 کیلومتر مربع است (شکل 1). لیقوان چای زهکش اصلی حوضه بوده که با جهت‌گیری جنوبی ـ شـمالی در ابتـدا بـا نـام باغچـا دره‌سی‌ چای، از نقاط و قلل مرتفع کوه سهند، مانند گیروه داغ (3596 متر)، کمال‌داغ (3236متر)، شرشر داغ و ... سرچشمه می‌گیرد. در طول مسیر آبراهه‌های متعددی چون، توله سرچای و بارالی‌چای به آن مـی‌پیونـدد. ایـن رودخانـه بـا گـذر از روستاهای سفیده‌خوان، لیقوان، هربی و دیزج عبدل، به سمت شمال تا شهر باسمنج ادامه می‌یابد و با عبور از شهر تبریـز، در حوالی فرودگاه، به رودخانة آجی‌چای می‌پیوندد (کرمی و همکاران، 1385: 128).

 

شکل (1) موقعیت منطقه مورد مطالعه در استان آذربایجان شرقی

مواد و روش­ها

در تحقیق حاضر از داده، مواد و ابزارهایی به شرح ذیل استفاده گردیده است:

از تصویر ماهواره‌ای لندست 8 (Operational Land Imager) برای استخراج نقشه کاربری و پوشش اراضی مربوط به 10 جولای سال 2013 استفاده گردید. از مدل رقومی ارتفاعی استر با اندازه پیکسل 30 متری برای استخراج حوضه و تصحیح توپوگرافیکی تصویر لندست استفاده گردید. از نقشه‌های زمین‌شناسی  1:100000 برای تولید رقومی لایه سنگ‌شناسی، نقشه‌های توپوگرافی 25000­:1 برای تصحیح هندسی تصویر ماهواره‌ای، استخراج لایه شبکه زهکشی و شیب، از دستگاه سیستم موقعیت جهانی و از سامانه گوگل ارث برای برداشت نقاط تعلیمی و از داده‌های ایستگاه‌های باران­سنجی، تبخیرسنجی و سینوپتیک اشاره شده در جدول (1)، برای تهیه لایه حداکثر بارش 24 ساعته استفاده گردید. لازم به ذکر است که لایه بارش تهیه شده با روش زمین آمار جبری و بر اساس تغییرات بارش برحسب تغییرات ارتفاع به­دست آمد. نرم‌افزارهای تخصصی ENVI4.8، IDRISI17  و ARC GIS10.1 در مطالعه حاضر مورد استفاده قرار گرفتند.

جدول (1) ایستگاه‌های مورد استفاده در تحقیق حاضر

مختصات UTM

سازمان مربوطه

ارتفاع از سطح دریا (متر)

نوع ایستگاه

نام ایستگاه

ردیف

X

Y

626681

4223677

وزارت نیرو

1500

تبخیر سنجی

سد نهند

1

612548

4215839

هواشناسی

1364

سینوپتیک

تبریز

2

611438

4189928

هواشناسی

2110

بارانسنجی

زینجناب

3

639111

4201450

هواشناسی

1950

بارانسنجی

سعید آباد

4

626141

4188470

وزارت نیرو

2200

تبخیر سنجی

لیقوان

5

630394

4197786

وزارت نیرو

1980

بارانسنجی

هروی

6

               

برای طبقه‌بندی تصویر و استخراج لایه کاربری ـ پوشش اراضی از روش نظارت شده ناپارامتریک ماشین بردار پشتیبان[9] استفاده گردید. این روش طبقه‌بندی­کننده به­صورت باینری عمل کرده و دو کلاس را با استفاده از یک فرا صفحه[10] از هم جدا می‌کند. به­طور کلی در این تحقیق از روش نظارت­شده‌ ناپارامتریک پیکسل پایه به­دلیل دقت زیاد آن در طبقه‌بندی (واپنیک[11]، 1999: 139؛ مانتروو همکاران[12]، 2005: 560؛ اندریانی، 1393: 110) استفاده گردید. برای برآورد ارتفاع رواناب از روش SCS و از رابطه (1) استفاده گردید:

رابطه (1)

 

در رابطه فوق Q: رواناب بر حسب میلی‌متر مکعب،P : میزان بارش بر حسب میلی‌متر (در تحقیق حاضر به­جای پارامتر بارش لایه استخراج شده حداکثر بارش 24 ساعته استفاده گردید)، S: عامل مربوط به نگهداشت آب در سطح زمین می‌باشد و با احتساب تلفات اولیه (2/0)، برای سیستم متریک از رابطه (2) محاسبه گردید (علیزاده، 1388: 520).

رابطه (2)

 

در رابطه فوق پارامتر [13]CN عبارت است از شماره منحنی و عددی بی‌بعد می‌باشد. دامنه عددی این پارامتر بین 100 – 0 متغیر می‌باشد، طوری که در عدد 100 تمام باران تبدیل به رواناب شده مثل مخزن آب و در عدد صفر روانابی تشکیل نمی‌شود. به­­عبارت، دیگر چنین چیزی با توجه به رابطه فوق تعریف نشده است. شماره منحنی بر اساس جدول ارائه شده از طرف سازمان حفاظت خاک امریکا که با مدنظر قرار دادن کاربری و پوشش زمین و گروه هیدرولوژی خاک به­دست آمده است، تعیین گردید (سازمان حفاظت خاک امریکا، 1973: 6). به­عبارت دیگر لایه کاربری و پوشش اراضی با گروه هیدرولوژیکی خاک منطقه مورد مطالعه در محیط ARC GIS همپوشانی شدند و مقادیر جدول )2( با توجه به اطلاعات هر ردیف پایگاه داده تشکیل گردید. لازم به ذکر است در تحقیق حاضر لایه گروه‌های هیدرولوژیکی خاک تولید شده توسط «مهندسین مشاور آب اندیشان آذر (1393)» مورد استفاده قرار گرفت.

جدول (2) مقادیر شماره منحنی برای گروه هیدرولوژیکی خاک ـپوشش موجود در منطقه

انواع کاربری و پوشش اراضی

زیر کلاس­های موجود

گروه هیدرولوژیکی خاک

B

C

D

اراضی باغی

پوشش متوسط

60

73

79

اراضی بایر

بدون پوشش و زمین لخت

86

91

94

اراضی دیم

غالباً غلات،با پوشش نسبتاً فقیر

76

84

88

مرتع

مرتع درجه یک

61

74

80

مرتع درجه دو

69

79

84

مرتع درجه سه

79

86

89

مسکونی

مناطق غیر قابل نفوذ و جاده‌ها

84

90

92

بعد از تولید لایه‌های ارتفاع، شیب، سنگ‌شناسی، تراکم زهکشی، ضریب گراویلیوس، ارتفاع رواناب، رتبه‌بندی استراهلر و کاربری و پوشش اراضی در هر یک از 13 زیر حوضه، اقدام به وزن­دهی سلسله مراتبی AHP برای فاکتورهای ذکر شده گردید. سپس زیر معیارهای طبقه‌بندی شده با توجه به اهمیت آن‌ها در تولید سیلاب با نظر کارشناسی فازی­سازی گردیده (طبق جدول 5) و سپس با روش ترکیب خطی وزن­دار مدل‌سازی و پهنه­بندی سیلاب در منطقه مورد مطالعه صورت گرفت. روش‌های ذکر شده از مبانی زیر تبعیت می‌کنند:

روش سلسله­مراتبی در سه گام اصلی: 1- تعریف معیار مورد بررسی برای تصمیم‌گیری 2- مقایسه دو به دوی معیارهای تعریف شده و 3- ارزش‌دهی معیارها با توجه به اهمیت وجودی آن‌ها صورت گرفت. البته لازم به ذکر است که اهمیت و ارزش معیارها بایستی به­صورت کارشناسی انتخاب شود. با وجود این­که ساعتی و ورگاس (1991) شاخصی را برای ناسازگاری تصمیم‌های اتخاذ شده، تعریف کرده‌اند، ولی با این وجود سلایق و دیدگاه‌های کارشناسان نیز متفاوت بوده و این مسئله تأثیر بسزایی در نتیجه تصمیم‌گیری دارد.

در  ماتریس A که معیارها دو به دو مقایسه می‌شوند به­صورت رابطه (3) قابل بیان است:

رابطه (3)

 

در رابطه (3) aijترجیح یک معیار به دیگری می‌باشد که می‌تواند به­صورت  و aij=1،nبرابر تعداد معیارها ‌باشد که در مطالعه حاضر n برابر 8 می‌باشد. wi و wj عناصر نرمال کردن وزن بردار ویژه W می‌باشد به­عبارتی جمع وزن‌های هر مقایسه به‌صورت رابطه 4 قابل بیان است: 

رابطه (4)

 

از مقایسه دو به دوی ماتریس A وزن بردار Wنیز طبق رابطه 5 بیان شد:

رابطه (5)

 

در رابطه (7)،  بزرگترین مقدار ویژه ماتریس A می‌باشد. چنانکه نتیجه ماتریس به‌صورت رابطه 6 باشد به عنوان ماتریس سازگار تلقی می‌شود.

رابطه (6)     

در این روش، نرخ ناسازگاری (CR) کمتر از 1/0 در نظر گرفته می‌شود و در صورت بیشتر شدن از این مقدار بایستی در مقایسه زوجی معیارها تجدیدنظر شود. برای به­دست آوردن نرخ ناسازگاری نیاز به محاسبه شاخص ناسازگاری (CI، در رابطه 7) و شاخص ناسازگاری تصادفی (RI) وجود دارد.

رابطه (7)

 

شاخص ناسازگاری تصادفی از جدول (3) بر اساس تعداد معیارها به دست آمد.

جدول (3) مقدار شاخص ناسازگاری تصادفی برای تعداد مختلف n (معیار) (ساعتی، 1977)

n

2

3

4

5

6

7

8

9

10

RI

0/0

58/0

9/0

12/1

1/1

32/1

41/1

45/1

49/1

نرخ ناسازگاری به‌صورت رابطه (8) بیان می‌شود:

رابطه (8)

 

وزن­دهی به معیارها از طریق روش تحلیل سلسله مراتبی و با استفاده از نظر کارشناسی صورت گرفت. بدین طریق که ابتدا در ماتریس 8 در 8 معیارها بر اساس میزان تأثیر آنها مرتب شده، سپس ما بین 1 تا 9 در مقایسه دو به دو معیارها، امتیاز به آنها تعلق گرفت.

روش ترکیب خطی وزن‌دار: در تحقیق حاضر برای تولید نقشه پهنه بندی مناطق مستعد سیلاب از روش ترکیب خطی وزن‌دار (Weighted Linear Combination) در محیط ARC GISاستفاده گردید. هدفاز ارزیابیچندمعیاری،انتخاببهترین جایگزین (Alternative)برمبنایرتبه‌بندیآنهاازطریقارزیابیچند معیاراصلیاست.در اینارزیابیبرایدستیابیبههدف، معیارهاتعریفوتعیینشدند. معیارهابه‌صورت عامل (Factor)دسته‌بندی گردیدند(مالچوفسکی[14]، 1999 :199). دراینروشتصمیم‌گیری،مقدارهر جایگزین بر اساس رابطه 9 قابل محاسبه است (پنتیوس و همکاران[15]، 2001 :195؛ اوروک[16]، 2003 :2451). 

رابطه (9)

 

در رابطه فوق: =SWمیزان مطلوبیت،  =Wiوزن هر معیار، =Siارزش استاندارد شده هر معیار می‌باشد. لازم به ذکر است که ارزش استاندارد شده هر معیار با روش فازی­سازی و نظر کارشناسی صورت گرفتبا توجه به این­که در این مجموعه‌ها مرز کاملاً مشخصی وجود ندارد، به عبارت دیگر تبدیل آنها بین عضویت و عدم عضویت در یک مجموعه به‌صورت تدریجی است، بنابراین مجموعه فازی به وسیله درجه عضویت فازی آن مشخص شده و دامنه‌ای بین 0 تا 1 به آنها تعلق گرفت که این دامنه نشانگر افزایش پیوسته از عدم عضویت به سمت عضویت کامل می‌باشد.

یافته­های تحقیق

عوامل مؤثر بر رخداد سیلاب

الف) ارتفاع رواناب

هیدروگراف هر رودخانه تحت تأثیر کاربری و پوشش گیاهی آن می‌باشد، به­طوری که اگر حوضه آبریز، فاقد پوشش گیاهی باشد، بارندگی بدون این­که جذب خاک گردد ابتدا چاله‌ها را پر کرده و آنگاه جریان سطحی آغاز می‌گردد. در صورتی که حوضه از پوشش گیاهی خوبی برخودار باشد مقداری از بارندگی توسط شاخه و برگ و ریشه گیاهان جذب می‌گردد. بنابراین CN حوضه که دربرگیرنده وضعیت پوشش گیاهی و نیز گروه‌های هیدرولوژیکی خاک از نقطه نظر سرعت نفوذ آب در آنها می‌باشد مدنظر قرار می‌گیرد. از طرفی میزان نگهداشت خاک نیز به تبعیت از شماره منحنی و در جهت عکس آن عدد مربوطه را به­خود می‌گیرد به عبارتی دیگر در مکان‌هایی که نفوذپذیری دارای مقدار کمی است میزان نگهداشت خاک آن مکان‌ها نیز کاهش می‌یابد که این مسئله با افزایش مقدار شماره منحنی قابل مشاهده است. شکل‌های (2 و 3) نقشه کاربری - پوشش اراضی، گروه هیدرولوژیکی خاک، شماره منحنی و عامل مربوط به نگهداشت خاک منطقه مورد مطالعه را نشان می‌دهد.

 

 

شکل (2) به ترتیب از راست به چپ نقشه کاربری - پوشش اراضی و گروه هیدرولوژیکی خاک

 

شکل (3) به­ترتیب از راست به چپ نقشه شماره منحنی و عامل مربوط به نگهداشت خاک

با توجه به این­که برای محاسبه ارتفاع رواناب نیاز به بارش حداکثر بارش 24 ساعته وجود دارد لذا در تحقیق حاضر بین شش ایستگاه مورد استفاده و ارتفاع آن‌ها رابطه خطی برقرار گردید. شکل (4) میزان ضریب تبیین و معادله مربوطه را نشان می‌دهد.

 

شکل (4) نمودار معادله خطی و ضریب تبیین بین ارتفاع و بارش ایستگاه‌های مورد استفاده

در معادله خطی به­دست آمده x ارتفاع هر پیکسل در لایه رقومی ارتفاعی و y میزان بارش به­دست آمده به ازای ارتفاع هر پیکسل می‌باشد. شکل (5) لایه بارش بر حسب میلی‌متر و ارتفاع رواناب به تبعیت از لایه بارش برحسب میلی‌متر را نشان می‌دهد.

 

شکل (5) به­ترتیب از راست به چپ میزان بارش 24 ساعته و ارتفاع رواناب حاصل از آن

ب) شیب حوضه

شیب حوضه شاید اصلی‌ترین عامل در بروز یا عدم بروز سیلاب در یک حوضه باشد و تأثیر عمده‌ای بر روی زمان تمرکز، مقدار جریان، نفوذ آب و نحوه تغییرات سطح آب زیرزمینی و میزان رطوبت خاک دارد. در حوضه‌هایی که از شیب بالایی برخوردار هستند، هیدروگراف مربوط به سیلاب دارای دامنه کم (مدت پایه هیدروگراف کم) و نقطه اوج بالایی است. یعنی حجم زیاد آب در مدت کوتاهی از یک مقطع عبور می‌کند. شکل (6) شیب منطقه، رتبه آبراهه‌ای منطقه به روش استراهلر[17] و دامنه و تعداد طبقات آن‌ها را نشان می‌دهد. لازم به ذکر است که رتبه‌بندی استراهلر به­صورت یک لایه جداگانه در کنار دیگر لایه‌ها در تصمیم‌گیری مورد استفاده قرار گرفت چرا که وجود رتبه بالاتر در هر یک از زیر حوضه‌ها می‌تواند زیر حوضه‌های تحت اشغال را تحت تأثیر خطر سیلاب قرار دهد.

 

شکل (6) به­ترتیب از راست به چپ نقشه شیب بر حسب درصد و رتبه‌بندی آبراهه‌ها به روش استراهلر در هر یک از زیرحوضه‌ها

ت ـ تراکم شبکه آبراهه‌ها

تراکم شبکه آبراهه‌ها در واحد سطح حوضه، طول کل آبراهه، شکل و ویژگی هندسی آبراهه‌ها و نحوه استقرار آنها در یک حوضه در ایجاد سیلاب نقش مهمی دارند. در حوضه‌ای که نسبت انشعاب آبراهه‌ها زیاد باشد طبیعتاً آب سطح حوضه توسط این کانال‌ها زودتر زهکشی می­شود. در منطقه مورد مطالعه تراکم آبراهه برای کل حوضه 07/1 به­دست آمد ولی برای دخیل دانستن معیار تراکم آبراهه در تک تک زیرحوضه، این عامل در 23 زیرحوضه محاسبه گشت (شکل 7). زیرحوضه‌ها از نظر تراکم زهکشی، دامنه عددی 77/0 الی 67/1 را به خود اختصاص داده‌اند که این امر می­تواند دلیلی بر اختلاف کم زیرحوضه‌ها از لحاظ شیب، لیتولوژی و مساحت باشد. برای دخالت دادن شکل حوضه نیز ضریب گراویلیوس مورد استفاده قرار گرفت. اگر فرم حوضه نزدیک به دایره باشد آبراهه‌ها کمترین مسیر را داشته و در نتیجه در چنین حوضه‌ای سیلاب به وقوع خواهد پیوست. در ضریب گراویلیوس پایین، تخلیه‌ی رواناب به سرعت انجام می‌گیرد و زمان تمرکز کاهش می‌یابد که این خود می‌تواند باعث افزایش خطر وقوع سیل در این مناطق‌ شود. شکل 7 تراکم آبراهه و ضریب گراویلیوس را برای زیرحوضه‌ها نشان می‌دهد. دامنه عددی این ضریب مابین 53/1 الی 95/2 می­باشد که این مقادیر نشان­دهنده دایره کامل و یا کشیدگی کامل زیرحوضه‌ها نمی‌باشد. به­عبارتی دیگر این ضریب نشان دهنده سیل­خیزی منطقه در حد متوسط می­باشد. لازم به ذکر است زیرحوضه‌هایی با طبقه یکسان در هم ادغام گردیده است.

 

شکل (7) به ترتیب از راست به چپ میزان تراکم زهکشی و ضریب گراویلیوس در هر یک از زیرحوضه‌ها

پ ـ سنگ‌شناسی منطقه و فاصله از آبراهه

نوع سنگ و پوشش خاک، ظرفیت نفوذ را تحت تأثیر قرار می‌دهد. خاک یا سنگ قابل نفوذ، شرایط نفوذ آب به داخل زمین را فراهم می‌کند و تخلیه‌ی آن را به داخل آبراهه‌ی اصلی به تأخیر می‌اندازد، از این رو رواناب سطحی کاهش می‌یابد. حوضه‌های با سنگ بستر یا خاک به نسبت غیرقابل نفوذ، حجم بالایی از رواناب سطحی ایجاد می‌کنند (گارد، به نقل از خیری­زاده و همکاران، 1391: 48).لذا در وزن­دهی و فازی­سازی میزان نفوذپذیری واحدها مدنظر قرار گرفته است. از آنجایی که منطقه مورد مطالعه در دامنه شمالی سهند واقع شده است لذا بیشتر واحدها از فعالیت آتشفشانی تأثیر پذیرفته‌اند طوری که  در منطقه مورد مطالعه چینه­شناسی پلیو – کواترنری با واحد PLQC(کنگلومرا با لایه‌هایی از لاهار) با مساحت 211 کیلومتر مربع بیشترین گسترش را در سطح منطقه دارد.به عبارتی 57 درصد از مساحت منطقه از واحد فوق‌الذکر تشکیل شده است. این واحد کنگلومرائی که خاستگاه آبرفتی دارد و در خشکی و یا در مناطق دریاچه‌ای رسوب‌گذاری کرده است دارای قطعاتی با بافت و جنس متفاوت است و گهگاه قطر قطعات آن به بیش از یک متر می‌رسد قطعات آن از سنگ‌های ولکانیکی فعالیت‌های قدیمی سهند تشکیل گردیده و سیمان آن‌ها غالباً سست بوده و از ماسه، رس و خاکسترهای آتشفشانی تشکیل شده است. با این توصیف می‌توان گفت در بیشتر مساحت منطقه نفوذپذیری سنگ‌ها کم بوده و در تولید رواناب با توجه به­وجود پوشش مرتعی ضعیف با درصد زیاد (مرتع درجه 3) بیشترین امتیاز را به­خود اختصاص می‌دهد. شکل (8) سنگ‌شناسی و واحدهای موجود در منطقه مورد مطالعه را نشان می‌دهد.

 

      شکل (8) سنگ‌شناسی منطقه مورد مطالعه                                شکل (9) لایه فاصله از آبراهه‌های منطقه

با توجه به اهمیت بالای آبراهه‌ در وقوع سیلاب، این لایه نیز به­صورت فاصله اقلیدسی در طبقاتی به مدل معرفی گردید. شکل (9) نقشه فاصله از آبراهه‌ها را نشان می‌دهد.

در مطالعه حاضر برای وزن دهی تک‌تک معیارها نسبت به هم به­صورت زوجی از روش سلسله­مراتبی با توجه به نظر کارشناسی استفاده گردید. جدول (4) نتیجه وزن­دهی به معیارهای موثر در بروز خطر سیلاب را نشان می‌دهد.

جدول (4) ماتریس امتیازدهی به معیارهای مؤثر در بروز خطر سیلاب

معیارها

1

2

3

4

5

6

7

وزن‌های مقدار ویژه

ناسازگاری

ارتفاع رواناب (1)

1

 

 

 

 

 

 

34/0

 

شیب (2)

5/0

1

 

 

 

 

 

23/0

 

فاصله از آبراهه (3)

33/0

5/0

1

 

 

 

 

18/0

 

لیتولوژی (4)

33/0

25/0

33/0

1

 

 

 

11/0

 

کاربری-پوشش اراضی(5)

14/0

2/0

25/0

25/0

1

 

 

05/0

 

تراکم آبراهه (6)

12/0

25/0

16/0

2/0

25/0

1

 

03/0

 

ضریب گراویلیوس (7)

11/0

2/0

14/0

16/0

33/0

5/0

1

02/0

 

 

 

 

 

 

 

 

 

 

08/0

جدول (4) و نحوه امتیازدهی نشان‌دهنده تأثیر بیشتر فاکتورهای شیب و سپس میزان ارتفاع رواناب بر رخداد سیلاب می‌باشد. مقدار ناسازگاری 08/ 0 نیز نشان‌دهنده امتیازدهی درست بر معیارها می‌باشد. بعد از امتیازدهی به معیارها، برای زیرمعیارها و یا طبقه‌های هر یک از معیارها (برای وارد شدن به مدل میانگین خطی وزن­دار) با استفاده از روش فازی­سازی و بر اساس میزان اهمیت آنها در رخداد سیلاب، دامنه عددی بین 0 و 1 با تابع تعریف کاربر در نظر گرفته شد. لازم به ذکر است نظر کارشناسی و پرسشنامه‌ای ملاک وزن­دهی در تحقیق حاضر بود. جدول (5) وزن هر یک از زیر معیارها را نشان می‌دهد که بر اساس این وزن‌ها طبقه‌بندی مجدد برای هر یک از معیارها صورت گرفت.


جدول (5) تعریف طبقات هر فاکتور و استانداردسازی فازی آنها

معیار

زیرمعیارها

شماره طبقه

وزن فازی­سازی (1-0)

فاصله از آبراهه (متر)

200

1

1/0

 

400 - 200

2

3/0

 

600 - 400

3

5/0

 

1000 - 600

4

8/0

 

2000 - 1000

5

1

تراکم آبراهه

< 95/0

1

1/0

 

1/1 – 95/0

2

3/0

 

25/1 – 1/1

3

5/0

 

45/1 – 25/1

4

8/0

 

45/1< 

5

1

شیب(درصد)

10-0

1

2/0

 

20-10

2

5/0

 

30-20

3

55/0

 

60 - 30

4

8/0

 

60 < 

5

1

سنگ‌شناسی

آبرفت­های جدید رودخانه‌ای

1

1/0

 

تراس­های جوان، رسوبات دشت و مخروط افکنه­ها

1

1/0

 

خاکستر آتشفشانی، لاهار

2

5/0

 

داسیت-آندزیت

4

8/0

 

داسیتیک آندزیت، کوارتز آندزیت

4

8/0

 

مارن با بین لایه­ای ماسه سنگی

3

6/0

 

کنگلومرا همراه با لایه­هایی از لاهار، توف

3

6/0

ضریب گراویلیوس

 7/1 >

1

1

 

9/1 – 7/1

2

6/0

 

1/2 – 9/1

3

4/0

 

3/2 – 1/2

4

2/0

 

3/2<

5

1/0

کاربری-پوشش اراضی

باغ و زراعت آبی

1

1/0

 

بایر

4

6/0

 

دیم

4

6/0

 

مرتع درجه 1

2

2/0

 

مرتع درجه 2

3

3/0

 

مرتع درجه 3

5

9/0

 

مسکونی

6

9/0

ارتفاع رواناب(میلی‌متر)

5 >

1

2/0

 

10 - 5

2

4/0

 

17 - 10

3

6/0

 

27 - 17

4

8/0

 

27<

5

1

با اعمال مدل تصمیم‌گیری چندمعیاره ترکیب خطی وزن­دار بر روی لایه‌های تولید شده و طبقه‌بندی مجدد شده منطقه مورد مطالعه به پنج طبقه بر اساس شکل (10) طبقه‌بندی گردید. شکل (11) درصد هر یک از پهنه‌های سیلابی در منطقه مورد مطالعه را نشان می‌دهد.

 

شکل (10) شناسایی حساسیت منطقه نسبت به رخداد سیل

 

شکل (11) درصد پهنه‌های سیلابی در منطقه مورد مطالعه

نتیجه‌گیری

شناسایی مناطق حساس به خطر سیل در حوضه‌های آبریز علی‌الخصوص در حوضه‌هایی که رودخانه دائمی داشته و از قطب‌های جمعیتی عبور می‌کنند، با اعمال روش‌های مدیریتی بهینه و کارآمد می‌تواند خسارت‌های وارده از طریق این مخاطره طبیعی را کاهش دهد. در مطالعه حاضر نیز با شناسایی عوامل موثر در رخداد سیلاب از جمله کاربری و پوشش اراضی، ارتفاع، حجم رواناب و ... با بهره‌گیری از تصاویر ماهواره‌ای، داده‌های زمینی و تکنیک­های سنجش از دور و سیستم اطلاعات جغرافیایی، مناطق حساس به سیل‌خیزی مورد پهنه بندی قرار گرفت. بدین منظور ابتدا با استفاده از تصویر ماهواره‌ای لندست 8، نقشه کاربری و پوشش اراضی مربوط به تیرماه سال 2013 استخراج گردید. که نتیجه نشان داد منطقه اکثراً تحت پوشش اراضی مرتعی به خصوص مرتع درجه 3 می‌باشد که این امر نشان‌دهنده سیل‌خیز بودن منطقه و نفوذناپذیری منطقه می‌باشد. از آن جایی که بسترهای گرانیتی نفوذناپذیر می‌باشد در نتیجه سنگ‌های این بستر جریان آب را بدون این که جذب زمین گردد به سیستم رودخانه‌ای منتقل می‌کند، در حالی که بسترهای آهکی و گچی نفوذپذیر بوده و آب را سریعاً جذب می‌کنند. منطقه مورد مطالعه به دلیل قرار گرفتن در دامنه شمالی سهند اکثر دارای لیتولوژی نفوذناپذیر می‌باشد. در ادامه برای تولید لایه مربوط به ارتفاع رواناب نیز از لایه گروه‌های هیدرولوژیکی خاک و کاربری و پوشش اراضی منطقه و ایستگاه‌های درون و نزدیک به منطقه استفاده گردید. در لایه گروه‌های هیدرولوژیکی منطقه نیز بیش از 100 کیلومترمربع جزو گروه‌های C و D می‌باشد که در این گروه­ها نفوذ آب به لایه‌های زیرین کمتر اتفاق افتاده و حجم بیشتری رواناب تولید می‌کند. لایه‌های دیگر نیز به­گونه‌ای سیل‌خیز بودن منطقه را نشان می‌دهند. به هر حال در مطالعه حاضر معیارها و زیر معیارهای آن‌ها با نظر کارشناسی به ترتیب با مدل سلسله مراتبی و فازی وزن­دهی شدند سپس با استفاده از مدل میانگین خطی وزن دار مورد پهنه بندی قرار گرفت. نتایج حاصله نشان‌دهنده سیل‌خیز بودن منطقه در پایین‌دست حوضه و ورودی شهر تبریز می‌باشد به طوری که بیش از 28 درصد از مساحت منطقه جزو مناطق با پتانسیل سیل‌خیزی از نوع زیاد و بسیار زیاد می‌باشد. با توجه به این که هر سال به دلیل وقوع سیل در مناطق ورودی شهر تبریز آسیب‌ها و خسارت‌های جانی و مالی زیادی وارد می‌شود لازم است اقدامات اساسی در بالادست حوضه به ویژه در مناطق با پهنه سیل‌خیزی بالا صورت گیرد این اقدامات را می‌توان در راستای اجرای برنامه‌های مدیریت جامع حوضه آبریز، مدیریت ریسک سیل و آبخیزداری برنامه‌ریزی نمود به گونه‌ای که در دراز مدت علاوه بر بهبود وضعیت پوشش گیاهی در سطح حوضه منجر به کاهش آسیب‌های جانی و مالی وارده گردد. در عین حال با پیگیری و اجرای مستمر قوانین موجود و تهیه و تنظیم قوانین جدید مانع تغییر کاربری اراضی در بستر و حریم رودخانه‌ها و مسیل‌های سیل‌خیز شد.



[1]-  Kusky

[2]- Rozalis et al.,

[3]- Bronstert

[4]- Christensen

[5]- Ho & Umistsu  

[6]- Shuttle Radar Topographic Mission

[7]- Zhaoli et al.,

[8]- Set pair analysis

[9]- Suport Vector Machine

[10]- Hyperplane

[11]-Vapnik

[12]-Mantero et al.,

[13]- Curve Number

[14]- Malczewski

[15]- Portius et al.,

[16]- Oruc

[17]- Sterahler

منابع
ـ آقاعلیخانی، مرضیه (1388)، پهنه­بندی پتانسیل سیلخیزی و سیل­گیری حوضه فرحزاد تهران با استفاده از مدل منطق فازی، پایان­نامه کارشناسی ارشد، دانشگاه خوارزمی تهران.
ـ امیدوار، کمال؛ آمنه کیانفر (1389)، پهنه­بندی پتانسیل سیل­خیزی حوضه آبریز کنجانچم، پژوهش­های جغرافیای طبیعی، شماره 72،صص-73-90.
ـ بهشتی­جاوید، ابراهیم (1390)، پهنه­بندی پتانسیل سیل­خیزی در حوضه رودخانه بالخلوچای، پایان­نامه کارشناسی ارشد، دانشگاه خوارزمی تهران، صص150-142.
ـ بهشتی، مسعود. سادات فیض­نیا. علی سلاجقه. حسن احمدی (1388)، بررسی کارایی پهنه­بندی زمین لغزش فاکتور اطمینان: مطالعه موردی حوضه آبخیز معلم کلایه، فصلنامه جغرافیایی طبیعی، سال دوم، شماره 5، ص 32-20.
ـ پورطاهری، مهدی (1389)، کاربرد روش­های تصمیم­گیری چند شاخصه در جغرافیا، چاپ اول، انتشارات سمت، تهران. 
ـ دادرسی سبزوار، احمد (1387)، مقایسهمدلمنطقفازیباسایرمدل­هایمفهومیسازگار با GIS درمکان‌یابیمناطقمستعدگسترش سیلابباکاربرداطلاعاتماهواره‌ایسنجنده ETM. همایش ژئوماتیک 87 ، 22 تا 23 اردیبهشت ‌ماه، سازمان نقشه‌برداری کشور ، تهران.
ـ حسین­زاده، سیدرضا؛ مهناز جهادی طرقی (1385)، تجزیه و تحلیل ژئومورفولوژیک سیلاب­های کاتاستروفیک رودخانه­ها در سو(جنگل گلستان)، جغرافیا و توسعه ناحیه­ای، شماره 7، ص 89-115.
ـ خیری­زاده، منصور؛ ملکی، جبرئیل؛ عمونیا، حمید (1392)، پهنهبندی پتانسیل خطر وقوع سیلاب در حوضهی آبریز مردق چای با استفاده از مدل ANP، فصلنامه پژوهش های ژئومورفولوژی کمی ، سال اول، شماره 3،صص 56-39.
ـ سوری، مهشید؛ جعفری، محمد؛ آذرنیوند، حسین؛ فرخ­زاده، بهنوش (1391)، تعیین مناطق مناسب اجرای پروژه پخش سیلاب با استفاده از فرآیند تحلیل سلسله­مراتبی و سیستم اطلاعات جغرافیایی (مطالعه موردی: حوضه آبخیز میخوران کرمانشاه)، پژوهش های آبخیزداری (پژوهش و سازندگی)، زمستان 1391، دوره 25، شماره 4 صفحه 92 تا صفحه 103.
ـ صیاد اصغری سرسکانرود، صالح اصغری سرکانرود، بتول زینالی (1392)، برآورد سیلاب حوضه آبخیز و مخاطرات آن در حوضه آلانق، دومین کنفرانس بین­المللی مخاطرات محیطی، تهران، ایران.
ـ کرم امیر، درخشان فرزانه (1391)، پهنه­بندی سیل خیزی، برآورد سیلاب و ارزیابی کارایی کانال­های دفع آب­های سطحی در حوضه­های شهری (مطالعه موردی: حوضه آبشوران در کرمانشاه)، جغرافیای طبیعی:  تابستان 1391 , دوره  5 , شماره  16 , از صفحه 37 تا صفحه 54.
ـ کرمی، حجت؛ اردشیز، عبدالله؛ حسینی، سیدهادی؛ میکائیلی، محمدعلی (1387)، پهنه­بندی سیلابدشت با تلفیق مدل هیدرولیکی و سامانه اطلاعات جغرافیایی، مقاله سمینار چهاردهمین کنفرانس دانشجویان مهندسی عمران سراسر کشور.
ـ قاسمی، علی، سلاجقه، علی، ملکیان،آرش، اسمعلی، اباذر (1393)، بررسی سیل خیزی و تعیین عوامل مؤثر در آن در حوضه رودخانه بالقلی چای با استفاده از تکنیک GIS، RSو AHP، فصلنامه محیط­شناسی، سال چهلم، شماره 70، صص 4000-389.
ـ مؤمنی، محمد و شریفی سلیم، علی (1390)، مدل­ها و نرم­افزارهای تصمیم­گیری چندشاخصه، چاپ اول، نشر مؤلفین، تهران.
-Agha Alikhani, M. (2009), Flood zonation of Farahzad Tehran basin using fuzzy logic models, Thesis, University of Tehran Kharaizmi. In Persian
-Akgun, A. Kincal, C. (2010), Landslide risk as an environmental threat for Izmir (West Turkey) City and its urban expansion. EGU General Assembly Conference Abstracts p. 1249.
-Akgun, A. Türk, N. (2010), Landslide susceptibility mapping for Ayvalik (Western Turkey) and its vicinity by multicriteria decision analysis. Environ. Earth Sci. 61, 595–611.
-Akgun, A. Sezer, E.A. Nefeslioglu, H.A. Gokceoglu, C. Pradhan, B. (2012), An easy-to-use MATLAB program (MamLand) for the assessment of landslide susceptibility using a Mamdani fuzzy algorithm. Comput. Geosci. 38 (1), 23–34.
-Aleotti, P. Chowdhury, R. (1999), Landslide hazard assessment: summary review and new perspectives. Bull. Eng. Geol. Environ. 58, 21–44.
-Althuwaynee, O.F. Pradhan, B. Lee, S. (2012a), Application of an evidential belief function model in landslide susceptibility mapping. Comput. Geosci. 44, 120–135.
-Althuwaynee, O.F. Pradhan, B. Mahmud, A.R. Yusoff, Z.M. (2012b), Prediction of slope failures using bivariate statistical based index of entropy model, Humanities, Science and Engineering (CHUSER). IEEE Colloq. IEEE 362–367.
-Bates, P.D. (2004), Remote sensing and flood inundation modelling. Hydrol. Process. 18, 2593–2597.
-Beheshti.J,E. E. (2011), Zonation Balkhli chay flooding in the river basin, MSc Thesis, University of Tehran Kharaizmi. In Persian.
-Beheshti, M.Feiz nia, S. A. Salajegheh. Ahmad, H. (2009), Studied the efficacy of zoning landslide Confidence Factor: A Case Study of Moalem Kalaye watershed, physiographic Quarterly, Vol 2, No. 5, pp 20-32. In Persian.
-Dadrasi Sabzevar, A. (2008), Comparison of fuzzy logic model is consistent with conceptual models with GIS for locating flood-prone areas by using satellite data, Landsat ETM. Geomatics 87, 22 and 23 April, the topography of the country, Tehran. In Persian.
-Devkota, K. Regmi, A. Pourghasemi, H. Yoshida, K. Pradhan, B. Ryu, I. Dhital, M. Althuwaynee, O. (2013), Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling–Narayanghat road section in Nepal Himalaya. Nat. Hazards 65, 135–165.
-Elshorbagy, A. Corzo, G. Srinivasulu, S. Solomatine, D. (2010a), Experimental investigation of the predictive capabilities of data driven modeling techniques in hydrology-Part 1: concepts and methodology. Hydrol. Earth Syst. Sci. 14, 1931–1941.
-Feng LH, Lu J (2010), The practical research on flood forecasting based on artificial neural networks. Expert Syst Appl37:2974–2977
-Fenicia, F. Kavetski, D. Savenije, H.H. Clark, M.P. Schoups, G. Pfister, L. Freer, J. (2013), Catchment properties, function, and conceptual model representation: is, there a correspondence Hydrol. Process... http://dx.doi.org/10.1002/hyp.9726.79-85.
-Hagen E.J.F. Shrodr Jr. X.X.Lu, John F.Teufert. (2010), Reverse engineered flood hazard mapping in Afghanistan: A parsimonious flood map model for developing Countries. Quaternary International, Pages 1-10.
-Khu, S.T. Liong, S.Y. Babovic, V. Madsen, H. Muttil, N. (2001), Genetic programming and its application in real-time runoff forecasting1, Water Resour. Assoc. 37, pp 439–451.
-Lee, Seungbum, Walsh, Patrick, (2011), SWOT and AHP hybrid model for sport marketing outsourcing, Sport Management Review, Vol.14, No. 4, 361-369.
-Mugagga, F. Kakembo, V. Buyinza, M. (2012), Land use changes on the slopes of Mount Elgon and the implications for the occurrence of landslides. Catena 90, pp 39–46.
-Movoh, J.A. Zhu, J.z. (1998), Application of AHP/ANP to Unit Commitment in the Deregulated Power Industry, In: IEEE International Conference on System, Man and Cybernetics, Vol. 1 San Diego, pp. 817-822.
-Moore, I.D. Burch, G.J. (1986), Sediment transport capacity of sheet and rill flow: application of unit stream power theory.Water Resource 22, 1350-1360.
-Nefeslioglu,H.A. Duman, T.Y. Duemaz,S. (2008), Landslide susceptibility mapping for a part of tectonic Kelkit Valley (Eastern Black Sea region of Turkey). Geomorphology 94,pp 401–418.
-Oh, H.J. Pradhan, B. (2011), Application of a neuro-fuzzy model to landslidesusceptibility mapping for shallow landslides in a tropical hilly area. Comput. Geosci. 37, pp 1264–1276.
-Pradhan, B. Shafiee, M. (2009), Flood hazrad assessment for cloud prone rainy areas in a typical tropical environment. Disaster Adv. 2, pp 7–15.
-Pradhan, B. Mansor, S. Pirasteh, S. Buchroithner, M.F. (2011), Landslide hazard and risk analyses at a landslide prone catchment area using statistical based geospatial model. Int. J. Remote Sens. 32 (14), pp 4075–4087.
-Pradhan, B. (2011), Use of GIS-based fuzzy logic relations and its cross application to produce landslide susceptibility maps in three test areas in Malaysia. Environ. Earth Sci. 63, pp 329–349.
-Pradhan, B. Hagemann, U. Shafapour Tehrany, M. Prechtel, N. (2014), An easy to use ArcMap based texture analysis program for extraction of flooded areas from TerraSAR-X satellite image. Comput. Geosci. 63, pp 34–43.
-Shafapour Tehrany, M., Moung-Jin Lee ., Pradhan, B., Neamah Jebur,M., Lee,S. (2014), Flood susceptibility mapping using integrated bivariate and multivariate statistical models, Environ Earth Sci., DOI 10.1007/s12665-014-3289-3.
-Sahu. R.K, S.K.Mishra, T.I.Eldho. (2010), ComparativeevaluationofSCS-CN-inspiredmodelsinapplicationstoclassifieddatasets. Agricultural Water Management 97(2010) pp 749–756.
-Schmidt, F. Persson, A. (2003), Comparison of DEM data capture and topographic wetness indices. Precision Agriculture 4, pp 179–192.
-Smith, K. Ward, R. (1998), Floods: Physical Processes and Human Impacts. Wiley, Chichester, pp. 382.
-Sorensen, R. Zinko, U. Seibert, J. (2006), On the calculation of the topographic wetness index: evaluation of different methods based on field observations. Hydrology and Earth System Sciences Discussions 2, pp 1807–1834.
-Wanders, N. Bierkens, M.F. de Jong, S.M. de Roo, A. Karssenberg, D. (2013), The benefits of using remotely sensed soil moisture in parameter identification of large-scale hydrological models. In: EGU General Assembly Conference Abstracts, vol. 15, pp. 10215.
-Zou, Q. Zhou, J. Zhou, C. Song, L. Guo, J. (2013), Comprehensive flood risk assessment based on set pair analysis-variable fuzzy sets model and fuzzy AHP. Stoch. Environ. Res. Risk Assess. 27, 525–546.