Giới thiệu dự án

Sự gia tăng nhanh chóng của tốc độ đô thị hóa, mật độ phương tiện cơ giới và hoạt động xây dựng hạ tầng kỹ thuật tại các đô thị miền núi đang đặt ra thách thức nghiêm trọng đối với công tác kiểm soát ô nhiễm không khí. Tại thành phố Sơn La (tỉnh Sơn La) – trung tâm kinh tế, chính trị của khu vực Tây Bắc với đặc trưng địa hình lòng chảo karst chia cắt phức tạp (độ cao trung bình 700 - 800 m so với mực nước biển), sự gia tăng lưu lượng giao thông trên trục Quốc lộ 6 cùng hoạt động khai thác khoáng sản, sản xuất vật liệu xây dựng đã làm suy giảm đáng kể chất lượng môi trường không khí xung quanh.

Công tác quan trắc môi trường truyền thống chủ yếu dựa trên các điểm đo mặt đất riêng lẻ, tạo ra các tập dữ liệu rời rạc, thiếu tính liên tục không gian và đòi hỏi chi phí đầu tư mạng lưới trạm cố định rất tốn kém. Nhằm giải quyết triệt để bài toán này, đề tài “Ứng dụng GIS đánh giá chất lượng môi trường không khí tại thành phố Sơn La, tỉnh Sơn La” do sinh viên Lã Đức Cường thực hiện dưới sự hướng dẫn của PGS.TS Trần Quang Bảo (Trường Đại học Lâm nghiệp Việt Nam) đã tích hợp Hệ thống Thông tin Địa lý (GIS) và phương pháp nội suy không gian để số hóa, mô hình hóa và trực quan hóa phân bố các chất ô nhiễm.

       [Mạng lưới quan trắc mặt đất]
        (6 vị trí đại diện KK1 - KK6)
  [Phân tích phòng thí nghiệm & Hiện trường]
  - Khí: SO2, NO2, CO, Cl2, H2S, TSP
  - Vi khí hậu & Tiếng ồn (Leq)
       [Hệ thống Thông tin Địa lý (GIS)]
      (ArcGIS 10.8 / ArcPy / PostGIS 3.1)
  - Thuật toán nội suy IDW (Power p=2)
  - Chuyển đổi Vector -> Raster Surface
     [Bản đồ phân vùng ô nhiễm không gian]
  (Đối sánh QCVN 05:2013, QCVN 26:2010)

Mục tiêu dự án

  1. Đánh giá thực trạng ô nhiễm: Định lượng nồng độ các chất ô nhiễm vi khí hậu, khí độc ($\text{SO}_2$, $\text{NO}_2$, $\text{CO}$, $\text{H}_2\text{S}$, $\text{Cl}2$), tổng bụi lơ lửng (TSP) và tiếng ồn ($L{eq}$) tại các nút giao thông, khu dân cư, điểm du lịch và dịch vụ trọng điểm.
  2. Xây dựng cơ sở dữ liệu không gian: Chuẩn hóa tập dữ liệu quan trắc địa lý theo hệ quy chiếu VN-2000 kinh tuyến trục $104^\circ00'$, múi chiếu $3^\circ$.
  3. Mô hình hóa biến động: Ứng dụng kỹ thuật nội suy không gian để thành lập bản đồ chuyên đề về mức độ lan truyền ô nhiễm không khí và tiếng ồn qua các giai đoạn 2014 và 2018.
  4. Đề xuất giải pháp quản lý: Thiết lập luận cứ khoa học hỗ trợ các cơ quan quản lý nhà nước (Sở Tài nguyên và Môi trường tỉnh Sơn La) quy hoạch không gian đô thị và kiểm soát khí thải giao thông hiệu quả.

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

  • Phạm vi không gian: Toàn bộ địa giới hành chính 7 phường trung tâm thuộc thành phố Sơn La gồm: Chiềng Lề, Chiềng An, Tô Hiệu, Quyết Thắng, Quyết Tâm, Chiềng Cơi và Chiềng Sinh.
  • Phạm vi thông số: Quan trắc và phân tích 6 điểm đặc trưng (ký hiệu KK1 đến KK6) với các chỉ tiêu: Nhiệt độ, độ ẩm, tốc độ gió, áp suất, tiếng ồn ($L_{eq}$), $\text{CO}$, $\text{SO}_2$, $\text{NO}_2$, $\text{H}_2\text{S}$, $\text{Cl}_2$ và tổng bụi lơ lửng (TSP).
  • Giới hạn kỹ thuật: Dữ liệu quan trắc hiện trường được lấy vào mùa khô (tháng 03/2018), kết hợp số liệu đối chứng tháng 02/2014; chưa tích hợp mô hình phân tán động lực học khí quyển dạng 3D thời gian thực.

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

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

