Giới thiệu dự án

Ô nhiễm không khí dạng hạt và các hợp chất hữu cơ độc hại bền vững đang là thách thức môi trường nghiêm trọng tại các siêu đô thị đang phát triển. Theo báo cáo của Tổ chức Y tế Thế giới (WHO), ô nhiễm không khí xung quanh là nguyên nhân dẫn đến hơn 4,2 triệu ca tử vong sớm mỗi năm trên toàn cầu do các bệnh lý tim mạch, đột quỵ, ung thư phổi và nhiễm trùng đường hô hấp cấp tính. Tại Việt Nam, Thành phố Hồ Chí Minh (TP.HCM) nói chung và Thành phố Thủ Đức nói riêng là trung tâm kinh tế - công nghiệp trọng điểm với mật độ dân cư và lưu lượng phương tiện giao thông cơ giới dày đặc. Nồng độ tổng bụi lơ lửng (Total Suspended Particles - TSP) tại các trục giao thông chính thường xuyên vượt ngưỡng khuyến cáo, dao động từ $50 - 90\ \mu\text{g/m}^3$.

Tuy nhiên, phần lớn các chương trình quan trắc hiện nay chỉ dừng lại ở việc đo lường khối lượng bụi tổng hoặc bụi $PM_{2.5}$ đơn thuần mà chưa phân tích sâu thành phần hóa học độc hại hấp phụ trên các phân đoạn kích thước hạt khác nhau, đặc biệt là nhóm Hydrocacbon thơm đa vòng (Polycyclic Aromatic Hydrocarbons - PAHs). PAHs là nhóm hợp chất hữu cơ bán bay hơi có từ 2 đến nhiều vòng benzen ngưng tụ, sinh ra chủ yếu từ quá trình nhiệt phân và đốt cháy không hoàn toàn nhiên liệu hóa thạch, sinh khối và rác thải đô thị. Việc thiếu hụt dữ liệu phân bố kích thước đa tầng (multi-fraction size distribution) và động lực học phân pha Khí - Hạt (Gas-Particle Partitioning) của PAHs tạo nên rào cản lớn trong việc đánh giá chính xác rủi ro ung thư đối với sức khỏe cộng đồng.

       [Phát thải Đô thị]                 [Động lực học Khí quyển]              [Tác động Sức khỏe]
 (Giao thông, Đốt sinh khối)               (Phân chia Pha & Kích thước)          (Hít thở & Tích tụ)

Đồ án tốt nghiệp chuyên ngành Công nghệ Kỹ thuật Môi trường với đề tài "Đánh giá sự phân bố Hydrocacbon thơm đa vòng (PAHs) trong không khí tại Thành phố Thủ Đức" được thực hiện nhằm giải quyết trực tiếp khoảng trống nghiên cứu trên thông qua 4 mục tiêu cụ thể:

  1. Quan trắc và định lượng phân bố nồng độ bụi theo 6 cấp kích thước hạt ($\le 0,1\ \mu\text{m}$; $0,1 - \le 0,5\ \mu\text{m}$; $0,5 - \le 1,0\ \mu\text{m}$; $1,0 - \le 2,5\ \mu\text{m}$; $2,5 - \le 10\ \mu\text{m}$; $> 10\ \mu\text{m}$) trong chu kỳ 1 năm bao gồm cả mùa mưa và mùa khô.
  2. Định danh và phân tích sự phân pha của 27 hợp chất PAHs (bao gồm 16 hợp chất theo chuẩn USEPA và 15 hợp chất theo chuẩn EU) giữa pha khí và pha bụi.
  3. Truy xuất và nhận diện nguồn gốc phát thải PAHs bằng phương pháp Tỷ lệ chẩn đoán (Diagnostic Ratios - DR), Phân tích thành phần chính (Principal Component Analysis - PCA) kết hợp Mô hình quỹ đạo khối khí ngược HYSPLIT (Hybrid Single Particle Lagrangian Integrated Trajectory).
  4. Đánh giá định lượng rủi ro sức khỏe cộng đồng (độc tính tương đương $\text{BaP-TEQ}$ và rủi ro ung thư tích lũy trọn đời - Incremental Lifetime Cancer Risk - ILCR) theo từng nhóm tuổi phơi nhiễm.

