چرا تدقیق دادهٔ اولیه مهم‌ترین گام است؟
مدل هیدرولوژیکی، آستانهٔ هشدار و سنجش دقت پیش‌بینی همه بر دادهٔ ایستگاه‌های زمینی تکیه دارند. یک خطای ساعت ثبت، یک جابه‌جایی صفر اشل، یا منحنی دبی–اشل قدیمی می‌تواند اوج سیلاب را چند ده درصد جابه‌جا کند. در حوضهٔ گرگان‌رود این خطر بیشتر است: سیلاب‌های تند با رسوب لسی بالا بستر را تغییر می‌دهند، ایستگاه‌ها در سیلاب آسیب می‌بینند، و سدهای گلستان و وشمگیر جریان طبیعی را تغییر می‌دهند.
۴۵ ایستگاه هیدرومتری
داخل حوضه (فایل ملی waterse.ir)
۷۶ باران‌سنج
داخل حوضه؛ بیشتر غیرثبات و روزانه
۵۰ ایستگاه اقلیم‌شناسی
داخل حوضه
۱ ایستگاه سینوپتیک
داخل حوضه؛ ساعتی
زنجیرهٔ کیفیت داده (Data Levels) — مطابق رویهٔ WMO
Level 0 · Rawدادهٔ خام همان‌طور که از نهاد رسیده؛ هرگز تغییر نمی‌کند.
Level 1 · Automated checksآزمون‌های خودکار و پرچم‌گذاری؛ هیچ مقداری حذف نمی‌شود.
Level 2 · Expert reviewبررسی کارشناس، تصمیم دربارهٔ مقادیر مشکوک، اصلاح با ثبت دلیل.
Level 3 · Corrected / filledاصلاح سیستماتیک، سرشکنی خطا و پر کردن خلأ با برآورد عدم قطعیت.
Level 4 · Model-readyسری‌های همگن و تأییدشده برای واسنجی مدل و آستانه.
پرچم‌های کیفیت (Quality flags) — به هر مقدار یک پرچم
G · Good
S · Suspect
B · Bad / rejected
C · Corrected
E · Estimated (gap-filled)
M · Missing
انواع خطا در داده‌های زمینی
| نوع خطا | بارندگی (Rainfall) | هیدرومتری (Stage / Discharge) |
|---|
Gross errors
خطای فاحش | رقم اشتباه، جابه‌جایی ممیز، ثبت در روز اشتباه، واحد اشتباه | خواندن غلط اشل، تبدیل اشتباه واحد، جابه‌جایی ستون‌ها |
Systematic errors
خطای سیستماتیک | کم‌برآورد باد (wind undercatch)، تبخیر از ظرف، تغییر محل یا نوع دستگاه | جابه‌جایی صفر اشل (datum shift)، منحنی دبی–اشل قدیمی، رانش حسگر (sensor drift) |
Random errors
خطای تصادفی | خطای خواندن، پاشش، ریزش ناقص | نوسان سطح آب، خطای اندازه‌گیری سرعت در دبی‌سنجی |
Timing errors
خطای زمان | ساعت محلی/UTC، ثبت بارش ۲۴ ساعته در روز بعد | ساعت دیتالاگر، تأخیر انتقال تله‌متری |
Event-specific
ویژهٔ سیلاب | سرریز ظرف در رگبار شدید، آسیب ایستگاه | تخریب ایستگاه، رسوب و گرفتگی، برون‌یابی منحنی در دبی‌های بالا، اثر پس‌زدگی (backwater) |
آزمون‌های رهیابی خطا (Error detection)
بارندگی (Rainfall)
- Range / limit check: مقدار منفی یا بیشتر از حد اقلیمی (مثلاً بیشینهٔ ثبت‌شدهٔ منطقه × ضریب ایمنی).
- Internal consistency: جمع ساعتی = روزانه؛ روزانه = ماهانه.
- Temporal consistency: مقدار ثابت تکراری (flatline)، جهش نامعقول، روزهای صفر طولانی در فصل بارش.
- Spatial consistency: مقایسه با ایستگاه‌های همسایه (spatial regression test، برآورد IDW از همسایه‌ها و آزمون انحراف).
- Double-mass curve: شکست شیب منحنی تجمعی ایستگاه در برابر میانگین همسایه‌ها = تغییر سیستماتیک.
- Homogeneity tests: SNHT (Alexandersson)، Pettitt، Buishand range، von Neumann ratio.
- Outlier tests: Grubbs، Tukey IQR، روی مقادیر غیرصفر.
- Independent cross-check: مقایسه با GPM IMERG، رادار هواشناسی و ERA5 (فقط برای آشکارسازی، نه جایگزینی).
هیدرومتری (Hydrometric)
- Rating-curve check: باقیماندهٔ دبی‌سنجی‌ها نسبت به منحنی؛ آشکارسازی rating shift پس از هر سیل.
- Rate-of-change limit: سرعت بالا/پایین رفتن اشل بیش از حد فیزیکی.
- Continuity / mass balance: دبی پایین‌دست ≈ بالادست + سرشاخه‌ها + ورودی جانبی، با تأخیر روندیابی (routing lag).
- Rainfall–runoff plausibility: ضریب رواناب رویداد بین ۰ و ۱؛ حجم سیلاب بیشتر از حجم باران ممکن نیست.
- Reservoir water balance: تغییر حجم سدهای گلستان، وشمگیر و بوستان = ورودی − خروجی − تبخیر.
- Flatline / spike detection: گیر کردن حسگر یا قطع ارتباط تله‌متری.
- Datum check: مقایسهٔ صفر اشل با نقشه‌برداری و ارتفاع DEM.
- Hydrograph shape: زمان رسیدن اوج در ایستگاه‌های پشت سر هم باید با فاصله و شیب سازگار باشد.
سرشکن کردن و اصلاح خطا (Error distribution & correction)
سرشکنی کمترین مربعات (Least-squares adjustment / Data reconciliation)
دبی ایستگاه‌های یک شبکهٔ رودخانه‌ای باید در پیوستگی جرم صدق کنند. باقیماندهٔ پیوستگی به نسبت عدم قطعیت هر ایستگاه بین آن‌ها تقسیم می‌شود: ایستگاه کم‌دقت‌تر سهم بیشتری از اصلاح می‌گیرد.
min Σ (Q̂ᵢ − Qᵢ)² / σᵢ² subject to Q̂_down − Σ Q̂_up − Q_lateral = 0
اصلاح سیستماتیک با Double-mass
اگر شیب منحنی تجمعی ایستگاه پس از تاریخ t تغییر کرده، دادهٔ دورهٔ نامعتبر با نسبت شیب‌ها به دورهٔ مرجع برگردانده می‌شود.
P_corrected = P_observed × (b_reference / b_observed)
بازواسنجی منحنی دبی–اشل (Rating-curve recalibration)
برازش دوباره با دبی‌سنجی‌های جدید و تفکیک دوره‌ها پس از هر جابه‌جایی بستر؛ روش‌های بیزی مانند BaRatin باند عدم قطعیت دبی را هم می‌دهند، به‌ویژه برای برون‌یابی در دبی‌های سیلابی.
Q = a · (h − h₀)ᵇ (fitted per stable period, with uncertainty band)
اصلاح کم‌برآورد باران‌سنج (Undercatch correction)
ضریب اصلاح بر پایهٔ سرعت باد و نوع باران‌سنج (رویکرد WMO intercomparison)؛ برای ایستگاه‌های بادخیز دشت و ساحل.
پر کردن خلأ داده (Gap filling)
| روش | کاربرد | یادداشت |
|---|
| Normal ratio method | بارش روزانه از ایستگاه‌های همسایه | Pₓ = (1/n) Σ (Nₓ/Nᵢ) Pᵢ — N: میانگین سالانه |
| Inverse Distance Weighting (IDW) | بارش با همسایه‌های نزدیک | ساده؛ اثر ارتفاع را نمی‌بیند |
| Regression / MOVE | تکمیل و تطویل سری دبی از ایستگاه همبسته | واریانس سری حفظ می‌شود (MOVE.1) |
| Kriging with external drift | بارش با ارتفاع (DEM) به‌عنوان متغیر کمکی | مناسب کوهپایه‌های البرز |
| Machine learning (Random Forest, LSTM) | خلأهای بلند؛ ترکیب ایستگاه، ماهواره و ERA5 | فقط با اعتبارسنجی متقاطع (cross-validation) |
| Gauge–satellite merging | نقشهٔ بارش ساعتی پیوسته | conditional merging، quantile mapping برای تصحیح IMERG |
قاعده: مقدار پرشده هرگز جای دادهٔ اصلی نمی‌نشیند؛ با پرچم E و برآورد عدم قطعیت ذخیره می‌شود.
بهینه‌سازی (Optimization)
بهینه‌سازی شبکهٔ ایستگاه‌ها
- مقایسه با حداقل تراکم توصیه‌شدهٔ WMO برای مناطق کوهستانی و دشتی.
- Kriging variance reduction: جایی که ایستگاه جدید بیشترین کاهش عدم قطعیت را می‌دهد.
- Entropy-based network design: حذف ایستگاه‌های تکراری و افزودن ایستگاه پرارزش.
- اولویت: ایستگاه‌های ثبات و تله‌متری در سرشاخه‌های تندابی (مادرسو، اوغان، چهل‌چای).
بهینه‌سازی فرایند (Automation)
- اجرای خودکار آزمون‌های Level 1 با هر دادهٔ جدید.
- صف بررسی کارشناس فقط برای مقادیر Suspect.
- ثبت هر اصلاح با دلیل، تاریخ و کاربر (audit trail).
- داشبورد کیفیت: درصد G/S/B/M هر ایستگاه.
بهینه‌سازی برای مدل (Fit for purpose)
- تمرکز QA/QC روی پنجره‌های ۱۶ رویداد تاریخی و ایستگاه‌های کلیدی آستانه.
- کمی‌سازی عدم قطعیت ورودی و انتقال آن به پیش‌بینی (uncertainty propagation).
- بازخورد مدل: باقیماندهٔ بزرگ مدل در یک ایستگاه = نشانهٔ احتمالی خطای داده.
خروجی‌های QA/QC در سامانه
QA report
گزارش هر ایستگاه: کامل بودن داده، پرچم‌ها، اصلاحات.
Station metadata
مختصات، ارتفاع، نوع دستگاه، تاریخچهٔ تغییرات.
Corrected series
سری‌های Level 3 و 4 با پرچم و عدم قطعیت.
QA layer on WebGIS
رنگ ایستگاه بر اساس کیفیت داده روی نقشه.
مراجع روش: WMO-No. 168 Guide to Hydrological Practices؛ WMO-No. 8 Guide to Instruments and Methods of Observation؛ WMO Manual on Stream Gauging (WMO-No. 1044). این روش در ماژول M0 (هاب داده و کنترل کیفیت) پیاده می‌شود؛ نقشهٔ اکوسیستم را ببینید.