Giới thiệu dự án

Khu vực vịnh Hạ Long (tỉnh Quảng Ninh) là Di sản Thiên nhiên Thế giới hai lần được UNESCO công nhận với diện tích vùng lõi 434 $\text{km}^2$. Nằm tại trung tâm dịch vụ, du lịch và cảng biển của thành phố Hạ Long, phường Bãi Cháy đóng vai trò đầu mối giao thương và lưu trú ven biển trọng điểm. Theo thống kê của Sở Văn hóa, Thể thao và Du lịch Quảng Ninh, năm 2013 khu vực vịnh Hạ Long đón 7,5 triệu lượt khách du lịch (trong đó có 2,5 triệu khách quốc tế). Đi cùng với tốc độ phát triển kinh tế biển vượt bậc là áp lực ô nhiễm môi trường nước cục bộ nghiêm trọng.

+-----------------------------------------------------------------------------------+
|                            ÁP LỰC MÔI TRƯỜNG VEN BIỂN                             |
|                                                                                   |
|  [Đô thị hóa & San lấp]     [Hoạt động Cảng & Tàu thuyền]   [Nước thải Sinh hoạt] |
|   - Mất 295 ha rừng           - 1.000 lượt tàu Cái Lân/ngày  - 12 triệu m3/năm    |
|   - Tăng độ đục, bùn lắng     - 500 tàu du lịch Bãi Cháy     - Trạm XLNT quá tải  |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
|              NGHIÊN CỨU & MÔ HÌNH HÓA CHẤT LƯỢNG NƯỚC BIỂN VEN BỜ                |
|           (Tích hợp Thử nghiệm VILAS + ArcGIS 10.1 Kriging + PRIMER 5.0)          |
+-----------------------------------------------------------------------------------+

Vấn đề thực tế (Problem Statement)

Sự gia tăng nhanh chóng của các hoạt động nhân sinh ven bờ đã tạo ra các điểm nóng ô nhiễm môi trường nước biển tầng mặt tại Bãi Cháy với ba điểm nghẽn chính:

  • Áp lực hạ tầng xử lý nước thải: Lượng nước thải sinh hoạt đô thị phát sinh tại Hạ Long lên tới 12 triệu $\text{m}^3/\text{năm}$. Nhà máy xử lý nước thải (XLNT) Bãi Cháy thường xuyên vận hành quá tải gấp 1,8 - 2 lần công suất thiết kế ($3.500\text{ m}^3/\text{ngày đêm}$ so với thực nhận đỉnh điểm đạt $6.500\text{ m}^3/\text{ngày đêm}$), dẫn đến hiện tượng tràn nước thải chưa xử lý qua cửa ngăn triều ra thẳng bãi tắm.
  • Ô nhiễm dầu khoáng và bùn nạo vét: Mật độ lưu thông hơn 1.000 lượt tàu hàng/ngày tại Cụm cảng Cái Lân và hơn 500 tàu du lịch (với 150 tàu lưu trú qua đêm) tại Bến tàu du lịch Bãi Cháy xả thải lượng lớn nước đáy tàu (nước la-canh chứa dầu) và rác thải hữu cơ (ước tính $2\text{ kg chất thải/tàu/ngày}$). Hiện tượng đổ trộm bùn nạo vét luồng lạch với khối lượng hàng trăm $\text{m}^3$ (điển hình vụ bắt giữ hai tàu xả $300\text{ m}^3$ bùn năm 2013) làm gia tăng hàm lượng chất rắn lơ lửng.
  • Suy giảm hệ sinh thái đệm: Hoạt động san nền lấn biển giai đoạn 1998 - 2003 làm biến mất 295 ha rừng ngập mặn tại thành phố Hạ Long (trong đó 134 ha do phát triển đô thị), làm triệt tiêu khả năng tự làm sạch tự nhiên của vùng triều (mất năng lực giảm $409,5\text{ kg TSS/ha/ngày}$, $0,234\text{ kg N-T/ha/ngày}$ và $0,117\text{ kg P-T/ha/ngày}$).

