Giới thiệu dự án

Bối cảnh và thực trạng ngành

Trong tiến trình công nghiệp hóa và hiện đại hóa, ngành công nghiệp khai thác khoáng sản và năng lượng nhiệt điện giữ vai trò then chốt trong việc bảo đảm an ninh năng lượng quốc gia. Theo số liệu thống kê từ Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (VINACOMIN), hiệu quả khai thác và vận hành nhà máy phụ thuộc trực tiếp vào chất lượng và sự ổn định của lực lượng lao động. Tại tỉnh Quảng Nam, quy hoạch phát triển nhân lực giai đoạn 2016–2020 đặt mục tiêu đào tạo 239.000 lao động; tuy nhiên, khu vực này đối mặt với làn sóng dịch chuyển nhân sự chất lượng cao về các đô thị lớn như Đà Nẵng và TP. Hồ Chí Minh.

Công ty Cổ phần Than Điện Nông Sơn - VINACOMIN (địa bàn xã Quế Trung, huyện Nông Sơn, tỉnh Quảng Nam) hoạt động trong lĩnh vực khai thác, sàng tuyển than và vận hành Nhà máy Nhiệt điện Nông Sơn. Với quy mô tổng tài sản đạt 567,270 tỷ đồng (năm 2014) và 495 cán bộ công nhân viên (trong đó lao động trực tiếp chiếm 87,88%, lao động nam chiếm 94,95%), doanh nghiệp phải đối mặt với áp lực duy trì năng suất và giữ chân người lao động trong điều kiện làm việc nặng nhọc, độc hại.

                  CƠ CẤU LAO ĐỘNG VÀ VỐN TẠI VINACOMIN NÔNG SƠN

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

Mặc dù doanh thu năm 2014 đạt 751,7 tỷ đồng và lợi nhuận sau thuế tăng trưởng 87,98% so với năm 2013, công tác quản trị nguồn nhân lực tại CTCP Than Điện Nông Sơn còn tồn tại nhiều điểm nghẽn:

  1. Chính sách phân phối thu nhập chưa tối ưu: Chế độ lương - thưởng chưa phản ánh chính xác mức độ hao phí lao động và mức độ rủi ro trong môi trường hầm mỏ, nhiệt điện.
  2. Áp lực môi trường làm việc và an toàn lao động: Tính chất công việc nặng nhọc, nguy cơ bệnh nghề nghiệp cao, điều kiện cơ sở vật chất vùng sâu gây tâm lý bất an.
  3. Cơ chế đào tạo và thăng tiến thiếu tính minh bạch: 78,99% lao động là phổ thông nhưng thiếu lộ trình nâng bậc thợ rõ ràng, dẫn đến tình trạng suy giảm động lực cống hiến dài hạn.

Mục tiêu dự án

  1. Hệ thống hóa cơ sở lý luận về tạo động lực lao động thông qua tích hợp các học thuyết quản trị nhân sự kinh điển (Maslow, Herzberg, Vroom, Adams, Locke, Skinner).
  2. Đo lường và đánh giá thực trạng mức độ thỏa mãn của nhân viên đối với 7 nhóm nhân tố cấu thành động lực làm việc thông qua phương pháp định lượng thực nghiệm ($N = 135$).
  3. Xác định mô hình hồi quy tuyến tính bội phản ánh mức độ tác động của từng nhân tố đến động lực làm việc và lòng trung thành của nhân viên.
  4. Xây dựng hệ thống giải pháp khả thi, đồng bộ nhằm tối ưu hóa năng suất lao động và giữ chân nhân tài tại doanh nghiệp.

Phương pháp tiếp cận và kết quả kỳ vọng

