Tổng quan nghiên cứu

Rừng giữ vai trò nền tảng trong việc bảo tồn đa dạng sinh học, điều tiết nguồn nước, phòng hộ đầu nguồn và giảm phát thải khí nhà kính nhằm ứng phó với biến đổi khí hậu toàn cầu. Tại tỉnh Bình Dương, tiến trình công nghiệp hóa và đô thị hóa diễn ra nhanh chóng với tốc độ tăng trưởng kinh tế huyện Phú Giáo đạt bình quân 12,44% mỗi năm, tạo ra áp lực rất lớn lên công tác quản lý tài nguyên thiên nhiên. Theo số liệu thống kê ngành lâm nghiệp, tổng diện tích đất có rừng của huyện Phú Giáo đạt 4.651,60 ha, trong đó rừng tự nhiên chỉ chiếm 376,01 ha (tương đương 8,08%) và rừng trồng chiếm 4.275,59 ha (chiếm 91,92%), bên cạnh 881,79 ha đất có rừng trồng chưa thành rừng thuộc quy hoạch.

Phương pháp kiểm kê rừng truyền thống ngoài thực địa đòi hỏi nguồn kinh phí lớn, tốn nhiều nhân lực và thường kéo dài chu kỳ từ 3 đến 5 năm, dẫn đến sự chậm trễ trong việc cập nhật biến động tài nguyên. Trong khi đó, việc ứng dụng ảnh viễn thám quang học đơn lẻ thường bị gián đoạn do mây mù nhiệt đới che phủ, còn ảnh radar đơn lẻ lại gặp trở ngại bởi hiện tượng nhiễu hạt speckle phức tạp. Luận văn tập trung giải quyết bài toán định lượng trữ lượng rừng bằng cách tích hợp dữ liệu viễn thám mở Sentinel-1 (Radar khẩu độ tổng hợp băng C) và Sentinel-2 (Quang học đa phổ) tại huyện Phú Giáo, tỉnh Bình Dương trong giai đoạn 2015-2017.

Nghiên cứu mang ý nghĩa thực tiễn to lớn khi thiết lập quy trình giám sát rừng với chi phí giảm hơn 60% so với phương pháp ngoại nghiệp truyền thống, rút ngắn chu kỳ cập nhật dữ liệu xuống còn 5 đến 12 ngày một lần, cung cấp cơ sở dữ liệu số chuẩn xác phục vụ quy hoạch bảo vệ rừng và chi trả dịch vụ môi trường rừng.

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 2 lý thuyết viễn thám không gian chủ đạo:

  • Lý thuyết phản xạ phổ quang học: Dựa trên đặc tính quang phổ của thảm thực vật, diệp lục trong lá cây hấp thụ mạnh bức xạ mặt trời ở dải sóng đỏ (Red, bước sóng 665 nm) và tán xạ mạnh mẽ ở dải sóng cận hồng ngoại (NIR, bước sóng 842 nm). Tỷ số giữa 2 dải sóng này cho phép xác định Chỉ số khác biệt thực vật (NDVI) phản ánh mật độ sinh khối xanh.
  • Lý thuyết tán xạ ngược vi ba Radar khẩu độ tổng hợp (SAR): Dựa trên sự tương tác của sóng vi ba băng C (bước sóng 5,6 cm) với tán cây và thân cây. Các chế độ phân cực kép VV (phát dọc - thu dọc) và VH (phát dọc - thu ngang) cung cấp thông tin chuyên sâu về độ nhám bề mặt, độ ẩm và cấu trúc hình thái không gian của lâm phần.

Khung phân tích vận dụng 4 khái niệm then chốt: Trữ lượng lâm phần (MGO, tính bằng m³/ha), Hệ số tán xạ ngược Sigma Nought (đo bằng thang decibel - dB), Chỉ số thực vật chuẩn hóa (NDVI từ -1,0 đến +1,0) và Sai số trung phương (RMSE) dùng để kiểm định độ chính xác mô hình.

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

Nguồn dữ liệu sử dụng gồm ảnh vệ tinh Sentinel-1A chụp ngày 22/04/2016 (chế độ IW GRDH) và ảnh Sentinel-2A chụp ngày 04/04/2016 (mức xử lý L1C) với trường phủ rộng 290 km, độ phân giải không gian 10m x 10m. Dữ liệu thực địa gồm 102 ô tiêu chuẩn hình chữ nhật diện tích 1.000 m² (kích thước 33,3m x 30m), định vị bằng máy GPS chuyên dụng có độ chính xác từ 3m đến 5m.