Phạm vi nghiên cứu được triển khai tại điểm quan trắc đô thị đặc trưng (sân thượng tòa nhà trung tâm, Trường Đại học Sư phạm Kỹ thuật TP.HCM, TP. Thủ Đức) từ tháng 07/2023 đến tháng 06/2024, tuân thủ nghiêm ngặt phương pháp chuẩn TO-13A của Cơ quan Bảo vệ Môi trường Hoa Kỳ (USEPA).

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

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

Các phương pháp quan trắc chất lượng không khí truyền thống tại Việt Nam hiện bộc lộ nhiều hạn chế về mặt phân giải hóa học và độc học:

Tiêu chí so sánh Trạm quan trắc tự động liên tục Lấy mẫu lọc bụi tổng (High-Vol TSP/PM10) Phương pháp Phân tầng Đa cấp + GC-MS (Đồ án)
Độ phân giải kích thước hạt Chỉ đo $PM_{2.5}$ hoặc $PM_{10}$ tổng 1 cấp hạt duy nhất trên màng lọc 6 cấp phân đoạn ($\le 0,1\ \mu\text{m}$ đến $> 10\ \mu\text{m}$)
Phân tích pha khí Không hỗ trợ hợp chất hữu cơ vết Không thu giữ được pha khí Thu hồi triệt để bằng XAD-2 / PUF
Định danh hợp chất độc hại Chỉ đo khí vô cơ ($SO_2, NO_x, CO, O_3$) Giới hạn ở kim loại nặng hoặc ion vô cơ Định lượng 27 đồng phân PAHs độc tính cao
Đánh giá rủi ro ung thư Dựa trên chỉ số AQI tổng thể Ước lượng thô, thiếu chính xác Chi tiết hóa $\text{BaP-TEQ}$ và ILCR theo nhóm tuổi
Truy xuất nguồn gốc Định tính theo hoa gió tức thời Khó phân tách nguồn giao thông / sinh khối Kết hợp đa mô hình DR + PCA + HYSPLIT

Ưu tiên hóa yêu cầu kỹ thuật của giải pháp quan trắc theo ma trận MoSCoW:

  • Must have: Thu thập mẫu bụi 6 tầng kích thước; thu mẫu pha khí bằng vật liệu hấp phụ polyme xốp XAD-2/PUF; phân tích sắc ký khí khối phổ (GC-MS) đạt độ thu hồi $70 - 130%$ theo TO-13A.
  • Should have: Tính toán hệ số phân pha thực nghiệm $\log K_p$; mô phỏng quỹ đạo ngược HYSPLIT 72 giờ ở các độ cao 100m, 500m, 1000m.
  • Could have: Đối chiếu rủi ro sức khỏe riêng biệt giữa tiêu chuẩn USEPA và tiêu chuẩn Liên minh Châu Âu (EU).
  • Won't have (lần này): Quan trắc trực tuyến thời gian thực (Online real-time MS) do giới hạn chi phí thiết bị.

Thiết kế hệ thống

Quy trình nghiên cứu và kiến trúc xử lý dữ liệu được thiết kế thành một chuỗi khép kín từ hiện trường đến phòng thí nghiệm và mô hình hóa tính toán:

flowchart TD
    subgraph FIELD["1. THU MẪU HIỆN TRƯỜNG (Chu kỳ 1 năm)"]
        A["Không khí xung quanh (TP. Thủ Đức)"] --> B["Cascade Impactor 6 tầng (Pha bụi: <=0.1 đến >10 um)"]
        A --> C["Vật liệu hấp phụ XAD-2 / PUF (Pha khí)"]
        D["Trạm khí tượng (Nhiệt độ, Độ ẩm, Tốc độ gió, Áp suất)"]
    end

    subgraph LAB["2. TIỀN XỬ LÝ & PHÂN TÍCH HÓA HỌC (ISO/USEPA TO-13A)"]
        B --> E["Cân vi lượng 5 số (Xác định dM/dLogDp)"]
        E --> F["Chiết Soxhlet (DCM:Hexane 1:1, 16h)"]
        C --> F
        F --> G["Cô quay chân không & Thổi khí N2 đến 1.0 mL"]
        G --> H["Hệ thống GC-MS (Chế độ SIM, Cột DB-5MS)"]
    end

    subgraph CHEMOMETRICS["3. XỬ LÝ DỮ LIỆU & MÔ HÌNH HÓA"]
        H --> I["Hệ số đáp ứng tương đối (RRF) & Đường chuẩn"]
        I --> J["Định lượng 27 hợp chất PAHs"]
        J --> K["Hệ số phân pha Log Kp vs Áp suất bay hơi Pv"]
        J --> L["Tỷ lệ chẩn đoán (DR) & Ma trận tương quan Spearman"]
        J --> M["Phân tích thành phần chính (PCA Biplot)"]
        D --> N["Mô hình quỹ đạo khí ngược HYSPLIT (NOAA GDAS)"]
    end

    subgraph RISK["4. ĐÁNH GIÁ RỦI RO SỨC KHỎE"]
        K & L & M & N --> O["Độc tính tương đương BaP-TEQ"]
        O --> P["Mô hình ILCR (Trẻ em, Vị thành niên, Người lớn)"]
    end

Cấu hình công nghệ và thông số vận hành thiết bị

  • Thiết bị lấy mẫu phân tầng: Bộ lấy mẫu va đập nhiều tầng (Cascade Impactor) vận hành ở lưu lượng $28,3\ \text{L/phút}$, màng lọc sợi thủy tinh/thạch anh gia nhiệt $450^\circ\text{C}$ trong 4 giờ trước khi dùng.
  • Hệ thống Sắc ký khí Khối phổ: GC-MS (Agilent Technologies) trang bị cột mao quản phân cực trung bình DB-5MS ($30\ \text{m} \times 0,25\ \text{mm} \times 0,25\ \mu\text{m}$).
    • Khí mang: Heli siêu sạch ($99,999%$), tốc độ dòng không đổi $1,0\ \text{mL/phút}$.
    • Chương trình nhiệt độ lò: Giữ $60^\circ\text{C}$ (1 phút) $\rightarrow$ tăng $10^\circ\text{C/phút}$ lên $200^\circ\text{C}$ $\rightarrow$ tăng $5^\circ\text{C/phút}$ lên $300^\circ\text{C}$ (giữ 10 phút).
    • Chế độ ion hóa: Bắn phá điện tử (EI, $70\ \text{eV}$), quét ion chọn lọc (Selected Ion Monitoring - SIM) theo các mảnh $m/z$ định lượng đặc trưng cho từng đồng phân.
  • Công cụ tính toán & Thống kê: Python v3.11 (scipy, scikit-learn, pandas), R v4.3.2 cho phân tích PCA biplot, và NOAA HYSPLIT v5.2 với nguồn cấp dữ liệu khí tượng toàn cầu GDAS $1^\circ \times 1^\circ$.

Methodology

Nghiên cứu áp dụng quy trình kiểm soát chất lượng QA/QC nghiêm ngặt theo khuyến cáo USEPA TO-13A:

  • Kiểm soát mẫu trắng (Blank controls): Mẫu trắng hiện trường (Field blanks) và mẫu trắng phòng thí nghiệm (Method blanks) được phân tích song song cho mỗi đợt 10 mẫu để trừ nền nhiễu.
  • Hiệu chuẩn độ thu hồi (Surrogate Recovery): Thêm chất chuẩn đánh dấu vào tất cả các mẫu trước khi chiết Soxhlet; giới hạn chấp nhận độ thu hồi đạt $78,4% - 112,6%$.
  • Độ tuyến tính đường chuẩn: Hệ số tương quan $R^2 \ge 0,995$ trên khoảng nồng độ 5 điểm; hệ số biến thiên tương đối của hệ số đáp ứng ($\text{RSD of RRF}$) $< 15%$.
  • Giới hạn phát hiện phương pháp (MDL): Dao động từ $0,005 - 0,08\ \text{ng/m}^3$ đối với 27 hợp chất PAHs.

