Tổng quan nghiên cứu

Thành phố Hồ Chí Minh là trung tâm kinh tế, khoa học và công nghệ lớn nhất cả nước với diện tích tự nhiên khoảng 2.056,5 km² (chiếm 0,76% diện tích toàn quốc) cùng quy mô dân số thực tế đạt gần 14 triệu người. Tốc độ đô thị hóa nhanh chóng và mật độ xây dựng dày đặc đã làm suy giảm nghiêm trọng không gian xanh, thay thế bề mặt tự nhiên bằng các lớp phủ nhân tạo không thấm. Quá trình này kích hoạt hiệu ứng đảo nhiệt đô thị (Urban Heat Island - UHI), khiến nền nhiệt nội đô tăng cao hơn các vùng ngoại ô lân cận từ 3°C đến 10°C, thường xuyên xuất hiện những đợt nắng nóng gay gắt với mức nhiệt ngoài trời vượt ngưỡng 37°C đến 40°C trong mùa khô.

Tình trạng căng thẳng nhiệt đô thị không chỉ làm gia tăng mức tiêu thụ năng lượng làm mát mà còn ảnh hưởng trực tiếp đến sức khỏe cộng đồng, gia tăng nguy cơ sốc nhiệt và suy giảm năng suất lao động. Để giải quyết bài toán giám sát nhiệt độ trên diện rộng mà phương pháp quan trắc trạm khí tượng truyền thống bị giới hạn bởi tính cục bộ, luận văn thạc sĩ chuyên ngành Quản lý Tài nguyên và Môi trường tập trung vào mục tiêu: Ứng dụng kỹ thuật xử lý ảnh viễn thám vệ tinh kết hợp mô hình thống kê để xác định các tham số vi khí hậu gồm nhiệt độ không khí và độ ẩm tương đối, từ đó mô phỏng phân bố không gian chỉ số bất tiện nghi nhiệt (Discomfort Index - DI).

Phạm vi không gian nghiên cứu bao gồm 19 quận nội thành và 4 huyện ngoại thành phía Bắc Thành phố Hồ Chí Minh (không tính huyện Cần Giờ) dựa trên nguồn dữ liệu ảnh vệ tinh quang học và nhiệt Landsat 8 thu nhận vào ngày 19/01/2019. Kết quả đề tài mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng phục vụ công tác quy hoạch đô thị bền vững, thích ứng với biến đổi khí hậu.

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 truyền xạ bức xạ nhiệt và định luật vật lý Planck, mô tả mối quan hệ giữa nhiệt độ động lực bề mặt đất (Land Surface Temperature - Ts) với năng lượng bức xạ phát ra từ vật thể ở các kênh phổ hồng ngoại nhiệt (Thermal Infrared - TIR). Để hiệu chỉnh độ phát xạ bề mặt, nghiên cứu vận dụng định luật Kirchhoff và mô hình tương quan độ phát xạ thực nghiệm thông qua tỷ lệ che phủ của thực vật do Valor & Caselles phát triển.

Bên cạnh đó, luận văn áp dụng tiêu chuẩn tiện nghi nhiệt ASHRAE 55 của Hiệp hội Kỹ sư Nhiệt lạnh Hoa Kỳ và mô hình chỉ số bất tiện nghi nhiệt DI do Thom đề xuất năm 1959. Chỉ số DI được xác định thông qua hàm tuyến tính tích hợp giữa nhiệt độ không khí (Air Temperature - Ta) và độ ẩm tương đối (Relative Humidity - RH). Các ngưỡng tiện nghi nhiệt được đối chiếu khoa học với Quy chuẩn kỹ thuật quốc gia QCVN 26:2016/BYT về vi khí hậu tại nơi làm việc và Luật An toàn, vệ sinh lao động số 84/2015/QH13 nhằm lượng hóa mức độ căng thẳng nhiệt của cơ thể người.

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

