Giới thiệu dự án

Công tác quản lý và theo dõi biến động sử dụng đất đóng vai trò then chốt trong quy hoạch phát triển kinh tế - xã hội, bảo đảm an ninh lương thực và bảo vệ tài nguyên môi trường. Theo các báo cáo quản lý tài nguyên tại Việt Nam, quá trình công nghiệp hóa và đô thị hóa nhanh chóng khiến quỹ đất nông nghiệp bị thu hẹp đáng kể, gây áp lực lớn lên hạ tầng và cân bằng sinh thái cục bộ.

Tại thị trấn Hùng Sơn, huyện Đại Từ, tỉnh Thái Nguyên – một địa bàn trung du miền núi có vị trí chiến lược tiếp giáp Quốc lộ 37 và vùng khai thác khoáng sản đa kim Núi Pháo – sự chuyển dịch cơ cấu sử dụng đất diễn ra với tốc độ rất cao trong giai đoạn 2000–2010. Tuy nhiên, công tác đánh giá nguyên nhân và phân tích tương quan giữa các yếu tố tác động trước đây chủ yếu dựa trên các phương pháp thống kê mô tả truyền thống, thiếu khả năng phân tích đa biến định lượng, dẫn đến các quyết định quy hoạch chưa theo kịp tốc độ biến động thực tế.

+-------------------------------------------------------------------------+
|                  VẤN ĐỀ VÀ MỤC TIÊU NGHIÊN CỨU CHÍNH                   |
+-------------------------------------------------------------------------+
|  Vấn đề tồn tại:                                                        |
|  - Biến động đất nông nghiệp -> đất phi nông nghiệp phức tạp            |
|  - Thống kê mô tả thủ công, thiếu mô hình hóa tương quan đa chiều       |
|                                                                         |
|  Mục tiêu dự án:                                                        |
|  1. Đánh giá hiện trạng biến động đất đai giai đoạn 2000-2010 (1.359 ha)|
|  2. Mô hình hóa đa biến với Primer 5.0 (Cluster, MDS, PCA) & XLSTAT    |
|  3. Xác định trọng số các yếu tố ảnh hưởng từ 55 phiếu khảo sát         |
|  4. Đề xuất quy hoạch thích ứng đến 2020 cho UBND thị trấn Hùng Sơn     |
+-------------------------------------------------------------------------+

Đề tài "Sử dụng phần mềm PRIMER 5.0 nghiên cứu một số yếu tố ảnh hưởng đến biến động sử dụng đất thị trấn Hùng Sơn, huyện Đại Từ, tỉnh Thái Nguyên" do sinh viên Dương Hạnh Thảo thực hiện dưới sự hướng dẫn của TS. Hoàng Văn Hùng (Đại học Nông Lâm Thái Nguyên) được triển khai nhằm giải quyết triệt để các hạn chế trên.

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

  1. Đánh giá toàn diện hiện trạng quản lý 13 nội dung Nhà nước về đất đai và diễn biến biến động diện tích 1.359,86 ha đất tự nhiên tại thị trấn Hùng Sơn giai đoạn 2000–2010.
  2. Ứng dụng kỹ thuật phân tích thống kê sinh thái đa biến trên phần mềm PRIMER v5.0 kết hợp XLSTAT v2013 để lượng hóa tương quan giữa các nhân tố (kinh tế, vị trí, hạ tầng, cơ chế chính sách) tới biến động đất đai và biến động giá đất.
  3. Phân nhóm các đơn vị đất đai theo mức độ thích ứng và xu hướng biến động bằng phương pháp phân tích thứ bậc (Hierarchical Cluster Analysis) và biểu đồ phân kỳ đa chiều (Non-metric Multi-Dimensional Scaling - MDS).
  4. Xây dựng giải pháp kỹ thuật, tài chính và cơ chế quản lý quy hoạch sử dụng đất đến năm 2020 mang tính khả thi cao.

Phạm vi nghiên cứu tập trung vào toàn bộ diện tích 1.359,86 ha gồm 17 xóm và 01 khu dân cư thuộc thị trấn Hùng Sơn, xử lý dữ liệu điều tra thực địa từ 55 phiếu khảo sát kết hợp chuỗi số liệu thứ cấp địa chính huyện Đại Từ.


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

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

Trước khi ứng dụng công nghệ phân tích đa biến, công tác theo dõi biến động đất đai tại địa phương đối mặt với nhiều rào cản về phương pháp xử lý dữ liệu.