Trước khi triển khai mô hình GIS, các giải pháp quan trắc và đánh giá chất lượng môi trường không khí thường đối mặt với sự đánh đổi lớn giữa chi phí và độ phủ không gian:

Phương pháp Ưu điểm Nhược điểm Chi phí triển khai Độ phủ không gian
Trạm quan trắc tự động liên tục Dữ liệu thời gian thực ($24/7$), độ chính xác cao, chuẩn hóa tự động. Chi phí đầu tư và bảo dưỡng rất lớn ($>2$ tỷ VNĐ/trạm), bán kính đại diện hạn chế. Rất cao Điểm rời rạc (Point-based)
Viễn thám thuần túy (Satellite AOT) Phủ vùng rộng lớn, quan sát đa thời gian (Landsat 8, Sentinel-5P, MODIS). Sai số cao ở cấp độ vi đô thị, bị ảnh hưởng nặng bởi mây che và địa hình núi cao. Thấp - Trung bình Vùng vĩ mô ($1\times1\text{ km}$)
GIS kết hợp quan trắc điểm (IDW/Kriging) Chi phí hợp lý, trực quan hóa bề mặt liên tục, tích hợp linh hoạt đa chỉ tiêu. Phụ thuộc vào mật độ trạm quan trắc mặt đất và cấu hình thuật toán nội suy. Tối ưu Bề mặt đô thị liên tục ($10\times10\text{ m}$)

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must have (Bắt buộc): Module nhập liệu tọa độ trạm quan trắc (WGS84/VN-2000); thuật toán nội suy nghịch đảo khoảng cách có trọng số (IDW); công cụ phân lớp raster theo ngưỡng Quy chuẩn Kỹ thuật Quốc gia (QCVN).
  • Should have (Nên có): Bảng phân tích tương quan biến động chuỗi thời gian 2014 - 2018; module tính toán diện tích vượt quy chuẩn tự động.
  • Could have (Có thể có): Khả năng xuất bản đồ tương thích chuẩn WebGIS qua GeoServer; module nhập liệu tự động từ file Excel/CSV chuẩn hóa.
  • Won't have (Chưa hỗ trợ): Mô phỏng lan truyền khí động lực học 3D (CFD) theo hình thái kiến trúc vi mô; dự báo tự động bằng Machine Learning.
graph TD
    A["Dữ liệu Quan trắc Hiện trường<br/>(KK1 - KK6)"] --> B["Phòng Thí nghiệm Phân tích<br/>(TCVN/QCVN)"]
    C["Dữ liệu Không gian Nền<br/>(Bản đồ Hành chính Sơn La)"] --> D["Hệ cơ sở dữ liệu GIS<br/>(PostGIS / Shapefile)"]
    B --> D
    D --> E["Module Xử lý Không gian<br/>(ArcGIS 10.8 / ArcPy Engine)"]
    E --> F["Thuật toán Nội suy IDW<br/>(Spatial Analyst)"]
    F --> G["Raster Mặt phẳng Ô nhiễm<br/>(Cell size: 10m x 10m)"]
    G --> H["Phân vùng Cảnh báo<br/>(So sánh QCVN 05 & 26)"]
    H --> I["Bản đồ Chuyên đề & Báo cáo Ra quyết định"]