Phương pháp chọn mẫu phân tầng ngẫu nhiên có hệ thống được triển khai trên 3 trạng thái rừng điển hình: rừng trồng keo, rừng trồng cao su và rừng thường xanh phục hồi. Toàn bộ 102 ô mẫu được chia thành 2 tập dữ liệu độc lập: 73 ô mẫu (chiếm 71,57%) dùng để thiết lập mô hình hồi quy và 29 ô mẫu (chiếm 28,43%) dùng để kiểm định độ chính xác mô hình.

Quy trình phân tích dữ liệu được thực hiện trên phần mềm SNAP của Cơ quan Không gian Châu Âu (ESA), Excel và ArcGIS. Ảnh Sentinel-1 được định chuẩn bức xạ về giá trị Sigma Nought và lọc nhiễu bằng thuật toán Sigma Lee với các kích thước cửa sổ 5x5 và 7x7. Ảnh Sentinel-2 được tính toán chỉ số NDVI. Lý do lựa chọn phân tích hồi quy tích hợp đa biến quang học và radar là nhằm khai thác tối đa ưu thế tương hỗ: viễn thám quang học nhạy cảm với diệp lục tầng tán trên, trong khi sóng radar băng C đâm xuyên ghi nhận cấu trúc thân cành, qua đó triệt tiêu hiện tượng bão hòa tín hiệu quang học ở các lâm phần có mật độ tán dày.

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

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

  • Đặc trưng phản xạ radar: Sau khi định chuẩn và lọc nhiễu bằng thuật toán Sigma Lee trên phần mềm SNAP, giá trị tán xạ ngược Sentinel-1 đạt ngưỡng từ -25,5119 dB đến 28,4213 dB ở phân cực VV, và từ -28,5090 dB đến 10,9877 dB ở phân cực VH. Phân cực chéo VH thể hiện độ nhạy cao hơn 18,4% so với phân cực VV trong việc phân biệt cấu trúc tán và mật độ cành nhánh cây rừng.
  • Phân bố chỉ số thực vật: Kênh ảnh NDVI chiết xuất từ Sentinel-2 dao động từ -0,3046 đến 1,0000. Trong đó, thảm cỏ và đất sản xuất nông nghiệp đạt giá trị từ 0,30 đến 0,60; các trạng thái rừng trồng cao su, keo và rừng tự nhiên thường xanh phục hồi có chỉ số NDVI đạt trên 0,60, phản ánh độ che phủ tầng tán rất tốt.
  • Hiệu quả mô hình tích hợp: Mô hình đơn biến sử dụng riêng NDVI chỉ đạt hệ số xác định R² khoảng 0,512, mô hình dùng riêng phân cực VH đạt R² khoảng 0,465. Khi kết hợp đa biến giữa NDVI và tỷ số phân cực radar (Sig_VV/Sig_VH), hệ số xác định R² tăng lên mức 0,748, tương đương mức cải thiện khả năng giải thích phương sai tăng hơn 32,5% so với mô hình quang học đơn lẻ.
  • Kiểm định độ tin cậy: Đánh giá độc lập trên 29 ô tiêu chuẩn cho thấy sai số trung phương RMSE của mô hình tích hợp giảm xuống còn 17,85 m³/ha, thấp hơn rõ rệt so với mức sai số 28,40 m³/ha khi chỉ dùng ảnh quang học đơn lẻ (sai số giảm 37,15%).

Thảo luận kết quả

Sự vượt trội của mô hình tích hợp bắt nguồn từ cơ chế tương tác vật lý bổ trợ giữa hai loại cảm biến. Bước sóng quang học Sentinel-2 phản ánh màu sắc và hoạt tính quang hợp của lớp lá mặt ngoài nhưng nhanh chóng chạm ngưỡng bão hòa khi trữ lượng lâm phần vượt quá 120 m³/ha. Ngược lại, sóng vi ba Sentinel-1 có khả năng xuyên thấu tán lá, tương tác với các cành nhánh và thân cây, giúp lượng hóa cấu trúc thể tích thực của rừng.

Kết quả này hoàn toàn nhất quán với các nghiên cứu quốc tế tiêu biểu như công trình của Trung tâm Viễn thám Malaysia khi kết hợp ảnh radar ERS-1 và quang học Landsat TM làm tăng độ chính xác từ 69,3% lên 74,6%. Tại Việt Nam, việc làm chủ nguồn dữ liệu mở Sentinel cho thấy tiềm năng thay thế các nguồn tư liệu ảnh thương mại đắt tiền như SPOT-5 hay ALOS PALSAR.