Implementation và kết quả

Development process

Quá trình triển khai thực địa và xử lý dữ liệu được tóm tắt qua các mốc giải thuật chính:

1. Phương trình động lực học phân pha Khí - Hạt

Hệ số phân bố khí - hạt thực nghiệm ($K_p, \text{m}^3/\mu\text{g}$) được xác định thông qua nồng độ PAHs trong pha hạt ($C_p, \text{ng/m}^3$), pha khí ($C_g, \text{ng/m}^3$) và nồng độ bụi tổng ($\text{TSP}, \mu\text{g/m}^3$):

$$K_p = \frac{C_p / \text{TSP}}{C_g} \iff \log K_p = \log\left(\frac{C_p}{C_g}\right) - \log(\text{TSP})$$

Mối quan hệ giữa $\log K_p$ và áp suất hơi bão hòa dưới dạng lỏng lỏng ($\log P_L^\circ$) thể hiện cơ chế hấp phụ bề mặt hoặc hòa tan hữu cơ: $\log K_p = m \log P_L^\circ + b$.

2. Tính toán Đương lượng độc tính ($\text{BaP-TEQ}$) và Rủi ro ung thư ($\text{ILCR}$)

Độ độc tương đương Benzo[a]pyrene ($\text{BaP-TEQ}$) được tính theo hệ số độc tính tương đương ($\text{TEF}_i$):

$$\text{BaP-TEQ} = \sum_{i=1}^{27} \left( C_i \times \text{TEF}_i \right)$$

Rủi ro ung thư trọn đời do phơi nhiễm qua đường hô hấp ($\text{ILCR}_{\text{inh}}$):

$$\text{ILCR}{\text{inh}} = \frac{\text{BaP-TEQ} \times \text{IR} \times \text{EF} \times \text{ED} \times \text{CSF}{\text{inh}}}{\text{BW} \times \text{AT} \times 10^6}$$

Trong đó: $\text{IR}$ là tốc độ hô hấp ($\text{m}^3/\text{ngày}$); $\text{EF}$ là tần suất phơi nhiễm ($350\ \text{ngày/năm}$); $\text{ED}$ là thời gian phơi nhiễm ($6$ năm đối với trẻ em, $30$ năm đối với người lớn); $\text{BW}$ là trọng lượng cơ thể ($15\ \text{kg}$ cho trẻ em, $60\ \text{kg}$ cho người lớn); $\text{AT}$ là thời gian trung bình đời sống ($70 \times 365\ \text{ngày}$); $\text{CSF}_{\text{inh}} = 3,85\ (\text{mg/kg-ngày})^{-1}$ là hệ số độ dốc gây ung thư đường hô hấp theo USEPA.

import numpy as np
import pandas as pd

def calculate_pah_metrics(df_concentrations: pd.DataFrame, tsp_ug_m3: float) -> pd.DataFrame:
    """
    Tính toán BaP-TEQ, Log Kp và Rủi ro Ung thư (ILCR) cho mẫu quan trắc PAHs.
    df_concentrations: DataFrame chứa cột ['PAH_Name', 'C_gas_ng_m3', 'C_part_ng_m3', 'TEF']
    """
    results = df_concentrations.copy()
    
    # 1. Tính tổng nồng độ và tỷ lệ phân pha
    results['C_total_ng_m3'] = results['C_gas_ng_m3'] + results['C_part_ng_m3']
    results['Particulate_Ratio_%'] = (results['C_part_ng_m3'] / results['C_total_ng_m3']) * 100
    
    # 2. Tính hệ số phân pha Log Kp (tránh chia cho 0)
    results['Log_Kp'] = np.where(
        (results['C_gas_ng_m3'] > 0) & (results['C_part_ng_m3'] > 0),
        np.log10(results['C_part_ng_m3'] / results['C_gas_ng_m3']) - np.log10(tsp_ug_m3),
        np.nan
    )
    
    # 3. Tính Đương lượng độc tính BaP-TEQ (ng/m3)
    results['TEQ_gas'] = results['C_gas_ng_m3'] * results['TEF']
    results['TEQ_part'] = results['C_part_ng_m3'] * results['TEF']
    results['TEQ_total'] = results['C_total_ng_m3'] * results['TEF']
    
    return results

