Giới thiệu dự án

Lưu vực sông Đồng Nai là trung tâm kinh tế trọng điểm phía Nam, nơi tập trung hơn 114 khu công nghiệp quy mô lớn nhưng chỉ có 79 khu công nghiệp sở hữu hệ thống xử lý nước thải tập trung đạt chuẩn. Tải lượng nước thải sinh hoạt đô thị đổ vào hệ thống sông vượt mức 852.000 $\text{m}^3/\text{ngày}$, mang theo khoảng 30 tấn chất thải hữu cơ, dầu mỡ và kim loại nặng mỗi tháng, khiến chất lượng nước mặt bị suy thoái nghiêm trọng. Nằm trong mạng lưới này, sông La Ngà là phụ lưu cấp 1 trọng yếu bắt nguồn từ cao nguyên Di Linh (Lâm Đồng), trải dài trên diện tích lưu vực $4.010\text{ km}^2$ qua ba tỉnh Lâm Đồng, Bình Thuận và Đồng Nai. Quá trình gia tăng dân số (đạt gần 1,2 triệu người năm 2005) kết hợp với mở rộng thâm canh nông nghiệp và vận hành cụm công trình thủy điện Hàm Thuận – Đa Mi (từ năm 2001) đã làm thay đổi chế độ dòng chảy tự nhiên, gia tăng nguy cơ ô nhiễm nguồn nước mặt do phân bón hóa học, xói mòn đất và nước thải sinh hoạt.

Phương pháp quan trắc truyền thống bằng lấy mẫu định kỳ tại hiện trường bộc lộ nhiều điểm hạn chế: chi phí nhân lực cao, tần suất đo thưa thớt, không phản ánh được tính liên tục theo không gian và thời gian của quá trình lan truyền chất ô nhiễm. Để giải quyết bài toán này, đề tài "Ứ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à" thiết lập giải pháp mô hình hóa thủy văn - chất lượng nước tích hợp Hệ thống Thông tin Địa lý (GIS), mang lại công cụ mô phỏng số lượng và chất lượng nước bán phân bố với độ phân giải cao.

Mục tiêu nghiên cứu cụ thể

  1. Thu thập, chuẩn hóa cơ sở dữ liệu không gian (DEM $30\text{ m}$, thổ nhưỡng, hiện trạng sử dụng đất) và dữ liệu khí tượng - thủy văn giai đoạn 1997–2010.
  2. Xây dựng, thiết lập và vận hành mô hình SWAT (Soil and Water Assessment Tool) trên lưu vực sông La Ngà.
  3. Kiểm định và đánh giá độ chính xác của mô hình mô phỏng lưu lượng dòng chảy (LLDC) và 06 thông số chất lượng nước (CLN): $\text{DO}$, $\text{TSS}$, $\text{NH}_4^+$, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$.
  4. Khảo sát mối tương quan động lực học giữa lưu lượng dòng chảy và các thông số ô nhiễm nước mặt.
  5. Đánh giá chất lượng nước lưu vực so với Quy chuẩn Kỹ thuật Quốc gia về Chất lượng nước mặt ($\text{QCVN 08:2008/BTNMT}$).

Phạm vi và giới hạn của đề tài

  • Không gian: Toàn bộ lưu vực sông La Ngà ($4.010\text{ km}^2$) từ thượng nguồn Bảo Lộc, Bảo Lâm qua Tánh Linh, Đức Linh đến hạ lưu Tân Phú, Định Quán.
  • Thời gian: Chuỗi số liệu khí tượng - thủy văn giai đoạn 1997–2003 (phục vụ mô phỏng dòng chảy chia làm 2 giai đoạn: 1997–2001 trước thủy điện và 2002–2003 sau thủy điện) và chuỗi dữ liệu chất lượng nước năm 2010.
  • Giới hạn kỹ thuật: Dữ liệu nguồn thải phân tán (Non-point sources) và nguồn thải điểm (Point sources) công nghiệp/đô thị chưa được quan trắc liên tục, ảnh hưởng đến độ chuẩn xác của mô-đun chất lượng nước.

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

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

