Tổng quan nghiên cứu

Sau khi mở rộng địa giới hành chính theo Nghị quyết số 15/2008/QH12 ngày 29/5/2008 của Quốc hội khóa XII, thành phố Hà Nội bước vào thời kỳ đô thị hóa và phát triển hạ tầng kỹ thuật với tốc độ nhanh chóng, nâng tổng diện tích tự nhiên toàn thành phố lên 335.859 ha vào năm 2018. Khu vực phía Tây thành phố, trọng tâm là 3 huyện Phúc Thọ, Thạch Thất và Quốc Oai, giữ vai trò cửa ngõ giao thương chiến lược với các trục hành lang kinh tế trọng điểm như Quốc lộ 21, Quốc lộ 32 và Đại lộ Thăng Long (đường Láng - Hòa Lạc).

Sự hình thành của các khu đô thị vệ tinh, khu công nghiệp và đặc biệt là Khu công nghệ cao Hòa Lạc đã tạo ra áp lực chuyển dịch mục đích sử dụng đất vô cùng mạnh mẽ. Quỹ đất nông nghiệp truyền thống nhanh chóng bị thu hẹp, nhường chỗ cho đất xây dựng cơ sở hạ tầng, khu dân cư và các dịch vụ phi nông nghiệp. Để giải quyết bài toán giám sát tài nguyên khách quan, luận văn thạc sĩ ngành Quản lý đất đai (Mã số: 8850103) của tác giả Khuất Thị Thu, dưới sự hướng dẫn khoa học của TS. Nguyễn Thị Thu Hiền tại Học viện Nông nghiệp Việt Nam, đã ứng dụng tích hợp công nghệ Viễn thám (Remote Sensing - RS) và Hệ thống thông tin địa lý (GIS).

Nghiên cứu được triển khai trong phạm vi không gian 3 huyện phía Tây Hà Nội với mốc thời gian 8 năm (2010 - 2018). Đóng góp nổi bật của luận văn là xây dựng quy trình chiết xuất thông tin biến động đạt độ chính xác phân loại tổng thể 96,11% (năm 2010) và 89,44% (năm 2018), cung cấp cơ sở dữ liệu số hóa tin cậy hỗ trợ các nhà hoạch định chính sách cân đối giữa mục tiêu tăng trưởng kinh tế và bảo vệ môi trường sinh thái bền vữ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 dựa trên nền tảng của hai lý thuyết trụ cột: Lý thuyết biến động lớp phủ và sử dụng đất (LUCC) của Turner và Lambin, kết hợp với Lý thuyết phân tích véc tơ thay đổi phổ (Change Vector Analysis - CVA) trong công nghệ viễn thám quang học. Mô hình nghiên cứu tích hợp quy trình xử lý ảnh vệ tinh đa thời gian và kỹ thuật phân loại thống kê có kiểm định theo thuật toán Xác suất cực đại (Maximum Likelihood Classification).

Bốn khái niệm khoa học cốt lõi được chuẩn hóa trong nghiên cứu bao gồm:

  1. Hiện trạng sử dụng đất: Biểu thị trạng thái khai thác bề mặt lãnh thổ thông qua hoạt động kinh tế, kỹ thuật của con người theo định nghĩa của FAO (1976) và Docuchaev.
  2. Biến động sử dụng đất: Sự thay đổi về diện tích, ranh giới và tính chất lớp phủ mặt đất qua thời gian dưới tác động tổng hợp của các nhân tố tự nhiên, kinh tế - xã hội và quy hoạch (Muller, 2003; Ellis, 2010).
  3. Công nghệ viễn thám: Phương pháp thu nhận thông tin bức xạ điện từ của các đối tượng mặt đất từ xa qua các kênh phổ phản xạ (dải sóng nhìn thấy 0,52 - 0,69 µm, cận hồng ngoại NIR 0,90 µm và hồng ngoại sóng ngắn SWIR 1,75 - 2,35 µm).
  4. Hệ thống thông tin địa lý (GIS): Hệ thống tích hợp phần cứng, phần mềm và dữ liệu không gian phục vụ thu thập, lưu trữ, chồng xếp bản đồ và mô hình hóa biến động không gian.

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

Nguồn dữ liệu của đề tài bao gồm ảnh vệ tinh Landsat 7 (bộ cảm ETM+) chụp năm 2010 và Landsat 8 (bộ cảm OLI/TIRS) chụp năm 2018 với độ phân giải không gian 30 m, kết hợp số liệu kiểm kê đất đai định kỳ của Bộ Tài nguyên và Môi trường, niên giám thống kê thành phố Hà Nội và bản đồ địa chính cơ sở.

