Giới thiệu dự án

Quá trình công nghiệp hóa, hiện đại hóa (CNH-HĐH) và đô thị hóa tại Việt Nam đang diễn ra với tốc độ nhanh chóng, kéo theo nhu cầu chuyển đổi mục đích sử dụng đất quy mô lớn. Tại tỉnh Thái Bình – một trong những trọng điểm nông nghiệp thuộc Đồng bằng sông Hồng với diện tích tự nhiên đạt $158.461\text{ ha}$, trong đó đất sản xuất nông nghiệp chiếm tới $57,39%$ ($90.941\text{ ha}$) và cơ cấu lao động nông nghiệp chiếm $27%$ năm 2022 – công tác thu hồi đất nông nghiệp (THĐNN) phục vụ phát triển các khu công nghiệp, cụm công nghiệp và hạ tầng giao thông đô thị đang diễn ra mạnh mẽ. Năm 2022, GRDP của tỉnh Thái Bình đạt $63.000\text{ tỷ đồng}$ (tăng trưởng $9,52%$), khu vực Công nghiệp - Xây dựng tăng trưởng $14,09%$, đặt ra sức ép lớn lên quỹ đất nông nghiệp truyền thống.

+-----------------------------------------------------------------------------------+
|                           BỐI CẢNH KINH TẾ THÁI BÌNH                              |
|  - Tổng diện tích đất tự nhiên: 158.461 ha                                       |
|  - Đất sản xuất nông nghiệp: 90.941 ha (57,39%)                                   |
|  - Tăng trưởng GRDP (2022): 9,52% (Quy mô: 63.000 tỷ VND)                        |
|  - Tỷ trọng Công nghiệp - Xây dựng: 42,9% GRDP | Lao động nông nghiệp: 27%       |
+-----------------------------------------------------------------------------------+

Vấn đề nghiên cứu (Problem Statement)

Việc chuyển đổi tư liệu sản xuất cốt lõi (đất đai) sang vốn tài chính bồi thường một lần đặt ra các rủi ro dài hạn đối với sinh kế và phúc lợi bền vững của nông hộ:

  • Mất nguồn sinh kế truyền thống: Người nông dân mất tư liệu sản xuất chính, dẫn đến suy giảm nguồn thu nhập ổn định từ trồng trọt, chăn nuôi.
  • Rào cản chuyển đổi nghề nghiệp: Tỷ lệ lao động nông thôn có trình độ chuyên môn kỹ thuật hạn chế, gặp khó khăn khi gia nhập thị trường lao động phi nông nghiệp hoặc làm việc tại các khu công nghiệp địa phương.
  • Rủi ro phân bổ dòng tiền bồi thường: Nguồn tiền đền bù chủ yếu được sử dụng cho tiêu dùng ngắn hạn, xây dựng nhà ở thay vì tái đầu tư sinh kế, làm gia tăng nguy cơ tái nghèo và bất bình đẳng kinh tế.
  • Tổn thương mạng lưới an sinh: Sự suy giảm phúc lợi phi vật chất (môi trường sống, ô nhiễm không khí, áp lực hạ tầng và biến động cấu trúc văn hóa - xã hội làng xã).

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

  1. Hệ thống hóa cơ sở lý luận: Làm rõ khung phân tích tác động của chính sách THĐNN đến phúc lợi hộ gia đình theo cách tiếp cận đa chiều (vật chất, việc làm, an sinh xã hội, môi trường).
  2. Khảo sát và lượng hóa thực trạng: Thu thập số liệu sơ cấp từ $N = 372$ hộ gia đình bị thu hồi đất tại 4 địa bàn trọng điểm của tỉnh Thái Bình (Đông Hưng, Hưng Hà, TP. Thái Bình, Tiền Hải).
  3. Mô hình hóa kinh tế lượng: Ứng dụng mô hình hồi quy Logit thứ bậc (Ordered Logit Model) nhằm định lượng xác suất thay đổi mức thu nhập và các nhân tố quyết định đến phúc lợi hộ dân sau khi bị thu hồi đất.
  4. Đề xuất hệ thống giải pháp chính sách: Kiến nghị nhóm giải pháp khả thi về đào tạo nghề, hỗ trợ tiếp cận tín dụng, hoàn thiện cơ chế bồi thường - hỗ trợ - tái định cư và nâng cao năng lực quản lý đất đai cấp cơ sở.

Phương pháp tiếp cận giải pháp

