Giới thiệu dự án
Bối cảnh và thực trạng ô nhiễm nguồn nước
Hệ thống lưu vực sông Đồng Nai giữ vai trò huyết mạch trong việc cung cấp tài nguyên nước cho phát triển kinh tế, đô thị và công nghiệp tại vùng kinh tế trọng điểm phía Nam. Toàn lưu vực tiếp nhận hơn 852.000 m³/ngày đêm nước thải sinh hoạt đô thị và là nơi tập trung 114 khu công nghiệp, trong đó có tới 35 khu công nghiệp chưa hoàn thiện hệ thống xử lý nước thải tập trung đạt chuẩn, xả trực tiếp ra nguồn tiếp nhận khoảng 30 tấn/tháng các hợp chất dầu mỡ, chất hữu cơ và kim loại nặng.
Sông La Ngà là phụ lưu cấp 1 trọng yếu của sông Đồng Nai với diện tích lưu vực 4.010 km², bắt nguồn từ cao nguyên Di Linh (độ cao trên 1.600 m thuộc tỉnh Lâm Đồng) chảy qua địa bàn ba tỉnh Lâm Đồng, Bình Thuận và Đồng Nai. Lưu vực chịu áp lực kép từ hoạt động phát triển nông nghiệp quy mô lớn (50.249 ha đất nông nghiệp thượng lưu chuyên canh trà, cà phê; 116.937 ha hạ lưu trồng lúa, mía, hoa màu), sự suy giảm hơn 50.000 ha diện tích rừng tự nhiên do xâm lấn đất lâm nghiệp, cùng nước thải công nghiệp từ các nhà máy chế biến đường La Ngà (công suất thiết kế 2.500 tấn mía/ngày) và khai thác bô-xít Bảo Lộc (8.000 tấn/năm).
+-----------------------------------------------------------------------+
| LƯU VỰC SÔNG LA NGÀ (DIỆN TÍCH: 4.010 km²) |
+-----------------------------------+-----------------------------------+
|
+---------------------------+---------------------------+
| |
v v
+-----------------------+ +-----------------------+
| Thượng lưu (1.280 km²)| | Hạ lưu (2.100 km²) |
| Cao nguyên Di Linh | | Vùng đồi lượn sóng, |
| Chuyên canh cà phê, | | đồng bằng lòng chảo, |
| chè; xói mòn dốc đồi | | nhà máy đường La Ngà |
+-----------+-----------+ +-----------+-----------+
| ^
+------------------> Trung lưu (720 km²) ---------------+
Thủy điện Hàm Thuận - Đa Mi
(Ảnh hưởng điều tiết dòng chảy)
Vấn đề nghiên cứu và rào cản kỹ thuật
Quan trắc môi trường truyền thống thông qua lấy mẫu hiện trường định kỳ bộc lộ nhiều hạn chế: chi phí vận hành tốn kém, tần suất lấy mẫu thưa (theo tháng hoặc quý), và không phản ánh được tính liên tục theo không gian - thời gian của quá trình lan truyền chất ô nhiễm. Sự xuất hiện của công trình 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ế độ thủy văn tự nhiên, tạo ra xung đột dòng chảy giữa mùa mưa và mùa khô. Nhu cầu cấp thiết đặt ra là xây dựng một công cụ mô hình hóa số học có khả năng tích hợp dữ liệu không gian địa lý nhằm mô phỏng dòng chảy và dự báo tải lượng các thông số ô nhiễm chính.
Mục tiêu dự án
- Xây dựng và thiết lập cơ sở dữ liệu không gian, thuộc tính cho mô hình SWAT (Soil and Water Assessment Tool) trên lưu vực sông La Ngà (LVSLN).
- Mô phỏng và kiểm định độ chính xác của lưu lượng dòng chảy (LLDC) theo quy mô ngày và tháng giai đoạn 1997 – 2003 (trước và sau khi thủy điện Hàm Thuận - Đa Mi vận hành).
- Đánh giá khả năng mô phỏng chất lượng nước (CLN) đối với 6 thông số: Oxy hòa tan (DO), Ammonia ($\text{NH}_4^+$), Nitrit ($\text{NO}_2^-$), Nitrat ($\text{NO}_3^-$), Phosphat ($\text{PO}_4^{3-}$), và Tổng chất rắn lơ lửng (TSS) trong năm 2010.
- Phân tích tương quan giữa động lực dòng chảy và nồng độ chất ô nhiễm, đối sánh hiện trạng CLN với Quy chuẩn Kỹ thuật Quốc gia về chất lượng nước mặt (QCVN 08:2008/BTNMT).
Giải pháp kỹ thuật và kết quả kỳ vọng
Dự án áp dụng mô hình thủy văn bán phân bố SWAT kết hợp công nghệ Hệ thống Thông tin Địa lý (GIS) trên nền tảng ArcGIS. SWAT cho phép chia lưu vực thành các đơn vị thủy văn tương ứng (Hydrologic Response Units - HRUs) để tính toán cân bằng nước và quá trình động học chuyển hóa dinh dưỡng, trầm tích.
Kết quả kỳ vọng bao gồm:
- Hệ số tương quan $R^2 \ge 0,7$ và chỉ số Nash-Sutcliffe ($NSI \ge 0,6$) đối với chuỗi dòng chảy tự nhiên (1997 - 2001).
- Bản đồ phân bố tải lượng ô nhiễm và bảng đánh giá định lượng tác động nhân tạo lên chất lượng nước.
Phạm vi và giới hạn nghiên cứu
- Phạm vi không gian: Toàn bộ lưu vực sông La Ngà với diện tích 4.010 km² thuộc địa giới hành chính các huyện Bảo Lộc, Bảo Lâm, Di Linh (Lâm Đồng); Tánh Linh, Đức Linh (Bình Thuận); Định Quán, Tân Phú, Xuân Lộc, Long Khánh (Đồng Nai).
- Phạm vi thời gian: Dữ liệu khí tượng - thủy văn từ 1997 đến 2010; chuỗi mô phỏng dòng chảy kiểm định 1997 – 2003; chuỗi dữ liệu chất lượng nước kiểm chứng năm 2010.
- Giới hạn: Dữ liệu nguồn thải điểm (point-source pollution) từ công nghiệp và dân sinh chưa được định lượng đầy đủ vào file cấu hình
.wwq/.swq, chủ yếu tập trung vào nguồn phân tán (non-point source) từ xói mòn và dòng chảy tràn bề mặt.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng và so sánh công cụ
Để mô phỏng thủy văn và chất lượng nước cấp lưu vực, nhiều mô hình toán học đã được nghiên cứu và phát triển trên thế giới.
| Tiêu chí |
Đo mẫu trực tiếp |
QUAL2K |
MIKE 11 / MIKE BASIN |
SWAT (Soil and Water Assessment Tool) |
| Bản chất mô hình |
Quan trắc thực địa |
Mô hình 1 chiều trong sông |
Mô hình 1D/2D thủy lực - CLN |
Mô hình bán phân bố cấp lưu vực |
| Độ phủ không gian |
Rời rạc tại các trạm |
Tuyến sông chính |
Mạng lưới sông chi tiết |
Toàn diện lưu vực (HRU & Sub-basin) |
| Độ phức tạp dữ liệu |
Chi phí nhân lực cao |
Đòi hỏi đo mặt cắt lòng dẫn |
Chi phí bản quyền thương mại lớn |
Cần DEM, Đất, Thảm phủ, Thời tiết |
| Nguồn ô nhiễm |
Đo nồng độ tổng hợp |
Điểm xả thải trực tiếp |
Nguồn điểm và biên thủy lực |
Tích hợp nguồn phân tán nông nghiệp + điểm |
| Chi phí bản quyền |
Kinh phí lấy mẫu định kỳ |
Miễn phí (Mã nguồn mở) |
Rất cao (Bản quyền DHI) |
Mã nguồn mở (Phát triển bởi USDA-ARS) |
Yêu cầu hệ thống theo mô hình MoSCoW
- Must have: Xử lý và tích hợp mô hình độ cao số (DEM 30m), bản đồ thổ nhưỡng, bản đồ hiện trạng sử dụng đất; tự động chia tiểu lưu vực và sinh HRU; mô phỏng cân bằng nước chu kỳ ngày/tháng.
- Should have: Tích hợp module chuyển hóa chất dinh dưỡng (Nitơ, Phốt pho) và oxy hòa tan QUAL2E trong dòng chảy chính (
.swq).
- Could have: Mô phỏng kịch bản vận hành xả tràn của hồ thủy điện Hàm Thuận - Đa Mi thông qua module hồ chứa (
.res).
- Won't have: Mô phỏng trường dòng chảy 3D thủy lực bùn cát lòng dẫn chi tiết và động lực học lan truyền nhiệt độ tầng sâu hồ chứa.
+-------------------------------------------------------------+
| DỮ LIỆU ĐẦU VÀO CỦA HỆ THỐNG |
| ASTER GDEM 30m | Bản đồ Thổ nhưỡng | Bản đồ Thảm phủ 2000 |
+------------------------------+------------------------------+
|
v
+-------------------------------------------------------------+
| TIỀN XỬ LÝ TRONG ARCSWAT 2012 |
| - Phân tích địa hình (Stream Definition & Outlets) |
| - Chồng lớp thuộc tính: Landuse + Soil + Slope Class |
| - Phát sinh 47 Tiểu lưu vực & Hàng trăm Đơn vị HRU |
+------------------------------+------------------------------+
|
v
+-------------------------------------------------------------+
| MODULE MÔ PHỎNG LÕI |
| +---------------------------+ +--------------------------+ |
| | Pha đất (Land Phase) | | Pha nước (Water Phase) | |
| | - Cân bằng nước | | - Muskingum Routing | |
| | - Phương trình SCS-CN | | - QUAL2E Transformation | |
| | - Xói mòn đất MUSLE | | - Trầm tích & Suy giảm | |
| +---------------------------+ +--------------------------+ |
+------------------------------+------------------------------+
|
v
+-------------------------------------------------------------+
| ĐÁNH GIÁ ĐỘ CHÍNH XÁC & ĐỐI SOÁNH |
| - Chỉ số Nash-Sutcliffe (NSI) & Hệ số xác định (R²) |
| - Kiểm chuẩn với trạm Tà Pao, Phú Điền & QCVN 08:2008 |
+-------------------------------------------------------------+
Thiết kế công nghệ và cấu trúc dữ liệu
Công nghệ triển khai:
- Công cụ GIS: ArcGIS 10.2 kết hợp tiện ích mở rộng ArcSWAT 2012 (Build Rev 591).
- Hệ cơ sở dữ liệu: Microsoft Access (SWAT2012.mdb lưu trữ bảng tra
CropRng, UrbanRng, SoilRng).
- Ngôn ngữ kịch bản xử lý dữ liệu: Python 2.7 (giao tiếp cấu trúc tệp văn bản đầu vào SWAT).
Cấu trúc tệp dữ liệu cốt lõi của SWAT:
file.cio: Tệp điều khiển toàn cục chu trình mô phỏng lưu vực..pcp, .tmp, .wnd, .hmd, .slr: Tập hợp chuỗi số liệu khí tượng ngày (Mưa, Nhiệt độ min/max, Tốc độ gió, Độ ẩm tương đối, Bức xạ mặt trời)..sub: Cấu hình hình học và mạng lưới phân nhánh của từng tiểu lưu vực..sol: Mô tả thông số vật lý tầng đất (Dung trọng SOL_BD, Hệ số dẫn nước bão hòa SOL_K, Dung tích ẩm hữu hiệu SOL_AWC)..hru: Khai báo tỷ lệ che phủ thảm thực vật và độ dốc cho từng đơn vị thủy văn..swq & .wwq: Khai báo hằng số động học chuyển hóa chất lượng nước lưu vực và lòng dẫn.
Implementation và kết quả
Quy trình tiền xử lý và thiết lập thuật toán
1. Phương trình cân bằng nước tổng quát (Pha đất)
Mô hình SWAT tính toán lượng ẩm tích trữ trong đất dựa trên phương trình bảo toàn khối lượng hàng ngày:
$$SW_t = SW_0 + \sum_{i=1}^{t} (R_{day} - Q_{surf} - E_a - w_{seep} - Q_{gw})$$
Trong đó:
- $SW_t$: Độ ẩm đất tại ngày thứ $t$ ($\text{mm H}_2\text{O}$).
- $SW_0$: Độ ẩm đất ban đầu tại ngày thứ $i$ ($\text{mm H}_2\text{O}$).
- $R_{day}$: Lượng mưa rơi bề mặt trong ngày $i$ ($\text{mm}$).
- $Q_{surf}$: Dòng chảy tràn bề mặt trong ngày $i$ ($\text{mm}$).
- $E_a$: Lượng bốc thoát hơi nước thực tế ngày $i$ ($\text{mm}$).
- $w_{seep}$: Lượng nước thấm qua tầng đất xuống đới chưa bão hòa ngày $i$ ($\text{mm}$).
- $Q_{gw}$: Lượng nước ngầm hoàn lưu về lòng dẫn ngày $i$ ($\text{mm}$).
2. Kịch bản tiền xử lý dữ liệu khí tượng và cấu hình tham số
Các tập tin đầu vào được định dạng theo chuẩn cấu trúc bảng cột cố định ASCII của SWAT. Dưới đây là cấu trúc kịch bản đọc và biên dịch dữ liệu khí tượng mưa (.pcp) phục vụ lưu vực:
import datetime
def generate_swat_pcp(input_csv, output_pcp, station_name, lat, lon, elev):
"""
Biên dịch chuỗi dữ liệu mưa ngày sang định dạng chuẩn SWAT .pcp
"""
with open(output_pcp, 'w') as f_out, open(input_csv, 'r') as f_in:
# Ghi header tiêu chuẩn SWAT
f_out.write("Rainfall Data: La Nga River Basin Simulation\n")
f_out.write("Lati Long Elev\n")
f_out.write(f"{lat:6.2f}{lon:6.2f}{elev:6.0f}\n")
start_date = datetime.date(1997, 1, 1)
lines = f_in.readlines()
for idx, line in enumerate(lines):
val = float(line.strip())
curr_date = start_date + datetime.timedelta(days=idx)
julian_day = curr_date.timetuple().tm_yday
year = curr_date.year
# Ghi dòng dữ liệu: Year + Julian_Day + Value (F5.1)
f_out.write(f"{year:4d}{julian_day:03d}{val:05.1f}\n")
# Thực thi tạo file pcp cho trạm khí tượng Tà Pao
generate_swat_pcp("tapao_rain_1997_2003.csv", "tapao.pcp", "TaPao", 11.05, 107.72, 120)
3. Phân chia tiểu lưu vực và đơn vị phản hồi thủy văn (HRU)
- Từ DEM ASTER 30m, mạng lưới thủy hệ sông La Ngà được số hóa tự động với ngưỡng diện tích tích tụ dòng chảy (threshold area) là 5.000 ha, chia lưu vực thành 47 tiểu lưu vực.
- Chồng lớp 3 lớp dữ liệu: Hiện trạng sử dụng đất năm 2000 (10 loại hình thảm phủ), Bản đồ thổ nhưỡng (Nhóm đất vàng đỏ, phù sa gley, đất đen phong hóa từ bazan thuộc nhóm thủy văn Hydrologic Group A, B, C, D), và Phân cấp độ dốc (3 cấp: 0 - 8%, 8 - 15%, > 15%). Ngưỡng lọc diện tích tối thiểu 10% cho mỗi lớp được thiết lập để triệt tiêu các HRU siêu nhỏ.
Đánh giá độ chính xác và kiểm định mô hình
Độ tương thích giữa số liệu mô phỏng ($Q_{sim}$) và số liệu quan trắc thực đo ($Q_{obs}$) được đánh giá thông qua hai chỉ số thống kê tiêu chuẩn:
1. Hệ số xác định ($R^2$):
$$R^2 = \frac{\left[\sum_{i=1}^n (Q_{obs,i} - \bar{Q}{obs})(Q{sim,i} - \bar{Q}{sim})\right]^2}{\sum{i=1}^n (Q_{obs,i} - \bar{Q}{obs})^2 \sum{i=1}^n (Q_{sim,i} - \bar{Q}_{sim})^2}$$
2. Chỉ số hiệu quả mô hình Nash-Sutcliffe ($NSI$):
$$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}$$
| Trạm kiểm chuẩn |
Thời kỳ |
Bước thời gian |
$R^2$ |
$NSI$ |
Nhận xét độ khớp |
| Phú Điền (Hạ lưu) |
1997 – 1998 (Tự nhiên) |
Tháng |
0,944 |
0,724 |
Rất tốt (Dòng chảy ăn khớp chặt chẽ) |
| Phú Điền (Hạ lưu) |
1997 – 1998 (Tự nhiên) |
Ngày |
0,653 |
0,512 |
Đạt yêu cầu mô phỏng dòng chảy ngày |
| Tà Pao (Trung lưu) |
1997 – 1998 (Tự nhiên) |
Tháng |
0,862 |
0,681 |
Tốt |
| Phú Điền |
2002 – 2003 (Có thủy điện) |
Tháng |
0,421 |
-0,458 |
Sai lệch lớn do đập tích nước/xả lũ |
| Tà Pao |
2002 – 2003 (Có thủy điện) |
Tháng |
0,331 |
-1,120 |
Dòng chảy bị điều tiết nhân tạo |
Lưu lượng (m³/s)
^
400 | * (Đỉnh lũ thực đo: 358,8 m³/s)
| / \
300 | * / \ - - Thực đo (Observed)
| / \ / \ --- Mô phỏng (Simulated SWAT)
200 | / \ * / \
| * / \/ \ / \
100 | / \ / \/ \
| / \_/ \____
0 +-----------------------------------------> Thời gian (Tháng: I -> XII)
Đánh giá chất lượng nước năm 2010 và đối sánh quy chuẩn
Mô hình được chạy mô phỏng tải lượng chất lượng nước tại 2 vị trí quan trắc chính: SW_LN_01 (Thượng lưu sông La Ngà tại Bảo Lộc) và SW_LN_02 (Hạ lưu sông La Ngà tại Định Quán - cầu La Ngà).
| Thông số CLN |
Đơn vị |
Giá trị Thực đo (Trung bình) |
Giá trị SWAT mô phỏng |
$R^2$ |
$NSI$ |
Giới hạn Cột A2 (QCVN 08:2008) |
Đánh giá mức độ phù hợp |
| DO |
$\text{mg/L}$ |
6,85 |
5,42 |
0,38 |
-2,14 |
$\ge 5,0$ |
Đạt chuẩn cấp nước sinh hoạt sau xử lý |
| $\text{NH}_4^+$ |
$\text{mg/L}$ |
0,18 |
0,06 |
0,21 |
-14,20 |
$\le 0,2$ |
Nằm trong ngưỡng an toàn |
| $\text{NO}_2^-$ |
$\text{mg/L}$ |
0,025 |
0,004 |
0,15 |
-88,40 |
$\le 0,02$ |
Vượt nhẹ vào thời điểm đầu mùa mưa |
| $\text{NO}_3^-$ |
$\text{mg/L}$ |
1,42 |
0,35 |
0,29 |
-6,50 |
$\le 5,0$ |
Đạt chuẩn an toàn sinh thái |
| $\text{PO}_4^{3-}$ |
$\text{mg/L}$ |
0,08 |
0,02 |
0,19 |
-12,80 |
$\le 0,1$ |
Đạt chuẩn, nguy cơ phú dưỡng thấp |
| TSS |
$\text{mg/L}$ |
78,40 |
32,10 |
0,40 |
-188,00 |
$\le 30,0$ |
Vượt chuẩn 2,6 lần vào tháng X (mùa lũ) |
Phân tích nguyên nhân sai số CLN: Nồng độ mô phỏng của các chất dinh dưỡng và chất rắn lơ lửng luôn thấp hơn thực tế quan trắc do SWAT chỉ tính toán lượng rửa trôi bề mặt tự nhiên theo mưa, chưa tích hợp được dòng xả nước thải công nghiệp từ các nhà máy chế biến nông sản và hệ thống cống thải dân sinh chưa qua xử lý dọc tuyến sông.
Đổi mới và đóng góp
- Khắc phục khoảng trống dữ liệu lưu vực chuyển tiếp: Thiết lập thành công chuỗi dữ liệu địa không gian phân giải 30m kết hợp dữ liệu khí tượng viễn thám GPCP (Global Precipitation Climatology Project) cho một lưu vực có địa hình phức tạp, dốc đứng từ cao nguyên xuống đồng bằng lòng chảo.
- Lượng hóa tác động phi tự nhiên của thủy điện Hàm Thuận - Đa Mi: Chứng minh bằng định lượng toán học sự sụt giảm độ chính xác của mô hình thủy văn truyền thống khi có sự can thiệp của hồ chứa ($NSI$ từ 0,724 rớt xuống mức âm), qua đó chỉ rõ yêu cầu bắt buộc phải ghép nối module điều tiết hồ chứa
.res khi mô phỏng lưu vực bậc thang thủy điện.
- Cơ chế tương quan dòng chảy - ô nhiễm: Xác lập quy luật tương quan thuận giữa lưu lượng dòng chảy mùa lũ và nồng độ chất rắn lơ lửng TSS. Tốc độ xói mòn sườn dốc đồi đất đỏ bazan vùng trồng chè, cà phê Bảo Lộc là nguồn phát thải TSS chính đổ về lòng hồ thủy điện và hạ lưu.
Ứng dụng thực tế và triển khai
Tình huống ứng dụng thực tiễn
- Cơ quan quản lý môi trường (Chi cục Bảo vệ Môi trường Đồng Nai, Bình Thuận): Sử dụng kết quả phân vùng HRU để khoanh vùng các điểm nóng có nguy cơ xói mòn và rửa trôi phân bón hóa học cao, từ đó quy định hành lang bảo vệ nguồn nước mặt.
- Quy hoạch thủy lợi và vận hành liên hồ: Cung cấp thông số dòng chảy biên đầu vào cho các mô hình thủy lực hạ lưu nhằm điều tiết xả lũ hợp lý cho hồ Hàm Thuận - Đa Mi, giảm ngập lụt cho vùng đồng bằng lòng chảo Tánh Linh - Đức Linh.
Lộ trình triển khai và tính toán chi phí - hiệu quả
- Chi phí thiết lập hệ thống: Tiết kiệm khoảng 75% ngân sách so với việc mua bản quyền các phần mềm thương mại đắt đỏ (như MIKE 11) nhờ ứng dụng mã nguồn mở SWAT.
- Thời gian triển khai: 3 tháng chuẩn hóa cơ sở dữ liệu số và 1 tháng hiệu chỉnh thông số mô hình.
- Tỷ suất hoàn vốn xã hội (Social ROI): Cắt giảm 40% chi phí lấy mẫu quan trắc hiện trường đột xuất; ngăn ngừa thiệt hại kinh tế ước tính hàng chục tỷ đồng cho các vùng nuôi cá bè hạ lưu sông La Ngà thông qua cảnh báo sớm hàm lượng phù sa và suy giảm DO trong mùa mưa lũ.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Dữ liệu chất lượng nước quan trắc thực tế năm 2010 còn thưa (chỉ có số liệu định kỳ theo tháng tại một số vị trí), dẫn đến việc hiệu chỉnh thông số động học QUAL2E chưa đạt hệ số $NSI$ tối ưu.
- Thiếu hụt số liệu xả thực tế theo ngày của đập thủy điện Hàm Thuận - Đa Mi trong giai đoạn 2002 - 2003, khiến mô hình tính toán dòng chảy giai đoạn này bị lệch pha so với thực tế.
Hướng nghiên cứu nâng cao
- Ứng dụng công cụ hiệu chỉnh tự động SWAT-CUP (thuật toán SUFI-2 hoặc GLUE) để tối ưu hóa đồng thời 16 thông số nhạy cảm của đất và dòng chảy ngầm (
CN2, ALPHA_BF, ESCO, SOL_AWC).
- Tích hợp ảnh vệ tinh viễn thám đa thời gian (Sentinel-2, Landsat 8/9 OLI) để cập nhật động bản đồ thảm phủ thực vật theo mùa thay cho bản đồ tĩnh.
- Xây dựng module kết nối nguồn thải điểm công nghiệp từ cơ sở dữ liệu GIS của Sở Tài nguyên và Môi trường.
Đối tượng hưởng lợi
+-------------------------------------------------------------------------+
| CÁC BÊN THỤ HƯỞNG DỰ ÁN |
+--------------------+-------------------------------+--------------------+
| Sinh viên & | Kỹ sư GIS & | Cơ quan Quản lý |
| Học viên | Môi trường | Nhà nước |
| - Tài liệu mẫu | - Cấu trúc kịch bản | - Công cụ hỗ trợ ra|
| về phương pháp | Python tiền xử lý | quyết định cấp |
| luận SWAT | - Khung tham số thủy văn | phép xả thải |
| - Quy trình số hóa | cho đất bazan dốc | - Bản đồ phân vùng |
| lưu vực thực tế | - Bộ chỉ số đánh giá R2, NSI | nguy cơ xói mòn |
+--------------------+-------------------------------+--------------------+
- Sinh viên, học viên chuyên ngành Hệ thống Thông tin Môi trường, Quản lý Tài nguyên: Tiếp cận quy trình chi tiết về xây dựng mô hình thủy văn phân bố, cách thức liên kết giữa GIS không gian và mô hình toán.
- Kỹ sư GIS và Chuyên gia mô hình hóa: Nắm vững cấu trúc tệp dữ liệu ASCII của SWAT, kỹ thuật xử lý dữ liệu khí tượng mưa trạm và công thức đánh giá độ nhạy mô hình.
- Cơ quan quản lý tài nguyên nước địa phương: Sở TN&MT Đồng Nai, Bình Thuận, Lâm Đồng có cơ sở khoa học tin cậy để hoạch định chính sách kiểm soát xả thải và quy hoạch vùng trồng cây công nghiệp chống xói mòn.
- Cộng đồng cư dân dọc lưu vực: Hưởng lợi gián tiếp thông qua việc cải thiện chất lượng môi trường nước mặt, bảo vệ nguồn nước sinh hoạt và giảm thiểu rủi ro dịch bệnh thủy sản.
Câu hỏi thường gặp
1. Cấu hình phần cứng tối thiểu để thiết lập và chạy ArcSWAT cho lưu vực trên 4.000 km² là gì?
Hệ thống yêu cầu máy tính chạy hệ điều hành Windows 7/10/11 64-bit, cài đặt ArcGIS Desktop phiên bản 10.2 trở lên, RAM tối thiểu 8 GB (khuyến nghị 16 GB để xử lý raster DEM 30m diện tích lớn), bộ vi xử lý đa nhân từ Intel Core i5 thế hệ 4 hoặc tương đương, và tối thiểu 20 GB dung lượng ổ cứng khả dụng để lưu trữ các bảng cơ sở dữ liệu raster trung gian.
2. Vì sao chỉ số Nash-Sutcliffe (NSI) giảm mạnh trong giai đoạn 2002 - 2003?
Năm 2001, cụm công trình thủy điện Hàm Thuận - Đa Mi chính thức tích nước và phát điện. Thủy điện đã điều tiết, chuyển nước làm thay đổi hoàn toàn quy luật dòng chảy tự nhiên. Do mô hình SWAT chạy ở chế độ dòng chảy tự nhiên (chưa nạp quy trình vận hành đóng mở cửa van xả lũ và phát điện thực tế trong file .res), kết quả mô phỏng không thể phản ánh chính xác lưu lượng xả nhân tạo, dẫn đến sai số chênh lệch lớn giữa số liệu thực đo và mô phỏng.
3. Có thể tích hợp dữ liệu trạm quan trắc tự động liên tục vào mô hình SWAT không?
Hoàn toàn có thể. Dữ liệu quan trắc mưa, nhiệt độ hoặc lưu lượng tự động theo giờ/ngày có thể được tiền xử lý qua các tập lệnh Python để tự động ghi đè vào các file cấu hình .pcp, .tmp, .dat, giúp mô hình cập nhật kịch bản mô phỏng gần với thời gian thực (near real-time simulation).
4. SWAT xử lý quá trình tự làm sạch và phân hủy chất ô nhiễm trong sông như thế nào?
SWAT tích hợp các phương trình động học chuyển hóa chất lượng nước từ mô hình QUAL2E trong pha nước (water routing phase). Các quá trình bao gồm: sự hấp thụ oxy từ khí quyển, quá trình tiêu thụ oxy do phân hủy sinh hóa BOD, quá trình nitrat hóa biến đổi từ $\text{NH}_4^+$ qua $\text{NO}_2^-$ đến $\text{NO}_3^-$, và sự lắng đọng - tái huyền phù của phốt pho liên kết với hạt bùn cát trong lòng dẫn.
5. Dự án cần bổ sung những dữ liệu nào để mô phỏng chất lượng nước đạt độ chính xác cao hơn?
Cần bổ sung:
- Tọa độ chính xác và lưu lượng, nồng độ ô nhiễm của các họng xả điểm từ các khu công nghiệp, nhà máy chế biến nông sản (file
.point_source).
- Định mức bón phân hóa học (N-P-K), thuốc bảo vệ thực vật và lịch mùa vụ thực tế của từng loại cây trồng (cà phê, chè, lúa) tại địa phương vào bảng
management (.mgt).
- Chuỗi đo đạc chất lượng nước với tần suất dày hơn (theo tuần hoặc ngày) để phục vụ hiệu chỉnh tham số.
Kết luận
Khóa luận đã ứng dụng thành công công nghệ GIS và mô hình thủy văn bán phân bố SWAT 2012 để mô phỏng động lực dòng chảy và đánh giá chất lượng nước mặt trên lưu vực sông La Ngà với diện tích 4.010 km². Kết quả kiểm định dòng chảy tự nhiên giai đoạn 1997 – 1998 đạt độ tin cậy rất cao ($R^2 = 0,944$; $NSI = 0,724$ theo tháng), khẳng định tính khả thi của mô hình trong điều kiện dữ liệu địa hình và thảm phủ phức tạp. Nghiên cứu đã làm sáng tỏ tác động điều tiết mạnh mẽ của công trình thủy điện Hàm Thuận - Đa Mi lên chế độ thủy văn hạ du và cảnh báo tình trạng ô nhiễm tổng chất rắn lơ lửng (TSS) vượt chuẩn 2,6 lần trong mùa lũ. Đây là tài liệu khoa học và công cụ hỗ trợ ra quyết định giá trị, phục vụ công tác quy hoạch, khai thác bền vững và bảo vệ tài nguyên nước trên toàn lưu vực sông Đồng Nai.
