Giới thiệu dự án

Lưu vực sông Vàm Cỏ Đông là trục thủy đạo chiến lược thuộc hệ thống sông Đồng Nai, giữ vai trò quyết định đối với an ninh nguồn nước, sinh hoạt và phát triển kinh tế của hai tỉnh Tây Ninh và Long An. Với tổng diện tích lưu vực khoảng 8.500 km² (diện tích lưu vực kín gần 6.000 km²), chiều dài đoạn chảy qua Tây Ninh là 151 km và Long An là 145 km, dòng sông cung cấp nguồn nước tưới tiêu cho hơn 70% diện tích đất nông nghiệp của vùng và tiếp nhận nước thải của hơn 1,27 triệu cư dân cùng hàng loạt khu - cụm công nghiệp tập trung.

Tuy nhiên, áp lực kép từ biến đổi khí hậu, xâm nhập mặn theo chế độ bán nhật triều không đều và tốc độ công nghiệp hóa - đô thị hóa nhanh chóng đang đẩy nguồn nước mặt sông Vàm Cỏ Đông vào tình trạng suy thoái chất lượng nghiêm trọng. Việc kiểm soát ô nhiễm gặp nhiều thách thức lớn:

  • Tải lượng nước thải sinh hoạt đô thị chưa qua xử lý đạt xấp xỉ $127.468 \text{ m}^3/\text{ngày}$ (tính theo định mức $100 \text{ lít/người/ngày}$, hệ số thải 0,8 của WHO).
  • Nước thải chăn nuôi gia súc tại Tây Ninh phát sinh $15.000 \text{ m}^3/\text{ngày}$, tải lượng ô nhiễm hữu cơ lên tới $84 \text{ tấn BOD}_5/\text{ngày}$, $253 \text{ tấn TSS/ngày}$, $19,6 \text{ tấn Tổng N/ngày}$ và $5,8 \text{ tấn Tổng P/ngày}$.
  • Nước thải công nghiệp từ 18 khu công nghiệp (KCN) và 42 cụm công nghiệp (CCN) theo quy hoạch Long An (lưu lượng xả thải đơn lẻ của các nhà máy mía đường, dệt nhuộm lên đến $7.000 \text{ m}^3/\text{ngày}$) và các KCN trọng điểm tại Tây Ninh như Trảng Bàng ($12.800 \text{ m}^3/\text{ngày}$), Bourbon An Hòa ($19.200 \text{ m}^3/\text{ngày}$), Phước Đông ($19.200 \text{ m}^3/\text{ngày}$).
  • Hiện tượng xả phèn nội đồng từ hơn 50% diện tích lưu vực là đất phèn hoạt động trong hai giai đoạn đầu mùa mưa (tháng 4 - 7) và cuối mùa mưa (tháng 11 - 1), kéo độ pH của nước sông xuống dưới ngưỡng sinh thái.
  • Dữ liệu quan trắc phân tán, lưu trữ dạng bảng biểu tĩnh, thiếu công cụ trực quan hóa không gian đa thời gian phục vụ ra quyết định điều hành liên tỉnh.
flowchart TD
    A[Áp lực ô nhiễm lưu vực Sông Vàm Cỏ Đông] --> B[Nguồn thải sinh hoạt: 127.468 m³/ngày]
    A --> C[Nguồn thải công nghiệp: >50.000 m³/ngày]
    A --> D[Nguồn thải nông nghiệp: 15.000 m³/ngày, 84 tấn BOD₅/ngày]
    A --> E[Hiện tượng tự nhiên: Rửa phèn & Xâm nhập mặn 12,6‰]
    
    B & C & D & E --> F[Hệ thống tích hợp: WQI + GIS IDW + Lohani Matrix]
    F --> G[Bản đồ hóa không gian chất lượng nước ArcGis]
    F --> H[Bảng xếp hạng ưu tiên can thiệp môi trường]
    G & H --> I[Quyết định quản lý & Giải pháp công nghệ FAST/Biogas]

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

  1. Thu thập, phân tích và đánh giá biến động của 10 thông số lý hóa - vi sinh quan trắc chất lượng nước mặt tại 10 trạm trọng điểm dọc toàn tuyến sông Vàm Cỏ Đông trong hai năm 2013 và 2014 theo hai mùa tương phản (mùa khô tháng 3 và mùa mưa tháng 9).
  2. Chuẩn hóa và áp dụng thuật toán tính toán Chỉ số Chất lượng Nước (Water Quality Index - WQI) kết hợp kỹ thuật nội suy không gian nghịch đảo trọng số khoảng cách (Inverse Distance Weighting - IDW) trên nền tảng ArcGIS để xây dựng hệ thống bản đồ số chuyên đề về diễn biến ô nhiễm.
  3. Ứng dụng ma trận phân tích đa tiêu chí Lohani để định lượng, xếp hạng mức độ ưu tiên của các áp lực môi trường trên lưu vực.
  4. Đề xuất đồng bộ gói giải pháp quản lý lưu vực liên tỉnh và giải pháp công nghệ kỹ thuật (xử lý nước thải phi tập trung FAST, Biogas công nghiệp, máy trục vớt lục bình cơ giới).

