Giới thiệu dự án

Hệ thống thủy lợi (HTTL) Liễn Sơn là công trình thủy nông huyết mạch của tỉnh Vĩnh Phúc, giữ vai trò cung cấp nguồn nước tưới tiêu cho 24.876 ha đất canh tác nông nghiệp, 1.709 ha cây lâu năm, 1.788 ha nuôi trồng thủy sản trên tổng diện tích tự nhiên 47.481 ha. Đồng thời, hệ thống còn đảm nhiệm cấp nước cho phường Bạch Hạc (thành phố Việt Trì, tỉnh Phú Thọ) và huyện Mê Linh (Hà Nội). Tuy nhiên, dưới áp lực đô thị hóa và phát triển công nghiệp nhanh chóng, HTTL Liễn Sơn đang chịu tải lượng xả thải lên tới 41.351,66 $\text{m}^3/\text{ngày đêm}$, trong đó có đến 51,00% (21.056,06 $\text{m}^3/\text{ngày đêm}$) chưa qua xử lý xả trực tiếp vào nguồn nước.

+-------------------------------------------------------------------------------+
|                       HỆ THỐNG THỦY LỢI LIỄN SƠN (VĨNH PHÚC)                   |
|       (Diện tích tự nhiên: 47.481 ha | Diện tích tưới: 24.876 ha đất NN)      |
+-------------------------------------------------------------------------------+
                                        |
       +--------------------------------+-------------------------------+
       |                                                                |
       v                                                                v
[NGUỒN CẤP NƯỚC CHÍNH]                                        [ÁP LỰC Ô NHIỄM NƯỚC THẢI]
- Sông Phó Đáy qua Đập Liễn Sơn (17 m3/s)                     Tổng thải: 41.351,66 m3/ngày-đêm
- TB Bạch Hạc (sông Lô: 11,2 m3/s)                            - Chưa xử lý: 51,00% (21.056 m3/ng-đ)
- TB Đại Định (sông Hồng: 11,2 m3/s)                          - Nước thải sinh hoạt: 55,17%
- Trục tiêu & tiếp nước: Sông Phan (>800 km2)                 - KCN / CCN: 22,94% | Cơ sở SX: 13,11%
                                                              - Chăn nuôi: 5,89%  | Y tế: 1,23%

Vấn đề thực tiễn (Problem Statement)

Chất lượng nguồn nước mặt trên hệ thống Liễn Sơn, đặc biệt là lưu vực sông Phan, đang suy thoái nghiêm trọng trong mùa kiệt (tháng 1 đến tháng 3). Nồng độ các chất ô nhiễm hữu cơ và dinh dưỡng vượt quy chuẩn kỹ thuật quốc gia:

  • $\text{BOD}_5$ vượt từ 1,9 đến 2,2 lần so với QCVN 08-MT:2015/BTNMT (Cột B1).
  • $\text{COD}$ vượt từ 2,6 đến 2,8 lần.
  • $\text{N-NH}_4^+$ vượt từ 1,2 đến 1,6 lần.
  • Tình trạng bèo tây, rác thải làng nghề tái chế nhựa, nước thải chăn nuôi chưa qua xử lý gây bồi lắng lòng dẫn và cản trở dòng chảy.

Mục tiêu dự án

  1. Đánh giá toàn diện hiện trạng chất lượng nước (CLN) trên HTTL Liễn Sơn giai đoạn 2018–2020 theo không gian (15 vị trí quan trắc) và thời gian (19 đợt quan trắc mùa khô và mùa mưa).
  2. Kiểm kê, phân loại và định lượng các nguồn thải chính đổ vào hệ thống (sinh hoạt, công nghiệp, làng nghề, chăn nuôi, y tế).
  3. Phân tích diễn biến chỉ số chất lượng nước (WQI) theo Quyết định 1460/QĐ-TCMT và đối sánh theo QCVN 08-MT:2015/BTNMT (Cột B1).
  4. Dự báo xu thế ô nhiễm nguồn nước đến năm 2025 theo 3 kịch bản phát triển kinh tế - xã hội của tỉnh Vĩnh Phúc.
  5. Đề xuất giải pháp công trình và phi công trình nhằm bảo vệ nguồn nước, kiểm soát xả thải và nâng cao khả năng tự làm sạch của hệ thống.

