Giới thiệu dự án

Bối cảnh và thực trạng quản lý tài nguyên nước

Huyện Thủy Nguyên là trung tâm kinh tế công nghiệp, dịch vụ trọng điểm phía Bắc thành phố Hải Phòng với tốc độ đô thị hóa nhanh chóng (dự kiến trở thành thành phố trực thuộc Hải Phòng). Địa bàn có diện tích tự nhiên $24.279,9\text{ ha}$, dân số trên $341.399\text{ người}$, được bao bọc bởi 4 hệ thống sông lớn: sông Kinh Thầy, sông Cửa Cấm, sông Hàn Mấu và sông Bạch Đằng. Sự bùng nổ công nghiệp hóa – đô thị hóa cùng với biến đổi khí hậu đang tạo áp lực gay gắt lên hệ thống thủy lợi:

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

Hệ thống thủy lợi Thủy Nguyên quản lý $74\text{ km}$ đê (trên $42\text{ km}$ đê cấp 3), $79$ cống dưới đê, gần $200$ trạm bơm điện và $229$ tuyến kênh chính/cấp 1 dài trên $300\text{ km}$. Tuy nhiên, công tác vận hành đang gặp các thách thức nghiêm trọng:

  1. Xâm nhập mặn mùa kiệt: Độ mặn vượt ngưỡng $1,0‰$ lấn sâu $30 - 50\text{ km}$ từ biển vào sông Cửa Cấm và Bạch Đằng, vô hiệu hóa phần lớn các cống lấy nước tự chảy nhỏ ngoại trừ $02$ cống đầu mối An Sơn II và Phi Liệt.
  2. Áp lực ô nhiễm đa nguồn: Hơn $3.098$ cơ sở công nghiệp xả thải trực tiếp vào hệ thống kênh mương cùng với nước thải từ bãi rác Da Lợn (xã Minh Tân), nước thải sinh hoạt đô thị Núi Đèo - Minh Đức, nước thải làng nghề (đúc Mỹ Đồng, hương Kiền Bái) và dư lượng hóa chất bảo vệ thực vật từ nông nghiệp.
  3. Mâu thuẫn nhu cầu sử dụng nước: Nhu cầu cấp nước thô hàng năm đạt gần $200\text{ triệu m}^3$, tiêu thoát trên $230\text{ triệu m}^3$, trong khi cơ sở hạ tầng bị chia cắt do các dự án giao thông, đô thị mới.

Mục tiêu dự án

  1. Đánh giá toàn diện hiện trạng thủy văn, diễn biến xâm nhập mặn và chất lượng nguồn nước mặt theo bộ quy chuẩn QCVN 08-MT:2015/BTNMT (Cột B1).
  2. Phân tích cân bằng nước và xác định khả năng đáp ứng của $02$ công trình đầu mối chính (Cống An Sơn II và Cống Phi Liệt).
  3. Đề xuất khung kiến trúc công nghệ và giải pháp kỹ thuật - điều hành giám sát tài nguyên nước tự động phục vụ đa mục tiêu (nông nghiệp, công nghiệp, sinh hoạt).

Phạm vi và giới hạn (Scope & Limitations)

  • Không gian: Toàn bộ lưu vực thủy nông Thủy Nguyên thuộc tọa độ $20^\circ 52'$ – $21^\circ 01'$ Vĩ Bắc và $106^\circ 31'$ – $106^\circ 46'$ Kinh Đông.
  • Thời gian: Chuỗi số liệu khí tượng - thủy văn quan trắc giai đoạn 1958 – 2023 từ trạm Phù Liễn, Bến Bình, Hòn Dấu và mẫu phân tích chất lượng nước thực địa tháng 5/2023.
  • Giới hạn kỹ thuật: Nghiên cứu tập trung vào nguồn nước mặt hở; nước dưới đất chỉ xét ở mức độ hiện trạng sử dụng cục bộ.

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

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