Phạm vi và giới hạn

  • Phạm vi không gian: Trải dài 296 km dọc dòng chính sông Vàm Cỏ Đông, từ trạm VCD1 (vùng biên giới Việt Nam - Campuchia tại xã Biên Giới, huyện Châu Thành, Tây Ninh) đến trạm VCD10 (khu vực ngã ba cửa sông Vàm Cỏ Đông hợp lưu với sông Vàm Cỏ Tây tại Tân Trụ, Long An).
  • Phạm vi thời gian: Chuỗi số liệu phân tích chuyên sâu mùa khô (tháng 3) và mùa mưa (tháng 9) trong 2 năm liên tiếp (2013 - 2014), đối chiếu quy chuẩn kỹ thuật quốc gia QCVN 08:2008/BTNMT (cột A2, B1, B2).

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

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

Trước khi triển khai phương pháp tích hợp WQI - GIS, công tác quản trị dữ liệu chất lượng nước tại các địa phương đối mặt với nhiều bất cập lớn:

Phương pháp đánh giá Ưu điểm Nhược điểm Khả năng ứng dụng thực tế
So sánh đơn thông số với TCVN/QCVN Đơn giản, tra cứu trực tiếp từng chỉ tiêu vi phạm. Rời rạc, không phản ánh được bức tranh tổng thể nguồn nước; khó truyền thông cho nhà quản lý và cộng đồng. Phù hợp thanh tra, xử phạt từng nguồn thải cụ thể.
Mô phỏng thủy lực 1D/2D (HEC-RAS, MIKE 11) Dự báo chi tiết vận tốc dòng chảy, quá trình lan truyền chất ô nhiễm. Đòi hỏi dữ liệu địa hình mặt cắt đáy sông, khí tượng và biên thủy văn cực kỳ chi tiết; chi phí tính toán cao. Phù hợp nghiên cứu quy hoạch dài hạn cấp vùng.
Tích hợp WQI + GIS Interpolation (Đề tài ứng dụng) Thu gọn đa thông số về chỉ số vô đinh lượng (1 - 100 hoặc thang 10 cấp); hiển thị trực quan không gian thực tế qua màu sắc. Phụ thuộc vào mật độ trạm quan trắc và hàm phân bố trọng số. Tối ưu nhất cho giám sát hiện trạng, cảnh báo sớm và ra quyết định quản lý môi trường.

Ma trận MoSCoW xác định yêu cầu hệ thống

  • Must have: Thuật toán tính chỉ số WQI theo bộ trọng số chuẩn Tôn Thất Lãng (2006); Bản đồ nội suy IDW đẳng trị WQI và các thông số đơn lẻ ($pH, DO, BOD_5, TSS, N, P, Fe, Coliform$); Ma trận xếp hạng Lohani.
  • Should have: Bảng đối chiếu vi phạm vượt ngưỡng QCVN 08:2008/BTNMT cột A2/B1/B2; Báo cáo phân tích tương quan thủy văn triều bán nhật.
  • Could have: Module hồi quy tuyến tính dự báo xu thế chất lượng nước giai đoạn 2015 - 2020; Mô phỏng phân tán ô nhiễm cục bộ bằng ENVIMWQ.
  • Won't have (ở giai đoạn này): Cảm biến truyền dẫn dữ liệu thời gian thực IoT trực tiếp từ hiện trường.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu môi trường không gian bao gồm 4 tầng chức năng:

[Tầng 1: Thu thập & Phân tích mẫu chuẩn] 
   -> Quan trắc hiện trường (APHA, TCVN) tại 10 trạm (VCD1 - VCD10)
   -> 10 thông số: pH, Độ mặn, TSS, DO, BOD5, Tổng N, Tổng P, Tổng Fe, Coliform, E.Coli
        |
[Tầng 2: Chuẩn hóa & Tính toán chỉ số chuyên sâu]
   -> Biến đổi thông số sang chỉ số phụ Pi (Thang điểm 1 - 10)
   -> Áp dụng bộ trọng số Ti: BOD5 (0.23), DO (0.18), TSS (0.16), pH (0.15), Tổng N (0.15), Coliform (0.13)
   -> Công thức: WQI = Σ (Pi * Ti)
        |
[Tầng 3: Mô hình hóa Không gian GIS (ArcGIS 10.3 / ArcPy)]
   -> Dữ liệu nền địa lý: MicroStation V8i (*.dgn), MapInfo Pro (*.tab) sang ESRI Shapefile (*.shp)
   -> Raster Spatial Interpolation: IDW (Power p = 2, Search Radius Variable)
   -> Phân vùng chất lượng nước 6 cấp màu theo chuẩn quản lý
        |