Thiết kế hệ thống và ngăn xếp công nghệ

  • Hệ thống Xử lý Không gian: ESRI ArcGIS Desktop 10.8 (ArcMap, ArcCatalog) kết hợp Extension Spatial Analyst phục vụ nội suy bề mặt.
  • Ngôn ngữ Lập trình & Scripting: Python 3.8 với thư viện ArcPy, GeoPandas 0.10.2, GDAL/OGR 3.4.1, NumPy 1.21.0 để tự động hóa đường ống xử lý dữ liệu raster.
  • Hệ quản trị Cơ sở dữ liệu Không gian: PostgreSQL 13 tích hợp không gian PostGIS 3.1.
  • Hệ quy chiếu chuẩn: VN-2000 UTM Zone 48N (Kinh tuyến trục $104^\circ$, múi chiếu $3^\circ$, Ellipsoid WGS84 biến đổi).
  • Thiết bị đo đạc & phương pháp phân tích:
    • Tiếng ồn: Thiết bị đo mức âm tích phân theo tiêu chuẩn TCVN 5964:1995 và QCVN 26:2010/BTNMT.
    • Bụi TSP: Phương pháp khối lượng sử dụng màng lọc sợi thủy tinh theo TCVN 5067:1995.
    • Khí $\text{SO}_2, \text{NO}_2, \text{CO}$: Phương pháp hấp thụ hóa chất kết hợp quang phổ so màu (TCVN 5971:1995, TCVN 6137:2009, TCVN 5972:1995).

Implementation và kết quả

Quy trình phân tích và thuật toán nội suy không gian

Thuật toán Nghịch đảo Khoảng cách có Trọng số (Inverse Distance Weighting - IDW) được áp dụng để tính toán giá trị nồng độ ô nhiễm tại các điểm chưa quan trắc dựa trên giả định mức độ ảnh hưởng của nguồn phát thải giảm dần theo khoảng cách địa lý.

Công thức toán học của thuật toán IDW:

$$\hat{Z}(s_0) = \sum_{i=1}^{N} \lambda_i Z(s_i)$$

Trong đó trọng số $\lambda_i$ được xác định bởi:

$$\lambda_i = \frac{d_{i0}^{-p}}{\sum_{j=1}^{N} d_{j0}^{-p}} \quad \text{với} \quad \sum_{i=1}^{N} \lambda_i = 1$$

  • $\hat{Z}(s_0)$: Giá trị nội suy nồng độ chất ô nhiễm tại tọa độ $s_0$.
  • $Z(s_i)$: Giá trị nồng độ quan trắc thực tế tại trạm $s_i$.
  • $d_{i0}$: Khoảng cách Euclid giữa điểm quan trắc $s_i$ và điểm cần tính $s_0$.
  • $p$: Bậc lũy thừa trọng số khoảng cách (lựa chọn $p = 2$ cho lan truyền không khí đô thị).

Đoạn mã tự động hóa quy trình nội suy và phân lớp chất lượng không khí bằng ArcPy:

import arcpy
from arcpy.sa import *

# Cấu hình môi trường làm việc
arcpy.env.workspace = "C:/GIS_SonLa_AirQuality/Workspace.gdb"
arcpy.env.overwriteOutput = True
arcpy.CheckOutExtension("Spatial")

input_points = "Monitoring_Stations_2018"
cell_size = 10.0  # Kích thước pixel 10m x 10m
power_val = 2     # Bậc khoảng cách p = 2

pollutants = {
    "TSP": {"field": "TSP_ug_m3", "qcvn": 300.0},
    "NO2": {"field": "NO2_ug_m3", "qcvn": 200.0},
    "SO2": {"field": "SO2_ug_m3", "qcvn": 350.0},
    "Noise": {"field": "Leq_dBA", "qcvn": 70.0}
}

for name, params in pollutants.items():
    # 1. Thực hiện nội suy IDW
    idw_raster = Idw(
        in_point_features=input_points,
        z_value_field=params["field"],
        cell_size=cell_size,
        power=power_val,
        search_radius=RadiusVariable(number_of_points=6)
    )
    idw_raster.save(f"Raster_IDW_{name}")
    
    # 2. Phân vùng vượt ngưỡng quy chuẩn kỹ thuật quốc gia
    exceed_raster = Con(Raster(f"Raster_IDW_{name}") > params["qcvn"], 1, 0)
    exceed_raster.save(f"Raster_Exceed_{name}")
    
    print(f"[SUCCESS] Hoàn thành nội suy và phân vùng cảnh báo cho: {name}")

