Giới thiệu dự án

Ngành dệt may Việt Nam giữ vị trí trụ cột trong cơ cấu xuất khẩu quốc gia, đạt kim ngạch 35 tỷ USD vào năm 2018 với tốc độ tăng trưởng kép hàng năm (CAGR) giai đoạn 2011–2018 đạt 12%. Trong đó, sản phẩm may mặc chiếm tỷ trọng chi phối 80%, với thị trường xuất khẩu trọng điểm là Hoa Kỳ (chiếm 47% tổng kim ngạch toàn ngành). Mặc dù hưởng lợi lớn từ làn sóng dịch chuyển chuỗi cung ứng toàn cầu sau căng thẳng thương mại Mỹ - Trung, phần lớn doanh nghiệp dệt may trong nước vẫn hoạt động ở phân khúc gia công cắt - may - hoàn thiện (CMT) và FOB cấp 1 với biên lợi nhuận mỏng. Trước áp lực gia tăng chi phí nhân công và cạnh tranh gay gắt từ các quốc gia như Bangladesh và Campuchia, tối ưu hóa năng suất lao động (NSLĐ) trở thành bài toán sống còn.

Công ty Scavi Huế (thành viên của Tập đoàn Corele International - Pháp, Top 2 tại Pháp và Top 7 thế giới về dịch vụ outsourcing thời trang nội y) sở hữu quy mô 3 nhà máy tại Khu công nghiệp Phong Điền với 6.310 lao động (trong đó 5.896 lao động trực tiếp, chiếm 93,44% tổng nhân sự năm 2018). Mặc dù doanh thu năm 2018 đạt 2.583,51 tỷ VNĐ và lợi nhuận sau thuế đạt 227,42 tỷ VNĐ, tốc độ tăng trưởng lợi nhuận đã chững lại ở mức 4,50% (so với mức tăng 97,32% của năm 2017), không đạt kỳ vọng tăng trưởng 25%/năm của Tập đoàn.

Vấn đề cốt lõi đặt ra là sự thiếu hụt mô hình lượng hóa chuẩn xác về mức độ tác động của các nhân tố vi mô (tâm lý, môi trường, kỹ thuật, quản trị) đến năng suất thực tế của công nhân trực tiếp, dẫn đến tình trạng quản lý định tính và phân bổ nguồn lực chưa tối ưu.

+-----------------------------------------------------------------------------------+
|                        MỤC TIÊU VÀ PHẠM VI NGHIÊN CỨU                             |
+-----------------------------------------------------------------------------------+
| 1. Hệ thống hóa cơ sở lý luận về năng suất lao động trong công nghiệp may mặc.     |
| 2. Đo lường định lượng mức độ ảnh hưởng của 5 nhóm nhân tố cốt lõi đến NSLĐ.       |
| 3. Xây dựng phương trình hồi quy tuyến tính chuẩn hóa làm căn cứ ra quyết định.  |
| 4. Đề xuất gói giải pháp công thái học, kỹ thuật và quản trị nâng cao NSLĐ.      |
+-----------------------------------------------------------------------------------+

Dự án áp dụng phương pháp tiếp cận định lượng hiện đại dựa trên lý thuyết tiếp cận toàn diện của Ủy ban Năng suất châu Âu (EPA) và mô hình hồi quy tuyến tính đa biến. Nghiên cứu giới hạn phạm vi khảo sát tại Nhà máy 1 và Nhà máy 2 của Scavi Huế trong giai đoạn quý IV/2019, sử dụng dữ liệu sản xuất kinh doanh thứ cấp chuỗi thời gian 2016–2018.

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

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

Tại các doanh nghiệp dệt may quy mô lớn, các phương pháp đánh giá năng suất truyền thống thường tập trung vào phương pháp đo đạc thời gian bấm giờ (time-study) hoặc định mức kỹ thuật đơn thuần, bỏ qua các biến số hành vi và điều kiện môi trường làm việc.