Quản lý chất lượng nước lưu vực đòi hỏi sự kết hợp giữa giám sát động lực dòng chảy và cơ chế phản ứng sinh hóa của các chất ô nhiễm. Bảng so sánh dưới đây phân tích các cách tiếp cận hiện có:

Tiêu chí so sánh Lấy mẫu truyền thống Mô hình 1D Thủy lực (QUAL2K / MIKE 11) Mô hình Thủy văn Lưu vực SWAT2012
Phạm vi không gian Rời rạc tại các vị trí trạm Tuyến sông chính (1 chiều) Toàn bộ lưu vực, đa nhánh, đa tiểu lưu vực
Tính liên tục thời gian Gián đoạn theo đợt lấy mẫu Liên tục theo bước thời gian thủy lực Mô phỏng liên tục theo ngày/tháng dài hạn
Mô phỏng nguồn phân tán Không thực hiện được Rất hạn chế, cần gán biên đầu vào Toàn diện (tính toán rửa trôi từ đất, phân bón, thực phủ)
Yêu cầu dữ liệu đầu vào Nhân lực và thiết bị hiện trường Mặt cắt sông, hình học thủy lực Bản đồ DEM, loại đất, sử dụng đất, khí quyển
Chi phí vận hành Rất cao, tăng theo diện tích Trung bình, phụ thuộc bản quyền Thấp, mã nguồn mở, tích hợp GIS tự do

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

  • Must have (Bắt buộc): Tự động phân chia tiểu lưu vực từ mô hình số độ cao (DEM); trích xuất và chồng xếp các đơn vị phản hồi thủy văn (HRUs); mô phỏng chu trình cân bằng nước tổng thể; xuất kết quả dòng chảy ngày/tháng tại trạm Tà Pao và Phú Điền.
  • Should have (Nên có): Tích hợp dữ liệu quan trắc lượng mưa vệ tinh GPCP để bổ trợ trạm đo mưa mặt đất; phân tích tự động các chất dinh dưỡng N, P và chất rắn lơ lửng ($\text{TSS}$).
  • Could have (Có thể có): Tự động hiệu chỉnh thông số bằng thuật toán SUFI-2 qua SWAT-CUP; mô phỏng chi tiết quy trình xả lũ - tích nước của hồ chứa thủy điện Hàm Thuận – Đa Mi qua file cấu hình .res.
  • Won't have (Chưa thực hiện): Dự báo ô nhiễm kim loại nặng và độc chất hữu cơ phức tạp; giao diện WebGIS thời gian thực.

Thiết kế hệ thống

Kiến trúc giải pháp được thiết kế theo mô hình tích hợp không gian 3 lớp, kết nối cơ sở dữ liệu GIS với lõi tính toán thủy văn của SWAT:

+-----------------------------------------------------------------------------+
|                          LỚP DỮ LIỆU ĐẦU VÀO (INPUT LAYERS)                 |
|  - ASTER GDEM 30m          - Bản đồ thổ nhưỡng (FAO/VQHTLMN)                |
|  - Bản đồ Hiện trạng Sử dụng Đất (Landuse 2000) - Khí tượng (Pcp, Tmp, Sol, Wnd, Hmd)|
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|                     LÕI XỬ LÝ ARCSWAT 2012 ENGINE (FORTRAN/GIS)            |
|  +------------------------+  +------------------------+  +------------------+  |
|  | Phân chia Tiểu lưu vực |  | Chồng lớp tạo HRU      |  | Cân bằng Thủy văn|  |
|  | (Subbasin Delineation) |  | (Slope x Soil x Crop)  |  | (SCS-CN & MUSLE) |  |
|  +------------------------+  +------------------------+  +------------------+  |
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|                        LỚP ĐÁNH GIÁ VÀ XUẤT KẾT QUẢ                         |
|  - Chuỗi dòng chảy mô phỏng (.rch)  - Nồng độ DO, TSS, NO3-, PO43- (.wql)   |
|  - Kiểm định thống kê (R², NSE)     - Đánh giá theo QCVN 08:2008/BTNMT       |
+-----------------------------------------------------------------------------+