[Tầng 4: Phân tích quyết định ưu tiên & Đề xuất giải pháp]
   -> Ma trận Lohani: Ui = Pi * Ri * Ci
   -> Định hình gói can thiệp: Ban quản trị lưu vực + FAST + Biogas + Máy vớt cơ giới

Technology Stack và Đặc tả kỹ thuật

  • GIS Platform: ESRI ArcGIS Desktop 10.3 (ArcMap, ArcCatalog, Spatial Analyst Extension)
  • Spatial Processing Core: Python 2.7 / ArcPy API, NumPy, SciPy
  • CAD/GIS Data Conversion: Bentley MicroStation V8i, Pitney Bowes MapInfo Professional 12.5
  • Chuẩn phân tích phòng thí nghiệm:
    • $pH$: TCVN 6492:2000 (Đo điện thế)
    • Độ mặn: APHA 2520, 20th Ed. (Độ dẫn điện)
    • $TSS$: TCVN 4560:1988 (Khối lượng sấy khô $105^\circ\text{C}$)
    • $DO$: APHA 4500-O G, 20th Ed. (Điện cực màng quang học / Winkler)
    • $BOD_5$: TCVN 6001:1995 (Ủ 5 ngày ở $20^\circ\text{C}$)
    • Tổng N: TCVN 6624-1:2000 (Vô cơ hóa Kjeldahl)
    • Tổng P: TCVN 6202:1996 (Trắc phổ Amoni Molypdat)
    • Tổng Fe: TCVN 6177:1996 (Trắc quang Phenanthroline)
    • Vi sinh: APHA 9221B-9222C (Coliform MPN), SMEWW 9221F ($E. coli$)

Cơ sở dữ liệu thuộc tính và không gian (GeoDatabase Schema)

Bảng dữ liệu quan trắc được thiết kế chuẩn hóa quan hệ (PostGIS / Geodatabase Feature Class):

  • station_id (VARCHAR(10), Primary Key): Định danh trạm (VCD1 - VCD10)
  • geom (GEOMETRY(Point, 4326)): Tọa độ trạm (WGS 84)
  • sample_date (DATE): Thời gian lấy mẫu (03/2013, 09/2013, 03/2014, 09/2014)
  • season (VARCHAR(10)): Mùa quan trắc (Kho/Mua)
  • ph_val, do_val, bod5_val, tss_val, tot_n_val, coliform_val: Dữ liệu đo thực nghiệm
  • sub_ph, sub_do, sub_bod5, sub_tss, sub_tot_n, sub_coliform: Chỉ số phụ chuẩn hóa $P_i$ ($1.0 - 10.0$)
  • wqi_score (NUMERIC(4,2)): Điểm WQI tổng hợp
  • water_quality_class (VARCHAR(50)): Phân loại cấp nước (Xanh dương, Lục, Vàng, Cam, Đỏ, Nâu)

Implementation và kết quả

Development Process & Key Algorithms

1. Thuật toán tính toán Chỉ số Chất lượng Nước (WQI)

Chỉ số WQI được tính toán thông qua cấu trúc tính trung bình cộng có trọng số của các chỉ số phụ:

$$WQI = \sum_{i=1}^{n} (P_i \times T_i)$$

Trong đó:

  • $P_i$: Chỉ số phụ của thông số chất lượng nước thứ $i$ (thang điểm 1 đến 10 dựa trên hàm chuyển đổi nội suy tuyến tính từng đoạn theo nồng độ chất ô nhiễm).
  • $T_i$: Trọng số tương ứng của thông số thứ $i$ theo quy chuẩn phân bổ:
Thông số Trọng số ($T_i$) Rationale môi trường
$BOD_5$ 0,23 Đại diện chính cho ô nhiễm hữu cơ và tiêu hao oxy sinh hóa
$DO$ 0,18 Điều kiện tiên quyết duy trì hệ sinh thái thủy sinh
$TSS$ 0,16 Đánh giá hàm lượng phù sa, bùn cát lơ lửng và độ đục
$pH$ 0,15 Đánh giá độ chua phèn và cân bằng ion hydro
Tổng N 0,15 Đánh giá mức độ phú dưỡng nguồn nước
Coliform 0,13 Chỉ thị vi sinh và nguy cơ lây truyền mầm bệnh đường ruột

2. Kỹ thuật nội suy không gian nghịch đảo khoảng cách (IDW)

Giá trị chất lượng nước tại điểm chưa lấy mẫu $Z(s_0)$ trên mặt lưới raster sông được ước lượng thông qua công thức:

$$Z(s_0) = \frac{\sum_{i=1}^{N} w_i(s_0) Z(s_i)}{\sum_{i=1}^{N} w_i(s_0)} \quad \text{với} \quad w_i(s_0) = \frac{1}{d(s_0, s_i)^p}$$

