Tổng quan nghiên cứu

Quá trình đô thị hóa diễn ra mạnh mẽ tại Việt Nam trong hơn hai thập kỷ qua đã tạo ra áp lực to lớn đối với tài nguyên đất đai. Thành phố Hải Phòng, với vị thế là 1 trong 3 trung tâm đô thị lớn nhất cả nước, đang chứng kiến tốc độ chuyển dịch đất đai diễn ra nhanh chóng, đặc biệt tại các khu vực ven đô. Huyện Kiến Thụy là địa bàn ngoại thành điển hình với quá trình mở rộng khu dân cư và công nghiệp, khiến quỹ đất nông nghiệp bị thu hẹp đáng kể qua các chu kỳ kiểm kê đất đai định kỳ 5 năm một lần. Thực tế này đòi hỏi cơ quan quản lý phải có phương pháp giám sát liên tục, trực quan và có độ tin cậy cao thay vì chỉ dựa vào các phương pháp thống kê thủ công tốn kém chi phí.

Mục tiêu cốt lõi của nghiên cứu là xây dựng quy trình tích hợp tư liệu viễn thám đa thời gian, hệ thống thông tin địa lý và mô hình kết hợp chuỗi Markov cùng mạng tế bào tự động nhằm đánh giá biến động lớp phủ mặt đất huyện Kiến Thụy giai đoạn 2000 - 2010, đồng thời mô phỏng dự báo xu thế biến đổi không gian đến năm 2015 và 2020. Phạm vi không gian nghiên cứu bao trùm toàn bộ ranh giới hành chính huyện Kiến Thụy với 18 đơn vị cấp xã, tập trung vào 4 nhóm đối tượng chính gồm đất nông nghiệp, đất ở, đất bằng chưa sử dụng và mặt nước. Nghiên cứu mang ý nghĩa thực tiễn sâu sắc, cung cấp cơ sở khoa học giúp nâng cao độ chính xác trong công tác lập quy hoạch sử dụng đất với mức độ tương thích không gian đạt trên 85%, hỗ trợ đắc lực cho mục tiêu phát triển bền vững của địa phươ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 hệ thống lý thuyết chuẩn mực về khoa học địa lý và mô hình hóa không gian:

Khái niệm lớp phủ mặt đất được tiếp cận theo định nghĩa của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc, xác định đây là tập hợp tất cả các đối tượng vật chất nhân tạo hoặc tự nhiên quan sát được trên bề mặt Trái Đất, phân biệt rõ ràng với khái niệm mục đích sử dụng đất.

Lý thuyết chuỗi Markov do nhà toán học Andrei Markov đề xuất năm 1907 và được Kolmogorov mở rộng năm 1936, đóng vai trò xác định xác suất chuyển dịch trạng thái giữa các loại hình lớp phủ theo thời gian rời rạc dựa trên giả định trạng thái tương lai chỉ phụ thuộc vào trạng thái hiện tại.

Mô hình mạng tế bào tự động được khởi xướng bởi Von Neumann năm 1966 và ứng dụng vào không gian địa lý bởi Tobler năm 1979, cho phép mô phỏng sự lan truyền không gian dựa trên trạng thái của các ô lân cận theo cấu trúc lân cận Von Neumann 4 ô hoặc Moore 8 ô.

Phương pháp đánh giá đa chỉ tiêu được tích hợp nhằm xác định các trọng số và ngưỡng thích hợp cho từng loại hình biến động đất đai dưới tác động của các yếu tố kinh tế, xã hội và hạ tầng giao thông.

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

Nghiên cứu sử dụng nguồn tư liệu ảnh vệ tinh đa thời gian bao gồm Landsat MSS, Landsat TM với độ phân giải không gian 30m và ảnh SPOT 3, SPOT 5 có độ phân giải từ 15m đến 20m, thu nhận qua các mốc thời gian năm 2000, 2005 và 2010. Bộ cảm biến quang học ghi nhận năng lượng bức xạ điện từ trong dải sóng nhìn thấy và cận hồng ngoại từ 400nm đến 25cm với 256 mức xám.