Nghiên cứu kết hợp phương pháp nghiên cứu định tính (phỏng vấn bán cấu trúc, tổng quan tài liệu thứ cấp) và định lượng chuyên sâu (thống kê mô tả, kiểm định phi tham số, mô hình kinh tế lượng Ordered Logit). Cách tiếp cận này khắc phục hạn chế của các phân tích thống kê thuần túy, cho phép xác định rõ độ co giãn và mức độ tác động biên (marginal effects) của từng yếu tố chính sách.

Kết quả kỳ vọng và chỉ số đo lường

  • Bộ dữ liệu vi mô chuẩn hóa gồm 372 quan sát với hơn 20 biến số kinh tế - xã hội.
  • Mô hình định lượng đạt mức ý nghĩa thống kê ($p < 0,05$), kiểm định thỏa mãn giả định đường hồi quy song song (Parallel Regression Assumption).
  • Xác định chính xác trọng số ảnh hưởng của: tỷ lệ thu nhập phi nông nghiệp, quy mô lao động, thời gian thu hồi đất, trình độ đào tạo và khả năng tiếp cận vốn tín dụng đến phúc lợi hộ gia đình.

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

  • Không gian: 4 đơn vị hành chính cấp huyện đại diện cho các vùng kinh tế của Thái Bình: Huyện Đông Hưng ($82\text{ hộ}$ - $22,04%$), Huyện Hưng Hà ($85\text{ hộ}$ - $22,85%$), TP. Thái Bình ($101\text{ hộ}$ - $27,15%$), Huyện Tiền Hải ($104\text{ hộ}$ - $27,96%$).
  • Thời gian: Dữ liệu khảo sát sơ cấp thu thập trong giai đoạn 2022 - 2023, đánh giá biến động trước và sau thu hồi đất.
  • Nội dung: Tập trung vào phúc lợi hộ gia đình nông thôn bị thu hồi đất nông nghiệp phục vụ các dự án phát triển kinh tế - xã hội vì lợi ích quốc gia, công cộng theo Luật Đất đai 2013.

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

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

Các nghiên cứu trước đây tại Việt Nam và quốc tế đã tiếp cận vấn đề thu hồi đất dưới nhiều góc độ khác nhau:

Phương pháp tiếp cận Ưu điểm Nhược điểm / Khoảng trống
Thống kê mô tả & So sánh trước - sau (Nguyễn Thị Diễn, 2012; Bùi Thị Minh Nguyệt, 2017) Đơn giản, trực quan, dễ dàng theo dõi biến động cơ cấu nghề nghiệp và thu nhập trung bình. Không kiểm soát được các yếu tố nhiễu đồng thời; không đo lường được tác động biên của từng biến độc lập.
Hồi quy OLS tuyến tính cổ điển (Phan Văn Phúc et al., 2023) Định lượng được tương quan tuyến tính giữa diện tích đất mất đi và mức giảm thu nhập. Vi phạm giả định khi biến phụ thuộc (thu nhập theo nhóm/mức độ thỏa dụng phúc lợi) mang bản chất thứ bậc, phân nhóm rời rạc.
Mô hình Hồi quy Logit thứ bậc (Cameron & Trivedi, 2005 - Phương pháp đề xuất) Xử lý chính xác dữ liệu phụ thuộc dạng thứ bậc ($1 \to 5$), tối ưu hóa hàm hợp lý cực đại (MLE), giải thích rõ xác suất tích lũy. Yêu cầu kiểm định chặt chẽ giả định tỷ lệ tỷ số chênh không đổi (Proportional Odds Assumption).

Yêu cầu nghiên cứu theo ma trận MoSCoW

  • Must have (Bắt buộc): Bộ dữ liệu khảo sát 372 hộ thực tế; mô hình Ordered Logit ước lượng bằng phần mềm STATA; phân tích so sánh chi tiết cơ cấu nghề nghiệp, thu nhập, tín dụng và bồi thường.
  • Should have (Cần có): Đo lường mức độ hài lòng đa khía cạnh (hạ tầng, giao thông, y tế, giáo dục, môi trường, an ninh trật tự).
  • Could have (Có thể có): So sánh đối chiếu với các mô hình thu hồi đất tại các tỉnh lân cận thuộc Đồng bằng sông Hồng (Hải Dương, Hưng Yên).
  • Won't have (Chưa thực hiện kỳ này): Dữ liệu bảng dài hạn (Longitudinal Panel Data) theo dõi hộ gia đình trong chu kỳ 10-15 năm.

Thiết kế hệ thống mô hình hóa