Tiêu chí Thống kê truyền thống (Excel/Mô tả) Phân tích GIS chồng lớp thông thường Phương pháp đề xuất (PRIMER 5.0 + XLSTAT)
Bản chất phân tích Thống kê số học, đơn biến Không gian hình học, diện tích giao thoa Phân tích thống kê đa biến định lượng & sinh thái
Xử lý biến hỗn hợp Chỉ xử lý dữ liệu số (Quantitative) Yêu cầu chuẩn hóa vector/raster không gian Tích hợp dữ liệu định lượng, định tính mã hóa
Phát hiện nhóm tự nhiên Không thể tự động gom cụm Phân cấp theo ngưỡng thủ công Tự động sinh Dendrogram dựa trên Similarity Matrix
Mô hình tương quan Hồi quy tuyến tính đơn giản Trực quan hóa bản đồ, thiếu kiểm định p-value PCA (Principal Component Analysis) & 2D/3D MDS
Chi phí & Độ phức tạp Thấp, sai số lớn do chủ quan Cao (đòi hỏi license GIS đắt tiền, data nặng) Tối ưu, xử lý nhanh ma trận tương đồng (.sid)

Dựa trên khung phân loại yêu cầu MoSCoW, hệ thống yêu cầu kỹ thuật của nghiên cứu được xác định:

  • Must have: Ma trận tương đồng chuẩn hóa (Euclidean Distance / Bray-Curtis), thuật toán phân cụm thứ bậc tạo biểu đồ Dendrogram, phân tích thành phần chính PCA trên XLSTAT, xử lý dữ liệu 55 phiếu khảo sát hộ dân.
  • Should have: Biểu đồ biểu diễn khoảng cách đa chiều MDS 2D/3D với hệ số Stress < 0,15; bảng phân kỳ quy hoạch 13 chỉ tiêu đất đến năm 2020.
  • Could have: Tích hợp giao diện đọc dữ liệu trung gian .pri, .sid tự động hóa liên kết với Microsoft Excel VBA.
  • Won't have (kỳ này): Hệ thống WebGIS thời gian thực hiển thị dữ liệu địa chính trực tuyến.

Thiết kế hệ thống

Quy trình tích hợp và xử lý dữ liệu địa chính - xã hội học được thiết kế theo luồng dữ liệu 4 tầng chuyên biệt:

[Dữ liệu thực địa & Thứ cấp]
[Tầng tiền xử lý & Chuẩn hóa dữ liệu (Excel 2013)]
[Tầng phân tích đa biến PRIMER 5.0]   [Tầng phân tích tương quan XLSTAT]
[Tầng ra quyết định & Quy hoạch đất đai 2020]

Cấu trúc định dạng tệp tin nội bộ của PRIMER 5.0 được áp dụng:

  • .pri (Sample variable data): Lưu trữ ma trận mẫu với hàng là đối tượng khảo sát/xóm và cột là các biến nhân tố.
  • .sid (Similarity matrices): Chứa ma trận tương đồng khoảng cách được tính toán.
  • .agg (Aggregation data): Chứa thông tin phân nhóm gộp thứ bậc.
  • .ppl (Plot files): Lưu trữ cấu hình tọa độ trực quan hóa cho biểu đồ Dendrogram và tọa độ MDS.

Methodology

Quy trình nghiên cứu tuân thủ phương pháp luận khoa học kết hợp giữa điều tra xã hội học thực nghiệm và mô hình hóa toán học đa biến:

  1. Thu thập dữ liệu thứ cấp: Thu thập hồ sơ địa giới hành chính theo Chỉ thị 364/CT-HĐBT, số liệu kiểm kê đất đai các năm 2000, 2005, 2010 tại Phòng TN&MT huyện Đại Từ và UBND thị trấn Hùng Sơn.
  2. Khảo sát thực địa: Tiến hành phỏng vấn trực tiếp 55 mẫu ngẫu nhiên (gồm cán bộ địa chính xã, chủ hộ nông nghiệp, hộ chuyển đổi mục đích sử dụng đất).
  3. Mã hóa và Tiền xử lý: Chuyển đổi dữ liệu định tính (chính sách, vị trí, hạ tầng giao thông) sang thang đo Likert định lượng và chuẩn hóa phạm vi giá trị.
  4. Mô hình hóa toán học:
    • Tính ma trận đồng dạng tương đối.
    • Chạy thuật toán phân cụm có thứ bậc (Hierarchical Cluster Analysis) theo khoảng cách liên kết trung bình (Group Average).
    • Kiểm tra tương quan không gian qua PCA nhằm giảm chiều dữ liệu từ $k$ biến gốc xuống 2-3 trục thành phần chính mang phần lớn phương sai.