Tiêu chí Giải pháp quản lý thủ công (Truyền thống) Giải pháp vận hành theo quy trình tĩnh Hệ thống giám sát & Điều hành tích hợp (Đề xuất)
Cơ chế đo độ mặn Lấy mẫu cơ học tại trạm theo ca Đo định kỳ 2 lần/ngày Cảm biến IoT EC/Độ mặn thời gian thực trực tuyến
Độ trễ vận hành cống $2 - 4\text{ giờ}$ sau đỉnh triều Theo lịch con nước cố định Đóng/mở tự động kích hoạt bởi ngưỡng mặn $< 1,0‰$
Giám sát xả thải Thanh tra định kỳ tại chỗ Báo cáo kiểm toán môi trường năm Bản đồ GIS Web tích hợp dữ liệu $3.098$ điểm xả thải
Khả năng ứng phó sự cố Bị động, dễ nhiễm mặn lòng kênh Phản ứng sau khi nước ngọt bị nhiễm Cảnh báo sớm tự động qua mô hình cân bằng dòng chảy

Phân loại yêu cầu hệ thống theo mô hình MoSCoW:

  • Must have: Tự động hóa kiểm soát đóng mở cống An Sơn II và Phi Liệt khi độ mặn $> 1,0‰$; module tính toán cân bằng nước theo lưu lượng thực tế $Q_{\text{Bến Bình}} = 779\text{ m}^3/\text{s}$.
  • Should have: Tích hợp module GIS giám sát $3.098$ nguồn xả thải công nghiệp; hệ thống cảnh báo ô nhiễm vượt chuẩn QCVN 08-MT:2015/BTNMT.
  • Could have: Tích hợp thuật toán dự báo xâm nhập mặn theo triều Hòn Dấu bằng Machine Learning (LSTM).
  • Won't have: Điều khiển tự động cấp nước đến tận bể lắng hộ gia đình trong giai đoạn này.

Thiết kế hệ thống

Ngăn xếp công nghệ (Technology Stack)

  • IoT & SCADA Edge: Modbus RTU over RS-485, ESP32/PLC S7-1200 firmware v4.5, giao thức truyền tin MQTT v5.0.
  • Backend & Database: Python 3.11, FastAPI 0.104.1, PostgreSQL 15 kết hợp extension PostGIS 3.3 lưu trữ không gian mạng lưới kênh mương.
  • Engine tính toán thủy lực: Tích hợp mô hình toán thủy văn trên NumPy 1.26 và GeoPandas 0.14.
  • Frontend Dashboard: React 18.2, TypeScript 5.2, TailwindCSS 3.3, Leaflet 1.9.4 thể hiện bản đồ dòng chảy và vị trí quan trắc.

Phương pháp nghiên cứu và quản lý (Methodology)

Nghiên cứu áp dụng quy trình đánh giá hệ số phát thải kết hợp mô hình cân bằng nước theo phương trình tổng quát:

$$W_{\text{dự trữ}}(t) = W(t-1) + \left( Q_{\text{An Sơn II}} + Q_{\text{Phi Liệt}} + P_{\text{mưa}} \cdot A_{\text{lưu vực}} \right) - \left( E_{\text{bốc hơi}} + Q_{\text{tưới}} + Q_{\text{sinh hoạt}} + Q_{\text{công nghiệp}} + Q_{\text{thau chua}} \right)$$

Trong đó tốc độ bốc hơi mặt nước được hiệu chỉnh theo công thức ống Piche:

$$E_{\text{nước}} = K \cdot E_{\text{Piche}}$$

Với hệ số $K$ xác định dựa trên chuỗi số liệu 1958–2019 trạm Phù Liễn ($E_{\text{năm}} = 713,7\text{ mm}$).

Implementation và kết quả

Quy trình phân tích dữ liệu và thuật toán điều khiển

Hệ thống triển khai thuật toán kiểm soát đóng/mở công trình lấy nước tự động nhằm ngăn chặn hiện tượng "đầu mặn, cuối chua" tại khu vực Đông Nam huyện:

"""
Module: water_gate_controller.py
Chức năng: Thuật toán quyết định vận hành cống An Sơn II và Phi Liệt 
dựa trên quan trắc độ mặn, mực nước triều và ngưỡng tiêu chuẩn QCVN 08-MT:2015.
"""
from dataclasses import dataclass
from datetime import datetime
from enum import Enum