Phạm vi và giới hạn nghiên cứu

  • Phạm vi không gian: Toàn bộ hệ thống trải dài trên 7 huyện/thành phố/thị xã của Vĩnh Phúc (Lập Thạch, Tam Dương, Vĩnh Tường, Vĩnh Yên, Yên Lạc, Bình Xuyên, Phúc Yên). Mạng lưới quan trắc gồm 8 điểm trên kênh chính và 7 điểm trên hệ thống kênh tiêu/sông Phan.
  • Phạm vi thời gian: Thu thập chuỗi số liệu quan trắc 3 năm (2018: 5 đợt; 2019: 5 đợt; 2020: 9 đợt), phân chia theo mùa khô (tháng 1, 2, 3) và mùa mưa (tháng 5, 6, 7).

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

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

Hiện nay, việc đánh giá CLN tại các hệ thống thủy nông Việt Nam thường gặp phải tình trạng thiếu đồng bộ về tần suất lấy mẫu, phương pháp tính toán cũ (Quyết định 879/QĐ-TCMT) chưa phản ánh hết các nhóm thông số kim loại và vi sinh, hoặc thiếu mô hình hóa tải lượng ô nhiễm từ lưu vực.

Phương pháp tiếp cận Ưu điểm Nhược điểm Khả năng ứng dụng tại HTTL Liễn Sơn
Đánh giá đơn chỉ tiêu (QCVN) Rõ ràng, dễ xác định vi phạm theo từng thông số Không đưa ra bức tranh tổng thể đa chỉ tiêu Áp dụng bắt buộc để kiểm định nguồn tưới B1
Chỉ số WQI cũ (QĐ 879/QĐ-TCMT) Tính toán đơn giản, 5 mức phân loại Thiếu các hệ số trọng số mới, bỏ sót chỉ tiêu đặc thù Lỗi thời, cần nâng cấp
Chỉ số WQI cải tiến (QĐ 1460/QĐ-TCMT) 6 mức thang điểm chi tiết, tích hợp 15 thông số hóa lý - vi sinh Yêu cầu số liệu đầu vào nghiêm ngặt và chính xác cao Giải pháp trọng tâm của đề tài

Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)

  • Must have: Phân tích đủ 15 thông số hóa lý/vi sinh/kim loại nặng; chuẩn hóa thuật toán tính $WQI_{SI}$, $WQI_{DO}$, $WQI_{pH}$; xác định chính xác 3 vùng ô nhiễm trọng điểm trên sông Phan.
  • Should have: Dự báo tải lượng chất ô nhiễm theo 3 kịch bản tăng trưởng KT-XH 2025; tính toán cân bằng oxy hòa tan ($\text{DO}$).
  • Could have: Tích hợp bản đồ GIS phân vùng chất lượng nước; sơ đồ công nghệ xử lý sinh học BASTAF và Mương oxy hóa cải tiến.
  • Won't have: Mô hình hóa thủy lực 2D dòng chảy biến đổi liên tục theo thời gian thực (IoT tự động).

Thiết kế hệ thống mạng lưới quan trắc và công nghệ phân tích

+---------------------------------------------------------------------------------+
|                    KIẾN TRÚC MẠNG LƯỚI QUAN TRẮC & ĐÁNH GIÁ CLN                 |
+---------------------------------------------------------------------------------+
| [8 Điểm Kênh Chính]                             [7 Điểm Kênh Tiêu / Sông Phan]  |
| LS1: Đập Liễn Sơn (Đồng Tĩnh)                   LS9: Cầu Vàng (Hoàng Lâu)       |
| LS2: TB Bạch Hạc (Việt Xuân)                    LS10: Cầu Vĩnh Tường (Yên Lập)  |
| LS3: TB Đại Định (Cao Đại)                      LS11: Cầu Xuân Lai (Vũ Di)      |
| LS4: Cống Vân Tập (Vân Hội)                     LS12: Đập Vĩnh Sơn (Vĩnh Sơn)   |
| LS5: Cống Lương Lập (Xuân Lai)                  LS13: Đập Lạc Ý (Đồng Tâm)      |
| LS6: Cống Nguyệt Đức (Thanh Lãng)               LS14: Cầu Lò Cang (Hương Canh)  |
| LS7: Cầu Cơ Khí (Thanh Lãng)                    LS15: Sông Cà Lỗ Cụt (Tiên Châu)|
| LS8: Cầu Đất (Tân Phong)                                                        |
+---------------------------------------------------------------------------------+
                                        |
                                        v
