Tổng quan nghiên cứu

Ô nhiễm không khí đang là thách thức môi trường cấp bách hàng đầu tại các đô thị lớn trên toàn cầu và Việt Nam. Theo báo cáo từ Cơ quan Năng lượng Quốc tế (IEA), ô nhiễm không khí là nguyên nhân gây ra khoảng 6,5 triệu ca tử vong sớm mỗi năm trên thế giới. Đáng chú ý, ước tính 87% dân số toàn cầu và tới 98% người dân Việt Nam đang sống trong môi trường có nồng độ bụi mịn PM2.5 vượt ngưỡng khuyến cáo của Tổ chức Y tế Thế giới (WHO). Tại Thành phố Hồ Chí Minh, áp lực môi trường gia tăng nhanh chóng cùng tốc độ đô thị hóa, nơi đang quản lý hơn 7,97 triệu phương tiện giao thông (trong đó xe mô tô chiếm đến 92%), gây thiệt hại kinh tế ước tính trên 6,1 tỷ USD mỗi năm, tương đương 13,4% tổng sản phẩm trên địa bàn (GRDP).

Trong đánh giá chất lượng không khí, việc xác định quy luật lan truyền và khuếch tán chất ô nhiễm phụ thuộc chặt chẽ vào ba nhóm yếu tố: nguồn phát thải, địa hình và điều kiện khí tượng. Tuy nhiên, mạng lưới quan trắc khí quyển bề mặt tại Thành phố Hồ Chí Minh còn thưa thớt, các trạm đo tự động chỉ cung cấp thông tin tức thời tại một số nút giao thông trọng điểm như An Sương hay Điện Biên Phủ, chưa đáp ứng được yêu cầu chuỗi dữ liệu khí tượng cao không và dữ liệu mặt đất liên tục theo từng giờ.

Xuất phát từ thực tiễn đó, luận văn thạc sĩ chuyên ngành Quản lý Tài nguyên và Môi trường đã thực hiện đề tài ứng dụng mô hình Nghiên cứu và Dự báo Thời tiết (Weather Research and Forecasting - WRF) nhằm xây dựng bộ dữ liệu khí tượng chuẩn hóa, độ phân giải cao. Nghiên cứu tập trung vào không gian Thành phố Hồ Chí Minh, thực nghiệm chi tiết tại Phường 2, Quận 10 và huyện Nhà Bè trong năm 2017. Kết quả cung cấp bộ thông số khí tượng đầu vào chuẩn xác cho mô hình lan truyền AERMOD, mở ra giải pháp công nghệ hiện đại giúp tối ưu hóa công tác giám sát, đánh giá và quản lý chất lượng môi trường không khí đô thị.

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 trên nền tảng lý thuyết động lực học khí quyển phi thủy tĩnh và cơ chế khuếch tán chất ô nhiễm trong lớp biên hành tinh (PBL). Luận văn tích hợp hai khối lý thuyết trọng tâm:

  1. Lý thuyết động lực học khí quyển số trị (Mô hình WRF): Hệ thống mô hình WRF phiên bản nghiên cứu nâng cao (ARW v3.6) giải hệ phương trình thông lượng Euler đầy đủ trong hệ tọa độ áp suất thủy tĩnh thẳng đứng theo địa hình ($\eta$). Mô hình tích hợp phương pháp tích phân thời gian Runge-Kutta bậc 3 cùng mạng lưới so le không gian Arakawa-C, cho phép tái hiện chân thực các hoàn lưu nhiệt, gió địa phương và biến động xáo trộn nhiệt động lực học từ quy mô khu vực đến vi mô.
  2. Lý thuyết khuếch tán chất ô nhiễm lớp biên (Mô hình AERMOD): Khung lý thuyết của Cơ quan Bảo vệ Môi trường Hoa Kỳ (US EPA) mô tả sự lan truyền các dòng khí thải ổn định trong lớp biên khí quyển, kết hợp giữa mô hình phân tán Gauss cho điều kiện khí quyển ổn định và phân bố phi Gauss cho điều kiện đối lưu bất ổn định.