class GateState(Enum):
    OPEN = "MO_CONG_LAY_NGOT"
    CLOSE = "DONG_CONG_NGAN_MAN"
    DRAIN = "MO_CONG_XUAT_TIEU"

@dataclass
class HydrologicalData:
    timestamp: datetime
    gate_name: str
    salinity_permil: float      # Độ mặn (‰)
    water_level_outer: float     # Mực nước ngoài sông (m)
    water_level_inner: float     # Mực nước trong đồng (m)
    do_mg_l: float               # Nồng độ oxy hòa tan (mg/L)
    cod_mg_l: float              # COD (mg/L)

class WaterQualityDecisionEngine:
    SALINITY_THRESHOLD_PERMIL: float = 1.0   # Ngưỡng cho phép tưới nông nghiệp
    MIN_DO_THRESHOLD: float = 4.0            # QCVN 08 Cột B1 (DO >= 4 mg/L)
    MAX_COD_THRESHOLD: float = 30.0          # QCVN 08 Cột B1 (COD <= 30 mg/L)

    def evaluate_gate_action(self, data: HydrologicalData) -> dict:
        # 1. Kiểm tra an toàn độ mặn
        if data.salinity_permil >= self.SALINITY_THRESHOLD_PERMIL:
            return {
                "decision": GateState.CLOSE.value,
                "reason": f"Độ mặn {data.salinity_permil}‰ vượt ngưỡng an toàn ({self.SALINITY_THRESHOLD_PERMIL}‰)",
                "alert_level": "HIGH"
            }
        
        # 2. Kiểm tra chất lượng nước đầu vào
        if data.do_mg_l < self.MIN_DO_THRESHOLD or data.cod_mg_l > self.MAX_COD_THRESHOLD:
            return {
                "decision": GateState.CLOSE.value,
                "reason": "Chất lượng nước sông bị ô nhiễm hữu cơ vượt ngưỡng",
                "alert_level": "CRITICAL"
            }

        # 3. Đánh giá chênh lệch mực nước để quyết định lấy nước hoặc tiêu thoát
        if data.water_level_outer > data.water_level_inner + 0.15:
            return {
                "decision": GateState.OPEN.value,
                "reason": "Thỏa mãn điều kiện nước ngọt tự chảy vào hệ thống",
                "alert_level": "NORMAL"
            }
        elif data.water_level_inner > data.water_level_outer + 0.30:
            return {
                "decision": GateState.DRAIN.value,
                "reason": "Chênh lệch áp lực phù hợp mở cống tiêu úng nội đồng",
                "alert_level": "NORMAL"
            }
        
        return {
            "decision": GateState.CLOSE.value,
            "reason": "Chênh lệch mực nước không đủ điều kiện thủy lực",
            "alert_level": "LOW"
        }

-- Schema cơ sở dữ liệu PostGIS phục vụ giám sát nguồn xả thải huyện Thủy Nguyên
CREATE EXTENSION IF NOT EXISTS postgis;

CREATE TABLE water_monitoring_station (
    station_id VARCHAR(20) PRIMARY KEY,
    station_name VARCHAR(100) NOT NULL,
    waterbody_name VARCHAR(50), -- Sông Kinh Thầy, Bạch Đằng, Giá, Hòn Ngọc
    location GEOMETRY(Point, 4326),
    is_active BOOLEAN DEFAULT TRUE
);

CREATE TABLE discharge_sources (
    source_id VARCHAR(30) PRIMARY KEY,
    enterprise_name VARCHAR(255) NOT NULL,
    industry_type VARCHAR(100), -- Đúc kim loại, xi măng, dệt may, chăn nuôi
    daily_discharge_volume_m3 NUMERIC(10, 2),
    location GEOMETRY(Point, 4326),
    receives_treatment BOOLEAN DEFAULT FALSE,
    outfall_waterway VARCHAR(100)
);

