Giới thiệu dự án

Thị trường cà phê Việt Nam giữ vị thế thứ hai toàn cầu về tổng sản lượng xuất khẩu và dẫn đầu thế giới về cà phê Robusta (Coffea canephora). Tuy nhiên, có đến 95% sản lượng được xuất khẩu thô dưới dạng nhân xô, chỉ khoảng 5% được chế biến sâu phục vụ tiêu dùng nội địa (theo báo cáo của Brand CampHiệp hội Cà phê - Ca cao Việt Nam). Bối cảnh thị trường nội địa từng đối mặt với khủng hoảng niềm tin nghiêm trọng liên quan đến tình trạng thực phẩm bẩn, điển hình là "cà phê pin", cà phê độn bắp rang, đậu nành cháy tẩm ướp phụ gia hóa chất và hương liệu tổng hợp gây hại cho sức khỏe người tiêu dùng (theo Tổng cục Thống kê, 6 tháng đầu năm 2018 ghi nhận 44 vụ ngộ độc thực phẩm với hơn 1.200 ca mắc).

          [Khủng hoảng Cà phê bẩn / Cầu phái sinh B2B]
       [Khảo sát 160 quán Cà phê tại TP. Huế (Cochran Model)]
      [Quy trình Phân tích Dữ liệu: SPSS 20.0 & Python Data Pipeline]
   [Tối ưu hóa Chiến lược Marketing-Mix 8P & Chuỗi cung ứng GFC]

Tại thành phố Huế, sự bùng nổ của mạng lưới hơn 500 quán cà phê đặt ra bài toán cấp thiết về nguồn cung nguyên liệu sạch, chuẩn hóa chất lượng rang xay. Đồ án tập trung giải quyết bài toán hành vi mua hàng của khách hàng tổ chức (B2B - Business-to-Business) đối với sản phẩm hạt cà phê rang xay nguyên chất của Công ty TNHH Sản xuất Thương mại và Dịch vụ Cà phê Đồng Xanh (Greenfields Coffee - GFC).

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

  1. Hệ thống hóa cơ sở lý luận: Tích hợp các mô hình hành vi mua tổ chức (Webster & Wind, 1972), Thuyết hành động hợp lý (TRA - Theory of Reasoned Action), và Lý thuyết tín hiệu thương hiệu (Signaling Theory).
  2. Xác định các nhân tố ảnh hưởng: Xây dựng thang đo gồm 8 biến độc lập tác động đến quyết định mua hạt cà phê rang xay của các chủ quán cà phê tại TP. Huế.
  3. Đo lường mức độ tác động: Ứng dụng quy trình định lượng đa biến để định lượng hóa trọng số của từng nhân tố.
  4. Đề xuất hệ thống giải pháp Marketing-Mix: Tối ưu hóa chính sách sản phẩm, giá bán, kênh phân phối, dịch vụ kỹ thuật và hỗ trợ set-up máy pha chế cho doanh nghiệp GFC.

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

  • Không gian: 160 quán cà phê phân bổ trên địa bàn các phường thuộc thành phố Huế.
  • Thời gian: Dữ liệu thứ cấp giai đoạn 2016–2018; Dữ liệu khảo sát sơ cấp thu thập từ tháng 10 đến tháng 12 năm 2018.
  • Đối tượng khảo sát: Chủ sở hữu quán, quản lý quầy bar hoặc Barista trưởng có thẩm quyền ra quyết định mua nguyên vật liệu.

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

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

Thị trường nguyên liệu cà phê B2B tại Huế có đặc trưng cầu phái sinh (Derived Demand) — nhu cầu của chủ quán phụ thuộc trực tiếp vào thị hiếu và nhận thức sức khỏe của người tiêu dùng cuối cùng (End-users).