Các khái niệm then chốt được vận dụng gồm: cấu trúc lưới lồng đa tỷ lệ (nested domain), tham số hóa lớp biên hành tinh theo sơ đồ YSU (Yonsei University), mô hình truyền nhiệt - ẩm mặt đất Unified Noah 4 lớp, và chương trình giao tiếp dữ liệu khí tượng trung gian MMIF (Mesoscale Model Interface Program).

   +--------------------------------------------------------+
   |            Dữ liệu toàn cầu NCEP FNL (1°x1°)           |
   |              (6 giờ/lần, 62 lớp, 352 biến)             |
   +---------------------------+----------------------------+
                               |
                               v
   +--------------------------------------------------------+
   |             Hệ thống tiền xử lý WPS & WRF              |
   |           Lưới lồng 4 cấp: D01 (30km) -> D04 (1.11km)  |
   +---------------------------+----------------------------+
                               |
                               v
   +--------------------------------------------------------+
   |                 Chương trình MMIF                      |
   |    (Chuyển đổi dữ liệu NetCDF sang định dạng AERMOD)    |
   +---------------------------+----------------------------+
                               |
                               v
   +--------------------------------------------------------+
   |              Mô hình lan truyền AERMOD                 |
   |   (Tính toán nồng độ ô nhiễm CO theo từng giờ & 24h)   |
   +--------------------------------------------------------+

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

Luận văn kết hợp đồng bộ các phương pháp nghiên cứu định lượng tiên tiến:

  • Nguồn dữ liệu đầu vào: Dữ liệu điều kiện biên khí quyển toàn cầu FNL từ Trung tâm Dự báo Môi trường Quốc gia Hoa Kỳ (NCEP) với độ phân giải không gian 1° x 1°, cập nhật định kỳ 6 giờ/lần gồm 62 lớp dữ liệu và 352 trường biến số. Dữ liệu địa hình và lớp phủ bề mặt thu thập từ cơ sở dữ liệu SRTM độ phân giải cao.
  • Quy mô và phương pháp chọn mẫu kiểm định: Cỡ mẫu kiểm định bao gồm 48 bản tin khí tượng METAR liên tục trong ngày (tần suất thu thập 30 phút/lần) phát từ trạm quan trắc sân bay quốc tế Tân Sơn Nhất theo tiêu chuẩn QCVN 16:2008/BTNMT. Dữ liệu phát thải giao thông được thu thập qua phương pháp đếm phương tiện mặt cắt ngang trực tiếp tại 5 trục đường huyết mạch thuộc Phường 2, Quận 10 (Hùng Vương, Nguyễn Chí Thanh, Sư Vạn Hạnh, Ngô Gia Tự, Lê Hồng Phong).
  • Lý do lựa chọn phương pháp phân tích: Thiết lập miền tính 4 lưới lồng đồng tâm (D01: 30 km; D02: 10 km; D03: 3,33 km; D04: 1,11 km). Cấu hình này giúp chuyển tiếp trơn tru năng lượng từ hoàn lưu quy mô lớn trên toàn lãnh thổ Việt Nam xuống vùng đô thị Thành phố Hồ Chí Minh. Việc sử dụng công cụ MMIF xử lý file đầu ra NetCDF của WRF thành bộ tham số định dạng bề mặt và cao không của AERMET giúp loại bỏ hoàn toàn các sai số chủ quan khi thiếu trạm thám không tại địa phương.

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

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

  1. Khả năng tái lập trường nhiệt độ bề mặt xuất sắc: Kết quả chạy mô hình WRF cho thấy diễn biến nhiệt độ 2m (T2) bám sát xu thế đo đạc thực tế tại trạm Tân Sơn Nhất. Chỉ số sai số trung bình (Bias) giữa nhiệt độ mô phỏng và số liệu quan trắc thực tế duy trì ở mức thấp, dao động trong khoảng từ 0,5°C đến 1,4°C trong suốt chu kỳ 24 giờ.
  2. Mô phỏng trường gió và hoàn lưu đô thị chân thực: WRF tái hiện chính xác trường vận tốc và hướng gió mùa Tây - Tây Nam đặc trưng của Thành phố Hồ Chí Minh trong mùa mưa (tháng 8/2017). Tốc độ gió mô phỏng ghi nhận mức trung bình từ 2,1 m/s đến 5,3 m/s, chênh lệch vận tốc gió so với quan trắc đạt dưới 0,8 m/s, bắt trọn các pha đổi hướng gió giữa ngày và đêm do hiệu ứng đảo nhiệt đô thị.
  3. Phân bố nồng độ CO phản ánh sát thực tế lưu lượng xe: Kết quả tích hợp chuỗi số liệu khí tượng vào mô hình AERMOD để tính toán nồng độ khí CO tại Phường 2, Quận 10 vào ngày 07/08/2017 cho thấy: nồng độ CO trung bình 1 giờ cực đại đạt khoảng 12,8 mg/m3 tại các nút giao cắt lớn, 100% điểm tính toán đều nằm trong giới hạn cho phép của QCVN 05:2013/BTNMT (ngưỡng 30 mg/m3). Nồng độ CO trung bình 24 giờ dao động trong khoảng 1,5 – 4,2 mg/m3, phân bố tập trung cao dọc theo các tuyến hẻm sâu và trục đường có mật độ trên 8.000 lượt xe/giờ.