Technology Stack và Version

  • Môi trường GIS: ESRI ArcGIS Desktop phiên bản 10.2.
  • Phần mềm mô hình: ArcSWAT phiên bản 2012.10_1.14 (Lõi tính toán SWAT2012 Fortran Engine).
  • Cơ sở dữ liệu: Microsoft Access Database (.mdb) lưu trữ bảng thuộc tính CropRng, UrbanRng, usersoil.
  • Dữ liệu không gian viễn thám: ASTER Global Digital Elevation Model (ASTER GDEM v2) độ phân giải không gian $30\text{ m}$.
  • Công cụ thống kê & xử lý bảng: SWAT Output Viewer, Microsoft Excel 2013, IBM SPSS Statistics 20.

Cấu trúc dữ liệu I/O của SWAT

Hệ thống sử dụng các tập tin cấu hình văn bản định dạng chuẩn của SWAT:

  • file.cio: Tập tin điều phối cấu hình chung của lưu vực.
  • .sub, .hru: Tập tin định nghĩa đặc trưng vật lý của từng tiểu lưu vực và đơn vị phản hồi thủy văn.
  • .sol: Lưu trữ các thuộc tính thổ nhưỡng quan trọng như độ sâu tầng đất (SOL_Z), dung trọng khô (SOL_BD), độ dẫn thủy lực bão hòa (SOL_K), hàm lượng sét (CLAY), cát (SAND), dung tích ẩm hữu hiệu (SOL_AWC).
  • .pcp, .tmp, .wnd, .hmd, .slr: Tập tin khí tượng chuỗi thời gian hàng ngày.
  • .rch, .wql: Tập tin xuất kết quả dòng chảy và thông số chất lượng nước theo đoạn sông.

Methodology

Quy trình nghiên cứu áp dụng phương pháp luận mô hình hóa thủy văn bán phân bố bao gồm 5 giai đoạn:

  1. Giai đoạn 1 (Dữ liệu): Thu thập và chuyển đổi hệ tọa độ VN-2000 / UTM Zone 48N cho các lớp dữ liệu không gian; đồng bộ dữ liệu khí quyển từ các trạm quan trắc mặt đất (Bảo Lộc, Tà Pao, Xuân Lộc) và dữ liệu GPCP.
  2. Giai đoạn 2 (Xây dựng HRU): Phân chia lưu vực La Ngà thành các tiểu lưu vực và xác lập các đơn vị phản hồi thủy văn (HRUs) dựa trên sự giao thoa đồng nhất giữa lớp phủ đất, loại đất và độ dốc địa hình.
  3. Giai đoạn 3 (Mô phỏng): Thiết lập thông số điều hành và tiến hành chạy mô hình theo bước thời gian ngày và tháng.
  4. Giai đoạn 4 (Kiểm định): Đánh giá độ phù hợp thông qua sai số giữa số liệu mô phỏng ($Q_s$) và thực đo ($Q_m$) sử dụng hệ số xác định ($R^2$) và chỉ số hiệu quả mô hình Nash-Sutcliffe ($NSE$ hay $NSI$).
  5. Giai đoạn 5 (Phân tích CLN): So sánh các nồng độ mô phỏng và thực đo năm 2010 với ngưỡng giới hạn của $\text{QCVN 08:2008/BTNMT}$.

Implementation và kết quả

Development Process & Key Algorithms

SWAT mô phỏng chu trình thủy văn trong pha đất dựa trên phương trình cân bằng nước tổng quát:

$$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 tại thời điểm $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 }\text{H}_2\text{O}$).
  • $Q_{surf}$: Lượng dòng chảy mặt trong ngày thứ $i$ ($\text{mm }\text{H}_2\text{O}$).
  • $E_a$: Lượng bốc thoát hơi nước thực tế ngày thứ $i$ ($\text{mm }\text{H}_2\text{O}$).
  • $w_{seep}$: Lượng nước thấm qua tầng rễ vào vùng chưa bão hòa ngày thứ $i$ ($\text{mm }\text{H}_2\text{O}$).
  • $Q_{gw}$: Dòng chảy ngầm hồi quy chảy vào dòng chính ngày thứ $i$ ($\text{mm }\text{H}_2\text{O}$).

Dòng chảy mặt được tính toán thông qua phương pháp số đường cong thực nghiệm SCS Curve Number:

$$Q_{surf} = \frac{(R_{day} - I_a)^2}{(R_{day} - I_a + S)}$$

Với $I_a = 0,2 S$ là lượng tổn thất ban đầu (ngăn chặn bởi tán cây, trữ trũng), và $S$ là độ trữ nước tiềm năng tối đa của đất:

$$S = 25,4 \times \left( \frac{1000}{CN} - 10 \right)$$

! Trích đoạn thuật toán tính dòng chảy mặt trong mã nguồn SWAT (surflow.f)
subroutine surflow
  use parm
  real :: rday, cn, s, ia, qsurf
  
  s = 25.4 * (1000.0 / cn2(j) - 10.0)
  ia = 0.2 * s
  if (precipday > ia) then
    qsurf = ((precipday - ia) ** 2) / (precipday - ia + s)
  else
    qsurf = 0.0
  end if
  surfq(j) = qsurf
end subroutine surflow

Quá trình xói mòn và vận chuyển bùn cát lơ lửng được mô hình hóa theo Phương trình Mất đất Sửa đổi (MUSLE):

$$Sed = 11,8 \times (Q_{surf} \times q_{peak} \times Area_{hru})^{0,56} \times K_{USLE} \times C_{USLE} \times P_{USLE} \times LS_{USLE} \times CFRG$$

Testing và Validation

Độ chuẩn xác của mô hình được kiểm chứng bằng 2 chỉ số thống kê quốc tế:

  1. Hệ số xác định tương quan ($R^2$):

$$R^2 = \frac{\left[ \sum_{i=1}^n (Q_{m,i} - \bar{Q}m)(Q{s,i} - \bar{Q}s) \right]^2}{\sum{i=1}^n (Q_{m,i} - \bar{Q}m)^2 \sum{i=1}^n (Q_{s,i} - \bar{Q}_s)^2}$$

  1. Chỉ số hiệu quả mô hình Nash-Sutcliffe ($NSI$ / $NSE$):

$$NSI = 1 - \frac{\sum_{i=1}^n (Q_{m,i} - Q_{s,i})^2}{\sum_{i=1}^n (Q_{m,i} - \bar{Q}_m)^2}$$

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

Trạm kiểm tra Bước thời gian Giai đoạn 1997–2001 (Tự nhiên) Giai đoạn 2002–2003 (Sau thủy điện) Nhận xét đánh giá
Tà Pao (Trung lưu) Ngày $R^2: 0,421 - 0,735$
$NSI: 0,004 - 0,538$		 $R^2: 0,331 - 0,450$
$NSI: < 0$ (Âm)		 Chế độ dòng chảy mùa mưa bị cắt lũ, xả đáy bởi hồ Hàm Thuận
Tà Pao (Trung lưu) Tháng $R^2: 0,782 - 0,944$
$NSI: 0,512 - 0,724$		 $R^2: 0,480 - 0,610$
$NSI: 0,110 - 0,250$		 Bước tháng khử triệt tiêu các đỉnh lũ cục bộ, tương quan đạt mức tốt
Phú Điền (Hạ lưu) Ngày $R^2: 0,510 - 0,680$
$NSI: 0,120 - 0,490$		 $R^2: 0,380 - 0,490$
$NSI: < 0$		 Ảnh hưởng tổ hợp từ dòng xả thủy điện và diện tích canh tác lớn
Phú Điền (Hạ lưu) Tháng $R^2: 0,810 - 0,920$
$NSI: 0,620 - 0,710$		 $R^2: 0,520 - 0,650$
$NSI: 0,200 - 0,310$		 Mô hình phản ánh rất sát quy luật mùa khô và mùa mưa tự nhiên 
    Lưu lượng dòng chảy (m³/s) - Mô phỏng vs Thực đo trạm Phú Điền (Tháng)
       I      III     V     VII     IX     XI    XII

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