Nguồn dữ liệu chính của đề tài là cảnh ảnh vệ tinh Landsat 8 (tích hợp bộ cảm biến OLI và TIRS) thu nhận ngày 19/01/2019 với độ phân giải không gian kênh quang học 30m và kênh nhiệt 100m, được cung cấp bởi Cục Khảo sát Địa chất Hoa Kỳ (USGS). Đồng thời, đề tài sử dụng chuỗi số liệu quan trắc thực tế về nhiệt độ và độ ẩm thu thập tại các trạm khí tượng thủy văn trên địa bàn Thành phố Hồ Chí Minh. Quy trình phân tích được thực hiện từ tháng 02/2019 đến tháng 06/2019 theo các bước phương pháp luận chặt chẽ:

  1. Xử lý ảnh số và hiệu chỉnh bức xạ: Chuyển đổi giá trị số nguyên DN của kênh nhiệt sang độ chói phổ thông qua các hệ số khuếch đại (Gain) và độ lệch (Bias), chuyển đổi về nhiệt độ chói Tb, sau đó tích hợp chỉ số thực vật khác biệt chuẩn hóa (NDVI) để tính toán độ phát xạ bề mặt và suy ra nhiệt độ bề mặt đất Ts.
  2. Xây dựng mô hình hồi quy thống kê: Phân tích tương quan tuyến tính và phi tuyến giữa Ts với nhiệt độ không khí Ta, cũng như giữa NDVI với độ ẩm tương đối RH tại các vị trí trạm quan trắc mặt đất. Lý do lựa chọn mô hình hồi quy là nhằm mở rộng dữ liệu điểm quan trắc thành bề mặt trường thông tin không gian liên tục trên toàn đô thị.
  3. Phân loại lớp phủ và kiểm định độ chính xác: Phân chia 4 kiểu bề mặt (mặt nước, đô thị, đất trống, thực vật) theo ngưỡng NDVI. Cỡ mẫu kiểm định được xác định ngẫu nhiên phân tầng với 200 điểm mẫu phân bố đều trên 23 quận huyện, đối chiếu trực tiếp qua ảnh độ phân giải cao Google Earth và khảo sát thực địa để tính toán ma trận nhầm lẫn (confusion matrix).

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

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

Quá trình tiền xử lý và nắn chỉnh hình học trên nền bản đồ địa hình tỷ lệ 1/25.000 (hệ tọa độ UTM, WGS 84, múi chiếu 48N) đạt độ chính xác cao với sai số toàn phương trung bình RMSE nhỏ hơn 0,5 pixel. Kết quả phân loại lớp phủ bề mặt đạt độ chính xác toàn cục 97,0334% và hệ số Kappa vượt 0,94. Các phát hiện cốt lõi từ dữ liệu phân tích bao gồm:

  1. Phân bố lớp phủ bề mặt: Ngưỡng chỉ số NDVI phân định rõ nét 4 kiểu lớp phủ: mặt nước (NDVI dưới -0,26), đô thị và bề mặt không thấm (NDVI từ -0,26 đến 0,1), đất trống (NDVI từ 0,1 đến 0,2) và thực vật (NDVI trên 0,2). Khu vực 13 quận nội thành cũ có tỷ lệ bề mặt không thấm bê tông hóa chiếm trên 70% diện tích.
  2. Căng thẳng nhiệt tập trung cao độ tại vùng lõi đô thị: Kết quả mô phỏng chỉ số DI cho thấy khu vực đô thị có tỷ lệ diện tích thuộc vùng không tiện nghi nhiệt lên tới 76%, trong đó vùng rất không tiện nghi chiếm tới 87%. Nhiệt độ không khí trung bình tại các quận trung tâm như Quận 1, Quận 3, Tân Bình, Bình Thạnh dao động từ 31°C đến 34°C vào thời điểm khảo sát.
  3. Hiệu ứng làm mát vượt trội từ thảm thực vật: Tại các khu vực ngoại thành có mật độ cây xanh và mặt nước lớn (như Củ Chi, Hóc Môn), tỷ lệ diện tích đạt mức tiện nghi nhiệt lên đến 91%. Vùng thực vật chỉ ghi nhận 11,5% diện tích không tiện nghi và mức rất không tiện nghi chỉ chiếm khoảng 15%.
  4. Tương quan vi khí hậu giữa các phân vùng: Độ ẩm tương đối không khí tại vùng nội thành giảm thấp đáng kể (dao động từ 55% đến 65%) so với các vùng nông thôn ngoại thành (thường trên 75%), tạo ra sự tương phản nhiệt ẩm gay gắt giữa vùng trung tâm và vùng ven.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự phân hóa nhiệt độ môi trường là do cấu trúc bề mặt đô thị chứa nhiều vật liệu nhân tạo như bê tông, nhựa đường và các khối nhà cao tầng. Những vật liệu này có dung lượng nhiệt lớn, hấp thụ bức xạ sóng ngắn mặt trời vào ban ngày và phát xạ mạnh mẽ trở lại khí quyển lớp biên. Đồng thời, sự thiếu vắng thảm thực vật làm suy giảm quá trình thoát hơi nước làm mát tự nhiên.

