Tổng quan nghiên cứu

Quá trình công nghiệp hóa và đô thị hóa nhanh chóng tại các vùng kinh tế trọng điểm phía Nam, đặc biệt là Thành phố Hồ Chí Minh, đã tạo ra áp lực môi trường gay gắt khi hơn 70% lượng khí thải độc hại phát sinh từ các cơ sở sản xuất và mạng lưới giao thông. Vấn đề cốt lõi đặt ra cho công tác quản lý đô thị là sự thiếu hụt các công cụ tin học tự động hóa có khả năng mô phỏng và dự báo chính xác động lực học lan truyền chất ô nhiễm theo không gian và thời gian. Cơ chế cấp phép xả thải truyền thống dựa trên các chỉ số trung bình hóa tĩnh không thể ứng phó kịp thời trước những biến động đột xuất của khí tượng thủy văn hoặc các tai biến môi trường cục bộ.

Mục tiêu trọng tâm của nghiên cứu là thiết kế, xây dựng và hoàn thiện phần mềm mô phỏng sự nhiễm bẩn không khí (phần mềm CAP) có khả năng dự báo trường nồng độ chất ô nhiễm từ nhiều nguồn thải công nghiệp trên một vùng lãnh thổ xác định. Phạm vi nghiên cứu bao quát từ quy mô vi mô địa phương (diện tích dưới 10 km) đến quy mô trung bình của vùng đô thị (100 km đến 1.000 km), với thang thời gian dự báo linh hoạt từ vài phút tức thời đến chu kỳ 24 giờ và dài hạn một năm. Nghiên cứu mang ý nghĩa chiến lược sâu sắc trong việc hỗ trợ các nhà quản lý đưa ra quyết định nhanh chóng, giảm thiểu sai số dự báo phát tán xuống dưới 15% và nâng cao hiệu suất xử lý dữ liệu khí tượng môi trường gấp 3 đến 5 lần so với các phương pháp tính toán thủ công trước đây.

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 xây dựng trên nền tảng thủy nhiệt động lực học khí quyển kết hợp hệ phương trình vi phân đạo hàm riêng bảo toàn khối lượng ba chiều cho các chất ô nhiễm. Khung lý thuyết tích hợp ba mô hình nền tảng: lý thuyết khuếch tán rối bán thực nghiệm của Monin-Obukhov cho lớp biên khí quyển gần mặt đất, quy luật vệt nâng dòng khí nóng hai phần ba của Briggs kết hợp mô hình vệt khói Gauss, và hệ phương trình lan truyền chất ô nhiễm Navier-Stokes mở rộng theo trường phái của Viện sĩ Marchuk.

Hệ thống khái niệm chính được chuẩn hóa trong nghiên cứu bao gồm:

  • Độ cao hữu dụng của nguồn thải: Tổng độ cao hình học của ống khói và độ nâng nhiệt động của vệt khói phụ thuộc vào hệ số nổi và vận tốc gió.
  • Tenxơ hệ số trao đổi rối: Các hệ số khuếch tán không gian theo phương ngang và phương đứng (thứ nguyên mét vuông trên giây), phản ánh sự xáo trộn của các xoáy khí quyển.
  • Độ ổn định tầng kết khí quyển: Phân loại trạng thái khí quyển từ bất ổn định mạnh đến ổn định cao theo thang tham số nổi và gradien nhiệt độ.
  • Nghịch đảo nhiệt độ: Hiện tượng nhiệt độ tăng theo độ cao làm suy giảm sự đối lưu thẳng đứng, giữ bẫy chất ô nhiễm tại tầng sát mặt đất ở độ cao dưới 300 m.

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

Nghiên cứu sử dụng cơ sở dữ liệu không gian lưới 3 chiều với quy mô tính toán tiêu chuẩn gồm 20 x 20 x 10 điểm (tương đương 4.000 nút lưới) và mở rộng đến 4.420 điểm lưới đối với các bài toán vệt khói mở rộng. Chuỗi số liệu đầu vào tích hợp dữ liệu trắc quan thực địa về vận tốc gió nền (từ 1 m/s đến 12 m/s), nhiệt độ khí thải từ 80°C đến 180°C, cùng nồng độ đo đạc của các chất chỉ thị như lưu huỳnh đioxit, oxit cacbon và bụi lơ lửng.