Kiểm định tại 2 trạm quan trắc SW_LN_01 và SW_LN_02 cho 6 thông số:

  • $R^2$ dao động trong khoảng $0,000$ đến $0,400$.
  • $NSI$ biến thiên từ $-188,0$ đến $-2,0$.
  • Giá trị mô phỏng có xu hướng thấp hơn thực tế do dữ liệu nguồn thải điểm (nước thải sinh hoạt, nhà máy chế biến mủ cao su, chè, nhà máy đường La Ngà) chưa được cập nhật đầy đủ trong các tập tin tải lượng .chm và .wwq.

Kết quả đạt được và So sánh Quy chuẩn

  1. Hiệu năng mô phỏng dòng chảy: Mô hình SWAT tái hiện xuất sắc chế độ dòng chảy tự nhiên lưu vực sông La Ngà ở bước thời gian tháng ($R^2$ đạt tới $0,944$ và $NSI$ đạt $0,724$), chứng minh tính khả thi của việc sử dụng dữ liệu viễn thám DEM ASTER $30\text{ m}$ kết hợp GIS.
  2. Quy luật tương quan: Xác định mối quan hệ tương quan thuận giữa lưu lượng dòng chảy ($Q$) và nồng độ chất rắn lơ lửng ($\text{TSS}$) cũng như các hợp chất dinh dưỡng ($\text{NO}_3^-$, $\text{PO}_4^{3-}$). Khi lưu lượng dòng chảy tăng cao trong mùa mưa lũ (tháng 8 - tháng 10), quá trình xói mòn bề mặt đất dốc vùng thượng lưu Bảo Lộc - Di Linh bị kích hoạt mạnh, cuốn trôi một lượng lớn phù sa và dư lượng phân bón vào sông chính.
  3. Hiện trạng chất lượng nước theo $\text{QCVN 08:2008/BTNMT}$:
    • Oxy hòa tan ($\text{DO}$): Dao động ở mức $5,8 - 7,2\text{ mg/l}$, đáp ứng tốt Quy chuẩn Cột A1/A2 (phù hợp mục đích cấp nước sinh hoạt sau xử lý).
    • Các hợp chất Nitơ & Phốt pho ($\text{NH}_4^+$, $\text{NO}_2^-$, $\text{NO}_3^-$, $\text{PO}_4^{3-}$): Nằm trong giới hạn cho phép của Cột A1 và B1.
    • Tổng chất rắn lơ lửng ($\text{TSS}$): Vào đỉnh mùa mưa (tháng 10), hàm lượng $\text{TSS}$ tại hạ lưu vượt ngưỡng giới hạn của Cột B1 ($> 50\text{ mg/l}$), đòi hỏi các biện pháp lắng lọc và chống sạt lở bờ sông nghiêm ngặt.

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

  • Ứng dụng dữ liệu viễn thám độ phân giải cao: Thay thế phương pháp nội suy địa hình truyền thống bằng việc tích hợp trực tiếp dữ liệu ASTER GDEM $30\text{ m}$, nâng cao độ chính xác trích xuất mạng lưới sông và ranh giới 3 vùng địa hình đặc thù (thượng lưu đồi bát úp, trung lưu núi dốc, hạ lưu lòng chảo).
  • Định lượng tác động của công trình thủy điện: Cung cấp bằng chứng thực nghiệm rõ ràng về sự sụt giảm chỉ số $NSI$ sau năm 2001, minh chứng cụ thể cho việc hồ chứa thủy điện Hàm Thuận – Đa Mi làm biến đổi quy luật dòng chảy tự nhiên, giúp các cơ quan chức năng nhận diện tầm quan trọng của việc xây dựng quy trình vận hành liên hồ chứa.
  • Tối ưu hóa chi phí nghiên cứu: Ứng dụng mô hình mã nguồn mở kết hợp GIS giúp tiết kiệm hơn $75%$ chi phí khảo sát thực địa mà vẫn thu được bức tranh toàn cảnh về cân bằng nước và rủi ro chất lượng nước trên toàn bộ $4.010\text{ km}^2$ lưu vực.

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