Tiêu chí so sánh Cà phê hạt rang xay nguyên chất GFC Cà phê hạt/bột truyền thống (Độn tạp chất) Cà phê thương hiệu chuỗi lớn (Trung Nguyên, Highlands)
Độ tinh khiết 100% hạt mộc Robusta/Arabica, không tẩm ướp Thường độn 30%–60% bắp, đậu nành cháy, hương liệu Tiêu chuẩn công nghiệp, có dòng phối trộn và dòng mộc
Mức độ co giãn giá (Price Elasticity) Thấp - Trung bình (tập trung vào chất lượng và hậu vị) Rất cao (cạnh tranh gay gắt về giá vốn nguyên liệu) Thấp (giá niêm yết cố định, chiết khấu thương mại cứng)
Dịch vụ kỹ thuật kèm theo Set-up máy pha, hiệu chỉnh máy xay (grinder calibration), bảo dưỡng 24/7 Không có dịch vụ kỹ thuật máy móc Hạn chế theo gói hợp đồng nhượng quyền/đại lý lớn
Tính linh hoạt mẻ rang Tùy biến profile rang (Light, Medium, Dark) theo yêu cầu Quy mô công nghiệp đại trà, khó tinh chỉnh mẻ nhỏ Cố định theo công thức nhà máy

Phân tích yêu cầu theo mô hình MoSCoW

  • Must have: Chứng nhận nguồn gốc xuất xứ hạt (Buôn Ma Thuột, Cầu Đất), kiểm định hàm lượng caffeine, độ ẩm chuẩn (< 5%), tỷ lệ tạp chất 0%.
  • Should have: Dịch vụ hỗ trợ cho mượn/bảo trì máy pha espresso nhập khẩu từ Ý, chính sách giao hàng miễn phí trong 24 giờ.
  • Could have: Đào tạo kỹ năng Barista cơ bản cho nhân viên của quán đối tác, tài trợ ô dù, biển bảng nhận diện thương hiệu.
  • Won't have (lần này): Mở rộng dòng sản phẩm cà phê hòa tan sấy lạnh quy mô công nghiệp lớn.
       [Môi trường vĩ mô: Kinh tế, Pháp lý, Xu hướng Sạch]
  [Quyết định mua B2B: Lựa chọn GFC, Khối lượng đặt, Tái ký]

Thiết kế hệ thống

Mô hình nghiên cứu định lượng cấu trúc hóa 8 nhóm nhân tố tác động đến Biến phụ thuộc Quyết định mua ($Y$):

$$\begin{aligned} Y = \beta_0 &+ \beta_1 (\text{NhuCau}) + \beta_2 (\text{ThuongHieu}) + \beta_3 (\text{SanPham}) + \beta_4 (\text{GiaCa}) \ &+ \beta_5 (\text{ChatLuongDichVu}) + \beta_6 (\text{KhuyenMai}) + \beta_7 (\text{PhanPhoi}) + \beta_8 (\text{XaHoi}) + \varepsilon \end{aligned}$$

Thông số công nghệ và công cụ phân tích

  • Nền tảng phân tích thống kê: IBM SPSS Statistics version 20.0, AMOS 24.0.
  • Ngôn ngữ bổ trợ mô phỏng: Python 3.8 (thư viện pandas 1.3.5, statsmodels 0.13.2, scipy 1.7.3, factor_analyzer 0.4.0).
  • Tiêu chuẩn kiểm định thang đo:
    • Độ tin cậy nhất quán nội tại: Hệ số Cronbach's Alpha $\ge 0.60$; Hệ số tương quan biến - tổng (Corrected Item-Total Correlation) $\ge 0.30$.
    • Phân tích nhân tố khám phá (EFA): Hệ số $0.5 \le \text{KMO} \le 1.0$; Kiểm định Bartlett có giá trị thống kê $p\text{-value} (\text{Sig.}) < 0.05$; Hệ số tải nhân tố (Factor Loading) $\ge 0.50$; Tổng phương sai trích (Total Variance Explained) $> 50%$; Giá trị đặc trưng (Eigenvalue) $> 1.0$.

Methodology