Về mặt trực quan hóa, toàn bộ dữ liệu ước tính được thể hiện dưới dạng bản đồ raster phân giải 10m x 10m (100 m² mỗi điểm ảnh) trên nền tảng ArcGIS. Phân bố trữ lượng được hiển thị rõ nét thông qua biểu đồ phân tán hồi quy (Scatter Plot) kết hợp bảng phân cấp 5 mức trữ lượng: rừng nghèo dưới 50 m³/ha, rừng trung bình từ 50 đến 100 m³/ha, rừng giàu từ 100 đến 150 m³/ha và rừng rất giàu trên 150 m³/ha, hỗ trợ trực quan cho công tác quy hoạch lâm nghiệp tại địa phương.

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

Nhằm chuyển giao kết quả nghiên cứu vào thực tiễn quản lý rừng tại huyện Phú Giáo và tỉnh Bình Dương, 4 giải pháp trọng tâm được kiến nghị thực hiện:

  • Tự động hóa quy trình xử lý ảnh viễn thám đa nguồn: Chi cục Kiểm lâm tỉnh Bình Dương chủ trì phối hợp với Viện Sinh thái rừng và Môi trường xây dựng quy trình tiền xử lý ảnh bán tự động trên nền tảng mã nguồn mở. Mục tiêu giảm 70% thời gian xử lý dữ liệu từ 15 ngày xuống dưới 48 giờ mỗi chu kỳ ảnh, duy trì sai số RMSE dưới 16,5 m³/ha. Thời gian thực hiện từ quý I/2024 đến quý IV/2024.
  • Xây dựng cơ sở dữ liệu số và bản đồ biến động trữ lượng định kỳ: Ủy ban nhân dân huyện Phú Giáo phối hợp với Phòng Nông nghiệp và PTNT tổ chức cập nhật định kỳ 6 tháng một lần bản đồ trữ lượng cho toàn bộ 4.651,60 ha đất có rừng. Mục tiêu giám sát 100% diện tích rừng tự nhiên (376,01 ha) và rừng trồng (4.275,59 ha), phát hiện kịp thời các biến động suy giảm rừng với độ chính xác không gian trên 92%. Thời gian thực hiện trong giai đoạn 2024 - 2026.
  • Chuẩn hóa mạng lưới ô tiêu chuẩn định vị vĩnh viễn ngoài thực địa: Ban Quản lý rừng và các doanh nghiệp cao su, lâm nghiệp địa phương nâng cấp từ 102 ô tiêu chuẩn hiện có lên hệ thống 160 ô tiêu chuẩn định vị GPS cố định. Mục tiêu cung cấp dữ liệu kiểm chứng mặt đất liên tục, nâng độ chính xác hiệu chuẩn mô hình đạt trên 88% trong giai đoạn 2025 - 2027.
  • Đào tạo chuyển giao kỹ thuật công nghệ viễn thám tiên tiến: Trường Đại học Lâm nghiệp phối hợp với Sở Tài nguyên và Môi trường tỉnh Bình Dương tổ chức các khóa tập huấn chuyên sâu về phần mềm SNAP và ArcGIS. Mục tiêu bồi dưỡng cho ít nhất 45 lượt cán bộ kỹ thuật kiểm lâm cấp huyện, đảm bảo 100% đơn vị cơ sở có khả năng tự chủ vận hành mô hình trước năm 2028.

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

Công trình nghiên cứu mang lại giá trị học thuật và giải pháp ứng dụng thiết thực cho 4 nhóm đối tượng:

  • Cơ quan quản lý nhà nước về lâm nghiệp và tài nguyên môi trường (Cục Lâm nghiệp, Chi cục Kiểm lâm các tỉnh): Ứng dụng khung kỹ thuật tích hợp Sentinel-1 và Sentinel-2 để thực hiện điều tra kiểm kê rừng định kỳ, theo dõi biến động diện tích và chất lượng rừng với chi phí ngân sách tiết kiệm hơn 60% so với phương pháp đo đếm thủ công.
  • Các viện nghiên cứu và trường đại học chuyên ngành (Lâm nghiệp, Khoa học Môi trường, Địa tin học): Sử dụng quy trình định chuẩn Sigma Nought, lọc nhiễu Sigma Lee và tích hợp chỉ số NDVI làm tài liệu tham khảo cho các nghiên cứu mô hình hóa sinh khối rừng nhiệt đới bằng viễn thám đa nguồn.
  • Doanh nghiệp lâm nghiệp và các nông trường cao su quy mô lớn (điển hình như các đơn vị quản lý hơn 36.249 ha cao su tại Phú Giáo): Khai thác mô hình toán học để định lượng định kỳ trữ lượng gỗ, dự báo sinh khối khai thác và tối ưu hóa kế hoạch luân canh rừng trồng.
  • Các tổ chức phát triển dự án các-bon rừng và dịch vụ môi trường rừng (REDD+, thị trường tín chỉ các-bon): Sử dụng bản đồ raster trữ lượng 10m x 10m để thiết lập đường phát thải tham chiếu (FREL/FRL), phục vụ công tác đo đạc, báo cáo và thẩm định (MRV) độc lập với độ tin cậy vượt trên 85%.

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

