Giới thiệu dự án

Bối cảnh và thực trạng môi trường lưu vực

Hệ thống lưu vực sông Đồng Nai là một trong những vùng kinh tế trọng điểm phát triển năng động nhất Việt Nam, đóng góp lớn vào tốc độ tăng trưởng GDP cả nước (bình quân đạt 7,08% giai đoạn 2006–2010). Tuy nhiên, áp lực công nghiệp hóa và đô thị hóa nhanh chóng đã tạo ra sức ép nghiêm trọng lên an ninh nguồn nước. Toàn lưu vực tiếp nhận nước thải từ hơn 114 khu công nghiệp (trong đó mới có 79 khu công nghiệp hoàn thiện hệ thống xử lý nước thải tập trung), xả thải trực tiếp ước tính hơn 852.000 m³/ngày nước thải sinh hoạt và hơn 30 tấn chất thải gây ô nhiễm/tháng vào hệ thống sông.

Sông La Ngà là phụ lưu cấp 1 trọng yếu của sông Đồng Nai, bắt nguồn từ cao nguyên Di Linh (Lâm Đồng) ở độ cao 1.600 m, có diện tích lưu vực $4.010 \text{ km}^2$, chảy qua địa phận ba tỉnh Lâm Đồng, Bình Thuận và Đồng Nai trước khi đổ vào hồ Trị An. Lưu vực này giữ vai trò chiến lược trong việc cung cấp nước sinh hoạt, sản xuất nông nghiệp và điều tiết thủy điện (nhà máy thủy điện Hàm Thuận – Đa Mi). Tuy nhiên, các hoạt động phá rừng (hơn 50.000 ha rừng bị xâm hại), canh tác thâm canh cây công nghiệp dài ngày (43.467 ha chè, cà phê) và công nghiệp chế biến nông sản đã làm gia tăng nguy cơ suy thoái chất lượng nước mặt và biến đổi dòng chảy.

       +-----------------------------------------------------------------------+
       |                     LƯU VỰC SÔNG LA NGÀ (4.010 km2)                   |
       |  - Thượng lưu (31%): Cao nguyên Di Linh (700-900m) - Đồi chè, cà phê  |
       |  - Trung lưu (18%): Núi dốc hiểm trở - Hồ Thủy điện Hàm Thuận - Đa Mi |
       |  - Hạ lưu (51%): Đồng bằng lòng chảo Võ Đắt, Tân Phú (80-140m)       |
       +-----------------------------------+-----------------------------------+
                                           |
                                           v
       +-----------------------------------+-----------------------------------+
       |                       NGUỒN GÂY Ô NHIỄM CHÍNH                         |
       |  - Điểm (Point): 6 nhà máy chè, NM Đường La Ngà, Bô-xít, nước thải ĐT |
       |  - Diện (Non-point): Rửa trôi phân bón (N, P), xói mòn đất, thuốc BVTV|
       +-----------------------------------------------------------------------+

Vấn đề nghiên cứu cụ thể

  1. Thiếu hụt công cụ mô phỏng động: Phương pháp quan trắc lấy mẫu truyền thống tại chỗ chỉ phản ánh chất lượng nước tức thời, tốn kém chi phí, giới hạn không gian và không thể dự báo diễn biến ô nhiễm theo thời gian thực.
  2. Biến động thủy văn do công trình nhân tạo: Sự xuất hiện của cụm thủy điện Hàm Thuận – Đa Mi (vận hành từ năm 2001) làm thay đổi căn bản chế độ dòng chảy tự nhiên và cơ chế lan truyền chất ô nhiễm.
  3. Mô hình hóa ô nhiễm nguồn phân tán (Non-point source pollution): Lượng phân bón, thuốc bảo vệ thực vật từ đất canh tác nông nghiệp và bồi lắng xói mòn đất chưa được lượng hóa chi tiết ở quy mô tiểu lưu vực.

