Tổng quan nghiên cứu

Ô nhiễm không khí do ôzôn tầng mặt (O3) đang nổi lên như một thách thức môi trường cấp bách tại các quốc gia đang phát triển. Khác với tầng bình lưu nơi ôzôn đóng vai trò tấm khiên ngăn chặn bức xạ cực tím, ôzôn tại tầng đối lưu là chất ô nhiễm thứ cấp nguy hại và là khí nhà kính mạnh có tiềm năng gây ấm lên toàn cầu cao gấp khoảng 1.000 lần so với khí CO2. Tại Việt Nam, sự tăng trưởng vượt bậc của các phương tiện giao thông cùng các hoạt động công nghiệp hóa đã làm gia tăng đáng kể nồng độ các chất tiền chất như NOx và VOCs, thúc đẩy các phản ứng quang hóa hình thành O3. Tiếp xúc ngắn hạn với nồng độ ôzôn trong khoảng 200 - 500 µg/m³ có thể gây viêm phế quản, tổn thương mô biểu mô phổi và làm gia tăng tỷ lệ nhập viện do bệnh hô hấp từ 3% đến 5%.

Mục tiêu trọng tâm của nghiên cứu là xác định quy luật biến động nồng độ ôzôn tầng mặt theo không gian tại 3 vùng địa hình đặc trưng (đồng bằng, miền núi và ven biển) và theo thời gian (chu trình giờ trong ngày, các ngày trong tuần, các tháng và các mùa trong năm). Phạm vi nghiên cứu được thực hiện tại 3 trạm quan trắc không khí tự động liên tục đại diện cho miền Bắc Việt Nam gồm Hà Nội, Phú Thọ và Quảng Ninh với chuỗi dữ liệu hoàn chỉnh trong suốt 365 ngày của năm 2014.

Kết quả của đề tài mang ý nghĩa khoa học và thực tiễn sâu sắc khi lần đầu tiên làm sáng tỏ mối liên hệ định lượng giữa ôzôn với bức xạ mặt trời, tỷ số NO2/NO và SO2 tại miền Bắc. Nghiên cứu cung cấp cơ sở dữ liệu xác thực phục vụ kiểm soát ô nhiễm, cảnh báo rủi ro sức khỏe cộng đồng và đề xuất giải pháp tích hợp thông số O3 vào mạng lưới quan trắc môi trường quốc gia.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên 3 nền tảng lý thuyết cốt lõi của hóa học khí quyển và môi trường:

  • Lý thuyết sương khói quang hóa (Photochemical Smog Theory): Giải thích sự hình thành ôzôn tầng đối lưu (khoảng 9% tổng lượng ôzôn khí quyển) từ phản ứng quang hóa giữa các chất tiền thân dưới tác động của năng lượng bức xạ mặt trời.
  • Chu trình phản ứng quang hóa NOx và gốc tự do ROx: Mô tả cơ chế quang phân NO2 ở bước sóng ánh sáng dưới 424 nm tạo ra nguyên tử oxy hoạt tính (O3P), kết hợp với O2 tạo thành O3. Chu trình này liên tục tái tạo NO và NO2 thông qua sự tham gia của các gốc oxy hóa trung gian như HO2 và RO2 từ các hợp chất hữu cơ dễ bay hơi VOCs.
  • Lý thuyết động lực học lớp biên hành tinh (Planetary Boundary Layer): Làm rõ sự phân tán, tích tụ và trao đổi ôzôn theo chiều thẳng đứng dưới tác động của hoàn lưu gió mùa, hiện tượng nghịch nhiệt và các yếu tố khí tượng vi mô.

Các khái niệm chuyên ngành then chốt được chuẩn hóa bao gồm: ôzôn tầng mặt (tropospheric ozone), tiền chất quang hóa (photochemical precursors), tỷ số chuyển hóa NO2/NO, và nồng độ giới hạn bảo vệ sức khỏe theo quy chuẩn QCVN 05:2013/BTNMT cùng khuyến cáo của Tổ chức Y tế Thế giới (WHO).

Phương pháp nghiên cứu