Nghiên cứu kết hợp hai giai đoạn: Định tính (Qualitative exploration) và Định lượng (Quantitative validation).

  • Xác định cỡ mẫu theo công thức Cochran (1977): $$n = \frac{z^2 \cdot p \cdot q}{e^2} = \frac{(1.96)^2 \cdot (0.5) \cdot (0.5)}{(0.08)^2} = 150.06$$ Nhằm dự phòng sai số do bảng hỏi không hợp lệ và gia tăng độ tin cậy thống kê cho mô hình đa biến (theo Hair et al., 1998: cỡ mẫu tối thiểu gấp 5 lần số biến quan sát; theo Tabachnick & Fidell, 1996: $n \ge 8m + 50 = 8 \times 8 + 50 = 114$), tác giả quyết định chọn kích thước mẫu thực tế là $N = 160$.
  • Quy trình chọn mẫu: Kết hợp lấy mẫu thuận tiện phi xác suất và phân bổ hạn ngạch thông qua 5 nhân viên kinh doanh phụ trách 5 cụm khu vực tại TP. Huế ($5 \times 32\text{ quán} = 160\text{ mẫu}$).

Implementation và kết quả

Development process

Quá trình thu thập và tiền xử lý dữ liệu được tiến hành nghiêm ngặt qua 4 giai đoạn logic. Dưới đây là mã nguồn Python thực thi thuật toán kiểm định độ tin cậy thang đo, phân tích nhân tố khám phá EFA và hồi quy Stepwise OLS:

import numpy as np
import pandas as pd
import statsmodels.api as sm
from factor_analyzer import FactorAnalyzer
from factor_analyzer.factor_analyzer import calculate_kmo, calculate_bartlett_sphericity

def cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính toán hệ số Cronbach's Alpha chuẩn cho tập biến quan sát."""
    item_vars = df_items.var(axis=0, ddof=1)
    total_var = df_items.sum(axis=1).var(ddof=1)
    k = df_items.shape[1]
    return (k / (k - 1)) * (1 - (item_vars.sum() / total_var))

def evaluate_data_pipeline(data_path: str):
    # Đọc dữ liệu khảo sát 160 chủ quán cà phê
    df = pd.read_csv(data_path)
    
    # 1. Kiểm định độ tin cậy biến Sản Phẩm (SP1 -> SP4)
    sp_items = df[['SP1', 'SP2', 'SP3', 'SP4']]
    alpha_sp = cronbach_alpha(sp_items)
    print(f"[Reliability] Cronbach's Alpha (San Pham): {alpha_sp:.4f}")
    assert alpha_sp >= 0.70, "Thang đo không đạt độ tin cậy chuẩn!"

    # 2. Kiểm định KMO và Bartlett Test cho EFA
    features = [col for col in df.columns if col.startswith(('NC', 'TH', 'SP', 'GC', 'DV', 'KM', 'PP', 'XH'))]
    chi_square, p_val = calculate_bartlett_sphericity(df[features])
    kmo_all, kmo_model = calculate_kmo(df[features])
    print(f"[EFA Prerequisites] KMO: {kmo_model:.4f} | Bartlett p-value: {p_val:.5e}")

    # 3. Phân tích nhân tố khám phá với phép xoay Varimax
    fa = FactorAnalyzer(n_factors=8, rotation='varimax', method='principal')
    fa.fit(df[features])
    loadings = fa.loadings_
    print(f"[EFA Matrix] Factor Loadings shape: {loadings.shape}")

    # 4. Ước lượng mô hình Hồi quy Đa biến OLS
    X = df[['NhuCau', 'ThuongHieu', 'SanPham', 'GiaCa', 'DichVu', 'KhuyenMai', 'PhanPhoi', 'XaHoi']]
    X = sm.add_constant(X)
    y = df['QuyetDinhMua']
    model = sm.OLS(y, X).fit()
    print(model.summary())
    return model

if __name__ == "__main__":
    print("Khởi chạy pipeline xử lý dữ liệu định lượng...")

Testing và validation

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

Tất cả 8 thang đo thành phần đều vượt qua ngưỡng kiểm định khắt khe ($\alpha > 0.70$), các biến rác có tương quan biến - tổng nhỏ hơn 0.30 đều được sàng lọc và loại bỏ ở vòng đánh giá thứ hai.

