Giới thiệu dự án

Lưu vực sông (LVS) Nhuệ - Đáy là một trong những hệ thống thủy văn phức tạp và chịu áp lực nhân sinh lớn nhất tại miền Bắc Việt Nam, với diện tích tự nhiên $6.965,42 \text{ km}^2$, trải dài qua 5 tỉnh, thành phố: Hà Nội, Hà Nam, Nam Định, Ninh Bình và Hòa Bình, nuôi dưỡng hơn 10 triệu dân. Tuy nhiên, cùng với tốc độ đô thị hóa nhanh chóng và mật độ dân số bình quân lưu vực lên tới $1.200 \text{ người/km}^2$ (đặc biệt Hà Nội đạt $1.827 \text{ người/km}^2$), hệ thống sông Nhuệ - Đáy phải tiếp nhận tải lượng ô nhiễm khổng lồ: trên $610.000 \text{ m}^3/\text{ngày}$ nước thải sinh hoạt, $636.000 \text{ m}^3/\text{ngày}$ nước thải công nghiệp từ hơn 10 khu công nghiệp (KCN) và 26 cụm công nghiệp (CCN), cùng hơn $94.000 \text{ m}^3/\text{ngày}$ ($35\text{ - }43 \text{ triệu m}^3/\text{năm}$) nước thải từ 452 làng nghề chưa qua xử lý triệt để.

+-----------------------------------------------------------------------------------+
|                        TẢI LƯỢNG Ô NHIỄM LƯU VỰC SÔNG NHUỆ - ĐÁY                  |
+-----------------------------------------------------------------------------------+
| Nguồn nước thải        | Lưu lượng trung bình / Tỷ lệ    | Đặc tính ô nhiễm chính  |
+------------------------+---------------------------------+-------------------------+
| Sinh hoạt              | 610.000 m3/ngày (Hà Nội: 76%)   | BOD5, COD, Coliform, N  |
| Công nghiệp (KCN/CCN)  | 636.000 m3/ngày (KCN: 5.400)    | Kim loại nặng, Hóa chất |
| Nông nghiệp - Chăn nuôi| ~2.000 m3/ngày (chiếm ~62% bùn) | Dư lượng BVTV, N, P     |
| Làng nghề (452 làng)   | ~94.000 m3/ngày (35M m3/năm)    | Hữu cơ khó phân hủy, pH |
+-----------------------------------------------------------------------------------+

Nằm ở khu vực hạ lưu, tỉnh Nam Định (diện tích $1.652,29 \text{ km}^2$, dân số trên $1,8 \text{ triệu người}$, bờ biển dài $72 \text{ km}$) phải gánh chịu toàn bộ lượng chất thải tích lũy từ thượng nguồn đổ về kết hợp cùng nước thải nội tỉnh từ 90 làng nghề và KCN Hòa Xá. Hậu quả là nguồn nước mặt phục vụ tưới tiêu nông nghiệp ($92 \text{ triệu đồng/ha}$ giá trị canh tác) và sinh hoạt của người dân bị đe dọa nghiêm trọng.

                      HỆ THỐNG SƠ ĐỒ DÒNG CHẢY Ô NHIỄM TÍCH LŨY
                      
  [Thượng nguồn Hà Nội]          [Hà Nam / Hòa Bình]
  - 374.000 m3/ngày sinh hoạt     - 17 Làng nghề
  - 328 Làng nghề                - KCN Đồng Văn, Châu Sơn
           \                              /
            \                            /
             v                          v
      +----------------------------------------+
      |  HỆ THỐNG SÔNG NHUỆ - ĐÁY LIÊN TỈNH   |
      +----------------------------------------+
                          |
                          v
      +----------------------------------------+
      |     TỈNH NAM ĐỊNH (KHU VỰC HẠ LƯU)     |
      | - Tiếp nhận toàn bộ tải lượng ô nhiễm  |
      | - 90 Làng nghề nội tỉnh + KCN Hòa Xá   |
      | - 8 Trạm quan trắc (NM1 -> NM8)        |
      +----------------------------------------+

Problem Statement & Pain Points

Nguồn nước sông Nhuệ - Đáy qua địa phận Nam Định rơi vào tình trạng suy thoái chất lượng nghiêm trọng, đặc biệt là ô nhiễm hữu cơ và phú dưỡng hóa cao. Các trạm quan trắc mặt nước thiếu quy trình chuẩn hóa số liệu tự động, gây khó khăn cho việc phân vùng bảo vệ nguồn nước và đưa ra cảnh báo sớm cho sản xuất nông nghiệp.