Thảo luận kết quả

Độ chính xác cao của bộ dữ liệu khí tượng tạo ra từ WRF bắt nguồn từ việc lựa chọn tối ưu tổ hợp các tham số vật lý: sơ đồ vi vật lý WSM3 kết hợp mô hình lớp biên hành tinh YSU và sơ đồ bức xạ sóng dài RRTM. Cấu hình này giúp mô tả chính xác quá trình trao đổi nhiệt hiện và nhiệt ẩn tại bề mặt đô thị nhiệt đới. Khi so sánh với các nghiên cứu trước đây tại Việt Nam vốn chỉ dừng lại ở bước lưới thô từ 6 km đến 27,5 km, bước lưới 1,11 km của miền D04 trong luận văn này đã tạo ra bước đột phá về độ chi tiết không gian.

Thông số khí tượng / Môi trường Giá trị mô phỏng WRF / AERMOD Giá trị quan trắc thực tế Đánh giá sai số & Quy chuẩn
Nhiệt độ cực đại ngày 07/08 33,2 °C 32,0 °C Độ lệch Bias: +1,2 °C
Vận tốc gió trung bình 3,4 m/s 3,1 m/s Sai số: 0,3 m/s (tương quan > 85%)
Nồng độ CO cực đại 1 giờ 12,8 mg/m³ 11,5 – 13,2 mg/m³ Đạt QCVN 05:2013/BTNMT (30 mg/m³)
Nồng độ CO trung bình 24 giờ 3,8 mg/m³ 3,2 – 4,0 mg/m³ Đạt QCVN 05:2013/BTNMT (5 mg/m³)

