Giới thiệu dự án

Lưu vực sông Đồng Nai là một trong những vùng kinh tế trọng điểm phía Nam, nơi tiếp nhận tải lượng nước thải sinh hoạt và công nghiệp rất lớn. Toàn lưu vực hiện có hơn 114 khu công nghiệp nhưng chỉ khoảng 70% có trạm xử lý nước thải tập trung đạt chuẩn, xả thải trực tiếp ra môi trường khoảng 30 tấn chất ô nhiễm/tháng (gồm dầu mỡ, chất hữu cơ, kim loại nặng). Trong đó, lưu vực sông La Ngà (diện tích $4.010\text{ km}^2$) – phụ lưu cấp I quan trọng bắt nguồn từ cao nguyên Di Linh (Lâm Đồng) chảy qua Bình Thuận và Đồng Nai – đang đối mặt với áp lực ô nhiễm nghiêm trọng do khai thác đất nông nghiệp thâm canh, chế biến nông sản (chè, cà phê, dâu tằm), công nghiệp mía đường và khai thác bô-xít.

                  +----------------------------------------------+
                  |  Cao nguyên Di Linh (Lâm Đồng) - 1.280 km²   |
                  |  - Nông nghiệp: Cà phê, Trà, Dâu tằm (86.5%)  |
                  |  - Công nghiệp chế biến chè, bô-xít          |
                  +-----------------------+----------------------+
                                          |
                                          v
                  +----------------------------------------------+
                  |  Vùng trung lưu (Tà Pao) - 720 km²           |
                  |  - Thủy điện Hàm Thuận - Đa Mi (từ 2001)     |
                  |  - Rừng núi dốc, địa hình hiểm trở           |
                  +-----------------------+----------------------+
                                          |
                                          v
                  +----------------------------------------------+
                  |  Vùng hạ lưu (Đồng Nai) - 2.100 km²          |
                  |  - Vùng lúa, cây ngắn ngày, Nhà máy đường LN |
                  |  - Đổ vào hồ Trị An / Hệ thống Sông Đồng Nai |
                  +----------------------------------------------+

Vấn đề thực tế (Problem Statement)

Phương pháp đánh giá chất lượng nước truyền thống phụ thuộc vào việc lấy mẫu định kỳ tại các vị trí cố định. Cách tiếp cận này bộc lộ ba điểm nghẽn lớn:

  • Chi phí nhân lực và hóa chất phân tích cao, chu kỳ lấy mẫu ngắt quãng (thường 1 - 3 tháng/lần) không phản ánh được tính biến thiên liên tục của dòng chảy.
  • Không có khả năng dự báo lan truyền ô nhiễm theo không gian và thời gian thực dưới tác động của biến đổi khí hậu và chuyển đổi sử dụng đất.
  • Thiếu công cụ định lượng đóng góp của nguồn ô nhiễm diện (nông nghiệp, xói mòn đất, thuốc bảo vệ thực vật) so với nguồn ô nhiễm điểm (nước thải công nghiệp, sinh hoạt đô thị).

Mục tiêu dự án

  1. Xây dựng cơ sở dữ liệu không gian và thuộc tính tích hợp GIS cho lưu vực sông La Ngà (độ cao số DEM 30m, thổ nhưỡng, hiện trạng sử dụng đất LULC năm 2000, khí tượng 1997–2010).
  2. Thiết lập và mô phỏng dòng chảy mặt và tải lượng chất lượng nước bằng mô hình bán phân bố SWAT (Soil and Water Assessment Tool).
  3. Đánh giá độ chính xác của mô hình thủy văn theo chỉ số Nash-Sutcliffe ($NSI$) và hệ số xác định ($R^2$) trong hai giai đoạn: tiền thủy điện (1997–2001) và hậu thủy điện Hàm Thuận – Đa Mi (2002–2003).
  4. Phân tích tương quan giữa lưu lượng dòng chảy (LLDC) và 6 thông số chất lượng nước (CLN): Oxy hòa tan (DO), Tổng chất rắn lơ lửng (TSS), Ammonia ($\text{NH}_4^+$), Nitrit ($\text{NO}_2^-$), Nitrat ($\text{NO}_3^-$), Phosphat ($\text{PO}_4^{3-}$).
  5. Đánh giá mức độ phù hợp của chất lượng nước mặt lưu vực theo Quy chuẩn Kỹ thuật Quốc gia QCVN 08:2008/BTNMT.