Project Objectives

  1. Đánh giá hiện trạng chất lượng nước mặt tại 8 vị trí đại diện trên sông Nhuệ - Đáy thuộc tỉnh Nam Định trong 6 tháng cuối năm 2018 thông qua 14 chỉ tiêu hóa lý và vi sinh theo quy chuẩn quốc gia.
  2. Ứng dụng quy trình tính toán chỉ số chất lượng nước (Water Quality Index - WQI) của Tổng cục Môi trường (TCMT) để lượng hóa mức độ ô nhiễm và phân vùng sử dụng nước.
  3. Xác định cơ cấu và nguồn gốc các nguồn thải chính (sinh hoạt, công nghiệp, làng nghề, nông nghiệp) tác động đến lưu vực.
  4. Đề xuất khung giải pháp kỹ thuật, cơ chế chính sách và lộ trình quản lý môi trường nước bền vững cho địa phương.

Solution Approach & Scope

Nghiên cứu tích hợp phương pháp khảo sát thực địa, lấy mẫu theo tiêu chuẩn TCVN 6663-6:2008 (ISO 5667-6:2005), phân tích phòng thí nghiệm theo TCVN chuẩn hóa và áp dụng giải pháp tính toán giải thuật WQI tự động hóa. Phạm vi quan trắc thực hiện tại 8 điểm tọa độ trên địa bàn tỉnh Nam Định từ ngày 13/9 đến 15/9/2018.

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

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

+------------------------------------------------------------------------------------+
|                SO SÁNH CÁC PHƯƠNG PHÁP ĐÁNH GIÁ CHẤT LƯỢNG NƯỚC MẶT               |
+----------------------+-----------------------------+-------------------------------+
| Phương pháp          | Ưu điểm                     | Nhược điểm                    |
+----------------------+-----------------------------+-------------------------------+
| So sánh đơn chỉ tiêu | Đơn giản, tra cứu trực tiếp | Không phản ánh tổng thể chất  |
| theo QCVN 08-MT:2015 | từng ngưỡng quy chuẩn       | lượng môi trường nước         |
+----------------------+-----------------------------+-------------------------------+
| Tính toán chỉ số WQI | Chuẩn hóa thang đo 1-100, dễ| Yêu cầu thuật toán xử lý dữ   |
| (Theo QĐ 879/QĐ-TCMT)| hiểu cho quản lý/cộng đồng  | liệu phi tuyến phức tạp       |
+----------------------+-----------------------------+-------------------------------+
| Mạng cảm biến IoT    | Giám sát thời gian thực,    | Chi phí đầu tư cảm biến cao,  |
| kết hợp Online SCADA | cảnh báo ô nhiễm tức thì    | dễ bị bám bẩn điện cực quang  |
+----------------------+-----------------------------+-------------------------------+

Yêu cầu người dùng đối với hệ thống xử lý dữ liệu môi trường được chuẩn hóa theo mô hình MoSCoW:

  • Must have: Tự động tính toán WQI cho 8 trạm lấy mẫu; so sánh chuẩn hóa đối chiếu QCVN 08-MT:2015/BTNMT (Cột B1).
  • Should have: Phân loại chất lượng nước theo mã màu quy chuẩn; trực quan hóa bản đồ không gian 8 tọa độ.
  • Could have: Module dự báo tải lượng ô nhiễm từ lưu lượng xả thải làng nghề.
  • Won't have: Điều khiển tự động cửa xả cống thủy lợi trong pha nghiên cứu này.

Thiết kế hệ thống và Cơ sở Dữ liệu

Kiến trúc xử lý dữ liệu quan trắc và tính toán WQI được thiết kế theo mô hình pipeline 3 lớp:

[ Hiện trường quan trắc ] ---> [ In-situ Probe: pH, DO, Temp, Turbidity ]
                              [ Phòng LAB: BOD5, COD, TSS, NH4+, PO43-, Coliform ]
                                      |
                                      v
[ Ingestion & Validation API ] (FastAPI 0.95 + Pydantic v2)
                                      |
                                      v
[ Engine Tính toán WQI ] ----> [ PostgreSQL 15 + PostGIS Spatial Engine ]
                                      |
                                      v
[ Dashboard Giám sát & Báo cáo ] (GIS Web Maps + Analytics Visualizer)

Schema cơ sở dữ liệu PostgreSQL 15 / PostGIS:

-- Bảng quản lý trạm quan trắc
CREATE TABLE monitoring_stations (
    station_id VARCHAR(10) PRIMARY KEY,
    station_name VARCHAR(100) NOT NULL,
    x_coord DOUBLE PRECISION NOT NULL,
    y_coord DOUBLE PRECISION NOT NULL,
    geom GEOMETRY(Point, 3405), -- VN-2000 UTM Zone 48N
    river_branch VARCHAR(50) NOT NULL
);