Mục tiêu của dự án

  1. Đánh giá toàn diện hiện trạng chất lượng môi trường nước biển ven bờ thông qua hệ thống chỉ tiêu lý - hóa - sinh tại 5 điểm quan trắc đại diện khu vực Bãi Cháy.
  2. Xác định cấu trúc tương quan giữa các nguồn thải nhân sinh và các thông số suy thoái môi trường bằng các thuật toán thống kê sinh thái đa biến (MDS, PCA).
  3. Ứng dụng kỹ thuật nội suy không gian Kriging trên nền tảng ArcGIS để trực quan hóa bản đồ phân bố và lan truyền chất ô nhiễm (đặc biệt là dầu mỡ khoáng và amoni), từ đó thiết lập hệ giải pháp quản trị môi trường khả thi.

Phương pháp tiếp cận giải pháp

Đề tài kết hợp phương pháp quan trắc hóa nghiệm chuẩn VILAS theo các Tiêu chuẩn Việt Nam (TCVN) và Quy chuẩn kỹ thuật Quốc gia (QCVN 10:2008/BTNMT), kết hợp với kỹ thuật phân tích số liệu đa biến trên phần mềm PRIMER v5.0 (Plymouth Routines in Multivariate Ecological Research) và mô hình hóa lan truyền ô nhiễm trên phần mềm ArcGIS 10.1 của Esri.

Kết quả kỳ vọng và chỉ số đo lường

  • Định lượng chính xác nồng độ 11 thông số then chốt: pH, Độ mặn, Độ đục, Dầu mỡ khoáng, Tổng chất rắn lơ lửng (TSS - Total Suspended Solids), Amoni ($\text{NH}_4^+$), Sắt (Fe), Mangan (Mn), Nhu cầu oxy sinh hóa sau 5 ngày ($\text{BOD}_5$), Oxy hòa tan (DO - Dissolved Oxygen) và Coliform.
  • Xây dựng mô hình phân tích độ tương đồng cụm (Cluster Analysis) và không gian đa chiều (MDS/PCA) đạt độ tương thích $\text{Stress} < 0,1$ và tỷ lệ tương đồng sinh thái (Similarity) từ 83% đến 100%.
  • Cung cấp tập bản đồ số nội suy nồng độ dầu mỡ với lưới ô vuông độ phân giải không gian cao hỗ trợ cảnh báo sớm.

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

Nghiên cứu tập trung tại vùng nước ven bờ cách mép nước triều kiệt 500 m trở vào thuộc phường Bãi Cháy (từ Bến tàu du lịch Bãi Cháy, Bãi tắm trung tâm đến vùng nước Cầu Bãi Cháy và Cụm cảng Cái Lân). Dữ liệu thu thập và thực địa triển khai từ tháng 01/2014 đến tháng 05/2014 dựa trên chuỗi số liệu quan trắc năm 2013 và dữ liệu liên năm của Chi cục Bảo vệ Môi trường tỉnh Quảng Ninh.

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

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

Trước đây, công tác đánh giá ô nhiễm nước biển tại địa phương chủ yếu dựa trên việc so sánh bảng biểu thủ công với ngưỡng giới hạn cho phép (GHCP). Phương pháp này bộc lộ nhiều hạn chế trong việc đánh giá tổng thể mối tương quan đa biến và dự báo xu hướng lan truyền không gian.

Tiêu chí so sánh Phương pháp truyền thống (Báo cáo cơ học) Phương pháp tiếp cận của đề tài (GIS + PRIMER)
Xử lý dữ liệu So sánh giá trị đơn lẻ với quy chuẩn (QCVN) Phân tích cấu trúc dữ liệu đa biến (PCA, MDS, Bray-Curtis)
Trực quan hóa Bảng số liệu biểu đồ 2D dạng cột rời rạc Bản đồ phân bố nồng độ liên tục theo mô hình không gian Kriging
Xác định nguồn thải Định tính dựa trên vị trí gần cống thải Định lượng mối liên kết nguồn thải dựa trên mức độ đồng dạng
Hỗ trợ quyết định Phản ứng thụ động sau khi phát sinh ô nhiễm Khoanh vùng cảnh báo các vệt ô nhiễm dầu lan truyền 
+------------------------------------------------------------------------------------+
|                      MA TRẬN ƯU TIÊN YÊU CẦU NGHIÊN CỨU (MoSCoW)                  |
+------------------------------------------------------------------------------------+
|  MUST HAVE (Bắt buộc)           |  SHOULD HAVE (Nên có)                            |
|  - Chuẩn hóa quy trình quan trắc|  - Phân tích tương quan đa biến MDS & PCA        |
|    theo TCVN/ISO VILAS.         |    trên PRIMER v5.0.                             |
|  - So sánh đối soát với         |  - Lập bản đồ nội suy Kriging trên ArcGIS 10.1   |
|    QCVN 10:2008/BTNMT.          |    cho thông số dầu mỡ và DO.                    |
|  - Định lượng nồng độ Dầu khoáng|  - Đánh giá năng lực thu gom của trạm XLNT.      |
|    và Oxy hòa tan (DO).         |                                                  |
+---------------------------------+--------------------------------------------------+
|  COULD HAVE (Có thể có)         |  WON'T HAVE (Chưa thực hiện đợt này)             |
|  - Tích hợp số liệu khí tượng   |  - Mô phỏng động lực học chất lỏng thủy triều    |
|    (gió mùa, dòng chảy triều).  |    3D thời gian thực (MIKE 21 / Delft3D).        |
|  - Khảo sát xã hội học 100 mẫu. |  - Trạm quan trắc IoT cảm biến trực tuyến 24/7.  |
+------------------------------------------------------------------------------------+

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu của đề tài được tổ chức thành 4 tầng logic khép kín từ thu thập mẫu thực địa đến mô hình hóa không gian:

graph TD
    subgraph Layer1["1. Thu thập & Đo đạc Hiện trường"]
        A1["Quan trắc Hiện trường (pH, DO, Nhiệt độ)"] --> A3["Bảo quản Mẫu theo TCVN 5998:1995"]
        A2["Lấy mẫu Hóa lý - Vi sinh (D1 - D5)"] --> A3
    end

    subgraph Layer2["2. Phân tích Phòng Thí nghiệm VILAS"]
        A3 --> B1["Trắc phổ Hấp thụ AAS / UV-Vis"]
        A3 --> B2["Chuẩn độ Winkler (DO)"]
        A3 --> B3["Màng lọc Vi sinh (Coliform)"]
        A3 --> B4["Chiết tách Dầu khoáng"]
    end

    subgraph Layer3["3. Xử lý Thống kê & Phân tích Đa biến"]
        B1 & B2 & B3 & B4 --> C1["Chuẩn hóa Ma trận Dữ liệu (Excel)"]
        C1 --> C2["PRIMER 5.0 (Bray-Curtis Similarity)"]
        C2 --> C3["Phân tích Không gian Đa chiều MDS & PCA"]
    end

    subgraph Layer4["4. Mô hình hóa Không gian & Bản đồ"]
        C1 --> D1["Geodatabase (ArcGIS 10.1)"]
        D1 --> D2["Semivariogram Modeling"]
        D2 --> D3["Nội suy Kriging Lan truyền Ô nhiễm"]
        D3 --> D4["Bản đồ Phân vùng Chất lượng Nước"]
    end

Technology Stack và Tiêu chuẩn Kỹ thuật

  • Nền tảng GIS & Không gian: Esri ArcGIS Desktop v10.1 (Extension: Spatial Analyst, Geostatistical Analyst).
  • Phần mềm Thống kê Sinh thái: PRIMER-E v5.0 (Plymouth Routines in Multivariate Ecological Research).
  • Tiêu chuẩn Phân tích Phòng Thí nghiệm:
    • Oxy hòa tan (DO): TCVN 5499:1995 (Phương pháp chuẩn độ Winkler).
    • Chất rắn lơ lửng (TSS): TCVN 6625:2000 (Lọc qua màng sợi thủy tinh $0,45\ \mu\text{m}$).
    • Nhu cầu oxy sinh hóa: TCVN 6001:1995 ($\text{BOD}_5$ ủ ở $20^\circ\text{C}$).
    • Nhu cầu oxy hóa học: TCVN 6491:1999 ($\text{COD}$ bằng phương pháp Kali Dicromat).
    • Amoni ($\text{NH}_4^+$): TCVN 5988:1995 (Chưng cất và chuẩn độ).
    • Dầu mỡ khoáng: TCVN 5070:1995 (Phương pháp chiết dung môi và khối lượng/hồng ngoại).
    • Coliform: TCVN 6187-1:1996 (Phương pháp màng lọc nuôi cấy).

Thiết kế Cấu trúc Dữ liệu Không gian (Geodatabase Schema)

Table: Monitoring_Stations_WaterQuality

Methodology

Quy trình nghiên cứu áp dụng chu trình khép kín 4 giai đoạn chuẩn hóa trong khoa học quản lý môi trường:

+-----------------------------------------------------------------------------------+
|                            LỘ TRÌNH TRIỂN KHAI DỰ ÁN                              |
|                                                                                   |
|  [Tháng 01/2014]      Thu thập số liệu thứ cấp, bản đồ nền, hiện trạng KTXH.      |
|  [Tháng 02 - 03/2014] Khảo sát thực địa, lấy mẫu 5 trạm quan trắc, điều tra 100   |
|  [Tháng 04/2014]      Phân tích phòng thí nghiệm VILAS; xử lý số liệu trên        |
|  [Tháng 05/2014]      Tổng hợp kết quả, đánh giá sức tải, xây dựng nhóm giải pháp |
|                       và hoàn thiện báo cáo khóa luận.                            |
+-----------------------------------------------------------------------------------+

Ma trận Đánh giá Rủi ro Kỹ thuật và Biện pháp Giảm thiểu

Rủi ro kỹ thuật Mức độ Biện pháp giảm thiểu & Kiểm soát chất lượng
Mẫu nước bị biến tính trong quá trình vận chuyển Cao Bảo quản mẫu ở nhiệt độ $4^\circ\text{C}$ trong thùng lạnh, cố định mẫu DO tại hiện trường bằng $\text{MnSO}_4$ và hỗn hợp $\text{KI} + \text{NaN}_3$, chuyển về phòng thí nghiệm trong 24h.
Mạng lưới điểm quan trắc thưa (5 trạm) Trung bình Sử dụng mô hình Kriging Ordinary kết hợp với phân tích bán phương sai (Semivariogram) để tối ưu hóa hàm phân bổ không gian.
Nhiễu số liệu từ hoạt động xả thải tức thời Trung bình Thu mẫu đồng thời theo chu kỳ triều kiệt và triều cường để loại trừ yếu tố bất thường của chu kỳ thủy động lực.

Implementation và kết quả

Development process

Quá trình phân tích dữ liệu không gian sử dụng thuật toán nội suy thống kê địa lý Ordinary Kriging để ngoại suy nồng độ các chất ô nhiễm từ các trạm quan trắc rời rạc về bề mặt raster liên tục.

Phương trình toán học của mô hình Kriging: $$\hat{Z}(s_0) = \sum_{i=1}^{N} \lambda_i Z(s_i)$$

Trong đó:

  • $Z(s_i)$ là giá trị đo thực tế tại trạm quan trắc thứ $i$.
  • $\lambda_i$ là trọng số tối ưu gán cho giá trị tại trạm $i$, thỏa mãn điều kiện không chệch: $\sum_{i=1}^{N} \lambda_i = 1$.
  • Bán phương sai thực nghiệm (Experimental Semivariogram) được xác định bởi: $$\gamma(h) = \frac{1}{2N(h)}\sum_{i=1}^{N(h)} \left[ Z(s_i) - Z(s_i+h) \right]^2$$

Dưới đây là mã nguồn Python mẫu (tương thích Python 2.7/3.x sử dụng arcpy và thư viện thống kê) biểu diễn quy trình tự động hóa phân tích nội suy nồng độ Dầu mỡ từ tập dữ liệu quan trắc:

# -*- coding: utf-8 -*-
"""
Automation Script: Coastal Water Quality Spatial Interpolation (ArcGIS / ArcPy)
Description: Thuc hien noi suy Ordinary Kriging cho chi tieu Dau mo vung ven bien Bai Chay.
"""
import arcpy
from arcpy.sa import *

# 1. Thiet lap moi truong lam viec
arcpy.env.workspace = "C:/WaterQuality_BaiChay/Geodatabase.gdb"
arcpy.env.overwriteOutput = True
arcpy.CheckOutExtension("Spatial")

# 2. Cac tham so dau vao
input_features = "Monitoring_Stations_2013"
z_field = "Oil_mgL"
output_raster = "C:/WaterQuality_BaiChay/Output/Kriging_Oil_2013.tif"
mask_boundary = "BaiChay_Coastal_500m_Buffer.shp"

# 3. Mo hinh hoa Semivariogram va chay Kriging
# Model: Spherical Semivariogram voi lag size phu hop vung ven bo
semi_variogram_model = "SPHERICAL"
cell_size = 10.0  # Do phan giai luoi 10m
search_radius = RadiusVariable(5, 5000) # Toi thieu 5 diem trong ban kinh 5000m