Phương pháp tiếp cận và Phạm vi nghiên cứu

Dự án áp dụng mô hình toán bán phân bố SWAT tích hợp GIS để mô phỏng chu trình thủy văn và chất lượng nước. Không gian nghiên cứu bao quát $4.010\text{ km}^2$ thuộc địa bàn 3 tỉnh Lâm Đồng, Bình Thuận và Đồng Nai, chia thành 3 tiểu vùng: thượng lưu ($1.280\text{ km}^2$), trung lưu ($720\text{ km}^2$) và hạ lưu ($2.100\text{ km}^2$). Chuỗi dữ liệu mô phỏng thủy văn vận hành từ năm 1997 đến 2003 và hiệu chỉnh chất lượng nước với chuỗi quan trắc thực đo năm 2010.

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

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

Hiện nay có nhiều mô hình toán thủy lực và chất lượng nước đang được ứng dụng. Bảng so sánh dưới đây phân tích ưu - nhược điểm của các mô hình tiêu biểu:

Tiêu chí SWAT (Soil and Water Assessment Tool) QUAL2K MIKE 11 / MIKE BASIN WASP7 (EPA)
Loại mô hình Bán phân bố (Semi-distributed), liên tục theo thời gian 1D dừng/bán động trên trục sông 1D thủy lực động lực học đầy đủ Động lực học chất lượng nước (1D, 2D, 3D)
Nguồn ô nhiễm diện Mô phỏng chi tiết theo từng loại đất, cây trồng, phân bón Hạn chế, chỉ nhập tải lượng biên đơn giản Cần ghép nối thêm mô hình thu nước mặt Không tích hợp sẵn tính toán dòng chảy mặt
Tích hợp GIS Tự động hóa qua ArcGIS/QGIS (ArcSWAT, QSWAT) Nhập thủ công qua bảng tính Excel Tích hợp sâu vào bộ phần mềm DHI Ghép nối ngoài thông qua file cấu hình
Chi phí bản quyền Mã nguồn mở, Miễn phí Miễn phí Bản quyền thương mại rất cao Miễn phí (EPA)
Khả năng ứng dụng quy mô lớn Tối ưu cho lưu vực lớn ($>1.000\text{ km}^2$) Chỉ phù hợp đoạn sông cụ thể Phù hợp thủy lực kênh nhánh phức tạp Phù hợp vùng nước đứng, hồ chứa, cửa sông

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

  • Must have: Phân chia lưu vực tự động từ DEM 30m; tạo lập các đơn vị phản ứng thủy văn (HRUs); mô phỏng cân bằng nước hàng ngày/tháng; tính toán chỉ số thống kê $R^2$ và $NSI$.
  • Should have: Mô phỏng sự biến thiên nồng độ 6 chỉ tiêu CLN (DO, TSS, $\text{NH}_4^+$, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$); tự động gán dữ liệu khí tượng từ các trạm lân cận.
  • Could have: Mô phỏng biến kịch bản điều tiết hồ chứa thủy điện Hàm Thuận – Đa Mi; đánh giá rủi ro xói mòn đất theo công thức MUSLE.
  • Won't have (trong giai đoạn này): Mô phỏng lan truyền thủy lực 2D dòng xoáy ven bờ; tích hợp dữ liệu xả thải SCADA thời gian thực từ các nhà máy.

Thiết kế hệ thống và Kiến trúc dữ liệu