Kết quả này hoàn toàn nhất quán với các công trình nghiên cứu về chỉ số DI tại thành phố Hàng Châu (Trung Quốc) của Wang Weiwu và nghiên cứu đảo nhiệt tại thành phố Akure (Nigeria) của Ahmed A. Balogun. Về mặt trình bày khoa học, các trường dữ liệu vi khí hậu được trực quan hóa sinh động qua hệ thống bản đồ nhiệt độ không khí Ta, bản đồ độ ẩm tương đối RH và bản đồ phân vùng tiện nghi nhiệt DI tỷ lệ 1/25.000. Dữ liệu diện tích và tỷ lệ % mức độ tiện nghi theo từng đơn vị hành chính được tổng hợp chi tiết qua các bảng ma trận thống kê và biểu đồ cột so sánh trực quan, minh chứng rõ rệt cho quy luật: mật độ xây dựng càng tăng thì chỉ số bất tiện nghi nhiệt càng mở rộng.

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

Dựa trên kết quả thực nghiệm, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm giảm thiểu tác động của căng thẳng nhiệt đô thị cho Thành phố Hồ Chí Minh:

  1. Phát triển và phục hồi mạng lưới mảng xanh đô thị: Ủy ban nhân dân Thành phố Hồ Chí Minh cùng Sở Xây dựng cần ưu tiên quỹ đất để nâng chỉ tiêu cây xanh đô thị đạt tối thiểu 3 đến 5 m²/người vào năm 2025 và hướng tới 7 m²/người vào năm 2030 tại 13 quận nội thành. Thiết lập các dải cây xanh cách ly có tán rộng quanh khu công nghiệp và phát triển công viên công cộng nhằm giảm nhiệt độ bức xạ bề mặt cục bộ từ 1,5°C đến 3°C.
  2. Quy chuẩn hóa vật liệu xây dựng và kiến trúc giảm nhiệt: Sở Quy hoạch - Kiến trúc cần ban hành quy định bắt buộc các dự án phát triển mới phải sử dụng vật liệu lát có độ phản xạ nhiệt cao (suất phản chiếu albedo trên 0,4), tăng cường giải pháp mái nhà xanh (green roofs) và vỉa hè tự thấm nước, đảm bảo tỷ lệ bề mặt thấm nước đạt tối thiểu 30% tổng diện tích khuôn viên dự án trong giai đoạn 2024 - 2028.
  3. Thiết lập hệ thống quan trắc vi khí hậu tự động bằng GIS và viễn thám: Sở Tài nguyên và Môi trường cần xây dựng cơ sở dữ liệu giám sát đảo nhiệt và lập bản đồ phân vùng chỉ số DI định kỳ 6 tháng một lần dựa trên ảnh vệ tinh đa thời gian Landsat và Sentinel, tích hợp vào cổng thông tin điều hành đô thị thông minh để cảnh báo sớm các điểm nóng nhiệt môi trường.
  4. Tăng cường kiểm soát an toàn vi khí hậu nơi làm việc: Liên đoàn Lao động Thành phố và ngành Y tế cần phối hợp triển khai giám sát nghiêm ngặt việc tuân thủ QCVN 26:2016/BYT và Luật An toàn, vệ sinh lao động số 84/2015/QH13. Doanh nghiệp cần bố trí ca làm việc hợp lý cho công nhân ngoài trời, trang bị hệ thống phun sương làm mát và thông gió nhà xưởng khi chỉ số DI vượt ngưỡng 28°C để bảo vệ sức khỏe người lao động trước năm 2025.

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

Nghiên cứu cung cấp nguồn tư liệu học thuật và giá trị ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  1. Cơ quan quản lý nhà nước và nhà hoạch định chính sách: Sở Tài nguyên và Môi trường, Sở Xây dựng, Sở Quy hoạch - Kiến trúc sử dụng dữ liệu bản đồ phân vùng nhiệt DI làm căn cứ pháp lý và khoa học để điều chỉnh quy hoạch phân khu, quy hoạch không gian xanh và chiến lược ứng phó biến đổi khí hậu cấp đô thị.
  2. Kiến trúc sư và kỹ sư quy hoạch đô thị: Ứng dụng mô hình tương quan giữa lớp phủ bề mặt và vi khí hậu để lựa chọn giải pháp vật liệu xây dựng sinh thái, thiết kế hướng đón gió tự nhiên và phân bổ hệ thống hồ điều hòa hợp lý cho các khu đô thị mới.
  3. Cộng đồng nghiên cứu khoa học, giảng viên và sinh viên: Khai thác quy trình công nghệ tích hợp viễn thám và GIS, công thức chuyển đổi nhiệt độ bề mặt Ts sang nhiệt độ không khí Ta và phương pháp tính toán chỉ số DI để mở rộng nghiên cứu sang các tỉnh thành lân cận như Bình Dương, Đồng Nai.
  4. Chuyên gia bảo hộ lao động và y tế công cộng: Ứng dụng các ngưỡng phân cấp tiện nghi nhiệt để đánh giá rủi ro nghề nghiệp, xây dựng cẩm nang phòng chống sốc nhiệt và cải thiện điều kiện vi khí hậu tại các nhà máy, công trường xây dựng.

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