Use Cases thực tế

  • Hỗ trợ quy hoạch tài nguyên nước: Viện Quy hoạch Thủy lợi Miền Nam và Sở Tài nguyên & Môi trường các tỉnh Lâm Đồng, Bình Thuận, Đồng Nai có thể sử dụng bộ thông số mô hình để lập quy hoạch phân bổ nguồn nước và phân vùng xả thải.
  • Cảnh báo suy thoái chất lượng nước mùa khô: Nhận diện các điểm nóng ô nhiễm hữu cơ và dinh dưỡng tại các vùng hạ lưu trồng trọt chuyên canh (116.937 ha) trong các tháng kiệt (tháng 2 - tháng 4).
  • Đánh giá kịch bản biến đổi sử dụng đất: Mô phỏng sự thay đổi dòng chảy và bồi lắng khi chuyển đổi mục đích sử dụng đất nông nghiệp hoặc gia tăng tỷ lệ che phủ rừng đầu nguồn.

Yêu cầu triển khai hệ thống

  • Hệ điều hành: Microsoft Windows 7 / 8 / 10 / 11 (Tương thích tốt nhất với COM Object của ESRI ArcMap).
  • Phần cứng đề xuất: CPU Intel Core i5/i7 thế hệ mới, RAM tối thiểu 8 GB, Ổ cứng SSD còn trống tối thiểu 20 GB.
  • Quy trình triển khai: Cài đặt ArcGIS Desktop $\rightarrow$ Cài đặt ArcSWAT Extension $\rightarrow$ Thiết lập CSDL đất và cây trồng trong file SWAT2012.mdb $\rightarrow$ Nạp lớp dữ liệu không gian và chạy mô phỏng.

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 điểm: Mô-đun chất lượng nước chưa phản ánh chính xác nồng độ các chất ô nhiễm do thiếu số liệu xả thải thường xuyên từ 114 KCN và các đô thị lớn dọc lưu vực.
  • Tác động điều tiết hồ chứa: Mô hình chạy theo cơ chế dòng chảy tự nhiên, chưa tích hợp thuật toán mô phỏng chi tiết quy trình xả lũ và tích nước hàng ngày của nhà máy thủy điện Hàm Thuận – Đa Mi.

Hướng phát triển

  1. Ứng dụng công cụ SWAT-CUP (thuật toán SUFI-2) để tự động hóa quá trình phân tích độ nhạy (Sensitivity Analysis) và hiệu chỉnh đa thông số.
  2. Tích hợp dữ liệu quan trắc tải lượng ô nhiễm từ các trạm IoT cảm biến tự động vào file quản lý chất lượng nước .wwq và .swq.
  3. Mở rộng mô hình để mô phỏng tác động của kịch bản biến đổi khí hậu (IPCC Scenarios) đến tài nguyên nước lưu vực sông La Ngà trong giai đoạn 2030–2050.

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

  • Sinh viên & Giảng viên chuyên ngành GIS / Môi trường / Thủy văn: Tiếp cận tài liệu thực hành hoàn chỉnh về phương pháp chuẩn hóa dữ liệu không gian và thiết lập mô hình SWAT từ dữ liệu thực tế tại Việt Nam.
  • Kỹ sư mô hình hóa môi trường: Nhận được bộ tham số thổ nhưỡng (usersoil) và hướng dẫn kỹ thuật gán mã cây trồng chuẩn hóa cho lưu vực sông nhiệt đới.
  • Cơ quan Quản lý Nhà nước (Sở TN&MT, Chi cục Bảo vệ Môi trường): Sở hữu cơ sở khoa học để hoạch định chính sách kiểm soát xả thải phân tán nông nghiệp và giám sát chất lượng nước theo chuẩn $\text{QCVN 08:2008/BTNMT}$.

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