graph TD
    A[ASTER GDEM 30m] -->|Khai thác mạng thủy văn| E[Phân chia Tiểu lưu vực - Subbasins]
    B[Bản đồ Thổ nhưỡng FAO/VQHTLMN] --> F[Chồng lớp không gian]
    C[Bản đồ LULC 2000] --> F
    D[Bản đồ Độ dốc Slope] --> F
    F -->|Ngưỡng diện tích 10%/10%/10%| G[Đơn vị phản ứng thủy văn - HRUs]
    E --> G
    H[Dữ liệu khí tượng 1997-2010: Mưa, Nhiệt độ, Gió, Bức xạ, Độ ẩm] --> I[SWAT Core Engine - Water Balance]
    G --> I
    I --> J[Pha đất: Chu trình nước, xói mòn MUSLE, dinh dưỡng]
    J --> K[Pha nước: Dẫn truyền kênh chính QUAL2E]
    K --> L[Đầu ra: Lưu lượng dòng chảy & Tải lượng CLN]
    L --> M[Hiệu chỉnh / Kiểm định: R², NSI với Dữ liệu thực đo]

Technology Stack

  • GIS Platform: ESRI ArcGIS 10.2 / ArcSWAT 2012 Engine.
  • Mô hình cốt lõi: SWAT 2012 Revision 622 (viết bằng Fortran/C, tích hợp thuật toán QUAL2E để tính toán lan truyền chất lượng nước trong lòng kênh).
  • Dữ liệu không gian: ASTER Global Digital Elevation Model (GDEM v2) độ phân giải 30m; Bản đồ hiện trạng sử dụng đất tỷ lệ 1:100.000 (VQHTLMN, 2000); Bản đồ đất phân loại theo FAO/UNESCO.
  • Công cụ thống kê & tiền xử lý: SWAT-CUP 5.1.6.2, Python 2.7 (ArcPy, NumPy, Pandas), Microsoft Access Driver (SWAT2012.mdb).

Cấu trúc dữ liệu đầu vào và đầu ra SWAT

Mô hình quản lý cấu trúc file văn bản theo định dạng cố định:

  • .cio: Master control file quy định toàn bộ tiến trình đọc/ghi của lưu vực.
  • .pcp, .tmp, .slr, .wnd, .hmd: Tập hợp chuỗi dữ liệu chuỗi thời gian khí tượng theo ngày của các trạm Bảo Lộc, Tà Pao, Xuân Lộc.
  • .sol: Bảng thuộc tính vật lý của 12 loại đất (Dung trọng SOL_BD, Độ dẫn thủy lực bão hòa SOL_K, Suất phản chiếu SOL_ALB, Hàm lượng sét CLAY, cát SAND).
  • .sub, .hru: Định nghĩa thuộc tính hình học và tham số quản lý canh tác của từng đơn vị thủy văn.
  • .rch, .wql: Tệp kết quả đầu ra mô phỏng lưu lượng dòng chảy và chất lượng nước tại các cửa xả lưu vực theo từng bước thời gian.

Implementation và kết quả

Quy trình tính toán và Thuật toán cốt lõi

1. Phương trình cân bằng nước tổng thể (Pha đất)

Toàn bộ chu trình thủy lực trong mô hình SWAT được thiết lập dựa trên nguyên lý bảo toàn khối lượng hàng 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 của đất tại ngày thứ $t$ ($\text{mm }\text{H}_2\text{O}$).
  • $SW_0$: Độ ẩm ban đầu của đất tại ngày thứ $i$ ($\text{mm }\text{H}_2\text{O}$).
  • $R_{day}$: Lượng mưa rơi trong ngày thứ $i$ ($\text{mm}$).
  • $Q_{surf}$: Lượng dòng chảy mặt sinh ra trong ngày thứ $i$ ($\text{mm}$), tính theo phương pháp SCS Curve Number (CN-II):

$$Q_{surf} = \frac{(R_{day} - I_a)^2}{(R_{day} - I_a + S)} = \frac{(R_{day} - 0.2S)^2}{(R_{day} + 0.8S)} \quad (\text{với } R_{day} > 0.2S)$$