+-----------------------------------------------------------------------------------------+
|                               QUY TRÌNH MÔ HÌNH HÓA ĐỊNH LƯỢNG                          |
|                                                                                         |
|  [Thu thập mẫu sơ cấp]      [Làm sạch & Mã hóa biến]       [Ước lượng Ordered Logit]     |
|   (N = 372 phiếu hỏi)   -->  - STATA 17.0 / Python     -->  - Ước lượng hợp lý cực đại   |
|   - 4 huyện trọng điểm       - Kiểm tra Missing, Outliers   - Tọa độ cắt (Cut points)   |
|                                                                                         |
|                                      |                                                  |
|                                      v                                                  |
|  [Đề xuất chính sách 4 nhóm] <-- [Đánh giá & Kiểm định]                                 |
|   1. Đào tạo chuyển đổi nghề       - Brant Test (Kiểm định Odds)                        |
|   2. Tín dụng sinh kế              - Pseudo R-squared, Log-Likelihood                   |
|   3. Bồi thường & Tái định cư      - Kiểm tra đa cộng tuyến (VIF)                       |
|   4. Quản lý đất đai cơ sở                                                              |
+-----------------------------------------------------------------------------------------+

Công nghệ và công cụ sử dụng

  • STATA 17.0 MP: Công cụ xử lý kinh tế lượng vi mô, ước lượng mô hình hồi quy ologit, tính toán hiệu ứng biên margins, và thực hiện kiểm định brant.
  • Python 3.10 (Pandas, Statsmodels, Scipy, Seaborn): Dùng để pipeline hóa quá trình tiền xử lý dữ liệu, trực quan hóa biểu đồ và kiểm tra phân phối phân bố mẫu.
  • LaTeX / Excel Data Engine: Xây dựng bảng biểu và chuẩn hóa hệ thống công thức toán kinh tế.

Thiết kế cơ sở dữ liệu và cấu trúc biến (Data Schema)

Biến phụ thuộc ($Y$ - Thunhap) phản ánh mức thu nhập hộ gia đình sau thu hồi đất được lượng hóa thành thang đo thứ tự 5 bậc:

  • $Y = 1$: Thu nhập từ $0 - 3\text{ triệu VNĐ/tháng}$
  • $Y = 2$: Thu nhập từ $3 - 5\text{ triệu VNĐ/tháng}$
  • $Y = 3$: Thu nhập từ $5 - 10\text{ triệu VNĐ/tháng}$
  • $Y = 4$: Thu nhập từ $10 - 15\text{ triệu VNĐ/tháng}$
  • $Y = 5$: Thu nhập trên $20\text{ triệu VNĐ/tháng}$
Tên biến Ký hiệu mô hình Loại biến Thang đo / Đơn vị Kỳ vọng dấu
Tuổi của chủ hộ Tuoi Liên tục Năm tuổi $(-)$
Thời gian thu hồi đất Thoi_gian Liên tục Năm $(+)$
Tham gia hoạt động XH Thamgia Giả (Dummy) $1 = \text{Có}, 0 = \text{Không}$ $(+)$
Tỷ lệ thu nhập phi nông nghiệp Phi_Nongnghiep Liên tục $%$ trong tổng thu nhập $(+)$
Tình trạng sức khỏe Suc_khoe Định danh $1 = \text{Tốt}, 2 = \text{Bình thường}, 3 = \text{Yếu}, 4 = \text{Bệnh}$ $(-)$
Số lao động trong hộ Laodong Liên tục Người $(+)$
Thuê thêm đất canh tác Thue_dat Giả (Dummy) $1 = \text{Có}, 0 = \text{Không}$ $(+)$
Đánh giá hỗ trợ TĐC Danh_gia Định danh/Thứ bậc Mức độ đánh giá chính sách $(+)$
Hộ có vay vốn tín dụng Vayvon Giả (Dummy) $1 = \text{Có}, 0 = \text{Không}$ $(+)$
Quy mô tiền đền bù Tiendenbu Liên tục Triệu VNĐ $(+)$

Phương pháp luận (Methodology)

Mô hình Logit thứ bậc (Ordered Logit) giả định tồn tại một biến tiềm ẩn liên tục $Y^*$ (Continuous Latent Variable) chưa quan sát được:

$$Y_i^* = X_i \beta + \varepsilon_i, \quad \varepsilon_i \sim \text{Logistic}(0, \pi^2/3)$$

Mối liên hệ giữa biến tiềm ẩn $Y_i^*$ và biến thứ bậc quan sát được $Y_i$ thông qua các điểm cắt $\alpha_1 < \alpha_2 < \alpha_3 < \alpha_4$:

$$Y_i = \begin{cases} 1 & \text{nếu } Y_i^* \le \alpha_1 \ 2 & \text{nếu } \alpha_1 < Y_i^* \le \alpha_2 \ 3 & \text{nếu } \alpha_2 < Y_i^* \le \alpha_3 \ 4 & \text{nếu } \alpha_3 < Y_i^* \le \alpha_4 \ 5 & \text{nếu } Y_i^* > \alpha_4 \end{cases}$$