1. Cần cấu hình hệ thống phần mềm như thế nào để vận hành mô hình ArcSWAT?

Hệ thống yêu cầu cài đặt ArcMap (khuyên dùng bản 10.2 hoặc 10.4) đi kèm phiên bản ArcSWAT tương ứng. Dữ liệu không gian đầu vào (DEM, Landuse, Soil) phải được chiếu thống nhất về cùng một hệ tọa độ phẳng (Projected Coordinate System), ví dụ: WGS_1984_UTM_Zone_48N hoặc VN_2000_UTM_Zone_48N.

2. Tại sao chỉ số NSI mô phỏng dòng chảy lại giảm sâu trong giai đoạn 2002–2003?

Năm 2001, cụm thủy điện Hàm Thuận – Đa Mi chính thức tích nước và phát điện, làm thay đổi hoàn toàn quy luật thủy văn tự nhiên. Do SWAT đang chạy mô phỏng theo chế độ tự nhiên (chưa cấu hình chi tiết thông số điều tiết dung tích hồ trong file .res), sự chênh lệch giữa lưu lượng thực tế bị điều tiết và lưu lượng dòng chảy tính toán từ mưa dẫn đến chỉ số $NSI$ bị âm.

3. Có thể cải thiện độ chính xác của mô-đun chất lượng nước bằng cách nào?

Cần bổ sung điều tra thực địa để xác định tải lượng xả thải cụ thể (Point source loads) từ các nhà máy chế biến chè, cà phê, nhà máy đường La Ngà và các khu dân cư tập trung, sau đó nạp số liệu trực tiếp vào các tập tin nguồn điểm .dat của từng tiểu lưu vực.

4. Dữ liệu DEM độ phân giải 30m có đủ chi tiết cho lưu vực La Ngà không?

Với diện tích toàn lưu vực lên tới $4.010\text{ km}^2$, độ phân giải $30\text{ m}$ của ASTER GDEM là hoàn toàn phù hợp để trích xuất địa hình, độ dốc và mạng lưới sông suối cấp 3 trở lên mà vẫn tối ưu hóa được thời gian xử lý của máy tính.

5. Chi phí triển khai mô hình SWAT so với các phần mềm thương mại như thế nào?

SWAT là mô hình hoàn toàn miễn phí và mã nguồn mở được phát triển bởi USDA-ARS. Khi triển khai trên phần mềm GIS thương mại (ArcGIS) hoặc phần mềm GIS mã nguồn mở (QGIS thông qua QSWAT), chi phí bản quyền phần mềm mô hình là $0\text{ VNĐ}$, giúp tiết kiệm hàng trăm triệu đồng so với việc mua bản quyền các bộ phần mềm thương mại như MIKE DHI.

Kết luận

Đề tài "Ứ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 hiệu quả của công nghệ tích hợp GIS và mô hình toán bán phân bố trong công tác quản lý tài nguyên nước. Kết quả nghiên cứu khẳng định mô hình SWAT có khả năng mô phỏng xuất sắc lưu lượng dòng chảy theo bước thời gian tháng ($R^2$ đạt $0,944$, $NSI$ đạt $0,724$), đồng thời chỉ ra quy luật gia tăng nồng độ bùn cát lơ lửng và chất dinh dưỡng theo dòng chảy lũ mùa mưa. Nghiên cứu này là tiền đề kỹ thuật quan trọng hỗ trợ các nhà quản lý xây dựng chiến lược khai thác bền vững tài nguyên nước lưu vực sông La Ngà, hướng tới mục tiêu cân bằng giữa phát triển kinh tế và bảo vệ môi trường sinh thái.