Cỡ mẫu giải đoán được thiết lập với hơn 300 điểm mẫu phân bố đồng đều trên địa bàn huyện, áp dụng phương pháp chọn mẫu ngẫu nhiên phân tầng kết hợp đối soát tọa độ thực địa bằng thiết bị định vị toàn cầu và bản đồ hiện trạng sử dụng đất các năm 2005, 2010.

Phương pháp phân tích chính là phân loại có kiểm định theo thuật toán xác suất cực đại kết hợp kỹ thuật so sánh sau phân loại trên nền tảng phần mềm Idrisi Andes 15 và mô-đun Land Change Modeler. Lý do lựa chọn tổ hợp phương pháp này là vì thuật toán xác suất cực đại khai thác tối đa đặc trưng phản xạ phổ của từng đối tượng tự nhiên, trong khi sự kết hợp giữa mô-đun Markov và mạng tế bào tự động giúp khắc phục triệt để nhược điểm thiếu yếu tố vị trí không gian của mô hình Markov thuần túy, mang lại bức tranh dự báo trực quan và định lượng chính xác theo từng pixel raster.

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

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

Kết quả giải đoán và phân tích dữ liệu giai đoạn 2000 - 2010 đã làm rõ quy luật biến động lớp phủ mặt đất trên địa bàn huyện Kiến Thụy:

Thứ nhất, giai đoạn 2000 - 2005 ghi nhận diện tích đất nông nghiệp, đặc biệt là đất trồng lúa và hoa màu, chiếm ưu thế tuyệt đối với tỷ lệ trên 68% diện tích tự nhiên, nhưng bắt đầu xu hướng giảm dần khoảng 2,5% do sự xuất hiện của các công trình công cộng và mở rộng đất ở nông thôn.

Thứ hai, giai đoạn 2005 - 2010 chứng kiến tốc độ chuyển dịch đất đai diễn ra mạnh mẽ nhất, diện tích đất nông nghiệp giảm hơn 4,8%, nhường chỗ cho đất chuyên dùng, khu dân cư mới và diện tích nuôi trồng thủy sản nước lợ tại các vùng trũng.

Thứ ba, kết quả mô phỏng dự báo đến năm 2020 chỉ ra rằng diện tích đất ở và đất phi nông nghiệp sẽ tăng trưởng khoảng 15% đến 20% so với năm 2010, tập trung chủ yếu dọc theo các tuyến đường giao thông liên xã và khu vực giáp ranh nội thành Hải Phòng.

Thứ tư, dải rừng ngập mặn phòng hộ ven biển tại xã Đại Hợp duy trì diện tích tương đối ổn định với độ biến động dưới 1,2%, đóng vai trò quan trọng trong việc bảo vệ đê điều và cân bằng sinh thái địa phương.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự suy giảm liên tục của đất canh tác nông nghiệp xuất phát từ tác động kép của việc triển khai Quy hoạch phát triển kinh tế - xã hội thành phố Hải Phòng đến năm 2020 và nhu cầu nâng cấp hạ tầng dân sinh. Quá trình đô thị hóa tự phát kết hợp với chuyển dịch cơ cấu mùa vụ đã thúc đẩy các hộ dân chuyển đổi từ trồng lúa sang đào ao thả tôm, cá tại các vùng ven sông, ven biển.

Các dữ liệu biến động trong nghiên cứu được thể hiện trực quan qua bảng ma trận chuyển dịch trạng thái diện tích hai chiều, đồng thời trực quan hóa bằng chuỗi bản đồ raster phân giải 30m. Bản đồ phân ngưỡng mức độ thích hợp được xây dựng chi tiết cho từng loại đối tượng như giao thông, thổ cư và mặt nước. Khi đối chiếu với các đồ án quy hoạch tổng thể của địa phương, mô hình mô phỏng đã phản ánh chính xác xu hướng mở rộng không gian đô thị hóa tập trung. Kết quả kiểm chứng độc lập giữa bản đồ mô phỏng năm 2010 và bản đồ phân loại thực tế năm 2010 cho thấy độ tương đồng đạt chỉ số Kappa trên 0,86, khẳng định độ tin cậy vượt trội của thuật toán tích hợp so với các phương pháp dự báo xu thế thống kê đơn biến truyền thống.

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

