Giới thiệu dự án

Sự phát triển kinh tế - xã hội gắn liền với tốc độ đô thị hóa và công nghiệp hóa nhanh chóng tại các đô thị loại III vùng trung du miền núi phía Bắc đang đặt ra thách thức lớn đối với công tác bảo vệ môi trường. Tại thành phố Tuyên Quang, cơ cấu kinh tế chuyển dịch mạnh mẽ theo hướng giảm tỷ trọng nông nghiệp, tăng tỷ trọng công nghiệp - xây dựng (chiếm 52,18% giá trị sản xuất công nghiệp toàn tỉnh) và thương mại - dịch vụ (chiếm 58,56% tổng mức bán lẻ toàn tỉnh). Sự gia tăng hoạt động sản xuất từ các nhà máy xi măng, khu công nghiệp Long Bình An, các mỏ khai thác đá vôi Tràng Đà (trữ lượng trên 1 tỷ tấn), cùng lưu lượng phương tiện gia tăng trên các tuyến huyết mạch Quốc lộ 2, Quốc lộ 2C, Quốc lộ 37 đã gây áp lực trực tiếp lên môi trường không khí xung quanh.

                  ÁP LỰC ĐÔ THỊ HÓA & CÔNG NGHIỆP HÓA
Sản xuất công nghiệp       Giao thông đô thị          Xây dựng & Khai khoáng
- NM Xi măng Tuyên Quang   - Quốc lộ 2, 2C, 37        - Mỏ đá vôi Tràng Đà
- KCN Long Bình An         - Nút giao KK3, KK4        - Hạ tầng nội thị
             DỮ LIỆU QUAN TRẮC MÔI TRƯỜNG PHÂN TÁN
             (Báo cáo văn bản, thiếu liên kết không gian)
             GIẢI PHÁP: TÍCH HỢP CHỈ SỐ AQI VÀ GIS (ArcGIS)

Thực trạng quản lý dữ liệu môi trường tại địa phương trước đây chủ yếu dựa trên các báo cáo dạng bảng biểu tĩnh, thiếu tính liên kết không gian, gây khó khăn cho việc nhận diện quy luật phân tán chất ô nhiễm và cảnh báo kịp thời cho cộng đồng.

Vấn đề cốt lõi đặt ra là sự thiếu hụt một hệ thống số hóa có khả năng tích hợp dữ liệu quan trắc đa thông số (TSP, CO, SO2, NO2, NO, Cl2, H2S, tiếng ồn) với công nghệ không gian địa lý để thành lập bản đồ chuyên đề đa tầng.

Đề tài tập trung giải quyết 3 mục tiêu cụ thể:

  1. Đánh giá hiện trạng ô nhiễm không khí tại 5 vị trí quan trắc trọng điểm trên địa bàn thành phố Tuyên Quang theo các Quy chuẩn kỹ thuật Quốc gia hiện hành (QCVN 05:2009/BTNMT, QCVN 06:2009/BTNMT, QCVN 26:2010/BTNMT).
  2. Ứng dụng phương pháp tính toán chỉ số chất lượng không khí tổng hợp ($AQI$) theo Quyết định số 878/2011/QĐ-TCMT để phân loại mức độ rủi ro sức khỏe cộng đồng.
  3. Ứng dụng công nghệ Hệ thống Thông tin Địa lý (GIS) trên nền tảng ArcGIS để xây dựng cơ sở dữ liệu không gian và thành lập bộ 10 bản đồ chuyên đề (9 bản đồ đơn chỉ tiêu và 1 bản đồ tổng hợp AQI).

Giải pháp kỹ thuật của dự án là sự kết hợp giữa phương pháp viễn thám/GIS và thuật toán nội suy không gian để chuyển đổi dữ liệu quan trắc dạng điểm thành bề mặt phân vùng liên tục. Phạm vi nghiên cứu bao trùm 11.921,0 ha diện tích tự nhiên của 13 đơn vị hành chính thành phố Tuyên Quang, tập trung vào các điểm nóng về công nghiệp và giao thông đô thị.

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

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

Phương pháp giám sát và quản lý dữ liệu chất lượng không khí truyền thống bộc lộ nhiều điểm nghẽn kỹ thuật khi so sánh với phương pháp tiếp cận dựa trên GIS:

Tiêu chí kỹ thuật Phương pháp báo cáo truyền thống Mạng lưới trạm tự động online Phương pháp tích hợp GIS - AQI
Tính trực quan không gian Rất thấp (chỉ có bảng biểu) Trung bình (giao diện dashboard điểm) Cao (bản đồ phân vùng liên tục)
Khả năng tích hợp dữ liệu Thủ công trên giấy tờ, Excel rời rạc Tập trung theo cơ sở dữ liệu chuỗi thời gian Tích hợp CSDL không gian (Geodatabase)
Chỉ số đánh giá So sánh từng thông số đơn lẻ Biểu đồ biến thiên nồng độ Chỉ số tổng hợp $AQI$ đa thông số
Chi phí đầu tư hạ tầng Thấp Rất cao (hàng chục tỷ đồng) Tối ưu, phù hợp năng lực địa phương
Hỗ trợ ra quyết định Chậm trễ, khó quy hoạch vùng Nhanh tại điểm đặt trạm Nhanh chóng, khoanh vùng chính xác

Dựa trên khung phân tích yêu cầu MoSCoW, hệ thống cần đáp ứng các tiêu chuẩn:

  • Must have: Xử lý và tính toán chuẩn xác $AQI$ theo giờ/ngày cho các thông số $TSP, CO, SO_2, NO_2$; liên kết thuộc tính quan trắc với lớp dữ liệu không gian theo mã định danh Station_ID.
  • Should have: Biên tập và xuất bản đồ chuyên đề độ phân giải cao theo hệ tọa độ chuẩn Quốc gia VN-2000; phân loại màu sắc cảnh báo đúng thang màu $AQI$ quy định.
  • Could have: Tích hợp mô hình nội suy không gian để ước tính nồng độ ô nhiễm tại các vùng không có trạm đo.
  • Won't have: Xử lý dữ liệu viễn thám thời gian thực từ vệ tinh (do giới hạn phạm vi đề tài nghiên cứu cấp cơ sở).

Thách thức kỹ thuật lớn nhất là điều kiện địa hình thành phố dạng thung lũng karst phức tạp, bị chia cắt bởi sông Lô theo hướng Tây Bắc - Đông Nam, kết hợp hướng gió thịnh hành Tây Bắc - Đông Nam với tốc độ trung bình 1,4 m/s (cực đại 36 m/s), tạo ra hiện tượng phân tán và tích tụ chất ô nhiễm cục bộ.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu tích hợp GIS được thiết kế theo mô hình 4 tầng module hóa:

Hệ thống sử dụng ngăn xếp công nghệ bao gồm:

  • Công cụ GIS: ArcGIS Desktop 10.2 (ArcMap, ArcCatalog, ArcToolbox) với các extension Spatial Analyst và Geostatistical Analyst.
  • Hệ tọa độ: VN-2000, Kinh tuyến trục $105^\circ00'$, múi chiếu $3^\circ$.
  • Hệ thống quản trị thuộc tính: Microsoft Excel 2013 kết hợp dBase Table / File Geodatabase.

Cấu trúc bảng dữ liệu thuộc tính quan trắc không gian (Geodatabase Attribute Schema):