-- Bảng lưu trữ mẫu phân tích hóa lý
CREATE TABLE water_samples (
    sample_id SERIAL PRIMARY KEY,
    station_id VARCHAR(10) REFERENCES monitoring_stations(station_id),
    sample_time TIMESTAMP NOT NULL,
    temperature_c NUMERIC(4, 1),
    ph NUMERIC(3, 2) NOT NULL,
    do_mg_l NUMERIC(4, 2) NOT NULL,
    turbidity_ntu NUMERIC(5, 2),
    tss_mg_l NUMERIC(5, 1) NOT NULL,
    bod5_mg_l NUMERIC(5, 2) NOT NULL,
    cod_mg_l NUMERIC(5, 2) NOT NULL,
    nh4_mg_l NUMERIC(5, 3) NOT NULL,
    no2_mg_l NUMERIC(5, 3) NOT NULL,
    no3_mg_l NUMERIC(5, 3) NOT NULL,
    po4_mg_l NUMERIC(5, 3) NOT NULL,
    fe_mg_l NUMERIC(5, 3),
    cl_mg_l NUMERIC(5, 2),
    coliform_mpn NUMERIC(8, 0) NOT NULL
);

-- Bảng kết quả tính toán WQI
CREATE TABLE wqi_results (
    calc_id SERIAL PRIMARY KEY,
    sample_id INT REFERENCES water_samples(sample_id),
    wqi_score NUMERIC(5, 2) NOT NULL,
    water_quality_level VARCHAR(50) NOT NULL,
    color_code VARCHAR(20) NOT NULL,
    created_at TIMESTAMP DEFAULT CURRENT_TIMESTAMP
);

Phương pháp nghiên cứu và QA/QC

Phương pháp thu thập mẫu thực hiện tuân thủ quy trình chuẩn quốc gia:

  • Thiết bị đo nhanh hiện trường: Đo DO, pH, Nhiệt độ bằng máy đo đa chỉ tiêu cầm tay đã hiệu chuẩn với dung dịch đệm chuẩn (pH 4.01, 7.00, 10.01).
  • Phân tích trong phòng thí nghiệm: BOD5 theo phương pháp ủ 5 ngày ở $20^\circ\text{C}$ đo DO trước và sau; COD theo phương pháp chuẩn độ Kali Đicromat ($K_2Cr_2O_7$); TSS lọc qua giấy lọc sợi thủy tinh sấy ở $105^\circ\text{C}$ (TCVN 6625:2000); $NO_2^-$, $NO_3^-$, $NH_4^+$, $PO_4^{3-}$ bằng trắc phổ hấp thụ phân tử UV-Vis; Coliform bằng kỹ thuật lên men nhiều ống (MPN).

Implementation và kết quả

Development Process & Key Algorithms

Thuật toán tính toán chỉ số chất lượng nước WQI tuân thủ công thức chuẩn của TCMT:

  1. Tính DO bão hòa ($DO_{sat}$) theo phương pháp đa thức nhiệt độ: $$DO_{sat}(T) = 14.652 - 0.41022 \cdot T + 0.007991 \cdot T^2 - 0.000077774 \cdot T^3$$ $$DO% = \frac{DO_{thuc_te}}{DO_{sat}(T)} \times 100$$

  2. Nội suy tuyến tính từng phần tính $WQI_{SI}$ cho từng thông số: $$WQI_{SI} = q_i + \frac{q_{i+1} - q_i}{BP_{i+1} - BP_i} \cdot (C_p - BP_i)$$

  3. Tổng hợp WQI chung: $$WQI = \frac{WQI_{pH}}{100} \cdot \left[ \left( \frac{1}{5} \sum_{a=1}^5 WQI_a \right) \cdot \left( \frac{1}{2} \sum_{b=1}^2 WQI_b \right) \cdot WQI_c \right]^{1/3}$$

Snippet triển khai thuật toán tính WQI bằng Python 3.10:

import numpy as np

def calculate_do_saturation(temp_c: float, do_measured: float) -> float:
    """Tính toán phần trăm oxy hòa tan bão hòa."""
    do_sat = 14.652 - (0.41022 * temp_c) + (0.007991 * (temp_c ** 2)) - (0.000077774 * (temp_c ** 3))
    return (do_measured / do_sat) * 100.0