+---------------------------------------------------------------------------------+
|                         GIAO THỨC PHÂN TÍCH CHUẨN TCVN                          |
|  - Độ đục: Turb-350 LR (TCVN 6184:2008)       - TSS: DR5000 630nm (TCVN 6625)   |
|  - DO, pH: Máy TOA WQC-22A (TCVN 7321, 6492)  - COD: K2Cr2O7 (TCVN 6491-1999)  |
|  - BOD5: YSI 5000 (TCVN 6001-1:2008)          - N-NH4+: Indophenol (TCVN 5988)  |
|  - P-PO4: UV-Vis Molypdat (TCVN 6202:2008)    - Cu, Pb, Cd: AAS (TCVN 6193)    |
+---------------------------------------------------------------------------------+

Phương pháp luận (Methodology)

  1. Thu thập và kế thừa số liệu: Chuỗi dữ liệu KT-XH, khí tượng thủy văn (Trạm Vĩnh Yên, Tam Đảo), mạng lưới 230 công trình trên sông chính.
  2. Quy trình lấy và bảo quản mẫu: Tuân thủ TCVN 6663-1:2011 và TCVN 6663-3:2003. Mẫu lấy ở độ sâu $0,5 \text{ m}$ tại tâm dòng chảy ổn định, bảo quản lạnh $4^\circ\text{C}$ tránh ánh sáng và chuyển về Phòng Thí nghiệm Tổng hợp - Viện Nước, Tưới tiêu và Môi trường.
  3. Đánh giá theo Quy chuẩn và Chỉ số:
    • QCVN 08-MT:2015/BTNMT (Cột B1 - dùng cho mục đích tưới tiêu nông nghiệp).
    • Quyết định 1460/QĐ-TCMT phân hạng CLN thành 6 mức: Rất tốt (91-100), Tốt (76-90), Trung bình (51-75), Kém (26-50), Ô nhiễm nặng (10-25), Ô nhiễm rất nặng (<10).

Triển khai và kết quả thực nghiệm

Quy trình tính toán và thuật toán xử lý

Thuật toán nội suy tuyến tính giá trị thông số $WQI_{SI}$ cho từng thông số quy chuẩn theo Quyết định 1460/QĐ-TCMT:

$$WQI_{SI} = \frac{q_{i+1} - q_i}{BP_{i+1} - BP_i} \times (C_p - BP_i) + q_i$$

Trong đó:

  • $C_p$: Nồng độ thực tế của thông số quan trắc.
  • $BP_i, BP_{i+1}$: Nồng độ giới hạn dưới và trên tại khoảng giá trị tương ứng mức $i$ và $i+1$.
  • $q_i, q_{i+1}$: Điểm WQI chuẩn tương ứng tại các mức $BP_i, BP_{i+1}$.

Đối với thông số Oxy hòa tan ($\text{DO}$), nồng độ bão hòa được tính theo hàm phi tuyến bậc 3 của nhiệt độ môi trường nước $T$ ($^\circ\text{C}$):

$$\text{DO}_{\text{bão hòa}} = 14.652 - 0.41022 \cdot T + 0.007991 \cdot T^2 - 0.000077774 \cdot T^3$$

$$\text{DO}{% \text{bão hòa}} = \left( \frac{\text{DO}{\text{quan trắc}}}{\text{DO}_{\text{bão hòa}}} \right) \times 100%$$

def calculate_do_saturation(temp_c: float, do_measured: float) -> float:
    """Tính toán tỷ lệ phần trăm oxy bão hòa DO% theo nhiệt độ nước."""
    do_sat_val = 14.652 - (0.41022 * temp_c) + (0.007991 * (temp_c ** 2)) - (0.000077774 * (temp_c ** 3))
    do_percent = (do_measured / do_sat_val) * 100.0
    return do_percent