$$S = 25.4 \left( \frac{1000}{CN} - 10 \right)$$

  • $E_a$: Lượng bốc thoát hơi nước thực tế ($\text{mm}$).
  • $w_{seep}$: Lượng nước thấm qua tầng rễ xuống vùng chưa bão hòa ($\text{mm}$).
  • $Q_{gw}$: Lượng nước ngầm chảy bổ cập vào dòng chính ($\text{mm}$).

2. Chỉ số đánh giá độ tin cậy mô hình

Độ chính xác giữa giá trị mô phỏng ($S_i$) và thực đo ($O_i$) được đánh giá qua hai chỉ số:

$$R^2 = \frac{\left[ \sum_{i=1}^n (O_i - \bar{O})(S_i - \bar{S}) \right]^2}{\sum_{i=1}^n (O_i - \bar{O})^2 \sum_{i=1}^n (S_i - \bar{S})^2}$$

$$NSI = 1 - \frac{\sum_{i=1}^n (O_i - S_i)^2}{\sum_{i=1}^n (O_i - \bar{O})^2}$$

import numpy as np

def calculate_model_performance(observed, simulated):
    """
    Tinh toan chi so R2 va Nash-Sutcliffe Efficiency (NSI) giua thuc do va mo phong.
    """
    obs = np.array(observed)
    sim = np.array(simulated)
    
    # Loai bo cac diem du lieu Missing (-99)
    mask = (obs != -99) & (sim != -99)
    obs, sim = obs[mask], sim[mask]
    
    mean_obs = np.mean(obs)
    mean_sim = np.mean(sim)
    
    # He so xac dinh R2
    numerator_r2 = np.sum((obs - mean_obs) * (sim - mean_sim)) ** 2
    denominator_r2 = np.sum((obs - mean_obs) ** 2) * np.sum((sim - mean_sim) ** 2)
    r2 = numerator_r2 / denominator_r2 if denominator_r2 != 0 else 0
    
    # Chi so Nash-Sutcliffe (NSI)
    numerator_nsi = np.sum((obs - sim) ** 2)
    denominator_nsi = np.sum((obs - mean_obs) ** 2)
    nsi = 1 - (numerator_nsi / denominator_nsi) if denominator_nsi != 0 else -np.inf
    
    return {"R2": round(r2, 4), "NSI": round(nsi, 4)}

# Vi du du lieu dong chay thuc do vs mo phong tai Tram Ta Pao (m3/s)
obs_flow = [24.81, 15.14, 12.56, 18.11, 30.77, 64.93, 112.62, 182.96, 181.47, 165.38, 95.01, 46.21]
sim_flow = [26.10, 14.80, 11.20, 19.50, 35.40, 72.10, 125.80, 198.30, 192.50, 174.10, 88.30, 42.10]

metrics = calculate_model_performance(obs_flow, sim_flow)
print(f"Chi so Kiem dinh Thuy van: R2 = {metrics['R2']}, NSI = {metrics['NSI']}")

Kết quả kiểm định và Đánh giá mô hình

1. Mô phỏng Lưu lượng dòng chảy (LLDC) giai đoạn 1997 – 2003

Mô hình được chạy mô phỏng tại 2 trạm thủy văn then chốt: Tà Pao (Trung lưu) và Phú Điền (Hạ lưu).

Trạm thủy văn Thời kỳ Bước thời gian $R^2$ $NSI$ Nhận xét độ tin cậy
Tà Pao 1997 – 1998 (Tiền thủy điện) Tháng 0,782 – 0,944 0,542 – 0,724 Rất tốt; Dòng chảy mô phỏng bám sát thực đo
Tà Pao 1997 – 2001 Ngày 0,412 – 0,685 0,215 – 0,480 Khá; Mô phỏng ngày chịu ảnh hưởng cục bộ của mưa dông
Phú Điền 1997 – 1998 Tháng 0,810 – 0,915 0,601 – 0,688 Tốt; Cân bằng nước chuẩn xác ở quy mô hạ lưu
Tà Pao & Phú Điền 2002 – 2003 (Hậu thủy điện) Tháng 0,331 – 0,512 -1,42 đến 0,08 Kém; Thủy điện Hàm Thuận - Đa Mi tích nước làm đảo lộn dòng chảy tự nhiên

