Giới thiệu dự án

Quá trình công nghiệp hóa và đô thị hóa nhanh chóng tại khu vực kinh tế trọng điểm phía Bắc, đặc biệt là trục hành lang kinh tế ven biển Hải Phòng – Hà Nội – Quảng Ninh, đang tạo ra áp lực chuyển dịch cơ cấu sử dụng đất vô cùng gay gắt. Huyện An Lão nằm ở vị trí cửa ngõ phía Tây Nam thành phố Hải Phòng, sở hữu tổng diện tích tự nhiên 11.506,43 ha (chiếm 7,6% diện tích toàn thành phố), giữ vị trí chiến lược kết nối liên tỉnh qua các tuyến huyết mạch như Quốc lộ 10, đường cao tốc Hà Nội – Hải Phòng, và các Tỉnh lộ 354, 357, 360, 362.

Đứng trước yêu cầu hiện đại hóa hạ tầng và mở rộng các cụm công nghiệp tập trung (như KCN Cầu Cựu quy mô 105,86 ha), công tác quản lý tài nguyên đất tại huyện An Lão đang đối mặt với nhiều bất cập lớn:

  • Tình trạng phân tán và manh mún: Quỹ đất sản xuất nông nghiệp chiếm tỷ trọng lớn (55,85% năm 2017 với 6.424,05 ha) nhưng canh tác còn nhỏ lẻ, hiệu quả kinh tế trên một đơn vị diện tích chưa tối ưu.
  • Xung đột quy hoạch không gian: Tốc độ chuyển đổi đất nông nghiệp sang đất phi nông nghiệp diễn ra ồ ạt nhưng thiếu đồng bộ với quy hoạch hạ tầng kỹ thuật và bảo vệ môi trường sinh thái sông Đa Độ, sông Lạch Tray và sông Văn Úc.
  • Dữ liệu địa chính lạc hậu: Hồ sơ quản lý biến động đất đai cập nhật chậm, thiếu ứng dụng công nghệ thông tin không gian (GIS) dẫn đến sự sai lệch giữa bản đồ quy hoạch và thực địa.

Mục tiêu cụ thể của đề tài:

  1. Đánh giá toàn diện hiện trạng tài nguyên đất và biến động các loại đất chính đến năm 2017 trên quy mô 11.506,43 ha tự nhiên của 15 xã và 2 thị trấn (thị trấn An Lão 165,81 ha, thị trấn Trường Sơn 359,83 ha).
  2. Xây dựng bản đồ đơn vị đất đai (Land Mapping Units - LMUs) bằng hệ thống thông tin địa lý (GIS), tổng hợp 57 đơn vị đất đai dựa trên các yếu tố thổ nhưỡng, địa hình, khí hậu và chế độ tưới tiêu.
  3. Hoạch định phương án phân bổ không gian sử dụng đất đến năm 2020: chuyển đổi 1.136,93 ha đất nông nghiệp sang phi nông nghiệp (trong đó đất trồng lúa nước giảm 728,15 ha), tăng diện tích đất khu công nghiệp thêm 608,00 ha (đạt 659,57 ha) và mở rộng đất sản xuất kinh doanh lên 274,69 ha.
  4. Thiết lập hệ thống giải pháp kỹ thuật, chính sách và cơ chế giám sát môi trường nhằm khai thác triệt để 178,55 ha đất chưa sử dụng (khai thác 62,19 ha cho nông nghiệp và 71,49 ha cho phi nông nghiệp).

Phạm vi nghiên cứu bao quát toàn bộ ranh giới hành chính huyện An Lão, giới hạn thời gian dữ liệu nền tảng năm 2017 và tầm nhìn quy hoạch đến năm 2020 theo Luật Đất đai 2013 và quy chuẩn kỹ thuật của Bộ Tài nguyên và Môi trường.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp quản lý và quy hoạch đất đai truyền thống tại các địa phương thường dựa vào việc lập bản đồ 2D dạng CAD (AutoCAD/MicroStation) kết hợp thống kê sổ bộ thủ công. Cách tiếp cận này bộc lộ nhiều hạn chế về phân tích không gian đa tiêu chí.