Trong đó $d(s_0, s_i)$ là khoảng cách Euclidean giữa điểm dự báo $s_0$ và điểm trạm quan trắc $s_i$. Bậc lũy thừa khoảng cách được xác lập tối ưu là $p = 2$, đảm bảo mức độ suy giảm ảnh hưởng của trạm ở xa diễn ra mượt mà dọc theo hình thái uốn khúc của lòng sông.

3. Kịch bản ArcPy xử lý tự động tính toán và xuất bản đồ

Đoạn mã Python tích hợp ArcPy thực hiện tính toán chỉ số phụ, WQI và tự động nội suy Raster:

# -*- coding: utf-8 -*-
"""
Tự động hóa tính toán WQI và nội suy không gian IDW cho sông Vàm Cỏ Đông
Platform: ArcGIS 10.3 / Python 2.7 (ArcPy)
"""
import arcpy
from arcpy.sa import *

# Cấu hình môi trường làm việc
arcpy.env.workspace = "C:/VamCoDong_GIS/Geodatabase.gdb"
arcpy.env.overwriteOutput = True
arcpy.CheckOutExtension("Spatial")

# Trọng số chuẩn Tôn Thất Lãng (2006)
WEIGHTS = {
    'BOD5': 0.23,
    'DO': 0.18,
    'TSS': 0.16,
    'pH': 0.15,
    'TotN': 0.15,
    'Coliform': 0.13
}

def calculate_wqi(sub_indices):
    """Tính điểm WQI tổng hợp từ dict chứa các chỉ số phụ Pi"""
    wqi_score = sum(sub_indices[param] * WEIGHTS[param] for param in WEIGHTS)
    return round(wqi_score, 2)

def execute_idw_interpolation(input_features, z_field, mask_river_poly, output_raster_name):
    """Thực thi nội suy không gian IDW với p=2 cắt theo ranh giới lòng sông"""
    try:
        cell_size = 25.0  # Độ phân giải không gian 25m
        power_p = 2
        search_radius = RadiusVariable(10, 50000) # Lấy tối đa 10 điểm trong bán kính 50km
        
        # Thực hiện IDW
        idw_out = Idw(input_features, z_field, cell_size, power_p, search_radius)
        
        # Mask cắt theo lòng sông Vàm Cỏ Đông
        clipped_idw = ExtractByMask(idw_out, mask_river_poly)
        clipped_idw.save(output_raster_name)
        arcpy.AddMessage(u"Đã tạo thành công Raster: {0}".format(output_raster_name))
    except Exception as e:
        arcpy.AddError(u"Lỗi trong quá trình nội suy: {0}".format(str(e)))

# Thực thi mẫu cho đợt quan trắc Tháng 03/2014
if __name__ == "__main__":
    input_points = "Monitoring_Stations_2014_Dry"
    river_boundary = "VamCoDong_River_Polygon"
    out_raster = "Raster_WQI_Dry_2014"
    execute_idw_interpolation(input_points, "WQI_Score", river_boundary, out_raster)

4. Đánh giá ưu tiên môi trường theo ma trận Lohani

Mức độ ưu tiên của từng vấn đề môi trường ($U_i$) được xác định theo mô hình toán học đa biến:

$$U_i = P_i \times R_i \times C_i$$

Trong đó:

  • $R_i$ (Range index - Chỉ số địa lý): Cấp 1 (Cấp tỉnh), Cấp 2 (Vùng sinh thái liên huyện), Cấp 3 (Toàn bộ lưu vực).
  • $P_i$ (Persistence index - Chỉ số đối kháng/thời gian giải quyết): 1 (đến năm 2015), 2 (đến năm 2020), 3 (đến năm 2025), 4 (sau năm 2025).
  • $C_i$ (Complexity index - Chỉ số phức hợp): Tổng số thành phần môi trường chịu tác động trong 9 yếu tố sinh thái - nhân văn - tài nguyên (biến thiên từ 1 đến 9).

Testing và Validation

Dữ liệu quan trắc thực tế tại 10 trạm qua 4 đợt khảo sát (Tháng 3/2013, Tháng 9/2013, Tháng 3/2014, Tháng 9/2014) được thẩm định nghiêm ngặt qua các chỉ tiêu lý hóa:

Bảng số liệu phân tích diễn biến các thông số quan trắc then chốt