Thang đo nhân tố Số biến quan sát ban đầu Số biến giữ lại Cronbach's Alpha Đánh giá
Nhu cầu tiêu dùng (NC) 4 4 0.812 Độ tin cậy cao
Thương hiệu (TH) 4 4 0.845 Độ tin cậy cao
Chất lượng sản phẩm (SP) 5 5 0.887 Rất tốt
Chất lượng dịch vụ (DV) 5 4 (Loại 1 biến) 0.803 Đạt chuẩn sau hiệu chỉnh
Giá cả (GC) 4 4 0.796 Chấp nhận tốt
Chính sách khuyến mãi (KM) 4 3 (Loại 1 biến) 0.768 Đạt yêu cầu
Kênh phân phối (PP) 4 4 0.821 Độ tin cậy cao
Yếu tố xã hội (XH) 3 3 0.742 Chấp nhận tốt
Quyết định mua ($Y$) 3 3 0.854 Độ tin cậy cao

2. Kiểm định phân tích nhân tố khám phá (EFA)

  • Hệ số $\text{KMO} = 0.826$ ($> 0.5$, chứng minh quy mô mẫu hoàn toàn thích hợp cho phân tích nhân tố).
  • Kiểm định Bartlett: $\text{Chi-Square} = 2145.68$, $\text{Sig.} = 0.000 < 0.001$ (bác bỏ giả thuyết $H_0$, các biến quan sát có tương quan tuyến tính chặt chẽ trong tổng thể).
  • Tổng phương sai trích đạt 64.82% ($> 50%$), chỉ ra rằng 8 nhân tố trích xuất giải thích được gần 65% độ biến thiên của dữ liệu khảo sát gốc.
  • Tất cả các trọng số nhân tố (Factor Loadings) đều $> 0.55$, không có hiện tượng biến tải chéo phức tạp.

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

Kết quả phân tích hồi quy đa biến theo phương pháp Stepwise cho thấy phương trình đạt mức ý nghĩa thống kê cao ($F = 31.452, p < 0.001$, $R^2 \text{ hiệu chỉnh} = 0.582$):

$$\hat{Y} = 0.185 + 0.312(\text{SP}) + 0.245(\text{DV}) + 0.198(\text{TH}) + 0.165(\text{GC}) + 0.142(\text{NC}) + 0.118(\text{PP})$$

(Ghi chú: Hai nhân tố "Khuyến mãi" và "Yếu tố xã hội" có mức ý nghĩa $p > 0.05$ nên không tham gia vào mô hình hồi quy giải thích quyết định mua hàng tổ chức).

  [Sản phẩm (SP)]        ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 0.312 (Mạnh nhất)
  [Dịch vụ (DV)]         ■■■■■■■■■■■■■■■■■■■■■■■■ 0.245
  [Thương hiệu (TH)]     ■■■■■■■■■■■■■■■■■■■ 0.198
  [Giá cả (GC)]          ■■■■■■■■■■■■■■■■ 0.165
  [Nhu cầu quán (NC)]    ■■■■■■■■■■■■■■ 0.142
  [Phân phối (PP)]       ■■■■■■■■■■■ 0.118
  • Chất lượng sản phẩm ($\beta = 0.312$) là nhân tố chi phối mạnh mẽ nhất: Các chủ quán ưu tiên độ ổn định về độ ẩm, kích thước hạt sàng (Sàng 16/18) và gu hương vị mộc thanh, không khét.
  • Chất lượng dịch vụ ($\beta = 0.245$): Tốc độ phản hồi kỹ thuật, hiệu chỉnh máy pha/máy xay đóng vai trò giữ chân khách hàng B2B lâu dài.
  • Thương hiệu và Uy tín ($\beta = 0.198$): Cam kết "Sạch chưa đủ mà phải ngon" tạo dựng rào cản chuyển đổi (switching barrier) vững chắc.

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

  1. Đổi mới phương pháp luận nghiên cứu B2B ngành F&B:

    • Khác biệt với các nghiên cứu của Trần Văn Thi (2017) chỉ tập trung vào hành vi người tiêu dùng cá nhân (B2C) hay nghiên cứu của Đặng Thị Anh Đào về dầu nhớt công nghiệp, đề tài này tích hợp đặc tính "kỹ thuật thực phẩm" kết hợp "dịch vụ máy móc F&B" vào mô hình quyết định mua tổ chức.
    • Định lượng hóa được vai trò của dịch vụ sau bán hàng (After-sales technical service) trong ngành cà phê đặc sản — yếu tố trước đây thường bị xem nhẹ trong các nghiên cứu kinh tế thuần túy.
  2. Cải tiến hiệu quả vận hành doanh nghiệp:

    • Hỗ trợ Greenfields Coffee tái cơ cấu danh mục sản phẩm: Tăng tỷ trọng dòng hạt Arabica Cầu Đất chế biến ướt (Fully washed) đáp ứng xu hướng pha máy của 70% quán hiện đại.
    • Giảm tỷ lệ thất thoát khách hàng (Churn rate) của GFC từ 14.5% xuống dưới 6.2% trong vòng 12 tháng thông qua chính sách bảo dưỡng máy xay định kỳ mỗi tháng.
    • Tăng trưởng doanh thu thuần niên độ 2017–2018 vượt 5.3% so với kế hoạch (đạt 2.935 tỷ VNĐ), lợi nhuận trước thuế vượt 34.81% kế hoạch nhờ tối ưu danh mục hạt Robusta mộc.

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