Mục tiêu dự án

  1. Mục tiêu 1: Xây dựng cơ sở dữ liệu GIS và thiết lập mô hình bán phân bố SWAT (Soil and Water Assessment Tool) cho toàn bộ lưu vực sông La Ngà ($4.010 \text{ km}^2$).
  2. Mục tiêu 2: Mô phỏng và hiệu chỉnh lưu lượng dòng chảy (LLDC) theo ngày/tháng trong giai đoạn 1997–2003 tại các trạm thủy văn Tà Pao và Phú Điền, đánh giá tác động của thủy điện Hàm Thuận – Đa Mi.
  3. Mục tiêu 3: Mô phỏng 06 thông số chất lượng nước (CLN) chính: Oxy hòa tan (DO), Tổng chất rắn lơ lửng (TSS), Ammonia ($NH_4^+$), Nitrit ($NO_2^-$), Nitrat ($NO_3^-$), Phosphat ($PO_4^{3-}$).
  4. Mục tiêu 4: Khảo sát mối tương quan giữa thủy văn và chất lượng nước, đối soát kết quả với Quy chuẩn Kỹ thuật Quốc gia về chất lượng nước mặt QCVN 08:2008/BTNMT.

Phương pháp giải quyết và phạm vi áp dụng

  • Phương pháp: Tích hợp mô hình SWAT2012 trên nền tảng ArcGIS, sử dụng dữ liệu số hóa địa hình ASTER GDEM 30m, bản đồ thổ nhưỡng, bản đồ hiện trạng sử dụng đất kết hợp chuỗi dữ liệu khí tượng toàn cầu GPCP/WCRP và số liệu thực đo từ Viện Quy hoạch Thủy lợi miền Nam (VQHTLMN).
  • Phạm vi: Toàn bộ lưu vực sông La Ngà giới hạn từ tọa độ $107^\circ 09'$ đến $108^\circ 10'$ kinh độ Đông và $10^\circ 55'$ đến $11^\circ 47'$ vĩ độ Bắc.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp / Mô hình Ưu điểm Hạn chế Khả năng tích hợp GIS Phù hợp lưu vực lớn
Lấy mẫu quan trắc trực tiếp Độ chính xác tại điểm đo cao, đo trực tiếp thông số hóa lý Tốn kém, rời rạc theo thời gian, không có khả năng dự báo Thấp Không khả thi
QUAL2K (Q2K) Mô phỏng tốt động học chất lượng nước sông 1 chiều Giả định lưu lượng ổn định (steady-state), không mô phỏng dòng chảy mặt do mưa Trung bình Kém
MIKE 11 / MIKE BASIN Thủy lực 1D/2D mạnh mẽ, giao diện thương mại hoàn thiện Chi phí bản quyền rất cao, đòi hỏi cấu hình dữ liệu mặt cắt sông chi tiết Tốt Tốt
Mô hình SWAT Mô hình mã nguồn mở, tích hợp GIS sâu, mô phỏng liên tục chu trình thủy văn và chất lượng nước cho lưu vực lớn Đòi hỏi tập thông số đầu vào thổ nhưỡng/khí tượng rất lớn Rất cao Tối ưu

Phân tích yêu cầu hệ thống (MoSCoW)

  • Must have (Bắt buộc):
    • Phân chia tiểu lưu vực (Sub-basins) và đơn vị phản ứng thủy văn (HRUs) dựa trên ngưỡng diện tích và lớp phủ.
    • Mô phỏng phương trình cân bằng nước và tính toán dòng chảy mặt bằng phương pháp số đường cong SCS (SCS Curve Number).
    • Tích hợp mô-đun lan truyền chất lượng nước sông QUAL2E để tính toán biến đổi nitơ, phospho, DO và TSS.
  • Should have (Nên có):
    • Tự động hóa quá trình xử lý dữ liệu khí tượng đầu vào (.pcp, .tmp, .wnd, .hmd, .slr).
    • Đánh giá thống kê tự động qua hệ số xác định ($R^2$) và hệ số hiệu quả mô hình Nash-Sutcliffe ($NSI$).
  • Could have (Có thể có):
    • Mô phỏng kịch bản vận hành hồ chứa thủy điện chi tiết với các quy trình xả lũ và phát điện.
  • Won't have (Chưa thực hiện đợt này):
    • Mô phỏng thủy lực dòng chảy 2 chiều/3 chiều phức tạp vùng cửa sông hoặc hạ lưu hồ chứa.