Mã trạm Vị trí quan trắc pH (Mùa khô / Mưa 2014) DO [mg/L] (Khô / Mưa 2014) TSS [mg/L] (Khô / Mưa 2014) Độ mặn [‰] (Khô / Mưa 2014) Đặc điểm nguồn tiếp nhận chính
VCD1 Xã Biên Giới (Campuchia) 6,6 / 6,7 5,5 / 5,7 24,0 / 42,5 0,0 / 0,1 Kiểm soát biên giới, ít xả thải cục bộ
VCD2 Khu vực Cảng Bến Kéo 6,3 / 6,4 3,8 / 4,2 52,0 / 61,0 0,1 / 0,0 Giao thông thủy, bốc dỡ hàng hóa, dân cư
VCD3 Khu vực Cầu Gò Dầu 6,1 / 6,3 3,5 / 4,0 48,0 / 58,5 0,3 / 0,0 Đô thị Gò Dầu, làng nghề bánh tráng, KCN
VCD4 Thượng lưu NM đường Hiệp Hòa 5,4 / 5,6 1,2 / 2,1 35,0 / 49,0 0,8 / 0,0 Nước thải chế biến đường, mật rỉ, hữu cơ cao
VCD5 Bến đò Xuân Khánh - Hậu Nghĩa 5,3 / 5,5 2,4 / 3,1 42,0 / 68,0 1,2 / 0,0 Bến đò, xói lở bờ, nông nghiệp thoát thủy
VCD6 Kênh Trà Cú (cách ngã 3 sông 1km) 5,0 / 5,1 2,8 / 3,6 19,6 / 31,5 1,8 / 0,0 Đất phèn hoạt động mạnh, nước chua xả nội đồng
VCD7 Kênh Xáng (gần cống xả Đ.Ấn Độ) 5,2 / 5,4 2,1 / 2,9 45,0 / 62,0 2,5 / 0,0 Tiếp nhận nước thải công nghiệp mía đường
VCD8 Kênh Thủ Thừa 5,2 / 5,5 3,1 / 3,8 28,0 / 44,0 3,4 / 0,0 Nước tiêu nội đồng, rửa phèn Đồng Tháp Mười
VCD9 Khu vực Cầu Bến Lức 5,8 / 6,0 3,6 / 4,3 38,0 / 54,0 5,8 / 0,2 Đô thị hóa cao, KCN Thuận Đạo, Nhựt Chánh
VCD10 Cửa sông Vàm Cỏ Đông (Tân Trụ) 6,6 / 6,5 4,2 / 4,8 58,0 / 65,0 12,5 / 0,3 Hỗn hợp bán nhật triều, xâm nhập mặn biển Đông
QCVN QCVN 08:2008/BTNMT (A2) 6,0 - 8,5 $\ge$ 5,0 $\le$ 30 - Ngưỡng an toàn cấp nước sinh hoạt 
BIỂU ĐỒ DIỄN BIẾN OXY HÒA TAN (DO) MÙA KHÔ NĂM 2014 SO VỚI QUY CHUẨN (mg/L)
[QCVN 08 Cột A2: >= 5.0 mg/L] ========================================== (5.0)

6.0 |  * (5.5)
5.0 | 
4.0 |         * (3.8)  * (3.5)                          * (3.6)   * (4.2)
3.0 |                                   * (2.8)  * (3.1)
2.0 |                          * (2.4)          * (2.1)
1.0 |                 * (1.2)
0.0 +-----+--------+--------+--------+--------+--------+--------+--------+--------+--------+
    VCD1  VCD2     VCD3     VCD4     VCD5     VCD6     VCD7     VCD8     VCD9     VCD10

Kết quả đạt được

  1. Phân hóa pH và hiện tượng axit hóa nguồn nước:

    • pH toàn lưu vực dao động từ 5,0 đến 6,7. Đoạn trung lưu từ trạm VCD4 đến VCD8 hoàn toàn không đạt QCVN 08:2008/BTNMT cột A2 ($6,0 - 8,5$).
    • Trạm VCD6 (Kênh Trà Cú) có pH thấp nhất đạt $5,0$ do lượng phèn hoạt động từ đất sulphate chua xả tràn vào mùa mưa và tích tụ đậm đặc vào mùa khô. Nước chỉ kiềm nhẹ trở lại tại VCD1 ($6,6 - 6,7$) và vùng cửa sông VCD10 ($6,5 - 6,6$) nhờ đệm carbonate từ nước biển.

  2. Ô nhiễm hữu cơ và tình trạng suy kiệt Oxy hòa tan ($DO - BOD_5$):

    • DO tại 9/10 trạm (trừ VCD1) đều không đạt chuẩn A2 ($\ge 5 \text{ mg/L}$) vào mùa khô. Điểm đen nghiêm trọng nhất là VCD4 (thượng lưu cống xả nhà máy đường Hiệp Hòa) với DO chạm đáy $1,2 \text{ mg/L}$ (năm 2013) và $1,5 \text{ mg/L}$ (năm 2014), thấp hơn cả tiêu chuẩn B2 ($\ge 2 \text{ mg/L}$ dùng cho giao thông thủy), chứng tỏ vi sinh vật dị dưỡng kỵ khí đã tiêu thụ sạch oxy để phân hủy dòng thải hữu cơ cực lớn.
    • Nồng độ $BOD_5$ và các chất dinh dưỡng (Tổng N, Tổng P) tăng đột biến tại các đoạn sông tiếp nhận nước thải sinh hoạt đô thị và công nghiệp chế biến như Gò Dầu, Bến Lức.

  3. Hiện tượng xói lở và hàm lượng chất rắn lơ lửng ($TSS$):

    • Mùa mưa, 100% trạm quan trắc vượt quy chuẩn A2 ($30 \text{ mg/L}$), đạt đỉnh tại VCD5 ($68 \text{ mg/L}$) do dòng chảy tràn bề mặt gây xói mòn vách bờ sông dốc (tốc độ xâm thực bờ đạt $0,2 - 0,4 \text{ m/năm}$) kết hợp bùn cát từ kênh tiêu nội đồng đổ về.

  4. Động thái xâm nhập mặn:

    • Độ mặn biến thiên cực đoan từ $0,0‰$ ở thượng nguồn lên đến $12,5 - 12,6‰$ tại cửa sông VCD10 trong tháng 3 (mùa khô). Tuy nhiên, dòng nước ngọt bổ sung từ thượng nguồn và lưu lượng xả điều tiết từ hồ Dầu Tiếng ($10 - 12 \text{ m}^3/\text{s}$) đã đẩy mặn lùi về hạ lưu, giúp mùa mưa độ mặn toàn tuyến giảm xuống dưới $0,3‰$.

  5. Tổng hợp xếp hạng chỉ số WQI và không gian hóa:

    • Giá trị WQI tính toán phân bố từ 2,8 (Mức 5: Đỏ - Ô nhiễm nặng tại VCD4 và VCD7) đến 7,9 (Mức 2: Lục - Ô nhiễm rất nhẹ tại VCD1).
    • Bản đồ ArcGis chỉ ra rõ rệt: chất lượng nước năm 2014 có xu hướng cải thiện nhẹ so với 2013 (WQI trung bình tăng từ $4,2$ lên $4,8$) nhờ các biện pháp siết chặt kiểm tra cống xả các nhà máy đường, tinh bột mì của Sở TN&MT Tây Ninh và Long An.