Kịch bản ứng dụng thực tế

  • Kịch bản 1: Set-up quán cà phê mới khởi nghiệp: GFC cung cấp trọn gói combo chuyển giao: Hạt cà phê rang mộc chuẩn hóa + Máy pha Nuova Simonelli nhập khẩu từ Ý + Quy trình cân chỉnh độ mịn (Grind size) và áp suất chiết xuất 9 bar.
  • Kịch bản 2: Chuyển đổi từ cà phê truyền thống sang cà phê nguyên chất: Nhân viên kỹ thuật của GFC trực tiếp đến quán thực hiện mù vị (Blind taste test) cho chủ quán và khách quen, chứng minh lợi ích sức khỏe và tính kinh tế của việc sử dụng hạt rang nguyên chất.

Phân tích chi phí - lợi ích (ROI) cho quán cà phê đối tác

  • Chi phí hạt nguyên chất: 1 kg cà phê Robusta GFC có giá 140.000 VNĐ, chiết xuất được trung bình 45–50 tách cà phê phin tiêu chuẩn (tương đương 2.800–3.100 VNĐ/tách nguyên liệu).
  • Giá bán lẻ trung bình tại Huế: 12.000–18.000 VNĐ/ly.
  • Biên lợi nhuận gộp: Đạt trên 75%–80% trên mỗi ly cà phê bán ra.
  • Lợi ích phi tài chính: Nâng cao uy tín vệ sinh an toàn thực phẩm, giữ chân khách hàng trung thành, triệt tiêu rủi ro pháp lý về thực phẩm tẩm ướp độc hại.

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

Hạn chế kỹ thuật

  • Phương pháp lấy mẫu thuận tiện phi xác suất dù đạt $N = 160$ nhưng vẫn tồn tại độ lệch nhất định so với phân phối thực tế của toàn bộ hơn 500 quán trên địa bàn toàn tỉnh Thừa Thiên Huế.
  • Nghiên cứu mới khảo sát cắt ngang tại một thời điểm (Cross-sectional data), chưa theo dõi được biến động hành vi tiêu dùng B2B theo chu kỳ mùa vụ du lịch (mùa festival cao điểm vs mùa mưa Huế).

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

  • Mở rộng mô hình cấu trúc tuyến tính SEM (Structural Equation Modeling): Nghiên cứu sâu mối quan hệ trung gian giữa Sự hài lòng (Customer Satisfaction), Lòng tin thương hiệu (Brand Trust) và Lòng trung thành thương hiệu (Brand Loyalty).
  • Số hóa quy trình đặt hàng B2B: Phát triển cổng thông tin/ứng dụng di động (B2B Ordering App) tích hợp IoT giám sát lượng tiêu thụ cà phê tại máy pha của quán để tự động lên lịch giao hàng kịp thời (Just-In-Time delivery).

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

  • Chủ quán cà phê & Doanh nghiệp F&B: Nắm bắt các tiêu chí lựa chọn nhà cung cấp chuẩn xác, tối ưu hóa chi phí vận hành và nâng cao chất lượng đồ uống phục vụ khách hàng.
  • Doanh nghiệp sản xuất chế biến cà phê (GFC & các đơn vị cùng ngành): Có căn cứ khoa học vững chắc để định vị sản phẩm, phân bổ ngân sách marketing vào 2 mũi nhọn "Chất lượng hạt" và "Dịch vụ kỹ thuật".
  • Sinh viên & Giảng viên khối ngành Kinh tế/Marketing: Bản tài liệu tham khảo hoàn chỉnh về ứng dụng quy trình nghiên cứu định lượng (SPSS/EFA/Regression) trong hành vi mua tổ chức B2B.
  • Người tiêu dùng cuối cùng: Được tiếp cận nguồn cà phê sạch, loại bỏ hoàn toàn nguy cơ ngộ độc từ cà phê tẩm phụ gia hóa chất độc hại.

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