[!IMPORTANT] Sự suy giảm đột ngột của hệ số $NSI$ trong giai đoạn 2002–2003 (xuống mức âm) phản ánh chính xác tác động can thiệp nhân tạo của công trình thủy điện Hàm Thuận – Đa Mi (vận hành từ năm 2001). Hồ chứa đã cắt giảm đỉnh lũ mùa mưa và xả nước tăng cường vào mùa khô, làm chuỗi số liệu tự nhiên không còn tương thích nếu không có mô-đun điều tiết hồ chứa (.res).

2. Mô phỏng Chất lượng nước (CLN) năm 2010

Mô phỏng nồng độ chất lượng nước được kiểm chứng tại 2 vị trí quan trắc SW_LN_01 (Tà Pao) và SW_LN_02 (Phú Điền) cho 6 thông số:

Thông số CLN $R^2$ $NSI$ Nồng độ mô phỏng trung bình Nồng độ thực đo trung bình So sánh QCVN 08:2008 (Cột A1/B1)
DO (mg/L) 0,38 -2,14 5,82 6,45 Đạt chuẩn Cột A1 ($\ge 6\text{ mg/L}$)
TSS (mg/L) 0,41 -12,85 32,4 68,2 Vượt chuẩn vào tháng X ($128\text{ mg/L} > 50\text{ mg/L}$)
$\text{NH}_4^+$ (mg/L) 0,15 -45,10 0,08 0,22 Đạt mức A2/B1 ($\le 0,5\text{ mg/L}$)
$\text{NO}_2^-$ (mg/L) 0,08 -88,20 0,005 0,018 Nằm trong ngưỡng cho phép
$\text{NO}_3^-$ (mg/L) 0,22 -15,40 0,45 1,85 Đạt chuẩn Cột A1 ($\le 2\text{ mg/L}$)
$\text{PO}_4^{3-}$ (mg/L) 0,12 -188,40 0,02 0,11 Đạt chuẩn Cột A2 ($\le 0,2\text{ mg/L}$)

Nhận định sai số chất lượng nước: Giá trị mô phỏng luôn thấp hơn thực tế do dữ liệu đầu vào chưa tích hợp đầy đủ các điểm xả thải tập trung (Point Sources) từ nhà máy chế biến nông sản và hệ số rửa trôi phân bón hóa học trong nông nghiệp còn bị ước lượng thấp.

Đổi mới và đóng góp

Điểm cải tiến kỹ thuật

  • Phân giải không gian nâng cao: Xây dựng mô hình số độ cao DEM 30m chuẩn hóa từ dữ liệu ASTER GDEM, tối ưu hóa việc phân chia lưu vực La Ngà thành 37 tiểu lưu vực và 245 đơn vị phản ứng thủy văn (HRUs) đồng nhất về loại đất, độ dốc và thảm phủ.
  • Định lượng tác động của hồ chứa: Cung cấp bằng chứng định lượng rõ nét về sự thay đổi thủy văn trước và sau khi xây dựng công trình thủy điện Hàm Thuận – Đa Mi thông qua sự biến thiên các chỉ số hiệu quả $NSI$ và $R^2$.
  • Mô hình hóa mối quan hệ Dòng chảy – Ô nhiễm: Xác lập được hàm tương quan thuận giữa lưu lượng dòng chảy và tải lượng chất rắn lơ lửng ($TSS$), nitrat và phosphat, chứng minh cơ chế rửa trôi bề mặt do dòng chảy tràn đóng vai trò chủ đạo trong việc tải chất ô nhiễm vào mùa mưa lũ.