Nguồn dữ liệu sơ cấp được trích xuất từ 3 trạm quan trắc không khí tự động, liên tục do Trung tâm Quan trắc môi trường - Tổng cục Môi trường quản lý và vận hành. Các trạm đại diện cho 3 hình thái không gian: Trạm Nguyễn Văn Cừ (Hà Nội) đặc trưng cho đô thị đồng bằng với mật độ giao thông khoảng 8.000 xe/giờ; Trạm Vân Cơ (Phú Thọ) đại diện cho vùng trung du bán sơn địa; Trạm Hồng Hà (Quảng Ninh) đại diện cho khu vực ven biển gắn liền hoạt động công nghiệp khai khoáng và nhiệt điện.

Tổng cỡ mẫu nghiên cứu đạt hơn 315.000 bản ghi dữ liệu đo lường liên tục với tần suất lấy mẫu chuẩn 5 phút/lần trong toàn bộ 365 ngày của năm 2014. Tỷ lệ thiếu hụt dữ liệu được kiểm soát chặt chẽ ở mức thấp (Hà Nội thiếu 11%, Quảng Ninh thiếu 13%, Phú Thọ thiếu 9%), đảm bảo độ tin cậy vượt tiêu chuẩn 80% theo yêu cầu quan trắc quốc tế. Phương pháp chọn mẫu là quan trắc tự động liên tục tại các điểm cố định theo phương pháp trắc quang tia cực tím.

Quy trình phân tích dữ liệu được thực hiện đa tầng nhằm tối ưu hóa độ chính xác:

  • Ngôn ngữ lập trình R và gói package Openair: Được lựa chọn nhờ khả năng xử lý chuỗi dữ liệu thời gian lớn, phân tích xu hướng biến động và vẽ biểu đồ hoa gió quang hóa chuyên sâu.
  • Phần mềm Microsoft Excel: Quản lý, chuẩn hóa định dạng dữ liệu bảng .CSV và phân loại dữ liệu theo các khung thời gian trung bình giờ, ngày, tuần, tháng.
  • Phần mềm IBM SPSS: Thực hiện các phép kiểm định thống kê mô tả, ma trận tương quan Pearson, mô hình hồi quy tuyến tính đa biến và phân tích thành phần chính (PCA) nhằm cô lập các cụm nhân tố chi phối nồng độ O3.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phân tích chuỗi số liệu quan trắc năm 2014 đã chỉ ra 4 phát hiện quan trọng về hiện trạng ô nhiễm ôzôn tầng mặt:

Thứ nhất, nồng độ O3 có sự phân hóa không gian sâu sắc. Nồng độ trung bình 24 giờ tại Hà Nội đạt mức cao nhất với 53,5 µg/m³, cao gấp 1,96 lần so với Quảng Ninh (27,3 µg/m³) và cao gấp 2,72 lần so với Phú Thọ (19,7 µg/m³). Nồng độ cực đại theo giờ tại Hà Nội chạm mức 490,4 µg/m³, Quảng Ninh ghi nhận 306,5 µg/m³, trong khi Phú Thọ chỉ đạt mức cao nhất là 143,6 µg/m³.

Thứ hai, chu trình biến thiên ngày đêm tuân theo quy luật quang hóa 4 giai đoạn rõ rệt. Nồng độ O3 ban ngày (07:00 - 18:00) dao động từ 28,4 đến 73,3 µg/m³, cao gấp 2 đến 3 lần so với ban đêm (11,0 - 33,8 µg/m³). Giá trị cực tiểu xuất hiện vào khoảng 05:00 - 06:00 sáng, sau đó tăng vọt và đạt đỉnh cực đại vào thời điểm 12:00 - 14:00 chiều. Tại mốc 12:00 trưa, nồng độ O3 tại Hà Nội đạt 115,4 µg/m³, cao hơn 109% so với Quảng Ninh (55,1 µg/m³) và cao hơn 211% so với Phú Thọ (37,1 µg/m³).