def compute_ilcr(total_bap_teq_ng_m3: float, age_group='adult') -> float:
    # Tham số sinh học theo USEPA
    params = {
        'child': {'IR': 7.6, 'BW': 15.0, 'ED': 6.0},
        'adult': {'IR': 16.0, 'BW': 60.0, 'ED': 30.0}
    }
    p = params[age_group]
    EF, AT, CSF = 350.0, 70.0 * 365.0, 3.85
    
    # Đổi TEQ từ ng/m3 sang mg/m3 (* 1e-6)
    teq_mg_m3 = total_bap_teq_ng_m3 * 1e-6
    ilcr = (teq_mg_m3 * p['IR'] * EF * p['ED'] * CSF) / (p['BW'] * AT)
    return ilcr

Testing và validation

Hiệu năng của phương pháp phân tích và độ tin cậy của tập dữ liệu thực nghiệm được kiểm chứng qua các tiêu chuẩn QA/QC:

Thông số kiểm soát Mục tiêu tiêu chuẩn (TO-13A) Kết quả thực tế đạt được Đánh giá
Độ thu hồi mẫu thêm chuẩn (Spike %) $70% - 130%$ $84,2% - 106,8%$ Đạt yêu cầu
Độ lặp lại phân tích sắc ký (RSD %) $< 10%$ $2,1% - 5,8%$ Độ chụm cao
Độ tuyến tính đường chuẩn ($R^2$) $\ge 0,990$ $0,996 - 0,999$ Tuyến tính xuất sắc
Độ biến thiên RRF ($\text{RSD}_{\text{RRF}}$) $< 15%$ $4,3% - 9,7%$ Ổn định ion hóa
Nhiễu nền mẫu trắng hiện trường $< \text{MDL}$ Không phát hiện / $< 0,01\ \text{ng/m}^3$ Không nhiễm bẩn chéo

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

Nghiên cứu cung cấp bức tranh định lượng toàn diện về ô nhiễm bụi và PAHs tại TP. Thủ Đức qua 2 mùa:

Chỉ số quan trắc Mùa mưa (Mean ± SD) Mùa khô (Mean ± SD) Xu hướng biến thiên
Nồng độ TSP ($\mu\text{g/m}^3$) $26,15 \pm 5,75$ (dao động $20,4 - 31,9$) $52,25 \pm 15,35$ (dao động $36,9 - 67,6$) Tăng $100%$ trong mùa khô
$\sum_{27}\text{PAHs}$ tổng ($\text{ng/m}^3$) $180,1 \pm 38,55$ $351,5 \pm 71,69$ Tăng $95,2%$ trong mùa khô
$\sum_{16}\text{USEPA PAHs}$ ($\text{ng/m}^3$) $145,4 \pm 34,82$ ($81%$ tổng PAHs) $293,8 \pm 60,70$ ($84%$ tổng PAHs) Chiếm ưu thế tuyệt đối
$\sum_{15}\text{EU PAHs}$ ($\text{ng/m}^3$) $17,3 \pm 6,31$ $31,4 \pm 14,50$ Tăng $81,5%$ trong mùa khô
$\text{BaP-TEQ}$ tổng ($\text{ng/m}^3$) $24,4 \pm 6,10$ $47,4 \pm 15,70$ Vượt 24 - 47 lần chuẩn EU ($1\ \text{ng/m}^3$) 
    Phân bố kích thước hạt (dM/dLogDp)               Đóng góp Độc tính TEQ theo Pha
    Bụi mịn (0.2 - 2.5 um) : [████████████] 58%       Pha hạt (Particulate): [████████████████] 92%
    Bụi thô (2.5 - 10 um)  : [███████     ] 34%       Pha khí (Gas phase)  : [██              ]  8%
    Bụi siêu mịn (<0.1 um) : [██          ]  8%       (Dù nồng độ khối lượng pha khí chiếm >70%)