Tiêu chí phân tích Phương pháp truyền thống (Định mức kỹ thuật) Mô hình kinh tế lượng EFA & OLS (Đề xuất)
Bản chất đo lường Đơn biến, thuần cơ học thao tác Đa biến, kết hợp tâm lý học lao động & quản trị
Xác định nguyên nhân Chỉ phản ánh kết quả sản lượng đầu ra Lượng hóa chính xác trọng số từng yếu tố
Độ tin cậy dữ liệu Dễ sai lệch do hiệu ứng Hawthorne Kiểm định chặt chẽ (Cronbach's $\alpha$, KMO, Bartlett)
Khả năng dự báo Kém, mang tính phản ứng thụ động Dự báo độ co giãn của NSLĐ theo từng biến số

Áp dụng khung ưu tiên hóa yêu cầu phân tích theo tiêu chuẩn MoSCoW:

  • Must-have: Xác định cấu trúc nhân tố thông qua EFA; kiểm định đa cộng tuyến (VIF < 2.0); kiểm định tự tương quan Durbin-Watson.
  • Should-have: Đánh giá thực trạng 24 biến quan sát theo thang đo Likert 5 mức độ; phân tích ANOVA sai khác nhân khẩu học.
  • Could-have: Tích hợp thuật toán hồi quy tự động hóa trên nền tảng ngôn ngữ lập trình thống kê.
  • Won't-have (giai đoạn này): Phân tích chuỗi thời gian ARIMA cho biến động năng suất theo mùa vụ.

Thiết kế hệ thống

Kiến trúc phân tích dữ liệu và mô hình hóa kinh tế lượng được thiết kế theo quy trình luồng chuẩn hóa:

Cấu trúc biến số trong mô hình nghiên cứu bao gồm:

  • Biến phụ thuộc ($Y$ - NSLĐ): Đo lường qua mức độ hoàn thành kế hoạch, chất lượng thành phẩm, tỷ lệ sai hỏng và sự hăng hái sản xuất (3 biến quan sát).
  • Biến độc lập ($X_1 \dots X_5$): Gồm 21 biến quan sát chia thành 5 nhóm:
    • $X_1$: Bản thân người lao động (Trình độ, tay nghề, kỷ luật, sức khỏe)
    • $X_2$: Sự quản lý và phân công lao động của cấp trên (Giao việc, hỗ trợ kỹ thuật, khen thưởng)
    • $X_3$: Điều kiện làm việc (Ánh sáng, vi khí hậu, độ rung, an toàn vệ sinh lao động)
    • $X_4$: Sự cải tiến trong sản xuất (Kaizen, thiết bị tự động hóa, cữ gá dưỡng, bố trí chuyền Lean)
    • $X_5$: Môi trường làm việc (Văn hóa doanh nghiệp, quan hệ đồng nghiệp, áp lực công việc)

Methodology

  • Phương pháp thu thập dữ liệu: Khảo sát trực tiếp công nhân tại chuyền may và khu vực hoàn thiện trong các ca sản xuất. Quy mô mẫu $N = 120$ được xác định dựa trên hai tiêu chuẩn kiểm nghiệm:
    1. Theo Hair et al. (1998): $N \ge 5 \times k = 5 \times 24 = 120$ (với $k = 24$ biến quan sát).
    2. Theo Tabachnick & Fidell (1991): $N \ge 8m + 50 = 8 \times 5 + 50 = 90$ (với $m = 5$ biến độc lập). Kích thước mẫu $N = 120$ thỏa mãn đồng thời cả hai điều kiện toán học.
  • Công cụ và nền tảng xử lý:
    • IBM SPSS Statistics version 20.0 (phân tích nhân tố EFA, độ tin cậy, hồi quy OLS).
    • Python version 3.8.10 kết hợp thư viện Pandas 1.2.4 và Statsmodels 0.12.2 (xây dựng pipeline kiểm tra giả định).
  • Tiêu chuẩn kiểm soát chất lượng: Loại bỏ mẫu lỗi (missing values, straight-lining), chuẩn hóa dữ liệu Z-score trước khi đưa vào ma trận xoay.

Implementation và kết quả

Development process

Quá trình phân tích kinh tế lượng được thực thi qua 4 giai đoạn logic với các tiêu chuẩn toán học nghiêm ngặt:

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

def validate_econometric_model(df, independent_cols, dependent_col):
    """
    Pipeline thực thi kiểm định giả định và hồi quy tuyến tính OLS
    cho mô hình năng suất lao động Scavi Huế.
    """
    X = df[independent_cols]
    y = df[dependent_col]
    
    # 1. Bổ sung hệ số chặn (Intercept / Constant)
    X_with_const = sm.add_constant(X)
    
    # 2. Ước lượng mô hình Ordinary Least Squares (OLS)
    model = sm.OLS(y, X_with_const).fit()
    
    # 3. Tính toán Variance Inflation Factor (VIF) kiểm tra đa cộng tuyến
    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))]
    
    # 4. Trích xuất chỉ số Durbin-Watson kiểm tra tự tương quan
    dw_stat = sm.stats.durbin_watson(model.resid)
    
    return {
        "summary": model.summary(),
        "vif": vif_data,
        "durbin_watson": dw_stat,
        "params": model.params,
        "pvalues": model.pvalues,
        "rsquared_adj": model.rsquared_adj
    }

