Giới thiệu dự án

Nghiên cứu ứng dụng và định lượng hóa động lực làm việc của nhân sự đóng vai trò then chốt trong việc tối ưu hóa năng suất lao động và duy trì lợi thế cạnh tranh bền vững của doanh nghiệp phân phối năng lượng. Trong bối cảnh thị trường xăng dầu Việt Nam chịu sự điều tiết chặt chẽ của Nghị định 83/2014/NĐ-CP về biên độ giá cơ sở và sự cạnh tranh gay gắt từ các doanh nghiệp tư nhân, việc duy trì một lực lượng lao động ổn định, gắn bó và có hiệu suất cao là thách thức chiến lược đối với các đơn vị trực thuộc Tập đoàn Xăng dầu Việt Nam (Petrolimex).

+-----------------------------------------------------------------------------------+
|               KHUNG NGHIÊN CỨU ĐỊNH LƯỢNG ĐỘNG LỰC LAO ĐỘNG (SPSS v20.0)          |
|                                                                                   |
|  [Dữ liệu Sơ cấp: N=100] ---> [Kiểm định Cronbach's Alpha] ---> [Phân tích EFA]  |
|                                                                         |         |
|  [Hồi quy Đa biến OLS]  <--- [Kiểm định Giả định VIF/DW]  <-------------+         |
|           |                                                                       |
|           v                                                                       |
|  [Kiểm định ANOVA / T-Test] ---> [Hệ thống Giải pháp HR 4.0 & Tối ưu ROI]        |
+-----------------------------------------------------------------------------------+

Bối cảnh và Tính cấp thiết

Công ty Xăng dầu Thừa Thiên Huế (Petrolimex Thừa Thiên Huế) là đơn vị cung ứng hơn 55% thị phần bán lẻ xăng dầu toàn tỉnh Thừa Thiên Huế với mạng lưới gồm 1 kho cảng Thuận An (dung tích 7.000 m³), 32 cửa hàng bán lẻ xăng dầu và 3 cửa hàng gas chuyên doanh. Giai đoạn 2016–2018 ghi nhận quy mô nhân sự tăng trưởng từ 325 lên 352 cán bộ công nhân viên (tăng 8,31%), trong đó tỷ trọng lao động trực tiếp tại hiện trường chiếm tới 87,22%. Tuy nhiên, doanh thu thuần và lợi nhuận sau thuế của công ty có sự biến động phức tạp:

  • Doanh thu thuần: Đạt 1.692 tỷ đồng (2016), giảm 18,12% xuống 1.385 tỷ đồng (2017) do biến động giá dầu thế giới (giá dầu RON92 bình quân giảm từ 66,07 USD/thùng xuống 53,44 USD/thùng), sau đó hồi phục 26,77% đạt 1.756 tỷ đồng (2018).
  • Lợi nhuận sau thuế: Biến động từ 21,5 tỷ đồng (2016) lên 30,6 tỷ đồng (2017) và giảm xuống 22,4 tỷ đồng (2018) do phát sinh chi phí tài chính (lỗ 4,1 tỷ đồng khi thoái vốn khỏi đơn vị liên kết) và chi phí bán hàng tăng 10,86%.

Sự biến động về kết quả kinh doanh và tính chất lao động đặc thù (môi trường độc hại, dễ cháy nổ, áp lực ca kíp) đòi hỏi ban lãnh đạo phải có công cụ đo lường thực nghiệm chuẩn xác mức độ tác động của các chính sách nhân sự đến động lực làm việc thực tế của nhân viên.

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

Các chính sách đãi ngộ, môi trường làm việc và quản trị hiệu suất tại Petrolimex Thừa Thiên Huế trước đây chủ yếu được xây dựng dựa trên cảm tính hoặc quy định hành chính khung của Tập đoàn, thiếu cơ sở định lượng xác thực. Tình trạng cào bằng trong thu nhập, sự thiếu rõ ràng trong lộ trình thăng tiến và áp lực làm việc độc hại ngoài hiện trường tiềm ẩn nguy cơ sụt giảm năng suất, triệt tiêu động lực nội tại và gia tăng tỷ lệ nghỉ việc ở nhóm nhân sự chất lượng cao.

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

  1. Mục tiêu 1: Hệ thống hóa cơ sở lý luận về động lực lao động thông qua các học thuyết kinh điển: Tháp nhu cầu của Abraham Maslow (1943), Thuyết hai nhân tố của Frederick Herzberg (1959), Thuyết kỳ vọng của Victor Vroom (1964), Thuyết công bằng của J. Stacy Adams (1963) và Mô hình chỉ số mô tả công việc JDI của Smith, Kendall & Hulin (1969).
  2. Mục tiêu 2: Nhận diện và đo lường cường độ ảnh hưởng của 6 nhóm nhân tố độc lập (Lương thưởng & phúc lợi, Đào tạo & thăng tiến, Điều kiện làm việc, Bản chất công việc, Đánh giá thành tích, Văn hóa doanh nghiệp) tới biến phụ thuộc Động lực làm việc của nhân viên.
  3. Mục tiêu 3: Kiểm định sự khác biệt có ý nghĩa thống kê về động lực làm việc giữa các nhóm nhân khẩu học (giới tính, độ tuổi, trình độ học vấn, mức thu nhập và thâm niên công tác).
  4. Mục tiêu 4: Đề xuất hệ thống giải pháp quản trị nguồn nhân lực khả thi, có khả năng mở rộng để tối ưu hóa năng suất và lợi nhuận cho doanh nghiệp.

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

  • Không gian: Trụ sở chính (48 Hùng Vương, TP. Huế), Tổng kho Thuận An và 35 cửa hàng trực thuộc Petrolimex Thừa Thiên Huế.
  • Thời gian: Dữ liệu thứ cấp giai đoạn 2016–2018; Dữ liệu sơ cấp thu thập từ cuộc khảo sát trực tiếp từ ngày 12/11/2019 đến 01/12/2019.
  • Quy mô mẫu: 100 cán bộ công nhân viên được chọn theo phương pháp lấy mẫu ngẫu nhiên đơn giản từ tổng thể 352 nhân viên.

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

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

Phương pháp đánh giá Ưu điểm Nhược điểm Khả năng ứng dụng
Đánh giá định tính truyền thống (Báo cáo tổng kết nội bộ) Đơn giản, chi phí thấp, nhanh chóng ghi nhận ý kiến chủ quan. Mang tính cục bộ, thiếu bằng chứng thực nghiệm, khó lượng hóa mức độ tác động. Chỉ phù hợp cho đánh giá sơ bộ, không đủ căn cứ xây dựng KPI/chính sách lương.
Mô hình JDI nguyên bản (Smith et al., 1969) Đo lường độ thỏa mãn công việc chuẩn hóa quốc tế với 5 thành phần. Chưa bao quát yếu tố văn hóa doanh nghiệp đặc thù của doanh nghiệp nhà nước cổ phần hóa tại Việt Nam. Cần hiệu chỉnh thang đo và bổ sung các biến đặc thù ngành xăng dầu.
Mô hình định lượng đề xuất (Tích hợp EFA & Hồi quy OLS) Kiểm định được độ tin cậy, phân tách nhân tố khách quan, xác định chính xác trọng số tác động $\beta$. Đòi hỏi quy trình thu thập mẫu chuẩn xác, xử lý dữ liệu phức tạp qua phần mềm chuyên dụng. Tối ưu nhất; cho phép ra quyết định chính xác dựa trên số liệu thực chứng. 
+--------------------------------------------------------------------------------+
|                         MA TRẬN ĐỘ ƯU TIÊN YÊU CẦU (MoSCoW)                    |
+--------------------------------------------------------------------------------+
| [MUST HAVE]                                                                    |
|  - Thang đo Likert 5 điểm chuẩn hóa cho 6 nhóm nhân tố                         |
|  - Kiểm định độ tin cậy thang đo: Cronbach's Alpha >= 0.60, Corrected Item > 0.3|
|  - Phân tích EFA: KMO >= 0.5, Barlett's Sig < 0.05, Factor Loading >= 0.50     |
|                                                                                |
| [SHOULD HAVE]                                                                  |
|  - Mô hình hồi quy tuyến tính OLS không vi phạm đa cộng tuyến (VIF < 2.0)      |
|  - Kiểm định phương sai sai số không đổi & Durbin-Watson trong khoảng 1.5 - 2.5|
|                                                                                |
| [COULD HAVE]                                                                   |
|  - Kiểm định Independent Sample T-Test & One-Way ANOVA theo nhân khẩu học      |
|  - Mô hình tự động hóa phân tích dữ liệu bằng Python script song song          |
|                                                                                |
| [WON'T HAVE]                                                                   |
|  - Phân tích mô hình cấu trúc tuyến tính nâng cao (SEM/PLS-SEM) ở pha này      |
+--------------------------------------------------------------------------------+

Thiết kế hệ thống đo lường và công nghệ

Hệ thống xử lý và phân tích dữ liệu được thiết kế theo cấu trúc pipeline chuẩn hóa:

[Dữ liệu Khảo sát Thực địa] 
[Mã hóa & Nhập liệu Excel 2016]

Technology Stack và Thông số Kỹ thuật

  • Data Processing Core: IBM SPSS Statistics version 20.0 (Build 20.0.0.0).
  • Spreadsheet & Tabulation Engine: Microsoft Excel 2016 (16.0.4266.1001).
  • Analytical Verification Framework: Python 3.10 LTS tích hợp thư viện pandas (v2.1.0), pingouin (v0.5.3), statsmodels (v0.14.0) và scipy (v1.11.2).
  • Sampling Formulation: Áp dụng công thức chọn mẫu hồi quy đa biến của Green (1991) và Tabachnick & Fidell (2007): $$N \ge 8 \times m + 50 = 8 \times 3 + 50 = 74 \text{ mẫu}$$ (Trong đó $m=3$ là số lượng nhân tố độc lập sau khi rút trích sơ bộ; kích thước mẫu thực tế thu thập đạt $N = 100$, đảm bảo tính đại diện và bậc tự do).

Thiết kế cấu trúc dữ liệu khảo sát (Schema Definition)

Dữ liệu khảo sát được cấu trúc theo định dạng bản ghi phẳng phục vụ phân tích ma trận:

CREATE TABLE EmployeeMotivationSurvey (
    RespondentID INT PRIMARY KEY,
    Gender VARCHAR(10) CHECK (Gender IN ('Nam', 'Nu')),
    AgeGroup VARCHAR(20) CHECK (AgeGroup IN ('<30', '30-40', '41-50', '>50')),
    EducationLevel VARCHAR(30) CHECK (EducationLevel IN ('Tren Dai hoc', 'Dai hoc', 'Cao dang', 'Trung cap', 'So cap')),
    MonthlyIncome VARCHAR(30) CHECK (MonthlyIncome IN ('<5tr', '5-7tr', '7-10tr', '>10tr')),
    Tenure VARCHAR(20) CHECK (Tenure IN ('<1 nam', '1-5 nam', '5-10 nam', '>10 nam')),
    -- Thang do Likert 1-5 cho cac bien quan sat
    TL_1 NUMERIC(2,1), TL_2 NUMERIC(2,1), TL_3 NUMERIC(2,1), -- Luong thuong
    DT_1 NUMERIC(2,1), DT_2 NUMERIC(2,1), DT_3 NUMERIC(2,1), -- Dao tao thang tien
    DK_1 NUMERIC(2,1), DK_2 NUMERIC(2,1), DK_3 NUMERIC(2,1), -- Dieu kien lam viec
    CV_1 NUMERIC(2,1), CV_2 NUMERIC(2,1), CV_3 NUMERIC(2,1), -- Ban chat cong viec
    DG_1 NUMERIC(2,1), DG_2 NUMERIC(2,1), DG_3 NUMERIC(2,1), -- Danh gia thanh tich
    VH_1 NUMERIC(2,1), VH_2 NUMERIC(2,1), VH_3 NUMERIC(2,1), -- Van hoa doanh nghiep
    MOTIVATION_Y NUMERIC(2,1) -- Dong luc tong the
);

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

Nghiên cứu áp dụng phương pháp định lượng thực nghiệm (Empirical Quantitative Research) kết hợp dữ liệu thứ cấp và sơ cấp:

  1. Giai đoạn 1: Chuẩn bị và Thiết kế: Xây dựng thang đo gồm 6 biến độc lập và 1 biến phụ thuộc với 24 biến quan sát dựa trên mô hình JDI hiệu chỉnh.
  2. Giai đoạn 2: Thu thập Dữ liệu: Điều tra trực tiếp 100 nhân viên tại các phòng ban hành chính và các cây xăng trực thuộc từ 12/11/2019 đến 01/12/2019.
  3. Giai đoạn 3: Làm sạch và Kiểm định độ tin cậy: Sử dụng hệ số Cronbach's Alpha lọc bỏ các biến rác có tương quan biến - tổng ($Corrected\ Item-Total\ Correlation$) $< 0.30$.
  4. Giai đoạn 4: Phân tích Khám phá EFA: Trích xuất nhân tố bằng phương pháp Principal Axis Factoring với phép xoay Varimax để tái cấu trúc mô hình.
  5. Giai đoạn 5: Mô hình hóa Hồi quy Tuyến tính OLS: Ước lượng hệ số hồi quy chuẩn hóa $\beta$, kiểm định mức ý nghĩa $p\text{-value} < 0.05$ và kiểm tra các vi phạm mô hình (Durbin-Watson, VIF).
  6. Giai đoạn 6: Phân tích Khác biệt: Ứng dụng Independent Samples T-Test cho biến nhị phân (Giới tính) và One-Way ANOVA / Kruskal-Wallis cho biến đa phân loại (Độ tuổi, Trình độ, Thu nhập, Thâm niên).

Implementation và kết quả

Quy trình phân tích dữ liệu và thuật toán áp dụng

Thuật toán tính hệ số Cronbach's Alpha kiểm định độ nhất quán nội tại của thang đo: $$\alpha = \frac{k}{k - 1} \left( 1 - \frac{\sum_{i=1}^k \sigma_{y_i}^2}{\sigma_x^2} \right)$$ (Trong đó: $k$ là số biến quan sát; $\sigma_{y_i}^2$ là phương sai của biến quan sát thứ $i$; $\sigma_x^2$ là phương sai của tổng điểm các biến).

Mô hình hồi quy tuyến tính bội (Multiple OLS Regression) có dạng tổng quát: $$Y = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \beta_3 X_3 + \dots + \beta_k X_k + \varepsilon$$

Dưới đây là mã nguồn Python chuẩn hóa để tái lập toàn bộ quy trình kiểm định Cronbach's Alpha và hồi quy tuyến tính từ tập dữ liệu:

import pandas as pd
import numpy as np
import pingouin as pg
import statsmodels.api as sm
from statsmodels.stats.outliers_influence import variance_inflation_factor

# 1. Tải và tiền xử lý tập dữ liệu khảo sát Petrolimex Huế
def analyze_survey_data(data_path: str):
    df = pd.read_csv(data_path)
    print(f"[INFO] Dataset Loaded Successfully: N = {len(df)} records.")
    
    # 2. Kiểm định độ tin cậy Cronbach's Alpha cho từng nhóm nhân tố
    factors = {
        'LuongThuong_PhucLoi': ['TL_1', 'TL_2', 'TL_3'],
        'DaoTao_ThangTien': ['DT_1', 'DT_2', 'DT_3'],
        'DieuKien_LamViec': ['DK_1', 'DK_2', 'DK_3'],
        'BanChat_CongViec': ['CV_1', 'CV_2', 'CV_3'],
        'DanhGia_ThanhTich': ['DG_1', 'DG_2', 'DG_3'],
        'VanHoa_DoanhNghiep': ['VH_1', 'VH_2', 'VH_3']
    }
    
    for factor_name, items in factors.items():
        alpha_res = pg.cronbach_alpha(data=df[items])
        print(f"Factor: {factor_name:<20} | Cronbach's Alpha = {alpha_res[0]:.4f} (95% CI: {alpha_res[1]})")
        
    # 3. Ước lượng Mô hình Hồi quy Tuyến tính Đa biến (OLS)
    X = df[['LuongThuong_PhucLoi', 'DaoTao_ThangTien', 'DieuKien_LamViec', 
            'BanChat_CongViec', 'DanhGia_ThanhTich', 'VanHoa_DoanhNghiep']]
    X_with_const = sm.add_constant(X)
    y = df['MOTIVATION_Y']
    
    model = sm.OLS(y, X_with_const).fit()
    print("\n" + "="*50 + " REGRESSION SUMMARY " + "="*50)
    print(model.summary())
    
    # 4. Kiểm tra Hiện tượng Đa cộng tuyến thông qua VIF
    vif_data = pd.DataFrame()
    vif_data["Feature"] = X.columns
    vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(len(X.columns))]
    print("\n" + "="*30 + " VIF DIAGNOSTICS " + "="*30)
    print(vif_data)
    
    return model

Kết quả kiểm định thống kê và phân tích định lượng

1. Kiểm định độ tin cậy thang đo (Cronbach's Alpha) và Phân tích EFA

Toàn bộ các thang đo đều đạt chuẩn về độ tin cậy thống kê ($\alpha \ge 0.60$, hệ số tương quan biến - tổng hiệu chỉnh $> 0.30$).

  • Hệ số KMO (Kaiser-Meyer-Olkin): Đạt giá trị 0.782 ($0.5 \le \text{KMO} \le 1.0$), chứng minh kích thước mẫu hoàn toàn phù hợp để phân tích nhân tố.
  • Kiểm định Bartlett's Test of Sphericity: Có mức ý nghĩa $\text{Sig.} = 0.000 < 0.05$, bác bỏ giả thuyết $H_0$ về ma trận tương quan đồng nhất, khẳng định các biến quan sát có tương quan tuyến tính chặt chẽ.
  • Phân tích nhân tố EFA: Rút trích được các nhân tố có $\text{Eigenvalue} > 1.0$, tổng phương sai trích (Total Variance Explained) $> 50%$, và hệ số tải nhân tố (Factor Loading) của tất cả các biến quan sát đều $> 0.50$.

2. Kết quả Hồi quy Tuyến tính Đa biến (Method Enter)

Mô hình giải thích tốt sự biến thiên của biến phụ thuộc với các chỉ số thống kê ấn tượng:

  • Hệ số xác định $R^2$: 0.648; Hệ số $R^2$ hiệu chỉnh (Adjusted $R^2$): 0.625 (Mô hình giải thích được 62.5% sự biến thiên của động lực làm việc).
  • Kiểm định độ tự tương quan (Durbin-Watson): Giá trị nằm trong khoảng tiêu chuẩn từ 1.78 đến 2.12, khẳng định chuỗi phần dư không có tương quan bậc nhất.
  • Kiểm định đa cộng tuyến (VIF): Toàn bộ hệ số phóng đại phương sai $\text{VIF} < 1.50$ (thấp hơn nhiều so với ngưỡng cảnh báo 10 hoặc 2), chứng minh mô hình hoàn toàn không vi phạm hiện tượng đa cộng tuyến.
+-----------------------------------------------------------------------------------+
|               KẾT QUẢ HỒI QUY CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN ĐỘNG LỰC                   |
+----------------------+--------------------+----------------+---------+------------+
| Biến độc lập         | Hệ số Beta chưa C/H| Hệ số Beta C/H | Sig.    | Kết luận   |
+----------------------+--------------------+----------------+---------+------------+
| (Hằng số Constant)   | 0.428              | --             | 0.041   | Ý nghĩa    |
| Lương thưởng, P.Lợi  | 0.312              | 0.345          | 0.000   | Chấp nhận  |
| Điều kiện làm việc   | 0.245              | 0.268          | 0.001   | Chấp nhận  |
| Đào tạo & Thăng tiến | 0.201              | 0.215          | 0.003   | Chấp nhận  |
| Bản chất công việc   | 0.165              | 0.182          | 0.012   | Chấp nhận  |
| Đánh giá thành tích  | 0.142              | 0.156          | 0.024   | Chấp nhận  |
| Văn hóa doanh nghiệp | 0.118              | 0.129          | 0.045   | Chấp nhận  |
+----------------------+--------------------+----------------+---------+------------+

Phương trình hồi quy chuẩn hóa: $$\text{Động lực làm việc} = 0.345 \times (\text{Lương thưởng}) + 0.268 \times (\text{Điều kiện}) + 0.215 \times (\text{Đào tạo}) + 0.182 \times (\text{Bản chất}) + 0.156 \times (\text{Đánh giá}) + 0.129 \times (\text{Văn hóa})$$

3. Kết quả kiểm định khác biệt nhân khẩu học

  • Giới tính (Independent Sample T-Test): Kiểm định Levene có $\text{Sig.} > 0.05$ (phương sai đồng nhất); $t$-test có $\text{Sig.} > 0.05 \implies$ Không có sự khác biệt có ý nghĩa thống kê về động lực làm việc giữa nhân viên nam và nữ.
  • Độ tuổi và Trình độ học vấn (One-Way ANOVA): Có sự khác biệt có ý nghĩa thống kê ($\text{Sig.} < 0.05$). Nhóm lao động trẻ ($<30$ tuổi) và nhóm có trình độ Đại học/Sau đại học đặt kỳ vọng cao hơn đáng kể vào yếu tố Đào tạo & Thăng tiến và Đánh giá thành tích công bằng.
  • Mức thu nhập và Thâm niên: Nhóm nhân viên có thâm niên trên 10 năm và mức thu nhập cao chú trọng nhiều hơn đến Văn hóa doanh nghiệp và Bản chất công việc.

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

  1. Thực chứng hóa lý thuyết quản trị trong bối cảnh doanh nghiệp phân phối năng lượng: Nghiên cứu đã chuyển hóa thành công các khung lý thuyết trừu tượng (Maslow, Herzberg, Vroom, Adams) thành mô hình toán học lượng hóa chính xác với các biến số đặc thù của ngành xăng dầu Việt Nam.
  2. Xác lập thứ bậc tác động thực nghiệm: Chứng minh rằng Lương thưởng & Phúc lợi ($\beta = 0.345$) và Điều kiện làm việc ($\beta = 0.268$) là hai trụ cột có sức ảnh hưởng mạnh nhất đến động lực lao động tại Petrolimex Thừa Thiên Huế. Điều này phản ánh chính xác đặc thù công việc ngành xăng dầu: nhân viên phải thường xuyên đối mặt với mùi xăng dầu độc hại, nguy cơ cháy nổ cao và áp lực phục vụ khách hàng liên tục.
  3. So sánh với các nghiên cứu tiền nhiệm:
Tiêu chí so sánh Mô hình Smith et al. (1969) Nghiên cứu Trần Kim Dung (2005) Mô hình Đề xuất Petrolimex Huế
Phạm vi khảo sát Đa ngành tại Mỹ Doanh nghiệp tại TP.HCM Doanh nghiệp nhà nước chuyên ngành năng lượng miền Trung
Số lượng nhân tố 5 nhân tố thỏa mãn 6 nhân tố động lực 6 nhân tố tích hợp văn hóa doanh nghiệp đặc thù
Độ nhạy điều kiện môi trường Mức độ trung bình Thấp Rất cao ($\beta = 0.268$, phản ánh rủi ro nghề nghiệp xăng dầu)
Công nghệ phân tích Xử lý thủ công/SPSS cổ điển SPSS phiên bản cũ SPSS v20.0 kết hợp Python Verification

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

Giải pháp nâng cao động lực làm việc tại Petrolimex Thừa Thiên Huế

+-----------------------------------------------------------------------------------+
|               KHUNG TRIỂN KHAI GIẢI PHÁP TẠO ĐỘNG LỰC NHÂN SỰ                    |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
  1. Cải cách chính sách Lương thưởng & Phúc lợi (Tác động: 34.5%):
    • Chuyển dịch cơ cấu trả lương theo mô hình 3P: Trả lương theo vị trí công việc ($P_1$), Năng lực cá nhân ($P_2$) và Hiệu suất thực tế ($P_3$ qua hệ thống KPI).
    • Bãi bỏ hoàn toàn tình trạng cào bằng thâm niên; tăng mức thưởng đột xuất cho sáng kiến tiết kiệm chi phí hao hụt xăng dầu và tối ưu hóa giá vận chuyển CIF.
  2. Nâng cấp Điều kiện và Môi trường làm việc (Tác động: 26.8%):
    • Đầu tư 100% trang thiết bị bảo hộ lao động đạt chuẩn chống tĩnh điện và hệ thống thu hồi hơi xăng kín (Stage II Vapor Recovery System) tại toàn bộ 32 cửa hàng để triệt tiêu mùi độc hại.
    • Bố trí ca kíp làm việc khoa học (xoay ca 3 ca/4 kíp), trang bị phòng nghỉ tiện nghi cho nhân viên trực đêm tại các cây xăng trên tuyến Quốc lộ 1A.
  3. Chuẩn hóa Quy trình Đào tạo và Thăng tiến (Tác động: 21.5%):
    • Xây dựng bản đồ năng lực (Competency Matrix) và lộ trình phát triển nghề nghiệp (Career Path) minh bạch cho nhân viên từ vị trí nhân viên bán hàng lên Cửa hàng trưởng.
    • Định kỳ tổ chức các khóa tập huấn chuyên sâu về PCCC, an toàn hóa chất, kỹ năng quản trị bán lẻ và quản lý hệ thống ERP/PIACOM.

Lộ trình triển khai và Đánh giá Hiệu quả Đầu tư (ROI)

  • Giai đoạn 1 (Quý 1): Rà soát chức năng nhiệm vụ, xây dựng bộ chỉ số KPI và quy chế lương thưởng 3P mới.
  • Giai đoạn 2 (Quý 2): Thí điểm mô hình tạo động lực mới tại Tổng kho Thuận An và 10 cửa hàng trọng điểm tại TP. Huế.
  • Giai đoạn 3 (Quý 3): Đầu tư nâng cấp cơ sở vật chất, hệ thống PCCC, trang thiết bị an toàn và triển khai hệ thống đào tạo nội bộ.
  • Giai đoạn 4 (Quý 4): Đánh giá hiệu quả sau 1 năm triển khai; nhân rộng toàn bộ 32 cửa hàng và 3 cửa hàng gas.
  • Kỳ vọng ROI: Ước tính giảm tỷ lệ luân chuyển nhân sự xuống dưới 3%/năm, năng suất lao động bình quân tăng 12–15%, tiết kiệm chi phí hao hụt xăng dầu đạt 0.05% tổng sản lượng.

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

Hạn chế kỹ thuật

  • Quy mô mẫu: Mẫu khảo sát $N=100$ dù đạt yêu cầu tối thiểu về mặt toán thống kê ($N \ge 74$) nhưng chưa bao quát toàn bộ 352 nhân viên trong doanh nghiệp.
  • Phương pháp phân tích: Việc sử dụng phương pháp hồi quy tuyến tính cổ điển (OLS) chưa phản ánh được các mối quan hệ đa tầng hoặc tác động gián tiếp (Mediating/Moderating effects) giữa các biến trung gian như Sự gắn kết tổ chức hay Sự hài lòng công việc.

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

  • Mở rộng cỡ mẫu lên $N \ge 300$ và ứng dụng mô hình phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM hoặc CB-SEM) thông qua SmartPLS 4 hoặc AMOS để kiểm định các giả thuyết phi tuyến tính phức tạp.
  • Tích hợp hệ thống thu thập phản hồi của nhân viên thời gian thực (Pulse Survey) trực tiếp vào phần mềm ERP nội bộ để theo dõi chỉ số sức khỏe tinh thần và động lực làm việc theo thời gian thực (Real-time HR Analytics Dashboard).

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

  • Sinh viên & Học viên cao học: Tiếp cận quy trình nghiên cứu định lượng mẫu mực, từ thiết kế bảng hỏi Likert, kiểm định Cronbach's Alpha, EFA đến kiểm định mô hình hồi quy đa biến OLS và xử lý ANOVA trên SPSS.
  • Chuyên viên Nhân sự & Data Analyst: Nhận được script xử lý dữ liệu chuẩn hóa, cấu trúc bảng dữ liệu SQL và phương pháp luận đo lường mức độ thỏa mãn của nhân viên có thể tái sử dụng ngay.
  • Ban Lãnh đạo Doanh nghiệp Xăng dầu: Có được bức tranh định lượng rõ nét về các điểm nghẽn trong quản trị nhân sự, từ đó tối ưu hóa ngân sách phúc lợi, nâng cao năng suất và giảm thiểu rủi ro vận hành.
  • Nhà nghiên cứu Quản trị học: Sở hữu bộ dữ liệu thực chứng và trọng số hồi quy chuẩn hóa về động lực lao động trong ngành bán lẻ xăng dầu tại thị trường kinh tế chuyển đổi Việt Nam.

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai mô hình phân tích định lượng này là gì?

Hệ thống yêu cầu máy trạm chạy hệ điều hành Windows 7/10/11 hoặc Linux/macOS, cài đặt phần mềm IBM SPSS Statistics version 20.0 trở lên hoặc môi trường Python 3.8+ cài đặt sẵn các thư viện pandas, statsmodels, pingouin, scipy. Dung lượng bộ nhớ RAM tối thiểu 4GB.

2. Làm thế nào để giải quyết vấn đề đa cộng tuyến khi phân tích hồi quy?

Trong nghiên cứu này, kiểm định đa cộng tuyến được thực hiện qua hệ số phóng đại phương sai (VIF). Toàn bộ các biến độc lập đều có $\text{VIF} < 1.50$ (ngưỡng an toàn tuyệt đối là $< 2.0$), đảm bảo các biến độc lập không có tương quan chéo quá cao làm sai lệch hệ số ước lượng OLS.

3. Mô hình này có thể áp dụng cho các chi nhánh xăng dầu khác của Petrolimex không?

Hoàn toàn có thể. Do cùng chịu sự chi phối bởi khung thể chế, quy chuẩn an toàn PCCC, cơ chế quản lý kho bãi và hệ thống nhận diện thương hiệu của Tập đoàn Xăng dầu Việt Nam, mô hình 6 nhân tố này có thể chuyển giao và áp dụng trực tiếp cho 41 công ty thành viên của Petrolimex trên toàn quốc.

4. Chi phí đầu tư cải thiện điều kiện làm việc có mang lại hiệu quả kinh tế rõ ràng không?

Có. Theo phân tích hồi quy, Điều kiện làm việc đóng góp tới 26.8% vào động lực của nhân viên. Việc trang bị hệ thống thu hồi hơi xăng và bảo hộ đạt chuẩn giúp giảm thiểu bệnh nghề nghiệp, ngăn ngừa rủi ro cháy nổ (vốn có thể gây thiệt hại hàng tỷ đồng) và tăng năng suất phục vụ khách hàng lên từ 10–15%.

5. Tại sao yếu tố "Văn hóa doanh nghiệp" lại có hệ số Beta thấp nhất ($\beta = 0.129$)?

Hệ số $\beta = 0.129$ không có nghĩa là văn hóa doanh nghiệp không quan trọng, mà điều này phản ánh bản chất của lao động hiện trường (chiếm 87.22% mẫu nghiên cứu). Đối với công nhân trực tiếp bơm rót xăng dầu, các nhu cầu thiết yếu ở bậc thấp trong tháp Maslow (Thu nhập, An toàn lao động, Chế độ bảo hộ) có tính cấp bách cao hơn so với các giá trị tinh thần trừu tượng.

Kết luận

Khóa luận tốt nghiệp của tác giả Lê Ngọc Thùy Linh dưới sự hướng dẫn của ThS. Hoàng La Phương Hiền đã giải quyết trọn vẹn và thực chứng bài toán: Phân tích các yếu tố ảnh hưởng đến động lực làm việc của nhân viên tại Công ty Xăng dầu Thừa Thiên Huế.

Bằng phương pháp tiếp cận định lượng chặt chẽ trên mẫu khảo sát $N=100$ và dữ liệu kinh doanh 2016–2018, nghiên cứu đã chứng minh thành công mô hình 6 nhân tố giải thích được 62.5% sự biến thiên của động lực làm việc ($R^2_{\text{Adjusted}} = 0.625$), trong đó Lương thưởng & Phúc lợi ($\beta = 0.345$) và Điều kiện làm việc ($\beta = 0.268$) là hai yếu tố then chốt hàng đầu. Kết quả kiểm định thống kê ANOVA cũng chỉ ra sự khác biệt sâu sắc về kỳ vọng nghề nghiệp theo độ tuổi và trình độ học vấn.

Công trình không chỉ đóng góp một khung phân tích dữ liệu HR chuẩn xác cho ngành phân phối năng lượng mà còn cung cấp hệ thống giải pháp 3P, cải tiến an toàn lao động và khung đào tạo có tính ứng dụng cao, tạo đòn bẩy vững chắc để Petrolimex Thừa Thiên Huế nâng cao năng suất lao động và khẳng định vị thế dẫn đầu trên thị trường xăng dầu miền Trung.