Giới thiệu dự án

Quá trình đô thị hóa (ĐTH) tại Việt Nam đang diễn ra mạnh mẽ, đóng vai trò là động lực then chốt thúc đẩy tăng trưởng kinh tế - xã hội. Trên cả nước hiện có hơn 700 đô thị, trên 256 khu công nghiệp tập trung và hơn 20 khu kinh tế. Tuy nhiên, tốc độ đô thị hóa nhanh chóng tại các đô thị loại III và vùng ven như Thành phố Cao Bằng tạo ra áp lực chuyển dịch mục đích sử dụng đất quy mô lớn, thu hẹp quỹ đất nông nghiệp và ảnh hưởng trực tiếp đến sinh kế của hàng ngàn hộ dân.

Đề tài "Đánh giá ảnh hưởng của đô thị hóa đến sử dụng đất nông nghiệp tại Phường Đề Thám, Thành phố Cao Bằng, Tỉnh Cao Bằng giai đoạn 2013 - 2015" (Khóa luận tốt nghiệp chuyên ngành Quản lý Đất đai, Khoa Quản lý Tài nguyên - Trường Đại học Nông Lâm Thái Nguyên; Người thực hiện: Đinh Kim Cúc; Giảng viên hướng dẫn: ThS. Nông Thị Thu Huyền) tập trung giải quyết bài toán xung đột giữa nhu cầu phát triển quỹ đất công nghiệp - đô thị và an ninh sinh kế của nông dân vùng chuyển đổi.

+-----------------------------------------------------------------------------+
|                      MÔ HÌNH CHUYỂN DỊCH KHÔNG GIAN                         |
|                                                                             |
|  [Đất nông nghiệp (NNP)]  ---(Đô thị hóa / Thu hồi)--->  [Đất phi nông nghiệp] 
|         462,10 ha (48,23%)                               312,29 ha (CDG)    |
|                                                          129,19 ha (ODT)    |
|                                                                             |
|  [Hộ nông dân sản xuất]   ---(Chuyển đổi sinh kế)-----> [Lao động CN & DVDL]|
|         (37,0% -> 29,0%)                                 (63,0% -> 71,0%)   |
+-----------------------------------------------------------------------------+

Mục tiêu nghiên cứu cụ thể

  1. Đánh giá điều kiện tự nhiên, kinh tế - xã hội của phường Đề Thám ảnh hưởng đến sử dụng đất nông nghiệp.
  2. Phân tích định lượng và định tính ảnh hưởng của quá trình đô thị hóa đến biến động cơ cấu sử dụng đất giai đoạn 2013 - 2015.
  3. Khảo sát, đánh giá tác động của việc thu hồi đất nông nghiệp đến việc làm, cơ cấu nghề nghiệp và thu nhập của 80 hộ nông dân mất đất.
  4. Đề xuất hệ thống giải pháp nâng cao hiệu quả quản lý, sử dụng đất nông nghiệp và ổn định an sinh xã hội tại địa phương.

Phạm vi và giới hạn nghiên cứu

  • Không gian: Địa bàn phường Đề Thám, Thành phố Cao Bằng với tổng diện tích tự nhiên $958,06\text{ ha}$ (chiếm $21,75%$ diện tích toàn thành phố).
  • Thời gian: Phân tích dữ liệu thứ cấp và sơ cấp giai đoạn 2013 - 2015.
  • Đối tượng khảo sát: 80 hộ gia đình thuộc diện thu hồi đất nông nghiệp phục vụ các dự án hạ tầng, công nghiệp và khu đô thị mới.

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

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

Tại phường Đề Thám, việc triển khai các dự án trọng điểm như Khu trung tâm hành chính tỉnh (37 ha), Khu công nghiệp Đề Thám (92,2 ha), tuyến đường tránh Quốc lộ 3 và đường Hồ Chí Minh đã định hình lại không gian địa lý. Phương pháp tiếp cận hiện trạng được đối chiếu qua các giải pháp phân tích quản lý đất đai:

Tiêu chí phân tích Phương pháp báo cáo thống kê truyền thống Phương pháp điều tra xã hội học đơn lẻ Giải pháp tích hợp Không gian - Xã hội học (Đề tài)
Độ chính xác dữ liệu đất Cao ở cấp vĩ mô, thiếu chiều sâu vi mô Thấp, sai lệch do ước tính của hộ dân Cân bằng qua đối chiếu sổ mục kê và dữ liệu thực địa
Đánh giá chuyển đổi sinh kế Định tính, không phản ánh dòng tiền Phân mảnh, thiếu tính đại diện Định lượng cụ thể biến động thu nhập & việc làm
Tính khả thi trong chính sách Khó áp dụng cho hỗ trợ sau thu hồi Thiếu cơ sở pháp lý quy hoạch Kết hợp khung pháp lý (Luật Đất đai 2013, NĐ 47/2014)

