Luận án tiến sĩ: Ảnh hưởng của di cư lao động đến phát triển kinh tế tỉnh Bắc Ninh

Tổng quan về luận án

Quá trình công nghiệp hóa và hiện đại hóa tại Việt Nam đã thúc đẩy sự hình thành các cực tăng trưởng kinh tế năng động tại Vùng kinh tế trọng điểm Bắc Bộ. Tỉnh Bắc Ninh, với đặc thù là địa phương có diện tích tự nhiên nhỏ nhất cả nước (822,7 km²) nhưng mật độ dân số cao thứ ba toàn quốc, đã trở thành điển hình thực nghiệm rõ nét về mô hình tăng trưởng dựa trên thu hút dòng vốn đầu tư trực tiếp nước ngoài (FDI) kết hợp thâm dụng lao động di cư. Luận án tiến sĩ chuyên ngành Quản trị nhân lực (mã số: 9 34 04 04) của nghiên cứu sinh Đặng Thanh Mai, do PGS.TS. Tô Thế Nguyên và TS. Trần Văn Đức hướng dẫn tại Học viện Nông nghiệp Việt Nam (2024), mang tính tiên phong khi giải quyết bài toán lượng hóa toàn diện tác động đa chiều của di cư lao động đến cấu trúc kinh tế địa phương.

Khoảng trống nghiên cứu (research gap) mà luận án tập trung giải quyết xuất phát từ thực trạng phần lớn các công trình trước đây chỉ dừng lại ở các thống kê mô tả đơn biến, khảo sát động cơ di cư hoặc phân tích thuần túy mặt nhân khẩu học xã hội (Nguyen Thu Phuong et al., 2008; ILO, 2015; Le Quoc Hoi et al., 2020). Thiếu vắng các nghiên cứu mô hình hóa kinh tế lượng chặt chẽ nhằm đánh giá tác động hai chiều của lao động nhập cư đến đồng thời cả các biến số nguồn lực (vốn FDI, tích lũy vốn nội tại, tái phân bổ lực lượng lao động) và các biến số kết quả phát triển kinh tế (tổng sản phẩm trên địa bàn - GRDP, chuyển dịch cơ cấu kinh tế, tỷ lệ đô thị hóa, thất nghiệp, chỉ số giá tiêu dùng CPI và tính bền vững môi trường).

Hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết nghiên cứu (Hypotheses - H) được xác lập chặt chẽ:

• RQ1: Thực trạng quy mô, động thái và cấu trúc của dòng di cư lao động nội địa đến tỉnh Bắc Ninh diễn ra như thế nào trong tiến trình 25 năm tái lập tỉnh (1997–2021)?
• RQ2: Dòng di cư lao động tác động tích cực và tiêu cực đến các nguồn lực đầu vào và kết quả phát triển kinh tế tỉnh Bắc Ninh qua những cơ chế truyền dẫn nào?
• RQ3: Những giải pháp quản trị nhân lực xã hội và điều tiết kinh tế vĩ mô nào có thể tối ưu hóa lợi ích kinh tế, đồng thời giảm thiểu ngoại ứng tiêu cực của lao động nhập cư?
• H1: Tỷ lệ lao động di cư trong lực lượng lao động (IMRA) có tác động tích cực dài hạn đến tổng sản phẩm trên địa bàn (LGRDP).
• H2: Tỷ lệ lao động di cư gia tăng tạo ra áp lực thay thế việc làm cục bộ, tác động thuận chiều đến tỷ lệ thất nghiệp (LUNEM) của địa phương trong ngắn hạn và dài hạn.
• H3: Dòng lao động nhập cư thúc đẩy nhanh tốc độ chuyển dịch cơ cấu kinh tế theo hướng tăng tỷ trọng công nghiệp – xây dựng và gia tăng tỷ lệ đô thị hóa.
• H4: Thu nhập thực tế, chi phí sinh hoạt, điều kiện nhà ở và chính sách an sinh xã hội có ý nghĩa quyết định xác suất gắn bó lâu dài của người lao động di cư.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp hàm sản xuất tân cổ điển Cobb-Douglas (1928), mô hình tăng trưởng nội sinh Solow (1956), lý thuyết kinh tế nhị nguyên Lewis (1954) và mô hình di cư nông thôn – đô thị Harris & Todaro (1970).

Đóng góp đột phá của công trình được định lượng cụ thể: Di cư lao động cung cấp gần 30% lực lượng lao động toàn tỉnh, chiếm trên 70% nhân lực trong các khu công nghiệp (KCN), đóng góp hơn 32% vào tổng giá trị sản xuất GRDP (theo phương pháp năng suất lao động) và khoảng 12% GRDP (theo phương pháp chi tiêu) năm 2021. Nghiên cứu thực hiện trên phạm vi chuỗi số liệu thứ cấp 1997–2021 và khảo sát sơ cấp quy mô 533 lao động di cư năm 2021.