Phương pháp chọn mẫu được tiến hành theo mạng lưới không gian đa tỷ lệ đại diện cho ba loại nguồn chính: nguồn điểm (ống khói nhà máy riêng lẻ), nguồn đường (trục giao thông chính) và nguồn diện (khu công nghiệp tập trung). Mật độ điểm mẫu được tối ưu hóa dày hơn tại các khu vực có gradien nồng độ cao hoặc địa hình phức tạp.

Lý do lựa chọn phương pháp sai phân hữu hạn phân rã (splitting method theo sơ đồ Yanenko) kết hợp phương pháp hạt lớn (Particle-in-cell) và kỹ thuật hàm liên hợp xuất phát từ yêu cầu khắc phục triệt để hiện tượng khuếch tán nhân tạo do sai số xấp xỉ số hạng bình lưu (khử độ nhớt giả). Cách tiếp cận này giúp tiết kiệm hơn 40% dung lượng bộ nhớ máy tính, đảm bảo tính ổn định số học và cho phép xử lý bài toán trên máy tính cá nhân trong khoảng thời gian từ năm 1996 đến năm 1998 dưới sự chỉ đạo của Chương trình Công nghệ Thông tin Thành phố Hồ Chí Minh.

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

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

Nghiên cứu đã đạt được các kết quả định lượng cụ thể thông qua việc mô phỏng thực nghiệm:

  • Tác động của nghịch đảo nhiệt: Khi xuất hiện lớp nghịch đảo nhiệt ở độ cao dưới 300 m, nồng độ cực đại của chất ô nhiễm tại mặt đất tăng vọt lên gấp 2,0 lần so với điều kiện khí quyển trung hòa, đồng thời vùng ảnh hưởng nguy hiểm bị nén lại trong bán kính từ 1,5 km đến 3,0 km quanh nguồn thải.
  • Ảnh hưởng của trường gió và độ nhám địa hình: Khi vận tốc gió mặt đất giảm xuống dưới 1,5 m/s, nồng độ tích tụ chất ô nhiễm tại các vùng trũng đô thị tăng thêm khoảng 35% đến 48% do hiệu ứng suy giảm đối lưu ngang.
  • Tối ưu hóa thuật toán số trị: Việc ứng dụng tham số hiệu chỉnh thời điểm phân bố nồng độ bậc một và bậc hai giúp triệt tiêu hơn 85% sai số nhớt giả trong sơ đồ sai phân, duy trì định luật bảo toàn khối lượng khi bước thời gian biến thiên từ 1 phút đến 30 phút.
  • Khả năng tích hợp đa nguồn: Phần mềm CAP mô phỏng đồng thời trên 50 nguồn thải công nghiệp xen kẽ trong khu vực đô thị, tách bạch mức độ đóng góp nồng độ của từng ống khói với độ tương đồng đạt trên 90% so với dữ liệu quan trắc thực tế.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự tập trung nồng độ cao tại mặt đất là do sự tương tác giữa hiệu ứng đảo nhiệt đô thị và các luồng hoàn lưu xoáy cục bộ dọc theo các tuyến đường phố. Về mặt biểu diễn dữ liệu, toàn bộ trường nồng độ ô nhiễm được trực quan hóa thông qua hệ thống bản đồ số với các đường đồng mức nồng độ 2D và bề mặt phân bố 3D liên kết trực tiếp với hệ cơ sở dữ liệu địa lý môi trường GIS. Các bảng tổng hợp ma trận nhạy cảm tham số cho phép so sánh định lượng mức độ tác động của từng yếu tố khí tượng (nhiệt độ, gió, độ ẩm) đến sai số nồng độ.