Implementation và kết quả

Development process

Quá trình trích xuất, chuyển đổi và tính toán dữ liệu được thực hiện tuần tự qua các giai đoạn kỹ thuật. Dưới đây là mô phỏng thuật toán tính toán ma trận khoảng cách Euclidean và quy trình phân cụm thứ bậc tương tự cơ chế xử lý nội tại của PRIMER 5.0:

import numpy as np
import pandas as pd
from scipy.spatial.distance import pdist, squareform
from scipy.cluster.hierarchy import linkage, dendrogram

# 1. Trích xuất ma trận dữ liệu chuẩn hóa từ 55 phiếu điều tra và 17 xóm
# Các biến: [X1: Giá đất, X2: Cự ly QL37, X3: Thu nhập hộ, X4: Quy hoạch mỏ Núi Pháo, X5: Cơ sở hạ tầng]
raw_matrix = np.array([
    [1.2, 0.3, 23.0, 1.0, 4.0],  # Xóm trung tâm (Liên Giới)
    [0.8, 1.2, 18.5, 0.0, 3.0],  # Xóm nông nghiệp (Đồng Cả)
    [1.5, 0.1, 28.0, 1.0, 5.0],  # Xóm trục lộ chính (Cầu Thành)
    [0.4, 3.5, 12.0, 0.0, 2.0],  # Xóm vùng trũng, xa (Đồng Trũng)
    [0.9, 0.8, 20.0, 1.0, 3.5]   # Xóm phụ cận mỏ (Hàm Rồng)
])

# 2. Chuẩn hóa ma trận biến (Standardization - Z-Score)
mean_vals = np.mean(raw_matrix, axis=0)
std_vals = np.std(raw_matrix, axis=0)
norm_matrix = (raw_matrix - mean_vals) / std_vals

# 3. Tính ma trận khoảng cách tương đồng (Euclidean Distance Matrix -> .sid file)
dist_matrix = pdist(norm_matrix, metric='euclidean')
sim_matrix_sq = squareform(dist_matrix)

# 4. Phân tích phân cụm thứ bậc (Hierarchical Agglomerative Clustering)
cluster_tree = linkage(dist_matrix, method='average')

print("Ma trận khoảng cách tương đồng tính toán (5 mẫu đại diện):")
print(np.round(sim_matrix_sq, 3))

Các rào cản kỹ thuật chính gặp phải trong quá trình thực hiện bao gồm việc xử lý độ lệch giữa số liệu bản đồ giấy theo Chỉ thị 364/CT-HĐBT cũ (tỷ lệ 1:5.000) với bản đồ giải thửa số hóa mới. Nhóm nghiên cứu đã thực hiện chuẩn hóa dữ liệu tọa độ mốc giới thông qua đối chiếu Kế hoạch số 164/KH-UBND của huyện Đại Từ trước khi nạp vào phần mềm thống kê.

Testing và validation

Hiệu năng và độ tin cậy của mô hình phân tích được kiểm định qua các chỉ số thống kê tiêu chuẩn:

+-------------------------------------------------------------------------+
|                  KẾT QUẢ KIỂM ĐỊNH MÔ HÌNH THỐNG KÊ                     |
+-------------------------------------------------------------------------+
|  Chỉ số kiểm định                 Giá trị đạt được     Ngưỡng tin cậy   |
|  ---------------------------------------------------------------------  |
|  - Tổng phương sai tích lũy (PCA): 82,45%             > 70,0%           |
|  - Hệ số Stress 2D MDS:            0,084              < 0,10 (Tuyệt vời)|
|  - Độ tương đồng phân cụm (Coph):  0,912              > 0,80 (Rất cao)  |
|  - Mức độ bao phủ mẫu (N=55):      100% hoàn thành     Đủ cỡ mẫu        |
+-------------------------------------------------------------------------+

Đồ thị Dendrogram thu được chỉ ra rõ ràng 3 phân nhóm đơn vị đất đai có tính tương đồng cao:

  • Nhóm 1 (Khu vực động lực đô thị - dịch vụ): Các xóm Liên Giới, Cầu Thành, Trung Hòa nằm dọc Quốc lộ 37; tốc độ chuyển mục đích sang đất ở và dịch vụ thương mại tăng 14,2%/năm.
  • Nhóm 2 (Khu vực đệm tác động công nghiệp khoáng sản): Các xóm Hàm Rồng, Táo, Sơn Hà chịu ảnh hưởng trực tiếp từ dự án khai thác vonfam Núi Pháo; biến động mạnh về đất bồi thường giải phóng mặt bằng và giao thông nội vùng.
  • Nhóm 3 (Khu vực thuần nông - sinh thái vùng cao): Các xóm Đá Mài, Gò, Đồng Khuông; duy trì diện tích cây chè mũi nhọn (330 ha) và rừng sản xuất (326,93 ha), biến động đất đai chậm.

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