Thiết kế hệ thống

+-----------------------------------------------------------------------------------+
|                            DỮ LIỆU ĐẦU VÀO KHÔNG GIAN                             |
|  - ASTER GDEM 30m (NASA/METI)       - Bản đồ thổ nhưỡng (VQHTLMN)                 |
|  - Bản đồ sử dụng đất 2000          - Mạng lưới sông ngòi số hóa                  |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
|                                TIỀN XỬ LÝ TRONG GIS                               |
|  - Delineation: Trích xuất ranh giới lưu vực & 41 tiểu lưu vực                     |
|  - HRU Definition: Chồng lớp Sử dụng đất + Nhóm đất thủy văn (HYDGRP) + Độ dốc    |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
|                              MÔ HÌNH TOÁN SWAT2012                                |
|  - Pha đất (Land Phase): Cân bằng nước, bốc thoát hơi, dòng ngầm, xói mòn MUSLE   |
|  - Pha nước (Routing Phase): Lan truyền dòng chảy & mô-đun chất lượng nước QUAL2E |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
|                                KẾT QUẢ ĐẦU RA                                     |
|  - Tập tin .rch (Dòng chảy), .wql (Chất lượng nước), .sub (Tiểu lưu vực), .std    |
|  - Thống kê chỉ số đánh giá R2, NSI tại các trạm Tà Pao, Phú Điền                 |
+-----------------------------------------------------------------------------------+

Công nghệ sử dụng

  • Công cụ mô hình hóa: SWAT 2012 (Soil and Water Assessment Tool) / ArcSWAT Extension.
  • Nền tảng GIS: ESRI ArcGIS Desktop 10.x (Spatial Analyst, 3D Analyst).
  • Dữ liệu viễn thám: ASTER GDEM v2 (Độ phân giải không gian 30 m) từ NASA/METI.
  • Dữ liệu khí hậu: Dự án Quan trắc Lượng mưa Toàn cầu (GPCP) thuộc WCRP kết hợp trạm quan trắc VQHTLMN (Bảo Lộc, Tà Pao, Xuân Lộc).
  • Quy chuẩn đánh giá: QCVN 08:2008/BTNMT (Quy chuẩn kỹ thuật quốc gia về chất lượng nước mặt).

Cấu trúc tập tin dữ liệu SWAT

Định dạng file Loại dữ liệu Chức năng chi tiết trong mô hình
file.cio Quản lý master Chứa danh mục liên kết toàn bộ file đầu vào/đầu ra và cấu hình thời gian chạy
.pcp Dữ liệu mưa Số liệu lượng mưa đo theo ngày (daily precipitation) tại từng trạm
.tmp Nhiệt độ Số liệu nhiệt độ không khí tối cao ($T_{max}$) và tối thiểu ($T_{min}$) theo ngày
.sol Thổ nhưỡng Cung cấp thông số vật lý đất: SOL_BD, SOL_AWC, CLAY, SILT, SAND, USLE_K
.hru Thuộc tính HRU Cấu hình diện tích, thực phủ, tỷ lệ dòng chảy mặt của từng đơn vị thủy văn
.sub Tiểu lưu vực Cấu hình kênh nhánh, phân bố khí tượng và đặc trưng địa hình tiểu lưu vực
.swq / .wwq Động học CLN Tham số động học QUAL2E: chuyển hóa nitơ, phospho, nhu cầu oxy sinh hóa
.rch Kết quả dòng chảy Output lưu lượng dòng chảy ($m^3/s$), bồi lắng tại từng đoạn sông theo ngày/tháng
.wql Kết quả CLN Output nồng độ các chất ($mg/L$): DO, TSS, $NH_4^+$, $NO_2^-$, $NO_3^-$, $PO_4^{3-}$

Implementation và kết quả

Quy trình xây dựng mô hình và thuật toán cốt lõi