Trong biểu đồ hoa gió (wind rose) phân tích đối sánh, tần suất xuất hiện hướng gió Tây Nam từ mô hình đạt tỷ lệ tương đồng trên 88% so với dữ liệu trạm quan trắc Tân Sơn Nhất. Các bản đồ đẳng nồng độ CO dạng 2D (isopleth maps) trích xuất từ AERMOD minh chứng rõ nét ảnh hưởng của hiệu ứng hẻm phố đô thị (street canyon), nơi các khối nhà cao tầng làm giảm vận tốc gió cục bộ xuống dưới 1,0 m/s, khiến chất ô nhiễm tích tụ cục bộ vào các khung giờ cao điểm sáng (7h00 - 8h30) và chiều (17h00 - 18h30).

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

  1. Chuẩn hóa quy trình vận hành chuỗi mô hình WRF-MMIF-AERMOD: Sở Tài nguyên và Môi trường cùng các đơn vị tư vấn cần đưa chuỗi mô hình này vào quy chuẩn thẩm định Đánh giá tác động môi trường (ĐTM). Mục tiêu đạt 100% báo cáo ĐTM cho các khu công nghiệp và dự án giao thông lớn áp dụng mô phỏng khí tượng độ phân giải dưới 2 km, hoàn thành lộ trình ban hành hướng dẫn kỹ thuật trong giai đoạn 2024 – 2025.
  2. Nâng cấp sơ đồ tán che đô thị (Urban Canopy Model - UCM): Đội ngũ nghiên cứu cần tích hợp mô hình phân lớp nhà đô thị (BEP/BEM) vào miền tính D04 của WRF. Mục tiêu giảm sai số nhiệt độ bề mặt xuống dưới 0,8°C và cải thiện độ chính xác trường gió bề mặt thêm 15% trong vòng 12 tháng nghiên cứu tiếp theo.
  3. Thiết lập hệ thống dự báo chất lượng không khí tự động: Ủy ban nhân dân Thành phố Hồ Chí Minh cần cấp ngân sách xây dựng trung tâm tính toán hiệu năng cao để chạy mô hình WRF-CMAQ/AERMOD hàng ngày. Mục tiêu cung cấp bản tin dự báo nồng độ bụi PM2.5, NO2, CO trước 48 giờ cho 22 quận, huyện và thành phố Thủ Đức, đưa vào thử nghiệm từ quý IV/2025.
  4. Tăng cường kiểm soát phát thải nguồn giao thông: Sở Giao thông Vận tải cần đẩy nhanh lộ trình kiểm định khí thải xe mô tô, xe gắn máy đang lưu hành, kết hợp mở rộng mạng lưới giao thông công cộng chạy điện. Mục tiêu giảm thiểu 20 – 25% lượng phát thải khí CO và bụi PM10 từ nguồn giao thông đường bộ tại các quận nội thành đến năm 2030.

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

  • 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 Thành phố Hồ Chí Minh sử dụng tài liệu làm cơ sở khoa học để xây dựng chính sách quản lý chất lượng không khí xung quanh và quy hoạch mạng lưới trạm quan trắc tự động.
  • Chuyên viên tư vấn môi trường và lập báo cáo ĐTM: Kỹ sư môi trường tại các công ty tư vấn có thể trực tiếp ứng dụng quy trình thiết lập lưới lồng WRF và trích xuất tệp dữ liệu khí tượng qua MMIF để mô phỏng phát tán khí thải ống khói công nghiệp hoặc giao thông trên phần mềm AERMOD.
  • Học viên cao học, nghiên cứu sinh và giảng viên: Nhóm nghiên cứu thuộc các trường đại học chuyên ngành Môi trường, Khí tượng Thủy văn, Biến đổi Khí hậu sử dụng luận văn như tài liệu học thuật mẫu mực về phương pháp số trị mô hình hóa khí quyển và kỹ thuật xử lý dữ liệu NetCDF, GRIB2.
  • Kỹ sư khoa học dữ liệu và phần mềm môi trường: Các chuyên gia công nghệ thông tin phát triển nền tảng đô thị thông minh (Smart City) có thể khai thác cấu trúc dữ liệu tự động hóa để tích hợp các chỉ số khí tượng dự báo vào hệ thống điều hành giao thông và cảnh báo sức khỏe cộng đồng.

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

1. Tại sao cần sử dụng mô hình WRF thay vì thu thập trực tiếp số liệu từ các trạm quan trắc mặt đất?
Các trạm quan trắc mặt đất tại Thành phố Hồ Chí Minh có số lượng hạn chế, chi phí đầu tư lớn và chỉ cung cấp thông tin tại một điểm cố định. Trong khi đó, mô hình lan truyền như AERMOD đòi hỏi số liệu khí tượng bề mặt và cao không liên tục 24 giờ/ngày. Mô hình WRF cung cấp dữ liệu phủ kín không gian với độ phân giải cao tới 1,11 km, giải quyết triệt để bài toán thiếu hụt chuỗi dữ liệu đầu vào.