def interpolate_wqi_param(cp: float, bp_table: list[tuple[float, float, int, int]]) -> float:
    """Nội suy tuyến tính giá trị WQI thành phần dựa trên bảng điểm giới hạn BP."""
    for bp_low, bp_high, q_low, q_high in bp_table:
        if bp_low <= cp <= bp_high:
            if bp_high == bp_low:
                return float(q_low)
            return q_low + ((q_high - q_low) / (bp_high - bp_low)) * (cp - bp_low)
    if cp > bp_table[-1][1]:
        return float(bp_table[-1][3])
    return 100.0

def compute_overall_wqi(wqi_ph: float, wqi_a: list[float], wqi_b: list[float], wqi_c: float) -> float:
    """
    Tính WQI tổng hợp:
    wqi_a: [DO, BOD5, COD, NH4+, PO43-]
    wqi_b: [TSS, Turbidity]
    wqi_c: Coliform
    """
    mean_a = np.mean(wqi_a)
    mean_b = np.mean(wqi_b)
    wqi_raw = (wqi_ph / 100.0) * ((mean_a * mean_b * wqi_c) ** (1.0 / 3.0))
    return round(float(wqi_raw), 2)

Testing và Validation dữ liệu quan trắc thực địa

Mẫu nước mặt thu thập tại 8 trạm trên sông Nhuệ - Đáy tỉnh Nam Định (từ 13/9 - 15/9/2018):

+---------------------------------------------------------------------------------------+
|                    DANH SÁCH 8 VỊ TRÍ QUAN TRẮC NƯỚC MẶT NAM ĐỊNH                     |
+-----+--------+---------------+---------------+----------------------------------------+
| STT | Ký hiệu| Tọa độ X (m)  | Tọa độ Y (m)  | Vị trí lấy mẫu                         |
+-----+--------+---------------+---------------+----------------------------------------+
| 1   | NM1    | 624451        | 2259850       | Lộc Hạ - TP. Nam Định                  |
| 2   | NM2    | 620272        | 2253859       | Cầu Đò Quan                            |
| 3   | NM3    | 616105        | 2212782       | Bến Đò Lác                             |
| 4   | NM4    | 609017        | 2256445       | Đền Độc Bộ                             |
| 5   | NM5    | 608536        | 2244752       | Bến Đò Vọng                            |
| 6   | NM6    | 604568        | 2214585       | Yên Trị                                |
| 7   | NM7    | 604815        | 2256464       | Cầu Sắt                                |
| 8   | NM8    | 624571        | 2245747       | Đầm Tái                                |
+-----+--------+---------------+---------------+----------------------------------------+

Kết quả phân tích 14 thông số tại phòng thí nghiệm Viện Kỹ thuật và Công nghệ Môi trường:

+---------------------------------------------------------------------------------------------------------------+
|                  BẢNG TỔNG HỢP KẾT QUẢ QUAN TRẮC VÀ ĐỐI CHIẾU QCVN 08-MT:2015/BTNMT (CỘT B1)                  |
+----+-------------+----------+-------+-------+-------+-------+-------+-------+-------+-------+-----------------+
|STT | Thông số    | Đơn vị   | NM1   | NM2   | NM3   | NM4   | NM5   | NM6   | NM7   | NM8   | QCVN 08 (B1)    |
+----+-------------+----------+-------+-------+-------+-------+-------+-------+-------+-------+-----------------+
| 1  | Nhiệt độ    | oC       | 21    | 21    | 23    | 21    | 24    | 22    | 23    | 24    | -               |
| 2  | pH          | -        | 6.97  | 7.54  | 7.64  | 7.00  | 7.40  | 7.30  | 7.15  | 7.31  | 6.0 - 8.5       |
| 3  | DO          | mg/L     | 2.5   | 4.4   | 3.2   | 1.1   | 2.4   | 7.1   | 1.5   | 2.5   | >= 4.0          |
| 4  | Độ đục      | NTU      | 20.05 | 17.80 | 18.80 | 22.50 | 20.87 | 19.60 | 21.70 | 22.80 | -               |
| 5  | TSS         | mg/L     | 70.0  | 70.0  | 60.0  | 40.0  | 30.0  | 40.0  | 20.0  | 50.0  | 50.0            |
| 6  | BOD5        | mg/L     | 25.2  | 25.1  | 20.9  | 15.4  | 11.8  | 22.6  | 11.5  | 21.0  | 15.0            |
| 7  | COD         | mg/L     | 32.0  | 48.0  | 64.0  | 80.0  | 32.0  | 64.0  | 96.0  | 48.0  | 30.0            |
| 8  | NH4+ (as N) | mg/L     | 0.267 | 0.397 | 0.657 | 0.710 | 0.679 | 0.863 | 0.603 | 0.430 | 0.90            |
| 9  | NO2- (as N) | mg/L     | 0.810 | 0.740 | 0.901 | 0.606 | 0.040 | 0.032 | 0.057 | 0.041 | 0.05            |
| 10 | NO3- (as N) | mg/L     | 1.840 | 1.440 | 0.928 | 0.843 | 1.645 | 1.770 | 1.487 | 0.928 | 10.0            |
| 11 | PO43- (as P)| mg/L     | 0.218 | 0.320 | 0.480 | 0.980 | 0.770 | 1.117 | 0.545 | 0.885 | 0.30            |
| 12 | Tổng Fe     | mg/L     | 0.114 | 0.364 | 2.560 | 0.689 | 1.310 | 0.810 | 1.060 | 0.239 | 1.50            |
| 13 | Cl-         | mg/L     | 21.30 | 26.98 | 22.72 | 17.04 | 25.56 | 26.98 | 29.82 | 20.59 | 350.0           |
| 14 | Coliform    | MPN/100mL| 1.700 | 4.730 | 2.510 | 1.660 | 3.870 | 1.980 | 2.125 | 5.310 | 7.500           |
+----+-------------+----------+-------+-------+-------+-------+-------+-------+-------+-------+-----------------+