def calculate_wqi_do(do_percent: float) -> float:
    """Tính toán chỉ số WQI_DO theo QĐ 1460/QĐ-TCMT."""
    if do_percent <= 20.0 or do_percent >= 200.0:
        return 1.0
    elif 88.0 <= do_percent <= 112.0:
        return 100.0
    elif 20.0 < do_percent < 88.0:
        # Nội suy bậc 1 theo thang quy chuẩn
        bp_low, bp_high, q_low, q_high = 20.0, 88.0, 1.0, 100.0
        return ((q_high - q_low) / (bp_high - bp_low)) * (do_percent - bp_low) + q_low
    else: # 112.0 < do_percent < 200.0
        bp_low, bp_high, q_low, q_high = 112.0, 200.0, 100.0, 1.0
        return ((q_high - q_low) / (bp_high - bp_low)) * (do_percent - bp_low) + q_low

Kiểm chứng thực nghiệm và phân tích dữ liệu quan trắc (2018–2020)

Phân bố tải lượng xả thải theo ngành

Tổng lượng nước thải phát sinh xả vào HTTL Liễn Sơn đạt 41.351,66 $\text{m}^3/\text{ngày đêm}$:

  • Nước thải sinh hoạt: $22.811,18 \text{ m}^3/\text{ngày đêm}$ (chiếm 55,17%). Thành phố Vĩnh Yên chiếm tỷ trọng lớn nhất với mật độ dân cư tập trung.
  • Khu/Cụm công nghiệp (KCN/CCN): $9.486,78 \text{ m}^3/\text{ngày đêm}$ (22,94%). KCN Khai Quang chiếm tới 91,84% lượng thải công nghiệp toàn hệ thống.
  • Cơ sở sản xuất ngoài KCN: $5.422,32 \text{ m}^3/\text{ngày đêm}$ (13,11%).
  • Chăn nuôi gia súc, gia cầm: $2.433,69 \text{ m}^3/\text{ngày đêm}$ (5,89%). Đáng báo động khi 80,28% (1.953,69 $\text{m}^3/\text{ngày đêm}$) chưa qua xử lý, phát sinh từ tổng đàn 51.594 con với 8,84 tấn chất thải rắn/ngày.
  • Cơ sở y tế: $507,97 \text{ m}^3/\text{ngày đêm}$ (1,23%), tập trung 88,56% tại TP. Vĩnh Yên.
  • Làng nghề: $1.016,5 \text{ m}^3/\text{ngày đêm}$ từ 12 làng nghề (tái chế nhựa Yên Đồng, mộc Lũng Hạ, chế biến thực phẩm Hòn Loan).
TỶ LỆ CÁC NGUỒN NƯỚC THẢI XẢ VÀO HTTL LIỄN SƠN
============================================================
Sinh hoạt      [████████████████████████████] 55,17% (22.811 m3/ng-đ)
KCN / CCN      [███████████] 22,94% (9.487 m3/ng-đ)
Ngoài KCN/CCN  [██████] 13,11% (5.422 m3/ng-đ)
Chăn nuôi      [███] 5,89% (2.434 m3/ng-đ)
Y tế & Làng nghề [█] 2,89% (1.524 m3/ng-đ)
============================================================

Đánh giá diễn biến các thông số hóa lý quan trọng