CREATE TABLE water_quality_log (
    log_id BIGSERIAL PRIMARY KEY,
    station_id VARCHAR(20) REFERENCES water_monitoring_station(station_id),
    recorded_at TIMESTAMP WITH TIME ZONE DEFAULT CURRENT_TIMESTAMP,
    temperature_c NUMERIC(4, 2),
    ph NUMERIC(3, 2),
    do_mg_l NUMERIC(4, 2),
    tss_mg_l NUMERIC(6, 2),
    cod_mg_l NUMERIC(6, 2),
    ammonia_n_mg_l NUMERIC(4, 3),
    chloride_mg_l NUMERIC(6, 2),
    coliform_mpn NUMERIC(8, 2)
);

Kết quả thử nghiệm và kiểm chứng chất lượng nước

Kết quả quan trắc thực địa tại $04$ điểm nhạy cảm đầu mối (ngày 14/05/2023) được đối chiếu theo QCVN 08-MT:2015/BTNMT (Cột B1):

Thông số Đơn vị N1 (Thượng lưu An Sơn II) N2 (Hạ lưu An Sơn II) N3 (Thượng lưu Phi Liệt) N4 (Hạ lưu Phi Liệt) QCVN 08 Cột B1
Nhiệt độ $^\circ\text{C}$ 24,7 25,1 24,8 25,2 —
pH — 7,35 7,24 7,19 7,22 $5,5 - 9,0$
DO $\text{mg/L}$ 5,98 6,13 6,09 6,11 $\ge 4,0$
TSS $\text{mg/L}$ 20,0 14,0 18,0 13,0 $\le 50$
COD $\text{mg/L}$ 15,0 10,0 12,0 11,0 $\le 30$
Amoni ($NH_4^+ - N$) $\text{mg/L}$ 0,10 0,11 0,11 0,09 $\le 0,9$
Clorua ($Cl^-$) $\text{mg/L}$ 3,3 3,1 3,2 2,8 $\le 350$
Sắt tổng số ($Fe$) $\text{mg/L}$ 0,07 0,08 0,07 0,09 $\le 1,5$
Coliform $\text{VK/100ml}$ 2.500 3.400 3.100 2.900 $\le 7.500$

Đánh giá lưu lượng và phân phối dòng chảy

  • Lưu lượng bình quân nhiều năm tại trạm thủy văn Bến Bình (thượng lưu Kinh Thầy) đạt $Q_{\text{TB}} = 779\text{ m}^3/\text{s}$.
  • Mùa lũ (tháng VI - X) chiếm $71%$ tổng dòng chảy ($8,5\text{ tỷ m}^3$), dòng chảy tháng VIII đạt đỉnh chiếm $20%$.
  • Mùa cạn (tháng XI - V năm sau) chỉ chiếm $29%$ ($3,5\text{ tỷ m}^3$), thấp nhất vào tháng I ($221\text{ m}^3/\text{s}$) và tháng III ($3%$).
  • Việc điều hành lấy nước qua cống An Sơn II và Phi Liệt đã cấp ổn định $200\text{ triệu m}^3/\text{năm}$, bảo đảm tưới tiêu cho hơn $4.000\text{ ha}$ diện tích sản xuất nông nghiệp và cung cấp $15,7\text{ triệu m}^3/\text{năm}$ nước sinh hoạt thô từ hồ Sông Giá.

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

Các điểm cải tiến kỹ thuật chính

  1. Chuẩn hóa quy trình vận hành theo ngưỡng triều động: Thay vì đóng mở cống cố định theo giờ, hệ thống đề xuất điều tiết mở cống đón đỉnh triều lấy ngọt chỉ khi độ mặn $< 1,0‰$ (tại sông Kinh Thầy từ trạm Cao Kênh trở lên và sông Bạch Đằng từ trạm Đồn Sơn trở lên).
  2. Số hóa mạng lưới $3.098$ điểm nguồn thải: Thiết lập cơ sở dữ liệu GIS tích hợp toàn bộ các nguồn xả thải công nghiệp, nông thôn và làng nghề vào bản đồ quản lý không gian mạng lưới kênh mương.
  3. Cơ chế rửa mặn - thau chua đa tầng: Tối ưu hóa chu trình dẫn nước từ cống An Sơn II qua hệ thống sông Sau - sông Si - kênh Hòn Ngọc để bổ cập cho sông Giá, giúp đẩy mặn và giảm tải ô nhiễm hữu cơ nội đồng, tiêu thoát an toàn $> 230\text{ triệu m}^3/\text{năm}$.