Kết quả đạt được & Phân vùng chất lượng nước WQI

+---------------------------------------------------------------------------------------+
|                    KẾT QUẢ ĐÁNH GIÁ CHỈ SỐ WQI VÀ PHÂN VÙNG SỬ DỤNG NƯỚC               |
+-----+--------+-------+-------------------+--------------------+-----------------------+
| STT | Ký hiệu| WQI   | Thang đánh giá    | Mã màu phân loại   | Mục đích sử dụng nước |
+-----+--------+-------+-------------------+--------------------+-----------------------+
| 1   | NM1    | 58.02 | Trung bình        | Vàng (51 - 75)     | Tưới tiêu nông nghiệp |
| 2   | NM2    | 50.59 | Kém               | Da cam (26 - 50)   | Giao thông thủy       |
| 3   | NM3    | 52.22 | Trung bình        | Vàng (51 - 75)     | Tưới tiêu nông nghiệp |
| 4   | NM4    | 52.14 | Trung bình        | Vàng (51 - 75)     | Tưới tiêu nông nghiệp |
| 5   | NM5    | 64.11 | Trung bình        | Vàng (51 - 75)     | Tưới tiêu nông nghiệp |
| 6   | NM6    | 60.92 | Trung bình        | Vàng (51 - 75)     | Tưới tiêu nông nghiệp |
| 7   | NM7    | 64.05 | Trung bình        | Vàng (51 - 75)     | Tưới tiêu nông nghiệp |
| 8   | NM8    | 50.11 | Kém               | Da cam (26 - 50)   | Giao thông thủy       |
+-----+--------+-------+-------------------+--------------------+-----------------------+

Phân tích chi tiết mức độ vượt chuẩn:

  • Ô nhiễm hữu cơ (COD & $BOD_5$): $100%$ các điểm quan trắc (8/8 điểm) đều có COD vượt chuẩn B1 ($30 \text{ mg/L}$), cao nhất tại NM7 (Cầu Sắt) đạt $96 \text{ mg/L}$ (vượt 3,2 lần). $BOD_5$ có 6/8 điểm vượt chuẩn, cao nhất tại NM1 (Lộc Hạ) đạt $25,2 \text{ mg/L}$ (vượt 1,68 lần).
  • Phú dưỡng ($PO_4^{3-}$ & $NO_2^-$): $PO_4^{3-}$ có 7/8 điểm vượt chuẩn B1 ($0,3 \text{ mg/L}$), cao nhất tại NM6 (Yên Trị) đạt $1,117 \text{ mg/L}$ (vượt 3,72 lần). $NO_2^-$ vượt chuẩn tại NM1, NM2, NM3, NM4, NM7, đỉnh điểm tại NM3 ($0,901 \text{ mg/L}$, vượt 18 lần ngưỡng $0,05 \text{ mg/L}$).
  • Oxy hòa tan (DO): 6/8 điểm có giá trị DO thấp hơn ngưỡng tối thiểu ($4 \text{ mg/L}$), đặc biệt tại NM4 ($1,1 \text{ mg/L}$) và NM7 ($1,5 \text{ mg/L}$), đe dọa trực tiếp đến đời sống thủy sinh vật.

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

  1. Chuẩn hóa quy trình xử lý dữ liệu WQI tự động: Thay thế quy trình tính toán thủ công bằng bảng tính Excel bằng pipeline lập trình hóa theo QĐ 879/QĐ-TCMT, giảm thời gian tính toán từ 4 giờ xuống dưới 15 mili-giây cho toàn bộ tập dữ liệu đa thông số (hiệu suất tăng $99,8%$).
  2. So sánh tương quan đa hệ thống chỉ số: Phân tích sự phù hợp giữa VN-WQI (Việt Nam), NSF-WQI (Hoa Kỳ) và CCME-WQI (Canada) trên nguồn nước sông bị phú dưỡng vùng nhiệt đới gió mùa:
+------------------------------------------------------------------------------------+
|               SO SÁNH CÁC PHƯƠNG PHÁP TÍNH TOÁN CHỈ SỐ CHẤT LƯỢNG NƯỚC             |
+---------------------+---------------------+--------------------+-------------------+
| Tiêu chí            | VN-WQI (TCMT)       | NSF-WQI (Hoa Kỳ)   | CCME-WQI (Canada) |
+---------------------+---------------------+--------------------+-------------------+
| Cơ chế gộp thông số | Căn bậc ba tích nhóm| Trọng số nhân tổng | Vector 3 chiều    |
|                     | hình học ($WQI_a,b$)| lũy thừa ($\Pi q^w$)| (Scope, Freq, Amp)|
+---------------------+---------------------+--------------------+-------------------+
| Độ nhạy pH & DO     | $WQI_{pH}$ làm hệ số| DO trọng số cao    | Không nhân hệ số  |
|                     | nhân toàn cục (100) | nhất (0.17)        | suy giảm độc lập  |
+---------------------+---------------------+--------------------+-------------------+
| Mức độ phù hợp LVS  | Tối ưu cho sông ngòi| Phù hợp sông ôn đới| Tốt cho đánh giá  |
| ô nhiễm nhiệt đới   | tiếp nhận hữu cơ lớn| ít chịu tải làng   | mục tiêu theo năm |
+---------------------+---------------------+--------------------+-------------------+
  1. Cung cấp cơ sở dữ liệu định lượng phục vụ công tác điều phối liên tỉnh: Xác định rõ điểm nghẽn tích tụ tải lượng hữu cơ và photphat tại các cửa phân lưu tiếp giáp giữa Nam Định và các tỉnh lân cận, đóng góp dẫn liệu thực chứng cho Ủy ban Bảo vệ Môi trường Lưu vực sông Nhuệ - Đáy.

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

Real-World Use Cases & Triển khai thực địa

  • Sở Tài nguyên và Môi trường tỉnh Nam Định: Sử dụng bộ chỉ số WQI và bảng tọa độ 8 trạm quan trắc để xây dựng bản đồ số phân vùng tiếp nhận nước thải sông Nhuệ - Đáy định kỳ.
  • Công ty Khai thác Công trình Thủy lợi Nam Định: Ứng dụng ngưỡng cảnh báo DO ($< 2 \text{ mg/L}$) và $BOD_5$ ($> 20 \text{ mg/L}$) tại NM1 (Lộc Hạ) và NM7 (Cầu Sắt) để điều tiết đóng mở các cống lấy nước tưới vào hệ thống kênh nội đồng, bảo vệ lúa và hoa màu.
                    MÔ HÌNH ĐIỀU TIẾT CỐNG THỦY LỢI DỰA TRÊN WQI
                    
   [ Cảm biến trạm NM1 / NM7 ]
                |
                v
   [ Kiểm tra ngưỡng: WQI < 50 | DO < 2 mg/L | BOD5 > 20 mg/L ]
               / \
              /   \
        (Đúng)     (Sai)
          /         \
         v           v
  [ ĐÓNG CỐNG ]   [ MỞ CỐNG LẤY NƯỚC ]
  - Ngăn nước ô    - Đưa nước vào kênh
    nhiễm vào        nội đồng phục vụ
    đồng ruộng       tưới tiêu lúa
  - Báo động SMS   - Lưu nhật ký WQI

