آموزش استخراج حوضه آبریز و شبکه آبراهه از DEM در ArcGIS Pro و QGIS (گام‌به‌گام با داده ایران)

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

در این راهنما، کل مسیر را «از DEM خام تا لایه حوضه آماده مدل هیدرولیکی» با هم می‌رویم: انتخاب DEM مناسب برای ایران، آماده‌سازی و سیستم مختصات، اجرای ابزارها در ArcGIS Pro و QGIS، تعیین آستانه آبراهه با عدد واقعی، رتبه‌بندی استرالر، محاسبه پارامترهای فیزیوگرافی و چک‌لیست خطاهای رایج دانشجویان.

حوضه آبریز چیست و چرا استخراج دقیق آن مهم است؟

حوضه آبریز (Watershed یا Catchment) سطحی از زمین است که همه رواناب ناشی از بارش در آن، از طریق شبکه آبراهه به یک نقطه مشخص در پایین‌دست می‌رسد. این نقطه را خروجی حوضه یا Pour Point می‌نامیم؛ معمولاً محل ایستگاه هیدرومتری، محل سد، پل یا تلاقی دو رودخانه است. مرز حوضه از روی خط‌الرأس‌ها (ارتفاعات جداکننده آب) عبور می‌کند.

در GIS حوضه را با منطق «هر سلول آبش را به کدام همسایه می‌دهد؟» محاسبه می‌کنیم. وقتی برای همه سلول‌های DEM جهت جریان مشخص شد، کافی است از نقطه خروجی به سمت بالادست حرکت کنیم و همه سلول‌هایی را که در نهایت به آن می‌رسند علامت بزنیم؛ این مجموعه همان حوضه است.

تفاوت Basin و Watershed در ArcGIS Pro

این دو ابزار اغلب با هم اشتباه گرفته می‌شوند:

• Basin: بدون نیاز به نقطه خروجی، کل نقشه را به حوضه‌هایی تقسیم می‌کند که به لبه‌ی نقشه یا یک چاله ختم می‌شوند. برای دید کلی و تقسیم‌بندی منطقه مفید است.
• Watershed: حوضه‌ی بالادست یک یا چند نقطه خروجی مشخص را استخراج می‌کند. برای پروژه‌های مهندسی (سد، پل، ایستگاه) همیشه از این ابزار استفاده کنید.

انتخاب DEM مناسب برای استخراج حوضه آبریز در ایران

کیفیت مدل رقومی ارتفاع (DEM) مستقیم روی دقت مرز و شبکه آبراهه اثر می‌گذارد. برای بیشتر حوضه‌های ایران گزینه‌های رایگان زیر در دسترس است:

اگر با DEM پهپادی کار می‌کنید، حتماً کیفیت آن را با نقاط کنترل بسنجید؛ راهنمای نقاط کنترل و چک‌پوینت در نقشه‌برداری پهپاد این ارزیابی را توضیح می‌دهد. درباره روش‌های تولید DEM هم مطلب روش‌های تولید DEM در GIS را ببینید.

آماده‌سازی داده پیش از اجرای ابزارهای هیدرولوژی

1. دانلود با حاشیه: محدوده‌ای بزرگ‌تر از حوضه تخمینی (حداقل ۲ تا ۵ کیلومتر حاشیه) دانلود کنید تا خط‌الرأس واقعی بیرون از کادر نماند.
2. موزاییک و برش: اگر چند کاشی دارید، ابتدا با Mosaic To New Raster یکی کنید (نوع پیکسل 32 bit float) و سپس با Clip برش بزنید.
3. پروجکت با بازنمونه Bilinear: ابزار Project Raster را با روش Bilinear و اندازه سلول ثابت (مثلاً ۳۰ متر) اجرا کنید. روش Nearest برای DEM مناسب نیست و پله‌پله می‌شود.
4. بررسی NoData: سلول‌های بدون مقدار در دره‌ها مسیر جریان را قطع می‌کنند. با Fill حل نمی‌شوند؛ باید با Focal Statistics یا درون‌یابی پر شوند.
5. تنظیم محیط: در Environments، گزینه‌های Snap Raster و Cell Size را برابر DEM قرار دهید تا همه خروجی‌ها هم‌تراز باشند. جزئیات در هم‌ترازی رسترها در GIS.

استخراج حوضه آبریز در ArcGIS Pro (گام‌به‌گام)

همه ابزارهای این بخش در Spatial Analyst Tools ← Hydrology قرار دارند و به اکستنشن Spatial Analyst نیاز دارند.

گام ۱: پر کردن چاله‌ها با Fill

چاله (Sink) سلولی است که از همه همسایه‌هایش پایین‌تر است و آب در آن گیر می‌افتد. بیشتر چاله‌ها خطای DEM هستند، نه واقعیت. ابزار Fill آن‌ها را تا ارتفاع کمترین همسایه بالا می‌آورد.

• Input surface raster: DEM پروجکت‌شده
• Z limit: خالی (همه چاله‌ها پر شوند) مگر اینکه چاله واقعی مثل دریاچه یا کارستی‌ها داشته باشید.