Các phát hiện khoa học then chốt:

  1. Quy luật phân bố kích thước hạt: Hạt bụi tồn tại chủ yếu ở 2 chế độ (Bimodal distribution): Chế độ tích tụ hạt mịn ($0,2 - 2,5\ \mu\text{m}$) và chế độ hạt thô ($2,5 - 10\ \mu\text{m}$).
  2. Quy luật phân pha theo khối lượng phân tử:
    • LMW-PAHs (2–3 vòng: NAP, ACY, ACE, FLU, PHE, ANT) có áp suất hơi cao, phân bố chủ yếu trong pha khí ($> 85%$).
    • MMW-PAHs (4 vòng: FL, PYR, BaA, CHR) chuyển tiếp phân bố đồng đều giữa hai pha.
    • HMW-PAHs (5–6 vòng: BbF, BkF, BaP, IcdP, DahA, BghiP) kỵ nước và áp suất hơi cực thấp, liên kết chặt chẽ và tập trung $> 90%$ trong tầng bụi mịn ($\le 2,5\ \mu\text{m}$).
  3. Nhận diện nguồn thải (Diagnostic Ratios & PCA):
    • Tỷ lệ $\text{IcdP} / (\text{IcdP} + \text{BghiP}) = 0,38 - 0,46$ kết hợp $\text{Flu} / (\text{Flu} + \text{Pyr}) = 0,42 - 0,55$ chỉ ra nguồn phát thải hỗn hợp từ động cơ giao thông (xăng & diesel) và hoạt động đốt sinh khối rơm rạ/nhiên liệu rắn.
    • PCA trích xuất 3 thành phần chính giải thích $> 82%$ tổng phương sai, xác nhận đóng góp lần lượt từ giao thông nội đô ($48%$), đốt nông nghiệp/sinh khối ngoại vi ($24%$) và vận chuyển liên vùng qua HYSPLIT ($10%$).

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

  • Đột phá về phân giải kích thước hạt tại TP.HCM: Đây là công trình đầu tiên tại khu vực đô thị TP. Thủ Đức phân tích đồng thời 27 hợp chất PAHs trên chuỗi 6 tầng kích thước hạt bụi chi tiết ($\le 0,1\ \mu\text{m}$ đến $> 10\ \mu\text{m}$) kết hợp pha khí liên tục trong vòng 1 năm.
  • Chứng minh nghịch lý "Nồng độ - Độc tính": Làm sáng tỏ bằng chứng thực nghiệm: Mặc dù pha khí chiếm ưu thế về khối lượng phát thải ($72% - 88%$), nhưng pha hạt mịn chiếm tới $86% - 98%$ tổng đương lượng độc tính gây ung thư ($\text{BaP-TEQ}$) do sự tích tụ chọn lọc của các hợp chất vòng lớn (5–6 vòng).
  • So sánh với các nghiên cứu tiền nhiệm:
Nghiên cứu Địa bàn quan trắc Số lượng PAHs Phân tầng hạt bụi Pha khí Đánh giá phân pha Log Kp
Hien et al. (2007) Hà Nội / TP.HCM 16 USEPA Không (Chỉ đo $PM_{2.5}$) Không Không
Dat et al. (2017) Đông Á / Đông Nam Á 16–20 PAHs 1–2 tầng thô/mịn Có (gián đoạn) Thô sơ
Đồ án này (2024) TP. Thủ Đức, TP.HCM 27 PAHs (EPA + EU) 6 tầng phân đoạn Có (Đồng thời 100%) Chi tiết theo $\log P_L^\circ$
  • Đóng góp học thuật và chính sách: Cung cấp bộ cơ sở dữ liệu hóa học khí quyển chính xác phục vụ xây dựng Tiêu chuẩn Khí thải Đô thị và Quy hoạch Giao thông Xanh tại TP. Thủ Đức giai đoạn 2025–2030.

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