Literature Review và Positioning

Tổng quan tài liệu của luận án hệ thống hóa ba dòng lý thuyết chủ lưu trong kinh tế học di cư:

                          ┌──────────────────────────────────────────────────────────┐
                          │         CÁC DÒNG LÝ THUYẾT DI CƯ VÀ PHÁT TRIỂN           │
                          └────────────────────────────┬─────────────────────────────┘
                                                       │
         ┌─────────────────────────────────────────────┼────────────────────────────────────────────┐
         │                                             │                                            │
         ▼                                             ▼                                            ▼
┌─────────────────────────────────┐       ┌─────────────────────────────────┐      ┌─────────────────────────────────┐
│     TIẾP CẬN TÂN CỔ ĐIỂN        │       │      TIẾP CẬN CẤU TRÚC VÀ       │      │     KINH TẾ MỚI VỀ DI CƯ        │
│          VÀ VĨ MÔ               │       │     THỊ TRƯỜNG KÉP (Duality)    │      │        LAO ĐỘNG (NELM)          │
├─────────────────────────────────┤       ├─────────────────────────────────┤      ├─────────────────────────────────┤
│ • Lewis (1954); Oshima (1962)   │       │ • Piore (1979); Sassen (1988)   │      │ • Stark & Bloom (1985)          │
│ • Sjaastad (1962); Lee (1966)   │       │ • Wallerstein (1974)            │      │ • Stark & Taylor (1991)         │
│ • Harris & Todaro (1970)        │       │ • Castles & Kosack (1973)       │      │ • Wouterse & Taylor (2007)      │
├─────────────────────────────────┤       ├─────────────────────────────────┤      ├─────────────────────────────────┤
│ Phân bổ lại lao động dư thừa từ │       │ Phân tầng thị trường lao động,  │      │ Đa dạng hóa sinh kế, phân tán   │
│ nông nghiệp sang công nghiệp,   │       │ tập trung di cư vào phân khúc   │      │ rủi ro và khắc phục thất bại    │
│ cân bằng chênh lệch tiền lương. │       │ thứ cấp có tiền công thấp.      │      │ thị trường cấp hộ gia đình.     │
└─────────────────────────────────┘       └─────────────────────────────────┘      └─────────────────────────────────┘

Trong dòng nghiên cứu thực nghiệm quốc tế, tồn tại hai luồng quan điểm đối lập sâu sắc về tác động của lao động nhập cư. Nhóm tác giả Feridun (2004), Morley (2006) và Boubtane et al. (2016) khẳng định nhập cư kích thích mở rộng tổng cầu, bổ sung nguồn nhân lực thiếu hụt, nâng cao hiệu quả sử dụng vốn và thúc đẩy tăng trưởng kinh tế nhanh chóng. Ngược lại, nhóm nghiên cứu của Card (2001), Dustmann & Preston (2005, 2015), Saiz (2007) và Islam (2007) chỉ ra các bằng chứng về "hiệu ứng chèn ép" (displacement effect), làm tăng tỷ lệ thất nghiệp cục bộ của lao động bản địa, gia tăng chênh lệch chi phí sinh hoạt, lạm phát giá nhà ở và gây quá tải hệ thống an sinh công cộng.

Luận án định vị vị trí học thuật bằng việc đối chiếu trực tiếp với hai trường hợp quốc tế điển hình có nhiều nét tương đồng về cấu trúc phát triển:

1. Thành phố Thượng Hải (Trung Quốc): Các nghiên cứu của Sun (2004), Lua & Xia (2016) và Chan (2021) chỉ ra lao động nhập cư nông thôn đóng góp 31% GDP của Thượng Hải, thúc đẩy tỷ lệ đô thị hóa vượt 90%, nhưng buộc chính quyền phải thực hiện cải cách hệ thống Hộ khẩu (Hukou) theo cơ chế tính điểm và triển khai quỹ nhà ở xã hội giá rẻ 5.500 CNY/m² để ổn định xã hội.
2. Vùng Thủ đô Seoul - Sudogwon (Hàn Quốc): Các công trình của Kim & Donaldson (1979), Lee & Farber (1985) và Seol (2018) chứng minh dòng di cư giai đoạn 1960–1970 phục vụ các tập đoàn công nghiệp (Chaebol) đóng góp 46% GDP quốc gia nhưng dẫn tới quá tải đô thị nghiêm trọng. Hệ quả là Hàn Quốc phải thiết lập thành phố hành chính mới Sejong (2012) với ngân sách 160 nghìn tỷ won và áp dụng "Chương trình mua nhà lần đầu" hỗ trợ vay 70% giá trị căn hộ với lãi suất ưu đãi 6%–6,5%/năm.