Testing và validation

Toàn bộ thang đo trải qua kiểm định hai bước: kiểm định tính nhất quán nội tại (Cronbach's Alpha) và kiểm định giá trị hội tụ/phân biệt (EFA).

+--------------------------------------------------------------------------------------+
|                           BẢNG KIỂM ĐỊNH ĐỘ TIN CẬY THANG ĐO                         |
+------------------------------------+-----------------+------------------+------------+
| Thang đo nhân tố                   | Số biến quan sát| Cronbach's Alpha | Đánh giá   |
+------------------------------------+-----------------+------------------+------------+
| Bản thân người lao động (NL)       | 5               | 0.814            | Rất tốt    |
| Quản lý và phân công (QL)          | 4               | 0.792            | Tốt        |
| Điều kiện làm việc (DK)            | 4               | 0.768            | Tốt        |
| Cải tiến trong sản xuất (CT)       | 4               | 0.825            | Rất tốt    |
| Môi trường làm việc (MT)           | 4               | 0.781            | Tốt        |
| Năng suất lao động (NS)            | 3               | 0.803            | Rất tốt    |
+------------------------------------+-----------------+------------------+------------+

Tất cả các biến quan sát đều có hệ số tương quan biến - tổng (Corrected Item-Total Correlation) $> 0.35$, không có biến nào bị loại.

Kết quả phân tích nhân tố khám phá EFA (biến độc lập):

  • Hệ số kiểm định Kaiser-Meyer-Olkin (KMO) = $0.782$ ($> 0.50$), chứng minh kích thước mẫu hoàn toàn phù hợp với phân tích nhân tố.
  • Kiểm định Bartlett's Test of Sphericity: Chi-Square = $842.16$, mức ý nghĩa $\text{Sig.} = 0.000 < 0.001$, khẳng định các biến quan sát có tương quan tuyến tính có ý nghĩa thống kê.
  • Giá trị Eigenvalue của nhân tố thứ 5 = $1.214 > 1.0$; tổng phương sai trích (Total Variance Explained) đạt $62.45% > 50%$, chứng minh 5 nhân tố giải thích được 62,45% độ biến thiên của bộ dữ liệu gốc.
  • Ma trận xoay nhân tố Varimax cho thấy hệ số tải nhân tố (Factor Loading) của tất cả 21 biến quan sát đều $> 0.55$.

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

Phương trình hồi quy tuyến tính đa biến chuẩn hóa sau khi kiểm định không vi phạm các giả định thống kê ($1.5 < \text{Durbin-Watson} = 1.874 < 2.5$; hệ số $\text{VIF} < 1.45$; phần dư chuẩn hóa có phân phối chuẩn):

$$\text{NSLĐ} = \beta_0 + 0.312 \cdot \text{CT} + 0.264 \cdot \text{QL} + 0.218 \cdot \text{MT} + 0.185 \cdot \text{NL} + 0.142 \cdot \text{DK} + e_i$$

+------------------------------------------------------------------------------------+
|                         KẾT QUẢ PHÂN TÍCH HỒI QUY ĐA BIẾN                          |
+------------------------------+--------------------+---------------+----------------+
| Biến độc lập                 | Hệ số Beta chuẩn hóa| T-statistic   | P-value (Sig.) |
+------------------------------+--------------------+---------------+----------------+
| Cải tiến trong sản xuất (CT) | 0.312              | 4.128         | 0.000          |
| Quản lý và phân công (QL)    | 0.264              | 3.491         | 0.001          |
| Môi trường làm việc (MT)     | 0.218              | 2.915         | 0.004          |
| Bản thân người lao động (NL) | 0.185              | 2.462         | 0.015          |
| Điều kiện làm việc (DK)      | 0.142              | 2.018         | 0.046          |
+------------------------------+--------------------+---------------+----------------+
| R-Square: 0.584 | Adjusted R-Square: 0.566 | F-Statistic: 32.042 (Sig. = 0.000)   |
+------------------------------------------------------------------------------------+

Hệ số $R^2$ hiệu chỉnh đạt $0.566$, chỉ ra rằng 5 nhóm nhân tố đưa vào mô hình giải thích được $56.6%$ sự biến động trong năng suất lao động của công nhân may trực tiếp tại Scavi Huế. Trong đó, yếu tố "Cải tiến trong sản xuất" ($\beta = 0.312$) và "Quản lý và phân công lao động" ($\beta = 0.264$) có tác động mạnh nhất.

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

Nghiên cứu mang lại những đóng góp kỹ thuật và phương pháp luận rõ nét:

  1. Lượng hóa mối tương quan đa chiều trong môi trường Lean Garment: Khác với nghiên cứu của Thomas P. Triebs & Subal C. Kumbhakar (2012) tại Ấn Độ vốn chỉ dừng lại ở phân tích định tính về sự đánh đổi giữa kỹ thuật và quản lý, công trình này cung cấp bộ trọng số cụ thể cho ngành nội y cao cấp.
  2. Chứng minh vai trò vượt trội của cải tiến vi mô (Kaizen/Gá dưỡng): Kết quả chứng minh cải tiến kỹ thuật thao tác ($\beta = 0.312$) mang lại hiệu quả gia tăng năng suất cao gấp 2,2 lần so với việc chỉ tác động đơn thuần vào điều kiện làm việc vật lý ($\beta = 0.142$).
  3. Mô hình hóa dữ liệu nhân khẩu học đặc thù: Bóc tách chi tiết cơ cấu 79,65% lao động nữ và 95,7% trình độ phổ thông, giải thích cơ chế tâm lý xã hội tác động trực tiếp đến tính ổn định của đường cong học tập (learning curve) trên chuyền may.
Tiêu chí Mô hình APO (2000) Mô hình Nguyễn Thị Phương Minh (2014) Mô hình Nghiên cứu Scavi Huế
Đối tượng Tổng quan SME dệt may Mã hàng NANO tại Dệt may 29-3 Dây chuyền sản xuất nội y Scavi Huế
Phương pháp Khung lý thuyết định tính Bấm giờ thao tác & Định mức Khảo sát Likert + EFA + Hồi quy OLS
Biến số quản trị Chung chung, cấp cao Cá nhân & Tiền lương Tách biệt quản lý trực tiếp & Môi trường
Tính ứng dụng Thấp cho chuyền may Hẹp cho một mã hàng Cao, áp dụng toàn bộ cụm nhà máy

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

Dựa trên trọng số hồi quy thực nghiệm, lộ trình triển khai nâng cao năng suất tại Công ty Scavi Huế được xây dựng qua 3 giai đoạn:

Đánh giá hiệu quả kinh tế (Cost-Benefit & ROI)

  • Chi phí đầu tư dự kiến: 850 triệu VNĐ (mua sắm bổ sung cữ gá bán tự động, đèn LED chống lóa chuyên dụng, phần mềm theo dõi sản lượng real-time trên 40 chuyền may Nhà máy 3).
  • Lợi ích kinh tế lượng hóa:
    • Tăng năng suất bình quân của 5.896 công nhân trực tiếp thêm $3.5%$.
    • Với doanh thu 2.583,51 tỷ VNĐ/năm, giá trị sản lượng gia tăng ước đạt: $2.583,51 \times 3.5% \approx 90,42$ tỷ VNĐ.
    • Tiết kiệm chi phí thời gian chờ đợi và sai hỏng vải: ước đạt 1,2 tỷ VNĐ/năm.
    • Thời gian hoàn vốn (Payback Period): Dưới 3 tháng kể từ khi vận hành đồng bộ.

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

  • Hạn chế mẫu và thời gian: Nghiên cứu khảo sát mẫu thuận tiện $N = 120$ tại Nhà máy 1 và Nhà máy 2; chưa bao quát toàn diện toàn bộ 6.310 công nhân và chưa phân tầng sâu theo từng dòng sản phẩm chuyên biệt (như Sport Bra phức tạp vs Quần lót Brief cơ bản).
  • Giới hạn biến số: Mô hình giải thích được $56.6%$ biến thiên của NSLĐ ($R^2_{\text{adj}} = 0.566$); còn $43.4%$ chịu ảnh hưởng bởi các yếu tố vĩ mô chưa được đưa vào mô hình như: gián đoạn chuỗi cung ứng nguyên phụ liệu quốc tế, biến động tỷ giá hối đoái, và chính sách làm thêm giờ.
  • Hướng nghiên cứu tiếp theo:
    1. Ứng dụng mô hình Cấu trúc Tuyến tính SEM (Structural Equation Modeling) nhằm phân tích các biến trung gian (Job Satisfaction, Work Stress).
    2. Xây dựng hệ thống IoT đo lường thời gian chu kỳ sản xuất (Cycle Time) tự động theo từng đường may để nạp dữ liệu huấn luyện các mô hình Machine Learning dự báo điểm nghẽn (Bottleneck Detection).

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

+-----------------------------------------------------------------------------------+
|                            GIÁ TRỊ MANG LẠI CHO CÁC BÊN                           |
+--------------------+--------------------------------------------------------------+
| Nhóm đối tượng     | Giá trị thực tiễn và định lượng                              |
+--------------------+--------------------------------------------------------------+
| Sinh viên & Học viên| Bộ khung mẫu chuẩn mực về phương pháp luận nghiên cứu định  |
| Kinh tế / QTKD     | lượng EFA - OLS ứng dụng trong môi trường công nghiệp thực tế.|
+--------------------+--------------------------------------------------------------+
| Quản trị viên & IE  | Bộ trọng số tác động ($\beta$) giúp ưu tiên đầu tư gá dưỡng,|
| (Kỹ sư công nghiệp)| thiết kế trạm thao tác chuẩn công thái học để tăng yield rate.|
+--------------------+--------------------------------------------------------------+
| Ban Lãnh đạo Scavi | Cơ sở dữ liệu khoa học để điều chỉnh kế hoạch sản xuất, cải  |
| và Doanh nghiệp May| tiến cơ chế khoán sản phẩm và đạt mục tiêu tăng trưởng 40%.   |
+--------------------+--------------------------------------------------------------+
| Giới nghiên cứu    | Khẳng định tính phù hợp của thang đo hành vi tổ chức trong   |
| Kinh tế ứng dụng   | bối cảnh lực lượng lao động phổ thông nữ tại Đông Nam Á.     |
+--------------------+--------------------------------------------------------------+

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

1. Yêu cầu phần cứng và phần mềm tối thiểu để chạy mô hình phân tích dữ liệu của đề tài là gì?

Hệ thống yêu cầu máy tính văn phòng tiêu chuẩn: CPU lõi kép (Intel Core i3 thế hệ 4 trở lên hoặc tương đương), RAM tối thiểu 4GB, ổ cứng trống 2GB. Về mặt phần mềm, cần cài đặt hệ điều hành Windows 7/10/11 hoặc Linux/macOS, tích hợp IBM SPSS Statistics version 20.0 (hoặc mới hơn) hoặc môi trường Python 3.8+ kèm các gói pandas, statsmodels, scipy.

2. Vì sao biến "Cải tiến trong sản xuất" lại có tác động lớn nhất đến năng suất lao động ($\beta = 0.312$)?

Trong ngành may mặc thời trang nội y (Lingerie), các chi tiết vải ren, voan, chun có tính co giãn cực lớn và đường may nhỏ, đòi hỏi độ chính xác tính bằng milimet. Việc áp dụng cữ gá dưỡng, máy may bán tự động chém mép và lập trình mẫu may giúp giảm phụ thuộc vào kỹ năng khéo léo cá nhân, giảm chu kỳ may (Cycle time) từ $15 - 25%$ cho mỗi công đoạn phức tạp.

3. Làm thế nào để tích hợp kết quả nghiên cứu vào hệ thống ERP / MES hiện tại của nhà máy?

Dữ liệu trọng số hồi quy có thể được sử dụng làm hệ số điều chỉnh (Weight coefficients) trong module cân bằng chuyền (Line Balancing Module) của hệ thống MES (Manufacturing Execution System). Khi hệ thống phân công lệnh sản xuất, thuật toán sẽ tự động tính toán năng suất kỳ vọng dựa trên mức độ tương thích tay nghề công nhân và mức độ hỗ trợ của thiết bị gá lắp.

4. Chi phí triển khai các giải pháp có vượt quá ngân sách hoạt động của doanh nghiệp vừa và nhỏ không?

Không. Nghiên cứu chỉ ra rằng các giải pháp quản trị phi tài chính (như phân công lao động theo tay nghề, cải thiện giao tiếp của tổ trưởng chuyền, sắp xếp ánh sáng chống chói) chiếm hơn $48%$ tổng mức độ tác động nhưng gần như không phát sinh chi phí đầu tư lớn. Các giải pháp kỹ thuật có thể triển khai dạng cuốn chiếu theo từng chuyền thử nghiệm (Pilot Line).

5. Khung thời gian hoàn vốn (ROI) và rủi ro chính khi áp dụng giải pháp là gì?

Thời gian hoàn vốn dự kiến dưới 3 tháng dựa trên mức tăng sản lượng $3.5%$. Rủi ro lớn nhất là tâm lý kháng cự thay đổi của công nhân lâu năm khi chuyển đổi thao tác may mới hoặc áp dụng bảng kiểm tra chất lượng tại chỗ (Inline AQL). Rủi ro này được triệt tiêu thông qua chính sách thưởng nóng Kaizen và đào tạo mẫu tại trạm thao tác chuẩn.

Kết luận

Đề tài "Nghiên cứu các yếu tố ảnh hưởng đến năng suất lao động của bộ phận sản xuất trực tiếp tại Công ty Scavi Huế" đã hoàn thành toàn diện các mục tiêu nghiên cứu đặt ra. Thông qua việc kết hợp chặt chẽ giữa khảo sát thực địa ($N = 120$), kiểm định thang đo Cronbach's Alpha, phân tích nhân tố khám phá EFA và xây dựng mô hình hồi quy tuyến tính đa biến OLS, nghiên cứu đã chứng minh rằng năng suất lao động chịu tác động thuận chiều của 5 nhân tố với thứ tự ưu tiên: Cải tiến trong sản xuất ($\beta = 0.312$) > Sự quản lý và phân công của cấp trên ($\beta = 0.264$) > Môi trường làm việc ($\beta = 0.218$) > Bản thân người lao động ($\beta = 0.185$) > Điều kiện làm việc ($\beta = 0.142$).

Kết quả nghiên cứu cung cấp luận cứ khoa học thực tiễn vững chắc cho Ban Giám đốc Công ty Scavi Huế và các doanh nghiệp dệt may trong việc tái cơ cấu dây chuyền, tối ưu hóa công thái học và nâng cao năng lực cạnh tranh quốc tế trong kỷ nguyên chuyển dịch chuỗi cung ứng toàn cầu.