Vị trí quan trắc $\text{DO}$ mùa khô ($\text{mg/L}$) $\text{COD}$ ($\text{mg/L}$) $\text{BOD}_5$ ($\text{mg/L}$) $\text{N-NH}_4^+$ ($\text{mg/L}$) Phân cấp WQI Đánh giá chất lượng nước
LS1 (Đập Liễn Sơn) 5,8 – 6,5 12 – 18 6 – 9 0,15 – 0,35 78 – 86 (Tốt) Đạt chuẩn tưới tiêu B1
LS4 (Cống Vân Tập) 4,9 – 5,4 22 – 28 11 – 14 0,45 – 0,60 62 – 71 (Trung bình) Phù hợp tưới tiêu
LS7 (Cầu Cơ Khí) 1,8 – 2,5 68 – 85 32 – 45 2,10 – 3,40 18 – 24 (Ô nhiễm nặng) Không đạt, nguy hại cho cây trồng
LS8 (Cầu Đất) 1,2 – 2,1 75 – 92 38 – 52 2,80 – 4,10 12 – 19 (Ô nhiễm nặng) Bị ô nhiễm nghiêm trọng bởi nước thải
LS9 (Sông Phan - Cầu Vàng) 3,1 – 3,8 45 – 58 21 – 29 1,20 – 1,85 35 – 44 (Kém) Hạn chế sử dụng, cần pha loãng
LS11 (Cầu Xuân Lai) 1,5 – 2,2 70 – 88 35 – 48 2,40 – 3,60 15 – 22 (Ô nhiễm nặng) Vùng ô nhiễm tập trung
LS13 (Đập Lạc Ý) 0,8 – 1,6 82 – 105 42 – 60 3,50 – 5,20 8 – 14 (Rất nặng) Suy thoái sinh thái trầm trọng 
BIẾN THIÊN NỒNG ĐỘ COD (mg/L) TRÊN HỆ THỐNG MÙA KHÔ
======================================================================
QCVN 08 (B1: 30 mg/L) |--------------------|
LS1 (Đập Liễn Sơn)    |========= 18
LS4 (Cống Vân Tập)    |============= 28
LS7 (Cầu Cơ Khí)      |========================================= 85
LS8 (Cầu Đất)         |============================================= 92
LS13 (Đập Lạc Ý)      |=================================================== 105
======================================================================

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

  1. Đổi mới phương pháp đánh giá tải lượng và phân vùng ô nhiễm:

    • Áp dụng đầy đủ thuật toán Quyết định 1460/QĐ-TCMT trên toàn mạng lưới thủy nông, phân hóa chính xác 3 vùng ô nhiễm trọng điểm trên sông Phan:
      • Vùng 1: Đoạn Cầu Cơ Khí đến Cầu Đất (nhận nước thải làng nghề gỗ, sinh hoạt thị trấn Thanh Lãng).
      • Vùng 2: Trục sông Phan qua TP. Vĩnh Yên (tiếp nhận KCN Khai Quang, nước thải đô thị).
      • Vùng 3: Từ Trạm bơm Bạch Hạc qua Cầu Xuân Lai đến TB Đại Định.
    • Xác lập vùng an toàn cấp nước tưới tự chảy từ Đập dâng Liễn Sơn đến Cống Vân Tập ($WQI > 75$).

  2. So sánh với các nghiên cứu tiền nhiệm tại các đại hệ thống thủy nông:

Tiêu chí phân tích HTTL Bắc Đuống (Viện Nước 2018) HTTL Bắc Hưng Hải (Lê Xuân Quang 2019) HTTL Liễn Sơn (Nghiên cứu này, 2021)
Quy mô lưu vực 41.000 ha tưới / 53.000 ha tiêu 150.000 ha tưới / 185.000 ha tiêu 24.876 ha tưới / 40.279 ha tiêu
Tỷ lệ trạm ô nhiễm 65% số điểm vượt QCVN 80% (12/15 điểm ô nhiễm nặng) 46,7% (7/15 điểm ô nhiễm nặng/kém)
Tác nhân chính 62 làng nghề sông Ngũ Huyện Khê Nước thải sinh hoạt Hà Nội, Hưng Yên Nước thải sinh hoạt (55,17%), KCN (22,94%)
Công cụ đánh giá Đơn thông số, WQI cũ (QĐ 879) WQI cũ, tần suất định kỳ WQI 1460/QĐ-TCMT + Mô hình 3 kịch bản 2025
Mức cải thiện đề xuất Nạo vét cục bộ Gạn tháo nước đệm Mương oxy hóa + BASTAF bãi lọc ngầm + Điều tiết Đập
  1. Đề xuất giải pháp công nghệ sinh thái thích ứng địa phương:
    • Khu vực đô thị/thị trấn: Mô hình Mương oxy hóa cải tiến (Carrousel / Modified Oxidation Ditch) với hiệu suất xử lý $\text{BOD}_5 > 90%$, $\text{COD} > 85%$, khử Nitơ/Photpho thông qua kiểm soát chu kỳ thiếu khí - hiếu khí xen kẽ.
    • Khu vực nông thôn và làng nghề: Dây chuyền Module BASTAF (Bể tự hoại cải tiến có vách ngăn mỏng và ngăn lọc kỵ khí hướng lên) kết hợp Bãi lọc ngầm trồng cây dòng chảy ngầm ngang/đứng (Subsurface Flow Constructed Wetland) sử dụng sậy (Phragmites australis) và cỏ Vetiver. Giảm chi phí vận hành 65 - 75% so với công nghệ cơ điện truyền thống.

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