Dựa trên kết quả dự báo biến động lớp phủ đến năm 2020, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm tối ưu hóa công tác quản lý tài nguyên đất:

Thiết lập ranh giới bảo vệ nghiêm ngặt tối thiểu 40% diện tích đất nông nghiệp chuyên canh lúa nước có độ phì nhiêu cao, nghiêm cấm việc chuyển đổi mục đích sử dụng đất tràn lan sang đất ở và đất công nghiệp; cơ quan chủ trì thực hiện là Ủy ban nhân dân huyện Kiến Thụy phối hợp cùng Sở Tài nguyên và Môi trường thành phố Hải Phòng với lộ trình hoàn thiện trước năm 2015.

Xây dựng hệ thống cơ sở dữ liệu số và ứng dụng công nghệ viễn thám định kỳ theo chu kỳ 1 đến 2 năm một lần nhằm giám sát biến động lớp phủ, thay thế cho phương thức điều tra thống kê thủ công 5 năm một lần; mục tiêu giảm 50% chi phí điều tra thực địa và rút ngắn 70% thời gian cập nhật bản đồ; đơn vị triển khai là Phòng Tài nguyên và Môi trường huyện trong giai đoạn 2014 - 2017.

Quy hoạch bảo tồn và phục hồi dải rừng ngập mặn ven biển tại khu vực xã Đại Hợp và vùng cửa sông, duy trì diện tích rừng phòng hộ trên 250 hecta nhằm tăng cường khả năng chống chịu bão lũ, hạn chế xói mòn bờ biển và giảm thiểu rủi ro biến đổi khí hậu; cơ quan thực hiện là Ban Quản lý rừng phòng hộ phối hợp với chính quyền các xã ven biển từ năm 2014 đến năm 2020.

Nâng cao năng lực chuyên môn cho 100% cán bộ địa chính cấp xã về kỹ năng vận hành hệ thống thông tin địa lý và khai thác ảnh vệ tinh trong quản lý hồ sơ địa chính; đơn vị chủ quản là Sở Khoa học và Công nghệ phối hợp cùng các trường đại học chuyên ngành tổ chức đào tạo trong giai đoạn 2015 - 2018.

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

Nghiên cứu này là tài liệu tham khảo giá trị cho nhiều nhóm đối tượng hoạt động trong lĩnh vực quản lý và nghiên cứu:

Cơ quan quản lý nhà nước về đất đai và môi trường: Cán bộ thuộc Sở Tài nguyên và Môi trường, Phòng Tài nguyên và Môi trường cấp huyện có thể sử dụng quy trình xử lý ảnh và ma trận chuyển dịch để xây dựng phương án điều chỉnh quy hoạch sử dụng đất sát với thực tiễn.

Giảng viên, nghiên cứu sinh và học viên cao học: Những người nghiên cứu chuyên ngành Địa chính, Quản lý đất đai, Hệ thống thông tin địa lý và Viễn thám có thể tham khảo khung phương pháp kết hợp mô hình Markov và mạng tế bào tự động làm nền tảng phát triển các đề tài chuyên sâu.

Các đơn vị tư vấn quy hoạch đô thị và nông thôn: Chuyên gia quy hoạch có thể ứng dụng trực tiếp các bản đồ dự báo không gian để mô phỏng sự phát triển của hệ thống giao thông, điểm dân cư nông thôn và phân bổ mạng lưới hạ tầng kỹ thuật.

Tổ chức nghiên cứu bảo tồn sinh thái và biến đổi khí hậu: Các tổ chức phi chính phủ và viện nghiên cứu môi trường có thể khai thác dữ liệu biến động dải rừng ngập mặn và mặt nước để đánh giá tính dễ bị tổn thương sinh thái vùng duyên hải.

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

Mô hình kết hợp Markov và mạng tế bào tự động có ưu thế gì so với mô hình chuỗi Markov đơn thuần? Mô hình chuỗi Markov đơn thuần chỉ tính toán được xác suất và tổng lượng diện tích biến động giữa các trạng thái qua thời gian mà không định vị được vị trí không gian cụ thể. Việc tích hợp mạng tế bào tự động bổ sung yếu tố tương tác không gian lân cận và các điều kiện ràng buộc địa lý, giúp xác định chính xác vị trí chuyển đổi của từng pixel trên bản đồ.