Đề tài áp dụng quy trình nghiên cứu kết hợp: Nghiên cứu định tính (phỏng vấn sâu $n=10$) để hiệu chỉnh thang đo, kết hợp nghiên cứu định lượng ($N=135$) xử lý bằng phần mềm SPSS 20.0. Kết quả kỳ vọng cung cấp phương trình hồi quy chuẩn hóa làm căn cứ định lượng giúp ban lãnh đạo tái cấu trúc quỹ lương, cải thiện điều kiện an toàn, nâng cao tỷ lệ giữ chân lao động kỹ thuật lên trên 95%.

Phạm vi và giới hạn

  • Không gian: Toàn bộ các phòng ban, phân xưởng khai thác, vận tải, sàng tuyển, cơ điện và Ban quản lý dự án NMNĐ Nông Sơn tại Quảng Nam.
  • Thời gian: Dữ liệu thứ cấp giai đoạn 2012–2014; dữ liệu sơ cấp khảo sát từ tháng 02/2015 đến tháng 05/2015.
  • Giới hạn: Khảo sát tập trung vào khối lao động chính thức; không bao gồm lao động hợp đồng thời vụ dưới 3 tháng.

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

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

Tiêu chí Quản trị nhân sự truyền thống Mô hình định lượng tạo động lực
Cơ chế đãi ngộ Cào bằng theo thâm niên, thiếu minh bạch Đãi ngộ theo hiệu suất (3P: Position, Person, Performance)
Đo lường sự hài lòng Định tính, cảm tính của nhà quản lý Thang đo Likert 5 mức độ, kiểm định Cronbach's Alpha & EFA
Phân tích nguyên nhân Giải quyết sự vụ khi nhân viên nghỉ việc Dự báo trước qua các hệ số hồi quy đa biến $\beta$
Tính bền vững Thấp, dễ phát sinh xung đột nội bộ Cao, gắn kết lợi ích người lao động với mục tiêu doanh nghiệp

Phân loại nhu cầu nhân viên theo mô hình MoSCoW

                  MA TRẬN YÊU CẦU QUẢN TRỊ THEO NGUYÊN TẮC MOSCOW

Thiết kế hệ thống nghiên cứu

Mô hình nghiên cứu lý thuyết bao gồm 6 biến độc lập ($X_1$ đến $X_6$) và 1 biến phụ thuộc ($Y$):

$$\text{Động lực làm việc } (Y) = f(X_1, X_2, X_3, X_4, X_5, X_6)$$

Trong đó:

  • $X_1$: Lương thưởng và phúc lợi (Compensation & Benefits)
  • $X_2$: Môi trường làm việc (Work Environment)
  • $X_3$: An toàn công việc (Job Safety & Security)
  • $X_4$: Mối quan hệ giữa đồng nghiệp và cấp trên (Workplace Relationships)
  • $X_5$: Chính sách đào tạo và thăng tiến (Training & Promotion)
  • $X_6$: Sự hấp dẫn của bản thân công việc (Job Characteristics)
  • $Y$: Động lực làm việc và sự trung thành với tổ chức (Employee Motivation & Loyalty)
                            MÔ HÌNH NGHIÊN CỨU EFA VÀ HỒI QUY

Công cụ và ngăn xếp công nghệ phân tích (Technology Stack)

  • Hệ điều hành: Windows 7 / Linux Ubuntu 14.04 LTS
  • Phần mềm xử lý thống kê: IBM SPSS Statistics version 20.0
  • Phần mềm bảng tính & xử lý dữ liệu thô: Microsoft Excel 2013
  • Thang đo: Likert 5 điểm (1: Rất không đồng ý; 2: Không đồng ý; 3: Trung lập; 4: Đồng ý; 5: Rất đồng ý)