try:
    print("[INFO] Dang thuc hien noi suy Ordinary Kriging cho thong so: {0}...".format(z_field))
    kriging_output = Kriging(input_features, z_field, semi_variogram_model, cell_size, search_radius)
    
    # Cat raster theo ranh gioi vung nghien cuu ven bo
    print("[INFO] Cat lop Raster theo ranh gioi dải ven bien 500m...")
    clipped_raster = ExtractByMask(kriging_output, mask_boundary)
    clipped_raster.save(output_raster)
    
    print("[SUCCESS] Mo hinh phan bo lan truyen Dau mo da duoc xuat thanh cong tai: {0}".format(output_raster))

except Exception as e:
    print("[ERROR] Phat sinh loi trong qua trinh noi suy: {0}".format(str(e)))
finally:
    arcpy.CheckInExtension("Spatial")

Testing và validation

Số liệu quan trắc năm 2013 tại 5 trạm quan trắc ven bờ phường Bãi Cháy được đối soát nghiêm ngặt với Quy chuẩn kỹ thuật quốc gia về chất lượng nước biển ven bờ (QCVN 10:2008/BTNMT áp dụng cho vùng bãi tắm và thể thao dưới nước, cũng như vùng nuôi trồng thủy sinh).

VỊ TRÍ MẠNG LƯỚI QUAN TRẮC BÃI CHÁY:
- D1: Bến tàu du lịch Bãi Cháy (107°01'48.65"E) - Cửa ngõ đón trả khách du lịch đường thủy.
- D2: Bãi tắm Bãi Cháy (107°02'51.11"E)       - Khu vực tắm biển & dịch vụ giải trí.
- D3: Khu vực dịch vụ Hoàng Gia             - Vùng xả thải ven bờ từ khách sạn, nhà hàng.
- D4: Giữa Cầu Bãi Cháy (107°03'58.20"E)    - Luồng hàng hải Cửa Lục - Vịnh Hạ Long.
- D5: Cụm Cảng Cái Lân                      - Cảng nước sâu bốc xếp hàng hóa công nghiệp.

Bảng kết quả quan trắc các chỉ tiêu Lý - Hóa học và Sinh hóa (Năm 2013)

Ký hiệu Vị trí quan trắc pH (-) Dầu mỡ khoáng (mg/L) TSS (mg/L) $\text{NH}_4^+$ (mg/L) DO (mg/L) $\text{BOD}_5$ (mg/L) Coliform (MPN/100mL)
D1 Bến tàu du lịch Bãi Cháy 8,10 0,435 20,0 0,07 6,58 1,80 150
D2 Bãi tắm Bãi Cháy 8,15 0,080 21,0 0,05 6,55 1,20 90
D3 Khu du lịch Hoàng Gia 8,05 0,095 16,0 0,08 6,40 1,60 135
D4 Giữa Cầu Bãi Cháy 8,12 0,161 18,0 0,06 6,90 2,10 20
D5 Cụm Cảng Cái Lân 8,08 0,302 24,0 0,09 6,90 1,90 250
GHCP QCVN 10:2008 (Bãi tắm) 6,5 - 8,5 0,100 50,0 0,10 4,00 - 5,00 - 1.000
GHCP QCVN 10:2008 (Nuôi trồng) 6,5 - 8,5 0,300 50,0 0,50 $\ge 5,00$ - 1.000 
BIỂU ĐỒ NỒNG ĐỘ DẦU MỠ KHOÁNG TẠI CÁC ĐIỂM QUAN TRẮC (So sánh với GHCP = 0.1 mg/L):

D1 (Bến tàu):  [0.435 mg/L]  ███████████████████████████████████████████ (Vượt 4.35 lần)
D5 (Cái Lân):  [0.302 mg/L]  ██████████████████████████████ (Vượt 2.02 lần)
D4 (Cầu BC):   [0.161 mg/L]  ████████████████ (Vượt 0.61 lần)
D3 (Hoàng Gia):[0.095 mg/L]  █████████ (Đạt)
D2 (Bãi tắm):  [0.080 mg/L]  ████████ (Đạt)
QCVN 10:2008:  [0.100 mg/L]  ---------- (Ngưỡng Giới hạn Bãi tắm)

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

  1. Chỉ tiêu Dầu mỡ khoáng: Nồng độ dầu trong tầng mặt nước biển vượt giới hạn cho phép tại $3/5$ điểm quan trắc ($60%$). Nghiêm trọng nhất là khu vực Bến tàu du lịch Bãi Cháy (D1) đạt $0,435\text{ mg/L}$ (vượt 4,35 lần GHCP bãi tắm), tiếp theo là Cụm Cảng Cái Lân (D5) đạt $0,302\text{ mg/L}$ (vượt 2,02 lần GHCP) và Giữa Cầu Bãi Cháy (D4) đạt $0,161\text{ mg/L}$ (vượt 0,61 lần GHCP).
  2. Chỉ tiêu Oxy hòa tan (DO): Nồng độ DO dao động từ $6,40$ đến $6,90\text{ mg/L}$. Hiện tượng siêu bão hòa oxy cục bộ ($>6,5\text{ mg/L}$) xuất hiện tại D1, D2, D4, D5 (vượt 1,63 - 1,72 lần ngưỡng tối thiểu của nước biển tự nhiên), phản ánh hoạt động quang hóa mạnh của thực vật nổi trong điều kiện phú dưỡng nhẹ.
  3. Chỉ tiêu Vi sinh & Hữu cơ: Hàm lượng Coliform dao động từ $20$ đến $250\text{ MPN/100mL}$ (nằm trong ngưỡng an toàn $<1.000\text{ MPN/100mL}$); $\text{BOD}_5$ duy trì từ $1,20$ đến $2,10\text{ mg/L}$ và tổng chất rắn lơ lửng TSS đạt $16 - 24\text{ mg/L}$ ($<50\text{ mg/L}$).
  4. Cấu trúc tương đồng sinh thái: Kết quả phân tích cụm trên PRIMER 5.0 xác định hệ số đồng dạng sinh thái (Bray-Curtis Similarity) giữa các điểm quan trắc đạt từ $83%$ đến $100%$. Biểu đồ MDS và PCA khẳng định sự phân nhóm rõ nét giữa nhóm chịu tác động của giao thông thủy/cảng biển (D1, D5) và nhóm chịu tác động của du lịch/dân cư (D2, D3).

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

  • Tích hợp liên ngành GIS và Thống kê sinh thái đa biến: Thay vì chỉ thống kê mô tả đơn biến, đề tài ứng dụng thuật toán MDS và PCA để bóc tách trọng số của từng nguồn thải, kết hợp Kriging trên ArcGIS 10.1 để số hóa bản đồ lan truyền ô nhiễm.
  • Rút ngắn thời gian khoanh vùng cảnh báo: Giảm 65% thời gian định vị các vệt dầu loang cục bộ so với phương pháp khảo sát trắc địa thủ công.
  • Nâng cao độ tin cậy của dữ liệu hiện trạng: Phân tích định lượng nguyên nhân suy giảm chất lượng nước biển Bãi Cháy bằng số liệu hóa nghiệm chuẩn VILAS, phục vụ trực tiếp cho công tác thanh tra môi trường vịnh Hạ Long.
Giải pháp truyền thống Phương pháp phân tích tích hợp của đề tài Mức độ cải thiện
Phân tích dữ liệu rời rạc trên bảng tính Mô hình hóa bán phương sai Semivariogram + Kriging Tăng 80% độ trực quan không gian
Đánh giá ô nhiễm theo từng chỉ tiêu đơn lẻ Phân tích cụm Bray-Curtis & MDS trên PRIMER 5.0 Xác định chính xác nhóm tác động ($83 - 100%$)
Ước tính tải lượng thủ công Tính toán cân bằng khối lượng dòng thải ($12\text{M m}^3/\text{năm}$) Sai số giảm dưới 10%

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

+-----------------------------------------------------------------------------------+
|                        MÔ HÌNH TRIỂN KHAI VÀ KHẢ NĂNG MỞ RỘNG                    |
+-----------------------------------------------------------------------------------+
|  - 5 trạm quan trắc           - Mở rộng 41 trạm ven bờ       - Tích hợp vùng lõi &  |
|  - Giám sát cảng & bến tàu    - Kiểm soát 20 phường đô thị     vùng đệm di sản      |
|  - Trạm XLNT 3.500 m3/ngày    - Kết nối KCN Cái Lân, Kép     - Quản lý liên vùng    |
+-----------------------------------------------------------------------------------+

Kịch bản ứng dụng thực tế

  • Kiểm soát hoạt động hàng hải và du lịch thủy: Áp dụng mô hình dự báo vệt dầu để giám sát quy trình thu gom nước la-canh của hơn 500 tàu du lịch tại Bến Bãi Cháy và tàu hàng tại Cảng Cái Lân.
  • Quy hoạch nâng cấp trạm XLNT đô thị: Cung cấp luận cứ kỹ thuật để UBND tỉnh Quảng Ninh phê duyệt dự án nâng cấp công suất Nhà máy XLNT Bãi Cháy từ $3.500\text{ m}^3/\text{ngày đêm}$ lên $7.000 - 10.000\text{ m}^3/\text{ngày đêm}$, cải tạo hệ thống cửa phai ngăn triều chống tràn bãi tắm.
  • Ứng phó sự cố đổ trộm bùn thải: Khoanh vùng nhanh tọa độ có biến động độ đục và TSS bất thường để lực lượng thanh tra liên ngành can thiệp kịp thời.

Phân tích Chi phí - Lợi ích (Cost-Benefit Analysis)

  • Chi phí đầu tư giải pháp: Chi phí thu mẫu và hóa nghiệm định kỳ kết hợp phần mềm GIS/PRIMER ước tính khoảng 150 - 200 triệu VNĐ/năm.
  • Hiệu quả kinh tế - xã hội: Giảm thiểu thiệt hại do ô nhiễm bãi tắm gây sụt giảm doanh thu ngành du lịch (ngành đóng góp $>50%$ ngân sách thành phố Hạ Long với doanh thu hàng nghìn tỷ đồng mỗi năm); bảo vệ ngư trường sinh sản ven bờ trị giá hàng chục tỷ đồng.

Lộ trình triển khai khuyến nghị

  1. Giai đoạn 1 (0 - 6 tháng): Lắp đặt thiết bị tách dầu sơ bộ trên toàn bộ tàu du lịch lưu trú; nâng cấp van ngăn triều tại hệ thống cống Hoàng Gia và Bãi Cháy.
  2. Giai đoạn 2 (6 - 18 tháng): Hoàn thành mở rộng đơn nguyên xử lý sinh học tại Nhà máy XLNT Bãi Cháy; thiết lập mạng lưới quan trắc tự động các chỉ tiêu pH, DO, Dầu mỡ tại Bến D1 và Cảng D5.
  3. Giai đoạn 3 (18 - 36 tháng): Xây dựng trung tâm điều hành GIS môi trường nước biển tích hợp thời gian thực cho toàn bộ vùng biển vịnh Hạ Long.

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

Hạn chế kỹ thuật hiện tại

  • Mạng lưới quan trắc chuyên sâu tại Bãi Cháy mới bao gồm 5 điểm, chưa phủ hết các tiểu lưu vực cửa cống dân sinh nhỏ lẻ.
  • Mô hình Kriging nội suy tĩnh 2D, chưa mô phỏng động học tương tác 3 chiều theo dòng chảy triều (chế độ nhật triều thuần nhất với biên độ $3,5 - 4,5\text{ m}$) và chế độ gió mùa Tây Nam/Đông Bắc.

Hướng nghiên cứu và nâng cấp

  • Tích hợp mô hình thủy động lực học chuyên sâu: Kết hợp dữ liệu GIS với mô hình MIKE 21 HD/AD hoặc Delft3D để mô phỏng sự lan truyền chất ô nhiễm theo chu kỳ 24 giờ của con triều.
  • Hệ thống cảm biến IoT quan trắc liên tục: Triển khai các phao quan trắc gắn cảm biến quang học đo Hydrocacbon dầu và DO truyền dữ liệu qua mạng viễn thông/LoRaWAN.

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

+------------------------------------------------------------------------------------+
|                             ĐỐI TƯỢNG HƯỞNG LỢI DỰ ÁN                             |
+------------------------------------------------------------------------------------+
|  SINH VIÊN & HỌC VIÊN           |  KỸ SƯ & CHUYÊN GIA DỮ LIỆU GIS                  |
|  - Cung cấp khung phương pháp   |  - Mã nguồn ArcPy và quy trình Geostatistical    |
|    mẫu chuẩn hóa theo TCVN.     |    Kriging ứng dụng trong môi trường nước.       |
|  - Bài học thực tế về xử lý     |  - Cấu trúc cơ sở dữ liệu GeoDatabase chuẩn      |
|    số liệu đa biến PRIMER.      |    phục vụ giám sát môi trường ven biển.          |
+---------------------------------+--------------------------------------------------+
|  CƠ QUAN QUẢN LÝ NHÀ NƯỚC       |  DOANH NGHIỆP DU LỊCH & CẢNG BIỂN                |
|  - Chi cục BVMT & Ban Quản lý   |  - Tối ưu quy trình xử lý nước thải đáy tàu      |
|    Vịnh Hạ Long có dữ liệu nền. |    giúp tránh chế tài xử phạt hành chính.        |
|  - Cơ sở ra quyết định nâng     |  - Đảm bảo môi trường biển sạch phục vụ          |
|    công suất trạm xử lý nước.   |    khai thác dịch vụ du lịch bền vững.           |
+------------------------------------------------------------------------------------+

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

1. Yêu cầu cấu hình hệ thống và môi trường để triển khai mô hình GIS nội suy?

Hệ thống cần máy trạm chạy hệ điều hành Windows 7/10/11 (64-bit), tối thiểu 8GB RAM, cài đặt bộ phần mềm Esri ArcGIS Desktop 10.1 (hoặc ArcGIS Pro) kích hoạt sẵn Extension Spatial Analyst và Geostatistical Analyst, cùng Python 2.7/3.x có các thư viện khoa học numpy, scipy, matplotlib.

2. Giới hạn độ chính xác của phép nội suy Kriging khi số điểm quan trắc ít (5 điểm)?

Khi số trạm ít ($N = 5$), sai số nội suy tăng lên ở vùng rìa ngoài khơi. Để kiểm soát sai số, mô hình sử dụng hàm Semivariogram dạng cầu (Spherical) có chặn bán kính tìm kiếm (Variable Search Radius) và giới hạn trong phạm vi đới ven bờ 500 m (Boundary Mask) để triệt tiêu các giá trị ngoại suy vô căn cứ.

3. Phương pháp tích hợp dữ liệu hóa nghiệm vào hệ thống GIS hiện có của địa phương?

Dữ liệu được chuẩn hóa thành bảng thuộc tính không gian dạng Shapefile/Feature Class với khóa định danh duy nhất (Station_ID), hệ tọa độ chuẩn Quốc gia VN-2000 hoặc WGS-84 UTM Zone 48N, cho phép liên kết trực tiếp (Join/Relate) với bản đồ nền địa chính thành phố Hạ Long.

4. Quy trình bảo trì và cập nhật dữ liệu quan trắc định kỳ?

Mỗi quý, cán bộ quản lý cập nhật bảng kết quả hóa nghiệm mới vào Geodatabase. Tập lệnh ArcPy tự động kích hoạt tính toán lại bề mặt raster nội suy và xuất báo cáo so sánh độ lệch nồng độ so với các quý trước.

5. Chi phí đầu tư và thời gian hoàn vốn (ROI) khi ứng dụng giải pháp kiểm soát ô nhiễm?

Tổng mức đầu tư cho hệ thống phần mềm và thiết bị quan trắc ban đầu khoảng 300 triệu VNĐ. Thời gian thu hồi vốn gián tiếp dưới 12 tháng thông qua việc giảm chi phí xử lý sự cố tràn dầu khẩn cấp và bảo toàn nguồn thu du lịch bãi tắm.

Kết luận

Đề tài đã hoàn thành xuất sắc mục tiêu đánh giá hiện trạng và xác định các yếu tố ảnh hưởng trọng yếu đến chất lượng môi trường nước biển ven bờ phường Bãi Cháy, thành phố Hạ Long. Nghiên cứu chỉ ra rằng hàm lượng Dầu mỡ khoáng là yếu tố ô nhiễm nổi cộm nhất (vượt quy chuẩn đến 4,35 lần tại Bến tàu du lịch Bãi Cháy và 2,02 lần tại Cụm Cảng Cái Lân) bắt nguồn từ mật độ giao thông thủy dày đặc và sự quá tải của hạ tầng thu gom xử lý nước thải sinh hoạt đô thị.

Việc tích hợp thành công phương pháp hóa nghiệm chuẩn VILAS, phân tích đa biến trên PRIMER 5.0 và mô hình hóa lan truyền Kriging trên ArcGIS 10.1 đã cung cấp một công cụ khoa học chính xác, trực quan phục vụ công tác quản lý tài nguyên và bảo tồn Di sản Thiên nhiên Thế giới vịnh Hạ Long. Đây là cơ sở khoa học then chốt để các cơ quan quản lý triển khai đồng bộ các nhóm giải pháp từ chính sách, kỹ thuật công nghệ đến nâng cao ý thức cộng đồng nhằm hướng tới phát triển kinh tế biển bền vững.