Xác suất tích lũy để hộ gia đình $i$ rơi vào nhóm thu nhập $j$ hoặc thấp hơn:

$$\Pr(Y_i \le j \mid X_i) = \frac{\exp(\alpha_j - X_i \beta)}{1 + \exp(\alpha_j - X_i \beta)} = \Lambda(\alpha_j - X_i \beta), \quad j = 1, 2, 3, 4$$

Xác suất rơi vào từng phân loại cụ thể:

$$\Pr(Y_i = j \mid X_i) = \Pr(Y_i \le j \mid X_i) - \Pr(Y_i \le j-1 \mid X_i)$$

Thực thi và kết quả

Quy trình triển khai xử lý dữ liệu

Quá trình phân tích thực nghiệm được cấu trúc hóa thông qua quy trình xử lý dữ liệu chuẩn trên STATA và Python:

* ====================================================================
* KHÓA LUẬN TỐT NGHIỆP: ƯỚC LƯỢNG MÔ HÌNH ORDERED LOGIT CHO PHÚC LỢI NÔNG HỘ
* Địa bàn: Tỉnh Thái Bình (N = 372)
* ====================================================================

clear all
set more off

* 1. Nhập và kiểm tra dữ liệu sơ cấp
use "thaibinh_land_welfare_2023.dta", clear
describe
summarize Thunhap Tuoi Thoi_gian Thamgia Phi_Nongnghiep Suc_khoe Laodong Thue_dat Danh_gia Vayvon Tiendenbu

* 2. Kiểm tra đa cộng tuyến giữa các biến độc lập
regress Thunhap Tuoi Thoi_gian Thamgia Phi_Nongnghiep Suc_khoe Laodong Thue_dat Danh_gia Vayvon Tiendenbu
estat vif

* 3. Ước lượng mô hình Hồi quy Logit thứ bậc (Ordered Logit)
ologit Thunhap Tuoi Thoi_gian Thamgia Phi_Nongnghiep i.Suc_khoe Laodong Thue_dat Danh_gia Vayvon Tiendenbu, robust

* 4. Tính toán Tỷ số chênh (Odds Ratios) và Hiệu ứng biên (Marginal Effects)
ologit, or
margins, dydx(*) predict(outcome(5))
margins, dydx(*) predict(outcome(1))

* 5. Kiểm định giả định hồi quy song song (Brant Test)
brant, detail

Triển khai song song mô hình trên Python (statsmodels.miscmodels.ordinal_model) để kiểm tra độ tin cậy của tham số:

import numpy as np
import pandas as pd
from statsmodels.miscmodels.ordinal_model import OrderedModel

# Load empirical survey dataset
df = pd.read_csv("thaibinh_survey_processed.csv")

# Define dependent and independent variables
formula = "Thunhap ~ Tuoi + Thoi_gian + Thamgia + Phi_Nongnghiep + C(Suc_khoe) + Laodong + Thue_dat + Danh_gia + Vayvon + Tiendenbu"

# Fit Ordered Logit Model
mod_ologit = OrderedModel.from_formula(formula, data=df, distr='logit')
res_ologit = mod_ologit.fit(method='bfgs', maxiter=1000)

print(res_ologit.summary())

Kiểm định và đánh giá độ tin cậy mô hình

  1. Kiểm tra hiện tượng đa cộng tuyến (Multicollinearity Test): Hệ số phóng đại phương sai $\text{VIF} < 2,15$ cho toàn bộ các biến độc lập (trung bình $\text{Mean VIF} = 1,42 < 5$), khẳng định mô hình không xảy ra hiện tượng đa cộng tuyến nghiêm trọng.
  2. Kiểm định giả định đường hồi quy song song (Brant Test): Giá trị $\chi^2$ toàn bộ mô hình không có ý nghĩa thống kê ở mức $5%$ ($p > 0,10$), chứng minh giả định Proportional Odds thỏa mãn hoàn toàn.
  3. Độ phù hợp của mô hình: Chỉ số $\text{Pseudo } R^2 = 0,218$, giá trị $\text{LR } \chi^2(10) = 142,65$ ($p = 0,0000$), chứng minh các biến độc lập giải thích có ý nghĩa cho sự biến động mức thu nhập của hộ gia đình.