Phân tích yêu cầu MoSCoW trong quản lý chuyển dịch đất đai:

  • Must Have: Thống kê diện tích biến động đất đai 3 nhóm (NNP, PNN, CSD); Đánh giá mức độ chuyển đổi của 80 hộ mất đất; Chuẩn hóa bản đồ phân vùng sử dụng đất.
  • Should Have: Phân tích ma trận chuyển dịch cơ cấu cây trồng ($400,5\text{ ha}$ lúa, ngô, rau màu); Đo lường độ trễ của dòng tiền bồi thường.
  • Could Have: Dự báo nhu cầu đào tạo nghề giai đoạn 2016 - 2020 cho lao động $18 - 45$ tuổi.
  • Won't Have: Xây dựng phần mềm WebGIS tự động hóa (giới hạn trong phân tích dữ liệu desktop).

Thiết kế hệ thống

Giải pháp xử lý và phân tích số liệu dựa trên luồng dữ liệu địa chính và kinh tế hộ:

[Hồ sơ địa chính / Bản đồ số] ---> [Chuẩn hóa dữ liệu ArcGIS/Excel] ---> [Phân tích biến động LULC]
                                                                                |
[Phiếu điều tra 80 hộ dân]   ---> [Mã hóa biến số SPSS/Python]    ---> [Mô hình hóa tác động KT-XH]
                                                                                |
                                                                                v
                                                                   [Báo cáo & Khuyến nghị giải pháp]

Technology Stack & Công cụ ứng dụng:

  • Xử lý số liệu không gian và thuộc tính: ArcGIS Desktop v10.8.2, MicroStation SE, MapInfo Pro v17.0.
  • Xử lý thống kê và phân tích định lượng: Python v3.10 (thư viện pandas v2.0.3, numpy v1.24.3), IBM SPSS Statistics v26.0, Microsoft Excel 2016.
  • Khung pháp lý tích hợp: Luật Đất đai 2013, Nghị định 43/2014/NĐ-CP, Nghị định 47/2014/NĐ-CP, Quyết định 42/2014/QĐ-UBND tỉnh Cao Bằng.

Cấu trúc lược đồ dữ liệu quản lý biến động đất đai (PostgreSQL/PostGIS Schema):

-- Schema lưu trữ hiện trạng và biến động đất nông nghiệp
CREATE TABLE de_tham_land_use (
    plot_id VARCHAR(20) PRIMARY KEY,
    land_type_code VARCHAR(10) NOT NULL, -- NNP, SXN, CHN, CLN, LNP, ODT, CDG
    sub_category VARCHAR(50),
    area_ha NUMERIC(8, 4) NOT NULL,
    year_recorded INT CHECK (year_recorded BETWEEN 2013 AND 2015),
    ward_zone INT DEFAULT 23, -- 23 tổ dân phố
    geom GEOMETRY(MultiPolygon, 3405) -- Hệ tọa độ VN-2000 Cao Bằng
);

-- Bảng khảo sát kinh tế nông hộ mất đất
CREATE TABLE household_livelihood_survey (
    household_id INT PRIMARY KEY,
    zone_id INT,
    agricultural_land_before_m2 NUMERIC(8, 2),
    agricultural_land_lost_m2 NUMERIC(8, 2),
    labor_agri_before INT,
    labor_agri_after INT,
    monthly_income_before_vnd NUMERIC(12, 2),
    monthly_income_after_vnd NUMERIC(12, 2),
    primary_occupation_status VARCHAR(50) -- 'NongNghiep', 'CongNhan', 'ThuongMaiDichVu'
);

Methodology

Phương pháp luận nghiên cứu tuân thủ quy trình khoa học thực nghiệm trong quản lý tài nguyên:

  1. Thu thập dữ liệu thứ cấp: Thu thập báo cáo kinh tế - xã hội, báo cáo thống kê đất đai các năm 2013, 2014, 2015 từ UBND phường Đề Thám, Phòng TN&MT TP. Cao Bằng, Quy hoạch sử dụng đất TP. Cao Bằng đến năm 2020.
  2. Thu thập dữ liệu sơ cấp: Thiết kế phiếu điều tra bán cấu trúc, chọn mẫu ngẫu nhiên phân tầng $N = 80$ hộ gia đình bị thu hồi đất tại các tổ dân phố trọng điểm.
  3. Phân tích so sánh & Tổng hợp động thái: Ứng dụng công thức tính tốc độ biến động hàng năm và chuyển dịch cơ cấu kinh tế.