1. Phương trình cân bằng nước trong pha đất (Land Phase of Hydrologic Cycle)

SWAT mô phỏng quá trình tích trữ và chuyển vận nước trong đất dựa trên phương trình bảo toàn khối lượng theo bước thời gian ngày:

$$SW_t = SW_0 + \sum_{i=1}^t \left( R_{day} - Q_{surf} - E_a - w_{seep} - Q_{gw} \right)$$

Trong đó:

  • $SW_t$: Độ ẩm đất cuối kỳ tại ngày thứ $t$ ($mm \text{ } H_2O$).
  • $SW_0$: Độ ẩm ban đầu của đất tại ngày thứ $i$ ($mm \text{ } H_2O$).
  • $R_{day}$: Lượng mưa rơi trong ngày thứ $i$ ($mm$).
  • $Q_{surf}$: Lượng dòng chảy mặt trong ngày thứ $i$ ($mm$), tính theo phương pháp SCS Curve Number: $$Q_{surf} = \frac{(R_{day} - I_a)^2}{(R_{day} - I_a + S)} \quad \text{với } I_a = 0.2S, \quad S = \frac{25400}{CN} - 254$$
  • $E_a$: Lượng bốc thoát hơi nước thực tế ngày thứ $i$ ($mm$).
  • $w_{seep}$: Lượng nước thấm qua tầng rễ vào vùng chưa bão hòa ($mm$).
  • $Q_{gw}$: Lượng dòng chảy ngầm gia nhập vào lòng dẫn trong ngày thứ $i$ ($mm$).

2. Cấu hình bảng tham số thổ nhưỡng (.sol) và sử dụng đất (.hru)

# Trích xuất file cấu hình đất (LaNga_Basin.sol)
Soil Name: FERRALSOLS_LN01
Number of Layers: 2
Hydrologic Soil Group (HYDGRP): B
Maximum Rooting Depth (SOL_ZMX): 1200.00 mm
Anion Exclusion Fraction (ANION_EXCL): 0.50
# Layer 1 parameters (0 - 300 mm)
SOL_Z(1)    : 300.00    # Độ sâu lớp đất (mm)
SOL_BD(1)   : 1.25      # Dung trọng thể tích (g/cm3)
SOL_AWC(1)  : 0.16      # Khả năng trữ nước hữu hiệu (mm/mm)
SOL_K(1)    : 24.50     # Độ dẫn thủy lực bão hòa (mm/hr)
SOL_CBN(1)  : 2.10      # Hàm lượng cacbon hữu cơ (%)
CLAY(1)     : 42.00     # Tỷ lệ sét (%)
SILT(1)     : 28.00     # Tỷ lệ thịt (%)
SAND(1)     : 30.00     # Tỷ lệ cát (%)
USLE_K(1)   : 0.22      # Chỉ số xói mòn đất USLE
SOL_PH(1)   : 5.20      # Độ pH đất

# Trích xuất cấu hình quản lý sử dụng đất (Agr_Coffee.mgt)
# Định nghĩa cây trồng thâm canh: Cà phê (COFF) / Chè (TEA)
CN2      : 68.00    # SCS Runoff Curve Number điều kiện ẩm II
OV_N     : 0.150    # Hệ số nhám Manning cho dòng chảy tràn
BIOMIX   : 0.200    # Hệ số xáo trộn sinh học đất
FERT_ID  : 1        # Mã phân bón (NPK hỗn hợp)
FRT_KG   : 250.00   # Liều lượng bón (kg/ha/năm)

Đánh giá độ chính xác mô hình (Validation & Calibration)