arcpy.CheckInExtension("Spatial")

Kết quả phân tích hiện trạng và biến động chất lượng không khí

Dữ liệu quan trắc thực tế tại 6 vị trí trọng điểm trên địa bàn thành phố Sơn La (tháng 03/2018) được tổng hợp và đối chiếu trực tiếp với các Quy chuẩn kỹ thuật Quốc gia hiện hành:

Mã điểm Vị trí quan trắc Tọa độ địa lý ($N, E$) Tiếng ồn $L_{eq}$ ($\text{dBA}$) Bụi TSP ($\mu\text{g/m}^3$) $\text{SO}_2$ ($\mu\text{g/m}^3$) $\text{NO}_2$ ($\mu\text{g/m}^3$) $\text{CO}$ ($\mu\text{g/m}^3$)
KK1 Ngã tư Cầu Trắng $21^\circ19'37.50''\text{N}; 103^\circ54'52.09''\text{E}$ 76 915 43 43 14.800
KK2 Ngã ba Quyết Thắng $21^\circ19'13.11''\text{N}; 103^\circ55'17.26''\text{E}$ 75 240 67 56 14.050
KK3 Đồi Khau Cả $21^\circ19'45.19''\text{N}; 103^\circ54'33.82''\text{E}$ 61 305 35 32 14.050
KK4 Ngã tư Xe Khách $21^\circ19'30.76''\text{N}; 103^\circ55'14.48''\text{E}$ 74 315 62 68 14.050
KK5 Cổng BVĐK tỉnh Sơn La $21^\circ20'59.23''\text{N}; 103^\circ54'39.96''\text{E}$ 68 290 45 45 14.050
KK6 Khu vực Bến xe TP $21^\circ18'03.84''\text{N}; 103^\circ56'34.77''\text{E}$ 76 280 78 72 15.600
QCVN Ngưỡng giới hạn (1 giờ) — 70 (QCVN 26) 300 (QCVN 05) 350 (QCVN 05) 200 (QCVN 05) 30.000 (QCVN 05)

Đánh giá chi tiết theo từng chỉ tiêu

  1. Ô nhiễm tiếng ồn ($L_{eq}$):
    • 4/6 vị trí quan trắc vượt ngưỡng QCVN 26:2010/BTNMT đối với khu vực thông thường từ 6h - 21h ($70\text{ dBA}$).
    • Điểm KK1 (Ngã tư Cầu Trắng) và KK6 (Bến xe thành phố) ghi nhận mức ồn cao nhất đạt $76\text{ dBA}$ (vượt $8,57%$), điểm KK2 đạt $75\text{ dBA}$ và KK4 đạt $74\text{ dBA}$. Nguyên nhân xuất phát trực tiếp từ xung đột giao thông và mật độ phương tiện hạng nặng lưu thông qua trục đường chính.
  2. Tổng bụi lơ lửng (TSP):
    • Điểm KK1 ghi nhận nồng độ bụi cực đại lên tới $915,\mu\text{g/m}^3$, vượt 3,05 lần giới hạn cho phép ($300,\mu\text{g/m}^3$).
    • Các điểm KK3 ($305,\mu\text{g/m}^3$) và KK4 ($315,\mu\text{g/m}^3$) đều vượt ngưỡng từ $1,02$ đến $1,05$ lần do tập trung nhiều công trình cải tạo chỉnh trang đô thị và mật độ xe khách liên tỉnh.
  3. Các thông số khí quyển ($\text{SO}_2$, $\text{NO}_2$, $\text{CO}$, $\text{H}_2\text{S}$, $\text{Cl}_2$):
    • Toàn bộ nồng độ các chất độc hại đều nằm dưới ngưỡng giới hạn cho phép của QCVN 05:2013/BTNMT và QCVN 06:2009/BTNMT. Tuy nhiên, giá trị $\text{NO}_2$ và $\text{SO}_2$ tại bến xe thành phố (KK6) cao vượt trội (đạt lần lượt $72,\mu\text{g/m}^3$ và $78,\mu\text{g/m}^3$), phản ánh lượng phát thải khí xả động cơ diesel chưa qua xử lý triệt để.