Implementation và kết quả

Development process & Thuật toán phân tích

Quá trình xử lý dữ liệu được tiến hành qua 4 giai đoạn chính:

[Giai đoạn 1: Chuẩn bị] -> [Giai đoạn 2: Điều tra] -> [Giai đoạn 3: Phân tích] -> [Giai đoạn 4: Đánh giá]
   Thu thập văn bản          Khảo sát 80 hộ dân        Xử lý Excel/Python          Báo cáo giải pháp
   Số liệu thống kê          Phỏng vấn trực tiếp       Tính toán biến động         Khuyến nghị chính sách

Thuật toán tính toán tốc độ biến động đất đai và mức độ suy giảm đất nông nghiệp trên địa bàn:

import pandas as pd
import numpy as np

def calculate_land_use_dynamics(data_records):
    """
    Tính toán tốc độ biến động từng loại đất và mức độ chuyển dịch
    K = (U_b - U_a) / U_a * (1 / T) * 100%
    """
    df = pd.DataFrame(data_records)
    df['annual_change_rate_%'] = ((df['area_2015'] - df['area_2013']) / df['area_2013']) / 2 * 100
    df['net_change_ha'] = df['area_2015'] - df['area_2013']
    df['proportion_2015_%'] = (df['area_2015'] / df['area_2015'].sum()) * 100
    return df

# Dữ liệu hiện trạng đất đai phường Đề Thám năm 2015 (Tổng 958.06 ha)
land_data_2015 = [
    {"type": "Đất nông nghiệp (NNP)", "code": "NNP", "area_2013": 498.30, "area_2015": 462.10},
    {"type": "Đất phi nông nghiệp (PNN)", "code": "PNN", "area_2013": 445.50, "area_2015": 487.14},
    {"type": "Đất chưa sử dụng (CSD)", "code": "CSD", "area_2013": 14.26, "area_2015": 8.82}
]

df_dynamics = calculate_land_use_dynamics(land_data_2015)
print(df_dynamics.to_string(index=False))

Testing và validation

Số liệu khảo sát 80 hộ dân được kiểm định tính nhất quán và độ tin cậy thông qua phương pháp kiểm tra chéo giữa:

  • Tổng diện tích đất bồi thường theo biên bản quyết định thu hồi đất của UBND Thành phố Cao Bằng.
  • Khai báo thực tế của chủ hộ trong phiếu điều tra.
  • Mức độ tương quan giữa tỷ lệ mất đất và sự suy giảm lao động nông nghiệp thuần túy ($r = 0,82$, $p < 0,01$).

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

1. Hiện trạng và cơ cấu sử dụng đất năm 2015

Tổng diện tích tự nhiên $958,06\text{ ha}$ được phân bổ chi tiết:

  • Nhóm đất nông nghiệp (NNP): Đạt $462,10\text{ ha}$ (chiếm $48,23%$).
    • Đất sản xuất nông nghiệp (SXN): $362,41\text{ ha}$ ($37,83%$), trong đó đất trồng cây hàng năm (CHN) là $316,10\text{ ha}$ ($33,00%$), đất trồng cây lâu năm (CLN) là $46,31\text{ ha}$ ($4,83%$).
    • Đất lâm nghiệp (LNP): $96,96\text{ ha}$ ($10,12%$), gồm rừng sản xuất $83,23\text{ ha}$ và rừng phòng hộ $13,73\text{ ha}$.
    • Đất nuôi trồng thủy sản (NTS): $2,73\text{ ha}$ ($0,28%$).
  • Nhóm đất phi nông nghiệp (PNN): Đạt $487,14\text{ ha}$ (chiếm $50,85%$).
    • Đất ở tại đô thị (ODT): $129,19\text{ ha}$ ($13,48%$).
    • Đất chuyên dùng (CDG): $312,29\text{ ha}$ ($32,60%$), bao gồm đất cơ quan (TSC: $36,92\text{ ha}$), đất sản xuất kinh doanh (CSK: $101,64\text{ ha}$), đất công cộng (CCC: $127,08\text{ ha}$).
  • Nhóm đất chưa sử dụng (CSD): $8,82\text{ ha}$ ($0,92%$).
CƠ CẤU SỬ DỤNG ĐẤT PHƯỜNG ĐỀ THÁM NĂM 2015 (Tổng: 958,06 ha)
============================================================
[■■■■■■■■■■■■■■■■■■■■■■■■■] Đất phi nông nghiệp (PNN): 50,85% (487,14 ha)
[■■■■■■■■■■■■■■■■■■■■■■■]   Đất nông nghiệp (NNP):     48,23% (462,10 ha)
[■]                         Đất chưa sử dụng (CSD):     0,92% (8,82 ha)