Tại sao cần kết hợp ảnh Sentinel-1 và Sentinel-2 thay vì chỉ sử dụng một loại ảnh đơn lẻ? Ảnh quang học Sentinel-2 phản ánh mật độ diệp lục tầng tán qua chỉ số NDVI nhưng dễ bị bão hòa ở rừng rậm, trong khi Sentinel-1 (Radar SAR băng C) có khả năng xuyên tán ghi nhận cấu trúc thân cành qua phân cực VH và VV. Sự tích hợp bổ trợ giúp tăng hệ số xác định R² từ 0,512 lên 0,748 và giảm 37,15% sai số ước tính trữ lượng.

Sai số trung phương RMSE của mô hình đạt mức bao nhiêu và có đáp ứng tiêu chuẩn kiểm kê rừng không? Kết quả thử nghiệm trên 29 ô tiêu chuẩn kiểm chứng độc lập cho sai số trung phương RMSE đạt 17,85 m³/ha. Mức sai số này hoàn toàn đáp ứng tiêu chuẩn ngành lâm nghiệp Việt Nam đối với công tác điều tra kiểm kê tài nguyên rừng cấp huyện ở tỷ lệ bản đồ 1:25.000.

Dữ liệu viễn thám Sentinel mang lại lợi thế kinh tế gì so với các nguồn ảnh vệ tinh thương mại? Dữ liệu Sentinel-1 và Sentinel-2 được Cơ quan Không gian Châu Âu cung cấp hoàn toàn miễn phí với độ phân giải không gian 10m và chu kỳ chụp lặp lại 5 ngày. Điều này giúp địa phương tiết kiệm hàng tỷ đồng kinh phí mua ảnh thương mại như SPOT-5 hay ALOS PALSAR mà vẫn bảo đảm độ chính xác giám sát.

Quy mô 102 ô tiêu chuẩn thực địa có đủ độ tin cậy đại diện cho toàn huyện Phú Giáo không? Mẫu điều tra gồm 102 ô tiêu chuẩn diện tích 1.000 m² được bố trí phân tầng ngẫu nhiên trên toàn bộ 4.651,60 ha đất có rừng của huyện, bao gồm rừng trồng keo, cao su và rừng thường xanh phục hồi. Cỡ mẫu này đảm bảo tính đại diện thống kê với độ tin cậy đạt trên 90% theo quy chuẩn điều tra rừng.

Phương pháp trong luận văn có thể áp dụng cho các khu vực đồi núi có độ dốc lớn không? Phương pháp hoàn toàn có thể mở rộng cho các vùng địa hình phức tạp, nhưng cần tích hợp thêm bước hiệu chỉnh bức xạ địa hình (Terrain Flattening) bằng mô hình số độ cao DEM độ phân giải cao nhằm loại bỏ ảnh hưởng bóng râm và biến dạng góc chiếu của sóng radar trên các sườn dốc trên 15 độ.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp xác định trữ lượng rừng thông qua việc tích hợp nguồn dữ liệu viễn thám mở Sentinel-1 và Sentinel-2 tại huyện Phú Giáo, tỉnh Bình Dương.
  • Thiết lập mô hình hồi quy tối ưu kết hợp giữa chỉ số NDVI và tỷ số tán xạ ngược radar (Sig_VV/Sig_VH) đạt hệ số R² = 0,748, giảm sai số RMSE xuống còn 17,85 m³/ha trên 29 ô mẫu kiểm chứng độc lập.
  • Thành lập bản đồ phân bố trữ lượng rừng số hóa tỷ lệ 1:25.000 với độ phân giải điểm ảnh 10m x 10m trên nền GIS, phục vụ quản lý 4.651,60 ha đất có rừng toàn huyện.
  • Hoàn thiện quy trình công nghệ chuẩn từ tiền xử lý SNAP, lọc nhiễu Sigma Lee đến raster hóa không gian, giúp tiết kiệm trên 60% kinh phí so với phương pháp điều tra thực địa truyền thống.
  • Kế hoạch giai đoạn 2024 - 2026 sẽ mở rộng thuật toán trên nền tảng Google Earth Engine để giám sát biến động sinh khối cho toàn bộ 36.249 ha diện tích cao su trên địa bàn tỉnh Bình Dương. Các nhà quản lý và nghiên cứu lâm nghiệp nên áp dụng ngay quy trình viễn thám Sentinel tích hợp để hiện đại hóa công tác quản trị rừng và sẵn sàng tham gia thị trường tín chỉ các-bon quốc tế.