Đóng góp cho ngành quản lý môi trường và tài nguyên nước

  • Mặt khoa học: Cung cấp bức tranh thực trạng thủy hóa, thủy triều và động lực học dòng chảy tại khu vực bán sơn địa chuyển tiếp châu thổ sông Hồng.
  • Mặt thực tiễn: Cung cấp cơ sở khoa học trực tiếp giúp Công ty TNHH MTV Khai thác công trình thủy lợi Thủy Nguyên quản lý an toàn nguồn vốn đầu tư nâng cấp công trình (trên $165\text{ tỷ đồng}$ trong 10 năm qua cho 31 cống đập và 11 cầu thông kênh).

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

Kịch bản vận hành thực địa (Use Cases)

Yêu cầu triển khai hệ thống (Deployment Requirements)

  • Hạ tầng phần cứng:
    • $04$ trạm quan trắc tự động chuyên dụng (đo liên tục: pH, DO, Độ mặn/EC, TSS, Mực nước) tại cống An Sơn II, Phi Liệt, Bính Động và trạm Đồn Sơn.
    • Máy chủ Server tối thiểu $8\text{ Core CPU}$, $32\text{ GB RAM}$, lưu trữ $1\text{ TB SSD RAID-1}$ đặt tại trung tâm điều hành.
  • Môi trường phần mềm:
    • Hệ điều hành: Ubuntu Server 22.04 LTS.
    • Cơ sở dữ liệu: PostgreSQL 15 + PostGIS 3.3.
    • Nền tảng GIS Server: GeoServer 2.24 phục vụ xuất bản bản đồ WMS/WFS.

Hiệu quả kinh tế - xã hội và lộ trình (Roadmap)

  • Hiệu quả kinh tế: Tiết kiệm khoảng $15 - 20%$ điện năng vận hành của gần $200$ trạm bơm nhờ tận dụng tối đa chu kỳ triều tự chảy; giảm thiểu thiệt hại do mặn gây chết cây trồng trên $4.000\text{ ha}$ đất canh tác.
  • Lộ trình 3 giai đoạn:
    • Giai đoạn 1 (06 tháng): Lắp đặt mạng lưới cảm biến và số hóa dữ liệu $3.098$ nguồn xả thải.
    • Giai đoạn 2 (12 tháng): Triển khai phần mềm Web GIS và kết nối tự động với hệ thống điều khiển cống.
    • Giai đoạn 3 (18 tháng): Tích hợp mô hình dự báo ngập úng và xâm nhập mặn toàn huyện thời gian thực.

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

Hạn chế kỹ thuật

  • Chuỗi số liệu phân tích chất lượng nước chuyên sâu mới tập trung vào đợt đo tháng 5/2023, cần mở rộng đo đạc liên tục qua nhiều mùa mưa/mùa cạn để đánh giá xu thế biến đổi dài hạn.
  • Cơ sở hạ tầng kênh mương cấp 2, cấp 3 tại một số xã vùng núi phía Bắc còn phân tán, chưa được kiên cố hóa hoàn chỉnh gây khó khăn cho việc kiểm soát thất thoát.

Hướng nâng cấp đề xuất

  • Ứng dụng công nghệ ảnh viễn thám (Sentinel-2) để tự động nhận diện vùng nước mặt bị bèo rác xâm lấn làm tắc nghẽn dòng chảy.
  • Tích hợp mô hình trí tuệ nhân tạo (AI/Deep Learning) dự báo nồng độ ô nhiễm và lưu lượng dòng chảy Bến Bình trước 48 giờ dựa trên dữ liệu khí tượng vệ tinh.

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

  • Công ty Khai thác Công trình Thủy lợi: Sở hữu bộ công cụ tự động hóa quyết định vận hành, bảo vệ an toàn công trình trị giá hàng trăm tỷ đồng.
  • Người dân và Doanh nghiệp trên địa bàn: Đảm bảo an ninh cấp nước cho sản xuất công nghiệp và đời sống dân sinh, loại trừ nguy cơ nước máy bị nhiễm mặn.
  • Cộng đồng Kỹ sư & Nghiên cứu Môi trường: Tiếp cận phương pháp luận tích hợp mô hình thủy lực, số liệu khí tượng - thủy hải văn thực tế và giải pháp IoT GIS trong bảo vệ môi trường nước.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai hệ thống quản lý nước thông minh tại Thủy Nguyên là gì?