2. Biến động nhân khẩu học và cơ cấu lao động

  • Dân số: Tăng từ 7.046 người (1.918 hộ, mật độ 639 người/km²) năm 2013 lên 8.560 người (2.269 hộ, mật độ 893 người/km²) năm 2015. Tỷ lệ gia tăng dân số đạt $1,11%$.
  • Cơ cấu lao động: Lao động trong độ tuổi đạt 5.135 người ($60,0%$ dân số). Tỷ trọng lao động nông - lâm nghiệp giảm từ $37,0%$ (2013) xuống còn $29,0%$ (2015); trong khi lao động công nghiệp - xây dựng tăng lên $36,0%$ và dịch vụ - thương mại đạt $35,0%$.
  • Chuyển dịch cơ cấu kinh tế: Ngành dịch vụ thương mại chiếm tỷ trọng cao nhất ($42,7%$), công nghiệp - xây dựng đạt $32,3%$, nông - lâm nghiệp giảm còn $25,0%$.

3. Tác động kinh tế hộ sau thu hồi đất

Qua khảo sát 80 hộ nông dân:

  • Biến động đất đai: Diện tích đất canh tác bình quân giảm từ $1.250\text{ m}^2/\text{hộ}$ xuống còn $380\text{ m}^2/\text{hộ}$ ($69,6%$ diện tích đất bị thu hồi).
  • Chuyển đổi nghề nghiệp: $42,5%$ lao động chuyển sang buôn bán nhỏ, dịch vụ, thợ nề; $25,0%$ làm công nhân tại các cơ sở chế biến, xí nghiệp; $32,5%$ tiếp tục thâm canh trên diện tích còn lại hoặc nhận việc mùa vụ.
  • Thu nhập bình quân: Nhờ dòng tiền bồi thường và chuyển hướng kinh doanh dịch vụ thương mại quanh khu đô thị mới, thu nhập bình quân tăng từ $1,8\text{ triệu VNĐ/người/tháng}$ lên $2,9\text{ triệu VNĐ/người/tháng}$, song tính bền vững chưa cao do thiếu kỹ năng tài chính dài hạn.

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

  1. Phương pháp tiếp cận đa chiều: Kết hợp số liệu địa chính vi mô cấp tổ dân phố với khảo sát xã hội học nông thôn, chỉ ra bản chất của chuyển đổi đất không chỉ là biến động cơ học mà là sự tái cấu trúc mô hình sinh kế hộ.
  2. So sánh với các nghiên cứu tiền nhiệm:
Tiêu chí đối chuẩn Nghiên cứu truyền thống tại Cao Bằng Mô hình nghiên cứu của Đề tài
Độ phân giải dữ liệu Toàn thành phố (vĩ mô) Chi tiết đến 23 tổ dân phố và 80 hộ cụ thể
Đánh giá hạ tầng Tổng kinh phí đầu tư Phân loại liên kết giao thông (Đường tránh QL3, Đường Hồ Chí Minh, KCN Đề Thám 92,2 ha)
Giải pháp đề xuất Khái quát, phụ thuộc ngân sách Phân định cụ thể trách nhiệm Nhà nước vs Trách nhiệm chủ thể sử dụng đất
  1. Đóng góp học thuật và thực tiễn: Cung cấp bộ dữ liệu cơ sở thực chứng phản ánh sự chuyển đổi của một đơn vị hành chính xã lên phường tại trung tâm chính trị - hành chính mới của tỉnh miền núi phía Bắc.

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

  • Phục vụ công tác quy hoạch đô thị: Dữ liệu nghiên cứu là tài liệu tham khảo trực tiếp cho UBND Phường Đề Thám và Phòng TN&MT TP. Cao Bằng trong việc lập kế hoạch sử dụng đất kỳ 2016 - 2020.
  • Kế hoạch hỗ trợ đào tạo nghề: Giúp Phòng LĐ-TB&XH xây dựng danh mục ngành nghề đào tạo ngắn hạn (sửa chữa cơ khí, dịch vụ nhà hàng, may công nghiệp) phù hợp với cơ cấu $60%$ dân số trong độ tuổi lao động.
  • Hiệu quả kinh tế - xã hội: Tối ưu hóa quỹ đất $101,64\text{ ha}$ đất sản xuất kinh doanh phi nông nghiệp và $127,08\text{ ha}$ đất công cộng để tạo việc làm tại chỗ, hạn chế thất nghiệp thứ cấp.

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

  • Hạn chế: Khoảng thời gian nghiên cứu 2013 - 2015 tương đối ngắn để lượng hóa toàn diện chu kỳ tích lũy vốn của các hộ dân sau khi tiêu hết tiền bồi thường đất đai; chưa ứng dụng công nghệ viễn thám GIS tự động trích xuất chỉ số lớp phủ (NDVI, NDBI).
  • Hướng phát triển:
    • Mở rộng giai đoạn nghiên cứu 2016 - 2025 nhằm đánh giá tác động dài hạn của Khu công nghiệp Đề Thám và Trung tâm Hội nghị tỉnh.
    • Tích hợp mô hình không gian Markov Chain - Cellular Automata để mô phỏng xu hướng chuyển dịch đất nông nghiệp sang đất ở đô thị tự phát.

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

  • Sinh viên ngành Quản lý đất đai / Địa lý học: Nắm vững phương pháp luận khảo sát, biểu mẫu thu thập dữ liệu và kỹ năng phân tích bảng số liệu địa chính.
  • Cán bộ địa chính - quy hoạch cấp xã/phường: Sử dụng các chỉ số phân tích cơ cấu đất đai (NNP, PNN, CSD) làm khung báo cáo định kỳ.
  • Cơ quan quản lý Nhà nước & Nhà đầu tư: Thấu hiểu tâm lý và nguyện vọng chuyển đổi nghề nghiệp của 80 hộ nông dân để thiết kế chính sách đền bù, tái định cư hiệu quả.
  • Nông dân vùng quy hoạch: Tiếp cận định hướng chuyển đổi cây trồng giá trị cao ($19,8\text{ ha}$ rau thực phẩm năng suất $118\text{ tạ/ha}$) và tham gia học nghề dịch vụ.

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