Thứ ba, mức độ vượt quy chuẩn QCVN 05:2013/BTNMT (giới hạn 1 giờ là 200 µg/m³) diễn ra đáng lo ngại tại khu vực đô thị. Hà Nội ghi nhận tới 231 lần vượt ngưỡng (chiếm 2,36% tổng số phép đo trong năm, phân bổ trên 38 ngày), tập trung vào các tháng 1, tháng 10 và tháng 12. Quảng Ninh ghi nhận 87 lần vượt ngưỡng (0,99% số phép đo trên 24 ngày). Ngược lại, Phú Thọ có 0 lần vượt chuẩn (0%).

Thứ tư, hiệu ứng ngày trong tuần thể hiện rõ nét. Nồng độ các chất ô nhiễm và O3 vào các ngày đầu tuần (thứ Hai đến thứ Năm) cao hơn hẳn cuối tuần. Nồng độ O3 trung bình ngày thứ Hai tại Hà Nội đạt đỉnh 94,1 µg/m³, trong khi ngày Chủ nhật giảm 46,5% xuống còn 50,3 µg/m³ do lưu lượng phương tiện giảm sút.

Thảo luận kết quả

Sự gia tăng nồng độ O3 vào đầu giờ chiều bắt nguồn trực tiếp từ cường độ bức xạ mặt trời cực đại, thúc đẩy tốc độ quang phân NO2. Kết quả hồi quy tuyến tính cho thấy tương quan thuận chặt chẽ giữa nồng độ O3 và bức xạ nhiệt tại cả 3 trạm quan trắc.

Khi xem xét tương quan tiền chất, tỷ số nồng độ NO2/NO giải thích được tới 57% sự biến thiên của O3 tại Hà Nội (R² = 0,57) và 45% tại Phú Thọ (R² = 0,45). Điều này phản ánh rõ cơ chế: khi bức xạ tăng, NO2 phân giải giải phóng O3; khi ánh sáng yếu, O3 phản ứng ngược lại với NO tạo thành NO2. Ngược lại, tại Quảng Ninh, tỷ số NO2/NO chỉ giải thích được 10% biến động O3, nhưng nồng độ SO2 lại giải thích được tới 63% sự thay đổi O3 (R² = 0,63). Đây là hiện tượng đặc thù phản ánh tác động tổng hợp của các nguồn thải công nghiệp nhiệt điện và đốt than quy mô lớn tại địa phương.

Về mặt trực quan hóa, toàn bộ mối quan hệ này được chứng minh qua các đồ thị phân tán hồi quy (Scatter Plot) và biểu đồ chuỗi thời gian (Time-series) biểu diễn sự dịch chuyển nồng độ theo chuỗi NO → NO2 → O3. So sánh với quốc tế, mức nồng độ đỉnh 490,4 µg/m³ tại Hà Nội đã vượt gấp 2,45 lần quy chuẩn kỹ thuật quốc gia và cao hơn gấp 2 lần mức cảnh báo nguy hại cho sức khỏe của WHO (240 µg/m³), phản ánh tình trạng ô nhiễm quang hóa tương đương với nhiều siêu đô thị tại châu Á như Bangkok hay Bắc Kinh trong các giai đoạn cao điểm.

Đề xuất và khuyến nghị