So sánh với các công trình nghiên cứu trước đây

Đặc điểm nghiên cứu Nghiên cứu của Trần Xuân Lộc (2012) Nghiên cứu của Nguyễn Hà Trang (2009) Nghiên cứu này (Nguyễn Đỗ Ngọc Uyên, 2014)
Đối tượng Hồ Cầu Mới (Quy mô nhỏ) Toàn bộ lưu vực sông Đồng Nai Lưu vực sông La Ngà ($4.010\text{ km}^2$)
Chỉ tiêu chất lượng nước Bồi lắng trầm tích Đánh giá sơ bộ BOD, COD Chi tiết 6 chỉ tiêu: DO, TSS, $\text{NH}_4^+$, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$
Kiểm định mô hình Hiệu chỉnh hạn chế Sai số lớn, chưa kiểm định sâu Kiểm định riêng biệt 2 giai đoạn (trước/sau thủy điện)
Độ tin cậy thủy văn $R^2 \approx 0,65$ $R^2 \approx 0,58$ $R^2$ tháng đạt tới 0,944 (giai đoạn 1997–1998)

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

Kịch bản ứng dụng trong quản lý tài nguyên

  1. Quy hoạch phân vùng xả thải: Giúp Sở Tài nguyên và Môi trường tỉnh Đồng Nai, Bình Thuận xác định sức chịu tải môi trường nước mặt trên từng tiểu lưu vực, từ đó cấp phép lưu lượng xả thải tối đa cho các nhà máy công nghiệp.
  2. Cảnh báo suy thoái đất và xói mòn: Nhận diện các khu vực có độ dốc $>20%$ tại thượng lưu Bảo Lộc bị mất rừng để ưu tiên triển khai các dự án trồng rừng phòng hộ và nông - lâm kết hợp.
  3. Vận hành liên hồ chứa thông minh: Cung cấp chuỗi lưu lượng dòng chảy biên đầu vào cho các kịch bản xả lũ và điều tiết tưới tiêu của đập Hàm Thuận – Đa Mi và đập dâng Tà Pao.
                   +-----------------------------------------------+
                   |     SWAT Cloud & GIS Spatial Server           |
                   |  (Cơ sở dữ liệu DEM, Khí tượng, Sử dụng đất)  |
                   +-----------------------+-----------------------+
                                           |
                    +----------------------+----------------------+
                    |                                             |
                    v                                             v
+---------------------------------------+     +---------------------------------------+
|   Sở TN&MT Đồng Nai / Bình Thuận      |     |  Ban Quản lý Thủy điện / Thủy lợi     |
|   - Giám sát chỉ số DO, TSS, Dinh dưỡng |     |  - Dự báo lưu lượng đến hồ            |
|   - Cấp phép và thanh tra nguồn xả    |     |  - Điều tiết lưu lượng xả dòng chảy   |
+---------------------------------------+     +---------------------------------------+

Yêu cầu cấu hình triển khai

  • Hệ điều hành: Windows Server / Linux Ubuntu 20.04 LTS (chạy SWAT core binary).
  • Phần cứng đề xuất: CPU Intel Xeon hoặc Core i7 đa nhân (tối thiểu 8 cores để chạy mô phỏng Monte Carlo đa kịch bản); RAM 16GB; Ổ cứng SSD NVMe 256GB.
  • Môi trường phần mềm: ArcGIS Desktop 10.x kèm ArcSWAT extension hoặc QGIS 3.x với plugin QSWAT, Python 3.8+ hỗ trợ tự động hóa tiền/hậu xử lý dữ liệu.

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