+-----------------------------------------------------------------------------------+
|                        KẾT QUẢ KIỂM ĐỊNH MÔ HÌNH KINH TẾ LƯỢNG                    |
|  - Tổng số quan sát (Sample Size): N = 372                                        |
|  - Log-Likelihood = -342.18                                                       |
|  - LR chi2(10) = 142,65 (p-value = 0.0000 ***)                                    |
|  - Pseudo R2 (McFadden) = 0,2184                                                  |
|  - Mean VIF = 1.42 (Không có đa cộng tuyến)                                       |
|  - Brant Test p-value = 0.184 (Thỏa mãn giả định Proportional Odds)               |
+-----------------------------------------------------------------------------------+

Kết quả phân tích thực nghiệm

1. Chuyển dịch cơ cấu nghề nghiệp và việc làm

Khảo sát chỉ ra bước ngoặt chuyển dịch lao động: Tỷ lệ lao động thuần nông nghiệp giảm mạnh từ $68,4%$ trước thu hồi xuống còn $21,5%$ sau thu hồi. Tỷ lệ lao động chuyển sang làm việc trong các khu công nghiệp, xưởng may mặc cơ khí tăng lên $38,2%$, lao động kinh doanh dịch vụ buôn bán nhỏ chiếm $24,1%$, và lao động tự do/thời vụ chiếm $16,2%$.

2. Kết quả ước lượng kinh tế lượng Ordered Logit

  • Tỷ lệ thu nhập phi nông nghiệp (Phi_Nongnghiep): Có hệ số dương ($\beta = 0,0342$, $p < 0,01$, $\text{OR} = 1,035$). Hộ có tỷ lệ thu nhập từ thương mại, dịch vụ, công nghiệp càng cao thì khả năng chuyển dịch lên các bậc thu nhập cao ($> 15 - 20\text{ triệu/tháng}$) tăng tương ứng $3,5%$ với mỗi $1%$ tăng thêm của tỷ trọng thu nhập phi nông.
  • Số lượng lao động (Laodong): Tác động tích cực đáng kể ($\beta = 0,4215$, $p < 0,01$). Hộ gia đình có thêm 1 lao động trong độ tuổi làm việc làm tăng tỷ số chênh chuyển lên nhóm thu nhập cao hơn gấp $1,52$ lần ($\text{OR} = 1,524$).
  • Vay vốn tín dụng (Vayvon): Tác động thuận chiều ở mức ý nghĩa $5%$ ($\beta = 0,6120$, $p = 0,021$). Các hộ tiếp cận vốn ngân hàng chính sách hoặc thương mại để mở rộng xưởng cơ khí, dịch vụ ăn uống có xác suất đạt thu nhập trên 10 triệu/tháng cao hơn $84,4%$ so với hộ không vay vốn.
  • Tuổi của chủ hộ (Tuoi): Tác động ngược chiều ($\beta = -0,0281$, $p < 0,05$). Chủ hộ lớn tuổi gặp nhiều rào cản thích ứng với công việc mới, khiến xác suất rơi vào nhóm thu nhập thấp ($< 3 - 5\text{ triệu/tháng}$) tăng lên.
  • Quy mô tiền đền bù (Tiendenbu): Có hệ số dương nhưng tác động biên tương đối thấp ($\beta = 0,0012$, $p < 0,10$), cho thấy tiền đền bù nếu không đi kèm phương án chuyển đổi sinh kế bền vững chỉ mang lại hiệu ứng thu nhập tức thời trong ngắn hạn.

3. Đo lường mức độ hài lòng phúc lợi đa chiều

Tiêu chí phúc lợi Rất hài lòng / Hài lòng Bình thường Không hài lòng / Bất bình
Cơ sở hạ tầng đường sá giao thông $76,34%$ $18,28%$ $5,38%$
Cơ sở vật chất trường học, trạm y tế $68,82%$ $24,46%$ $6,72%$
Đời sống vật chất tổng thể $48,12%$ $34,41%$ $17,47%$
Đời sống tinh thần, giao lưu cộng đồng $42,20%$ $38,71%$ $19,09%$
Chất lượng môi trường không khí & tiếng ồn $28,49%$ $31,18%$ $40,33%$
Mức giá đền bù và chính sách hỗ trợ $31,18%$ $33,60%$ $35,22%$

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

  1. Đổi mới phương pháp luận kinh tế lượng vi mô: Nghiên cứu vượt ra khỏi cách tiếp cận mô tả đơn thuần bằng việc áp dụng chuẩn mực mô hình Hồi quy Logit thứ bậc (Ordered Logit) với đầy đủ kiểm định chẩn đoán (Brant test, VIF, Marginal effects), giúp lượng hóa xác suất chuyển dịch giữa 5 tầng mức thu nhập của hộ nông dân.
  2. Tiếp cận phúc lợi toàn diện (Multidimensional Welfare Framework): Không chỉ đánh giá phúc lợi dựa trên thu nhập tiền tệ, đề tài tích hợp đánh giá phúc lợi xã hội theo tiêu chuẩn HDI và ILO, bao gồm chất lượng hệ thống y tế - giáo dục, cơ sở hạ tầng nông thôn mới, an ninh trật tự và biến đổi chất lượng môi trường tự nhiên.
  3. Định lượng tác động chuyển đổi sinh kế: Nghiên cứu chứng minh bằng thực nghiệm: Thu nhập phi nông nghiệp đóng vai trò là "chất đệm" giúp hộ gia đình tăng trưởng thu nhập bền vững sau khi mất đất sản xuất (hệ số tác động $\text{OR} = 1,035$, $p < 0,01$).
  4. Đối chiếu thực nghiệm với các công trình tiêu biểu:
Công trình đối chiếu Địa bàn Phương pháp Kết quả khác biệt được bổ sung từ khóa luận
Lê Ngọc Phương Quý et al. (2021) Đô Lương, Nghệ An Thống kê so sánh mô tả Bổ sung mô hình kinh tế lượng đo lường hiệu ứng biên của quy mô lao động và tiếp cận tín dụng đến mức độ phục hồi thu nhập.
Nguyễn Thị Diễn et al. (2012) Hưng Yên Khung sinh kế bền vững DFID Lượng hóa cụ thể tác động của cấu trúc tuổi chủ hộ và kiểm định giả định phân phối thứ bậc trong điều kiện Đồng bằng sông Hồng.
Brahima Coulibaly et al. (2020) Sebougou, Mali Khảo sát định tính ven đô Chứng minh trong bối cảnh Việt Nam, đầu tư hạ tầng công cộng đi kèm dự án công nghiệp bù đắp một phần mức độ suy giảm phúc lợi đất đai.

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

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

  • Ủy ban Nhân dân cấp Huyện (Ban Bồi thường & GPMB): Sử dụng các biến số định lượng để xây dựng phương án phân nhóm đối tượng thụ hưởng chính sách: Hộ có chủ hộ trên 50 tuổi cần chế độ trợ cấp ổn định đời sống dài hạn; hộ có lao động trẻ cần gói tài trợ học nghề tại các cơ sở đào tạo liên kết với khu công nghiệp.
  • Ngân hàng Chính sách Xã hội Chi nhánh Thái Bình: Thiết lập gói tín dụng ưu đãi chuyên biệt phục vụ các hộ bị thu hồi từ $30%$ diện tích đất nông nghiệp trở lên, giảm thiểu yêu cầu tài sản thế chấp truyền thống (do đất đã bị thu hồi).
+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH TRIỂN KHAI CHÍNH SÁCH ĐỀ XUẤT                     |
|                                                                                   |
|  [Giai đoạn 1: Chuẩn bị & Điều tra] (Tháng 1 - 3)                                 |
|  - Rà soát hồ sơ địa chính, phân tầng nhóm tuổi & cấu trúc lao động hộ dân.      |
|  - Khảo sát nhu cầu đào tạo nghề nghiệp và nhu cầu vay vốn sản xuất.             |
|                                                                                   |
|  [Giai đoạn 2: Thực thi bồi thường & Hỗ trợ chuyển đổi] (Tháng 4 - 9)             |
|  - Chi trả bồi thường minh bạch theo khung giá sát thị trường (Đ.114 Luật ĐĐ).   |
|  - Mở lớp đào tạo kỹ thuật ngắn hạn, kết nối tuyển dụng tại KCN Tiền Hải/Đông Hưng|
|  - Giải ngân vốn vay ưu đãi phát triển thương mại - dịch vụ làng nghề.             |
|                                                                                   |
|  [Giai đoạn 3: Giám sát & Tái thiết lập an sinh] (Tháng 10 - 24)                  |
|  - Quan trắc chất lượng môi trường không khí, nguồn nước khu vực dự án.          |
|  - Đánh giá định kỳ chỉ số phục hồi thu nhập và bảo trợ xã hội nhóm yếu thế.      |
+-----------------------------------------------------------------------------------+

Đánh giá hiệu quả kinh tế - xã hội (ROI & Cost-Benefit Analysis)

  • Hiệu quả ngân sách: Việc chuyển dịch một hộ sang thu nhập phi nông nghiệp ổn định giúp tăng đóng góp thuế nội địa cấp cơ sở từ $1,5 - 2,2\text{ lần}$ so với sản xuất lúa thuần túy.
  • Giảm chi phí an sinh: Phòng ngừa tình trạng khiếu kiện đất đai kéo dài giúp tiết kiệm $25 - $30%$ chi phí hành chính và đẩy nhanh tiến độ giải phóng mặt bằng các dự án trọng điểm từ $6 - 8\text{ tháng}$.

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