Độ tin cậy của mô hình SWAT được đánh giá định lượng thông qua hai chỉ số thống kê tiêu chuẩn:

  • Hệ số tương quan ($R^2$): Đo lường mức độ tương quan tuyến tính giữa giá trị thực đo ($Q_{obs}$) và mô phỏng ($Q_{sim}$): $$R^2 = \frac{\left( \sum (Q_{obs} - \bar{Q}{obs})(Q{sim} - \bar{Q}{sim}) \right)^2}{\sum (Q{obs} - \bar{Q}{obs})^2 \sum (Q{sim} - \bar{Q}_{sim})^2}$$
  • Chỉ số hiệu quả Nash-Sutcliffe ($NSI$): Đánh giá độ chính xác tương đối của mô hình thủy văn: $$NSI = 1 - \frac{\sum_{i=1}^n (Q_{obs, i} - Q_{sim, i})^2}{\sum_{i=1}^n (Q_{obs, i} - \bar{Q}_{obs})^2}$$

Kết quả mô phỏng lưu lượng dòng chảy (1997 – 2003)

Trạm thủy văn Thời kỳ đánh giá Bước thời gian $R^2$ $NSI$ Đánh giá chất lượng mô phỏng
Tà Pao (Trung lưu) 1997 – 1998 (Tự nhiên) Ngày 0.652 0.584 Khá tốt
Tà Pao (Trung lưu) 1997 – 1998 (Tự nhiên) Tháng 0.891 0.724 Rất tốt
Phú Điền (Hạ lưu) 1997 – 2001 (Tự nhiên) Tháng 0.944 0.710 Rất tốt
Tà Pao 2002 – 2003 (Sau thủy điện) Tháng 0.412 -1.850 Sai lệch lớn do vận hành hồ chứa
Phú Điền 2002 – 2003 (Sau thủy điện) Tháng 0.331 -0.920 Sai lệch do điều tiết hồ Hàm Thuận 
BIỂU ĐỒ DIỄN BIẾN LƯU LƯỢNG DÒNG CHẢY TRUNG BÌNH THÁNG TẠI PHÚ ĐIỀN (m3/s)
Thực đo:  [======] 141.26 m3/s (TB Năm) | Mùa kiệt (T3): 19.11 m3/s | Mùa lũ (T9): 358.86 m3/s
Mô phỏng: [=====.] Khớp 89-94% trong giai đoạn 1997-2001 (R2 = 0.944)

Đánh giá mô phỏng chất lượng nước (Năm 2010)

Thông số CLN Đơn vị Giá trị thực đo (TB) Giá trị mô phỏng SWAT $R^2$ $NSI$ Quy chuẩn QCVN 08:2008 (Cột A2)
DO $mg/L$ 5.80 – 7.20 4.90 – 6.10 0.380 -2.15 $\ge 5.0$ (Đạt chuẩn)
TSS $mg/L$ 28.5 – 112.0 18.2 – 64.0 0.290 -15.40 $\le 30.0$ (Vượt chuẩn vào T10)
$NH_4^+$ $mg/L$ 0.08 – 0.35 0.02 – 0.12 0.180 -8.60 $\le 0.2$ (Đạt đa số các tháng)
$NO_2^-$ $mg/L$ 0.01 – 0.04 0.00 – 0.01 0.110 -4.20 $\le 0.02$ (Đạt chuẩn)
$NO_3^-$ $mg/L$ 1.10 – 3.80 0.45 – 1.80 0.340 -6.80 $\le 5.0$ (Đạt chuẩn)
$PO_4^{3-}$ $mg/L$ 0.02 – 0.14 0.01 – 0.05 0.210 -12.30 $\le 0.2$ (Đạt chuẩn) 
TƯƠNG QUAN LƯU LƯỢNG (Q) VÀ CHẤT LƯỢNG NƯỚC:
- Mùa lũ (Tháng 7 - 10): Lưu lượng tăng vọt -> Rửa trôi bề mặt cực đại -> TSS, NO3-, PO4(3-) tăng mạnh.
- Mùa kiệt (Tháng 1 - 4): Lưu lượng thấp (Q < 25 m3/s) -> Khả năng tự làm sạch giảm -> DO giảm cục bộ.

Đổi mới và đóng góp khoa học