Phương pháp luận nghiên cứu (Methodology)

  1. Thiết kế chọn mẫu: Phương pháp chọn mẫu ngẫu nhiên phân tầng (Stratified Random Sampling). Cỡ mẫu tối thiểu tính theo quy tắc Hair et al. ($N \ge 5 \times k$, với $k=27$ biến quan sát $\Rightarrow N_{min} = 135$).
  2. Quy trình kiểm định chất lượng thang đo:
    • Bước 1: Đánh giá độ tin cậy thang đo bằng hệ số Cronbach's Alpha (điều kiện: $\alpha \ge 0,6$; hệ số tương quan biến - tổng $\ge 0,3$).
    • Bước 2: Phân tích nhân tố khám phá (EFA) với phép trích Principal Axis Factoring / Principal Components và phép xoay Varimax (điều kiện: $0,5 \le \text{KMO} \le 1$; kiểm định Bartlett có Sig. $< 0,05$; Eigenvalue $> 1$; tổng phương sai trích $> 50%$; hệ số tải Factor Loading $> 0,5$).
    • Bước 3: Phân tích tương quan Pearson và ước lượng mô hình hồi quy OLS (Ordinary Least Squares).

Implementation và kết quả

Quy trình phân tích dữ liệu thực nghiệm

Dưới đây là đoạn mã thực thi phân tích hồi quy và kiểm định độ tin cậy tương đương quy trình xử lý trên SPSS Syntax / Python statsmodels:

import numpy as np
import pandas as pd
import statsmodels.api as sm
from factor_analyzer import FactorAnalyzer
from statsmodels.stats.outliers_influence import variance_inflation_factor

# 1. Nạp tập dữ liệu sơ cấp N = 135 tại VINACOMIN Nông Sơn
df = pd.read_csv("vinacomin_motivation_survey.csv")

# 2. Kiểm định độ tin cậy thang đo Cronbach's Alpha
def cronbach_alpha(df_subset):
    item_vars = df_subset.var(axis=0, ddof=1)
    total_var = df_subset.sum(axis=1).var(ddof=1)
    k = df_subset.shape[1]
    return (k / (k - 1)) * (1 - item_vars.sum() / total_var)

# 3. Phân tích hồi quy tuyến tính bội OLS
X = df[['X1_LuongThuong', 'X2_MoiTruong', 'X3_AnToan', 'X4_QuanHe', 'X5_DaoTao', 'X6_HapDan']]
y = df['Y_DongLuc']
X_const = sm.add_constant(X)

ols_model = sm.OLS(y, X_const).fit()
print(ols_model.summary())

# 4. Loại biến không có ý nghĩa thống kê (p-value > 0.05: X4) và tái ước lượng
X_optimized = X.drop(columns=['X4_QuanHe'])
X_opt_const = sm.add_constant(X_optimized)
final_model = sm.OLS(y, X_opt_const).fit()
print(final_model.summary())

Kiểm định và xác thực mô hình

1. Kiểm định độ tin cậy thang đo Cronbach's Alpha

Kết quả kiểm định lần 1 và lần 2 cho thấy toàn bộ các biến quan sát thuộc 6 nhóm nhân tố độc lập và 1 nhóm biến phụ thuộc đều đạt yêu cầu ($\alpha > 0,7$, hệ số tương quan biến - tổng $> 0,35$).

Thang đo nhân tố Số biến quan sát Cronbach's Alpha Đánh giá
Lương thưởng và phúc lợi ($X_1$) 4 0.812 Đạt độ tin cậy cao
Môi trường làm việc ($X_2$) 4 0.764 Đạt độ tin cậy tốt
An toàn công việc ($X_3$) 3 0.795 Đạt độ tin cậy tốt
Mối quan hệ đồng nghiệp & cấp trên ($X_4$) 4 0.743 Đạt độ tin cậy tốt
Chính sách đào tạo và thăng tiến ($X_5$) 4 0.826 Đạt độ tin cậy cao
Sự hấp dẫn của công việc ($X_6$) 3 0.781 Đạt độ tin cậy tốt
Động lực & Lòng trung thành ($Y$) 3 0.845 Đạt độ tin cậy cao