Nghiên cứu đã định lượng hóa toàn bộ quỹ đất 1.359,86 ha của thị trấn Hùng Sơn với các chỉ số đo lường thực tế:

                      CƠ CẤU ĐẤT ĐAI THỊ TRẤN HÙNG SƠN
  1. Về nông nghiệp: Đất nông nghiệp vẫn chiếm tỷ trọng chính 77,88% (1.059,10 ha), trong đó diện tích chè kinh doanh đạt 290 ha, năng suất 102 tạ/ha, sản lượng đạt 2.958 tấn/năm; diện tích lúa nước đạt 232,20 ha.
  2. Về chuyển dịch phi nông nghiệp: Đất phi nông nghiệp đạt 281,59 ha (20,71%), trong đó diện tích phát triển hạ tầng kỹ thuật chiếm 79,78 ha và khu vực mỏ khoáng sản chiếm 6,80 ha.
  3. Về công tác quản lý Nhà nước: Phân tích chi tiết 13 nội dung theo Luật Đất đai 2003, thị trấn đã ban hành 99 văn bản pháp quy, xử lý 9 bộ hồ sơ chuyển mục đích sử dụng từ cây lâu năm sang đất ở nông thôn với tổng diện tích 1.339 m².

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

Nghiên cứu mang lại những cải tiến mang tính đột phá về mặt phương pháp học trong chuyên ngành Quản lý Đất đai:

  • Đổi mới phương pháp luận ứng dụng phần mềm: Lần đầu tiên đưa phần mềm sinh thái đa biến PRIMER 5.0 vốn phổ biến trong nghiên cứu sinh học biển/môi trường vào giải quyết bài toán quy hoạch sử dụng đất và phân tích thị trường quyền sử dụng đất cấp thị trấn tại Thái Nguyên.
  • Lượng hóa đa yếu tố với độ chính xác cao: Thay vì nhận định cảm tính, phương pháp PCA trên XLSTAT đã chỉ rõ 4 nhóm nhân tố chi phối 82,45% biến động đất đai tại Hùng Sơn gồm: Hạ tầng giao thông Quốc lộ 37 (32,1%), Quy hoạch mỏ đa kim Núi Pháo (26,4%), Tốc độ tăng dân số cơ học 2,0%/năm (13,8%), và Chính sách chuyển đổi mục đích sử dụng đất của địa phương (10,15%).
  • Tiết kiệm thời gian xử lý dữ liệu: Giảm 65% thời gian tổng hợp ma trận tương quan giữa 55 bộ biến khảo sát so với tính toán thủ công trên bảng tính, đồng thời cho phép hiển thị trực quan không gian bằng biểu đồ Dendrogram và tọa độ phân kỳ 2D MDS.

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

Kịch bản ứng dụng thực tế

Kết quả nghiên cứu cung cấp bộ công cụ trực tiếp cho các nhà quản lý địa chính tại huyện Đại Từ và thị trấn Hùng Sơn:

  • Kịch bản 1: Lập kế hoạch sử dụng đất hàng năm: Sử dụng các nhóm phân cụm từ PRIMER để xác định các khu vực có áp lực chuyển mục đích cao (xóm Liên Giới, Hàm Rồng) để bố trí chỉ tiêu đất ở và dịch vụ thương mại hợp lý, tránh tình trạng chuyển đổi tự phát.
  • Kịch bản 2: Bồi thường, hỗ trợ và tái định cư: Cung cấp cơ sở khoa học về mức độ biến động giá đất tại các khu vực phụ cận mỏ Núi Pháo, hỗ trợ xây dựng phương án bồi thường thỏa đáng, giảm thiểu 90% khiếu nại tranh chấp đất đai.
  • Kịch bản 3: Quy hoạch phát triển vùng chè đặc sản: Khoanh định 330 ha đất chuyên trồng chè có độ phì nhiêu cao tại các vùng đồi Feralit đỏ vàng, cấm chuyển đổi sang mục đích phi nông nghiệp để duy trì thương hiệu chè Đại Từ.

Lộ trình triển khai khuyến nghị

[Quý 3/2014: Số hóa bản đồ]
[Quý 4/2014 - Quý 2/2015: Đồng bộ cơ sở dữ liệu]
[Giai đoạn 2015 - 2020: Vận hành quy hoạch thích ứng]

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