PHÂN VÙNG CHẤT LƯỢNG NƯỚC SÔNG VÀM CỎ ĐÔNG QUA CHỈ SỐ WQI (MÙA KHÔ 2014)
=======================================================================================
Vị trí:   [Thượng nguồn VCD1] ====> [Trung lưu VCD4 - VCD7] ====> [Hạ lưu VCD9 - VCD10]
WQI:             7,8                      3,1 - 3,8                      5,2 - 6,1
Màu sắc:     XANH LỤC                        ĐỎ / CAM                       VÀNG
Chất lượng: Ô nhiễm rất nhẹ               Ô nhiễm nặng                 Ô nhiễm nhẹ
Mục đích:   Cấp sinh hoạt (xử lý)         Giao thông thủy              Tưới tiêu nông nghiệp
=======================================================================================

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

  1. Chuẩn hóa khung phương pháp luận tích hợp WQI - GIS IDW - Lohani: Khắc phục hoàn toàn tình trạng "dữ liệu nằm trên giấy" bằng việc chuyển đổi toàn bộ chuỗi số liệu giám sát môi trường phức tạp thành bản đồ đẳng trị trực quan số hóa. Giúp các cấp lãnh đạo không chuyên về hóa môi trường có thể nhận diện ngay các "điểm nóng" cần can thiệp khẩn cấp.

  2. So sánh cải tiến kỹ thuật so với các công trình trước đây:

Tiêu chí Tiếp cận truyền thống (Sở TN&MT 2007 - 2012) Mô hình số trị chuyên sâu (ENVIMWQ, MIKE) Giải pháp tích hợp của đề tài (WQI + GIS + Lohani)
Tính trực quan Kém (Báo cáo bảng biểu số liệu tĩnh) Cao (Biểu đồ 2D/3D dạng số hóa) Rất cao (Bản đồ phân vùng chuyên đề chuẩn GIS)
Chi phí & Thời gian Thấp nhưng hiệu quả chỉ đạo thấp Cực kỳ tốn kém, cần hiệu chỉnh biên phức tạp Tối ưu, xử lý nhanh chóng trong vài giờ
Độ phủ không gian Cục bộ tại từng vị trí cống xả Giới hạn theo lưới tính toán Bao quát liên tục 296 km toàn tuyến liên tỉnh
Hỗ trợ quyết định Thủ công, định tính Mang tính nghiên cứu học thuật Định lượng thứ tự ưu tiên bằng số qua ma trận Lohani
  1. Đóng góp vào cơ sở dữ liệu môi trường vùng Đông Nam Bộ: Cung cấp bộ dữ liệu nền phân tích chi tiết $pH$, $DO$, $BOD_5$, $TSS$, vi sinh và phèn sắt ($Fe$) chuẩn hóa theo hệ tọa độ VN-2000/WGS-84 cho Chi cục Bảo vệ Môi trường Tây Ninh, Long An và các viện nghiên cứu tài nguyên nước.

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