2. Phân tích nhân tố khám phá (EFA)

  • Hệ số $\text{KMO} = 0,782$ ($> 0,5$), thỏa mãn điều kiện phân tích nhân tố.
  • Kiểm định Bartlett's Test of Sphericity: $\chi^2 = 1428,65$ với $\text{Sig.} = 0,000 < 0,001$ (các biến quan sát có tương quan tuyến tính có ý nghĩa).
  • Tổng phương sai trích (Cumulative Variance Explained) đạt $63,48% > 50%$, giá trị Eigenvalue của nhân tố thứ 6 đạt $1,142 > 1,0$.
                              BẢNG KẾT QUẢ PHÂN TÍCH EFA

3. Phân tích hồi quy tuyến tính bội OLS

Khi đưa toàn bộ 6 biến vào mô hình theo phương pháp Enter, biến $X_4$ (Mối quan hệ đồng nghiệp và cấp trên) có mức ý nghĩa $\text{Sig.} = 0,312 > 0,05$, không có ý nghĩa thống kê đối với biến phụ thuộc trong tập mẫu. Sau khi loại biến $X_4$, kết quả ước lượng mô hình tối ưu đạt được:

Biến độc lập Hệ số chưa chuẩn hóa ($B$) Sai số chuẩn (Std. Error) Hệ số chuẩn hóa ($\beta$) Giá trị $t$ Mức ý nghĩa ($\text{Sig.}$) Thống kê đa cộng tuyến (VIF)
Hằng số (Constant) 0.412 0.185 - 2.227 0.028 -
$X_1$ (Lương & Phúc lợi) 0.328 0.048 0.342 6.833 0.000 1.245
$X_3$ (An toàn công việc) 0.245 0.045 0.258 5.444 0.000 1.189
$X_6$ (Bản chất công việc) 0.198 0.042 0.205 4.714 0.000 1.312
$X_5$ (Đào tạo & Thăng tiến) 0.174 0.041 0.182 4.243 0.000 1.278
$X_2$ (Môi trường làm việc) 0.136 0.039 0.141 3.487 0.001 1.156
  • Hệ số xác định $R^2 = 0,648$; $R^2$ hiệu chỉnh $= 0,634$ (Mô hình giải thích được $63,4%$ sự biến thiên của động lực làm việc).
  • Kiểm định ANOVA: $F(5, 129) = 47,42$ với $\text{Sig.} = 0,000$ (Mô hình hồi quy hoàn toàn phù hợp).
  • Hệ số phóng đại phương sai $\text{VIF} < 2,0$ khẳng định không xảy ra hiện tượng đa cộng tuyến.

Phương trình hồi quy chuẩn hóa

$$Y = 0,342 \cdot X_1 + 0,258 \cdot X_3 + 0,205 \cdot X_6 + 0,182 \cdot X_5 + 0,141 \cdot X_2$$

                           MỨC ĐỘ TÁC ĐỘNG CỦA CÁC NHÂN TỐ (BETA)
  X1: Lương & Phúc lợi      ████████████████████████████████████ 0.342
  X3: An toàn công việc     ███████████████████████████ 0.258
  X6: Bản chất công việc    █████████████████████ 0.205
  X5: Đào tạo & Thăng tiến  ███████████████████ 0.182
  X2: Môi trường làm việc   ███████████████ 0.141

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

Điểm mới trong nghiên cứu

  1. Bản địa hóa mô hình nghiên cứu cho doanh nghiệp khai khoáng đặc thù: Khác với các nghiên cứu tổng quát tại khối dịch vụ hoặc sản xuất nhẹ, đề tài bổ sung và chuẩn hóa thang đo "An toàn công việc" ($X_3$) với hệ số tác động xếp thứ 2 ($\beta = 0,258$), phản ánh chính xác tâm lý người lao động ngành than - nhiệt điện.
  2. Khám phá thực nghiệm về biến loại trừ: Chứng minh biến $X_4$ (Quan hệ đồng nghiệp và cấp trên) không có tác động trực tiếp mang tính quyết định đến động lực làm việc trong môi trường công nghiệp nặng VINACOMIN, giúp doanh nghiệp tập trung nguồn lực tài chính vào các giải pháp trực tiếp hơn.