Kịch bản diễn biến tải lượng ô nhiễm đến năm 2025

Nghiên cứu xây dựng 3 kịch bản phát triển KT-XH tỉnh Vĩnh Phúc đến năm 2025 để dự báo lưu lượng nước thải xả vào HTTL Liễn Sơn:

DỰ BÁO TẢI LƯỢNG NƯỚC THẢI HTTL LIỄN SƠN ĐẾN NĂM 2025
+-------------------------------------------------------------------------------+
| Kịch bản 1: Phát triển tự nhiên, hạ tầng xử lý chậm tiến độ                   |
| -> Tổng nước thải: 58.450 m3/ng-đ | Chưa xử lý: 64,2% -> WQI sông Phan < 15   |
+-------------------------------------------------------------------------------+
| Kịch bản 2: Đầu tư theo Quy hoạch thoát nước đô thị Vĩnh Phúc                 |
| -> Tổng nước thải: 62.100 m3/ng-đ | Đã xử lý: 68,5%   -> WQI phục hồi 45 - 55 |
+-------------------------------------------------------------------------------+
| Kịch bản 3: Xử lý tập trung KCN 100% + BASTAF nông thôn đạt 80% (Tối ưu)      |
| -> Tổng nước thải: 65.800 m3/ng-đ | Đã xử lý: 88,2%   -> WQI toàn tuyến > 65  |
+-------------------------------------------------------------------------------+

DỰ BÁO TỔNG TẢI LƯỢNG NƯỚC THẢI ĐƯỢC XỬ LÝ NĂM 2025 (m3/ngày-đêm)
======================================================================
Kịch bản 1 (35,8%) |================== 20.925 / 58.450
Kịch bản 2 (68,5%) |================================== 42.538 / 62.100
Kịch bản 3 (88,2%) |============================================ 58.035 / 65.800
======================================================================

Lộ trình triển khai giải pháp công trình và vận hành liên hồ chứa

LỘ TRÌNH TRIỂN KHAI BẢO VỆ NGUỒN NƯỚC HTTL LIỄN SƠN (2021 - 2025)
--------------------------------------------------------------------------------
Giai đoạn 1 (2021 - 2022): Kiểm soát & Tối ưu vận hành
* Kiểm toán 100% cơ sở xả thải trên lưu vực sông Phan và KCN Khai Quang.
* Thiết lập quy trình vận hành đóng mở liên động Đập Liễn Sơn và cống Vân Tập để tăng cường lưu lượng xả rửa đáy mùa kiệt (Q_xa duy trì > 5 m3/s).

Giai đoạn 2 (2022 - 2024): Xây dựng hạ tầng xử lý phân tán
* Xây dựng trạm XLNT tập trung tại TT. Thanh Lãng (công suất 2.500 m3/ng-đ) và TT. Thổ Tang (3.000 m3/ng-đ).
* Lắp đặt thí điểm 15 cụm BASTAF kết hợp bãi lọc ngầm tại các điểm nóng chăn nuôi xã Yên Đồng và Hải Lựu.

Giai đoạn 3 (2024 - 2025): Hoàn thiện mạng lưới quan trắc tự động
* Lắp đặt 5 trạm quan trắc tự động liên tục (Online Telemetry) tại Đập Liễn Sơn, Cầu Đất, Cầu Cơ Khí, Đập Lạc Ý, Cống Nguyệt Đức truyền dữ liệu về Chi cục Thủy lợi Vĩnh Phúc.
--------------------------------------------------------------------------------

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