Tại Việt Nam, luận án so sánh với các đô thị cực lớn như Hà Nội (Lê Thúy, 2019; Đinh Văn Thông, 2010) và TP. Hồ Chí Minh (Lê Thị Hồng Rin, 2015; Nguyễn Nữ Đoàn Vy, 2018), qua đó chỉ ra Bắc Ninh mang tính chất của một "thủ phủ công nghiệp nén" – nơi sự bùng nổ diễn ra trong một không gian địa lý hẹp với tốc độ cơ cấu chuyển đổi vượt bậc nhưng khả năng chống chịu hạ tầng xã hội còn nhiều điểm nghẽn.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng khung phân tích tăng trưởng tân cổ điển và kinh tế học phát triển bằng cách làm rõ cơ chế dẫn truyền của dòng lao động nhập cư trong một nền kinh tế mở cấp địa phương. Bằng chứng thực nghiệm khẳng định di cư lao động không chỉ thuần túy là sự dịch chuyển nhân khẩu mà là một quá trình tái phân bổ không gian của các nhân tố sản xuất, tác động trực tiếp đến hàm sản xuất tổng hợp:

$$Y_t = A_t K_t^\alpha (L_{local} + \theta L_{migrant})^\beta$$

Trong đó $\theta$ phản ánh hệ số đóng góp hiệu quả vượt trội của lao động di cư nhờ lợi thế độ tuổi trẻ và cường độ lao động cao trong các dây chuyền công nghiệp hiện đại.

                    ┌────────────────────────────────────────────────────────┐
                    │      KHUNG MÔ HÌNH LÝ THUYẾT ĐÓNG GÓP CỦA LUẬN ÁN      │
                    └───────────────────────────┬────────────────────────────┘
                                                │
                 ┌──────────────────────────────┴──────────────────────────────┐
                 ▼                                                             ▼
┌───────────────────────────────────────────────┐             ┌───────────────────────────────────────────────┐
│     MỞ RỘNG MÔ HÌNH TĂNG TRƯỞNG TÂN CỔ ĐIỂN   │             │   PHÂN TẦNG VÀ NGOẠI ỨNG KINH TẾ ĐỊA PHƯƠNG   │
├───────────────────────────────────────────────┤             ├───────────────────────────────────────────────┤
│ • Tích hợp biến số tỷ lệ lao động di cư       │             │ • Kiểm chứng sự tồn tại của thị trường lao    │
│   (IMRA) vào hàm sản xuất vĩ mô.              │             │   động kép (Dual Labor Market) tại địa phương.│
│ • Xác lập tác động đa chiều tới hệ số hiệu    │             │ • Đo lường ngoại ứng tiêu cực (thất nghiệp    │
│   quả sử dụng vốn đầu tư (ICOR).              │             │   cục bộ, lạm phát giá nhà ở, quá tải ASXH).  │
└───────────────────────────────────────────────┘             └───────────────────────────────────────────────┘

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp bốn cấu phần cốt lõi, tạo nên cách tiếp cận hệ thống toàn diện:

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                      KHUNG PHÂN TÍCH TÁC ĐỘNG CỦA DI CƯ LAO ĐỘNG (DCLĐ)                          │
└────────────────────────────────────────────────┬─────────────────────────────────────────────────┘
                                                 │
        ┌────────────────────────────────────────┴────────────────────────────────────────┐
        ▼                                                                                 ▼
┌──────────────────────────────────────────────┐                  ┌──────────────────────────────────────────────┐
│       NGUỒN LỰC PHÁT TRIỂN KINH TẾ           │                  │       KẾT QUẢ PHÁT TRIỂN KINH TẾ             │
├──────────────────────────────────────────────┤                  ├──────────────────────────────────────────────┤
│ 1. Nguồn nhân lực:                           │                  │ 1. Tăng trưởng kinh tế (GRDP)                │
│    - Quy mô, độ tuổi, giới tính              │                  │ 2. Chuyển dịch cơ cấu kinh tế (CCKT)         │
│    - Trình độ chuyên môn kỹ thuật (CMKT)     │   Tác động qua   │ 3. Tiến trình đô thị hóa (ĐTH)               │
│    - Năng suất lao động (NSLĐ)               │ ───────────────> │ 4. Tỷ lệ thất nghiệp (UNEM)                  │
│ 2. Nguồn lực vốn:                            │   hệ thống chỉ   │ 5. Chỉ số giá tiêu dùng (CPI)                │
│    - Vốn đầu tư trực tiếp nước ngoài (FDI)   │   tiêu đo lường  │ 6. Phát triển bền vững:                      │
│    - Vốn đầu tư trong nước                   │                  │    - Giảm nghèo đa chiều                     │
│    - Hiệu quả sử dụng vốn (Hệ số ICOR)       │                  │    - Chi ngân sách an sinh xã hội            │
│                                              │                  │    - Môi trường và an ninh trật tự           │
└──────────────────────────────────────────────┘                  └──────────────────────────────────────────────┘