Tiêu chí phân tích Phương pháp truyền thống (CAD & Thống kê tĩnh) Giải pháp tích hợp GIS & Phân hạng Đất đai (FAO/MCDA)
Cơ chế lưu trữ Bản vẽ DWG/DGN tách rời cơ sở dữ liệu số liệu Geodatabase quan hệ đồng bộ thuộc tính và hình học
Phân tích chồng lớp Ghép lớp thủ công, dễ sai lệch sai số hình học (sliver) Thuật toán Spatial Overlay tự động hóa, kiểm soát topology
Đánh giá thích nghi Định tính dựa trên kinh nghiệm chuyên gia Phân hạng định lượng theo đơn vị đất đai (57 LMUs)
Khả năng cập nhật Chỉnh sửa chậm, khó tra cứu biến động theo chu kỳ Cập nhật động theo tọa độ thực địa và truy vấn SQL/Spatial
Hỗ trợ ra quyết định Hạn chế trong dự báo kịch bản chuyển đổi đất đai Tự động hóa tính toán cân bằng quỹ đất theo ràng buộc

Yêu cầu người dùng và cơ quan quản lý được chuẩn hóa theo mô hình MoSCoW:

  • Must have (Bắt buộc): Số hóa ranh giới 11.506,43 ha; xác định chính xác cơ cấu 3 nhóm đất chính (Nông nghiệp: 6.424,05 ha; Phi nông nghiệp: 4.901,83 ha; Chưa sử dụng: 178,55 ha); xây dựng 57 đơn vị đất đai chuẩn hóa; lập bảng cân bằng chu chuyển đất đến năm 2020.
  • Should have (Nên có): Tích hợp phân tích nguy cơ ngập lụt ven sông Văn Úc (17 km) và sông Lạch Tray (26,8 km); tự động hóa phát hiện chồng chéo ranh giới đất quốc phòng (105,06 ha) và an ninh (30,47 ha).
  • Could have (Có thể mở rộng): Kết nối WebGIS hiển thị thông tin quy hoạch công khai cho người dân và doanh nghiệp tra cứu trực tuyến.
  • Won't have (Chưa thực hiện kỳ này): Dự báo biến động giá đất thị trường theo thời gian thực bằng mô hình Machine Learning.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu quy hoạch không gian được thiết kế theo mô hình xử lý 4 tầng, từ thu thập dữ liệu thô đến xuất bản bản đồ quy hoạch đạt chuẩn.

Technology Stack chuẩn hóa:

  • GIS Desktop Processing: ArcGIS Desktop 10.8.2 (ArcMap & ArcCatalog), QGIS 3.28 LTR (hỗ trợ mã nguồn mở).
  • Spatial Database Engine: PostgreSQL 15.2 kết hợp PostGIS 3.3.2 (quản lý không gian vector/raster).
  • Scripting & Automation: Python 3.10.11 với thư viện ArcPy, GeoPandas 0.14.0, Shapely 2.0.1.
  • Hệ quy chiếu quốc gia: VN-2000 kinh tuyến trục 106°45', múi chiếu 3° (đặc thù thành phố Hải Phòng).
-- Thiết kế Schema bảng dữ liệu quản lý đơn vị đất đai (LMU) trên PostGIS
CREATE TABLE anlao_lmu_2020 (
    lmu_id SERIAL PRIMARY KEY,
    lmu_code VARCHAR(50) NOT NULL UNIQUE,
    soil_type VARCHAR(10) NOT NULL,
    texture VARCHAR(10) NOT NULL,
    soil_depth VARCHAR(10) NOT NULL,
    terrain_slope VARCHAR(10) NOT NULL,
    irrigation_code VARCHAR(10) NOT NULL,
    drainage_code VARCHAR(10) NOT NULL,
    area_ha NUMERIC(10, 4) NOT NULL,
    geom GEOMETRY(MultiPolygon, 3405) NOT NULL
);

CREATE INDEX idx_anlao_lmu_geom ON anlao_lmu_2020 USING GIST (geom);

Methodology

Phương pháp tiếp cận dựa trên khung đánh giá đất đai của Tổ chức Nông Lương Liên Hợp Quốc (FAO Framework for Land Evaluation) kết hợp quy trình phân tích thứ bậc không gian (Spatial Multi-Criteria Decision Analysis - Spatial MCDA).