So sánh với các công trình nghiên cứu trước

Tiêu chí Đề tài Nhà máy Nông sản Sơn Long (2013) Đề tài Ngân hàng ACB - CN Huế (Nguyễn Khắc Hoàn, 2010) Đề tài Tổng Công ty Xi măng VN (Vũ Thị Uyên, 2010) Nghiên cứu tại CTCP Than Điện Nông Sơn (2015)
Đối tượng Lao động chế biến nông sản Nhân viên ngân hàng (Dịch vụ) Toàn thể nhân viên ngành xi măng Lao động than - nhiệt điện ($87,88%$ trực tiếp)
Cỡ mẫu ($N$) Mẫu chọn điểm nhỏ $N = 32$ Khảo sát diện rộng $N = 135$ (Phân tầng ngẫu nhiên)
Nhân tố then chốt Tiền lương, Máy móc thiết bị Môi trường, Hứng thú công việc Đánh giá thực hiện công việc Lương thưởng ($\beta=0,342$), An toàn ($\beta=0,258$)
Hạn chế Thiếu phân tích định lượng EFA sâu Cỡ mẫu quá nhỏ ($N=32$) Không phân tách đặc thù phân xưởng Kiểm định đầy đủ EFA, Cronbach's Alpha, OLS

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

Hệ thống giải pháp can thiệp chiến lược

                    HỆ THỐNG GIẢI PHÁP TẠO ĐỘNG LỰC TẠI NÔNG SƠN

Lộ trình triển khai thực tế (Roadmap)

                            LỘ TRÌNH TRIỂN KHAI 4 GIAI ĐOẠN

Phân tích hiệu quả kinh tế và lợi ích đầu tư (Cost-Benefit Analysis)

  • Chi phí dự kiến: Trích $3,5 - 5%$ lợi nhuận sau thuế (khoảng 3,3 - 4,7 tỷ đồng từ nguồn LNST năm 2014) vào quỹ phát triển nhân sự và nâng cấp an toàn mỏ.
  • Lợi ích định lượng:
    • Giảm tỷ lệ luân chuyển lao động kỹ thuật cao từ $8,5%$ xuống dưới $2,5%$/năm.
    • Tăng năng suất lao động bình quân phân xưởng khai thác và sàng tuyển lên $12 - 15%$.
    • Tiết kiệm chi phí tuyển dụng và đào tạo lại ước tính 1,2 tỷ đồng mỗi năm.
    • Tối thiểu hóa rủi ro tai nạn lao động, bảo đảm vận hành liên tục cho Nhà máy Nhiệt điện Nông Sơn.

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

Hạn chế kỹ thuật

  1. Phạm vi không gian: Dữ liệu thu thập tại một đơn vị thành viên (Than Điện Nông Sơn), việc tổng quát hóa cho toàn bộ Tập đoàn VINACOMIN cần bổ sung thêm mẫu từ các mỏ than tại Quảng Ninh và Tây Nguyên.
  2. Dữ liệu cắt ngang (Cross-sectional Data): Khảo sát được thực hiện trong thời điểm từ tháng 02 đến tháng 05/2015; chưa theo dõi được sự biến động động lực theo chu kỳ kinh doanh dài hạn qua các năm tiếp theo.

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

  • Mở rộng ứng dụng mô hình cấu trúc tuyến tính SEM (Structural Equation Modeling) nhằm phân tích tác động trung gian của sự thỏa mãn công việc đến hiệu quả tài chính doanh nghiệp.
  • Xây dựng phần mềm HRM tích hợp thuật toán Machine Learning dự báo nguy cơ nghỉ việc (Churn Prediction) dựa trên biến thiên định kỳ của các chỉ số $X_1, X_3, X_5$.

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

                       CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI

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