Tình huống ứng dụng thực tế

  1. Hệ thống Cảnh báo Sớm Rủi ro Sức khỏe Môi trường: Tích hợp bộ trọng số $\text{BaP-TEQ}$ thực nghiệm vào các trạm quan trắc $PM_{2.5}$ hiện có của Sở Tài nguyên & Môi trường để quy đổi trực tiếp từ nồng độ bụi ra nguy cơ gây ung thư hô hấp theo thời gian thực.
  2. Kiểm soát Ô nhiễm Giao thông Cục bộ: Định vị chính xác các điểm nóng giao thông có nồng độ HMW-PAHs cao để áp dụng chính sách hạn chế phương tiện sử dụng dầu diesel cũ hoặc điều tiết luồng xe tải nặng trong giờ cao điểm.
       [Mạng lưới Cảm biến PM2.5/PM10]

Chiến lược nhân rộng và Phân tích Chi phí - Lợi ích

  • Khả năng mở rộng: Mô hình phân tích có thể mở rộng áp dụng cho toàn bộ 21 quận/huyện còn lại của TP.HCM và các vùng kinh tế trọng điểm phía Nam (Bình Dương, Đồng Nai).
  • Hiệu quả kinh tế (ROI Sức khỏe): Giảm thiểu $10\ \mu\text{g/m}^3$ nồng độ bụi mịn chứa PAHs giúp giảm $3,51%$ tỷ lệ nhập viện do bệnh hô hấp ở trẻ em dưới 5 tuổi (theo nghiên cứu Luong et al., 2020), tiết kiệm hàng chục tỷ đồng chi phí điều trị y tế công hàng năm.

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

Hạn chế kỹ thuật

  • Tần suất lấy mẫu: Tần suất quan trắc 2 đợt/tháng do chi phí màng lọc thạch anh và vật liệu hấp phụ siêu sạch XAD-2 còn cao.
  • Mạng lưới trạm: Dữ liệu mới chỉ đại diện cho một điểm đô thị đặc thưng tại trung tâm TP. Thủ Đức, chưa bao phủ hết các phân khu công nghiệp chuyên biệt như Khu chế xuất Linh Trung hay Khu công nghệ cao.

Hướng nghiên cứu mở rộng

  • Ứng dụng các thiết bị khối phổ trực tuyến thời gian thực (như SP-AMS hoặc PTR-TOF-MS) để ghi nhận động lực học hình thành PAHs thứ cấp theo từng giờ.
  • Kết hợp Machine Learning (Random Forest / XGBoost) để dự báo nồng độ $\text{BaP-TEQ}$ dựa trên tập dữ liệu khí tượng vệ tinh và mật độ giao thông tức thời.

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

  • Sinh viên & Học viên cao học: Tiếp cận tài liệu mẫu mực về quy trình chuẩn hóa ISO/USEPA TO-13A, kỹ thuật chiết Soxhlet và phương pháp giải phổ GC-MS SIM.
  • Kỹ sư Môi trường & Lập trình viên: Nắm vững thuật toán phân pha khí quyển, công thức tính $\text{BaP-TEQ}$, mã nguồn Python tự động hóa tính toán độc học và kỹ năng chạy mô hình NOAA HYSPLIT.
  • Cơ quan Quản lý Nhà nước (Sở TN&MT): Cơ sở dữ liệu khoa học để hoạch định vùng phát thải thấp (Low Emission Zones) và nâng cấp quy chuẩn chất lượng không khí.
  • Viện nghiên cứu & Nhà khoa học: Dữ liệu chuẩn mực về 27 hợp chất PAHs tại vùng khí hậu nhiệt đới gió mùa Đông Nam Á để so sánh toàn cầu.

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