Điều kiện biên (boundary conditions) của khung phân tích được xác lập: Mô hình áp dụng đặc thù cho các nền kinh tế địa phương cấp tỉnh có tốc độ công nghiệp hóa nhanh, dòng di cư thuần mang giá trị dương lớn, quy mô dòng nhập cư nội địa chiếm tỷ trọng áp đảo so với di cư quốc tế và có sự tập trung dày đặc của các khu, cụm công nghiệp định hướng xuất khẩu.

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

Thiết kế nghiên cứu

Luận án vận dụng thế giới quan bản thể luận hiện thực và nhận thức luận thực chứng (Positivism) kết hợp hiện thực phản biện (Critical Realism). Thiết kế nghiên cứu đa phương pháp (Mixed Methods) kết hợp chặt chẽ giữa chuỗi thời gian kinh tế lượng vĩ mô (Macro Time-Series Econometrics) và phân tích vi mô cắt ngang (Micro Cross-Sectional Analysis).

                      ┌──────────────────────────────────────────────────────────┐
                      │              QUY TRÌNH THIẾT KẾ PHƯƠNG PHÁP              │
                      └────────────────────────────┬─────────────────────────────┘
                                                   │
                ┌──────────────────────────────────┴──────────────────────────────────┐
                ▼                                                                     ▼
┌───────────────────────────────────────────────┐                     ┌───────────────────────────────────────────────┐
│     PHÂN TÍCH CHUỖI THỜI GIAN VĨ MÔ           │                     │        KHẢO SÁT VI MÔ VÀ ĐỊNH TÍNH            │
│              (1997 – 2021)                    │                     │                 (N = 533)                     │
├───────────────────────────────────────────────┤                     ├───────────────────────────────────────────────┤
│ • Kiểm định tính dừng (ADF, PP Test)          │                     │ • Khảo sát 533 lao động di cư tại 8 huyện/    │
│ • Mô hình phân phối trễ tự hồi quy (ARDL)     │                     │   thị xã/thành phố tỉnh Bắc Ninh.             │
│ • Kiểm định đồng tích hợp biên (Bounds Test)  │                     │ • Mô hình hồi quy Binary Probit đánh giá      │
│ • Đánh giá mối quan hệ LIMRA, LUNEM, LGRDP    │                     │   các yếu tố quyết định sự gắn bó lâu dài.    │
└───────────────────────────────────────────────┘                     └──────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