So sánh với các mô hình hộp đơn giản (Box model) vốn giả định chất ô nhiễm pha trộn đồng nhất tức thời trên toàn thành phố, mô hình trong luận văn thể hiện sự vượt trội khi khắc phục được tình trạng ước tính thấp hơn thực tế tới 50% ở các khu vực lân cận nguồn thải lớn. Kết quả nghiên cứu hoàn toàn tương thích với các công trình thực nghiệm quốc tế của Briggs và Marchuk nhưng được tinh chỉnh phù hợp với đặc thù khí hậu nhiệt đới gió mùa tại Việt Nam.

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

  • Tích hợp module tính toán của phần mềm CAP vào mạng lưới trạm quan trắc tự động liên tục của các cơ quan quản lý môi trường, đặt mục tiêu nâng tỷ lệ tự động hóa giám sát khí thải lên 100% tại các khu công nghiệp trọng điểm trong lộ trình 12 tháng, do các Sở Tài nguyên và Môi trường chủ trì thực hiện.
  • Chuẩn hóa hệ thống cơ sở dữ liệu địa lý môi trường GIS ở tỷ lệ 1:10.000 đến 1:50.000 cho hơn 500 cơ sở sản xuất có lưu lượng khí thải lớn, hoàn thành trong vòng 18 tháng dưới sự phối hợp giữa Viện Cơ học Ứng dụng và Chi cục Bảo vệ Môi trường.
  • Ban hành cơ chế cấp phép xả thải động dựa trên dự báo khí tượng thời gian thực thay thế định ngạch cố định, nhằm giảm ít nhất 25% tổng lượng phát thải công nghiệp vào những ngày dự báo có nghịch đảo nhiệt cao, áp dụng thí điểm trong 6 tháng tại các khu chế xuất lớn.
  • Nâng cấp cấu trúc thuật toán của phần mềm CAP lên phiên bản mới nhằm bổ sung module tính toán phản ứng hóa học quang hóa và sol khí mịn, rút ngắn thời gian xử lý mô phỏng 3D xuống dưới 3 phút cho mỗi kịch bản dự báo, triển khai trong giai đoạn 24 tháng bởi nhóm nghiên cứu tin học môi trường.

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

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường: Các chuyên viên tại Bộ và Sở Tài nguyên và Môi trường có thể ứng dụng luận văn làm cơ sở khoa học để thẩm định báo cáo đánh giá tác động môi trường, quy hoạch không gian các cụm công nghiệp và thiết lập kịch bản ứng phó sự cố rò rỉ khí độc.
  • Kỹ sư môi trường và chuyên viên tư vấn công nghệ: Các kỹ sư sử dụng trực tiếp các công thức toán lý và mô hình vệt khói trong tài liệu để tính toán chiều cao tối ưu cho ống khói nhà máy, thiết kế hệ thống xử lý khí thải đạt chuẩn xả thải môi trường.
  • Nhà nghiên cứu và giảng viên chuyên ngành Tin học Môi trường: Các nhà khoa học thuộc lĩnh vực cơ học chất lỏng, khí tượng và mô hình hóa tính toán có thể kế thừa các thuật toán sai phân phân rã, phương pháp hạt lớn và bài toán liên hợp để phát triển các phần mềm mô phỏng chất lượng không khí thế hệ mới.
  • Ban quản lý khu chế xuất và khu công nghiệp: Các nhà quản lý hạ tầng công nghiệp áp dụng mô hình để giám sát tổng lượng phát thải nội khu, điều phối lịch vận hành của các nhà máy thành viên theo điều kiện khí tượng nhằm hạn chế ô nhiễm không khí cục bộ.

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

Phần mềm CAP giải quyết bài toán lan truyền chất ô nhiễm dựa trên khung mô hình toán học nào? Phần mềm kết hợp mô hình vệt khói Gauss cho các nguồn điểm đơn lẻ và hệ phương trình vi phân bảo toàn khối lượng 3 chiều cho quy mô toàn vùng. Nhờ áp dụng phương pháp phân rã toán tử theo sơ đồ Yanenko trên lưới 20 x 20 x 10 điểm, hệ thống mô phỏng chuẩn xác quá trình bình lưu và khuếch tán rối, loại bỏ hiệu quả sai số khuếch tán nhân tạo.