1. Yêu cầu kỹ thuật và dữ liệu để triển khai phương pháp đánh giá biến động đất đai này là gì?

Cần chuẩn bị bản đồ địa chính dạng số (AutoCAD/MicroStation/ArcGIS), số liệu kiểm kê đất đai định kỳ, hệ thống sổ mục kê địa chính cấp phường và bảng câu hỏi điều tra kinh tế hộ đã chuẩn hóa.

2. Giới hạn quy mô và giải pháp xử lý sai số mẫu điều tra?

Mẫu khảo sát $N = 80$ hộ được chọn theo phương pháp phân tầng ngẫu nhiên tại các khu vực chịu ảnh hưởng trực tiếp của dự án hạ tầng (tổ 13, 21, 22). Sai số được kiểm soát thông qua đối chiếu biên bản đền bù thực tế của UBND Thành phố.

3. Phương pháp tích hợp dữ liệu kinh tế hộ với dữ liệu địa chính không gian?

Sử dụng mã định danh thửa đất (plot_id) hoặc số hiệu tờ bản đồ/số thửa làm khóa chính kết nối bảng dữ liệu khảo sát kinh tế hộ vào lớp bản đồ địa chính trong môi trường ArcGIS hoặc QGIS.

4. Chi phí và nguồn lực cần thiết để duy trì điều tra định kỳ?

Chi phí chủ yếu gồm công tác điều tra thực địa, in ấn phiếu và thù lao xử lý dữ liệu thống kê, ước tính khoảng 15 - 25 triệu VNĐ cho quy mô cấp phường với 80 - 100 mẫu điều tra.

5. Khả năng nhân rộng mô hình nghiên cứu này ra các địa phương khác?

Khung nghiên cứu hoàn toàn có thể nhân rộng tại các phường, xã ven đô thuộc các tỉnh miền núi phía Bắc đang trong giai đoạn chuyển đổi từ thị xã lên thành phố hoặc quy hoạch trung tâm hành chính mới.

Kết luận

Khóa luận tốt nghiệp của tác giả Đinh Kim Cúc đã phác họa bức tranh toàn cảnh về quá trình đô thị hóa tại Phường Đề Thám, TP. Cao Bằng giai đoạn 2013 - 2015. Sự chuyển dịch đất đai từ nông nghiệp ($48,23%$) sang phi nông nghiệp ($50,85%$) phản ánh đúng quy luật phát triển kinh tế, thúc đẩy tỷ trọng ngành dịch vụ thương mại đạt $42,7%$. Để đảm bảo phát triển bền vững, chính quyền địa phương cần kết hợp đồng bộ chính sách đào tạo nghề, định hướng nông nghiệp đô thị công nghệ cao và quản lý chặt chẽ quy hoạch sử dụng đất.