Phân tích biến động môi trường giai đoạn 2014 – 2018

  • Tiếng ồn: Mức độ ồn trung bình toàn đô thị tăng từ $2\text{ dBA}$ đến $7\text{ dBA}$ so với năm 2014. Tại ngã tư Cầu Trắng, mức ồn tăng từ $71\text{ dBA}$ (2014) lên $76\text{ dBA}$ (2018).
  • Bụi TSP: Sự gia tăng cơ học của phương tiện ô tô cá nhân và xe tải xây dựng khiến nồng độ TSP tại điểm KK1 tăng đột biến từ $305,\mu\text{g/m}^3$ (2014) lên $915,\mu\text{g/m}^3$ (2018), tương đương mức tăng $200%$.

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

  1. Phương pháp luận số hóa tích hợp: Khóa luận đã thiết lập quy trình chuẩn hóa kết hợp giữa điều tra xã hội học nhanh nông thôn (PRA với 30 hộ dân), số liệu quan trắc lý hóa phòng thí nghiệm đạt chuẩn TCVN và mô hình xử lý không gian GIS.
  2. Khắc phục điểm mù dữ liệu quan trắc rời rạc: Chuyển đổi dữ liệu bảng biểu 1D thành các bản đồ phân bố trường nồng độ 2D bề mặt liên tục với độ phân giải không gian $10\times10\text{ m}$, cho phép nhận diện ngay các "điểm nóng" (hotspots) về ô nhiễm bụi và tiếng ồn tại thành phố Sơn La.
  3. Tối ưu hóa chi phí giám sát môi trường: Mô hình cung cấp giải pháp quy hoạch mạng lưới lấy mẫu tối ưu, giúp giảm thiểu đến $65%$ kinh phí vận hành so với phương án lắp đặt trạm đo cố định dày đặc tại địa hình miền núi hiểm trở.
Tiêu chí so sánh Phương pháp truyền thống Phương pháp đề tài (GIS-IDW) Hệ thống quan trắc ngoại vi IoT
Độ bao phủ không gian Cực bộ tại điểm đo Bề mặt không gian liên tục toàn thành phố Mạng lưới điểm lưới dày đặc
Chi phí triển khai Rất cao nếu tăng số trạm Tối ưu, tận dụng trạm hiện hữu Trung bình - Cao
Thời gian mô hình hóa Thủ công (vài ngày) Tự động hóa bằng Python/ArcPy ($<5$ phút) Thời gian thực (Real-time)
Khả năng trực quan hóa Báo cáo bảng số liệu tĩnh Bản đồ phân lớp theo chuẩn QCVN Dashboard Web/Mobile

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

Kịch bản ứng dụng trong quản lý đô thị

  • Quy hoạch mạng lưới phân luồng giao thông: Phân bố tải lượng xe tải hạng nặng và xe khách liên tỉnh tránh trục đường Chu Văn Thịnh, Tô Hiệu và ngã tư Cầu Trắng trong các khung giờ cao điểm ($6\text{h}30 - 8\text{h}30$ và $16\text{h}30 - 18\text{h}30$).
  • Cảnh báo và bảo vệ các khu vực nhạy cảm: Thiết lập hành lang cách ly tiếng ồn xung quanh Bệnh viện Đa khoa tỉnh Sơn La (KK5) và Di tích Lịch sử Nhà tù Sơn La trên Đồi Khau Cả (KK3).
  • Hỗ trợ công tác thẩm định ĐTM: Cung cấp lớp dữ liệu nền tảng phục vụ đánh giá tác động môi trường cho các dự án mở rộng hạ tầng giao thông đô thị loại II của thành phố Sơn La.
       [Giai đoạn 1: Chuẩn bị]
       - Khảo sát & đo đạc 6 trạm quan trắc (KK1-KK6)
       - Xử lý số liệu lý hóa theo TCVN
       [Giai đoạn 2: Số hóa & Xử lý]
       - Xây dựng CSDL không gian GIS (VN-2000)
       - Chạy Script ArcPy nội suy IDW (10m x 10m)
       [Giai đoạn 3: Phân tích & Tích hợp]
       - Đối soát tự động QCVN 05:2013 & QCVN 26:2010
       - Phân tích chuỗi biến động 2014 - 2018
       [Giai đoạn 4: Đóng gói & Chuyển giao]
       - Xuất bản tập bản đồ chuyên đề chất lượng không khí
       - Chuyển giao Sở TN&MT Sơn La phục vụ quy hoạch