Mô hình xử lý các yếu tố khí tượng phức tạp và địa hình bề mặt bằng cách nào? Mô hình vận dụng lý thuyết đồng dạng Monin-Obukhov để xác định profin gió thẳng đứng và công thức thực nghiệm Briggs tính toán độ nâng vệt khói hữu dụng. Tác động của địa hình và đặc tính lớp phủ mặt đất được số hóa thông qua tham số độ nhám bề mặt và góc nghiêng tọa độ, giúp mô tả chính xác sự tích tụ chất ô nhiễm sát mặt đất.

Ưu điểm nổi bật của phương pháp số trị trong luận văn so với mô hình hộp truyền thống là gì? Khác với mô hình hộp giả định pha trộn đều khiến nồng độ đỉnh bị giảm ước tính tới 50%, phương pháp sai phân phân rã kết hợp hạt lớn cho phép theo dõi chi tiết gradien nồng độ theo 3 chiều không gian. Thuật toán giúp tiết kiệm hơn 40% bộ nhớ máy tính và cho phép thay đổi bước thời gian linh hoạt từ 1 đến 30 phút.

Phần mềm CAP có thể ứng dụng cho những loại nguồn thải và chất ô nhiễm cụ thể nào? Hệ thống xử lý hiệu quả các chất ô nhiễm công nghiệp phổ biến như lưu huỳnh đioxit, oxit cacbon, nitơ đioxit và bụi lắng trọng trường. Nguồn phát thải được phân loại linh hoạt gồm nguồn điểm (ống khói đơn lẻ), nguồn đường (giao thông vận tải) và nguồn diện (khu công nghiệp), cho phép tính toán đồng thời hơn 50 nguồn thải xen kẽ.

Dữ liệu đầu vào bắt buộc để vận hành mô hình dự báo lan truyền không khí gồm những gì? Người dùng cần cung cấp ba nhóm thông tin cơ bản: thông số nguồn thải (tọa độ vị trí, lưu lượng khí thoát, nhiệt độ phát thải, đường kính ống), thông số khí tượng quan trắc (vận tốc gió nền, hướng gió, độ ổn định khí quyển) và bản đồ số hóa địa hình khu vực. Cơ sở dữ liệu của CAP sẽ tự động đồng bộ để tính toán và xuất bản đồ đường đồng mức.

Kết luận

  • Xây dựng thành công hệ thống mô hình toán lý 3 chiều toàn diện mô tả quá trình bình lưu và khuếch tán chất ô nhiễm trong lớp biên khí quyển đô thị.
  • Hoàn thiện phần mềm CAP tích hợp trực tiếp với cơ sở dữ liệu chuyên dụng và công nghệ thông tin địa lý GIS để xuất kết quả trực quan dạng bản đồ đường đồng mức.
  • Giải quyết triệt để bài toán khử độ nhớt giả trong các sơ đồ sai phân hữu hạn, đảm bảo độ chính xác tính toán đạt trên 90% so với thực nghiệm.
  • Ứng dụng mô phỏng thành công cho các kịch bản phát thải công nghiệp thực tế tại Thành phố Hồ Chí Minh và khu vực kinh tế trọng điểm phía Nam.
  • Khẳng định tính tiên phong trong việc kết hợp công nghệ thông tin và khoa học môi trường, đặt nền móng cho các hệ thống hỗ trợ ra quyết định tại Việt Nam.

Đóng góp lớn nhất của công trình là cung cấp một công cụ phần mềm tự động hóa phục vụ công tác đánh giá hiện trạng và dự báo chất lượng không khí mang tầm chiến lược. Trong lộ trình 12 tháng tới, nhóm nghiên cứu hướng tới mở rộng mô hình phản ứng quang hóa và hoàn thiện giải pháp dự báo thời gian thực trong 24 tháng tiếp theo. Các cơ quan quản lý đô thị và doanh nghiệp cần đẩy mạnh việc ứng dụng các giải pháp số hóa như phần mềm CAP vào thực tiễn quy hoạch để kiến tạo môi trường sống trong lành và phát triển bền vững.