Hạn chế nghiên cứu

  1. Tính chất dữ liệu cắt ngang (Cross-sectional Data): Dữ liệu thu thập tại một thời điểm (năm 2022 - 2023) phản ánh các biến số dựa trên hồi tưởng của người dân đối với thời kỳ trước thu hồi, có thể tồn tại sai số nhớ lại (Recall Bias).
  2. Phạm vi địa lý: Mặc dù chọn 4 huyện tiêu biểu đại diện cho 4 vùng kinh tế, mẫu khảo sát vẫn chưa bao phủ toàn bộ 8 huyện/thành phố của tỉnh Thái Bình (chưa bao gồm Quỳnh Phụ, Thái Thụy, Kiến Xương, Vũ Thư).
  3. Môi trường tự nhiên: Đo lường tác động môi trường chủ yếu dựa trên cảm nhận chủ quan qua thang đo Likert, chưa tích hợp các chỉ số quan trắc kỹ thuật lý - hóa (nồng độ bụi mịn PM2.5, chỉ số BOD/COD nguồn nước mặt quanh khu công nghiệp).

Hướng nghiên cứu tiếp theo

  • Xây dựng mô hình dữ liệu bảng (Panel Data) theo dõi dài hạn mẫu 372 hộ trong giai đoạn 3 năm, 5 năm và 10 năm sau thu hồi đất.
  • Tích hợp mô hình kinh tế lượng không gian (Spatial Econometrics) để phân tích tác động lan tỏa (Spillover Effects) từ các khu công nghiệp đến các xã vùng ven lân cận.
  • Ứng dụng viễn thám GIS và chỉ số thực vật NDVI kết hợp dữ liệu kinh tế lượng vi mô để giám sát biến động quỹ đất nông nghiệp tự động.

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

+-----------------------------------------------------------------------------------------+
|                                    ĐỐI TƯỢNG HƯỞNG LỢI                                  |
|                                                                                         |
|  [Sinh viên & Giảng viên]       [Chuyên viên Phân tích]        [Nhà Quản lý Chính sách]  |
|  - Khung lý luận kinh tế        - Source code STATA/Python     - Căn cứ thực tiễn ban   |
|    chính trị đất đai chuẩn      - Phương pháp Ordered Logit      hành Quyết định bồi    |
|  - Dữ liệu khảo sát chuẩn        hoàn chỉnh & Brant test         thường, quy hoạch TĐC   |
+-----------------------------------------------------------------------------------------+
  • Sinh viên & Học viên cao học: Tài liệu tham khảo toàn diện về phương pháp luận nghiên cứu kinh tế chính trị ứng dụng, kỹ thuật thiết kế bảng hỏi bán cấu trúc và phương pháp định lượng nâng cao.
  • Chuyên viên phân tích dữ liệu & Kinh tế lượng: Tiếp cận trọn vẹn quy trình kiểm định mô hình Ordered Logit trên STATA và Python, bao gồm các phương pháp kiểm tra vi phạm giả định và tính toán Marginal Effects.
  • Cơ quan quản lý Nhà nước (Sở TN&MT, UBND tỉnh Thái Bình): Bằng chứng thực nghiệm định lượng rõ ràng phục vụ công tác sửa đổi bảng giá đất, điều chỉnh chính sách đào tạo nghề và lập quy hoạch tái định cư đồng bộ.
  • Các tổ chức nghiên cứu chính sách & Phát triển nông thôn: Mô hình tham chiếu tin cậy để triển khai các chương trình an sinh xã hội tại các địa phương có tốc độ đô thị hóa nhanh tại Việt Nam.

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

1. Tại sao đề tài sử dụng mô hình Ordered Logit thay vì mô hình hồi quy tuyến tính OLS thông thường?

Biến phụ thuộc "Thu nhập của hộ gia đình" sau thu hồi đất được mã hóa dưới dạng các nhóm thứ bậc ($1 \to 5$), phản ánh các khoảng thu nhập rời rạc và có thứ tự nhưng khoảng cách giữa các bậc không hoàn toàn bằng nhau. Sử dụng hồi quy tuyến tính OLS cổ điển trong trường hợp này sẽ dẫn đến các ước lượng chệch (biased), không hiệu quả và có thể dự báo xác suất nằm ngoài khoảng $[0, 1]$. Mô hình Ordered Logit tối ưu hóa hàm hợp lý cực đại (MLE) và giả định phân phối Logistic, đảm bảo tính nhất quán và vững của các ước lượng tham số.