Hạn chế kỹ thuật

  • Thiếu hụt dữ liệu nguồn thải phân tán: Số liệu phân bón sử dụng trong canh tác nông nghiệp và nước thải sinh hoạt nông thôn chưa được điều tra chi tiết, dẫn đến việc mô hình hóa chất lượng nước có sai số âm tính ($NSI < 0$).
  • Gián đoạn chuỗi khí tượng: Một số trạm đo mưa bị gián đoạn, phải bù số liệu thông qua dự án GPCP toàn cầu có sai số nhất định so với vi khí hậu địa phương.
  • Chưa tích hợp mô-đun vận hành hồ chứa: Chưa đưa trực tiếp biểu đồ xả nước hàng ngày của nhà máy thủy điện Hàm Thuận – Đa Mi vào tệp .res, làm giảm tính chính xác của mô hình sau năm 2001.

Hướng phát triển tiếp theo

  • Hiệu chỉnh tự động nâng cao: Ứng dụng thuật toán SUFI-2 trong phần mềm SWAT-CUP để phân tích độ nhạy (Sensitivity Analysis) và hiệu chỉnh tự động 16 tham số dòng chảy (CN2, ALPHA_BF, GW_DELAY, ESCO, v.v.).
  • Tích hợp ảnh viễn thám độ phân giải cao: Sử dụng ảnh vệ tinh Sentinel-2/Landsat 8 để cập nhật biến động lớp phủ đất (LULC) theo từng mùa vụ thay vì dùng bản đồ tĩnh.
  • Ghép nối SWAT-MODFLOW: Mở rộng mô phỏng tương tác 2 chiều giữa tầng nước mặt sông La Ngà và mực nước ngầm khu vực thềm phù sa cổ Định Quán – Tân Phú.

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

                                  ĐỐI TƯỢNG HƯỞNG LỢI
                                           |
         +-----------------+---------------+----------------+----------------+
         |                 |                                |                |
         v                 v                                v                v
   [SINH VIÊN]        [KỸ SƯ GIS]                     [CƠ QUAN QUẢN LÝ]  [CỘNG ĐỒNG & DN]
   - Mẫu khóa luận    - Pipeline xử lý DEM/LULC        - Dữ liệu quy hoạch- Bảo vệ nguồn nước
   - Quy trình SWAT   - Script Python kiểm định NSI    - Báo cáo QCVN 08  - Giảm thiệt hại lũ
  • Sinh viên & Học viên cao học: Tiếp cận toàn bộ phương pháp luận chuẩn hóa về mô hình hóa thủy văn, quy trình xử lý dữ liệu viễn thám/GIS và script tính toán chỉ số thống kê $R^2$, $NSI$.
  • Kỹ sư GIS & Lập trình viên Môi trường: Nắm bắt cấu trúc tệp dữ liệu SWAT, logic kết nối CSDL không gian và giải thuật tự động hóa mô phỏng.
  • Cơ quan quản lý nhà nước (Sở TN&MT, Viện Quy hoạch Thủy lợi): Sở hữu bộ công cụ khoa học phục vụ công tác giám sát, đánh giá tải lượng ô nhiễm và lập quy hoạch tài nguyên nước mặt lưu vực sông La Ngà.
  • Doanh nghiệp nông nghiệp & thủy điện: Dự báo được quy luật phân phối dòng chảy kiệt và dòng chảy lũ, tối ưu hóa kế hoạch trữ nước và vận hành phát điện.

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

1. Yêu cầu kỹ thuật tối thiểu để chạy mô hình SWAT cho lưu vực quy mô $4.000\text{ km}^2$?

Để thiết lập mô hình ArcSWAT cho lưu vực khoảng $4.000\text{ km}^2$ với độ phân giải DEM 30m (khoảng 37 tiểu lưu vực, 250 HRUs), hệ thống cần tối thiểu CPU 4 nhân (2.5 GHz trở lên), 8GB RAM và phần mềm ArcGIS Desktop (bản 10.2 - 10.8) kèm ArcSWAT tương thích. Thời gian chạy một kịch bản 10 năm mất khoảng 30–60 giây.

2. Vì sao chỉ số $NSI$ của mô hình dòng chảy lại giảm sâu sau năm 2001?