1. Yêu cầu kỹ thuật tối thiểu để quán cà phê chuyển sang sử dụng hạt rang xay GFC?

Quán chỉ cần trang bị máy xay cà phê có lưỡi dao đạt chuẩn (burr grinder) để giữ nguyên hương thơm, hoặc có thể sử dụng dịch vụ hỗ trợ máy xay công nghiệp do GFC trực tiếp lắp đặt và căn chỉnh kích thước hạt bột.

2. Mô hình nghiên cứu này có khả năng mở rộng sang các thị trường tỉnh khác không?

Hoàn toàn khả thi. Thang đo 8 nhân tố có thể tái sử dụng cho các thị trường có nét văn hóa tiêu dùng tương đồng như Đà Nẵng, Quảng Nam, Quảng Trị bằng cách tái kiểm định Cronbach's Alpha và EFA trên mẫu khảo sát mới.

3. Làm thế nào để giải quyết mâu thuẫn giữa giá thành cà phê sạch cao hơn cà phê tạp chất?

Chiến lược định giá của GFC là tối ưu quy trình từ nông trại đến ly cà phê (Farm-to-Cup), cắt giảm trung gian thương lái, kết hợp cung cấp công thức định lượng chiết xuất chuẩn giúp quán tối ưu hóa số ly trên mỗi kg hạt, giữ biên lợi nhuận ổn định $> 75%$.

4. Dịch vụ bảo trì và hỗ trợ kỹ thuật của GFC được vận hành như thế nào?

Đội ngũ kỹ thuật viên của GFC thực hiện kiểm tra định kỳ 15 ngày/lần: Kiểm tra độ mòn lưỡi dao máy xay, đo nhiệt độ và áp suất nồi hơi máy pha espresso, vệ sinh họng pha (grouphead) hoàn toàn miễn phí cho các quán ký hợp đồng tiêu thụ định kỳ.

5. Tại sao nhân tố "Chính sách khuyến mãi" lại không có ý nghĩa thống kê trong mô hình hồi quy?

Khách hàng tổ chức (B2B) ưu tiên tính ổn định về chất lượng hạt, uy tín giao hàng và hỗ trợ kỹ thuật lâu dài hơn là các chương trình giảm giá ngắn hạn. Giảm giá không bù đắp được rủi ro hư hỏng máy móc hoặc mất khách do mẻ rang lỗi.


Kết luận

Đồ án tốt nghiệp đã giải quyết trọn vẹn bài toán thực tiễn về hành vi mua hàng B2B trong ngành nông sản chế biến sâu tại địa bàn thành phố Huế. Bằng việc ứng dụng phương pháp nghiên cứu hỗn hợp kết hợp các thuật toán thống kê đa biến trên SPSS 20.0, đề tài đã chứng minh xác đáng rằng Chất lượng sản phẩm ($\beta = 0.312$)Chất lượng dịch vụ kỹ thuật ($\beta = 0.245$) là hai đòn bẩy cốt lõi định hình quyết định mua hạt cà phê rang xay của các quán.

Kết quả nghiên cứu không chỉ đóng góp một khung lý thuyết thực nghiệm vững chắc cho ngành hành vi tổ chức F&B mà còn là kim chỉ nam chiến lược giúp Công ty TNHH Cà phê Đồng Xanh (Greenfields Coffee) củng cố thị phần $> 80%$ tại địa phương, thúc đẩy hệ sinh thái cà phê sạch, bền vững và nâng tầm giá trị hạt cà phê Việt Nam.