Phương pháp chọn mẫu phân tầng ngẫu nhiên được áp dụng để thiết lập tập mẫu giải đoán gồm 6 lớp đối tượng: đất trồng lúa, đất trồng cây hàng năm khác, đất trồng rừng - cây lâu năm, đất xây dựng, đất mặt nước và đất chưa sử dụng. Cỡ mẫu khảo sát thực địa gồm 120 điểm định vị không gian bằng hệ thống GPS cầm tay phân bố đều trên toàn địa bàn 3 huyện, đảm bảo tính đại diện cho từng loại hình thảm phủ.

Lý do lựa chọn thuật toán Maximum Likelihood trong phần mềm ENVI là nhờ khả năng tính toán xác suất tiên nghiệm dựa trên phân bố chuẩn của các pixel mẫu, giúp giảm thiểu tối đa sai số hỗn hợp phổ giữa các đối tượng phức tạp. Ảnh vệ tinh sau khi nắn chỉnh hình học (sai số RMS < 0,5 pixel), hiệu chỉnh bức xạ khí quyển và cắt theo ranh giới hành chính được phân loại độc lập cho từng năm. Tiếp đó, dữ liệu vector được chuyển nạp vào phần mềm ArcGIS để tiến hành kỹ thuật chồng xếp không gian (Overlay), thành lập bản đồ biến động sử dụng đất tỷ lệ 1:50.000 trong giai đoạn nghiên cứu 2010 - 2018.

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

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

Quá trình phân tích không gian và thống kê diện tích đa thời gian đã chỉ ra 4 phát hiện quan trọng về thực trạng biến động đất đai tại 3 huyện Phúc Thọ, Thạch Thất và Quốc Oai:

  1. Độ chính xác phân loại ảnh số đạt độ tin cậy vượt trội: Kết quả đánh giá bằng ma trận sai số cho thấy ảnh Landsat năm 2010 đạt độ chính xác tổng thể 96,11% với hệ số Kappa tương ứng κ = 0,95. Ảnh Landsat 8 năm 2018 đạt độ chính xác tổng thể 89,44% với hệ số Kappa κ = 0,88, chứng minh tính chuẩn xác của quy trình giải đoán viễn thám.
  2. Diện tích đất trồng lúa sụt giảm sâu: Trong giai đoạn 2010 - 2018, đất trồng lúa ghi nhận mức suy giảm mạnh nhất với diện tích giảm 1.228,48 ha, tương đương tỷ lệ giảm khoảng 8,5% so với tổng quỹ đất lúa ban đầu tại khu vực.
  3. Đất xây dựng mở rộng nhanh chóng: Quá trình đô thị hóa và nâng cấp hạ tầng đã đẩy diện tích đất xây dựng năm 2018 tăng thêm 556,56 ha so với năm 2010, thể hiện tốc độ gia tăng trên 12% ở cả 3 huyện ngoại thành.
  4. Diện tích đất trồng rừng và cây lâu năm tăng trưởng mạnh mẽ: Nhóm đất rừng và cây lâu năm ghi nhận mức tăng ròng lên tới 1.423,08 ha sau 8 năm, phản ánh làn sóng tái cơ cấu nội ngành nông nghiệp rất rõ nét.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự thu hẹp 1.228,48 ha đất trồng lúa và sự gia tăng 556,56 ha đất xây dựng bắt nguồn từ việc thực hiện Quy hoạch chung xây dựng Thủ đô Hà Nội. Việc triển khai xây dựng hạ tầng các cụm công nghiệp, quy hoạch điểm dân cư nông thôn và mở rộng các tuyến giao thông liên vùng đã chuyển hóa trực tiếp các dải đất nông nghiệp màu mỡ thành mặt bằng phi nông nghiệp. Đồng thời, xu hướng chuyển đổi cơ cấu cây trồng từ độc canh lúa có năng suất thấp sang mô hình trang trại trồng cây ăn quả chất lượng cao và khoanh nuôi cây lâu năm đã thúc đẩy diện tích đất trồng rừng - cây lâu năm tăng 1.423,08 ha, mang lại hiệu quả kinh tế gấp 3 đến 5 lần cho người dân địa phương.