Yêu cầu triển khai hệ thống

  • Hạ tầng phần cứng: Máy trạm chuyên dụng xử lý GIS (CPU 8 Cores, RAM 16GB, GPU rời 4GB, ổ cứng SSD NVMe 512GB).
  • Môi trường phần mềm: Hệ điều hành Windows 10/11 Pro hoặc Ubuntu 20.04 LTS; ArcGIS Desktop 10.8 / ArcGIS Pro 3.x hoặc QGIS 3.22 LTR; PostgreSQL 13/PostGIS 3.1.
  • Dữ liệu đầu vào định kỳ: Tập tin bảng tính Excel/CSV chứa tọa độ kinh vi độ (WGS84/VN-2000) và nồng độ các thông số môi trường đã kiểm chuẩn QA/QC.

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

Hạn chế kỹ thuật

  1. Mật độ điểm quan trắc: Do nguồn lực hạn chế, nghiên cứu mới chỉ phân tích tại 6 trạm mặt đất đại diện; khoảng cách giữa các trạm tương đối thưa, có thể tạo ra sai số biên tại các khu vực ngoại vi đô thị.
  2. Tính chu kỳ của số liệu: Dữ liệu phản ánh trạng thái môi trường điển hình vào mùa khô tháng 3; chưa tích hợp chuỗi số liệu quan trắc 4 mùa liên tục trong năm để đánh giá hiện tượng nghịch nhiệt mùa đông.
  3. Mô hình tĩnh: Chưa kết hợp dữ liệu khí tượng bề mặt động (trường gió tức thời, mô hình AERMOD/CALPUFF) để mô phỏng sự khuếch tán chất ô nhiễm theo địa hình thực tế.

Hướng phát triển nâng cao

  • Tích hợp cảm biến IoT: Xây dựng mạng lưới cảm biến không khí chi phí thấp (Low-cost Air Quality Sensors) truyền dữ liệu liên tục qua giao thức MQTT/LoRaWAN về máy chủ trung tâm.
  • Nâng cấp nền tảng WebGIS: Xây dựng cổng thông tin giám sát chất lượng không khí trực tuyến trên nền tảng Leaflet/Mapbox, giúp người dân theo dõi chỉ số AQI (Air Quality Index) theo thời gian thực.
  • Tích hợp ảnh vệ tinh viễn thám độ phân giải cao: Kết hợp dữ liệu ảnh Sentinel-5P (đo mật độ cột khí $\text{NO}_2$, $\text{SO}_2$, $\text{CO}$) và thuật toán Random Forest/XGBoost để nâng cao độ chính xác nội suy không gian.

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

  • Sinh viên ngành Môi trường & Quản lý Đất đai: Tiếp cận quy trình hoàn chỉnh từ khâu thiết kế mạng lưới lấy mẫu, phân tích số liệu theo quy chuẩn TCVN/QCVN đến kỹ thuật xử lý không gian trong GIS.
  • Kỹ sư GIS & Lập trình viên Môi trường: Tham khảo cấu trúc mã nguồn ArcPy tự động hóa đường ống xử lý dữ liệu địa lý từ bảng tính sang raster chuyên đề.
  • Cơ quan Quản lý Nhà nước (Sở Tài nguyên và Môi trường Sơn La): Nhận bàn giao bộ cơ sở dữ liệu số hóa và các lớp bản đồ chuyên đề hỗ trợ xây dựng báo cáo hiện trạng môi trường định kỳ và quy hoạch phân vùng đô thị.
  • Cộng đồng cư dân đô thị: Nâng cao nhận thức về các điểm nóng ô nhiễm bụi và tiếng ồn, từ đó chủ động thực hiện các biện pháp bảo vệ sức khỏe đường hô hấp cá nhân.

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