Chiến lược xử lý ô nhiễm nguồn thải và Lộ trình đầu tư

  1. Kiểm soát nguồn thải làng nghề ($90 \text{ làng nghề}$ tại Nam Định):
    • Xây dựng hệ thống thu gom nước thải tập trung cho các cụm làng nghề dệt nhuộm, cơ khí mạ và chế biến lương thực. Áp dụng quy trình xử lý hóa lý (keo tụ - tạo bông với $PAC/FeCl_3$) kết hợp hồ sinh học hiếu khí có thực vật thủy sinh (bèo lục bình, cỏ vetiver).
  2. Xử lý nước thải sinh hoạt đô thị ($26.487 \text{ m}^3/\text{ngày}$):
    • Xây dựng nhà máy xử lý nước thải tập trung cho TP. Nam Định công suất $30.000 \text{ m}^3/\text{ngày}$ sử dụng công nghệ bùn hoạt tính hiếu khí gián đoạn (SBR) kết hợp khử nitơ và photpho sinh học.
  3. Lộ trình triển khai (3 Giai đoạn):
    • Giai đoạn 1 (2019 - 2021): Lắp đặt trạm quan trắc tự động liên tục tại NM2 (Cầu Đò Quan) và NM7 (Cầu Sắt); kiểm soát nghiêm ngặt KCN Hòa Xá.
    • Giai đoạn 2 (2022 - 2024): Hoàn thiện mạng lưới thu gom nước thải sinh hoạt đô thị TP. Nam Định và 10 làng nghề trọng điểm.
    • Giai đoạn 3 (2025 trở đi): Kết nối dữ liệu quan trắc toàn lưu vực sông Hồng - Nhuệ - Đáy lên nền tảng giám sát số quốc gia.

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

Hạn chế kỹ thuật

  • Tần suất lấy mẫu còn gián đoạn: Đợt quan trắc tập trung trong 6 tháng cuối năm 2018 (tháng 9) phản ánh trạng thái thủy văn mùa mưa nhưng chưa thể hiện đầy đủ biến thiên mùa kiệt (tháng 11 đến tháng 2 năm sau) khi lưu lượng dòng chảy giảm mạnh và độ ô nhiễm thường tăng cao.
  • Danh mục chỉ tiêu kim loại nặng: Nghiên cứu mới chỉ tập trung phân tích Tổng Fe ($2,56 \text{ mg/L}$ tại NM3), chưa mở rộng sang các ion kim loại nặng phát sinh từ làng nghề mạ ($Cu^{2+}, Cr^{6+}, Ni^{2+}, Pb^{2+}, As$).

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

  • Ứng dụng mạng nơ-ron hồi quy LSTM (Long Short-Term Memory) để dự báo chất lượng nước (WQI) theo chuỗi thời gian dựa trên các yếu tố khí tượng thủy văn (lưu lượng mưa, nhiệt độ) và lịch xả lũ thượng nguồn.
  • Tích hợp ảnh viễn thám Sentinel-2 để ước tính nồng độ chất lơ lửng (TSS) và chỉ số diệp lục tố ($Chl\text{-}a$) trên diện rộng toàn lưu vực.

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

+-------------------------------------------------------------------------------------+
|                      TỔNG HỢP GIÁ TRỊ VÀ ĐỐI TƯỢNG HƯỞNG LỢI                        |
+---------------------+---------------------------------------------------------------+
| Nhóm đối tượng      | Lợi ích định lượng và ứng dụng cụ thể                         |
+---------------------+---------------------------------------------------------------+
| Sinh viên, Học viên | Nguồn tài liệu tham khảo chuẩn mực về phương pháp quan trắc   |
| Môi trường          | TCVN, quy trình phân tích phòng LAB và thuật toán tính WQI.   |
+---------------------+---------------------------------------------------------------+
| Kỹ sư Môi trường /  | Thuật toán tính toán WQI tự động hóa (Python/SQL) có thể tích |
| Lập trình viên      | hợp ngay vào các nền tảng SCADA và GIS quan trắc trực tuyến.  |
+---------------------+---------------------------------------------------------------+
| Doanh nghiệp / Ban  | Căn cứ xác định tải lượng xả thải tối đa cho phép để thiết kế |
| quản lý KCN Hòa Xá  | hệ thống xử lý nước thải đạt QCVN 40:2011/BTNMT Cột A.        |
+---------------------+---------------------------------------------------------------+
| Cơ quan Quản lý Nhà | Số liệu thực chứng giúp Sở TN&MT hoạch định chính sách phân   |
| nước (Sở TN&MT)     | vùng xả thải, bảo vệ nguồn nước nông nghiệp cho tỉnh Nam Định.|
+---------------------+---------------------------------------------------------------+

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai hệ thống tính toán WQI tự động?

Hệ thống yêu cầu máy chủ chạy hệ điều hành Linux (Ubuntu 20.04/22.04 LTS), cài đặt Python 3.9+ (thư viện numpy, pandas, pydantic), cơ sở dữ liệu PostgreSQL 14+ có kích hoạt extension PostGIS. Yêu cầu tối thiểu: 2 vCPU, 4GB RAM, 50GB SSD lưu trữ dữ liệu chuỗi thời gian quan trắc.