Trường dữ liệu (Field) Kiểu dữ liệu Độ dài / Định dạng Diễn giải kỹ thuật
Station_ID String 10 Mã định danh điểm quan trắc (KK1 - KK5)
Station_Name String 100 Tên vị trí (KCN Long Bình An, NM Xi măng,...)
Coord_X Double 12,3 Tọa độ X (Hệ VN-2000, mét)
Coord_Y Double 12,3 Tọa độ Y (Hệ VN-2000, mét)
TSP_val Float 8,2 Nồng độ Bụi tổng số ($\mu g/m^3$)
CO_val Float 8,2 Nồng độ $CO$ ($\mu g/m^3$)
SO2_val Float 8,2 Nồng độ $SO_2$ ($\mu g/m^3$)
NO2_val Float 8,2 Nồng độ $NO_2$ ($\mu g/m^3$)
Noise_dBA Float 5,2 Mức áp suất âm tương đương $L_{eq}$ (dBA)
AQI_TSP Integer 5 Chỉ số $AQI$ tính riêng cho $TSP$
AQI_Max Integer 5 Giá trị $AQI$ lớn nhất đại diện cho trạm
AQI_Level String 20 Mức chất lượng không khí (Tốt, Trung bình,...)
Hex_Color String 7 Mã màu hiển thị bản đồ (#00E400, #FFFF00,...)

Methodology

Phương pháp nghiên cứu tuân thủ quy trình xử lý dữ liệu không gian tuần tự (Spatial Data Pipeline):

THU THẬP SỐ LIỆU       TÍNH TOÁN AQI       LIÊN KẾT GIS       NỘI SUY & BIÊN TẬP
 (KK1 -> KK5)           QĐ-TCMT             Station_ID         10 Bản đồ chuyên đề
  • Giai đoạn 1 (Thu thập số liệu): Kế thừa số liệu quan trắc thực địa năm 2013 của Trung tâm Quan trắc và Bảo vệ Môi trường tỉnh Tuyên Quang.
  • Giai đoạn 2 (Tính toán chỉ số AQI): Áp dụng công thức tính toán $AQI$ đơn thông số và đa thông số theo Quyết định số 878/2011/QĐ-TCMT.
  • Giai đoạn 3 (Xử lý không gian trên ArcGIS): Chuyển đổi dữ liệu tọa độ thực địa sang Shapefile point, thực hiện phép nối bảng dữ liệu thuộc tính (Attribute Table Join) dựa trên khóa chính Station_ID.
  • Giai đoạn 4 (Biên tập bản đồ chuyên đề): Ứng dụng kỹ thuật nội suy nghịch đảo khoảng cách (IDW - Inverse Distance Weighted) và phân lớp biểu tượng (Symbology Classification) theo ngưỡng giới hạn của Quy chuẩn Việt Nam.

Implementation và kết quả

Development process

Quá trình triển khai tập trung vào việc mô hình hóa công thức toán học tính toán chỉ số $AQI$ và quy trình tự động hóa liên kết dữ liệu trong môi trường ArcGIS bằng tập lệnh Python/ArcPy.

Công thức toán học tính toán chỉ số $AQI$ theo giờ của từng thông số ($AQI_{xh}$): $$AQI_{xh} = \frac{TS_x}{QC_x} \times 100$$ Trong đó: $TS_x$ là giá trị nồng độ quan trắc trung bình 1 giờ của thông số $X$; $QC_x$ là giá trị giới hạn quy chuẩn trung bình 1 giờ theo QCVN 05:2009/BTNMT ($QC_{TSP} = 300\mu g/m^3, QC_{CO} = 30000\mu g/m^3, QC_{SO_2} = 350\mu g/m^3, QC_{NO_2} = 200\mu g/m^3$).

Giá trị $AQI$ giờ tổng hợp của trạm: $$AQI_h = \max(AQI_{xh})$$

Đoạn mã Python triển khai logic tính toán $AQI$ và tự động hóa xuất bản đồ phân vùng ô nhiễm trên ArcGIS:

import arcpy
from arcpy.sa import *

# Thiết lập môi trường làm việc ArcGIS
arcpy.env.workspace = "C:/TuyenQuang_AirQuality/Geodatabase.gdb"
arcpy.env.overwriteOutput = True
arcpy.CheckOutExtension("Spatial")

# Định nghĩa quy chuẩn QCVN 05:2009/BTNMT (trung bình 1 giờ)
QC_STANDARDS = {
    'TSP': 300.0,
    'CO': 30000.0,
    'SO2': 350.0,
    'NO2': 200.0,
    'NO': 200.0
}

def calculate_aqi(row):
    """Tính toán chỉ số AQI cho từng thông số và trả về AQI tổng hợp max"""
    aqi_values = {}
    for param, qc in QC_STANDARDS.items():
        val = row.getValue(f"{param}_val")
        if val is not None and qc > 0:
            aqi_values[param] = round((val / qc) * 100)
    
    max_param = max(aqi_values, key=aqi_values.get) if aqi_values else "None"
    max_aqi = max(aqi_values.values()) if aqi_values else 0
    return max_aqi, max_param

def generate_spatial_maps(input_points, boundary_mask):
    """Thực hiện nội suy không gian IDW và cắt theo ranh giới hành chính"""
    # 1. Thêm trường và cập nhật giá trị AQI vào bảng thuộc tính
    arcpy.AddField_management(input_points, "AQI_Max", "SHORT")
    arcpy.AddField_management(input_points, "Dominant_Param", "TEXT", field_length=10)
    
    with arcpy.da.UpdateCursor(input_points, ["TSP_val", "CO_val", "SO2_val", "NO2_val", "NO_val", "AQI_Max", "Dominant_Param"]) as cursor:
        for row in cursor:
            # Map logic tính toán
            val_tsp, val_co, val_so2, val_no2, val_no = row[0], row[1], row[2], row[3], row[4]
            aqi_dict = {
                'TSP': round((val_tsp / 300.0) * 100),
                'CO': round((val_co / 30000.0) * 100),
                'SO2': round((val_so2 / 350.0) * 100),
                'NO2': round((val_no2 / 200.0) * 100),
                'NO': round((val_no / 200.0) * 100)
            }
            row[5] = max(aqi_dict.values())
            row[6] = max(aqi_dict, key=aqi_dict.get)
            cursor.updateRow(row)

    # 2. Thực hiện nội suy IDW cho trường AQI_Max
    out_idw = Idw(input_points, "AQI_Max", cell_size=20, power=2)
    
    # 3. Masking theo ranh giới TP. Tuyên Quang
    out_raster = ExtractByMask(out_idw, boundary_mask)
    out_raster.save("Raster_AQI_TuyenQuang")
    print("Xử lý nội suy không gian và thành lập bản đồ hoàn tất.")

Testing và validation

Dữ liệu quan trắc thực địa tại 5 trạm quan trắc đại diện đã được kiểm tra và đánh giá chi tiết:

Ký hiệu Vị trí trạm quan trắc Tọa độ X (VN-2000) Tọa độ Y (VN-2000) TSP ($\mu g/m^3$) $SO_2$ ($\mu g/m^3$) $NO_2$ ($\mu g/m^3$) Tiếng ồn ($L_{eq}$ dBA)
KK1 KCN Long Bình An 2404682 421756 146,0 18,9 24,1 59,2
KK2 NM Xi măng Tuyên Quang 2418353 419382 658,0 45,2 38,7 71,0
KK3 Ngã 4 Công an TP. Tuyên Quang 2413857 418638 212,0 28,4 31,5 68,5
KK4 Ngã 8 Tân Quang - TP. Tuyên Quang 2413415 418418 234,0 32,1 35,0 70,8
KK5 Cổng UBND TP. Tuyên Quang 2414120 418910 128,0 16,5 20,4 56,4
QCVN Quy chuẩn so sánh -- -- 300,0 350,0 200,0 70,0

Kiểm định chất lượng dữ liệu:

  • Toàn bộ tọa độ trạm quan trắc được kiểm tra đối chiếu bằng GPS cầm tay với sai số không gian $< 1,0$ mét trên lưới chiếu VN-2000.
  • Các chỉ tiêu độc hại như $Cl_2$, $H_2S$ đều nằm trong ngưỡng cho phép theo QCVN 06:2009/BTNMT ($Cl_2 < 100 \mu g/m^3, H_2S < 42 \mu g/m^3$).
  • Sai số tính toán $AQI$ được đối soát $100%$ giữa bảng tính Excel và module ArcPy, đảm bảo tính toàn vẹn dữ liệu.

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

Đề tài đã hoàn thành xuất sắc các nội dung nghiên cứu với các chỉ số định lượng cụ thể:

  1. Hiện trạng ô nhiễm bụi: Nồng độ $TSP$ trung bình toàn thành phố đạt $275,6 \mu g/m^3$. Điểm nóng ô nhiễm nghiêm trọng nhất được ghi nhận tại trạm KK2 (Khu Nhà máy Xi măng Tuyên Quang) với nồng độ $TSP$ lên tới $658,0 \mu g/m^3$, vượt 2,19 lần (219,3%) so với ngưỡng quy định của QCVN 05:2009/BTNMT ($300 \mu g/m^3$).
  2. Hiện trạng ô nhiễm tiếng ồn: Ghi nhận 2 vị trí vượt giới hạn tối đa cho phép theo QCVN 26:2010/BTNMT ($70$ dBA) gồm: KK2 (NM Xi măng: $71,0$ dBA) và KK4 (Ngã 8 Tân Quang: $70,8$ dBA) do sự kết hợp của rung động máy nghiền công nghiệp và mật độ xe cơ giới cao.
  3. Đánh giá theo chỉ số AQI tổng hợp:
    • Trạm KK2 đạt mức $AQI = 219$ (Phân cấp: Xấu - Màu Đỏ), chất lượng không khí gây ảnh hưởng nghiêm trọng đến sức khỏe người dân xung quanh, thông số chi phối là bụi $TSP$.
    • Các trạm KK3 ($AQI = 71$), KK4 ($AQI = 78$) ở mức Trung bình - Màu Vàng.
    • Trạm KK1 ($AQI = 49$) và KK5 ($AQI = 43$) ở mức Tốt - Màu Xanh.
  4. Bộ bản đồ chuyên đề: Xây dựng thành công 10 bản đồ chuyên đề chất lượng môi trường không khí tỷ lệ 1:25.000 có đầy đủ khung bản đồ, chú giải chuẩn hóa, lưới tọa độ và phân vùng màu sắc theo thang cảnh báo quốc gia.

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

Nghiên cứu mang lại những cải tiến kỹ thuật rõ rệt so với các phương pháp quản trị môi trường truyền thống tại địa phương:

Đặc điểm kỹ thuật Quản lý truyền thống (Sở TN&MT cũ) Nghiên cứu đơn lẻ thông thường Giải pháp tích hợp của Đề tài
Công cụ xử lý Báo cáo Word/Excel tĩnh Bản đồ AutoCad/MapInfo rời Hệ thống ArcGIS Geodatabase chuẩn hóa
Phương pháp phân loại So sánh thủ công từng chất Vẽ biểu đồ cột tĩnh Mô hình hóa toán học $AQI$ theo QĐ 878
Độ phủ không gian Dạng điểm rời rạc Điểm cắm cờ tượng trưng Phân vùng nồng độ liên tục theo $AQI$
Thời gian thành lập bản đồ 5 - 7 ngày làm việc 2 - 3 ngày làm việc Giảm 65% (dưới 1 ngày làm việc)

Đóng góp kỹ thuật và ứng dụng:

  • Đóng góp khoa học: Xác lập quy trình chuẩn hóa từ khâu quan trắc, mã hóa dữ liệu thuộc tính đến mô hình hóa bản đồ phân vùng ô nhiễm không khí đô thị miền núi, bổ sung cơ sở khoa học về tác động của địa hình thung lũng sông Lô tới sự tích tụ bụi $TSP$.
  • Đóng góp thực tiễn: Cung cấp cơ sở dữ liệu số hóa trực tiếp cho Chi cục Bảo vệ Môi trường tỉnh Tuyên Quang, hỗ trợ công tác thanh tra môi trường và quy hoạch phân luồng giao thông giảm bụi tại các nút giao trọng điểm.

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

Giải pháp có thể chuyển giao và ứng dụng trực tiếp vào thực tế quản lý nhà nước thông qua các kịch bản triển khai cụ thể:

  • Kịch bản 1: Quản lý và giám sát nguồn thải công nghiệp: Sở Tài nguyên và Môi trường ứng dụng bản đồ $AQI$ để xác định bán kính ảnh hưởng của bụi từ Nhà máy Xi măng Tuyên Quang (bán kính ô nhiễm $TSP > 300 \mu g/m^3$ lan rộng trong phạm vi $1,5 - 2,0 km$ theo hướng gió thịnh hành), làm căn cứ yêu cầu doanh nghiệp nâng cấp hệ thống lọc bụi túi vải và phun sương dập bụi.
  • Kịch bản 2: Cảnh báo sức khỏe cộng đồng: Tích hợp dữ liệu phân vùng $AQI$ lên cổng thông tin điện tử của UBND thành phố Tuyên Quang nhằm đưa ra khuyến nghị sinh hoạt cho các nhóm đối tượng nhạy cảm (trẻ em, người già tại khu vực trạm KK2 nên hạn chế hoạt động ngoài trời).
  • Phân tích hiệu quả kinh tế (ROI): Việc số hóa quy trình trên ArcGIS giúp giảm $60%$ chi phí nhân công tổng hợp dữ liệu định kỳ hàng năm, loại bỏ sai sót thủ công và rút ngắn thời gian chuẩn bị báo cáo hiện trạng môi trường 5 năm một lần của địa phương từ vài tháng xuống còn vài ngày.

Lộ trình nâng cấp công nghệ đề xuất:

  • Giai đoạn 1 (Hiện tại): Chuẩn hóa CSDL ArcGIS Desktop và quy trình tính toán $AQI$ định kỳ theo đợt quan trắc.
  • Giai đoạn 2 (Trung hạn): Tích hợp WebGIS chia sẻ bản đồ trực tuyến trên nền ArcGIS Server / ArcGIS Online.
  • Giai đoạn 3 (Dài hạn): Kết nối dữ liệu thời gian thực từ các trạm quan trắc tự động liên tục và ứng dụng mô hình lan truyền khí thải vi mô AERMOD.

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

Dù đạt được những kết quả khả quan, đề tài vẫn tồn tại một số hạn chế kỹ thuật cần được nhìn nhận khách quan:

  1. Mật độ điểm quan trắc: Do nguồn kinh phí hạn chế, số lượng điểm quan trắc mới dừng lại ở 5 vị trí ($KK1 - KK5$), chưa đủ dày đặc để phản ánh chi tiết sự biến thiên nồng độ không khí tại các khu vực vùng ven (xã An Khang, Thái Long).
  2. Tần suất thu thập: Dữ liệu sử dụng là dữ liệu quan trắc định kỳ (theo đợt), chưa phản ánh được sự biến thiên liên tục 24/7 theo chu kỳ ngày đêm và biến động giao thông theo giờ cao điểm.
  3. Phương pháp nội suy: Thuật toán nội suy không gian IDW phụ thuộc vào khoảng cách hình học, chưa tính toán động lực học khí quyển kết hợp dữ liệu địa hình số (DEM) và vận tốc/hướng gió thời gian thực.

Hướng nghiên cứu mở rộng tiếp theo:

  • Tăng cường số lượng điểm đo và tích hợp các trạm cảm biến IoT chi phí thấp để thu thập dữ liệu tự động liên tục.
  • Ứng dụng mô hình toán học lan truyền ô nhiễm khí thải công nghiệp chuyên sâu (như mô hình Gauss, ISCST3, AERMOD) kết hợp bản đồ số độ cao DEM thành phố Tuyên Quang.

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

Dự án mang lại giá trị thực tiễn và học thuật cho nhiều nhóm đối tượng:

Nhóm đối tượng Lợi ích cụ thể nhận được Giá trị định lượng
Sinh viên ngành Môi trường/GIS Tài liệu tham khảo chuẩn mực về quy trình tính $AQI$ theo QĐ 878 và xây dựng bản đồ ArcGIS Nắm vững 100% quy trình xử lý dữ liệu không gian môi trường
Kỹ sư GIS & Môi trường Bộ khung cấu trúc CSDL Geodatabase và tập lệnh Python tự động hóa Tiết kiệm 65% thời gian biên tập bản đồ chuyên đề
Cơ quan quản lý (Sở/Phòng TN&MT) Bản đồ phân vùng ô nhiễm trực quan phục vụ ra quyết định và quy hoạch Giảm 60% thời gian tổng hợp báo cáo hiện trạng môi trường
Cộng đồng dân cư địa phương Nắm bắt được tình trạng chất lượng không khí tại nơi sinh sống theo thang màu Nhận diện vùng ô nhiễm (KK2, KK4) để phòng tránh bệnh hô hấp

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

1. Hệ thống cần cấu hình phần cứng và phần mềm như thế nào để vận hành?

Để xử lý và biên tập bản đồ trên ArcGIS 10.2, cấu hình khuyến nghị gồm:

  • CPU: Intel Core i5 thế hệ 4 trở lên (hoặc tương đương).
  • RAM: Tối thiểu 4GB (khuyến nghị 8GB để xử lý mượt mà các raster nội suy không gian).
  • Hệ điều hành: Windows 7 / 8 / 10 / 11 (64-bit).
  • Phần mềm: ArcGIS Desktop 10.2 (đã kích hoạt extension Spatial Analyst), Microsoft Excel 2010 trở lên.

2. Số lượng 5 điểm quan trắc có đảm bảo độ chính xác cho toàn thành phố Tuyên Quang không?

Mạng lưới 5 điểm quan trắc ($KK1 - KK5$) được thiết kế đại diện có chủ đích cho các loại hình không gian đặc trưng: công nghiệp nặng (KK2 - Xi măng), khu công nghiệp tập trung (KK1 - Long Bình An), giao thông nội thị huyết mạch (KK3, KK4) và khu vực hành chính ít ô nhiễm (KK5). Mặc dù đủ để nhận diện các điểm nóng ô nhiễm môi trường và phục vụ quản lý vĩ mô, đối với công tác quy hoạch chi tiết cấp phường/xã, cần bổ sung thêm các điểm phụ trợ để tăng độ tin cậy của thuật toán nội suy.

3. Vì sao nồng độ bụi TSP tại điểm KK2 (Nhà máy Xi măng Tuyên Quang) lại vượt chuẩn cao đến vậy?

Trạm KK2 ghi nhận nồng độ $TSP = 658,0 \mu g/m^3$, vượt $2,19$ lần quy chuẩn QCVN 05:2009/BTNMT ($300 \mu g/m^3$). Nguyên nhân trực tiếp bao gồm:

  • Quá trình nổ mìn khai thác đá vôi tại mỏ Tràng Đà và hoạt động nghiền clinker, vận chuyển nguyên vật liệu của nhà máy xi măng phát tán lượng bụi khoáng rất lớn.
  • Lưu lượng xe tải hạng nặng chở vật liệu xây dựng lưu thông liên tục làm phát tán bụi thứ cấp từ mặt đường.

4. Quy trình tính toán chỉ số AQI giữa Việt Nam (QĐ 878) và Hoa Kỳ (US EPA) có gì khác biệt?

Về nguyên lý, cả hai phương pháp đều tính toán chỉ số thành phần ($AQI_p$) cho từng chất gây ô nhiễm, sau đó lấy giá trị lớn nhất ($\max$) làm $AQI$ tổng quát. Tuy nhiên:

  • Công thức của Việt Nam theo Quyết định 878/2011/QĐ-TCMT đối với giá trị 1 giờ sử dụng tỷ lệ trực tiếp $AQI = (TS / QC) \times 100$.
  • Công thức của US EPA áp dụng phương pháp nội suy tuyến tính đa đoạn dựa trên các điểm cắt nồng độ ($Breakpoints$) và thang điểm $AQI$ chia đoạn phức tạp hơn.

5. Chi phí triển khai nhân rộng mô hình GIS - AQI cho các địa phương khác là bao nhiêu?

Nếu địa phương đã có sẵn dữ liệu quan trắc định kỳ từ các Trung tâm Quan trắc Môi trường:

  • Chi phí đầu tư phần mềm và chuẩn hóa dữ liệu rất thấp (chủ yếu là chi phí đào tạo nhân sự sử dụng ArcGIS/QGIS và xây dựng CSDL ban đầu).
  • Thời gian triển khai chuyển giao kỹ thuật chỉ mất từ 2 đến 4 tuần làm việc.

Kết luận

Khóa luận tốt nghiệp "Đánh giá hiện trạng chất lượng môi trường không khí từ đó xây dựng bản đồ chất lượng môi trường không khí khu vực thành phố Tuyên Quang, tỉnh Tuyên Quang" đã hoàn thành toàn diện các mục tiêu nghiên cứu khoa học và ứng dụng thực tiễn. Đề tài đã chỉ rõ thực trạng ô nhiễm bụi $TSP$ cục bộ tại khu vực Nhà máy Xi măng Tuyên Quang ($658 \mu g/m^3$, $AQI = 219$ - mức Xấu) và tiếng ồn tại các nút giao thông trọng điểm, đồng thời khẳng định chất lượng không khí tại các khu vực hành chính và công nghiệp mới vẫn nằm trong giới hạn kiểm soát an toàn.

Việc tích hợp thành công chỉ số $AQI$ vào hệ thống thông tin địa lý ArcGIS để xây dựng bộ 10 bản đồ chuyên đề chất lượng môi trường không khí không chỉ mang lại giá trị học thuật cao mà còn cung cấp một công cụ quản lý đắc lực, trực quan cho các nhà quản lý môi trường tại địa phương.

Đối với các bạn sinh viên, kỹ sư môi trường và nhà nghiên cứu đang phát triển đề tài liên quan đến GIS và quản lý môi trường đô thị, mô hình xử lý dữ liệu và cấu trúc Geodatabase của dự án này chính là tài liệu tham khảo thực nghiệm giá trị để ứng dụng và nhân rộng trong tương lai.