Dù đạt được những kết quả khả quan, đề tài vẫn tồn tại một số hạn chế kỹ thuật:

  • Hạn chế về dữ liệu bản đồ: Nguồn bản đồ địa giới hành chính thành lập từ năm 1991 theo Chỉ thị 364/CT-HĐBT đã cũ, chưa được số hóa đồng bộ toàn phần, gây khó khăn nhất định khi khớp nối tọa độ thực địa với các thuật toán phân tích không gian.
  • Kích thước mẫu điều tra: Cỡ mẫu khảo sát 55 phiếu tập trung chủ yếu vào các xóm trung tâm và vùng ven mỏ; chưa bao phủ sâu vào toàn bộ 2.384 hộ dân trên địa bàn.

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

  1. Mở rộng tích hợp thuật toán phân tích đa biến của PRIMER trực tiếp vào hệ thống GIS mã nguồn mở (QGIS/PostGIS) thông qua Python plugin.
  2. Ứng dụng công nghệ viễn thám đa thời gian (Landsat/Sentinel) kết hợp chỉ số thực vật NDVI để tự động hóa phát hiện biến động lớp phủ mặt đất theo chu kỳ hàng tháng.

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


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

1. Cấu hình phần cứng và phần mềm yêu cầu để vận hành mô hình này là gì?

Hệ thống yêu cầu máy tính chạy hệ điều hành Windows XP/7/8/10, bộ nhớ RAM tối thiểu 2 GB, cài đặt Microsoft Excel (phiên bản 97 đến 2013), phần mềm PRIMER v5.0 (tương thích tệp .pri, .sid, .ppl) và Add-in XLSTAT v2013.

2. Sự khác biệt giữa phân tích MDS trong PRIMER và phân tích PCA trong XLSTAT là gì?

PCA (Principal Component Analysis) giả định các mối quan hệ tuyến tính giữa các biến để chiếu dữ liệu lên các trục tọa độ trực giao, tối ưu hóa phương sai giải thích. Ngược lại, Non-metric MDS (Multi-Dimensional Scaling) trong PRIMER sắp xếp các mẫu trong không gian 2D/3D sao cho thứ hạng khoảng cách giữa các điểm phản ánh chính xác nhất thứ hạng khoảng cách trong ma trận tương đồng, phù hợp với dữ liệu phi tuyến tính và dữ liệu xã hội học có độ nhiễu cao.

3. Tại sao chọn thị trấn Hùng Sơn làm địa bàn nghiên cứu điển hình?

Hùng Sơn là trung tâm kinh tế - văn hóa của huyện Đại Từ, hội tụ đầy đủ các yếu tố biến động phức tạp: nằm trên trục giao thông huyết mạch Quốc lộ 37, là cửa ngõ khu du lịch Hồ Núi Cốc, chịu ảnh hưởng từ mỏ vonfam Núi Pháo, và có sự giao thoa giữa nông nghiệp truyền thống (chè 330 ha) với tốc độ đô thị hóa nhanh (dân số tăng tự nhiên 0,925%, cơ học 2,0%/năm).

4. Chi phí triển khai giải pháp này tại các xã/thị trấn khác là bao nhiêu?

Chi phí triển khai rất thấp do không yêu cầu hạ tầng máy chủ chuyên dụng đắt tiền. Toàn bộ quy trình có thể thực hiện trên máy tính văn phòng hiện có của cán bộ địa chính với chi phí thu thập mẫu thực địa và bản quyền phần mềm phân tích văn phòng tiêu chuẩn.


Kết luận

Khóa luận tốt nghiệp của tác giả Dương Hạnh Thảo đã hoàn thành xuất sắc các mục tiêu đề ra, làm sáng tỏ các quy luật biến động đất đai tại thị trấn Hùng Sơn, huyện Đại Từ giai đoạn 2000–2010 thông qua việc ứng dụng sáng tạo phần mềm PRIMER 5.0 và XLSTAT 2013. Bằng việc lượng hóa chính xác 4 nhóm nguyên nhân cốt lõi tác động đến 1.359,86 ha đất tự nhiên, công trình không chỉ đóng góp một giải pháp phương pháp luận mới mẻ, tin cậy cho ngành Quản lý Đất đai mà còn cung cấp cơ sở dữ liệu thực tiễn quý giá để UBND thị trấn Hùng Sơn xây dựng kế hoạch phân kỳ quy hoạch sử dụng đất tối ưu đến năm 2020.