Dựa trên các bằng chứng khoa học thu được, 4 nhóm giải pháp chiến lược được đề xuất nhằm kiểm soát ô nhiễm ôzôn tầng mặt:

  • Mở rộng và hiện đại hóa mạng lưới trạm quan trắc tự động: Cục Kiểm soát ô nhiễm môi trường phối hợp cùng các địa phương triển khai lắp đặt thêm tối thiểu 15 trạm quan trắc không khí tự động đồng bộ tại các vùng kinh tế trọng điểm Bắc Bộ trước năm 2025. Mục tiêu nâng tỷ lệ thu nhận dữ liệu đạt chuẩn trên 95%, tích hợp đồng thời cảm biến đo bức xạ tử ngoại UV và nồng độ VOCs để phục vụ mô hình hóa dự báo.
  • Kiểm soát nghiêm ngặt nguồn phát thải tiền chất từ giao thông đô thị: Bộ Giao thông Vận tải và Ủy ban nhân dân thành phố Hà Nội cần đẩy mạnh thực thi quy chuẩn khí thải mức Euro 5 cho 100% xe cơ giới sản xuất mới từ năm 2024, đồng thời thiết lập các vùng phát thải thấp (LEZ) tại khu vực nội đô nhằm cắt giảm tối thiểu 25% lượng phát thải NOx trong khung giờ cao điểm từ 07:00 đến 09:00 sáng.
  • Tăng cường kiểm soát và xử lý khí thải công nghiệp, nhiệt điện: Sở Tài nguyên và Môi trường tỉnh Quảng Ninh chủ trì thanh kiểm tra định kỳ, yêu cầu 100% nhà máy nhiệt điện than và cơ sở đốt than lắp đặt hệ thống khử lưu huỳnh (FGD) và khử nitơ oxit (SCR) trước năm 2026, hướng tới mục tiêu giảm 30% phát thải khí SO2 và NOx cục bộ.
  • Thiết lập hệ thống cảnh báo sớm rủi ro quang hóa cho cộng đồng: Trung tâm Quan trắc môi trường quốc gia hoàn thiện thuật toán dự báo chỉ số AQI ôzôn theo thời gian thực trước năm 2025, tự động phát tín hiệu khuyến cáo người dân hạn chế các hoạt động thể chất ngoài trời trong khung giờ 12:00 - 15:00 khi nồng độ O3 dự báo vượt ngưỡng 160 µg/m³.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn đa dạng cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước về môi trường: Cán bộ tại Cục Kiểm soát ô nhiễm môi trường, Sở Tài nguyên và Môi trường các tỉnh Hà Nội, Quảng Ninh, Phú Thọ có thể sử dụng số liệu 231 lần vượt ngưỡng và mô hình tương quan tiền chất làm căn cứ khoa học để xây dựng Kế hoạch quản lý chất lượng không khí cấp tỉnh giai đoạn 2025 - 2030.
  • Nhà nghiên cứu và giảng viên chuyên ngành môi trường: Các nhà khoa học khí quyển có thể kế thừa bộ cơ sở dữ liệu hơn 315.000 bản ghi cùng quy trình phân tích dữ liệu chuyên sâu bằng gói Openair trên nền tảng R để phát triển các mô hình dự báo ô nhiễm thứ cấp mở rộng.
  • Kỹ sư quy hoạch giao thông và đô thị: Các chuyên gia quy hoạch có thể tham khảo mối liên hệ giữa mật độ 8.000 phương tiện/giờ với nồng độ phát thải NOx để thiết kế các phương án phân luồng giao thông hợp lý, giảm thiểu tắc nghẽn vào các khung giờ phát sinh phản ứng quang hóa mạnh.
  • Tổ chức y tế công cộng và chuyên gia dịch tễ học: Các bác sĩ và nhà nghiên cứu y tế có thể sử dụng các mốc vượt ngưỡng phơi nhiễm ôzôn trên 100 - 200 µg/m³ để đánh giá mối tương quan với các đợt bùng phát bệnh lý hô hấp ở trẻ em và người cao tuổi, từ đó thiết lập quy trình can thiệp y tế dự phòng hiệu quả.

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

1. Tại sao ôzôn ở tầng bình lưu có lợi nhưng ôzôn tầng mặt lại gây độc hại?
Ở độ cao 20 - 50 km, tầng ôzôn hấp thụ tia cực tím độc hại từ mặt trời để bảo vệ sự sống. Tuy nhiên, ở tầng đối lưu sát mặt đất (0 - 10 km), O3 là chất oxy hóa cực mạnh. Khi tiếp xúc với nồng độ trên 200 µg/m³, O3 gây viêm niêm mạc hô hấp, suy giảm chức năng phổi và là tác nhân gây hiệu ứng nhà kính mạnh gấp khoảng 1.000 lần CO2.