Kết quả này hoàn toàn tương đồng với xu thế biến động đất đai trên phạm vi cả nước giai đoạn 2010 - 2018, khi đất trồng cây lâu năm toàn quốc tăng 850,2 nghìn ha và đất phi nông nghiệp tăng 44,7 nghìn ha theo báo cáo của Bộ Tài nguyên và Môi trường. Dữ liệu nghiên cứu được biểu diễn trực quan qua ma trận sai số (Error Matrix), bảng thống kê so sánh diện tích chi tiết và sơ đồ biến động sử dụng đất tỷ lệ 1:50.000. Trên các biểu đồ phân bố và biểu đồ Sankey dòng luân chuyển đất đai, người xem có thể quan sát rõ ràng sự chuyển đổi không gian của hàng nghìn hecta giữa các nhóm đất nông nghiệp, phi nông nghiệp và đất chưa sử dụng qua từng giai đoạn.

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

Dựa trên các phân tích khoa học và thực tiễn biến động, luận văn đưa ra 4 nhóm giải pháp cụ thể:

  1. Thiết lập quy trình giám sát tài nguyên đất định kỳ bằng viễn thám: Sở Tài nguyên và Môi trường Hà Nội chủ trì phối hợp với các viện nghiên cứu ứng dụng ảnh vệ tinh Sentinel-2 (độ phân giải 10 m) và Landsat 8/9 để quét hiện trạng định kỳ 6 tháng một lần, nhằm phát hiện sớm các vi phạm lấn chiếm đất công với độ chính xác đạt trên 95%, hoàn thành khung cơ chế vận hành từ quý IV/2024.
  2. Số hóa và đồng bộ hóa cơ sở dữ liệu GIS địa chính: Ủy ban nhân dân các huyện Phúc Thọ, Thạch Thất, Quốc Oai chỉ đạo Phòng Tài nguyên và Môi trường hoàn thành việc chuẩn hóa, tích hợp 100% bản đồ biến động và hồ sơ địa chính lên nền tảng GIS dùng chung toàn thành phố, đảm bảo sai số vị trí ranh giới thửa đất dưới 0,5 m trong lộ trình 2024 - 2026.
  3. Quy hoạch khoanh vùng bảo vệ đất lúa và kiểm soát chuyển đổi cây trồng: Sở Nông nghiệp và Phát triển nông thôn phối hợp với chính quyền địa phương khoanh định ranh giới bảo vệ nghiêm ngặt tối thiểu 80% diện tích đất chuyên trồng lúa nước chất lượng cao; đồng thời hướng dẫn chuyển đổi có kiểm soát khoảng 1.500 ha đất canh tác kém hiệu quả sang cây ăn quả đặc sản tại các vùng gò đồi trong giai đoạn 2025 - 2030.
  4. Kiểm soát tốc độ bê tông hóa và tăng chỉ tiêu không gian xanh: Ban Quản lý Khu công nghệ cao Hòa Lạc và Ủy ban nhân dân các huyện siết chặt quy chuẩn cấp phép xây dựng dự án, bắt buộc duy trì mật độ che phủ cây xanh và mặt nước chuyên dùng tối thiểu 20% trên tổng diện tích các khu đô thị, cụm công nghiệp mới từ năm 2024 nhằm giảm áp lực suy giảm tầng nước ngầm và biến đổi vi khí hậu.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng cao cho 4 nhóm đối tượng trọng tâm:

  1. Cán bộ quản lý nhà nước và chuyên viên địa chính cấp huyện, xã: Cung cấp phương pháp luận và quy trình số hóa bản đồ hiện trạng giúp nâng cao chất lượng công tác lập quy hoạch, kế hoạch sử dụng đất hàng năm và giám sát biến động thực địa nhanh chóng.
  2. Các nhà khoa học, giảng viên và học viên cao học ngành Quản lý đất đai, Trắc địa - Viễn thám: Sử dụng tài liệu như một nghiên cứu điển hình về quy trình phân loại ảnh đa phổ (ENVI), thuật toán Maximum Likelihood và đánh giá độ chính xác phân loại bằng hệ số Kappa (κ = 0,88 - 0,95).
  3. Các đơn vị tư vấn quy hoạch đô thị và doanh nghiệp phát triển hạ tầng: Khai thác số liệu định lượng về tốc độ mở rộng 556,56 ha đất xây dựng và xu hướng dịch chuyển không gian để đánh giá tiềm năng thị trường, lập báo cáo tiền khả thi cho các dự án khu đô thị và cụm công nghiệp tại phía Tây Hà Nội.
  4. Sinh viên các ngành Khoa học Môi trường, Địa lý học và Hệ thống thông tin: Nắm bắt phương pháp xử lý dữ liệu không gian, kỹ năng chồng xếp lớp bản đồ trên ArcGIS và cách thức liên kết số liệu thống kê với các mô hình biến động lãnh thổ.

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