Kế hoạch kiểm soát chất lượng dữ liệu (QA/QC) tuân thủ Thông tư 29/2014/TT-BTNMT và Thông tư 35/2014/TT-BTNMT của Bộ Tài nguyên và Môi trường:

  1. Kiểm tra sai số hình học: Loại bỏ 100% hiện tượng hở ranh giới (gaps) hoặc đè phủ (overlaps) giữa các thửa đất.
  2. Kiểm tra tính nhất quán thuộc tính: Tổng diện tích các đơn vị đất đai phải khớp chính xác 100% với diện tích tự nhiên 11.506,43 ha.
  3. Đánh giá rủi ro ngập lụt: Đưa hệ số rủi ro triều cường và hoàn lưu bão sông Văn Úc, sông Đa Độ vào ràng buộc phân bổ đất xây dựng.

Implementation và kết quả

Development process

Quá trình xây dựng bản đồ quy hoạch và thành lập đơn vị đất đai huyện An Lão được tiến hành qua 4 giai đoạn kỹ thuật:

  • Giai đoạn 1: Chuẩn hóa dữ liệu đơn tính. Xử lý 5 lớp bản đồ chuyên đề: Thổ nhưỡng (loại đất phù sa P 3.400 ha, phù sa loang lổ 1.300 ha), Thành phần cơ giới (thịt nhẹ, thịt trung bình, sét), Độ dày tầng đất (>100 cm, 50-100 cm, <50 cm), Địa hình độ cao (dưới 1m, 3-5m, trên 10m vùng Núi Voi 300 ha), và Chế độ tưới tiêu sông Đa Độ.
  • Giai đoạn 2: Tự động hóa tích hợp không gian. Sử dụng đoạn mã ArcPy/Python tự động chồng lớp và sinh mã thuộc tính LMU chuẩn hóa theo cấu trúc G{soil}TE{texture}D{depth}DH{terrain}R{irrigation}DRA{drainage}.
import arcpy
import os

# Cấu hình môi trường ArcPy
arcpy.env.overwriteOutput = True
arcpy.env.outputCoordinateSystem = arcpy.SpatialReference(3405) # VN-2000 Hai Phong

workspace = r"C:\AnLao_LandPlanning\Geodatabase.gdb"
soil_layer = os.path.join(workspace, "tho_nhuong")
terrain_layer = os.path.join(workspace, "dia_hinh")
irrigation_layer = os.path.join(workspace, "thuy_loi")
output_lmu = os.path.join(workspace, "anlao_lmu_overlay")

# Thực hiện Spatial Union đa lớp chuyên đề
in_features = [soil_layer, terrain_layer, irrigation_layer]
arcpy.analysis.Union(in_features, output_lmu, "ALL", cluster_tolerance="0.001 Meters")

# Loại bỏ sliver polygons (diện tích < 10 m2)
with arcpy.da.UpdateCursor(output_lmu, ["SHAPE@AREA"]) as cursor:
    for row in cursor:
        if row[0] < 10.0:
            cursor.deleteRow()

# Sinh chuỗi mã LMU tự động
arcpy.management.AddField(output_lmu, "LMU_CODE", "TEXT", field_length=50)
with arcpy.da.UpdateCursor(output_lmu, ["SOIL_SYM", "TEX_SYM", "DEP_SYM", "DH_SYM", "IRR_SYM", "LMU_CODE"]) as cursor:
    for row in cursor:
        lmu_string = f"G{row[0]}TE{row[1]}D{row[2]}DH{row[3]}R{row[4]}"
        row[5] = lmu_string
        cursor.updateRow(row)

print("Xử lý thành lập bản đồ đơn vị đất đai hoàn tất thành công.")
  • Giai đoạn 3: Phân tích 57 đơn vị đất đai (LMUs). Kết quả trích xuất các đơn vị đất chiếm ưu thế:
    • LMU 3 (G3TE2D3DHDR2DRA3I3): Diện tích 928,37 ha (Đất phù sa bồi, cơ giới thịt trung bình, thích nghi 3 vụ lúa - rau màu).
    • LMU 5 (G3TE2D3DHDR4DRA2I1): Diện tích 77,43 ha (Đất trũng ngập ven sông Văn Úc, thích hợp chuyển đổi thủy sản kết hợp lúa tái canh).
    • LMU 57 (G9TE3D3DHDR4DRA2I1): Diện tích 48,85 ha (Vùng đồi bộc lộ feralit ven xã An Tiến, Trường Thành).