Cần trang bị các trạm đo độ mặn và mực nước tự động tại các vị trí xung yếu (An Sơn II, Phi Liệt, trạm Cao Kênh, Đồn Sơn) hỗ trợ giao thức Modbus/RS485 hoặc 4G/NB-IoT, truyền dữ liệu thời gian thực về máy chủ trung tâm chạy PostgreSQL/PostGIS.

2. Giới hạn độ mặn cho phép lấy nước phục vụ nông nghiệp là bao nhiêu?

Theo quy chuẩn kỹ thuật và thực tiễn vận hành hệ thống thủy lợi Thủy Nguyên, ngưỡng giới hạn độ mặn an toàn tối đa để lấy nước tưới nông nghiệp là dưới $1,0‰$. Khi độ mặn ngoài sông vượt quá $1,0‰$, các cống đầu mối phải đóng hoàn toàn để ngăn mặn xâm nhập nội đồng.

3. Hệ thống xử lý thế nào khi xảy ra mâu thuẫn giữa nhu cầu tưới nông nghiệp và cấp nước thô cho công nghiệp/sinh hoạt?

Thuật toán ưu tiên hàng đầu cho việc tích trữ nước ngọt tại dung tích hồ chứa sông Giá (nguồn cấp sinh hoạt $15,7\text{ triệu m}^3/\text{năm}$). Sau đó, nguồn nước từ kênh Hòn Ngọc và sông Si sẽ được điều tiết phân phối cho các vùng sản xuất nông nghiệp thông qua việc đóng mở luân phiên hệ thống cống nội đồng.

4. Chi phí đầu tư và thời gian thu hồi vốn (ROI) của giải pháp?

Ước tính chi phí đầu tư hạ tầng quan trắc tự động và phần mềm điều hành khoảng $1,5 - 2,5\text{ tỷ đồng}$. Thời gian hoàn vốn đạt từ $2 - 3\text{ năm}$ thông qua việc tiết kiệm điện năng bơm tưới, giảm thiểu chi phí nạo vét xử lý sự cố ô nhiễm và bảo vệ năng suất mùa vụ nông nghiệp.

5. Khả năng tương thích và tích hợp với hệ thống đô thị thông minh Hải Phòng?

Hệ thống được thiết kế theo kiến trúc Microservices với chuẩn giao tiếp RESTful API / GeoJSON, cho phép tích hợp trực tiếp vào Trung tâm Giám sát Điều hành Đô thị thông minh (IOC) thành phố Hải Phòng và thành phố Thủy Nguyên trong tương lai.

Kết luận

Đồ án đã phân tích sâu sắc thực trạng hệ thống thủy nông Thủy Nguyên – công trình có quy mô $74\text{ km}$ đê, $79$ cống, phục vụ cấp nước $200\text{ triệu m}^3$ và tiêu úng $230\text{ triệu m}^3$ hàng năm. Kết quả nghiên cứu làm sáng tỏ bức tranh thủy lực, cơ chế xâm nhập mặn mùa kiệt và phân tích chất lượng nước đạt chuẩn QCVN 08-MT:2015/BTNMT (Cột B1) tại $02$ công trình đầu mối huyết mạch cống An Sơn II và Phi Liệt.

Việc tích hợp giải pháp quan trắc thông minh, cơ sở dữ liệu GIS giám sát $3.098$ nguồn thải và quy trình vận hành tự động theo ngưỡng độ mặn $1,0‰$ là bước chuyển dịch tất yếu nhằm đảm bảo an ninh nguồn nước, phục vụ đắc lực cho tiến trình đưa huyện Thủy Nguyên trở thành thành phố hiện đại, bền vững.