1. Yêu cầu trang thiết bị tối thiểu để triển khai quy trình phân tích PAHs theo đồ án là gì?

Cần phòng thí nghiệm đạt chuẩn an toàn hóa chất trang bị: Bộ lấy mẫu Cascade Impactor 6 tầng ($28,3\ \text{L/phút}$), hệ chiết Soxhlet cách thủy, máy cô quay chân không, cân vi lượng 5 số lẻ ($0,01\ \text{mg}$), nguồn khí $N_2$ siêu sạch và hệ thống GC-MS có khả năng vận hành chế độ SIM với cột DB-5MS.

2. Tại sao nồng độ PAHs trong mùa khô lại cao gấp gần 2 lần so với mùa mưa?

Do mùa khô có hiện tượng nghịch nhiệt tầng đối lưu làm giảm độ dày lớp biên khí quyển, kết hợp với bức xạ mặt trời gay gắt thúc đẩy quá trình bay hơi từ mặt đường và lượng mưa rửa trôi hạt bụi thấp hơn nhiều so với mùa mưa.

3. Làm thế nào để phân biệt PAHs sinh ra từ xăng xe hay từ việc đốt rác/sinh khối?

Sử dụng các tỷ lệ chẩn đoán đồng phân (Diagnostic Ratios): Tỷ lệ $\text{IcdP}/(\text{IcdP}+\text{BghiP}) < 0,20$ đặc trưng cho xăng dầu chưa cháy; từ $0,20 - 0,50$ là khí xả động cơ đốt trong (xăng/diesel); và $> 0,50$ là dấu hiệu đặc trưng của quá trình đốt than, gỗ và sinh khối nông nghiệp.

4. Chi phí vận hành một đợt quan trắc phân tầng PAHs ước tính bao nhiêu?

Chi phí vật tư tiêu hao (màng lọc sợi thạch anh, dung môi HPLC Grade DCM/Hexane, ống hấp phụ XAD-2, chuẩn nội $16\ \text{PAHs}$ đồng vị bền) dao động từ 15 – 25 triệu VNĐ cho một chuỗi mẫu phân tầng hoàn chỉnh.

5. Ý nghĩa của chỉ số $\text{BaP-TEQ} = 47,4\ \text{ng/m}^3$ trong mùa khô đối với sức khỏe là gì?

Giá trị này vượt hơn 47 lần ngưỡng khuyến nghị của Tiêu chuẩn Không khí Châu Âu ($1,0\ \text{ng BaP/m}^3$). Điều này chỉ ra rằng nguy cơ rủi ro tích lũy ung thư ($\text{ILCR}$) đối với người dân sinh sống lâu năm tại khu vực nghiên cứu vượt ngưỡng an toàn ($10^{-6}$), đòi hỏi các biện pháp bảo vệ hô hấp cá nhân và can thiệp chính sách giảm phát thải cấp bách.

Kết luận

Đồ án tốt nghiệp "Đánh giá sự phân bố Hydrocacbon thơm đa vòng (PAHs) trong không khí tại Thành phố Thủ Đức" đã hoàn thành xuất sắc các mục tiêu nghiên cứu với độ tin cậy khoa học cao:

  • Thiết lập hệ thống dữ liệu thực nghiệm chuẩn mực về nồng độ bụi 6 tầng kích thước và 27 hợp chất PAHs trong suốt 1 năm.
  • Chứng minh rõ ràng quy luật: Pha khí chiếm khối lượng lớn nhưng pha bụi mịn ($\le 2,5\ \mu\text{m}$) nắm giữ hơn $90%$ độc tính gây ung thư ($\text{BaP-TEQ}$).
  • Nhận diện thành công 2 nguồn phát thải chính là khí xả giao thông và đốt sinh khối, đặt nền móng kỹ thuật vững chắc cho công tác quan trắc môi trường và bảo vệ sức khỏe cộng đồng tại các đô thị Việt Nam.