1. Cỡ mẫu $N = 135$ có đảm bảo tính đại diện cho 495 nhân viên không?

Có. Theo công thức tính cỡ mẫu của Hair et al. ($N \ge 5 \times 27 = 135$) và công thức Yamane với mức ý nghĩa $e = 7,3%$, mẫu $N = 135$ phân bổ theo phương pháp ngẫu nhiên phân tầng bao phủ đầy đủ các phân xưởng khai thác, vận tải, sàng tuyển, cơ điện và khối phòng ban chức năng.

2. Tại sao nhân tố "Mối quan hệ đồng nghiệp và cấp trên" ($X_4$) bị loại khỏi mô hình?

Trong phân tích hồi quy, biến $X_4$ có $\text{Sig.} = 0,312 > 0,05$. Điều này phản ánh rằng trong môi trường công nghiệp nặng, người lao động xem sự hỗ trợ của đồng nghiệp là điều kiện hiển nhiên của công việc tổ đội, trong khi các yếu tố vật chất (thu nhập), sự an toàn tính mạng và cơ hội thăng tiến mới là động lực thúc đẩy nỗ lực vượt trội.

3. Làm thế nào để áp dụng kết quả nghiên cứu vào doanh nghiệp khác?

Các doanh nghiệp thuộc khối công nghiệp nặng (khai khoáng, xi măng, luyện kim, nhiệt điện) có thể sử dụng trực tiếp thang đo 27 biến quan sát của đề tài, thực hiện khảo sát nội bộ và chạy lại mô hình OLS để tìm ra trọng số $\beta$ đặc thù cho đơn vị mình.

4. Doanh nghiệp cần bao nhiêu ngân sách để triển khai gói giải pháp?

Doanh nghiệp nên trích từ $3 - 5%$ lợi nhuận sau thuế hàng năm. Khoản đầu tư này được hoàn vốn (ROI) trong vòng 12–18 tháng thông qua việc giảm chi phí tuyển dụng, giảm tỷ lệ nghỉ việc và gia tăng sản lượng khai thác.

5. Thang đo Likert 5 mức độ trong nghiên cứu được kiểm soát sai số như thế nào?

Bảng câu hỏi được chuẩn hóa qua bước phỏng vấn thử nghiệm ($n=10$) để loại bỏ câu hỏi đa nghĩa; kiểm tra hệ số Cronbach's Alpha loại bỏ các biến có tương quan biến - tổng $< 0,3$, và kiểm tra VIF $< 2,0$ để triệt tiêu hiện tượng đa cộng tuyến.

Kết luận

Khóa luận tốt nghiệp của tác giả Lương Thị Lệ Thủy dưới sự hướng dẫn của ThS. Hồ Khánh Ngọc Bích đã giải quyết trọn vẹn bài toán nâng cao động lực làm việc cho người lao động tại CTCP Than Điện Nông Sơn - VINACOMIN. Thông qua việc kết hợp các học thuyết quản trị nhân sự với mô hình phân tích định lượng chuyên sâu (Cronbach's Alpha, EFA, OLS Regression trên SPSS 20.0), đề tài chứng minh được 5 nhân tố then chốt tác động theo thứ tự ưu tiên: Lương thưởng ($\beta = 0,342$), An toàn lao động ($\beta = 0,258$), Bản chất công việc ($\beta = 0,205$), Đào tạo & thăng tiến ($\beta = 0,182$) và Môi trường làm việc ($\beta = 0,141$). Kết quả nghiên cứu không chỉ mang ý nghĩa học thuật mà còn là tài liệu tham khảo giá trị cho các nhà quản lý doanh nghiệp khai khoáng và nhiệt điện trong việc hoạch định chính sách nhân sự bền vững.