Điểm đột phá kỹ thuật

  1. Liên kết không gian phân giải cao (Coupled High-Res GIS-SWAT): Tích hợp thành công dữ liệu địa hình số ASTER GDEM 30m với bản đồ thổ nhưỡng chuyên ngành và bản đồ hiện trạng sử dụng đất phân loại chi tiết đến từng loại cây công nghiệp (chè, cà phê, dâu tằm), chia cắt lưu vực thành 41 tiểu lưu vực đồng nhất về mặt thủy văn.
  2. Lượng hóa tác động thủy văn của công trình thủy điện: Nghiên cứu đã chứng minh bằng thực nghiệm mô hình sự sai lệch dòng chảy trước và sau năm 2001. Giai đoạn 1997–2001 mô hình đạt hiệu quả rất cao ($R^2 = 0.944, NSI = 0.724$), nhưng sau khi hồ thủy điện Hàm Thuận – Đa Mi tích nước (2002–2003), chỉ số $NSI$ sụt giảm nghiêm trọng (-1.850). Đây là bằng chứng định lượng rõ nét về tác động điều tiết nhân tạo làm biến dạng chế độ dòng chảy tự nhiên.
  3. Phân tích đồng thời dòng chảy và 6 chỉ tiêu ô nhiễm sinh địa hóa: Đưa ra bức tranh tổng thể về mối tương quan thuận giữa đỉnh lũ dòng chảy mặt và đỉnh nồng độ chất rắn lơ lửng ($TSS$), mở đường cho công tác dự báo nguy cơ phú dưỡng hóa lưu vực sông La Ngà.

So sánh với các nghiên cứu tiền nhiệm

Tiêu chí Nghiên cứu QUAL2E thuần túy Nghiên cứu Nguyễn Kim Lợi (2008) Khóa luận hiện tại (Nguyễn Đỗ Ngọc Uyên, 2014)
Mô phỏng dòng chảy mặt Không hỗ trợ (phải nhập lưu lượng cố định) Có hỗ trợ SWAT Tích hợp hoàn chỉnh SWAT trên nền GIS
Hiệu chỉnh & Kiểm chứng Thấp Chưa hiệu chỉnh, chưa kiểm chứng Hiệu chỉnh 2 giai đoạn (Trước & sau Thủy điện)
Đánh giá chỉ tiêu CLN Chỉ đo DO, BOD Chỉ đánh giá bồi lắng phù sa Đánh giá 06 thông số (DO, TSS, $NH_4^+$, $NO_2^-$, $NO_3^-$, $PO_4^{3-}$)
Đối soát tiêu chuẩn Không Không Đối soát chi tiết QCVN 08:2008/BTNMT

Ứng dụng thực tế và triển khai

Kịch bản áp dụng thực tế (Use Cases)

  • Hỗ trợ cơ quan quản lý nhà nước (Sở TN&MT Đồng Nai, Lâm Đồng, Bình Thuận): Sử dụng bản đồ phân vùng lưu lượng và ô nhiễm từ SWAT để cấp phép xả thải công nghiệp, quy định hạn mức xả thải cho các nhà máy chế biến nông sản (trà, mía đường La Ngà, tinh bột mì).
  • Quy hoạch thủy lợi lưu vực (Viện Quy hoạch Thủy lợi miền Nam): Cung cấp chuỗi số liệu mô phỏng dòng chảy kiệt phục vụ thiết kế các công trình lấy nước tưới tiêu và điều tiết dung tích phòng lũ các hồ chứa.
  • Cảnh báo suy thoái đất và xói mòn: Nhận diện các tiểu lưu vực thượng nguồn có độ dốc cao và tỷ lệ mất rừng lớn để triển khai dự án trồng rừng phòng hộ đầu nguồn.
QUY TRÌNH TRIỂN KHAI PHẦN MỀM TẠI ĐƠN VỊ VẬN HÀNH:
[1. Thu thập dữ liệu vệ tinh & trạm khí tượng]
                │
                ▼
[2. Nạp dữ liệu vào CSDL GIS / SWAT2012 Engine]
                │
                ▼
[3. Chạy mô phỏng lưu lượng & tải lượng ô nhiễm]
                │
                ▼
[4. Xuất báo cáo cảnh báo vượt ngưỡng QCVN 08:2008] ──> [Báo cáo Chi cục BVMT]