Giải pháp Quản lý và Tổ chức

  1. Thành lập Ban điều phối Quản lý Tổng hợp Lưu vực Sông Vàm Cỏ Đông: Cơ quan thường trực liên tỉnh (Tây Ninh - Long An - TP.HCM) nhằm thống nhất quy chế cấp phép xả thải, giám sát lưu lượng xả thải tự động của các doanh nghiệp vượt quá $500 \text{ m}^3/\text{ngày}$.
  2. Quy hoạch phân vùng chức năng tiếp nhận nguồn nước:
    • Vùng thượng lưu (từ biên giới đến Gò Dầu): Vùng bảo tồn nghiêm ngặt để cấp nước sinh hoạt.
    • Vùng trung lưu (từ Hiệp Hòa đến Bến Lức): Vùng kiểm soát tiếp nhận nước thải công nghiệp có điều kiện xả thải loại A theo QCVN 40:2011/BTNMT.
    • Vùng hạ lưu (từ Bến Lức đến cửa sông): Phân vùng kiểm soát xâm nhập mặn phục vụ cấp nước thủy sản và giao thông thủy.
graph LR
    subgraph Ban Quản Lý Lưu Vực Liên Tỉnh
        A[Sở TN&MT Tây Ninh] --- C[Trung Tâm Giám Sát Tự Động]
        B[Sở TN&MT Long An] --- C
    end
    C --> D[Cảnh báo sớm ô nhiễm & Điều tiết nước Hồ Dầu Tiếng]
    C --> E[Thanh tra cống xả tự động KCN]

Giải pháp Công nghệ và Kỹ thuật

SƠ ĐỒ CÔNG NGHỆ XỬ LÝ NƯỚC THẢI KHU DÂN CƯ TẬP TRUNG (FAST - CÔNG SUẤT 500 - 2.000 m³/ngày)
========================================================================================
[Nước thải sinh hoạt] 
      │
      ▼
┌──────────────┐     ┌──────────────┐     ┌───────────────────────┐
│ Song chắn rác│ ──> │   Bể lắng    │ ──> │   Bể sinh học hiếu khí│
│   thô & tinh │     │   cát sơ cấp │     │   giá thể dính bám    │
└──────────────┘     └──────────────┘     │   (FAST Media Module) │
                                          └───────────────────────┘
                                                      │
                                                      ▼
┌──────────────┐     ┌──────────────┐     ┌───────────────────────┐
│ Nguồn tiếp   │ <── │ Bể khử trùng │ <── │      Bể lắng bùn      │
│ nhận sông VCĐ│     │   (Chlorine) │     │       thứ cấp         │
└──────────────┘     └──────────────┘     └───────────────────────┘
  1. Công nghệ xử lý nước thải sinh hoạt phân tán FAST (Fixed Activated Sludge Treatment):
    • Ứng dụng cụm module màng vi sinh ngập nước cố định cho các cụm dân cư ven sông (Hậu Nghĩa, Bến Lức, Gò Dầu), hiệu suất khử BOD đạt $>90%$, TSS đạt $>85%$, không phát sinh mùi hôi, diện tích chiếm đất chỉ bằng $1/3$ hồ sinh học truyền thống.
  2. Hệ thống xử lý nước thải chăn nuôi heo khép kín:
    • Bể Biogas yếm khí phủ bạt HDPE kết hợp hồ sinh học trồng thực vật thủy sinh (bèo lục bình, cỏ vetiver), chuyển hóa $15.000 \text{ m}^3/\text{ngày}$ nước thải phân heo thành khí sinh học phát điện, giảm $95%$ tải lượng BOD xả thẳng ra sông.
  3. Cơ giới hóa thu gom xử lý vấn nạn lục bình:
    • Vận hành cụm tàu vớt cơ giới công suất 50 tấn lục bình/giờ tại khúc sông nghẽn (thị trấn Gò Dầu), kết hợp ủ sinh học lục bình làm phân compost hữu cơ phục vụ nông nghiệp sạch.

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

Hạn chế kỹ thuật

  • Tần suất lấy mẫu chỉ thực hiện 2 đợt/năm (tháng 3 và tháng 9), chưa phản ánh được các biến động đột xuất theo con nước triều cường ngày rằm/mồng một hay các sự cố xả lén ban đêm.
  • Mô hình IDW thuần túy tính toán khoảng cách Euclidean đường thẳng, có thể tạo sai số nhỏ tại các khúc quanh uốn gấp của sông nếu không có mask chặn đường bờ chi tiết.

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

  • Nâng cấp mô hình nội suy IDW sang thuật toán Geostatistical Kriging kết hợp thủy lực dòng chảy 2D (ENVIMWQ / MIKE 21).
  • Tích hợp mạng lưới trạm phao quan trắc tự động liên tục (IoT Water Stations) truyền dữ liệu qua giao thức MQTT/LoRaWAN về trung tâm dữ liệu GIS Server.
  • Xây dựng mô hình trí tuệ nhân tạo (LSTM / Random Forest) dự báo sớm rủi ro cá chết và xâm nhập mặn trước 72 giờ phục vụ người dân lấy nước ngọt.

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

  • Sinh viên & Học viên cao học: Tiếp cận quy trình hoàn chỉnh từ khảo sát hiện trường chuẩn TCVN/APHA đến xử lý dữ liệu nâng cao trên phần mềm GIS và phương pháp toán môi trường.
  • Kỹ sư Môi trường & Chuyên viên GIS: Tham khảo bộ mã nguồn ArcPy và cấu trúc bảng trọng số để tùy biến cho các lưu vực sông khác (sông Sài Gòn, sông Đồng Nai, sông Tiền).
  • Cơ quan Quản lý Nhà nước (Sở TN&MT Tây Ninh, Long An): Sở hữu cơ sở dữ liệu số và bản đồ chuyên đề phục vụ thẩm định báo cáo ĐTM, cấp phép xả thải và quy hoạch các phân khu kinh tế đến năm 2030.
  • Doanh nghiệp & Cộng đồng ven sông: Định hướng cải tiến công nghệ tiền xử lý nước cấp sản xuất, chủ động bảo vệ nguồn nước nuôi trồng thủy sản và giảm thiểu thiệt hại do xâm nhập mặn.

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