2. Các yêu cầu kỹ thuật để tái lập (reproduce) kết quả mô hình này là gì?

  • Môi trường STATA: Phiên bản STATA 15.0 trở lên (khuyến nghị STATA 17.0/18.0), cài đặt bổ sung package brant thông qua lệnh ssc install spost9_legacy hoặc spost13.
  • Môi trường Python: Python $\ge 3.8$, cài đặt các thư viện lõi: pandas>=1.3.0, statsmodels>=0.13.0, numpy>=1.20.0, seaborn>=0.11.0.
  • Cấu hình phần cứng: RAM tối thiểu 4GB, vi xử lý 2 nhân (Core i3 trở lên), dung lượng ổ cứng khả dụng 500MB cho việc lưu trữ các tệp log và kết quả trung gian.

3. Làm thế nào để xử lý khi mô hình Ordered Logit vi phạm giả định hồi quy song song (Parallel Regression Assumption)?

Khi kiểm định Brant trả về giá trị $p < 0,05$ (vi phạm giả định tỷ lệ tỷ số chênh không đổi), nhà nghiên cứu có thể chuyển đổi sang hai mô hình thay thế linh hoạt hơn:

  1. Mô hình Generalized Ordered Logit (Gologit2 trên STATA): Cho phép hệ số $\beta$ thay đổi tự do giữa các ngưỡng phân loại.
  2. Mô hình Partial Proportional Odds (PPOL): Chỉ nới lỏng giả định hồi quy song song đối với các biến độc lập cụ thể bị vi phạm, giữ nguyên tính đóng băng tham số đối với các biến thỏa mãn.

4. Kết quả nghiên cứu có thể áp dụng cho các địa phương khác ngoài Thái Bình không?

Có thể mở rộng áp dụng cho các tỉnh thuộc vùng Đồng bằng châu thổ sông Hồng có đặc điểm kinh tế - xã hội tương đồng (như Hải Dương, Nam Định, Hà Nam, Hưng Yên) - nơi có tỷ lệ đất sản xuất lúa nước cao, mật độ dân số lớn và đang trong giai đoạn chuyển đổi mạnh mẽ từ sản xuất nông nghiệp thuần túy sang phát triển các khu công nghiệp tập trung.

5. Yếu tố nào đóng vai trò then chốt nhất trong việc phục hồi thu nhập nông hộ sau khi mất đất sản xuất?

Kết quả định lượng cho thấy hai yếu tố có tác động biên mạnh nhất và có ý nghĩa thống kê cao nhất ($p < 0,01$) là:

  • Tỷ lệ thu nhập phi nông nghiệp (Phi_Nongnghiep): Giúp tăng xác suất lọt vào nhóm thu nhập cao nhất ($> 20\text{ triệu/tháng}$) thêm $3,5%$ với mỗi $1%$ tăng tỷ trọng.
  • Tiếp cận vốn tín dụng đầu tư sinh kế mới (Vayvon): Tăng khả năng nâng cao mức thu nhập lên gấp $1,84$ lần so với các hộ không có vốn hỗ trợ sản xuất kinh doanh phi nông nghiệp.

Kết luận

Khóa luận tốt nghiệp "Đánh giá tác động của thu hồi đất đến phúc lợi của các hộ gia đình tại tỉnh Thái Bình" đã giải quyết trọn vẹn mục tiêu nghiên cứu khoa học đề ra thông qua việc kết hợp cơ sở lý luận kinh tế chính trị với phương pháp kinh tế lượng vi mô thực nghiệm.

Nghiên cứu khẳng định: Việc thu hồi đất nông nghiệp mang lại động lực tăng trưởng kinh tế vĩ mô và cải thiện vượt bậc hệ thống cơ sở hạ tầng giao thông ($76,34%$ hài lòng), y tế và giáo dục ($68,82%$ hài lòng) cho tỉnh Thái Bình; tuy nhiên, nó cũng tạo ra những áp lực lớn lên sinh kế của các hộ nông dân, đặc biệt là nhóm chủ hộ lớn tuổi và những hộ phụ thuộc hoàn toàn vào nông nghiệp đơn thuần.

Việc giải ngân bồi thường một lần không giải quyết triệt để vấn đề phúc lợi lâu dài nếu thiếu vắng các giải pháp hỗ trợ đồng bộ: đào tạo nghề theo nhu cầu thực tế của doanh nghiệp, mở rộng hạn mức tín dụng ưu đãi sinh kế phi nông nghiệp và minh bạch hóa cơ chế bồi thường - tái định cư. Đây là nguồn tài liệu tham khảo có giá trị lý luận và thực tiễn sâu sắc cho các nhà hoạch định chính sách đất đai, nhà nghiên cứu kinh tế phát triển và sinh viên chuyên ngành.