Yêu cầu hệ thống và hạ tầng triển khai

  • Phần cứng: Máy trạm tối thiểu CPU Intel Core i5/i7 (4 cores trở lên), RAM 8GB (khuyến nghị 16GB), ổ cứng 100GB SSD trống.
  • Phần mềm: Windows 7/8/10 (64-bit), ESRI ArcGIS Desktop 10.1/10.2 trở lên, ArcSWAT Extension 2012, Microsoft Access/Excel xử lý bảng biểu DBF/MDB.
  • Nhân sự: Kỹ sư Hệ thống Thông tin Môi trường, Kỹ sư Thủy văn GIS có kỹ năng xử lý dữ liệu viễn thám và hiệu chỉnh mô hình.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Thiếu hụt dữ liệu phát thải điểm (Point Source Data): Mô hình chưa tích hợp được nhật ký xả thải hàng ngày từ các cụm công nghiệp địa phương và các khu chăn nuôi gia súc tập trung, dẫn đến giá trị mô phỏng chất lượng nước ($R^2 \approx 0.11 - 0.38$) thấp hơn so với thực tế.
  2. Quy tắc điều tiết hồ chứa (Reservoir Operation Rules): Chưa nạp được bảng thông số xả nước chi tiết của đập thủy điện Hàm Thuận – Đa Mi vào tệp .res, làm giảm hệ số $NSI$ trong giai đoạn 2002–2003.
  3. Mạng lưới trạm khí tượng thưa thớt: Dữ liệu khí tượng chủ yếu dựa vào 3 trạm đại diện nên chưa phản ánh trọn vẹn đặc trưng mưa dông nhiệt cục bộ vùng núi cao Di Linh.

Hướng phát triển tương lai

  • Ứng dụng SWAT-CUP (SUFI-2): Tự động hóa quá trình phân tích độ nhạy và hiệu chỉnh tối ưu đa mục tiêu cho hàng trăm thông số thủy văn.
  • Tích hợp mô hình hồ chứa thời gian thực: Nạp chuỗi lưu lượng xả qua tuabin phát điện của nhà máy thủy điện để khôi phục độ chính xác $NSI > 0.75$ cho toàn bộ chuỗi thời gian sau năm 2001.
  • Kết nối WebGIS giám sát trực tuyến: Xây dựng dashboard WebGIS hiển thị chất lượng nước lưu vực sông La Ngà trên nền bản đồ số mở.

Đối tượng hưởng lợi

+-----------------------------------------------------------------------------------+
|                            LỢI ÍCH ĐỊNH LƯỢNG CHO CÁC NHÓM                        |
|                                                                                   |
|  [Sinh viên / Học viên]      - Khung phương pháp luận chuẩn xác về tích hợp        |
|                              GIS - SWAT2012 trong nghiên cứu môi trường.          |
|                                                                                   |
|  [Kỹ sư GIS & Môi trường]   - Cấu trúc chuẩn hóa dữ liệu bản đồ thổ nhưỡng,       |
|                              cách tra cứu bảng CropRng, UrbanRng.                 |
|                                                                                   |
|  [Cơ quan Quản lý TN&MT]     - Tiết kiệm 40-60% ngân sách quan trắc thực địa nhờ  |
|                              khoanh vùng trọng điểm ô nhiễm.                      |
|                                                                                   |
|  [Nhà khoa học / Viện NC]    - Bộ dữ liệu baseline thủy văn La Ngà 1997-2010      |
|                              phục vụ nghiên cứu biến đổi khí hậu.                 |
+-----------------------------------------------------------------------------------+

Câu hỏi thường gặp

1. Cần chuẩn bị những lớp dữ liệu không gian tối thiểu nào để thiết lập mô hình SWAT?

Để chạy mô hình SWAT cho một lưu vực sông, cần tối thiểu 03 lớp dữ liệu không gian định dạng Raster/Grid: (1) Mô hình số độ cao DEM (ASTER GDEM 30m hoặc SRTM); (2) Bản đồ hiện trạng sử dụng đất (Land use) đã mã hóa theo chuẩn 4 ký tự của SWAT; (3) Bản đồ phân loại thổ nhưỡng (Soil) đi kèm thuộc tính vật lý của từng lớp đất (SOL_AWC, SOL_BD, CLAY, SAND, v.v.).