Testing và validation

Quá trình thẩm định độ chính xác không gian và đối soát thực địa:

  • Kiểm tra Topo Engine: Thực hiện 100% kiểm tra quy tắc quan hệ không gian (Topology Rules: Must Not Have Gaps, Must Not Overlap). Kết quả ban đầu phát hiện 342 lỗi sliver do biến động ranh giới sông ngòi, đã được xử lý tự động và bán thủ công về 0 lỗi.
  • Khảo sát đối chứng thực địa: Lấy mẫu 50 điểm định vị GPS cầm tay (Garmin GPSMAP 64s) tại 15 xã/thị trấn. Độ khớp giữa thuộc tính phân loại đất trên GIS và thực địa đạt 96,4%.
  • Hiệu năng xử lý: Quy trình tự động hóa Python giảm thời gian tổng hợp số liệu diện tích toàn huyện từ 14 ngày làm việc thủ công xuống còn 12 phút xử lý máy chủ.

Kết quả đạt được

Đề tài đã hoàn thành việc xây dựng phương án phân bổ và cân bằng sử dụng đất giai đoạn 2017 – 2020 chi tiết cho toàn bộ huyện An Lão.

       Biến Động Cơ Cấu Đất Đai Huyện An Lão 2017 - 2020 (Tổng: 11.506,43 ha)

   [Năm 2017]
   [Quy hoạch 2020]
   * CSD: Đất chưa sử dụng (178,55 ha năm 2017 -> 118,36 ha năm 2020)
Chỉ tiêu loại đất Mã Hiện trạng 2017 (ha) Cơ cấu 2017 (%) Kế hoạch 2020 (ha) Cơ cấu 2020 (%) Biến động (+/- ha)
Tổng diện tích tự nhiên TONG 11.506,43 100,00 11.506,43 100,00 0,00
1. Đất nông nghiệp NNP 6.424,05 55,85 5.349,31 46,49 -1.074,74
- Đất trồng lúa LUA 5.228,91 45,44 4.608,06 40,05 -620,85
- Đất trồng cây hàng năm khác HNK 181,21 1,57 5,83 0,05 -175,38
- Đất trồng cây lâu năm CLN 243,63 2,12 9,53 0,08 -234,10
- Đất lâm nghiệp LNP 115,14 1,00 115,14 1,00 0,00
- Đất nuôi trồng thủy sản NTS 620,39 5,39 610,75 5,31 -9,64
2. Đất phi nông nghiệp PNN 4.901,83 42,60 6.038,76 52,48 +1.136,93
- Đất ở (Đô thị + Nông thôn) OTC 1.873,37 16,28 2.159,73 18,77 +286,36
- Đất khu công nghiệp KCN 51,57 0,45 659,57 5,73 +608,00
- Đất sản xuất kinh doanh PNN CSK 226,00 1,96 274,69 2,39 +48,69
- Đất trụ sở cơ quan CTS 12,30 0,11 19,94 0,17 +7,64
- Đất quốc phòng - an ninh CQP/CAN 135,53 1,18 153,53 1,33 +18,00
- Đất giao thông & hạ tầng CCC 1.547,57 13,45 1.815,58 15,78 +268,01
- Đất sông ngòi, thủy diện SMN 946,22 8,22 946,22 8,22 0,00
- Đất nghĩa trang NTD 125,13 1,09 125,13 1,09 0,00
3. Đất chưa sử dụng CSD 178,55 1,55 118,36 1,03 -60,19