Tại sao nghiên cứu lại ưu tiên sử dụng ảnh Landsat và SPOT thay vì ảnh vệ tinh siêu cao như QuickBird hay Ikonos? Ảnh vệ tinh Landsat và SPOT có độ phân giải không gian từ 15m đến 30m, hoàn toàn đáp ứng độ chính xác cho việc thành lập bản đồ lớp phủ cấp huyện tỷ lệ 1:25.000 đến 1:50.000. Đồng thời, các tư liệu này có chu kỳ chụp liên tục qua nhiều thập kỷ với chi phí hợp lý, phù hợp cho phân tích chuỗi thời gian dài từ năm 2000 đến 2010.

Độ tin cậy của kết quả dự báo biến động lớp phủ mặt đất đến năm 2020 được kiểm chứng như thế nào? Quy trình dự báo được kiểm chứng thông qua bước mô phỏng ngược hiện trạng năm 2010 dựa trên chuỗi dữ liệu 2000 - 2005. Kết quả mô phỏng năm 2010 sau đó được so sánh chéo với bản đồ phân loại thực tế năm 2010 bằng chỉ số Kappa đạt trên 0,86, chứng minh mô hình có độ tin cậy toán học vững chắc trước khi tiến hành dự báo đến năm 2020.

Phương pháp đánh giá đa chỉ tiêu đóng vai trò gì trong cấu trúc mô hình dự báo? Phương pháp đánh giá đa chỉ tiêu đóng vai trò chuẩn hóa các yếu tố tác động tự nhiên, kinh tế, xã hội như khoảng cách tới đường giao thông, trung tâm dân cư và nguồn nước thành các bản đồ phân ngưỡng thích hợp. Các bản đồ này trở thành trọng số ràng buộc không gian, đảm bảo quá trình chuyển dịch tế bào tự động diễn ra đúng quy luật thực tế.

Kết quả của đề tài có thể nhân rộng sang các địa phương khác ngoài huyện Kiến Thụy không? Quy trình công nghệ và cấu trúc thuật toán được xây dựng trong nghiên cứu có tính phổ quát cao, hoàn toàn có thể áp dụng cho các huyện ngoại thành, thị xã hoặc các đô thị ven biển khác tại Việt Nam có tốc độ đô thị hóa nhanh và đặc điểm phân hóa sinh thái tương đồng.

Kết luận

Đề tài đã hoàn thành xuất sắc các mục tiêu nghiên cứu và mang lại những đóng góp khoa học nổi bật:

  • Hệ thống hóa toàn diện cơ sở lý luận về ứng dụng tích hợp công nghệ viễn thám, hệ thống thông tin địa lý và mô hình Markov - Cellular Automata trong nghiên cứu biến động bề mặt đất.
  • Xây dựng thành công chuỗi bản đồ hiện trạng lớp phủ huyện Kiến Thụy các năm 2000, 2005, 2010 với độ chính xác phân loại đạt trên 85%.
  • Lập ma trận chuyển dịch định lượng phản ánh rõ xu hướng giảm diện tích đất nông nghiệp và sự gia tăng nhanh chóng của đất ở, đất chuyên dùng.
  • Hoàn thiện mô hình dự báo không gian giai đoạn 2010 - 2020, làm rõ quy luật phát triển đô thị hóa tập trung dọc theo các tuyến hành lang giao thông chính.
  • Đề xuất hệ thống 4 nhóm giải pháp quy hoạch khả thi gắn liền với lộ trình thực hiện cụ thể cho chính quyền địa phương.

Đóng góp lớn nhất của công trình là cung cấp một công cụ hỗ trợ ra quyết định hiện đại, thay thế phương pháp thống kê truyền thống bằng mô hình hóa không gian số. Trong giai đoạn tiếp theo từ năm 2015 đến năm 2020, các nhà quản lý và nghiên cứu cần mở rộng tích hợp thêm các kịch bản nước biển dâng và dữ liệu viễn thám radar để hoàn thiện hệ thống dự báo. Hãy khai thác ngay khung phương pháp này để nâng cao hiệu quả quản trị tài nguyên đất đai tại địa phương của bạn.