Chỉ số bất tiện nghi nhiệt (DI) được tính toán từ ảnh vệ tinh như thế nào? Chỉ số DI được tính toán thông qua công thức kinh nghiệm của Thom (1959), kết hợp giữa nhiệt độ không khí Ta và độ ẩm tương đối RH. Cả hai thông số này được trích xuất gián tiếp từ kênh nhiệt Landsat 8 (thông qua nhiệt độ bề mặt Ts) và kênh quang học (thông qua chỉ số thực vật NDVI) nhờ các hàm hồi quy tương quan đa biến.

Tại sao nghiên cứu chỉ thực hiện ở phía Bắc TP.HCM mà không bao gồm huyện Cần Giờ? Huyện Cần Giờ có đặc thù sinh thái riêng biệt với rừng ngập mặn chiếm phần lớn diện tích tự nhiên và được bao bọc bởi mạng lưới sông ngòi ven biển. Việc loại trừ Cần Giờ nhằm tập trung phân tích chuyên sâu vào khu vực có tốc độ đô thị hóa nhanh và mật độ dân cư tập trung đông đúc tại 19 quận và 4 huyện phía Bắc.

Độ tin cậy của phương pháp phân loại lớp phủ bề mặt trong luận văn đạt bao nhiêu? Độ chính xác toàn cục của quá trình phân loại lớp phủ đạt mức 97,0334% và hệ số Kappa đạt trên 0,94. Kết quả này được kiểm định nghiêm ngặt thông qua ma trận nhầm lẫn với 200 điểm mẫu khảo sát thực địa kết hợp đối chiếu ảnh không gian độ phân giải cao trên Google Earth.

Hiệu ứng đảo nhiệt đô thị ảnh hưởng ra sao đến các quận nội thành TP.HCM? Nghiên cứu chỉ ra rằng khu vực nội thành có tỷ lệ diện tích không tiện nghi nhiệt lên đến 76% và mức rất không tiện nghi lên đến 87%. Hiện tượng đảo nhiệt khiến nhiệt độ bề mặt và nhiệt độ không khí tại các quận trung tâm cao hơn các vùng ngoại thành từ 3°C đến 5°C, gây oi bức kéo dài.

Các quy chuẩn pháp lý nào được sử dụng để đối chiếu đánh giá tiện nghi nhiệt? Nghiên cứu áp dụng và đối chiếu trực tiếp với Tiêu chuẩn tiện nghi nhiệt quốc tế ASHRAE 55, Quy chuẩn kỹ thuật quốc gia QCVN 26:2016/BYT về giá trị cho phép vi khí hậu tại nơi làm việc và các quy định của Luật An toàn, vệ sinh lao động số 84/2015/QH13 nhằm lượng hóa rủi ro sức khỏe cộng đồng.

Kết luận

  • Luận văn đã xây dựng thành công quy trình công nghệ ứng dụng ảnh vệ tinh Landsat 8 để trích xuất các tham số khí tượng Ta, RH và mô phỏng không gian chỉ số tiện nghi nhiệt DI với độ chính xác tin cậy.
  • Đánh giá chính xác hiện trạng lớp phủ bề mặt với độ chính xác toàn cục 97,0334%, chỉ ra tỷ lệ bê tông hóa cao tại 13 quận nội thành là nguyên nhân chính gây căng thẳng nhiệt.
  • Định lượng rõ nét vai trò của thảm thực vật khi vùng có độ phủ xanh cao đạt tỷ lệ tiện nghi nhiệt lên tới 91%, đối lập hoàn toàn với vùng đô thị có mức rất không tiện nghi chiếm tới 87%.
  • Đóng góp luận cứ khoa học thực chứng hỗ trợ tích cực cho công tác quy hoạch đô thị sinh thái, quản lý môi trường và bảo vệ sức khỏe người lao động theo tiêu chuẩn QCVN 26:2016/BYT.
  • Đề xuất lộ trình triển khai giám sát vi khí hậu định kỳ giai đoạn 2024 - 2030, mở ra hướng ứng dụng viễn thám đa thời gian phục vụ phát triển đô thị bền vững và thích ứng biến đổi khí hậu.