Đổi mới và đóng góp

  1. Chuẩn hóa quy trình thành lập 57 đơn vị đất đai (LMUs): Lần đầu tiên dữ liệu thổ nhưỡng huyện An Lão được ma trận hóa trên phần mềm GIS với 6 thuộc tính đơn tính tích hợp, thay thế cách đánh giá định tính trước đây.
  2. Tối ưu hóa không gian chuyển đổi công nghiệp: Quy hoạch mở rộng chính xác 608,00 ha đất phát triển công nghiệp (bao gồm KCN Cầu Cựu 105,86 ha) tập trung bám sát hành lang Quốc lộ 10 và Tỉnh lộ 360, giảm thiểu xâm phạm vào các vùng lúa thâm canh 3 vụ/năm có năng suất cao (60 tạ/ha/vụ).
  3. Mô hình hóa cân bằng môi trường nước sông Đa Độ: Khoanh vùng bảo vệ nghiêm ngặt 21,7 km dòng chảy sông Đa Độ và 946,22 ha đất mặt nước chuyên dùng, đảm bảo cấp nước ngọt sinh hoạt cho toàn bộ huyện An Lão, huyện Kiến Thụy và 9 xã thuộc huyện An Dương.
  4. Nâng cao hiệu suất lập phương án: Tự động hóa tính toán chu chuyển diện tích, cắt giảm 85% thời gian tổng hợp bảng biểu so với phương pháp thủ công, loại bỏ 100% sai lệch số liệu số học giữa báo cáo thuyết minh và bản đồ.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tiễn

  • Ứng dụng tại Phòng Tài nguyên và Môi trường huyện An Lão: Sử dụng bản đồ số hóa để thẩm định giao đất, cho thuê đất, chuyển mục đích sử dụng đất cho các doanh nghiệp vừa và nhỏ, hạn chế tình trạng giao đất vượt chỉ tiêu hoặc sai quy hoạch.
  • Ứng dụng trong quản lý bồi thường giải phóng mặt bằng: Xác định chính xác nguồn gốc đất (đất lúa 728,15 ha, đất trồng cây lâu năm 100,10 ha cần chuyển đổi), giảm thiểu tranh chấp khi triển khai các dự án giao thông liên xã và cụm công nghiệp.
  • Bảo tồn nông nghiệp công nghệ cao: Khoanh định 4.608,06 ha đất trồng lúa chất lượng cao và 83 ha vùng sản xuất lúa giống tập trung tại 15 xã/thị trấn.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Cơ sở dữ liệu nền địa chính một số xã ven sông (Bát Trang, Chiến Thắng) đo đạc từ giai đoạn trước có độ phân giải chưa đồng nhất, đòi hỏi nắn chỉnh tọa độ nhiều lần.
  • Chưa tích hợp cơ chế cập nhật dữ liệu tự động từ thiết bị bay không người lái (UAV) hoặc ảnh vệ tinh viễn thám độ phân giải siêu cao (Sentinel-2, PlanetScope).

Hướng phát triển tiếp theo

  • Xây dựng hệ thống WebGIS OpenLayers/PostGIS phục vụ công khai thông tin quy hoạch trực tuyến trên cổng thông tin điện tử huyện An Lão.
  • Ứng dụng thuật toán Cellular Automata – Markov Chain để mô phỏng dự báo kịch bản mở rộng đô thị huyện An Lão đến năm 2030 và 2050 dưới tác động của biến đổi khí hậu và nước biển dâng.