1. Cần cấu hình hệ thống như thế nào để vận hành quy trình nội suy GIS này?

Hệ thống yêu cầu máy tính chạy hệ điều hành Windows 10/11 64-bit hoặc Linux, vi xử lý tối thiểu 4 nhân, 8GB RAM (khuyến nghị 16GB để xử lý raster mượt mà), cài đặt sẵn bộ phần mềm ArcGIS Desktop 10.8 (có giấy phép Spatial Analyst) hoặc nền tảng mã nguồn mở QGIS 3.22+ kèm môi trường Python 3.8.

2. Thuật toán IDW có ưu điểm gì hơn so với phương pháp Kriging trong đề tài này?

IDW có tính toán trực quan, tốc độ xử lý nhanh và không yêu cầu giả định phân phối chuẩn nghiêm ngặt của tập dữ liệu như Kriging. Với tập mẫu nhỏ (6 điểm quan trắc tập trung tại đô thị), IDW giúp tái tạo bề mặt suy giảm nồng độ ô nhiễm theo khoảng cách một cách ổn định và tránh được hiện tượng ngoại suy sai lệch.

3. Quy trình chuyển đổi tọa độ điểm quan trắc sang hệ quy chiếu VN-2000 được thực hiện ra sao?

Tọa độ thu thập tại hiện trường bằng thiết bị GPS cầm tay (hệ tọa độ WGS84) được chuẩn hóa sang định dạng Decimal Degrees. Sau đó, sử dụng công cụ Project trong ArcGIS với tham số chuyển đổi 7 thông số (Datum transformation) tương ứng với khu vực Sơn La để chuyển sang hệ tọa độ phẳng VN-2000 UTM Zone 48N (Kinh tuyến trục $104^\circ00'$).

4. Tại sao nồng độ bụi TSP tại Ngã tư Cầu Trắng (KK1) lại vượt quy chuẩn cao gấp 3 lần?

Điểm KK1 nằm tại nút giao thông huyết mạch trên trục Quốc lộ 6 cắt ngang trung tâm thành phố. Mật độ phương tiện cơ giới dày đặc, nhiều xe tải chở vật liệu xây dựng lưu thông kết hợp với hoạt động cải tạo mặt đường trong giai đoạn khảo sát là nguyên nhân chính dẫn đến nồng độ bụi TSP đo được lên tới $915,\mu\text{g/m}^3$.

5. Khóa luận có thể mở rộng để giám sát các đô thị khác không?

Quy trình công nghệ và cấu trúc mã nguồn ArcPy được thiết kế hoàn toàn theo dạng module hóa. Người dùng chỉ cần cập nhật lớp ranh giới hành chính và bảng dữ liệu quan trắc mới là có thể tự động xuất bản đồ phân vùng ô nhiễm cho bất kỳ đô thị hay khu công nghiệp nào trên toàn quốc.

Kết luận

Khóa luận tốt nghiệp “Ứng dụng GIS đánh giá chất lượng môi trường không khí tại thành phố Sơn La, tỉnh Sơn La” đã giải quyết trọn vẹn bài toán tích hợp công nghệ địa không gian vào quản lý tài nguyên môi trường. Bằng việc lượng hóa chính xác thực trạng ô nhiễm bụi TSP và tiếng ồn vượt chuẩn tại các nút giao thông trọng điểm, công trình không chỉ cung cấp bức tranh toàn cảnh về sự biến động chất lượng không khí giai đoạn 2014 – 2018 mà còn cung cấp khung công cụ số hóa thiết thực, giúp các nhà quản lý địa phương tối ưu hóa chính sách quy hoạch đô thị hướng tới phát triển bền vững.