Hạn chế kỹ thuật

  • Phương pháp lấy mẫu định kỳ thủ công (Grab sampling) chưa phản ánh liên tục các biến động xả thải trộm vào ban đêm tại các làng nghề phôi giấy Hợp Thịnh và tái chế nhựa.
  • Nguồn ô nhiễm diện (phân bón hóa học, thuốc BVTV từ đồng ruộng) chưa được mô hình hóa chi tiết bằng các công cụ chuyên sâu như SWAT (Soil and Water Assessment Tool).

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

  • Tích hợp mô hình toán thủy lực - chất lượng nước 1D/2D (MIKE 11 / Delft3D) kết hợp với GIS để mô phỏng sự lan truyền chất ô nhiễm $\text{BOD}_5, \text{NH}_4^+$ theo thời gian thực khi vận hành đóng mở cống tiêu.
  • Xây dựng hệ thống cảnh báo sớm ô nhiễm nước tưới tiêu ứng dụng cảm biến IoT năng lượng mặt trời kết nối mạng truyền thông không dây LoRaWAN/4G.

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

  • Cơ quan quản lý nhà nước & Doanh nghiệp khai thác công trình thủy lợi (Công ty TNHH MTV Thủy lợi Liễn Sơn, Chi cục Thủy lợi Vĩnh Phúc): Cung cấp cơ sở khoa học và bộ số liệu quan trắc thực chứng để xây dựng quy trình vận hành điều tiết nước, cấp phép xả thải và bảo vệ nguồn nước theo Luật Thủy lợi 2017.
  • Hộ nông dân & Hợp tác xã nông nghiệp: Đảm bảo nguồn nước tưới đạt chuẩn QCVN 08-MT:2015/BTNMT (Cột B1), ngăn ngừa hiện tượng cháy lá, ngộ độc hữu cơ và tồn dư kim loại nặng trên lúa và rau màu.
  • Cộng đồng dân cư ven sông Phan: Giảm thiểu mùi hôi thối, ô nhiễm vi sinh gây dịch bệnh tiêu hóa, phục hồi cảnh quan sinh thái mặt nước tại Đầm Vạc và hạ lưu sông Cà Lỗ.
  • Các nhà nghiên cứu, sinh viên chuyên ngành Kỹ thuật Tài nguyên nước, Môi trường: Cung cấp tài liệu tham khảo thực tế về phương pháp tích hợp WQI theo chuẩn QĐ 1460/QĐ-TCMT và thiết kế hệ sinh thái xử lý nước thải phân tán.

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

1. Tại sao chất lượng nước sông Phan lại suy thoái nghiêm trọng hơn nhiều so với kênh chính Tả Ngạn Liễn Sơn?

Kênh chính Tả Ngạn nhận trực tiếp nước bổ cập từ sông Phó Đáy qua Đập Liễn Sơn với lưu lượng thiết kế $17 \text{ m}^3/\text{s}$ và được tiếp nguồn bởi trạm bơm Bạch Hạc ($11,2 \text{ m}^3/\text{s}$), dòng chảy liên tục nên khả năng pha loãng cao. Ngược lại, sông Phan đóng vai trò trục tiêu thoát, lưu thông kém vào mùa khô (dòng chảy tù hãm), nhưng lại tiếp nhận hơn 80% tổng lượng nước thải sinh hoạt đô thị, nước thải KCN Khai Quang, nước thải làng nghề và chăn nuôi chưa qua xử lý.

2. Phương pháp tính toán WQI theo QĐ 1460/QĐ-TCMT có điểm gì ưu việt hơn QĐ 879/QĐ-TCMT trước đây?

Quyết định 1460/QĐ-TCMT chia thang đánh giá thành 6 cấp độ chi tiết (thay vì 5 cấp), phân nhóm rõ ràng các thông số: Nhóm I (pH), Nhóm II (Thông số kim loại: Cu, Pb, Cd...), Nhóm III (Thông số hữu cơ và dinh dưỡng: $\text{DO}, \text{BOD}_5, \text{COD}, \text{NH}_4^+, \text{NO}_3^-$), Nhóm IV (Độ đục, TSS), Nhóm V (Vi sinh - Coliform). Quy định công thức phi tuyến chặt chẽ cho $\text{DO}$ và phân tách ngưỡng nồng độ chính xác hơn, tránh hiện tượng che lấp rủi ro ô nhiễm độc học.