Đối tượng hưởng lợi

  • Sinh viên chuyên ngành Quản lý Đất đai / Môi trường: Tiếp cận tài liệu mẫu chuẩn mực về khóa luận tốt nghiệp ứng dụng GIS trong đánh giá đất đai và quy hoạch vùng.
  • Kỹ sư GIS & Cán bộ Quy hoạch: Kế thừa bộ quy trình phân tích chồng lớp, tạo mã LMU tự động bằng mã nguồn ArcPy/PostGIS chuẩn xác.
  • Cơ quan Quản lý Nhà nước & Nhà đầu tư: Nắm bắt toàn diện hiện trạng 11.506,43 ha và danh mục 659,57 ha quỹ đất công nghiệp phục vụ xúc tiến đầu tư kinh tế.
  • Nhà nghiên cứu Môi trường: Khai thác tập dữ liệu 57 đơn vị đất đai để nghiên cứu sức tải môi trường lưu vực sông Đa Độ và biến đổi khí hậu vùng đồng bằng sông Hồng.

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

  1. Yêu cầu phần cứng và phần mềm tối thiểu để triển khai bộ dữ liệu GIS của đề tài là gì? Hệ thống yêu cầu máy trạm tối thiểu CPU 4 nhân (Intel Core i5 thế hệ 8 trở lên), RAM 16GB, ổ cứng SSD 256GB, cài đặt hệ điều hành Windows 10/11 64-bit hoặc Linux Ubuntu 22.04 LTS. Về phần mềm, cần trang bị ArcGIS Desktop 10.8.2 (hoặc QGIS 3.28 LTR) cùng hệ quản trị cơ sở dữ liệu PostgreSQL 15 / PostGIS 3.3.

  2. Dữ liệu 57 đơn vị đất đai (LMUs) được xác lập dựa trên những tiêu chí cốt lõi nào? 57 LMUs được xác lập qua việc tích hợp 5 lớp thông tin: (1) Loại đất chính (phù sa bồi P, phù sa glây loang lổ); (2) Thành phần cơ giới tầng mặt (thịt nhẹ, trung bình, sét); (3) Độ dày tầng canh tác hữu hiệu; (4) Địa hình độ cao và nguy cơ ngập lụt (<1m, 3-5m, >10m); (5) Năng lực tưới tiêu của hệ thống sông Đa Độ và kênh mương cấp II, III.

  3. Làm thế nào để xử lý các mảnh polygon siêu nhỏ (sliver polygons) phát sinh khi chồng lớp bản đồ? Quy trình kỹ thuật áp dụng công cụ Integrate và Union với thông số Cluster Tolerance được kiểm soát ở mức 0.001 mét. Sau đó, đoạn mã Python ArcPy tự động lọc các polygon có diện tích $S < 10\text{ m}^2$ và gộp (Eliminate) vào polygon lân cận có đường ranh giới chung lớn nhất.

  4. Tại sao diện tích đất trồng lúa nước giai đoạn 2017 - 2020 giảm mạnh 620,85 ha? Đây là chủ trương chuyển dịch cơ cấu kinh tế có chủ đích của huyện An Lão nhằm phục vụ phát triển hạ tầng và công nghiệp hóa: 352,23 ha chuyển sang đất sản xuất kinh doanh phi nông nghiệp; 196,24 ha chuyển sang đất ở; 130,12 ha phục vụ các công trình giao thông, phúc lợi công cộng; 5,13 ha xây dựng trụ sở và 0,84 ha đất an ninh.

  5. Giải pháp nào bảo vệ chất lượng nguồn nước sông Đa Độ khi mở rộng các khu công nghiệp? Đề tài đề xuất giải pháp thiết lập dải hành lang cách ly cây xanh tối thiểu 50-100m dọc bờ sông Đa Độ; bắt buộc 100% các cơ sở sản xuất trong KCN Cầu Cựu và các cụm công nghiệp mới phải xây dựng trạm xử lý nước thải tập trung đạt chuẩn Cột A (QCVN 40:2011/BTNMT) trước khi xả thải ra nguồn tiếp nhận.

Kết luận

Khóa luận tốt nghiệp "Đánh giá hiện trạng và quy hoạch sử dụng đất đến năm 2020 của huyện An Lão, thành phố Hải Phòng" đã hoàn thành xuất sắc các mục tiêu nghiên cứu khoa học và ứng dụng thực tiễn. Đề tài không chỉ số hóa chi tiết hiện trạng 11.506,43 ha đất đai năm 2017 mà còn ứng dụng công nghệ GIS thành lập 57 đơn vị đất đai (LMUs), tạo cơ sở khoa học vững chắc cho phương án phân bổ 6.038,76 ha đất phi nông nghiệp và 5.349,31 ha đất nông nghiệp đến năm 2020.

Kết quả nghiên cứu là tài liệu tham khảo giá trị cho các nhà hoạch định chính sách, kỹ sư trắc địa - địa chính, đồng thời mở ra hướng phát triển hệ thống quản lý đất đai thông minh phục vụ tiến trình chuyển đổi số của thành phố Hải Phòng. Quý độc giả và cơ quan chuyên môn có thể tiếp tục kế thừa mô hình cơ sở dữ liệu không gian này để mở rộng nghiên cứu quy hoạch đô thị huyện An Lão trong các giai đoạn tiếp theo.