2. Độ chính xác của dữ liệu khí tượng do mô hình WRF mô phỏng đạt mức độ nào?
Dữ liệu mô phỏng từ WRF đã được kiểm định chéo với 48 bản tin METAR thực tế từ trạm khí tượng sân bay Tân Sơn Nhất. Kết quả cho thấy độ lệch nhiệt độ trung bình chỉ dao động từ 0,5°C đến 1,4°C, hướng gió và tốc độ gió bám sát hoàn lưu thực tế với hệ số tương quan trên 85%, đảm bảo độ tin cậy tuyệt đối cho mô hình hóa phân tán ô nhiễm.

3. Chương trình MMIF đóng vai trò gì trong quy trình liên kết WRF và AERMOD?
Mô hình WRF xuất dữ liệu ở định dạng khoa học dạng mảng NetCDF, trong khi AERMOD yêu cầu các tệp dữ liệu khí tượng định dạng bảng chuyên biệt (.SFC và .PFL). Chương trình MMIF thực hiện nhiệm vụ trích xuất, nội suy và chuyển đổi các trường động lực học từ NetCDF sang định dạng tương thích với AERMOD mà không làm mất mát độ chính xác vật lý.

4. Mô hình AERMOD tính toán nồng độ CO dựa trên những nguồn phát thải nào trong nghiên cứu?
Nghiên cứu tập trung vào nguồn thải giao thông đường bộ tại Phường 2, Quận 10. Dữ liệu được tính toán dựa trên số liệu đếm lưu lượng phương tiện thực tế (xe máy, ô tô con, xe buýt, xe tải) nhân với hệ số phát thải riêng của từng loại phương tiện theo định mức tiêu thụ nhiên liệu và vận tốc di chuyển trung bình.

5. Kết quả của nghiên cứu có thể mở rộng ứng dụng cho các chất ô nhiễm khác không?
Hoàn toàn có thể. Bộ dữ liệu khí tượng chuẩn hóa tạo ra từ mô hình WRF là trường khí quyển chung, có thể sử dụng trực tiếp để mô phỏng phát tán cho nhiều chất ô nhiễm khác như bụi mịn PM2.5, PM10, khí SO2, NO2 hoặc ứng dụng vào các mô hình hóa học quang hóa phức tạp hơn như CMAQ và CAMx.

Kết luận

  • Xây dựng thành công chuỗi công nghệ liên hoàn: Luận văn đã ứng dụng xuất sắc mô hình số trị WRF v3.6 với cấu hình 4 lưới lồng đạt độ phân giải cực cao 1,11 km tại Thành phố Hồ Chí Minh.
  • Chuẩn hóa dữ liệu trung gian qua MMIF: Tạo lập thành công bộ dữ liệu khí tượng bề mặt và cao không theo từng giờ, tương thích hoàn toàn với yêu cầu đầu vào nghiêm ngặt của mô hình phân tán AERMOD.
  • Kiểm định thực nghiệm chuẩn xác: Độ tin cậy của mô hình được chứng minh qua số liệu quan trắc sân bay Tân Sơn Nhất với sai số nhiệt độ thấp dưới 1,4°C và tái hiện hoàn hảo hướng gió mùa Tây Nam.
  • Mô phỏng thành công lan truyền ô nhiễm CO đô thị: Đánh giá chi tiết sự phân bố nồng độ CO tại Phường 2, Quận 10, khẳng định vai trò quyết định của khí tượng và địa hình cục bộ đối với sự tích tụ ô nhiễm giao thông.
  • Mở ra hướng nghiên cứu mở rộng: Đặt nền móng vững chắc cho việc ứng dụng chuỗi mô hình WRF-AERMOD trong quản lý môi trường, dự báo ô nhiễm thời gian thực và quy hoạch hạ tầng đô thị bền vững.

Để phát huy tối đa giá trị thực tiễn của công trình, các cơ quan quản lý và viện nghiên cứu cần tiếp tục hợp tác triển khai tự động hóa quy trình dự báo khí tượng - chất lượng không khí trên nền tảng điện toán đám mây trong giai đoạn 2024 – 2026. Hãy khai thác ngay bộ khung phương pháp luận của luận văn này để nâng cao chất lượng các báo cáo đánh giá tác động môi trường và xây dựng các giải pháp kiểm soát ô nhiễm không khí hiệu quả cho đô thị của bạn.