2. Tại sao nồng độ ôzôn tầng mặt luôn đạt đỉnh vào khoảng 12:00 đến 14:00 chiều?
Nồng độ O3 phụ thuộc trực tiếp vào cường độ bức xạ mặt trời để kích hoạt phản ứng quang phân NO2 ở bước sóng dưới 424 nm. Vào thời điểm trưa từ 12:00 đến 14:00, bức xạ tổng đạt mức cao nhất trong ngày, thúc đẩy quá trình tạo ôzôn diễn ra mãnh liệt nhất, đưa nồng độ O3 tại Hà Nội lên mức đỉnh 115,4 µg/m³.

3. Nguyên nhân khiến nồng độ ôzôn tại Hà Nội cao vượt trội so với Quảng Ninh và Phú Thọ là gì?
Hà Nội có mật độ giao thông rất cao đạt khoảng 8.000 xe/giờ, phát thải lượng lớn các tiền chất như NO, NO2 và CO (nồng độ NOx trung bình ngày đầu tuần lên tới 96,5 µg/m³). Sự tích tụ đậm đặc của các tiền chất này kết hợp với hiệu ứng đô thị đã tạo điều kiện tối ưu cho các phản ứng quang hóa tạo O3 diễn ra mạnh mẽ hơn 2 đến 3 lần các khu vực khác.

4. Hiệu ứng ngày trong tuần tác động thế nào đến nồng độ ôzôn?
Các hoạt động kinh tế, giao thông và công nghiệp diễn ra sôi động từ thứ Hai đến thứ Sáu làm tăng mạnh lượng phát thải tiền chất. Kết quả quan trắc cho thấy nồng độ O3 trung bình ngày thứ Hai tại Hà Nội đạt 94,1 µg/m³, cao hơn tới 87,1% so với ngày Chủ nhật (chỉ đạt 50,3 µg/m³) khi lưu lượng xe cộ giảm rõ rệt.

5. Khí SO2 đóng vai trò gì trong sự hình thành ôzôn tại khu vực Quảng Ninh?
Tại Quảng Ninh, nồng độ SO2 duy trì ở mức cao trung bình khoảng 31,2 µg/m³ do các hoạt động đốt than và nhiệt điện. Phân tích hồi quy cho thấy SO2 giải thích tới 63% sự biến thiên nồng độ O3 tại đây, chứng minh ảnh hưởng phức hợp của các nguồn thải công nghiệp đặc thù vùng mỏ đến hóa học khí quyển địa phương.

Kết luận

  • Luận văn đã thiết lập bức tranh toàn diện đầu tiên về quy luật biến động nồng độ O3 tầng mặt tại 3 vùng địa hình đặc trưng miền Bắc thông qua chuỗi dữ liệu chuẩn hóa hơn 315.000 bản ghi trong năm 2014.
  • Xác định sự phân hóa không gian rõ rệt, trong đó Hà Nội là điểm nóng ô nhiễm với nồng độ đỉnh 490,4 µg/m³ và 231 lần vượt quy chuẩn QCVN 05:2013/BTNMT, cao vượt trội so với Quảng Ninh (87 lần) và Phú Thọ (0 lần).
  • Làm sáng tỏ cơ chế quang hóa tự nhiên với chu trình đỉnh - đáy ngày đêm và chứng minh định lượng vai trò của tỷ số NO2/NO tại Hà Nội (R² = 0,57) cùng tác động đặc thù của khí SO2 tại Quảng Ninh (R² = 0,63).
  • Cung cấp cơ sở khoa học chuẩn xác để hoàn thiện quy chuẩn kỹ thuật quốc gia và đề xuất các giải pháp kỹ thuật kiểm soát phát thải tiền chất giao thông, công nghiệp với lộ trình cụ thể đến năm 2026.
  • Mở ra hướng tiếp cận phân tích dữ liệu môi trường hiện đại bằng ngôn ngữ lập trình R và SPSS, đóng góp quan trọng vào tiến trình lồng ghép thông số O3 vào hệ thống quan trắc tự động quốc gia.

Các cơ quan quản lý và nhà nghiên cứu cần khẩn trương phối hợp ứng dụng các kết quả này vào thực tiễn quản lý đô thị trước năm 2025 nhằm bảo vệ sức khỏe cộng đồng và xây dựng môi trường sống bền vững.