Tại sao nghiên cứu lại kết hợp ảnh Landsat 7 và Landsat 8 để đánh giá giai đoạn 2010 - 2018? Ảnh vệ tinh Landsat cung cấp nguồn tư liệu đa phổ hoàn toàn miễn phí, có tính kế thừa liên tục với độ phân giải không gian 30 m ở các kênh quang học. Sự kết hợp giữa bộ cảm ETM+ (năm 2010) và OLI/TIRS (năm 2018) đảm bảo tính đồng nhất về mặt hình học và quang phổ, giúp việc so sánh biến động đạt độ tin cậy khoa học cao nhất.

Độ tin cậy của thuật toán Maximum Likelihood trong phân loại ảnh được đảm bảo như thế nào? Thuật toán phân loại dựa trên hàm mật độ xác suất Bayes gán từng pixel vào lớp có xác suất cao nhất. Độ tin cậy được kiểm chứng qua ma trận sai số với 120 điểm mẫu thực địa, mang lại độ chính xác tổng thể 96,11% (năm 2010) và 89,44% (năm 2018), với chỉ số Kappa đều vượt ngưỡng 0,88, thể hiện sự trùng khớp rất cao giữa kết quả giải đoán và thực địa.

Diện tích đất trồng lúa giảm 1.228,48 ha đã chuyển dịch chủ yếu sang các loại hình đất nào? Phần lớn diện tích đất lúa giảm đã được chuyển hóa sang đất xây dựng (tăng 556,56 ha) nhằm phục vụ các dự án khu dân cư, hạ tầng giao thông kết nối và chuyển sang nhóm đất trồng cây lâu năm, cây ăn quả (tăng 1.423,08 ha) do người dân chủ động chuyển đổi mô hình canh tác nâng cao thu nhập.

Công nghệ Viễn thám và GIS mang lại lợi ích gì vượt trội so với phương pháp đo đạc truyền thống? Phương pháp viễn thám kết hợp GIS giúp cắt giảm hơn 60% thời gian và chi phí khảo sát trên diện tích lớn hàng chục nghìn hecta. Hệ thống cho phép cập nhật thông tin đa thời gian có tính chu kỳ, lưu trữ dữ liệu không gian nhất quán và hỗ trợ truy vấn, trích xuất dữ liệu tự động phục vụ ra quyết định tức thì.

Kết quả nghiên cứu của luận văn có ý nghĩa gì đối với công tác quy hoạch sử dụng đất hiện nay? Dữ liệu của luận văn cung cấp bức tranh định lượng chính xác về quy mô và tốc độ dịch chuyển quỹ đất giữa các ngành. Đây là căn cứ khoa học then chốt giúp các cấp chính quyền cân đối chỉ tiêu sử dụng đất, hạn chế việc phát triển đô thị tràn lan xâm lấn đất nông nghiệp màu mỡ và duy trì cân bằng sinh thái.

Kết luận

  • Luận văn đã định lượng chính xác bức tranh biến động sử dụng đất tại 3 huyện Phúc Thọ, Thạch Thất và Quốc Oai giai đoạn 2010 - 2018: đất lúa giảm 1.228,48 ha, đất cây lâu năm và rừng tăng 1.423,08 ha, đất xây dựng mở rộng 556,56 ha.
  • Khẳng định tính hiệu quả và độ tin cậy vượt trội của việc tích hợp công nghệ Viễn thám và GIS, đạt độ chính xác phân loại tổng thể 96,11% (năm 2010) và 89,44% (năm 2018) với hệ số Kappa tương ứng là 0,95 và 0,88.
  • Cung cấp hệ thống sản phẩm khoa học hoàn chỉnh gồm bản đồ hiện trạng, bản đồ biến động sử dụng đất tỷ lệ 1:50.000 và hệ thống bảng biểu chi tiết, phục vụ đắc lực cho công tác quản trị tài nguyên đất đai.
  • Lộ trình giai đoạn 2024 - 2030 cần tiếp tục mở rộng tích hợp nguồn ảnh vệ tinh độ phân giải siêu cao (Sentinel-2, SPOT) và thuật toán học máy (Machine Learning) để tự động hóa công tác giám sát biến động theo thời gian thực.
  • Các cơ quan quản lý nhà nước và viện nghiên cứu cần sớm thể chế hóa quy trình ứng dụng viễn thám - GIS vào công tác lập quy hoạch đất đai định kỳ, hướng tới xây dựng hệ sinh thái đô thị thông minh và phát triển kinh tế bền vững.