Năm 2001, nhà máy thủy điện Hàm Thuận – Đa Mi chính thức tích nước và phát điện. Việc hồ chứa điều tiết dòng chảy tự nhiên (cắt lũ mùa mưa, xả đáy mùa khô) làm thay đổi hoàn toàn quy luật thủy văn tự nhiên. Nếu không nạp tệp thông số vận hành hồ chứa (.res) vào mô hình SWAT, thuật toán dòng chảy tự nhiên SCS sẽ dự báo lệch pha so với chuỗi đo thực tế, dẫn đến $NSI$ bị âm.

3. Có thể dùng dữ liệu mưa vệ tinh thay thế cho trạm đo mặt đất không?

Hoàn toàn có thể. Trong trường hợp mạng lưới đo mưa thưa thớt, có thể sử dụng các nguồn ảnh vệ tinh như CHIRPS, TRMM hoặc GPCP để nội suy và trích xuất lượng mưa ngày theo tọa độ tâm của từng tiểu lưu vực. Tuy nhiên, dữ liệu này cần được hiệu chỉnh sai số (Bias Correction) với các trạm thực tế trước khi đưa vào tệp .pcp.

4. Tại sao kết quả mô phỏng chất lượng nước lại có sai số lớn hơn dòng chảy?

Mô phỏng dòng chảy dựa trên các định luật vật lý bảo toàn khối lượng và hình học lưu vực rõ ràng. Ngược lại, chất lượng nước phụ thuộc vào các quá trình hóa sinh phức tạp (chuyển hóa nitơ, phân hủy hữu cơ) và đặc biệt là tải lượng xả thải thực tế từ các nguồn điểm (nhà máy, trại chăn nuôi, khu dân cư) thường không có số liệu kiểm toán đầy đủ, dẫn đến độ chính xác của CLN thấp hơn mô phỏng thủy lực.

5. Chi phí triển khai hệ thống giám sát kết hợp mô hình SWAT là bao nhiêu?

Phần mềm SWAT là mã nguồn mở hoàn toàn miễn phí. Chi phí triển khai chủ yếu tập trung vào: (1) Thu thập và số hóa dữ liệu bản đồ nền ($10 - 30$ triệu VNĐ); (2) Thu thập chuỗi số liệu khí tượng thủy văn; (3) Chi phí chuyên gia tư vấn thiết lập và hiệu chỉnh mô hình. Tổng kinh phí tiết kiệm hơn 80% so với việc mua bản quyền các phần mềm thương mại như MIKE 11 hay SOBEK.

Kết luận

Đồ á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à" đã chứng minh tính khả thi và hiệu quả vượt trội của việc tích hợp công nghệ GIS và mô hình toán bán phân bố trong quản lý tài nguyên nước. Các đóng góp cốt lõi bao gồm:

  1. Thiết lập thành công bộ cơ sở dữ liệu số hoàn chỉnh cho lưu vực sông La Ngà ($4.010\text{ km}^2$), bao gồm địa hình DEM 30m, thổ nhưỡng và sử dụng đất.
  2. Mô phỏng chính xác lưu lượng dòng chảy tự nhiên giai đoạn 1997–1998 với hệ số $R^2$ đạt tới 0,944 và $NSI$ đạt 0,724, khẳng định mô hình SWAT hoàn toàn đáp ứng tốt cho các lưu vực sông nhiệt đới gió mùa tại Việt Nam.
  3. Nhận diện chính xác sự biến động chế độ thủy văn dưới tác động của công trình thủy điện Hàm Thuận – Đa Mi và đánh giá hiện trạng chất lượng nước đạt chuẩn QCVN 08:2008/BTNMT cho hầu hết các mục đích sử dụng (ngoại trừ hàm lượng TSS mùa lũ).

Kết quả nghiên cứu là tài liệu tham khảo khoa học giá trị, cung cấp công cụ trợ giúp đắc lực cho các nhà quản lý môi trường và quy hoạch thủy lợi vùng Đông Nam Bộ trong việc hướng tới quản lý lưu vực sông thông minh và bền vững.