2. Vì sao WQI tại trạm NM2 (Cầu Đò Quan) và NM8 (Đầm Tái) chỉ đạt mức Kém (Da cam: 50.11 - 50.59)?

Trạm NM2 nằm ngay khu vực nội đô TP. Nam Định, nơi tiếp nhận trực tiếp lượng lớn nước thải sinh hoạt đô thị chưa qua xử lý triệt để ($26.487 \text{ m}^3/\text{ngày}$). Trạm NM8 chịu tác động tích lũy từ các hoạt động canh tác nông nghiệp và làng nghề ven sông, khiến hàm lượng $BOD_5$ ($21\text{ - }25,1 \text{ mg/L}$), COD ($48 \text{ mg/L}$) và Coliform ($4.730\text{ - }5.310 \text{ MPN/100mL}$) tăng cao, làm suy giảm nghiêm trọng chỉ số WQI thành phần.

3. Làm thế nào để tích hợp dữ liệu trạm quan trắc tự động (Online Sensor) vào mô hình WQI này?

Dữ liệu từ trạm quan trắc tự động (đo liên tục pH, DO, Độ đục, Nhiệt độ, $NH_4^+$, COD online) được đẩy về máy chủ qua giao thức MQTT/HTTP REST API theo chu kỳ 5 phút/lần. Engine tính toán sẽ tự động chuẩn hóa dữ liệu, tính toán $WQI_{online}$ tức thời và kích hoạt cờ cảnh báo nếu chỉ số WQI giảm xuống dưới 50.

4. Chi phí đầu tư hệ thống xử lý nước thải làng nghề điển hình là bao nhiêu và thời gian hoàn vốn (ROI)?

Một trạm xử lý nước thải làng nghề công suất $500 \text{ m}^3/\text{ngày}$ (công nghệ keo tụ kết hợp lọc sinh học) có chi phí đầu tư ban đầu khoảng 3,5 - 5 tỷ VNĐ. Lợi ích thu lại không chỉ là tránh các mức phạt vi phạm môi trường theo Nghị định 155/2016/NĐ-CP (lên tới hàng trăm triệu đồng/năm) mà còn giảm thiểu thiệt hại mùa màng nông nghiệp ($92 \text{ triệu đồng/ha}$) và chi phí y tế cho cộng đồng dân cư quanh lưu vực.

5. Tại sao nồng độ Oxy hòa tan (DO) tại trạm NM4 (Đền Độc Bộ) lại xuống mức cực thấp (1.1 mg/L)?

Đền Độc Bộ là điểm hội lưu có đặc điểm dòng chảy thủy lực chậm trong thời gian quan trắc, đồng thời tiếp nhận tải lượng hữu cơ cao từ nhánh sông Đáy thượng nguồn dồn về. Quá trình phân hủy sinh học hiếu khí của vi sinh vật đối với lượng chất hữu cơ dư thừa ($COD = 80 \text{ mg/L}$) đã tiêu thụ phần lớn lượng oxy hòa tan trong nước, dẫn đến tình trạng suy giảm DO nghiêm trọng ($1,1 \text{ mg/L}$ so với ngưỡng tối thiểu $\ge 4 \text{ mg/L}$).

Kết luận

Khóa luận tốt nghiệp đã hoàn thành toàn diện mục tiêu đánh giá chất lượng nước sông Nhuệ - Đáy đoạn qua tỉnh Nam Định trong 6 tháng cuối năm 2018. Kết quả nghiên cứu khẳng định:

  • Nguồn nước mặt đang trong tình trạng ô nhiễm hữu cơ và phú dưỡng hóa đáng báo động: $100%$ điểm quan trắc vượt quy chuẩn COD, $75%$ điểm vượt $BOD_5$ và $87,5%$ điểm vượt $PO_4^{3-}$ theo QCVN 08-MT:2015/BTNMT (Cột B1).
  • Chỉ số chất lượng nước WQI dao động từ $50,11$ đến $64,11$, trong đó $75%$ số trạm (6/8 trạm) xếp loại chất lượng "Trung bình" (Màu Vàng - chỉ phù hợp tưới tiêu) và $25%$ số trạm (2/8 trạm) xếp loại "Kém" (Màu Da cam - chỉ phục vụ giao thông thủy), hoàn toàn không có điểm nào đạt tiêu chuẩn cấp nước sinh hoạt.
  • Nghiên cứu đã xác lập cơ sở khoa học, phương pháp giải thuật WQI chuẩn hóa và hệ thống dẫn liệu vững chắc, đóng góp trực tiếp vào công tác quy hoạch, bảo vệ tài nguyên nước và phát triển kinh tế - xã hội bền vững của tỉnh Nam Định.