1. Yêu cầu cấu hình phần cứng và phần mềm tối thiểu để chạy mô hình GIS nội suy này là gì?

Hệ thống yêu cầu máy tính chạy Windows 7/10 (64-bit), tối thiểu 8GB RAM, bộ vi xử lý Intel Core i5 thế hệ 4 trở lên, cài đặt phần mềm ArcGIS Desktop 10.2 - 10.8 (đã kích hoạt Extension Spatial Analyst) cùng môi trường Python 2.7 đi kèm ArcPy.

2. Tại sao trọng số $T_i$ của $BOD_5$ lại được đặt cao nhất ($0,23$) trong công thức WQI?

Do sông Vàm Cỏ Đông tiếp nhận tải lượng cực lớn từ nước thải sinh hoạt đô thị, nước thải chăn nuôi ($84 \text{ tấn BOD}_5/\text{ngày}$) và các nhà máy chế biến nông sản thực phẩm (mía đường, mì). $BOD_5$ là thước đo trực tiếp phản ánh mức độ suy giảm oxy sinh hóa, quyết định sự sống còn của hệ sinh thái thủy sinh trên sông.

3. Phương pháp nội suy IDW có điểm gì ưu việt hơn Kriging hay Spline trong bài toán này?

IDW bảo toàn giá trị cực trị tại các trạm quan trắc thực tế, thuật toán tính toán nhanh, trực quan, không tạo ra các giá trị ngoại lai âm (như Spline thường gặp) và không đòi hỏi giả định phân phối chuẩn bán biến thiên phức tạp như Geostatistical Kriging trên chuỗi dữ liệu 10 điểm trạm.

4. Tại sao giá trị pH ở đoạn kênh Trà Cú (trạm VCD6) luôn ở mức thấp nhất ($pH \approx 5,0$)?

Khu vực Kênh Trà Cú nằm trọn trong vùng trũng phèn Đồng Tháp Mười mở rộng với hơn 50% diện tích là đất phèn tiềm tàng và hoạt động. Nước mưa rửa trôi các hợp chất pyrite ($FeS_2$) bị oxy hóa từ đồng ruộng vào kênh rạch, tạo thành axit sunfuric ($H_2SO_4$) tự do làm nguồn nước bị chua trầm trọng.

5. Chi phí đầu tư hệ thống xử lý phân tán FAST có khả thi cho các khu dân cư ven sông không?

Công nghệ FAST có chi phí vận hành thấp hơn $40%$ so với bùn hoạt tính truyền thống nhờ không cần hệ thống cào bùn hay tuần hoàn phức tạp. Thời gian hoàn vốn đầu tư xã hội hóa và cải thiện giá trị hệ sinh thái sinh hoạt cho người dân ước tính từ 3 - 5 năm.

Kết luận

Bằng phương pháp tiếp cận khoa học liên ngành, đồ án đã giải quyết trọn vẹn bài toán đánh giá diễn biến chất lượng nước mặt trên toàn tuyến 296 km sông Vàm Cỏ Đông. Việc ứng dụng thành công Chỉ số WQI kết hợp công nghệ GIS nội suy không gian và ma trận đánh giá ưu tiên Lohani đã làm sáng tỏ các quy luật biến động lý hóa: tính axit hóa do đất phèn ở trung lưu, ô nhiễm hữu cơ kiệt quệ DO cục bộ tại các cống xả công nghiệp và xâm nhập mặn $12,6‰$ ở hạ lưu mùa khô. Kết quả nghiên cứu không chỉ mang giá trị học thuật sâu sắc mà còn là luận cứ khoa học thực tiễn vững chắc giúp chính quyền hai tỉnh Tây Ninh và Long An đưa ra các chính sách quản lý lưu vực thông minh, cân bằng giữa phát triển kinh tế công nghiệp và bảo vệ môi trường sinh thái bền vững.