برای اطمینان، قبل از Fill ابزار Sink را اجرا کنید و تعداد و عمق چاله‌ها را ببینید. اگر در DEM ۳۰ متری چاله‌ای عمیق‌تر از ۳۰ متر دیدید، احتمالاً معدن روباز، چاله برداشت شن یا خطای داده است و باید جداگانه بررسی شود.

گام ۲: جهت جریان با Flow Direction

این ابزار برای هر سلول مشخص می‌کند آب به کدام یک از هشت همسایه می‌رود. در روش D8 کدها ۱، ۲، ۴، ۸، ۱۶، ۳۲، ۶۴ و ۱۲۸ هستند (شرق، جنوب‌شرق، جنوب، جنوب‌غرب، غرب، شمال‌غرب، شمال و شمال‌شرق).

• Input: خروجی Fill
• Flow direction type: D8 (ابزار Watershed و Stream Order فقط با D8 کار می‌کنند)
• Force all edge cells to flow outward: برای حوضه‌هایی که تا لبه کادر می‌رسند فعال کنید.

گام ۳: تجمع جریان با Flow Accumulation

خروجی این ابزار برای هر سلول نشان می‌دهد چند سلول بالادست به آن آب می‌دهند. سلول‌های با مقدار زیاد بستر آبراهه‌اند و سلول‌های صفر روی خط‌الرأس قرار دارند. برای دیدن بهتر، Symbology را روی Classify با روش Geometrical Interval یا Stretch با کشش لگاریتمی بگذارید.

اگر بخواهید به‌جای «تعداد سلول» مساحت بالادست را داشته باشید، کافی است خروجی را در مساحت یک سلول ضرب کنید. برای DEM ۳۰ متری هر سلول ۹۰۰ متر مربع است؛ پس مقدار ۱۱۱۱ سلول معادل یک کیلومتر مربع است.

گام ۴: ساخت و جانمایی نقطه خروجی (Pour Point)

1. یک Feature Class نقطه‌ای با همان سیستم مختصات DEM بسازید.
2. با زوم کافی روی لایه Flow Accumulation، نقطه را دقیقاً روی سلول‌های پرمقدار آبراهه، در محل ایستگاه یا سد، قرار دهید.
3. ابزار Snap Pour Point را با Snap distance برابر ۲ تا ۵ برابر اندازه سلول (مثلاً ۶۰ تا ۱۵۰ متر برای DEM ۳۰ متری) اجرا کنید.

گام ۵: اجرای Watershed و تبدیل به پلیگون

• Input D8 flow direction raster: خروجی گام ۲
• Input raster or feature pour point data: خروجی Snap Pour Point

خروجی رستری را با Raster to Polygon (گزینه Simplify را غیرفعال کنید تا مرز با سلول‌ها منطبق بماند) به پلیگون تبدیل کنید. سپس با Calculate Geometry و Field Calculator مساحت (کیلومتر مربع) و محیط (کیلومتر) را محاسبه کنید.

استخراج آبراهه از DEM و تعیین آستانه مناسب

شبکه آبراهه را با یک شرط ساده روی Flow Accumulation می‌سازیم. در Raster Calculator:

Con("FlowAcc" >= 1111, 1)

عدد ۱۱۱۱ یعنی آبراهه از جایی شروع شود که حداقل یک کیلومتر مربع سطح بالادست داشته باشد (برای DEM ۳۰ متری). انتخاب آستانه مهم‌ترین تصمیم این مرحله است:

بهترین روش برای انتخاب علمی آستانه، مقایسه چند خروجی با آبراهه‌های نقشه توپوگرافی یا تصویر ماهواره‌ای و محاسبه تراکم زهکشی (طول کل آبراهه تقسیم بر مساحت) است. نقطه‌ای که با افزایش آستانه، تراکم زهکشی ناگهان تغییر شیب می‌دهد معمولاً انتخاب منطقی است.

رتبه‌بندی استرالر و تبدیل آبراهه به خط

1. ابزار Stream Order را با ورودی رستر آبراهه و Flow Direction و روش STRAHLER اجرا کنید.
2. با Stream to Feature خروجی را به لایه خطی تبدیل کنید؛ فیلد grid_code رتبه هر شاخه است.
3. در Symbology ضخامت خط را بر اساس رتبه تنظیم کنید تا نقشه خوانا شود؛ برای چیدمان نهایی از راهنمای طراحی لی‌اوت نقشه کمک بگیرید.

استخراج حوضه آبریز در QGIS (رایگان)

در QGIS دو مسیر اصلی دارید: الگوریتم‌های GRASS (r.watershed و r.water.outlet) و SAGA. مسیر GRASS سریع‌تر و پایدارتر است:

1. DEM پروجکت‌شده را باز کنید. از Processing Toolbox، r.fill.dir را برای پر کردن چاله‌ها اجرا کنید (یا مستقیم به مرحله بعد بروید؛ r.watershed چاله‌ها را مدیریت می‌کند).
2. r.watershed را اجرا کنید. Minimum size of exterior watershed basin را بر حسب تعداد سلول وارد کنید (مثلاً ۱۱۱۱). خروجی‌های Drainage direction، Number of cells (accumulation) و Stream segments را تیک بزنید.
3. مختصات خروجی را از روی لایه accumulation بردارید (کلیک راست روی نقشه ← Copy Coordinate).
4. r.water.outlet را با ورودی Drainage direction و مختصات خروجی اجرا کنید؛ رستر حوضه ساخته می‌شود.
5. با Raster ← Conversion ← Polygonize به وکتور تبدیل و با Field Calculator مساحت را با $area / 1000000 محاسبه کنید.

مقایسه ArcGIS Pro و QGIS در استخراج حوضه

مثال عددی: حوضه‌ای در دامنه زاگرس

فرض کنید DEM Copernicus ۳۰ متری یک حوضه کوهستانی در استان فارس (UTM زون ۳۹ شمالی) را پردازش کرده‌ایم. نتایج نمونه پس از اجرای مراحل بالا:

ضریب گراویلیوس ۱٫۶۲ نشان می‌دهد حوضه کشیده است و احتمالاً دبی اوج تیزی ندارد؛ این همان اطلاعاتی است که در فصل دوم پایان‌نامه یا گزارش مطالعات آبخیزداری نیاز دارید. برای شیب متوسط، ابتدا نقشه شیب را بسازید و سپس با Zonal Statistics میانگین را داخل مرز حوضه بگیرید.

محاسبه زمان تمرکز با فرمول کرپیچ

با داشتن طول آبراهه اصلی (L بر حسب متر) و اختلاف ارتفاع بین بالاترین نقطه آبراهه و خروجی (H بر حسب متر)، زمان تمرکز به دقیقه از رابطه Tc = 0.0195 × L^0.77 × S^-0.385 به دست می‌آید که در آن S شیب آبراهه (H/L) است. در مثال ما با L برابر ۳۱٬۴۰۰ متر و H برابر ۱۱۸۰ متر، شیب ۰٫۰۳۸ و زمان تمرکز حدود ۱۹۳ دقیقه یعنی نزدیک ۳٫۲ ساعت می‌شود. این عدد ورودی مستقیم روش استدلالی و هیدروگراف واحد است.

خطاهای رایج در استخراج حوضه آبریز و راه حل

آماده‌سازی خروجی برای HEC-RAS و مدل‌های هیدرولوژی

لایه حوضه و آبراهه‌ای که ساختید ورودی مستقیم مدل‌هایی مثل HEC-HMS (هیدرولوژی) و HEC-RAS (هیدرولیک) است. برای انتقال بدون دردسر:

• همه لایه‌ها را در یک سیستم تصویر متریک و ترجیحاً در GeoPackage یا Geodatabase تحویل دهید.
• DEM را با همان مرز حوضه و یک بافر ۵۰۰ متری Clip کنید و با فرمت GeoTIFF و نوع Float32 ذخیره کنید.
• برای HEC-HMS زیرحوضه‌ها را با چند Pour Point روی محل تلاقی‌ها بسازید و برای هر کدام مساحت، طول آبراهه، شیب و زمان تمرکز را در جدول ثبت کنید.
• برای HEC-RAS خط مرکزی رودخانه (River Centerline) را از لایه آبراهه رتبه بالا جدا کنید. معرفی نرم‌افزار را در آموزش HEC-RAS ببینید و مرحله بعد، یعنی پهنه‌بندی سیلاب با GIS، را دنبال کنید.

برای دیدن کاربردهای گسترده‌تر این لایه‌ها در هیدرولوژی، سلسله مطالب کاربرد GIS در هیدرولوژی و دوره GIS در هیدرولوژی و مدیریت منابع آب را ببینید. مستندات رسمی ابزار Flow Accumulation در ArcGIS Pro و راهنمای r.watershed در GRASS هم مرجع معتبری برای جزئیات پارامترهاست.

چک‌لیست تحویل لایه حوضه آبریز در پروژه یا پایان‌نامه

• نام و منبع DEM، تاریخ دانلود، تفکیک و سیستم مختصات در متادیتا ثبت شده است.
• تعداد و عمق چاله‌های پرشده گزارش شده است.
• آستانه آبراهه با دلیل (مقایسه با نقشه یا تراکم زهکشی) انتخاب شده است.
• محل Pour Point با مختصات و دلیل انتخاب (ایستگاه، سد، پل) مشخص است.
• مرز حوضه با تصویر ماهواره‌ای یا نقشه ۱:۵۰۰۰۰ مقایسه و اختلاف‌ها توضیح داده شده است.
• جدول پارامترهای فیزیوگرافی (مساحت، محیط، ضریب شکل، طول آبراهه، شیب، زمان تمرکز) تکمیل است.
• نقشه نهایی دارای جهت شمال، مقیاس، راهنما و منبع داده است.

پرسش‌های متداول درباره استخراج حوضه آبریز

مطالب مرتبط