3. Nước thải chăn nuôi ảnh hưởng cụ thể như thế nào đến nguồn nước tưới tiêu của hệ thống?

Với $2.433,69 \text{ m}^3/\text{ngày đêm}$ nước thải chăn nuôi (80,28% chưa qua xử lý) cùng 8,84 tấn chất thải rắn/ngày, nguồn thải này đưa vào lòng kênh lượng cực lớn hợp chất hữu cơ, vi khuẩn Salmonella, Coliform, trứng giun sán và dư lượng kháng sinh. Điều này làm cạn kiệt oxy hòa tan ($\text{DO} < 1,5 \text{ mg/L}$), bùng phát hiện tượng phú dưỡng hóa, gây ngộ độc rễ cây trồng và lây lan dịch bệnh lở mồm long móng, cúm gia cầm.

4. Công nghệ BASTAF kết hợp bãi lọc ngầm trồng cây có ưu thế gì đối với khu vực nông thôn Vĩnh Phúc?

BASTAF kết hợp bãi lọc ngầm trồng sậy/cỏ Vetiver hoạt động hoàn toàn dựa trên quá trình lắng, lọc kỵ khí và hấp phụ sinh học tự nhiên của rễ cây, không cần tiêu tốn điện năng liên tục cho sục khí, không sử dụng hóa chất đắt tiền. Chi phí đầu tư chỉ bằng 30–40% và chi phí vận hành bảo dưỡng chỉ bằng 15–20% so với các công nghệ hiếu khí nhân tạo, rất phù hợp với năng lực tài chính và trình độ quản lý của cấp xã/thôn.

5. Biện pháp vận hành thủy lực nào có thể cải thiện tức thời chất lượng nước vào mùa kiệt?

Áp dụng quy trình xả nước đệm và pha loãng cưỡng bức: Mở cống đầu mối Đập Liễn Sơn kết hợp vận hành tối đa công suất Trạm bơm Bạch Hạc ($Q = 11,2 \text{ m}^3/\text{s}$) và Trạm bơm Đại Định ($Q = 11,2 \text{ m}^3/\text{s}$) vào đầu các đợt đổ ải vụ Đông Xuân. Nguồn nước sạch từ sông Hồng và sông Lô sẽ đẩy trôi dòng nước tù đọng ô nhiễm trên sông Phan ra cống Nam Viêm để tiêu thoát vào sông Cà Lỗ, làm tăng nồng độ $\text{DO}$ và hạ thấp tức thời hàm lượng $\text{COD}, \text{NH}_4^+$.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán cấp bách về bảo vệ môi trường nước tại Hệ thống Thủy lợi Liễn Sơn thông qua:

  1. Xác lập cơ sở dữ liệu quan trắc tin cậy qua 19 đợt khảo sát thực địa (2018–2020) với 15 thông số hóa lý, vi sinh và kim loại nặng trên 15 trạm trọng điểm.
  2. Định lượng chính xác bức tranh phát sinh chất thải $41.351,66 \text{ m}^3/\text{ngày đêm}$, chứng minh nguồn ô nhiễm chính đến từ nước thải sinh hoạt ($55,17%$) và công nghiệp ($22,94%$), phân vùng rõ 3 điểm nóng ô nhiễm trên sông Phan ($WQI < 25$).
  3. Ứng dụng thành công thuật toán WQI theo Quyết định 1460/QĐ-TCMT và hoạch định 3 kịch bản tải lượng đến năm 2025 gắn liền với quy hoạch phát triển bền vững của tỉnh Vĩnh Phúc.
  4. Đề xuất đồng bộ chuỗi giải pháp kỹ thuật công nghệ sinh thái (Mương oxy hóa, BASTAF bãi lọc ngầm) kết hợp quy trình điều tiết thủy lực liên hồ, cung cấp luận cứ khoa học thực tiễn cho công tác quản lý tài nguyên nước và phát triển nông nghiệp sinh thái hiện đại.