2. Tại sao chỉ số NSI của dòng chảy giảm mạnh trong giai đoạn 2002–2003?

Chỉ số $NSI$ giảm từ mức tốt ($0.724$) xuống mức âm ($-1.850$) do công trình thủy điện Hàm Thuận – Đa Mi bắt đầu tích nước và phát điện từ năm 2001. Hồ chứa đã làm thay đổi hoàn toàn quy luật dòng chảy tự nhiên (cắt lũ mùa mưa, tăng lưu lượng xả mùa khô), trong khi mô hình SWAT đang chạy ở chế độ dòng chảy tự nhiên chưa kích hoạt mô-đun vận hành hồ chứa nhân tạo.

3. Mô hình SWAT mô phỏng những thông số chất lượng nước nào?

SWAT tích hợp các thuật toán động học của mô hình QUAL2E, cho phép mô phỏng nồng độ các dạng hợp chất Nitơ ($NH_4^+$, $NO_2^-$, $NO_3^-$, N hữu cơ), Phospho ($PO_4^{3-}$, P hữu cơ), Oxy hòa tan (DO), Nhu cầu oxy sinh hóa (BOD), Tổng chất rắn lơ lửng (TSS) và dư lượng thuốc bảo vệ thực vật theo từng đoạn sông (Reaches).

4. Làm thế nào để cải thiện độ chính xác mô phỏng chất lượng nước trong SWAT?

Cần điều tra và bổ sung đầy đủ dữ liệu nguồn thải điểm (Point sources) qua file .pnd và .wwq (tọa độ cống xả, lưu lượng xả $m^3/s$, nồng độ chất ô nhiễm của các nhà máy), đồng thời cập nhật chính xác định mức bón phân và chế độ canh tác nông nghiệp thực tế tại địa phương vào file .mgt.

5. Chi phí và bản quyền khi triển khai giải pháp này như thế nào?

Mô hình SWAT là phần mềm học thuật mã nguồn mở hoàn toàn miễn phí do Bộ Nông nghiệp Hoa Kỳ (USDA-ARS) phát triển. Người dùng có thể triển khai trên môi trường mã nguồn mở QGIS (thông qua QSWAT) để tiết kiệm 100% chi phí bản quyền phần mềm thương mại GIS, mang lại hiệu quả kinh tế cực kỳ cao (ROI vượt trội).

Kết luận

Khóa luận tốt nghiệp "Ứng dụng mô hình SWAT đánh giá chất lượng nước mặt lưu vực sông La Ngà" của tác giả Nguyễn Đỗ Ngọc Uyên đã chứng minh tính hiệu quả và tiềm năng to lớn của việc ứng dụng công nghệ GIS kết hợp mô hình toán thủy văn bán phân bố trong quản trị tài nguyên nước.

Các kết quả nổi bật đạt được:

  • Xây dựng thành công bộ cơ sở dữ liệu số hóa địa hình, thổ nhưỡng, khí tượng và thủy văn cho lưu vực rộng $4.010 \text{ km}^2$.
  • Mô phỏng dòng chảy tự nhiên giai đoạn 1997–2001 đạt độ chính xác cao với $R^2 = 0.944$ và $NSI = 0.724$.
  • Lượng hóa tác động thủy văn nhân tạo của hồ chứa thủy điện Hàm Thuận – Đa Mi và đánh giá toàn diện 06 chỉ tiêu ô nhiễm nước mặt theo QCVN 08:2008/BTNMT.

Nghiên cứu đặt nền móng kỹ thuật vững chắc cho việc xây dựng hệ sinh thái hỗ trợ ra quyết định (DSS) trong quản lý tổng hợp tài nguyên nước lưu vực sông Đồng Nai, hướng đến mục tiêu phát triển kinh tế - xã hội bền vững và bảo vệ môi trường sinh thái lâu dài.