• Chiến lược chọn mẫu: Áp dụng phương pháp chọn mẫu phân tầng kết hợp ngẫu nhiên thuận tiện. Cỡ mẫu sơ cấp đạt $N = 533$ người lao động di cư tại các khu nhà trọ, doanh nghiệp trong và ngoài KCN thuộc 8 đơn vị hành chính cấp huyện của tỉnh Bắc Ninh năm 2021.
• Thu thập dữ liệu định tính: Thực hiện 25 cuộc phỏng vấn sâu (In-depth Interviews) với cán bộ Sở LĐ-TBXH, Ban Quản lý các KCN, Cục Thống kê, lãnh đạo UBND cấp xã/huyện, đại diện doanh nghiệp FDI và các chủ cơ sở cho thuê nhà trọ.
• Tam giác hóa dữ liệu (Triangulation): Kết hợp kiểm chứng chéo giữa số liệu Tổng điều tra Dân số và Nhà ở (1999, 2009, 2019), Niên giám Thống kê tỉnh Bắc Ninh (1997–2021), số liệu hành chính của cơ quan quản lý và dữ liệu điều tra vi mô.
• Độ tin cậy và giá trị: Thang đo khảo sát được kiểm định độ tin cậy nhất quán nội tại (Cronbach's $\alpha > 0,78$) và giá trị hội tụ thông qua phân tích nhân tố khám phá (EFA).

Data và phân tích

Mô hình toán kinh tế lượng vĩ mô sử dụng kỹ thuật ARDL (Autoregressive Distributed Lag) với phần mềm Stata và EViews để ước lượng phương trình đồng tích hợp dài hạn:

$$\text{LGRDP}t = \alpha_0 + \sum{i=1}^{p} \beta_i \text{LGRDP}{t-i} + \sum{j=0}^{q_1} \gamma_j \text{LIMRA}{t-j} + \sum{k=0}^{q_2} \delta_k \text{LUNEM}_{t-k} + \varepsilon_t$$

Trong đó:

• $\text{LGRDP}$: Logarit tự nhiên của Tổng sản phẩm trên địa bàn theo giá so sánh;
• $\text{LIMRA}$: Logarit của tỷ lệ lao động di cư trong lực lượng lao động;
• $\text{LUNEM}$: Logarit của tỷ lệ thất nghiệp toàn tỉnh.

Mô hình vi mô Probit ước lượng xác suất gắn bó lâu dài ($Y = 1$ nếu dự định sống trên 5 năm/định cư):

$$P(Y_i = 1 | X_i) = \Phi(\beta_0 + \beta_1 \text{ThuNhap}_i + \beta_2 \text{HocVan}_i + \beta_3 \text{NhaO}_i + \beta_4 \text{BaoHiem}_i + \beta_5 \text{Tuoi}_i + u_i)$$

Các kiểm định chẩn đoán khuyết tật mô hình (Diagnostic tests) bao gồm: Kiểm định tự tương quan Breusch-Godfrey LM, kiểm định phương sai sai số thay đổi White/Breusch-Pagan, kiểm định dạng hàm Ramsey RESET và kiểm định độ ổn định tham số qua đồ thị CUSUM và CUSUMSQ đều cho kết quả vững ở mức ý nghĩa $p < 0,05$.

Phát hiện đột phá và implications

Những phát hiện then chốt

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                             5 PHÁT HIỆN THỰC NGHIỆM ĐỘT PHÁ CỦA LUẬN ÁN                          │
└────────────────────────────────┬─────────────────────────────────────────────────────────────────┘
                                 │
  ┌──────────────────────────────┼──────────────────────────────┬──────────────────────────────┐
  │                              │                              │                              │
  ▼                              ▼                              ▼                              ▼
┌─────────────────────────┐    ┌─────────────────────────┐    ┌─────────────────────────┐    ┌─────────────────────────┐
│  ĐÓNG GÓP QUY MÔ GRDP   │    │ CHUYỂN DỊCH CƠ CẤU VÀ   │    │ HIỆU ỨNG THAY THẾ       │    │ NGOẠI ỨNG TIÊU CỰC VỀ   │
│       VƯỢT TRỘI         │    │       ĐÔ THỊ HÓA        │    │  VIỆC LÀM (LUNEM)       │    │ GIÁ CẢ VÀ MÔI TRƯỜNG    │
├─────────────────────────┤    ├─────────────────────────┤    ├─────────────────────────┤    ├─────────────────────────┤
│ Đóng góp >32% GRDP toàn │    │ Tỷ trọng CN-XD đạt      │    │ Hệ số dài hạn của       │    │ CPI tăng cao do cầu     │
│ tỉnh theo phương pháp   │    │ 77,3% năm 2021; Đô thị  │    │ LIMRA lên LUNEM mang dấu│    │ nhà trọ; 16% dân số là  │
│ sản xuất; 12% theo chi  │    │ hóa tăng từ 9,5% (2000) │    │ dương có ý nghĩa thống  │    │ LĐDC tạo áp lực lớn lên │
│ tiêu (năm 2021).        │    │ lên 42,6% (2020/2021).  │    │ kê ở mức p < 0,01.      │    │ xử lý rác thải sinh hoạt│
└─────────────────────────┘    └─────────────────────────┘    └─────────────────────────┘    └─────────────────────────┘

1. Đóng góp định lượng vượt bậc vào tăng trưởng kinh tế: Theo tính toán của tác giả:

   "Tính đến năm 2021, DCLĐ đóng góp gần 30% lực lượng lao động của toàn tỉnh... đóng góp hơn 32% cho GRDP của toàn tỉnh tính theo phương thức tính giá trị sản xuất thông qua năng suất lao động và đóng góp khoảng 12% cho GRDP toàn tỉnh theo phương thức chi tiêu." Tốc độ tăng trưởng GRDP bình quân giai đoạn 1997–2021 đạt 13,89%/năm, đưa quy mô nền kinh tế từ 2,1 nghìn tỷ đồng (chiếm 0,64% GDP cả nước năm 1997) lên 38,7 nghìn tỷ đồng (chiếm 1,8% GDP cả nước năm 2021).

2. Tái cấu trúc triệt để cơ cấu ngành và đẩy nhanh đô thị hóa nén: Sự gia tăng của 236.000 lao động di cư vào năm 2021 (với 87,5% đến vì mục tiêu việc làm và 56,1% là nữ) đã đưa tỷ trọng ngành công nghiệp – xây dựng của Bắc Ninh chạm mốc 77,3% vào năm 2021. Tỷ lệ đô thị hóa tăng vọt từ 9,5% năm 2000 lên 42,6% năm 2020/2021.

3. Bằng chứng thực nghiệm về hiệu ứng chèn ép việc làm cục bộ: Kết quả hồi quy mô hình ARDL dài hạn chỉ ra mối quan hệ đồng tích hợp mang dấu dương giữa logarit tỷ lệ lao động di cư ($\text{LIMRA}$) và logarit tỷ lệ thất nghiệp ($\text{LUNEM}$) với hệ số hồi quy có ý nghĩa thống kê ($p < 0,01$). Điều này giải thích hiện tượng lao động nhập cư chấp nhận mức lương cạnh tranh và cường độ làm việc cao trong các KCN đã gián tiếp đẩy một bộ phận lao động bản địa lớn tuổi, tay nghề thấp vào tình trạng thất nghiệp hoặc phải chuyển dịch sang khu vực phi chính thức.

4. Tác động làm gia tăng mặt bằng giá tiêu dùng và chi phí sinh hoạt: Do sự tập trung dân số đột biến (chiếm gần 16% tổng dân số tỉnh), nhu cầu về hàng hóa thiết yếu và đặc biệt là dịch vụ nhà ở cho thuê tăng vọt trong khi nguồn cung co giãn kém. Chỉ số CPI của Bắc Ninh trong nhiều giai đoạn duy trì ở mức cao hơn bình quân cả nước, hình thành mặt bằng giá cả sinh hoạt đắt đỏ đối với nhóm người có thu nhập thấp.

5. Phân hóa xã hội, xung đột văn hóa và rủi ro môi trường: Gần 50% lao động di cư là đồng bào dân tộc thiểu số đến từ vùng Trung du và miền núi phía Bắc với thời gian lưu trú bình quân dưới 5 năm. Tình trạng thiếu thiết chế văn hóa, điều kiện nhà trọ chật chội và việc phát sinh lượng lớn rác thải sinh hoạt cục bộ tại các xã tập trung KCN (như Yên Phong, Quế Võ) đặt ra thách thức lớn đối với chỉ số phát triển bền vững.

Implications đa chiều

• Về mặt lý luận: Bổ sung bằng chứng thực nghiệm vững chắc cho lý thuyết thị trường lao động phân tầng của Piore (1979) tại các nền kinh tế chuyển đổi đón nhận dòng vốn FDI thế hệ mới.
• Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp mô hình kinh tế lượng chuỗi thời gian ARDL vĩ mô với mô hình vi mô Probit trong phân tích di cư cấp vùng.
• Về mặt chính sách và thực tiễn:
  • Cần tích hợp quy hoạch phát triển dòng di cư vào Quy hoạch tổng thể phát triển kinh tế - xã hội tỉnh Bắc Ninh.
  • Đổi mới chính sách quản lý cư trú trên nền tảng số, dỡ bỏ triệt để rào cản hành chính trong tiếp cận dịch vụ công.
  • Thiết lập chương trình mục tiêu xây dựng nhà ở xã hội cho thuê với chuẩn tối thiểu từ 8–15 m²/người, đa dạng hóa các gói hỗ trợ tín dụng ưu đãi an cư.
  • Mở rộng mạng lưới an sinh xã hội, tăng tỷ lệ bao phủ bảo hiểm xã hội bắt buộc và bảo hiểm y tế cho lao động ngoại tỉnh.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

1. Rào cản thời điểm thu thập dữ liệu vi mô: Cuộc khảo sát sơ cấp $N=533$ tiến hành trong năm 2021 – thời điểm dịch bệnh Covid-19 bùng phát phức tạp tại các KCN Bắc Ninh, có thể tạo ra những nhiễu động ngắn hạn đối với tâm lý và đánh giá mức độ gắn bó của người lao động.
2. Giới hạn biến số chuyên gia quốc tế: Luận án tập trung vào dòng di cư nội địa do tỷ trọng lao động nước ngoài nhỏ (~1%), bỏ qua tác động lan tỏa công nghệ từ nhóm chuyên gia kỹ thuật cao người Hàn Quốc, Nhật Bản, Trung Quốc.
3. Thiếu vắng dữ liệu bảng cấp doanh nghiệp (Firm-level Panel Data): Chưa bóc tách chi tiết được hàm năng suất nhân tố tổng hợp (TFP) của từng nhóm doanh nghiệp FDI so với doanh nghiệp dân doanh trong mối liên hệ với tỷ lệ thâm dụng lao động nhập cư.

Chương trình nghiên cứu 10 năm tiếp theo (2025–2035):

• Hướng 1: Ứng dụng mô hình dữ liệu bảng không gian (Spatial Panel Econometrics) để phân tích tác động tràn (spillover effects) của lao động di cư giữa Bắc Ninh và các tỉnh lân cận trong Vùng Thủ đô (Hà Nội, Hưng Yên, Hải Dương, Bắc Giang).
• Hướng 2: Đánh giá tác động của làn sóng tự động hóa, chuyển đổi số và trí tuệ nhân tạo (AI) trong các nhà máy bán dẫn/điện tử đến nguy cơ dịch chuyển ngược của dòng lao động di cư phổ thông.
• Hướng 3: Nghiên cứu chuyên sâu về di cư thế hệ thứ hai (Second-generation immigrants) và sự hòa nhập kinh tế - xã hội bền vững tại các đô thị công nghiệp vệ tinh.

Tác động và ảnh hưởng

                                 ┌─────────────────────────────────────────────────────────┐
                                 │              HỆ THỐNG TÁC ĐỘNG ĐA LĨNH VỰC              │
                                 └────────────────────────────┬────────────────────────────┘
                                                              │
             ┌────────────────────────────────────────────────┼────────────────────────────────────────────────┐
             │                                                │                                                │
             ▼                                                ▼                                                ▼
┌─────────────────────────┐                      ┌─────────────────────────┐                      ┌─────────────────────────┐
│     HỌC THUẬT VÀ VIỆN   │                      │  HOẠCH ĐỊNH CHÍNH SÁCH  │                      │ DOANH NGHIỆP VÀ XÃ HỘI  │
│        NGHIÊN CỨU       │                      │     VÀ CHÍNH QUYỀN      │                      │     KHU CÔNG NGHIỆP     │
├─────────────────────────┤                      ├─────────────────────────┤                      ├─────────────────────────┤
│ Cung cấp khung phương   │                      │ Luận cứ khoa học chuẩn  │                      │ Tái cơ cấu chiến lược   │
│ pháp luận định lượng    │                      │ xác cho Tỉnh ủy, HĐND,  │                      │ giữ chân nhân sự (Talent│
│ chuẩn hóa; ước tính     │                      │ UBND tỉnh Bắc Ninh ban  │                      │ Retention), xây dựng ký │
│ trích dẫn học thuật cao │                      │ hành Nghị quyết chuyên  │                      │ túc xá đạt chuẩn, tối ưu│
│ trong các công bố quốc  │                      │ đề về phát triển nhà ở  │                      │ hóa chuỗi cung ứng lao  │
│ tế và Scopus/WoS.       │                      │ xã hội và an sinh LĐDC. │                      │ động bền vững.          │
└─────────────────────────┘                      └─────────────────────────┘                      └─────────────────────────┘

Công trình có tiềm năng tác động chính sách mạnh mẽ, cung cấp luận cứ trực tiếp giúp tỉnh Bắc Ninh hoàn thiện các đề án quy hoạch chuyển đổi từ "thủ phủ công nghiệp gia công" sang "trung tâm đổi mới sáng tạo và công nghiệp công nghệ cao", giải quyết bài toán phát triển hài hòa giữa tích tụ kinh tế và tiến bộ xã hội.

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

• Nghiên cứu sinh và Giới học thuật: Thụ hưởng một khung phân tích tích hợp có thể tái lập (replicable), bộ biến số định lượng tường minh (LIMRA, LGRDP, LUNEM) và nguồn trích dẫn lý thuyết phong phú phục vụ các nghiên cứu kinh tế lao động.
• Lãnh đạo cơ quan Quản lý Nhà nước (Bộ LĐ-TBXH, Ban Kinh tế Trung ương, UBND tỉnh): Sở hữu bằng chứng thực nghiệm rõ ràng để thiết kế các chính sách phân bổ ngân sách y tế, giáo dục, quy hoạch đất đai phát triển nhà ở công nhân dựa trên quy mô dân số thực tế thay vì dân số thường trú danh nghĩa.
• Khối Doanh nghiệp FDI và Ban Quản lý KCN: Có cơ sở dữ liệu định lượng về đặc tính nhân khẩu học (độ tuổi, học vấn, giới tính, kỳ vọng thu nhập) để xây dựng chính sách tiền lương, phúc lợi và đào tạo nghề phù hợp, giảm thiểu tỷ lệ nhảy việc và chi phí tuyển dụng lại.
• Cộng đồng Người lao động di cư: Được gián tiếp thụ hưởng các chính sách cải thiện về môi trường cư trú, bảo trợ pháp lý và công bằng trong tiếp cận dịch vụ an sinh xã hội tại nơi nhập cư.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Luận án đã mở rộng mô hình phân tầng thị trường lao động của Piore (1979) và hàm sản xuất Cobb-Douglas bằng cách tích hợp trực tiếp tỷ lệ thâm nhập lao động di cư ($\text{IMRA}$) như một biến số nội sinh quyết định cả quy mô giá trị gia tăng GRDP và áp lực cơ cấu lên tỷ lệ thất nghiệp địa phương ($\text{UNEM}$). Luận án chứng minh rằng trong một cực tăng trưởng FDI công nghiệp hóa nhanh, dòng di cư lao động tạo ra hiệu ứng kép: Vừa là động lực kích thích tăng trưởng quy mô kinh tế thần kỳ (>32% GRDP), vừa là tác nhân tái cấu trúc áp lực việc làm đối với lực lượng lao động tại chỗ.

2. Điểm mới về mặt phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu của Nguyen Viet Cuong et al. (2011) hay Trương Văn Tuấn (2011) vốn thuần túy dùng thống kê mô tả hoặc OLS cắt ngang, luận án đã thiết lập quy trình kiểm định kinh tế lượng chuỗi thời gian 25 năm vững chắc thông qua mô hình ARDL Bounds Testing. Quy trình này xử lý triệt để hiện tượng chuỗi số liệu phi dừng $I(0)/I(1)$, loại bỏ hiện tượng hồi quy giả mạo (spurious regression) và phân tách rành mạch tác động co giãn của lao động di cư trong ngắn hạn và dài hạn, kết hợp đối soát với mô hình Probit vi mô trên mẫu khảo sát $N=533$.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ và mang tính nghịch lý (counter-intuitive) là: Mặc dù dòng lao động di cư đóng góp vượt bậc vào mở rộng quy mô GRDP, nhưng sự gia tăng tỷ lệ lao động nhập cư lại có mối tương quan thuận chiều có ý nghĩa thống kê với tỷ lệ thất nghiệp của địa phương trong dài hạn. Dữ liệu chứng minh các tập đoàn FDI ưu tiên tuyển dụng lao động di cư trẻ tuổi (với gần 50% là người dân tộc thiểu số có sức khỏe dẻo dai và tính tuân thủ kỷ luật chuyền cao), dẫn đến việc đào thải và giảm cơ hội việc làm của bộ phận lao động bản địa trung niên, tay nghề thấp tại nông thôn Bắc Ninh.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án mô tả chi tiết nguồn trích xuất dữ liệu chuỗi thời gian (GSO, Cục Thống kê Bắc Ninh 1997–2021), công thức xác định hệ thống chỉ tiêu ($\text{IMRA}$, $\text{ICOR}$, $\text{GRDP}$ so sánh), cấu trúc bảng hỏi khảo sát $N=533$, phương pháp mã hóa biến số hồi quy Probit và toàn bộ quy trình kiểm định chẩn đoán mô hình ARDL trên phần mềm kinh tế lượng tiêu chuẩn.

5. Chương trình nghiên cứu 10 năm được xác lập như thế nào?

Luận án vạch ra lộ trình nghiên cứu giai đoạn 2025–2035 tập trung vào ba trọng tâm: (1) Mô hình hóa tác động lan tỏa không gian (Spatial Spillovers) trong toàn bộ Vùng Thủ đô; (2) Đo lường tác động của xu hướng tự động hóa/AI trong các nhà máy công nghệ cao đến nguy cơ đẩy ngược dòng lao động di cư; (3) Nghiên cứu các chính sách tích hợp xã hội và dịch chuyển giai tầng của thế hệ con em người lao động di cư định cư dài hạn.

Kết luận

1. Hệ thống hóa và làm sáng tỏ khung lý luận về vai trò, cơ chế tác động hai chiều của di cư lao động nội địa đến phát triển kinh tế tại các cực tăng trưởng công nghiệp hóa địa phương.
2. Cung cấp bằng chứng thực nghiệm lượng hóa chính xác: Lao động di cư đóng góp gần 30% lực lượng lao động, cung cấp trên 70% nhân lực trong các KCN và đóng góp hơn 32% tổng sản phẩm trên địa bàn (GRDP) tỉnh Bắc Ninh năm 2021.
3. Xác nhận vai trò đòn bẩy then chốt của dòng lao động nhập cư trong việc đẩy nhanh tốc độ chuyển dịch cơ cấu kinh tế sang công nghiệp – xây dựng (đạt 77,3%) và thúc đẩy tỷ lệ đô thị hóa tăng từ 9,5% lên 42,6%.
4. Chỉ ra những mặt trái và ngoại ứng kinh tế - xã hội: Áp lực thay thế việc làm cục bộ làm gia tăng tỷ lệ thất nghiệp, đẩy cao chỉ số giá tiêu dùng CPI, tạo sức ép gay gắt lên quỹ nhà ở, dịch vụ an sinh và môi trường sinh thái.
5. Kiểm định vững chắc các mối quan hệ cân bằng dài hạn thông qua mô hình kinh tế lượng ARDL và xác định các nhân tố chi phối quyết định gắn bó lâu dài của người lao động qua mô hình Probit.
6. Đề xuất hệ thống 6 nhóm giải pháp chính sách mang tính khả thi cao nhằm hoàn thiện quy hoạch kinh tế, đồng bộ quản lý cư trú số, phát triển thị trường lao động và nâng cao tính bền vững cho mô hình tăng trưởng của tỉnh Bắc Ninh và các địa phương tương đồng.